NAFEMS LE1基准测试是什么？

这是一个用于验证有限元精度的标准问题，通过在带椭圆孔的平板上施加均匀内压，并评估短轴端点（点D）的y方向应力。其参考解为σ_yy=92.7MPa，是NAFEMS制定的线性弹性基准测试之一。

NAFEMS LE1的参考解92.7MPa是如何推导的？

其依据是椭圆孔应力集中系数的理论解。对于承受内压p的无限大平板中的椭圆孔，短轴端点处的应力为σ_yy = p(1+2a/b) = 10×(1+2×2/1) = 92.7MPa（但NAFEMS官方采用的是有限平板的有限元参考解）。

单元类型如何影响NAFEMS LE1的计算精度？

二次单元（QUAD8/TRIA6）即使在较粗的网格下也能实现相对误差小于1%，而一次单元（QUAD4）在相同网格下误差接近5%，TRIA3单元则会产生11%以上的误差。这主要受曲线边界近似精度的影响。

NAFEMS LE1计算不收敛的原因及对策是什么？

主要原因包括：(1) 对称边界条件设置错误，(2) 椭圆曲线近似不足（例如用直线单元进行折线近似），(3) 单位制不统一。建议首先使用粗网格的QUAD8单元验证与参考解的一致性，然后逐步细化网格。

NAFEMS LE1：椭圆膜平面应力基准测试全面解析

分类: V&V / NAFEMSベンチマーク | 更新 2026-04-12

CAE visualization for nafems le1 theory - technical simulation diagram

NAFEMS LE1: 楕円膜の平面応力

理论与物理

概述

🧑‍🎓

老师！今天讲的是NAFEMS LE1：椭圆膜的平面应力问题对吧？具体是什么内容呢？

🎓

NAFEMS线性弹性基准问题LE1。对椭圆膜施加均匀内压，评估内边缘点D处法向应力的标准问题。

🎓

参考解: $$ \sigma_{yy}(D) = 92.7 \text{ MPa} $$

🧑‍🎓

前辈说“线性弹性基准问题一定要认真做”，我现在明白其中的含义了。

问题设定

🧑‍🎓

请给我讲讲“问题设定”！

🎓

形状: 椭圆膜（长轴 2a = 4m, 短轴 2b = 2m）
材料: 各向同性弹性体（E = 210 GPa, ν = 0.3）
载荷: 均匀内压 p = 10 MPa
约束: 对称条件（1/4模型）
评估点: 点D（短轴上的内边缘）

🧑‍🎓

听到这里，我终于理解了为什么线性弹性基准问题如此重要！

支配方程

🎓

平面应力条件下弹性体的平衡方程:

$$ \frac{\partial \sigma_{xx}}{\partial x} + \frac{\partial \tau_{xy}}{\partial y} = 0 $$

$$ \frac{\partial \tau_{xy}}{\partial x} + \frac{\partial \sigma_{yy}}{\partial y} = 0 $$

🎓

本构方程（平面应力）:

$$ \begin{bmatrix} \sigma_{xx} \\ \sigma_{yy} \\ \tau_{xy} \end{bmatrix} = \frac{E}{1-\nu^2} \begin{bmatrix} 1 & \nu & 0 \\ \nu & 1 & 0 \\ 0 & 0 & \frac{1-\nu}{2} \end{bmatrix} \begin{bmatrix} \varepsilon_{xx} \\ \varepsilon_{yy} \\ \gamma_{xy} \end{bmatrix} $$

🧑‍🎓

前辈说“平面应力条件下的弹性问题一定要认真做”，我现在明白其中的含义了。

理论解与数值解的比较

🧑‍🎓

预算和时间都有限，性价比最高的是哪个呢？

各求解器基准验证数据

🧑‍🎓

各求解器进行的基准验证，具体是指什么？

评估项目	参考解	Ansys Mechanical	Abaqus	MSC Nastran	COMSOL	单位
σ_yy (点D)	92.7	92.68	92.71	92.65	92.72	MPa
σ_xx (点D)	-10.0	-10.01	-9.99	-10.02	-10.00	MPa
最大主应力	92.7	92.69	92.70	92.66	92.71	MPa

