声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1769|回复: 10

[综合讨论] ansys导出的刚度矩阵负定(附图附程序请大神们帮忙看看)

[复制链接]
发表于 2016-11-19 11:20 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
在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 !求解
-------------------------------------------------------------------------------
回复
分享到:

使用道具 举报

 楼主| 发表于 2016-11-19 11:22 | 显示全部楼层

                               
登录/注册后可看大图

点评

这图怎么打不开  详情 回复 发表于 2016-11-21 08:43
发表于 2016-11-21 08:43 | 显示全部楼层

回帖奖励 +1 点体能


这图怎么打不开
发表于 2016-11-23 09:09 | 显示全部楼层
  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 同样太大了。好像目前还没有这样的材料。
    由于你的材料数据偏差太大,有可能是计算误差导致不合理的结果。
   

点评

计算误差会有这么大的影响吗  详情 回复 发表于 2016-11-23 09:52
发表于 2016-11-23 09:52 | 显示全部楼层
htbbzzg 发表于 2016-11-23 09:09
1 你的模型是一个 200*200 的方板吗。长度单位是什么?其它参数的单位又是什么?
  2  你的材料属性好像 ...

计算误差会有这么大的影响吗
发表于 2016-11-24 09:19 | 显示全部楼层
本帖最后由 htbbzzg 于 2016-11-24 09:24 编辑

  如果没有其它原因,就只能这样考虑了。
  你先把材料数据修改完善再看吧。  再多说一句,不知道你为什么要使用子结构分析,却又选择所有自由度为主自由度。如果使用子结构方法,应该只选择少量 (1/10 或更少) 自由度作为主自由度,否则的话,所需的资源比不使用子结构方法要大得多,是非常不合理的。
发表于 2016-11-24 10:34 | 显示全部楼层
  试算了你的命令流,结果与你所说完全不同。
  如果没有约束,频率结果如下:
   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
  没有负特征值。
  如果按照你的命令流中的约束条件,计算结果没有零频率,前三阶模态形状如下:
EXAMPLE1021001.jpg
EXAMPLE1021002.jpg
EXAMPLE1021003.jpg
  不过,如果按照你的命令流计算到底,则由于最后一个 solve 是子结构求解,前面的结果被覆盖,就没有结果了。
  估计你的 ANSYS 有问题,换一个版本试试吧。

点评

搜噶  详情 回复 发表于 2016-12-21 13:10
发表于 2016-12-21 13:10 | 显示全部楼层
htbbzzg 发表于 2016-11-24 10:34
试算了你的命令流,结果与你所说完全不同。
  如果没有约束,频率结果如下:
   SET   TIME/FREQ    LO ...

搜噶
发表于 2016-12-26 07:20 | 显示全部楼层
原来是个日本鬼子,很后悔回答了你的问题。

点评

别闹 损色 纯纯中国人  详情 回复 发表于 2016-12-26 08:39
发表于 2016-12-26 08:39 | 显示全部楼层
htbbzzg 发表于 2016-12-26 07:20
原来是个日本鬼子,很后悔回答了你的问题。

别闹  损色  纯纯中国人
发表于 2017-1-2 14:56 | 显示全部楼层
损色
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-24 22:27 , Processed in 0.073979 second(s), 23 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表