!/output,OUTPUT,lis
finish			! von Präprozessor-Level auf Begin-Level wechseln
/clear			! eventuell vorhandene Datenbasis löschen
/graphics,full		! wichtig, sonst nur Spannungen-Resultate für Oberflächenknoten
/NERR,0,99999999	! Programmabbruch von ANSYS nach 99999999 Warnungen/Fehlermeldungen
/UIS,msgpop,3		! nur Fehlermeldungen anzeigen, keine Warnungen
/cmap
/PNUM,LINE,1  		! Einschalten der Liniennummerierung und Farbkodierung im Plot
/PNUM,KPOI,1  		! Einschalten der Keypointnummerierung und Farbkodierung im Plot
/PNUM,area,1  		! Einschalten der Flächennummerierung und Farbkodierung im Plot
/PNUM,volu,1  		! Einschalten der Volumennummerierung und Farbkodierung im Plot
/PLOPTS,INFO,auto	! Auto-legend-Modus, es folgen Legendeneinstellungen
/PLOPTS,LEG1,on
/PLOPTS,LEG2,on
/PLOPTS,LEG3,on
/PLOPTS,FRAME,on 	! Rahmen um Darstellungsfenster
/PLOPTS,TITLE,on 	! Titel an
/PLOPTS,MINM,on  	! Minima- und Maximasymbole ein
/PLOPTS,FILE,on  	! Dateiname in Legende wiedergeben
/PLOPTS,WINS,on  	! Darstellung automatisch skalieren und an Rahmen anpassen
/PLOPTS,WP,0		! Darstellung der Arbeitsebenen aus
/PLOPTS,DATE,2  	! Zeige Datum und Zeit in Legende
/VSCALE,1,1,0		! Vektoren, z. B. Koordinatenkreuz (triad), skalieren
/TRIAD,ORIG 		! Globales Koordinatenkreuz im Ursprung anzeigen
/pbc,forc,,2		! Numerischen Wert der Kraft in Darstellung wiedergeben
/REPLOT
/units,si		! SI-Einheiten (1N/mm^2=1e6N/m^2)

*abbr,Platte-Zylinder,/input,Platte-Zylinder.inp		!Abkürzung anlegen
/title,Kontakt Zylinder-Ebene


!Variablendeklarierung
!=====================

PI=2*ASIN(1)


!Elementgeometrien

Rz=30			![mm]	Radius Zylinder
Bz=15			![mm]	Breite Zylinder

Bp=20			![mm]	Plattenbreite
Lp=40			![mm]	Plattenlänge
Dp=30			![mm]	Plattendicke

r=50			![mm]	Radius Balligkeit !=0 wenn keine Balligkeit 

!Werkstoffdaten


Emodulz=2.1e5		![N/mm²] E-Modul Zylinder
Nuehz=0.3		!	 Querkontraktion Zylinder
Emodulp=2.1e5		![N/mm²] E-Modul Platte
Nuehp=0.3		!	 Querkontraktion Platte

f=0.1			!	Reibzahl im Kontakt


!Vernetzungsdaten

Netzdichtez=3		!	Netzdichte Zylinder Breitenrichtung
Netzdichtezb=5		!	Netzdichte Zylinder Breitenrichtung

Netzdichtep=5		!	Netzdichte Platte

!Randbedingungen

F=1000			![N] 	Kraft auf Zylinder in y-Richtung


!Variablenberechnung

*if,r,gt,0
	l=r-sqrt((r**2)-((Bz**2)/4))
*endif

*if,r,eq,0,then
	l=0
*endif

Verschiebung=Rz/100

!Präprozessor
!============

/prep7

!Elementtyp und Materialzuweisung


ET,1,SOLID65		!Elementtyp für Zylinder
ET,2,SOLID65		!Elementtyp für Platte


!Geometrie

Bz=Bz/2			!Halbierung, da symmetrie ausgenutzt wird
Bp=Bp/2

mp,ex,1,Emodulz 	!E-Modul Zylinder  
mp,prxy,1,Nuehz		!Querkontraktion Zylinder

