/clear
finish
/filname,byb
x1=83.29
x2=84.42
x3=260
r=56
h=8
r1=64
r2=82
degree=8
width=682
length=1000
vel_p=-18800
vtime0=0
vtime1=2
vuprotate0=-vel_p/r1
vuprotate1=vuprotate0
vdownrotate0=vel_p/r2
vdownrotate1=vdownrotate0

/PREP7
ET,1,SHELL163
ET,2,SHELL163
ET,3,SHELL163
R, 1, 6/6, 4, h
R,2,6/6, 4,0.01
R,3,6/6, 4,0.01
keyopt,1,1,10
keyopt,2,1,10
keyopt,3,1,10
MP,DENS,1,7.8e-6
MP,EX,1,2.06e8
MP,NUXY,1,0.3
TB,BISO,1,,,,
TBDAT,1,2.45e5
TBDAT,2,5.8e5
EDMP,RIGI,2,7,5
MP,DENS,2,7.8e-6
MP,EX,2,2.06e8
MP,NUXY,2,0.3
EDMP,RIGI,2,7,5
MPCOPY, ,2,3
TBCOPY,ALL,2,3

*afun,deg

k,1,0,h/2
k,2,x1,h/2
k,3,0,-h/2
k,4,x2,-h/2
k,5,x2+x3*cos(degree),x3*sin(degree)-h/2
k,6,0,h/2+r1
k,7,x1,h/2+r1
k,8,0,-h/2-r2
k,9,x2,-h/2-r2
lstr,1,2
lstr,3,4
lstr,4,5
lfillt,2,3,r
arotat,1,,,,,,6,7
arotat,2,3,4,,,,8,9

wprota,,-90,0
rectng,0,width/2,-x3*sin(degree)-r2-20,-x3*sin(degree)-r2-20-length
wprota,,90,0

ESIZE,10

TYPE, 1
MAT, 1
REAL, 1
MSHAPE,0,2D
MSHKEY,1
AMESH,17

TYPE, 2
MAT, 2
REAL, 2
MSHAPE,0,2D
MSHKEY,1
ASEL,S,,,1,4,1
AMESH,ALL

TYPE,3
MAT,3
REAL,3
MSHAPE,0,2D
MSHKEY,1
ASEL,S,,,5,16,1
AMESH,ALL

EDPART,CREATE

EDCGEN,ASTS, 1, 2,_Z2,0.15,0,0,0, , , , ,0,10000000,0,0
EDCGEN,ASTS, 1, 3,_Z2,0.15,0,0,0, , , , ,0,10000000,0,0

finish
/solu

ASEL,S,,,17
NSLA,S,1
CM,slab,NODE
EDVE,VELO,SLAB,0,0,vel_p,0,0,0, , , , , ,

NSEL,R,LOC,X,0
D,ALL, , , , , ,UX,UY, , , ,
ASEL,S,,,17
NSLA,S,1
NSEL,R,LOC,X,0,60
NSEL,R,LOC,Z,-(-x3*sin(degree)-r2-20),-(-x3*sin(degree)-r2-20-20)
D,ALL,UY,0

ALLSEL,ALL

*afun,rad
*dim,vtime,array,2
*dim, vuprotate, array, 2
*dim, vdownrotate, array, 2
vtime(1)=vtime0,vtime1
vuprotate(1)=vuprotate0,vuprotate1
vdownrotate(1)=vdownrotate0,vdownrotate1

EDLOAD,ADD,RBRX,0, 2,VTIME,VUPROTATE, 0, , , , ,
EDLOAD,ADD,RBRX,0, 3,VTIME,VDOWNROTATE, 0, , , , ,

time,0.1
EDENERGY,1,1,1,1
EDOPT,ADD,blank,BOTH
EDRST,200,
EDHTIME,1000,
EDDUMP,1,
EDOUT,RCFORC

solve
finish

效果圖: