Next Article in Journal
Dynamic Behavior of Track Bridge on Serviced Long-Span Road/Railway Bridge
Previous Article in Journal
Substitution Effects on the Optoelectronic Properties of Coumarin Derivatives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Partial Frequency Assignment for Torsional Vibration Control of Complex Marine Propulsion Shafting Systems

1
College of Power and Energy Engineering, Harbin Engineering University, Harbin 150001, China
2
School of Engineering, University of Liverpool, Liverpool L69 3GH, UK
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(1), 147; https://doi.org/10.3390/app10010147
Submission received: 22 November 2019 / Revised: 15 December 2019 / Accepted: 16 December 2019 / Published: 23 December 2019
(This article belongs to the Section Acoustics and Vibrations)

Abstract

:
With the large-scale and complexity of ship propulsion shafting, it is more difficult to analyze and control the torsional vibration of shafting. Therefore, an effective control method for the torsional vibration of shafting is of great significance in the field of ship engineering. The main strategy of torsional vibration control adopted in this paper is to keep the natural frequency of a shaft system away from the excitation frequency through structural modifications. In addition, because the basic parameters of much of the equipment in engineering applications cannot be changed, this restriction cannot be ignored when seeking solutions related to structural modifications. This paper studies the partial eigenvalue assignment for the torsional vibration control of complex ship propulsion shafting using the gradient flow method, which can shift a “dangerous” natural frequency to a safe value, while satisfying complex physical constraints. The models of a ship propulsion system and a diesel generator set are established to demonstrate several different desired modification schemes and constraint conditions in practice. In particular, close frequencies are shifted. The numerical simulation results demonstrate that it is effective and feasible to make a partial frequency assignment of torsional vibration, which provides a reliable approach for the control of torsional vibration for complex shaft systems in practical engineering.

1. Introduction

In recent years, with the improvement of ship power and the widespread use of large-scale and complex shafting systems, the vibration problem of shafting systems becomes more and more serious. The torsional vibration of a shafting system is one of the main causes of shafting fracture and other faults, which has attracted wide attention. The dynamic performance of torsional vibration is one of the main aspects to guarantee the productivity and safety of machines and structures, and thus has great economic and social significance. However, there are many situations in which the structural design cannot meet the desired structural dynamic characteristics so that some accidents happen in serious cases, especially shaft breakage. Structural modifications (SMs) are known to be a very important technique to improve the dynamic characteristics of a structure in the form of fundamental structural parameters (such as mass, stiffness, and damping) to avoid resonance or the creation of a node on the system at a certain frequency. Additionally, SMs are also a procedure aimed at identifying the changes needed and determining the values of the physical parameters of a structure to achieve desirable dynamic characteristics (usually modal properties such as natural frequencies, mode shapes, and frequency response). Several theories and many algorithms have been established to solve SMs as a direct problem and as an inverse problem. Early studies about direct structural modification are mainly based on small modifications, localized modifications, and modal approximation [1]. The first useful method for structural modification was formulated and broaden by many researchers [2], which was based on a matrix perturbation approach to get an approximate solution for the modal properties of the modified structure. However, the limitation of the perturbation approach is that the modification should be small in order to get an accurate solution. Other methods such as the analytical approach of structural modification by adding different mass and stiffness attachments [3,4,5], the Rayleigh quotient method [6], and the sensitivity analysis method [7,8] were applied as useful tools for structural modification.
The inverse problem, being more complex, tries to identify, in the frame of a given set of possible locations, the most appropriate modifications required to obtain the desired dynamic behavior [9]. One major objective of doing that is to keep the frequency of excitation away from a natural frequency so that a large structural response can be prevented. The methods for inverse structural modifications are based on physical model parameters [10,11,12] (such as the mass matrix, stiffness matrix, and damping matrix) and experimental modal analysis [13,14,15,16] (e.g., natural frequencies and modal shapes derived from modal testing or directly by the frequency response function) by either passive procedues, active procedures, or combining the above two procedures in a hybrid approach. There are large numbers of published papers in the literature describing the study of inverse problems, and therefore only a brief introduction to some major ones is given as follows.
Hu et al. [17] studied the problem of natural frequency assignment for mass-chain systems with inertias to determine whether an arbitrary set of positive numbers can be assigned as the natural frequencies of a chain of n masses. The solution and solvability condition for the inverse problem of linear undamped gyroscopic systems were presented by Yuan and Dai [18], and the results were applied to model updating. A geared rotor-bearing system was designed and tested by Tsai et al. [19] to validate the method for solving the inverse structural modification using the measured receptances. Liu et al. presented the strategy for the torsional vibration control of shaft systems based on the measured torsional receptances of the system [20].
All of these above approaches can effectively change some natural frequencies of a structure to targeted values. However, a problem has received less attention: the undesired changes of the other natural frequencies of the original structure as a result of the modifications, which is known as the frequency “spill-over” phenomenon. In order to overcome this problem, many scholars have studied this challenging research topic using the partial eigenvalue assignment problem, which aims to assign some natural frequencies to desired locations and keep all the other frequencies unchanged.
Partial frequency assignment can be achieved fairly easily by means of active control. Many researchers [21,22,23,24,25] studied this problem based on the stiffness matrix perturbation, the receptance measurements, or a combination of these two. Bai et al. [26,27,28] presented hybrid optimization methods for the robust and minimum norm partial quadratic eigenvalue assignment problems, which are suitable for systems with and without time delays. Based on [28], Liu et al. [29] proposed a modified conjugate gradient method to solve the problem only using receptance. Compared with the hybrid method given by Ram et al. [22] for the single-input systems, Bai [30] solved the multi-input partial quadratic eigenvalue assignment problem using a multi-step hybrid method combining the system matrices and the receptance measurements. Ram and Mottershead [31] described a new theory that was applicable to both the single-input-multiple-output and multiple-input-multiple-output vibration control of practical engineering structures using the receptance method. Zhang et al. [32,33,34] presented some algorithms to solve the partial eigenvalue assignment for high-order systems by single-input control or multi-input control. Many researchers [35,36,37] also studied this problem based on various feedback signals, such as the sound pressure signal, displacement feedback signal, and acceleration feedback signal.
In fact, considering the low cost and maintenance required for system stability, it is expected to achieve a partial eigenvalues (or natural frequencies) assignment, depending on passive control or passive structural modification. However, this is a far more difficult task and challenge; as a result, few research papers have been introduced to solve this problem.
In recent years, some researchers have tried to study partial frequency assignments via passive control, which is more challenging but more useful than active control. Passive structural modifications of mass-spring systems for the partial assignment of natural frequencies were addressed by Ouyang and Zhang [38], and two different methods were presented. Compared with the work in [38], Belotti et al. [39] proposed a new method to assign a subset of natural frequencies of general structures with low spill-over, and optimal structural modifications were determined through a three-step procedure considering both the prescribed eigenvalues and the feasibility constraints.
The method mentioned above in [38] was shown to work for only simple systems and small modifications. In this paper, the main objective is a further theoretical development of the passive partial frequency assignment methods proposed by the research group of one of the authors (Ouyang and Zhang [38]). Its practical application to realistic large systems, which involves more stringent physical constrains, and its feasibility to bigger frequency modifications are investigated. In addition, it is for the first time demonstrated that the assignment of close eigenvalues can be achieved through the theory of structural modification.
In this paper, a simplified model of a “real” marine diesel engine propulsion system and a complex branched shaft system of a diesel generator set for torsional vibration control are investigated to achieve the partial assignment of natural frequencies. Various operation conditions and practical restrictions are considered. Additionally, modifications for assigning close eigenvalues and bigger frequencies are designed to show the strong applicability of the method. The outline of this paper is as follows. In Section 2, the passive control method is presented. The desired modifications and the achieved results of the simplified “real” propulsion system with different physical constraints at various operating conditions are analyzed in Section 3. Section 4 deals with the assignment of close eigenvalues of a complex branched shaft system. Finally, some conclusions are drawn in Section 5.

