Left

CivilFEM Online Help

Official Websites: CivilFEM | CivilFEM Licenses | Ingeciber
CivilFEM powered by Marc: Installation | Documentation | Python Guide | Theory Manual | Verification Manual | PDF | Class 3
CivilFEM for Ansys: Installation | Documentation | Verification Manual | Instalación | Documentación | PDF | PDF (ES) | Class 3
CivilFEM Apps for Ansys Workbench: Installation | Licensing | FAQ | Reinforced Concrete | Structural Steel | Theory Manual | PDF | Class 3
Right

CFVR3203 - Slope stability with water pressure

Slope stability analysis, with water pressure and an impermeable layer.

A plane model as shown in the figure is used to check the slope stability analysis. Fellenius, Bishop, Janbu and Modified Janbu methods are used.
The model is made off four layers of materials with the following properties (from top to bottom):

  • Specific weight 1 = 12600 N/m2; Porosity 1 = 0.10
  • Specific weight 2 = 14400 N/m2
  • Specific weight 3 = 11200 N/m2; Porosity 3 = 0.40
  • Specific weight 4 = 9000 N/m2; Porosity 4 = 0.10
The terrain has a water table defined as a polygonal line. The second layer is considered impermeable, no water pore pressure will act on it (RUSI=0).
The only action on the structure is its self weight.

Element types used in the model: PLANE82

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

The example cannot be launched on an educational license.
KEYWORDS
Geotechnics
Slope Stability


Model Statistics
Number of elements 197
Number of nodes 650
Number of civil materials 4
Number of cross sections 0
Number of shell vertices 0

Log file: CFVR3203.DAT


  FINISH
  ~CFCLEAR,,1
  ~CFACTIV,GETC,Y
  NomFile='CFVR3203'
  /TITLE, %NomFile%, Slope stability with water pressure

! 1. Initial data
! --------------------------------------------------
! CivilFEM SETUP
  ~UNITS,SI

/PREP7
! 2. Materials
! --------------------------------------------------
~CFMP,1,LIB,SOIL,,OH
~CFMP,2,LIB,SOIL,,CH
~CFMP,3,LIB,SOIL,,SC
~CFMP,4,LIB,SOIL,,PH
~CFMP,1,SOIL,GAMD,,12600 ! OH
~CFMP,2,SOIL,GAMD,,14400 ! CH
~CFMP,3,SOIL,GAMD,,11200 ! SC
~CFMP,4,SOIL,GAMD,,9000  ! PH
*DO,II,1,4
  ~CFMP,II,SOIL,RUSI,,1
*ENDDO
~CFMP,2,SOIL,RUSI,,0

! 3. Elements
! --------------------------------------------------
ET,1,PLANE82

! 4. Model
! --------------------------------------------------
K, 1,-12,28
K, 2,  1,29
K, 3, 10,37
K, 4, 19,38
K, 5, 35,28
K, 6, 50,27
K, 7,  6,33
K, 8, 20,29
K, 9, 26,34
K,10, 18,24
K,11, 50,21
K,12,-12,19
K,13, 50,16
K,14,-12, 6
K,15, 50,10

A,7,8,9,4,3
A,2,7,8,9,5,6,11,10
A,1,2,10,11,13,12
A,15,13,12,14

AESIZE,ALL,3

*DO,I,1,4
  MAT,I
  AMESH,I
*ENDDO
/NUMBER,1
/PNUM,MAT,1
EPLOT
/IMAGE, SAVE, %NomFile%,BMP

! 5. Capturing Model
! --------------------------------------------------
pp1=node(kx(1),ky(1),kz(1))
pp2=node(kx(6),ky(6),kz(6))
pp3=node(kx(5),ky(5),kz(5))
~SLPIN,PP1,PP2,PP3

! 6. Circles Centres Grid Input
! --------------------------------------------------
K,16,21,41
K,17,35,32
K,18,42,43
k,19,28,52
~SLPCIRK,16,19,17,5,6

! 7. Tangents Input
! --------------------------------------------------
K,20,34,27
K,21,34,11
K,22, 0,11
K,23,18,37
~SLPTANK,20,23,22,21,10

