Analytical Solutions of Functionally Graded Curved Beams under an Arbitrarily Directed Single Force

Zhang, Liangliang; Shang, Lange; Gao, Yang

doi:https://doi.org/10.1155/2019/4925379

Mathematical Problems in Engineering

On this page

Abstract Introduction Discussion Data Availability Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2019 | Article ID 4925379 | https://doi.org/10.1155/2019/4925379

Analytical Solutions of Functionally Graded Curved Beams under an Arbitrarily Directed Single Force

Liangliang Zhang,^1,2Lange Shang,¹and Yang Gao¹

Academic Editor: Gaetano Giunta

Received29 Nov 2018

Accepted07 Feb 2019

Published27 Feb 2019

Abstract

A functionally graded curved beam subjected to a shear tension force as well as a concentrated force at the free end is solved based on the inverse method, and a general two-dimensional solution is presented. The explicit expressions are derived by assuming that the elastic properties within curved beams vary in the radial direction according to a power law, i.e., E = E₀rⁿ, but are constant across the depth. After degenerating it into the isotropic homogeneous elastic cases, the results are in good consistency with existing analytical solutions. The stresses and displacements are firstly observed in different forms in terms of the different power function exponent . These results will be useful as a guide for designing devices or as benchmark to assess other approximate methodologies.

1. Introduction

Functionally graded materials (FGMs) are inhomogeneous materials in which the volume fractions of two or more materials are varied continuously as a function of position along certain direction(s) of the structure to achieve required properties. Structures made of FGMs offer great application potential in engineering and industrial fields such as aerospace, automotive, petrochemical, and civil engineering fields [1–5]. Due to the increasing interest in FGMs in recent years, lots of studies have been performed on the static and dynamic behavior of functionally graded structures [6–12].

Curved beams are ones of frequently used components in frame structures, whose mechanical properties have been obtained by many researchers [13–16]. Many approximate methods, such as classical beam theory [17], Timoshenko’s beam theory [18–21], and the high order beam theory [22], have been used to study the static displacement and stress distribution problems. For curved beams made of homogeneous materials, Timoshenko and Goodier [23] gave explicit solutions for curved beam under pure bending and curved cantilever beam under free end loads by using the Airy stress function method. Kilic and Aktas [24] got the elasticity solution for curved cantilever beams subjected to a single arbitrarily directed force at the free end. Bagci [25] studied the curved beams and rings of exponential thickness by using polar elasticity and state of plane stress. For functionally graded curved beams, Lekhnitskii [26] derived the solutions with elastic modulus being a product of a periodical function of the tangential coordinate and a power law function along the radial coordinate. Dryden [27] investigated the stress distribution of a functionally graded curved beam under pure bending where Young’s modulus was in a more general form with proportional to and Poisson’s ratio was held constant. However, in his study, only pure bending problem is considered. Based on the Airy stress function method, Wang and Liu [28] extended Dryden’s study [27] to the ones of the orthotropic functionally graded curved beams subjected to a uniform load on the outer surface and having various constraints or/and loads at ends. In his research, the beams are composed of n layers and each layer is orthotropic and functionally graded. Pydah and Sabale [29] analyzed the flexure of bidirectional functionally graded circular beams by using the kinematical assumptions of the Euler-Bernoulli theory. In the paper, the Young’s modulus was graded along both radial direction and tangential direction. However, due to the mathematical complexity, they failed to carry out the explicit formulation of the mechanical field. They also gave out a review of literature for functionally graded curved beams. By using the laminated model and the Taylor series expansion method, Shi el al. [30, 31] obtained the exact solutions of bending problem of a piezoelectric layered curved beam. Besides analytical studies, finite element methods play a significant role in the studied of functionally graded curved beams. De Miguel et al. [32] presented a new class of refined curved beam elements for the accurate stress analysis of composite structures. In their study, the kinematics of the beam element based on layer-wise and component-wise models and the locking-free curved beams elements based on the mixed interpolation of tensorial components method are provided. De Pietro et al. [33] investigated the mechanical behavior of slender and thick curved beams based on Navier-type analytical solution and a Carrera’s unified formulation based weak form finite element solution.

In this paper, the elastic analysis for a functionally graded curved beam with Young’s modulus varies as E(r) = E₀rⁿ is carried out. The variation of exponent is considered, and the explicit expressions are derived based on the Airy stress function method. The stresses and the displacements are observed in different forms in terms of the different power function exponent . These results will be useful as a guide for designing functionally graded devices or as benchmark to assess other approximate methodologies.

2. Basic Equations

