众所周知，在进行非线性有限元分析时有TL和UL两种列式，膜结构非线性有限元分析通常采用UL列式。
UL列式在求解过程中单元的节点坐标通常是在每一增量步结束更新一次，使其比较地接近变形状态。带动坐标系的迭代法实质上是一种U. L.法，在切线刚度阵中可略去大位移矩阵Ka，算例表明这样做较常规的U. L.法收敛快。在带有动坐标的的迭代法中，每迭代一次单元节点坐标就更新一次，使其紧紧地跟随单元变形状态。因此在计算不平衡力时其几何矩阵B可用其线性部分BL代替，（吕和祥等:《大转动梁的几何非线性分析讨论》，计算结构力学及其应用）。
对于某些讲膜结构的书中，通常在每一迭代步完后进行坐标更新，而不是在每一增量步结束更新一次，（比如杨庆山《张拉膜结构分析与设计》），因此其用的是拖动坐标法。那么在进行应力更新时，是不是用节点位移增量乘以线性BL就可以得到其应力增量。

进行应力更新时，是用节点位移增量乘以线性BL就可以得到其应力增量。
按照杨老师《张拉膜结构分析与设计》这本书的思路做即可.

我想说的是:杨老师的做法是不合适的,我曾经编程发现,
这会导致结果在很小的地方跳动,就是无法达到非常
小的容差,因此我建议你还是采用标准的非线性处理
方法,计算应力增量时候考虑非线性部分应变

无法达到非常 小的容差,我也遇到了这个问题，也许可能是找形时没有进行应力更新的缘故，我在找找原因。
欢迎多多交流！

怎么会不收敛呢?我好像没有这个问题,你所说的容差到底是多少?