Element types used in the model: SHELL_181

Needed CivilFEM Modules:
Geotechnical Module
Bridges and Civil Non Linearities Module
Advanced Prestressed Concrete Module

   FINISH
  ~cfclear,,1
   NomFile='CFVR0995'
  /TITLE, %NomFile%, Shell cracking reinforcement check according to EC2

! ------------------------------------------------------------------------------
! Model definition and solve
! ------------------------------------------------------------------------------
	!Code
  ~CODESEL,,EC2-08
  	!Units
  ~UNITS,SI 
  
    /PREP7
  !Material
  ~CFMP,1,LIB,CONCRETE,EC2,C35/45 
  ~CFMP,2,LIB,REINF,EC2,S500  
  
  ! Element types
   ET,1,181
   
   ! Shell Vertex
   t=0.3        ! Thickness
   r=0.04        ! Cover 
  ~SHLRNF,1,t,1,2,r,,,,,0,0.000,45.00,0.000E+00,25,1/3,25,0.000E+00,25,0.000E+00,25
  ~SHLMDF,1,NAME,,,Shell Vertex 1   

  !Shell properties
  ~BMSHPRO,1,SHELL,1,1,1,1, 181,,, ,Shell 1   
  ~BMSHOFF, 1,0,, ,
  
  !Geometry
  RECTNG,0,1,0,1, 
  
  !Mesh
  AESIZE,1,1,  
  AMESH,1
  
  !BC
  D,1,ALL
  D,4,ALL
  F,2,MY,-50000
  F,3,MY,-50000
  F,2,FX,-50000
  F,3,FX,-50000
  
  FINISH
  /SOLU
   SOLVE
   
   FINISH
  /POST1
   ~CFSET,,1
   ~CHKCON,SHELL,CRACK, ,0.0003,0.8,0.4,3.4,0.425, , , ,0,0
   ~CFSET,1
 !~PLSHCON,SIG_X_B,0,0,2
 !~PLSHCON,WK_X_B,0,0,2
 !~PLSHCON,CRTTOT_X,0,0,2
!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp = 4
  NComp_ch = 0

! Matrix dim.
  *DIM,LABEL,CHAR,Ncomp,1
  *DIM,LABEL_CH,CHAR,Ncomp_ch,1
  *DIM,VALUE,,Ncomp,3
  *DIM,VALUE_CH,CHAR,Ncomp_ch,3
  *DIM,TOLER,,Ncomp,2

! Labels
!--------------------------------------------------------------------------------------
  LABEL(1) = 'SIG_X_B'
  LABEL(2) = 'WK_X_B'
  LABEL(3) = 'CRTTOT_X'

! Correct values
  VALUE( 1,1)= 248549379
  VALUE( 2,1)= 0.0002635
  VALUE( 3,1)= 0.8784

! Comparison
 ~CFSET,1
 ~CFGET,VALUE(1,2),ELEMENT,1,RESULT,SIG_X_B,I
 ~CFGET,VALUE(2,2),ELEMENT,1,RESULT,WK_X_B,I
 ~CFGET,VALUE(3,2),ELEMENT,1,RESULT,CRTTOT_X,I

  ! Warning and error tolerances
  TOLER(1,1)= 0.04*VALUE(1,1) $ TOLER(1,2)= 0.04*VALUE(1,1)
  TOLER(2,1)= 0.06*VALUE(2,1) $ TOLER(2,2)= 0.06*VALUE(2,1) 
  TOLER(3,1)= 0.06*VALUE(3,1) $ TOLER(3,2)= 0.06*VALUE(3,1) 
  
!--------------------------------------------------------------------------------------
! Results comparison
!--------------------------------------------------------------------------------------
  COMPARA.MAC