层合板分析的无网格局部Petrov-Galerkin方法

基于Kirchhoff均匀各向异性板控制方程的等效积分弱形式和对挠度函数采用移动最小二乘近似函数进行插值, 进一步研究无网格局部Petrov-Galerkin方法在纤维增强对称层合板弯曲问题中的应用。该方法不需要任何形式的网格划分, 所有的积分都在规则形状的子域及其边界上进行,其问题的本质边界条件采用罚因子法来施加。通过数值算例和与其他方法的结果比较, 表明无网格局部Petrov-Galerkin法求解层合薄板弯曲问题具有解的精度高、收敛性好等一系列优点。      

薄板屈曲分析的局部Petrov-Galerkin方法

利用薄板控制微分方程的等效积分对称弱形式和对变量(挠度)采用移动最小二乘近似函数进行插值,研究了无网格局部Petrov-Galerkin方法在薄板屈曲问题中的应用。它不需要任何形式的网格划分,所有的积分都在规则形状的子域及其边界上进行,并用罚因子法施加本质边界条件。数值算例表明,无网格局部Petrov-Galerkin法不但能够求解弹性静力学问题,而且在求解弹性稳定性问题时仍具有收敛快,稳定性好,精度高的特点。     

基于滑动Kriging插值的无网格MLPG法求解结构动力问题

利用基于滑动Kriging插值的无网格局部Petrov-Galerkin (MLPG) 法来求解二维结构动力问题,Heaviside分段函数作为局部弱形式的权函数并采用精细积分法来离散时间域。基于滑动Kriging插值构造的形函数满足Kronecker Delta性质,因此可以直接施加本质边界条件。刚度矩阵形成过程中只涉及到边界积分,而没有涉及到区域积分和奇异积分。计算结果表明：基于滑动Kriging插值的MLPG法具有模拟简单、计算精度高等优点。     

摩擦接触裂纹问题的扩展有限元法

扩展有限元法(XFEM)是一种在常规有限元框架内求解强和弱不连续问题的新型数值方法。扩展有限元法分析闭合型裂纹时,必须考虑裂纹面间的接触问题。已有文献均采用迭代法求解裂纹面的接触问题。该文建立了闭合型摩擦裂纹问题的扩展有限元线性互补模型,将裂纹面非线性摩擦接触转化为一个线性互补问题求解,不需要迭代求解。算例分析说明了该方法的正确性和有效性,同时表明扩展有限元法结合线性互补法求解接触问题具有较好的前景。     

高分子复合材料接触力学的有限元研究

采用有限元方法，并结合边界元方法，研究了高分子复合材料的力学问题，用可能接触区接触节点对方法，比例加载方法和库仑摩擦定律来模拟该问题。在受载过程中，随载的增加而逐渐改变接触边界条件。由于边界条件的不定性，导致了该问题的非线性，故了数学规划方法求解高度不定性的弹性接触方程组，进而求得接触区的接触位移和接触压力。     

热弹性动力学耦合问题的插值型移动最小二乘无网格法研究

该文基于插值型移动最小二乘法,将无网格局部Petrov-Galerkin（MLPG）法用于二维耦合热弹性动力学问题的求解。修正的Fourier热传导方程和弹性动力控制方程通过加权余量法来离散,Heaviside分段函数作为局部弱形式的权函数,从而得到描述热耦合问题的二阶常微分方程组。然后利用微分代数方法,温度和位移作为辅助变量,将上述二阶常微分方程组转换成常微分代数系统,采用Newmark逐步积分法进行求解。该方法无需Laplace变换可直接得到温度场和位移场数值结果,同时插值型移动最小二乘法构造的形函数由于满足Kroneckerdelta特性,因此能直接施加本质边界条件。最后通过两个数值算例来验证该方法的有效性。     

线弹性断裂力学问题的扩展自然单元法

为了更加有效地求解线弹性断裂问题,提出了扩展自然单元法。该方法基于单位分解的思想,在自然单元法的位移模式中加入扩展项表征不连续位移场和裂纹尖端奇异场。通过水平集方法确定裂纹面和裂纹尖端区域,并基于虚位移原理推导了平衡方程的离散线性方程。由于自然单元法的形函数满足Kronecker delta函数性质,本质边界条件易于施加。混合模式裂纹的应力强度因子由相互作用能量积分方法计算。数值算例结果表明扩展自然单元法可以方便地求解线弹性断裂力学问题。     

有摩擦弹性接触问题边界元分析的一种新方法

