Variational Principles for Bending and Vibration of Partially Composite Timoshenko Beams

Abstract

Variational principles are established for the partially composite Timoshenko beam using the semi-inverse method. The principles are derived directly from governing differential equations for bending and vibration of the beam considered. It is concluded that the semi-inverse method is a powerful tool for searching for variational principles directly from the governing equations. Comparison between our results and the results reported in literature is given.

(ProQuest: ... denotes non-US-ASCII text omitted.)

Academic Editor:Robertt A. Valente

Division of Mechanics, Mechanical Engineering Department, Yildiz Technical University (YTU), Yildiz, Besiktas, 34349 Istanbul, Turkey

Received 19 April 2014; Accepted 29 June 2014; 13 July 2014

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

Composite beams composed of different elastic materials have been widely used in many engineering applications. The individual beam components of the composite beam are combined by using the shear connectors. Therefore, the overall behavior of the composite beam depends on the stiffness of connectors. Connector having infinite stiffness eliminates any interlayer shear slip between the individual beam components, which leads to the full interaction connection. However, the stiffness of connector has a finite value and the interlayer slip between the individual components occurs. This type of connection is called partial-interaction connection. Therefore, analysis of the partial-interaction composite beams requires the consideration of the interlayer slip between the beam components. The Euler-Bernoulli beam theory has been extensively used in bending, vibration, and buckling analyses. Ecsedi and Baksa [1] analyzed the static behavior of elastic two-layer beams with interlayer slip and developed closed-form solutions for displacements and interlayer slips. Girhammar and Pan [2] presented general solutions for the deflection and internal actions for partially composite Euler-Bernoulli beams and beam-columns. Ranzi et al. [3] presented an analytical formulation for the analysis of two-layered composite beams with longitudinal and vertical partial-interaction. Their formulation is based on the principle of virtual work expressed in terms of the vertical and axial displacements of the two layers. The model was presented in both its weak and its strong forms. Xu and Wu [4] developed a new plane stress model of composite beams with interlayer slips using the one-dimensional theory. They concluded that the shear force produced by the shear connectors increases with the increase in rigidity of shear connectors.

However, the effect of transverse shear deformation was neglected in the Euler-Bernoulli beam theory. When the beam is thick, the effect of shear deformation becomes significant and cannot be neglected for a valid analysis. The most widely used and fundamentally simpler theory was developed by Timoshenko [5]. Sousa and da Silva [6] studied the behavior of the general case of multilayered composite beams with interlayer slip, under Euler-Bernoulli as well as Timoshenko beam theory (TBT) assumptions. Xu and Wang [7] formulated the principle of virtual work and reciprocal theorem of work for the partial-interaction composite beams using the kinematic assumptions of Timoshenko's beam theory. The variational principles for the frequency of free vibration and critical load of buckling were also deduced. Xu and Wang [8] derived the relationships of solutions between single-span Euler-Bernoulli and Timoshenko partial-interaction composite beams.

Variational formulations provide the basis for a number of approximate and numerical methods. Recently, two significant variational methods are proposed by He; one is the semi-inverse method [9, 10] and the other method is the variational iteration method [11]. The semi-inverse method is used to establish variational principles directly from the governing differential equations. His second method, the variational iteration method, depends on constructing a correction functional by a general Lagrange multiplier. Then, the optimal value of the Lagrange multiplier is identified by using the stationary conditions [12, 13]. However, in the semi-inverse method, the term involving the Lagrange multiplier is replaced by an unknown function F . The semi-inverse method eliminates two important variational crises; one is that the Lagrange multiplier is equal to zero and the other crisis is that making the Lagrangian stationary leads to only some parts of Euler equations. In this study, we will apply the semi-inverse method to establish variational principles directly from the governing differential equations defining the bending and vibration of Timoshenko composite beam with partial-interaction. The variational formulations were obtained by following the rules of the calculus of variations.

2. Timoshenko Composite Beam with Partial-Interaction

Before applying the semi-inverse method, the problem is briefly discussed in Figure 1. Figure 1 shows a partial-interaction composite beam that is composed of two-layer beams with different materials.

Figure 1: A two-layer composite beam.

[figure omitted; refer to PDF]

In Figure 1, _Ei , _Gi , _Ai , _Ii , and _ρi (i=1,2) denote the elasticity modulus, shear modulus, cross-sectional area, and moment of inertia of two beam components, respectively. L is the beam length, H is the beam height, and h is the distance between the centroids of two beam sections. q and m denote the distributed load and distributed bending moment, respectively. As seen in Figure 1, shear connectors are used to connect the beam members of the composite beam. Figure 2 shows geometrical relationship among the interlayer slip, rotary angle (the rotation of the cross section), and longitudinal displacements.

Figure 2: Kinematic model of a two-layer partially composite Timoshenko beam [7].

[figure omitted; refer to PDF]

