3 Cylindrical Hole in an Infinite Mohr-Coulomb Medium
|
|
- Meagan King
- 6 years ago
- Views:
Transcription
1 Cylindrical Hole in an Infinite Mohr-Coulomb Medium Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3.1 Problem Statement The problem concerns the determination of stresses and displacements for the case of a cylindrical hole in an infinite elasto-plastic medium subjected to in-situ stresses. The medium is assumed to be linearly elastic, perfectly plastic, with a failure surface defined by the Mohr-Coulomb criterion with both the associated (dilatancy = friction angle) and non-associated (dilatancy = 0) flow rules. This problem tests the Mohr-Coulomb plasticity model, the plane-strain condition and axisymmetric geometry in FLAC. The Mohr-Coulomb material is assigned the following properties: density (ρ) 2500 kg/m 3 shear modulus (G) 2.8 GPa bulk modulus (K) 3.9 GPa cohesion (c) 3.45 MPa friction angle (φ) 30 dilation angle (ψ) 0 and 30 An isotropic in-situ stress state exists with stresses equal to 30 MPa (tension positive). It is assumed that the problem is symmetric about both the horizontal and vertical axes. The radius of the hole is 1 m, and is assumed to be small compared to the length of the cylinder. This permits the use of the plane-strain condition.
2 3-2 Verification Problems 3.2 Closed-Form Solution The yield zone radius, R o, is given analytically by a theoretical model based on the solution of Salencon (1969): R o [ 2 = a K p + 1 P o + P i + q ] 1/(Kp 1) K p 1 q K p 1 (3.1) where: a = radius of hole K p = 1+sin φ 1 sin φ q =2 c tan (45 + φ/2) P o = initial in-situ stress magnitude and P i = internal pressure. The radial stress at the elastic/plastic interface is σ re = 1 K p + 1 (2P o q) (3.2) The stresses in the plastic zone are: σ r = σ θ = q K p 1 (P i + q K p 1 K p(p i + q K p 1 ) ( r a )K p 1 q K p 1 ) ( r a )K p 1 (3.3) where r = distance to the center of the hole. The stresses in the elastic zone are: σ r = P o + (P o σ re ) ( R o r )2 (3.4) σ θ = P o (P o σ re ) ( R o r )2
3 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-3 The displacements in the elastic and plastic regions are given by Salencon (1969). For the elastic region, ( u r = P o ( 2P o q K p + 1 ) ) ( R o 2G )(R o ) r (3.5) and for the plastic region: u r = r 2G χ χ = (2ν 1) ( P o + ) q K p 1 ( ) (1 ν)(k 2 ( )( ) p 1) q (Kp 1) ( Ro Ro + P i + K p + K ps K p 1 a r ( + (1 ν) (K )( ) pk ps + 1) (K p + K ps ) ν q ( r ) (Kp 1) P i + K p 1 a ) (Kps +1) (3.6) where: K ps = 1+sin ψ 1 sin ψ ψ = dilation angle ν = Poisson s ratio and G = shear modulus.
4 3-4 Verification Problems 3.3 FLAC Model The problem is first modeled as a two-dimensional plane-strain calculation using quarter-symmetry. The boundary conditions applied to the model are shown in Figure 3.1. The area representing the problem is divided into finite-difference zones, as shown in Figures 3.2 and 3.3, using the FISH data file HOLE.FIS (see Section 3 in the FISH volume). The outer boundary is located 10 m (five hole diameters) from the hole center. The model contains 900 rectangular zones oriented in a radial pattern, as indicated in Figures 3.2 and 3.3. This pattern minimizes the influence of the grid on localization effects (as discussed in Section in the User s Guide). P o P o Figure 3.1 Model for FLAC analysis quarter-symmetry
5 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-5 JOB TITLE :. (*10^1) May-08 9:46 step E+00 <x< 1.167E E+00 <y< 1.167E Grid plot 0 2E (*10^1) Figure 3.2 FLAC grid quarter-symmetry JOB TITLE :. 16-May-08 9:46 step E-01 <x< 2.000E E-01 <y< 2.000E+00 Grid plot E Figure 3.3 Zone geometry in region around hole
6 3-6 Verification Problems The problem is also modeled using axisymmetric geometry. Figure 3.4 shows the boundary conditions, and Figure 3.5 shows the zone geometry created for this calculation. A height of two zones is used to facilitate stress contour plots. Both models are subjected to an isotropic compressive stress of 30 MPa. The initial stress state is applied throughout each model first then the hole is removed. The data file MHOLE.DAT in Section 3.6 contains the FLAC commands for both the quartersymmetry grid model and the axisymmetric model. The comparison of FLAC results to the analytical solution is performed using FISH functions also contained in this data file. The theoretical relations are solved in FISH function theor. Then, the FLAC values for stress and displacement are compared to the solution values. The error distribution in radial and tangential stress is calculated in the function evals, and the error distribution in displacement is calculated in the function evald. Axis of Symmetry P o Figure 3.4 Model for FLAC analysis axisymmetric geometry
7 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-7 JOB TITLE : May-08 9:47 step E-01 <x< 1.050E E+00 <y< 5.500E Grid plot 0 2E (*10^1) Figure 3.5 FLAC grid axisymmetric geometry 3.4 Results and Discussion Figures 3.6, 3.7 and 3.8 show a direct comparison between FLAC results and the analytical solution along a radial line for the quarter-symmetry case. Normalized stresses, σ r /P o and σ θ /P o, are plotted versus normalized radius, r/a, infigure 3.6 normalized displacement, u r /a, is plotted versus normalized radius in Figure 3.7 for the associated flow case, and in Figure 3.8 for the nonassociated flow case. (Note that the stress states are identical for both plastic flow rules.) The error distribution throughout the entire FLAC grid is also presented. Figures 3.9 through 3.11 show contour plots of tangential stress error, radial stress error and radial displacement error, respectively, for the non-associated flow case.
8 3-8 Verification Problems JOB TITLE :. 16-May-08 9:50 step 3262 Table Plot Table 10 Table 11 Table 20 Table Figure 3.6 Stress solution comparison tangential stress: analytical (Table 10) vs FLAC (Table 11) radial stress: analytical (Table 20) vs FLAC (Table 21) JOB TITLE :. 16-May-08 9:50 step 3262 Table Plot Table 30 Table (10 ) Figure 3.7 Radial displacement solution comparison associated flow: analytical (Table 30) vs FLAC (Table 31)
9 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-9 JOB TITLE :. -02 (10 ) 16-May-08 9:51 step Table Plot Table 30 Table Figure 3.8 Radial displacement solution comparison non-associated flow: analytical (Table 30) vs FLAC (Table 31) JOB TITLE :. (*10^1) May-08 10:59 step E+00 <x< 1.167E E+00 <y< 1.167E Boundary plot 0 2E 0 EX_ 6 Contours -3.50E E E E E E E E E E Contour interval= 5.00E-01 Extrap. by averaging (*10^1) Figure 3.9 Error distributions for tangential stress
10 3-10 Verification Problems JOB TITLE :. (*10^1) May-08 11:01 step E+00 <x< 1.167E E+00 <y< 1.167E Boundary plot 0 2E 0 EX_ 7 Contours -2.00E E E E E E E E E Contour interval= 2.50E-01 Extrap. by averaging (*10^1) Figure 3.10 Error distributions for radial stress JOB TITLE :. (*10^1) May-08 9:53 step E+00 <x< 1.167E E+00 <y< 1.167E Boundary plot 0 2E 0 EX_ 8 Contours -5.00E E E E E E E E E Contour interval= 5.00E (*10^1) Figure 3.11 Error distributions for displacements
11 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-11 For the axisymmetric geometry, Figures 3.12, 3.13 and 3.14 show a direct comparison between FLAC results and the analytical solution. Figures 3.15, 3.16 and 3.17 show the error distributions for tangential stress, radial stress and radial displacement. JOB TITLE : Hole in Mohr-Coulomb Medium, Axisymmetric Case 16-May-08 9:56 step 1902 Table Plot Table 10 Table 11 Table 20 Table Figure 3.12 Stress solution comparison axisymmetric case tangential stress: analytical (Table 10) vs FLAC (Table 11) radial stress: analytical (Table 20) vs FLAC (Table 21)
12 3-12 Verification Problems JOB TITLE : Hole in Mohr-Coulomb Medium, Axisymmetric Case 16-May-08 9:56 step 1902 Table Plot Table 30 Table (10 ) Figure 3.13 Radial displacement solution comparison axisymmetric case (associated flow): analytical (Table 30) vs FLAC (Table 31) JOB TITLE : Hole in Mohr-Coulomb Medium, Axisymmetric Case 16-May-08 9:57 step 2127 Table Plot Table 30 Table (10 ) Figure 3.14 Radial displacement solution comparison axisymmetric case (non-associated flow): analytical (Table 30) vs FLAC (Table 31)
13 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-13 JOB TITLE : Hole in Mohr-Coulomb medium, Axisymmetric Case May-08 11:04 step E-01 <x< 1.050E E+00 <y< 5.500E Boundary plot 0 2E 0 EX_ 6 Contours -3.00E E E E E E E E E+00 Contour interval= 5.00E-01 Extrap. by averaging (*10^1) Figure 3.15 Error distributions in tangential stress axisymmetric case JOB TITLE : Hole in Mohr-Coulomb medium, Axisymmetric Case May-08 11:04 step E-01 <x< 1.050E E+00 <y< 5.500E Boundary plot 0 2E 0 EX_ 7 Contours -1.75E E E E E E E E+00 Contour interval= 2.50E-01 Extrap. by averaging (*10^1) Figure 3.16 Error distributions in radial stress axisymmetric case
14 3-14 Verification Problems JOB TITLE : Hole in Mohr-Coulomb medium, Axisymmetric Case May-08 11:05 step E-01 <x< 1.050E E+00 <y< 5.500E Boundary plot 0 2E 0 EX_ 8 Contours 2.50E E E E E E E E E-01 Contour interval= 2.50E (*10^1) Figure 3.17 Error distribution in radial displacements axisymmetric case In both the quarter-symmetry model and the axisymmetric model, the errors are small the average stress error is less than 2%, while the average displacement error is less than 3%. Additional plots indicating the accuracy of the FLAC results are presented in Figures 3.18 through 3.20.
15 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-15 JOB TITLE :. (*10^1) May-08 10:05 step E+00 <x< 1.167E E+00 <y< 1.167E Plasticity Indicator * at yield in shear or vol. X elastic, at yield in past Boundary plot E (*10^1) Figure 3.18 Plasticity state indicators quarter-symmetry JOB TITLE :. 07 (10 ) 16-May-08 10:04 step E+00 <x< 1.167E E+00 <y< 1.167E Linear Profile Y-axis : YY-stress X-axis : Distance From ( 0.00E+00, 0.00E+00) To ( 1.00E+01, 0.00E+00) Figure 3.19 Linear profile of tangential stress quarter-symmetry
16 3-16 Verification Problems JOB TITLE :. 16-May-08 10:05 step 3262 Failure Surface Plot Major Prin. Stress vs. Minor Prin. Stress Zone Stress States Mohr-Coulomb Fail. Surf. Friction = E+01 Cohesion = E+06 Tension = E (10 ) (10 ) Figure 3.20 Zone stress states quarter-symmetry 3.5 Reference Salencon, J. Contraction Quasi-Statique D une Cavité a Symétrie Sphérique Ou Cylindrique Dans Un Milieu Elastoplastique, Annales Des Ponts Et Chaussées, 4, (1969).
17 Cylindrical Hole in an Infinite Mohr-Coulomb Medium Data File MHOLE.DAT Project Record Tree export *** Branch: Plane strain **** new... State: m3a.sav... title HOLE IN MOHR-COULOMB MEDIUM (ASSOCIATED FLOW) config extra=8 g mo mo call hole.fis set rmin=1 rmul=10 gratio=1.1 hole prop shear=2.8e9 bulk=3.9e9 dens=2500 coh=3.45e6 fric=30 ten=1e10 prop dil 0.0 non-associated flow prop dil 30.0 associated flow ini sxx=-30e6 syy=-30e6 szz=-30e6 fix y j 1 fix x j 31 app sxx=-30e6 syy=-30e6 i 31 hist unbal hist xd i 1 j 1 hist xv i 1 j 1 hist sxx i 1 j 1 hist syy i 1 j 1 hist szz i 1 j 1 solve save m3a.sav... State: m3a v.sav... rest m3aa.sav ****************** define the constants ********************* def parm p0=30e6 p1=0.0 rmin=1.0 s=cohesion(1,1) fi=friction(1,1)*degrad dil=dilation(1,1)*degrad bm=bulk mod(1,1) sm=shear mod(1,1) nu=(3.0*bm-2.0*sm)/(6.0*bm+2.0*sm)
18 3-18 Verification Problems kp=(1.0+sin(fi))/(1.0-sin(fi)) kps=(1.0+sin(dil))/(1.0-sin(dil)) q=2*s*cos(fi)/(1.0-sin(fi)) rp=rmin*(2.0/(kp+1.0)*(p0+q/(kp-1.0))/(p1+q/(kp-1.0)))ˆ(1.0/(kp-1)) sre=(2.0*p0-q)/(kp+1.0) parm *********** calculate the theoretical results ************** the theoretical results are stored in the following arrays tangential stress... EX 1 radial stress... EX 2 x displacements... EX 3 y displacements... EX 4 displacements magnitude... EX 5 def theor loop i (1,izones) loop j (1,jzones) xc=.25*(x(i,j)+x(i,j+1)+x(i+1,j+1)+x(i+1,j)) yc=.25*(y(i,j)+y(i,j+1)+y(i+1,j+1)+y(i+1,j)) rz=sqrt(xcˆ2+ycˆ2) if rz<=rp then ex 1(i,j)=q/(kp-1)-kp*(p1+q/(kp-1.0))*(rz/rmin)ˆ(kp-1.0) ex 2(i,j)=q/(kp-1)-(p1+q/(kp-1.0))*(rz/rmin)ˆ(kp-1.0) else if rz=0.0 then ex 1(i,j) = 0 ex 2(i,j) = 0 else ex 1(i,j)=-p0-(p0-sre)*(rp/rz)ˆ2 ex 2(i,j)=-p0+(p0-sre)*(rp/rz)ˆ2 if if loop i (1,igp) loop j (1,jgp) ro=sqrt(x(i,j)ˆ2+y(i,j)ˆ2) if ro#0 then if ro<=rp then d1=(2.0*nu-1.0)*(p0+q/(kp-1.0)) d2a=(1.0-nu)*(kpˆ2-1.0)/(kp+kps) d2b=(p1+q/(kp-1.0))*(rp/rmin)ˆ(kp-1.0)*(rp/ro)ˆ(kps+1.0) d3a=(1.0-nu)*(kp*kps+1.0)/(kp+kps)-nu
19 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-19 d3b=(p1+q/(kp-1.0))*(ro/rmin)ˆ(kp-1.0) dd=ro*(d1+d2a*d2b+d3a*d3b)/(2.0*sm) else dd=(p0-sre)*rp/(2.0*sm)*(rp/ro) if ex 3(i,j)=-dd*x(i,j)/ro ex 4(i,j)=-dd*y(i,j)/ro ex 5(i,j)=-dd if theor ************** evaluate the error in stresses ********************** the errors in stress calculations are evaluated for each zone and are stored in the following arrays: the total average errors are calculated and stored in: average average def evals ert=0 err=0 loop i (1,izones) loop j (1,jzones) temp1=.5*(sxx(i,j)+syy(i,j)) temp2=sqrt(sxy(i,j)ˆ2+.25*(sxx(i,j)-syy(i,j))ˆ2) stm=temp1-temp2 ex 6(i,j)=100*(stm-ex 1(i,j))/p0 ert=ert+ex 6(i,j)ˆ2 srm=temp1+temp2 ex 7(i,j)=100*(srm-ex 2(i,j))/p0 err=err+ex 7(i,j)ˆ2 ert=sqrt(ert/(izones*jzones)) err=sqrt(err/(izones*jzones)) evals ************** evaluate the error in displacements **************** the errors in displacement calculations are evaluated for each zone and are stored in the following array: the total average error is calculated and stored in: average
20 3-20 Verification Problems def evald erd=0 loop i (1,igp) loop j (1,jgp) temp3=100*sqrt((xdisp(i,j)-ex 3(i,j))ˆ2+(ydisp(i,j)-ex 4(i,j))ˆ2) if ex 5(1,1)#0 then ex 8(i,j)=temp3/ex 5(1,1) else ex 8(i,j)=0.0 if erd=erd+ex 8(i,j) erd=erd/(igp*jgp) evald ************** compare results in tables **************** def tabm The theoretical and numerical results are stored in terms of normalized radius in the following tables: table 10: tangential stress theoretical table 11: tangential stress numerical table 20: radial stress theoretical table 21: radial stress numerical table 30: radial displacement theoretical table 31: radial displacement numerical loop i(1,izones) j = 1 xc=.25*(x(i,j)+x(i,j+1)+x(i+1,j+1)+x(i+1,j)) yc=.25*(y(i,j)+y(i,j+1)+y(i+1,j+1)+y(i+1,j)) rz=sqrt(xcˆ2+ycˆ2) temp1=.5*(sxx(i,j)+syy(i,j)) temp2=sqrt(sxy(i,j)ˆ2+.25*(sxx(i,j)-syy(i,j))ˆ2) stm=temp1-temp2 srm=temp1+temp2 xtable(10,i) = rz ytable(10,i) = -ex 1(i,1)/p0 xtable(11,i) = rz ytable(11,i) = -stm/p0 xtable(20,i) = rz ytable(20,i) = -ex 2(i,1)/p0 xtable(21,i) = rz ytable(21,i) = -srm/p0
21 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-21 loop i (1,igp) j = 1 rg = sqrt(x(i,j)ˆ2+y(i,j)ˆ2)/rmin temp1=sqrt(xdisp(i,j)ˆ2+ydisp(i,j)ˆ2)/rmin xtable(30,i) = rg ytable(30,i) = -ex 5(i,j) / rmin xtable(31,i) = rg ytable(31,i) = temp1 tabm ****************** create plots ******************************* scline scline save m3a v.sav *** Branch: Plane strain-non **** new... State: m3a-non.sav... title HOLE IN MOHR-COULOMB MEDIUM (ASSOCIATED FLOW) config extra=8 g mo mo call hole.fis set rmin=1 rmul=10 gratio=1.1 hole prop shear=2.8e9 bulk=3.9e9 dens=2500 coh=3.45e6 fric=30 ten=1e10 prop dil 0.0 non-associated flow prop dil 30.0 associated flow ini sxx=-30e6 syy=-30e6 szz=-30e6 fix y j 1 fix x j 31 app sxx=-30e6 syy=-30e6 i 31 hist unbal hist xd i 1 j 1 hist xv i 1 j 1 hist sxx i 1 j 1 hist syy i 1 j 1 hist szz i 1 j 1 solve save m3a-non.sav... State: m3a-non v.sav...
22 3-22 Verification Problems rest m3aa.sav ****************** define the constants ********************* def parm p0=30e6 p1=0.0 rmin=1.0 s=cohesion(1,1) fi=friction(1,1)*degrad dil=dilation(1,1)*degrad bm=bulk mod(1,1) sm=shear mod(1,1) nu=(3.0*bm-2.0*sm)/(6.0*bm+2.0*sm) kp=(1.0+sin(fi))/(1.0-sin(fi)) kps=(1.0+sin(dil))/(1.0-sin(dil)) q=2*s*cos(fi)/(1.0-sin(fi)) rp=rmin*(2.0/(kp+1.0)*(p0+q/(kp-1.0))/(p1+q/(kp-1.0)))ˆ(1.0/(kp-1)) sre=(2.0*p0-q)/(kp+1.0) parm *********** calculate the theoretical results ************** the theoretical results are stored in the following arrays tangential stress... EX 1 radial stress... EX 2 x displacements... EX 3 y displacements... EX 4 displacements magnitude... EX 5 def theor loop i (1,izones) loop j (1,jzones) xc=.25*(x(i,j)+x(i,j+1)+x(i+1,j+1)+x(i+1,j)) yc=.25*(y(i,j)+y(i,j+1)+y(i+1,j+1)+y(i+1,j)) rz=sqrt(xcˆ2+ycˆ2) if rz<=rp then ex 1(i,j)=q/(kp-1)-kp*(p1+q/(kp-1.0))*(rz/rmin)ˆ(kp-1.0) ex 2(i,j)=q/(kp-1)-(p1+q/(kp-1.0))*(rz/rmin)ˆ(kp-1.0) else if rz=0.0 then ex 1(i,j) = 0 ex 2(i,j) = 0 else ex 1(i,j)=-p0-(p0-sre)*(rp/rz)ˆ2 ex 2(i,j)=-p0+(p0-sre)*(rp/rz)ˆ2 if
23 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-23 if loop i (1,igp) loop j (1,jgp) ro=sqrt(x(i,j)ˆ2+y(i,j)ˆ2) if ro#0 then if ro<=rp then d1=(2.0*nu-1.0)*(p0+q/(kp-1.0)) d2a=(1.0-nu)*(kpˆ2-1.0)/(kp+kps) d2b=(p1+q/(kp-1.0))*(rp/rmin)ˆ(kp-1.0)*(rp/ro)ˆ(kps+1.0) d3a=(1.0-nu)*(kp*kps+1.0)/(kp+kps)-nu d3b=(p1+q/(kp-1.0))*(ro/rmin)ˆ(kp-1.0) dd=ro*(d1+d2a*d2b+d3a*d3b)/(2.0*sm) else dd=(p0-sre)*rp/(2.0*sm)*(rp/ro) if ex 3(i,j)=-dd*x(i,j)/ro ex 4(i,j)=-dd*y(i,j)/ro ex 5(i,j)=-dd if theor ************** evaluate the error in stresses ********************** the errors in stress calculations are evaluated for each zone and are stored in the following arrays: the total average errors are calculated and stored in: average average def evals ert=0 err=0 loop i (1,izones) loop j (1,jzones) temp1=.5*(sxx(i,j)+syy(i,j)) temp2=sqrt(sxy(i,j)ˆ2+.25*(sxx(i,j)-syy(i,j))ˆ2) stm=temp1-temp2 ex 6(i,j)=100*(stm-ex 1(i,j))/p0 ert=ert+ex 6(i,j)ˆ2 srm=temp1+temp2 ex 7(i,j)=100*(srm-ex 2(i,j))/p0 err=err+ex 7(i,j)ˆ2
24 3-24 Verification Problems ert=sqrt(ert/(izones*jzones)) err=sqrt(err/(izones*jzones)) evals ************** evaluate the error in displacements **************** the errors in displacement calculations are evaluated for each zone and are stored in the following array: the total average error is calculated and stored in: average def evald erd=0 loop i (1,igp) loop j (1,jgp) temp3=100*sqrt((xdisp(i,j)-ex 3(i,j))ˆ2+(ydisp(i,j)-ex 4(i,j))ˆ2) if ex 5(1,1)#0 then ex 8(i,j)=temp3/ex 5(1,1) else ex 8(i,j)=0.0 if erd=erd+ex 8(i,j) erd=erd/(igp*jgp) evald ************** compare results in tables **************** def tabm The theoretical and numerical results are stored in terms of normalized radius in the following tables: table 10: tangential stress theoretical table 11: tangential stress numerical table 20: radial stress theoretical table 21: radial stress numerical table 30: radial displacement theoretical table 31: radial displacement numerical loop i(1,izones) j = 1 xc=.25*(x(i,j)+x(i,j+1)+x(i+1,j+1)+x(i+1,j)) yc=.25*(y(i,j)+y(i,j+1)+y(i+1,j+1)+y(i+1,j)) rz=sqrt(xcˆ2+ycˆ2) temp1=.5*(sxx(i,j)+syy(i,j))
25 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-25 temp2=sqrt(sxy(i,j)ˆ2+.25*(sxx(i,j)-syy(i,j))ˆ2) stm=temp1-temp2 srm=temp1+temp2 xtable(10,i) = rz ytable(10,i) = -ex 1(i,1)/p0 xtable(11,i) = rz ytable(11,i) = -stm/p0 xtable(20,i) = rz ytable(20,i) = -ex 2(i,1)/p0 xtable(21,i) = rz ytable(21,i) = -srm/p0 loop i (1,igp) j = 1 rg = sqrt(x(i,j)ˆ2+y(i,j)ˆ2)/rmin temp1=sqrt(xdisp(i,j)ˆ2+ydisp(i,j)ˆ2)/rmin xtable(30,i) = rg ytable(30,i) = -ex 5(i,j) / rmin xtable(31,i) = rg ytable(31,i) = temp1 tabm ****************** create plots ******************************* scline scline save m3a-non v.sav *** Branch: axisymmetric **** new... State: m3b.sav... title HOLE IN MOHR-COULOMB MEDIUM (AXISYMMETRY, ASSOCIATED FLOW RULE) config ax extra=8 g 31 2 mo mo gen rat prop shear=2.8e9 bulk=3.9e9 dens=2500 coh=3.45e6 fric=30 ten=6e6 prop dil 0.0 non-associated flow prop dil 30.0 associated flow ini sxx=-30e6 syy=-30e6 szz=-30e6 fix y j 1 fix y j 3 app press 30e6 i 32 hist unbal
26 3-26 Verification Problems hist xd i 1 j 1 hist xv i 1 j 1 hist sxx i 1 j 1 hist syy i 1 j 1 hist sig1 i 1 j 1 hist sig2 i 1 j 1 hist szz i 1 j 1 solve save m3b.sav... State: m3b v.sav... rest m3bb.sav ****************** define the constants ********************* def parm p0=30e6 p1=0.0 rmin=1.0 s=cohesion(1,1) fi=friction(1,1)*degrad dil=dilation(1,1)*degrad bm=bulk mod(1,1) sm=shear mod(1,1) nu=(3.0*bm-2.0*sm)/(6.0*bm+2.0*sm) kp=(1.0+sin(fi))/(1.0-sin(fi)) kps=(1.0+sin(dil))/(1.0-sin(dil)) q=2*s*cos(fi)/(1.0-sin(fi)) rp=rmin*(2.0/(kp+1.0)*(p0+q/(kp-1.0))/(p1+q/(kp-1.0)))ˆ(1.0/(kp-1)) sre=(2.0*p0-q)/(kp+1.0) parm *********** calculate the theoretical results ************** the theoretical results are stored in the following arrays tangential stress... EX 1 radial stress... EX 2 x displacements... EX 3 y displacements... EX 4 displacements magnitude... EX 5 def theor loop i (1,izones) loop j (1,jzones) rz=.25*(x(i,j)+x(i,j+1)+x(i+1,j+1)+x(i+1,j)) if rz<=rp then ex 1(i,j)=q/(kp-1)-kp*(p1+q/(kp-1.0))*(rz/rmin)ˆ(kp-1.0)
27 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-27 ex 2(i,j)=q/(kp-1)-(p1+q/(kp-1.0))*(rz/rmin)ˆ(kp-1.0) else ex 1(i,j)=-p0-(p0-sre)*(rp/rz)ˆ2 ex 2(i,j)=-p0+(p0-sre)*(rp/rz)ˆ2 if loop i (1,igp) loop j (1,jgp) ro=x(i,j) if ro<=rp then d1=(2.0*nu-1.0)*(p0+q/(kp-1.0)) d2a=(1.0-nu)*(kpˆ2-1.0)/(kp+kps) d2b=(p1+q/(kp-1.0))*(rp/rmin)ˆ(kp-1.0)*(rp/ro)ˆ(kps+1.0) d3a=(1.0-nu)*(kp*kps+1.0)/(kp+kps)-nu d3b=(p1+q/(kp-1.0))*(ro/rmin)ˆ(kp-1.0) dd=ro*(d1+d2a*d2b+d3a*d3b)/(2.0*sm) else dd=(p0-sre)*rp/(2.0*sm)*(rp/ro) if ex 3(i,j)=-dd ex 4(i,j)=0.0 ex 5(i,j)=-dd theor ************** evaluate the error in stresses ********************** the errors in stress calculations are evaluated for each zone and are stored in the following arrays: the total average errors are calculated and stored in: average average def evals ert=0 err=0 loop i (1,izones) loop j (1,jzones) ex 6(i,j)=100*(szz(i,j)-ex 1(i,j))/p0 ert=ert+ex 6(i,j)ˆ2 ex 7(i,j)=100*(sxx(i,j)-ex 2(i,j))/p0 err=err+ex 7(i,j)ˆ2 ert=sqrt(ert/(izones*jzones))
28 3-28 Verification Problems err=sqrt(err/(izones*jzones)) evals ************** evaluate the error in displacements **************** the errors in displacement calculations are evaluated for each zone and are stored in the following array: the total average error is calculated and stored in: average def evald erd=0 loop i (1,igp) loop j (1,jgp) temp3=100*sqrt((xdisp(i,j)-ex 3(i,j))ˆ2+(ydisp(i,j)-ex 4(i,j))ˆ2) if ex 8(i,j)#0.0 then ex 8(i,j)=temp3/ex 5(1,1) else ex 8(i,j)=temp3 if erd=erd+ex 8(i,j) erd=erd/(igp*jgp) evald ************** compare results in tables **************** def tabm The theoretical and numerical results are stored in terms of normalized radius in the following tables: table 10: tangential stress theoretical table 11: tangential stress numerical table 20: radial stress theoretical table 21: radial stress numerical table 30: radial displacement theoretical table 31: radial displacement numerical loop i(1,izones) j = 1 rz=.25*(x(i,j)+x(i,j+1)+x(i+1,j+1)+x(i+1,j)) xtable(10,i) = rz ytable(10,i) = -ex 1(i,1)/p0 xtable(11,i) = rz ytable(11,i) = -szz(i,j)/p0 xtable(20,i) = rz
29 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-29 ytable(20,i) = -ex 2(i,1)/p0 xtable(21,i) = rz ytable(21,i) = -sxx(i,j)/p0 loop i (1,igp) j = 1 rg = sqrt(x(i,j)ˆ2+y(i,j)ˆ2)/rmin temp1=sqrt(xdisp(i,j)ˆ2+ydisp(i,j)ˆ2)/rmin xtable(30,i) = rg ytable(30,i) = -ex 5(i,j) / rmin xtable(31,i) = rg ytable(31,i) = temp1 tabm ****************** create plots ******************************* scline save m3b v.sav *** Branch: axisymmetric-non **** new... State: m3b-non.sav... title HOLE IN MOHR-COULOMB MEDIUM (AXISYMMETRY, ASSOCIATED FLOW RULE) config ax extra=8 g 31 2 mo mo gen rat prop shear=2.8e9 bulk=3.9e9 dens=2500 coh=3.45e6 fric=30 ten=6e6 prop dil 0.0 non-associated flow prop dil 30.0 associated flow ini sxx=-30e6 syy=-30e6 szz=-30e6 fix y j 1 fix y j 3 app press 30e6 i 32 hist unbal hist xd i 1 j 1 hist xv i 1 j 1 hist sxx i 1 j 1 hist syy i 1 j 1 hist sig1 i 1 j 1 hist sig2 i 1 j 1 hist szz i 1 j 1 solve save m3b-non.sav
30 3-30 Verification Problems... State: m3b-non v.sav... rest m3bb.sav ****************** define the constants ********************* def parm p0=30e6 p1=0.0 rmin=1.0 s=cohesion(1,1) fi=friction(1,1)*degrad dil=dilation(1,1)*degrad bm=bulk mod(1,1) sm=shear mod(1,1) nu=(3.0*bm-2.0*sm)/(6.0*bm+2.0*sm) kp=(1.0+sin(fi))/(1.0-sin(fi)) kps=(1.0+sin(dil))/(1.0-sin(dil)) q=2*s*cos(fi)/(1.0-sin(fi)) rp=rmin*(2.0/(kp+1.0)*(p0+q/(kp-1.0))/(p1+q/(kp-1.0)))ˆ(1.0/(kp-1)) sre=(2.0*p0-q)/(kp+1.0) parm *********** calculate the theoretical results ************** the theoretical results are stored in the following arrays tangential stress... EX 1 radial stress... EX 2 x displacements... EX 3 y displacements... EX 4 displacements magnitude... EX 5 def theor loop i (1,izones) loop j (1,jzones) rz=.25*(x(i,j)+x(i,j+1)+x(i+1,j+1)+x(i+1,j)) if rz<=rp then ex 1(i,j)=q/(kp-1)-kp*(p1+q/(kp-1.0))*(rz/rmin)ˆ(kp-1.0) ex 2(i,j)=q/(kp-1)-(p1+q/(kp-1.0))*(rz/rmin)ˆ(kp-1.0) else ex 1(i,j)=-p0-(p0-sre)*(rp/rz)ˆ2 ex 2(i,j)=-p0+(p0-sre)*(rp/rz)ˆ2 if loop i (1,igp) loop j (1,jgp)
31 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-31 ro=x(i,j) if ro<=rp then d1=(2.0*nu-1.0)*(p0+q/(kp-1.0)) d2a=(1.0-nu)*(kpˆ2-1.0)/(kp+kps) d2b=(p1+q/(kp-1.0))*(rp/rmin)ˆ(kp-1.0)*(rp/ro)ˆ(kps+1.0) d3a=(1.0-nu)*(kp*kps+1.0)/(kp+kps)-nu d3b=(p1+q/(kp-1.0))*(ro/rmin)ˆ(kp-1.0) dd=ro*(d1+d2a*d2b+d3a*d3b)/(2.0*sm) else dd=(p0-sre)*rp/(2.0*sm)*(rp/ro) if ex 3(i,j)=-dd ex 4(i,j)=0.0 ex 5(i,j)=-dd theor ************** evaluate the error in stresses ********************** the errors in stress calculations are evaluated for each zone and are stored in the following arrays: the total average errors are calculated and stored in: average average def evals ert=0 err=0 loop i (1,izones) loop j (1,jzones) ex 6(i,j)=100*(szz(i,j)-ex 1(i,j))/p0 ert=ert+ex 6(i,j)ˆ2 ex 7(i,j)=100*(sxx(i,j)-ex 2(i,j))/p0 err=err+ex 7(i,j)ˆ2 ert=sqrt(ert/(izones*jzones)) err=sqrt(err/(izones*jzones)) evals ************** evaluate the error in displacements **************** the errors in displacement calculations are evaluated for each zone and are stored in the following array: the total average error is calculated and stored in:
32 3-32 Verification Problems average def evald erd=0 loop i (1,igp) loop j (1,jgp) temp3=100*sqrt((xdisp(i,j)-ex 3(i,j))ˆ2+(ydisp(i,j)-ex 4(i,j))ˆ2) if ex 8(i,j)#0.0 then ex 8(i,j)=temp3/ex 5(1,1) else ex 8(i,j)=temp3 if erd=erd+ex 8(i,j) erd=erd/(igp*jgp) evald ************** compare results in tables **************** def tabm The theoretical and numerical results are stored in terms of normalized radius in the following tables: table 10: tangential stress theoretical table 11: tangential stress numerical table 20: radial stress theoretical table 21: radial stress numerical table 30: radial displacement theoretical table 31: radial displacement numerical loop i(1,izones) j = 1 rz=.25*(x(i,j)+x(i,j+1)+x(i+1,j+1)+x(i+1,j)) xtable(10,i) = rz ytable(10,i) = -ex 1(i,1)/p0 xtable(11,i) = rz ytable(11,i) = -szz(i,j)/p0 xtable(20,i) = rz ytable(20,i) = -ex 2(i,1)/p0 xtable(21,i) = rz ytable(21,i) = -sxx(i,j)/p0 loop i (1,igp) j = 1 rg = sqrt(x(i,j)ˆ2+y(i,j)ˆ2)/rmin temp1=sqrt(xdisp(i,j)ˆ2+ydisp(i,j)ˆ2)/rmin xtable(30,i) = rg
33 Cylindrical Hole in an Infinite Mohr-Coulomb Medium 3-33 ytable(30,i) = -ex 5(i,j) / rmin xtable(31,i) = rg ytable(31,i) = temp1 tabm ****************** create plots ******************************* scline save m3b-non v.sav *** plot commands **** plot name: grid plot hold grid plot name: Stresses comparison plot hold table 21 cross 20 line 11 cross 10 line plot name: Rad. disp. comp. plot hold table 31 cross 30 line plot name: Error - tang. stress plot hold bound ex 6 zone fill plot name: Error - radial stress plot hold bound ex 7 zone fill plot name: Error - disp plot hold bound ex 8 fill plot name: Plasticity indicator plot hold plasticity bound plot name: Zone stress state set pltc pltf 30.0 pltt 0.0 plot hold fail principal plot name: Profile of syy plot hold syy line (0.0,0.0) (10.0,0.0) 31
34 3-34 Verification Problems
5 Spherical Cavity in an Infinite Elastic Medium
Spherical Cavity in an Infinite Elastic Medium 5-1 5 Spherical Cavity in an Infinite Elastic Medium 5.1 Problem Statement Stresses and displacements are determined for the case of a spherical cavity in
More informationExample-3. Title. Description. Cylindrical Hole in an Infinite Mohr-Coulomb Medium
Example-3 Title Cylindrical Hole in an Infinite Mohr-Coulomb Medium Description The problem concerns the determination of stresses and displacements for the case of a cylindrical hole in an infinite elasto-plastic
More information4 Cylindrical Hole in an Infinite Hoek-Brown Medium
Cylindrical Hole in an Infinite Hoek-Brown Medium 4-1 4 Cylindrical Hole in an Infinite Hoek-Brown Medium 4.1 Problem Statement Stresses and displacements are calculated for the case of a cylindrical hole
More information8 Lined Circular Tunnel in an Elastic Medium with Anisotropic Stresses
Lined Circular Tunnel in an Elastic Medium with Anisotropic Stresses 8-1 8 Lined Circular Tunnel in an Elastic Medium with Anisotropic Stresses 8.1 Problem Statement This problem concerns the analysis
More information8 Displacements near the Face of an Advancing Shaft
Displacements near the Face of an Advancing Shaft 8-1 8 Displacements near the Face of an Advancing Shaft 8.1 Problem Statement A circular shaft is excavated in chalk and lined with monolithic precast
More information7 Uniaxial Compressive Strength of a Jointed Rock Sample
Uniaxial Compressive Strength of a Jointed Rock Sample 7-1 7 Uniaxial Compressive Strength of a Jointed Rock Sample 7.1 Problem Statement The uniaxial compressive strength of a jointed rock sample is a
More information1 Slope Stability for a Cohesive and Frictional Soil
Slope Stability for a Cohesive and Frictional Soil 1-1 1 Slope Stability for a Cohesive and Frictional Soil 1.1 Problem Statement A common problem encountered in engineering soil mechanics is the stability
More informationCylindrical Hole in an Infinite Mohr-Coulomb Medium
Benchmark Example No. 29 Cylindrical Hole in an Infinite Mohr-Coulomb Medium SOFiSTiK 2018 VERiFiCATiON MANUAL BE29: Cylindrical Hole in an Infinite Mohr-Coulomb Medium VERiFiCATiON MANUAL, Version 2018-7
More information10 Slope Stability Analysis of a Rock Slope
Slope Stability Analysis of a Rock Slope 10-1 10 Slope Stability Analysis of a Rock Slope 10.1 Problem Statement Limit equilibrium methods are commonly used to evaluate the stability of slopes in rock
More informationA circular tunnel in a Mohr-Coulomb medium with an overlying fault
MAP3D VERIFICATION EXAMPLE 9 A circular tunnel in a Mohr-Coulomb medium with an overlying fault 1 Description This example involves calculating the stresses and displacements on a fault overlying a 5 m
More information13 Plastic Flow in a Punch Problem
Plastic Flow in a Punch Problem 13-1 13 Plastic Flow in a Punch Problem 13.1 Problem Statement Difficulties are sometimes reported in the modeling of plastic flow where large velocity gradients exist.
More information4 Undrained Cylindrical Cavity Expansion in a Cam-Clay Medium
Undrained Cylindrical Cavity Expansion in a Cam-Clay Medium 4-1 4 Undrained Cylindrical Cavity Expansion in a Cam-Clay Medium 4.1 Problem Statement The stress and pore pressure changes due to the expansion
More information1 Slope Stability for a Cohesive and Frictional Soil
Slope Stability for a Cohesive and Frictional Soil 1-1 1 Slope Stability for a Cohesive and Frictional Soil 1.1 Problem Statement A common problem encountered in engineering soil mechanics is the stability
More information15 Drained and Undrained Triaxial Compression Test on a Cam-Clay Sample
Drained and Undrained Triaxial Compression Test on a Cam-Clay Sample 15-1 15 Drained and Undrained Triaxial Compression Test on a Cam-Clay Sample 15.1 Problem Statement Conventional drained and undrained
More informationNumerical Modeling in Geotechnical Engineering October 17,
1 Numerical Modeling in Geotechnical Engineering October 17, 18 2014 PRESENTED BY: Siavash Zamiran Ph.D. Student, Research and Teaching Assistant, Department of Civil Engineering, Southern Illinois University
More informationFigure A7 22: Micrograph of PX (a) before and (b) after heat treatment to 50ºC
Figure A7 22: Micrograph of PX (a) before and (b) after heat treatment to 50ºC 293 Figure A7 23: Micrograph of PX after heat treatment to 100ºC and 140ºC 294 Figure A7 24: Micrograph of VTA (a) before
More informationEarth Pressure Theory
Lateral Earth Pressure Page 1 Earth Pressure Theory Examples of Retaining Walls Lateral Earth Pressure Page 2 At-Rest, Active and Passive Earth Pressure Wednesday, August 17, 2011 12:45 PM At-rest condition
More informationCVEEN Table of Contents
CVEEN 7330 Table of Contents Table of Contents... 1 Objectives:... 2 FLAC Input:... 3 Model Inputs... 3 Create Grid and Assign Properties... 4 Create Foundation... 4 Create Forcing Function... 4 Create
More information13 Step-Path Failure of Rock Slopes
Step-Path Failure of Rock Slopes 13-1 13 Step-Path Failure of Rock Slopes 13.1 Problem Statement While UDEC represents a jointed rock structure as a system of discrete blocks by default, it is also a straightforward
More information2 Open Stoping Using Vertical Retreat
Open Stoping Using Vertical Retreat 2-1 2 Open Stoping Using Vertical Retreat 2.1 Problem Statement A distinct element simulation of a large blasthole, open stoping operation is shown to demonstrate the
More information10 Steady-State Fluid Flow with a Free Surface
Steady-State Fluid Flow with a Free Surface 10-1 10 Steady-State Fluid Flow with a Free Surface 10.1 roblem Statement This numerical simulation analyzes the steady-state seepage flow through a homogeneous
More information2 CONSTITUTIVE MODELS: THEORY AND IMPLEMENTATION
CONSTITUTIVE MODELS: THEORY AND IMPLEMENTATION 2-1 2 CONSTITUTIVE MODELS: THEORY AND IMPLEMENTATION 2.1 Introduction There are twelve basic constitutive models provided in, arranged into null, elastic
More information16 Rainfall on a Slope
Rainfall on a Slope 16-1 16 Rainfall on a Slope 16.1 Problem Statement In this example, the stability of a generic slope is analyzed for two successive rainfall events of increasing intensity and decreasing
More informationPractice Final Examination. Please initial the statement below to show that you have read it
EN175: Advanced Mechanics of Solids Practice Final Examination School of Engineering Brown University NAME: General Instructions No collaboration of any kind is permitted on this examination. You may use
More information13 Dewatered Construction of a Braced Excavation
Dewatered Construction of a Braced Excavation 13-1 13 Dewatered Construction of a Braced Excavation 13.1 Problem Statement A braced excavation is constructed in saturated ground. The excavation is dewatered
More informationEffect of embedment depth and stress anisotropy on expansion and contraction of cylindrical cavities
Effect of embedment depth and stress anisotropy on expansion and contraction of cylindrical cavities Hany El Naggar, Ph.D., P. Eng. and M. Hesham El Naggar, Ph.D., P. Eng. Department of Civil Engineering
More information3 Thick-Walled Cylinder Subject to Internal Pressure
Thick-Walled Cylinder Subject to Internal Pressure 3-1 3 Thick-Walled Cylinder Subject to Internal Pressure 3.1 Problem Statement This problem involves the plane strain elastic-plastic analysis of a thick-walled
More informationVerification Manual GT
Verification Manual GT Written by: The SoilVision Systems Ltd. Team Last Updated: Tuesday, February 20, 2018 SoilVision Systems Ltd. Saskatoon, Saskatchewan, Canada Software License The software described
More informationCVEEN Table of Contents
CVEEN 7330 Table of Contents Table of Contents... 1 Objectives... 2 ackground... 2 Model Setup and Analysis Input... 2 Data for Interfaces... 5 Required Outputs from FLAC... 10 Required Calculations and
More informationu a 1 Cyclic Loading of a Specimen with a Slipping Crack Cyclic Loading of a Specimen with a Slipping Crack Problem Statement
Cyclic Loading of a Specimen with a Slipping Crack 1-1 1 Cyclic Loading of a Specimen with a Slipping Crack 1.1 Problem Statement This problem concerns an elastic block with an inclined internal closed
More informationModule 6: Stresses around underground openings. 6.2 STRESSES AROUND UNDERGROUND OPENING contd.
LECTURE 0 6. STRESSES AROUND UNDERGROUND OPENING contd. CASE : When σ x = 0 For σ x = 0, the maximum tangential stress is three times the applied stress and occurs at the boundary on the X-axis that is
More informationThermal Analysis Contents - 1
Thermal Analysis Contents - 1 TABLE OF CONTENTS 1 THERMAL ANALYSIS 1.1 Introduction... 1-1 1.2 Mathematical Model Description... 1-3 1.2.1 Conventions and Definitions... 1-3 1.2.2 Conduction... 1-4 1.2.2.1
More informationPlane and axisymmetric models in Mentat & MARC. Tutorial with some Background
Plane and axisymmetric models in Mentat & MARC Tutorial with some Background Eindhoven University of Technology Department of Mechanical Engineering Piet J.G. Schreurs Lambèrt C.A. van Breemen March 6,
More informationLab Practical - Discontinuum Analysis & Distinct Element Method
Lab Practical - Discontinuum Analysis & Distinct Element Method Part A The Basics The Universal Distinct Element Code (UDEC) is a two-dimensional numerical program based on the distinct element method
More information2D Embankment and Slope Analysis (Numerical)
2D Embankment and Slope Analysis (Numerical) Page 1 2D Embankment and Slope Analysis (Numerical) Sunday, August 14, 2011 Reading Assignment Lecture Notes Other Materials FLAC Manual 1. 2. Homework Assignment
More informationChapter 3. Load and Stress Analysis
Chapter 3 Load and Stress Analysis 2 Shear Force and Bending Moments in Beams Internal shear force V & bending moment M must ensure equilibrium Fig. 3 2 Sign Conventions for Bending and Shear Fig. 3 3
More informationSettlement and Bearing Capacity of a Strip Footing. Nonlinear Analyses
Settlement and Bearing Capacity of a Strip Footing Nonlinear Analyses Outline 1 Description 2 Nonlinear Drained Analysis 2.1 Overview 2.2 Properties 2.3 Loads 2.4 Analysis Commands 2.5 Results 3 Nonlinear
More informationRecent Research on EPS Geofoam Seismic Buffers. Richard J. Bathurst and Saman Zarnani GeoEngineering Centre at Queen s-rmc Canada
Recent Research on EPS Geofoam Seismic Buffers Richard J. Bathurst and Saman Zarnani GeoEngineering Centre at Queen s-rmc Canada What is a wall (SEISMIC) buffer? A compressible inclusion placed between
More informationStresses Analysis of Petroleum Pipe Finite Element under Internal Pressure
ISSN : 48-96, Vol. 6, Issue 8, ( Part -4 August 06, pp.3-38 RESEARCH ARTICLE Stresses Analysis of Petroleum Pipe Finite Element under Internal Pressure Dr.Ragbe.M.Abdusslam Eng. Khaled.S.Bagar ABSTRACT
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting Lectures & 3, 9/31 Aug 017 www.geosc.psu.edu/courses/geosc508 Discussion of Handin, JGR, 1969 and Chapter 1 Scholz, 00. Stress analysis and Mohr Circles Coulomb Failure
More informationLab Practical - Discontinuum Analysis & Distinct Element Method
Lab Practical - Discontinuum Analysis & Distinct Element Method Part A The Basics The Universal Distinct Element Code (UDEC) is a twodimensional numerical program based on the distinct element method for
More informationExercise: concepts from chapter 8
Reading: Fundamentals of Structural Geology, Ch 8 1) The following exercises explore elementary concepts associated with a linear elastic material that is isotropic and homogeneous with respect to elastic
More informationVerification of the Hyperbolic Soil Model by Triaxial Test Simulations
1 Introduction Verification of the Hyperbolic Soil Model by Triaxial Test Simulations This example simulates a series of triaxial tests that can be used to verify that the Hyperbolic constitutive model
More informationLateral Earth Pressure
1 of 11 6/2/2012 4:28 AM Lateral Earth Pressure The magnitude of lateral earth pressure depends on: 1. Shear strength characteristics of soil 2. Lateral strain condition 3. Pore water pressure 4. State
More informationUSER S MANUAL 1D Seismic Site Response Analysis Example University of California: San Diego August 30, 2017
USER S MANUAL 1D Seismic Site Response Analysis Example http://www.soilquake.net/ucsdsoilmodels/ University of California: San Diego August 30, 2017 Table of Contents USER'S MANUAL TABLE OF CONTENTS Page
More informationDynamic Analysis Contents - 1
Dynamic Analysis Contents - 1 TABLE OF CONTENTS 1 DYNAMIC ANALYSIS 1.1 Overview... 1-1 1.2 Relation to Equivalent-Linear Methods... 1-2 1.2.1 Characteristics of the Equivalent-Linear Method... 1-2 1.2.2
More informationEffect of the intermediate principal stress on fault strike and dip - theoretical analysis and experimental verification
Effect of the intermediate principal stress on fault strike and dip - theoretical analysis and experimental verification B. Haimson University of Wisconsin, USA J. Rudnicki Northwestern University, USA
More informationStress distribution in inclined backfilled stopes
MINEFILL 2007 Montreal, Qc, Canada Stress distribution in inclined backfilled stopes L. Li 1, M. Aubertin 1,3, A. Shirazi 1,T. Belem 2, R. Simon 1 1 École Polytechnique de Montréal 2 Université du Québec
More informationROTATING RING. Volume of small element = Rdθbt if weight density of ring = ρ weight of small element = ρrbtdθ. Figure 1 Rotating ring
ROTATIONAL STRESSES INTRODUCTION High centrifugal forces are developed in machine components rotating at a high angular speed of the order of 100 to 500 revolutions per second (rps). High centrifugal force
More informationELASTICITY (MDM 10203)
ELASTICITY () Lecture Module 3: Fundamental Stress and Strain University Tun Hussein Onn Malaysia Normal Stress inconstant stress distribution σ= dp da P = da A dimensional Area of σ and A σ A 3 dimensional
More informationSEMM Mechanics PhD Preliminary Exam Spring Consider a two-dimensional rigid motion, whose displacement field is given by
SEMM Mechanics PhD Preliminary Exam Spring 2014 1. Consider a two-dimensional rigid motion, whose displacement field is given by u(x) = [cos(β)x 1 + sin(β)x 2 X 1 ]e 1 + [ sin(β)x 1 + cos(β)x 2 X 2 ]e
More informationON THE FACE STABILITY OF TUNNELS IN WEAK ROCKS
33 rd 33 Annual rd Annual General General Conference conference of the Canadian of the Canadian Society for Society Civil Engineering for Civil Engineering 33 e Congrès général annuel de la Société canadienne
More informationMechanics of Earthquakes and Faulting
Mechanics of Earthquakes and Faulting www.geosc.psu.edu/courses/geosc508 Surface and body forces Tensors, Mohr circles. Theoretical strength of materials Defects Stress concentrations Griffith failure
More informationAdvanced model for soft soils. Modified Cam-Clay (MCC)
Advanced model for soft soils. Modified Cam-Clay (MCC) c ZACE Services Ltd August 2011 1 / 62 2 / 62 MCC: Yield surface F (σ,p c ) = q 2 + M 2 c r 2 (θ) p (p p c ) = 0 Compression meridian Θ = +π/6 -σ
More informationUsing the Timoshenko Beam Bond Model: Example Problem
Using the Timoshenko Beam Bond Model: Example Problem Authors: Nick J. BROWN John P. MORRISSEY Jin Y. OOI School of Engineering, University of Edinburgh Jian-Fei CHEN School of Planning, Architecture and
More informationThe effect of stope inclination and wall rock roughness on backfill free face stability
The effect of stope inclination and wall rock roughness on backfill free face stability Dirige, A. P. E., McNearny, R. L., and Thompson, D. S. Montana Tech of the University of Montana, Butte, Montana,
More informationUSER S MANUAL. 1D Seismic Site Response Analysis Example. University of California: San Diego.
USER S MANUAL 1D Seismic Site Response Analysis Example http://www.soilquake.net/ucsdsoilmodels/ University of California: San Diego August 2, 2017 Table of Contents USER'S MANUAL TABLE OF CONTENTS Page
More informationBearing Capacity, Comparison of Results from FEM and DS/EN DK NA 2013
NGM 2016 Reykjavik Proceedings of the 17 th Nordic Geotechnical Meeting Challenges in Nordic Geotechnic 25 th 28 th of May Bearing Capacity, Comparison of Results from FEM and DS/EN 1997-1 DK NA 2013 Bjørn
More informationMATERIAL ELASTIC ANISOTROPIC command
MATERIAL ELASTIC ANISOTROPIC command.. Synopsis The MATERIAL ELASTIC ANISOTROPIC command is used to specify the parameters associated with an anisotropic linear elastic material idealization. Syntax The
More informationModule 7: Micromechanics Lecture 29: Background of Concentric Cylinder Assemblage Model. Introduction. The Lecture Contains
Introduction In this lecture we are going to introduce a new micromechanics model to determine the fibrous composite effective properties in terms of properties of its individual phases. In this model
More informationDYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION
October 1-17,, Beijing, China DYNAMIC ANALYSIS OF PILES IN SAND BASED ON SOIL-PILE INTERACTION Mohammad M. Ahmadi 1 and Mahdi Ehsani 1 Assistant Professor, Dept. of Civil Engineering, Geotechnical Group,
More informationCavity Expansion Methods in Geomechanics
Cavity Expansion Methods in Geomechanics by Hai-Sui Yu School of Civil Engineering, University of Nottingham, U. K. KLUWER ACADEMIC PUBLISHERS DORDRECHT / BOSTON / LONDON TABLE OF CONTENTS Foreword Preface
More informationComputational models of diamond anvil cell compression
UDC 519.6 Computational models of diamond anvil cell compression A. I. Kondrat yev Independent Researcher, 5944 St. Alban Road, Pensacola, Florida 32503, USA Abstract. Diamond anvil cells (DAC) are extensively
More informationEMA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading
MA 3702 Mechanics & Materials Science (Mechanics of Materials) Chapter 2 Stress & Strain - Axial Loading MA 3702 Mechanics & Materials Science Zhe Cheng (2018) 2 Stress & Strain - Axial Loading Statics
More informationElastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount
Elastoplastic Deformation in a Wedge-Shaped Plate Caused By a Subducting Seamount Min Ding* 1, and Jian Lin 2 1 MIT/WHOI Joint Program, 2 Woods Hole Oceanographic Institution *Woods Hole Oceanographic
More informationMEMORANDUM SUBJECT: CERTIFICATE IN ROCK MECHANICS PAPER 1 : THEORY SUBJECT CODE: COMRMC MODERATOR: H YILMAZ EXAMINATION DATE: OCTOBER 2017 TIME:
MEMORANDUM SUBJECT: CERTIFICATE IN ROCK MECHANICS PAPER 1 : THEORY EXAMINER: WM BESTER SUBJECT CODE: COMRMC EXAMINATION DATE: OCTOBER 2017 TIME: MODERATOR: H YILMAZ TOTAL MARKS: [100] PASS MARK: (60%)
More informationElements of Rock Mechanics
Elements of Rock Mechanics Stress and strain Creep Constitutive equation Hooke's law Empirical relations Effects of porosity and fluids Anelasticity and viscoelasticity Reading: Shearer, 3 Stress Consider
More informationDiscrete Element Modelling of a Reinforced Concrete Structure
Discrete Element Modelling of a Reinforced Concrete Structure S. Hentz, L. Daudeville, F.-V. Donzé Laboratoire Sols, Solides, Structures, Domaine Universitaire, BP 38041 Grenoble Cedex 9 France sebastian.hentz@inpg.fr
More informationStrength of Material. Shear Strain. Dr. Attaullah Shah
Strength of Material Shear Strain Dr. Attaullah Shah Shear Strain TRIAXIAL DEFORMATION Poisson's Ratio Relationship Between E, G, and ν BIAXIAL DEFORMATION Bulk Modulus of Elasticity or Modulus of Volume
More information20. Rheology & Linear Elasticity
I Main Topics A Rheology: Macroscopic deformation behavior B Linear elasticity for homogeneous isotropic materials 10/29/18 GG303 1 Viscous (fluid) Behavior http://manoa.hawaii.edu/graduate/content/slide-lava
More informationChapter 3. Load and Stress Analysis. Lecture Slides
Lecture Slides Chapter 3 Load and Stress Analysis 2015 by McGraw Hill Education. This is proprietary material solely for authorized instructor use. Not authorized for sale or distribution in any manner.
More informationTentamen/Examination TMHL61
Avd Hållfasthetslära, IKP, Linköpings Universitet Tentamen/Examination TMHL61 Tentamen i Skademekanik och livslängdsanalys TMHL61 lördagen den 14/10 2000, kl 8-12 Solid Mechanics, IKP, Linköping University
More informationSTANDARD SAMPLE. Reduced section " Diameter. Diameter. 2" Gauge length. Radius
MATERIAL PROPERTIES TENSILE MEASUREMENT F l l 0 A 0 F STANDARD SAMPLE Reduced section 2 " 1 4 0.505" Diameter 3 4 " Diameter 2" Gauge length 3 8 " Radius TYPICAL APPARATUS Load cell Extensometer Specimen
More informationExample Application 12. Embankment Loading on a Cam-Clay Foundation
Example Application 12 Embankment Loading on a Cam-Clay Foundation 1 Foundation Soil Properties 2 Modeling Procedure Step 1 Create the foundation model grid, set boundary conditions, assign material model
More informationGeology 229 Engineering Geology. Lecture 5. Engineering Properties of Rocks (West, Ch. 6)
Geology 229 Engineering Geology Lecture 5 Engineering Properties of Rocks (West, Ch. 6) Common mechanic properties: Density; Elastic properties: - elastic modulii Outline of this Lecture 1. Uniaxial rock
More informationProf. B V S Viswanadham, Department of Civil Engineering, IIT Bombay
51 Module 4: Lecture 2 on Stress-strain relationship and Shear strength of soils Contents Stress state, Mohr s circle analysis and Pole, Principal stressspace, Stress pathsin p-q space; Mohr-coulomb failure
More information1.5 STRESS-PATH METHOD OF SETTLEMENT CALCULATION 1.5 STRESS-PATH METHOD OF SETTLEMENT CALCULATION
Module 6 Lecture 40 Evaluation of Soil Settlement - 6 Topics 1.5 STRESS-PATH METHOD OF SETTLEMENT CALCULATION 1.5.1 Definition of Stress Path 1.5. Stress and Strain Path for Consolidated Undrained Undrained
More informationMaterial is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f
Material is perfectly elastic until it undergoes brittle fracture when applied stress reaches σ f Material undergoes plastic deformation when stress exceeds yield stress σ 0 Permanent strain results from
More informationEFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES
EFFECT OF ELLIPTIC OR CIRCULAR HOLES ON THE STRESS DISTRIBUTION IN PLATES OF WOOD OR PLYWOOD CONSIDERED AS ORTHOTROPIC MATERIALS Information Revied and Reaffirmed March 1956 No. 1510 EFFECT OF ELLIPTIC
More informationThe Frictional Regime
The Frictional Regime Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 1/25/2016 10:08 AM We Discuss The Frictional Regime Processes of Brittle Deformation
More informationExamination in Damage Mechanics and Life Analysis (TMHL61) LiTH Part 1
Part 1 1. (1p) Define the Kronecker delta and explain its use. The Kronecker delta δ ij is defined as δ ij = 0 if i j 1 if i = j and it is used in tensor equations to include (δ ij = 1) or "sort out" (δ
More informationCHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES
CHAPTER THREE SYMMETRIC BENDING OF CIRCLE PLATES * Governing equations in beam and plate bending ** Solution by superposition 1.1 From Beam Bending to Plate Bending 1.2 Governing Equations For Symmetric
More informationStructural Analysis I Chapter 4 - Torsion TORSION
ORSION orsional stress results from the action of torsional or twisting moments acting about the longitudinal axis of a shaft. he effect of the application of a torsional moment, combined with appropriate
More informationSHEAR STRENGTH OF SOIL
Soil Failure Criteria SHEAR STRENGTH OF SOIL Knowledge about the shear strength of soil important for the analysis of: Bearing capacity of foundations, Slope stability, Lateral pressure on retaining structures,
More information8.1. What is meant by the shear strength of soils? Solution 8.1 Shear strength of a soil is its internal resistance to shearing stresses.
8.1. What is meant by the shear strength of soils? Solution 8.1 Shear strength of a soil is its internal resistance to shearing stresses. 8.2. Some soils show a peak shear strength. Why and what type(s)
More informationPillar strength estimates for foliated and inclined pillars in schistose material
Pillar strength estimates for foliated and inclined pillars in schistose material L.J. Lorig Itasca Consulting Group, Inc., Minneapolis, MN, USA A. Cabrera Itasca S.A., Santiago, Chile ABSTRACT: Pillar
More informationForce and Stress. Processes in Structural Geology & Tectonics. Ben van der Pluijm. WW Norton+Authors, unless noted otherwise 1/9/ :35 PM
Force and Stress Processes in Structural Geology & Tectonics Ben van der Pluijm WW Norton+Authors, unless noted otherwise 1/9/2017 12:35 PM We Discuss Force and Stress Force and Units (Trigonometry) Newtonian
More information3D ANALYSIS OF STRESSES AROUND AN UNLINED TUNNEL IN ROCK SUBJECTED TO HIGH HORIZONTAL STRESSES
3D ANALYSIS OF STRESSES AROUND AN UNLINED TUNNEL IN ROCK SUBJECTED TO HIGH HORIZONTAL STRESSES Abdel Meguid, M. Graduate Student, Department of Civil Engineering, University of Western Ontario, London,
More informationI hereby declare that, except where specifically indicated, the work submitted herein is my own original work.
Finite Element Studies on the Mechanical Stability of Arbitrarily Oriented and Inclined Wellbores Using a Vertical 3D Wellbore Model by Di Zien Low (F) Fourth-Year Undergraduate Project in Group D, 2010/2011
More informationTectonics. Lecture 12 Earthquake Faulting GNH7/GG09/GEOL4002 EARTHQUAKE SEISMOLOGY AND EARTHQUAKE HAZARD
Tectonics Lecture 12 Earthquake Faulting Plane strain 3 Strain occurs only in a plane. In the third direction strain is zero. 1 ε 2 = 0 3 2 Assumption of plane strain for faulting e.g., reverse fault:
More informationStress Integration for the Drucker-Prager Material Model Without Hardening Using the Incremental Plasticity Theory
Journal of the Serbian Society for Computational Mechanics / Vol. / No., 008 / pp. 80-89 UDC: 59.74:004.0 Stress Integration for the Drucker-rager Material Model Without Hardening Using the Incremental
More informationCVEEN 7330 Modeling Exercise 2c
CVEEN 7330 Modeling Exercise 2c Table of Contents Table of Contents... 1 Objectives:... 2 FLAC Input:... 2 DEEPSOIL INPUTS:... 5 Required Outputs from FLAC:... 6 Required Output from DEEPSOIL:... 6 Additional
More informationA non-linear elastic/perfectly plastic analysis for plane strain undrained expansion tests
Bolton, M. D. & Whittle, R. W. (999). GeÂotechnique 49, No., 33±4 TECHNICAL NOTE A non-linear elastic/perfectly plastic analysis for plane strain undrained expansion tests M. D. BOLTON and R. W. WHITTLE{
More informationPILE-SUPPORTED RAFT FOUNDATION SYSTEM
PILE-SUPPORTED RAFT FOUNDATION SYSTEM Emre Biringen, Bechtel Power Corporation, Frederick, Maryland, USA Mohab Sabry, Bechtel Power Corporation, Frederick, Maryland, USA Over the past decades, there has
More informationWhen you are standing on a flat surface, what is the normal stress you exert on the ground? What is the shear stress?
When you are standing on a flat surface, what is the normal stress you exert on the ground? What is the shear stress? How could you exert a non-zero shear stress on the ground? Hydrostatic Pressure (fluids)
More informationLecture #7: Basic Notions of Fracture Mechanics Ductile Fracture
Lecture #7: Basic Notions of Fracture Mechanics Ductile Fracture by Dirk Mohr ETH Zurich, Department of Mechanical and Process Engineering, Chair of Computational Modeling of Materials in Manufacturing
More informationME 7502 Lecture 2 Effective Properties of Particulate and Unidirectional Composites
ME 75 Lecture Effective Properties of Particulate and Unidirectional Composites Concepts from Elasticit Theor Statistical Homogeneit, Representative Volume Element, Composite Material Effective Stress-
More informationUNIT I SIMPLE STRESSES AND STRAINS
Subject with Code : SM-1(15A01303) Year & Sem: II-B.Tech & I-Sem SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road 517583 QUESTION BANK (DESCRIPTIVE) UNIT I SIMPLE STRESSES
More informationSupplementary Figures
Fracture Strength (GPa) Supplementary Figures a b 10 R=0.88 mm 1 0.1 Gordon et al Zhu et al Tang et al im et al 5 7 6 4 This work 5 50 500 Si Nanowire Diameter (nm) Supplementary Figure 1: (a) TEM image
More informationBrittle Deformation. Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm
Lecture 6 Brittle Deformation Earth Structure (2 nd Edition), 2004 W.W. Norton & Co, New York Slide show by Ben van der Pluijm WW Norton, unless noted otherwise Brittle deformation EarthStructure (2 nd
More information