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 Beton [N/m^2]
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.20			!Verschiebung
z=4			!Anazhl Betonelemente in y-Richtung
y=2			!Anzahl Stahlelemente in y-Richtung
x=9			!Anazhl Elemente in Kreisrichtung
Pi=acos(-1)			!PI
a=(pi*r)/(2)		!Länge der Platte


ElGross=a/10 
! --- GEOMETRIE --------------------------------------------------


/prep7

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

k,5,-a,r,0
k,6,a,r,0
k,7,a,r,l
k,8,-a,r,l

a,1,2,3,4
a,5,6,7,8



!-----äußerer Stahlring-----

cylind,r-tr,r,0,l,  0, 90
cylind,r-tr,r,0,l, 90,180
cylind,r-tr,r,0,l,180,270
cylind,r-tr,r,0,l,270,360
vglue,1,2,3,4


!-----Betonring-----
cylind,r1,r-tr,0,l,  0, 90
cylind,r1,r-tr,0,l, 90,180
cylind,r1,r-tr,0,l,180,270
cylind,r1,r-tr,0,l,270,360
vglue,2,3,4,8

!-----innerer Stahlring-----
cylind,r1-tr,r1,0,l,  0, 90
cylind,r1-tr,r1,0,l, 90,180
cylind,r1-tr,r1,0,l,180,270
cylind,r1-tr,r1,0,l,270,360
vglue,3,4,8,12

!-----Richtiges Vernetzen----


lsel,s,,,10  
lsel,a,,,15  
lsel,a,,,22  
lsel,a,,,25  
lsel,a,,,27  
lsel,a,,,30  
lsel,a,,,45  
lsel,a,,,47  
lsel,a,,,58  
lsel,a,,,60  
lsel,a,,,61 
lsel,a,,,64  
lsel,a,,,65 
lsel,a,,,68  
lsel,a,,,89 
lsel,a,,,90
lsel,a,,,91
lsel,a,,,92
lsel,a,,,93 
lsel,a,,,94
lsel,a,,,95
lsel,a,,,96
lsel,a,,,97
lsel,a,,,98
lsel,a,,,99
lsel,a,,,100
lsel,a,,,121 
lsel,a,,,123 
lsel,a,,,126 
lsel,a,,,127
lsel,a,,,130 
lsel,a,,,131

LESIZE,all, , ,x
lsel,all

lsel,s,,,9   
lsel,a,,,11  
lsel,a,,,14  
lsel,a,,,16  
lsel,a,,,23  
lsel,a,,,26  
lsel,a,,,35  
lsel,a,,,38  
lsel,a,,,39 
lsel,a,,,44  
lsel,a,,,48  
lsel,a,,,52  
lsel,a,,,74  
lsel,a,,,78  
lsel,a,,,87  
lsel,a,,,102 
LESIZE,all, , ,y
lsel,all



lsel,s,,,21  
lsel,a,,,24  
lsel,a,,,28  
lsel,a,,,32  
lsel,a,,,42  
lsel,a,,,46  
lsel,a,,,55  
lsel,a,,,70 
LESIZE,all, , ,z
lsel,all
 

!------Material Stahl-----
vsel,s,,,1
vsel,a,,,5
vsel,a,,,6
vsel,a,,,7
vsel,a,,,3
vsel,a,,,13
vsel,a,,,14
vsel,a,,,15


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
vsel,a,,,9
vsel,a,,,10
vsel,a,,,11

!et,1,solid185
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

!esize,ElGross		!Zuweisung Elementgröße


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





! --- KONTAKT ---------------------------------------------

