!********************************************************************************!
! TITLE: CFVR6103 - Pile Wailing
! SUBTITLE: Rectangular wailing sustained by a mesh of piles
!
! DESCRIPTION: A rectangular wailing (4.9x2.4x2.5 m),is sustained by 14 piles
! DESCRIPTION: separated 0.7 meters within each other.
! DESCRIPTION: The piles are located to 0.35 meters in relation to the slab.
! DESCRIPTION: The structure is submitted to:
! DESCRIPTION:  - Axial force : 2123 Mp
! DESCRIPTION:  - X axis moment:193 Mp
! DESCRIPTION:  - Y axis moment:53.6 Mp.
! DESCRIPTION: The aim of this example is to calculate
! DESCRIPTION: the maximum axial force supported by the piles.
!
! ELEMENT TYPE: BEAM4, SHELL43, SURF154, LINK8
! MODULES: GETC
! UNITS: User
! KEYWORD1: Materials
! KEYWORD2: Geotechnics
! KEYWORD3: Piles
!
!*********************************************************************************!
FINISH
~CFCLEAR,,1
 NomFile='CFVR6103'
/TITLE, %NomFile%,Rectangular Pile Wailing

! -----------------------------------------------------------------------------------
! CivilFEM Setup:Codes,units and modules
! -----------------------------------------------------------------------------------
~CODESEL,EC3-92,EHE-98,EHE-98,,EC8-04    ! Set EHE codes
~UNITS,,FORC,MP
~CFACTIV,GETC,Y

! -----------------------------------------------------------------------------------
! Parameters
! -----------------------------------------------------------------------------------
Loc_Pil_X=1.05
Loc_Pil_Y=-2.1
Loc_Pil_Z=1.25

! ------------------------------------------------------------------------------------
! Setup
! ------------------------------------------------------------------------------------
/PREP7

! Materials
~CFMP,1,LIB,CONCRETE,EHE,HA-30
~CFMP,2,LIB,REINF,EHE,B500S

~CFMP,3,LIB,SOIL,,SW
~CFMP,3,SOIL,SPT,,20
~CFMP,3,SOIL,PHIMCEFF,,33

~CFMP,4,LIB,SOIL,,SW
~CFMP,4,SOIL,SPT,,50
~CFMP,4,SOIL,PHIMCEFF,,41

! terrain
~TERDEF,1,NEW,2,,Z,,,,,,sand
~TERDEF,1,LAYER,1,3,6 ,-1 ,CLMB
~TERDEF,1,LAYER,2,4,25 ,-1 ,CLMB


~GTPD,CFVR6103,rwCF

!---------------------------------------------------------------------------------------
! DATA CHECK
!---------------------------------------------------------------------------------------
! Data comparison number
  NComp = 1
  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) = 'FX'

! Correct values
  VALUE(1,1)= -86.47

! Comparison
 /POST1
 ~CFSET,,1
 ~PLLSFOR,F,X,-1,
 NSEL,S,LOC,Z,Loc_Pil_Z
 NSEL,R,LOC,X,Loc_Pil_X
 NSEL,R,LOC,Y,Loc_Pil_Y
 ESLN,S
 ESEL,R,TYPE,,EPIL_
 *GET,ELE_AXIL,ELEM,,NUM,MAX
 ~CFGET,VALUE(1,2),ELEMENT,ELE_AXIL,FORCE,FX,J

! Warning and error tolerances
 TOLER( 1, 1)= 0.1E+01 $ TOLER( 1, 2)= 0.1E+01
! -----------------------------------------------------------------------------------
! Results Comparison
! -----------------------------------------------------------------------------------
COMPARA.MAC