Consider a functionally graded curved beam with unit depth referred to a polar coordinate system (r, θ). As shown in Figure 1, the curved beam is bounded by inner and outer radii equal to and with the origin at the center of this curved beam and let θ = 0 and and θ = α be the free end and the fixed end, respectively. Industrial loading applied on a curved beam could be decomposed to a moment and an arbitrarily directed concentrated force. As for pure bending problem, it was studied by Dryden [27]. In this paper, an arbitrarily directed single force was considered, which can be decomposed to a shear tension force and a normal concentrated force applied on the free end as shown in Figure 1.

To start, the constitutive equations arewhere E(r) = E₀rⁿ is the varying function of the Young’s modulus, is the Young’s modulus at r = a of the curved beam, and Poisson’s ratio is constant. Also, in (1), the shear modulus can be considered by letting E(r) = 2(1+)G(r). The power function is either monotonically increasing or decreasing depending upon the index of the exponent . As in previous treatments dealing with the analysis of inhomogeneous solids, Poisson’s ratio is considered to be constant. Since the main interest is to find how the stiffness variation influences the flexural stress and displacement, this is not a serious deficiency.

The strain-displacement equations are

The strain compatibility equations are

With a stress function φ, the stress expressions are

The boundary conditions are as follows:

The objective of this work is to present analytical solutions of functionally graded curved beams, assuming that the elastic properties vary in the radial direction according to a power law, i.e., E = E₀rⁿ, but are constant across the depth. In Sections 3 and 4, a shear tension force and a normal concentrated force applied on the free end of functionally graded curved beams are analytically formulated in terms of Airy stress function. In Section 5, numerical results and discussion are provided for a specific example.

3. Problem with a Shear Tension Force Applied

Under a shear tension force Q, the stress function will relate to both and θ. Here we assume that the stress function is φ(r, θ) = R(r) with the boundary conditions (5), where Q ≠ 0 and P = 0. Substituting φ(r, θ) into (1) then into (3), we obtainwhich satisfies the Euler equation with solution of R = Dr^k. Substituting the solution into (6), we obtainThen the characteristic roots of the aforementioned equation are = 1 and = 1 + n, and the other two roots and must satisfy Due to the values of which are –1 < < 1/2, the roots and are different from each other. Considering the particular solutions of n, this problem can be divided into two cases, namely, Case 1: n ≠ –4/(1–) and Case 2: n = –4/(1–). In Case 1, all roots are distinct to each other; in the Case 2, there are two pairs of repeated roots.

Case 1 (n ≠ –4/(1–)). In this case, the stress function can be assumed as where can be ignored for it has no effect on the stresses and displacements. In this paper, we take i, j = 1, 2, 3. The summation convention has been used throughout this paper: the Einstein summation over repeated lower case indices is applied, while upper case indices take on the same numbers as the corresponding lower case ones but are not summed. By substituting φ(r, θ) into (4), the stresses can be derived From the boundary conditions (5), we determinedwhere ε_ijk is Levi-Civita symbol. The Levi-Civita symbol is defined asSubstituting (10) into (1) and after integration, we can get the displacements aswhere f(θ) and (r) are expressed as where , , and are arbitrary constants which can be determined from the constraint conditions of the curved beams. The corresponding displacements of f(θ) and (r) are and , as follows:In (15), it is shown that there are rigid body displacements in and , which is in and in . These two items are related to the translation and rotation of the rigid body, which can be balanced by applying a corresponding constraint. Thus, it can be ignored and then the displacements of the curved beam can be carried out as

Case 2 (n = –4/(1–)). In this case, we assumewhere are constants to be determined. Similarly, the stresses and displacements are summarized as follows:If we set n = 0, then these results of Cases 1 and 2 can be degenerated to isotropic homogeneous elastic cases. According to (6), the characteristic roots , , , and are 0, 0, –1, and 3, respectively. As a result, we can assume the stress functions as Thus in the similar way of Cases 1 and 2 by using this stress function, the solutions deduced here coincide well with existing solutions [17].

Assume the curved beam is expanded to be a curved ring by letting α = 2π as shown in Figure 2. Considering in (16) and (19), there is an edge dislocation δ which is caused by the curved ring itself. To avoid this edge dislocation, we can add an additional shear tension force on both the inner and the outer faces. By substituting into (11) and (20), the additional shear tension force can be obtained for the two cases, respectively,

4. Problem with a Concentrated Force Applied

