/PREP7

!Materialdaten Femur

!Spongiosa

mptemp,1,-0.78,-0.28,0.22,0.72,1.22,1.72              !BMD-Werte von 0-2.4


et,1,solid45
mp,ex,1,1157,2059.46,803.19,-58.9,17.96     !Wirtz transversal 
mp,nuxy,1,0.12


!Kortikalis

et,2,solid45
mp,ex,2,2314,3632.98,1035.4,148.4,53.06      !Wirtz transversal
mp,nuxy,2,0.3


!Materialdaten Prothese

ET, 3, solid45
MP, nuxy,3,0.3
MP, ex, 3, 295e+3 


!Material für Rotation

ET,10,MPC184 
KEYOPT,10,1,1		  !rigid beam

!.
!.
!. Geometrieerzeugung
!.
!.



/solu

!!1.Lastschritt 

solcontrol,on      !automatische Wahl der Zeitschritte
antype,static,new  !Berechnungstyp statisch
outres,nsol,all    !alle Ergebnisse speichern fuer PostProcessing 
!nlgeom,on          !grosse Verformung einbeziehen
kbc,0              !Last linear steigern in jedem SubStep 
time,t             !Zeit fuer Lastschritt
nsubst,5,10,2      !Vorgabe SubstepAnzahl (5 aktuell, 10 maximal, 2 minimal)
neqit,100          !maximale Anzahl gleicher Iterationsschritte 
alls

nrop,unsymm
solve    


!2.Lastschritt 

kbc,0              !Last linear steigern in jedem SubStep 
esel,s,mat,,3
nsle,s,1
ddele,all,all,,,on
alls
time,2*t           !Zeit fuer Lastschritt
nsubst,5,10,2      !Vorgabe SubstepAnzahl (5 aktuell, 10 maximal, 2 minimal)
neqit,100          !maximale Anzahl gleicher Iterationsschritte 
alls

nrop,unsymm
solve              !Berechnung der Loesung



/prep7


!Beamelemente erzeugen

n,100000,-2.828,-0.7,-2.828	          !erzeugt Knoten im Prothesenzentrum
esel,s,mat,,3
nsle,s,1
cmsel,s,Momentknoten,node
nplot

type,10
mat,10
*get,nanz,node,,count      !Anzahl der Knoten

narb=0   	           !Arbeits-Knoten
*do,i,1,nanz		   !Schleife mit Nanz durchläufen  

narb=ndnext(narb)        	
e,100000,narb
*enddo

alls


/solu

antype,0,rest,2,4,continue


!Lastschritt Gangzyklus

!1.!!!!!!!!!!!!!!!!!!

kbc,0              
F,knoten2,FZ,-FZ1
F,knoten2,FX,-FX1
F,100000,Mz,M1
time,3*t
nsubst,5,10,2      
neqit,100          
alls

nrop,unsymm
solve