A new parallel overlapped domain decomposition method for nonlinear dynamic finite element analysis

In this paper a new parallel algorithm for nonlinear transient dynamic analysis of large structures has been presented. An unconditionally stable Newmark-beta method (constant average acceleration technique) has been employed for time integration. The proposed parallel algorithm has been devised within the broad framework of domain decomposition techniques. However, unlike most of the existing parallel algorithms (devised for structural dynamic applications) which are basically derived using nonoverlapped domains, the proposed algorithm uses overlapped domains. The parallel overlapped domain decomposition algorithm proposed in this paper has been formulated by splitting the mass, damping and stiffness matrices arises out of finite element discretisation of a given structure. A predictor-corrector scheme has been formulated for iteratively improving the solution in each step. A computer program based on the proposed algorithm has been developed and implemented with message passing interface as software development environment. PARAM-10000 MIMD parallel computer has been used to evaluate the performances. Numerical experiments have been conducted to validate as well as to evaluate the performance of the proposed parallel algorithm. Comparisons have been made with the conventional nonoverlapped domain decomposition algorithms. Numerical studies indicate that the proposed algorithm is superior in performance to the conventional domain decomposition algorithms. (C) 2003 Elsevier Ltd. All rights reserved.

Stability of the flow of a fluid through a flexible tube at high Reynolds number

