finish
/clear
/title, Knoten verschieben
!
!----------
!Einheiten:	
!----------
		    !Weg/Länge	mm
		    !Kraft	N
		    !Moment	Nmm
		    !E-Modul	N/mm²
		    !Druck	MPa [N/mm²]
!
!---------------
!Geometriedaten:
!---------------
a=100		!Länge in x-Richtung
b=10		!Höhe in y-Richtung
c=60		!Tiefe in z-Richtung
!
!--------------
!Materialdaten:
!--------------
e1=210000	!E-Modul
e2=210000
v1=0.3		!Querkontraktionszahl
v2=0.3
!
!*******************************************************************************
!		                               Preprocessor
!*******************************************************************************
/PREP7		!Preprocessor starten
!
!-------------------
!Materialfestlegung:
!-------------------
mp,ex,1,e1	!E-Modul Material 1
mp,nuxy,1,v1
!
mp,ex,2,e2	!E-Modul Material 2
mp,nuxy,2,v2
!
!---------------------
!Contact-Reibungszahl:
!---------------------
mp,mu,3,0.1
!
!------------------
!Elementdefinition:
!------------------
et,1,solid185
et,2,targe170,,1
et,3,conta174
!keyopt,3,10,4
!
r,1		!real constants definieren
r,2
!
!--------------------
!Geometrieerstellung:
!--------------------
!Quader 1, Gegenlauffläche (Welle):
block,0,a,0,b,0,c
!
!Quader 2, Dichtkörper:
block,a/4,(a/4)*3,b,b*8,c/6,c-c/6
!
!------------------
!Volumen vernetzen:
!------------------
type,1    !Elementtype
mat,1     !Material
real,1    !real constant
mshkey,1  !free or mapped meshing-1=mapped
esize,5   !Elementgröße
vmesh,1   !Volumen 1 vernetzen
eplot     !Produces an element display
!
!------------------
!Volumen vernetzen:
!------------------
type,1    !Elementtype
mat,2     !Material
real,1    !real constant
mshkey,1  !free or mapped meshing-1=mapped
esize,5   !Elementgröße
vmesh,2   !Volumen 1 vernetzen
eplot
!
!-------------------
!Kontakte erstellen:
!-------------------
!Nach cad.de müssen die Kontakte doppelt definiert werden, aber mit vertauschtem Contact und Target um eine Übertragung in beide Richtungen gleichmäßig zu haben. Dabei bekommt jedes Paar eine eigene Realkonstante.
!A4 ist die obere Fläche des unteren Quaders 1
!A9 ist die untere Fläche des oberen Quaders 2
!
!Quader 1 A4 Target:
asel,s,,,4
nsla,s,1
type,2      !2=Target,3=Contact
mat,3
real,2
esurf
!
!Quader 2 A9 Contact:
asel,s,,,9
nsla,s,1
type,3      !2=Target,3=Contact
esurf
allsel
!
!Quader 1 A4 Contact:
asel,s,,,4
nsla,s,1
type,3      !2=Target,3=Contact
mat,3
real,3
esurf
!
!Quader 1 A9 Target:
asel,s,,,9
nsla,s,1
type,2      !2=Target,3=Contact
esurf
allsel
!
!close gap/reduce penetration-Einstellung 
/COM, CONTACT PAIR PROPERTIES - START   
RMODIF,2,7, 
RMODIF,2,8, 
KEYOPT,3,5,3
/COM, CONTACT PAIR PROPERTIES - END
!
!*******************************************************************************
!		                               Solution
!*******************************************************************************
/solu
!------------------------
!nichtlineare Berechnung:
!------------------------
antype,static,new   !statische Analyse
eqslv,sparse        !equation solver
nlgeom,on           !nichtlineare Geometrie
pred,-1             !Quelle cad.de, weiß auch nicht, was die Angabe bringt
outres,nsol,last    !controls the solution data written to the database
nsubst,1,10,1       !nsubst,NSUBST,NSBSTP,NSBMX,NSBMN,Carry
!
!----------------
!Randbedingungen:
!----------------
asel,s,,,3
nsla,s,1
d,all,all
allsel
!
!-------------------------------------------------------------------------------
!1.Lastfall, Kraft/Druck aufbringen
!-------------------------------------------------------------------------------
time,1
NSEL,S,LOC,Y,80 !alle Knoten mit Y=80 auswählen
F,ALL,FY,-500  !Kraft FY=-1000N auf die oben ausgewählten Knoten aufbringen
allsel
solve
!-------------------------------------------------------------------------------
!2.Lastfall, Knoten verschieben
!-------------------------------------------------------------------------------
time,2
nsel,s,loc,y,80
ddele,all,all
allsel
solve
/Prep7
MODM,nocheck
asel,s,,,9
nsla,s,1
NMODIF,all,,11,
/post1
plnsol, s,eqv