!**************************************************************************************!
! TITLE: CFVR1252 - Shear checking according to ITER Structural Design Code for Buildings
! SUBTITLE: Check a reinforced concrete beam subjected to shear force according to ITER Structural Design Code for Buildings
!
! DESCRIPTION: A 5m length beam which is simply supported, has a T concrete reinforced section
! DESCRIPTION: and is subjected to a uniform distributed load of 84 kN/m.
! DESCRIPTION:
! DESCRIPTION: The aim of the example is to obtain the check criterion according to ITER Structural Design Code for Buildings.
!
! ELEMENT TYPE: BEAM3
! MODULES:
! UNITS: User
! KEYWORD1: Beam
! KEYWORD2: Cross Section
! KEYWORD3: Code Checking/Design
! KEYWORD4: Reinforced Concrete
! KEYWORD5: ITER
!
!**************************************************************************************!
   FINISH
  ~CFCLEAR,,1
   NomFile='CFVR1252'
  /TITLE,%NomFile%, Shear checking according to ITER Structural Design Code for Buildings

! -------------------------------------------------------------------------------------
! Initial data
! -------------------------------------------------------------------------------------
! CivilFEM SETUP
  ~UNITS,,LENG,M
  ~UNITS,,TIME,S
  ~UNITS,,FORC,KN
  ~CODESEL,EC3-92,ITER,EC2-08,,EC8-04

! Preprocessor
  /PREP7

! Parameters
  L=5          ! Length between supports, in m
  Q=84         ! Uniform distributed load in kN/m

! Material Definition
 ~CFMP,1,LIB,CONCRETE,EC2,C25/30
 ~CFMP,2,LIB,REINF,EC2,S500

! Element Type Definition
  ET,1,Beam3                           ! Element Type 1: 2D Beam

! Cross-section definition
 ~CSECDMS,1,T,1,0.375,0.15,0.75,0.15   ! T Section
! Reinforcement definition
 ~RNFDEF,1,1,2,1,1,0.02,22E-4          ! Longitudinal reinforcement
 ~RNFMDF,1,1,UFC,2
 ~RNFDEF,1,SHEAR,2,,,10.94E-4          ! Shear reinforcement
! Code properties definition
 ~SECMDF,1,EC2PROP ,REC   ,,0.05       ! 5 cm Cover
 ~SECMDF,1,EC2PROP ,KEYAST,,0
 ~SECMDF,1,EC2PROP ,RHO1  ,,0.045
! Beam properties
  ~BMSHPRO,1,BEAM,1,1,,,3,1,,,

! Model Construction
  N,1                       ! Nodes input
  N,26,L
  FILL,1,26

  E,1,2                     ! Elements input
  EGEN,25,1,1,1

  D,1, , , , , ,UX,UY       ! Constraint conditions
  D,26, , , , , ,UY

/SOLU
  SFBEAM,ALL,1,PRES,Q
  SOLVE               ! Calculation

/POST1
! Postprocessor
  ~CFSET,,1,1       ! Read first load step

! Shear checking
  ~CHKCON,SHT,SHEAR

!--------------------------------------------------------------------------------------
! 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,1) = 'VRDC'
  LABEL(2,1) = 'VRDS'
  LABEL(3,1) = 'VRDMAX'
  LABEL(4,1) = 'CRT_TOT'

! Correct values
  VALUE(1,1)=38.46
  VALUE(2,1)=139.134
  VALUE(3,1)=215.777
  VALUE(4,1)=1.51

! Comparison
 ~PLLSCON,VRDC
 *GET,VALUE(1,2),ELEM,1,ETAB,CFETAB_I
 ~PLLSCON,VRDS
  *GET,VALUE(2,2),ELEM,1,ETAB,CFETAB_I
 ~PLLSCON,VRDMAX
  *GET,VALUE(3,2),ELEM,1,ETAB,CFETAB_I
 ~PLLSCON,CRT_TOT
  *GET,VALUE(4,2),ELEM,1,ETAB,CFETAB_I

! Warning and error tolerances
  TOLER( 1, 1)= 1E-02 $ TOLER( 1, 2)= 1E-02
  TOLER( 2, 1)= 1E-02 $ TOLER( 2, 2)= 1E-02
  TOLER( 3, 1)= 1E-02 $ TOLER( 3, 2)= 1E-02
  TOLER( 4, 1)= 1E-02 $ TOLER( 4, 2)= 1E-02

!--------------------------------------------------------------------------------------
! Results comparison
!--------------------------------------------------------------------------------------
  COMPARA.MAC
