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

第 19 卷 　第 2 期 应 用 力 学 学 报 Vol. 19 　No. 2 2002 年 6 月 CHINESE JOURNAL OF APPL IED MECHANICS J un. 2002 文章编号 (2002) 大变形中摩擦接触问题的数值模拟及应用 蔡中义　李明哲　陈庆敏 ( 吉林大学辊锻研究所　长春　130025) 摘要 :大变形的摩擦接触是复杂的非线性问题 ,本文介绍了一种处理摩擦接触问题的数值方法。采 用接触单元技术模拟接触界面 ,基于弹塑性理论形式的非经典 Coulomb 摩擦定律及罚函数方法建 立了摩擦接触的增量本构关系。结合大变形的增量分析格式给出了积分摩擦接触本构方程的回映 方法。这种处理摩擦接触问题的方法计算简单、使用方便。给出的计算实例及应用实例说明了方 法的精度与稳定性。 关键词 :摩擦 ;接触 ;Coulomb 定律 ;大变形 ;有限元 中图分类号: TG302 　　　文献标识码 : 　A 微观位移 ,因此也可在一定程度上克服经典摩擦定 1 　引　言 律的缺陷。 从纯数学的角度来讲 ,求解摩擦接触问题实际 摩擦接触界面的处理是大变形数值分析中一个 上就是求解受变分不等式约束控制的物理方程组的 十分重要的问题。目前摩擦问题有限元模拟使用较 问题,约束边界就是接触区域。这类问题通常采用 多的仍是 Coulomb 在 1781 年提出的经典干摩擦定 Lagrange 乘 子 法[5 ,6 ] 及 罚 函 数 方 法 ( Penalty 律 ,即处于接触的两物体之间的相对滑移 ,只有当接 method) [7 ,8 ]等技术进行求解。罚函数方法不增加 触面的切向力达到临界值时才能发生 , 临界值的大 结构自由度 ,数学处理也比较方便。本文将根据弹 小与法向接触压力成比例 , 比例常数即为摩擦系数。 塑性理论形式的非经典摩擦理论 ,采用罚函数技术 , 当切向力小于临界值时 ,处于纯粘着状态 ,接触面的 结合大变形问题的增量分析过程 ,给出处理摩擦接 相对滑移量为零。一旦产生滑移 ,滑移量的大小不 触问题的数值模拟方法。 影响摩擦力。 事实上 ,对于任何小于临界值的摩擦力都可产 2 　摩擦接触约束的描述 生一定的微小相对位移。为克服经典 Coulomb 摩 擦定律的缺陷 ,人们提出了非经典摩擦理论。比较 2 1 　Coulomb 摩擦的 KuhnTucker 约束 有代表性的是 Oden 和 Pires[1 ,2 ] 提出的以切向相对 由经典Coulomb 摩擦定律定义滑动函数 Φ,则 滑移 为 函 数 的 摩 擦 理 论 以 及 Frericksson[3 ] 、 接触点处切向约束表示成 KuhnTucker 条件为 Curnier [4 ]等提出的类似于弹塑性理论形式的摩擦 Φ( μ ) μ t T , , tN ∶= ‖t T ‖- | tN | ≤0 理论。前者能够反映摩擦问题的非线性特征及非局 Φ d g T + d λ = 0 部特征 ,理论形式比较完备 ,但所涉及的参数不易确