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     

Characteristic Galerkin method for convection-diffusion equations and implicit algorithm using precise integration

This paper presents a finite element procedure for solving transient, multidimensional convection-diffusion equations. The procedure is based on the characteristic Galerkin method with an implicit algorithm using precise integration method. With the operator splitting procedure, the precise integration method is introduced to determine the material derivative in the convection-diffusion equation, consequently, the physical quantities of material points. An implicit algorithm with a combination of both the precise and the traditional numerical integration procedures in time domain in the Lagrange coordinates for the characteristic Galerkin method is formulated. The stability analysis of the algorithm shows that the unconditional stability of present implicit algorithm is enhanced as compared with that of the traditional implicit numerical integration procedure. The numerical results validate the presented method in solving convection-diffusion equations. As compared with SUPG method and explicit characteristic Galerkin method, the present method gives the results with higher accuracy and better stability. The project sponsored by the State Scientific and Technological Commission of China through “China State Key Project: the Theory and Methodology for Scientific and Engineering Computations with Large Scale”, the National Natural Science Foundation of China and the European Commission Research Project CI1^*CT94-0014.     

A new simple method of implicit time integration for dynamic problems of engineering structures

This paper presents a new simple method of implicit time integration with two control parameters for solving initial-value problems of dynamics such that its accuracy is at least of order two along with the conditional and unconditional stability regions of the parameters. When the control parameters in the method are optimally taken in their regions, the accuracy may be improved to reach of order three. It is found that the new scheme can achieve lower numerical amplitude dissipation and period dispersion than some of the existing methods, e.g. the Newmark method and Zhai’s approach, when the same time step size is used. The region of time step dependent on the parameters in the new scheme is explicitly obtained. Finally, some examples of dynamic problems are given to show the accuracy and efficiency of the proposed scheme applied in dynamic systems. The project supported by the National Key Basic Research and Development Foundation of the Ministry of Science and Technology of China (G2000048702, 2003CB716707), the National Science Fund for Distinguished Young Scholars (10025208), the National Natural Science Foundation of China (Key Program) (10532040), the Research Fund for Oversea Chinese (10228028). The English text was polished by Keren Wang.     

A finite element model for numerical simulation of thermo-mechanical frictional contact problems

Two kinds of variational principles for numerical simulation of heat transfer and contact analysis are respectively presented. A finite element model for numerical simulation of the thermal contact problems is developed with a pressure dependent heat transfer constitutive model across the contact surface. The numerical algorithm for the finite element analysis of the thermomechanical contact problems is thus developed. Numerical examples are computed and the results demonstrate the validity of the model and algorithm developed. The project supported by the National Key Basic Research Special Foundation (G1999032805), the National Natural Science Foundation of China (50178016, 10225212) and the Foundation for University Key Teacher by the Ministry of Education of China     

关于动力分析精细积分算法精度的讨论

对动力问题分析的精细积分算法的精度问题进行深入研究,并在此基础上提出对原有的算法的改进策略,改进后的算法可以较好地克服算法精度对积分时间步长的依赖性问题。     

A note for analysis of thermo-mechanical contact problems

A discussion about the bifurcation and non-uniqueness of solutions in the analysis of thermo-mechanical contact problems with initial gap is given. Without loss of generality, a mechanical contact problem coupled with steady heat transfer is studied and an example of non-uniqueness of solutions caused by the thermo-mechanical mechanism is presented. The important work is that the non-uniqueness of solutions, which is different from that found in the analysis of the traditional frictional contact problems, is studied in detail. The possible oscillation and non-convergence problems in the iteraction process of the numerical computation are discussed, and an enhanced algorithm is put forward to overcome the difficulties. Project sypported by the National Natural Science Foundation of China (Nos. 50178016, 10225212 and 19872016), the National Key Basic Research Special Foundation (No. G1999032805) and the Foundation for University Key Teacher by the Ministry of Education.     

Mixed energy method for solution of quadratic programming problems and elastic-plastic analysis of truss structures

