CFVR3207 - Slope stability
Slope stability analysis, comparing results of having punctual loads or surface loads on the slope.
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 a single material with a dry specific weight of 12600 N/m2. The terrain is considered dry. No water table is defined.
Apart from self weight, two different load scenarios have been considered and compared:
- Horizontal punctual force of -106 N and vertical punctual force of -106 N on the middle point of the surface of the slope.
- On the adjacent elements to the node on which the previous punctual forces where applied, a surface load is applied with a value equivalent to the punctual forces. The value of the surface load is 5*105 N/m2, since the two adjacent elements have 1 m sized sides.
Element types used in the model: PLANE82 Needed CivilFEM Modules: |
|
| Model Statistics | |
| Number of elements | 273 |
| Number of nodes | 906 |
| Number of civil materials | 1 |
| Number of cross sections | 0 |
| Number of shell vertices | 0 |
Log file: CFVR3207.DAT
FINISH ~CFCLEAR,,1 ~CFACTIV,GETC,Y NomFile='CFVR3207' /TITLE, %NomFile%, Slope stability (pressures vs forces) ! ------------------------------------------------------------------------------------- ! 1. Initial data ! ------------------------------------------------------------------------------------- ! CivilFEM SETUP ~UNITS,SI /PREP7 ! 2. Materials ! -------------------------------------------------- ~CFMP,1,LIB,SOIL,,OH ~CFMP,1,Soil,GAMD,,12600 ! OH ! 3. Elements ! -------------------------------------------------- ET,1,PLANE82 ! 4. Model ! -------------------------------------------------- K, 1, 0, 0 K, 2,150, 0 K, 3,150,30 K, 4, 70,30 K, 5, 51,49 K, 6, 50,50 K, 7, 49,51 K, 8, 30,70 K, 9, 0,70 A,1,2,3,4,5,6,7,8,9 MAT,1 ESIZE,5 AMESH,1 ! 5. Capturing Model ! -------------------------------------------------- pp1=node(kx(3),ky(3),kz(3)) pp2=node(kx(9),ky(9),kz(9)) pp3=node(kx(6),ky(6),kz(6)) pp5=node(kx(5),ky(5),kz(5)) pp7=node(kx(7),ky(7),kz(7)) Force = 1e6 F,pp3,FX,-Force F,pp3,FY,-Force ! 6. Circles Centres Grid Input ! -------------------------------------------------- K,21,75,77 K,22,77,75 K,23,77,79 ! 7. Tangents Input ! -------------------------------------------------- K,30,37,39 K,31,39,37 K,32, 0,11 K,33,18,27 ! 8. Solve ! -------------------------------------------------- ~SLPIN,PP1,PP2,PP3 ~SLPCIRK,21,22,23,3,1 ~SLPTANK,30,31,32,33,1 !-------------------------------------------------------------------------------------- ! 9. Data Check !-------------------------------------------------------------------------------------- ! Data comparison number NComp = 12 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 !-------------------------------------------------------------------------------------- /POST1 ~SLPSOL,0 LABEL( 1) = 'FEL 1' $ ~CFGET,VALUE( 1,1),SLOPE,,RES,F,,2,1 LABEL( 2) = 'FEL 2' $ ~CFGET,VALUE( 2,1),SLOPE,,RES,F,,2,2 LABEL( 3) = 'FEL 3' $ ~CFGET,VALUE( 3,1),SLOPE,,RES,F,,2,3 ~SLPSOL,1 LABEL( 4) = 'BIS 1' $ ~CFGET,VALUE( 4,1),SLOPE,,RES,F,,2,1 LABEL( 5) = 'BIS 2' $ ~CFGET,VALUE( 5,1),SLOPE,,RES,F,,2,2 LABEL( 6) = 'BIS 3' $ ~CFGET,VALUE( 6,1),SLOPE,,RES,F,,2,3 ~SLPSOL,2 LABEL( 7) = 'JAN 1' $ ~CFGET,VALUE( 7,1),SLOPE,,RES,F,,2,1 LABEL( 8) = 'JAN 2' $ ~CFGET,VALUE( 8,1),SLOPE,,RES,F,,2,2 LABEL( 9) = 'JAN 3' $ ~CFGET,VALUE( 9,1),SLOPE,,RES,F,,2,3 ~SLPSOL,3 LABEL( 10) = 'JMD 1' $ ~CFGET,VALUE(10,1),SLOPE,,RES,F,,2,1 LABEL( 11) = 'JMD 2' $ ~CFGET,VALUE(11,1),SLOPE,,RES,F,,2,2 LABEL( 12) = 'JMD 3' $ ~CFGET,VALUE(12,1),SLOPE,,RES,F,,2,3 /PREP7 FDELE,ALL,ALL NWPLAN,-1,PP3,PP7,144 CSYS,4 NSEL,S,LOC,Y,0 NSEL,R,LOC,X,-1*SQRT(2),1*SQRT(2) CSYS,0 SF,ALL,PRESS,FORCE/2,FORCE/2 ALLSEL /PSF,PRES,NORM,2,0,1 EPLOT /IMAGE, SAVE, %NomFile%,BMP ~SLPIN,PP1,PP2,PP3 ~SLPCIRK,21,22,23,3,1 ~SLPTANK,30,31,32,33,1 /POST1 ~SLPSOL,0 ~CFGET,VALUE( 1,2),SLOPE,,RES,F,,2, 1 ~CFGET,VALUE( 2,2),SLOPE,,RES,F,,2, 2 ~CFGET,VALUE( 3,2),SLOPE,,RES,F,,2, 3 ~SLPSOL,1 ~CFGET,VALUE( 4,2),SLOPE,,RES,F,,2, 1 ~CFGET,VALUE( 5,2),SLOPE,,RES,F,,2, 2 ~CFGET,VALUE( 6,2),SLOPE,,RES,F,,2, 3 ~SLPSOL,2 ~CFGET,VALUE( 7,2),SLOPE,,RES,F,,2, 1 ~CFGET,VALUE( 8,2),SLOPE,,RES,F,,2, 2 ~CFGET,VALUE( 9,2),SLOPE,,RES,F,,2, 3 ~SLPSOL,3 ~CFGET,VALUE( 10,2),SLOPE,,RES,F,,2, 1 ~CFGET,VALUE( 11,2),SLOPE,,RES,F,,2, 2 ~CFGET,VALUE( 12,2),SLOPE,,RES,F,,2, 3 ! Warning and error tolerances *DO,II,1,NComp TOLER(II,1) = VALUE(II,1)*0.006 TOLER(II,2) = VALUE(II,1)*0.006 *ENDDO !-------------------------------------------------------------------------------------- ! 11. Results comparison !-------------------------------------------------------------------------------------- COMPARA.MAC |
Results
| Label | Target | CivilFEM | Ratio | Tolerance |
| FEL 1 | 0.71125 | 0.70784 | 1.005 | 0.004268 |
| FEL 2 | 0.72469 | 0.72575 | 0.999 | 0.004348 |
| FEL 3 | 0.75316 | 0.75129 | 1.002 | 0.004519 |
| BIS 1 | 0.82141 | 0.81702 | 1.005 | 0.004928 |
| BIS 2 | 0.84579 | 0.84715 | 0.998 | 0.005075 |
| BIS 3 | 0.88963 | 0.88735 | 1.003 | 0.005338 |
| JAN 1 | 0.70148 | 0.69912 | 1.003 | 0.004209 |
| JAN 2 | 0.71799 | 0.71871 | 0.999 | 0.004308 |
| JAN 3 | 0.74162 | 0.74031 | 1.002 | 0.00445 |
| JMD 1 | 0.75999 | 0.75743 | 1.003 | 0.00456 |
| JMD 2 | 0.77897 | 0.77975 | 0.999 | 0.004674 |
| JMD 3 | 0.80563 | 0.80422 | 1.002 | 0.004834 |
Contains proprietary and confidential information of Ingeciber, S.A.