刚看到一个做热模态的例子,不知道这APDL格式解释解释
6、 两端固结的梁,长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,-90, , , //温度
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,-90, , ,
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
不考虑预应力的结果 考虑预应力的结果
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
|