Abstract

The dynamic effects of cable attachment on a cantilever beam with tip mass are investigated by an improved Chebyshev spectral element method. The cabled beam is modeled as a double-beam system connected by springs at several discrete locations. By utilizing high order Chebyshev polynomials as basis functions and meshing the system at the locations of connections, precise numerical results of the natural frequencies and mode shapes can be obtained using only a few elements. The accuracy of this method is validated through comparing the results of finite element method and those of spectral element method in literature. The validated method is implemented to investigate the effects of parameters, including spring stiffness, number of connections, density, and Young’s modulus of cable. The results show that the mode shapes of the cabled beam system can be classified into two types: beam mode shapes and cable mode shapes, according to their main deformation. Their corresponding natural frequencies change in very different ways with the variation of system parameters. This work can be applied to optimize the dynamic characteristics of precise spacecraft structures with cable attachments.

1. Introduction

In spacecraft applications, modeling and evaluating dynamic behavior accurately during launch and in orbit are extremely important. Cables as power and signal transfers are used in almost all satellites. Previous tests have shown that cables have remarkable influences on precision space structures, not only changing the mass distribution but also adding damping and in general changing the entire dynamic response [1, 2]. The influence of cables on space structures as a structural element is a growing but relatively new concern, and as such, there are only a few published studies that investigate modeling space cables specifically [3].

Currently, the most commonly used method for modeling cabled structure is nonstructural mass method, including lumped mass method which adds lumped mass at attachment locations and distributed mass method which adds cable as a distributed mass by changing the density of the host structure. But both of these two methods have been proved insufficient to modeling precise cabled structures connected at discrete locations because they ignore the dynamic interaction between the cable and its host structure [4, 5]. In space mission, an imprecise dynamic model might degrade mission performance and destabilize structural control loops [6]. In order to improve modeling precision, the double-beam system was proposed for modeling cabled beams, in which cable is modeled as a beam to contain its bending stiffness and shear deformation, and springs are used to represent discrete connections [7]. Although the dynamic characteristics of a single beam can be obtained exactly by many ways, when a secondary beam is attached to the main beam by springs at several places, its solution becomes complicated. Previous studies mostly focus on beams with continuous elastica connection such as sandwich beam and composite beam [8, 9], or beams with string connected at their ends [10]. Only a few researchers have performed dynamic analysis for the double-beam system connected by springs at several discrete locations. Babuška et al. [11] established a finite element model for a cabled beam using software Nastran and got its frequency response functions. But the preprocess of finite element method (FEM) is usually tedious, especially when the number or location of connections is changed. Spak et al. developed a method to calculate cable properties that are suitable to use for modeling purpose from direct physical measurements of cables [3, 12], compared different damping models for space flight cables [13], and used the distributed transfer function method to model cables and a simple cabled structure [5]. Choi [14, 15] studied the natural frequencies and damping of a cabled-harnessed structure using spectral element method (SEM). The effects of the number of connections, the spring stiffness of interconnections, and the mass portion of an attached cable had also been analyzed. However, this method is actually a semianalytical and seminumerical method, needing to resolve analytical solution of dispersion relations for every type of model. Besides, because the SEM is a time-domain method, it is inconvenient to obtain the mode shapes of cabled structures. Remedia et al. [16] investigated the influence of harness on spacecraft honeycomb panel through an experimental campaign. Zhang et al. [17] presented a general analytical approach for studying transverse vibrations of double-beam systems, made of two Timoshenko beams connected by discrete springs and coupled with multiple discontinuities.

However, mostly previous researches focus on beams with free-free boundary, and the changes of mode shapes due to attached cables are also rarely discussed. Cantilever beam with a tip mass as shown in Figure 1 is a typical model for many beam-like space structures equipped with instrument at its tip. Cables are usually essential for signal and power transmission between the hub and tip equipment. If the tip equipment is a very sophisticated machine and has a high request to the dynamics characteristics, the effects of the attached cable must be taken into account. This is just the goal of this research. To that end, the effects of cables property parameters, such as density and Young’s modulus, need to be investigated and the influences of connections’ number, stiffness, and locations also should be studied. In order to get the precise numerical solutions of both natural frequencies and mode shapes, our approach used here is a Chebyshev spectral element method (CSEM) improved from spectral-Tchebychev technique which was developed by Yagci for solving linear and nonlinear beam equations [18, 19]. This technique uses orthogonal Chebyshev polynomials as basis functions and applies Galerkin’s method to obtain the solution. It has high accuracy and exponential convergence, and only a small number of polynomials are sufficient to obtain the machine-precision accuracy. It also requires fewer elements than the FEM approach, not needing remeshing after the connections’ locations changed. Using projection matrices, any linear boundary conditions can be incorporated. Different from the meaning of SEM in literature [14] where Fourier spectral decomposition is used for time-domain analysis, “spectral” in CSEM means to perform spatial discretization in elements using Chebyshev spectral method.

