!Flaechenmodell der Peppe

fini
/clear
/filna,Peppe_flae2

!Modellgenerierung

/prep7

parres,,Jogi182,txt !Holt sich die Variablen aus der Datei "Jogi182.txt"

!Knotenpunkte fuer Ursprungsgeometrie

k,1,x1,y1 $k,2,x2,y2 $k,3,x3,y3 $k,4,x4,y4 $k,5,x5,y5 $k,6,x6,y6 $k,7,x7,y7 $k,8,x8,y8 $k,9,x9,y9 $k,10,x10,y10 $k,11,x11,y11 $k,12,x12,y12 $k,13,x13,y13 $k,14,x14,y14

!Knotenpunkte - Zusaetzlich, wegen Aufteilung der Flaeche fuer die Vernetzung

k,15,x10,y1 $k,16,x4+2,y4 $k,17,x4+2,y7 $k,18,x4,y7 $k,19,x14,y3 $k,20,x10,y3 $k,21,x13,y3 $k,22,x13,y1


!Flaechen

a,11,12,13,10 $a,13,21,20,10 $a,21,22,15,20 $a,13,14,19,21 $a,21,19,1,22 $a,15,2,3,20 $a,3,20,10,9 $a,3,9,8,4 $a,8,7,18,4 $a,4,18,17,16 $a,17,16,5,6

!Linien fuer die Vernetzung vorbereiten

lesize,1,,,8 $lesize,2,,,3 $lesize,3,,,8 $lesize,4,,,3 $lesize,5,,,2 $lesize,6,,,8 $lesize,7,,,2 $lesize,8,,,7 $lesize,9,,,8 $lesize,10,,,7
lesize,11,,,1 $lesize,12,,,2 $lesize,13,,,1 $lesize,14,,,7 $lesize,15,,,1 $lesize,16,,,3 $lesize,17,,,7 $lesize,18,,,3 $lesize,19,,,3 $lesize,20,,,2
lesize,21,,,10 $lesize,22,,,2 $lesize,23,,,10 $lesize,24,,,5 $lesize,25,,,2 $lesize,26,,,5 $lesize,27,,,10 $lesize,28,,,5 $lesize,29,,,10 $lesize,30,,,8
lesize,31,,,5 $lesize,32,,,8


!Materialdaten und Vernetzung

et,1,182,1,,1,,,1
TB,hyper,1,,3,MOONEY $TBDATA,,c10,c01,c11,dd               !Hier Mooney-Rivlin Gesetz
local,11,0,,,,,, $esys,11 $type,1 $mat,1 $mshape,0 $amesh,all $nrotat,all


!Berechnung

/solu
antype,static,new	!Neue statische Berechnung

!Randbedingungen
nsel,s,loc,y,0 $d,all,all $allsel $nsel,s,loc,y,y1 $*get,master,node,,num,min $cp,1,uy,all $allsel $nsel,s,loc,y,y1 $sf,all,pres,last2 $allsel

nlgeom,on	!große Verformungen an
autots,off	!Automatische Schrittvorgabe aus, um nach der Newton-Raphson Methode (bevorzugt für statische Instabilitaet) zu berechnen
nsubst,1000,10000,10	!zwischenschritte für einen Lastschritt werden spezifiziert
outres,all,10	!Loesungen die mitgeschrieben werden
arclen,on,25,0.0001	!Methode fuer Newton-Raphson, autots muss aus sein
solve