A new algorithm for the solution of quadratic programming problems is put forward in terms of the mixed energy theory and is further used for the incremental solution of elastic-plastic truss structures. The method proposed is different from the traditional one, for which the unknown variables are selected just in one class such as displacements or stresses. The present method selects the variables in the mixed form with both displacement and stress. As the method is established in the hybrid space, the information found in the previous incremental step can be used for the solution of the present step, making the algorithm highly efficient in the numerical solution process of quadratic programming problems. The results obained in the examples of the elastic-plastic solution of the truss structures verify what has been predicted in the theoretical analysis. Project supported by the National Natural Science Foundation of China (No. 50178916, No. 19732020 and No. 19872016), the National Key Basic Research Special Foundation (No. G1999032805), the Special Funds for Major State Basic Research Projects and the Foundation for University Key Teachers by the Ministry of Education of China.     

AN ANALYSIS OF 3D FINITE CRACKED BODIES

The new-type traction boundary integral equations developed by Hu and with no hypersingular integral are applied to analysis of 3D finite cracked bodies. A numerical algorithm for general 3D problems and a semi-analytical one for axisymmetric problems are presented. Some examples of thick plates and cylindrical columns including penny-shaped crack(s), and rectangular plates including an elliptical crack normal to the surface are analyzed. The comparison between present results and those in literature shows the high accuracy and effectiveness of the present method. Project supported by the National Natural Science Foundation of China (No. 19972060) and the Foundation of Key Laboratory of the Ministry of Education, Tongji University.     

Solving frictional contact problems by two aggregate-function-based algorithms

Three dimensional frictional contact problems are formulated as linear complementarity problems based on the parametric variational principle. Two aggregate-function-based algorithms for solving complementarity problems are proposed. One is called the self-adjusting interior point algorithm, the other is called the aggregate function smoothing algorithm. Numerical experiment shows the efficiency of the proposed two algorithms. The project supported by the National Natural Science Foundation of China(10225212, 50178016, 10302007), the National Key Basic Research Special Foundation and the Ministry of Education of China The English text was polished by Ron Marshall.     

?????????????????????Krylov?????????

提出了一种新的精细时程积分法来求解大型动力系统. 结合Krylov子空间法、培德级数 近似以及一般载荷的维数扩展法,进一步提高精细时程积分法的计算效率. 利用维数扩展法 避免计算微分方程特解,并可处理任意载荷. 对于大型动力系统,通过Krylov子空间的降维 分析将问题转化到一个子空间,计算效率得到极大提高. 对于迭代次数$N$ 的选择作了详细讨论,进一步提高了计算效率.     

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     

无条件稳定的更新精细积分方法

提出将Pade逼近与精细积分方法中的指数矩阵运算技巧结合起来,建立了精细积分法的更新形式及计算过程,对该更新精细积分方法的稳定性进行了论证与探讨.结果表明,该更新精细积分方法是无条件稳定的,整个积分方法的精度取决于所取Pade逼近的阶数与高斯积分点的数量.数值例题也显示了该方法的高效率及其可行性.     

A method for the analysis of dynamic response of structure containing non-smooth contactable interfaces

A novel single-step method is proposed for the analysis of dynamic response of visco-elastic structures containing non-smooth contactable interfaces. In the method, a two-level algorithm is employed for dealing with a nonlinear boundary condition caused by the dynamic contact of interfaces. At the first level, and explicit method is adopted to calculate nodal displacements of global viscoelastic system without considering the effect of dynamic contact of interfaces and at the second level, by introducing contact conditions of interfaces, a group of equations of lower order is derived to calculate dynamic contact normal and shear forces on the interfaces. The method is convenient and efficient for the analysis of problems of dynamic contact. The accuracy of the method is of the second order and the numerical stability condition is wider than that of other explicit methods. The project supported by the National Natural Science Foundation of China (59578032) and the Key Project of the Ninth Five-Year Plan (96221030202)     

Time discontinuous Galerkin finite element method for dynamic analyses in saturated poro-elasto-plastic medium

A time-discontinuous Galerkin finite element method for dynamic analyses in saturated poro-elasto-plastic medium is proposed. As compared with the existing discontinuous Galerkin finite element methods, the distinct feature of the proposed method is that the continuity of the displacement vector at each discrete time instant is automatically ensured, whereas the discontinuity of the velocity vector at the discrete time levels still remains. The computational cost is then obviously reduced, particularly, for material non-linear problems. Both the implicit and explicit algorithms to solve the derived formulations for material non-linear problems are developed. Numerical results show a good performance of the present method in eliminating spurious numerical oscillations and providing with much more accurate solutions over the traditional Galerkin finite element method using the Newmark algorithm in the time domain. The project supported by the National Natural Science Foundation of China (19832010, 50278012, 10272027) and the National Key Basic Research and Development Program (973 Program, 2002CB412709)     