/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,3 
REAL,3  
ET,3,170
ET,4,174
KEYOPT,4,9,0
KEYOPT,4,10,2   
R,3,
RMORE,  
RMORE,,0
RMORE,0 
! Generate the target surface   
ASEL,S,,,2  
CM,_TARGET,AREA 
AATT,-1,3,3,-1  
TYPE,3  
AMESH,ALL   
! Generate the contact surface  
ASEL,S,,,5  
ASEL,A,,,33 
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   
/VIEW,1,1,2,3   
/ANG,1  
/REP,FAST   
/VIEW,1,1,1,1   
/ANG,1  
/REP,FAST   
/VIEW,1,1,2,3   
/ANG,1  
/REP,FAST   
!*  
/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,,,1  
CM,_TARGET,AREA 
AATT,-1,4,5,-1  
TYPE,5  
AMESH,ALL   
! Generate the contact surface  
ASEL,S,,,29 
ASEL,A,,,37 
CM,_CONTACT,AREA
TYPE,6  
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
*SET,_REALID,4  
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,5  
ESEL,A,TYPE,,6  
ESEL,R,REAL,,4  
ASEL,S,REAL,,4  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
! Reverse target normals
CM,_Y,AREA  
ASEL, , , ,       1 
CM,_YLN,LINE
CM,_YEL,ELEM
CM,_YND,NODE
LSLA,S,1
NSLA,S,1
ESLN,S,1
ESEL,R,REAL,,_REALID
ESURF,,REVERSE  
CMSEL,S,_Y  
CMSEL,S,_YLN
CMSEL,S,_YEL
CMSEL,S,_YND
CMDELE,_Y   
CMDELE,_YLN 
CMDELE,_YEL 
CMDELE,_YND 
/REPLOT 
!*  
ESEL,ALL
ESEL,S,TYPE,,5  
ESEL,A,TYPE,,6  
ESEL,R,REAL,,4  
ASEL,S,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  
ASEL,S,REAL,,4  
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,5 
REAL,5  
ET,7,170
ET,8,174
KEYOPT,8,9,0
KEYOPT,8,10,2   
R,5,
RMORE,  
RMORE,,0
RMORE,0 
! Generate the target surface   
ASEL,S,,,3  
ASEL,A,,,4  
ASEL,A,,,5  
ASEL,A,,,6  
ASEL,A,,,27 
ASEL,A,,,28 
ASEL,A,,,29 
ASEL,A,,,30 
ASEL,A,,,31 
ASEL,A,,,32 
ASEL,A,,,33 
ASEL,A,,,34 
ASEL,A,,,35 
ASEL,A,,,36 
ASEL,A,,,37 
ASEL,A,,,38 
CM,_TARGET,AREA 
TYPE,7  
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,,,9  
ASEL,A,,,10 
ASEL,A,,,11 
ASEL,A,,,12 
ASEL,A,,,47 
ASEL,A,,,48 
ASEL,A,,,49 
ASEL,A,,,50 
ASEL,A,,,51 
ASEL,A,,,52 
ASEL,A,,,53 
ASEL,A,,,54 
ASEL,A,,,55 
ASEL,A,,,56 
ASEL,A,,,57 
ASEL,A,,,58 
CM,_CONTACT,AREA
TYPE,8  
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
!* Create Companion Pair - Start
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 
TYPE,9  
ESEL,S,TYPE,,8  
NSLE,S  
ESLN,S,0
ESURF   
TYPE,10 
ESEL,S,TYPE,,7  
NSLE,S  
ESLN,S,0
ESURF   
!* Create Companion Pair - End  
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,7  
ESEL,A,TYPE,,8  
ESEL,R,REAL,,5  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,7  
ESEL,A,TYPE,,8  
ESEL,R,REAL,,5  
ESEL,A,TYPE,,9  
ESEL,A,TYPE,,10 
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,7 
REAL,7  
ET,11,170   
ET,12,174   
KEYOPT,12,9,0   
KEYOPT,12,10,2  
R,7,
RMORE,  
RMORE,,0
RMORE,0 
! Generate the target surface   
ASEL,S,,,9  
ASEL,A,,,10 
ASEL,A,,,11 
ASEL,A,,,12 
ASEL,A,,,47 
ASEL,A,,,48 
ASEL,A,,,49 
ASEL,A,,,50 
ASEL,A,,,51 
ASEL,A,,,52 
ASEL,A,,,53 
ASEL,A,,,54 
ASEL,A,,,55 
ASEL,A,,,56 
ASEL,A,,,57 
ASEL,A,,,58 
CM,_TARGET,AREA 
TYPE,11 
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,,,16 
ASEL,A,,,17 
ASEL,A,,,18 
ASEL,A,,,19 
ASEL,A,,,67 
ASEL,A,,,68 
ASEL,A,,,69 
ASEL,A,,,70 
ASEL,A,,,71 
ASEL,A,,,72 
ASEL,A,,,73 
ASEL,A,,,74 
ASEL,A,,,75 
ASEL,A,,,76 
ASEL,A,,,77 
ASEL,A,,,78 
CM,_CONTACT,AREA
TYPE,12 
NSLA,S,1
ESLN,S,0
NSLE,A,CT2 ! CZMESH patch (fsk qt-40109 8/2008) 
ESURF   
!* Create Companion Pair - Start
R,8 
REAL,8  
ET,13,170   
ET,14,174   
KEYOPT,14,9,0   
KEYOPT,14,10,2  
R,8,
RMORE,  
RMORE,,0
RMORE,0 
TYPE,13 
ESEL,S,TYPE,,12 
NSLE,S  
ESLN,S,0
ESURF   
TYPE,14 
ESEL,S,TYPE,,11 
NSLE,S  
ESLN,S,0
ESURF   
!* Create Companion Pair - End  
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,11 
ESEL,A,TYPE,,12 
ESEL,R,REAL,,7  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,11 
ESEL,A,TYPE,,12 
ESEL,R,REAL,,7  
ESEL,A,TYPE,,13 
ESEL,A,TYPE,,14 
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   
/MREP,EPLOT 

/COM, CONTACT PAIR PROPERTIES - START   
KEYOPT,4,4,1
/COM, CONTACT PAIR PROPERTIES - END 
/COM, CONTACT PAIR PROPERTIES - START   
KEYOPT,6,4,1
/COM, CONTACT PAIR PROPERTIES - END 
/MREP,EPLOT 

!-----------Lagerung---------------------------------------------
ksel,s,loc,x,0
ksel,r,loc,y,r
ksel,r,loc,z,l
dk,all,ux,0
dk,all,uz,0
ksel,all


ksel,s,loc,x,0
ksel,r,loc,y,r1
ksel,r,loc,z,l
dk,all,ux,0
dk,all,uz,0
ksel,all


ksel,s,loc,x,r
ksel,r,loc,y,0
ksel,r,loc,z,l
dk,all,uz,0
ksel,all

lsel,s,,,41
nsll,s,1
d,all,ux,0
lsel,all
nsel,all

lsel,s,,,43
nsll,s,1
d,all,ux,0
lsel,all
nsel,all

! --- BELASTUNG ---------------------------------------------------


/solu

!nlgeom,on
!solcon,on
CNVTOL,U,,0.05 		! 5% Tolerance
!CNVTOL,U,1, ,0, , 	!Konvergenzkreterium  
!neqit,50		!50 Gleichgewichtsiterationen
!eqslv,front		!lösungsmethode "frontsolver"
nsub,50,200		! erlaube bis zu 200 Sub-Iterationen
!outres,all,all		!gib alle resultate aus


*do,i,1,s
asel,s,loc,y,r
	da,all,uy,-(i/s)*v
  asel,all
  time,i
  solve
*enddo