!*******************************************************************************
! TITLE:  CFVR0994 - CCRACKING SHELLS PURE BENDING 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='CFVR0994'
  /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 
   ! As = 15.61e-4 ;numBar = As/PI*(0.025^2)/4 ; Sep = 1/numBar
   ~SHLRNF,1,t,1,2,r,,,,,0,0.000,45.00,0.000E+00,25,0.314461,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
  
  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)= 267050000
  VALUE( 2,1)= 0.0003001
  VALUE( 3,1)= 1.003

! 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.02*VALUE(1,1) $ TOLER(1,2)= 0.02*VALUE(1,1)
  TOLER(2,1)= 0.01*VALUE(2,1) $ TOLER(2,2)= 0.01*VALUE(2,1) 
  TOLER(3,1)= 0.01*VALUE(3,1) $ TOLER(3,2)= 0.01*VALUE(3,1) 
  
!--------------------------------------------------------------------------------------
! Results comparison
!--------------------------------------------------------------------------------------
  COMPARA.MAC