finish
/CLEAR,nostart

/prep7

d=0.1

!Geometrie
t=2
r=50.0
l=65.0
r2=r+t

!Mesh
meshdens=1
nip=3

!=============== define element types
ET,2,SOLID185  

MP,EX,1,70000        ! MPa
MP,nuXY,1,0.30
MP,dens,1,2.78e-09  !g/cm³
MP,mu,1,0.5

!************* Inner cylinder
local,11,1,0,0,0
MSHKEY,1	!Specifies whether free meshing (0) or mapped meshing (1) should be used to mesh a model. (2) Use mapped meshing if possible; otherwise, use free meshing
CYLIND,r,r+t+0.1,0,l,0,180
csys,11
esys,11
LESIZE,1,,, 1.5
ESIZE,meshdens
VATT,1,,2
vmesh,1
esel,all
cm,vol1_elm,elem
nsel,all
cm,vol1_nodes,node
esel,none
nsel,none

! Symmetry Nodes

csys,0
cmsel,s,vol1_nodes
nsel,r,loc,x,r-t,r+t+0.2
nsel,r,loc,y,0,0.1
cm,sym_nodes1,node
nsel,none

cmsel,s,vol1_nodes
nsel,r,loc,x,-r+t,-r-t-0.2
nsel,r,loc,y,0,0.1
cm,sym_nodes2,node
nsel,none

! Radial Nodes

csys,11
cmsel,s,vol1_nodes
cmsel,u,sym_nodes1
cmsel,u,sym_nodes2
nsel,r,loc,x,r+t,r+t+0.1
cm,rad_inner_cyl_outer_nodes,node
nsel,none

cmsel,s,vol1_nodes
cmsel,u,sym_nodes1
cmsel,u,sym_nodes2
nsel,r,loc,x,r,r+0.1
cm,rad_inner_cyl_inner_nodes,node
nsel,none

! Edge Nodes

cmsel,s,sym_nodes1
nsel,r,loc,x,r,r+0.1
cm,edge_nodes1,node
nsel,none

cmsel,s,sym_nodes1
nsel,r,loc,x,r+t,r+t+0.1
cm,edge_nodes2,node
nsel,none

cmsel,s,sym_nodes2
nsel,r,loc,x,r,r+0.1
cm,edge_nodes3,node
nsel,none

cmsel,s,sym_nodes2
nsel,r,loc,x,r+t,r+t+0.1
cm,edge_nodes4,node
nsel,none

nsel,all
csys,11
nrotat,all

!************* Outer cylinder

MSHKEY,1	!Specifies whether free meshing (0) or mapped meshing (1) should be used to mesh a model. (2) Use mapped meshing if possible; otherwise, use free meshing
CYLIND,r2,r2+t,0,l+30,0,180
csys,11
esys,11
LESIZE,13,,, 1.5
ESIZE,meshdens
VATT,1,,2
vmesh,2
esel,all
cmsel,u,vol1_elm
cm,vol2_elm,elem
nsel,all
cmsel,u,vol1_nodes
cm,vol2_nodes,node
esel,none
nsel,none

! Symmetry Nodes

csys,0
cmsel,s,vol2_nodes
nsel,r,loc,x,r2-t,r2+t+0.2
nsel,r,loc,y,0,0.1
cm,sym_nodes3,node
nsel,none

cmsel,s,vol2_nodes
nsel,r,loc,x,-r2+t,-r2-t-0.2
nsel,r,loc,y,0,0.1
cm,sym_nodes4,node
nsel,none

! Fixation Nodes

cmsel,s,vol2_nodes
nsel,r,loc,z,l+30,l+29.90
cm,fixation,nodes
nsel,none

! Radial Nodes

csys,11
cmsel,s,vol2_nodes
cmsel,u,sym_nodes3
cmsel,u,sym_nodes4
nsel,r,loc,x,r2+t,r2+t+0.1
cm,rad_outer_cyl_outer_nodes,node
nsel,none

cmsel,s,vol2_nodes
cmsel,u,sym_nodes3
cmsel,u,sym_nodes4
nsel,r,loc,x,r2,r2+0.1
cm,rad_outer_cyl_inner_nodes,node
nsel,none

! Edge Nodes

cmsel,s,sym_nodes3
nsel,r,loc,x,r2,r2+0.1
cm,edge_nodes5,node
nsel,none

cmsel,s,sym_nodes3
nsel,r,loc,x,r2+t,r2+t+0.1
cm,edge_nodes6,node
nsel,none

cmsel,s,sym_nodes4
nsel,r,loc,x,r2,r2+0.1
cm,edge_nodes7,node
nsel,none

cmsel,s,sym_nodes4
nsel,r,loc,x,r2+t,r2+t+0.1
cm,edge_nodes8,node
nsel,none

! Middle Nodes for BCs

csys,0
nsel,s,loc,x,-0.3,0.3
cm,middle_nodes,node
nsel,none

! Contact Nodes

