发现一个很严重的问题

rp_no

最近要用comsol做结构分析，涉及到接触。弄了很久，发现一个比接触收敛性更头疼的问题。
当添加接触后，geometric nonlinearity是默认的，所有其他区域全部按geometric nonlinear算法计算。而力学（不光是线弹性力学部分）从engineering strain变为Green-Lagrange strain，问题就出来了，gl strain理论上导致材料刚度变小，在-0.4的时候会导致刚度消失，于是就超级大变形，残差无法收敛。comsol blog里提到，当solid.ep3最小主应变小到-0.3的时候，基本计算就差不多要崩溃了。主要，这不是说大面积出现才崩溃，而是只要有一个点出现，就完蛋

blog还说，一些时候必须要使用geometric nonlinearity，比如hyperelastic, contact这样，另外某些几何结构导致singularity也会大应变。这些现象并不是不可接受，或者只是正常计算过程中出现的一些暂时的误差，但是当使用geometric nonlinearity时，却导致无法收敛。

调了半个月模型了，一直以为是接触不对，但现在发现，其实更大的问题是这个。接触通过重写接触应力公式，好歹一般可以保证收敛性。