Research Article | Open Access
Variational Principles for Bending and Vibration of Partially Composite Timoshenko Beams
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.
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  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  presented general solutions for the deflection and internal actions for partially composite Euler-Bernoulli beams and beam-columns. Ranzi et al.  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  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 . Sousa and da Silva  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  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  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 . 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 . 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.
In Figure 1, , , , , and denote the elasticity modulus, shear modulus, cross-sectional area, and moment of inertia of two beam components, respectively. is the beam length, is the beam height, and is the distance between the centroids of two beam sections. and 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.
In Figure 2, is the rotary angle and is the interlayer slip between two beam layers. and 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: The bending moment, shear force, and interlayer shear force are given, respectively, as [7, 8]:where denotes the deflection of the composite beam in the -direction (see Figure 1) and denotes the rigidity of the shear connectors. The other quantities used in (2a), (2b), and (2c) are defined as follows: in which and are the shear correction factors of the Timoshenko beam. is the flexural stiffness of the composite beam in full interaction, is the flexural stiffness of the composite beam without shear connection, is the effective axial stiffness, and is the shear rigidity of the whole cross sections. Deflection of the composite beam is then obtained using the relation below: where is the deflection of the full interaction composite beam. and 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 Using the semi-inverse method, a trial variational principle can be constructed as follows [9, 10]: where 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: where is an unknown function of 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 , one can get the following trial Euler equation for : where the operator is called a variational operator and is the first order variation of . is called He’s variational derivative with respect to , which is defined as We search for such an so that (8) is equivalent to (5b). From (9), the unknown function can be determined as By adding the above relation, the trial Lagrangian can be renewed as follows: It can be easily proved that the stationary condition of the above Lagrangian with respect to satisfies (5b). In (11), is a newly introduced undetermined function of and/or its derivatives and is free from the variables and . Making the new trial Lagrangian (11) stationary with respect to results in the relation below: which is the last trial Euler equation. The second term on the left is the variational derivative with respect to and reads from which the unknown can be determined as Substituting into (11) and rearranging lead to the necessary variational principle as which is the total potential energy of partial-interaction composite Timoshenko beam subjected to uniformly distributed load and bending moment (see ) and yields the minimum potential energy principle by letting .
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 where , , and are unknown constants, which can be determined from the following stationary conditions: 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 :where denotes time, , and . We can construct the following trial variational principle using the semi-inverse method [9, 10]: 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: with being an unknown function of 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 , we obtain the following trial Euler equation: where is defined as From the above relation, we can identify in the form Then, the Lagrangian (20) is further updated as follows: It is obvious that making the renewed trial functional stationary with respect to satisfies (18b). In (24), is a new undetermined function of 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 . Finally, we can easily arrive at the required variational principle: The above functional is the same as that reported in  and yields Hamilton’s principle by letting . The fifth term inside the braces is the kinetic energyof the beam componentsand reads
By following the same procedures performed for beam bending, the approximate solutions are obtained conveniently for beam vibrating by the Ritz method.
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.
- I. Ecsedi and A. Baksa, “Static analysis of composite beams with weak shear connection,” Applied Mathematical Modelling, vol. 35, no. 4, pp. 1739–1750, 2011.
- U. A. Girhammar and 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.
- G. Ranzi, F. Gara, and 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.
- R. Xu and 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.
- 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.
- J. B. M. Sousa and A. R. da Silva, “Analytical and numerical analysis of multilayered beams with interlayer slip,” Engineering Structures, vol. 32, no. 6, pp. 1671–1680, 2010.
- R. Xu and 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.
- R. Xu and 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- 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.
- J. H. He, “On the semi-inverse method and variational principle,” Thermal Science, vol. 17, no. 5, pp. 1565–1568, 2013.
- I. Kucuk, I. S. Sadek, and S. Adali, “Variational principles for multiwalled carbon nanotubes undergoing vibrations based on nonlocal timoshenko beam theory,” Journal of Nanomaterials, vol. 2010, Article ID 461252, 7 pages, 2010.
Copyright © 2014 Halil Özer. 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.