2. Problem Description

Consider a cantilever beam with tip mass and cable attachment in Figure 1 which can be modeled as a double-beam system as shown in Figure 2. Both the host beam and the cable attachment are considered as Timoshenko beams to include the effects of shear deformation and rotatory inertia. They are connected by translational and torsional springs at several discrete points; and are the translational and torsional stiffness, respectively. The rotational inertia of tip mass is ignored. Because the two ends of the cable are usually slack, free-free boundary conditions are adopted. The excitation levels are thought to be in the linear region. Nonlinear behavior, such as slack cables slapping against the host beam, is not considered here. Based on this model, the Chebyshev spectral element model will be established for this cabled beam system, and numerical results of natural frequencies and mode shapes will be compared with the results of FEM and SEM. Then this validated model will be used to investigate the changes of dynamic characteristics of a cabled beam with variation of cable parameters, such as density and Young’s modulus. The influences of the connection number and springs’ stiffness will also be discussed.

3. Chebyshev Spectral Element Method

The discrete connections divide the beam and the cable into several parts. The CSEM used herein meshes the beam and the cable in the places of their connections and then assemblies them like FEM. It is assumed that there are connections, cable elements, and beam elements. Adding the tip mass element, there are totally elements. But unlike FEM which use low order polynomials in “small” elements, this method uses high order Chebyshev polynomials to approximate the deformation in “big” elements. The main advantage of the CSEM is that no remeshing is needed after the locations or number of the connections changed. Besides, only needing to change the number of polynomials, a good balance between accuracy and efficiency can be obtained. It can be considered as a combination of the Chebyshev spectral method (CSM) and the finite element method, combining the advantages of both methods. Furthermore, in order to deal with different boundary conditions in the same way, a projection matrix method will be used. In this section, we firstly review the properties of the Chebyshev polynomials and then give the Chebyshev spectral element model of this double-beam system.

3.1. Chebyshev Polynomials

The Chebyshev polynomials are a group of recursive orthogonal polynomials which can be given by [18]where and . For an arbitrary interval , a mappingcan be defined, and then scaled Chebyshev polynomials will be used. The Chebyshev polynomials form a complete set, which means any square-integrable function which is infinitely differentiable in the interval can be approximated by a truncated Chebyshev series expansion as where is the number of polynomials used for the truncated expansion.

Adopting Gauss-Lobatto sampling method,The values of function at spatial sampling points form a vector . There exists a one-to-one mapping between the expansion coefficients vector and sampling value vector : where is an forward transformation matrix and is the inverse matrix of .

The derivative of each Chebyshev polynomial can be recursively expressed in terms of lower order Chebyshev polynomials. Therefore, the approximation of the spatial derivative can be written as where coefficient vector can be got from through a different matrix which just relates to the number of Chebyshev polynomials and the definition domain of .

Similarly, the expansion coefficient vector of the th spatial derivative of can be given asand then the th spatial derivative of can be obtained aswhere is called the th derivative matrix.

The definite integral of can be described aswhere is the definite integral vector.

The inner product of functions and can also be expressed by Chebyshev polynomials as follows:where is the inner-product matrix.

3.2. Timoshenko Beam Element

The governing equation of a uniform Timoshenko beam is represented by [13]where and are the transverse displacement and the slope due to bending, respectively. is the material density, is the cross-sectional area, is the shear modulus, is Young’s modulus, and is the second moment of inertia. Damping for the displacement and rotation is included through coefficients and . and are the applied force and moment, respectively.

Assume that the approximate solution of (12) in the spectral form isHere, Chebyshev polynomials are chosen as trial functions. Due to the fact that (13) is an approximation, submitting (13) into (12) will produce an error or residual. According to Galerkin’s method, the weighted residuals must be minimized; namely,where and are the weighting functions coinciding with the trial functions. Using the inner-product expression from (11), the residuals can be rewritten as Since this equation should be satisfied for arbitrary and , there must beand this equation can be expressed in matrix form aswhereHere, is the element displacement vector, is a identity matrix, and is the th Chebyshev derivative matrix defined in (9). This is the discrete equation of motion of a Timoshenko viscous beam without imposing boundary conditions.

3.3. Beam-Cable Coupled System

Now, consider the double-beam system in Figure 2 which has elements. For the tip mass, at there areDefine where subscripts and represent beam and cable, respectively; ; is the global displacement vector, whose components and are elements’ displacement vectors defined in (19); , , and are matrices; , , , , , and are matrices defined in (20) and (21), and . is a column vector whose jth element is 1 and all other elements are zeros. are matrix block with all 0 elements. For clarity, in (25), the number of rows and columns is indicated by subscript.

