UNSTEADY／STEADY NUMERICAL SIMULATION OF THREE-DIMENSIONAL INCOMPRESSIBLE NAVIER-STOKES EQUATIONS ON ARTIFICIAL COMPRESSIBILITY

A mixed algorithm of central and upwind difference scheme for the solution of steady/unsteady incompressible Navier-Stokes equations is presented. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The numerical flux of semi-discrete equations is computed by using the Roe approximation. Time accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The algebraic turbulence model of Baldwin-Lomax is ulsed in this work. As examples, the solutions of flow through two dimensional flat, airfoil, prolate spheroid and cerebral aneurysm are computed and the results are compared with experimental data. The results show that the coefficient of pressure and skin friction are agreement with experimental data, the largest discrepancy occur in the separation region where the lagebraic turbulence model of Baldwin-Lomax could not exactly predict the flow.     

AN EXACT ANALYSIS OF THERMAL BUCKLING OF PIEZOELECTRIC LAMINATED PLATES

In this paper,the governing differential equations of elastic stability problems in ther-mopiezoelectric media are deduced.The solutions of the thermal buckling problems for piezoelectriclaminated plates are presented in the context of the mathematical theory of elasticity.Owing to thecomplexity of the eigenvalue problem involved,the critical temperature values of thermal bucklingmust be solved numerically.The numerical results for piezoelectric/non-piezoelectric laminated platesare presented and the influence of piezoelectricity upon thermal buckling temperature is discussed.     

ANALYTICAL RELATIONS BETWEEN EIGENVALUES OF CIRCULAR PLATE BASED ON VARIOUS PLATE THEORIES

Based on the mathematical similarity of the axisymmetric eigenvalue problems of a circular plate between the classical plate theory(CPT), the first-order shear deformation plate theory(FPT) and the Reddy's third-order shear deformation plate theory (RPT), analytical relations between the eigenvalues of circular plate based on various plate theories are investigated. In the present paper, the eigenvalue problem is transformed to solve an algebra equation. Analytical relationships that are expressed explicitly between various theories are presented. Therefore, from these relationships one can easily obtain the exact RPT and FPT solutions of critical buckling load and natural frequency for a circular plate with CPT solutions. The relationships are useful for engineering application, and can be used to check the validity, convergence and accuracy of numerical results for the eigenvalue problem of plates.     

EIGENVALUE PROBLEM OF A LARGE SCALE INDEFINITE GYROSCOPIC DYNAMIC SYSTEM

Gyroscopic dynamic system can be introduced to Hamiltonian system.Based on an adjoint symplectic subspace iteration method of Hamiltonian gyroscopic system, an adjoint symplectic subspace iteration method of indefinite Hamiltonian function gy- roscopic system was proposed to solve the eigenvalue problem of indefinite Hamiltonian function gyroscopic system.The character that the eigenvalues of Hamiltonian gyroscopic system are only pure imaginary or zero was used.The eigenvalues that Hamiltonian function is negative can be separated so that the eigenvalue problem of positive definite Hamiltonian function system was presented,and an adjoint symplectic subspace iteration method of positive definite Hamiltonian function system was used to solve the separated eigenvalue problem.Therefore,the eigenvalue problem of indefinite Hamiltonian function gyroscopic system was solved,and two numerical examples were given to demonstrate that the eigensolutions converge exactly.     

AN INVESTIGATION ON A NEW INTEGRAL TRANSFORM-BOUNDARY ELEMENT METHOD FOR ELASTODYNAMICS

The numerical methods of Fourier eigen transform FET and its inversion are discussedand applied to the boundary element method for elastodynamics. The program for solving elastody-namic problems with the boundary element method is developed and some examples are given. Fromthe numerical results of the examples, we know the method can increase the computing speed 5～10times and the accuracy is guaranteed.     

An adaptive reanalysis method for genetic algorithm with application to fast truss optimization

Although the genetic algorithm （GA） for structural optimization is very robust, it is very computationally intensive and hence slower than optimality criteria and mathematical programming methods. To speed up the design process, the authors present an adaptive reanalysis method for GA and its applications in the optimal design of trusses. This reanalysis technique is primarily derived from the Kirsch＇s combined approximations method. An iteration scheme is adopted to adaptively determine the number of basis vectors at every generation. In order to illustrate this method, three classical examples of optimal truss design are used to validate the proposed reanalysis-based design procedure. The presented numerical results demonstrate that the adaptive reanalysis technique affects very slightly the accuracy of the optimal solutions and does accelerate the design process, especially for large-scale structures.     

GENERALIZED FINITE SPECTRAL METHOD FOR 1D BURGERS AND KDV EQUATIONS

A generalized finite spectral method is proposed. The method is of high-order accuracy. To attain high accuracy in time discretization, the fourth-order Adams-Bashforth-Moulton predictor and corrector scheme was used. To avoid numerical oscillations caused by the dispersion term in the KdV equation, two numerical techniques were introduced to improve the numerical stability. The Legendre, Chebyshev and Her-mite polynomials were used as the basis functions. The proposed numerical scheme is validated by applications to the Burgers equation (nonlinear convection- diffusion problem) and KdV equation (single solitary and 2-solitary wave problems), where analytical solutions are available for comparison. Numerical results agree very well with the corresponding analytical solutions in all cases.     