In Figure 2, ψ is the rotary angle and _us is the interlayer slip between two beam layers. _u1 and _u2 are the longitudinal displacements at the centroids of beams 1 and 2, respectively. From Figure 2, the kinematic relationship among the interlayer slip, rotary angle, and longitudinal displacements can be written as follows: [figure omitted; refer to PDF] The bending moment, shear force, and interlayer shear force are given, respectively, as [7, 8]: [figure omitted; refer to PDF] where w denotes the deflection of the composite beam in the z -direction (see Figure 1) and _ks denotes the rigidity of the shear connectors. The other quantities used in (2a), (2b), and (2c) are defined as follows: [figure omitted; refer to PDF] in which _k1 and _k2 are the shear correction factors of the Timoshenko beam. D¯ is the flexural stiffness of the composite beam in full interaction, D is the flexural stiffness of the composite beam without shear connection, EA¯ is the effective axial stiffness, and C is the shear rigidity of the whole cross sections. Deflection of the composite beam is then obtained using the relation below: [figure omitted; refer to PDF] where _w0 is the deflection of the full interaction composite beam. _wslip and _wshear are the additional deflections due to the interlayer slip and transverse shear deformation, respectively.

In the next sections, using the semi-inverse method, we will illustrate how to establish variational principles directly from the governing differential equations for bending and vibration of the partially composite Timoshenko beams.

3. Derivation of Variational Principle for Bending of the Composite Beam

Consider the governing differential equations for bending of the partial-interaction composite Timoshenko beam under uniformly distributed load and bending moment [7] [figure omitted; refer to PDF] Using the semi-inverse method, a trial variational principle can be constructed as follows [9, 10]: [figure omitted; refer to PDF] where L is a trial Lagrangian. There are many approaches for constructing the trial Lagrangian; see [14-17]. We search for such a trial Lagrangian, so that its trial Euler equation gives one of the governing equations, say (5a). Referring to (5a), an energy-like trial Lagrangian can be constructed as follows: [figure omitted; refer to PDF] where _F1 is an unknown function of w and/or its derivatives. The advantage of the above trial Lagrangian lies in the fact that the stationary condition with respect to ψ results in (5a). Now by making (7) stationary with respect to w , one can get the following trial Euler equation for δw : [figure omitted; refer to PDF] where the operator δ is called a variational operator and δw is the first order variation of w . δ_F1 /δw is called He's variational derivative with respect to w , which is defined as [figure omitted; refer to PDF] We search for such an _F1 so that (8) is equivalent to (5b). From (9), the unknown function _F1 can be determined as [figure omitted; refer to PDF] By adding the above relation, the trial Lagrangian can be renewed as follows: [figure omitted; refer to PDF] It can be easily proved that the stationary condition of the above Lagrangian with respect to w satisfies (5b). In (11), _F2 is a newly introduced undetermined function of _us and/or its derivatives and is free from the variables ψ and w . Making the new trial Lagrangian (11) stationary with respect to _us results in the relation below: [figure omitted; refer to PDF] which is the last trial Euler equation. The second term on the left is the variational derivative with respect to _us and reads [figure omitted; refer to PDF] from which the unknown _F2 can be determined as [figure omitted; refer to PDF] Substituting _F2 into (11) and rearranging lead to the necessary variational principle as [figure omitted; refer to PDF] which is the total potential energy of partial-interaction composite Timoshenko beam subjected to uniformly distributed load and bending moment (see [7]) and yields the minimum potential energy principle by letting δJ=0 .

Proof.

Making the above functional (15) stationary with respect to ψ , w , and _us , the Euler equations turn out to be (5a)-(5c), respectively.

The Ritz method can be used to obtain an approximate analytical solution of the problem. We can write the one-term trial functions which satisfy the boundary conditions as [figure omitted; refer to PDF] where _w0 , _ψ0 , and _u0 are unknown constants, which can be determined from the following stationary conditions: [figure omitted; refer to PDF] By solving the system of (17) simultaneously, the unknown constants can then be obtained.

4. Derivation of Variational Principle for Free Vibration of the Composite Beam

Differential equations of motion for partial-interaction composite members under uniformly distributed load and bending moment can be written as [7]: [figure omitted; refer to PDF] where t denotes time, ρ_A0 =_ρ1A1 +_ρ2A2 , and ρ_I0 =_ρ1I1 +_ρ2I2 . We can construct the following trial variational principle using the semi-inverse method [9, 10]: [figure omitted; refer to PDF] Similarly, referring to (18a) and making some modifications so that the stationary condition with respect to ψ can identify (18a) lead us to the following trial Lagrangian: [figure omitted; refer to PDF] with _F3 being an unknown function of w and/or its derivatives. As can be seen easily, the stationary condition of the above Lagrangian with respect to ψ results in (18a). Now making (20) stationary with respect to w , we obtain the following trial Euler equation: [figure omitted; refer to PDF] where δ_F3 /δw is defined as [figure omitted; refer to PDF] From the above relation, we can identify _F3 in the form [figure omitted; refer to PDF] Then, the Lagrangian (20) is further updated as follows: [figure omitted; refer to PDF] It is obvious that making the renewed trial functional stationary with respect to w satisfies (18b). In (24), _F4 is a new undetermined function of _us and/or its derivatives. It must be noted that (18c) has the same form as (5c). Therefore, by following the same steps as before (see (12)-(13)), it is easily seen that _F4 =_F2 . Finally, we can easily arrive at the required variational principle: [figure omitted; refer to PDF] The above functional is the same as that reported in [7] and yields Hamilton's principle by letting δJ=0 . The fifth term inside the braces is the kinetic energyof the beam componentsand reads [figure omitted; refer to PDF]

