冲击载荷作用下结构动态响应分析的罚函数有限元法

基于罚函数法，根据接触－冲击问题的边界条件，运用双线性参数法推导出用有限元模拟三维接触－冲击问题的具体算法，并对杆和梁在冲击载荷下的动态响应进行了仿真分析，取得较好计算结果。     

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

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

正交各向异性弹塑性摩擦接触问题的数值求解

采用正交各向异性摩擦定律对三维弹塑性摩擦接触问题进行分析,基于参变量变分原理,经过有限元离散,将问题化为线性互补问题模型,之后给出一个求解互补问题的非内点光滑化算法．对三维接触问题,滑动方向的确定一直是个难点,为此,该文采用作者提出的组合规划法和迭代法对各向异性摩擦本构模型进行分析,数值结果说明了模型与算法的正确性。     

基于有限元法的轮轨蠕滑特性研究

轮轨滚动接触蠕滑率／力理论是轮轨相互作用一系列问题研究的基础．现有的几种计算蠕滑力的理论模型均建立在Hertz接触条件和半空间假设的基础上，已经无法完成对复杂的接触问题的进一步研究．使用有限元参数二次规划法来求解轮轨的三维弹性／弹塑性接触问题，得出了在不同的轴重、牵引力矩、摩擦系数、踏面形状、横向力条件下的轮轨接触力．提出用轮轨接触的轮周位移计算蠕滑率的新方法，并对在各种参数下所得到的蠕滑率进行了分析比较．     

三维接触问题的非光滑算法

给出了一种非光滑算法直接用于求解三维摩擦接触问题的不可微非线性互补模型,不再对模型进行光滑化处理,使算法更加简单。文中对非光滑算法的收敛性给出了严格的数学证明,数值实验表明该算法列式简单,但与光滑化算法同样有效。     

厚钣金件压弯翻边与回弹的数值研究

在独立研制的板材冲压成形性分析软件KMAS中，组入可有效模拟厚钣金件压弯成形与回弹变形的三维退化曲壳单元模型，简述了考虑凹，凸模间隔可变的基于有限元网格双面接触算法以及回弹的处理技术，数值研究了铁路客车牵引架厚钣金件压弯翻边与回弹变形过程，并将压弯过程模拟结果与国外基于单面接触算法的商品化软件的模拟结果以有实验结果进行比较，最后，给出了该类问题回弹量的预示方法。     

一种新的三维动力接触约束单元

本文应用广义变分原理，利用拉氏乘子法和罚函数法计入接触约束条件修正，建立了分析三维动力接触问题的一般有限元分析的理论模式；推导了处理这类问题的一种新的接触约束单元的接触刚度矩阵；采用增量求解模型解，研编了实施程序３ＤＤＣＦ；处理了方程中的病态问题、碰撞及能量释放条件问题。进行了算例考核，还首次给出了一个三维动力接触问题的算例结果。     

气体钻井中接触非线性问题的数值算法研究

根据动量守恒和钻柱力学的基本理论,建立了气体钻井三维弯曲井眼钻柱动力学模型和钻柱与井壁接触碰撞模型;提出了一种求解钻柱与井壁接触非线性有限元方程的元胞自动机算法。该方法具有编程简单、收敛速度快等优点,以实际工况为例进行仿真模拟,选取有代表性节点的仿真计算结果并加以分析。仿真结果表明:该方法较好地解决了三维弯曲井眼中钻柱与井壁接触碰撞的复杂边界条件问题,在解决气体钻井钻柱与井壁接触非线性问题的研究中具有良好的应用前景。     

岩质圆形隧洞围岩应力场弹塑性新解

针对动态接触问题的有限元并行计算，提出了一种新的接触算法. 新算法引入局部拉氏 乘子技术来计算接触力. 由于同时考虑了无穿透的接触约束条件和相邻接触对的相互影响， 较之广泛使用的罚参数法，新算法使接触约束条件和系统平衡方程得到更充分的满足. 虽然 为提高接触计算精度而在局部采用了迭代技术，但算法仍然具有较高的效率，且与显式时间 积分方案完全相容. 此外，通过构造专门的区域分解方案，实现了将现有为串行程序开发的 搜索算法平滑移植到并行环境的目标. 数值算例表明，所提出的接触算法具有很好的并行性， 在保证了接触问题并行计算精度的同时，取得了满意的并行效率.     

沉管隧道地震响应分析中的三维接触模型与算法研究

针对土体与埋置在土体中的管段以及沉管隧道各个柔性接头处不同构件之间的动力相互作用,采用基于超级计算机的三维接触模型对上海某沉管隧道进行了地震响应分析,并进行了三维接触算法的研究。整个模型单元数超过了120万,节点数超过了100万。研究证明,应用文中提出的接触模型和算法可将接触处理的CPU时间大大缩短。地震波载荷选用调幅后的唐山地震波,应用通用非线性动力分析有限元程序LS-DYNA进行求解,计算结果可为隧道的抗震设计提供有力依据。     

