膜褶皱（起皱）分析
膜褶皱（起皱）的理论基础
膜褶皱是什么
老师，膜结构的褶皱（wrinkling）是什么现象？
膜结构弯曲刚度几乎为零。产生压缩应力时就会形成褶皱。在太空的太阳能电池板、安全气囊、帐篷结构中会成为问题。
褶皱的力学
膜上产生压缩应力时：
1. 膜无法承受压缩 — 由于没有弯曲刚度，屈曲=褶皱
2. 褶皱的方向 — 垂直于压缩方向形成褶皱
3. 褶皱的波长 — 取决于膜的张力、板厚、曲率
FEM中的建模
两种方法：
1. 壳单元（薄板厚） — 直接模拟褶皱形状。NLGEOM=YES + 初始缺陷导致屈曲→褶皱
2. 膜单元（弯曲刚度为零）+ 褶皱模型 — 将压缩应力设为零的“张力场理论”
不能用膜单元直接表现褶皱吗？
膜单元因为没有弯曲刚度，受压时不会出现褶皱的“形状”。取而代之的是将压缩应力设为零，求得“产生褶皱状态下的应力场”。例如Abaqus的*NO COMPRESSION和膜褶皱算法。
总结
储罐液体晃动与膜褶皱的起源
膜褶皱理论源于拉伸与压缩的差异。区分不会产生褶皱的“张力膜”和存在压缩应力的“褶皱膜”的是Stein和Hedgepeth（1961年，NASA）。他们确立了在产生褶皱的区域将主压缩方向的应力设为0的“松弛主应力理论”。当前FEM褶皱分析的理论基础就源于这篇1961年的论文。
数值解法与实现
壳单元方法（求褶皱形状）
```
*SHELL SECTION
0.025, 5 $ 板厚0.025mm（膜）
*STEP, NLGEOM=YES
*STATIC, RIKS $ 褶皱是屈曲的一种
```
给定初始缺陷（一阶屈曲模态形状）以诱发褶皱图案。
膜单元方法（张力场）
```
*MEMBRANE SECTION
0.025
*NO COMPRESSION $ 将压缩应力设为零
```
不会出现褶皱形状，但能得到褶皱区域的应力场。
总结
褶皱有限元法：修正材料特性法
用FEM处理膜褶皱的方法有①将主应力设为零截断的“松弛刚度法”②仅计算无褶皱区域的“分析追踪法”③用精细网格再现实际几何褶皱的“屈曲分析法”。实用上修正材料特性法（松弛刚度法）最稳健，Abaqus的膜单元M3D4R与褶皱判定子程序的组合在产业界被广泛使用。
膜褶皱（起皱）膜褶皱（起皱）实践指南
褶皱分析的实务
在太空结构（太阳能电池板、膜结构天线）中最为重要。褶皱会劣化光学面的精度。
实务检查清单
太空太阳能发电卫星的展开膜分析
太空太阳能发电卫星（SSPS）的薄膜太阳能电池板（厚度0.01mm）在展开后可能产生热变形和褶皱。JAXA从2010年代开始进行聚酰亚胺薄膜（厚度12.5μm）的热-结构耦合褶皱分析，确定了太阳光入射压力与热膨胀组合产生的3~5mm波长褶皱对发电效率的影响，并反映在膜张力设计中。
膜褶皱（起皱）软件与求解器比较
褶皱分析的工具
Ansys Mechanical薄膜分析的特殊设置
Ansys Mechanical的Shell181或Membrane单元（SHELL181, KEYOPT(1)=1）可作为面外刚度为零的纯膜单元使用。褶皱判定与“仅允许 principal stress ≥ 0”的材料输入（USERFLD + 松弛刚度子程序）组合实现是标准方法。ESA在设计欧洲航天器的膜型大型天线（直径15m）时使用了此方法，以±3%的精度预测了轨道上的褶皱形状。
尖端技术
褶皱的尖端研究
褶皱尖端的奇异场与有效膜厚
在褶皱扩展的膜中，实际有效厚度变为“膜厚÷褶皱数”，等效弯曲刚度大幅降低。用分子动力学模拟计算厚度0.1mm铝箔的褶皱尖端时发现，局部应力在2020年代被证实可达体材料屈服应力的3~5倍。在折叠太空展开结构的寿命设计中，这种局部应力的大小是决定裂纹萌生寿命的重要参数。