2. Theoretical Development

2.1. The Problem Description

Suppose that a general n-degree-of-freedom vibration system consists of mass matrix M 0 and stiffness matrix K 0 . The damping matrix is neglected, as the damping of many practical systems is small, such as a shafting that undergoes torsional vibration concerned in this paper.
The undamped vibration system is described by a second-order ordinary differential equation as follows:
M 0 q ¨ + K 0 q = 0 .
In general, M 0 is a symmetric and positive definite matrix, and K 0 is a symmetric and positive-semidefinite matrix. q G n is the displacement vector, and q ¨ G n is the acceleration vector.
Assuming the displacement vector q t = x e i ω t and λ = ω 2 , and substituting them in Equation (1) yields
K 0 λ M 0 x = 0 .
According to the Cholesky decomposition, mass matrix can be defined as M 0 = D 0 2 . And then supposing u = D 0 x , Equation (2) can be turned into:
J 0 λ I u = 0
where
J 0 = D 0 1 K 0 D 0 1
so that J 0 has the same eigenvalues as the original system M 0 , K 0 .

2.2. Construction of a Real Smmetric Matrix J s with All Expected Eigenvalues

In order to achieve the goal of partial assignment of the natural frequencies of the system, the n system eigenvalues are denoted by λ 1 , λ 2 , , λ n . It is expected that a set of p < n original system eigenvalues Λ 1 = diag λ 1 , λ 2 , , λ p are replaced by the desired eigenvalues μ 1 = diag μ 1 , μ 2 , , μ p after accomplishing some related system mass and stiffness modifications. Furthermore, Λ 2 = diag λ p + 1 , λ p + 2 , , λ n contains all the eigenvalues that remain unchanged. X 1 = x 1 , x 2 , , x p is the mass-normalized eigenvector matrix of Equation (2) corresponding to Λ 1 .
From the spectral decomposition theorem of symmetric matrices, a new matrix J s can be expressed as follows:
J s = U 1 μ 1 U 1 T + U 2 Λ 2 U 2 T
which contains the desired eigenvalues of the modified system and has the same eigenvectors as J 0 . U 1 and U 2 are the normalized eigenvector matrices of Equation (3), corresponding to Λ 1 and Λ 2 , respectively.
Based on the relationships of U 1 = D 0 X 1 and U 2 Λ 2 U 2 T = J 0 U 1 Λ 1 U 1 T , Equation (5) can be changed to:
J s = J 0 + D 0 X 1 μ 1 Λ 1 X 1 T D 0 .
However, it should be noted that J s is not easy to determine in practice, although it contains all the desired eigenvalues; that is, it is difficult to find the relationship between J s and the mass and stiffness matrices. Then, it is useful to reconstruct a matrix J , which is in the same matrix form as J 0 so that it can be easily used to obtain the system matrices of new mass-spring systems. Besides, many machines have fixed or complex structural configurations, which are not easy to change or even not allowed to be changed. So, there are many special physical limitations to structural modifications in engineering applications. In order to solve this problem, it is critical to find a method that can not only reconstruct a matrix J , but also satisfy special constrains. In this paper, the method introduced in Ref. [40] will be applied to solve this problem, and then a partial frequency assignment for the torsional vibration control of shafting systems will be achieved.

2.3. The Gradient Flow Method

The gradient flow method described in [40] is briefly presented first and then applied to complex and realistic cases.
Let ξ G n × n denote a real symmetric matrix with distinct eigenvalues γ 1 , γ 2 , , γ n . In this paper, ξ = J s . According to the orthogonal similarity transforms, all matrices that have the same eigenvalues as ξ can be defined by the following matrix set
Q ξ = V ξ V T
where V is an n × n orthogonal matrix.
Some matrix entry locations are defined as L = i τ , j τ τ = 1 , and the corresponding desired values are set to δ = δ 1 ,   δ 2 ,   ,   δ . Then, all the matrices with the prescribed entries at the desired locations can be expressed
N L , δ = A R n × n A i τ j τ = δ τ , τ = 1 , , .
Any given matrix Z in Q ξ can be expressed as the sum
Z = Z L + Z L s
where the entries in Z L are the same as those in Z , except those at the matrix entry locations that do not belong to L are set identically at zero; and L s represents all locations that do not belong to L . Then, the projection P Z of any matrix Z onto the affine subspace N L , δ is given by
P Z = C L + Z L s
where C L is a constant matrix in N L , δ with zero entries at all locations corresponding to L s and known entries at all other locations.
A function can be defined to minimize the distance between Z and N L , δ by
f Z = 1 2 Z P Z , Z P Z
where Z P Z , Z P Z denotes the Frobenius inner product and Z P Z = Z L C L .
Since Z belongs to Q ξ , Z can be expressed as Z = V ξ V T . Then, this minimization with the objective function f Z can be rewritten in terms of V as follows:
h V = 1 2 V ξ V T P V ξ V T , V ξ V T P V ξ V T .
The gradient h of objective function h is given as:
h = V ξ V T P V ξ V T V ξ V ξ V T T V ξ V T P V ξ V T V .
Then, a gradient flow h V can be defined by
d V d t = Z T , Z P Z V
where Z P Z , Z T denotes the Lie bracket commutator, i.e., B , G = B G G B .
The system of ordinary differential Equation (14) can be readily integrated from a starting point, say, V 0 = I (the identity matrix). The integration termination criterion of Equation (14) can be chosen as follows:
min Z t k P Z t k F , Z t k T , Z t k P Z t k F 10 8
where     .   F denotes the Frobenius norm of a matrix. It should be noted that if a solution does not exist, a least-squares solution can be obtained from the formulation.
Combining Equations (14) and (15), orthogonal matrix V and matrix Z can be obtained. Then, the matrix can be obtained by Z = J according to the definition of Equation (7).

2.4. The Mass and Stiffness Matrix after Modification

In Section 2.3, the method is described to reconstruct a matrix J aiming at determining the mass and stiffness matrices after modification. Then, in this section, the process of obtaining the modified mass and stiffness matrices are shown below.
Since the mass matrix for the torsional vibration of a shafting system is generally simplified as a positive and diagonal matrix and the stiffness matrix may be different based on different connected forms, then the typical forms of M and K are as follows:
M = m 1 0 0 0 0 m 2 0 0 0 0 0 m n .
K = [ k 12 + k 13 + + k 1 n k 12 k 13 k 1 n k 12 k 12 + k 23 + + k 2 n k 23 k 2 n k 13 k 23 k 13 + k 23 + + k 3 , n k 3 , n k 1 , n k 2 , n k 3 , n k 1 , n + k 2 , n + + k 3 , n ] .
According to the connected form of the torsional vibration of the shafting system, let q 1 = 1 , 1 , , 1 T ; then, it can be easy to obtained as follows:
K q 1 = 0 , 0 , , 0 T .
Since K = D J D , Equation (18) becomes
D J D q 1 = D J D 1 , 1 , , 1 T = 0 , 0 , , 0
where D = M 1 / 2 = diag m 1 , m 2 , , m n . Then, Equation (19) is rewritten as
J m 1 , m 2 , , m n T = 0 , 0 , , 0 T .
For the torsional vibration of a shafting, the final problem becomes finding a special solution of Homogeneous Ordinary Differential Equation (20), which is denoted by q = q 1 , q 2 , q 3 , q n . Since such a solution remains a solution of Equation (20) when multiplied by a non-zero factor, the solution of Equation (20) can be expressed as
m 1 , m 2 , , m n T = e q
where e is a proportionality constant.
As discussed in Section 2.3, the matrix J = Z can be determined from the gradient flow method, which meets the requirement of containing the desired eigenvalues. Then, it can be substituted into Equation (20) to get q, and mi (i = 1, 2, …., n) can be found from Equation (21) based on q and e. Factor e is apparently not unique and depends on how a system mass or system masses are specified. For example, it can be assumed that m 1 is kept, and it can be known that e = m 1 / q 1 . Then, all the other masses of the modified mass matrix M can be computed by substituting e into Equation (21). Finally, the stiffness matrix can be obtained from equation K = D J D , in which D = diag m 1 , m 2 , , m n . So, the modified mass matrix M and stiffness matrix K are both obtained, which form a new system with those desired eigenvalues and those original eigenvalues.

3. Torsional Vibration Control and Analysis of a Marine Diesel Engine Propulsion System

As an important part of a marine power plant, the reliability and stability of ship propulsion shafting is very important to ship operation. A diesel engine propulsion system is the most common and typical shafting system in the shipping industry. During its operation, it is subject to various periodic loads, such as the excitation of the main engine and the excitation of the propeller, which can cause serious torsional vibration of the shafting and even lead to serious shaft failure. Therefore, the effective control of torsional vibration of the diesel engine propulsion system is of great engineering significance. In this chapter, a typical marine diesel engine propulsion system is taken as an example to study the partial eigenvalue assignment for the torsional vibration control of the shaft propulsion system at different operation conditions.

3.1. The Model Setup

The model studied in this paper is based on a typical marine diesel engine propulsion system mainly consisting of a vibration absorber, engine cylinders, a flywheel, a transmission shaft, a coupling, and a propeller, as shown in Figure 1.
A lumped parameter model is widely used in the calculation of torsional vibration of shafting, which includes three basic components: the homogeneous rigid disk element, the inertia-less damping element, and the inertia-less torsional spring element. In this paper, the ship propulsion shafting system is modeled by 12 polar moments of inertia (J1, J2, J3, J4, J5, J6, J7, J8, J9, J10, J11, J12) and 11 torsional springs (K12, K23, K34, K45, K56, K67, K78,K89, K9,10, K10,11, K11,12) with damping neglected, as shown in Figure 2. The model parameter values (stiffness values and inertial values) are calculated based on a real shafting system of a ship and are given in Table 1, and the natural frequencies of the original system are given in Table 2. Damping is small and thus is neglected.

3.2. Desired Modifications under Two Different Operating Conditions

