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

CFVR3210 - Slope stability with multiple water pressure lines

Slope stability analysis, with multiple water pressure lines 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 multiple polygonal lines to define several pressure levels, not following an hydrostatic pressure law. 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: CFVR3210.DAT


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

! -------------------------------------------------------------------------------------
! 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
~SLPPWP,7,10,-12,17, 01,17, 16,25, 41,15, 50,15
~SLPPWP,8,25,-12, 2, 01, 2, 16,10, 41, 0, 50, 0

! 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) =  2.187
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.992
VALUE( 14,1) =  1.669
VALUE( 15,1) =  1.532
VALUE( 16,1) =  1.585
VALUE( 17,1) =  1.964
VALUE( 18,1) =  2.702
VALUE( 19,1) =  1.868
VALUE( 20,1) =  1.604
VALUE( 21,1) =  1.475
VALUE( 22,1) =  1.624
VALUE( 23,1) =  2.074
VALUE( 24,1) =  3.119
VALUE( 25,1) =  1.781
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) =  2.766
VALUE( 34,1) =  2.605
VALUE( 35,1) =  2.570
VALUE( 36,1) =  3.363
VALUE( 37,1) =  2.309
VALUE( 38,1) =  1.888
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.677
VALUE( 47,1) =  2.160
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) =  2.662
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.746352432
VALUE( 92,1) =  2.113819272
VALUE( 93,1) =  1.863033364
VALUE( 94,1) =  2.771605293
VALUE( 95,1) =  2.695323123
VALUE( 96,1) =  3.379178768
VALUE( 97,1) =  2.32136147
VALUE( 98,1) =  1.856680598
VALUE( 99,1) =  1.696272308
VALUE(100,1) =  1.693979842
VALUE(101,1) =  2.024759364
VALUE(102,1) =  2.776590768
VALUE(103,1) =  2.093720093
VALUE(104,1) =  1.733334959
VALUE(105,1) =  1.574949787
VALUE(106,1) =  1.66013212
VALUE(107,1) =  1.858052621
VALUE(108,1) =  2.82942852
VALUE(109,1) =  1.950316335
VALUE(110,1) =  1.659562513
VALUE(111,1) =  1.511967284
VALUE(112,1) =  1.685219813
VALUE(113,1) =  2.181414864
VALUE(114,1) =  3.187496398
VALUE(115,1) =  1.849349508
VALUE(116,1) =  1.608915474
VALUE(117,1) =  1.500471577
VALUE(118,1) =  1.766556967
VALUE(119,1) =  2.367809689
VALUE(120,1) =  4.841717173

! Warning and error tolerances
  *DO,II,1,NComp
    TOLER(II,1) = VALUE(II,1)*0.05
    TOLER(II,2) = VALUE(II,1)*0.05
  *ENDDO
! Known exceptions
  I=  2 $ TOLER(  I,1) = VALUE(I,1)*0.40  $  TOLER(  I,2) = VALUE(I,1)*0.40
  I=  3 $ TOLER(  I,1) = VALUE(I,1)*0.40  $  TOLER(  I,2) = VALUE(I,1)*0.40
  I=  7 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I=  8 $ TOLER(  I,1) = VALUE(I,1)*0.30  $  TOLER(  I,2) = VALUE(I,1)*0.30
  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= 13 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 14 $ TOLER(  I,1) = VALUE(I,1)*0.30  $  TOLER(  I,2) = VALUE(I,1)*0.30
  I= 15 $ TOLER(  I,1) = VALUE(I,1)*0.30  $  TOLER(  I,2) = VALUE(I,1)*0.30
  I= 16 $ TOLER(  I,1) = VALUE(I,1)*0.20  $  TOLER(  I,2) = VALUE(I,1)*0.20
  I= 17 $ TOLER(  I,1) = VALUE(I,1)*0.20  $  TOLER(  I,2) = VALUE(I,1)*0.20
  I= 19 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 20 $ TOLER(  I,1) = VALUE(I,1)*0.20  $  TOLER(  I,2) = VALUE(I,1)*0.20
  I= 25 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  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.60  $  TOLER(  I,2) = VALUE(I,1)*0.60
  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= 47 $ 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= 63 $ TOLER(  I,1) = 100  $  TOLER(  I,2) = 100 ! Ignored
  I= 69 $ 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= 92 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 93 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I= 99 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15

!--------------------------------------------------------------------------------------
! 12. Results comparison
!--------------------------------------------------------------------------------------
  COMPARA.MAC

Results

