Shock and Vibration

Shock and Vibration / 2016 / Article

Research Article | Open Access

Volume 2016 |Article ID 3978424 | https://doi.org/10.1155/2016/3978424

Runze Zhang, Yipeng Cao, Wenping Zhang, Hongbo Li, Xiangmei Li, "An Improved Fourier Series Method for the Free Vibration Analysis of the Three-Dimensional Coupled Beams", Shock and Vibration, vol. 2016, Article ID 3978424, 18 pages, 2016. https://doi.org/10.1155/2016/3978424

An Improved Fourier Series Method for the Free Vibration Analysis of the Three-Dimensional Coupled Beams

Academic Editor: Brij N. Agrawal
Received07 Jul 2016
Accepted17 Oct 2016
Published20 Dec 2016

Abstract

This paper presents a free vibration analysis of three-dimensional coupled beams with arbitrary coupling angle using an improved Fourier method. The displacement and rotation of the coupled beams are represented by the improved Fourier series which consisted of Fourier cosine series and closed-form auxiliary functions. The coupling and boundary conditions are accomplished by setting coupling and boundary springs and assigning corresponding stiffness values to the springs. Modal parameters are determined through the application of Rayleigh-Ritz procedure to the system energy formulation. The accuracy and convergence of the present method are demonstrated by finite element method (FEM) result. Investigation on vibration of the propulsion shafting structure shows the extensive applicability of present method. The studies on the vibration suppression devices are also reported.

1. Introduction

For the propeller running in the nonuniform wake fluid field, the longitudinal and lateral vibrations of the propulsion shafting will be stimulated. In the following, the shafting longitudinal and lateral vibration energy transfers from shafting to bearing pedestal, thereby to the structure which will cause structure underwater radiated noise. Therefore, it is necessary to find a method that can not only calculate the vibration characteristics of the propulsion shafting but also consider the shell structure in the following step.

There is a vast amount of information available on the longitudinal or lateral vibration of propulsion shafting. Zhang and Zhao [1, 2] obtained the longitudinal vibration modes of the propulsion shafting using one-dimension elastic wave approach based on its continuous model. The coupled longitudinal-transverse dynamics of the marine propulsion shafting under primary and internal responses was investigated by Zou et al. [3]. The coupled longitudinal-transverse dynamic model is established by using the extended Hamilton principle and discretized by Galerkin method. Huang et al. [4] researched the coupled torsional-longitudinal vibration of the propulsion shafting. The changes of displacement and rotation angle caused by coupled torsional-longitudinal vibration of the propulsion shafting are studied by the finite element method. Zhou et al. [5] applied the finite element method to calculate the ship propulsion shafting vibration responses under the propeller bearing forces. Li et al. [6] studied the dynamic and acoustic characteristics of the coupled propeller and shafting system by transfer matrix method. Kane and McGoldrick [7] studied the dynamic behavior of the propulsion shafting caused by longitudinal excitation. Pan et al. [8] used a series of mass spring damper system to simulate the propulsion shafting and performed experiment analysis, they simplified propeller as lumped mass, and the fluid load used the added mass to be simulated.

Among those papers, the shafting is assumed to be a uniform cross-section beam. A wide spectrum of techniques has been developed for the vibrations of beams. Some authors have investigated the dynamic analyses of multispan beams which used modern numerical methods such as the finite element method.

Bishop and Johnson [9] presented a modal receptance approach to the case single and coupled beam structure. Farag and Pan [10] applied the receptance approach technique to the dynamic response and power flower flow of the three-dimensional coupled beam structure. Lee [11] used the assumed mode method to research the transverse vibration of a Timoshenko beam which is acted on by accelerating mass. Wu and Hsu [12] studied free vibration analysis of simply supported beams carrying multiple point mass and spring-mass system. Li [1316] proposed an improved Fourier series method to establish the multispan beam model with the elastic support and calculated the natural characteristics and dynamic response under moving load. An analytical solution consisting of the expansion of the displacement function into a Fourier series, the calculation of the kinetic energy and the elastic strain energy, and then solving the Lagrange equation is presented by Heidebrecht in [17].

