R_CENTER=10                                       
R_RADIUS=2                                        
C_OFFSET=1.0                                      
/prep7                                            
et,1,182                                          
keyopt,1,3,1                                      
!keyopt,1,6,1                                     
et,2,169                                          
et,3,171                                          
r,1                                               
r,2                                               
r,3                                               
r,4                                               
rmodif,2,3,2                                      
rmodif,3,3,2                                      
rmodif,4,3,2                                      
mp,mu ,1,0.1                                      
tb,hyper,1,1,3,ogden                              
tbtemp,0                                          
tbdata,1,6.3  ,1.3                                
tbdata,3,0.012,5.0                                
tbdata,5,-0.1 ,-2.0                               
tbdata,7,2e-4                                     
k,1,R_CENTER+R_RADIUS,         -C_OFFSET          
k,2,R_CENTER+R_RADIUS,         2*R_RADIUS+C_OFFSET
k,3,R_CENTER+R_RADIUS+C_OFFSET,2*R_RADIUS         
k,4,R_CENTER-R_RADIUS-C_OFFSET,2*R_RADIUS         
k,5,R_CENTER-R_RADIUS,         2*R_RADIUS+C_OFFSET
k,6,R_CENTER-R_RADIUS,         -C_OFFSET          
l,2,1                                             
latt,1,2,2                                        
lsel,none                                         
l,4,3                                             
latt,1,3,2                                        
lsel,none                                         
l,6,5                                             
latt,1,4,2                                        
wpoffs,R_CENTER,R_RADIUS                          
pcirc ,R_RADIUS,,0,90                             
aatt,1,1,1                                        
esize,0.3                                         
mshape,0                                          
mshkey,1                                          
amesh,all                                         
lsel,s,ext                                        
csys,4                                            
lsel,u,loc,y,0                                    
lsel,u,loc,x,0                                    
!lrefine,all                                      
!lrefine,all                                      
arsym,x,all                                       
arsym,y,all                                       
nummrg,node                                       
nummrg,kp                                         
lsel,r,loc,x,0,R_RADIUS                           
type,3                                            
real,2                                            
mat,1                                             
nsll,s,1                                          
esurf                                             
lsel,s,ext                                        
lsel,r,loc,y,0,R_RADIUS                           
type,3$real,3$mat,1                               
nsll,s,1                                          
esurf                                             
lsel,s,ext                                        
lsel,r,loc,x,0,-R_RADIUS                          
type,3$real,4$mat,1                               
nsll,s,1                                          
esurf                                             
csys,0                                            
lsel,all                                          
lsla,u                                            
lmesh,all                                         
ksel,s,kp,,1                                      
tshap,pilo                                        
katt,1,2,2                                        
kmesh,all                                         
allsel,all                                        
finish                                            
/solu                                             
antype,static                                     
nlgeom,on                                         
rescon,define,none                                
outres,all,all                                    
nsubst,50,1e5,20,off                              
ksel,s,kp,,1                                      
nslk                                              
d,all,ux,-R_RADIUS/3                              
allsel,all                                        
solve                                             
!  nsubst,50,1e5,10,off                           
lsel,s,line,,10,13,3                              
sfl,all,pres,20                                   
allsel,all                                        
solve                                             
finish                                            
wpsty,,,,,,,,0                                    
/post1                                            
set,last                                          
ples,nl,hpres                                     
finish                                            
                                                  
eplot                                             