CFVR3201 - Slope stability

Slope stability analysis, without water pressure.

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 dry specific weight (from top to bottom):

12600 N/m²
14400 N/m²
11200 N/m²
9000 N/m²

The terrain is considered dry. No water table is defined.
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: CFVR3201.DAT

  FINISH
  ~CFCLEAR,,1
  ~CFACTIV,GETC,Y
  NomFile='CFVR3201'
  /TITLE, %NomFile%, Slope stability

! 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

! 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. Capture the 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. 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

! 9. 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.731
VALUE(  2,1) =  2.075
VALUE(  3,1) =  2.704
VALUE(  4,1) =  2.504
VALUE(  5,1) =  2.486
VALUE(  6,1) =  3.172
VALUE(  7,1) =  2.308
VALUE(  8,1) =  1.854
VALUE(  9,1) =  1.675
VALUE( 10,1) =  1.669
VALUE( 11,1) =  1.859
VALUE( 12,1) =  2.785
VALUE( 13,1) =  2.104
VALUE( 14,1) =  1.738
VALUE( 15,1) =  1.560
VALUE( 16,1) =  1.590
VALUE( 17,1) =  1.964
VALUE( 18,1) =  2.702
VALUE( 19,1) =  1.968
VALUE( 20,1) =  1.664
VALUE( 21,1) =  1.500
VALUE( 22,1) =  1.628
VALUE( 23,1) =  2.082
VALUE( 24,1) =  3.119
VALUE( 25,1) =  1.870
VALUE( 26,1) =  1.613
VALUE( 27,1) =  1.486
VALUE( 28,1) =  1.716
VALUE( 29,1) =  2.279
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.903
VALUE( 32,1) =  2.277
VALUE( 33,1) =  2.766
VALUE( 34,1) =  2.605
VALUE( 35,1) =  2.570
VALUE( 36,1) =  3.362
VALUE( 37,1) =  2.428
VALUE( 38,1) =  1.967
VALUE( 39,1) =  1.787
VALUE( 40,1) =  1.787
VALUE( 41,1) =  2.058
VALUE( 42,1) =  3.194
VALUE( 43,1) =  2.182
VALUE( 44,1) =  1.812
VALUE( 45,1) =  1.637
VALUE( 46,1) =  1.681
VALUE( 47,1) =  2.160
VALUE( 48,1) =  3.203
VALUE( 49,1) =  2.022
VALUE( 50,1) =  1.717
VALUE( 51,1) =  1.557
VALUE( 52,1) =  1.724
VALUE( 53,1) =  2.313
VALUE( 54,1) =  3.494
VALUE( 55,1) =  1.910
VALUE( 56,1) =  1.653
VALUE( 57,1) =  1.534
VALUE( 58,1) =  1.821
VALUE( 59,1) =  2.485
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.629
VALUE( 62,1) =  1.986
VALUE( 63,1) =  2.662
VALUE( 64,1) =  2.553
VALUE( 65,1) =  2.485
VALUE( 66,1) =  3.079
VALUE( 67,1) =  2.262
VALUE( 68,1) =  1.775
VALUE( 69,1) =  1.577
VALUE( 70,1) =  1.577
VALUE( 71,1) =  1.835
VALUE( 72,1) =  2.641
VALUE( 73,1) =  2.060
VALUE( 74,1) =  1.672
VALUE( 75,1) =  1.484
VALUE( 76,1) =  1.536
VALUE( 77,1) =  1.913
VALUE( 78,1) =  2.575
VALUE( 79,1) =  1.925
VALUE( 80,1) =  1.607
VALUE( 81,1) =  1.438
VALUE( 82,1) =  1.572
VALUE( 83,1) =  1.990
VALUE( 84,1) =  2.944
VALUE( 85,1) =  1.831
VALUE( 86,1) =  1.563
VALUE( 87,1) =  1.432
VALUE( 88,1) =  1.651
VALUE( 89,1) =  2.191
VALUE( 90,1) =  4.504

