!***************************************************************************************************!
! TITLE: CFVR2816 -  Shear design according to IS456
! SUBTITLE: Design a reinforced concrete beam subjected to Shear force according to IS456
!
! DESCRIPTION: A beam which is fixed at the right and simply supported at the left,
! DESCRIPTION: has a rectangular reinforced concrete section
! DESCRIPTION: and it is subjected to a uniform distributed load of 7.5 Mp/m.
! DESCRIPTION:
! DESCRIPTION: The aim of the example is to obtain the necessary steel amount according to IS456.
!
! ELEMENT TYPE: BEAM3
! MODULES:
! UNITS: User
! KEYWORD1: Beam
! KEYWORD2: Cross Section
! KEYWORD3: Code Checking/Design
! KEYWORD4: Reinforced Concrete
! KEYWORD5: IS 456
!
!***************************************************************************************************!
  FINISH
  ~CFCLEAR,,1
  NomFile='CFVR2816'
 /TITLE,%NomFile%, SHEAR DESIGN ACCORDING TO IS456

! -------------------------------------------------------------------------------------
! Initial data
! -------------------------------------------------------------------------------------
! CivilFEM SETUP
  ~UNITS,,LENG,M
  ~UNITS,,TIME,S
  ~UNITS,,FORC,MP
  ~CODESEL,EC3-92,IS456,,,EC8-94

! Preprocessor
  /PREP7

! Parameters
  L = 10     ! Length between supports, in m
  b = 0.3    ! Cross-section width in m.
  h = 0.5    ! Depth in m.
  Q = 7.5    ! Uniform distributed load in Mp/m

! Material Definition
  ~CFMP,1,LIB,CONCRETE,IS456,M35
  ~CFMP,2,LIB,REINF,IS1786,Fe500

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

! Cross-section definition
  ~CSECDMS,1,REC,1,h,b

! Reinforcement definition
  ~RNFDEF,1,SHEAR,2,0.0,0.0,1E-4,0.0   ! Shear reinf.

! Code properties definition
  ~SECMDF,1,IS456P,REC     ,,5e-002

! Beam & Shell property
  ~BMSHPRO,1,BEAM,1,1,,,3,1,0,,Beam property

! Model
  N,1
  N,21,L
  FILL,1,21
  E,1,2
  EGEN,20,1,1,1

! Boundary conditions
  D, 1,UX,0
  D, 1,UY,0
  D,21,UY,0

/SOLU
! Solution
  SFBEAM,ALL,1,PRES,Q
  SOLVE               ! Calculation

/POST1
! Postprocessor
  ~CFSET,,1           ! Read first load step

! Shear in Y design
  ~DIMCON,SHT,SHEAR

!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparision number
  NComp = 3
  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) = 'VUS'
  LABEL(2,1) = 'VUC'
  LABEL(3,1) = 'ASSH'

! Correct Values
  VALUE(1,1) = 29.758
  VALUE(2,1) =  3.992
  VALUE(3,1) =  0.180E-2

! Comparison
 ~PLLSCON,VUS
 *GET,VALUE(1,2),ELEM,19,ETAB,CFETAB_J
 ~PLLSCON,VUC
 *GET,VALUE(2,2),ELEM,19,ETAB,CFETAB_J
 ~PLLSCON,ASSH
 *GET,VALUE(3,2),ELEM,19,ETAB,CFETAB_J

! 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-05 $ TOLER( 3, 2)= 1E-05
!--------------------------------------------------------------------------------------
! Results Comparison
!--------------------------------------------------------------------------------------
  COMPARA.MAC