! 8. Hydrostatic Pressure
! --------------------------------------------------
~SLPPWP,5,0,-12,27, 01,27, 16,35, 41,25, 50,25

! 9. Data Check
! --------------------------------------------------
! Data comparison number
  NComp    = 30*4
  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

! 10. Labels, Correct data, obtained data
! --------------------------------------------------
 *DIM,AUXVECT,ARRAY,10
 *DIM,AUXVEC2,ARRAY,10
/POST1
~SLPSOL,0 ! Fellenius
*DO,JJ,1,30
  LABEL(JJ) = 'FEL' %JJ%
  *DO,II,1,10
    ~CFGET,AUXVECT(II),SLOPE,,RES,F,,2,(JJ-1)*10 + II
  *ENDDO
  *VOPER,AUXVEC2,AUXVECT,MULT,AUXVECT
  *VFUN,AUXVEC2,SQRT,AUXVEC2
  *VFUN,AUXVEC2,ASORT,AUXVEC2
  VALUE(JJ,2) = AUXVEC2(1)
*ENDDO

VALUE(  1,1) =  2.592
VALUE(  2,1) =  1.971
VALUE(  3,1) =  1.000
VALUE(  4,1) =  2.504
VALUE(  5,1) =  2.486
VALUE(  6,1) =  3.172
VALUE(  7,1) =  1.963
VALUE(  8,1) =  1.773
VALUE(  9,1) =  1.609
VALUE( 10,1) =  1.634
VALUE( 11,1) =  1.857
VALUE( 12,1) =  2.785
VALUE( 13,1) =  1.720
VALUE( 14,1) =  1.283
VALUE( 15,1) =  1.415
VALUE( 16,1) =  1.279
VALUE( 17,1) =  1.646
VALUE( 18,1) =  2.702
VALUE( 19,1) =  1.598
VALUE( 20,1) =  1.282
VALUE( 21,1) =  1.419
VALUE( 22,1) =  1.624
VALUE( 23,1) =  2.074
VALUE( 24,1) =  3.119
VALUE( 25,1) =  1.538
VALUE( 26,1) =  1.558
VALUE( 27,1) =  1.463
VALUE( 28,1) =  1.711
VALUE( 29,1) =  2.280
VALUE( 30,1) =  4.802

~SLPSOL,1 ! Bishop
*DO,JJ,1,30
  LABEL(JJ+30) = 'BIS' %JJ+30%
  *DO,II,1,10
     ~CFGET,AUXVECT(II),SLOPE,,RES,F,,2,(JJ-1)*10 + II
  *ENDDO
  *VOPER,AUXVEC2,AUXVECT,MULT,AUXVECT
  *VFUN,AUXVEC2,SQRT,AUXVEC2
  *VFUN,AUXVEC2,ASORT,AUXVEC2
  VALUE(JJ+30,2) = AUXVEC2(1)
*ENDDO

VALUE( 31,1) =  2.765
VALUE( 32,1) =  2.174
VALUE( 33,1) =  1.500
VALUE( 34,1) =  2.605
VALUE( 35,1) =  2.570
VALUE( 36,1) =  3.363
VALUE( 37,1) =  2.309
VALUE( 38,1) =  1.853
VALUE( 39,1) =  1.758
VALUE( 40,1) =  1.784
VALUE( 41,1) =  2.056
VALUE( 42,1) =  3.194
VALUE( 43,1) =  2.072
VALUE( 44,1) =  1.745
VALUE( 45,1) =  1.610
VALUE( 46,1) =  1.674
VALUE( 47,1) =  2.033
VALUE( 48,1) =  3.203
VALUE( 49,1) =  1.924
VALUE( 50,1) =  1.658
VALUE( 51,1) =  1.533
VALUE( 52,1) =  1.720
VALUE( 53,1) =  2.313
VALUE( 54,1) =  3.494
VALUE( 55,1) =  1.821
VALUE( 56,1) =  1.599
VALUE( 57,1) =  1.511
VALUE( 58,1) =  1.817
VALUE( 59,1) =  2.486
VALUE( 60,1) =  5.074

