Research Article  Open Access
Jun Ren, Qianghao Zhang, "Structural Reanalysis Based on FRFs Using Sherman–Morrison–Woodbury Formula", Shock and Vibration, vol. 2020, Article ID 8212730, 12 pages, 2020. https://doi.org/10.1155/2020/8212730
Structural Reanalysis Based on FRFs Using Sherman–Morrison–Woodbury Formula
Abstract
Structural dynamic modification is a popular approach to obtain desire frequencies and dynamic characteristics. It has been observed that reanalyzing the modified structure usually involves complicated calculations when modifications are concerned with numerous degrees of freedom (DOFs), especially adding substructures to these DOFs. This paper proposed a method to reanalyze the frequency response functions (FRFs) of structures with multiple coordinates modifications. Two different cases are taken into consideration in the modifications, including adding (or decreasing) masses, stiffness, and damping, as well as adding springmass substructures, which makes the method more practical. This method is developed by employing Sherman–Morrison and Woodbury (SMW) formula based on the FRFs related to the modifications coordinates of the original system. The advantage of this method is that neither a physical model nor a modal model is required; instead, it needs only the FRFs, which can be directly measured by experimental modal testing. Another salient feature of this proposed strategy is that the FRFs of the modified structure can be calculated in only one step. Validation of this proposed method is demonstrated using various numerical examples. It is shown that the method is very effective and can be considered for real applications.
1. Introduction
Structural dynamic modification has been applied widely in practical engineering. In many engineering cases, structural dynamic modifications are used to obtain better dynamic characteristics of structures. Generally speaking, structural dynamic modification refers to a method to get certain structural dynamic characteristics by changing the local physical parameters (masses, stiffness, and damping), such as the need to avoid resonance or assignment of certain frequencies on desired locations. Structural dynamic modification is an economical and effective means of improving the dynamic characteristics of mechanical structures. This approach is widely used in aerospace, marine, automotive, civil engineering, bridge, and machinery industries.
The problems involved in structural dynamic modification can be divided into two categories: “forward problem” and “inverse problem.” The former aims to predict exact change to the structure’s dynamic properties when known modifications are made at a given location [1, 2]. The latter mainly studies how to modify an existing structure in order to achieve the expected dynamic characteristics (such as natural frequency and mode shape) [3–6]. The forward structural dynamic modification is also called structural reanalysis in early studies. In structural reanalysis, the dynamic effects of a modification on a structure are treated as an analysis problem involving the known dynamic properties of the original structure rather than a complete reanalysis of the modified structure [7]. This approach could avoid the blindness of the design, which reduces the design cost and has practical engineering application value.
Many works and literatures had conducted a comprehensive analysis of reanalysis methods. Reanalysis methods can be generally divided into two categories: direct methods and approximate methods. The approximate methods can be divided into three categories [8]: global approximation, local approximation, and combined approximation.
For the development of direct methods, early studies of structural reanalysis were reviewed and summarized by Baldwin and Hutton [7]. Several approaches such as Rayleigh quotient [9], sensitivity analysis [10], and perturbation approach [11] were used to address forward modification problems without a complete reanalysis of the whole structure.
This important issue has also been extensively discussed in recent years, and part of the relevant literature is summarized as follows. This problem was explored in [12] using developed successive matrix inverse method based on symmetry of corresponding stiffness matrix after constraint modification of the boundary. The numerical examples show that this method could quickly give accurate reanalysis results. In the same year, Liu et al. [13] proposed an approach for structural static reanalysis with unchanged number of degrees of freedom. This approach was based on a new preconditioner constructed by updating the Cholesky factorization of the original stiffness matrix, which could achieve fast convergence and accurate results. After two years, Song et al. [14] suggested using a direct reanalysis algorithm based on finding updated triangular factorization in sparse matrix solution to solve this problem. This algorithm is suitable for local modification, and the examples show that the algorithm improves reanalysis efficiency significantly, especially for highrank structural modification. Later on, a “cheap” algorithm, named independent coefficients method, was put forward to reanalyze structures with local modification, which leads to a lowrank change in the stiffness matrix [15]. Considering that previous work [13] suffers limitation of a structure as added degrees of freedom, a new and efficient reanalysis method [16] had been proposed by the same author. Another advantage of this method is that the Cholesky factorization of the stiffness matrix of the modified structure can be used as the initial information for reanalysis when the structure is further modified. In the work of [17], the issue was suggested using a nonlinear reanalysis method based on structural modification of residual incremental approximations. In contrast to other existing nonlinear reanalysis methods, which were based on the evaluation of changed stiffness matrices, only residual vectors need to be computed and stored. Kim and Eun [18] studied coupling and recoupling reanalysis methods. These methods were performed by using the concept of compatibility conditions at interface nodes between the substructures or between the original structure and the substructures. The most recent theory proposed in [19] offered a new method for freevibration reanalysis after structural topological modifications with added degrees of freedom. The implementation of this approach involved only LDL^{T} factorization of shifted substiffness matrices corresponding to the newly added DOFs, and the proposed method consists of matrixmatrix operations.
As for approximate methods, the combined approximation method is an effective solution method that combines the high efficiency of the local approximation method with the high quality of the global approximation method. The solution process of the method is based on the results of an accurate singlepoint analysis, and it is also a reanalysis method that has developed very fast in recent years. The original purpose of studying the CA method is to accelerate the optimization design. CA method has been widely used in many fields such as linear or nonlinear static analysis, dynamic analysis, modal analysis, and sensitivity analysis [8].
For the combined approximation method of reanalysis, Kirsch [20] first used forward and reverse substitution calculations to calculate the terms used as the basis vector binomial sequence in the CA method solving process. A new set of uncoupled basis vectors are generated and normalized by using the Gram–Schmidt orthogonalization process. This method can achieve an effective and accurate approximation for very large design changes. By 2006, Kirsch et al. [21] used combined approximation method to overcome the repeated eigenproblem solution of nonlinear dynamic reanalysis and solved its main problem. The method is based on the integration of several concepts and methods, including the basis of matrix factorization, series expansion, and reduction. In order to solve the frequencyconstrained structural optimization problem, Zuo et al. [22] proposed an adaptive eigenvalue reanalysis method based on genetic algorithm for structural optimization. The modified impulse analysis method is a combination approximation method from Kirsch, and it has a high level for repeated eigenvalue problems accuracy. Considered to integrate the Kirsch’s method into the result optimization process, a new adaptive method [23] that used the K condition number to determine the minimum number of basis vectors was proposed. Besides, on the reanalysis of sensitivity, Zuo et al. [24] proposes a new method for arbitrarily changing static displacement sensitive design variables. This method uses Taylor series expansion to approximate the current displacement of the modified sensitivity equation and then solves the direct sensitivity equation by a combined approximation method. One year later, the same author [25] conducted a sensitivity analysis of eigenvalues and eigenvectors using a combination approximation method, and the eigenvectors were solved by the Nelson method. This method can greatly improve the efficiency of sensitivity analysis and can accelerate the gradientbased structural optimization constraints with frequency and mode shape.
All these researches mentioned above are primarily based on the physical model, which requires the knowledge of mass, stiffness, and damping matrices. In practical engineering, however, these parameters matrices of vibration system structures are not easy to obtain. This is because the structure to be modified is usually a complex structure with multiple DOFs. Furthermore, in the process of reanalysis calculation, it is usually not very difficult to solve the singleelement change problem. When it comes to multipleelement change structures, however, the calculation becomes more complicated. Therefore, these problems limit the application of the above methods to some extent.
This paper proposed a method for reanalyzing FRFs of the modified structure. This method is developed by employing SMW formula [26, 27] based on the FRFs related to modifications coordinates of the original system. The advantage of this FRFsbased method is that the FRFs can be directly measured by model testing, without knowledge of the system matrices M, C, K, which are usually unavailable in practical engineering. Another salient feature of this proposed strategy is that the FRFs of a structure with multipleelements changing can be calculated in only one step, which improves the efficiency of reanalysis.
2. Theoretical Development
The equation of motion of a freevibration damping multidegreeoffreedom system can be expressed as
The dynamic stiffness matrix of original structure can be given bywhere Z is the dynamic stiffness matrix of original structure; represents the frequency variable and .
2.1. Adding Masses, Stiffness, and Damping to the Original Model Structure
It is assumed that the local modification of the structure involves n coordinate points, labeled 1, 2, …, n. The additional masses, stiffnesses, and damping at these points are denoted as Δm_{1}, Δm_{2}, …, Δm_{n}; Δk_{1}, Δk_{2}, …, Δk_{n}; Δc_{1}, Δc_{2}, …, Δc_{n}, respectively. The values of these additional stiffnesses, additional masses, and additional damping can be positive or negative. If the value is positive, then additional parameters are added to the original structure, while a negative value means that additional parameters are reduced from the original structure.
The above additional masses, stiffnesses, and damping can be expressed as the diagonal matrix of (3), (4), and (5), respectively.
After adding additional masses ΔM, additional stiffnesses ΔK, and additional damping ΔC,where Z is the structure dynamic stiffness matrix after adding additional parameters to original structure.
To make ΔZ more intuitive, it is expressed by the following formula:where U_{k} is a column vector of in which element of the k_{th} row is 1, and the other elements are zero. V_{k} represents a column vector of in which element of the k_{th} row is and the other elements are zero. Then the column vectors U and V can be given by
According to SMW formula [28], the dynamic stiffness Z^{−1} of modified structure with multielement change is simply written aswhere
Since the dynamic stiffness matrix of structure and the FRFs matrix are inverse matrices of each other, thenwhere H is FRFs matrix of original structure; H is FRFs matrix of the structure after adding additional masses, stiffnesses, and damping.
It is obvious that the FRFs matrix H can be calculated according to (12) when the FRFs H of original structure and added masses matrix ΔM, stiffnesses matrix ΔK, and damping matrix ΔC are known.
It should be noted that the FRFs matrices H and H mentioned above are displacement FRFs matrices. In the practical engineering cases, acceleration sensors are usually used to measure the response, so accelerance Ha is directly obtained. The relationship between H and Ha can be described as
Therefore, once accelerance matrix Ha of the original structure is obtained in practical application, the receptances H can be calculated according to (13), then substituting it into (10) and (11) to calculate receptances H of the modified structure, where H could be expressed as
Finally, the obtained H is substituted into (13) to calculate accelerances Ha of the modified structure.
2.2. Adding SpringMass Substructures to the Original Structure
In many engineering cases, the vibration system structures are designed to be immutable. When original structure does not satisfy the dynamic characteristics, adding springmass substructures to the original structural system is not a bad idea. The added springmass substructure is shown in Figure 1.
For the reason of adding a springmass substructure, the DOFs of original system has been changed, and one DOF is added to original system. The above calculated method is not suitable for this kind of condition anymore, and thus transforming the DOFs of modified structure to ones of original structure is needed [29].
The equation of motion of a linear undamped multidegreeoffreedom system can be expressed as
Assuming that the mass and stiffness modification of original system are ΔM, ΔK, then (15) can be written as
Assuming harmonic response and substituting it into (15), dynamic stiffness matrix Z and FRFs H of the system can be written as follows:
Then (16) yields
It is assumed that a springmass substructure is added at the i_{th} freedom of original structure. The mass and stiffness are Δm, Δk and the relative amplitude of vibration is Δ. Since one substructure is added to the original system structure, an extra freedom is added. The DOFs of original structure are changed from n to n + 1, and the matrices in (18) are enlarged by one row and column. Then, the equation of motion of the modified system is described by
The last row of above equation is
Substitute Δ with
Substituting (21) into (19), the right matrix of (19) can be written as follows:
Considering the i_{th} and (n + 1)_{th} row of (22), one can obtain
Then (19) can be written as follows:
It is easy to observe from (24) that when adding one springmass substructure at the i_{th} freedom of the original structure, element in the i_{th} row represented by V_{k} in the (8) is in this situation, and the other elements of V_{k} are zero; namely,
Without changing original structure of the system and adding one substructure to a certain coordinate of the system, after added one substructure, the DOFs of original system structure are changed from n to n + 1. According to a series of numerical operations, the DOFs of the system structure after adding one substructure are transformed from n + 1 to n, which could meet the calculation approach proposed in this paper.
3. Verification of the Method
In this section, three simulated examples are analyzed by the proposed method.
3.1. Numerical Experiment Setup
To verify the accuracy of the above method, a cantilever beam modal test model is constructed as shown in Figure 2. The physical parameters are shown in Table 1. The cantilever beam is discretely divided into 6 equal parts along the length direction, and 6 measuring points are evenly distributed.