LabelTargetCivilFEMRatioTolerance
'FEL'1 2.592 2.4851.041 0.1296
'FEL'2 1.971 1.35161.314 0.7884
'FEL'3 1 0.951161.049 0.4
'FEL'4 2.504 2.51640.995 0.1252
'FEL'5 2.486 2.50050.994 0.1243
'FEL'6 3.172 3.19750.992 0.1586
'FEL'7 2.187 1.97641.096 0.3281
'FEL'8 1.773 1.28261.277 0.5319
'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.09285
'FEL'12 2.785 2.68731.035 0.4178
'FEL'13 1.992 1.74061.126 0.2988
'FEL'14 1.669 1.26051.245 0.5007
'FEL'15 1.532 1.12981.263 0.4596
'FEL'16 1.585 1.28731.188 0.317
'FEL'17 1.964 1.66681.151 0.3928
'FEL'18 2.702 2.72870.990 0.1351
'FEL'19 1.868 1.61851.134 0.2802
'FEL'20 1.604 1.2971.191 0.3208
'FEL'21 1.475 1.42971.031 0.07375
'FEL'22 1.624 1.62830.997 0.0812
'FEL'23 2.074 2.08680.994 0.1037
'FEL'24 3.119 3.13440.995 0.156
'FEL'25 1.781 1.55651.126 0.2672
'FEL'26 1.558 1.55930.999 0.0779
'FEL'27 1.463 1.46361.000 0.07315
'FEL'28 1.711 1.71450.998 0.08555
'FEL'29 2.28 2.2890.996 0.114
'FEL'30 4.802 4.830.994 0.2401
'BIS'31 2.765 2.77560.996 0.1383
'BIS'32 2.174 2.12791.021 0.3261
'BIS'33 2.766 1.74821.368 1.66
'BIS'34 2.605 2.62230.993 0.1303
'BIS'35 2.57 2.5920.991 0.1285
'BIS'36 3.363 3.39630.990 0.1681
'BIS'37 2.309 2.31370.998 0.1155
'BIS'38 1.888 1.88251.003 0.2832
'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.1028
'BIS'42 3.194 3.23280.988 0.1597
'BIS'43 2.072 2.0760.998 0.1036
'BIS'44 1.745 1.75070.997 0.2617
'BIS'45 1.61 1.6011.006 0.2415
'BIS'46 1.677 1.6711.004 0.2516
'BIS'47 2.16 2.05241.050 0.324
'BIS'48 3.203 3.23560.990 0.1601
'BIS'49 1.924 1.92750.998 0.0962
'BIS'50 1.658 1.66290.997 0.2487
'BIS'51 1.533 1.54120.995 0.07665
'BIS'52 1.72 1.73070.994 0.086
'BIS'53 2.313 2.32640.994 0.1157
'BIS'54 3.494 3.51630.994 0.1747
'BIS'55 1.821 1.82420.998 0.09105
'BIS'56 1.599 1.6040.997 0.07995
'BIS'57 1.511 1.51790.996 0.07555
'BIS'58 1.817 1.82610.995 0.09085
'BIS'59 2.486 2.50010.994 0.1243
'BIS'60 5.074 5.11010.993 0.2537
'JAN'61 2.511 2.52330.995 0.1255
'JAN'62 1.907 1.88561.011 0.09535
'JAN'63 2.662 1.54161.421 100
'JAN'64 2.553 2.57850.990 0.1277
'JAN'65 2.485 2.50940.990 0.1243
'JAN'66 3.081 3.1230.986 0.154
'JAN'67 2.154 2.1560.999 0.1077
'JAN'68 1.704 1.7070.998 0.0852
'JAN'69 1.5 1.47611.016 0.225
'JAN'70 1.54 1.49011.032 0.077
'JAN'71 1.83 1.8690.979 0.0915
'JAN'72 2.642 2.58241.023 0.1321
'JAN'73 1.956 1.9590.999 0.0978
'JAN'74 1.608 1.60861.000 0.0804
'JAN'75 1.458 1.46080.998 0.0729
'JAN'76 1.463 1.47120.994 0.2195
'JAN'77 1.685 1.70530.988 0.08425
'JAN'78 2.571 2.59720.990 0.1286
'JAN'79 1.831 1.83380.999 0.09155
'JAN'80 1.549 1.54951.000 0.07745
'JAN'81 1.413 1.41361.000 0.07065
'JAN'82 1.567 1.57130.997 0.07835
'JAN'83 1.982 1.99940.991 0.0991
'JAN'84 2.945 2.95120.998 0.1472
'JAN'85 1.745 1.7461.000 0.08725
'JAN'86 1.51 1.51150.999 0.0755
'JAN'87 1.409 1.40931.000 0.07045
'JAN'88 1.646 1.65060.997 0.0823
'JAN'89 2.192 2.20530.994 0.1096
'JAN'90 4.505 4.53580.993 0.2253
'JMD'91 2.7464 2.71781.010 0.1373
'JMD'92 2.1138 1.98441.061 0.3171
'JMD'93 1.863 1.62991.125 0.2795
'JMD'94 2.7716 2.77131.000 0.1386
'JMD'95 2.6953 2.69571.000 0.1348
'JMD'96 3.3792 3.37931.000 0.169
'JMD'97 2.3214 2.29721.010 0.1161
'JMD'98 1.8567 1.80621.027 0.09283
'JMD'99 1.6963 1.55541.083 0.2544
'JMD'100 1.694 1.63121.037 0.0847
'JMD'101 2.0248 2.02521.000 0.1012
'JMD'102 2.7766 2.82470.983 0.1388
'JMD'103 2.0937 2.07151.011 0.1047
'JMD'104 1.7333 1.70871.014 0.08667
'JMD'105 1.5749 1.53861.023 0.07875
'JMD'106 1.6601 1.60331.034 0.08301
'JMD'107 1.8581 1.86140.998 0.0929
'JMD'108 2.8294 2.83011.000 0.1415
'JMD'109 1.9503 1.92841.011 0.09752
'JMD'110 1.6596 1.63721.014 0.08298
'JMD'111 1.512 1.49761.010 0.0756
'JMD'112 1.6852 1.68171.002 0.08426
'JMD'113 2.1814 2.16781.006 0.1091
'JMD'114 3.1875 3.18751.000 0.1594
'JMD'115 1.8493 1.82811.012 0.09247
'JMD'116 1.6089 1.5881.013 0.08045
'JMD'117 1.5005 1.48751.009 0.07502
'JMD'118 1.7666 1.76231.002 0.08833
'JMD'119 2.3678 2.36791.000 0.1184
'JMD'120 4.8417 4.84161.000 0.2421

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