The stability of the Hagen-Poiseuille flow of a Newtonian fluid in a tube of radius R surrounded by an incompressible viscoelastic medium of radius R < r < HR is analysed in the high Reynolds number regime. The dimensionless numbers that affect the fluid flow are the Reynolds number Re = (ρVR / η), the ratio of the viscosities of the wall and fluid ηr = (ηs/η), the ratio of radii H and the dimensionless velocity Γ = (ρV2/G)1/2. Here ρ is the density of the fluid, G is the coefficient of elasticity of the wall and Vis the maximum fluid velocity at the centre of the tube. In the high Reynolds number regime, an asymptotic expansion in the small parameter ε = (1/Re) is employed. In the leading approximation, the viscous effects are neglected and there is a balance between the inertial stresses in the fluid and the elastic stresses in the medium. There are multiple solutions for the leading-order growth rate do), all of which are imaginary, indicating that the fluctuations are neutrally stable, since there is no viscous dissipation of energy or transfer of energy from the mean flow to the fluctruations due to the Reynolds strees. There is an O(ε1/2) correction to the growth rate, s(1), due to the presence of a wall layer of thickness ε1/2R where the viscous stresses are O(ε1/2) smaller than the inertial stresses. An energy balance analysis indicates that the transfer of energy from the mean flow to the fluctuations due to the Reynolds stress in the wall layer is exactly cancelled by an opposite transfer of equal magnitude due to the deformation work done at the interface, and there is no net transfer from the mean flow to the fluctuations. Consequently, the fluctuations are stabilized by the viscous dissipation in the wall layer, and the real part of s(1) is negative. However, there are certain values of Γ and wavenumber k where s(l) = 0. At these points, the wail layer amplitude becomes zero because the tangential velocity boundary condition is identically satisfied by the inviscid flow solution. The real part of the O(ε) correction to the growth rate s(2) turns out to be negative at these points, indicating a small stabilizing effect due to the dissipation in the bulk of the fluid and the wall material. It is found that the minimum value of s(2) increases [is proportional to] (H − 1)−2 for (H − 1) [double less-than sign] 1 (thickness of wall much less than the tube radius), and decreases [is proportional to] (H−4 for H [dbl greater-than sign] 1. The damping rate for the inviscid modes is smaller than that for the viscous wall and centre modes in a rigid tube, which have been determined previously using a singular perturbation analysis. Therefore, these are the most unstable modes in the flow through a flexible tube.

Steady compressible flow of cohesionless granular materials through a wedge-shaped bunker

A continuum model based on the critical-state theory of soil mechanics is used to generate stress, density, and velocity profiles, and to compute discharge rates for the flow of granular material in a mass flow bunker. The bin–hopper transition region is idealized as a shock across which all the variables change discontinuously. Comparison with the work of Michalowski (1987) shows that his experimentally determined rupture layer lies between his prediction and that of the present theory. However, it resembles the former more closely. The conventional condition involving a traction-free surface at the hopper exit is abandoned in favour of an exit shock below which the material falls vertically with zero frictional stress. The basic equations, which are not classifiable under any of the standard types, require excessive computational time. This problem is alleviated by the introduction of the Mohr–Coulomb approximation (MCA). The stress, density, and velocity profiles obtained by integration of the MCA converge to asymptotic fields on moving down the hopper. Expressions for these fields are derived by a perturbation method. Computational difficulties are encountered for bunkers with wall angles θw [gt-or-equal, slanted] 15° these are overcome by altering the initial conditions. Predicted discharge rates lie significantly below the measured values of Nguyen et al. (1980), ranging from 38% at θw = 15° to 59% at θw = 32°. The poor prediction appears to be largely due to the exit condition used here. Paradoxically, incompressible discharge rates lie closer to the measured values. An approximate semi-analytical expression for the discharge rate is obtained, which predicts values within 9% of the exact (numerical) ones in the compressible case, and 11% in the incompressible case. The approximate analysis also suggests that inclusion of density variation decreases the discharge rate. This is borne out by the exact (numerical) results – for the parameter values investigated, the compressible discharge rate is about 10% lower than the incompressible value. A preliminary comparison of the predicted density profiles with the measurements of Fickie et al. (1989) shows that the material within the hopper dilates more strongly than predicted. Surprisingly, just below the exit slot, there is good agreement between theory and experiment.

Stress induced phase transition in a monomolecular perfluroalkylsilane film self assembled on aluminium surface

We control the stiffnesses of two dual double cantelevers placed in series to control penetration into a perflurooctyltrichlorosilane monolayer self assembled on aluminium and silicon substrates. The top cantilever which carries the probe is displaced with respect to the bottom cantilever which carries the substrate, the difference in displacement recorded using capacitors gives penetration. We further modulate the input displacement sinusoidally to deconvolute the viscoelastic properties of the monolayer. When the intervention is limited to the terminal end of the molecule there is a strong viscous response in consonance with the ability of the molecule to dissipate energy by the generation of gauche defects freely. When the intervention reaches the backbone, at a contact mean pressure of 0.2GPa the damping disappears abruptly and the molecule registers a steep rise in elastic modulus and relaxation time constant, with increasing contact pressure. We offer a physical explanation of the process and describe this change as due to a phase transition from a liquid like to a solid like state.

Axisymmetric turbulent boundary layer along a circular cylinder at constant pressure

A constant-pressure axisymmetric turbulent boundary layer along a circular cylinder of radius a is studied at large values of the frictional Reynolds number a+ (based upon a) with the boundary-layer thickness δ of order a. Using the equations of mean motion and the method of matched asymptotic expansions, it is shown that the flow can be described by the same two limit processes (inner and outer) as are used in two-dimensional flow. The condition that the two expansions match requires the existence, at the lowest order, of a log region in the usual two-dimensional co-ordinates (u+, y+). Examination of available experimental data shows that substantial log regions do in fact exist but that the intercept is possibly not a universal constant. Similarly, the solution in the outer layer leads to a defect law of the same form as in two-dimensional flow; experiment shows that the intercept in the defect law depends on δ/a. It is concluded that, except in those extreme situations where a+ is small (in which case the boundary layer may not anyway be in a fully developed turbulent state), the simplest analysis of axisymmetric flow will be to use the two-dimensional laws with parameters that now depend on a+ or δ/a as appropriate.

Capturing higher modes of vibration of micromachined resonators

MEMS resonators have potential applications in the areas of RF-MEMS, clock oscillators, ultrasound transducers, etc. The important characteristics of a resonator are its resonant frequency and Q-factor (a measure of damping). Usually large damping in macro structures makes it difficult to excite and measure their higher modes. In contrast, MEMS resonators seem amenable to excitation in higher modes. In this paper, 28 modes of vibration of an electrothermal actuator are experimentally captured–perhaps the highest number of modes experimentally captured so far. We verify these modes with FEM simulations and report that all the measured frequencies are within 5% of theoretically predicted values.

The Effect of Evacuated Backside Cavity on the Dynamic Characteristics of a Capacitive Micromachined Ultrasonic Transducer

In this paper, we investigate the effect of vacuum sealing the backside cavity of a Capacitive Micromachined Ultrasonic Transducer (CMUT). The presence or absence of air inside the cavity has a marked effect upon the system parameters, such as the natural frequency, damping, and the pull-in voltage. The presence of vacuum inside the cavity of the device causes a reduction in the effective gap height which leads to a reduction in the pull-in voltage. We carry out ANSYS simulations to quantify this reduction. The presence of vacuum inside the cavity of the device causes stress stiffening of the membrane, which changes the natural frequency of the device. A prestressed modal analysis is carried out to determine the change in natural frequency due to stress stiffening. The equivalent circuit method is used to evaluate the performance of the device in the receiver mode. The lumped parameters of the device are obtained and an equivalent circuit model of the device is constructed to determine the open circuit receiving sensitivity of the device. The effect of air in the cavity is included by incorporating an equivalent compliance and an equivalent resistance in the equivalent circuit.

Analytical Evaluation of Squeeze Film Forces in a CMUT With Sealed Air-Filled Cavity

The presence of vacuum inside the cavity of a capacitive micromachined ultrasonic transducer (CMUT) causes the membrane of the device (which is the main vibrating structural component) to deflect towards the substrate, thereby causing a reduction in the effective gap height. This reduction causes a drastic decrease in the pull-in voltage of the device limiting the DC bias at which the device can be operated for maximum efficiency. In addition, this initial deflection of the membrane due to atmospheric pressure, causes significant stress stiffening of the the membrane, changing the natural frequency of the device significantly from the design value. To circumvent the deleterious effects of vacuum in the sealed cavity, we investigate the possibility of using sealed CMUT cavities with air inside at ambient pressure. In order to estimate the transducer loss due to the presence of air in the sealed cavity, we evaluate the resulting damping and determine the forces acting on the vibrating membrane resulting from the compression of the trapped air film. We take into account the flexure of the top vibrating membrane instead of assuming the motion to be parallel-plate like. Towards this end, we solve the linearized Reynolds equation using the appropriate boundary conditions and show that, for a sealed CMUT cavity, the presence of air does not cause any squeeze film damping.

A Nonlinear Voltage Regulator With One Tunable Parameter for Multimachine Power Systems

This paper proposes a nonlinear voltage regulator with one tunable parameter for multimachine power systems. Based on output feedback linearization, this regulator can achieve simultaneous voltage regulation and small-signal performance objectives. Conventionally output feedback linearization has been used for voltage regulator design by taking infinite bus voltage as reference. Unfortunately, this controller has poor small-signal performance and cannot be applied to multimachine systems without the estimation of the equivalent external reactance seen from the generator. This paper proposes a voltage regulator design by redefining the rotor angle at each generator with respect to the secondary voltage of the step-up transformer as reference instead of a common synchronously rotating reference frame. Using synchronizing and damping torques analysis, we show that the proposed voltage regulator achieves simultaneous voltage regulation and damping performance over a range of system and operating conditions by controlling the relative angle between the generator internal voltage angle delta and the secondary voltage of the step up transformer. The performance of the proposed voltage regulator is evaluated on a single machine infinite bus system and two widely used multimachine test systems.

Viscoelastic phenomenology based structure assignment for closed-loop vibration control of a beam with sensors and actuators

In this paper we incorporate a novel approach to synthesize a class of closed-loop feedback control, based on the variational structure assignment. Properties of a viscoelastic system are used to design an active feedback controller for an undamped structural system with distributed sensor, actuator and controller. Wave dispersion properties of onedimensional beam system have been studied. Efficiency of the chosen viscoelastic model in enhancing damping and stability properties of one-dimensional viscoelastic bar have been analyzed. The variational structure is projected on a solution space of a closed-loop system involving a weakly damped structure with distributed sensor and actuator with controller. These assign the phenomenology based internal strain rate damping parameter of a viscoelastic system to the usual elastic structure but with active control. In the formulation a model of cantilever beam with non-collocated actuator and sensor has been considered. The formulation leads to the matrix identification problem of two dynamic stiffness matrices. The method has been simplified to obtain control system gains for the free vibration control of a cantilever beam system with collocated actuator-sensor, using quadratic optimal control and pole-placement methods.

Embedded Foundations Under Vertical Vibration

In this note, a simplified procedure based on energy consideration, has been developed, for the solution of steady-state vibration of a system with combined viscous and Coulomb friction damping, subjected to frequency in dependent and frequency dependent excitation, which yields results essentially same as the exact solution. The proposed method uses results essentially same as the exact solution. The proposed method uses equivalent damping which assumes that if the damping in a system is small, the total damping effect can be represented by that of an equivalent damper.

Relationship Between Physical Structure and Tribology of Single Soot Particles Generated by Burning Ethylene

Ethylene gas is burnt and the soot generated is sampled thermophoretically at different heights along the flame axis starting from a region close to the root of the flame. The morphology and crystallinity of the particle are recorded using high resolution transmission electron microscopes. The hardness of a single particle is measured using a nanoindenter. The frictional resistance and material removal of a particle are measured using an atomic force microscope. The particles present in the mid-flame region are found to have a crystalline shell. The ones at the flame root are found to be highly disordered and the ones at the flame tip and above have randomly distributed pockets of short range order. The physical state of a particle is found to relate, but not very strongly, with the mechanical and tribological properties of the particles.

Steady-State Stability of Current-Mode Active-Clamp ZVS DC-DC Converters

Active-clamp dc-dc converters are pulsewidth-modulated converters having two switches featuring zero-voltage switching at frequencies beyond 100 kHz. Generalized equivalent circuits valid for steady-state and dynamic performance have been proposed for the family of active-clamp converters. The active-clamp converter is analyzed for its dynamic behavior under current control in this paper. The steady-state stability analysis is presented. On account of the lossless damping inherent in the active-clamp converters, it appears that the stability region in the current-controlled active-clamp converters get extended for duty ratios, a little greater than 0.5 unlike in conventional hard-switched converters. The conventional graphical approach fails to assess the stability of current-controlled active-clamp converters, due to the coupling between the filter inductor current and resonant inductor current. An analysis that takes into account the presence of the resonant elements is presented to establish the condition for stability. This method correctly predicts the stability of the current-controlled active-clamp converters. A simple expression for the maximum duty cycle for subharmonic-free operation is obtained. The results are verified experimentally.

Localization induced base isolation in fractionally damped non linear system

In the present article we take up the study of nonlinear localization induced base isolation of a 3 degree of freedom system having cubic nonlinearities under sinusoidal base excitation. The damping forces in the system are described by functions of fractional derivative of the instantaneous displacements, typically linear and quadratic damping are considered here separately. Under the assumption of smallness of certain system parameters and nonlinear terms an approximate estimate of the response at each degree of freedom of the system is obtained by the Method of Multiple Scales approach. We then consider a similar system where the nonlinear terms and certain other parameters are no longer small. Direct numerical simulation is made use of to obtain the amplitude plot in the frequency domain for this case, which helps us to establish the efficacy of this method of base isolation for a broad class of systems. Base isolation obtained this way has no counterpart in the linear theory.

Response of Materials During Sliding on Various Surface Textures

In the present investigation, soft materials, such as Al-4Mg alloy, high-purity Al and pure Mg pins were slid against hard steel plates of various surface textures to study the response of materials during sliding. The experiments were conducted using an inclined pin-on-plate sliding apparatus under both dry and lubricated conditions in an ambient environment. Two kinds of frictional response, namely steady-state and stick-slip, were observed during sliding. In general, the response was dependent on material pair, normal load, lubrication, and surface texture of the harder material. More specifically, for the case of Al-4Mg alloy, the stick-slip response was absent under both dry and lubricated conditions. For Al, stick-slip was observed only under lubricated conditions. For the case of Mg, the stick-slip response was seen under both dry and lubricated conditions. Further, it was observed that the amplitude of stick-slip motion primarily depends on the plowing component of friction. The plowing component of friction was the highest for the surfaces that promoted plane strain conditions and was the lowest for the surfaces that promoted plane stress conditions near the surface.