CFVR1608 - Active Pressure acting on 2D solids

Element types used in the model: PLANE42

Needed CivilFEM Modules:
Geotechnical Module
Bridges and Civil Non Linearities Module
Advanced Prestressed Concrete Module

! DESCRUPTION: of 2.0 m.
  FINISH
  ~CFCLEAR,,1
  NomFile='CFVR1608'
  /TITLE, %NomFile%, Active Pressure on 2D Solids
  ~CFACTIV,GETC,Y ! Geotechnical module activation

! ---------------------------------------------------------------------------------
! Model definition
! ---------------------------------------------------------------------------------
! Civil Setup
~UNITS,,LENG,M
~UNITS,,TIME,S
~UNITS,,FORC,MP

/PREP7 ! Preprocessor
  ET,1,Plane42

  N,1,0,0
  N,2,0,1
  N,3,0.5774,2
  N,4,2*0.5774,1
  N,5,2*0.5774,0
  N,6,0.5774/2,1.5
  N,7,3*0.5774/2,1.5
  N,8,0.5774,0
  N,9,0.5774,1
  TYPE,1
  E,1,8,9,2
  E,2,9,6,6
  E,6,9,3,3
  E,3,9,7,7
  E,4,7,9,9
  E,8,5,4,9
  EPLOT
  /IMAGE, SAVE, %NomFile%,BMP

  DELTA= 15        ! Friction angle between the ground and wall, in degrees
  BETA=  10        ! Slope angle of the ground, in degrees
  Q=      0.1      ! Uniform surface load on the ground per unit of slope length
  EAKCN=  0        ! Reference number of the coordinate system in which the ground is defined
  EAH=    5.0      ! Ground level next to the wall
  WKEY=   1        ! Code identifying the influence of the water
  WH=     3.0      ! Water Table
  ~CFMP,1,LIB,SOIL,,SW
  ~CFMP,1,SOIL    ,N       ,,0.3
  ~CFMP,1,SOIL    ,GAMD    ,,2.0
  ~CFMP,1,SOIL    ,PHIMCEFF,,30
  ~TERDEF,1,NEW,2,EAKCN,Y,EAH,BETA, Q, WKEY, WH, Terrain 1
  ~TERDEF,1,LAYER,1,1,1.0,0,CLMB,DELTA
  ~TERDEF,1,LAYER,2,1,2.0,0,CLMB,DELTA

  NSEL,S,NODE,,1,7
  ~ETHSF,ACT,1      ! Active Pressure
  NSEL,ALL

/POST1
! ---------------------------------------------------------------------------------
! DATA CHECK
! ---------------------------------------------------------------------------------
! Data comparison number
  NComp = 20
  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) ='P1'
  LABEL( 2) ='P2'
  LABEL( 3) ='P3'
  LABEL( 4) ='P4'
  LABEL( 5) ='P5'
  LABEL( 6) ='P6'
  LABEL( 7) ='Fx1'
  LABEL( 8) ='Fy1'
  LABEL( 9) ='Fx2'
  LABEL(10) ='Fy2'
  LABEL(11) ='Fx3'
  LABEL(12) ='Fy3'
  LABEL(13) ='Fx4'
  LABEL(14) ='Fy4'
  LABEL(15) ='Fx5'
  LABEL(16) ='Fy5'
  LABEL(17) ='Fx6'
  LABEL(18) ='Fy6'
  LABEL(19) ='Fx7'
  LABEL(20) ='Fy7'

! Correct Values
  VALUE( 1,1)=  5.6517
  VALUE( 2,1)=  6.1115
  VALUE( 3,1)=  5.2126
  VALUE( 4,1)=  4.3137
  VALUE( 5,1)=  5.2126
  VALUE( 6,1)=  4.2208
  VALUE( 7,1)=  0.0
  VALUE( 8,1)= -0.3553
  VALUE( 9,1)= -0.1590
  VALUE(10,1)= -0.2975-.2754
  VALUE(11,1)=  0.0
  VALUE(12,1)= -0.4440
  VALUE(13,1)= -VALUE( 9,1)
  VALUE(14,1)=  VALUE(10,1)
  VALUE(15,1)= -VALUE( 7,1)
  VALUE(16,1)=  VALUE( 8,1)
  VALUE(17,1)= -.1436-.1436
  VALUE(18,1)= -.2487-.2487
  VALUE(19,1)= -VALUE(17,1)
  VALUE(20,1)=  VALUE(18,1)

! Obtained values
  *GET,VALUE( 1,2),ELEM,1,PRES,4
  *GET,VALUE( 2,2),ELEM,2,PRES,4
  *GET,VALUE( 3,2),ELEM,3,PRES,4
  *GET,VALUE( 4,2),ELEM,4,PRES,4
  *GET,VALUE( 5,2),ELEM,5,PRES,1
  *GET,VALUE( 6,2),ELEM,6,PRES,2
  *GET,VALUE( 7,2),NODE,1,F,FX
  *GET,VALUE( 8,2),NODE,1,F,FY
  *GET,VALUE( 9,2),NODE,2,F,FX
  *GET,VALUE(10,2),NODE,2,F,FY
  *GET,VALUE(11,2),NODE,3,F,FX
  *GET,VALUE(12,2),NODE,3,F,FY
  *GET,VALUE(13,2),NODE,4,F,FX
  *GET,VALUE(14,2),NODE,4,F,FY
  *GET,VALUE(15,2),NODE,5,F,FX
  *GET,VALUE(16,2),NODE,5,F,FY
  *GET,VALUE(17,2),NODE,6,F,FX
  *GET,VALUE(18,2),NODE,6,F,FY
  *GET,VALUE(19,2),NODE,7,F,FX
  *GET,VALUE(20,2),NODE,7,F,FY

! Warning and error tolerances
  *DO,I,1,NComp
    TOLER(I, 1)= 1E-4   $   TOLER(I, 2)= 1E-4
  *ENDDO

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