LENG_OVERALL = 115
LENG_CENTER  = 33
THICK      = 1.5
WID_CENTER   = 6
RAD_LARGE   = 25/2
RAD_SMALL   = 14
WID_OVERALL  = 2*RAD_LARGE
*afun,deg
/prep7
et,1,185
keyopt,1,6,1
r,1
TB,HYPE,1,1,3,YEOH  
TBTEMP,0
TBDATA,,0.6976,-0.24484,0.12629,,,
rectng,,LENG_CENTER/2,,WID_CENTER/2
asel,none
rectng,,LENG_OVERALL/2,,WID_OVERALL/2
!
wpoffs,LENG_CENTER/2,WID_CENTER/2+RAD_SMALL
pcirc,RAD_SMALL,,0,-90
cswpla,11,1
lsel,s,loc,x,RAD_SMALL
TEMPLINE1=lsnext(0)
!
wpcsys,-1,0
LENG_OFFSET=sqrt((RAD_SMALL+RAD_LARGE)**2-(RAD_SMALL+WID_CENTER/2)**2)
wpoffs,LENG_CENTER/2+LENG_OFFSET
pcirc,RAD_LARGE,,90,180
cswpla,12,1
lsel,s,loc,x,RAD_LARGE
TEMPLINE2=lsnext(0)
adele,all
lsel,all
lovlap,TEMPLINE1,TEMPLINE2
!
csys,11
lsel,s,loc,y,270,270+atan(LENG_OFFSET/(RAD_SMALL+WID_CENTER/2))
lsel,r,loc,x,RAD_SMALL
TEMPLINE1=lsnext(0)
!
csys,12
lsel,s,loc,y,90,180-atan((RAD_SMALL+WID_CENTER/2)/LENG_OFFSET)
lsel,r,loc,x,RAD_LARGE
TEMPLINE2=lsnext(0)
!
lsel,all
csys,0
asel,all
lsla,u
lsel,u,loc,x,LENG_OVERALL/2
lsel,u,line,,TEMPLINE1
lsel,u,line,,TEMPLINE2
ldele,all,,,1
lsel,all
nummrg,kp

a,2,6,7,12,15,3
wpcsys,-1,0
kwpave,6
rectng,,-RAD_LARGE,,WID_OVERALL/2
aovlap,all

mshape,0
mshkey,0
esize,WID_CENTER/5
vext,all,,,0,0,THICK/2
vsweep,all
asel,s,loc,x,0
da,all,ux
!
asel,s,loc,y,0
da,all,uy
!
asel,s,loc,z,0
da,all,uz
da,9,ux,0.58*LENG_CENTER
allsel,all
wpstyl,,,,,,,,0
rescon,define,none
finish

/solu
nlgeom,on
outres,all,all
nsubst,20,1e3,20

allsel,all
solve
finish