As shown in Figure 3. Haorilp are the accelerances curves of original structure, where Ha_{lp} are calculated accelerances relating points l and p.
(a)
(b)
(c)
Two kinds of examples are given to prove this approach proposed in this paper: one is to directly increase or decrease masses, stiffnesses, or damping on original structure. The other is to add a springmass substructure to a certain coordinate or multiple coordinates to the original structure. What should be pointed out is that the substructures are added in the vertical direction of the cantilever. For comparison purpose, exact accelerances corresponding with modified structure are also numerically calculated.
3.2. Adding Mass, Stiffness, and Damping to the Original Model Structure
As shown in Figure 4, the additional stiffness and damping are added at coordinates 2 and 4, stiffnesses Δk_{2} and Δk_{4} are 6000 N/m and 8000 N/m, respectively, and damping Δc_{2} and Δc_{4} are 30 Ns/m and 20 Ns/m, respectively. Additional masses are added at coordinates 2, 4, and 6, respectively, and the additional masses Δm_{2}, Δm_{4} and Δm_{6} are 0.38 Kg, 0.42 Kg, and 0.4 Kg, respectively. Response points are chosen at coordinates 2, 4, and 6. Hammer impact is moving sequentially from points 2, 4 to 6 for “measuring” accelerances . The goal of this example is to calculate the accelerances .
According to (3), (4), and (5), the additional mass matrix, stiffness matrix, and damping matrix in this example are
Expression of vector V obtained by (9) can be given as follows:
According to (13), receptances H of original structure can be obtained.
The accelerances H of modified structure can be calculated by (11) and (14).
It can be seen from the above that all the FRFs H (corresponding to original structure) required for calculating H (corresponding to modified structure) are related to the modification coordinates. Then, accelerances Ha of modified structure can be directly calculated by (13).
Although nine FRFs are calculated, only three of them are shown for the sake of brevity, as can be seen in Figure 5.
(a)
(b)
(c)
As shown in Figure 5, “exact” represents the exact accelerances curves obtained directly by numerical calculation. “Measured” represents the accelerances curves obtained by the method this paper proposed. “Original” represents the accelerances curves of original structure. It is obvious from Figure 5 that, after dynamic modification to the original structure of vibration system, the accelerances curves “Measured” as a whole are in quite good agreement with those of target accelerances curve “exact”, which indicates that the results of these two methods are consistent and accuracy of the method is verified.
3.3. Adding SpringMass Substructure to the Original Structure
One and multiple springmass substructures are, respectively, added to the original structure to calculate the FRFs of the modified structure for verification of the proposed method.
3.3.1. Adding One SpringMass Substructure to the Original Structure
As shown in Figure 6, a springmass substructure is added to the original structure at coordinate 1, where the stiffness Δk_{1} of the substructure is 6000 N/m, and the mass Δm_{1} of the substructure is 0.5 kg. For the accelerances of the structure after adding this springmass substructure, 4 natural frequencies are obtained within the range of 0–70 Hz. Response points are chosen at coordinates 1, 3, and 5. Hammer impact is moving sequentially from points 1, 3 to 5 for “measuring” accelerances . The goal of this example is to calculate the nine accelerances . Similarly, only three of them are shown for the sake of simplicity, as can be seen in Figure 7.
(a)
(b)
(c)
As shown in Figure 7, as one springmass substructure is added to the original structure, the exactly calculated accelerances curves of the modified structure and the ones calculated by the proposed method completely coincide. This result indicates that, in the condition of adding one springmass substructure, modified accelerances can be accurately calculated by the proposed method.
3.3.2. Adding Three SpringMass Substructures to the Original Structure
Figure 8 is a cantilever beam model after adding substructures. Three springmass substructures are added at coordinate points 2, 4, and 6, respectively. The masses Δm_{2}, Δm_{4}, and Δm_{6} of these substructures are 0.35 kg, 0.42 kg, and 0.39 kg, respectively, and stiffness Δk_{2}, Δk_{4} and Δk_{6} of these substructures are 3000 N/m, 4000 N/m, and 5000 N/m, respectively.
According to (24) and (25), the vector U and V can be written as follows.
Response points are chosen at coordinates 2, 4, and 6. Hammer impact is moving sequentially from points 2, 4 to 6. Then accelerances can be “measured”. The accelerances to be calculated are . Only the curves of , and are shown for the sake of brevity, as can be seen in Figure 9.
(a)
(b)
(c)
As shown in Figure 9, after adding three springmass substructures to the original structure, accelerances curves obtained by the proposed method coincide with the exactly values, which proves the accuracy of this method.
4. Conclusion
The work presented in this paper dealt with the problem of reanalyzing FRFs of a modified structure. A method is developed by employing SMW formula based on the FRFs related to the modification coordinates of original system. This method can solve the calculation problem of FRFs of modified structure in only one step. The accuracy and efficiency of this method are proved according to the simulated modal testing of cantilever beam model.
Two kinds of examples are used in verification. One is directly adding additional masses, stiffnesses, and damping to the cantilever beam. This kind of modification is mainly for the case that internal structural parameters matrices are changed. The other is adding springmass substructures to the cantilever beam. This condition is suitable for models whose original structure is immutable. Both numerical examples demonstrated good performance in the simulation verification. The results show that this approach is suitable for multipleelement change situation of structural modification and could improve the efficiency of structural reanalysis.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This work was supported by the Hubei Provincial Natural Science Foundation of China (2017CFB128), the Doctoral Scientific Research Foundation of Hubei University of Technology (BSQD14027), and Open fund of Hubei Key Laboratory of Modern Manufacturing Quality Engineering (KFJJ2020003).
References
 S.H. Tsai, H. Ouyang, and J.Y. Chang, “Inverse structural modifications of a geared rotorbearing system for frequency assignment using measured receptances,” Mechanical Systems and Signal Processing, vol. 110, pp. 59–72, 2018. View at: Publisher Site  Google Scholar
 D. J. Ewins, Modal Testing: Theory, Practice and Applications, Research Studies Press, England, UK, 2nd edition, 2000.
 J. He and Y. Li, “Relocation of antiresonances of a vibratory system by local structural changes,” International Journal of Analytical and Experimental Modal Analysis, vol. 10, no. 3, pp. 224–235, 1995. View at: Google Scholar
 K. T. Joseph, “Inverse eigenvalue problem in structural design,” AIAA Journal, vol. 30, no. 12, pp. 2890–2896, 1992. View at: Publisher Site  Google Scholar
 Y. M. Ram and S. G. Braun, “An inverse problem associated with modification of incomplete dynamic systems,” Journal of Applied Mechanics, vol. 58, no. 1, pp. 223–237, 1991. View at: Publisher Site  Google Scholar
 Y. G. Tsuei and E. K. L. Yee, “A method for modifying dynamic properties of undamped mechanical systems,” Journal of Dynamic Systems, Measurement, and Control, vol. 111, no. 3, pp. 403–408, 1989. View at: Publisher Site  Google Scholar
 J. F. Baldwin and S. G. Hutton, “Natural modes of modified structures,” AIAA Journal, vol. 23, no. 11, pp. 1937–1743, 1985. View at: Publisher Site  Google Scholar
 J.F. M. Barthelemy and R. T. Haftka, “Approximation concepts for optimum structural design—a review,” Structural Optimization, vol. 5, no. 3, pp. 129–144, 1993. View at: Publisher Site  Google Scholar
 B. P. Wang and W. D. Pilkey, “Eigenvalue reanalysis of locally modified structures using a generalized Rayleigh’s method,” AIAA Journal, vol. 24, no. 6, pp. 983–990, 1986. View at: Publisher Site  Google Scholar
 W. M. To and D. J. Ewins, “Structural modification analysis using Rayleigh quotient iteration,” International Journal of Mechanical Sciences, vol. 32, no. 3, pp. 169–179, 1990. View at: Publisher Site  Google Scholar
 S. H. Chen, X. W. Yang, and H. D. Lian, “Comparison of several eigenvalue reanalysis methods for modified structures,” Structural and Multidisciplinary Optimization, vol. 20, no. 4, pp. 253–259, 2000. View at: Publisher Site  Google Scholar
 Z. Li, B. S. Wu, and H. F. Liu, “Reanalysis method after structural boundary constraint modification,” Journal of Computational Mechanics, vol. 29, no. 5, pp. 699–703, 2012. View at: Google Scholar
 H. F. Liu, B. S. Wu, C. W. Lim, and Z. G. Li, “An approach for structural static reanalysis with unchanged number of degrees of freedom,” Structural and Multidisciplinary Optimization, vol. 45, no. 5, pp. 681–692, 2012. View at: Publisher Site  Google Scholar
 Q. Song, P. Chen, and S. Sun, “An exact reanalysis algorithm for local nontopological highrank structural modifications in finite element analysis,” Computers & Structures, vol. 143, pp. 60–72, 2014. View at: Publisher Site  Google Scholar
 G. Huang, H. Wang, and G. Li, “A reanalysis method for local modification and the application in largescale problems,” Structural and Multidisciplinary Optimization, vol. 49, no. 6, pp. 915–930, 2014. View at: Publisher Site  Google Scholar
 H. F. Liu, B. S. Wu, and Z. G. Li, “Method of updating the Cholesky factorization for structural reanalysis with added degrees of freedom,” Journal of Engineering Mechanics, vol. 140, no. 2, pp. 384–392, 2014. View at: Publisher Site  Google Scholar
 D. Materna and V. K. Kalpakides, “Nonlinear reanalysis for structural modifications based on residual increment approximations,” Computational Mechanics, vol. 57, no. 1, pp. 1–18, 2016. View at: Publisher Site  Google Scholar
 Y.S. Kim and H.C. Eun, “Seismic performance evaluation of footingtocircular RC pier connection reinforced by highmanganese steel bars (HMSBs),” Advances in Civil Engineering, vol. 2018, Article ID 3084078, 9 pages, 2018. View at: Publisher Site  Google Scholar
 S. P. Zheng, B. S. Wu, and Z. G. Li, “Free vibration reanalysis of structures with added degrees of freedom,” Computers & Structures, vol. 206, pp. 31–41, 2018. View at: Publisher Site  Google Scholar
 U. Krisch, “Combined approximations—a general reanalysis approach for structural optimization,” Structural and Multidisciplinary Optimization, vol. 20, no. 2, pp. 97–106, 2000. View at: Publisher Site  Google Scholar
 U. Kirsch, M. Bogomolni, and I. Sheinman, “Nonlinear dynamic reanalysis of structures by combined approximations,” Computer Methods in Applied Mechanics and Engineering, vol. 195, no. 33–36, pp. 4420–4432, 2006. View at: Publisher Site  Google Scholar
 W. Zuo, T. Xu, H. Zhang, and T. Xu, “Fast structural optimization with frequency constraints by genetic algorithm using adaptive eigenvalue reanalysis methods,” Structural and Multidisciplinary Optimization, vol. 43, no. 6, pp. 799–810, 2011. View at: Publisher Site  Google Scholar
 R. Sun, D. Liu, T. Xu, H. Zhang, and W. Zuo, “New adaptive technique of Kirsch method for structural reanalysis,” AIAA Journal, vol. 52, no. 3, pp. 486–495, 2014. View at: Publisher Site  Google Scholar
 W. Zuo, J. Bai, and J. Yu, “Sensitivity reanalysis of static displacement using Taylor series expansion and combined approximate method,” Structural and Multidisciplinary Optimization, vol. 53, no. 5, pp. 953–959, 2016. View at: Publisher Site  Google Scholar
 W. Zuo, K. Huang, J. Bai, and G. Guo, “Sensitivity reanalysis of vibration problem using combined approximations method,” Structural and Multidisciplinary Optimization, vol. 55, no. 4, pp. 1399–1405, 2017. View at: Publisher Site  Google Scholar
 J. Sherman and W. J. Morrison, “Adjustment of an inverse matrix corresponding to a change in one element of a given matrix,” The Annals of Mathematical Statistics, vol. 21, no. 1, pp. 124–127, 1950. View at: Publisher Site  Google Scholar
 M. A. Woodbury, Inverting Modified Matrices, Statistical Research Group, Princeton University, Princeton, NJ, USA, 1950.
 M. A. Akgün, J. H. Garcelon, and R. T. Haftka, “Fast exact linear and nonlinear structural reanalysis and the ShermanMorrisonWoodbury formulas,” International Journal for Numerical Methods in Engineering, vol. 50, no. 7, pp. 1587–1606, 2001. View at: Publisher Site  Google Scholar
 Z. Liu, W. Li, H. Ouyang, and D. Wang, “Eigen structure assignment in vibrating systems based on receptances,” Archive of Applied Mechanics, vol. 85, no. 6, pp. 713–724, 2015. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2020 Jun Ren and Qianghao Zhang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.