1
1
1
Abstract
This study is an attempt to solve the general governing differential equation for buckling of beams resting on elastic foundation for different types of beams (pin ended, fixed ended and cantilever beams) taking into account the effect of soil sub-grade reaction value (K), axial load, lateral load and relative end settlement on the deflection and bending moment of the beams for prismatic and nonprismatic beams. Approximate solutions are also presented using the extermization of the total potential energy equation (Rayleigh-Ritz method) using approximate shape functions for the deflection which satisfy the boundary conditions of the beam. The results show that one term series approximate solution gives acceptable results in comparison with the exact solution for practical case, the accuracy of the solution increases with increasing the number of terms up to 5 terms, beyond this limit the accuracy of the solution does not change. Buckling load increases linearly with increasing (K) while the deflection and the bending moment exponentially decrease with increasing (K). Keywords: Bending moment, Buckling, Cantilever beam, Deflection, Fixed ended beam, Non-prismatic beam, Pin ended beam, Prismatic beam, Settlement and Soil sub-grade reaction.
(axial load)
.(K)
(exponentially )
104
Vol.16
No.3
Aug.
2008
E : Modulus of elasticity. f(x) : Function of changing the moment of inertia. F : Tapered beam factor = F ( ) h1 : Smallest depth of the beam. h2 : Largest depth of the beam. I : Moment of inertia of the beam cross section. K : Soil sub-grade reaction. L : Span length of the beam. M : Bending moment. P : Axial load. Pcr : Euler buckling load. Po : Euler buckling load when (K=0). Q(x) : Lateral load. Qo : Intensity of the lateral load. w : Deflection of the beam. : Relative end settlement of the beam. : Total potential energy of the beam. : P/EI : K/EI : h2/h1 1
Introduction:
The critical buckling load for the beam under axial compression can be determined by considering the simultaneous action of compressive and bending moment by assuming an initial curvature or by considering the behavior of an ideal column, which is assumed initially to be perfectly straight, compressed by centrically applied load, the column to be perfectly elastic and the stress within the proportional limit. If the axial load (P Pcr); the bar remain straight and the column is in elastic stability condition, which means that if lateral force is applied, a small deflection is produced. The deflection disappears when the lateral load is removed and the bar retains its original form. If the axial load (P > Pcr); the bar becomes unstable and a small lateral force will produce a deflection which does not disappear when the lateral load is removed [1-3]. In beam on elastic foundation, the elastic media can be represented by an equally spaced elastic supports of equal rigidity. The reaction of the medium at any cross section of the bar is then proportional to the deflection at that section. The stability can be increased by using non-prismatic bars, these type of sections are very often used in steel structures. The solution of the differential equation of the gradually changing cross section including a truncated cone and pyramid is presented in ref. [1]. The main objectives of this study are to: 1- Extend the buckling problem of the beams on elastic foundation to include three types of beams (pin ended, fixed ended and cantilever) using exact solution of the governing differential equation and approximate method using (Rayleigh-Ritz) energy method. 2- Extend the analysis to include non-prismatic beams whose depths are varying linearly or parabolically. 3- Extend the analysis to include the effect of the relative end settlement on the deflection and the bending moment of the beam using both exact solution and approximate solution.
105
Hassan : Buckling Of Beams On Elastic Foundations Theory And Analysis: The total governing differential functional for the beam resting on elastic foundation subjected to axial compression force (P) and uniform lateral distributed load Q(x) can be written as [1-3]: = [EI/2 (w")2 + K/2 (w)2 P/2 (w')2 Q(x) w ] dx ------------(1) And the governing differential equation of the problem (EulerLagrange equation) for this functional can be written as below: ------------(2) EI wIV + P wII + K w = Q(x) I- Elementary buckling of prismatic beams [1,2, 4-8]: The governing differential equation for beams subjected to axial compressive force (P) and uniform lateral distributed load Q(x) for (K=0) becomes as: EI wIV + P wII = Q(x) ------------(3) The principle of superposition is not possible in equation (3) between the two conditions: 1- Q(x)=0; EI wIV + P wII = 0 ------------(4) 2- P = 0 ; EI wIV = Q(x) ------------(5) Because the effect of axial force on the deflection (w) will directly depends on the deflection that Q(x) induces, this means there is a non-linearity action in the behavior when (P) and Q(x) act together. The exact solution of equation (4) is given by: w = A sin(k1x) + B cos(k1x) + C x + D --------------(6) where: k12=P/EI The constants A, B, C and D are determined for different types of beams by applying the boundary conditions at the ends. 1- Pin ended beam: the deflection (w=0) and the bending moment (w''=0) at (x=0) and (x=L), apply these conditions to equation (6) and solve for the constants to obtain: B=C=D=0 and k1=n /L The deflection equation w=A sin(n x/L) and the buckling load Pcr= n2 2 EI/L2 2- Fixed ended beam: the deflection (w=0) and the slope (w'=0) at (x=0) and (x=L), the constants obtained from these equations: B+D=0 -----------(7-a) A k1 + C =0 -----------(7-b) A sin(k1L) + B cos(k1L) + C L + D = 0 -----------(7-c) A k1cos(k1L) B k1 sin(k1L) + C = 0 -----------(7-d) The characteristic equation is: sin(k1L/2) [k1L/2 cos(k1L/2) - sin(k1L/2)] = 0 ----------(8) one solution of this equation is sin(k1L/2)=0 and then (k1=2n /L) The deflection equation w=B[cos(2n x/L)1] and the buckling load Pcr= 4n2 2 EI/L2 3- Cantilever beam: the deflection (w=0) and the slope (w'=0) at (x=0) and the bending moment (w''=0) and the shear force (w''+k12w' =0) at (x=L). The deflection equation w=B [1cos ((2n1) x/2L)] and k1= (2n1) /2L and the buckling load Pcr = (2n1)2 2 EI/4L2 4- Fixed Hinged beam: the deflection (w=0) and the slope (w'=0) at (x=0) and the deflection (w=0) and the bending moment (w''=0) at (x=L) The constants obtained from these equations: B+D=0 -----------(9-a) A k1 + C =0 -----------(9-b) CL+D=0 -----------(9-c) A sin(k1L) + B cos(k1L) = 0 ----------(9-d) 106
Al-Rafidain Engineering
Vol.16
No.3
Aug.
2008
The characteristic equation is: tan(k1L) = k1L ---------(10) The smallest solution is k1L=4.493 and the buckling load is Pcr = 2 2 EI/L2 The buckling load can also be determined using stationary total potential energy method [3] as below: Pcr = EI (w'')2 dx / (w')2 dx ----------------(11) using the stress resultant displacement relation (M= EI w'' ), equation (11) becomes: -----------------(12) Pcr = M2/EI dx / (w')2 dx
II- Buckling of beams on elastic foundation [1-3]: 1- Pin ended beam: The general expression for the deflection curve of a bar with hinged ends can be represented by series of sine waves as shown in previous section: w = An sin n x/L ------------(13) Substitute in equation (1) and extremize the functional with respect to (An), the critical buckling load is: Pcr = 2 EI/L2 [n2 + KL4/(n2 4EI)] ----------------(14) or Pcr/Po = 1 + KL4/(n4 4EI)] ----------------(15) where n is an integer and Po is the Euler buckling load of the beam when (K=0). Equation (14) gives the critical load (Pcr) as a function of (n), which represents the number of half sine waves in which the bar is subdivided at buckling. When (K=0), the resulting critical load (Pcr = n2 2 EI/L2) is the exact form of Euler load. If (K) is very small and greater than zero, (n) should be taken as unity, thus for very flexible elastic foundation, the bar buckles without an intermediate inflection point. By gradually increasing (K), one finally conclude that at a condition (n=2) the critical load is smaller than at (n=1), but equation (14) should give the same result at (n=1) and (n=2) to give the limiting value of (K) that satisfies the condition of buckling with an intermediate inflection point. 1 + KL4/( 4 EI) = 4 + KL4/(4 4 EI) ----------------(16) 4 4 4 4 From which KL /( EI) = 4 and Ko = 4 EI/L when K Ko; the deflection curve of the buckled bar has no inflection point and (n=1). when K>Ko; the deflection curve of the buckled bar has an inflection point and (n=2) and the span is subdivided into two half sine waves. To generalize equation (16), the equation is to be written in the following form: n2 + KL4/(n2 4 EI) = (n+1)2 + KL4/((n+1)2 4 EI) ----------------(17) From which KL4/( 4 EI) = n2 (n+1)2 and Kon= n2 (n+1)2 4 EI/L4 Neglecting (1) as compared to (n), then: KL4/( 4 EI) = n4 or (L/n)4 = 4 EI/K and Kon= n4 4 EI/L4 Substituting in equation (14), to obtain: Pcr = 2 n2 2 EI/L2 or Pcr = 2 EI K The critical load (Pcr) of the beam on elastic foundation is twice the critical load for beam when (K=0) for the same properties. The governing differential equation (2) can be written in the following form: wIV + wII + w = Q(x)/EI ------------(18) where: = P/EI and = K/EI For the case (K>0) and ( /2)2 <( /2), the general solution of equation (18) is: w = A1 cos( 1x)+ A2 sin( 1x)+ A3 cos( 2x)+ A4 sin( 2x)+ Qo/K ------------(19) where: 107
( /2) ( /2)2 and 2= ( /2) + ( /2)2 Imposing the boundary condition of pin ended beam, A1=A3=0 and 1= 2= n /L and the deflection equation w=(A2+A4) sin(n x/L) or w=A sin(n x/L), this is the same shape function which is used in energy method and derived in previous section. The critical load equation is: Pcr/Po = 1 + KL4/(n4 4 EI) --------------(20) 2 2 where Po is the Euler load of the pin ended beam when (K=0) = n EI/L2 This equation is exactly similar to that obtained by energy (Rayleigh-Ritz) method. For the case ( /2)2 = or Pcr = 2 EI K , the general solution is: w = A1 cos( x)+ A2 sin( x)+ A3 x cos( x)+ A4 x sin( x)+ Qo/K ------------(21) 2 where: = P/(2EI) Imposing the boundary condition of hinged beam, and solve for the constants, one obtain: A1= Qo/K -----------(22-a) A4= Qo /(2K) ----------(22-b) A3= Qo /(2K) [cos( L)-1]/[ L cos( L) sin( L)] ----------(22-c) A2sin( L)=Qo/K[cos L+ L/2sin L1 L/2cos L(cos L-1)/( Lcos Lsin L)] -(22-d) For (Qo=0), A1= A3= A4=0 and A2 sin( L) = 0, then =n /L Also the same deflection equation is obtained w= A2 sin(n /L) and the buckling load Pcr = 2 n2 2 EI/L2 and K= n4 4 EI/L4 , these equations are similar to those obtained by the energy method. 2- Fixed ended beam: the deflection (w=0) and the slope (w'=0) at (x=0) and (x=L), apply these boundary conditions into equation (21) for case ( /2)2 = and solve for the constants: A1= Qo/K ------------(23-a) A2= Qo/K [ L + cos( L) sin( L)]/[sin( L)2 ( L)2] ------------(23-b) A3= Qo /K [ L + cos( L) sin( L)]/[sin( L)2 ( L)2] ------------(23-c) A4= Qo /K [sin( L)2]/[sin( L)2 ( L)2] ------------(23-d) The characteristic or transcendental equation is: 12 Z2 sinZ + Z sinZ sin4Z (2 sinZ + Z cosZ)(1-cos4Z) = 0 -----------(24) where: Z = L/2 One solution of the above equation is (Z=n ), then =2n /L= P/2EI The critical load in this case Pcr = 8 n2 2 EI/L2 But Po = 4 n2 2 EI/L2 , thus Pcr = 2 Po This is the same as that concluded in a pin ended beam, the critical load (Pcr) of the beam on elastic foundation is twice the critical load for beam when (K=0) for the same properties and for the case ( /2)2 = . The general expression for the deflection curve of a bar with fixed ends can be written as a series of cosine terms as shown in the previous section: w = An [1 cos(2n x/L)] ----------------(25) Substitute in equation (1) and extremis the functional with respect to (An) to obtain: Pcr = 4 2 EI/L2 [n2 + 3KL4/(16n2 4 EI)] ----------------(26) or Pcr/Po = 1 + 3KL4/(16n4 4 EI)] ----------------(27) when (K) is very small and greater than zero, (n) should be taken as unity, thus for very flexible elastic foundation, the bar buckle without an intermediate inflection point. By gradually increasing of (K), we finally arrive that at a condition for (n=2) the critical load is smaller than at (n=1), but equation (26) should gives the same result at (n=1) and (n=2). The limiting value of (K) that satisfy the condition of buckling with an intermediate inflection point. 1 + 3KL4/(16 4 EI) = 4 + 1/4 3KL4/(16 4 EI) ----------------(28) From which 3KL4/(16 4 EI) = 4 and Ko = 64/3 4 EI/L4 108
1=
Al-Rafidain Engineering
Vol.16
No.3
Aug.
2008
when K Ko, n=1; the deflection curve of the buckled bar has no inflection point. when K>Ko, n=2; the deflection curve of the buckled bar has an inflection point.and the span subdivided into two half sine waves. To generalize equation (28), the equation is to be written in the following form: n2 + 1/n2 3KL4/(16 4 EI) = (n+1)2 + 1/(n+1)2 3KL4/(16 4 EI) ---------------(29) 4 4 2 2 2 2 4 From which 3KL /(16 EI) = n (n+1) and Kon= n (n+1) 16 EI/(3L4) Neglecting (1) as compared to (n), then: 3KL4/(16 4 EI) = n4 or (L/n)4 = 16 4EI/(3K) and Kon= 16n4 4 EI/(3L4) Substituting in equation (26), to obtain: Pcr = 8 n2 2 EI/L2 or Pcr = 3 (2 EI K) This is exactly similar to the exact solution. Also (Pcr) for fixed ended beam is greater than pin ended beam by ( 3) with the same properties and for case ( /2)2 = . 3- Cantilever beam: the general approximate expression for the deflection curve that satisfies the boundary conditions is: w = An [1 cos(n x/2L)] ----------------(30) And following the same previous procedures, the critical buckling load is: Pcr/Po = 1 + 16[3n -8sin(n /2)]/(n ) KL4/(n4 4 EI) ----------------(31) For n=1, Pcr/Po = 1 + [16(3 -8)/ ] KL4/( 4 EI) ----------------(32) The limiting value of (K) that satisfy the condition of buckling with an intermediate inflection point is Ko = 0.55 4 EI/L4 Taking (n) terms series and following the same previous steps: 16[3n -8sin(n /2)]/(n ) KL4/( 4 EI) = n4 ----------------(33) Substituting in equation (31), to obtaine: Pcr = [3n -8sin(n /2)]/(n ) (2 EI K) --------(34) 5 5 4 and Kon = n EI/[16{3n -8sin(n /2)}L ] ---------(35) For n=1; Pcr = 2/3 (2 EI K) This means that (Pcr) for a cantilever beam is (2/3) of the pin ended beam with the same properties and for the case ( /2)2 = . Table (1) shows (Ko) and (Pcr) of the beams with respect to values of pin ended beam. Table (1): Comparison of (Ko and Pcr) for different beams Beam type Fixed ended beam Cantilever beam Ko/Ko(pin ended) 16/3 0.14 Pcr/Pcr(pin ended) 3 2/3 Equations (20, 27 and 31) can be written in one general form as: Pcr/Po = 1 + a KL4/( 4EI) or PcrL2/( 2EI) = a1 + a2 KL4/( 4EI) where (a, a1 and a2) are constants shown in Table (2)
-----------------(36) -----------------(37)
Table (2): Value of the constants (a, a1 and a2) Beam type a a1 a2 Pin ended 1 1 1 Fixed ended 3/16 4 3/4 Cantilever 16(3 -8)/ 1/4 4(3 -8)/
109
Figs.(1 and 2) show the plot of equations (36 and 37), as shown the value of buckling load of beams resting on elastic foundation increased linearly with increasing of the soil sub-grade reaction value.
III- Settlement of beams on elastic foundation: The governing differential equation of the problem (EulerLagrange equation) for beams resting on elastic foundation subjected to lateral distributed load can be written as: EI wIV + K w = Q(x) ------------(38) The general exact solution is: w = (k1 cosh x + k2 sinh x) cos x + (k3 cosh x + k4 sinh x) sin x + Qo/K ------ (39) where 4 = K/(4EI) and using the stress resultant displacement relation (M= EI w'' ), the bending moment equation is: M= 2EI 2 [(k1 sinh x + k2 cosh x) sin x + (k3 sinh x + k4 cosh x) cos x] -----(40) The constants k1, k2,k3 and k4 are determined for different types of beams by applying the boundary conditions at the ends. The relative end settlement is used as a boundary condition at one ends. 1- Pin ended beam: the deflection (w=0) and the bending moment (w''=0) at (x=0) and the deflection (w= ) and the bending moment (w''=0) at (x=L), apply these conditions on equation (39) and solve for the constants to obtain: k1 = Qo/K -----------(40-a) k2 = (k3 cot L + Qo/K) tanh L -----------(40-b) csch L sec L + Qo/K (coth L csch L sec L tanh L) ----------(40-c) k3 = ----------------------------------------------------------------(coth L tan L + tanh L cot L) k4 = 0 -----------(40-d) when ( =0), i.e. no settlement, the resulting equation of the deflection is exactly the same as that given in ref.[2]: cosh x cos (x-L) + cosh (x-L) cos x w = Qo/K [1 ------------------------------------------------ ] -----------(41) (cosh L + cos L) The bending moment equation becomes as: sinh x sin (x-L) + sinh (x-L) sin x 2 M = Qo/(2 ) [--------------------------------------------- ] -----------(42) (cosh L + cos L) Using the (Rayleigh. Ritz) energy method for the functional given in equation (1) and employing the (n) terms approximate function (w = An sin n x/L). A straight forward extremization of the functional with respect to (An) leads to the following approximation: 4QoL4/(n5 5EI) w = ------------------------------------- sin(n x/L) ------------(43) n=1 [1 + KL4/(n4 4EI) P/Po] where n=1,3,5, ..etc. and the bending moment: 4QoL2/(n3 3) M = ------------------------------------- sin(n x/L) ------------(44) n=1 4 4 4 [1 + KL /(n EI) P/Po]
110
Al-Rafidain Engineering
Vol.16
No.3
Aug.
2008
Including the relative end settlement into account, the approximate deflection equation becomes as below: w = An sin n x/L + x/L Following the same procedures of (Rayleigh-Ritz) method, the resulting equations of the deflection and the bending moment are: 2(2Qo K )L4/(n5 5EI) w = ------------------------------------ sin(n x/L) + x/L ------------(45) n=1 4 4 4 [1 + KL /(n EI) P/Po] and 2(2Qo K )L2/(n3 3) M = -------------------------------------- sin(n x/L) ------------(46) n=1 4 4 4 [1 + KL /(n EI) P/Po] The exact solution of pin ended beam using the elementary mechanics of material method [48]: w = QoL4/(24EI) [(x/L)4 2(x/L)2 + (x/L)] + x/L -----------(47) and M = QoL2/2 [(x/L) (x/L)2] ------------(48) The comparison of the results obtained from the exact solution and the (Rayleigh-Ritz) energy method is shown in Table (3) and Figs.(3 and 4). The results show that the energy method using one term series approximate solution (n=1) gives close results to the exact solution and practically accepted, the accuracy of the solution increases when two or more terms are used (n=3, 5, 7 and 9), beyond this limit the accuracy of the solution does not change, thus 5 terms (n=9) is quite enough to give results very close to the exact solution. Figs.(4-8) show the effect of the soil sub-grade reaction value (K) and the relative settlement ( ) on the deflection and bending moment, as shown increasing of (K) decreases the curvature of the beam, thus decreasing the deflection and bending moment. The effect of (K) is clearer in Fig.(9 and 10), the deflection and bending moment exponentially decreased as (K) is increased. The same behavior is observed for beams with settlements (12.5 and 25mm). Fig.(10) shows that negative moments are observed at (K/Ec=0.06%) for settlement (12.5mm) and (K/Ec=0.03%) for settlement (25mm), for this reason the concrete tie beams resting on elastic foundation is reinforced at the top and bottom to resist the positive and negative moments.
2- Fixed ended beam: the deflection (w=0) and the slope (w'=0) at (x=0) and the deflection (w= ) and the slope (w'=0) at (x=L), apply these conditions to equation (39) and solve for the constants to obtain: k1 = Qo/K -----------(49-a) k2 = k4/2 (coth L + cot L) + Qo/K (coth L cot L) -----------(49-b) k3 = k2 -----------(49-c)
2[Qo/K(cosh L cos L1)+ ] sinh L sin L+Qo/K(sinh L cos Lcos L sin L)2 k4 = -------------------------------------------------------------------------------------------------------2(sinh L sin L)2 + (sinh L cos L)2 (cosh L sin L)2
---(49-d)
Using the (Rayleigh-Ritz) energy method for the functional given in equation (1) and employing the (n) terms approximate function which satisfy the boundary conditions of the fixed ended beam: w = An [1 cos(2n x/L)] + [3 (x/L)2 2 (x/L)3] ---------(50) following the same procedures used in pin ended beam, the resulting equation of the deflection and the bending moment are shown below: (2Qo K )L4/(16n4 4EI) 111
112
Al-Rafidain Engineering
Vol.16
No.3
Aug.
2008
113
114
Al-Rafidain Engineering
Vol.16
No.3
Aug.
2008
115
116
Al-Rafidain Engineering
w=
n=1
Vol.16
No.3
Aug.
2008
and (2Qo K )L2/(4n2 2) M = ----------------------------------------- cos(2n x/L) + 6EI /L2(2 x/L 1) ---(52) n=1 [1 + 3KL4/(16n4 4EI) P/Po] The exact solution according to elementary mechanics of material method [4-8]: w = QoL4/(24EI) [(x/L)4 2(x/L)3 + (x/L)2] + [3(x/L)22(x/L)3] ----------(53) and M = QoL2/12 [6(x/L) 6(x/L)2 1] + 6EI /L2(2 x/L 1) ------------(54) Table (4) and Fig.(11) show that the (Rayleigh-Ritz) method gives close results to the exact method, the accuracy of the solution increased with number of terms up to 5 terms, beyond this limit the accuracy of the solution does not change. Rayleigh-Ritz method gives negative moment about (60 75%) of that obtained by exact method while the positive moment is greater by about (21 50%). Figs.(12-17) show the effect of (K) on the deflection and bending moment. The bending moment increased in one end and decreased in the other end due to the end settlement, because the bending moment developed due to the settlement (6EI /L2) is added algebraically to that developed due to the lateral load.
3- Cantilever beam: the deflection (w=0) and the slope (w'=0) at (x=0) and the bending moment (w''=0) and the shear force (w''+k12w' =0) at (x=L), apply these conditions to equation (39) and solve for the constants to obtain: k1 = Qo/K ------------(55-a) k2 = k3 -------------(55-b) k3 = (k1+k4)/2 tan L + (k1k4)/2 tanh L -------------(55-c) 2 2 (cosh L sin L) + (sinh L cos L) ---------------(55-d) k4 = Qo/K ----------------------------------------------------------2(cosh L cos L)2 + (cosh L sin L)2 (sinh L cos L)2 Using the (Rayleigh-Ritz) energy method for the functional given in equation (1) and employing the (n) terms approximate function which satisfy the boundary conditions of the fixed ended beam: w = An [1 cos(n x/2L)] ------------(56) Following the same procedures used in pin ended beam, the resulting equation of the deflection and the bending moment are shown below: 4 4 4 EI) 1QoL /(n w = --------------------------------------- [1 cos(n x/2L)] -----(57) n=1 4 4 4 [1 + 2KL /(n EI) P/Po] and 1QoL2/(4n2 2) M = ------------------------------------------ cos(n x/2L) -----(58) n=1 4 4 4 [1 + 2KL /(n EI) P/Po] where 1 = 32 ( 2)/ and 2 = 16 (3 8)/ The exact solution according to elementary strength of material method [4-8]: ------------(59) w = QoL4/(24EI) [(x/L)4 4(x/L)3 + 6(x/L)2] and M = QoL2/2 [2 (x/L) (x/L)2 1] ------------(60) 117
118
Al-Rafidain Engineering
Vol.16
No.3
Aug.
2008
119
The same conclusion and behavior are obtained as shown in Table (5) and Figs.(18-20). the bending moment obtained by (Rayleigh-Ritz) method is about (59 75%) of that obtained by the exact method. The existence of sub-grade soil supporting the cantilever beam has greater effect than pin ended and fixed ended beams, the reduction of the deflection and the bending moment to be less than (10%).
IV- Elastic buckling of non-prismatic beams on elastic foundation: Rayleigh-Ritz and Galerkin's methods are used to derive the general deflection equation and bending moment including the effect of lateral load Q(x), axial load (P), soil sub-grade reaction and variation of the cross section. Six types of linear and parabolic changing cross section are considered. Galerkin's method is another method for finding approximate solution to a differential equation. This method involves direct use of the differential equations. It does not require the existence of a functional, for this reason the method is used when the functional of the problem is not available, difficult or complicated. For the differential equation (Du = f), where D is the differential operator and u is the solution. If the approximate solution ( u = ai i ), where ( i) are shape functions and (ai) are constants. The integral of the out of balance of the differential equation can be written as below: (Du f) i dV = 0 --------------(61) Both methods (Rayleigh-Ritz and Galerkin's) are applied to the three types of beams (pin ended, fixed ended and cantilever) using the same deflection functions that used for prismatic beams in previous sections. 1- Pin ended beam: 4QoL4/(n5 5EI) w = ----------------------------------- sin(n x/L) ------------(62) n=1 4 4 4 [F + KL /(n EI) P/Po] and 4QoL2/(n3 3) f(x) M = ----------------------------------- sin(n x/L) ------------(63) n=1 4 4 4 [F + KL /(n EI) P/Po] where: F = 1 + a1 + a2 2 + a3 3 = h2/h1 1 h1 and h2 are the smallest and largest depth of the hunched beam respectively. f(x) is the function of changing the moment of inertia shown in Table (6). a1, a2, a3 are constants depend on the shape of the non-prismatic beam. The values of these constants are determined using Rayleigh-Ritz and Galerkin's methods and shown in Table (6). Equations (62 and 63) are similar to those derived for prismatic beams, only one additional factor (F), when (h1=h2, = 0 and F=1) the resulting form of equations (62 and 63) become similar to equations (43 and 44). 2- Fixed ended beam: QoL4/(8n4 4EI) w = -------------------------------------- [(1cos(2n x/L)] n=1 [F + 3KL4/(16n4 4EI) P/Po] and
-----(64)
120
Al-Rafidain Engineering
Vol.16
No.3
Aug.
2008
QoL2/(2n2 2) f(x) M = ----------------------------------------- cos(2n x/L) n=1 [F + 3KL4/(16n4 4EI) P/Po] Equations (64 and 65) are same as equations (51 and 52) when ( =0).
3- Cantilever beam: 4 4 4 EI) 1QoL /(n w = ----------------------------------------- [(1 cos(n x/2L)] n=1 [F + 2KL4/(n4 4EI) P/Po] And 1QoL2/(4n2 2) f(x) M = ---------------------------------------- cos(n x/2L) n=1 [F + 2KL4/(n4 4EI) P/Po] Equations (66 and 67) are similar to equations (57 and 58) when ( =0).
------(65)
-----(66)
-----(67)
Conclusions:
1- The critical load (Pcr) of the beam on elastic foundation is twice the critical load for beam when (K=0) for the same properties for the case ( /2)2 = . 2- The buckling load (Pcr) for fixed ended beam is greater than pin ended beam by ( 3) with the same properties and for the case ( /2)2 = . 3- The buckling load (Pcr) for cantilever beam is about (2/3) of the pin ended beam with the same properties and for the case ( /2)2 = . 4- The (Rayleigh-Ritz) method using one term series approximate solution (n=1) gives close results to the exact method, the accuracy of the solution increased when two or more terms are used up to five terms, beyond this limit the accuracy of the solution does not changed. 5- The increasing of soil sub-grade reaction (K) decreases the curvature of the beam, thus decreasing of the deflection and the bending moment exponentially as (K) is increased. 6- In fixed ended beams, Rayleigh-Ritz method gives negative moment about (60 75%) of that obtained by exact method while the positive moment is greater by about (21 50%). 7- In cantilever beams, the bending moment obtained by (Rayleigh-Ritz) method is about (59 75%) of that obtained by the exact method. The existence (K) reduces the deflection and the bending moment to less than (10%). 8- The general equation obtained for prismatic and non-prismatic beams are exactly similar for all types (pin ended, fixed ended and cantilever). The factor (F) is to be added to the equations which represent the effect of the variation of the cross section. This factor is unity in prismatic beams.
References: 1- Timoshenko, S.P. and Gere, J.M., "Theory of Elastic Stability", Second edition, McGraw-Hill International Book Company, 1985. 2- Dym, C.L. and Shames, I. H., "Solid Mechanics A Variational Approach", McGrawHill Kogakusha, Ltd., Tokyo, Japan,1970. 3- Washezu, K., "Variational Methods in Elasticity and Plasticity", 2nd edition, Pergamon Press, Oxford, England, 1975. 121
4- Gere, J. M. and Timoshenko, S.P., "Mechanics of Materials", 2nd edition, CBS Publishers and Distributors, New Delhi, India, 2004. 5- Gere, J. M., ""Mechanics of Materials", Brooks/Cole, 5th edition, USA, 2001. 6- Morrow, H.W., "Statics and Strength of Materials", 3rd edition, Prentice Hall, New Jersi, USA, 1998. 7- Spiegel, L. and Limbrunner, G.F., "Applied Statics and Strength of Materials", 4th edition, Pearson Prentice Hall, New Jersi, USA, 2004. 8- Raymond, P., "Solid Mechanics in Engineering", John Wiley and Sons, Ltd., England, 2001.
122