Shock and Vibration
Volume 2018 (2018), Article ID 5015807, 12 pages
https://doi.org/10.1155/2018/5015807

Research Article

Differential Equations of Motion for Naturally Curved and Twisted Composite Space Beams

Ying Hao, Wei He, and Yanke Shi

School of Civil Engineering and Communication, North China University of Water Resources and Electric Power, Zhengzhou 450045, China

Correspondence should be addressed to Ying Hao; yinghao_ncwu@163.com

Received 24 August 2017; Revised 30 March 2018; Accepted 16 April 2018; Published 15 May 2018

Academic Editor: Toshiaki Natsuki

Copyright © 2018 Ying Hao 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.

Abstract.

The differential equations of motion for naturally curved and twisted elastic space beams made of anisotropic materials with noncircular cross sections, being a coupled system consisting of 14 second-order partial differential equations with variable coefficients, are derived theoretically. The warping deformation of beam’s cross section, as a new design factor, is incorporated into the differential equations in addition to the anisotropy of material, the curvatures of the rod axis, the initial twist of the cross section, the rotary inertia, and the shear and axial deformations. Numerical examples show that the effect of warping deformation on the natural frequencies of the beam is significant under certain geometric and boundary conditions. This study focuses on improving and consummating the traditional theories to build a general curve beam theory, thereby providing new scientific research reference and design principle for curve beam designers.

1. Introduction

Composite space beams with noncircular cross sections have extensive applications in civil, mechanical, aircraft, and aerospace constructions given their outstanding engineering properties, such as specific strength, specific stiffness, weight, and fatigue life over metallic counterparts. The static and dynamic characteristics of composite beams have attracted considerable attention in recent decades.

Elastic couplings occur among all the forms of deformation given the anisotropy of the material. Moreover, in-plane and out-of-plane displacements may be coupled. Therefore, the differential equations of motion, obtained by the 3D elastic theory under isothermal conditions, become more complex by allowing initial curvature and torsion. Euler-Bernoulli and Timoshenko beam theories are unable to accurately describe the behavior of composite beams.

Many researchers have conducted profound and careful research on nonclassical influences, such as transverse shear deformations, torsion-related warping, and several elastic couplings that can develop into composite beams. The natural frequencies and modes of composite beams under various boundary conditions are obtained by using the first-order shear deformation theory [1]. Based on the third-order shear theory, the forced vibration response of composite beams is investigated [2]. Li et al. [3] discussed the influence of different support conditions on the natural frequency of composite beams using many numerical examples. In this study, the effect of shear deformation is noted. Ramalingeswara and Ganesan [4] and Zapfe and Lesieutre [5] delivered the natural frequencies of laminated beams by using higher order shear theory and iterative method.

Most of abovementioned studies [1–5] focus on straight beams. The research on the dynamic behavior is far less in composite curved beams than in straight beams considering the curvature and torsion. Yildirim [6] studied the free vibration of symmetrically orthogonal laminated curved beams, in which rotary inertia and axial and shear deformation effects are considered. Tseng et al. [7] studied the free vibration of laminated curved beams with variable curvature based on Timoshenko theory. However, a transverse shear strain is constant along the thickness direction, and the shear shape factor was used in this research to adjust the lateral shear stiffness. The accuracy of the results significantly depends on a selected shear shape factor. Matsunaga [8, 9] used 1D high-order theory to solve the natural frequencies and buckling loads of composite curved beams subjected to axial concentrated force considering the shear deformation and moment of inertia. Baba and Thoppul [10] explored the flexural strength and vibration characteristics of sandwich laminated curved beams determined the effects of laminated debonding and curvature on their natural frequencies. Yildirim [11] derived the differential equations of motion of spatial curved bars of anisotropic materials, but the examples in this research remain straight beams.

The abovementioned studies [1–11] included one or several items of the anisotropy of the rod material, the curvature and torsion of the rod axis, the effects of rotary inertia, and the shear and axial deformations. The effect of the warping deformation of the cross section on natural frequencies is disregarded. In recent years, many researchers have attempted to establish generalized curved beam theory with warping effect. Kim et al. [12] developed a finite element-based beam analysis for anisotropic beams with arbitrary shaped cross sections through asymptotic expansion method. Displacements at each order are split into fundamental and warping solutions in their research. The influence of the warping effect is analyzed through examples, thereby demonstrating the importance of including the warping degrees of freedom in the dynamic analysis [13]. Yu et al. [14] formulated a generalized finite element-based cross-sectional analysis for inhomogeneous, initial curved and twisted, anisotropic beams, in which the transverse shear deformation is classified as part of the warping field. Ghafari and Rezaeepazhand [15] discussed the free vibration of rotating composite beams with an arbitrary cross section on the basis of dimensional reduction method. These authors decomposed a 3D elasticity problem of the composite beam into a 2D cross-sectional and a 1D beam analysis. Recently, Sapountzakis and Tsiptsis [16–19] used the Analog Equation Method and the Isogeometric Tool to solve the static and dynamic analysis of the curved beams with arbitrary cross section taking into account nonuniform warping and shear deformation effects. Nonuniform warping distributions are taken into account by employing four independent warping parameters multiplying a shear warping function in each direction and two torsional warping functions in these papers.

