! Units mm, t, s .. N
! x-direction parallel railway track
! z-direction perpendicular

/prep7
youngs_concrete = 1!42300 + 4230                              ! Literature c 80/85 bspw. Schneider Tafel 5.36 or (schleuderbeton Quast 46500)           
dens_concrete = 1!2.498e-009                                  ! From Europoles plan  2.498 (consider eventually 0,5% mass loss)
 




ET,1,SOLID185                                               ! Volume elements for concrete
MP,EX,1,youngs_concrete 
MP,PRXY,1,0.2  
MP,DENS,1,dens_concrete

!-------------------------------------------------------------------------------
! Geometrical and material properties
!-------------------------------------------------------------------------------

! Meshing properties
y_1=0
y_2=100
x_1=0
x_2=100
z_1=0
z_2=100
!-------------------------------------------------------------------------------

!K, 1, 0, 0, 0
!K, 2, x_1, y_1, z_2
!K, 3, x_1, y_2, z_1

!K, 4, x_1,y_2,z_2
!K, 5, x_2,y_1,z_1
!K, 6, x_2,y_1,z_2

!K, 7, x_2,y_2,z_1
!K, 8, x_2,y_2,z_2


BLOCK,x_1,x_2,y_1,y_2,z_1,z_2,  
!A, 1,3,4,2
!A, 1,5,7,3
!A, 1,2,6,5
!A, 5,6,8,7
!A, 2,4,8,6
!A, 8,7,3,4
!allsel,all

TYPE,1
MAT,1
REAL,1
esize,20
VMES,all
NSEL,S,LOC,Y,0
D,all, , , , , ,UX,UZ,, , , 
NSEL,S,LOC,Z,0
D,all, , , , , ,UX,UY,, , , 
NSEL,S,LOC,X,0
D,all, , , , , ,UY,UZ,, , ,    

NSEL,S,LOC,X,x_2
D,all, , , , , ,UY,UZ,, , , 
NSEL,S,LOC,Y,y_2
D,all, , , , , ,UX,UZ,, , , 
NSEL,S,LOC,Z,z_2
D,all, , , , , ,UX,UY,, , , 


/SOLU
ANTYPE, STATIC
PSTRES,ON
SOLVE
FINISH
/SOL
ANTYPE,2
!*  
!*  
MODOPT,LANB,30 
EQSLV,SPAR  
MXPAND,30, , ,0 
LUMPM,0 
PSTRES,1
!*  
MODOPT,LANB,30,0,12, ,OFF   
/STATUS,SOLU
SOLVE  
/POST1  
SET,LIST  

/SOLU
ANTYPE,BUCKLE
BUCOPT,LANB,20    ! 10 ist die Anzahl der Eigenwerte in Ausgabe
EXPASS,ON
MXPAND,20
SOLVE
FINISH

!*  
!*  

/POST1
SET,LIST