张其林.索和膜结构.同济大学出版社
  
P137第十八章“空间曲面展开平面的等效杆单元有限元法”第137页[kj]矩阵中的L3从那里来？应该属于错误吧
  
本人认为该方法没有能够保证节点运动后三角形单元不重合！基于三角形边长变化最小收敛的
  
最后结果可能是一团相互重叠的三角形！

你说的这个问题，我也早就考虑过，
的确有这种可能性的，特别是曲率较大的薄膜展开，
因此本人提出了一种动态的展开方法，可以完全改善此中问题。
呵呵，不过正在等待发表，不便讨论

对于不可展曲面，尤其是大曲率面，仅仅基于边长变化最小也可能导致三角形角度变化很大，该方法缺少约束条件从而不可靠。
基于应变能变化最小的等效节点力动态展开可以解决该问题(如：五角锥展开视图.JPG)
影响展开应该还考虑膜材料的影响，就问题的优化而言，朋友是否可以谈谈考虑遗传算法的体会啊

上面图的单元应变能变化随迭代

再贴一张本人对展开/或者应力释放方法的一点思路的图
认为实际张拉得到的应力和形状和设计状态的差异可以通过2d  pattern shape得到改善