!*************************************************************************************!
! TITLE: CFVR6104 - Load Test Pile
! SUBTITLE: Circular pile subjected to an axial force
!
! DESCRIPTION: A circular pile is located in a terrain which is formed
! DESCRIPTION: by a six different layers.
! DESCRIPTION: The diameter of the pile is 0.85 meters and its length is 16 meters.
! DESCRIPTION: The structure is submitted to an axial force of 104 Mp.
! DESCRIPTION: The aim of this example is to calculate the displacements
! DESCRIPTION: of the pile in the Z axis as well as to verify
! DESCRIPTION: that the different skin friction contributions
! DESCRIPTION: is according to the Chapter 17-B.1.2 "Single Pile Load Capacity"
! DESCRIPTION: and its formulation.
!
! ELEMENT TYPE: BEAM4, LINK8
! MODULES: GETC
! UNITS: User
! KEYWORD1: Materials
! KEYWORD2: Geotechnics
! KEYWORD3: Piles
!
!*************************************************************************************!
FINISH
~CFCLEAR,,1
 NomFile='CFVR6104'
/TITLE, %NomFile%,Load Test Pile

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

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

! Materials
~CFMP,1,LIB,CONCRETE,EHE,HA-30
~CFMP,1,EHE_C,Ej,,3e11,7
~CFMP,2,LIB,REINF,EHE,B500S
~CFMP,3,LIB,SOIL,,CL
~CFMP,4,LIB,SOIL,,CL
~CFMP,5,LIB,SOIL,,SW
~CFMP,6,LIB,SOIL,,SP
~CFMP,7,LIB,SOIL,,SP
~CFMP,8,LIB,SOIL,,SP

~CFMP,3,SOIL,PHIMCEFF,,5
~CFMP,3,SOIL,SPT,,1
~CFMP,3,SOIL,QU,,0.001
~CFMP,3,SOIL,CMCEFF,,0

~CFMP,4,SOIL,SPT,,20
~CFMP,4,SOIL,PHIMCEFF,,33
~CFMP,4,SOIL,QU,,13*2
~CFMP,4,SOIL,CMCEFF,,4

~CFMP,5,SOIL,SPT,,6
~CFMP,5,SOIL,PHIMCEFF,,29
~CFMP,5,SOIL,CMCEFF,,0

~CFMP,6,SOIL,SPT,,19
~CFMP,6,SOIL,PHIMCEFF,,33
~CFMP,6,SOIL,CMCEFF,,0

~CFMP,7,SOIL,SPT,,36
~CFMP,7,SOIL,PHIMCEFF,,38
~CFMP,7,SOIL,CMCEFF,,0

~CFMP,8,SOIL,SPT,,15
~CFMP,8,SOIL,PHIMCEFF,,33
~CFMP,8,SOIL,CMCEFF,,0
~CFMP,8,SOIL,QU,,150000

! Terrains
~TERDEF,1,NEW,6,,Z,4,,,,,Load case
~TERDEF,1,LAYER,1,3,2.5,,DB
~TERDEF,1,LAYER,2,4,2.5,,DB
~TERDEF,1,LAYER,3,5,4.5,,DB
~TERDEF,1,LAYER,4,6,3.5,,DB
~TERDEF,1,LAYER,5,7,5,,DB
~TERDEF,1,LAYER,6,8,50,,DB


~GTPD,CFVR6104,ltCF


!---------------------------------------------------------------------------------------
! DATA CHECK
!---------------------------------------------------------------------------------------
! Data comparison number
  NComp = 7
  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) = 'Displacements_Z'
  LABEL(2,1) = 'Skin Friction'

! Correct values
  VALUE(1,1)= 0.0043
  VALUE(2,1)= 0.0/0.981
  VALUE(3,1)= 6.1/0.981
  VALUE(4,1)= 1.2/0.981
  VALUE(5,1)= 3.8/0.981
  VALUE(6,1)= 7.2/0.981
  VALUE(7,1)= 3.0/0.981

! Comparison
  Node1 = NODE(0,0,4)
  *GET,VALUE(1,2),NODE,Node1,U,Z
  VALUE(2,2)=FS_(1)
  VALUE(3,2)=FS_(2)
  VALUE(4,2)=FS_(3)
  VALUE(5,2)=FS_(4)
  VALUE(6,2)=FS_(5)
  VALUE(7,2)=FS_(6)

  ! Warning and error tolerances

  TOLER( 1, 1)= 0.5E-3 $ TOLER( 1, 2)= 0.5E-3
  TOLER( 2, 1)= 0.05   $ TOLER( 2, 2)= 0.05
  TOLER( 3, 1)= 0.05   $ TOLER( 3, 2)= 0.05
  TOLER( 4, 1)= 0.05   $ TOLER( 4, 2)= 0.05
  TOLER( 5, 1)= 0.01   $ TOLER( 5, 2)= 0.01
  TOLER( 6, 1)= 0.01   $ TOLER( 6, 2)= 0.01
  TOLER( 7, 1)= 0.01   $ TOLER( 7, 2)= 0.01


! -----------------------------------------------------------------------------------
! Results Comparison
! -----------------------------------------------------------------------------------
COMPARA.MAC