大位移接触摩擦单元在工程设计中应用.pdf

t／L 维普资讯 ， 华 东 交 通 大 学 学 报 Vo1．】lNo 1 1994年 3月 JournalofEastChinajiaotongUniversity Mar． 1994 帕 大位移接触摩擦单元在工程设计 中的应用 雷 晓 燕 (建筑 I程 系) 摘 要 一 种新 的大位移接触摩擦单元采用点一面接触方式 ，能够模拟大位移情况下两物 体间的张开 、闭合和摩擦滑动 。单元的基本未知量包括带点位移和接触应力。单元 的几何和静力约束是作为补充方程包含在接触摩擦单元的剐度方程 中。基于增量丹； 式的有 限元法 ．本文建立 7大位 秽接触 问题 的等效单元刚度矩 阵和荷载 向量．使之 与常规 的有限元一样 ，可在单元水平上进扦刚度和荷我 的遥加 。圆而具有算法简单 的优 点。 ， 关键词：盟 堡墼壁堡娅 工每 计l 0 引 言 J 接触摩擦是工程 中最常见和最重要的问题之一。许多研究人员致 力于这一领域的工作 早 期的研究集中在小位移问题方面．采用的是点一点接触模型，如节理单元、接触单元．罚函数 法和拉格朗日法 。点一点接触模型算法简单 ．但在应用上有局限性 ，不适宜大变形情况 。 近年来，随着计算机技术的发展 ，分析各种复杂大位移接触摩擦问题 已成为可能。目前． 研究人员的重点是寻求有效算法。这方面早期的工作要数Rothert 和HaiIquist ‘。自那 以后． 该领域的研究十分活跃 ．已有一系列的研究成果问世” 。由于在大位移接触问题中 ．点一点接 触模型不再适用，需要考虑接触运动。大位移问题中的接触条件是很复杂的．以致于出现丁 各种各样的理论 。目前的理论大多过于复杂．不便于实际应用。本文在文 [1]的基础上．采 用点一面接触摸型．提出了一种新的大位移接触摩擦单元。运用这一模型求解问题时 ．首先使 接触界面上的单元网格点点对应 随着荷载的增加．变形越来越大。在这一过程中．点一点接 触变成了点一面接触，单元必须采用更新的位形 。 接触摩擦单元定义于图 l。单元未知量包括节点位移和接触应力。单元的几何和静力约束 条件是作为补充方程包含在接触摩擦单元的刚度方程中 由此建立丁大位移接触 问题的等效 单元刚度矩阵和荷载 向量．可 以用常规的方法将接触摩攘单元的贡献迭加到总体有限元方程 中。 本文于 I993年 1l目291敢 到 一 — j—一 维普资讯 2 华 东 交 通 大 学 学 报 1 接触摩擦单元理论 在大位移问题中，物体的平衡方程必须建 立在当前的位形中。为此 ，在数值计算时需要 采用增量算法 。 考察两物体 A和B 的接触，A、B区域采 用普通有限元，接触界面采用接触元。接触初 始时，节点对应 。整个 问题酌描述采用笛卡尔 4 坐标系 z、 ，局部坐标系 n，s则用于描述接 触摩擦元 。 在接触摩擦分析中，随着外荷载的逐渐增 加 ，变形也越来越大 ，接触摩擦单元节 点对 由 初始时的点一点接触变成了点一面接触，如图2 圈l接触摩擦单元 所示 。这时 ，单元的 位形 1—2—3—4—5—6必须用新的位形 1—2—3—4一5一6更新。 在更新的位形 中，运用虚位移原理则有 妇 F 一 1( )互dr。 (1)