~SLPSOL,2 ! Janbu
*DO,JJ,1,30
  LABEL(JJ+60) = 'JAN' %JJ+60%
  *DO,II,1,10
    ~CFGET,AUXVECT(II),SLOPE,,RES,F,,2,(JJ-1)*10 + II
  *ENDDO
  *VOPER,AUXVEC2,AUXVECT,MULT,AUXVECT
  *VFUN,AUXVEC2,SQRT,AUXVEC2
  *VFUN,AUXVEC2,ASORT,AUXVEC2
  VALUE(JJ+60,2) = AUXVEC2(1)
*ENDDO

VALUE( 61,1) =  2.511
VALUE( 62,1) =  1.907
VALUE( 63,1) =  1.400
VALUE( 64,1) =  2.553
VALUE( 65,1) =  2.485
VALUE( 66,1) =  3.081
VALUE( 67,1) =  2.154
VALUE( 68,1) =  1.704
VALUE( 69,1) =  1.500
VALUE( 70,1) =  1.540
VALUE( 71,1) =  1.830
VALUE( 72,1) =  2.642
VALUE( 73,1) =  1.956
VALUE( 74,1) =  1.608
VALUE( 75,1) =  1.458
VALUE( 76,1) =  1.463
VALUE( 77,1) =  1.685
VALUE( 78,1) =  2.571
VALUE( 79,1) =  1.831
VALUE( 80,1) =  1.549
VALUE( 81,1) =  1.413
VALUE( 82,1) =  1.567
VALUE( 83,1) =  1.982
VALUE( 84,1) =  2.945
VALUE( 85,1) =  1.745
VALUE( 86,1) =  1.510
VALUE( 87,1) =  1.409
VALUE( 88,1) =  1.646
VALUE( 89,1) =  2.192
VALUE( 90,1) =  4.505

~SLPSOL,3 ! Modified Janbu
*DO,JJ,1,30
  LABEL(JJ+90) = 'JMD' %JJ+90%
  *DO,II,1,10
    ~CFGET,AUXVECT(II),SLOPE,,RES,F,,2,(JJ-1)*10 + II
  *ENDDO
  *VOPER,AUXVEC2,AUXVECT,MULT,AUXVECT
  *VFUN,AUXVEC2,SQRT,AUXVEC2
  *VFUN,AUXVEC2,ASORT,AUXVEC2
  VALUE(JJ+90,2) = AUXVEC2(1)
*ENDDO

VALUE( 91,1) =  2.72
VALUE( 92,1) =  2.98
VALUE( 93,1) =  1.63
VALUE( 94,1) =  2.77
VALUE( 95,1) =  2.69
VALUE( 96,1) =  3.38
VALUE( 97,1) =  2.30
VALUE( 98,1) =  1.81
VALUE( 99,1) =  1.55
VALUE(100,1) =  1.63
VALUE(101,1) =  2.02
VALUE(102,1) =  2.82
VALUE(103,1) =  2.07
VALUE(104,1) =  1.71
VALUE(105,1) =  1.54
VALUE(106,1) =  1.60
VALUE(107,1) =  1.86
VALUE(108,1) =  2.83
VALUE(109,1) =  1.93
VALUE(110,1) =  1.64
VALUE(111,1) =  1.50
VALUE(112,1) =  1.68
VALUE(113,1) =  2.17
VALUE(114,1) =  3.19
VALUE(115,1) =  1.83
VALUE(116,1) =  1.59
VALUE(117,1) =  1.49
VALUE(118,1) =  1.76
VALUE(119,1) =  2.37
VALUE(120,1) =  4.84

