刚度矩阵非正定

hejt2004

大家好！
最近遇到以下问题：
*** STIFFNESS MATRIX NOT POSITIVE DEFINITE ***
NODE=      2935 EQUATION=      8886  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      2976 EQUATION=      9018  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      2982 EQUATION=      9045  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      3020 EQUATION=      9168  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      3027 EQUATION=      9192  DOF= Z-ROTATION     PIVOT= -3.33786011E-06
NODE=      3076 EQUATION=      9354  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      3084 EQUATION=      9381  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      3120 EQUATION=      9498  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      3203 EQUATION=      9756  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      3244 EQUATION=      9888  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      3258 EQUATION=      9945  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      3296 EQUATION=     10068  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      3311 EQUATION=     10122  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      3352 EQUATION=     10254  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      3360 EQUATION=     10281  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
NODE=      3396 EQUATION=     10398  DOF= Z-ROTATION     PIVOT= -3.81469727E-06
  以上节点都是梁元和块元公用的（为了使两者耦合，梁元伸到块元里面去，在深度方向3个节点公用），请帮忙分析分析原因。
  祝好！

loveala

块体单元的节点是没有转角自由度的，梁单元是有转角自由度的，即使两者共用了节点也是不行的，两者的耦合要用约束方程。

hejt2004

谢谢loveala回复
不太会约束方程
上面的方法是我问adina技术人员的
应该可以的