Considering the relation between spring forces and element displacements at in Figure 2, the forces and moments acting on elements areNotice that all these forces are functions of displacements at connection locations. So, the global force vector due to spring connections can be written in terms of global displacement vector aswhere is a diagonal matrix whose th and th entries are −1 and and th entries are 1, . According to (17), we can assembly the global equation and getThis is the global discrete equation of the beam-cable coupled system without considering element compatibility and boundary conditions.

3.4. Compatibility and Boundary Conditions

Adjacent elements and in Figure 2 should meet displacement compatibility. That is, at the common points , the two sides are required to have equal displacements, slopes, forces, and moments. Mathematically, these compatibility conditions can be expressed asConsidering (9), their discrete form can be written asBesides, when spring stiffness is big enough, these connections in Figure 2 can be thought as rigid constraints. Instead of using global force vector in (28), we can handle them like compatibility conditions in (30). For instance, when linear spring stiffness is big enough, we can getSimilarly, generic linear boundary conditions can be written aswhere is the order of derivative, indicates th element, and denotes the particular boundary (right or left boundary of the element). Its discrete isFor example, for free-free boundary condition, there areTheir discrete form can be written as(31), (32), and (36) can be written in terms of global displacement vector aswhere is the compatibility and boundary condition matrix whose rows are determined by (31), (32), and (36) and is a column vector whose element is the temporal part of the boundary conditions. Here, is a zero vector.

An effective way of imposing the boundary is using projection matrices [19]; in this approach can be expressed as where and are and projection matrices, respectively, which can be determined using singular value of . And is a vector which satisfies homogeneous boundary conditions, while ensuring satisfies all the boundary conditions. And is the number of compatibility and boundary conditions.

Substituting (38) into (29) and left-multiplying both sides with result inwhere(39) is the discrete governing equation of the cabled beam system. Using state-space method, the numerical results of the natural frequencies and modal shapes can be obtained.

4. Method Evaluation

To evaluate the accuracy and convergence of CSEM, two cases were studied, including a double-beam system with free-free boundary and a cantilever cabled beam with tip mass. Their results were compared with the results of SEM in literature [14] and the results of FEM based on Nastran, respectively.

4.1. Free-Free Double-Beams System

Consider the free-free double-beam system in literature [14], whose material properties are shown in Table 1. Using CSEM, the calculated first 6 natural frequencies are summarized in Table 2, where 8 polynomials were used in every element. The results in literature [14] where SEM was used are also included in Table 2. It can be seen that the differences between CSEM and SEM are all less than 0.100%.

Furthermore, the mode shapes can also be obtained using CSEM, as shown in Figures 3 and 4. The blue lines represent the beam’s mode shapes, and the orange lines indicate the cable’s mode shapes. Comparing with the mode shapes of a bare beam, it is obvious that the mode shapes can be classified into two categories based on their main deformation component: mode shapes dominated by beam deformation (beam mode) and mode shapes dominated by cable deformation (cable mode). The corresponding natural frequencies of beam modes are underlined in Table 2. It should be noticed that the first four natural frequencies of the bare beam are 0 Hz, 11.1792 Hz, 30.8144 Hz, and 60.4046 Hz, respectively. They decrease slightly after being attached to a cable. It is also obvious that when the cable and beam are connected by 3 springs, the first flexible mode is no longer a beam mode. This indicates that the number of connections is a key factor to determine the first flexible mode shape of the cabled beam system.

4.2. Cantilever Beam with Tip Mass and Cable Attachment

Recall the cantilever beam with tip mass and cable attachment in Figure 1. The following parameters were used for the beam: length  m, radius of cross section  m, Young’s modulus  GPa, material density kg/m3, shear factor , Poisson ratio , and the tip mass  kg.

First, taking no account of the cable, we can get the natural frequencies and mode shapes of the cantilever beam with tip mass as shown in Figure 5 utilizing CSEM where 10 Chebyshev polynomials were used totally. Meanwhile, we use the finite element software MSC/NASTRAN to resolve this problem. A lumped point mass element (CONM2) and 50 beam elements (CBAR) were used to represent the tip mass and cantilever beam, respectively. The results we got are 11.626 Hz, 76.450 Hz, and 219.85 Hz. Both results agree well, but obviously CSEM is more efficient because it can get a high precise result only using a few polynomials.

And then, we consider the cable’s effects using CSEM. The parameters of the cable used here include length  m, radius  m, density  kg/m3, Young’s modulus of  GPa, shear factor , and Poisson ratio . The connection parameters include linear spring stiffness  N/m and torsional spring stiffness  Nm/rad, and number of connections is 5; that is, . The convergence of the first five natural frequencies is shown in Figure 6. The results show that, only using 6 polynomials, the relative error converges to . Table 3 lists both the results of FEM using different number of elements and the results of CSEM using different number of polynomials in each element. Figure 7 shows the first five mode shapes. We can see that the results from CSEM are fairly close to those of FEM.