The two-dimensional coupled beams model is used in the shafting vibration calculation. However, the effect of the torsional force caused by the propeller can not be expressed well. Therefore, in order to study the characteristic of shafting vibration with the numerical method better, the three-dimensional coupled beams are selected as the research object in this paper. The displacements and rotation components of the coupled beams are expressed as an improved Fourier series in which the closed-form auxiliary functions introduced to remove any potential discontinuities and ensure the convergence of the solution. All the coefficients will be solved directly from the Rayleigh-Ritz procedure. The convergence and accuracy of the present method are verified by numerical examples.

2. Theoretical Formulations

2.1. Theoretical Model

The three-dimensional coupled beams model is presented in Figure 1, together with the global and local coordinates employed in this paper. This model indicates ample information including general boundary condition, arbitrary coupling angles, and elastically coupled condition. In order to simulate the general boundary conditions and the elastically coupled condition, artificial spring technique is adopted here. Specifically, six sets of stiffness-like springs are used at the end of the beam with rectangular coordinates, including three sets of linear springs, respectively, along the -, -, -directions, and three sets of rotational springs around the -, -, -directions. The general boundary conditions can be achieved by varying the value of springs’ stiffness. The springs are assigned very high stiffness values for a clamped edge or zero for a free boundary condition. Similarly, another six sets of springs are used to link the adjacent beams. For a rigid connection, the stiffness of all the springs can be assigned an infinite valve. Through assigning the springs at proper stiffness, the isolator can be readily achieved. The form of springs is shown in Table 1. The counterclockwise is predefined as the positive direction of the coupling angle which is variable in the range from − to . The angle () of each beam was obtained by the angle between the local and the global coordinate. The arbitrary coupling angle () can be achieved by varying the value of the angle () of each beam.


The boundary springs at the right ends of beam The coupling springs at the junction of beams and The boundary springs at the left ends of beam

-direction spring
-direction spring
-direction spring
Around -axis rotational spring
Around -axis rotational spring
Around -axis rotational spring

2.2. Energy Functional of Coupled Beams Structure

The total energy function of the coupled beams is taken as the sum of the energy contributions from the beam components, coupling constraints, and boundary constraints. Thus, the Lagrangian energy function of the coupled structure can be expressed as

is the coefficient, and if , ; otherwise .

Also the strain energy and the kinetic energy for the coupled beams can be written as

The potential energy stored in the boundary and coupling springs can be described aswhere

is the transformation matrix to transform to the global coordinate system, and its form is as follows:

2.3. Solution for the Coupled Beams Structure

The displacements of the beam are described by , , , and in the -, -, -directions and rotation angle around -direction. And they can be written in the form of the improved Fourier series as follows:where is the angular frequency, represents time, , , , , and are the Fourier expanded coefficient, and , , , and are the expanded coefficients of the auxiliary functions. They can be solved by the Rayleigh-Ritz procedure. are the auxiliary functions which can remove any potential discontinuities of the original displacements and effectively enhance the convergence of the result. In this paper, these auxiliary functions are specially selected as

Substituting (2)-(3) and (6) together with the admissible functions into (1) and then performing the Rayleigh-Ritz procedure, the governing eigenvalue equation can be obtained and given aswhere the matrix is the stiffness matrix of the system and its form is as follows:

in which represents the total stiffness matrix of the th beam. The form is as follows:

The matrix is the coupling stiffness matrix caused by the coupling spring between the th and th beams, and its form is as follows:

The detail form of the above equation can be found in Appendix A.

The system mass matrix can be written aswhere is the mass matrix of the beam and can be expressed as

The detail form of the above equation can be found in Appendix B.

The definition of the coefficient vector in (8) iswhere

All the natural frequencies and mode shapes of the coupled beams can be obtained by solving (8), which is a standard characteristic equation.

3. Numerical Example and Discussions

In this section, several examples are presented to verify the accuracy and reliability of the present method. Firstly, the convergence of the present solution is checked and excellent accuracy is validated through comparison of the results obtained from the FEM. Then coupled beams with arbitrary coupling angle are studied. Influence of coupling springs is discussed as well. The properties for these coupled beams structures are summarized in Table 2.


ParameterUnitBeam 1Beam 2

