4節点四辺形要素(QUAD4)
理论与物理
Q4单元 — 二维FEM的主力
老师，4节点四边形单元（Q4）是二维FEM的基础吧。
形函数
Q4的形函数在自然坐标 $(\xi, \eta)$ 中是双线性的（bilinear）：
$$ N_i = \frac{1}{4}(1 + \xi_i \xi)(1 + \eta_i \eta) $$
Q4的形函数在自然坐标 $(\xi, \eta)$ 中是双线性的（bilinear）：
4个节点对应 $(\xi, \eta) = (\pm 1, \pm 1)$。
双线性意味着有 $1, \xi, \eta, \xi\eta$ 这4项。存在交叉项 $\xi\eta$。
正是由于这个交叉项，Q4比3节点三角形单元（CST）具有更高的弯曲表现能力。但要表现完整的弯曲变形，需要 $\xi^2, \eta^2$ 项，而Q4没有这些项。
剪切自锁问题
请给我讲讲Q4的剪切自锁。
使用完全积分（2×2 Gauss点）的Q4试图表现纯弯曲时，会产生寄生的剪切应变，导致单元“变硬”。这就是剪切自锁。
物理上会这样：
- 纯弯曲时上面拉伸，下面压缩
- Q4的双线性形函数在表现这种变形时，不可避免地会产生剪切应变
- 这种剪切应变会储存额外的能量，导致单元变得过于刚硬
对策呢？
有三种方法：
1. 减缩积分（1×1 Gauss点） — 用一个积分点评估剪切。自锁消失但会产生沙漏模式
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个自由度） | — |
| 曲面近似 | — | 仅直边（曲面用折线近似） |
总结
我来整理一下Q4的理论。
要点：
- 双线性形函数 — 包含交叉项 $\xi\eta$。精度高于CST
- 剪切自锁 — 完全积分时发生。对策：减缩积分或非协调模式
- Nastran的CQUAD4内置非协调模式 — FEM史上最成功的单元之一
- 二维实务标准 — 广泛用于平面应力、平面应变、轴对称问题
- 对扭曲形状敏感 — 管理长宽比和扭曲度很重要
在HEX8页面学到的自锁和沙漏的概念，在Q4上也完全适用呢。
Q4是HEX8的二维版，所以会出现完全相同的问题。在二维上理解了，扩展到三维就很容易。
Q4单元的等参公式化
Q4单元的等参公式化是1960年代由Irons和Ergatoudis等人在伯明翰大学开发的。通过使用自然坐标（ξ, η），可以将任意四边形形状映射到单位正方形，结合高斯积分实现了高精度的刚度矩阵计算。至今仍被用作商用代码的核心单元。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。你有过急刹车时身体被向前甩出去的经历吗？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动，而建筑物的质量“被落下”。静力分析中此项设为零，这是“缓慢加载所以加速度可忽略”的假设。对于冲击载荷或振动问题，此项绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时会感觉到“想恢复原状的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸得更长？当然是橡皮筋。这种“不易伸长性”就是杨氏模量 $E$，它决定了刚度。常见的误解：“刚度高=强度高”是不对的。刚度是“不易变形性”，强度是“不易破坏性”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（如重力）和表面力 $f_s$（如压力、接触力）。可以这样想——桥上卡车的重量是“作用在整个内容物上的力”（体积力），轮胎压路面的力是“只作用在表面的力”（表面力）。风压、水压、螺栓紧固力…全都是外力。这里容易犯的错误：弄错载荷方向。本想“拉伸”却变成了“压缩”——听起来像笑话，但在3D空间中坐标系旋转时确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。试着弹一下吉他的弦。声音会一直持续吗？不，会逐渐变小。因为振动能量通过空气阻力和弦的内部摩擦变成了热能。汽车的减震器也是同样原理——特意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃下去。实际上不会这样，所以设置适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入mm时，载荷·弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢：约210 GPa，铝：约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中为tonne/mm³（钢约为 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系用N，m系也用N统一 |
数值解法与实现
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 | 过大则精度下降 |
相关主题结构分析8節点四辺形要素(QUAD8)结构分析8節点六面体要素(HEX8)伝熱解析2次元定常熱伝導 — 数値解法用語集完全積分 — CAE用語解説伝熱解析3次元定常熱伝導 — 数値解法结构分析Plate Buckling
この記事の評価 ご回答ありがとうございます！ 参考に なった もっと 詳しく 誤りを 報告 |