Force-Balance and Mohr's Equations: Structural Geology
Force-Balance and Mohr's Equations: Structural Geology
Force-Balance and Mohr's Equations: Structural Geology
11 - Force-Balance
Consider the forces acting parallel to the front face of the triangular solid shown in Figure 11-1. These forces are acting in a plane so that this force-balance problem is two dimensional. If the area of the hypothenuse face is unity (i.e. 1), then the area of the left face is unity times sin whereas the area of the bottom face is unity times cos . Forces acting on the three faces are stresses multiplied by the area of the face [Force = (force/area) (area)]. Summing the forces in the horizontal and vertical directions (the forces add to zero in both directions), we obtain Fh = 0 = Fx - S sin - N cos Fv = 0 = Fy - N sin + S cos (11-1a) (11-1b)
Rewriting the force terms using components of stress, we obtain the normal and shear stresses acting on a plane an angle of to 1 in terms of the principal stresses. 1 cos n cos sin = 0 3 sin n sin + cos = 0 (11-2a) (11-2b)
In many geological applications the plane of interest is a fault plane. Now we can solve these two equations for n and n = 3 cos2 + 1 sin2 = (1 - 3) cos sin Remember that sin2 = 2 sin cos and also cos2 = 1/2(1 + cos2 ) and sin2 = 1/2(1 - cos2 ) We derive n = and =
1 1
(11-3a) (11-3b)
(11-4)
(11-5)
(11-6a)
/2(1 - 3)sin2 .
(11-6b)
Lect. 11 - Force-Balance
Note that in these equations is the angle between 1 and the normal to the fault plane. This convention leads to a positive sign in the middle of equation 11-6. Chapter 10 of Marshak and Mitra (1988) uses the convention that is the angle between the normal to the fault plane and 1. If the is the angle between 1 and the fault plane the sign convention equation 1-29a becomes n =
1
(11-7).
Unless each convention is memorized the reader will be confused when using these equations in the future. The general force-balance problem is shown in figure 11-2 where the back and side faces of the unit triangle are subject to both normal and shear stresses. In other words the coordinate system x-y is not parallel to the principal stress directions. The general case is given to illustrate a property of stress called invariance with respect to coordinate system. In two dimensions we define stress as a force per unit of line length, in contrast to the three dimensional situation where stress is a force per unit area (Fig. 11-1). Consider a coordinate system Oxy with an arbitrary line AB cutting the x and y axes such that the normal to the line AB makes an angle with the x-axis. This gives a right triangle AOB with sides OA (parallel to Ox) and OB (parallel to Oy) and a hypothenuse AB. Across the line AB a stress vector p can be applied making an angle with the x-axis. Remember that p = f/A when A 0, so a stress vector can represent stress at a point. Otherwise stress is defined on a line (2-D) or a surface (3-D). The stress vector p can be resolved into components parallel to the x and y axes: p = px + py. Even though P is called a vector it still has units of stress (force/length in two dimensions). Because the triangle ABO is in equilibrium the sum of the force-vectors on all sides must balance. In 2-D stress multiplied by line length will give a force vector. So pxAB = xOB + yxOA. (11-8)
Lect. 11 - Force-Balance
(Fig. 11-2)
If length a = AB, then acos = OB and asin = OA. If we divide by length a, then px = x cos + yx sin. Likewise py = y sin + xy cos (11-9b) (11-9a)
Now consider a normal stress n and shear stress across AB in terms of the components of the stress vector n = px cos + py sin The equality holds for = py cos - px sin. (11-9b) (11-9a)
Substituting for px and py and remembering that sin2 = 2 sin cos and also cos 2 = cos2 - sin2
Lect. 11 - Force-Balance n = x cos2 + 2xy sin cos + y sin2 and = 1/2(y - x) sin2 + xy cos2. (11-10b) (11-10a)
These are the general equations for a stress system in which the orientation of the principal stresses are unknown. The following exercise illustrates the case for the principal stresses being parallel to the coordinate axes.
(Fig. 11-3) Now we consider the effect of rotating the coordinate system on the values of stress (Fig. 11-3). To accomplish this coordinate rotation let n and be x' and x'y' where Ox'y' is rotated from Oxy. From a previous equation we have x' = x cos2 + 2xy sin cos + y sin2 (1111a)
Now we must find y' in the new coordinate system. This is accomplished by replacing by + 1/2 and we get y' = x sin2 - 2 xy sin cos + y cos2 If we add x' and y' and remember that sin2 + cos2 = 1, we get x' + y' = x + y. (11-12) (11-11b)
Lect. 11 - Force-Balance This shows that the sum of the normal stresses is invariant or unchanged by rotation of the coordinate system. Likewise x'y' = 0.5(y - x) sin2 + xy cos2. (11-13)
Principal stresses are found in planes containing no shear stress. If in these planes, x"y" = 0, then from the previous equation 2xy tan2 = ____ x - y (11-14)
where is the one angle between the coordinate system Ox"y" and Oxy where the shear stresses vanish along the directions Ox" and Oy". In this coordinate system the only stresses are the normal stresses x" and y". This coordinate system contains the principal stress axes and the components, x" and y" are known as the principal stresses.
(Fig. 11-4) Let us now call these principal stresses 1 and 2 (Fig. 10-4) and then we shall choose the coordinate system such that the x and y axes are in the direction of the principal stresses, 1 and 2. Now the normal n and shear stresses across a line whose normal is inclined at to 1 is n = 1.cos2 + 2.sin2 (11-15)
Lect. 11 - Force-Balance
This is the convention used in Chapter 10 of Marshak and Mitra (1988). Again remember that cos2 = 1/2(1 + cos2 ) and sin2 = 1/2(1 - cos2 ) we can rearrange the equation for the normal stress n = and shear stress =
1 1
(repeat of 11-6a)
From this last equation we see that shear stress is greatest when = /4 and 3/4. Using this same coordinate system where the axes are parallel to the principal stresses 1 and 2, we can look at the stress vectors px and py. They become px = 1.cos and py = 2.sin . (11-16)
(Fig. 11-5) Substituting into the equation sin2 + cos2 = 1, we can generate the equation for an ellipse where the px and py are on the ellipse called the stress ellipse.
The semi axes of the ellipse are 1 and 2. The Mohr's circle is a graphical method of representing the state of stress of a rock in two dimensions (Fig. 11-5). The equations used for the Mohr's circle representation are those derived above for the coordinate system with axes parallel to the principal stresses. The Mohr's circle may be used to derive the normal n and shear stresses on any plane whose normal is oriented at from 1. The coordinate system for the Mohr's circle representation is n along the horizontal axes with increasing compression to the right and along the vertical axes. Critical points along the n-axes are OP = 1, OQ = 2, and C = 1/2(1 + 2). The angle measured as PCA counter clockwise from OP is 2. Now we have n = OB = OC + CB = 1/2(1 + 2) + 1/2(1 - 2)cos2 and = AB = 1/2(1 - 2)sin2. (11-18b) (11-18a)