m0.240.8
m20.00090.0009
m4
Nm−2
Kgm−378007800

3.1. Convergence and Validation Study

Theoretically, there are infinite terms in the improved Fourier series solution. However, the series is numerically truncated and only finite terms are counted in actual calculations. The convergence of the present method will be checked. As an example, consider two beams connected together at 0° angle, as shown in Figure 2. Table 3 shows the first eight frequencies of the coupled beams with some classical boundary conditions, including C-C, F-F, and S-S, are derived by different truncation number . It is obvious that the natural frequencies converge quickly with the increasing truncated number. It can be also found that the results obtained by present method converge well at . Therefore, the truncated number will be uniformly selected as . And also, an excellent agreement of comparisons between the results and those from FEM shows the accuracy of present method. Some selected mode shapes for the C-C supported structure are given in Figure 3. It is obvious that the two sets of modes are essentially identical and the coupled beams are commonly coupled.


Boundary conditionsModeANSYSErrors/%

C-C184.5185.1585.4885.5085.6460.17
2235.15235.26235.61236.10237.190.46
3459.58461.71461.97462.17468.271.32
4773.66773.72773.98774.11781.460.95
51159.231160.051160.861160.881181.41.77
61641.331642.941643.211643.321674.11.87
72220.492221.012225.912226.402266.71.81
82917.792918.032923.762923.942967.71.50

S-S137.4037.5737.7037.7337.750.06
2150.26150.31150.58150.82151.560.49
3339.04339.12339.40339.53343.151.07
4603.43603.48603.54604.15615.401.86
5962.25962.49962.78962.88972.430.99
61384.011387.521387.631390.251419.702.12
71937.931938.411944.441945.211963.900.96
82569.582570.042583.402586.492613.101.03

F-F184.5485.3285.4685.4985.560.09
2234.93235.29235.64235.69236.790.47
3462.03462.03462.30462.36467.151.04
4763.66763.86764.72765.61779.021.75
51150.811151.171151.881152.571176.802.10
61633.811633.821634.621634.731666.101.92
72221.592221.782222.192222.432253.901.42
82924.932925.992926.582926.612948.100.73

3.2. Two Coupled Beams with Arbitrary Angle

In order to illustrate the application of the current method, the studies on free vibration of the coupled beams with a 45° and 90° coupling angle are carried out, as shown in Figure 4. The left end of beam 1 and the right end of beam 2 are clamped. Table 4 shows the comparison of nature frequencies obtained by present method and FEM. Some selected mode shapes are shown in Figures 5 and 6. It is seen that very good agreement of the result is obtained. The comparisons in Table 4 and Figures 5 and 6 indicated that the current analysis for coupled beams with various coupling angles is accurate and reliable.


AngleModel12345

= 90°Current91.35129.47192.02361.90396.75
ANSYS91.94129.91193.24365.96400.69
Error%0.640.340.641.120.99
Model678910
Current712.21738.021136.141161.201312.90
ANSYS727.85754.311158.001183.501343.50
Error%2.202.211.921.922.33

= 45°Current86.63128.99219.10362.03430.76
ANSYS86.15128.54217.32358.03425.28
Error%0.560.350.821.121.29
Model678910
Current718.88762.771106.601177.801191.00
ANSYS703.84755.181096.201154.531170.87
Error%2.141.000.952.021.72

3.3. Two Coupled Beams with Elastically Coupling Condition

In this section, the effects of the coupling spring on the mode shape of the coupled beams are investigated. Four types of coupling conditions as shown in Table 5 are considered in this part. The fourth mode shape of the coupled beams with various coupling conditions is present in Figure 7. Figure 7(a) shows that the mode shape is uncoupling and independent from each other when the coupling spring stiffness is . When the coupling spring stiffness changes from to , because of the elastic coupling, the displacement on the beams is no longer continuous across the junction. And there is a trend for the coupled beams from weak-coupling to strong-coupling. When the coupling spring stiffness is , the coupled beams realized the rigidity connection.


Four cases-direction spring-direction spring-direction springAround -axis rotational springAround -axis rotational springAround -axis rotational spring

Case 000000
Case
Case
Case

