接触问题的自适应无网格伽辽金方法

针对接触问题提出了基于应变能梯度的自适应无网格伽辽金方法。通过对自适应无网格伽辽金方法的原理讨论,给出了基于背景网格的误差估计方法,并对自适应求解流程进行了阐述。研究了无网格节点加密技术和自适应搜索半径等关键技术。编程求解了圆柱体接触问题,并与理论解和均匀密化数值解进行比较,给出了接触计算中合理的参数范围,对真实粗糙表面接触问题进行了计算,计算结果表明,在相当计算精度下,自适应加密计算耗时仅为同等整体加密计算的18.6%。     

二维自适应粗糙表面弹塑性接触模型

提出一种基于自适应粗糙表面描述的弹塑性接触模型。根据某一给定阈值,去除粗糙表面上对接触力学行为影响较小的结点以减少计算时间。采用在计算区域划分及结点排布上较为灵活性的无网格伽辽金—有限元耦合方法求解自适应粗糙表面弹塑性接触模型。算例表明,当接触压力相对误差约为5%时,自适应接触模型的计算时间相对于非自适应接触计算可减少约50%。分析不同阈值对接触压力分布、接触间隙及接触体弹塑性应力分布的影响。     

刚塑性无网格方法中体积闭锁问题的缓解算法

将无网格伽辽金方法引入金属塑性成形过程分析,在给出基于刚塑性理论的无网格伽辽金方法及其基本理论方程的基础上,重点研究应用无网格方法分析金属塑性成形问题时所面临体积闭锁问题的缓解算法,采用体积应变率映射方法,对能量速率泛函中的体积应变率进行修正处理,通过将速度场计算的体积应变率映射到低维空间,降低独立约束方程数目,从而建立体积闭锁现象的缓解算法。应用所建立的方法对典型锻造过程进行了无网格分析,结果表明,该算法能够有效解决体积闭锁问题。     

金属塑性成形过程无网格伽辽金法数值模拟技术研究

将无网格伽辽金法与刚(粘)塑性流动理论相结合,对无网格伽辽金法在金属塑性成形过程数值模拟中的应用技术进行了研究,分别采用混合变换法和罚函数法实现速度边界条件和体积不变条件,提出了基于刚(粘)塑性力学流动理论的无网格伽辽金法,给出了金属塑性成形无网格伽辽金法数值模拟的关键算法,拓展了无网格方法在金属成形领域的应用范围。典型算例的数值计算结果表明了方法的可行性。     

基于刚塑性流动理论的无网格伽辽金方法

为了对金属体积成形过程进行数值模拟，给出基于刚塑性流动理论的无网格伽辽金方法。由于移动最小二乘近似的非插值特征给无网格伽辽金方法的应用带来了一定的困难，将再生核质点法中的完全变换法与无网格伽辽金方法相结合，通过对速度场的移动最小二乘近似进行修正，实现速度边界条件的精确施加；采用罚函数法引入体积不可压缩条件，运用刚塑性材料的不完全广义变分原理，给出基于刚塑性流动理论的无网格伽辽金方法的离散控制方程以及收敛判据和收敛控制方法。数值计算结果表明方法的可行性及有效性。     

无网格伽辽金方法在线弹性断裂力学中的应用研究

通过对移动最小二乘形函数进行局部修正,将混合变换法应用于无网格伽辽金方法,给出分析线弹性断裂力学问题的有效的无网格伽辽金方法。这一方法克服了无网格伽辽金方法中常用的拉格朗日乘子法和罚函数法的缺点,实现了本质边界条件在节点处的精确施加。运用线弹性断裂力学理论,采用基于t-分布的新型权函数和部分扩展基函数,对有限板单边裂纹的应力强度因子和拉剪复合型裂纹的扩展进行分析。由于该方法仅需节点信息,而不需要节点的连接信息,从而避免了有限元方法中的网格重构,大大简化了裂纹扩展的分析过程。数值计算结果表明了方法的有效性。     

基于一致径向基函数的金属塑性成形过程无网格方法分析

径向基函数满足δ条件,可以直接施加本质边界条件.通过耦合多项式基函数可以精确重构任意多项式.将无网格伽辽金方法引入到塑性成形过程模拟,结合刚塑性材料的假设,提出基于一致径向基函数的无网格刚塑性伽辽金方法.采用体积应变率映射法处理体积不可压缩条件;采用反正切摩擦模型在局部坐标系下施加摩擦力边界条件,实现任意边界形状的塑性成形过程模拟.对典型塑性成形过程进行无网格方法分析,并与刚塑性有限元方法及试验数据进行比较,验证了建立方法的可行性.     