~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.83152206
VALUE( 92,1) =  2.15323452
VALUE( 93,1) =  2.3585
VALUE( 94,1) =  2.74424848
VALUE( 95,1) =  2.67011681
VALUE( 96,1) =  3.3793
VALUE( 97,1) =  2.41044791
VALUE( 98,1) =  1.90157239
VALUE( 99,1) =  1.69783065
VALUE(100,1) =  1.71645682
VALUE(101,1) =  1.98914166
VALUE(102,1) =  2.87643124
VALUE(103,1) =  2.178675
VALUE(104,1) =  1.77587309
VALUE(105,1) =  1.58308579
VALUE(106,1) =  1.64740663
VALUE(107,1) =  2.06495719
VALUE(108,1) =  2.80620208
VALUE(109,1) =  2.0245434
VALUE(110,1) =  1.69667833
VALUE(111,1) =  1.52352205
VALUE(112,1) =  1.68250515
VALUE(113,1) =  2.1582348
VALUE(114,1) =  3.18031197
VALUE(115,1) =  1.91728114
VALUE(116,1) =  1.6424169
VALUE(117,1) =  1.50344162
VALUE(118,1) =  1.76286865
VALUE(119,1) =  2.35274589
VALUE(120,1) =  4.807777881

! 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=  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.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

!--------------------------------------------------------------------------------------
! 10. Results comparison
!--------------------------------------------------------------------------------------
  COMPARA.MAC

Results

