Element types used in the model: BEAM3

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

   FINISH
  ~CFCLEAR,,1
  NomFile='CFVR2954'
  /TITLE,%NomFile%, SHEAR CHECKING ACCORDING TO SP 63.13330.2012

! -------------------------------------------------------------------------------------
! Initial data
! -------------------------------------------------------------------------------------
! CivilFEM SETUP
  ~UNITS,,LENG,M
  ~UNITS,,TIME,S
  ~UNITS,,FORC,KN
  ~CODESEL,,SP63133

! Preprocessor
  /PREP7

! Parameters
  L=5          ! Length between supports, in m
  Q=150        ! Uniform distributed load in kN/m

! Material Definition
~CFMP,1,LIB,CONCRETE,SP63133,B35
~CFMP,2,LIB,REINF,SP63133,A400  

! Element Type Definition
  ET,1,Beam3                               ! Element Type 1: 2D Beam

! Cross-section definition
~CSECDMS,1,REC,1,0.75,0.5
~SECMDF,1,SP63133P,REC,,0.05		! Mechanical cover for properties calculations   
~RNFDEF,1,SHEAR,2,,,0.5E-03     	! Shear Reinforcement=5cm2/m

! Beam properties
 ~BMSHPRO,1,BEAM,1,1,,,3,1,,,Beam Property

! 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     ! Boundary conditions
  D,26,UY

/SOLU
  SFBEAM,ALL,1,PRES,Q
  
  !CASE 1: 
  N=-875
  F,26,Fx,N
  SOLVE                   ! Calculation
  !CASE 2: 
  N=-3000
  F,26,Fx,N
  SOLVE                   ! Calculation
  !CASE 3: 
  N=-7350
  F,26,Fx,N
  SOLVE                   ! Calculation
  !CASE 4: 
  N=450
  F,26,Fx,N
  SOLVE                   ! Calculation

/POST1

! Postprocessor
  ~CFSET,,1,1       ! Read first load step
! Shear in Y checking for rectangular sections according to SP 63.13330.2012
  ~CHKCON,SHT,SHEAR

  ~CFSET,,2,1       ! Read first load step
! Shear in Y checking for rectangular sections according to SP 63.13330.2012
  ~CHKCON,SHT,SHEAR

  ~CFSET,,3,1       ! Read first load step
! Shear in Y checking for rectangular sections according to SP 63.13330.2012
  ~CHKCON,SHT,SHEAR
  
   ~CFSET,,4,1       ! Read first load step
! Shear in Y checking for rectangular sections according to SP 63.13330.2012
  ~CHKCON,SHT,SHEAR

!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp = 16
  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) = 'Q1_1'
  LABEL(2,1) = 'PHI_N_1'
  LABEL(3,1) = 'Q2_1'
  LABEL(4,1) = 'CRT_TOT_1'
  LABEL(5 ,1) = 'Q1_2'
  LABEL(6 ,1) = 'PHI_N_2'
  LABEL(7 ,1) = 'Q2_2'
  LABEL(8 ,1) = 'CRT_TOT_2'
  LABEL(9 ,1) = 'Q1_3'
  LABEL(10,1) = 'PHI_N_3'
  LABEL(11,1) = 'Q2_3'
  LABEL(12,1) = 'CRT_TOT_3'
  LABEL(13,1) = 'Q1_4'
  LABEL(14,1) = 'PHI_N_4'
  LABEL(15,1) = 'Q2_4'
  LABEL(16,1) = 'CRT_TOT_4'

! Correct values
  VALUE(1,1)=2304.57
  VALUE(2,1)=1.119
  VALUE(3,1)=363.553
  VALUE(4,1)=1.031
  VALUE(5,1)=2574.469
  VALUE(6 ,1)=1.25
  VALUE(7 ,1)=406.114
  VALUE(8 ,1)=0.923
  VALUE(9 ,1)=7.875
  VALUE(10,1)=3.82E-03
  VALUE(11,1)=1.241
  VALUE(12,1)=302.176
  VALUE(13,1)=1109.002
  VALUE(14,1)=0.538
  VALUE(15,1)=174.791
  VALUE(16,1)=2.145

! Comparison
*DO,I,0,3
 ~CFSET,1+I 
 ~PLLSCON,Q1
  *GET,VALUE(1+I*4,2),ELEM,1,ETAB,CFETAB_I
 ~PLLSCON,PHI_N
  *GET,VALUE(2+I*4,2),ELEM,1,ETAB,CFETAB_I
 ~PLLSCON,Q2
  *GET,VALUE(3+I*4,2),ELEM,1,ETAB,CFETAB_I
 ~PLLSCON,CRT_TOT
  *GET,VALUE(4+I*4,2),ELEM,1,ETAB,CFETAB_I
*ENDDO

! Warning and error tolerances
  TOLER( 1, 1) = VALUE( 1,  1)*0.01 $ TOLER(1  , 2)= VALUE( 1,  1)*0.01
  TOLER( 2, 1) = VALUE( 2,  1)*0.01 $ TOLER(2  , 2)= VALUE( 2,  1)*0.01
  TOLER( 3, 1) = VALUE( 3,  1)*0.01 $ TOLER(3  , 2)= VALUE( 3,  1)*0.01
  TOLER( 4, 1) = VALUE( 4,  1)*0.01 $ TOLER(4  , 2)= VALUE( 4,  1)*0.01
  TOLER( 5 , 1)= VALUE( 5 , 1)*0.01 $ TOLER(5  , 2)= VALUE( 5 , 1)*0.01
  TOLER( 6 , 1)= VALUE( 6 , 1)*0.01 $ TOLER(6  , 2)= VALUE( 6 , 1)*0.01
  TOLER( 7 , 1)= VALUE( 7 , 1)*0.01 $ TOLER(7  , 2)= VALUE( 7 , 1)*0.01
  TOLER( 8 , 1)= VALUE( 8 , 1)*0.01 $ TOLER(8  , 2)= VALUE( 8 , 1)*0.01
  TOLER( 9 , 1)= VALUE( 9 , 1)*0.03 $ TOLER(9  , 2)= VALUE( 9 , 1)*0.03
  TOLER( 10, 1)= VALUE( 10, 1)*0.03 $ TOLER(10 , 2)= VALUE( 10, 1)*0.03
  TOLER( 11, 1)= VALUE( 11, 1)*0.03 $ TOLER(11 , 2)= VALUE( 11, 1)*0.03
  TOLER( 12, 1)= VALUE(12,1)*0.03 $ TOLER( 12, 2)= VALUE(12,1)*0.03
  TOLER( 13, 1)= VALUE(13 , 1)*0.01 $ TOLER(13  , 2)= VALUE( 13 , 1)*0.01
  TOLER( 14, 1)= VALUE(14 , 1)*0.01 $ TOLER(14  , 2)= VALUE( 14 , 1)*0.01
  TOLER( 15, 1)= VALUE(15, 1)*0.01 $ TOLER(15 , 2)=VALUE( 15, 1)*0.01
  TOLER( 16, 1)= VALUE(16, 1)*0.01 $ TOLER(16 , 2)=VALUE( 16, 1)*0.01
  

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