mp,ex,2,Emodulz		!E-Modul Platte   
mp,prxy,2,Nuehz		!Querkontraktion Platte

local,100,0,0,E_Z1,0			! kartesisches Koordinatensystem
local,101,1,0,E_Z1,0			! Zylinderkoordinaten
local,102,0,0,-Interferenz		! kartesisches Koordinatensystem

csys,1

k,3,0.9*Rz,-52
k,4,0.9*Rz,-80
k,12,1.1*Rz,-54
k,13,1.1*Rz,-82

csys,0

k,1,0,0
k,2,0.9*Rz,0
k,5,0,-0.9*Rz
k,6,0.5*Rz,0
k,7,0.3*Rz,-0.3*Rz
k,8,0.18*Rz,-0.5*Rz
k,9,0,-0.5*Rz
k,10,0,-0.3*Rz
k,11,1.1*Rz,0
k,14,0,-1.1*Rz


Larc,2,3,7,0.9*Rz	!Kreisbögen innen	!L1
Larc,3,4,8,0.9*Rz				!L2
Larc,4,5,8,0.9*Rz				!L3
Larc,11,12,7,Rz		!Kreisbögen außen	!L4
Larc,12,13,8,Rz					!L5
Larc,13,14,9,Rz					!L6


L,6,7			!Verbindungslinien	!L7
L,7,8						!L8
L,8,9						!L9
L,5,9						!L10
L,9,10						!L11
L,10,1						!L12
L,1,6						!L13
L,6,2						!L14
L,3,7						!L15
L,4,8						!L16
L,10,7						!L17
L,2,11						!L18
L,3,12						!L19
L,4,13						!L20
L,5,14						!L21


al,1,15,7,14		!Definition der einzelnen Flächen
al,2,16,8,15
al,3,10,9,16
al,8,9,11,17
al,12,17,13,7
al,18,4,19,1
al,5,20,2,19
al,20,6,21,3


k,15,0,0,Bz		!Keypoints für rotierende Fläche
k,16,Rz,0
k,17,Rz-l,0,Bz

*if,r,eq,0,then

	L,16,17		!Verbindungslinie
*else

	Larc,16,17,15,r		!Kreisbogen
*endif

Arotat,22,,,,,,1,15,-90,1	!Rotation der Linie

*Do,i,1,8
	voffst,i,-Bz
*enddo	

aplot
	
!Sicht ändern

/VIEW,1,1,2,3   
/ANG,1  
/REP,FAST

FLST,2,3,6,ORDE,2   	!Trennen der Volumina durch Fläche
FITEM,2,6   
FITEM,2,-8  
VSBA,P51X,9 

FLST,2,3,6,ORDE,3   	!Löschen der äußeren Volumina
FITEM,2,9   
FITEM,2,11  
FITEM,2,13  
VDELE,P51X,,,1


allsel,all,all		!Verschmelzen der Linien Volumina und Flächen
nummrg,kp
numcmp,volu
numcmp,area
numcmp,line
numcmp,kp
allsel,all,all

lsel,s,length,,Bz
lesize,all,,,Netzdichtezb !Alle Linien in Breitenrichtung bekommen Netzdichte
alls

Nd=Netzdichtez

lesize,1,,,Nd
lesize,4,,,Nd
lesize,9,,,Nd
lesize,41,,,Nd

Nd=1.5*Netzdichtez

lesize,7,,,Nd
lesize,31,,,Nd

lesize,11,,,Nd
lesize,20,,,Nd

lesize,12,,,Nd
lesize,18,,,Nd

lesize,13,,,Nd
lesize,26,,,Nd

Nd=Netzdichtez

lesize,2,,,Nd
lesize,25,,,Nd

lesize,5,,,Nd
lesize,27,,,Nd

lesize,8,,,Nd
lesize,35,,,Nd

lesize,48,,,Nd
lesize,49,,,Nd



Nd=2*Netzdichtez

lesize,3,,,Nd
lesize,30,,,Nd

lesize,6,,,Nd
lesize,32,,,Nd

lesize,10,,,Nd
lesize,39,,,Nd

lesize,14,,,Nd
lesize,36,,,Nd