Retrospective time integral scheme and its applications to the advection equation

To put more information into a difference scheme of a differential equation for making an accurate prediction, a new kind of time integration scheme, known as the retrospective (RT) scheme, is proposed on the basis of the memorial dynamics. Stability criteria of the scheme for an advection equation in certain conditions are derived mathematically. The computations for the advection equation have been conducted with its RT scheme. It is shown that the accuracy of the scheme is much higher than that of the leapfrog (LF) difference scheme. The project supported by the National Key Program for Developing Basic Sciences (G1999043408 and G1998040901-1) and the National Natural Sciences Foundation of China (40175024 and 40035010)     

Application of MLPG in Large Deformation Analysis

Two-dimensional large deformation analysis of hyperelastic and elasto-plastic solids based on the Meshless Local Petrov–Galerkin method (MLPG) is presented. A material configuration based the nonlinear MLPG formulation is introduced for the large deformation analysis of both path-dependent and path-independent materials. The supports of the MLS approximation functions cover the same sets of nodes during material deformation, thus the shape function needs to be computed only in the initial stage. The multiplicative hyperelasto-plastic constitutive model is adopted to avoid objective time integration for stress update in large rotation. With this constitutive model, the computational formulations for path-dependent and path-independent materials become identical. Computational efficiency of the nonlinear MLPG method is discussed and optimized in several aspects to make the MLPG an O(N) algorithm. The numerical examples indicate that the MLPG method can solve large deformation problems accurately. Moreover, the MLPG computations enjoy better convergence rate than the FEM under very large particle distortion.The project supported by the National Natural Science Foundation of China (10472051). The English text was polished by Keren Wang.     

An improved level-set re-initialization solver

Re-initialization procedure in level-set interface capturing method were investigated. The algorithm accomplishes the re-initialization step through locking the interface positions. Better accuracy was obtained both on the interface positions and the total fluid volume keeping. Though one mare step of the interpolations is added in the procedure, there is no significant increase in total machine time spent.     

Formulations of the three-dimensional discontinuous deformation analysis method

This paper extends the original 2D discontinuous deformation analysis (DDA) method proposed by Shi to 3D cases, and presents the formulations of the 3D DDA. The formulations maintain the characteristics of the original 2D DDA approach. Contacts between the blocks are detected by using Common-Plane (C-P) approach and the non-smooth contact, such as of vertex-to-vertex, vertex-to-edge and edge-to-edge types, can be handled easily based on the C-P method. The matrices of equilibrium equations have been given in detail for programming purposes. TheC program codes for the 3D DDA are developed. The ability and accuracy of the formulations and the program are verified by the analytical solutions of several dynamic examples. The robustness and versatility of the algorithms presented in this paper are demonstrated with the aid of an example of scattering of densely packed cubes. Finally, implications and future extensions are discussed. The project supported by the National Natural Science Foundation of China (50139010)     

获得热传导问题“拟解析解”的精细积分算法

研究并讨论了用于热传导分析有限元解的精细积分算法，算法很好地克服了传统方法求解时的单调性问题，且对空间离散后所获得方程的解是解析的，因而算法的解将具有“拟解析解”的意义，论文证明了算法单调性．     

A SOLUTION PROCEDURE FOR LOWER BOUND LIMIT AND SHAKEDOWN ANALYSIS BY SGBEM

The symmetric Galerkin boundary element method (SGBEM) instead of the finite element method is used to perform lower bound limit and shakedown analysis of structures. The self-equilibrium stress fields are constructed by a linear combination of several basic self-equilibrium stress fields with parameters to be determined. These basic self-equilibrium stress fields are expressed as elastic responses of the body to imposed permanent strains and obtained through elastic-plastic incremental analysis. The complex method is used to solve nonlinear programming and determine the maximal load amplifier. The limit analysis is treated as a special case of shakedown analysis in which only the proportional loading is considered. The numerical results show that SGBEM is efficient and accurate for solving limit and shakedown analysis problems. Project supported by the National Natural Science Foundation of China (No. 19902007), the National Foundation for Excellent Doctorial Dissertation of China (No. 200025) and the Basic Research Foundation of Tsinghua University.