接触-碰撞.pptVIP

小结 举例 粗糙刚体离散的潜在影响 问题 单纯主从接触，平衡主从接触； 接触条件的运动学和动力学的表示，及其行为特征，单一接触条件的力学意义； 接触摩擦模型； Lagrange乘子和罚函数方法的物理意义； 接触问题弱形式和有限元格式； 显式程序实现。 接触界面应变不连续条件 相对于裂纹的强不连续，称其为弱不连续，如夹杂、损伤、剪切带等。 * * * 罚是拉回侵彻量的弹簧力 * 注意罚是等式，乘子是不等式 * * r是残数 ? -lamda * L乘子为刚塑性-罚为理想弹塑性 * L乘子为刚塑性-罚为理想弹塑性 * 对于Coulomb类型性能的屈服函数，对应于Coulomb摩擦条件： 在二维情况下，屈服函数包括以斜率为 的两条斜线。 在三维情况下，屈服函数是一个圆锥。 Mohr-Coulomb本构模型 界面本构方程 4 弱形式 对于Lagrangian网格，建立动量方程和接触界面条件的弱形式。当接触表面是作为L格式处理时，这一形式也适用于ALE网格。为了简单，从无摩擦接触开始，后面介绍切向面力的处理。将下面的公式限制在面力或者位移边界的情况，进而描述相应的所有面力或者速度分量。 接触表面既不是面力也不是位移边界。物体A的全部边界为 注意到 在接触问题中，奇异性发生在边界，如赫兹接触(1881)。然而与断裂力学不同的是，奇异性没有显示任何的工程意义，因为表面粗糙度的存在抵消了在应力中出现的更接近奇异性的行为。 自然变分原理是对原物理问题的微分方程和边界条件建立对应的泛函，再使泛函取驻值得到问题的解答，但是其未知场函数需要满足一定的附加条件。而广义变分原理（或称约束变分方程）不需要事先满足附加条件，采用Lagrange乘子法和罚函数法将附加条件引入泛函，重新构造一个修正泛函，将问题转化为求修正泛函的驻值。称为无附加条件的变分原理。 对于罚函数方法，将罚参数取正值，对修正泛函得到的近似解只是近似地满足附加条件，罚参数值越大，附加条件满足程度就越好，即鞍形条件。而在实际计算中，罚函数只能取有限值，所以利用罚函数求解只能得到近似解。 4 弱形式 Lagrange乘子弱形式 强加接触约束的通常方法是借助于Lagrange乘子。 令Lagrange乘子试函数为 相应的变分函数为 弱形式为 其中(1) 侵彻率变分为 这个弱形式不等式是等价于(2)动量方程、面力边界条件、内部连续条件（广义的动量平衡）和(1)接触界面条件：不可侵彻性、法向面力的动量平衡和无摩擦条件。Lagrange乘子场仅要求为C-1连续，因为它的导数并不出现在弱形式中。要求法向界面力是压力，这是对Lagrange乘子在试空间的一种限制。 和(2) 虚功率为 Lagrange乘子弱形式 比较Hu-Washizu变分原理，借助于Lagrange乘子，上面的方法是在弱形式中附加约束的标准方法。与Hu-W形式的唯一区别是约束为一个不等式。利用高斯积分关系，通过积分证明 从弱不等式可以推导出强形式，但是，必须考虑关于变分函数的不等式。一旦变分函数不受约束时，对于与变分函数相乘项的符号则没有限制，并且由密度原理，该项必须为零。 从上式中的前两个积分得到 即在物体A和B上，满足动量方程和自然边界条件。 Lagrange乘子强形式 从上式联立消去 得到关于法向面力的动量平衡条件： Lagrange乘子强形式 在接触表面被积函数的所有项中，除最后一项外，变分函数是任意的，因此获得了等式， 由上式被积函数最后一项的变分函数为负值 推论 必然是非正的，即弱不等式表示不可侵彻条件为 拉格朗日乘子法原理 拉格朗日乘子法的原理是首先不考虑接触，对物体的加速度、速度和位移进行更新计算，得到试探解答。然后，依据试探解进行接触判断，如果没有发生接触，试探解就是真实解；否则，接触发生并计算接触力，接触力的大小可以从接触的非穿透条件得到（一致性条件）。得到接触力后，修正试探解，使得该时间步最后得到的位移和速度等变量满足接触的非穿透条件。该算法的最大优点是总是能够保证接触界面节点满足非穿透条件。 在罚函数方法中，最终缝隙值将不是零，而是一个小量，取决于所选择的参数值。罚函数方法的优点是可以根据判别,能够给出一个可接受结果的实际需要的参数值，进行适当的调整。实际上，在一个复杂的问题中，这不是简单的工作，特别是对于涉及梁、板、或壳单元与实体单元之间的接触问题。 4 弱形式 罚函数方法 一种能够避免Lagrange乘子求解困难的方法是罚函数。罚方法是在接触节点之间附加一个刚度为 β的非线性弹簧。 罚函数方法原理 罚方法的原理是在每一时间步先