5. Cable Effects Analysis

There are many factors that might have influences on the dynamic characteristics of a cabled beam with tip mass. In this section we focus on the effects of connection parameters and cable parameters; both of them might change in a big range.

5.1. Effects of Connection Parameters

Using the same parameters in Section 4.2, the change of natural frequencies with the variation of connection spring stiffness is shown in Figure 8. The first natural frequency is nearly unchanged, and the other ones increase and finally converge to stable values with the increasing of spring stiffness . When  N/m, all of the first five natural frequencies stabilize. Observe the mode shapes and natural frequencies shown in Figure 7 in which  N/m; it is obvious that the 1st, 4th, and 5th modes are beam modes and the 2nd and 3rd modes are cable modes. So it is concluded that the frequencies of beam mode increase and converge faster than those of beam mode, and the lower frequencies converge faster than higher ones. It also can be seen that the changes of high order frequencies are much bigger than those of low order frequencies.

The change of natural frequencies with the variation of the number of connections is illustrated in Figure 9; here  N/m. As the number of connections increases, the natural frequencies of cable modes increase rapidly and those of beam modes, as shown in purple circle, keep nearly unchanged. This phenomenon results in mode shapes exchanges between adjacent modes. Remember that when the number of connections is 5, the 1st, 4th, and 5th modes are beam modes. Observe the mode shapes in Figure 10 where the number of connections changed to 6, the 1st, 2nd, and 5th modes become beam modes, and 3rd and 4th modes become cable modes. That is because the natural frequencies of cable modes changed from 60.360 Hz and 60.705 Hz to 75.821 Hz and 82.233 Hz, as shown in figure points and .

5.2. Effects of Cable Parameters

Modeling a spaceflight cable as a homogenous beam requires inputting various equivalent parameters, such as density and Young’s modulus. In engineering practice, there are many different kinds of cables with a wide range of material properties. The effects due to the changing of cable parameters, including density and Young’s modulus, are discussed in this section.

Figure 11 illustrates the variation of the natural frequencies with the cable’s density. It can be observed that the natural frequencies of cable modes decline gradually with the increasing of the cable’s density. The low order natural frequencies of beam modes keep nearly unchanged, and high order ones have a slight decline. As a result, adjacent frequency curves will cross at a series of special points, such as points and in Figure 11. Compare the mode shapes at points and , where the density, respectively, is 1500 kg/m3 and 3500 kg/m3, as shown in Figure 12. The natural frequency decreases to 213.12 Hz, and the deformation of the cable increases obviously, due to the increasing of density.

Figure 13 shows the change of natural frequencies with the increasing of Young’s modulus of cable material. It can be seen that all the natural frequencies of beam modes are nearly constant. But those of cable modes increase rapidly. Comparing the mode shapes at points and as shown in Figure 14, where the values of Young’s modulus are 10 GPa and 50 GPa, it can be seen that their natural frequencies almost have the same value. With the increase of Young’s modulus, the deformation curve of the cable mode becomes more close to that of the beam mode.

6. Conclusion

In this paper, a Chebyshev spectral element method (CSEM) is presented for exploring the dynamic effects of cable attachment on a cantilever beam with tip mass. The cabled beam is modeled as a double-beam system connected by several springs at discrete locations. Using high order Chebyshev polynomials as basis functions, precise numerical results of system’s natural frequencies and mode shapes can be obtained using fewer elements. This method divides the system according to the location of connections, not needing remeshing after the number of connections changed.

In addition, projection matrices are used to deal with differential boundary conditions. The results were compared with those of spectral element method in literature and those of finite element method based on software Nastran. It can be seen that the CSEM has a good convergence and accuracy.

Based on the validated CSEM, the effect of connection parameters, including their stiffness and number, was studied. The results show that the mode shapes can be classified into two types: beam mode shapes and cable mode shapes, according to their main deformation. With the increase of connection stiffness, the natural frequencies increase and finally converge to stable values. The low order frequencies increase and converge faster than the high order ones. The frequencies of beam modes increase and converge faster than those of cable modes. Changes of high order frequencies are much bigger than those of low order frequencies. As the number of connections increases, the natural frequencies of cable modes increase rapidly and the natural frequencies of beam modes keep nearly unchanged. The effects of cable’s parameters, such as density and Young’s modulus, are also investigated. It can be found that, with the increase of cable’s density, the natural frequencies of cable modes decrease far more rapidly than those of beam modes. With the increase of Young’s modulus of cable material, the natural frequencies of the beam modes keep nearly unchanged, and those of cable modes increase gradually. Meanwhile, the deformation curve of cable modes becomes more close to those of beam modes. These results can be applied to optimize the dynamic characteristics of precise cable-harnessed spacecraft structures.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work is supported by the National Basic Research Program of China (Grant no. 2013CB733004).