!***********************************************************************!
! TITLE: CFVR3002 - Rectangular Foundation Stiffness (SHELL99)
! SUBTITLE:
!
! DESCRIPTION:	This example tests the stiffness of a rectangular foudation, modelled as
! DESCRIPTION:  a rectangular shell.
! DESCRIPTION:
! DESCRIPTION:  For this purpose, it tests the stiffness at the following
! DESCRIPTION:  locations:
! DESCRIPTION:
! DESCRIPTION: 	- Corners
! DESCRIPTION:
! DESCRIPTION:  - Sides
! DESCRIPTION:
! DESCRIPTION:  - Center
! DESCRIPTION:
! DESCRIPTION:  - Maximum global value
!
! ELEMENT TYPE: SHELL99
! MODULES: GETC
! UNITS: SI
! KEYWORD1: Materials
! KEYWORD2: Geotechnics
! KEYWORD3: Material Models
!
!***********************************************************************!
  FINISH
  ~CFCLEAR,,1
  ~CFACTIV,GETC,Y

  NomFile='CFVR3002'
  /TITLE, %NomFile%, Rectangular foundation stiffness (SHELL99)

! -------------------------------------------------------------------------------------
! 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,99

! 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,,,10
  ! 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

!--------------------------------------------------------------------------------------
! DATA CHECK
!--------------------------------------------------------------------------------------
! Data comparison number
  NComp    = 7
  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)  = 'Corner1'
  LABEL( 2)  = 'Corner2'
  LABEL( 3)  = 'Corner3'
  LABEL( 4)  = 'Side'
  LABEL( 5)  = 'Center'
  LABEL( 6)  = 'Maximum'
  LABEL( 7)  = 'Corner Elem'

! Correct data, obtained data
!--------------------------------------------------------------------------------------
   *GET,NVal,PARM,FndN1,DIM,X
   VALUE( 1,1) = 298e3   ! Corner
   VALUE( 2,1) = 298e3   ! Corner
   VALUE( 3,1) = 298e3   ! Corner
   VALUE( 4,1) = 203e3   ! Side
   VALUE( 5,1) = 149e3   ! Center
   VALUE( 6,1) = 297e3   ! Maximum value
   VALUE( 7,1) = 240e3   ! Corner. Value on element
   ~CFGET,VALUE( 1,2),FOUNDAT,1,NODE,EFS,,1   ! Corner
   ~CFGET,VALUE( 2,2),FOUNDAT,1,NODE,EFS,,2   ! Corner
   ~CFGET,VALUE( 3,2),FOUNDAT,1,NODE,EFS,,22  ! Corner
   ~CFGET,VALUE( 4,2),FOUNDAT,1,NODE,EFS,,32  ! Side
   ~CFGET,VALUE( 5,2),FOUNDAT,1,NODE,EFS,,211 ! Center
   ~CFGET,VALUE( 6,2),FOUNDAT,1,EFS,MAX       ! Maximum value
   ~CFGET,VALUE( 7,2),FOUNDAT,1,ELEM,EFS,,1   ! Corner. Value on element

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

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