Label	Target	CivilFEM	Ratio	Tolerance
'FEL'1	2.731	2.7365	0.998	0.1366
'FEL'2	2.075	2.0866	0.994	0.1038
'FEL'3	2.704	2.3729	1.122	0.4056
'FEL'4	2.504	2.5164	0.995	0.1252
'FEL'5	2.486	2.5005	0.994	0.1243
'FEL'6	3.172	3.1975	0.992	0.1586
'FEL'7	2.308	2.3114	0.999	0.1154
'FEL'8	1.854	1.857	0.998	0.0927
'FEL'9	1.675	1.6811	0.996	0.08375
'FEL'10	1.669	1.6812	0.993	0.08345
'FEL'11	1.859	1.8753	0.991	0.09295
'FEL'12	2.785	2.6873	1.035	0.1393
'FEL'13	2.104	2.1057	0.999	0.1052
'FEL'14	1.738	1.7396	0.999	0.0869
'FEL'15	1.56	1.5622	0.999	0.078
'FEL'16	1.59	1.5957	0.996	0.0795
'FEL'17	1.964	1.9742	0.995	0.0982
'FEL'18	2.702	2.7287	0.990	0.1351
'FEL'19	1.968	1.9692	0.999	0.0984
'FEL'20	1.664	1.6656	0.999	0.0832
'FEL'21	1.5	1.5019	0.999	0.075
'FEL'22	1.628	1.6326	0.997	0.0814
'FEL'23	2.082	2.0967	0.993	0.1041
'FEL'24	3.119	3.1344	0.995	0.156
'FEL'25	1.87	1.8712	0.999	0.0935
'FEL'26	1.613	1.6136	1.000	0.08065
'FEL'27	1.486	1.4876	0.999	0.0743
'FEL'28	1.716	1.7197	0.998	0.0858
'FEL'29	2.279	2.289	0.996	0.114
'FEL'30	4.802	4.83	0.994	0.2401
'BIS'31	2.903	2.9114	0.997	0.1452
'BIS'32	2.277	2.2907	0.994	0.1139
'BIS'33	2.766	2.7875	0.992	0.1383
'BIS'34	2.605	2.6223	0.993	0.1303
'BIS'35	2.57	2.592	0.991	0.1285
'BIS'36	3.362	3.3963	0.990	0.1681
'BIS'37	2.428	2.4322	0.998	0.1214
'BIS'38	1.967	1.975	0.996	0.09835
'BIS'39	1.787	1.803	0.991	0.08935
'BIS'40	1.787	1.8045	0.990	0.08935
'BIS'41	2.058	2.0753	0.992	0.1029
'BIS'42	3.194	3.2328	0.988	0.1597
'BIS'43	2.182	2.1857	0.998	0.1091
'BIS'44	1.812	1.8183	0.997	0.0906
'BIS'45	1.637	1.6474	0.994	0.08185
'BIS'46	1.681	1.6942	0.992	0.08405
'BIS'47	2.16	2.1748	0.993	0.108
'BIS'48	3.203	3.2356	0.990	0.1601
'BIS'49	2.022	2.0251	0.998	0.1011
'BIS'50	1.717	1.7228	0.997	0.08585
'BIS'51	1.557	1.5658	0.994	0.07785
'BIS'52	1.724	1.7347	0.994	0.0862
'BIS'53	2.313	2.3264	0.994	0.1157
'BIS'54	3.494	3.5163	0.994	0.1747
'BIS'55	1.91	1.9127	0.999	0.0955
'BIS'56	1.653	1.6575	0.997	0.08265
'BIS'57	1.534	1.5413	0.995	0.0767
'BIS'58	1.821	1.8311	0.995	0.09105
'BIS'59	2.485	2.5001	0.994	0.1243
'BIS'60	5.074	5.1101	0.993	0.2537
'JAN'61	2.629	2.6387	0.996	0.1315
'JAN'62	1.986	2.0033	0.991	0.0993
'JAN'63	2.662	2.3007	1.136	0.3993
'JAN'64	2.553	2.5785	0.990	0.1277
'JAN'65	2.485	2.5094	0.990	0.1243
'JAN'66	3.079	3.123	0.986	0.154
'JAN'67	2.262	2.2633	0.999	0.1131
'JAN'68	1.775	1.7782	0.998	0.08875
'JAN'69	1.577	1.586	0.994	0.07885
'JAN'70	1.577	1.6214	0.972	0.07885
'JAN'71	1.835	1.869	0.982	0.09175
'JAN'72	2.641	2.5824	1.022	0.1321
'JAN'73	2.06	2.0615	0.999	0.103
'JAN'74	1.672	1.6723	1.000	0.0836
'JAN'75	1.484	1.4867	0.998	0.0742
'JAN'76	1.536	1.5484	0.992	0.0768
'JAN'77	1.913	1.9304	0.991	0.09565
'JAN'78	2.575	2.5972	0.991	0.1288
'JAN'79	1.925	1.9274	0.999	0.09625
'JAN'80	1.607	1.6091	0.999	0.08035
'JAN'81	1.438	1.4382	1.000	0.0719
'JAN'82	1.572	1.5757	0.998	0.0786
'JAN'83	1.99	2.0098	0.990	0.0995
'JAN'84	2.944	2.9512	0.998	0.1472
'JAN'85	1.831	1.8323	0.999	0.09155
'JAN'86	1.563	1.5646	0.999	0.07815
'JAN'87	1.432	1.4331	0.999	0.0716
'JAN'88	1.651	1.6558	0.997	0.08255
'JAN'89	2.191	2.2053	0.994	0.1095
'JAN'90	4.504	4.5358	0.993	0.2252
'JMD'91	2.8315	2.8421	0.996	0.1416
'JMD'92	2.1532	2.1713	0.992	0.1077
'JMD'93	2.3585	2.5177	0.933	0.3538
'JMD'94	2.7442	2.7713	0.990	0.1372
'JMD'95	2.6701	2.6957	0.990	0.1335
'JMD'96	3.3793	3.3793	1.000	0.169
'JMD'97	2.4104	2.4114	1.000	0.1205
'JMD'98	1.9016	1.9047	0.998	0.09508
'JMD'99	1.6978	1.707	0.995	0.08489
'JMD'100	1.7165	1.7648	0.972	0.08582
'JMD'101	1.9891	2.0252	0.982	0.09946
'JMD'102	2.8764	2.8247	1.018	0.1438
'JMD'103	2.1787	2.18	0.999	0.1089
'JMD'104	1.7759	1.7763	1.000	0.08879
'JMD'105	1.5831	1.5857	0.998	0.07915
'JMD'106	1.6474	1.6609	0.992	0.08237
'JMD'107	2.065	2.0834	0.991	0.1032
'JMD'108	2.8062	2.8301	0.992	0.1403
'JMD'109	2.0245	2.0268	0.999	0.1012
'JMD'110	1.6967	1.6987	0.999	0.08483
'JMD'111	1.5235	1.5237	1.000	0.07618
'JMD'112	1.6825	1.6864	0.998	0.08413
'JMD'113	2.1582	2.1791	0.990	0.1079
'JMD'114	3.1803	3.1875	0.998	0.159
'JMD'115	1.9173	1.9185	0.999	0.09586
'JMD'116	1.6424	1.6438	0.999	0.08212
'JMD'117	1.5034	1.5126	0.994	0.07517
'JMD'118	1.7629	1.7679	0.997	0.08814
'JMD'119	2.3527	2.3679	0.994	0.1176
'JMD'120	4.8078	4.8416	0.993	0.2404

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