Left

CivilFEM Online Help

Official Websites: CivilFEM | CivilFEM Licenses | Ingeciber | Current downloads
CivilFEM for Ansys: Install | Doc | Verif | Instala | Doc(ES) | PDF | PDF(ES) | Support | Soporte | Class3
CivilFEM Apps for Ansys Workbench: Install | Lic | FAQ | Concrete | Steel | Theory | PDF | Support | Soporte | Class3
CivilFEM powered by Marc: Install | Capabilities | Cap(ES) | Doc | Python | Theory | Verif | PDF | Support | Soporte | Class3
Right

CFVR1607 - Active Pressure acting on 3D solids

Active Pressure of a terrain acting on a 3D solid structure.

A terrain surrounds a structure made up of shell elements. The purpose of this example is to apply the equivalent loads of the terrain on the structural model.
The structure is a wall with inclined walls as shown in the figure, with a total height of 2.0 m.
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º
Active pressure is applied.

Element types used in the model: SOLID45

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 18
Number of civil materials 1
Number of cross sections 0
Number of shell vertices 0

Log file: CFVR1607.DAT


  FINISH
  ~CFCLEAR,,1
  NomFile='CFVR1607'
  /TITLE, %NomFile%, Active Pressure acting on 3D Solids
  ~CFACTIV,GETC,Y ! Geotechnical module activation

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

/PREP7  ! Preprocessor
  ET,1,Solid45

  N, 1,0,0,0
  N, 2,0,1,0
  N, 3,0.5774,2,0
  N, 4,0.5774*2,1,0
  N, 5,0.5774*2,0,0
  N, 6,0.5774/2,1.5,0
  N, 7,0.5774/2*3,1.5,0
  N, 8,0.5774,0,0
  N, 9,0.5774,1,0
  N,11,0,0,1
  N,12,0,1,1
  N,13,0.5774,2,1
  N,14,0.5774*2,1,1
  N,15,0.5774*2,0,1
  N,16,0.5774/2,1.5,1
  N,17,0.5774/2*3,1.5,1
  N,18,0.5774,0,1
  N,19,0.5774,1,1
  TYPE,1
  E,1,2,9,8,11,12,19,18
  E,8,9,4,5,18,19,14,15
  E,2,9,6,6,12,19,16,16
  E,6,9,3,3,16,19,13,13
  E,3,9,7,7,13,19,17,17
  E,4,7,9,9,14,17,19,19
  /VIEW,1,-0.8,0.9,1.1
  EPLOT
  /IMAGE, SAVE, %NomFile%,BMP

  DELTA= 15        ! Angle of friction 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        ! Ground level next to the wall
  WKEY=   1        ! Code identifying the influence of the water
  WH=     3        ! 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
  NSEL,A,NODE,,11,17
  ~ETHSF,ACT,1  ! Active pressure
  NSEL,ALL

! ---------------------------------------------------------------------------------
! DATA CHECK
! ---------------------------------------------------------------------------------
! Data comparison number
  NComp = 42
  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) ='Pn1'
  LABEL( 2) ='Pn2'
  LABEL( 3) ='Pn3'
  LABEL( 4) ='Pn4'
  LABEL( 5) ='Pn5'
  LABEL( 6) ='Pn6'
  LABEL( 7) ='Fx 1'
  LABEL( 8) ='Fy 1'
  LABEL( 9) ='Fx 2'
  LABEL(10) ='Fy 2'
  LABEL(11) ='Fx 3'
  LABEL(12) ='Fy 3'
  LABEL(13) ='Fx 4'
  LABEL(14) ='Fy 4'
  LABEL(15) ='Fx 5'
  LABEL(16) ='Fy 5'
  LABEL(17) ='Fx 6'
  LABEL(18) ='Fy 6'
  LABEL(19) ='Fx 7'
  LABEL(20) ='Fy 7'
  LABEL(21) ='Fx 8'
  LABEL(22) ='Fy 8'
  LABEL(23) ='Fx 9'
  LABEL(24) ='Fy 9'
  LABEL(25) ='Fx11'
  LABEL(26) ='Fy11'
  LABEL(27) ='Fx12'
  LABEL(28) ='Fy12'
  LABEL(29) ='Fx13'
  LABEL(30) ='Fy13'
  LABEL(31) ='Fx14'
  LABEL(32) ='Fy14'
  LABEL(33) ='Fx15'
  LABEL(34) ='Fy15'
  LABEL(35) ='Fx16'
  LABEL(36) ='Fy16'
  LABEL(37) ='Fx17'
  LABEL(38) ='Fy17'
  LABEL(39) ='Fx18'
  LABEL(40) ='Fy18'
  LABEL(41) ='Fx19'
  LABEL(42) ='Fy19'

