弹塑性接触问题的非光滑非线性方程组方法

将求解三维弹性摩擦接触问题的非光滑非线性方程组方法推广到弹塑性(Mises材料)情形，提出了两种应用方法：一种是将非光滑非线性方程组方法和求解弹塑性问题常用的Newton—Raphson迭代方法结合起来；另一种是将问题写成统一的非光滑非线性方程组，直接求解。数值算例验证了两种方法的有效性，并进行了结果比较。     

三维摩擦接触问题不动点算法

三维摩擦接触问题具有多重非线性性质,使求解变得比较困难,为了解决此问题,建立了三维摩擦接触问题的模型,介绍了一种非光滑混合不动点算法.该算法克服了在接触面上由于可能的滑动状态有无穷多个而难以确定的难点,算法未引入人工变量,计算量较小,计算结果精确满足接触状态条件,收敛性得到保证.根据此算法编制程序,将叠合悬臂梁算例数值计算结果与商用有限元软件进行比较,也表明了不动点算法的有效性.     

接触问题的非线性互补—接触柔度法

基于势能互补原理，不引入额外松弛变量，对库仑摩擦定律未做预先性化近似，提出了三维弹性摩擦接触问题的非线性互补－接触柔度法，收敛性和收敛速率得到了严格理论保证。数值实验结果表明方法可靠，有效，低存储量，高精度。     

非光滑方程组方法在求解非匹配网格接触问题中的应用

将非光滑方程组方法与Mortar StS接触模型(Mortar Segment-to-Segment)相结合,来求解接触面网格非匹配时的弹性接触问题.其中,非光滑方程组方法是求解弹性摩擦接触问题的有效方法,具有精确满足接触条件、迭代算法收敛性有理论保证的优点,但目前仅用于求解网格匹配的接触问题.Mortar StS接触模型可以较为方便地处理网格非匹配接触问题,其特点是不引入过多约束,满足接触分片检验条件,但目前大都采用“试验-误差”迭代方法求解控制方程,对于复杂接触问题,其收敛性不易保证.因此,将二者结合来处理网格非匹配接触问题,既可以提高求解精度,又能使得算法的收敛性得到理论保证.数值算例对接触分片检验和算法的计算精度进行了验证.     

接触问题的非线性互补－接触柔度法

基于势能互补原理，不引入额外松弛变量，对库仑摩擦定律未做预先线性化近似，提出了三维弹性摩擦接触问题的非线性互补－接触柔度法，收敛性和收敛速率得到了严格理论保证。数值实验结果表明方法可靠、有效、低存储量、高精度     

样条边界元法分析圆形基础板非线性非光滑接触

本文用样条边界元法分析轴对称圆形基础板的非线性非光滑接触。除事先给出外力与地基表面位移的对应关系外,地基的模型、板的厚薄(Reissner/Kirchhoff型)和外荷载沿径向的分布均可是任意的。文中将地基反力看作板的外载,对变量沿径向采用不等距祥条插值,并用无地基效应的板基本解。算例表明本文方法只要用少量自由度,便可得到令人满意的结果。     

三维摩擦接触问题的增量线性互补方法

借助罚因子建立了三维摩擦接触本构模型,接触条件被表示为类似于非关联流动的弹塑性本构关系的形式。采用增量描述,对Coulomb摩擦定律采用Taylor展开作线性近似,导出了接触问题的互补虚功方程,然后基于有限元离散建立了三维摩擦接触问题的增量线性互补方法。数值算例表明了本文方法的有效性。     

接触问题的互补变分原理及非线性互补模型

给出了接触问题的互补虚功原理，互补能量原理及变分不等式等变分提法，并讨论了三类提法的关系。其此建立了有限维非线性互补模型并给出解的存在唯一性定理。文中还阐明了有，无摩擦接触问题之间以及二维和三维有摩擦接触问题之间的本质区别。     

非光滑优化的双层规划模型及内点算法

对于包含接触约束的非光滑结构优化问题,其非光滑性体现在状态函数并不是处处可微的,针对含有应力约束及接触约束的非光滑结构优化问题,建立了一种双层规划模型,避免了求解时非光滑性所带来的问题,同时提出了一种迭代算法,用对偶内点二次规划进行分析,线性规划进行优化,算例表明这种方法十分有效。     

Block pivot methods for solving frictional contact problems

Based on elementary group theory, the block pivot methods for solving two-dimensional elastic frictional contact problems are presented in this paper. It is proved that the algorithms converge within a finite number of steps when the friction coefficient is “relative small”. Unlike most mathematical programming methods for contact problems, the block pivot methods permit multiple exchanges of basic and nonbasic variables. The project supported by the National Natural Science Foundation of China     

摩擦接触问题数学规划法的收敛性和算法的改进

对摩擦接触问题的数学规划法进行了收敛性分析，并对算法提出一种改进，对“较小”摩擦系数问题改进算法严格保证收敛。在“较大”摩擦系数情形，问题提法本身有可能发生解不唯一和不存在现象。     

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.     

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.     

F0RMULATION OF BOUNDARY ELEMENT-LINEAR COMPLEMENTARY EQUATION FOR THE FRICTIONAL ELASTIC CONTACT PROBLEMS

Boundary element-linear complementary equations are formulated to solve elasticcontact problems with Coulomb frictions.It is also a new attempt to solve free boundaryproblems in solid mechanics by means of boundary element-mathematical programmingtechniques.     

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     

Plane contact problem for periodic laminated composite involving frictional heating

Summary The stationary problem of a rigid thermally insulated punch sliding over the boundary surface of a periodic two-layered thermoelastic half-space is considered. The heat generated in the contact area is assumed to be caused by frictional forces. The problem is formulated within the framework of thermoelasticity with microlocal parameters, and it is reduced to a system of two integral equations, which is solved numerically. The effects connected with the composite structure are analyzed.