A combined programming and iteration algorithm for finite element analysis of three-dimensional contact problems

Comparing with two-dimensional contact problems, three-dimensional frictional contact problems are more difficult to deal with, because of the unknown slip direction of the tangential force and enormous computing time. In order to overcome these difficulties, a combined PQP (Parametric Quadratic Programming) and iteration method is derived in this paper. The iteration algorithm, which alleviates the difficulty of unknown slip direction, is used along with the PQP method to cut down computing costs. Numerical example is given to demonstrate the validity of the present algorithm. The project supported by the Machinary and Electronics Ministry of China     

An iterative parallel algorithm of finite element method

In this paper, a parallel algorithm with iterative form for solving finite element equation is presented. Based on the iterative solution of linear algebra equations, the parallel computational steps are introduced in this method. Also by using the weighted residual method and choosing the appropriate weighting functions, the finite element basic form of parallel algorithm is deduced. The program of this algorithm has been realized on the ELXSI-6400 parallel computer of Xi'an Jiaotong University. The computational results show the operational speed will be raised and the CPU time will be cut down effectively. So this method is one kind of effective parallel algorithm for solving the finite element equations of large-scale structures.     

一种新的并行自适应网格有限元算法

给出了一种新的适用于流体力学问题的并行自适应有限元算法。首先，基于初始稀网格上获得的事后误差估算值，应用反复谱对剖分方法对初网格进行划分，使各子域上总体误差近似相等，从而解决并行自适应计算中的负载平衡问题。然后在各处理器上独立地求解整体问题，并进行指定子域上的网格自适应处理。最后将各子域上的自适应网格组合成一个整体网格，应用基于粘接元技术的区域分裂法在该网格上获得最终解。文末给出了数值实验结果。     

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 parallel two-level finite element method for the Navier-Stokes equations

Based on domain decomposition, a parallel two-level finite element method for the stationary Navier-Stokes equations is proposed and analyzed. The basic idea of the method is first to solve the Navier-Stokes equations on a coarse grid, then to solve the resulted residual equations in parallel on a fine grid. This method has low communication complexity. It can be implemented easily. By local a priori error estimate for finite element discretizations, error bounds of the approximate solution are derived. Numerical results are also given to illustrate the high efficiency of the method.     

A contact searching algorithm for contact-impact problems

A new contact searching algorithm for contact-impact systems is proposed in this paper. In terms of the cell structure and the linked-list, this algorithm solves the problem of sorting and searching contacts in three dimensions by transforming it to a retrieving process from two one-dimensional arrays, and binary searching is no longer required. Using this algorithm, the cost of contact searching is reduced to the order ofO(N) instead ofO(Nlog₂ N) for traditional ones, whereN is the node number in the system. Moreover, this algorithm can handle contact systems with arbitrary mesh layouts. Due to the simplicity of this algorithm it can be easily implemented in a dynamic explicit finite element program. Our numerical experimental result shows that this algorithm is reliable and efficient for contact searching of three dimensional systems. The project supported by the National Natural Science Foundation of China (59875045), and the State Key Laboratory of Automobile Safety and Energy Saving (K9705)     

随机有限元分析的二级区域分解并行求解算法

采用蒙特卡罗模拟(MCS)和加权积分法对二维问题进行随机有限元分析。尽管MCS方法对任何有确定解的问题都具有求解精度高的优点,但由于求解所需的计算量巨大使其应用受到限制。利用并行求解技术可有效地处理这种密集型计算问题。基于有限元分裂对接法(FETI)的并行特性并利用预处理共轭梯度法(PCG)的求解高效性,结合整体子区域实现(GSI-PCG)和FETI法,提出二级求解算法,并在工作站机群上实现了数值算例。算例计算结果表明本文GSI(PCG)-FETI算法具有较高的并行加速比和并行效率,具有良好的性能,可有效地进行二维问题的随机有限元分析。     

非线性动力有限元重叠区域分裂的隐式并行算法

针对大规模结构非线性瞬态动力分析非常耗时,提出了相应的并行算法。该算法采用无条件稳定的Ne-wmark-β方法(平均加速技术)进行时间积分,并结合区域分裂技术进行分析。它不同于已有的采用非重叠区域的并行算法,而是采用重叠区域的并行算法。对给定结构有限元分析的质量、阻尼、刚度矩阵进行分裂可推出重叠区域分裂算法的计算公式。为改善每一步的求解,采用预估和校正子方案。编写了该算法的程序,在工作站机群上实现了数值算例,验证了算法的性能。计算结果表明该算法优于非重叠区域分裂算法。     

流体-结构耦合问题的有限元并行计算研究

流体-结构耦合问题广泛存在于各种工程领域,本文采用ALE显式有限元法求解该类问题,并对该方法的并行性进行讨论。同时根据流体-结构耦合问题与ALE显式有限元的计算特点,在坐标递归分区方法的基础上设计并程序实现了基于流体-结构耦合均衡的分区算法。通过与坐标递归分区方法的计算结果相比较,对于流体-结构耦合问题的求解,耦合均衡并行分区方法具有更好的加速比和并行效率。     

Mass redistribution method for finite element contact problems in elastodynamics

This paper is devoted to a new method dealing with the semi-discretized finite element unilateral contact problem in elastodynamics. This problem is ill-posed mainly because the nodes on the contact surface have their own inertia. We introduce a method based on an equivalent redistribution of the mass matrix such that there is no inertia on the contact boundary. This leads to a mathematically well-posed and energy conserving problem. Finally, some numerical tests are presented.     

一种新型并行化有限元结构模态分析集成系统

