!*************************************************************************************!
! TITLE: CFVR0988 - Shell shear reinforcement check according to EHE
! SUBTITLE: Shell shear reinforcement check according to the Spanish code EHE.
!
! DESCRIPTION: The model is made up of two square shell elements of 1.0m x 1.0m.
! DESCRIPTION: The nodes of one side are completely restrained and one of the nodes
! DESCRIPTION: Forces & Moments are taken from external file CFVR0988.FM
! DESCRIPTION: The reinforcement checking code is the Spanish code EHE.
! DESCRIPTION: Longitudinal Reinforcement is defined: 863.938E-06 m2/m
! DESCRIPTION: This example verifies D-02160, D-02162 (Class3 error)
!
! ELEMENT TYPE: Shell43
! MODULES:
! UNITS: User
! KEYWORD1: Shell
! KEYWORD2: Shell Vertex
! KEYWORD3: Reinforced Concrete
! KEYWORD4: EHE
!*************************************************************************************!
   FINISH
  ~cfclear,,1
   NomFile='CFVR0988'
  /TITLE, %NomFile%, Shell shear reinforcement check according to EHE
! ---------------------------------------------------------------------------------
! Model definition and solve
! ---------------------------------------------------------------------------------

! Setup
 ~UNITS,SI
 ~CODESEL,,EHE-08

  /PREP7
! Materials
~CFMP,1,LIB,CONCRETE,EHE,HA-35  
~CFMP,2,LIB,REINF,EHE,B 400 S  
 
! Element types
  ET,1,SHELL43
  
! Shell Vertex
  t=0.1    ! Thickness
  r=0.04   ! Cover

 ~SHLRNF,1,t,1,2,r,863.938E-06,863.938E-06,863.938E-06,863.938E-06,0,0.000,45.00  
 ~SHLSHR,1,2,,0,
 
! Beam & Shell property
 ~BMSHPRO,1,SHELL,1,1,1,1,43,,,,Shell property

! Nodes
  RECTNG,,1,,1,   
! Elements
  MAT,1 $ REAL,1 $ TYPE,1
  ESIZE,0.5
  
 AMESH,1
! Boundary conditions
  NSEL,S,LOC,X,0
  D,ALL,ALL
  NSEL,S,LOC,Y,1
  D,ALL, , , , , ,UX,UY,UZ, , ,  
  ALLSEL
! Solve
  /SOLU
  SOLVE

  /POST1
! ---------------------------------------------------------------------------------
! Check with CivilFEM
! ---------------------------------------------------------------------------------
 ~CFSET,,1
 
 ~FMREAD,CFVR0988,FM  ! Import the forces and moments results from an ASCII file to the RCV file.
 ~CHKCON,SHELL,SHEAR
 ~CFSET,1
 
!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp = 7
  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,1) = 'VU1'
  LABEL( 2,1) = 'CRTVU1'
  LABEL( 3,1) = 'VCU'
  LABEL( 4,1) = 'VU2'
  LABEL( 5,1) = 'CRTVU2'
  LABEL( 6,1) = 'CRT_TOT'


! Correct values
  VALUE( 1,1)=525.000E+03
  VALUE( 2,1)=0.122
  VALUE( 3,1)=116.191E+03
  VALUE( 4,1)=116.191E+03
  VALUE( 5,1)=0.552
  VALUE( 6,1)=0.552

! Comparison
 ~CFGET,VALUE(1,2), ELEMENT, 1 ,RESULT ,  VU1, I,
 ~CFGET,VALUE(2,2), ELEMENT, 1 ,RESULT ,  CRTVU1, I,
 ~CFGET,VALUE(3,2), ELEMENT, 1 ,RESULT ,  VCU, I,
 ~CFGET,VALUE(4,2), ELEMENT, 1 ,RESULT ,  VU2, I,
 ~CFGET,VALUE(5,2), ELEMENT, 1 ,RESULT ,  CRTVU2, I, 
 ~CFGET,VALUE(6,2), ELEMENT, 1 ,RESULT ,  CRT_TOT, I, 

! Warning and error tolerances
  TOLER( 1,1)= 1E0   $ TOLER( 1,2)= 1E0
  TOLER( 2,1)= 1E-03 $ TOLER( 2,2)= 1E-03
  TOLER( 3,1)= 1E1   $ TOLER( 3,2)= 1E1
  TOLER( 4,1)= 1E1   $ TOLER( 4,2)= 1E1
  TOLER( 5,1)= 1E-03 $ TOLER( 5,2)= 1E-03
  TOLER( 6,1)= 1E-03 $ TOLER( 6,2)= 1E-03


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