Natural frequencies of submerged piezoceramic hollow spheres

An exact 3D analysis of free vibration of a piezoceramic hollow sphere submerged in a compressible fluid is presented in this paper. A separation method is adopted to simplify the basic equations for spherically isotropic piezoelasticity. It is shown that there are two independent classes of vibration. The first one is independent of the fluid medium as well as the electric field, while the second is associated with both the fluid parameter and the piezoelectric effect. Exact frequency equations are derived and numerical results are obtained. The project supported by the National Natural Science Foundation of China (No. 19872060)     

Global analysis of secondary bifurcation of an elastic bar

In a three dimensional framework of finite deformation configurations, this paper investigates the secondary bifurcation of a uniform, isotropic and linearly elastic bar under compression in a large range of parameters. The governing differential equations and finite dimensional equations of this problem are discussed. It is found that, for a bar with two ends hinged, usually many secondary bifurcation points appear on the primary branches which correspond to the maximum bending stiffness. Results are shown on parameter charts. Secondary modes and branches are also calculated with numerical methods. The project supported in part by the National Natural Science Foundation of China     

Exact static analysis of a rotating piezoelectric spherical shell

This paper gives an exact analysis of a rotating piezoelectric spherical shell with arbitrary thickness. At first, three displacement functions are introduced to simplify the basic equations of a spherically isotropic, piezoelectric medium. By expanding the displacement functions as well as the electric potential in terms of spherical harmonics, the basic equations of equilibrium are converted to an uncoupled Euler type, second order ordinary differential equation and a coupled system of three second order ordinary differential equations. A general solution to the homogeneous equations of equilibrium is then derived. The static analysis of a rotating spherical shell is then performed and numerical example is presented. The project supported by the National Natural Science Foundation of China, the Zhejiang Provincial Natural Science Foundation, and the Japanese Committee of Culture, Education and Science.     

Effect of microstructure on the hardening and softening behaviors of polycrystalline shape memory alloys Part II: Numerical simulation under axisymmetrical loading

Based on the microstructure-based constitutive model established in Part I, a detailed numerical investigation on the role of each microstructure parameter in the kinematical and kinetic evolution of polycrystalline SMA under axisymmetrical tension loading is performed. Some macroscopic constitutive features of stress-induced martensite transformation are discussed. The subject supported by the Research Grant Committee (RGC) of Hong Kong SAR, the National Natural Science Foundation of China and the Provincial Natural Science Foundation of Jiangxi Province of China     

Analysis of interlaminar stress and nonlinear dynamic response for composite laminated plates with interfacial damage

By considering the effect of interfacial damage and using the variation principle, three-dimensional nonlinear dynamic governing equations of the laminated plates with interfacial damage are derived based on the general sixdegrees-of-freedom plate theory towards the accurate stress analysis. The solutions of interlaminar stress and nonlinear dynamic response for a simply supported laminated plate with interfacial damage are obtained by using the finite difference method, and the results are validated by comparison with the solution of nonlinear finite element method. In numerical calculations, the effects of interfacial damage on the stress in the interface and the nonlinear dynamic response of laminated plates are discussed.     

A combined parametric quadratic programming and precise integration method based dynamic analysis of elastic-plastic hardening/softening problems

The objective of the paper is to develop a new algorithm for numerical solution of dynamic elastic-plastic strain hardening/softening problems. The gradient dependent model is adopted in the numerical model to overcome the result mesh-sensitivity problem in the dynamic strain softening or strain localization analysis. The equations for the dynamic elastic-plastic problems are derived in terms of the parametric variational principle, which is valid for associated, non-associated and strain softening plastic constitutive models in the finite element analysis. The precise integration method, which has been widely used for discretization in time domain of the linear problems, is introduced for the solution of dynamic nonlinear equations. The new algorithm proposed is based on the combination of the parametric quadratic programming method and the precise integration method and has all the advantages in both of the algorithms. Results of numerical examples demonstrate not only the validity, but also the advantages of the algorithm proposed for the numerical solution of nonlinear dynamic problems. The project supported by the National Key Basic Research Special Foundation (G1999032805), the National Natural Science Foundation of China (19872016, 50178016, 19832010) and the Foundation for University Key Teacher by the Ministry of Education of China     

The solution for the generalized riccati algebraic equations of linear equality constraint system

Based on the dynamic equation, the performance functional and the system constraint equation of time-invariant discrete LQ control problem, the generalized Riccati equations of linear equality constraint system are obtained according to the minimum principle, then a deep discussion about the above equations is given, and finally numerical example is shown in this paper. Project supported by the National Natural Science Youth Foundation of China     

The three-point fifth-order accurate generalized compact scheme and its applications