This problem has the boundary conditions (5), where P ≠ 0, Q = 0. In the similar way of isotropic homogenous curved beams subjected to a concentrated force, we can assume the stress functions as φ(r, θ) = (r) + R(r), in which (r) is only related to . As a result, it is obvious to observe that (r) is the stress function of the pure bending problem which has been solved in our previous work [34] and R(r) is the stress function of the problem solved in Section 3. For pure bending problem, the characteristic equation is where the characteristic roots are = 0 and = 2 + n, and the other two roots and satisfyIn virtue of –1 < < 1/2, the roots and are different from each other. Considering the particular solutions of n, this problem can be divided into four cases, respectively, Case 3: n ≠ 1/, n ≠ –4/(1–), and n ≠ –2; all characteristic roots are distinct and ≠ 1 + n; Case 4: n = 1/; all characteristic roots are distinct; without loss of generality, we take = 1 + n and = 1; Case 5: n = –4/(1–); the characteristic roots are distinct and the roots have two pairs of repeated roots; Case 6: n = –2, the characteristic roots have a repeated root and the roots are distinct to each other. Combining with the cases which are given in Section 3, we can solve this problem in four different cases.

Case 3 (n ≠ 1/, n ≠ –4/(1–) and n ≠ –2). In this case, all the characteristic roots are distinct to each other and ≠ 1 + n. After the same manipulation as in Case 1, assume the stress functions as The stresses and the displacements without rigid body displacement in Case 3 can be stated as follows:

Case 4 (n = 1/). In this case, all the characteristic roots are distinct to each other. Without loss of generality, we take = 1 + n and = 1. The stress function has the same expression as in Case 3. We can get that the stresses are the same as in (26), and the displacements in a similar manipulation areC_i and D_i are the same as in (28).

Case 5 (n = –4/(1–ν)). In this case, the characteristic roots are distinct to each other but the roots have two pairs of repeated roots. As a result, we assume Then the stresses and the displacements without rigid body displacement can be obtained as

Case 6 (n = –2). In this case, the characteristic roots have a repeated root and the roots are distinct to each other. Thus, we assumeThen the stresses and the displacements without rigid body displacement can be obtained as

5. Numerical Results and Discussion

In this section, numerical studies and discussion will be presented. Before that, it is appropriate to verify the accuracy of the formulation. The validation studies are investigated by comparing the proposed solutions with finite element results for elastic and functionally graded curved beams. Specifically, we take SiC/C FGM as a specimen, whose is 400 GPa. The beam is set up as the same as Figure 1 (a = 6 mm, b = 12 mm) subjected to a tension force P = 100N at the free end (θ = 0). For convenience, we take Poisson’s ratio as 0.2. To conduct the finite element solution, commercial software ANSYS is used with the beam element BEAM188 and a mesh with 50 elements along the axis. The stress σ_θ at the section θ = π/3 is collected to make the comparisons which are shown in Figure 3 with elastic curved beams n = 0 (a) and functionally graded curved beams n = 2 (b). To avoid the size effect, we take ρ = (r – a)/a. It is worth mentioning that the present solutions show great agreement with the finite element results.

(a) ρ = (r – a)/a

(b) ρ = (r – a)/a

When it comes to numerical simulation, we consider the particular solutions of n; the graded exponent can be taken as –5, –2, 0, 2, and 5. Figure 4 shows the variation of displacements and u_θ and stresses , σ_θ, and τ_rθ at the section θ = π/3 with coordinate . To avoid the size effect, we take ρ = (r – a)/a.

(a)

(b) ρ = (r – a)/a

(c) ρ = (r – a)/a

(d) ρ = (r – a)/a

(e) ρ = (r – a)/a

It is observed that, in Figures 4(c) and 4(d), the displacement stays nearly uniform across the radius while the displacement u_θ varies nearly linearly with coordinate ρ. The distribution of stress σ_θ as shown in Figure 4(e) is nonlinear along the radius direction for FGM (n ≠ 0), different from the nearly linear distribution for a homogeneous material (n = 0). The variation of stress which is exhibited in Figure 4(a) is very similar to τ_rθ, except that the values of whose are opposite. It is worthy to notice that the values of σ_θ are much bigger compared to and τ_rθ at the inner and outer faces. From Figures 4(a) and 4(b), it is expected to observe that the stresses and τ_rθ increase from zero at the inner face (ρ = 0) to the peak and then decrease to zero at the outer face (ρ = 1). And the maximum values of and τ_rθ are shifted from the inner face to the outer face. The stress and displacement field of the curved beam is greatly influenced by graded exponent n, which means that the curved beam can be optimally designed according to given working conditions by adjusting the graded material properties.

6. Concluding Remarks

