!***********************************************************************!
! TITLE: CFVR3003 - Circular Foundation Stiffness (SURF154)
! SUBTITLE:
!
! DESCRIPTION:	This example tests the stiffness of a circular foundation, modelled as
! DESCRIPTION:  a circular surface.
! DESCRIPTION:
! DESCRIPTION:  For this purpose, it tests the stiffness at the following
! DESCRIPTION:  locations:
! DESCRIPTION:
! DESCRIPTION:  - Sides
! DESCRIPTION:
! DESCRIPTION:  - Center
!
! ELEMENT TYPE: SURF154
! MODULES: GETC
! UNITS: SI
! KEYWORD1: Materials
! KEYWORD2: Geotechnics
! KEYWORD3: Material Models
!
!***********************************************************************!
  FINISH
  ~CFCLEAR,,1
  ~CFACTIV,GETC,Y

  NomFile='CFVR3003'
  /TITLE, %NomFile%, Circular foundation stiffness (SURF154)

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

! Model
! --------------------------------------------------
  ! Keypoints
   k,1,0,0
  ! Lines
   circle,1,10,,,,4
   lesize,all,,,10
  ! Area
   a,2,3,4,5
  ! Hard Point
   hptcreate,AREA,1,10,COORD,0,0,0
  ! Elements
   amesh,1

! 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    = 2
  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)  = 'Center'
  LABEL( 2)  = 'Side'

! Correct data, obtained data
!--------------------------------------------------------------------------------------
   VALUE( 1,1) = 56980     ! Center
   VALUE( 2,1) = 89500     ! Side
   ~CFGET,VALUE( 1,2),FOUNDAT,1,NODE,EFS,,187 ! Center
   ~CFGET,VALUE( 2,2),FOUNDAT,1,NODE,EFS,,1   ! Side

! Warning and error tolerances
  TOLER(  1,1) = 1E2   $  TOLER(  1,2) = 1E2
  TOLER(  2,1) = 1E2   $  TOLER(  2,2) = 1E2

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