ansys导出的刚度矩阵负定(附图附程序请大神们帮忙看看)
在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 !求解
-------------------------------------------------------------------------------
桌面/QQ截图20161119111900.jpg superber 发表于 2016-11-19 11:22
这图怎么打不开 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 同样太大了。好像目前还没有这样的材料。
由于你的材料数据偏差太大，有可能是计算误差导致不合理的结果。
htbbzzg 发表于 2016-11-23 09:09
1 你的模型是一个 200*200 的方板吗。长度单位是什么？其它参数的单位又是什么？
2你的材料属性好像 ...
计算误差会有这么大的影响吗 本帖最后由 htbbzzg 于 2016-11-24 09:24 编辑
如果没有其它原因，就只能这样考虑了。
你先把材料数据修改完善再看吧。再多说一句，不知道你为什么要使用子结构分析，却又选择所有自由度为主自由度。如果使用子结构方法，应该只选择少量 (1/10 或更少) 自由度作为主自由度，否则的话，所需的资源比不使用子结构方法要大得多，是非常不合理的。
试算了你的命令流，结果与你所说完全不同。
如果没有约束，频率结果如下：
SET TIME/FREQ LOAD STEP SUBSTEPCUMULATIVE
1 0.0000 1 1 1
2 0.0000 1 2 2
30.80424E-04 1 3 3
40.40699E-03 1 4 4
50.98496E-03 1 5 5
60.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 发表于 2016-12-26 07:20
原来是个日本鬼子，很后悔回答了你的问题。
别闹损色纯纯中国人 损色{:{05}:}
页:
[1]