Proof.

Making the above functional (25) stationary with respect to ψ , w , and _us , the Euler equations correspond to (18a)-(18c), respectively.

By following the same procedures performed for beam bending, the approximate solutions are obtained conveniently for beam vibrating by the Ritz method.

5. Conclusion

We used the semi-inverse method to establish a set of variational principles directly from governing differential equations. By following the rules of the calculus of variations, we obtained necessary variational principles for bending and vibration of the Timoshenko composite beam with partial-interaction. The obtained variational principles have been compared with those reported in literature and proved to be correct. It is concluded that the semi-inverse method is a powerful tool for searching for variational principles directly from the governing equations. Moreover, introducing an unknown function instead of a Lagrange multiplier, additional variational principles can also be written by constraining the trial Lagrangian with the different boundary conditions, which may facilitate the implementation of complicated boundary conditions. The direct variational method such as the Ritz method can be used to obtain the approximate solutions of the problem.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

References

[1] I. Ecsedi, A. Baksa, "Static analysis of composite beams with weak shear connection," Applied Mathematical Modelling , vol. 35, no. 4, pp. 1739-1750, 2011.

[2] U. A. Girhammar, D. H. Pan, "Exact static analysis of partially composite beams and beam-columns," International Journal of Mechanical Sciences , vol. 49, no. 2, pp. 239-255, 2007.

[3] G. Ranzi, F. Gara, P. Ansourian, "General method of analysis for composite beams with longitudinal and transverse partial interaction," Computers & Structures , vol. 84, no. 31-32, pp. 2373-2384, 2006.

[4] R. Xu, Y. Wu, "Two-dimensional analytical solutions of simply supported composite beams with interlayer slips," International Journal of Solids and Structures , vol. 44, no. 1, pp. 165-175, 2007.

[5] S. P. Timoshenko, "On the correction for shear of the differential equation for transverse vibrations of prismatic bars," Philosophical Magazine , vol. 41, pp. 744-746, 1921.

[6] J. B. M. Sousa, A. R. da Silva, "Analytical and numerical analysis of multilayered beams with interlayer slip," Engineering Structures , vol. 32, no. 6, pp. 1671-1680, 2010.

[7] R. Xu, G. Wang, "Variational principle of partial-interaction composite beams using Timoshenkos beam theory," International Journal of Mechanical Sciences , vol. 60, no. 1, pp. 72-83, 2012.

[8] R. Xu, G. Wang, "Bending solutions of the Timoshenko partial-interaction composite beams using Euler-Bernoulli solutions," Journal of Engineering Mechanics , vol. 139, no. 12, pp. 1881-1885, 2013.

[9] J. H. He, "Semi-inverse method and generalized variational principles with multi-variables in elasticity," Applied Mathematics and Mechanics. Yingyong Shuxue he Lixue , vol. 21, no. 7, pp. 797-808, 2000.

[10] J. H. He, "Derivation of generalized variational principles without using lagrange multipliers, part I: applications to fluid mechanics," Mechanics, Automatic Control and Robotics , vol. 3, no. 11, pp. 5-12, 2001.

[11] J. H. He, "Variational iteration method--a kind of non-linear analytical technique: some examples," The International Journal of Non-Linear Mechanics , vol. 34, no. 4, pp. 699-708, 1999.

[12] H. Ozer, "Application of the variational iteration method to the boundary value problems with jump discontinuities arising in solid mechanics," International Journal of Nonlinear Sciences and Numerical Simulation , vol. 8, no. 4, pp. 513-518, 2007.

[13] H. Ozer, "Application of the variational iteration method to thin circular plates," International Journal of Nonlinear Sciences and Numerical Simulation , vol. 9, no. 1, pp. 25-30, 2008.

[14] J. H. He, "Variational principles for some nonlinear partial differential equations with variable coefficients," Chaos, Solitons & Fractals , vol. 19, no. 4, pp. 847-851, 2004.

[15] J. He, "A Lagrangian for von Karman equations of large deflection problem of thin circular plate," Applied Mathematics and Computation , vol. 143, no. 2-3, pp. 543-549, 2003.

[16] J. H. He, "On the semi-inverse method and variational principle," Thermal Science , vol. 17, no. 5, pp. 1565-1568, 2013.

[17] I. Kucuk, I. S. Sadek, S. Adali, "Variational principles for multiwalled carbon nanotubes undergoing vibrations based on nonlocal timoshenko beam theory," Journal of Nanomaterials , vol. 2010, 2010.