fini
/cle

nx=50		!Lokale Netzdichte in x definieren
ny=25		!Lokale Netzdichte in y definieren

e=210000
v=0.3
/prep7

!---Elementvereinbarung---
et,1,plane182   !Flächenelement für Elementtyp 1

keyopt,1,3,2 	!Plane Strain

et,2,targe169

et,3,conta171,,,2
keyopt,3,10,1

!---Materialvereinbarung----
mp,ex,1,e	!E-Modul für Material 1 zuweisen
mp,nuxy,1,v	!Querkontraktionszahl für Material 1 zuweisen 

mp,mu,2,0.0	!Zwecks Kontakt wird für Material 2 wird ein Reibungskoeffizient mit 0.0 definiert

tb,bkin,1
tbdata,1,1000,10000
!tbplot,bkin,1
!tblist,bkin,1

!---Realkonstantenvereinbarung---                   
r,2	!Realkonstantenset 1 für Elementtyp 2
r,3	!Realkonstantenset 2 für Elementtyp 3

!---Geometrievereinbarung----
rectng,0,10,0,5	!A1

k,10,10,4.5
k,11,14,4.5
k,12,14,8
k,13,10.5,8

a,10,11,12,13	!A2

!---Vernetzung---

type,1

lesize,1,,,nx	!Netzdichte für Linie 1
lesize,2,,,ny   !Netzdichte für Linie 2

mshkey,1
amesh,1

!---Kontaktvereinbarung----
!---Kontakt 1---

type,2
mat,2
real,2

lmesh,5,8

lsel,s,line,,5,8
nsll,s,1
esln,s,0
esurf
allsel

!---Kontakt 2---
type,3
mat,2
real,3

lsel,s,line,,2
nsll,s,1
esln,s,0
esurf

allsel
ncch
!---Randbedingung---
/solu
lsel,s,line,,5,8
nsll,s,1
d,all,ux,-0.5
allsel

lsel,s,line,,1
nsll,s,1
d,all,all,0
allsel

time,1
nlgeom,on
deltim,0.01,0.0001,0.05
autots,on
neqit,1000
outres,all,10
solve
!---Lösung---


/post1
set,1
pldisp,2