Element types used in the model: SHELL

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

FINISH
~cfclear,,1
 ~CODESEL,EC3-05,ACI_19,EC2-08,,EC8-04  
! ---------------------------------------------------------------------------------
! Model definition and solve
! ---------------------------------------------------------------------------------
 NomFile='CFVR9137'
 /TITLE, %NomFile%, check cracking shell according to ACI318 2019

/PREP7
ET,1,63
~CFMP,1,LIB,CONCRETE,EC2,C12/15 
~CFMP,2,LIB,REINF,EC2,S400

! Shell Vertex
  t=0.12     ! Thickness
  mc=0.021    ! Cover
  th = 0     ! Reinforcement angle

 ~SHLRNF,1,t,1,2,mc,0.0,0.0,0.0,0.0,0,th,45,0.20,12,0.20,12,0.20,12,0.20,12,
 ~BMSHPRO,1,SHELL,1,1,1,1,63,,,,Shell property
 
N,1 
N,2,1   
N,3,1,1 
N,4,,1  
TYPE,1 $ REAL,1
E,1,2,3,4 

*DO,II,1,2
	FINISH
	/FILNAME,dir_%II% ! 1->X ; 2->Y
	/SOLU
	FDELE,ALL,ALL 	$	DDELE,ALL,ALL 
	*IF,II,EQ,1,THEN
		! X
		 D,1,ALL
		 D,4,ALL
		 F,2,FX,500  
		 F,3,FX,500 
		 F,2,MY,-250   
		 F,3,MY,-250  
	*ELSE
		! Y
		 D,3,ALL
		 D,4,ALL
		 F,1,FY,-500  
		 F,2,FY,-500 
		 F,1,MX,-250   
		 F,2,MX,-250
	*ENDIF
	ALLSEL $ SOLVE,II
	 /POST1
	 FILE,dir_%II%
	 ~CFSET,,1   
	 ~CHKCON,SHELL,CRACK, ,0.5, , , , , , , ,0,0   
	 *IF,II,EQ,1,THEN
	    ~PLSHCON,S_X,0,0,2  
	   *GET,S_X,ELEM,1,ETAB,CF_ETAB
	   ~PLSHCON,FS_X,0,0,2  
	   *GET,FS_X,ELEM,1,ETAB,CF_ETAB
	   ~PLSHCON,SIGMA_X,0,0,2  
	   *GET,SIGMA_X,ELEM,1,ETAB,CF_ETAB
	   ~PLSHCON,CRTTOT_X,0,0,2  
	   *GET,CRTTOT_X,ELEM,1,ETAB,CF_ETAB
	 *ELSE
	   ~PLSHCON,S_Y,0,0,2  
	   *GET,S_Y,ELEM,1,ETAB,CF_ETAB
	   ~PLSHCON,FS_Y,0,0,2  
	   *GET,FS_Y,ELEM,1,ETAB,CF_ETAB
	   ~PLSHCON,SIGMA_Y,0,0,2  
	   *GET,SIGMA_Y,ELEM,1,ETAB,CF_ETAB
	   ~PLSHCON,CRTTOT_Y,0,0,2  
	   *GET,CRTTOT_Y,ELEM,1,ETAB,CF_ETAB
	 *ENDIF
*ENDDO

 !--------------------------------------------------------------------------------------
 ! DATA CHECK
 !--------------------------------------------------------------------------------------
 ! Data comparison number
   NComp    =  8
   NComp_ch =  0

 ! Marix 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,1) = 'S_X'
   LABEL(2,1) = 'FS_X'
   LABEL(3,1) = 'SIGMA_X'
   LABEL(4,1) = 'CRTTOT_X'
   LABEL(5,1) = 'S_Y'
   LABEL(6,1) = 'FS_Y'
   LABEL(7,1) = 'SIGMA_Y'
   LABEL(8,1) = 'CRTTOT_Y'
   
 ! Correct values
   VALUE(1,1)=7.651
   VALUE(2,1)=10.987e6
   VALUE(3,1)=1.01e6
   VALUE(4,1)=65.349e-3
   VALUE(5,1)=7.651
   VALUE(6,1)=10.987e6
   VALUE(7,1)=1.01e6
   VALUE(8,1)=65.349e-3

 ! Obtained values
   VALUE(1,2)=S_X
   VALUE(2,2)=FS_X
   VALUE(3,2)=SIGMA_X
   VALUE(4,2)=CRTTOT_X
   VALUE(5,2)=S_Y
   VALUE(6,2)=FS_Y
   VALUE(7,2)=SIGMA_Y
   VALUE(8,2)=CRTTOT_Y
   
 ! Warning and error tolerances
  TOLER(1,1)= 0.012*VALUE(1,1) $ TOLER(1,2)= 0.012*VALUE(1,1)
  TOLER(2,1)= 0.012*VALUE(2,1) $ TOLER(2,2)= 0.012*VALUE(2,1) 
  TOLER(3,1)= 0.012*VALUE(3,1) $ TOLER(3,2)= 0.012*VALUE(3,1) 
  TOLER(4,1)= 0.012*VALUE(4,1) $ TOLER(4,2)= 0.012*VALUE(4,1) 
  TOLER(5,1)= 0.012*VALUE(5,1) $ TOLER(5,2)= 0.012*VALUE(5,1) 
  TOLER(6,1)= 0.012*VALUE(6,1) $ TOLER(6,2)= 0.012*VALUE(6,1)
  TOLER(7,1)= 0.012*VALUE(7,1) $ TOLER(7,2)= 0.012*VALUE(7,1)
  TOLER(8,1)= 0.012*VALUE(8,1) $ TOLER(8,2)= 0.012*VALUE(8,1)

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