! Schalenelemente (Stütze)
	et,1,181
	keyopt,1,3,2
	keyopt,1,8,2
!
	sectype,1,shell                            				!Section definieren
	secdata,S,1,,3
	
! Schalenelemente (Palettenträger)
	et,2,181
!	keyopt,2,1,1
	keyopt,2,3,2
	keyopt,2,8,2
!
	sectype,2,shell                            				!Section definieren
	secdata,S1,1,,3
	
! Volumenelemente (Schrauben)
	et,10,185
	
! Kontaktelemente
	et,11,targe170
! Bonded Kontakt
	et,12,conta173
   	keyopt,12,5,3  
   	keyopt,12,11,1   
   	keyopt,12,12,3
  	et,13,conta175
  	keyopt,13,5,3
  	keyopt,13,11,1
  	keyopt,13,12,3
! Abhebender Kontakt
	et,15,conta173
   	keyopt,12,5,3  
   	keyopt,12,11,1   
   	keyopt,12,12,1

! Kontaktdefinition (für eine Schraube!)
!--------------------------------------------------------------------
!Kontaktpaare Schraubenschaft / Stütze
!--------------------------------------------------------------------
/prep7
mat,3   !Reibungszahl 0.2
!Generierung der Contactelemente auf der Stütze
    csys,100+i*100
    numstr,elem,60000+i*n_help
    type,13     	
    real,10+i*10
     r,10+i*10   
    cmsel,s,Stuetze,area
    lsla,s
    lsel,r,radius,,(d+delta_d)/2
    lsel,r,loc,x,-(d+delta_d)/2,(d+delta_d)/2
    lsel,r,loc,y,-(d+delta_d)/2,(d+delta_d)/2
    nsll,s,1
    esurf
    alls

!Generierung der Targetelemente auf dem Schraubenschaft  
!  csys,122+i*100
    type,11             
    real,10+i*10
    
    vsel,s,loc,z,-(kh+t),-(kh+t+20)
    vsel,r,loc,y,-d1/2,d1/2
    vsel,r,loc,x,-d1/2,d1/2

    aslv,s
    asel,u,loc,z,-(kh+t)
    asel,u,loc,z,-(kh+t+20)
    asel,u,loc,x,0
    asel,u,loc,y,0
    nsla,s,1
 
    esln,s
      esurf
    alls

!--------------------------------------------------------------------
!Kontaktpaare Schraubenschaft / Palettenträger    
!--------------------------------------------------------------------
numstr,elem,70000
!Generierung der Targetelemente auf dem Palettenträger
    type,13
    real,11+i
      r,11+i      
    cmsel,s,Palettenträger,area
    lsla,s
    lsel,r,radius,,(d+delta_d)/2
    lsel,r,loc,x,-(d+delta_d)/2,(d+delta_d)/2
    lsel,r,loc,y,-(d+delta_d)/2,(d+delta_d)/2
    nsll,s,1
      esurf
    alls 
!Generierung der Contactelemente auf dem Schraubenschaft
    type,11             
    real,11+i
    
    vsel,s,loc,z,-(kh+t),-(kh+t+20)
    vsel,r,loc,y,-d1/2,d1/2
    vsel,r,loc,x,-d1/2,d1/2
    aslv,s
    asel,u,loc,z,-(kh+t)
    asel,u,loc,z,-(kh+t+20)
    asel,u,loc,x,0
    asel,u,loc,y,0
    nsla,s,1  
    esln,s
      esurf
    alls

!--------------------------------------------------------------------
! Kontaktpaare Schraubenkopf / Palettenträger
!--------------------------------------------------------------------
mat,3   
numstr,elem,75000
! Generierung der Targetelemente auf dem Palettenträger
    type,11      	!Targetelement
    real,50+i
      r,50+i       !REAL-Werte auf default belassen!!!
    cmsel,s,Palettenträger,area
    nsla,s,1
    nsel,r,loc,z,-(kh+t+S1/2)
    nsel,r,loc,x,-w/2,w/2
    nsel,r,loc,y,-e/2,e/2
    esln,s
      esurf
      esurf,,reverse
    alls
! Generierung der Contactelemente auf dem Schraubenkopf   
    type,12             !Contactelement
    real,50+i
    vsel,s,loc,z,0,-(kh+t)
    vsel,r,loc,y,-e/2,e/2
    vsel,r,loc,x,-w/2,w/2
    aslv,s
    asel,r,loc,z,-(kh+t)
    nsla,s,1
    esln,s
    esurf
    alls

!--------------------------------------------------------------------
! Kontaktpaare Mutter / Stütze
!--------------------------------------------------------------------

mat,3  
! Generierung der Targetelemente auf der Stütze
numstr,elem,80000
  type,11      	!Targetelement
  real,60+i
    r,60+i      !REAL-Werte auf default belassen!!!
  cmsel,s,Stuetze,area
  nsla,s,1
  nsel,r,loc,z,-(kh+t+S1+S/2)
  nsel,r,loc,x,-sM/2,sM/2
  nsel,r,loc,y,-eM/2,eM/2
  esln,s
   esurf
   esurf,,reverse
  alls
! Generierung der Contactelemente auf der Mutter
  type,12             !Contactelement
  real,60+i
  vsel,s,loc,z,-(kh+t+S1+S+mh/2)
  vsel,r,loc,y,-eM/2,eM/2
  vsel,r,loc,x,-sM/2,sM/2
  aslv,s
  asel,r,loc,z,-(kh+t+S1+S)
  nsla,s,1
  esln,s
  esurf
  alls	

!--------------------------------------------------------------------
! Kontaktpaare Stütze / Palettenträger im Bereich der Schrauben
!--------------------------------------------------------------------	