A three-point fifth-order accurate generalized compact scheme (GC scheme) with a spectral-like resolution is constructed in a general way. The scheme satisfies the principle of stability and the principle about suppression of the oscillations, therefore numerical errors can decay automatically and no spurious oscillations are generated around shocks. The third-order TVD type Runge-Kutta method is employed for the time integration, thus making the GC scheme best suited for unsteady problems. Numerical results show that the GC scheme is shock-capturing. The time-dependent boundary conditions proposed by Thompson are well employed when the algorithm is applied to the Euler equations of gas dynamics. The project supported by the National Natural Science Foundation of China (19972038) and Foundation of the National CFD Laboratory of China     

GURTIN VARIATIONAL PRINCIPLE AND FINITE ELEMENT SIMULATION FOR DYNAMICAL PROBLEMS OF FLUID-SATURATED POROUS MEDIA

Based on the theory of porous media,a general Gurtin variational principle for theinitial boundary value problem of dynamical response of fluid-saturated elastic porous media isdeveloped by assuming infinitesimal deformation and incompressible constituents of the solid andfluid phase.The finite element formulation based on this variational principle is also derived.Asthe functional of the variational principle is a spatial integral of the convolution formulation,thegeneral finite element discretization in space results in symmetrical differential-integral equationsin the time domain.In some situations,the differential-integral equations can be reduced to sym-metrical differential equations and,as a numerical example,it is employed to analyze the reflectionof one-dimensional longitudinal wave in a fluid-saturated porous solid.The numerical results canprovide further understanding of the wave propagation in porous media.     

Dynamical Behavior of Viscoelastic Cylindrical Shells Under Axial Pressures

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｏｒｙｏｆｖｉｓｃｏｅｌａｓｔｉｃｉｔｙａｎｄｉｔｓａｐｐｌｉｃａｔｉｏｎｓｈａｖｅｂｅｅｎｂｅｃｏｍｉｎｇａｖｅｒｙａｔｔｒａｃｔｉｖｅｒｅｓｅａｒｃｈｆｉｅｌｄｓｉｎｃｅｖｉｓｃｏｅｌａｓｔｉｃｓｔｒｕｃｔｕｒｅｓｗｅｒｅｗｉｄｅｌｙａｐｐｌｉｅｄｉｎｔｏｖａｒｉｏｕｓｅｎｇｉｎｅｅｒｉｎｇ .Ｂｕｔｔｏｔｈｅａｕｔｈｏｒｓ’ｋｎｏｗｌｅｄｇｅ,ｔｈｅｒｅａｒｅｒａｒｅｆｏｒｐａｐｅｒｓｃｏｎｃｅｒｎｉｎｇｗｉｔｈｄｙｎａｍｉｃａｌｂｅｈａｖｉｏｒｏｆｖｉｓｃｏｅ…     

Chaotic attitude and reorientation maneuver for completely liquid-filled spacecraft with flexible appendage

The present paper investigates the chaotic attitude dynamics and reorientation maneuver for completely viscous liquid-filled spacecraft with flexible appendage. All of the equations of motion are derived by using Lagrangian mechanics and then transformed into a form consisting of an unperturbed part plus perturbed terms so that the system＇s nonlinear characteristics can be exploited in phase space. Emphases are laid on the chaotic attitude dynamics produced from certain sets of physical parameter values of the spacecraft when energy dissipation acts to derive the body from minor to major axis spin. Numerical solutions of these equations show that the attitude dynamics of liquid-filled flexible spacecraft possesses characteristics common to random, non- periodic solutions and chaos, and it is demonstrated that the desired reorientation maneuver is guaranteed by using a pair of thruster impulses. The control strategy for reorientation maneuver is designed and the numerical simulation results are presented for both the uncontrolled and controlled spins transition.     

Perturbation solution to 3-D nonlinear supercavitating flow problems

A 3-D nonlinear problem of supercavitating flow past an axisymmetric body at a small angle of attack is investigated by means of the perturbation method and Fourier-cosine-expansion method. The first three order perturbation equations are derived in detail and solved numerically using the boundary integral equation method and iterative techniques. Computational results of the hydrodynamic characteristics and cavity shapes of each order are presented for nonaxisymmetric supercavitating flow past cones with various apex-angles at different cavitation numbers. The numerical results are found in good agreement with experimental data. The project supported by the National Natural Science Foundation of China     

Interval parameter perturbation method for evaluating the bounds on natural frequencies of structures with interval parameters

For structural parameters with uncertainties, interval mathematics can, in the case where the probabilistic distribution density of uncertain variables is unavailable, deal with the influence of uncertainties in structural parameters on the response of structures. In order to evaluate the region containing natural frequencies of structures with interval parameters, the interval parameter perturbation method is presented in this paper. The advantage of the present method is its computational efficiency in evaluating the region containing natural frequencies. A numerical example is used to illustrate the efficiency of the method proposed. The project is supported by National Youth Natural Science Foundation of China and National Outstanding Youth Science Foundation of China.     

METHOD OF GREEN＇S FUNCTION OF CORRUGATED SHELLS

