 |
|
|
非常感谢civil fans的答复！！！
下面这个命令流是我根据你上次贴出的命令流改编的，我想对此进行模态分析，却出现了问题，请你帮我看看，毛病出在那里
谢谢！！！
/filname,a
/PREP7
ET,1,SHELL41
KEYOPT,1,1,2
ET,2,LINK10
R,1,0.001, , , , ,,
RMORE, ,
R,2,0.000314,0.9999999,
MP,EX,1,8e6
MP,PRXY,1,0.3
MP,alpx,1,1
mp,dens,1,1
UIMP,1,REFT,,,0
MP,EX,2,2e8
MP,PRXY,2,0.3
K,1,,,,
K,2,15,,,
K,3,7.5,6,,
K,4,7.5,-6,,
!!!定义结构几何模型
LSTR, 1, 3
LSTR, 3, 2
LSTR, 2, 4
LSTR, 4, 1
AL,all
AATT, 1, 1, 1, 0,
LATT,2,2,2, , , ,
LESIZE,ALL, , ,15, ,1, , ,1, !!!定义单元边界长度
MSHAPE,1,2D
MSHKEY,1
AMESH,all !!!三角划分面单元
!!!划分索单元
LMESH,ALL
FINISH
!!!设定边界并且第一次求解
/SOLU
DK,1,ALL
DK,2,ALL
DK,3,ALL
DK,4,ALL
NLGEOM,on
NSUBST,20,
AUTOTS,on
TIME,1
LUMPM,0
EQSLV, , ,0,
PRECISION,0
MSAVE,0
PIVCHECK,1
sstif,on
DK,1,UZ,5
DK,2,UZ,5
SOLVE
EPLOT
UPCOORD,1,OFF
DK,1,ALL
DK,2,ALL
FINISH
/PREP7
MP,EX,2,2E+11
R,2,0.000314,0.9999999e-3,
FINISH
/SOL
BFA,all,TEMP,-0.05
SOLVE
EPLOT
FINISH
/POST1
PLDISP,0
/VIEW, 1 ,,1
/REP,FAST
FINISH
!!!第一次循环更新坐标，释放不平衡应力
n=15
*do,i,1,n
/solu
antype,static
sstif,on
/solu
solve
UPCOORD,0.5,off,a,rst
*enddo
!!!施加真实应力
/SOL
BFA,all,TEMP,-2
SOLVE
FINISH
/POST1
AVPRIN,0, ,
PLESOL,S,EQV,0,1
FINISH
!!!第二次循环更新坐标，释放不平衡应力
n=5
*do,i,1,n
/solu
antype,static
sstif,on
/solu
solve
upcoord,0.5,off,a,rst
*enddo
!!!还原膜的真实弹性模量
/PREP7
MP,EX,1,8E+8
FINISH
/SOL
antype,static
pstress,on
BFA,all,TEMP,-2e-2
SOLVE
FINISH
/POST1
AVPRIN,0, ,
PLESOL,S,EQV,0,1
!!!第三次循环更新坐标，释放不平衡应力,找形结束
n=3
*do,i,1,n
/solu
antype,static
pstres,on
sstif,on
/solu
solve
upcoord,0.5,off,a,rst
*enddo
fini
/solu
antype,modal
pstres,on
modopt,subsp,3
mxpand,5
allsel,
outpr,all,all,,
/OUT,mode,vec
solve
fini
/post1
set,list,2
set,first
pldisp,2
fini
|
