Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
Elastic Bending of Tectonic Plates Two-Dimensional Bending of Flexure of Plates • Some loads cause the lithosphere to bend. • Flexure can be used to deduce elastic properties and thicknesses of plates. • A simple example of plate bending is cylindrical bending from heavy loads in which the width of the plate >> length. • The deflection can be determined by requiring it to be in equilibrium under the action of all the forces and torques exerted on it. • Forces: Q (downward force per unit area) V (net shear force per unit length) P (horizontal force per unit length) Because the net acceleration is zero, the net vertical force is zero. Hence: q(x)dx + dV =0 or q = -dV/dx • The net bending moment, M per unit length, is the integrated effect of moments exerted by normal forces on the cross-section (see figure next slide). • The moments M and M+dM give a net counterclockwise torque dM while V (shear) forces exert a torque Vdx in a clockwise sense. The P forces exert a counter clockwise torque (-Pdw). • The net torque is zero so dM –Pdw = Vdx (w is deflection of the plate) • Differentiating twice with respect to x, one obtains the following relationship: d2M/dx2 = -q + P* d2w/dx2 • M is inversely proportional to the local radius of curvature R (1/R = -d2w/dx2) -As the plate is deflected downward, the top part is contracted (σxx is positive) and the lower part is extended (σxx is negative). σxx is zero on the midplane at y=0 (neutral unstrained surface) -The bending moment is M= ∫(-h/2 to h/2)σxxy dy where h is the thickness of the plate. • The change in length (Δ l) = -y(l/R) where y is the distance from the midplane and R s the radius of curvature. • • • R Φ = dx, hence the subtending angle Φ = -d2w/dx2 * dx 1/R = Φ/l = Φ/dx = - d2w/dx2 = D M= D/R (flexural rigidity/ radius of curvature) • So… the general equation for the deflection of the plate is D d4w/dx4 = q(x) –P d2w/dx2 Bending of Plates under Applied Moments and Vertical Loads •Let us consider a plate embedded on only one end. •If we assume the plate is weightless, q=0, V=0 (no applied force), P=0. •Therefore, M is constant and equal to the applied torque. •If w=0 at x=0 and dw/dx =0 at x=0 we find that W=-Mqx2/2D (the plate has the shape of a parabola) • The bending of a plate embedded at one end with uniform loading: Q= constant, P=0 • D d4w/dx4 = q/D is the governing equation. • With the four boundary conditions, w=dw/dx =0 at x=0 and d2w/dx2 =0 at x=L, one obtains the following equation: W=qx2/D (x2/24 –Lx/6 + L2/4) • For this particular example, the maximum bending stress at x=0 is 3qL2/h2. The ratio of shear stress to maximum bending stress at x=0 is h/3L, which tells us that shear stresses in plates are small when compared to bending stresses. Buckling of a Plate under a Horizontal Load • Let us consider a plate pinned at both ends with an applied horizontal force, P. • The governing equation is D d4w/dx4 + P d2w/dx2 =0 Integrate twice: D d2w/dx2 + Pw = c1x + c2 Buckling of a Plate under a Horizontal Load, con’t • Solution of the differential equation gives the general solution of w=c1 sin (P/D)1/2 x • If c is not zero, sin (P/D)1/2L =0 so (P/D)1/2L = nπ, n= 1, 2, 3… • This tells us that P=(n2π2D)/L2, a series of values of which the smallest is Pc = (π2)D/L2 • This makes a half sine curve where w=c1 sin πx/L. c1, the amplitude, can not be determined by linear analysis. • Plate flexure theory in mountain belts is much more complicated due to layering. Deformation of Strata Overlying an Igneous Intrusion • A laccolith is a sill-like igneous intrusion with a round lense-like shape that forms when magma intrudes in layered rocks and magma raises and deforms the overburden. • The loading of the plate is the part of the upward pressure force in excess of the lithostatic pressure. Q= -p + ρgh • If we consider x=0 to be the center of the laccolith, x= +/- L/2 • With some algebra we obtain… W=- (p-ρgh)/24D (x4– L2x2/2 + L4/16) where the maximum deflection at x=0 is Wo = -(p – ρgh)L4/384D. Application to the Earth’s Lithosphere • We must remember to include the hydrostatic restoring force. For example, when oceanic lithosphere is depressed, water fills in space vacated by the mantle’s moved rocks. Continental crust is thickened by the amount that the Moho is depressed. • If water fills, the weight per unit area of a vertical column is ρwg(hw + w) + ρmgh • Deeper, where there is no plate deflection in the mantle, the weight per unit area is ρwghw + ρmg(h+w) • so there is an upward hydrostatic force of (ρm –ρw)gw which tries to restore the lithosphere to its original configuration • When the net force per unit area is q=qa- (ρm–ρw)gw, where qa is the applied load, the governing equation becomes: Dd4w/dx4 +Pd2w/dx2 + (ρm-ρw)gw = qa(x). Application to the Earth’s Lithosphere Periodic Loading (to answer the question of how mountains and valleys affect the lithosphere) • The elevation of the Earth’s topography can be approximated by using the following formula : h=ho sin 2 π x/λ where h is the topographic height and λ is the wavelength. • The load is qa(x) = ρcgho sin 2π x/λ, so the governing equation is: Dd4w/dx4 + (ρm – ρc)gw = ρcgho sin 2π x/λ • Because the loading is periodic, the deflection is also: Wo= ho/ ((ρm/ρc) -1 + D/ρg (2π/λ)4) • If the wavelength is short (λ<< 2π (D/ρcg) then (1/4wo << ho ) • Thus, short wavelength topography causes virtually no deformation of the lithosphere. Application to the Earth’s Lithosphere • If the wavelength is long (λ >> 2 π (D/ρcg)), 1/4w = ρcho/(ρm-ρc) • For long wavelength topography, the lithosphere has no rigidity, and the topography is in hydrostatic equilibrium. • The degree of compensation, C, of the topographic load is the ratio of the deflection of the lithosphere to its maximum or hydrostatic deflection. Stability of the Earth’s Lithosphere under an End Load (Buckling) • When the lithosphere is subjected to a horizontal force, P, the governing equation is: Dd4w/dx4 + Pd2w/dx2 + (ρm– ρw) gw = 0. • This equation can be satisfied by a sinusoidal deflection, giving: D(2π/λ)4 –P(2π/λ)2+ (ρm-ρw)g =0 • The solution provides, Pc, the minimum value of P for which the initially horizontal lithosphere will become unstable and accept the sinusoidal shape. Pc = 4Dg(ρm – ρw) ½ • The critical wavelength is λc = 2π (2D/ρc)1/2= 2π (Eh3/(12(1-v2)( ρm-ρw)g))1/4 Stability of the Earth’s Lithosphere under an End Load (Buckling) • Now, to determine whether buckling of the lithosphere can lead to the formation of a series of synclines and anticlines: • Let us consider an elastic lithosphere with a thickness of 50 km. Taking E=100 GPa, v= .25, ρm = 3300 kg/m3, and ρw = 1000 kg/m3, we find that the critical stress is 6.4 GPa. Therefore the lithosphere can support a horizontal compression of 6.4 GPa without buckling. Because the stress is very large, one can conclude that such buckling does not occur. • However, the buckling of thin elastic layers may contribute to the formation of folded structures. Bending of the Elastic Lithosphere under the Loads of Island Chains • Volcanic islands provide loads that cause the lithosphere to bend. An example is Hawaii, which has very deep seas surrounding the island (Hawaiian Deep). Bending of the Elastic Lithosphere under the Loads of Island Chains,con’t • The governing equation is: Dd4w/dx4 + (ρm –ρw)gw = 0 whose solution is: W= e(x/α) (c1 cos (x/α) + c2 sin (x/α)) + e(-x/α)(c3 cos (x/α) + c4 sin (x/α)) • α is the flexural parameter and is equal to (4D/ ((ρm –ρw)g)¼ • The above equation becomes w=c3 e(-x/α)(cos (x/α) + sin (w/α)), where x is positive and c3 is Voα3/8D, the maximum amplitude of deflection. Bending of the Elastic Lithosphere under the Loads of Island Chains,con’t • In the case of Hawaii, where volcanism may have weakened the lithosphere, we might also consider a case where the lithosphere is fractured along the line of the ridge. For the same load as the first case, the fractured model has twice the deflection. Bending of the Elastic Lithosphere at an Ocean Trench • Prior to subduction, considerable bending of the lithosphere occurs. To model this, we consider an elastic plate acted upon by an end load Vo, and a bending moment Mo. • Solution of the governing equation gives: w= e(-x/α) [c3 cos (x/α) + c4 sin (x/α)] where c4 = -Mo α2/2D and c3 = (Voα + Mo) α2/2D • Vo and Mo can not be determined directly, but the half-width of the forebulge is π/4 α, from which we can derive very complicated formulas to use. Bending of the Elastic Lithosphere at an Ocean Trench, con’t • The plot of (w/wb) vs. (x-xo)/(xb/xo) defines a universal flexure profile which is valid for any two-dimensional elastic flexure of the lithosphere under end loading. • A comparison of the universal flexure profile and subduction at the Marianas is excellent. Flexure and the Structure of Sedimentary Basins • Lithospheric flexure is also associated with the structure of many sedimentary basins. Sedimentary basins are areas of the Earth’s surface that have been depressed and the depressions filled with sediments. They are often important reservoirs of petroleum. • The loading of mountain belts can be used to explain some foreland basins, as in the Appalachian and Andes Mountains.