Instance Models
This section transforms the problem defined in the problem description into one which is expressed in mathematical terms. It uses concrete symbols defined in the data definitions to replace the abstract symbols in the models identified in theoretical models and general definitions.
The goals GS:Identify-Crit-and-FS, GS:Determine-Normal-Forces, and GS:Determine-Shear-Forces are met by the simultaneous solution of IM:fctSfty, IM:nrmShrFor, and IM:intsliceFs. The goal GS:Identify-Crit-and-FS is also contributed to by IM:crtSlpId.
The Morgenstern-Price method is a vertical slice, limit equilibrium slope stability analysis method. Analysis is performed by breaking the assumed slip surface into a series of vertical slices of mass. Static equilibrium analysis is performed, using two force equations and one moment equation as in TM:equilibrium. The problem is statically indeterminate with only these 3 equations and one constitutive equation (the Mohr Coulomb shear strength of TM:mcShrStrgth) so the assumption GD:normShrR and corresponding equation GD:normShrR are used. The force equilibrium equations can be modified to be expressed only in terms of known physical values, as done in GD:resShearWO and GD:mobShearWO.
Refname | IM:fctSfty |
---|---|
Label | Factor of safety |
Input | xslope, yslope, ywt, c′, φ′, γdry, γsat, γw, xslip, yslip, constf |
Output | FS |
Input Constraints | |
Output Constraints | |
Equation | FS=n−1∑i=1Rin−1∏v=iΨv+Rnn−1∑i=1Tin−1∏v=iΨv+Tn |
Description |
|
Notes |
|
Source | chen2005 and karchewski2012 |
RefBy | IM:nrmShrFor, IM:intsliceFs, IM:fctSfty, FR:Display-Interslice-Shear-Forces, FR:Display-Interslice-Normal-Forces, FR:Display-Factor-of-Safety, and FR:Determine-Critical-Slip-Surface |
The mobilized shear force defined in GD:bsShrFEq can be substituted into the definition of mobilized shear force based on the factor of safety, from GD:mobShr yielding Equation (1) below:
(Wi−Xi−1+Xi+Ug,icos(βi)+Qicos(ωi))sin(αi)−(−KcWi−Gi+Gi−1−Hi+Hi−1+Ug,isin(βi)+Qisin(ωi))cos(αi)=N′itan(φ′)+c′Lb,iFS
An expression for the effective normal forces, N′, can be derived by substituting the normal forces equilibrium from GD:normForcEq into the definition for effective normal forces from GD:resShearWO. This results in Equation (2):
N′i=(Wi−Xi−1+Xi+Ug,icos(βi)+Qicos(ωi))cos(αi)+(−KcWi−Gi+Gi−1−Hi+Hi−1+Ug,isin(βi)+Qisin(ωi))sin(αi)−Ub,i
Substituting Equation (2) into Equation (1) gives:
(Wi−Xi−1+Xi+Ug,icos(βi)+Qicos(ωi))sin(αi)−(−KcWi−Gi+Gi−1−Hi+Hi−1+Ug,isin(βi)+Qisin(ωi))cos(αi)=((Wi−Xi−1+Xi+Ug,icos(βi)+Qicos(ωi))cos(αi)+(−KcWi−Gi+Gi−1−Hi+Hi−1+Ug,isin(βi)+Qisin(ωi))sin(αi)−Ub,i)tan(φ′)+c′Lb,iFS
Since the interslice shear forces X and interslice normal forces G are unknown, they are separated from the other terms as follows:
(Wi+Ug,icos(βi)+Qicos(ωi))sin(αi)−(−KcWi−Hi+Hi−1+Ug,isin(βi)+Qisin(ωi))cos(αi)−(−Gi+Gi−1)cos(αi)+(−Xi−1+Xi)sin(αi)=((Wi+Ug,icos(βi)+Qicos(ωi))cos(αi)+(−KcWi−Hi+Hi−1+Ug,isin(βi)+Qisin(ωi))sin(αi)+(−Gi+Gi−1)sin(αi)+(−Xi−1+Xi)cos(αi)−Ub,i)tan(φ′)+c′Lb,iFS
Applying assumptions A:Seismic-Force and A:Surface-Load, which state that the seismic coefficient and the external forces, respectively, are zero, allows for further simplification as shown below:
(Wi+Ug,icos(βi))sin(αi)−(−Hi+Hi−1+Ug,isin(βi))cos(αi)−(−Gi+Gi−1)cos(αi)+(−Xi−1+Xi)sin(αi)=((Wi+Ug,icos(βi))cos(αi)+(−Hi+Hi−1+Ug,isin(βi))sin(αi)+(−Gi+Gi−1)sin(αi)+(−Xi−1+Xi)cos(αi)−Ub,i)tan(φ′)+c′Lb,iFS
The definitions of GD:resShearWO and GD:mobShearWO are present in this equation, and thus can be replaced by Ri and Ti, respectively:
Ti+(−Xi−1+Xi)sin(αi)−(−Gi+Gi−1)cos(αi)=Ri+((−Xi−1+Xi)cos(αi)+(−Gi+Gi−1)sin(αi))tan(φ′)FS
The interslice shear forces X can be expressed in terms of the interslice normal forces G using A:Interslice-Norm-Shear-Forces-Linear and GD:normShrR, resulting in:
Ti+(−λfi−1Gi−1+λfiGi)sin(αi)−(−Gi+Gi−1)cos(αi)=Ri+((−λfi−1Gi−1+λfiGi)cos(αi)+(−Gi+Gi−1)sin(αi))tan(φ′)FS
Rearranging yields the following:
Gi((λficos(αi)−sin(αi))tan(φ′)−(λfisin(αi)+cos(αi))FS)=Gi−1((λfi−1cos(αi)−sin(αi))tan(φ′)−(λfi−1sin(αi)+cos(αi))FS)+FSTi−Ri
The definitions for Φ and Ψ from DD:convertFunc1 and DD:convertFunc2 simplify the above to Equation (3):
GiΦi=Ψi−1Gi−1Φi−1+FSTi−Ri
Versions of Equation (3) instantiated for slices 1 to n are shown below:
G1Φ1=Ψ0G0Φ0+FST1−R1
G2Φ2=Ψ1G1Φ1+FST2−R2
G3Φ3=Ψ2G2Φ2+FST3−R3
…
Gn−2Φn−2=Ψn−3Gn−3Φn−3+FSTn−2−Rn−2
Gn−1Φn−1=Ψn−2Gn−2Φn−2+FSTn−1−Rn−1
GnΦn=Ψn−1Gn−1Φn−1+FSTn−Rn
Applying A:Edge-Slices, which says that G0 and Gn are zero, results in the following special cases: Equation (8) for the first slice:
G1Φ1=FST1−R1
and Equation (9) for the nth slice:
−(FSTn−RnΨn−1)=Gn−1Φn−1
Substituting Equation (8) into Equation (4) yields Equation (10):
G2Φ2=Ψ1(FST1−R1)+FST2−R2
which can be substituted into Equation (5) to get Equation (11):
G3Φ3=Ψ2(Ψ1(FST1−R1)+FST2−R2)+FST3−R3
and so on until Equation (12) is obtained from Equation (7):
Gn−1Φn−1=Ψn−2(Ψn−3(Ψ1(FST1−R1)+FST2−R2)+FSTn−2−Rn−2)+FSTn−1−Rn−1
Equation (9) can then be substituted into the left-hand side of Equation (12), resulting in:
−(FSTn−RnΨn−1)=Ψn−2(Ψn−3(Ψ1(FST1−R1)+FST2−R2)+FSTn−2−Rn−2)+FSTn−1−Rn−1
This can be rearranged by multiplying both sides by Ψn−1 and then distributing the multiplication of each Ψ over addition to obtain:
−(FSTn−Rn)=Ψn−1Ψn−2Ψ1(FST1−R1)+Ψn−1Ψn−2Ψ2(FST2−R2)+Ψn−1(FSTn−1−Rn−1)
The multiplication of the Ψ terms can be further distributed over the subtractions, resulting in the equation having terms that each either contain an R or a T. The equation can then be rearranged so terms containing an R are on one side of the equality, and terms containing a T are on the other. The multiplication by the factor of safety is common to all of the T terms, and thus can be factored out, resulting in:
FS(Ψn−1Ψn−2Ψ1T1+Ψn−1Ψn−2Ψ2T2+Ψn−1Tn−1+Tn)=Ψn−1Ψn−2Ψ1R1+Ψn−1Ψn−2Ψ2R2+Ψn−1Rn−1+Rn
Isolating the factor of safety on the left-hand side and using compact notation for the products and sums yields Equation (13), which can also be seen in IM:fctSfty:
FS=n−1∑i=1Rin−1∏v=iΨv+Rnn−1∑i=1Tin−1∏v=iΨv+Tn
FS depends on the unknowns λ (IM:nrmShrFor) and G (IM:intsliceFs).
Refname | IM:nrmShrFor |
---|---|
Label | Normal and shear force proportionality constant |
Input | xslope, yslope, ywt, γw, xslip, yslip, constf |
Output | λ |
Input Constraints | |
Output Constraints | |
Equation | λ=n∑i=1Cnum,in∑i=1Cden,i |
Description |
|
Notes |
|
Source | chen2005 and karchewski2012 |
RefBy | IM:nrmShrForNum, IM:nrmShrFor, IM:nrmShrForDen, IM:intsliceFs, IM:fctSfty, FR:Display-Interslice-Shear-Forces, FR:Display-Interslice-Normal-Forces, FR:Display-Factor-of-Safety, and FR:Determine-Critical-Slip-Surface |
From the moment equilibrium of GD:momentEql with the primary assumption for the Morgenstern-Price method of A:Interslice-Norm-Shear-Forces-Linear and associated definition GD:normShrR, Equation (14) can be derived:
0=−Gi(hz,i+bi2tan(αi))+Gi−1(hz,i−1−bi2tan(αi))−Hi(13hz,w,i+bi2tan(αi))+Hi−1(13hz,w,i−1−bi2tan(αi))+λbi2(Gifi+Gi−1fi−1)+−KcWihi2+Ug,isin(βi)hi+Qisin(ωi)hi
Rearranging the equation in terms of λ leads to Equation (15):
λ=−Gi(hz,i+bi2tan(αi))+Gi−1(hz,i−1−bi2tan(αi))−Hi(13hz,w,i+bi2tan(αi))+Hi−1(13hz,w,i−1−bi2tan(αi))+−KcWihi2+Ug,isin(βi)hi+Qisin(ωi)hi−bi2(Gifi+Gi−1fi−1)
This equation can be simplified by applying assumptions A:Seismic-Force and A:Surface-Load, which state that the seismic and external forces, respectively, are zero:
λ=−Gi(hz,i+bi2tan(αi))+Gi−1(hz,i−1−bi2tan(αi))−Hi(13hz,w,i+bi2tan(αi))+Hi−1(13hz,w,i−1−bi2tan(αi))+Ug,isin(βi)hi−bi2(Gifi+Gi−1fi−1)
Taking the summation of all slices, and applying A:Edge-Slices to set G0, Gn, H0, and Hn equal to zero, a general equation for the proportionality constant λ is developed in Equation (16), which combines IM:nrmShrFor, IM:nrmShrForNum, and IM:nrmShrForDen:
λ=n∑i=1bi(FxG+FxH)tan(αi)+hi−2Ug,isin(βi)n∑i=1bi(Gifi+Gi−1fi−1)
Equation (16) for λ is a function of the unknown interslice normal forces G (IM:intsliceFs) which itself depends on the unknown factor of safety FS (IM:fctSfty).
Refname | IM:nrmShrForNum |
---|---|
Label | Normal and shear force proportionality constant numerator |
Input | xslope, yslope, ywt, γw, xslip, yslip |
Output | Cnum |
Input Constraints | |
Output Constraints | |
Equation | Cnum,i={b1(G1+H1)tan(α1),i=1bi(FxG+FxH)tan(αi)+h−2Ug,isin(βi),2≤i≤n−1bn(Gn−1+Hn−1)tan(αn−1),i=n |
Description |
|
Notes |
|
Source | chen2005 and karchewski2012 |
RefBy | IM:nrmShrFor |
See IM:nrmShrFor for the derivation of Cnum.
Refname | IM:nrmShrForDen |
---|---|
Label | Normal and shear force proportionality constant denominator |
Input | xslip, constf |
Output | Cden |
Input Constraints | |
Output Constraints | |
Equation | Cden,i={b1f1G1,i=1bi(fiGi+fi−1Gi−1),2≤i≤n−1bnGn−1fn−1,i=n |
Description |
|
Notes |
|
Source | chen2005 and karchewski2012 |
RefBy | IM:nrmShrFor |
See IM:nrmShrFor for the derivation of Cden.
Refname | IM:intsliceFs |
---|---|
Label | Interslice normal forces |
Input | xslope, yslope, ywt, c′, φ′, γdry, γsat, γw, xslip, yslip, constf |
Output | G |
Input Constraints | |
Output Constraints | |
Equation | Gi={FST1−R1Φ1,i=1Ψi−1Gi−1+FSTi−RiΦi,2≤i≤n−10,i=0∨i=n |
Description |
|
Notes |
|
Source | chen2005 |
RefBy | IM:nrmShrFor, IM:intsliceFs, IM:fctSfty, FR:Display-Interslice-Shear-Forces, FR:Display-Interslice-Normal-Forces, FR:Display-Factor-of-Safety, and FR:Determine-Critical-Slip-Surface |
This derivation is identical to the derivation for IM:fctSfty up until Equation (3) shown again below:
GiΦi=Ψi−1Gi−1Φi−1+FSTi−Ri
A simple rearrangement of Equation (3) leads to Equation (17), also seen in IM:intsliceFs:
Gi=Ψi−1Gi−1+FSTi−RiΦi
The cases shown in IM:intsliceFs for when i=0, i=1, or i=n are derived by applying A:Edge-Slices, which says that G0 and Gn are zero, to Equation (17). G depends on the unknowns FS (IM:fctSfty) and λ (IM:nrmShrFor).
Refname | IM:crtSlpId |
---|---|
Label | Critical slip surface identification |
Input | xslope, yslope, xwt, ywt, c′, φ′, γdry, γsat, γw, constf |
Output | FSmin |
Input Constraints | |
Output Constraints | |
Equation | xslope=yslope |
Description |
|
Notes |
|
Source | li2010 |
RefBy | FR:Display-Graph and FR:Determine-Critical-Slip-Surface |