约束条件(Constraint) — CAE用语解说
约束条件(Constraint)是什么
约束条件就是"决定不动的地方"这个意思吗?
粗略说是的。在有限元法中,每个节点拥有最多6个自由度(DOF:平移3个方向+旋转3个方向),将不想动的方向固定为零是最基本的约束。例如,如果钢板的一端用螺栓固定到墙上,那么该端的全部DOF被约束——即所谓的完全固定(Encastre)。
"不动"以外还有各种各样的类型吗?
是的,大致可以分为3种。(1) SPC(单点约束):将特定节点的特定DOF固定为固定值(通常为零)。最简单。(2) MPC(多点约束):将多个节点的变位用数式相互关联。RBE2和RBE3是这个家族。(3) 对称、反对称约束:利用模型的对称性来减少计算量的约束。它们的区别仅在于对整体刚性方程 $[K]\{u\} = \{F\}$ 的约束条件施加方式不同。
DOF约束的基础
约束DOF在数学上是什么样的操作?
从整体刚性方程中实质上消除约束DOF的行和列的操作。例如,当指定节点 $i$ 的X方向变位 $u_i = 0$ 时,将 $[K]$ 的对应行、列设为零并在对角分量上放置大值,或者直接删除那一行和列。在Nastran中是 SPC 卡,在Abaqus中是边界条件的"ENCASTRE"或"PINNED"。
ENCASTRE和PINNED有什么区别?
ENCASTRE是约束所有6个DOF的完全固定。PINNED仅约束平移3个DOF,旋转是自由的——也就是说是销支座。在梁结构分析中,一端设为ENCASTRE(固定端),另一端设为PINNED(销支座)来对比悬臂梁和简支梁,这样的模式在实务中也是常见的。
刚体运动与约束不足
先生说"约束不足会因为刚体运动而计算失败",刚体运动是什么?
刚体运动(Rigid Body Motion)是指物体不发生变形而进行平移或旋转的运动。三维空间中有平移3个方向(X、Y、Z)和旋转3个方向(Rx、Ry、Rz)共6个刚体模式。如果不用约束来抑制这些,$[K]$ 会变成奇异矩阵,逆矩阵无法求出,求解器会报错。
那只要约束6个DOF就可以了吗?
原则上是的。但是"6个"的选择方式很重要。例如,如果一个点的全6个DOF被约束,该点的反力会集中,导致局部出现不现实的应力。实务中经常使用3点约束(1-2-3法)。第1个点抑制平移3个方向,第2个点抑制旋转2个方向,第3个点抑制剩余1个方向的旋转——这样约束反力分散,局部应力问题得以缓和。
具体会出现什么样的错误信息?
Nastran会显示"FATAL: MATRIX IS SINGULAR"或"MECHANISM DETECTED",Abaqus会显示"NUMERICAL SINGULARITY"或"ZERO PIVOT",这些是典型的。如果出现这样的错误,首先要审查约束条件。另外Nastran有 PARAM,BAILOUT,-1 选项可以在错误时也输出结果,查看变位的动画就能一目了然看出哪个方向发生刚体运动。
对称边界条件
我听说可以用对称边界条件把模型分成一半来计算,这是什么原理?
假设YZ平面是对称面。对对称面上的节点,约束垂直于对称面的平移($u_x = 0$)和对称面内的旋转($\theta_y = 0, \theta_z = 0$)。这样的话,以对称面为界存在"镜像模型"的效果就等价了。模型尺寸减半,计算时间就减少到1/4~1/8左右(内存也减半)。
还有反对称这种情况吗?
有的。反对称的情况下,约束平行于对称面的平移,垂直方向自由——也就是说和对称相反的组合。例如,左右对称的梁受到反对称荷载(如扭转等)时使用。但是只有在荷载也包括对称、反对称的情况下才能使用。也就是说,形状虽然对称但荷载不对称就无法使用——忽视这一点会导致结果异常。
MPC(多点约束)
MPC是"Multi-Point Constraint"吧。这是什么约束?
MPC是在多个节点之间定义数式关系的约束。一般形式为:
$$\sum_{j} A_j \cdot u_j = 0$$
其中 $A_j$ 是系数,$u_j$ 是各节点的DOF。例如"节点1和节点2的X方向变位相等"就是 $u_1 - u_2 = 0$,即 $A_1 = 1, A_2 = -1$。在实务中,用于连接不同网格密度的边界、螺栓联结部分的建模、刚体领域的表示等。应用范围非常广泛。
设定MPC时需要注意什么?
最常见的错误是"过约束"。当相同的DOF同时被SPC和MPC重复定义时,就会产生矛盾的约束,求解器会报错,最坏的情况下会无声地输出错误的解。另外,MPC的从属节点(dependent node)和独立节点(independent node)的选择也很重要。如果把有荷载作用的节点设为从属节点,该荷载可能无法正确传递。
RBE2和RBE3的区别
RBE2和RBE3的区别还是有点模糊……
这是FEM初级~中级用户几乎都会遇到的坑,让我们清楚地整理一下。
RBE2(Rigid Body Element, Type 2)是刚体结合。一个独立节点(independent node)对应多个从属节点(dependent nodes),从属节点完全追随独立节点的运动。独立节点移动,从属节点也同样移动。也就是说增加刚性。经常用于螺栓孔周围的联结或质量单元的固定。
RBE3是荷载分配单元。一个从属节点的变位等于多个独立节点变位的加权平均。不增加刚性是最大的要点。用于将面荷载集约到一个点或求测量点的代表值等。
等等,RBE2增加刚性,RBE3不增加?具体会有什么影响?
例如考虑圆孔周围的节点群。用RBE2将其连到中心的独立节点时,孔周围会完全刚体化,不会发生变形——这比实际的螺栓联结硬得多。而用RBE3将荷载输入中心时,荷载只是分配到周围节点,孔周围可以自由变形。
现场常见的失败是本应使用RBE3的地方却用了RBE2,导致局部引入过度刚性,应力值与实际完全不同。"想要传递力"还是"想要分配荷载"——根据这一点来选择是铁则。
实务要点
总结一下,设定约束条件时最重要的是什么?
有3个要点。
第一,物理上的妥当性。忠实再现实际结构物是如何支持的。螺栓固定到墙上就用ENCASTRE,轴承支座就只约束一个方向的平移——这样做。过度固定会导致应力出现偏低、评价在危险一侧,约束不足就根本得不到解。
第二,约束反力的确认。分析运行后一定要检查约束点的反力。全部反力的合计是否与外力平衡,特定约束点是否出现不现实的反力——通过这个可以检出大部分约束条件的错误。
第三,约束的影响范围。约束点附近容易产生应力集中,所以理想情况下应在要评价的领域充分远离的位置放置约束。根据圣维南原理,离约束代表尺寸的2~3倍距离处影响会衰减。
相关术语
- SPC(Single Point Constraint):将单个节点的特定DOF固定为固定值的单点约束
- MPC(Multi-Point Constraint):在多个节点间定义数式关系的多点约束
- RBE2:刚体结合单元。从属节点完全追随独立节点。增加刚性
- RBE3:荷载分配单元。从属节点的变位是独立节点群的加权平均。不增加刚性
- DOF(Degree of Freedom):自由度。节点拥有的独立运动分量
- 刚体运动(Rigid Body Motion):不伴随变形的平移、旋转运动。约束不足时产生
- 对称边界条件:约束垂直于对称面的平移和面内旋转来实现半模型的边界条件
- 圣维南原理:荷载或约束的局部详细特征在充分远处不产生影响的原理
正确理解CAE术语是团队内部交流的基础。— Project NovaSolver也在考虑实务工作者的学习支持。
请告诉我们在约束条件实务中感到的课题
Project NovaSolver致力于解决CAE工程师日常面临的课题——建模的复杂性、计算成本、结果的解释等。您的实务经验将成为开发更好工具的原动力。
咨询(筹备中)相关主题
更详细
错误