Element types used in the model: SHELL63

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

FINISH
~cfclear,,1
 NomFile='CFVR1274'
/TITLE,%NomFile%, Shell cracking checking according to EC2 2008
! ---------------------------------------------------------------------------------
! Model definition and solve
! ---------------------------------------------------------------------------------
/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,-2500   
		 F,3,MY,-2500  
	*ELSE
		! Y
		 D,3,ALL
		 D,4,ALL
		 F,1,FY,-500  
		 F,2,FY,-500 
		 F,1,MX,-2500   
		 F,2,MX,-2500
	*ENDIF
	ALLSEL $ SOLVE,II
	 /POST1
	 FILE,dir_%II%
	 ~CFSET,,1   
	 ~CHKCON,SHELL,CRACK, ,0.001, , , , , , , ,0,0   
	 *IF,II,EQ,1,THEN
	   ~PLSHCON,SIGMA_X,0,0,2  
	   *GET,SIGMA_X,ELEM,1,ETAB,CF_ETAB
	   ~PLSHCON,WK_X,0,0,2  
	   *GET,WK_X,ELEM,1,ETAB,CF_ETAB
	   ~PLSHCON,SRMAX_X,0,0,2  
	   *GET,SRMAX_X,ELEM,1,ETAB,CF_ETAB
	   ~PLSHCON,EM_X,0,0,2  
	   *GET,EM_X,ELEM,1,ETAB,CF_ETAB
	   ~PLSHCON,CRTTOT_X,0,0,2  
	   *GET,CRTTOT_X,ELEM,1,ETAB,CF_ETAB
	 *ELSE
	   ~PLSHCON,SIGMA_Y,0,0,2  
	   *GET,SIGMA_Y,ELEM,1,ETAB,CF_ETAB
	   ~PLSHCON,WK_Y,0,0,2  
	   *GET,WK_Y,ELEM,1,ETAB,CF_ETAB
	   ~PLSHCON,SRMAX_Y,0,0,2  
	   *GET,SRMAX_Y,ELEM,1,ETAB,CF_ETAB
	   ~PLSHCON,EM_Y,0,0,2  
	   *GET,EM_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    =  10
   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) = 'SIGMA_X'
   LABEL(2,1) = 'WK_X'
   LABEL(3,1) = 'SR_MAX_X'
   LABEL(4,1) = 'EM_X'
   LABEL(5,1) = 'CRTTOT_X'
   LABEL(6,1) = 'SIGMA_Y'
   LABEL(7,1) = 'WK_Y'
   LABEL(8,1) = 'SR_MAX_Y'
   LABEL(9,1) = 'EM_Y'
   LABEL(10,1) = 'CRTTOT_Y'
   
 ! Correct values
   VALUE(1,1)=105529484.9
   VALUE(2,1)=5.29E-05
   VALUE(3,1)=151E-3
   VALUE(4,1)=0.00035028
   VALUE(5,1)=5.29E-02
   VALUE(6,1)=105529484.9
   VALUE(7,1)=5.29E-05
   VALUE(8,1)=151E-3
   VALUE(9,1)=0.00035028
   VALUE(10,1)=5.29E-02

 ! Obtained values
   VALUE(1,2)=SIGMA_X
   VALUE(2,2)=WK_X
   VALUE(3,2)=SRMAX_X
   VALUE(4,2)=EM_X
   VALUE(5,2)=CRTTOT_X
   VALUE(6,2)=SIGMA_Y
   VALUE(7,2)=WK_Y
   VALUE(8,2)=SRMAX_Y
   VALUE(9,2)=EM_Y
   VALUE(10,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)
  TOLER(9,1)= 0.012*VALUE(9,1) $ TOLER(9,2)= 0.012*VALUE(9,1)
  TOLER(10,1)= 0.012*VALUE(10,1) $ TOLER(10,2)= 0.012*VALUE(10,1)
 !--------------------------------------------------------------------------------------
 ! Results comparison
 !--------------------------------------------------------------------------------------
   COMPARA.MAC