清华大学计算固体力学第八次单元技术.pptVIP

4 Q4和体积自锁 Q4－4节点四边形单元 考虑单元1，仅可能非零的节点速度是在3点 是任意值 单元 1 的所有其它节点速度必须为零以满足边界条件。 对于一个任意的运动，其膨胀率为 节点3的速度给出 ， 所以，膨胀率的常数项是非零的，除非 因此，一个等体积运动需要膨胀率为零，即 4 Q4和体积自锁 Q4－4节点四边形单元 当 其中 只有沿着直线 上式才为零！ 尽管单元的运动是一个常数体积运动，除了在该直线上，膨胀率是处处非零。为了满足在整个单元中的等体积运动，? 必须为零，并且节点3不能移动。 如果节点3不能移动，在单元2的左侧，则由节点2和3提供了刚性边界，并且对于单元2，重复这些讨论可以证明节点6是不能移动的。这一讨论则可以对网格中的所有单元重复，以证明所有节点的速度必须为零。即有限元模型自锁。这一讨论也适用于歪斜单元。 4 Q4和体积自锁 Q4-4节点四边形单元 另一种检验的方法是考虑单元的双线性速度场 膨胀率给出为 通过在整个单元域上积分膨胀率，计算一个单元面积的变化 线性分量 双线性项积分为零，其导数正交于常数场 证明对于任意的等体积速度场， 保持单元面积为常数的运动的膨胀率则是 是必要的 设上式结果为零，以反映等体积运动， 4 Q4和体积自锁 Q4－4节点四边形单元 尽管所有