finish
/clear,start
/title,niuju shi jia5
/prep7
cyl4,0,0,20,0,40,360,200
wpstyl,,5,0,100,,,0,1
wprota,,,90
vsbw,all
wprota,,90
vsbw,all
wpcsys,-1,0
wpstyl,,5,0,100,,,1,0
mp, ex,1,2e5
mp,prxy,1,0.3
et,1,solid185
et,2,mpc184
keyopt,2,1,1 !刚性梁
lesize,30,,,10
lesize,32,,,10
lesize,13,,,15
lesize,14,,,15
lesize,15,,,15
lesize,16,,,15
lesize,18,,,50
mshape,0,3d
mshkey,1
vmesh,all
nummrg,all,,,,low
numcmp,all
!.........................................
nsel,s,loc,z,200
n,100000,0,0,220
!*get,node_num,node,,count
!*status,parm
eplot
!............................生成梁单元
type,2 !千万不要忘记选择单元了
e,100000,2
*do,i,1,10
e,100000,51+i
*enddo
*do,i,1,15
e,100000,726+i
*enddo
*do,i,1,10
e,100000,1476+i
*enddo
*do,i,1,14
e,100000,1976+i
*enddo
*do,i,1,126
e,100000,2676+i
*enddo
*do,i,1,165
e,100000,8976+i
*enddo
*do,i,1,11
e,100000,17391+i
*enddo
*do,i,1,14
e,100000,17966+i
*enddo
*do,i,1,140
e,100000,18666+i
*enddo
*do,i,1,154
e,100000,25960+i
*enddo
allsel,all,all
!.........................................
/solu
nsel,s,loc,z,0
d,all,all
allsel
f,100000,mz,2000
!d,100000,rotz,0.7854 !施加转角位移,大小为45度,即0.7854弧度
nsubst,10
EQSLV,PCG,1E-8
nlgeom,on !打开大变形选项
solve