As an extension of [11, 20–22], the essential features and novel aspects of the present formulation compared with previous ones are summarized as follows.

(i)The improved beam theory is capable of spatial curved beams with arbitrary cross section, variable curvature, and composite material considering warping effect.(ii)The warping function of unidirectional composite beam with rectangular cross section can be obtained using the separation of variables method discussed in [23].(iii)Numerical examples are presented for two 3D curved beams with constant curvature and torsion. The results show that the effect of the warping deformation on natural frequencies of the beam is significant under certain geometric and boundary conditions.(iv)The influence of the warping deformation on the natural frequency is enhanced with the increase in the anisotropy of materials.

2. Geometrical Characteristics of a Space Anisotropic Beam

Considering a naturally curved and twisted anisotropic beam in space in Figure 1, indicate the tangential, normal, and binormal unit vectors of the beam axis, respectively. A second rectangular Cartesian frame is introduced to consider the initial twist of the cross section. The angle between the - and -axes is . For a smooth curve, the Frenet–Serret formulae are

where and are the arc coordinate, the curvature, and torsion of the curve, respectively. In order to take into account the initial twist of the cross section, a new Cartesian reference frame is introduced such that the -axis is in the direction of . If the unit vectors of and are represented by and , then

Figure 1: Geometry of a naturally curved and twisted anisotropic beam.

With simultaneous equations (1) and (2), we have

where , ,

3. Constitutive Equations

The following assumptions are used in the present study: the beam material is linearly elastic, homogeneous, and anisotropic. The pretwisting of the cross section is neglected . The generalized Hooke’s law for linearly elastic materials is

or

where , are the stiffness and flexibility matrices, correspondingly. Let

The following assumptions on stresses are valid in the classical rod theory.

can be represented as by solving (7) and (5) simultaneously. Then, the results are substituted to the residual expressions of (5), and the strain components are related to by

In (8),

So, is only related to but unrelated to when (8) is substituted to (4). Then, (4) is reduced to

The displacement components of an arbitrary point of the beam are given as [20–22, 24]

is not only a function of coordinates , but also a function of time under dynamic loads, in which

where , , are the three translations, , , are the three rotations, is a generalized warping coordinate which is used to describe the decay of deformation along beam length, and is the warping function of Saint-Venant’s torsion. Let , , be the axial and shear strain components, respectively,

in which , are shear correction factors, and . Based on the assumption that the initial curvature of the beam is gentle, then . Let

Let, for convenience, the resultant forces and moments on the cross section be expressed as , , , , , and shown as

where is the permutation tensor and ; ; .

The generalized warping moment is

The constitutive equations expressed by six displacement functions , and the generalized warping coordinate can be obtained by substituting (10), (13), and (14) into the above formulations (15)-(16).

4. Equations of Motion

The external force and moments per unit length of the beam axis are expressed as , . Also let , be the inertial force and the inertial moment per unit length of the rod axis. Owing to the displacements and the rotations that are assumed to be infinitesimal, the equations of motion can be written as [11, 24]

in which

The equation of motion in a scalar form can be obtained by substituting (12)–(14) into (17). The equations of , , , , , , can be obtained by solving the abovementioned equations simultaneously. The obtained results are substituted into the derivatives of the generalized warping moment to determine the final differential equations of motion of naturally curved and twisted elastic space beams of anisotropic materials.

5. Applications

5.1. Differential Equations of Motion for the Examples

A simple example is developed here for illustration. Two cylindrical helical springs with rectangular cross section can be regarded as 3D curved beams with constant curvature and torsion (Figure 2). The warping function of unidirectional composite beam with rectangular cross section (Figure 3) can be obtained using the separation of variables method [23]:

where , , , .

Figure 2: Cylindrical helical spring.

Figure 3: Unidirectional composite beam.

