!*******************************************************************************
! TITLE:  CFVR0993 - CCRACKING SHELLS PURE TENSION EC_02 - Shell cracking reinforcement check according to EC2
! SUBTITLE: Check a reinforced concrete shell by cracking according to Eurocode 2 (2008)
!
! ELEMENT TYPE: Shell_181
! MODULES:
! UNITS: User
! KEYWORD1: Shell
! KEYWORD2: Shell Vertex
! KEYWORD3: Reinforced Concrete
! KEYWORD3: Cracking check
! KEYWORD4: EC2
!*******************************************************************************
   FINISH
  ~cfclear,,1
   NomFile='CFVR0993'
  /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.05        ! Cover 
  ~SHLRNF,1,t,1,2,r,,,,,0,0.000,45.00,250.000E-03,25,250.000E-03,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
  DL,4, ,UX,  
  DL,4, ,UZ,  
  DL,2, ,UZ, 
  DL,1, ,UY, 
  SFL,2,PRES,-700000,  
 
   FINISH
  /SOLU
   SOLVE
   
   FINISH
  /POST1
   ~CFSET,,1
   ~CHKCON,SHELL,CRACK, ,0.0003,0.8,0.4,3.4,0.425, , , ,0,0
   ~CFSET,1
 
!--------------------------------------------------------------------------------------
! 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) = 'SIGMA_X'
  LABEL(2) = 'WK_X'
  LABEL(3) = 'SRMAX_X'
  LABEL(4) = 'CRTTOT_X'

! Correct values
  VALUE( 1,1)= 0.17826E9
  VALUE( 2,1)= 0.358E-3  
  VALUE( 3,1)= 0.668901
  VALUE( 4,1)= 1.193

! Comparison
 ~CFSET,1
 ~CFGET,VALUE(1,2),ELEMENT,1,RESULT,SIGMA_X,I
 ~CFGET,VALUE(2,2),ELEMENT,1,RESULT,WK_X,I
 ~CFGET,VALUE(3,2),ELEMENT,1,RESULT,SRMAX_X,I
 ~CFGET,VALUE(4,2),ELEMENT,1,RESULT,CRTTOT_X,I

  ! Warning and error tolerances
  TOLER(1,1)= 0.045*VALUE(1,1) $ TOLER(1,2)= 0.045*VALUE(1,1)
  TOLER(2,1)= 0.045*VALUE(2,1) $ TOLER(2,2)= 0.045*VALUE(2,1) 
  TOLER(3,1)= 0.022*VALUE(3,1) $ TOLER(3,2)= 0.022*VALUE(3,1) 
  TOLER(4,1)= 0.035*VALUE(4,1) $ TOLER(4,2)= 0.035*VALUE(4,1) 

!--------------------------------------------------------------------------------------
! Results comparison
!--------------------------------------------------------------------------------------
  COMPARA.MAC