mat,3   
! Symmetrischer Kontakt
! Generierung der Targetelemente auf der Stütze
  numstr,elem,85000
  type,11      	!Targetelement
  real,70+i
    r,70+i       !REAL-Werte auf default belassen!!!
  cmsel,s,Stuetze,area
  nsla,s,1
  nsel,r,loc,z,-(kh+t+S1+S/2)
  nsel,r,loc,x,-fak*(d+delta_d),fak*(d+delta_d)
  nsel,r,loc,y,-fak*(d+delta_d),fak*(d+delta_d)
  esln,s
   esurf
   esurf,,reverse
  alls
	
! Generierung der Contactelemente auf dem Palettenträger
  type,15             !Contactelement
  real,70+i
  cmsel,s,Palettenträger,area
  nsla,s,1
  nsel,r,loc,z,-(kh+t+S1/2)
  nsel,r,loc,x,-fak*(d+delta_d),fak*(d+delta_d)
  nsel,r,loc,y,-fak*(d+delta_d),fak*(d+delta_d)
  esln,s
  esurf
  alls

!--------------------------------------------------------------------
! Kontaktpaare Mutter / Schraubenschaft
!--------------------------------------------------------------------

mat,3 
! Generierung der Targetelemente auf der Mutter
numstr,elem,80000
  type,11      	!Targetelement
  real,80+i
    r,80+i     
  vsel,s,loc,z,-(kh+t+S1+S+mh/2)
  vsel,r,loc,y,-eM/2,eM/2
  vsel,r,loc,x,-sM/2,sM/2
  aslv,s
  asel,r,loc,x,-d1/2,d1/2
  asel,r,loc,y,-d1/2,d1/2
  nsla,s,1
  esln,s
   esurf
  alls
! Generierung der Contactelemente auf dem Schaft
  type,12             !Contactelement
  real,80+i
  vsel,s,loc,z,-(kh+t),-(kh+t+20)
  vsel,r,loc,y,-d1/2,d1/2
  vsel,r,loc,x,-d1/2,d1/2
  aslv,s
  asel,u,loc,z,-(kh+t)
  asel,u,loc,z,-(kh+t+20)
  asel,u,loc,x,0
  asel,u,loc,y,0
  nsla,s,1
  esln,s
  esurf
  alls


! Randbedingungen
! Fixierung aller Knoten oben und unten an der Stütze
  nsel,s,loc,y,L
  nsel,a,loc,y,0
  d,all,uy
  d,all,ux
  d,all,uz
  alls
	
! Schrauben in y- und x-Richtung halten
	i=1
	csys,100+i*100
	vsel,all
	nslv,s
	nsel,s,loc,z,0
	nsel,a,loc,z,-(kh+t+20)
	d,all,uy,0
	d,all,ux,0
	alls
	csys,0		
	
! Aufbringen der Pretension Elemente Schraube 1
	i=1
	csys,100+i*100
	vsel,s,loc,z,-(kh+t),-(kh+t+20)
  	vsel,r,loc,y,-d1/2,d1/2
  	vsel,r,loc,x,-d1/2,d1/2
	eslv,s
	psmesh,10,,,all,,0,z
	sload,10,PL01,lock,forc,50,1,2
	alls
	csys,0
	
! Aufbringen der Pretension Elemente Schraube 2	
	i=2
	csys,100+i*100
	vsel,s,loc,z,-(kh+t),-(kh+t+20)
 	vsel,r,loc,y,-d1/2,d1/2
  	vsel,r,loc,x,-d1/2,d1/2
	eslv,s
	psmesh,11,,,all,,0,z
	sload,11,PL01,lock,forc,50,1,2
	alls
	csys,0	
	
! Aufbringen der Pretension Elemente Schraube 3	
	i=1
	csys,500+i*100
	vsel,s,loc,z,-(kh+t),-(kh+t+20)
  	vsel,r,loc,y,-d1/2,d1/2
  	vsel,r,loc,x,-d1/2,d1/2
	eslv,s
	psmesh,12,,,all,,0,z
	sload,12,PL01,lock,forc,50,1,2
	alls
	csys,0

! Aufbringen der Pretension Elemente Schraube 4	
	i=2
	csys,500+i*100
	vsel,s,loc,z,-(kh+t),-(kh+t+20)
  	vsel,r,loc,y,-d1/2,d1/2
  	vsel,r,loc,x,-d1/2,d1/2
	eslv,s
	psmesh,13,,,all,,0,z
	sload,13,PL01,lock,forc,50,1,2
	alls
	csys,0		
	
finish
/solu
! Nichtlineare Berechnung
! 1. Load Step Vorspannung
	pred,off
	time,1 
	antype,static
	autots,on
	nlgeom,on
	nsubst,100,100,0
	eqslv,sparse
	pstres,on				
	kbc,0
	nropt,FULL

	outpr,all,all
	outres,all,all

	solve
	finish
		
! 2. Load Step Verschiebung
/solu
	pred,off
	time,2
	antype,static
	autots,on
	nlgeom,on
	nsubst,INCRE,INCRE,0
	eqslv,sparse
	pstres,on									
	kbc,0
	nropt,FULL

	outpr,all,all
	outres,all,all
	
! Lagerung Schrauben löschen
	i=1
	csys,100+i*100
	vsel,all
	nslv,s
	nsel,r,loc,z,0
	nsel,a,loc,z,-(kh+t+20)
	ddele,all,all
	alls
	csys,0

! Zwangsverschiebung aufbringen	
	csys,13
	nsel,s,loc,x,-(20+447.5-5),-(20+447.5+5)
	d,all,uy,uy
	alls
	csys,0

	solve
	finish