计算细长悬臂梁变形，慎用绑定

概述

机翼总长约为3000mm.机翼外蒙皮（T300复合材料）厚度从机身向翼展方向分别为* * * mm、* * mm、* mm，翼根螺栓连接局部区域厚度为* * * * mm，中间夹芯结构采用PVC泡沫。机翼通过铆钉与机翼固定座进行连接。

在机翼上方施加值为大小为 * N 的分布力载荷，计算变形量。

结果

注：实体单元泡沫为contact面，壳单元蒙皮为target面。

网格1，节点匹配度高

绑定，MPC算法(Multipoint constraint)
变形量251mm，与共节点的结果最为接近，误差5%内。

绑定，默认算法(Augmented Lagrangian)
变形量224mm

• 网格3，节点匹配度一般

绑定，MPC算法
变形量155mm

解决方法：

对于这种长薄梁结构，或推广至大面积薄壁夹芯结构，当模型包含1.壳单元与实体单元有接触关系，2.求解弯曲变形时，

1. 尽量将两者使用共节点进行描述。
2. 若无法共节点，使用MPC算法代替默认的Augmented Lagrangian或Pure Penalty算法。
3. 尽量提高两个体间的网格节点匹配度。（单通过加密网格可能无法解决这个问题）
4. 其它更好的方式：暂无。可能可以通过修改PINB改变结果，但难以知道合理的值。

来源：shengyutou