动力接触问题的间断galerkin时间有限元方法.pdfVIP

动力接触问题的间断Galerkin时间有限元方法 李建宇 沈晓阳 殷安琪 天津科技大学 机械工程学院，天津，300222 摘要：本文引入带时间参数的互补模型配合间断Galerkin时间有限元法研究动力接触问题的 数值方法，着重探讨了动力接触问题中冲击约束引起的空间域和时间域非光滑响应的模 拟。首先在空间域上给出了区别于静力接触互补模型的一类动力接触互补系统模型；然 后，引入间断Galerkin时间有限元方法求解在空间域上建立的互补系统模型。一个经典算例 验证了所提方法能够精确模拟由于接触冲击约束引起的接触-碰撞过程中的不连续响应。 关键词：动力接触；互补问题；非光滑响应；间断Galerkin时间有限元法 1 引言 工程实际应用的多数领域都会涉及接触现象，接触在结构和机械系统中 起着传递载荷和限制运动的作用，接触问题历来就是力学工作者研究的重要 课题。关于静态和准静态接触问题的数值方法研究已经比较成熟[1] ，然而， 动力情形下的接触分析要复杂的多。在数学模型上，动力接触边界体现为三 重非光滑性：第一，接触界面的不可穿透条件一般表述为不等式条件，法向 接触力与接触间隙之间需满足互补条件，这是空间上的不平滑性；第二，在 初始发生接触的瞬间时刻，接触边界上点的速度要发生突变，这是时间上的 不平滑性；第三，当接触面间存在摩擦时，例如Coulomb摩擦，摩擦力与接 触面切向相对速度之间的关系也是不平滑的。 [2] 在已有关于动力接触数值方法的研究中，一种策略如Hughes 等 ， [3] [4] Taylor和Papadopoulos 、Laursen等 给出的隐式算法，首先采用静力接触 分析中处理接触约束的办法获得一个预估解，然后通过某种补充条件修正给 出符合动力学性质的动力接触解。另一种常被使用办法是采用具有算法阻尼 的时间积分法，如耗散型Newmark法，广义- α法及单步Houbolt法等来剔除 [5] 由于接触条件不充分而导致的高频“假振荡” ，如Armero等 的工作，但这种 方法对于事先施加多大的算法阻尼才能完全剔除“假振荡”是未知的。显式算 [6] 法 采用罚函数法满足接触约束并采用中心差分法在时间域积分，罚函数法 具有某种正则化的作用，而中心差分法又仅用到位移和速度的信息，这些措 1 施在一定程度上放松了对动力接触问题解的光滑性的要求，从而缓减了接触 点速度的不连续性导致的“假振荡”现象。但是，罚因子的恰当选取至今仍是 一个公认的难题。 为了研究动力接触分析的更高效算法，本文继承了互补模型在静力接触 分析中精确建模的优点，在空间域内提出了动力接触分析的互补系统模型， 而在时间域内，考虑到接触约束引起动力响应的非光滑性，在动力接触分析 中引入不连续时间积分策略，并建立一类间断Galerkin时间有限元列式。 2 互补型动力接触界面约束条件 考虑如图1所示的两体接触问题，接触体记为 Ω Ω Ω , Ω 。假设 A , B 中接触边界上的点是一一对应 A B 的 ， 对 于 每 对 接 触 点 ， 定 义 其 位 移 向 量 T 1 2 3 1 2 3 T uc (ucA ,ucB ) (ucA ,ucA ,ucA ,ucB ,ucB ,ucB ) 。引入广义 (u1 ,u2 ,u3 ,u1 ,u2 ,u3 ) 坐标系 cA cA cA cB cB cB ，则广义坐标系中的一 组值就代表了接触点对的一个可能状态，各种可能 状态的全体构成接触点对的构型空间Q 。如图1所 示，接触边界上的法向非穿透条件可表达为关于接触点位移u 的不等式：