By using the fundamental equations of axisymmetric shallow shells of revolution, the nonlinear bending of a shallow corrugated shell with taper under arbitrary load has been investigated. The nonlinear boundary value problem of the corrugated shell was reduced to the nonlinear integral equations by using the method of Green＇s function. To solve the integral equations, expansion method was used to obtain Green＇s function. Then the integral equations were reduced to the form with degenerate core by expanding Green＇s function as series of characteristic function. Therefore, the integral equations become nonlinear algebraic equations. Newton＇ s iterative method was utilized to solve the nonlinear algebraic equations. To guarantee the convergence of the iterative method, deflection at center was taken as control parameter. Corresponding loads were obtained by increasing deflection one by one. As a numerical example,elastic characteristic of shallow corrugated shells with spherical taper was studied.Calculation results show that characteristic of corrugated shells changes remarkably. The snapping instability which is analogous to shallow spherical shells occurs with increasing load if the taper is relatively large. The solution is close to the experimental results.     

Dynamics of flexible multibody systems with tree topologies

The dynamic equations of flexible multibody systems with tree topological configuration are derived by using the Jourdain's principle. The independent joint coordinates are introduced to describe the large displacements of the bodies, and the modal coordinates are used to describe small deformations of flexible bodies based on the consistent mass finite element method and normal vibration mode analysis. The minimum differential equations are developed, which are compatible with the equations of multi-rigid body systems or structural dynamics. The stiff problem in the numerical integration is thus alleviated effectively. The method used in this paper can be extended to deal with systems with other topological configurations. Finally, the validity and feasibility of the presented mathematical model are demonstrated by a numerical example of a manipulator with two elastic links. The project supported by National Natural Science Foundation of China     

Numerical simulation of supersonic reacting mixing layer

In this paper, the supersonic chemically reacting mixing layer is simulated with the third order ENN scheme, based on the Navier-Stokes equations, containing transport equations of all species. The numerical results show that the thickness of mixing layer increases gradually along the flow direction, and that the Kelvin-Helmholtz instabilities may not exist in mixing layer flows. The project supported by the National Natural Science Foundation of China     

On the predictability of chaotic systems with respect to maximally effective computation time

The round-off error introduces uncertainty in the numerical solution. A computational uncertainty principle is explained and validated by using chaotic systems, such as the climatic model, the Rossler and super chaos system. Maximally effective computation time (MECT) and optimal stepsize (OS) are discussed and obtained via an optimal searching method. Under OS in solving nonlinear ordinary differential equations, the self-memorization equations of chaotic systems are set up, thus a new approach to numerical weather forecast is described. The project supported by the National Natural Science Foundation of China (40275031 and 40231006), the National Key Program for Developing Basic Sciences (G1999043408) and the Key Innovation Project of Chinese Academy of Sciences (K2CX1-10-07)     

An improved large eddy simulation of two-phase flows in a pump impeller

An improved large eddy simulation using a dynamic second-order sub-grid-scale (SGS) stress model has been developed to model the governing equations of dense turbulent particle-liquid two-phase flows in a rotating coordinate system, and continuity is conserved by a mass-weighted method to solve the filtered governing equations. In the current second-order SGS model, the SGS stress is a function of both the resolved strain-rate and rotation-rate tensors, and the model parameters are obtained from the dimensional consistency and the invariants of the strain-rate and the rotation-rate tensors. In the numerical calculation, the finite volume method is used to discretize the governing equations with a staggered grid system. The SIMPLEC algorithm is applied for the solution of the discretized governing equations. Body-fitted coordinates are used to simulate the two-phase flows in complex geometries. Finally the second-order dynamic SGS model is successfully applied to simulate the dense turbulent particle-liquid two-phase flows in a centrifugal impeller. The predicted pressure and velocity distributions are in good agreement with experimental results. The project supported by the National Natural Science Foundation of China (50779069 and 90510007), the Start-up Scientific Research Foundation of China Agricultural University (2006021) and the Beijing Natural Science Foundation (3071002).     

A computational method on process optimal control for finite deformation problems

On the principle of non-incremental algorithm, some basic ideas of process optimal control iterative algorithm, based on the Optimal Control Variational Principle in Mechanics, is proposed in this paper. Then the essential governing equations are presented. This work provides a new method to achieve the numerical solutions of the mechanic of finite deformation. Project supported by the National Natural Science Foundation of China (Grant No. 594305050).     

Localized analysis of thin-walled structure's buckling/initial post-buckling and its accuracy

A method of localization is proposed to lower the high order of equations in FEM calculation for the stability of a complex thin-walled structure. The localized analysis enables us to obtain both the upper and lower limits for the bifurcating point in a whole linear-elastic structural system, as well as an approximate solution to asymptotic post-buckling problem. Some numerical examples are included. Project supported by National Natural Science Foundation of China