4节点四边形单元(QUAD4)
4节点四边形单元(QUAD4)的理论基础
Q4单元 — 二维FEM的主力
老师，4节点四边形单元（Q4）是二维FEM的基础吧。
形状函数
Q4的形状函数在自然坐标 $(\xi, \eta)$ 中是双线性（bilinear）：
4个节点对应 $(\xi, \eta) = (\pm 1, \pm 1)$。
双线性意味着有 $1, \xi, \eta, \xi\eta$ 四项。包含交叉项 $\xi\eta$。
正是这个交叉项使Q4的弯曲表现能力比3节点三角形（CST）更强。但要完全表现弯曲变形需要 $\xi^2, \eta^2$ 项，Q4中没有这些。
剪切锁定的问题
请说明Q4的剪切锁定。
完全积分（2×2高斯点）的Q4在表现纯弯曲时，会产生寄生剪切应变，使单元"僵化"。这就是剪切锁定。
物理上是这样的：
- 纯弯曲时上面伸长，下面缩短
- Q4的双线性形状函数在表现这种变形时，不可避免地产生剪切应变
- 这个剪切应变的能量被多余地储存，使单元变硬
解决办法呢？
有三种方法：
1. 低减积分（1×1高斯点） — 只用一个积分点评估剪切。锁定消失但会出现沙漏模式
2. 非协调模式（Incompatible modes） — 添加对应 $\xi^2, \eta^2$ 的内部自由度。Nastran的CQUAD4、Abaqus的CPS4I
3. Assumed Natural Strain (ANS) — 单独假设剪切应变
Nastran的CQUAD4默认内置非协调模式吧。
是的。Nastran的CQUAD4被称为FEM历史上最成功的单元之一。MacNeal和Harder在1985年发表的改进版本，通过非协调模式大幅改善了弯曲精度，同时满足补丁测试。
Q4的优点和缺点
| 特性 | 优点 | 缺点 |
|---|---|---|
| 网格生成 | 映射网格容易生成 | 自由网格中混有三角形 |
| 精度 | 非协调模式下精度高 | 歪斜形状精度下降 |
| 计算成本 | 低（8个DOF） | — |
| 曲面近似 | — | 仅直边（曲面用折线近似） |
总结
让我整理一下Q4的理论。
要点：
- 双线性形状函数 — 包含交叉项 $\xi\eta$。精度高于CST
- 剪切锁定 — 在完全积分中产生。通过低减积分或非协调模式解决
- Nastran的CQUAD4内置非协调模式 — FEM史上最成功的单元之一
- 二维实务标准 — 广泛用于平面应力、平面应变、轴对称
- 对歪斜形状敏感 — 纵横比、歪斜度的管理很重要
在HEX8的页面中学到的锁定和沙漏概念，在Q4中也完全适用吧。
Q4是HEX8的二维版本，所以会发生完全相同的问题。在二维中理解它，扩展到三维就容易了。
Q4单元的等参数公式
Q4单元的等参数公式由Irons和Ergatoudis等人在1960年代于伯明翰大学开发。利用自然坐标（ξ, η）可以将任意四边形形状映射到单位正方形，与高斯积分组合实现高精度刚度矩阵计算。该方法至今仍被商用软件采用为核心单元。
4节点四边形单元(QUAD4)的数值计算方法
Q4的积分方案比较
请说明Q4各积分方案的特点。
沙漏有3个模式是说8个自由度中有3个是零能？
Q4的8个自由度减去刚体模式3个（平移2个+旋转1个）剩5个变形模式，低减积分只能评估其中2个。剩余3个就是沙漏模式。
求解器中的单元名称
| 变种 | Nastran | Abaqus | Ansys |
|---|---|---|---|
| 标准Q4 | CQUAD4 | CPS4（平面应力） | PLANE182 |
| 低减积分 | — | CPS4R | PLANE182(red.) |
| 非协调模式 | CQUAD4（默认） | CPS4I | PLANE182(EAS) |
| 混合型 | — | CPE4H（平面应变） | — |
Nastran的CQUAD4默认内置非协调模式，所以不需要特别设置吧。
Nastran用户生活在"使用CQUAD4就行"的简洁世界里。Abaqus和Ansys用户需要自己选择积分方案，这需要更多知识。
网格质量的影响
Q4对形状歪斜敏感吗？
非常敏感。理想的正方形Q4和歪斜平行四边形Q4的精度差异很大。
| 品质指标 | 理想值 | 允许范围 | 对精度的影响 |
|---|---|---|---|
| 纵横比 | 1.0 | < 5.0 | 大时精度下降 |
| 歪斜度 | 0° | < 45° | 大时精度下降 |
| 锥形比 | 0 | < 0.5 | 大时精度下降 |
| 内角 | 90° | 45°〜135° | 范围外精度下降 |
内角在45°以下的压扁四边形不行吧。
是的。压扁四边形使雅可比矩阵在不同位置差异很大，积分精度下降。即使是非协调模式单元，歪斜形状的精度也不能保证。Q4网格质量管理是二维FEM中最重要的工作。
补丁测试
补丁测试是什么？
对歪斜单元的集合（patch）施加均匀应变，检查是否能精确产生均匀应力的测试。这是FEM单元最低品质保证。
非协调模式单元有时不严格满足补丁测试（MacNeal-Harder的补丁测试通过，但一般歪斜patch有微小误差）。实际应用中不是问题，但理论上有人在意。
总结
Q4的数值方法，整理一下。
要点：
- Nastran的CQUAD4默认非协调模式 — 无需设置，高精度
- Abaqus中选择CPS4I（非协调模式） — CPS4R（低减积分）也可，但要注意沙漏
- 网格质量支配Q4的精度 — 严格管理纵横比、歪斜度
- 补丁测试 — FEM单元最低品质保证
- Q4+非协调模式是二维FEM最推荐 — 有疑问时用这个
Q4单元的2×2高斯积分
Q4单元通常使用2×2（4点）高斯积分。积分点坐标为±1/√3≈±0.5774，可以完全积分多项式次数3以下的函数。采用1点积分（低减积分）时计算速度约快4倍，但会出现称为沙漏模式的零能量模式，需要稳定化技术。
4节点四边形单元(QUAD4)的实务应用
Q4的实务应用
Q4在二维分析中的哪些地方使用？
Q4 vs. Q8 的使用选择
应该选择Q4还是8节点四边形（Q8）？
| 判断标准 | Q4 | Q8 |
|---|---|---|
| 应力集中评估 | 需要细网格 | 粗网格也精确 |
| 曲面近似 | 直边（折线近似） | 曲边（二次近似） |
| DOF | 8（4节点×2） | 16（8节点×2） |
| 计算成本 | 低 | Q4的2〜3倍 |
| 推荐情景 | 网格充分细致时 | 精度优先、网格粗时 |
网格充分细致的话Q4和Q8结果相同？
收敛的结果相同。但Q4收敛慢（应力收敛为 $O(h)$）而Q8快（$O(h^2)$），所以Q4需要Q8的2〜5倍网格。总DOF数上Q8通常更高效。
那就用Q8好了？
Q8虽然优越，但Q4的优势是与自动网格生成的兼容性好。Q4的自动网格很稳定，但Q8的自动网格在中间节点放置时有问题。而且显式法的标准是Q4（稳定时间增量大）。
网格生成的最佳实践
Q4网格生成的要点是？
O-grid是什么？
在孔周围采用放射状和同心圆状的Q4网格排列方式。这是精确捕捉孔应力集中的经典网格技巧。用HyperMesh或Abaqus/CAE的分割功能可以生成。
实务检查清单
请给Q4的检查清单。
为什么三角形退化这么差？
将Q4退化为三角形（两个节点在同一位置）时，雅可比矩阵在退化点变成零，积分精度严重下降。需要三角形时应使用专用的3节点三角形单元（CST/LST）。
Q4单元的纵横比管理
NASTRAN的QUAD4单元推荐纵横比在5以下。实务中对薄板结构进行面内应力求解时，将单元边长比控制在3以内既可以抑制Warning 3035又能改善精度。Altair HyperMesh的Quality Index在"Aspect Ratio"检查中以5.0为阈值设置已成为汽车部品分析的事实标准。
4节点四边形单元(QUAD4)的软件比较
Q4在各求解器中的特点
各求解器的Q4单元有区别吗？
"钻孔DOF"是什么？
Q4仅有面内的2个自由度（$u, v$），但在框架分析的Q4壳中有时需要面外的旋转自由度（钻孔旋转）。Nastran的CQUAD4具有钻孔DOF，因此与框架单元（梁单元）的连接自然。
选择指南
Q4的选择指南是？
Q4是"FEM基础的基础"。任何求解器都必须用到的单元吧。
正是这样。能否正确使用Q4是FEM工程师的基本技能。
各求解器Q4单元实现比较
Nastran的CQUAD4、Abaqus的S4R、Ansys的SHELL181都是4节点壳单元，但内部公式不同。S4R使用低减积分+沙漏控制，SHELL181基于Reissner-Mindlin，内置剪切锁定对策。LS-DYNA的Type 2壳默认采用1点积分，采用了针对高速碰撞分析的Belytschko-Lin-Tsay公式。
4节点四边形单元(QUAD4)的先端研究
Q4的先端研究
Q4还有改进空间吗？
已研究50多年的单元，仍在持续改进。
光滑有限元法 (S-FEM) 的二维版本
将Edge-based S-FEM (ES-FEM)应用到Q4，精度比普通Q4提高，接近Q8的结果。通过对边进行应变平滑化，隐式表现交叉项以上的变形模式。
多边形有限元法
不仅限于四边形，还可用五边形、六边形、N边形的二维单元。称为多边形单元法。利用Voronoi分割自然生成的网格，大大放松网格生成的约束。
Q4和Q8都超越了的通用多边形单元…。
用Talisca或Wachspress基函数为任意多边形构造形状函数。还在研究阶段，但未来有可能取代Q4/Q8。
基于物理的形状函数
通常Q4采用多项式形状函数，但也有使用支配方程的解作为形状函数的方法。例如以拉普拉斯方程的基本解作为形状函数的"Trefftz单元"。对特定问题（应力集中、裂纹尖端）精度比普通Q4高好几个数量级。
总结
Q4的先端研究，总结一下。
Q4有50多年的历史，但仍有改进空间。这表明"二维四边形单元"对结构力学有多么重要。
Q4单元的非协调模式扩展
1973年Wilson和Taylor等加州大学伯克利分校的研究人员为Q4单元添加了非协调位移模式（Wilson-Taylor单元），大幅改善了弯曲精度。增加的自由度在内部通过静力缩减，因此总体刚度矩阵尺寸不变。Abaqus的CPS4I和Nastran的CQUAD4（PSHELL）改进版继承了这一技术。
4节点四边形单元(QUAD4)的故障排除
Q4的故障
请说明Q4单元的常见故障。
二维FEM的基本问题都集中在Q4中。
位移过小（剪切锁定）
弯曲问题中位移比理论值小。
完全积分Q4中出现了剪切锁定。解决办法与三维HEX8相同：
- 改用CPS4I（非协调模式）
- 改用CPS4R（低减积分）
- 切换到Q8（CPS8R）
沙漏模式（位移呈锯齿状）
低减积分Q4位移呈波浪。
集中荷载激发了沙漏模式。解决办法是分散荷载或改用CPS4I（非协调模式）。
应力云图出现棋盘模式
应力呈棋盘花纹。
有两个可能：
1. 低减积分的沙漏 — 改用CPS4I
2. 网格歪斜 — 确认单元质量。是否存在纵横比 > 5 或歪斜度 > 45° 的情况
厚度设置遗漏
结果的力量级不对。
二维平面应力单元需要设置板厚。如果采用默认的1.0，实际板厚不同时力和挠度的量级就不对。在Abaqus的 *SOLID SECTION 中必须指定厚度。
边界条件过约束
我本想固定边界，却把所有节点的全部自由度都拘束了，结果应力不对。
由于泊松效应，横向膨胀变形被拘束了。必须明确"固定"的含义：
- 仅拘束面外位移 — 最常见
- 面内也拘束 — 泊松膨胀也被抑制。确认是否符合实际结构
总结
Q4的故障处理，整理一下。
Q4的故障与HEX8有相同的结构。在二维中理解，就可应用到三维。
Q4单元沙漏模式诊断
在1点积分Q4单元中产生沙漏（hourglass）模式时，相邻单元的位移会呈交替的正负振荡，在云图中显示为棋盘花纹。MSC Nastran可用DIA=16标记检测零能量模式。Abaqus默认采用人工沙漏控制刚度，HOURGLASS=ENHANCED选项据称精度最高。
相关主题
更多
报告