以成熟有限元软件的模态分析流程和大型稀疏矩阵特征值的并行求解为基础,开发出一种基于大规模并行机的新型有限元结构模态分析系统。通过对串行CAE软件的二次开发,将模态分析过程中计算量最大的特征值求解部分代之以并行计算。针对并行机特性以隐式重启动Lanczos算法为基础,编写了基于MPI的特征值并行求解程序,并通过实际算例验证了并行程序的加速比和扩展性;同时实现并行程序与其它串行分析步骤的无缝集成,使集成系统的界面友好,操作方便。本系统使结构模态分析的规模和速度大幅度提高,以大型CAE软件MSC／NASTRAN为并行化求解器开发平台,在“神威Ⅰ”超级计算机上验证了其可行性和高效性。     

A variational principle and finite element formulation for multi-physics PLZT ceramics

This research presents an extended variational principle and a finite element formulation for multi-physics analysis of PLZT (lanthanum zirconate titanate) ceramics by including the photovoltaic and optothermic effects. The photo-induced electrical and thermal, and mechanical fields, as well as the heat source effect due to light illumination are all considered in the formulation. A generalized variation approach for advanced multi-physics PLZT problems is developed and then the relevant finite element formulation is established.     

一类有限元结构分析的EBE计算方法

本文利用ＥＢＥ策略￣［１］的基本思想，给出一类有限元结构分析的ＥＢＥ计算方法，即ＥＢＥ共轭梯度法ＥＢＥ－ＣＧ和ＥＢＥ预处理共轭梯度法ＥＢＥ－ＰＣＧ，这类方法避免了传统有限元结构分析中总刚度阵的组集而可大大降低存储量要求。同时它们还特别适合在各种粒度下的多处理机系统上实现。初步数值试验结果表明：这类ＥＢＥ有限元结构分析方法对串行和并行计算都是很有效的。     

The random variational principle in finite deformation of elasticity and finite element method

ＴＨＥＲＡＮＤＯＭＶＡＲＩＡＴＩＯＮＡＬＰＲＩＮＣＩＰＬＥＩＮＦＩＮＩＴＥＤＥＦＯＲＭＡＴＩＯＮＯＦＥＬＡＳＴＩＣＩＴＹＡＮＤＦＩＮＩＴＥＥＬＥＭＥＮＴＭＥＴＨＯＤＧａｏＨａｎｇ－ｓｈａｎ（高行山）（ＮｏｒｔｈｗｅｓｔｅｎＰｏｌｙｔｅｃｈｎｉｃａｌＵ...     

The random variational principle and finite element method

In this paper, we introduced the random materials, geometrical shapes, force and displacement boundary condition directly into the functional variational formulations and developed a unified random variational principle and finite element method with the small parameter perturbation method. Numerical examples showed that the methods have the advantages of the simple and convenient program implementation, and are effective for the random mechanics problems.     

The boundary element method applied to elastic contact problems with friction

In the current paper the boundary integral equations (BIE) for elastic contact problems with friction are derived from the incremental virtual work principle. After introducing contact conditions of adhesion and slip into BIE all variants on boundary are made to discretize by quadratic isoparametric boundary element. In the current paper not only an auto-increment loading law is presented but also the iterative calculation laws for open, slip and adhesion condition are given. The results of numerical examples are satisfactory.     

参变量变分原理的提出、发展与应用

参变量变分原理及其参数二次规划算法是由钟万勰院士1985年针对弹性接触边界非线性问题首次提出来的,经过将近40年的不断发展,目前参变量变分原理已经成功应用于各个领域,其中包括弹塑性分析、接触问题、润滑力学、岩土力学、变刚度杆系结构、先进材料性能分析、材料的蠕变与损伤、柔性结构力学和LQ最优控制等各个工程领域。本文首先回顾了参变量变分原理的起源,介绍了参变量变分原理的基本概念,然后以弹塑性分析问题为例,阐明建立参变量变分原理的理论模型以及实现数值参数二次规划求解原理,最后详细回顾了参变量变分原理的基本理论与相应数值算法在各个领域的发展及其工程应用,展示了参变量变分原理在求解各类非线性问题的特色与优势。     

A version of the finite element method for frictional contact problems

We propose a method for solving frictional contact problems which is based on including the generalized coordinates of absolutely rigid bodies in the degrees of freedom of the system under study and on varying the functional of the variational problem with respect to these coordinates. As a result, one can include the generalized coordinates or the energy-conjugate generalized forces directly in the right-hand side of the resolving system of equations, which permits easily taking into account any laws of motion or loading of absolutely rigid bodies.     

The parametric variational principle for elastoplasticity

A parametric variational principle, the parametric minimum potential energy principle (abbreviated to PMPEP), is presented for the elastoplastic problems. The principle proposed is free from the restraint of Drucker's postulate and consequently suitable for solving the nonassociated plastic flow problems in rock, soil, concrete and other geomaterials.     

NON-INTERIOR SMOOTHING ALGORITHM FOR FRICTIONAL CONTACT PROBLEMS

A new algorithm for solving the three-dimensional elastic contact problem with friction is presented. The algorithm is a non-interior smoothing algorithm based on an NCP-function. The parametric variational principle and parametric quadratic programming method were applied to the analysis of three-dimensional frictional contact problem. The solution of the contact problem was finally reduced to a linear complementarity problem, which was reformulated as a system of nonsmooth equations via an NCP-function. A smoothing approximation to the nonsmooth equations was given by the aggregate function. A Newton method was used to solve the resulting smoothing nonlinear equations. The algorithm presented is easy to understand and implement. The reliability and efficiency of this algorithm are demonstrated both by the numerical experiments of LCP in mathematical way and the examples of contact problems in mechanics.