Left

CivilFEM Online Help

Right

CFVR1609 - Terrain weight acting on 2D solids

Weight of a terrain acting on a 2D solid structure.

A terrain surrounds a structure made up of shell elements. The purpose of this example is to apply the equivalent weight of the terrain on the structural model.
The structure is a wall with inclined walls as shown in the figure, with a total height
The terrain is defined by a single soil material with the following parameters:

  • Friction angle between the ground and wall: 15º
  • Slope angle of the ground: 10º
  • Uniform surface load on the ground per unit of slope length: 0.1 Tn/m2
  • Ground level: 5.0 m measured from the bottom of the structure
  • Water table level: 3.0 m measured from the bottom of the structure
  • Porosity: 0.3
  • Dry specific weight 2.0 Tn/m3
  • Internal friction angle: 30º

Element types used in the model: PLANE42

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

The example can be launched on an educational license.
KEYWORDS
Terrain


Model Statistics
Number of elements 6
Number of nodes 9
Number of civil materials 1
Number of cross sections 0
Number of shell vertices 0

Log file: CFVR1609.DAT


! DESCRUPTION: of 2.0 m.
  FINISH
  ~CFCLEAR,,1
  NomFile='CFVR1609'
  /TITLE, %NomFile%, Terrain Weight 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,WGH,1      ! Terrain weight
  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)=  3.0000
  VALUE( 2,1)=  3.6750
  VALUE( 3,1)=  3.0125
  VALUE( 4,1)=  2.3500
  VALUE( 5,1)=  3.0125
  VALUE( 6,1)=  2.0000
  VALUE( 7,1)=  0.0
  VALUE( 8,1)=  0.0
  VALUE( 9,1)= -0.4187
  VALUE(10,1)= -0.7253
  VALUE(11,1)=  0.0
  VALUE(12,1)= -1.1691
  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)= -0.7563
  VALUE(18,1)= -1.3099
  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-3   $   TOLER(I, 2)= 1E-3
  *ENDDO

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

Results

LabelTargetCivilFEMRatioTolerance
P1 3 31.000 0.001
P2 3.675 3.67521.000 0.001
P3 3.0125 3.01271.000 0.001
P4 2.35 2.35021.000 0.001
P5 3.0125 3.01271.000 0.001
P6 2 21.000 0.001
Fx1 0 00.000 0.001
Fy1 0-2.4493e-0160.000 0.001
Fx2 -0.4187 -0.41881.000 0.001
Fy2 -0.7253 -0.725331.000 0.001
Fx3 01.1102e-0160.000 0.001
Fy3 -1.1691 -1.16921.000 0.001
Fx4 0.4187 0.41881.000 0.001
Fy4 -0.7253 -0.725331.000 0.001
Fx5 0 00.000 0.001
Fy5 0-2.4493e-0160.000 0.001
Fx6 -0.7563 -0.756351.000 0.001
Fy6 -1.3099 -1.30991.000 0.001
Fx7 0.7563 0.756351.000 0.001
Fy7 -1.3099 -1.30991.000 0.001


Contains proprietary and confidential information of Ingeciber, S.A.