!*****************************************************************************************************!
! TITLE: CFVR1103 - EA (MV-103) check: Flexure checking.
! SUBTITLE: Flexure checking according to EA (MV-103) for a pinned-pinned beam.
!
! DESCRIPTION: A pinned-pinned steel beam is subjected to a compression force of 80 Mp and to a
! DESCRIPTION: Mz bending moment of 80*0.48 Mp*m.
! DESCRIPTION:
! DESCRIPTION: The beam has a length of 6 m.
! DESCRIPTION:
! DESCRIPTION: The cross section is a HE 240 B.
! DESCRIPTION:
! DESCRIPTION: The beam is made of A42 steel.
! DESCRIPTION:
! DESCRIPTION: The aim of the example is to obtain the flexure criterion
! DESCRIPTION: according to EA (MV-103).
!
! ELEMENT TYPE: BEAM3
! MODULES:
! UNITS: User
! KEYWORD1: Beam
! KEYWORD2: Cross Section
! KEYWORD3: Steel
! KEYWORD4: Code Checking/Design
! KEYWORD5: EA
!
!*************************************************************************************!
   FINISH
  ~CFCLEAR,,1
   NomFile='CFVR1103'
  /TITLE, %NomFile%, EA (MV-103) check: Flexure checking

! ---------------------------------------------------------------------------------
! Definition and solution of the model
! ---------------------------------------------------------------------------------
! CivilFEM Setup: Code & Units
  ~UNITS,,LENG,M
  ~UNITS,,TIME,S
  ~UNITS,,FORC,MP
  ~CODESEL,EA,EC2,,,EC8-94
/PREP7
! Preprocessor
! ---------------------------------------------------------------------------------
! Materials: A-42
  ~CFMP,1,LIB,STEEL,EA,A42
! Element Types
  ET,1,BEAM3  ! Type 1: 2D Beam
! Sections
  ~SSECLIB,2,1,6,8   !HE 240 B
! Member Properties
 ~MEMBPRO,2,EA,ALL,1,6,1,1
! Beam properties
 ~BMSHPRO,1,BEAM,2,,,,3,2,0,,

! Nodes
  N, 1
  N,11,6
  FILL, 1,11
! Elements
  REAL,1
  E,1,2
  EGEN,10,1,-1
  EPLOT

/SOLU
! Solution
! --------------------------------------------------
! Displacements
  D, 1,UX
  D, 1,UY
  D,11,UY
! Apply Load
  F,11,FX,-40.0
  F,11,MZ,-40*.48
! Solve
  SOLVE

/POST1
! Postprocessor
! --------------------------------------------------
! Read results
  ~CFSET,,1,1
! EA Bending Check
  ~CHKSTL,BENDING,

! ---------------------------------------------------------------------------------
! DATA CHECK
! ---------------------------------------------------------------------------------
! Data comparison number
  NComp = 5
  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) = 'SDYMZM '
  LABEL(2,1) = 'SCMAX  '
  LABEL(3,1) = 'SSMAX  '
  LABEL(4,1) = 'SCEQV  '
  LABEL(5,1) = 'CRTOT  '

! Correct values
! ---------------------------------------------------------------------------------
  VALUE(1,1)=  -3773.5849056604
  VALUE(2,1)= -24235.3966286
  VALUE(3,1)=    962.985254288294
  VALUE(4,1)=  24292.7246628326
  VALUE(5,1)=      0.9717

! Obtained values
! ---------------------------------------------------------------------------------
 ~PLLSSTL,SDYMZM
  *GET,VALUE(1,2),ELEM,10,ETAB,CFETAB_J
 ~PLLSSTL,SCMAX
  *GET,VALUE(2,2),ELEM,10,ETAB,CFETAB_J
 ~PLLSSTL,SSMAX
  *GET,VALUE(3,2),ELEM,10,ETAB,CFETAB_J
 ~PLLSSTL,SCEQV
  *GET,VALUE(4,2),ELEM,10,ETAB,CFETAB_J
 ~PLLSSTL,CRT_TOT
  *GET,VALUE(5,2),ELEM,10,ETAB,CFETAB_J

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

! ---------------------------------------------------------------------------------
! Results comparisson
! ---------------------------------------------------------------------------------
  COMPARA.MAC