ansys先静力分析后模态分析的程序
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
考虑热应力的模态分析 (2007-06-18 16:49:28)
标签： cae fea 模态有限元分析热应力
愿与大家分享，共同进步。
当然也欢迎大家把自己的想法与我交流。
考虑预应力影响的模态分析的方法如下：
1。
先进行静力分析，注意分析时打开预应力开关(PSTRES,ON)。
2。
改变分析类型，获取模态解。
(再用一次PSTRES,ON)
3。
特别注意：a.静力分析中和随后的模态分析中的集中质量选项(LUMPM)必须一致。
b.进行模态分析时，应保证静力分析的保存的.emat和.esav文件存在。
c.步骤1也可以是瞬态分析，但应当在需要的时间保存.emat和.esave文件。
下面是自己做的一个小例子。
两端固结
的梁，长10m，截面为0.1m*0.1m,材质为钢材。
对比降温90度前后的模态分析解。
!考虑温度预应力的模态分析
FINI
/CLEAR,NOSTART
/prep7
k,1,0,0,0
k,2,10,0,0
l,1,2
et,1,beam4
MP,EX,1,200e9
MP,NUXY,1,0.3
mp,alpx,1,0.000012
MP,dens,1,7800
R,1,0.01,8.3E-6,8.3E-6,0.1,0.1, , //定义常系数
LESIZE,ALL, , ,10, , , , ,1 //定义线上的单元数
LMESH,ALL //划分线生成线单元
/SOLU //进入求解器单元
ANTYPE,STATIC //定义分析类型
PSTRES,ON //是否计入预应力
lumpm,off
NSEL,S,LOC,X,0
D,all,all //施加约束ALLSEL,ALL //选择所有实体
NSEL,S,LOC,X,10
D,ALL,UY
D,all,all
esel,all
acel,,-9.8, //定义结构线性加速度BFE,all,TEMP,1,+70, , ,
ALLSEL,ALL
SOLVE
SAVE
FINI
/SOLU
ANTYPE,MODAL //选择分析类型
MODOPT,LANB,10 //模态分析选项
EQSLV,SPAR
MXPAND,0, , ,0
lumpm,off
PSTRES,1
MODOPT,LANB,10,0,500, ,OFF //模态分析选项和模态分析方法NSEL,S,LOC,X,0
D, all,all
ALLSEL,ALL
NSEL,S,LOC,X,10 //节点的选择
D, all,all
ALLSEL,ALL //选择所有实体
SOLVE
SAVE
/POST1
SET,LIST
!不考虑温度预应力
FINI
/CLEAR,NOSTART
/prep7
k,1,0,0,0
k,2,10,0,0
l,1,2
et,1,beam4
MP,EX,1,200e9
MP,NUXY,1,0.3
mp,alpx,1,0.000012
MP,dens,1,7800
R,1,0.01,8.3E-6,8.3E-6,0.1,0.1, , //定义常系数LESIZE,ALL, , ,10, , , , ,1 //定义线上的单元数LMESH,ALL //划分线生成线单元
/SOLU //进入求解器单元ANTYPE,STATIC //定义分析类型PSTRES,ON //是否计入预应力lumpm,off
NSEL,S,LOC,X,0
D,all,all //施加约束
ALLSEL,ALL //选择所有实体
NSEL,S,LOC,X,10
D,ALL,UY
D,all,all
esel,all
acel,,-9.8, //定义结构线性加速度
!BFE,all,TEMP,1, +70, , ,
ALLSEL,ALL
SOLVE
SAVE
FINI
/SOLU //进入求解器单元
ANTYPE,MODAL //分析类型
MODOPT,LANB,10 //模态分析选项
EQSLV,SPAR
MXPAND,0, , ,0
lumpm,off // Use the element-dependent default mass matrix formulation (default). PSTRES,1 //是否计入预应力
MODOPT,LANB,10,0,500, ,OFF //模态分析选项和模态分析方法
NSEL,S,LOC,X,0
D, all,all //施加约束ALLSEL,ALL //选择所有实体NSEL,S,LOC,X,10 //节点的选择
D, all,all //施加约束
ALLSEL,ALL //选择所有实体
SOLVE
SAVE
/POST1 //进入通用后处理器SET,LIST
不考虑预应力的结果考虑预应力的结果
SET TIME/FREQ SET TIME/FREQ
1 5.1946 1 10.446
2 5.1946 2 10.446
3 14.320 3 22.833
4 14.320 4 22.833
5 28.088 5 38.449
6 28.088 6 38.449
7 46.493 7 57.983
8 46.493 8 57.983
9 69.639 9 81.852
10 69.639 10 81.852。