4. Application of Present Method to the Propulsion Shafting System

In this section, the propulsion shafting system will be taken as an example to illustrate the application of present method in practice. The shafting is assumed to be the beam structure with arbitrary boundary conditions. Thrust bearing is simulated as a beam which can be both transmit force and bending moment. The propeller is considered to be the lumped mass and the middle bearing and stern bearing are considered to be the linear springs. The propulsion shafting is simplified in Figure 8 and the properties for shafting are summarized in Table 6.


ParameterUnitBeam 1Beam 2

m0.42.68
m2
m4
N/m2
kg/m378007800

4.1. Vibration Characteristic of the Propeller Shafting System

For propulsion shafting with a propeller, a lumped mass is adopted, it needs to consider the influence of the lumped mass to the mass matrix of coupling system, and the energy of lumped mass is as follows:

Through (1) to (8), we can get the vibration characteristics of propulsion shafting under multisupport bearings with propeller.

In Table 7, the nature frequencies of the propulsion shafting are compared to those obtained by FEM. From Table 7, we can see that the present results agree well with the referential data. The deviation is small and does not exceed 2.5% for the worst case, which validates the accuracy and feasibility of present method for propulsion shafting system. Some selected mode shapes are shown in Figure 9.


Model12345

Current26.4227.935579.0282.41
ANSYS26.5327.8654.7478.6681.87
Error%0.420.260.470.450.65

Model678910

Current134.89146.75175.84181.77271.76
ANSYS133.4145.8174.02179.71266.28
Error%1.110.641.031.142.02

4.2. Vibration Isolator

Longitudinal vibration isolator can be achieved by assigning the -direction springs between shafting and thrust bearing at proper stiffness value. In this section, the damping effect of the isolator with different stiffness valve is studied. The stiffness values were separately taken as 1E9 N/m, 5E8 N/m, 1E8 N/m, and 5E7 N/m. From Table 8, we can see that, with the decreasing of stiffness valve, the two natural frequencies of longitudinal vibration are gradually reduced. The simulation of the isolator can be achieved.


Longitudinal vibration frequency

11376.921374.991361.181336.97
22381.142380.012368.702346.83

4.3. Dynamic Vibration Absorber

Dynamic vibration absorbers are commonly designed and tuned to suppress vibrations of one vibration mode of a vibrating structure. As shown in Figure 10, the spring-mass systems which are attached to propulsion shafting are simulated as the dynamic vibration absorbers.

The mass is and the stiffness coefficient of the spring between the mass and the beam is :where is the mass of the shafting and is the frequency of the vibration mode which is suppressed.

The displacement of the mass has a form of

Accordingly, the kinetic energy and the potential energy stored in the spring-mass system can be described as

Substituting (18)–(20) into the Lagrangian energy function, we can get the vibration characteristics of propulsion shafting with dynamic vibration absorbers.

In Table 9, we can see that there are two natural frequencies on both sides of the absorbency frequency after the dynamic vibration absorber is attached, and the other natural frequencies remain unchanged. Therefore, the frequency characteristic of the dynamic vibration absorbers is met. Accordingly the spring-mass system can achieve the simulation of the vibration absorber.


ModelOrigin modelAbsorber model

126.4226.53
227.9327.81
30.44
35555.16
479.0279.05
582.4182.58
6134.89135.22
7146.75146.89
8175.84176.1
9181.77181.9
10271.76271.8

5. Conclusions

An improved Fourier series method has been developed for the free vibration analysis of three-dimensional coupled beams with arbitrary coupling angle. The displacement and rotation of the structure are expressed as the Fourier cosine series and closed-form auxiliary functions. On the basis of energy functional of the structure elements, the solutions are obtained using Rayleigh-Ritz procedure. By using the present method, the three-dimensional coupled beams with arbitrary coupling angle and boundary conditions can be solved in a unified form. The accuracy and convergence of the present method are demonstrated by the FEM results. Based on this, the methodology is expressed to be applicable to vibration analysis for the propulsion shafting system and proved accurate for the FEM result as reference and obtained an efficient solution for the vibration suppression devices.

Appendix

A. Stiffness Matrices

Various stiffness matrices are represented as