fini
/clear


l=100          	! Rohrlänge [mm]
tp=20           ! Plattenhöhe
tr=10 		! Stahlringdicke
r=150 		! Äußerer Radius außerer Stahlring
r1=30         	! Äußerer Radius innerer Stahlring !


E=210000
nu=0.3
Ecm=31900         ! E-Modul Beton [N/m^2]
fcm=33.0 	  ! einaxiale Druckfestigkeit
bfcm=3*fcm        ! biaxiale Druckfestigkeit Beton [N/m^2]
fctm=2.9          ! einaxiale Zugfestigkeit Beton [N/m^2]
nuc=0.2           ! Querdehn- bzw. -kontraktionszahl Beton [-]

s=20			!Schritte
v=0.02   		!Verschiebung

z=5			!Anazhl Betonelemente in y-Richtung
y=2 			!Anzahl Stahlelemente in y-Richtung
x=18			!Anazhl Elemente in Kreisrichtung

Pi=acos(-1) 		!PI
a=(pi*r)/(2) 		!Länge der Platte


ElGross=a/10 


/prep7
! --- GEOMETRIE --------------------------------------------------
!----- Platten-----
k,1,0,r,0
k,2,a,r,0
k,3,a,r,l
k,4,0,r,l
a,1,2,3,4

!-----äußerer Stahlring-----
cylind,r-tr,r,0,l/2,  0, 90

!-----Betonring-----
cylind,r1,r-tr,0,l/2,  0, 90

!-----innerer Stahlring-----
cylind,r1-tr,r1,0,l/2,  0, 90

!vglue,all

numcmp,volu

!V1,3: Ringe   V2:Beton 


!-----Richtiges Vernetzen----
csys,1
lsel,s,loc,y,45
LESIZE,all, , ,x
lsel,all

lsel,s,loc,x,(r-tr+r)/2
lsel,a,loc,x,(r1-tr+r1)/2
LESIZE,all, , ,y
lsel,all

lsel,s,loc,x,(r1+r-tr)/2   !Beton
LESIZE,all, , ,z
lsel,all
csys,0


!------Material Stahl-----
vsel,s,,,1
vsel,a,,,3

et,2,SOLID185 			!Element Stahl
mp,ex,2,E 			!E-Modul
mp,prxy,2,nu
!esize,ElGross
type,2 
mat,2

TB,BISO,2,1,2,  
TBTEMP,0
TBDATA,,235,23.5,,,,

vmesh,all
vsel,all


!------Material Beton-----
vsel,s,,,2

et,1,solid65 		!Element Beton
mp,ex,1,Ecm 		!E-Modul
mp,prxy,1,nuc 		!Querdehnzahl

tb,concr,1		
tbdata,3,fctm
tbdata,4,fcm


TB,miso,matid,1,22	
tbtemp,22.0         
TBPT,DEFI,0,0       

TBPT,DEFI,0.0001,3.19
TBPT,DEFI,0.0002,6.15
TBPT,DEFI,0.0003,9.12
TBPT,DEFI,0.0004,11.98
TBPT,DEFI,0.0005,14.69
TBPT,DEFI,0.0006,17.23
TBPT,DEFI,0.0007,19.58
TBPT,DEFI,0.0008,21.73
TBPT,DEFI,0.0009,23.67
TBPT,DEFI,0.0010,25.4
TBPT,DEFI,0.0011,26.92
TBPT,DEFI,0.0012,28.23
TBPT,DEFI,0.0013,29.36
TBPT,DEFI,0.0014,30.30
TBPT,DEFI,0.0015,31.07
TBPT,DEFI,0.0016,31.70
TBPT,DEFI,0.0017,32.18
TBPT,DEFI,0.0018,32.54
TBPT,DEFI,0.0019,32.79
TBPT,DEFI,0.0020,32.94
TBPT,DEFI,0.0021,33.0


type,1 			!Zuweisung Element
mat,1 			!Zuweisung Material
vmesh,all 		!Vernetzen
vsel,all

!-----------Lagerung---------------------------------------------
nsel,s,loc,x,0
d,all,ux,0

nsel,s,loc,y,0
d,all,uy,0

nsel,s,loc,z,0
d,all,uz,0
allsel

!/eof

!---- Kontakt---------------------------------------------------