刚(粘)塑性无网格伽辽金方法及其关键技术研究与应用

在国家杰出青年科学基金资助项目《塑性加工工艺及设备》(No.50425517)，以及国家自然科学基金资助项目《金属体积成形过程的刚(粘)塑性无网格伽辽金方法数值模拟理论及其关键技术研究》(No.50575125)支持下，开展刚(粘)塑性无网格伽辽金方法及其关键技术与应用的研究报告，对刚(粘)塑性无网格伽辽金方法的基础理论、数学模型建立方法、关键处理技术在金属塑性成形过程中应用研究的成果。 将无网格伽辽金方法引入塑性成形过程模拟，提出基于刚(粘)塑性理论的无网格伽辽金方法，推导刚度矩阵方程和求解列式。利用变换法施加本质边界条件，采用反正切摩擦模型描述摩擦力边界条件。对于模具边界任意的塑性成形过程，在局部坐标系下施加摩擦力边界条件，给出局部坐标系和整体坐标系的变换矩阵，解决了模具形状任意的二维塑性成形问题摩擦力边界条件的施加问题。采用直接迭代法获得初始速度场，利用Newton Raphson迭代方法求解刚度方程，给出模拟等温塑性成形问题的分析步骤。对于中高温条件下的塑性成形过程，推导出刚(粘)塑性无网格伽辽金方法热力耦合分析模型，给出热力     

二维无网格伽辽金-有限元耦合方法的研究

论述无网格伽辽金方法（element-free Galerkin method,EFG）的基本原理－－移动最小二乘原理（moving leastsquare,MLS）和EFG方法的位移近似函数，给出变分方程及离散方程。对无网格伽辽金－有限元耦合（EFG－FE coupling）方法进行详细阐述。编制相应程序，通过算例表明拉格朗日乘子对强制边界条件的作用及无网格伽辽金方法在权函数参数变化时的收敛特性，论证无网格伽辽金－有限元耦合的有效性。     

面向对象的无网格伽辽金法

采用Schmidt法对无网格伽辽金法(EFGM)形函数的基函数进行正交化,避免了形函数中的矩阵求逆运算,同时也减少了运算量并提高计算精度。对EFGM中的权函数、节点影响半径等关键问题进行了改进,并利用对象设计思想将无网格伽辽金法的对象抽象为独立的数据类,研究了无网格伽辽金法数据类的设计方法和组织方式并编制了相应的计算程序。数值算例结果表明,理论正确、方法有效。     

基于积分弱形式的无网格-有限元法耦合方法研究

无网格伽辽金法是基于移动最小二乘法的一种求解偏微分方程非常有效的方法,但是对于一些规模较大的问题,其计算费用较高,将无网格法与有限元耦合可以极大提高求解效率.文中根据D. Hegen提出的基于积分弱形式耦合方法,进一步明确处理耦合区域部分的数学、物理意义, 有效地引入一个综合控制因子反映权函数的紧支性, 并经算例验证方法的正确性及控制因子的有效性.     

非均质材料无网格Galerkin法

提出了采用无网格Galerkin方法来分析非均质材料。无网格伽辽金法（EFGM）是近些年发展起来的一种数值算法，它采用移动的最小二乘法构造形函数，从能量泛函的弱变分形式中得到控制方程，该法只需节点信息，不需将节点连成单元。在积分网格中，取高斯点的材料参数来模拟材料特性的变化。算例的结果显示用该法来分析非均质材料得到的结果精度高，表明该方法在非均质材料领域有着广阔的应用前景。     

B样条小波基自适应无网格迦辽金法应用于刚塑性成形模拟

无网格方法适合大变形计算和自适应分析,是进行刚塑性成形过程仿真的比较有潜力的方法。利用样条小波的多分辨分析特性,基于刚塑性变形的特点,通过实现刚塑性场量的两尺度分解技术,将场变量分解成低尺度成分和高尺度成分,利用高尺度成分判断求解高梯度区域;利用基于局部Delaunay三角化的细分方案实现高梯度区域节点的自适应细分,最终实现一个应用于刚塑性成形仿真的样条小波基自适应无网格迦辽金法,并用实例验证算法的稳定性和有效性。     

Adaptive multiscale method for two-dimensional nanoscale adhesive contacts