In order to better illustrate that this method can achieve different frequency assignment schemes satisfying complex physical constraints, three cases under two different operating speeds are proposed.
In the first case, it is assumed that the speed of the shafting system is 1500 rpm, so that 25 Hz is the main shafting excitation frequency. It is important to prevent the natural frequencies of the system from coinciding with this excitation frequency and its third harmonic (75 Hz). It is decided to shift the second frequency (24.96 Hz) and the fourth frequency (74.65 Hz) of the original system shown in Table 2 to 30 Hz and 90 Hz, while keeping others unchanged, as shown in the second column of Table 3. Additionally, considering an actual engineering restriction, it is assumed that the inertia and the stiffness of the diesel engine cannot be changed in structural modifications.
In the second case, the model is about an 888-rpm diesel engine system with a five-blade propeller. In this case, the excitation frequencies include the subharmonics and superharmonics of the rotation speed, and another important excitation frequency is propeller excitation (74 Hz). So, the third and fourth natural frequencies should be respectively modified to become 55 Hz and 80 Hz given in the third column of Table 3, which means that one modified frequency becomes bigger and the other becomes smaller. This case has the same physical constraint as the first case.
Compared with the second case, the third case requires that the third and fourth frequencies be shifted to smaller values (55 Hz and 70 Hz). In addition, the physical constraints are more complex compared with the above two cases. The results are listed in the fourth column of Table 3.
So far, three modification schemes for different trends of frequency changes have been given.

3.3. Results Obtained from the Gradient Flow Method

Based on the method presented in Section 2, the three cases are solved for the partial assignment of the desired natural frequencies of the shafting through modifying the value of the inertias and stiffness at specified locations aiming at meeting the requirements of practical engineering applications. Taking Case One as an example, the computational process of this method is further introduced. Firstly, according to the known inertia and stiffness parameters, Equations (1)–(6) are easy to be solved. Then, V is set to a unit matrix, and Z can be obtained from Equations (7) and (9). It is assumed that the inertia and stiffness parameters (J4, J5, J6, J7, J8 and K45, K56, K67, K78, respectively) of the engine are not changed. In other words, a subset of matrix entry locations L = {(4,5), (5,4), (5,5), (5,6), (6,5), (6,6), (6,7), (7,6), (7,7), (7,8), (8,7)} and their corresponding values δ = { K 45 J 4 J 5 ,   K 54 J 5 J 4 ,   K 55 J 5 , K 56 J 5 J 6 ,   K 65 J 6 J 5 ,   K 66 J 6 ,   K 67 J 6 J 7 ,   K 76 J 7 J 6 ,   K 77 J 7 ,   K 78 J 7 J 8 ,   K 87 J 8 J 7 }. Then, a gradient flow can be defined on the basis of Equations (10)–(14), and the ordinary differential equation can be solved to reconstruct a matrix J . Finally, based on Section 2.4, the modified mass and stiffness matrices can be obtained. The detailed results of the three cases obtained from the gradient flow algorithm are presented and analyzed in the next sections.

3.3.1. Case One (Desired Frequencies at 30 Hz and 90 Hz)

After solving the problem, the frequencies of the modified structure and the relative errors between them and the desired values are shown in Table 4, and the comparisons of the inertia values and stiffness values before and after modification are shown in Figure 3 and Figure 4. It can clearly be seen from the results shown in Table 4 that the method is excellent at achieving the desired frequencies for this shafting system: the second and third natural frequencies of the original shaft system are now shifted to 30 Hz and 90 Hz, and they make a perfect match with the desired results, which means that there is no spill-over. Furthermore, for case one, two non-adjacent natural frequencies are modified to show the universality of the method.
As shown in Figure 3, the differences between the inertia values before and after modification regarding the inertia of a diesel engine are small and can be considered to be almost identical. Figure 4 indicates that the stiffness parameters for the engine before and after modification are not exactly the same, which are similar to the inertia results; especially, stiffness K45 and stiffness K78 have some errors, but they are very small and acceptable. Therefore, the solution for modification is effective.

3.3.2. Case Two (Desired Frequencies at 55 Hz and 80 Hz)

In this case, the targets are that the third and fourth frequencies need to be changed to 55 Hz and 80 Hz while keeping the other frequencies unchanged. From Table 5, it can be shown that the third and fourth frequencies are successfully shifted to 55 Hz and 80 Hz, respectively. Furthermore, all the relative errors between the modified frequencies and the desired frequencies are equal to 0, except that the error of the 10th natural frequency is 0.11%. This means that all the modified natural frequencies are the same as the expected values.
In Figure 5 and Figure 6, the obtained results of modified inertias and stiffness are shown respectively. The arrows in the figures indicate that the values of the inertias and stiffness are almost the same before and after modification.

3.3.3. Case Three (Desired Frequencies at 55 Hz and 70 Hz)

Compared with the above two cases, the third and fourth frequencies of Case Three are required to become 55 Hz and 70 Hz. As shown in Table 6, the actually achieved values are 55.06 Hz and 70.04 Hz, respectively. Unlike the previous two cases, the physical constraints in the third case become more stringent, which requires the inertias and stiffness of both the engine and the propeller to remain unchanged. Figure 7 shows that the inertias that are not allowed to be modified are indeed kept unchanged. As shown in Figure 8, the stiffness parameters for the engine before and after modification are not exactly the same, which is similar to the inertia parameters; especially, stiffness K45, K56, and K67 have some errors. However, they are very small and acceptable. Therefore, the solution for modification is effective. This case further illustrates the universality of the method.
Through the above analysis and discussions, it is easy to see that the gradient flow method can achieve a partial assignment of natural frequencies by means of modifying the values of inertias and stiffness at certain locations while keeping the values of inertias and stiffness unchanged at other locations.

4. Torsional Vibration Control and Analysis of a Diesel Generator Set

Symmetric branched structures are widely used in engineering. In the marine engineering field, a symmetrical branched structure is widely used in a gear box system, a dual locomotive parallel system, and a diesel engine system with symmetrical auxiliary equipment. There are two main characteristics regarding the torsional vibration of symmetric branched systems: one is the phenomenon of close natural frequencies, and the other is the phenomenon of synchronous vibration and asynchronous vibration in two symmetric branched shafts. So, the torsional vibration of a branched shafting is more complex than that of a normal straight shafting. In order to better analyze and understand the torsional vibration characteristics of a branched shafting, many scholars [41,42,43,44,45] have studied the analysis methods for the torsional vibration of different branched systems. However, an effective control method for the torsional vibration of branched shafting has rarely been reported. In this section, from the point of view of structural modification, an effective solution is given to the problem of close frequencies of branched structures based on the method introduced in Section 3.