在将切向接触力与切向相对位移的关系表为带有惩罚因子的线性互补形式后,直接利用法向接触力与法向相对位移的互补关系,结合边界元技术,本文给出了一种新的求解有摩擦弹性接触问题的数学规划法。     

基于欧拉插值的最小二乘混合配点法在弹性力学平面问题中的应用

由于常规配点型无网格法存在求解不稳定、精度差和求解高阶导数等问题,提出了基于欧拉插值的最小二乘混合配点法。该方法同时以位移和应变作为未知量,通过欧拉插值将未知变量的导数表达出来,同时在插值中引入高斯权函数,并代入微分方程,从而形成以位移和应变为未知数的超定方程组,然后形成最小二乘意义下的法方程,法方程和相应的位移边界条件、应力边界条件一起形成定解体系。该方法不需要域积分,是一种真正的无网格法。一些典型的弹性力学平面问题表明本文方法具有良好的精度。     

轴对称动力学问题的无网格自然邻接点Petrov-Galerkin法

基于无网格自然邻接点Petrov-Galerkin法,提出了复杂轴对称动力学问题求解的一条新途径。几何形状和边界条件的轴对称特点,将原来的空间问题转化为平面问题求解。计算时仅仅需要横截面上离散节点的信息,无论积分还是插值都不需要网格。自然邻接点插值构造的试函数具有Kronecker delta函数性质,因此能够直接准确地施加本质边界条件。有限元三节点三角形单元的形函数作为权函数,可以减少域积分中被积函数的阶次,提高了计算效率。数值算例结果表明,本文提出的方法对求解轴对称动力学问题是行之有效的。     

A computational method for frictional contact problem using finite element method

Using the finite element method a numerical procedure is developed for the solution of the two-dimensional frictional contact problems with Coulomb's law of friction. The formulation for this procedure is reduced to a complementarity problem. The contact region is separated into stick and slip regions and the contact stress can be solved systematically by applying the solution technique of the complementarity problem. Several examples are given to demonstrate the validity of the present formulation.     

Finite element analysis of 3D elastic–plastic frictional contact problem for Cosserat materials

The objective of this paper is to develop a finite element model for 3D elastic–plastic frictional contact problem of Cosserat materials. Because 3D elastic–plastic frictional contact problems belong to the unspecified boundary problems with nonlinearities in both material and geometric forms, a large number of calculations are needed to obtain numerical results with high accuracy. Based on the parametric variational principle and the corresponding quadratic programming method for numerical simulation of frictional contact problems, a finite element model is developed for 3D elastic–plastic frictional contact analysis of Cosserat materials. The problems are finally reduced to linear complementarity problems (LCP). Numerical examples show the feasibility and importance of the developed model for analyzing the contact problems of structures with materials which have micro-polar characteristics.     

POST-BUCKLING ANALYSIS WITH FRICTIONAL CONTACTS COMBINING COMPLEMENTARITY RELATIONS AND AN ARC-LENGTH METHOD

A linear complementarity problem formulation combined with an arc-length method is presented for post-buckling analysis of geometrically non-linear structures with frictional contact constraints. The arc-length method with updated normal plane constraint is used to trace the equilibrium paths of the structures after limit points. Under the proportional loading assumption, the unknown load scale parameter used in the arc-length method is expressed in terms of contact forces, and eliminated to formulate as a linear complementarity problem. The unknown contact variables such as contact status and contact forces can be directly solved in this formulation without any ad hoc technique. Complicated non-linear buckling behaviours, such as snap-buckling, can be efficiently solved by the developed method, as shown by several buckling and post-buckling problems with frictional contact constraints.     

Parametric quadratic programming method for elastic contact fracture analysis

A solution procedure for elastic contact fracture mechanics has been proposed in this paper. The procedure is based on the quadratic programming and finite element method (FEM). In this paper, parametric quadratic programming method for two-dimensional contact mechanics analysis is applied to the crack problems involving the crack surfaces in frictional contact. Based on a linear complementary contact condition, the parametric variational principle and FEM, a linear complementary method is extended to analyze contact fracture mechanics. The near-tip fields are properly modeled in the analysis using special crack tip elements with quarter-point nodes. Stress intensity factor solutions are presented for some frictional contact fracture problems and are compared with known results where available.     

求解摩擦接触问题的IGA-B可微方程组方法

该文将等几何分析（Isogeometric analysis,简称IGA）和B可微方程组方法相结合,提出了求解弹性摩擦接触问题的IGA-B可微方程组方法。其中,接触边界的几何形状由非均匀有理B样条（NURBS）精确描述;接触条件则表示成B可微方程组的形式,可被严格满足,且在一定条件下求解该方程组的算法收敛性有理论保证。数值算例验证了该文方法用于求解弹性摩擦接触问题的有效性,计算精度较高,通过与ANSYS软件中的接触模型相比,计算自由度数可以大量节约。     

