刚度矩阵非正定
大家好！最近遇到以下问题：
*** STIFFNESS MATRIX NOT POSITIVE DEFINITE ***
NODE= 2935 EQUATION= 8886DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 2976 EQUATION= 9018DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 2982 EQUATION= 9045DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 3020 EQUATION= 9168DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 3027 EQUATION= 9192DOF= Z-ROTATION PIVOT= -3.33786011E-06
NODE= 3076 EQUATION= 9354DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 3084 EQUATION= 9381DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 3120 EQUATION= 9498DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 3203 EQUATION= 9756DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 3244 EQUATION= 9888DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 3258 EQUATION= 9945DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 3296 EQUATION= 10068DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 3311 EQUATION= 10122DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 3352 EQUATION= 10254DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 3360 EQUATION= 10281DOF= Z-ROTATION PIVOT= -3.81469727E-06
NODE= 3396 EQUATION= 10398DOF= Z-ROTATION PIVOT= -3.81469727E-06
以上节点都是梁元和块元公用的（为了使两者耦合，梁元伸到块元里面去，在深度方向3个节点公用），请帮忙分析分析原因。
祝好！ 块体单元的节点是没有转角自由度的，梁单元是有转角自由度的，即使两者共用了节点也是不行的，两者的耦合要用约束方程。 谢谢loveala回复
不太会约束方程
上面的方法是我问adina技术人员的
应该可以的
页:
[1]