csys,11
nsel,none
cmsel,s,rad_outer_cyl_inner_nodes
cmsel,a,edge_nodes5
cmsel,a,edge_nodes7
cm,outer_cyl_contact_nodes,node
nsel,none

cmsel,s,rad_inner_cyl_outer_nodes
cmsel,a,edge_nodes2
cmsel,a,edge_nodes4
cm,inner_cyl_contact_nodes,node
nsel,none

!************* contact cylinder

MSHKEY,1	!Specifies whether free meshing (0) or mapped meshing (1) should be used to mesh a model. (2) Use mapped meshing if possible; otherwise, use free meshing
CYLIND,r2+t,r2+2*t,0,l+30,0,180
csys,11
esys,11
LESIZE,25,,, 1.5
ESIZE,meshdens
VATT,1,,2
vmesh,3
esel,all
cmsel,u,vol1_elm
cmsel,u,vol2_elm
cm,vol3_elm,elem
nsel,all
cmsel,u,vol1_nodes
cmsel,u,vol2_nodes
cm,vol3_nodes,node
esel,none
nsel,none

! Radial Nodes

csys,11
cmsel,s,vol3_nodes
nsel,r,loc,x,r2+2*t,r2+2*t+0.1
cm,rad_cont_cyl_outer_nodes,node
nsel,none

cmsel,s,vol3_nodes
nsel,r,loc,x,r2+t,r2+t+0.1
cm,rad_cont_cyl_inner_nodes,node
nsel,none

nsel,all
csys,11
nrotat,all

/COM, CONTACT PAIR CREATION - START 
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
/GSAV,cwz,gsav,,temp
MP,MU,1,0   
MAT,1   
MP,EMIS,1,7.88860905221e-031
R,3 
REAL,3  
ET,3,170
ET,4,174
R,3,,,1.0,0.01,0,   
RMORE,,,1.0E20,0.0,1.0, 
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
RMORE,10.0  
KEYOPT,4,4,2
KEYOPT,4,5,0
KEYOPT,4,7,0
KEYOPT,4,8,0
KEYOPT,4,9,0
KEYOPT,4,10,2   
KEYOPT,4,11,0   
KEYOPT,4,12,6   
KEYOPT,4,2,2
KEYOPT,3,5,0
! Generate the target surface   
NSEL,S,,,RAD_OUTER_CYL_OUTER_NODES  
CM,_TARGET,NODE 
TYPE,3  
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM 
! Generate the contact surface  
NSEL,S,,,RAD_CONT_CYL_INNER_NODES   
CM,_CONTACT,NODE
TYPE,4  
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,3  
ESEL,A,TYPE,,4  
ESEL,R,REAL,,3  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,3  
ESEL,A,TYPE,,4  
ESEL,R,REAL,,3  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END 

! CONTACT

/COM, CONTACT PAIR CREATION - START 
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
/GSAV,cwz,gsav,,temp
MP,MU,1,0.5 
MAT,1   
MP,EMIS,1,7.88860905221e-031
R,3 
REAL,3  
ET,3,170
ET,4,174
R,3,,,1.0,0.001,0,   ! Penetration Tolerance
RMORE,,,1.0E20,0.0,1.0, 
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,4,4,0
KEYOPT,4,5,0
KEYOPT,4,7,0
KEYOPT,4,8,0
KEYOPT,4,9,2
KEYOPT,4,10,2   
KEYOPT,4,11,0   
KEYOPT,4,12,0   
KEYOPT,4,2,1
KEYOPT,3,5,0
! Generate the target surface   
NSEL,S,,,OUTER_CYL_CONTACT_NODES
CM,_TARGET,NODE 
TYPE,3  
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM 
! Generate the contact surface  
NSEL,S,,,INNER_CYL_CONTACT_NODES
CM,_CONTACT,NODE
TYPE,4  
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,3  
ESEL,A,TYPE,,4  
ESEL,R,REAL,,3  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,3  
ESEL,A,TYPE,,4  
ESEL,R,REAL,,3  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END
FINISH
/SOLU

! Static Analysis "Fügen"

ANTYPE,STATIC
NLGEOM,OFF
AUTOTS,ON
DELTIM,0.1,0.05,1.0
TIME,1
pstres,on

NLHIST,pair,Spalt,cont,pene,3
NLHIST,pair,Druck,cont,pres,3
NLHIST,pair,sliding,cont, slid,3

csys,11
nsel,none

cmsel,s,sym_nodes1
cmsel,a,sym_nodes2
cmsel,a,sym_nodes3
cmsel,a,sym_nodes4
DSYM,SYMM,Y,11

allsel

D,sym_nodes1,uz,0
D,sym_nodes2,uz,0
D,sym_nodes3,uz,0
D,sym_nodes4,uz,0

csys,0
D,fixation,ux,0

lswrite,01
*GET,wall_start,ACTIVE,0,TIME,WALL		! Start time of solution
lssolve,1,1,1
FINISH
*GET,wall_end,ACTIVE,0,TIME,WALL
solving_time = wall_end - wall_start