!***********************************************************************!
! TITLE: CFVR3004 - EFS application on a rectangular foundation (SHELL41)
! SUBTITLE:
!
! DESCRIPTION:	This example tests the EFS application on a rectangular foundation, modelled as
! DESCRIPTION:  a rectangular shell.
! DESCRIPTION:
! DESCRIPTION:  For this purpose, it applies the maximum ballast value to the foundation.
! DESCRIPTION:
! DESCRIPTION:  It test the following EFS values:
! DESCRIPTION:
! DESCRIPTION:   - Real value on element.
! DESCRIPTION:
! DESCRIPTION:   - EFS on element 1.
! DESCRIPTION:
! DESCRIPTION:   - Minimum value.
! DESCRIPTION:
! DESCRIPTION:   - Maximum value.
! DESCRIPTION:
! DESCRIPTION:   - Average value.
!
! ELEMENT TYPE: SHELL41
! MODULES: GETC
! UNITS: SI
! KEYWORD1: Materials
! KEYWORD2: Geotechnics
! KEYWORD3: Material Models
!
!***********************************************************************!
  FINISH
  ~CFCLEAR,,1
  ~CFACTIV,GETC,Y

  NomFile='CFVR3004'
  /TITLE, %NomFile%, EFS application on a rectangular foundation (SHELL41)

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

/PREP7
! Materiales
! --------------------------------------------------
! Soil
~CFMP,1,LIB,SOIL,,PH
~CFMP,1,SOIL,NUXYST,,0.35

! Elements
! --------------------------------------------------
   ET,1,41
   R,1,1,1,1,1

! Model
! --------------------------------------------------
  ! Keypoints
   k,1,0,0
   k,2,5,0
   k,3,5,10
   k,4,0,10
  ! Lines
   l,1,2
   l,2,3
   l,3,4
   l,4,1
   lesize,all,,,5
  ! Area
   a,1,2,3,4
  ! Elements
   amap,1,1,2,3,4

! Terrain
! --------------------------------------------------
  ~TERDEF, 1,NEW, 1, 0, Z
  ~TERDEF, 1,LAYER, 1, 1, 6, 0, DB

! Calculate foundation
! --------------------------------------------------
   ~efscalc,1   ,1

! 7. APPLY EFS
! --------------------------------------------------
!  ~EFSAPPL,UCIM,VALUE
   ~efsappl,1   ,2

!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp    = 5
  NComp_ch = 0

! Marix 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)  = 'Real val'
  LABEL( 2)  = 'EFS'
  LABEL( 3)  = 'Minimum'
  LABEL( 4)  = 'Maximum'
  LABEL( 5)  = 'Average'

! Correct data, obtained data
!--------------------------------------------------------------------------------------
   VALUE( 1,1) = 235e3    ! Real value on element 1.
   VALUE( 2,1) = 297e3    ! EFS on element 1.
   VALUE( 3,1) = 151e3    ! Minimum value.
   VALUE( 4,1) = 297e3    ! Maximum value.
   VALUE( 5,1) = 188e3    ! Average value.
   ~CFGET,VALUE( 1,2),FOUNDAT,1,ELEM,EFS,,1  ! Value on element 1.
   *GET,  VALUE( 2,2),RCON,   1,CONST,6      ! EFS on element 1.
   ~CFGET,VALUE( 3,2),FOUNDAT,1,EFS,Min      ! Minimum value.
   ~CFGET,VALUE( 4,2),FOUNDAT,1,EFS,Max      ! Maximum value.
   ~CFGET,VALUE( 5,2),FOUNDAT,1,EFS,Med      ! Average value.

! Warning and error tolerances
  *DO,II,1,5
    TOLER(II,1) = 1E3
    TOLER(II,2) = 1E3
  *ENDDO

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