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

CFVR3213 - Slope stability with water pressure (PLANE 42)

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.
This example is the same as CFVR3203 using PLANE42 elements instead of PLANE82 elements.

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
Geotechnics
Slope Stability


Model Statistics
Number of elements 198
Number of nodes 228
Number of civil materials 4
Number of cross sections 0
Number of shell vertices 0

Log file: CFVR3213.DAT


  FINISH
  ~CFCLEAR,,1
  ~CFACTIV,GETC,Y
  NomFile='CFVR3213'
  /TITLE, %NomFile%, Slope stability with water pressure (PLANE 42)

! -------------------------------------------------------------------------------------
! 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,PLANE42

! 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) =  2.000
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) =  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.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=  3 $ TOLER(  I,1) = VALUE(I,1)*0.10  $  TOLER(  I,2) = VALUE(I,1)*0.10
  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= 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= 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.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= 98 $ 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
  I=100 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15
  I=106 $ TOLER(  I,1) = VALUE(I,1)*0.15  $  TOLER(  I,2) = VALUE(I,1)*0.15

! 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.1
'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 2 1.74821.126 0.6
'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 2.662 1.54161.421 100
'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.7464 2.71781.010 0.08239
'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.08315
'JMD'95 2.6953 2.69571.000 0.08086
'JMD'96 3.3792 3.37931.000 0.1014
'JMD'97 2.3214 2.29721.010 0.06964
'JMD'98 1.8567 1.80621.027 0.2785
'JMD'99 1.6963 1.55541.083 0.2544
'JMD'100 1.694 1.63121.037 0.2541
'JMD'101 2.0248 2.02521.000 0.06074
'JMD'102 2.7766 2.82470.983 0.0833
'JMD'103 2.0937 2.07151.011 0.06281
'JMD'104 1.7333 1.70871.014 0.052
'JMD'105 1.5749 1.53861.023 0.04725
'JMD'106 1.6601 1.60331.034 0.249
'JMD'107 1.8581 1.86140.998 0.05574
'JMD'108 2.8294 2.83011.000 0.08488
'JMD'109 1.9503 1.92841.011 0.05851
'JMD'110 1.6596 1.63721.014 0.04979
'JMD'111 1.512 1.49761.010 0.04536
'JMD'112 1.6852 1.68171.002 0.05056
'JMD'113 2.1814 2.16781.006 0.06544
'JMD'114 3.1875 3.18751.000 0.09562
'JMD'115 1.8493 1.82811.012 0.05548
'JMD'116 1.6089 1.5881.013 0.04827
'JMD'117 1.5005 1.48751.009 0.04501
'JMD'118 1.7666 1.76231.002 0.053
'JMD'119 2.3678 2.36791.000 0.07103
'JMD'120 4.8417 4.84161.000 0.1453

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