4.1. The Model Setup

In this section, a diesel generator set is taken as the research object, which consists of a vibration absorber, two pumps, gears, engine cylinders, a flywheel, a transmission shaft, a coupling, and a generator, as shown in Figure 9.
A lumped parameter model for various parts of the shafting is established in Figure 10. It is composed of homogeneous rigid disc elements and massless torsion spring elements, with damping neglected. The model parameter values (stiffness values and inertial values) are calculated based on a real diesel generator set and are given in Table 7, and the natural frequencies of the original system are given in Table 8.

4.2. The Aim of Modification

From Table 8, it is found that the system has close eigenvalues, such as the 10th and 11th natural frequencies. The corresponding modes are given in Figure 11. From Figure 11, it can be concluded that the mode shape corresponding to the 10th natural frequency represents synchronous vibration, while the mode shape corresponding to the 11th natural frequency is asynchronous vibration. It is also clear that the modal amplitudes of inertias J9 and J11 are much larger than those of other inertias when the shafting is vibration in either of these two modes. As a result, once resonance occurs at these two frequencies, the equipment represented by these inertias can undergo large vibration and thus probably serious damage. Therefore, it is necessary to avoid these resonances.
In this case, it is assumed that the speed of the shafting system is 948 r/m, so that 15.8 Hz is the main shafting excitation frequency. It is important to prevent the excitation frequencies from coinciding with the natural frequencies of the system. So, the 17th harmonic frequency (268.6 Hz) needs to be avoided, and for this purpose, the 10th and 11th natural frequencies of the original system will be moved to 262.5 Hz and 263.0 Hz while keeping the others unchanged, as shown in Figure 12 corresponding to the 10th frequency (267.20 Hz) and the 11th frequency (267.77 Hz) shown in Table 8.

4.3. Results Obtained from the Gradient Flow Method

In this case, it is expected to change the original 10th natural frequency (267.20 Hz) and 11th natural frequency (267.77 Hz) to the corresponding desired frequencies of 262.50 Hz and 263.00 Hz, respectively. From the results given in Table 9, the 10th natural frequency and 11th natural frequency are shifted to 262.50 Hz and 262.99 Hz, respectively, and all the relative errors between the modified frequencies and the desired frequencies are equal to zero. In addition, according to the actual engineering requirements, it is impractical or unrealistic to change the parameters of diesel engines and generators after their designs have been finished or real machines have been produced. For this reason, it is assumed in this example that the relevant inertia (J3, J4, J5, J6, J8) and stiffness (K34, K45, K56, K78) parameters of the diesel engine and generator remain unchanged. From Table 10, it can be shown that the maximum changes of inertia and stiffness are 0.02 (kg∙m2) and 0.04 (105 N/m), respectively. These two small changes have almost no impact on the vibration characteristics; thus, they can be neglected in practice. Through the above analysis, it can be concluded that the torsional vibration of the shafting associated with the close frequencies is effectively controlled by shifting the relevant frequencies under the constraints of practical engineering operability. This example provides a useful scheme for the torsional vibration control of complex branched shafting.

5. Conclusions

With the development of the shipping industry, ship propulsion systems have become more and more complex, which makes the analysis and control of the torsional vibration of shafting systems more difficult. In order to achieve torsional vibration control, suitable structural modifications are an effective solution. However, in practical applications, some modifications of the actual shafting systems are limited by many physical constraints. In this paper, the gradient flow method is applied to determine the required modifications at any location with arbitrary eigenvalues, which are especially useful for structures or equipment that have locations that are not allowed to change, such as those in engineering applications.
In order to demonstrate the effectiveness of this method, two examples of complex ship engine systems are studied, which are a typical ship propulsion shaft system and a diesel generator set. In the first example, the natural frequencies of the original structure are shifted in three different scenarios under two different operating conditions and different physical constraints. In the second example, there are two very close frequencies due to the complex symmetric branches within the system. The partial assignment of these frequencies is achieved for the first time through structural modifications. The success of this method in solving these two examples shows its wide applicability and strength in controlling torsional vibration. It is expected to be applicable to linear passive vibration control for other machines and structures.

Author Contributions

Methodology, M.C. and H.O.; software and analysis, S.L. and D.W.; writing—original draft preparation, M.C.; writing—review and editing, M.C., H.O., and W.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research work is funded by the Fundamental Research Funds for Central Universities (No.3072019CFJ0303) and “National Special Support Program for High-Level Talents”.

Acknowledgments

Most of the research is carried out by the first author during her visit to the Universtiy of Liverpool under the supervision of the second author. Discussions wtih Shike Zhang of University of Liverpool have been helpful.

Conflicts of Interest

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

Notation

M 0 mass matrix V an n × n orthogonal matrix
K 0 stiffness matrix Q ξ all matrices that have the same eigenvalues as ξ
q G n displacement vector L some matrix entry locations
q ¨ G n acceleration vector δ the desired values corresponding to L
J 0 a real matrix with the same eigenvalues as the original system M 0 , K 0 N L , δ all matrices with the prescribed entries at the desired locations
J s a real symmetric matrix with all expected eigenvalues Z any given matrix in Q ξ
λ 1 , λ 2 , , λ n n system eigenvalues Z L the same as Z , except those at the matrix entry locations that do not belong to L are set identically zero
Λ 1 original system eigenvalues L s all locations that do not belong to L
μ 1 the desired eigenvalues P Z the projection of any matrix Z onto the affine subspace N L , δ
Λ 2 all the eigenvalues that remain unchanged C L a constant matrix in N L , δ with zero entries at all locations corresponding to L s and known entries at all other locations
X 1 the mass-normalized eigenvector matrix corresponding to Λ 1 f Z A function defined to minimize the distance between Z and N L , δ
U 1 and U 2 the normalized eigenvector matrices h V A function defined the minimization with objective function f Z in terms of V
ξ G n × n a real symmetric matrix h the gradient of objective function h