The differential equations of motion for cylindrical helical springs with rectangular cross section made of unidirectional composite materials are described as follows:

where

5.2. Effects of Material on the Natural Frequencies with Warping

The geometric parameters and finite element meshes of the two beams are exactly the same: the cylinder diameter  mm, the helix pitch angle , the number of active turns , and the rectangular cross section of  mm in the -direction and   in the -direction. The only difference between them is the material properties. Beam 1 is made of isotropic materials (the density is , the modulus of elasticity is , and Poisson’s ratio is ), while Beam 2 is made of unidirectional composite materials (the modulus of elasticity in the three main directions is , , the modulus of shearing of the three planes is , , the main Poisson’s ratio is , , and the density is ). We partitioned Beam 1 (Beam 2) into 720 Solid45 (Solid46) elements using ANSYS for analysis. The total number of nodal points is 1452. The first five natural frequencies of Beam 1 and Beam 2 are listed in Tables 1 and 2. The first five normal mode shapes of beams and the generalized warping coordinate are demonstrated in Figures 4–7, respectively.

Table 1: The first five frequencies of Beam 1 (made of isotropic material).

Mode number (Hz)
FEM
Present study (PS)

Warping ignored
Errors (%)
Warping included
Errors (%)

1
521.34
745.84
43.06
523.25
0.37

2
899.54
1009.30
12.20
897.71
0.20

3
915.97
1044.17
13.99
914.36
0.17

4
970.21
1122.22
15.67
973.88
0.38

5
1129.5
1338.65
18.52
1118.41
0.98

Table 2: The first five frequencies of Beam 2 (made of unidirectional composite material).

Mode number (Hz)
FEM
Present study (PS)

Warping ignored
Errors (%)
Warping included
Errors (%)

1
351.54
512.80
45.87
358.21
1.89

2
675.19
981.32
45.34
687.63
1.84

3
787.86
1117.87
41.89
805.74
2.27

4
880.95
1200.09
36.23
903.08
2.51

5
940.96
1357.22
44.24
958.56
1.87

Figure 4: First five mode shapes of Beam 1: (1) first mode (523.25 Hz); (2) second mode (897.71 Hz); (3) third mode (914.36 Hz); (4) fourth mode (973.88 Hz); (5) fifth mode (1118.41 Hz).

(a)

(b)

(c)

(d)

(e)

Figure 5: First five mode shapes of Beam 2 (made of unidirectional composite material): (a) first mode (358.21 Hz); (b) second mode (687.63 Hz); (c) third mode (805.74 Hz); (d) fourth mode (903.08 Hz); (e) fifth mode (958.56 Hz).

(a) First mode (523.25 Hz)

(b) Second mode (897.71 Hz)

(c) Third mode (914.36 Hz)

(d) Fourth mode (973.88 Hz)

(e) Fifth mode (1118.41 Hz)

Figure 6: First five normal mode shapes of the generalized warping coordinate of Beam 1.

(a) First mode (358.21 Hz)

(b) Second mode (687.63 Hz)

(c) Third mode (805.74 Hz)

(d) Fourth mode (903.08 Hz)

(e) Fifth mode (958.56 Hz)

Figure 7: First five normal mode shapes of the generalized warping coordinate of Beam 2.

6. Conclusion

The differential equations of motion for elastic space beams with noncircular cross sections and anisotropic materials have been derived systematically based on naturally curved and twisted beam theory. In the equations, the warping deformation of the cross section, the initial twist of the cross section, the shear and axial deformations, and the rotary inertia are studied theoretically.

The results show that the effect of warping deformations on the natural frequencies of the beam is very significant under certain geometric and boundary conditions. The errors can be excessively large when the warping effect is disregarded. The influence of warping on the natural frequency is enhanced with the increase in the anisotropy of materials.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work is partially supported by National Natural Science Foundation of China (11402090, 51779095, and 11502081).

References

K. Chandrashekhara, K. Krishnamurthy, and S. Roy, “Free vibration of composite beams including rotary inertia and shear deformation,” Composite Structures, vol. 14, no. 4, pp. 269–279, 1990.

M. P. Singh and A. S. Abdelnaser, “Random response of symmetric cross-ply composite beams with arbitrary boundary conditions,” AIAA Journal, vol. 30, no. 4, pp. 1081–1088, 1992.

J. Li, H. Hongxing, and S. Rongying, “Dynamic finite element method for generally laminated composite beams,” International Journal of Mechanical Sciences, vol. 50, no. 3, pp. 466–480, 2008.

