Page 310 - Physical Principles of Sedimentary Basin Analysis
P. 310

292                        Flexure of the lithosphere
                                             ∞

                                               w(x)dx = 2w max α.                   (9.39)
                                            −∞
                 The net area of the deflection is equal to two times the area of a rectangle with sides of
                 maximum deflection (w max ) and the length α.

                 Note 9.2 Equation (9.31) is solved for only a half-plate (x > 0) because the problem is
                 symmetric around x = 0. The half-plate formulation makes the load disappear from the
                 equation and it appears instead at the boundary x = 0. Equation (9.31) is first rewritten as

                                              1 d
                                                 4
                                                    + 1 w = 0                       (9.40)
                                              4 d ˆx 4
                                                x
                 using the dimensionless x-coordinate ˆ = x/α. The parameter α has unit length and it is
                 the characteristic length scale of the deflection problem. The solution of equation (9.40)is
                 done by solving the corresponding characteristic equation

                                                1  4
                                                 r + 1 = 0                          (9.41)
                                                4
                 which has four roots r k with k = 1, 2, 3, 4. The general solution of the fourth-order
                 equation (9.40) in terms of the roots r k is

                                                     4
                                                           x
                                                x
                                             w(ˆ) =    c k e r k ˆ                  (9.42)
                                                    k=1
                 which is seen by inserting the solution into equation (9.40). The characteristic equation is
                                                4
                                             r            iπ+i2πk

                                            √     =−1 = e                           (9.43)
                                              2
                 which gives
                                                 √
                                            r k =  2 e i(1/4+k/2)π                  (9.44)

                 for k = 1,..., 4. The four roots of the characteristic equation are r 1 = e iπ/4  = 1 + i,
                 r 2 = e i3π/4  =−1 + i, r 3 = e i5π/4  =−1 − i and r 4 = e i7π/4  = 1 − i. We recall that
                         a
                 e a+ib  = e (cos b +i sin b), and since we are solving the flexure problem for x > 0 we can
                 only use r 2 and r 3 , because r 1 and r 4 have positive real parts. Solution (9.42)is

                                    x
                                  w(ˆ) = e  −ˆ x  (c 2 + c 3 ) cos ˆx + i(c 2 − c 3 ) sin ˆx  (9.45)
                 when r 2 and r 3 are inserted. We want a solution with only real numbers, which is achieved
                 by replacing the coefficients c 2 and c 3 by real coefficients d 1 = c 2 +c 3 and d 2 = i(c 2 −c 3 ).
                 The solution in terms of the real coefficients d 1 and d 2 is then

                                          x
                                       w(ˆ) = e −ˆ x  d 1 cos ˆx + d 2 sin ˆx .     (9.46)
   305   306   307   308   309   310   311   312   313   314   315