!*************************************************************************************!
! TITLE: CFVR4055 - Reinforced concrete check according to GB 50010 2010 Chinese code
! SUBTITLE: Combined shear and torsion checking.
!
! DESCRIPTION: A concrete cantilever beam with rectangular shape section and
! DESCRIPTION: shear & torsion reinforcement is subjected to a vertical load and
! DESCRIPTION: torsion moment at free end.
! DESCRIPTION:
! DESCRIPTION: Determine the total criterion according to GB 50010-2010.
! DESCRIPTION:
!
! ELEMENT TYPE: BEAM4
! MODULES:
! UNITS: User
! KEYWORD1: Beam
! KEYWORD2: Cross Section
! KEYWORD3: Code Checking/Design
! KEYWORD4: Reinforced Concrete
! KEYWORD5: GB 50010
!
!*************************************************************************************!
FINISH
~CFCLEAR,,1

NomFile='CFVR4055'
/TITLE, %NomFile%, Combined shear and torsion checking according to GB50010 code (2010 version)

!-------------------------------------------------------------------------------
! initial data
!-------------------------------------------------------------------------------
~UNITS,,LENG,MM
~UNITS,,FORC,KN
~CODESEL,,GB50010N

/PREP7

! materials
~CFMP,1,LIB,CONCRETE,GB50010,C25
~CFMP,2,LIB,REINF,GB50010,HPB300
~CFMP,3,LIB,REINF,GB50010,HRB335

! element type
ET,1,BEAM4

! cross-section
~CSECDMS,1,REC,1,500.,250. 
~RNFDEF,1,SHEAR,2,,,0.153
~RNFDEF,1,TORSION,3,,,150.,8.,475.

! beam properties
~BMSHPRO,1,BEAM,1,1,,,4,1
~MEMBPRO,1,GB50010N,ALL,1,2,0.000  

! mesh
N,1,0
N,2,2000.
N,3,0.,0.,1.
E,1,2,3

! loads
V = 90.
T = 10000.
D,1,,,,,,ALL
F,2,FY,V
F,2,MX,T

! resolution & postprocessing
/SOLU
SOLVE
/POST1
~CFSET,,1,1
~CHKCON,SHT,BOTH

!-------------------------------------------------------------------------------
! data check
!-------------------------------------------------------------------------------

! arrays
NComp    = 13
NComp_ch = 0
*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)  = 'VRD1'
LABEL(2)  = 'CRVRD1'
LABEL(3)  = 'TRD1'
LABEL(4)  = 'CRTRD1'
LABEL(5)  = 'VRD2'
LABEL(6)  = 'CRVRD2'
LABEL(7)  = 'TRD2'
LABEL(8)  = 'CRTRD2'
LABEL(9)  = 'VRD3'
LABEL(10) = 'CRVRD3'
LABEL(11) = 'TRD3'
LABEL(12) = 'CRTRD3'
LABEL(13) = 'CRT_TOT'

! correct data
VALUE(1,1)  = 232.954
VALUE(2,1)  = 0.386
VALUE(3,1)  = 22948
VALUE(4,1)  = 0.4357
VALUE(5,1)  = 51.731
VALUE(6,1)  = 1.74
VALUE(7,1)  = 5753
VALUE(8,1)  = 1.74
VALUE(9,1)  = 70.74
VALUE(10,1) = 1.272
VALUE(11,1) = 1.52E04
VALUE(12,1) = 0.658
VALUE(13,1) = 1.272
     
! CivilFEM data
~PLLSCON,VRD1
*GET,VALUE(1,2),ELEM,1,ETAB,CFETAB_J
~PLLSCON,CRVRD1
*GET,VALUE(2,2),ELEM,1,ETAB,CFETAB_J
~PLLSCON,TRD1
*GET,VALUE(3,2),ELEM,1,ETAB,CFETAB_J
~PLLSCON,CRTRD1
*GET,VALUE(4,2),ELEM,1,ETAB,CFETAB_J
~PLLSCON,VRD2
*GET,VALUE(5,2),ELEM,1,ETAB,CFETAB_J
~PLLSCON,CRVRD2
*GET,VALUE(6,2),ELEM,1,ETAB,CFETAB_J
~PLLSCON,TRD2
*GET,VALUE(7,2),ELEM,1,ETAB,CFETAB_J
~PLLSCON,CRTRD2
*GET,VALUE(8,2),ELEM,1,ETAB,CFETAB_J
~PLLSCON,VRD3
*GET,VALUE(9,2),ELEM,1,ETAB,CFETAB_J
~PLLSCON,CRVRD3
*GET,VALUE(10,2),ELEM,1,ETAB,CFETAB_J
~PLLSCON,TRD3
*GET,VALUE(11,2),ELEM,1,ETAB,CFETAB_J
~PLLSCON,CRTRD3
*GET,VALUE(12,2),ELEM,1,ETAB,CFETAB_J
~PLLSCON,CRT_TOT
*GET,VALUE(13,2),ELEM,1,ETAB,CFETAB_J

! results comparison
EREL = 1.E-2
*DO,I,1,NCOMP
  TOLER(I,1) = EREL*ABS(VALUE(I,1))
  TOLER(I,2) = TOLER(I,1)
*ENDDO

COMPARA.MAC