CM,_CONTACT,AREA
!*  
/COM, CONTACT PAIR CREATION - START 
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
/GSAV,cwz,gsav,,temp
MP,MU,2,
MAT,2   
R,3 
REAL,3  
ET,3,170
ET,4,175
KEYOPT,4,9,0
KEYOPT,4,10,2   
R,3,
RMORE,  
RMORE,,0
RMORE,0 
! Generate the target surface   
ASEL,S,,,1  
CM,_TARGET,AREA 
AATT,-1,3,3,-1  
TYPE,3  
AMESH,ALL   
! Generate the contact surface  
ASEL,S,,,2  
ASEL,A,,,3  
ASEL,A,,,4  
ASEL,A,,,5  
ASEL,A,,,6  
ASEL,A,,,7  
CM,_CONTACT,AREA
TYPE,4  
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
*SET,_REALID,3  
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,3  
ESEL,A,TYPE,,4  
ESEL,R,REAL,,3  
ASEL,S,REAL,,3  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,3  
ESEL,A,TYPE,,4  
ESEL,R,REAL,,3  
ASEL,S,REAL,,3  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END   
/VIEW,1,,,1 
/ANG,1  
/REP,FAST   
CM,_TARGET,AREA 
!*  
CM,_CONTACT,AREA
!*  
/COM, CONTACT PAIR CREATION - START 
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
/GSAV,cwz,gsav,,temp
MP,MU,2,0   
MAT,2   
R,4 
REAL,4  
ET,5,170
ET,6,174
KEYOPT,6,9,0
KEYOPT,6,10,2   
R,4,
RMORE,  
RMORE,,0
RMORE,0 
! Generate the target surface   
ASEL,S,,,2  
ASEL,A,,,3  
ASEL,A,,,4  
ASEL,A,,,5  
ASEL,A,,,6  
ASEL,A,,,7  
CM,_TARGET,AREA 
TYPE,5  
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
NSLE,A,CT2  
ESURF   
CMSEL,S,_ELEMCM 
! Generate the contact surface  
ASEL,S,,,8  
ASEL,A,,,9  
ASEL,A,,,10 
ASEL,A,,,11 
ASEL,A,,,12 
ASEL,A,,,13 
CM,_CONTACT,AREA
TYPE,6  
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
!* Create Companion Pair - Start
R,5 
REAL,5  
ET,7,170
ET,8,174
KEYOPT,8,9,0
KEYOPT,8,10,2   
R,5,
RMORE,  
RMORE,,0
RMORE,0 
TYPE,7  
ESEL,S,TYPE,,6  
NSLE,S  
ESLN,S,0
ESURF   
TYPE,8  
ESEL,S,TYPE,,5  
NSLE,S  
ESLN,S,0
ESURF   
!* Create Companion Pair - End  
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,5  
ESEL,A,TYPE,,6  
ESEL,R,REAL,,4  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,5  
ESEL,A,TYPE,,6  
ESEL,R,REAL,,4  
ESEL,A,TYPE,,7  
ESEL,A,TYPE,,8  
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END   
CM,_TARGET,AREA 
!*  
CM,_CONTACT,AREA
!*  
/COM, CONTACT PAIR CREATION - START 
CM,_NODECM,NODE 
CM,_ELEMCM,ELEM 
CM,_KPCM,KP 
CM,_LINECM,LINE 
CM,_AREACM,AREA 
CM,_VOLUCM,VOLU 
/GSAV,cwz,gsav,,temp
MP,MU,2,0   
MAT,2   
R,6 
REAL,6  
ET,9,170
ET,10,174   
KEYOPT,10,9,0   
KEYOPT,10,10,2  
R,6,
RMORE,  
RMORE,,0
RMORE,0 
! Generate the target surface   
ASEL,S,,,8  
ASEL,A,,,9  
ASEL,A,,,10 
ASEL,A,,,11 
ASEL,A,,,12 
ASEL,A,,,13 
CM,_TARGET,AREA 
TYPE,9  
NSLA,S,1
ESLN,S,0
ESLL,U  
ESEL,U,ENAME,,188,189   
NSLE,A,CT2  
ESURF   
CMSEL,S,_ELEMCM 
! Generate the contact surface  
ASEL,S,,,14 
ASEL,A,,,15 
ASEL,A,,,16 
ASEL,A,,,17 
ASEL,A,,,18 
ASEL,A,,,19 
CM,_CONTACT,AREA
TYPE,10 
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
!* Create Companion Pair - Start
R,7 
REAL,7  
ET,11,170   
ET,12,174   
KEYOPT,12,9,0   
KEYOPT,12,10,2  
R,7,
RMORE,  
RMORE,,0
RMORE,0 
TYPE,11 
ESEL,S,TYPE,,10 
NSLE,S  
ESLN,S,0
ESURF   
TYPE,12 
ESEL,S,TYPE,,9  
NSLE,S  
ESLN,S,0
ESURF   
!* Create Companion Pair - End  
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,9  
ESEL,A,TYPE,,10 
ESEL,R,REAL,,6  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,9  
ESEL,A,TYPE,,10 
ESEL,R,REAL,,6  
ESEL,A,TYPE,,11 
ESEL,A,TYPE,,12 
CMSEL,A,_NODECM 
CMDEL,_NODECM   
CMSEL,A,_ELEMCM 
CMDEL,_ELEMCM   
CMSEL,S,_KPCM   
CMDEL,_KPCM 
CMSEL,S,_LINECM 
CMDEL,_LINECM   
CMSEL,S,_AREACM 
CMDEL,_AREACM   
CMSEL,S,_VOLUCM 
CMDEL,_VOLUCM   
/GRES,cwz,gsav  
CMDEL,_TARGET   
CMDEL,_CONTACT  
/COM, CONTACT PAIR CREATION - END 




! --- BELASTUNG ---------------------------------------------------
/solu
CNVTOL,U,,0.05 
outres,all,all !gib alle resultate aus
autots,on
nsubst,5,20,3


s=15
v=0.039

*do,i,1,s
nsel,s,loc,y,r
	d,all,uy,-(i/s)*v
  allsel
  time,i
  solve
*enddo

!---------------------------------------
/EXPAND,1,RECT,HALF,0,1e-005,0,1,RECT,HALF,1e-005,0,0, ,RECT,FULL,,,,  
/post1
/dscale,,10
pldi,2

plns,s,eqv

/type,,0
/DEVICE,VECTOR,1
plcr,0,0