Improved modal truncation method for eigensensitivity analysis of asymmetric matrix with repeated eigenvalues

An improved modal truncation method with arbitrarily high order accuracy is developed for calculating the second- and third-order eigenvalue derivatives and the first- and second-order eigenvector derivatives of an asymmetric and non-defective matrix with repeated eigenvalues. If the different eigenvalues λ1, λ2,, λrof the matrix satisfy |λ1| |λr| and |λs| |λs+1|(s r-1), then associated with any eigenvalue λi(i s), the errors of the eigenvalue and eigenvector derivatives obtained by the qth-order approximate method are proportional to |λi/λs+1|q+1, where the approximate method only uses the eigenpairs corresponding to λ1, λ2,, λs. A numerical example shows the validity of the approximate method. The numerical example also shows that in order to get the approximate solutions with the same order accuracy, a higher order method should be used for higher order eigenvalue and eigenvector derivatives.     

NONLINEAR VIBRATIONS OF ORTHOTROPIC SHALLOW SHELLS OF REVOLUTION

A set of nonlinearly coupled algebraic and differential eigenvalue equations of nonlinear axisymmetric free vibration of orthotropic shallow thin spherical and conical shells are formulated.following an assumed time-mode approach suggested in this paper. Analytic solutions are presented and an asymptotic relation for the amplitude-frequency response of the shells is derived. The effects of geometrical and material parameters on vibrations of the shells are investigated.     

A GENERALIZED VARIATIONAL PRINCIPLE OF COMPOSITE SHALLOW SHELLS AND ITS APPLICATION TO THE FOLDED SHELL

In this paper, a generalized variational principle of elastodynamics in compositeshallow shells with edge beams is presented,and its equivalence to corresponding basicequations. ridge conditions and boundary conditions is proved.Then this variationalprinciple is applied to the folded shell structure.By means of double series,the approximateanalytical solutions for statics and dynamics under common boundary conditions areobtained.The comparison of our results with FEM computations and experiments showsthe analytical solutions have good convergence and their accuracy is quite satisfactory.     

薄板振动特征值问题的一种改进积分方程解法

本文根据一种改进的边界元/有限元混合法求解薄板振动固有频率问题,既避开了标准的边界元法所导致的求解非代数特征值方程的困难,亦能够基本上消除通常的边界元/有限元混合法结果精度受区域内部单元划分影响较大的弊端。文中讨论了迭代算法的收敛问题,并用于薄板固有频率分析。数值结果表明,即便是在域内单元很粗疏划分的情况下,本文的方法仍能给出相当满意的结果。     

Time-harmonic BEM for 2-D piezoelectricity applied to eigenvalue problems

We will derive the fundamental generalized displacement solution, using the Radon transform, and present the direct formulation of the time-harmonic boundary element method (BEM) for the two-dimensional general piezoelectric solids. The fundamental solution consists of the static singular and the dynamics regular parts; the former, evaluated analytically, is the fundamental solution for the static problem and the latter is given by a line integral along the unit circle. The static BEM is a component of the time-harmonic BEM, which is formulated following the physical interpretation of Somigliana’s identity in terms of the fundamental generalized line force and dislocation solutions obtained through the Stroh–Lekhnitskii (SL) formalism. The time-harmonic BEM is obtained by adding the boundary integrals for the dynamic regular part which, from the original double integral representation over the boundary element and the unit circle, are reduced to simple line integrals along the unit circle.The BEM will be applied to the determination of the eigen frequencies of piezoelectric resonators. The eigenvalue problem deals with full non-symmetric complex-valued matrices whose components depend non-linearly on the frequency. A comparative study will be made of non-linear eigenvalue solvers: QZ algorithm and the implicitly restarted Arnoldi method (IRAM). The FEM results whose accuracy is well established serve as the basis of the comparison. It is found that the IRAM is faster and has more control over the solution procedure than the QZ algorithm. The use of the time-harmonic fundamental solution provides a clean boundary only formulation of the BEM and, when applied to the eigenvalue problems with IRAM, provides eigen frequencies accurate enough to be used for industrial applications. It supersedes the dual reciprocity BEM and challenges to replace the FEM designed for the eigenvalue problems for piezoelectricity.     

Axisymmetric viscoplastic deformation by the boundary element method