解三维摩擦接触问题的一个二阶锥线性互补法

针对三维摩擦接触问题的求解,给出了一种基于参变量变分原理的二阶锥线性互补法. 首先,基于三维Coulomb摩擦锥在数学表述上属于二阶锥的事实,利用二阶锥规划对偶理论,建立了三维Coulomb摩擦接触条件的参变量二阶锥线性互补模型,它是二维Coulomb摩擦接触条件参变量线性互补模型在三维情形下的自然推广;随后,利用参变量变分原理与有限元方法,建立了求解三维摩擦接触问题的二阶锥线性互补法. 较之于将三维Coulomb摩擦锥进行显式线性化的线性互补法,该方法无需对三维Coulomb摩擦锥进行线性化,因而在保证精度的前提下所解问题的规模要小很多. 最后通过算例展示了该方法的特点.      

微/纳尺度粘着滑动接触的多尺度分析

采用分子动力学与有限元耦合的多尺度方法,求解二维刚性圆柱表面压头与弹性平面的微/纳尺度粘着滑动接触问题,通过与全分子动力学模拟结果的比较验证了多尺度方法的有效性。对压头半径、滑动速度、下压深度以及是否考虑粘着效应等对滑动接触性能的影响进行了全面研究,通过不同条件下摩擦力及接触力分布的比较,揭示了上述各参数对粘着滑动接触...     

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

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

三维弹性接触问题的接触面单元法

基于虚功原理，推导了三维接触面单元的刚度矩阵，并引进预留单元的概念，避免了接触过程中由于接触面变化节点和单元需要重新编号的麻烦．文中算例验证了计算方法的有效性     

二维弹性接触问题的接触面单元法

本文基于虚功原理，推导了二维接触面单元的刚度矩阵，并引进预留单元的概念，避免了接触过程中由于接触面变化，节点和单元需要重新编号的麻烦。采用位移和应力联合控制的增量法控制加载过程，文中给出考题验证计算方法的有效性，并给出本文方法在汽检轮机中根与轮缘接触问题中的应用实例。     

表面形貌变形对塑性成形滑动接触界面摩擦的影响

为了更好地理解塑性成形滑动接触界面的摩擦行为,构建了一种新型的摩擦试验装置,运用表面纹理化技术制备了两类表面形貌的1050铝材试件,在不同的接触压力和滑动速度条件下进行一系列拉伸摩擦试验.对试验前后试件三维表面形貌进行了测量;提取真实接触面积比、封闭空体面积比和开放空体面积比等三维表面参数,来描述试件表面形貌的变化.试验发现:摩擦系数随名义接触压力和滑动速度增加而逐渐减小;试件初始表面形貌对摩擦有明显的影响;试件表面形貌和参数随接触条件出现了规律性变化.基于机械流变模型的分析表明:随着试件表面形貌变形,不同的机理决定界面摩擦行为,摩擦系数对名义接触压力和滑动速度的依赖性可分别归因于微观塑性流体动压润滑效应和入口区流体动压牵引效应.     

THREE DIMENSIONAL DYNAMIC ELASTICPLASTIC ANALYSIS OF THE INTERACTION OF A SET OF PERCUSSION BITS WITH ROCK MASSES

A three dimensional dynamic elastic-plastic finite element analysis of the interaction ofpercussion bits with rock roasses is presented.The dynamic contact problem was investigated with thenumerical method developed for high velocity impact calculations.The numerical analysis was carried out for three simultaneously loaded percussion bits arrangedon a straight line on the rock surface,the loading time was kept constant throughout this work.Thenumerical calculation for various strengths of the dynamic load on the percussion bits shows that a crit-ical minimum dynamic load exists for the generation of plastic deformations in the rock.For this criti-cal value of the load the maximum normal stress of the longitudinal wave in the percussion bit exceedsthe yield stress of the rock mass.     

PROGRAM-PATTERN MULTIPOLE BOUNDARY ELEMENT METHOD FOR FRICTIONAL CONTACT

A mathematical program is proposed for the highly nonlinear problem involving frictional contact. A program-pattern using the fast multipole boundary element method (FMBEM) is given for 3-D elastic contact with friction to replace the Monte Carlo method. A newoptimized generalized minimal residual (GMRES) algorithm is presented. Numerical examples demonstrate the validity of the program-pattern optimization model for node-to-surface contact with friction. The GMRES algorithm greatly improves the computational efficiency.     

大变形中摩擦接触问题的数值模拟及应用

大变形的摩擦接触是复杂的非线性问题 ,本文介绍了一种处理摩擦接触问题的数值方法。采用接触单元技术模拟接触界面 ,基于弹塑性理论形式的非经典Coulomb摩擦定律及罚函数方法建立了摩擦接触的增量本构关系。结合大变形的增量分析格式给出了积分摩擦接触本构方程的回映方法。这种处理摩擦接触问题的方法计算简单、使用方便。给出的计算实例及应用实例说明了方法的精度与稳定性     

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

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