R. S. Ramalingeswara and N. Ganesan, “Dynamic response of non-uniform composite beams,” Journal of Sound and Vibration, vol. 200, no. 5, pp. 563–577, 1997.

J. A. Zapfe and G. A. Lesieutre, “Vibration analysis of laminated beams using an iterative smeared laminate model,” Journal of Sound and Vibration, vol. 199, no. 2, pp. 275–284, 1997.

V. Yildirim, “Rotary inertia, axial and shear deformation effects on the in-plane natural frequencies of symmetric cross-ply laminated circular arches,” Journal of Sound and Vibration, vol. 224, no. 4, pp. 575–589, 1999.

Y. P. Tseng, C. S. Huang, and M. S. Kao, “In-plane vibration of laminated curved beams with variable curvature by dynamic stiffness analysis,” Composite Structures, vol. 50, no. 2, pp. 103–114, 2000.

H. Matsunaga, “In-plane vibration and stability of shallow circular arches subjected to axial forces,” International Journal of Solids and Structures, vol. 33, no. 4, pp. 469–482, 1996.

H. Matsunaga, “Effects of higher-order deformations on in-plane vibration and stability of thick circular rings,” Acta Mechanica, vol. 124, no. 1–4, pp. 47–61, 1997.

B. O. Baba and S. Thoppul, “Experimental evaluation of the vibration behavior of flat and curved sandwich composite beams with face/core debond,” Composite Structures, vol. 91, no. 1, pp. 110–119, 2009.

V. Yildirim, “Governing equations of initially twisted elastic space rods made of laminated composite materials,” International Journal of Engineering Science, vol. 37, no. 8, pp. 1007–1035, 1999.

J.-S. Kim, M. Cho, and E. C. Smith, “An asymptotic analysis of composite beams with kinematically corrected end effects,” International Journal of Solids and Structures, vol. 45, no. 7-8, pp. 1954–1977, 2008.

E. J. Sapountzakis and V. G. Mokos, “Dynamic analysis of 3-D beam elements including warping and shear deformation effects,” International Journal of Solids and Structures, vol. 43, no. 22-23, pp. 6707–6726, 2006.

W. Yu, D. H. Hodges, V. Volovoi, and C. E. S. Cesnik, “On Timoshenko-like modeling of initially curved and twisted composite beams,” International Journal of Solids and Structures, vol. 39, no. 19, pp. 5101–5121, 2002.

E. Ghafari and J. Rezaeepazhand, “Vibration analysis of rotating composite beams using polynomial based dimensional reduction method,” International Journal of Mechanical Sciences, vol. 115-116, pp. 93–104, 2016.

E. J. Sapountzakis and I. N. Tsiptsis, “B-splines in the Analog Equation Method for the generalized beam analysis including warping effects,” Computers & Structures, vol. 180, pp. 60–73, 2017.

I. N. Tsiptsis and E. J. Sapountzakis, “Generalized warping and distortional analysis of curved beams with isogeometric methods,” Computers & Structures, vol. 191, pp. 33–50, 2017.

E. J. Sapountzakis and I. N. Tsiptsis, “Generalized vibration analysis of beams including warping effects by isogeometric methods,” Journal of Vibration & Control, Article ID 107754631668567, 2017.

I. N. Tsiptsis and E. J. Sapountzakis, “Isogeometric analysis for the dynamic problem of curved structures including warping effects,” Mechanics Based Design of Structures and Machines, pp. 1–19, 2017.

A.-M. Yu and Y. Hao, “Warping effect in free vibration analysis of unidirectional composite non-cylindrical helical springs,” Meccanica, vol. 48, no. 10, pp. 2453–2465, 2013.

A. M. Yu and Y. Hao, “Effect of warping on natural frequencies of symmetrical cross-ply laminated composite non-cylindrical helical springs,” International Journal of Mechanical Sciences, vol. 74, pp. 65–72, 2013.

A. M. Yu and Y. Hao, “Free vibration analysis of cylindrical helical springs with noncircular cross-sections,” Journal of Sound and Vibration, vol. 330, no. 11, pp. 2628–2639, 2011.

R. Q. Xu, J. S. He, and W. Q. Chen, “Saint-Venant torsion of orthotropic bars with inhomogeneous rectangular cross section,” Composite Structures, vol. 92, no. 6, pp. 1449–1457, 2010.

Y. Hao and J. Guan, “Parametric study on the influence of warping deformation upon natural frequencies of die springs,” Shock and Vibration, vol. 2016, Article ID 6231676, 2016.