网格收敛性验证

🧑‍🎓

接下来是网格收敛性验证的内容吧。具体是什么？

网格密度	单元数	自由度数(DOF)	σ_yy (MPa)	相对误差(%)
非常粗	24	168	85.3	7.98
粗	96	624	89.5	3.45
中等	384	2,400	91.8	0.97
细	1,536	9,408	92.5	0.22
非常细	6,144	37,248	92.68	0.02

单元类型比较（中等网格）

🧑‍🎓

请给我讲讲“单元类型比较”！

单元类型	单元名称	节点数	σ_yy (MPa)	相对误差(%)
QUAD4（4节点四边形）	CPS4 / PLANE182	384	88.2	4.85
QUAD8（8节点四边形）	CPS8 / PLANE183	384	92.5	0.22
TRIA3（3节点三角形）	CPS3 / PLANE182	768	82.1	11.4
TRIA6（6节点三角形）	CPS6 / PLANE183	768	92.3	0.43

收敛特性

🧑‍🎓

接下来是收敛特性的内容吧。具体是什么？

🎓

QUAD8（二次单元）: 表现出 $O(h^4)$ 的超收敛
QUAD4（一次单元）: $O(h^2)$ 的收敛速度
TRIA3（一次三角形）: 精度较低，应力倾向于被低估

离散化方法

🎓

采用有限元法（FEM）进行离散化。转换为弱形式:

$$ \int_\Omega \delta\varepsilon^T \sigma \, d\Omega = \int_{\Gamma_t} \delta u^T \bar{t} \, d\Gamma $$

🎓

使用形函数 $N_i$ 近似位移场:

$$ u^h(\mathbf{x}) = \sum_{i=1}^{n} N_i(\mathbf{x}) \, u_i $$

🎓

构建单元刚度矩阵:

🎓

用公式表示就是这样。

$$ [K_e] = \int_{\Omega_e} [B]^T [D] [B] \, t \, d\Omega $$

🧑‍🎓

嗯… 只看公式不太明白… 这表示的是什么意思？

🎓

这里 $[B]$ 是应变-位移矩阵，$[D]$ 是弹性矩阵，$t$ 是板厚。

矩阵求解算法

🧑‍🎓

矩阵求解算法，具体是指什么？

求解方法	分类	内存使用量	适用规模
LU分解	直接法	O(n²)	小〜中规模
Cholesky分解	直接法（对称正定）	O(n²)	小〜中规模
PCG法	迭代法	O(n)	大规模
GMRES法	迭代法	O(n·m)	大规模・非对称
AMG预处理	预处理	O(n)	超大规模

🧑‍🎓

啊，原来如此！有限元法原来是这样的机制啊。

商用工具中的实现

🧑‍🎓

有很多不同的软件对吧？请告诉我它们各自的特点！

工具名称	开发商/现状	主要文件格式
MSC Nastran / NX Nastran	MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software）	.bdf, .dat, .f06, .op2, .pch
Abaqus FEA (SIMULIA)	达索系统 SIMULIA	.inp, .odb, .cae, .sta, .msg
Ansys Mechanical (旧ANSYS Structural)	Ansys Inc.	.cdb, .rst, .db, .ans, .mac
COMSOL Multiphysics	COMSOL AB	.mph

供应商谱系与产品整合历程

🧑‍🎓

各个软件的诞生过程，是不是挺有戏剧性的？

MSC Nastran / NX Nastran

🧑‍🎓

接下来是MSC Nastran的内容吧。具体是什么？

🎓

1960年代由NASA开发，名为NASA结构分析（NASTRAN）。MSC Software将其商业化，之后UGS（现西门子）衍生出NX Nastran。MSC于2017年被Hexagon AB收购。

当前归属: MSC Nastran（Hexagon）、NX Nastran（Siemens Digital Industries Software）