!**************************************************************************************!
! TITLE: CFVR1427 - Concrete predesing for rectangular sections
! SUBTITLE: Predesign a reinforced concrete rectangular sections by parameters
!
! DESCRIPTION: A beam fixed in its two ends and made up by different sections
! DESCRIPTION: is subjected to a load of 48000 N/m.
! DESCRIPTION:
! DESCRIPTION: The aim of this example is to obtain the reinforcement for the different cross sections
! DESCRIPTION: for each load state.
!
! ELEMENT TYPE: BEAM4
! MODULES:
! UNITS: SI
! KEYWORD1: Cross Section
! KEYWORD2: Reinforced Concrete
!
!**************************************************************************************!
   FINISH
   ~CFCLEAR,,1
   NomFile='CFVR1427'
  /TITLE,%NomFile%,Concrete predesing for rectangular sections

! -------------------------------------------------------------------------------------
! Initial data
! -------------------------------------------------------------------------------------
! CivilFEM SETUP
~UNITS,SI
~CODESEL,EC3-92,EHE-98,EC2-91,,EC8-94

! Parameters
  WIDTH=0.4  ! Rectangular section width
  HEIGHT=0.6  ! Rectangular section height
  L1=3       ! Length of Span 1
  L2=3       ! Length of Span 2
  L3=4       ! Length of Span 3
  NE1=10     ! Number of elements in span 1
  NE2=10     ! Number of elements in span 2
  NE3=5      ! Number of elements in span 3
  P = 48000  ! Pressure over the beam

/PREP7
! Preprocessor
! -------------------------------------------------------------------
! Materials
  ~CFMP,1,LIB,CONCRETE,EHE,HA-35
  ~CFMP,2,LIB,REINF,EHE,B500S
! Sections
  ~CSECDMS,1,REC,1,HEIGHT,WIDTH
  ~CSECDMS,2,REC,1,HEIGHT+0.4,WIDTH+0.2
  ~CSECDMS,3,CIRC,1,1
! Beam properties
  ~BMSHPRO,1,BEAM,1,1,,,4,1,0,,Beam 1
  ~BMSHPRO,2,BEAM,2,2,,,4,1,0,,Beam 2
  ~BMSHPRO,3,BEAM,3,3,,,4,1,0,,Beam 3

! Element type
  ET,1,BEAM4
! Model
  REAL,1 $ TYPE,1 $ MAT,1
  N,1,0,0
  N,NE1+1,L1,0
  FILL
  E,1,2
  EGEN,NE1,1,-1
  REAL,2 $ TYPE,1 $ MAT,1
  N,NE1+NE2+1,L1+L2,0
  FILL
  E,NE1+1,NE1+2
  EGEN,NE2,1,-1
  REAL,3 $ TYPE,1 $ MAT,1
  N,NE1+NE2+NE3+1,L1+L2+L3,0
  FILL
  E,NE1+NE2+1,NE1+NE2+2
  EGEN,NE3,1,-1
! Boundary conditions
  D,1,ALL,0
  D,NE1+NE2+NE3+1,ALL,0
  ALLSEL,ALL
! Loads
  SFBEAM,ALL,2,PRES,P

/SOLU
! Solution
  SOLVE

/POST1
! Postprocessor
  ~CFSET,,1

!--------------------------------------------------------------------------------------
! 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,6
  *DIM,VALUE_CH,CHAR,Ncomp_ch,3
  *DIM,TOLER,,Ncomp,6

! Labels
!--------------------------------------------------------------------------------------
  LABEL(1,1) = 'm(3,1)'
  LABEL(2,1) = 'm(4,1)'
  LABEL(3,1) = 'm(3,11)'
  LABEL(4,1) = 'm(4,11)'

 Dirkey=0
 McTop=0.05
 McBot=0.06
 RNFMAT=2
 KFILL=1
 MATRIX='matrix'
 FiTop=20
 FiBot=20

~CLPRD, Dirkey, McTop, McBot, RNFMAT, KFILL, MATRIX, FiTop, FiBot

 VALUE(1,1)=0.0002376
 VALUE(2,1)=0.00114915

 VALUE(1,2)=matrix(3,1)
 VALUE(2,2)=matrix(4,1)

 VALUE(3,1)=0.000792
 VALUE(4,1)=0.0002376

 VALUE(3,2)=matrix(3,11)
 VALUE(4,2)=matrix(4,11)

! Warning and error tolerances
  TOLER( 1, 1)= 1E-05 $ TOLER( 1, 2)= 1E-05
  TOLER( 2, 1)= 3E-05 $ TOLER( 2, 2)= 3E-05
  TOLER( 3, 1)= 1E-05 $ TOLER( 3, 2)= 1E-05
  TOLER( 4, 1)= 1E-05 $ TOLER( 4, 2)= 1E-05

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