! Correct Values
  VALUE( 1,1)=  5.6517
  VALUE( 2,1)=  4.2208
  VALUE( 3,1)=  5.2126
  VALUE( 4,1)=  4.3137
  VALUE( 5,1)=  5.2126
  VALUE( 6,1)=  6.1115
  VALUE( 7,1)=  0.0000
  VALUE( 8,1)= -0.1776
  VALUE( 9,1)= -0.0795
  VALUE(10,1)= -0.2865
  VALUE(11,1)=  0.0000
  VALUE(12,1)= -0.2220
  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.1436
  VALUE(18,1)= -0.2487
  VALUE(19,1)= -VALUE(17,1)
  VALUE(20,1)=  VALUE(18,1)
  VALUE(21,1)=  0.0000
  VALUE(22,1)=  0.0000
  VALUE(23,1)=  0.0000
  VALUE(24,1)=  0.0000
  VALUE(25,1)=  VALUE( 7,1)
  VALUE(26,1)=  VALUE( 8,1)
  VALUE(27,1)=  VALUE( 9,1)
  VALUE(28,1)=  VALUE(10,1)
  VALUE(29,1)=  VALUE(11,1)
  VALUE(30,1)=  VALUE(12,1)
  VALUE(31,1)=  VALUE(13,1)
  VALUE(32,1)=  VALUE(14,1)
  VALUE(33,1)=  VALUE(15,1)
  VALUE(34,1)=  VALUE(16,1)
  VALUE(35,1)=  VALUE(17,1)
  VALUE(36,1)=  VALUE(18,1)
  VALUE(37,1)=  VALUE(19,1)
  VALUE(38,1)=  VALUE(20,1)
  VALUE(39,1)=  VALUE(21,1)
  VALUE(40,1)=  VALUE(22,1)
  VALUE(41,1)=  VALUE(23,1)
  VALUE(42,1)=  VALUE(24,1)

! Obtained values
  *GET,VALUE( 1,2),ELEM,1,PRES,2
  *GET,VALUE( 2,2),ELEM,2,PRES,4
  *GET,VALUE( 3,2),ELEM,3,PRES,5
  *GET,VALUE( 4,2),ELEM,4,PRES,5
  *GET,VALUE( 5,2),ELEM,5,PRES,5
  *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
  *GET,VALUE(21,2),NODE, 8,F,FX
  *GET,VALUE(22,2),NODE, 8,F,FY
  *GET,VALUE(23,2),NODE, 9,F,FX
  *GET,VALUE(24,2),NODE, 9,F,FY
  *GET,VALUE(25,2),NODE,11,F,FX
  *GET,VALUE(26,2),NODE,11,F,FY
  *GET,VALUE(27,2),NODE,12,F,FX
  *GET,VALUE(28,2),NODE,12,F,FY
  *GET,VALUE(29,2),NODE,13,F,FX
  *GET,VALUE(30,2),NODE,13,F,FY
  *GET,VALUE(31,2),NODE,14,F,FX
  *GET,VALUE(32,2),NODE,14,F,FY
  *GET,VALUE(33,2),NODE,15,F,FX
  *GET,VALUE(34,2),NODE,15,F,FY
  *GET,VALUE(35,2),NODE,16,F,FX
  *GET,VALUE(36,2),NODE,16,F,FY
  *GET,VALUE(37,2),NODE,17,F,FX
  *GET,VALUE(38,2),NODE,17,F,FY
  *GET,VALUE(39,2),NODE,18,F,FX
  *GET,VALUE(40,2),NODE,18,F,FY
  *GET,VALUE(41,2),NODE,19,F,FX
  *GET,VALUE(42,2),NODE,19,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
Pn1 5.6517 5.65171.000 0.001
Pn2 4.2208 4.22081.000 0.001
Pn3 5.2126 5.21251.000 0.001
Pn4 4.3137 4.31371.000 0.001
Pn5 5.2126 5.21251.000 0.001
Pn6 6.1115 6.11141.000 0.001
Fx 1 0 00.000 0.001
Fy 1 -0.1776 -0.177631.000 0.001
Fx 2 -0.0795 -0.0795121.000 0.001
Fy 2 -0.2865 -0.286471.000 0.001
Fx 3 01.3878e-0170.000 0.001
Fy 3 -0.222 -0.221971.000 0.001
Fx 4 0.0795 0.0795121.000 0.001
Fy 4 -0.2865 -0.286471.000 0.001
Fx 5 0 00.000 0.001
Fy 5 -0.1776 -0.177631.000 0.001
Fx 6 -0.1436 -0.14361.000 0.001
Fy 6 -0.2487 -0.248691.000 0.001
Fx 7 0.1436 0.14361.000 0.001
Fy 7 -0.2487 -0.248691.000 0.001
Fx 8 0 00.000 0.001
Fy 8 0 00.000 0.001
Fx 9 0 00.000 0.001
Fy 9 0 00.000 0.001
Fx11 0 00.000 0.001
Fy11 -0.1776 -0.177631.000 0.001
Fx12 -0.0795 -0.0795121.000 0.001
Fy12 -0.2865 -0.286471.000 0.001
Fx13 01.3878e-0170.000 0.001
Fy13 -0.222 -0.221971.000 0.001
Fx14 0.0795 0.0795121.000 0.001
Fy14 -0.2865 -0.286471.000 0.001
Fx15 0 00.000 0.001
Fy15 -0.1776 -0.177631.000 0.001
Fx16 -0.1436 -0.14361.000 0.001
Fy16 -0.2487 -0.248691.000 0.001
Fx17 0.1436 0.14361.000 0.001
Fy17 -0.2487 -0.248691.000 0.001
Fx18 0 00.000 0.001
Fy18 0 00.000 0.001
Fx19 0 00.000 0.001
Fy19 0 00.000 0.001

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