A finite element solution method for contact problems with friction

A new finite element solution method for the analysis of frictional contact problems is presented. The contact problem is solved by imposing geometric constraints on the pseudo equilibrium configuration, defined as a configuration at which the compatibility conditions are violated. The algorithm does not require any a priori knowledge of the pairs of contactor nodes or segments. The contact condition of sticking, slipping, rolling or tension release is determined from the relative magnitudes of the normal and tangential global nodal forces. Contact iterations are in general found to converge within one or two iterations. The analysis method is applied to selected problems to illustrate the applicability of the solution procedure.     

Three‐dimensional quasi‐static frictional contact by using second‐order cone linear complementarity problem

A new formulation is presented for the three‐dimensional incremental quasi‐static problems with unilateral frictional contact. Under the assumptions of small rotations and small strains, a second‐order cone linear complementarity problem is formulated, which consists of complementarity conditions defined by bilinear functions and second‐order cone constraints. The equilibrium configurations are obtained by using a combined smoothing and regularization method for the second‐order cone complementarity problem. Copyright © 2005 John Wiley & Sons, Ltd.     

Contact Analysis for Solids Based on Linearly Conforming Radial Point Interpolation Method

To simulate the contact nonlinearity in 2D solid problems, a contact analysis approach is formulated using incremental form of the subdomain parametric variational principle (SPVP). The formulation is based on a linearly conforming radial point interpolation method (LC-RPIM) using nodal integration technique. Contact interface equations are also presented using a modified Coulomb frictional contact model and discretized by contact point-pairs. In the present approach, the global discretized system equations are transformed into a standard linear complementarity problem (LCP) that can be solved readily using the Lemke method. The present approach can simulate various contact behaviors including bonding/debonding, contacting/departing, and sticking/slipping. An intensive numerical study is performed to validate the proposed method via comparison with the ABAQUS^® and to investigate the effects of the various parameters used in computations. These parameters include normal and tangential adhesions, frictional coefficient, nodal density, the dimension of local nodal support domain, nodal irregularity, shape parameters used in the radial basis function and the external load. The numerical results have demonstrated that the present approach is accurate and stable for contact analysis of 2D solids.     

A contact analysis approach based on linear complementarity formulation using smoothed finite element methods

Based on the subdomain parametric variational principle (SPVP), a contact analysis approach is formulated in the incremental form for 2D solid mechanics problems discretized using only triangular elements. The present approach is implemented for the newly developed node-based smoothed finite element method (NS-FEM), the edge-based smoothed finite element method (ES-FEM) as well as standard FEM models. In the approach, the contact interface equations are discretized by contact point-pairs using a modified Coulomb frictional contact model. For strictly imposing the contact constraints, the global discretized system equations are transformed into a standard linear complementarity problem (LCP), which can be readily solved using the Lemke method. This approach can simulate different contact behaviors including bonding/debonding, contacting/departing, and sticking/slipping. An intensive numerical study is conducted to investigate the effects of various parameters and validate the proposed method. The numerical results have demonstrated the validity and efficiency of the present contact analysis approach as well as the good performance of the ES-FEM method, which provides solutions of about 10 times better accuracy and higher convergence rate than the FEM and NS-FEM methods. The results also indicate that the NS-FEM provides upper-bound solutions in energy norm, relative to the fact that FEM provides lower-bound solutions.     

A semi‐implicit time‐stepping model for frictional compliant contact problems

In this paper, we formulate a semi‐implicit time‐stepping model for multibody mechanical systems with frictional, distributed compliant contacts. Employing a polyhedral pyramid model for the friction law and a distributed, linear, viscoelastic model for the contact, we obtain mixed linear complementarity formulations for the discrete‐time, compliant contact problem. We establish the existence and finite multiplicity of solutions, demonstrating that such solutions can be computed by Lemke's algorithm. In addition, we obtain limiting results of the model as the contact stiffness tends to infinity. The limit analysis elucidates the convergence of the dynamic models with compliance to the corresponding dynamic models with rigid contacts within the computational time‐stepping framework. Finally, we report numerical simulation results with an example of a planar mechanical system with a frictional contact that is modelled using a distributed, linear viscoelastic model and Coulomb's frictional law, verifying empirically that the solution trajectories converge to those obtained by the more traditional rigid‐body dynamic model. Copyright © 2004 John Wiley Sons, Ltd.