我做的是一个2D隧道,初衬的内力计算模型,命令流如下,最后绘制出的内力图总觉得有问题,希望高手指点
*set,x1,-25
*set,y1,-24.3
*set,w,50
*set,h,60
*set,th,0.5
/prep7
NUMSTR,KP,1
NUMSTR,LINE,1
NUMSTR,AREA,1
NUMSTR,VOLU,1
CSYS,0
k, 1 ,5.700,.000,.000,
k, 2 ,3.189,4.725,.000,
k, 3 ,5.531,1.379,.000,
LARC,1 ,2 ,3
k, 4 ,-2.500,.000,.000,
k, 5 ,5.533,-1.645,.000,
k, 6 ,5.606,-1.237,.000,
LARC,5 ,1 ,6
k, 7 ,4.554,-1.444,.000,
k, 8 ,4.879,-2.390,.000,
k, 9 ,5.111,-2.275,.000,
LARC,8 ,5 ,9
k, 10 ,.000,11.794,.000,
k, 11 ,.000,-3.206,.000,
k, 12 ,1.241,-3.154,.000,
LARC,11 ,8 ,12
k, 13 ,2.500,.000,.000,
k, 14 ,-5.700,.000,.000,
k, 15 ,-5.533,-1.645,.000,
k, 16 ,-5.690,-.414,.000,
LARC,14 ,15 ,16
k, 17 ,-4.554,-1.444,.000,
k, 18 ,-4.879,-2.390,.000,
k, 19 ,-5.449,-1.890,.000,
LARC,15 ,18 ,19
k, 20 ,-3.689,-2.745,.000,
LARC,18 ,11 ,20
k, 21 ,.000,.000,.000,
k, 22 ,-3.189,4.725,.000,
k, 23 ,1.667,5.451,.000,
LARC,2 ,22 ,23
k, 24 ,-4.237,3.813,.000,
LARC,22 ,14 ,24
lsel,s,line,,1,2
lsel,a,line,,5,
lsel,a,line,,8,9
lcomb,all,,
alls
lsel,s,line,,3,4
lsel,a,line,,6,7
lcomb,all,,
alls
al,1,3
blc4,x1,y1,w,h
aovlap,2,1
et,1,beam3
et,2,plane42
keyopt,2,3,2
r,1,th,th*th*th/12,th, ,
mp,ex,1,2.6e7
mp,prxy,1,0.2
mp,dens,1,25
mp,ex,2,1.3e6
mp,prxy,2,0.32
tb,dp,2
tbdata,1,200,30,
mp,dens,2,22
mp,ex,3,1.3e6
mp,prxy,3,0.32
tb,dp,3
tbdata,1,200,30,
mp,dens,3,22
lsel,s,line,,1
lesize,all,,,50,
latt,1,1,1
lmesh,all
alls
lsel,s,line,,3
lesize,all,,,20,
latt,1,1,1
lmesh,all
alls
lsel,s,line,,5
lsel,a,line,,2
lesize,all,,,50
alls
lsel,s,line,,6
lsel,a,line,,4
lesize,all,,,40
alls
mopt,split,on
mopt,split,err
mopt,split,warn
mopt,qmesh,main
asel,s,area,,3
aatt,2,,2
amesh,all
alls
asel,s,area,,1
aatt,3,,2
amesh,all
alls
lsel,s,loc,x,x1
lsel,a,loc,x,x1+w
dl,all,,ux,0
alls
lsel,s,loc,y,y1
dl,all,,uy,0
alls
acel,,1
fini
/solu
antype,static
deltim,0.1,0.05,0.2
autots,on !使用自动时间步
pred,on !打开时间步长预测器
lnsrch,on !打开线性搜索
nlgeom,on !打开大位移效果
nropt,full !设定牛顿-拉普森选项
esel,s,type,,1
ekill,all
esel,all
esel,s,live
nsel,all
esel,all
solve
esel,s,mat,,3
ekill,all
esel,s,type,,1
ealive,all
nsle,s
ddele,all,all
esel,all
esel,s,live
nsle,s
nsel,inve
d,all,all,0
nsel,all
esel,all
solve
fini
/post1
etable,if,smisc,1
etable,jf,smisc,7
etable,im,smisc,6
etable,jm,smisc,12
etable,ism,nmisc,1
etable,jsm,nmisc,3
etable,isn,nmisc,2
etabel,jsn,nmisc,4 |