热-机械有限元分析方法研究及其在小型复合材料低温液氢贮箱案例中的演示

《Applied Mechanics》：An Investigation of Thermo-Mechanical Finite Element Analysis Methodologies and a Demonstration on the Case of a Small-Scale Composite Cryogenic Hydrogen Tank

　　

为了减少并最终消除CO₂排放，航空工业正在研究替代燃料，其中液氢（LH₂）是最有前景的选择之一。多个涉及液氢飞机子系统设计的研究项目已经启动，重点关注轻量化低温贮箱的开发。复合材料是低温航空应用的潜在候选材料，但其存在热应力问题，会导致各种类型的损伤，需对此进行深入研究，以实现可靠的复合材料低温结构。在此方向上，本研究针对特征几何结构构型，研究了替代性的有限元热-机械分析技术。通过将所研究分析方法得出的热-机械应力分析结果与已发表结果进行比较，验证了其正确性。对替代性建模技术的评估有助于将有限元建模开发策略调整为所需分析类型和预期结果，并可直接应用于低温部件的热-机械设计。研究成果以LH₂储存系统外贮箱的热-机械分析案例进行了演示，该分析采用壳单元与实体单元的最优组合完成。

**论文解读：复合材料低温液氢贮箱热-机械有限元分析方法研究**

**研究背景与现存问题**

为减少并最终消除CO₂排放，航空工业正积极开发替代燃料，液氢（LH₂）被公认为最具前景的选项之一。液氢飞机的子系统设计需要开发轻量化低温贮箱，复合材料凭借高比强度和高比刚度成为理想候选材料。然而，复合材料中纤维与基体的热膨胀系数（CTE）存在差异，温度变化时会产生两类热应力：层内由纤维-基体CTE失配引起的应力，以及层合板层面因各铺层CTE各向异性导致的层内与层间热应力。历史上X-33试验航天飞机贮箱爆炸即为热应力导致灾难性失效的典型案例。当前微观力学尺度上的热应力研究已较充分，但在层合板及全局结构尺度上，对有限元（FE）建模方法适用性的系统评估仍然匮乏。商业软件通常提供基于Reissner–Mindlin一阶剪切变形理论（FSDT）的壳单元与基于三维弹性理论的实体单元，二者精度与计算成本差异显著，如何按分析目标选择最优建模策略成为亟待解决的问题。

**研究内容与意义**

研究人员在欧盟H2ELIOS项目框架下，以ANSYS APDL为平台，系统比较了shell281分层壳单元与solid186实体单元两种建模技术在机械载荷、热载荷及热-机械耦合工况下的表现，覆盖平板、圆柱和球壳三类特征几何。通过四个文献基准算例验证数值模型准确性后，将最优策略应用于小型复合材料液氢贮箱的混合建模演示。研究填补了中尺度至全尺度结构FE建模方法评估的空白，为低温复合材料结构高保真热-机械分析奠定了基础。论文发表于《Applied Mechanics》。

**关键技术与方法**

研究采用shell281壳单元（基于Reissner–Mindlin FSDT，8节点，低内存需求）与solid186实体单元（基于三维弹性理论，20节点，高内存需求），对四个文献基准算例进行复现验证：Fotopoulos和Lampeas的机械载荷正交层合板、Sayyad等的热载荷层合板、Ehteram等的热载荷纤维缠绕管、Brischetto和Carrera的热载荷球壳夹芯板。材料参数取自各参考研究，应力结果按文献方法归一化以消除尺寸与载荷强度影响。演示模型采用IM7/977-2碳纤维复合材料及AISI 316不锈钢法兰，低温材料性能取自美国国家标准与技术研究院（NIST）数据库。

**研究结果**

**4. 机械载荷下的平板层合板**：实体单元层间剪应力与参考解最大偏差为0.5%–1.1%，壳单元则为2.9%–22.8%；两者面内正应力精度均较高（偏差不超过2%）。壳单元仅用625个单元完成全板建模，实体单元四分之一板即需25,600个单元。

**5. 热载荷下的平板层合板**：实体单元层间剪应力偏差为1.4%–6.4%，壳单元高达17.9%–160%，证实壳单元无法准确计算由内部热载荷引起的层间应力，其根源在于程序自动计算的剪切修正因子仅考虑外载。

**6. 热载荷下的复合材料圆柱管**：实体单元厚度方向剪应力rz偏差为21.6%，壳单元达169%；实体单元可计算厚度方向法向应力（偏差2.2%），壳单元则完全无法输出该分量。

**7. 热载荷下的球壳夹芯结构**：曲率增大使两种方法精度均下降；排除边界应力集中后实体单元层间剪应力偏差为4.9%，壳单元高达84.7%；壳单元厚度方向的位移结果仅能提供中面平均值。

**8. 小型LH₂贮箱演示**：在筒段与球冠等简单区域，壳单元与实体单元结果吻合良好（面内应力偏差2%–11%）；但在法兰连接区域，壳单元因无法考虑厚度方向应变及热收缩/膨胀，偏差高达85%–304%，验证了混合建模策略中关键区域必须采用实体单元的判断。

**讨论与结论**

讨论指出，壳单元计算效率极高，但仅适用于面内结果与薄壁简单结构；实体单元精度高但内存需求巨大——仅占模型16.2%的几何体却消耗了60.9%的单元数，故只宜用于关键局部区域。混合建模中壳-实体连接处需谨慎处理节点耦合，但仍可能出现间隙并导致应力集中。研究结论为：两种建模技术均经受住了机械与热载荷及三类几何的验证；壳单元适合计算面内热应力与位移，实体单元是热致厚度方向应力计算的唯一可靠选择，但计算成本高昂；曲率对两种方法精度均有负面影响；法兰等关键部件应建立局部实体模型。未来工作包括引入摩擦接触界面建模、垫片分离与泄漏预测，以及将复合材料视为纤维与基体两相分别建模的微观力学方法。