! Warning and error tolerances
  *DO,II,1,NComp
    TOLER(II,1) = VALUE(II,1)*0.03
    TOLER(II,2) = VALUE(II,1)*0.03
  *ENDDO
  ! Known exceptions
  I=  1 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I=  2 $ TOLER(  I,1) = VALUE(I,1)*0.35  $  TOLER(  I,2) = VALUE(I,1)*0.35
  I=  8 $ TOLER(  I,1) = VALUE(I,1)*0.40  $  TOLER(  I,2) = VALUE(I,1)*0.40
  I=  9 $ TOLER(  I,1) = VALUE(I,1)*0.40  $  TOLER(  I,2) = VALUE(I,1)*0.40
  I= 10 $ TOLER(  I,1) = VALUE(I,1)*0.40  $  TOLER(  I,2) = VALUE(I,1)*0.40
  I= 12 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 15 $ TOLER(  I,1) = VALUE(I,1)*0.25  $  TOLER(  I,2) = VALUE(I,1)*0.25
  I= 32 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 33 $ TOLER(  I,1) = VALUE(I,1)*0.30  $  TOLER(  I,2) = VALUE(I,1)*0.30
  I= 38 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 39 $ TOLER(  I,1) = VALUE(I,1)*0.20  $  TOLER(  I,2) = VALUE(I,1)*0.20
  I= 40 $ TOLER(  I,1) = VALUE(I,1)*0.20  $  TOLER(  I,2) = VALUE(I,1)*0.20
  I= 44 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 45 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 46 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 50 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 55 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 69 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 70 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 76 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 77 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 92 $ TOLER(  I,1) = VALUE(I,1)*0.35  $  TOLER(  I,2) = VALUE(I,1)*0.35
  I= 93 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I=  3 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 63 $ TOLER(  I,1) = VALUE(I,1)*0.30  $  TOLER(  I,2) = VALUE(I,1)*0.30

! 11. Results comparison
! --------------------------------------------------
  COMPARA.MAC

Results