The analytical elasticity solutions were derived for functionally graded curved beams by means of the inverse method. The mechanical properties of the material were assumed to be the same across the depth. Then a general two-dimensional solution was obtained for functionally graded curved beams subjected to a shear tension force as well as a concentrated force at the free end and was validated with finite element results, which showed great agreement with the proposed formulation. In this study, we observed the power function exponent had strong effect on stresses and displacements. The expressions of the solutions were in different forms when the exponent was assigned a characteristic value. The obtained solution is valid for graded variations as a power law, i.e., E = E₀rⁿ, so it could serve as a guide for designing or as a benchmark to assess other approximate methodologies.

Data Availability

The figure data (Figures 3(a)–3(e)) used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The work is supported by the National Natural Science Foundation of China (Nos. 11472299 and 51704015), China Agricultural University Education Foundation (No. 1101-240001).

References

K. H. Lee, Y. J. Lee, and S. B. Cho, “Characteristics of a transiently propagating crack in functionally graded materials,” Journal of Mechanical Science and Technology, vol. 23, no. 5, pp. 1306–1322, 2009.
View at: Publisher Site | Google Scholar
M. Z. Nejad and G. H. Rahimi, “Deformations and stresses in rotating FGM pressurized thick hollow cylinder under thermal load,” Scientific Research and Essays, vol. 4, no. 3, pp. 131–140, 2008.
View at: Google Scholar
R. S. Khabbaz, B. D. Manshadi, and A. Abedian, “Nonlinear analysis of FGM plates under pressure loads using the higher-order shear deformation theories,” Composite Structures, vol. 89, no. 3, pp. 333–344, 2009.
View at: Publisher Site | Google Scholar
X. Jin, L. Wu, Y. Sun, and L. Guo, “Microstructure and mechanical properties of ZrO2/NiCr functionally graded materials,” Materials Science and Engineering: A Structural Materials: Properties, Microstructure and Processing, vol. 509, no. 1-2, pp. 63–68, 2009.
View at: Publisher Site | Google Scholar
G. M. Genin and V. Birman, “Micromechanics and structural response of functionally graded, particulate-matrix, fiber-reinforced composites,” International Journal of Solids and Structures, vol. 46, no. 10, pp. 2136–2150, 2009.
View at: Publisher Site | Google Scholar
B. V. Sankar, “An elasticity solution for functionally graded beams,” Composites Science and Technology, vol. 61, no. 5, pp. 689–696, 2001.
View at: Publisher Site | Google Scholar
N. Tutuncu and O. Murat, “Exact solutions for stresses in functionally graded pressure vessels,” Composites Part B: Engineering, vol. 32, no. 8, pp. 683–686, 2001.
View at: Google Scholar
M. Jabbari, S. Sohrabpour, and M. R. Eslami, “Mechanical and thermal stresses in a functionally graded hollow cylinder due to radially symmetric loads,” International Journal of Pressure Vessels and Piping, vol. 79, no. 7, pp. 493–497, 2002.
View at: Publisher Site | Google Scholar
J. Ying, C. F. Lü, and W. Q. Chen, “Two-dimensional elasticity solutions for functionally graded beams resting on elastic foundations,” Composite Structures, vol. 84, no. 3, pp. 209–219, 2008.
View at: Publisher Site | Google Scholar
X.-F. Li, “A unified approach for analyzing static and dynamic behaviors of functionally graded Timoshenko and Euler-Bernoulli beams,” Journal of Sound and Vibration, vol. 318, no. 4-5, pp. 1210–1229, 2008.
View at: Publisher Site | Google Scholar
M. Şimşek and T. Kocatürk, “Free and forced vibration of a functionally graded beam subjected to a concentrated moving harmonic load,” Composite Structures, vol. 90, no. 4, pp. 465–473, 2009.
View at: Publisher Site | Google Scholar
F. Z. Jouneghani, R. Dimitri, M. Bacciocchi, and F. Tornabene, “Free vibration analysis of functionally graded porous doubly-curved shells based on the first-order shear deformation theory,” Applied Sciences, vol. 7, no. 12, p. 1252, 2017.
View at: Google Scholar
J. Sanchez-Hubert and E. Sanchez Palencia, “Statics of curved rods on account of torsion and flexion,” European Journal of Mechanics - A/Solids, vol. 18, no. 3, pp. 365–390, 1999.
View at: Publisher Site | Google Scholar | MathSciNet
N. Hu, B. Hu, B. Yan, H. Fukunaga, and H. Sekine, “Two kinds of C0-type elements for buckling analysis of thin-walled curved beams,” Computer Methods Applied Mechanics and Engineering, vol. 171, no. 1-2, pp. 87–108, 1999.
View at: Publisher Site | Google Scholar
P. Dadras, “Plane strain elastic-plastic bending of a strain-hardening curved beam,” International Journal of Mechanical Sciences, vol. 43, no. 1, pp. 39–56, 2001.
View at: Publisher Site | Google Scholar
C. P. Filipich and M. T. Piovan, “The dynamics of thick curved beams constructed with functionally graded materials,” Mechanics Research Communications, vol. 37, no. 6, pp. 565–570, 2010.
View at: Publisher Site | Google Scholar
J. T. Oden and A. R. Eugene, Mechanics of Elastic Structures, Hemisphere Publishing Corporation, New York, NY, USA, 1981.
S. P. Timoshenko, “LXVI. On the correction for shear of the differential equation for transverse vibrations of prismatic bars,” The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, vol. 41, no. 245, pp. 744–746, 1921.
View at: Publisher Site | Google Scholar
S. P. Timoshenko, “On the transverse vibrations of bars of uniform cross-section,” The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, vol. 43, no. 253, pp. 125–131, 1922.
View at: Google Scholar
C. W. Lim, C. M. Wang, and S. Kitipornchai, “Timoshenko curved beam bending solutions in terms of Euler-Bernoulli solutions,” Archive of Applied Mechanics, vol. 67, no. 3, pp. 179–190, 1997.
View at: Publisher Site | Google Scholar
J. Lee, “In-plane free vibration analysis of curved timoshenko beams by the pseudospectral method,” KSME International Journal, vol. 17, no. 8, pp. 1156–1163, 2003.
View at: Publisher Site | Google Scholar
A. Lyckegaard and O. T. Thomsen, “Nonlinear analysis of a curved sandwich beam joined with a straight sandwich beam,” Composites Part B: Engineering, vol. 37, no. 2-3, pp. 101–107, 2005.
View at: Publisher Site | Google Scholar
S. P. Timoshenko and J. N. Goodier, Theory of Elasticity, McGraw-Hill, New York, NY, USA, 1970.
View at: MathSciNet
O. Kiliç and A. Aktaş, “Determination of stress functions of a curved beam subjected to an arbitrarily directed single force at the free end,” Mathematical and Computational Applications, vol. 7, no. 2, pp. 181–188, 2002.
View at: Publisher Site | Google Scholar
C. Bagci, “Exact elasticity solutions for stresses and deflections in curved beams and rings of exponential and t-sections,” Journal of Mechanical Design, vol. 115, no. 3, pp. 346–358, 1993.
View at: Publisher Site | Google Scholar
S. G. Lekhnitskii, Theory of Elasticity of An Anisotropic Body, Mir Publishers, Moscow, Russia, 1981.
View at: MathSciNet
J. Dryden, “Bending of inhomogeneous curved bars,” International Journal of Solids and Structures, vol. 44, no. 11-12, pp. 4158–4166, 2007.
View at: Publisher Site | Google Scholar
M. Wang and Y. Liu, “Elasticity solutions for orthotropic functionally graded curved beams,” European Journal of Mechanics - A/Solids, vol. 37, pp. 8–16, 2013.
View at: Publisher Site | Google Scholar | MathSciNet
A. Pydah and A. Sabale, “Static analysis of bi-directional functionally graded curved beams,” Composite Structures, vol. 160, pp. 867–876, 2017.
View at: Publisher Site | Google Scholar
Z. Shi, “Bending behavior of piezoelectric curved actuator,” Smart Materials and Structures, vol. 14, no. 4, pp. 835–842, 2005.
View at: Publisher Site | Google Scholar
Z. Shi and T. Zhang, “Bending analysis of a piezoelectric curved actuator with a generally graded property for the piezoelectric parameter,” Smart Materials and Structures, vol. 17, no. 4, Article ID 045018, 2008.
View at: Google Scholar
A. G. de Miguel, G. De Pietro, E. Carrera, G. Giunta, and A. Pagani, “Locking-free curved elements with refined kinematics for the analysis of composite structures,” Computer Methods Applied Mechanics and Engineering, vol. 337, pp. 481–500, 2018.
View at: Publisher Site | Google Scholar | MathSciNet
G. De Pietro, A. G. de Miguel, E. Carrera, G. Giunta, S. Belouettar, and A. Pagani, “Strong and weak form solutions of curved beams via Carrera’s unified formulation,” Mechanics of Advanced Materials and Structures, pp. 1–12, 2018.
View at: Publisher Site | Google Scholar
L. Y. Yu and L. L. Zhang, “Bending solution of functionally graded curved-beam,” Journal of Engineering Mechanics, vol. 31, no. 12, pp. 4–10, 2014.
View at: Google Scholar

Copyright

Copyright © 2019 Liangliang Zhang et al. 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.

PDF Download Citation

Download other formats

Order printed copies