lesize,53,,,Nd,5
lesize,54,,,Nd,5


Nd=Netzdichtez

lesize,44,,,Nd
lesize,46,,,Nd

lesize,45,,,Nd
lesize,47,,,Nd

lesize,51,,,Nd
lesize,52,,,Nd

lesize,55,,,Nd,0.2
lesize,56,,,Nd,5



block,,Lp,-Rz,-Rz-Dp,,Bp



Nd=2*Netzdichtep

lesize,58,,,Nd,0.1
lesize,63,,,Nd,0.1
lesize,60,,,Nd,10
lesize,61,,,Nd,10

Nd=Netzdichtep

lesize,57,,,Nd,5
lesize,64,,,Nd,5
lesize,62,,,Nd,0.2
lesize,59,,,Nd,0.2

Nd=Netzdichtep


lesize,65,,,Nd
lesize,66,,,Nd
lesize,67,,,Nd
lesize,68,,,Nd
alls

!Vernetzung

MSHAPE,0,3D 
MSHKEY,1
!*  
FLST,5,9,6,ORDE,2   
FITEM,5,1   
FITEM,5,-9  
CM,_Y,VOLU  
VSEL, , , ,P51X 
CM,_Y1,VOLU 
CHKMSH,'VOLU'   
CMSEL,S,_Y  
!*  
VMESH,_Y1   
!*  
CMDELE,_Y   
CMDELE,_Y1  
CMDELE,_Y2  


!Kontaktbedingungen

CM,_TARGET,AREA 
!*  
CM,_CONTACT,AREA
!*  
/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,1000
MAT,1   
MP,EMIS,1,7.88860905221e-031
R,3 
REAL,3  
ET,3,170
ET,4,174
R,3,,,1.0,0.1,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
KEYOPT,4,4,0
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,0   
KEYOPT,4,2,0
KEYOPT,3,5,0

Keyopt,cid,4,3			!Mortar Algorithmus


! Generate the target surface   
ASEL,S,,,38 
ASEL,A,,,39 
ASEL,A,,,40 
ASEL,A,,,41 
ASEL,A,,,42 
ASEL,A,,,43 
CM,_TARGET,AREA 
TYPE,3  
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
NSLE,A,CT2  
ESURF   
CMSEL,S,_ELEMCM 
! Generate the contact surface  
ASEL,S,,,7  
ASEL,A,,,12 
ASEL,A,,,16 
ASEL,A,,,25 
ASEL,A,,,26 
ASEL,A,,,27 
ASEL,A,,,28 
ASEL,A,,,30 
ASEL,A,,,31 
ASEL,A,,,32 
ASEL,A,,,34 
ASEL,A,,,35 
ASEL,A,,,36 
ASEL,A,,,37 
CM,_CONTACT,AREA
TYPE,4  
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
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   
/MREP,EPLOT 


!Loesungsprozessor
!=================

/solu


!Randbedingungen

asel,s,loc,z,0		
nsla,s,1
d,all,uz
alls

asel,s,loc,x,0
nsla,s,1
d,all,ux
alls

asel,s,loc,y,-Rz-Dp
nsla,s,1
d,all,uy
alls

asel,s,loc,x,lp
nsla,s,1
d,all,ux
alls

NSEL,S,LOC,Y,0				
CP,1,UY,ALL				
*GET,NC,NODE,,NUM,MIN			
NSEL,ALL
allsel,all,all

D,NC,UY,-verschiebung

!asel,s,loc,y,0
!nsla,s,1
!d,all,uy,-Verschiebung
!alls

auto,on
nsubst,20,50
ncnv,0
neqit,100
nlgeom,on
pred,on
allsel,all,all
SOLVE		



nsel,s,loc,y,0
DDELE,all,UY		! DELETE IMPOSED DISPLACEMENT
F,all,FY,-F		! APPLY LOAD ON MODEL
alls
nsubst,8,8,8
AUTOTS,ON
NROPT,FULL,,OFF
!	ncnv,0
!	neqit,100
LNSRCH,auto
!	pred,on,,on
!	nlgeom,on
allsel,all,all
SOLVE                  	! SOLVE SECOND LOAD STEP