References

  1. Baldwin, J.F.; Hutton, S.G. Natural models of modified structures. AIAA 1985, 23, 1737–1743. [Google Scholar] [CrossRef]
  2. Chen, S.H.; Yang, X.W.; Lian, H.D. Comparison of several eigenvalue reanalysis methods for modified structures. Struct. Multidiscip. Optim. 2000, 20, 253–259. [Google Scholar] [CrossRef]
  3. Weissenburger, J.T. Effects of local modification on the vibration characteristics of linear systems. J. Appl. Mech. 1968, 35, 327–332. [Google Scholar] [CrossRef]
  4. Pomazal, R.J.; Snyder, V.W. Local modifications of damped linear systems. AIAA 1971, 9, 2216–2221. [Google Scholar] [CrossRef]
  5. Ram, Y.M.; Blech, J.J. The dynamic behaviour of a vibratory system after modification. J. Sound Vib. 1991, 150, 357–370. [Google Scholar] [CrossRef]
  6. Wang, B.P.; Pilkey, W.D. Eigenvalue reanalysis of locally modified structures using a generalized Rayleigh’s method. AIAA 1986, 24, 938–990. [Google Scholar] [CrossRef]
  7. To, W.M.; Ewins, D.J. Structural modification analysis using Rayleigh quotient iteration. Int. J. Mech. Sci. 1990, 32, 169–179. [Google Scholar] [CrossRef]
  8. Mottershead, J.E.; Link, M.; Friswell, M.I. The sensitivity method in finite element model updating: A tutorial. Mech. Syst. Signal Process. 2011, 25, 2275–2296. [Google Scholar] [CrossRef]
  9. Sestieri, A. Structural dynamic modification. Sadhana. Acad. Proc. Eng. Sci. 2000, 25, 247–260. [Google Scholar]
  10. Richiedei, D.; Trevisani, A.; Zanardo, G. A constrained convex approach to modal design optimization of vibrating systems. J. Mech. Design 2011, 133, 061011. [Google Scholar] [CrossRef]
  11. Farahani, K.; Bahai, H. An inverse strategy for relocation of eigenfrequencies in structural design. Part I: First order approximate solutions. J. Sound Vib. 2004, 274, 481–505. [Google Scholar] [CrossRef]
  12. Farahani, H.; Aryana, F. Design optimisation of structures vibration behaviour using first order approximation and local modification. Comput. Struct. 2002, 80, 1955–1964. [Google Scholar]
  13. Bucher, I.; Braun, S. The structural modification inverse problem: An exact solution. Mech. Syst. Signal Process. 1993, 7, 217–238. [Google Scholar] [CrossRef]
  14. Li, T.; He, J.M. Local structural modification using mass and stiffness changes. Eng. Struct. 1999, 21, 1028–1037. [Google Scholar] [CrossRef]
  15. Park, Y.H.; Park, Y.S. Structural modification based on measured frequency response functions: An exact eigenproperties reallocation. J. Sound Vib. 2000, 237, 411–426. [Google Scholar] [CrossRef]
  16. Mottershead, J.E. Structural modification for the assignment of zeros using measured receptances. J. Appl. Mech. 2001, 68, 791–798. [Google Scholar] [CrossRef]
  17. Hu, Y.L.; Michael, Z.Q.; Smith, M.C. Natural frequency assignment for mass-chain system with inerters. Mech. Syst. Signal Process. 2018, 108, 126–139. [Google Scholar] [CrossRef] [Green Version]
  18. Yuan, Y.X.; Dai, H. An inverse problem for undamped gyroscopic systems. J. Comput. Appl. Math. 2012, 236, 2574–2581. [Google Scholar] [CrossRef] [Green Version]
  19. Tsai, S.H.; Ouyang, H.; Chang, J.Y. Inverse structural modifications of a geared rotor-bearing system for frequency assignment using measured receptances. Mech. Syst. Signal Process. 2018, 110, 59–72. [Google Scholar] [CrossRef]
  20. Liu, Z.H.; Li, W.Y.; Ouyang, H. Structural modifications for torsional vibration control of shafting systems based on torsional receptances. Shock Vib. 2016. [Google Scholar] [CrossRef]
  21. Zhang, J.F.; Ye, J.P.; Ouyang, H.; Yin, X. An explicit formula of perturbating stiffness matrix for partial natural frequency assignment using static output feedback. J. Low Freq. Noise Vib. Act. 2018, 37, 1045–1052. [Google Scholar] [CrossRef] [Green Version]
  22. Singh, K.V.; Ling, X.X. Active control of viscoelastic systems by the method of receptance. J. Vib. Acoust. 2018, 140, 1–6. [Google Scholar] [CrossRef]
  23. Araujo, J.M.; Carlos, C.E.T.; Goncalves, L.M.G. Robustness of the quadratic partial eigenvalue assignment using spectrum sensitivities for state and derivative feedback designs. J. Low Freq. Noise Vib. Act. 2018, 37, 253–268. [Google Scholar] [CrossRef]
  24. Ram, Y.M.; Mottershead, J.E.; Tehrani, M.G. Partial pole placement with time delay in structures using the receptance and the system matrices. Linear Algebra Appl. 2011, 434, 1689–1696. [Google Scholar] [CrossRef] [Green Version]
  25. Wang, X.T.; Zhang, L. Partial eigenvalue assignment with time delay in high order system using the receptance. Linear Algebra Appl. 2017, 523, 335–345. [Google Scholar] [CrossRef]
  26. Bai, Z.J.; Chen, M.X.; Datta, B.N. Minimum norm partial quadratic eigenvalue assignment with time delay in vibrating structures using the receptance and the system matrices. J. Sound Vib. 2013, 332, 780–794. [Google Scholar] [CrossRef] [Green Version]
  27. Bai, Z.J.; Yang, J.K.; Datta, B.N. Robust partial quadratic eigenvalue assignment with time delay using the receptance and the system matrices. J. Sound Vib. 2016, 384, 1–14. [Google Scholar] [CrossRef]
  28. Bai, Z.J.; Lu, M.; Wan, Q.Y. Minimum norm partial quadratic eigenvalue assignment for vibrating structures using receptances and system matrices. Mech. Syst. Signal Process. 2018, 112, 265–279. [Google Scholar] [CrossRef]
  29. Liu, H.; He, B.X.; Chen, X.P. Minimum norm partial quadratic eigenvalue assignment for vibrating structures using receptance method. Mech. Syst. Signal Process. 2019, 123, 131–142. [Google Scholar] [CrossRef]
  30. Bai, Z.J.; Chen, M.X.; Yang, J.K. A multi-step hybrid method for multi-input partial quadratic eigenvalue assignment with time delay. Linear Algebra Appl. 2012, 437, 1658–1669. [Google Scholar] [CrossRef]
  31. Ram, Y.M.; Mottershead, J.E. Multiple-input active vibration control by partial pole placement using the method of receptances. Mech. Syst. Signal Process. 2013, 40, 727–735. [Google Scholar] [CrossRef]
  32. Zhang, L.; Wang, X.T. Partial eigenvalue assignment of high order systems by multi-input control. Mech. Syst. Signal Process. 2014, 42, 129–136. [Google Scholar] [CrossRef]
  33. Zhang, L. Multi-input partial eigenvalue assignment for high order control systems with time delay. Mech. Syst. Signal Process. 2016, 72–73, 376–382. [Google Scholar] [CrossRef]
  34. Zhang, L.; Yu, F.; Wang, X.T. An algorithm of partial eigenstructure assignment for high-order systems. Math. Methods Appl. Sci. 2018, 41, 6070–6079. [Google Scholar] [CrossRef]
  35. Li, S.; Yu, D.Z. Simulation study of active control of vibro-acoustic response by sound pressure feedback using modal pole assignment strategy. J. Low Freq. Noise Vib. Act. 2016, 35, 167–186. [Google Scholar] [CrossRef] [Green Version]
  36. Yuan, Y.X. Stiffness matrix modification with vibration test data by displacement feedback technique. J. Vibroeng. 2017, 19, 5149–5160. [Google Scholar] [CrossRef]
  37. Araujo, J.M.; Santos, T. A multiplicative eigenvalues perturbation and its application to natural frequency assignment in undamped second-order systems. Proc. Inst. Mech. Eng. Part I 2018, 232, 963–997. [Google Scholar] [CrossRef]
  38. Ouyang, H.; Zhang, J.F. Passive modifications for partial assignment of natural frequencies of mass-spring systems. Mech. Syst. Signal Process. 2015, 50–51, 214–226. [Google Scholar] [CrossRef] [Green Version]
  39. Belotti, R.; Ouyang, H.; Richiedei, D. A new method of passive modifications for partial frequency assignment of general structures. Mech. Syst. Signal Process. 2018, 99, 586–599. [Google Scholar] [CrossRef]
  40. Chu, M.T.; Diele, F.; Sgura, I. Gradient flow methods for matrix completion with prescribed eigenvalues. Linear Algebra Appl. 2004, 379, 85–112. [Google Scholar] [CrossRef] [Green Version]
  41. Huang, Y.M.; Hong, C.D. Extended transfer matrix method with complex numbers for branched torsional systems. J. Vib. Control 2001, 7, 155–166. [Google Scholar] [CrossRef]
  42. Zou, C.P.; Chen, D.S.; Hua, H.X. Torsional vibration analysis of complicated multi-branched shafting systems by modal synthesis method. J. Vib. Acoust. 2003, 125, 317–323. [Google Scholar] [CrossRef]
  43. Wu, J.S.; Chen, C.H. Torsional vibration analysis of gear-branched systems by finite element method. J. Sound Vib. 2001, 240, 159–182. [Google Scholar] [CrossRef]
  44. Liu, Z.S.; Chen, S.H.; Xu, T. Derivatives of eigenvalues for torsional vibration of geared shaft systems. J. Vib. Acoust. 1993, 115, 277–279. [Google Scholar]
  45. Crowther, A.R.; Zhang, N. Torsional finite elements and nonlinear numerical modelling in vehicle powertrain dynamics. J. Sound Vib. 2005, 284, 825–849. [Google Scholar] [CrossRef]