This paper presents a boundary element formulation and numerical implementation of the problem of small axisymmetric deformation of viscoplastic bodies. While the extension from planar to axisymmetric problems can be carried out fairly simply for the finite element method (FEM), this is far from true for the boundary element method (BEM). The primary reason for this fact is that the axisymmetric kernels in the integral equations of the BEM contain elliptic functions which cannot be integrated analytically even over boundary elements and internal cells of simple shape. Thus, special methods have to be developed for the efficient and accurate numerical integration of these singular and sensitive kernels over discrete elements. The accurate determination of stress rates by differentiation of the displacement rates presents another formidable challenge.A successful numerical implementation of the boundary element method with elementwise (called the Mixed approach) or pointwise (called the pure BEM or BEM approach) determination of stress rates has been carried out. A computer program has been developed for the solution of general axisymmetric viscoplasticity problems. Comparisons of numerical results from the BEM and FEM, for several illustrative problems, are presented and discussed in the paper. It is possible to get direct solutions for the simpler class of problems for cylinders of uniform cross-section, and these solutions are also compared with the BEM and FEM results for such cases.     

A generalized variational principle of composite shallow shells and its application to the folded shell

In this paper, a generalized variational principle of elastodynamics in composite shallow shells with edge beams is presented, and its equivalence to corresponding basic equations, ridge conditions and boundary conditions is proved. Then this variational principle is applied to the folded shell structure. By means of double series, the approximate analytical solutions for statics and dynamics under common boundary conditions are obtained. The comparison of our results with FEM computations and experiments shows the analytical solutions have good convergence and their accuracy is quite satisfactory.     

A DOMAIN DECOMPOSITION ALGORITHM WITH FINITE ELEMENT-BOUNDARY ELEMENT COUPLING

A domain decomposition algorithm coupling the finite element and the boundary element was presented. It essentially involves subdivision of the analyzed domain into sub-regions being independently modeled by two methods, i.e., the finite element method (FEM) and the boundary element method (BEM). The original problem was restored with continuity and equilibrium conditions being satisfied on the interface of the two sub-regions using an iterative algorithm. To speed up the convergence rate of the iterative algorithm, a dynamically changing relaxation parameter during iteration was introduced. An advantage of the proposed algorithm is that the locations of the nodes on the interface of the two sub-domains can be inconsistent. The validity of the algorithm is demonstrated by the consistence of the results of a numerical example obtained by the proposed method and those by the FEM, the BEM and a present finite element-boundary element (FE-BE) coupling method.     

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method (SBFEM) is a recently developed numerical method combining advantages of both finite element methods (FEM) and boundary element methods (BEM) and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors (SIFs) directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.The project supported by the National Natural Science Foundation of China (50579081) and the Australian Research Council (DP0452681)The English text was polished by Keren Wang.     

Method of fundamental solutions for multidimensional Stokes equations by the dual-potential formulation

A simple and accurate boundary-type meshless method of fundamental solutions (MFS) is applied to solve both 2D and 3D Stokes flows based on the dual-potential formulation of velocity potential and stream function vector. Using the dual-potential concept, the solutions of both 2D and 3D Stokes flows are obtained by combining the much simpler fundamental solutions of Laplace (potential) and bi-harmonic equations without using the complicated singular fundamental solutions such as Stokeslets and their derivatives as well as source doublet hypersingularity. The developed algorithm is used to test five numerical experiments for 2D flows: (1) circular cavity, (2) wave-shaped bottom cavity and (3) circular cavity with eccentric rotating cylinder; and for 3D flows: (4) a uniform flow passing a sphere and (5) a uniform flow passing a pair of spheres. Good results are obtained as comparing with solutions of analytical and numerical methods such as FEM, BEM and other meshfree schemes.     

A comparative study on algorithms of boundary element solution for free torsional vibration of bodies of revolution

In this paper the functionu=rsin in cylindrical coordinates (r,,z) is introduced into the equation for free torsional vibration of bodies of revolution (where=v / r represents the angle of twist). With the static fundamental solution (–1 /R) a mixed BEM / FEM equation is derived. The domain integral term in the equation is discretized by Serendipity elements instead of commonly used constant value finite elements in the literature. The equation is an algebraic eigenvalue one. The dynamic fundamental solution (e ^1R /R) is also used for deriving the other mixed BEM / FEM equation. An appropriate iterative solution procedure is described. An algebraic eigenvalue equation can be obtained and its solution accuracy is almost interior meshing independent. A number of examples are studied. The results show the good economy and high accuracy of the algorithms proposed.The Project is Supported by National Natural Science Foundation of China.     

层状弹性半空间轴对称动力问题的奇异解

本文利用Laplace-Hankel联合变换及传播矩阵技术导出了任意层数的层状弹性半空间轴对称动力问题时域奇异解的一般解析表达式,并给出了奇异解数值化实施的计算方法。文末的实例计算表明了本文给出解答的正确性以及数值化实施的可靠性,从而为进一步用边界元法直接解决由于层状介质而引起的非匀质动力问题开拓了一条潜在的途径。     

Vibration analysis of plates using the multivariable spline element method

The vibration analysis of plates using the multivariable spline element method is presented in this paper. The spline functions are applied to construct bending moments, twisting moments and transverse displacement field functions. The spline equations of eigenvalue problems with multiple variables of vibration of plates are derived based on the Hellinger-Reissner mixed variational principle. For simplicity, the boundary conditions which consist of three local spline points are amended to fit any specified boundary conditions. Several numerical solutions of plate vibration analysis are presented which illustrate the accuracy and convergence of the method.