ansys导出的刚度矩阵负定(附图附程序请大神们帮忙看看)

superber

在ansys中建模 施加约束发现刚度矩阵是负定的，请大神们帮忙看一下！！！
没约束时，总刚度矩阵有六个负特征值；施加了约束后还有一个负特征值；没想出来约束有什么问题！

这是ansys apdl程序：
--------------------------------------------------------------
finish
/clear
/cwd,E:\juzhen
/FILNAME, EXAMPLE1021
/prep7
density=7.15e-6
et,1,shell181
et,2,beam188
mp,ex,1,1.12e8
mp,prxy,1,0.25
mp,dens,1,density
r,1,0.8
k,1,0,0,0
k,2,200,0,0
k,3,200,200,0
k,4,0,200,0
a,1,2,3,4,1

sectype,1,shell
secdata,0.2
secoffset,mid
aatt,1,1,1,0,1
allsel
aesize,all,10
amesh,all
allsel

nsel,s,loc,x,0
d,all,,,,,,all,,,,,
allsel
nsel,s,loc,y,0
d,all,,,,,,all,,,,,
allsel
finish

/solu
antype,modal
modopt,lanb,5
mxpand,5
solve

/OUTPUT,juzhen,txt,,
/solu
antype,7 !substructuring分析类型
seopt,matname,2,1 !设置文件名称和刚度矩阵类型(刚度，质量，阻尼等)
nsel,all !选择所有节点
m,all,all !定义所有节点自由度为主自由度
solve !求解
-------------------------------------------------------------------------------

superber

                               
登录/注册后可看大图

怪咖先生

superber 发表于 2016-11-19 11:22

这图怎么打不开

htbbzzg

  1 你的模型是一个 200*200 的方板吗。长度单位是什么？其它参数的单位又是什么？
  2  你的材料属性好像有问题：
    如果你的长度单位是 m，那么弹性模量 1.12e8 Pa = 112 MPa，好像太小啦；而密度 7.15e-6 kg/m^3 同样太小了，而且是小了几个数量级。
    如果你的长度单位是 mm，那么弹性模量 1.12e8 MPa ，又太大了，而且是大了几个数量级；而密度 7.15e-6 T/mm^3 同样太大了。好像目前还没有这样的材料。
    由于你的材料数据偏差太大，有可能是计算误差导致不合理的结果。
   

think2015

htbbzzg 发表于 2016-11-23 09:09
1 你的模型是一个 200*200 的方板吗。长度单位是什么？其它参数的单位又是什么？
  2  你的材料属性好像 ...

计算误差会有这么大的影响吗

htbbzzg

  如果没有其它原因，就只能这样考虑了。
  你先把材料数据修改完善再看吧。  再多说一句，不知道你为什么要使用子结构分析，却又选择所有自由度为主自由度。如果使用子结构方法，应该只选择少量 (1/10 或更少) 自由度作为主自由度，否则的话，所需的资源比不使用子结构方法要大得多，是非常不合理的。

htbbzzg

  试算了你的命令流，结果与你所说完全不同。
  如果没有约束，频率结果如下：
   SET   TIME/FREQ    LOAD STEP   SUBSTEP  CUMULATIVE
     1      0.0000             1         1         1
     2      0.0000             1         2         2
     3  0.80424E-04         1         3         3
     4  0.40699E-03         1         4         4
     5  0.98496E-03         1         5         5
     6  0.10965E-02         1         6         6
     7      13.071             1         7         7
     8      18.963             1         8         8
     9      22.653             1         9         9
    10     33.542             1        10        10
  没有负特征值。
  如果按照你的命令流中的约束条件，计算结果没有零频率，前三阶模态形状如下：



  不过，如果按照你的命令流计算到底，则由于最后一个 solve 是子结构求解，前面的结果被覆盖，就没有结果了。
  估计你的 ANSYS 有问题，换一个版本试试吧。

野菊

htbbzzg 发表于 2016-11-24 10:34
试算了你的命令流，结果与你所说完全不同。
  如果没有约束，频率结果如下：
   SET   TIME/FREQ    LO ...

搜噶

htbbzzg

原来是个日本鬼子，很后悔回答了你的问题。

野菊

htbbzzg 发表于 2016-12-26 07:20
原来是个日本鬼子，很后悔回答了你的问题。

别闹  损色  纯纯中国人

realyyy

损色