There are two separate traditional approaches to model contact problems: continuum and atomistic theory. Continuum theory is successfully used in many domains, but when the scale of the model comes to nanometer, continuum approximation meets challenges. Atomistic theory can catch the detailed behaviors of an individual atom by using molecular dynamics (MD) or quantum mechanics, although accurately, it is usually time-consuming. A multiscale method coupled MD and finite element (FE) is presented. To mesh the FE region automatically, an adaptive method based on the strain energy gradient is introduced to the multiscale method to constitute an adaptive multiscale method. Utilizing the proposed method, adhesive contacts between a rigid cylinder and an elastic substrate are studied, and the results are compared with full MD simulations. The process of FE meshes refinement shows that adaptive multiscale method can make FE mesh generation more flexible. Comparison of the displacements of boundary atoms in the overlap region with the results from full MD simulations indicates that adaptive multiscale method can transfer displacements effectively. Displacements of atoms and FE nodes on the center line of the multiscale model agree well with that of atoms in full MD simulations, which shows the continuity in the overlap region. Furthermore, the Von Mises stress contours and contact force distributions in the contact region are almost same as full MD simulations. The method presented combines multiscale method and adaptive technique, and can provide a more effective way to multiscale method and to the investigation on nanoscale contact problems.     

Analysis of the free vibration of rectangular plates with central cut-outs using the discrete Ritz method

We present an analysis of the free vibration of plates with internal discontinuities due to central cut-outs. A numerical formulation for a basic L-shaped element which is divided into appropriate sub-domains that are dependent upon the location of the cut-out is used as the basic building element. Trial functions formed to satisfy certain boundary conditions are employed to define the transverse deflection of each sub-domain. Mathematical treatments in terms of the continuities in displacement, slope, moment, and higher derivatives between the adjacent sub-domains are enforced at the interconnecting edges. The energy functional results, from the proper assembly of the coupled strain and kinetic energy contributions of each sub-domain, are minimized via the Ritz procedure to extract the vibration frequencies and mode shapes of the plates. The procedures are demonstrated by considering plates with central cut-outs that are subjected to two types of boundary conditions.     

An element-free Galerkin forward solver for the complete-electrode model in electrical impedance tomography

An element-free Galerkin method (EFGM) is used for solving forward problem based on the complete electrode model (CEM) in electrical impedance tomography (EIT). The EFGM requires only nodal data and has the ability of providing mesh-independent solutions because no element connectivity is needed to be used in this method. However, direct imposition of Dirichlet boundary conditions for the EFGM is difficult because the shape functions employed in this method do not have the property of Dirac delta function. Solving the EIT forward problem based on the CEM by the EFGM, the effects of electrodes and contact impedances are taken into account and the complete solution of equations is provided without imposing Dirichlet boundary conditions. The numerical results are validated with experimental results obtained from a 2D circular homogeneous phantom, and the performance of the EFGM compared with the finite element method is also illustrated. Moreover, results obtained from the EFGM are presented for an inhomogeneous numerical phantom.     

Behaviour law for cutting process

To optimize the cutting conditions in machining, it is necessary to quantify energy parameters involved in the cutting process. These energy parameters are based on a hypothesis regarding the displacement of the material during the chip formation and the behaviour law of the material. However, a more refined model must be used to take into account complex phenomena in the cutting process. The strain gradient theory appears to be best in modelling rotation phenomena inside material. In order to solve this type of problem, the bases of the strain gradient theory are presented in this study.     

Integrated adaptive element free Galerkin model for elastoplastic contact of engineering surfaces

Abstract

This paper presents an h-adaptive element free Galerkin (EFG) model and its formulation based on stain energy gradient for solving elastoplastic contact problems. In this model, the element free Galerkin finite element (EFG–FE) coupling method uses the initial stiffness method, and an error estimation approach based on strain energy gradient and a local adaptive refinement strategy are combined. Two-dimensional elastic contact problem between a rigid cylinder and an elastic plane is analysed to validate the adaptive elastoplastic EFG model. The results indicate that the adaptively refined solutions are accurate as compared to the Hertz theoretical solution or the uniformly refined solutions, while the cost of CPU time used by the adaptive model is less than that by the uniformly refined model under the same condition. Furthermore, the elastoplastic contacts involving a rough surface of elastic perfectly plastic and elastoplastic properties are investigated.     

基于应变模态的模态应变能损伤识别方法

传统模态应变能计算需要完备模态振型信息,而模态振型信息中存在转角自由度难以准确获取的问题,为解决该问题,开展基于应变模态的模态应变能损伤识别研究,实现了结构损伤的定量识别。首先,通过基于应变与位移之间的联系,推导出应变模态与位移模态之间的转换矩阵;其次,利用应变模态代替位移模态计算单元模态应变能,建立基于灵敏度分析的损伤识别方程组;最后,根据奇异值截断法求解该方程组识别结构损伤。以一两端固支梁结构为对象,开展数值仿真和实验研究。结果表明,该方法可以有效识别出结构的损伤位置和损伤程度,相对于基于振型扩充的模态应变能损伤识别方法,具有更好精度和抗噪性能。