LabelTargetCivilFEMRatioTolerance
'FEL'1 2.592 2.4851.041 0.3888
'FEL'2 1.971 1.35161.314 0.6898
'FEL'3 1 0.951161.049 0.15
'FEL'4 2.504 2.51640.995 0.07512
'FEL'5 2.486 2.50050.994 0.07458
'FEL'6 3.172 3.19750.992 0.09516
'FEL'7 1.963 1.97640.993 0.05889
'FEL'8 1.773 1.28261.277 0.7092
'FEL'9 1.609 0.975831.394 0.6436
'FEL'10 1.634 1.08781.334 0.6536
'FEL'11 1.857 1.87530.990 0.05571
'FEL'12 2.785 2.68731.035 0.4178
'FEL'13 1.72 1.74060.988 0.0516
'FEL'14 1.283 1.26051.018 0.03849
'FEL'15 1.415 1.12981.202 0.3538
'FEL'16 1.279 1.28730.994 0.03837
'FEL'17 1.646 1.66680.987 0.04938
'FEL'18 2.702 2.72870.990 0.08106
'FEL'19 1.598 1.61850.987 0.04794
'FEL'20 1.282 1.2970.988 0.03846
'FEL'21 1.419 1.42970.993 0.04257
'FEL'22 1.624 1.62830.997 0.04872
'FEL'23 2.074 2.08680.994 0.06222
'FEL'24 3.119 3.13440.995 0.09357
'FEL'25 1.538 1.55650.988 0.04614
'FEL'26 1.558 1.55930.999 0.04674
'FEL'27 1.463 1.46361.000 0.04389
'FEL'28 1.711 1.71450.998 0.05133
'FEL'29 2.28 2.2890.996 0.0684
'FEL'30 4.802 4.830.994 0.1441
'BIS'31 2.765 2.77560.996 0.08295
'BIS'32 2.174 2.12791.021 0.3261
'BIS'33 1.5 1.74820.835 0.45
'BIS'34 2.605 2.62230.993 0.07815
'BIS'35 2.57 2.5920.991 0.0771
'BIS'36 3.363 3.39630.990 0.1009
'BIS'37 2.309 2.31370.998 0.06927
'BIS'38 1.853 1.88250.984 0.2779
'BIS'39 1.758 1.63171.072 0.3516
'BIS'40 1.784 1.65821.071 0.3568
'BIS'41 2.056 2.07530.991 0.06168
'BIS'42 3.194 3.23280.988 0.09582
'BIS'43 2.072 2.0760.998 0.06216
'BIS'44 1.745 1.75070.997 0.2617
'BIS'45 1.61 1.6011.006 0.2415
'BIS'46 1.674 1.6711.002 0.2511
'BIS'47 2.033 2.05240.991 0.06099
'BIS'48 3.203 3.23560.990 0.09609
'BIS'49 1.924 1.92750.998 0.05772
'BIS'50 1.658 1.66290.997 0.2487
'BIS'51 1.533 1.54120.995 0.04599
'BIS'52 1.72 1.73070.994 0.0516
'BIS'53 2.313 2.32640.994 0.06939
'BIS'54 3.494 3.51630.994 0.1048
'BIS'55 1.821 1.82420.998 0.2732
'BIS'56 1.599 1.6040.997 0.04797
'BIS'57 1.511 1.51790.996 0.04533
'BIS'58 1.817 1.82610.995 0.05451
'BIS'59 2.486 2.50010.994 0.07458
'BIS'60 5.074 5.11010.993 0.1522
'JAN'61 2.511 2.52330.995 0.07533
'JAN'62 1.907 1.88561.011 0.05721
'JAN'63 1.4 1.54160.899 0.42
'JAN'64 2.553 2.57850.990 0.07659
'JAN'65 2.485 2.50940.990 0.07455
'JAN'66 3.081 3.1230.986 0.09243
'JAN'67 2.154 2.1560.999 0.06462
'JAN'68 1.704 1.7070.998 0.05112
'JAN'69 1.5 1.47611.016 0.225
'JAN'70 1.54 1.49011.032 0.231
'JAN'71 1.83 1.8690.979 0.0549
'JAN'72 2.642 2.58241.023 0.07926
'JAN'73 1.956 1.9590.999 0.05868
'JAN'74 1.608 1.60861.000 0.04824
'JAN'75 1.458 1.46080.998 0.04374
'JAN'76 1.463 1.47120.994 0.2195
'JAN'77 1.685 1.70530.988 0.2527
'JAN'78 2.571 2.59720.990 0.07713
'JAN'79 1.831 1.83380.999 0.05493
'JAN'80 1.549 1.54951.000 0.04647
'JAN'81 1.413 1.41361.000 0.04239
'JAN'82 1.567 1.57130.997 0.04701
'JAN'83 1.982 1.99940.991 0.05946
'JAN'84 2.945 2.95120.998 0.08835
'JAN'85 1.745 1.7461.000 0.05235
'JAN'86 1.51 1.51150.999 0.0453
'JAN'87 1.409 1.40931.000 0.04227
'JAN'88 1.646 1.65060.997 0.04938
'JAN'89 2.192 2.20530.994 0.06576
'JAN'90 4.505 4.53580.993 0.1351
'JMD'91 2.72 2.71781.001 0.0816
'JMD'92 2.98 1.98441.334 1.043
'JMD'93 1.63 1.62991.000 0.2445
'JMD'94 2.77 2.77131.000 0.0831
'JMD'95 2.69 2.69570.998 0.0807
'JMD'96 3.38 3.37931.000 0.1014
'JMD'97 2.3 2.29721.001 0.069
'JMD'98 1.81 1.80621.002 0.0543
'JMD'99 1.55 1.55540.997 0.0465
'JMD'100 1.63 1.63120.999 0.0489
'JMD'101 2.02 2.02520.997 0.0606
'JMD'102 2.82 2.82470.998 0.0846
'JMD'103 2.07 2.07150.999 0.0621
'JMD'104 1.71 1.70871.001 0.0513
'JMD'105 1.54 1.53861.001 0.0462
'JMD'106 1.6 1.60330.998 0.048
'JMD'107 1.86 1.86140.999 0.0558
'JMD'108 2.83 2.83011.000 0.0849
'JMD'109 1.93 1.92841.001 0.0579
'JMD'110 1.64 1.63721.002 0.0492
'JMD'111 1.5 1.49761.002 0.045
'JMD'112 1.68 1.68170.999 0.0504
'JMD'113 2.17 2.16781.001 0.0651
'JMD'114 3.19 3.18751.001 0.0957
'JMD'115 1.83 1.82811.001 0.0549
'JMD'116 1.59 1.5881.001 0.0477
'JMD'117 1.49 1.48751.002 0.0447
'JMD'118 1.76 1.76230.999 0.0528
'JMD'119 2.37 2.36791.001 0.0711
'JMD'120 4.84 4.84161.000 0.1452

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