Figure 1. Schematic diagram of a marine diesel engine propulsion system.
Figure 1. Schematic diagram of a marine diesel engine propulsion system.
Applsci 10 00147 g001
Figure 2. The lumped parameter model of a marine diesel engine propulsion system.
Figure 2. The lumped parameter model of a marine diesel engine propulsion system.
Applsci 10 00147 g002
Figure 3. Comparisons of inertia values before and after modification in Case One where the desired frequencies are at 30 Hz and 90 Hz.
Figure 3. Comparisons of inertia values before and after modification in Case One where the desired frequencies are at 30 Hz and 90 Hz.
Applsci 10 00147 g003
Figure 4. Comparisons of stiffness values before and after modification in Case One where the desired frequencies are at 30 Hz and 90 Hz.
Figure 4. Comparisons of stiffness values before and after modification in Case One where the desired frequencies are at 30 Hz and 90 Hz.
Applsci 10 00147 g004
Figure 5. Comparisons of inertia values before and after modification in Case Two where the desired frequencies are at 55 Hz and 80 Hz.
Figure 5. Comparisons of inertia values before and after modification in Case Two where the desired frequencies are at 55 Hz and 80 Hz.
Applsci 10 00147 g005
Figure 6. Comparisons of stiffness values before and after modification in Case Two where the desired frequencies are at 55 Hz and 80 Hz.
Figure 6. Comparisons of stiffness values before and after modification in Case Two where the desired frequencies are at 55 Hz and 80 Hz.
Applsci 10 00147 g006
Figure 7. Comparisons of inertia values before and after modification in Case Three where the desired frequencies are at 55 Hz and 70 Hz.
Figure 7. Comparisons of inertia values before and after modification in Case Three where the desired frequencies are at 55 Hz and 70 Hz.
Applsci 10 00147 g007
Figure 8. Comparisons of stiffness values before and after modification in Case Three where the desired frequencies are at 55 Hz and 70 Hz.
Figure 8. Comparisons of stiffness values before and after modification in Case Three where the desired frequencies are at 55 Hz and 70 Hz.
Applsci 10 00147 g008
Figure 9. Schematic diagram of a diesel generator set.
Figure 9. Schematic diagram of a diesel generator set.
Applsci 10 00147 g009
Figure 10. The lumped parameter model of a diesel generator set.
Figure 10. The lumped parameter model of a diesel generator set.
Applsci 10 00147 g010
Figure 11. The mode shapes corresponding to the 10th order and 11th order frequencies.
Figure 11. The mode shapes corresponding to the 10th order and 11th order frequencies.
Applsci 10 00147 g011
Figure 12. Frequency response function with desired frequencies modification.
Figure 12. Frequency response function with desired frequencies modification.
Applsci 10 00147 g012
Table 1. System parameters.
Table 1. System parameters.
Parameter Value   [ kg m 2 ] Parameter Value   [ 10 5 N / m ]
J15K1210
J22K23200
J31K34150
J43K45100
J53K56100
J63K67100
J73K78100
J83K89200
J940K9,107
J105.5K10,114
J113.5K11,1250
J128
Table 2. Frequencies of modes of the original system.
Table 2. Frequencies of modes of the original system.
Mode number123456
f [Hz]024.9657.5274.65108.25232.69
Mode number789101112
f [Hz]234.95363.90467.27538.75577.991046.91
Table 3. Desired system frequencies and physical constraints.
Table 3. Desired system frequencies and physical constraints.
The First CaseThe Second CaseThe Third Case
Mode number233434
f [Hz]309055805570
Physical constraintsThe inertias and stiffness of the diesel engine remain unchangedThe inertias and stiffness of the diesel engine remain unchangedThe inertias and stiffness of the diesel engine and propeller remain unchanged
Table 4. Frequencies of modes of the modified system and the relative errors.
Table 4. Frequencies of modes of the modified system and the relative errors.
Mode number123456
f [Hz]029.9957.5390.08108.24232.68
Error (%)000000
Mode number789101112
f [Hz]234.86363.89467.46539.03578.251046.85
Error (%)000000
Table 5. Frequencies of modes of the modified system and the relative errors.
Table 5. Frequencies of modes of the modified system and the relative errors.
Mode number123456
f [Hz]024.9655.0080.00108.20232.69
Error (%)000000
Mode number789101112
f [Hz]234.81363.90467.69539.32578.531048.63
Error (%)0000.1100
Table 6. Frequencies of modes of the modified system and the relative errors.
Table 6. Frequencies of modes of the modified system and the relative errors.
Mode number123456
f [Hz]024.9755.0670.04108.56232.78
Error (%)000000
Mode number789101112
f [Hz]235.20364.20467.68539.18578.381051.33
Error (%)000000
Table 7. System parameters.
Table 7. System parameters.
Parameter Value   [ k g m 2 ] Parameter Value   [ 10 5 N / m ]
J115K12100
J22.8K2335
J39K3425
J49K4525
J59K5625
J69K67100
J7220K784
J8150K892
J90.6K9,1012.6
J102.8K10,112
J110.6K11,1212.6
J122.8
Table 8. Frequencies of modes of the original system.
Table 8. Frequencies of modes of the original system.
Mode number123456
f [Hz]010.0922.4236.6441.5474.87
Mode number789101112
f [Hz]120.97154.56195.72267.20267.77375.48
Table 9. Frequencies of modes of the modified system and the relative errors.
Table 9. Frequencies of modes of the modified system and the relative errors.
Mode number123456
f [Hz]010.0922.4236.6241.5374.86
Error (%)000000
Mode number789101112
f [Hz]120.94154.58195.71262.50262.99375.49
Error (%)000000
Table 10. Modified system parameter comparison.
Table 10. Modified system parameter comparison.
ParameterOriginal Value
[kg∙m2]
Obtained Value
[kg∙m2]
Changed
[kg∙m2]
ParameterOriginal Value
[105 N/m]
Obtained Value
[105 N/m]
Changed
[105 N/m]
J11514.80−0.2K12100109.299.29
J22.83.000.2K233533.78−1.22
J398.99−0.01K342524.99−0.01
J499.000K452525.010.01
J599.010.01K562525.040.04
J699.020.02K67100100.220.22
J7220220.140.14K7844.000
J8150150.020.02K8921.99−0.01
J90.60.600K9,1012.612.11−0.49
J102.82.77−0.03K10,1121.99−0.01
J110.60.600K11,1212.612.13−0.47
J122.82.78−0.02
(J3, J4, J5, J6, and J8 are unchanged, and K34, K45, K56, and K78 are unchanged).

Share and Cite

MDPI and ACS Style

Chen, M.; Ouyang, H.; Li, W.; Wang, D.; Liu, S. Partial Frequency Assignment for Torsional Vibration Control of Complex Marine Propulsion Shafting Systems. Appl. Sci. 2020, 10, 147. https://doi.org/10.3390/app10010147

AMA Style

Chen M, Ouyang H, Li W, Wang D, Liu S. Partial Frequency Assignment for Torsional Vibration Control of Complex Marine Propulsion Shafting Systems. Applied Sciences. 2020; 10(1):147. https://doi.org/10.3390/app10010147

Chicago/Turabian Style

Chen, Meilong, Huajiang Ouyang, Wanyou Li, Donghua Wang, and Siyuan Liu. 2020. "Partial Frequency Assignment for Torsional Vibration Control of Complex Marine Propulsion Shafting Systems" Applied Sciences 10, no. 1: 147. https://doi.org/10.3390/app10010147

APA Style

Chen, M., Ouyang, H., Li, W., Wang, D., & Liu, S. (2020). Partial Frequency Assignment for Torsional Vibration Control of Complex Marine Propulsion Shafting Systems. Applied Sciences, 10(1), 147. https://doi.org/10.3390/app10010147

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop