Next Article in Journal
LoRaWAN CubeSat with an Adaptive Data Rate: An Experimental Analysis of Path Loss Link Margin
Next Article in Special Issue
Space Robot On-Orbit Operation of Insertion and Extraction Impedance Control Based on Adaptive Neural Network
Previous Article in Journal
Modelling of Parachute Airborne Clusters Flight Dynamics and Parachute Interactions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis and Experimental Study of Excitation Force Transfer Path of Optical Satellite Refrigeration Module

1
State Key Laboratory of Robotics, Shenyang Institute of Automation, Chinese Academy of Sciences (CAS), Shenyang 110016, China
2
Institutes for Robotics and Intelligent Manufacturing, Chinese Academy of Sciences, Shenyang 110169, China
3
University of Chinese Academy of Sciences, Beijing 100049, China
4
School of Mechanical Engineering, Shenyang Ligong University, Shenyang 110159, China
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(1), 52; https://doi.org/10.3390/aerospace10010052
Submission received: 5 December 2022 / Revised: 21 December 2022 / Accepted: 29 December 2022 / Published: 4 January 2023
(This article belongs to the Special Issue Space Robotics and Mechatronics)

Abstract

:
When the problem of the reaction force under the receiver structure is complicated, the existing force-based transfer path analysis method is not suitable for its analysis. Therefore, an improved reaction force transfer path analysis method is proposed in this paper. In the transfer path analysis method based on the modal superposition method, the penalty function method is introduced, and the path contribution of the receiving end reaction force is evaluated by analyzing the displacement response of the connection between the transfer path and the receiving end. A simulation analysis of a satellite refrigeration module vibration isolation device designed by the research group identified the main factors affecting the vibration isolation efficiency of the vibration isolation ring, and provided guidance for subsequent further optimization. The device was tested and analyzed by building an experimental platform, and the obtained simulation data fit well with the test data. These conclusions further show that the proposed analysis method can more accurately and conveniently analyze the path contribution of the receiving end reaction force under specific working conditions. This method is suitable for analyzing the transfer path of the reaction force at the receiving end in the design stage.

1. Introduction

Vibration performance has always been the main aspect of evaluating the quality of machinery and equipment. Especially with the improvement of social productivity, the performance of various machinery and equipment has shown a trend of high speed, light weight and heavy load. Therefore, engineers need to continuously solve the problems of vibration control [1,2,3,4]. When the actual vibration mechanism is too complex, the transfer path analysis method (TPA) can be used to easily separate the source excitations from the structure/acoustic transmission characteristics, allowing engineers to eliminate the main path of vibration transmission, and then anticipate by making changes to either the source itself or the receiving structures that are connected to it [5]. Because of its convenient characteristics, TPA technology has been widely used in recent years. Based on the specific analysis of the transfer path in the vibration system, dynamic optimization technology can be guaranteed in the early stage of system design [6,7,8].
The traditional transfer path analysis method (TPA) is mainly based on the method of testing. It divides the mechanical system into three parts: the source, the transfer path, and the receiving end. The transfer function matrix between the source and the receiving end is obtained by isolating the vibration source, and then the source–path contribution information of the receiving end can be analyzed [9,10,11]. However, the isolation of the vibration source requires disassembly of the system, which requires a large workload, and the corresponding analysis accuracy is high [5,12].
In recent years, on the basis of traditional transfer path analysis methods, many scholars have done a lot of research to solve practical problems [1,5,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27]. In order to solve the problem of dismantling the system, a method of operation transfer path analysis is proposed [5,13,14,15,16,17,18]. De Sitter, G et al. used finite element simulation data to carry out operational transfer path analysis (OTPA) [5]. Zhang ZY et al. analyzed the operation transfer path of a heavy commercial truck cab installation system [13]. Deklerk and Ossipov applied OTPA based on singular value decomposition (SVD) to analyze tire noises [14]. Reising et al. applied OTPA based on SVD to effectively assess acoustic characteristics of the driven axle in automobiles [15]. Roozen et al. used a method for determining the number of effective transfer paths based on the singular value decomposition of the response matrix [16], and singular value decomposition (SVD) and principal component analysis (PCA) have been successfully used by Toome in his study of a source contribution prediction at low frequency [17] and by Strom in the OTPA of a high-speed train bogie [18]. Many scholars also analyze the transfer path from the perspective of energy flow [19,20,21,22,23,24]. Han X, et al. proposed an energy-based transfer path analysis method [19]. Kim and co-workers [20,21,22,23] studied vibration isolation and structure-borne noise transmission via multiple/multidimensional transfer paths. Ba-Leum K and Yoshida J et al. developed an improved transfer path analysis (MTPA) method to more accurately estimate the force of the main vibration source in a complex system with multiple vibration sources, basic excitations and multiple disturbances [24,25,26]. Lee identified the dominant vibration transfer path by experimentally measuring the multidimensional vibrational power flow through each isolator and obtained a 5 dB reduction in the sound pressure level in a vehicle by modifying the body structure accordingly [27].
However, most of the existing transfer path analysis methods use displacement, velocity, acceleration or sound pressure as the response for analysis, and rarely take reaction force as the research object. Maarten et al. introduced a method of reaction force transfer path analysis based on the traditional transfer path analysis method [5]. Wei Zhao et al. proposed the concept of the path transmission ratio, and derived the convenient formulae of reaction force based on the path transfer force analysis and the path disconnected method [1].
In the transfer path analysis methods of the reaction force introduced in the above documents, the reaction force is calculated by analyzing the response of the transfer path and the receiving end. In actual engineering applications, there are many cases where the force transmitted to the receiving end is used to measure the vibration performance of the system, and the receiving end is often much larger and more complex than the system. At this time, because the various system parameters of the receiving end are difficult to accurately obtain, the accuracy of the results obtained is difficult to guarantee by using the existing reaction force transfer path analysis method. At the same time, due to the complex structure of the receiving end, the amount of calculation in the analysis process is too large.
This article introduces in detail an improved reaction force transfer path analysis method, analyzes the existing reaction force transfer path analysis method, introduces the modal superposition method and the penalty function method, which will affect the support between the transfer path and the receiving end. The problem of solving the reaction force is transformed into a problem of solving the displacement between the two. Taking a satellite refrigeration module as an example, the reaction force between the transfer path and the receiving end is calculated, and the path contribution of the corresponding reaction force is analyzed. By setting up an experimental platform, the reaction force during the actual operation of the refrigeration module was measured and compared with the simulation data to verify the correctness of the proposed analysis method.

2. Theory of Force-Based Transfer Path Analysis

In engineering practice, many mechanical vibration systems have established their own transfer path models. For example, for the vehicle power transmission subsystem, the excitation of the running engine is transmitted to the vehicle body through a three-point suspension, thereby causing the vehicle body to vibrate. Normally, the transfer path model of a classic single-degree-of-freedom vibration system is shown in Figure 1.
In this system, the effects of stiffness, damping and mass are considered separately. Among them, m k , m 1 , m 2 , and m 3 are the equivalent masses of the vibration source, path1, path2, path3, and receiving end, respectively. k h g is the equivalent stiffness of the vibration source, and k l g is the equivalent stiffness of the receiver. k h 1 , k h 2 , and k h 3 are the equivalent stiffness between each path and the excitation source; k 11 , k 12 and k 13 are the equivalent stiffness between each path and the receiver. c h g and c 1 g are the equivalent damping of the vibration source and the receiver respectively. c h 1 , c h 2 and c h 3 are the equivalent damping of each path and vibration source respectively; c 11 , c 12 and c 13 are the equivalent damping of each path and receiver respectively. F ( t ) is the motivation of the system, which can be expressed as F ( t ) = F 0 e i ω t . x h , x 1 , x 2 , x 3 and x l are the displacements of the excitation end, each path, and the receiver respectively. F R 1 , F R 2 and F R 3 are the reaction forces between each path and the receiver respectively. Among these parameters, x 1 , x 2 , x 3 , F R 1 , F R 2 and F R 3 are all parameters to be solved.
The dynamic differential equation of the vibration system is
[ M ] { x ¨ } + [ C ] { x ˙ } + [ K ] { x } = { f ( t ) }
where [ M ] is the quality matrix of the system; [ K ] is the stiffness matrix; [ C ] is the damping matrix. { x ¨ } , { x ˙ } and { x } are the acceleration, velocity and displacement column vectors respectively. { f ( t ) } is the column vector of motivation.
In the system of Figure 1, mass, stiffness and damping are all known parameters, and the excitation force can be directly measured. The displacement, velocity and acceleration values of each position of the system can be obtained by solving the matrix.
The formula for calculating the reaction force between each path and the receiving end is
F Ri = k li [ x i ( t ) x l ( t ) ] + c li [ x ˙ i ( t ) x ˙ l ( t ) ]
Then the reaction force transfer function of the system can be expressed as
T Ri = | F Ri | | F 0 | ( i = 1 , 2 , 3 )
This analysis method has the following two defects:
(1) This method simplifies the relevant parameters of the system, and equates the mass, stiffness and damping parameters of each part to concentrated mass, concentrated stiffness and concentrated damping. Although this can improve the efficiency of calculation, it is difficult for these parameters to reflect the real boundary conditions of the system in the process of obtaining concentrated mass, concentrated stiffness and concentrated damping. Especially for complex systems, it is easy to get wrong results using this method.
(2) In actual engineering applications, there is a widespread situation in which a system A containing an excitation source is installed on another system B (system B is much larger than system A in terms of scale and complexity), as shown in Figure 2a. System A includes the excitation source and transfer path, and system B is the receiver. At this time, engineers often use the force of system A on system B as an index to measure the vibration performance of system A. When analyzing the transfer path of the reaction force between the two systems based on the above theory, engineers will spend a lot of energy and time in the process of obtaining the system parameters of system B, and it is difficult to ensure the accuracy of the system parameters of system B.
In engineering practice, in order to improve the efficiency of vibration analysis and ensure the accuracy of the vibration solution, the finite element analysis method is widely used in the early stage of mechanical design.
There are usually two methods to solve the dynamic response of structures by finite element: the modal superposition method and the stepwise integration method. The modal superposition method is widely used due to its high efficiency and high precision.
Combining with the solution of the vibration differential equation, it can be known that the displacement column vector can be expressed as the linear addition of each mode shape.
{ x } = α = 1 N q α { ψ α }
where { ψ α } is the α -th mode of the undamped system differential equation [ M ] { x ¨ } + [ K ] { x } = { f ( t ) } , q α is the mode participation factor of the α -th mode, and N is Modal order time.
Equation (4) is substituted into Equation (1) to get Equation (5)
[ M ] ( α = 1 N q ¨ α { ψ α } ) + [ C ] ( α = 1 N q ˙ α { ψ α } ) + [ K ] ( α = 1 N q α { ψ α } ) = { f ( t ) }
Both sides of Equation (5) are multiplied by { Ψ β } T to get Equation (6)
{ Ψ β } T [ M ] ( α = 1 N q ¨ α { ψ α } ) + { Ψ β } T [ C ] ( α = 1 N q ˙ α { ψ α } ) + { Ψ β } T [ K ] ( α = 1 N q α { ψ α } ) = { Ψ β } T { f ( t ) }
where { Ψ β } is the transpose matrix of the β -th order mode of the differential equation [ M ] { x ¨ } + [ K ] { x } = { f ( t ) } of the undamped system.
According to the orthogonality of the mode, Equation (7) is obtained.
{ { Ψ β } T [ M ] { ψ α } = { 0 , α β m α , α = β { Ψ β } T [ C ] { ψ α } = { 0 , α β c α , α = β { Ψ β } T [ K ] { ψ α } = { 0 , α β k α , α = β
where m α , c α and k α are the modal mass, modal damping and modal stiffness of the α -th mode of the system, respectively, and c α = am α + bk α ; a and b are proportionality constants.
Equation (6) can be simplified as
m α q ¨ α + c α q ˙ α + k α q α = { Ψ β } T { f ( t ) } = { Ψ α } T { f ( t ) }
The excitation load and response in Equation (8) are expressed in the complex form as Equation (9):
{ { f ( t ) } = { F } e j ω t q α = Q α e j ω t
Equation (9) is substituted into Equation (8) to get Equation (10)
( ω 2 m α + jc α + k α ) Q α e j ω t = { Ψ α } T { F } e j ω t
The expression of Q α is
Q α = { Ψ α } T { F } ω 2 m α + jc α + k α
Equation (9) and Equation (11) are substituted into Equation (4) to get Equation (12)
{ x } = α = 1 N q α { ψ α } = α = 1 N Q α e j ω t { ψ α } = α = 1 N { Ψ α } T { F } ω 2 m α + jc α + k α { ψ α } e j ω t
Thus, Equation (13) is obtained:
{ x } = α = 1 N { Ψ α } T { F } { ψ α } k α [ 1 ( ω 2 m α ) / k α + 2 j ζ α ( m α ω k α / m α ) / k α ] e j ω t
where ζ α is the damping ratio, ζ α = c α / c c = c α / ( 2 k α / m α ) , and c c is called the critical damping coefficient.
If the system only has an acting force F i at point i , Equation (14) is obtained.
{ F } = { 0 0   F i 0 0 } T
Equation (14) is substituted into Equation (13) to get Equation (15)
{ x } = α = 1 N Ψ i α F i { ψ α } k α [ 1 ( ω 2 m α ) / k α + 2 j ζ α ( m α ω k α / m α ) / k α ] e j ω t
where Ψ i α is the mode shape of the i -th degree of freedom of the α -th mode vector.
The response amplitude at point j of the system can be expressed as
{ x } = a = 1 N Ψ i a F i Ψ j a k a [ 1 ( ω 2 m a ) / k a + 2 j ζ a ( m a ω k a / m a ) / k α ]
where Ψ j α is the mode shape of the j -th degree of freedom of the α -th mode vector.
If the system is not affected by gravity and there is a branch reaction Force F Rj at point j of the system, then Equation (17) is obtained.
F R j = ( k α X j ) 2 + ( c α X ˙ j ) 2
If point j of the system is the constraint point, the penalty function method is adopted, that is, a spring with the stiffness of k s is connected at point j ( k s k α × 10 4 ), as shown in Figure 2b. Then Equation (18) is obtained.
F R j = [ ( k α + k s ) X j ] 2 + ( c α X ˙ j ) 2
Because of k s k α , Equation (19) is obtained.
F R j k s X j
Equation (16) and Equation (19) are combined to obtain Equation (20)
F R j = α = 1 N Ψ i α F i Ψ j α 1 ( ω 2 m α ) / k α + 2 j ζ α ( m α ω k α / m α ) / k α
and Equation (21) is obtained.
F R j F i = α = 1 N Ψ i α Ψ j α 1 ( ω 2 m α ) / k α + 2 j ζ α ( m α ω k α / m α ) / k α
Equation (21) defines the transfer function of the branch reaction force between the branch reaction force and the incentive force at the output point, which is expressed as
H Rji = F Rj F i

3. Simulation Analysis of Related Examples

For the vibration reduction device designed by the research group for an optical satellite refrigeration module, the design standard is that the reaction force in each direction of the output point is reduced by 90% compared with the excitation force in the corresponding direction. The vibration reduction device has been preliminarily designed.
The refrigeration module uses four refrigerators as excitation sources, which are labeled as I1, I2, I3, and I4. The centers of the four vibration isolation ring bases are the output points, which are marked as O1, O2, O3, and O4. The overall finite element model is shown in Figure 3.
The analysis method used in this paper is the force transfer path analysis method based on the modal superposition method. After the finite element modeling of the system is completed, the model is subjected to modal analysis. Generally speaking, the first few modes of the system have a large participation factor. Taking into account the speed of the finite element solution, this paper solves the first 10 modes of the system.

3.1. Measurement of Excitation Force

In engineering applications, the excitation force acquisition methods generally include the direct measurement method, inverse matrix method, dynamic stiffness method, etc. In this analysis, the excitation source is directly installed on the refrigeration module. Taking into account the convenience of measurement, the direct measurement method is adopted. At the same time, because the excitation source is not connected to the rest of the structure, the use of the direct measurement method will not affect the boundary conditions of the excitation source, which can also ensure that the obtained excitation force data has a higher accuracy.
In the process of measuring the excitation force, consider the four motors in this example as the excitation source; the excitation in the three directions of X, Y, and Z for each motor needs to be considered, and 12 excitation force data in total are measured. The frequency domain data obtained is shown in Figure 4.

3.2. Path Contribution Analysis

The analysis case in this article is a 12-input and 12-output model with a total of 144 transfer paths.
According to Equation (22), when the excitation force is only F i , that is { F } = { 0 0   F i 0 0 } T , Equation (23) is obtained.
F Rj = H Rji F i
If the excitation force is { F } = { F 1   F 2   F N } T , according to the linear superposition theorem, Equation (24) is obtained.
F R j = H R j 1 F 1 + H R j 2 F 2 + + H R j N F N = { H R j 1 H R j 2 H R j N } { F 1 F 2 F N } T
Then the reaction force of the system is expressed as:
{ F R } = { F R 1 F R 2 F R V } = [ H R 11 H R 12 H R I N H R 31 H R 22 H R 2 N H R N 1 H R N 2 H R N N ] { F 1 F 2 F N }
where [ H R ] is the transfer function matrix of the reaction force.
The extracted modal data is calculated to obtain the transfer function matrix of the system, and the extracted excitation force data is combined to calculate the reaction force data of the system. The data is shown in Figure 5.
Table 1 shows the maximum value of the reaction force and the damping efficiency data of each output point. It can be seen from the table that there are two places where the damping efficiency of the vibration isolation ring is lower than 90%, namely the Z-direction at the O1 point and the Z-direction at the O3 point. The path contributions of these two places are analyzed, and the results are shown in the Figure 6 and Figure 7.
As shown in Figure 6, the contribution of the Z-direction reaction force of O1 is analyzed, and it can be seen that the excitation points 1X, 1Y, 2X, and 2Y contribute more to the reaction force, but the contributions of the excitation points 2X and 2Y are both negative values, which can counteract the reaction force, so the paths I1:+X_O1:+Z and I1:+Y_O1:+Z are determined as the main reaction force transfer path.
As shown in Figure 7, the contribution of the Z-direction reaction force of the output point 3 is analyzed, and it can be seen that the contribution of the excitation point 4Z is the largest, and the contribution of the remaining paths is small. The decision path I4:+Z_O3:+Z is the main reaction force transfer path.
For the above three transfer paths, the transfer function of the full frequency range is shown in the following Figure 8.
The following two conclusions can be drawn by combining Figure 5 and Figure 8.
(1) For the two paths I1:+X_O1:+Z and I1:+Y_O1+Z where the contribution of output point O1 is larger, the excitation force amplitude is small. The reason for the larger reaction force at this point is that the path transfer characteristics of the two paths are not good. The subsequent optimization process is mainly aimed at optimizing the path transfer characteristics of the two transfer paths.
(2) For the transfer path I4:+Z_-O3:+Z, the transfer function at 64.7 Hz is small, the excitation force level is high, and the reason for the large reaction force in the 3Z direction of the output point is the large excitation force. The vibration reduction optimization is carried out around the excitation source.

4. Test Analysis of Reaction Force

The refrigeration module is applied to the satellite, and the analyzed working condition is the operating state of the satellite in orbit. In the process of building the experimental platform, a four-point suspension was used to simulate the non-gravity environment in space, and the suspension was suspended by a tension spring, which has a long spring length, which can ensure that the suspension hardly affects the accuracy of the experimental data. Figure 9 is an analysis test diagram of the reaction force transfer path of the refrigeration module.
In the simulation analysis, the vibration damping efficiency of the vibration isolation ring does not meet the design requirements, and the two reaction forces are in the Z direction at O1 point and the Z direction at O3 point. In the experiment, these two reaction forces are mainly measured and compared with the data obtained from the simulation analysis to verify the correctness of the simulation analysis.
Figure 10 and Figure 11 are comparison pictures of these two reaction forces.
It can be seen from the above two comparison figures that at the corresponding frequency, the fit of the experimental data and the simulated data can reach more than 90%, and the trend of the two data is basically fitted.
Compared with the simulated data, the experimental data has more peaks. This is because the data collected during the experiment contains noise, which is obvious in the low frequency range. At the same time, the experimental data is smaller. The reason is that the suspension mode is used to simulate the non-gravity environment in space under the test conditions. The suspension device bears a certain vibration input, which is also in line with expectations.

5. Results

Based on the traditional transfer path analysis, this paper adopts the modal superposition method and the penalty function method to establish an improved reaction force transfer path analysis method. This method is suitable for the vibration reduction analysis of the mechanical system in the design stage because there is no need to consider the structural parameters of the next system, and this method has a greater advantage in solving speed.
In this paper, by analyzing the transfer path of the force for the vibration damping design of a satellite refrigeration module by the research group, the output point positions with unsatisfactory vibration damping efficiency are identified, and the path contribution of these positions is analyzed. Through analysis, we can identify the main factors that affect the effect of vibration reduction, and clarify the direction for subsequent further optimization.
By building an experimental platform, the refrigeration module was tested and analyzed for the reaction force transfer path. Comparing the experimental data with the simulation data, it can be seen that the two types of data fit well, which proves the correctness of the established method for analyzing the reaction force transfer path.
Because this method is based on simulation, its solution accuracy depends on the accuracy of the finite element model. At the same time, idealizing the receiving end of the system during the analysis process will also bring certain errors to the analysis results. This is the next direction that needs improvement.

Author Contributions

Conceptualization, H.L. and Z.L.; methodology, H.L.; software, Z.L., C.F. and F.Z.; validation, H.L., Z.L. and C.F.; formal analysis, Z.L.; investigation, F.Z.; resources, H.L.; data curation, H.L.; writing—original draft preparation, Z.L.; writing—review and editing, Z.L.; visualization, C.F.; supervision, H.L.; project administration, H.L.; funding acquisition, H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was financially supported by the State Key Laboratory of Robotics (2022-Z01, 2022-Z09); National Natural Science Foundation of China (No. 51975567); Liao Ning Revitalization Talents Program (XLYC1907152); Natural Science Foundation of Liaoning Province (2019-MS-029, 2021-MS-029); Chinese Academy of Science Youth Innovation Promotion Association (No. 2018237) and Development Fund of Space Automation Technology Laboratory, SIA, CAS.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, W.; Zhang, Y.-M.; Feng, W.-Z. Path Contribution Analysis of Vibration Transfer Path Systems. Math. Probl. Eng. 2019, 2019, e1458149. [Google Scholar] [CrossRef]
  2. Vizzini, S.; Olsson, M.; Scattina, A. Component Mode Synthesis Methods for a Body-in-White Noise and Vibration Analysis. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2017, 231, 279–288. [Google Scholar] [CrossRef]
  3. Li, D.; Huang, Y. The Discomfort Model of the Micro Commercial Vehicles Interior Noise Based on the Sound Quality Analyses. Appl. Acoust. 2018, 132, 223–231. [Google Scholar] [CrossRef]
  4. Jafari, S.; Ioannou, P.; Rudd, L. Adaptive Feedback Suppression of Unknown Periodic Components of Acoustic Noises with Time-Varying Characteristics. J. Vib. Control 2017, 23, 526–538. [Google Scholar] [CrossRef]
  5. van der Seijs, M.V.; de Klerk, D.; Rixen, D.J. General Framework for Transfer Path Analysis: History, Theory and Classification of Techniques. Mech. Syst. Signal Process. 2016, 68–69, 217–244. [Google Scholar] [CrossRef]
  6. Cheng, W.; Zhang, Z.; Lee, S.; He, Z. Investigations of Denoising Source Separation Technique and Its Application to Source Separation and Identification of Mechanical Vibration Signals. J. Vib. Control 2014, 20, 2100–2117. [Google Scholar] [CrossRef]
  7. Cheng, W.; Lu, Y.; Zhang, Z. Tikhonov Regularization-Based Operational Transfer Path Analysis. Mech. Syst. Signal Process. 2016, 75, 494–514. [Google Scholar] [CrossRef]
  8. Cheng, W.; He, Z.; Zhang, Z. A Comprehensive Study of Vibration Signals for a Thin Shell Structure Using Enhanced Independent Component Analysis and Experimental Validation. J. Vib. Acoust. 2014, 136. [Google Scholar] [CrossRef]
  9. Plunt, J. Finding and Fixing Vehicle NVH Problems with Transfer Path Analysis. Sound Vib. 2005, 39, 12–17. [Google Scholar]
  10. Shin, T.; Kim, Y.S.; An, K.; Lee, S.K. Transfer Path Analysis of Rumbling Noise in a Passenger Car Based on In-Situ Blocked Force Measurement. Appl. Acoust. 2019, 149, 1–14. [Google Scholar] [CrossRef]
  11. Elliott, A.S.; Moorhouse, A.T.; Huntley, T.; Tate, S. In-Situ Source Path Contribution Analysis of Structure Borne Road Noise. J. Sound Vib. 2013, 332, 6276–6295. [Google Scholar] [CrossRef]
  12. De Sitter, G.; Devriendt, C.; Guillaume, P.; Pruyt, E. Operational Transfer Path Analysis. Mech. Syst. Signal Process. 2010, 24, 416–431. [Google Scholar] [CrossRef]
  13. Zhang, Z.; Pan, D.; Wu, W.; Huang, C. Vibration Source Identification of a Heavy Commercial Vehicle Cab Based on Operational Transfer Path Analysis. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2020, 234, 669–680. [Google Scholar] [CrossRef]
  14. de Klerk, D.; Ossipov, A. Operational Transfer Path Analysis: Theory, Guidelines and Tire Noise Application. Mech. Syst. Signal Process. 2010, 24, 1950–1962. [Google Scholar] [CrossRef]
  15. Reising, M.; Küçükay, F.; Eulert, S.; Decker, H.; Luehrs, G.-F. Hybrid Acoustic Synthesis–Assessment of Rear Final Drives in an Early Stage of Development, on the Basis of Transfer Path Analyses and Test Bench Measurements. Appl. Acoust. 2014, 80, 45–56. [Google Scholar] [CrossRef]
  16. Roozen, N.B.; Leclère, Q. On the Use of Artificial Excitation in Operational Transfer Path Analysis. Appl. Acoust. 2013, 74, 1167–1174. [Google Scholar] [CrossRef] [Green Version]
  17. Toome, M. Operational Transfer Path Analysis, a Study of Source Contribution Predictions at Low Frequency. Master’s Thesis, Chalmers University of Technology, Gothenburg, Sweden, 2012. [Google Scholar]
  18. Strom, R. Operational Transfer Path Analysis of Components of a High-Speed Train Bogie. Master’s Thesis, Chalmers University of Technology, Gothenburg, Sweden, 2014. [Google Scholar]
  19. Han, X.; Guo, Y.-J.; Zhao, Y.-E.; Lin, Z.-Q. The Application of Power-Based Transfer Path Analysis to Passenger Car Structure-Borne Noise. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2008, 222, 2011–2023. [Google Scholar] [CrossRef]
  20. Kim, S.; Singh, R. Vibration transmission through an isolator modelled by continuous system theory. J. Sound Vib. 2001, 248, 925–953. [Google Scholar] [CrossRef] [Green Version]
  21. Singh, R.; Kim, S. Examination of Multi-Dimensional Vibration Isolation Measures and Their Correlation to Sound Radiation over a Broad Frequency Range. J. Sound Vib. 2003, 262, 419–455. [Google Scholar] [CrossRef]
  22. Kim, S.; Inoue, A.; Singh, R. Experimental Study of Structure-Borne Noise Transfer Paths over the Mid-Frequency Regime; SAE Technical Technical Papers; SAE International: Warrendale, PA, USA, 2005. [Google Scholar] [CrossRef]
  23. Singh, R.; Kim, S. Prediction of structure-borne noise transmission through multiple/multi-dimensional transfer paths. In Proceedings of the International Congress on Noise and Vibration: Emerging Methods, Saint Raphael, France, 18–21 April 2005. [Google Scholar]
  24. Kim, B.-L.; Jung, J.-Y.; Oh, I.-K. Modified Transfer Path Analysis Considering Transmissibility Functions for Accurate Estimation of Vibration Source. J. Sound Vib. 2017, 398, 70–83. [Google Scholar] [CrossRef]
  25. Gajdatsy, P.; Janssens, K.; Desmet, W.; Van der Auweraer, H. Application of the Transmissibility Concept in Transfer Path Analysis. Mech. Syst. Signal Process. 2010, 24, 1963–1976. [Google Scholar] [CrossRef]
  26. Lennström, D.; Johnsson, R.; Ågren, A.; Nykänen, A. The Influence of the Acoustic Transfer Functions on the Estimated Interior Noise from an Electric Rear Axle Drive. SAE Int. J. Passeng. Cars-Mech. Syst. 2014, 7, 413–422. [Google Scholar] [CrossRef]
  27. Lee, S.K. Identification of a Vibration Transmission Path in a Vehicle by Measuring Vibrational Power Flow. Proc. Inst. Mech. Eng. Part D J. Automob. Eng. 2004, 218, 167–175. [Google Scholar] [CrossRef]
Figure 1. Transfer path model of single-degree-of-freedom vibration system.
Figure 1. Transfer path model of single-degree-of-freedom vibration system.
Aerospace 10 00052 g001
Figure 2. Schematic diagram of project example model: (a) Use rigid connection; (b) Use high-stiffness spring instead of connection.
Figure 2. Schematic diagram of project example model: (a) Use rigid connection; (b) Use high-stiffness spring instead of connection.
Aerospace 10 00052 g002
Figure 3. The overall finite element model of the satellite refrigeration module: (a) Positive; (b) Negative.
Figure 3. The overall finite element model of the satellite refrigeration module: (a) Positive; (b) Negative.
Aerospace 10 00052 g003
Figure 4. Frequency domain excitation force data of each motor within 5–200 Hz: (a) I1; (b) I2; (c) I3; (d) I4.
Figure 4. Frequency domain excitation force data of each motor within 5–200 Hz: (a) I1; (b) I2; (c) I3; (d) I4.
Aerospace 10 00052 g004
Figure 5. Reaction force data of each response point: (a) O1; (b) O2; (c) O3; (d) O4.
Figure 5. Reaction force data of each response point: (a) O1; (b) O2; (c) O3; (d) O4.
Aerospace 10 00052 g005
Figure 6. Analysis diagram of Z-direction path contribution of output point O1.
Figure 6. Analysis diagram of Z-direction path contribution of output point O1.
Aerospace 10 00052 g006
Figure 7. Analysis diagram of Z-direction path contribution of output point O3.
Figure 7. Analysis diagram of Z-direction path contribution of output point O3.
Aerospace 10 00052 g007
Figure 8. Analysis diagram of Z-direction path contribution of output point O3.
Figure 8. Analysis diagram of Z-direction path contribution of output point O3.
Aerospace 10 00052 g008
Figure 9. Transfer path analysis experiment of reaction force.
Figure 9. Transfer path analysis experiment of reaction force.
Aerospace 10 00052 g009
Figure 10. Test and simulation data comparison of Z-direction reaction force at O1 point.
Figure 10. Test and simulation data comparison of Z-direction reaction force at O1 point.
Aerospace 10 00052 g010
Figure 11. Test and simulation data comparison of Z-direction reaction force at O3 point.
Figure 11. Test and simulation data comparison of Z-direction reaction force at O3 point.
Aerospace 10 00052 g011
Table 1. Maximum reaction force of each output point and damping efficiency.
Table 1. Maximum reaction force of each output point and damping efficiency.
Output
Point
X-DirectionY-DirectionZ-Direction
Frequency
(Hz)
Amplitude
(N)
Damping Rate (%)Frequency
(Hz)
Amplitude
(N)
Damping Rate (%)Frequency
(Hz)
Amplitude
(N)
Damping Rate (%)
O152.90.18194.99.20.13794.09.20.49883.9
O252.90.12994.452.90.14793.052.90.27891.0
O352.90.24691.89.20.10895.364.70.32189.6
O452.90.18793.79.20.07097.052.90.22392.8
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Luo, H.; Liu, Z.; Fan, C.; Zhao, F. Analysis and Experimental Study of Excitation Force Transfer Path of Optical Satellite Refrigeration Module. Aerospace 2023, 10, 52. https://doi.org/10.3390/aerospace10010052

AMA Style

Luo H, Liu Z, Fan C, Zhao F. Analysis and Experimental Study of Excitation Force Transfer Path of Optical Satellite Refrigeration Module. Aerospace. 2023; 10(1):52. https://doi.org/10.3390/aerospace10010052

Chicago/Turabian Style

Luo, Haitao, Ziyang Liu, Chaohui Fan, and Fengqun Zhao. 2023. "Analysis and Experimental Study of Excitation Force Transfer Path of Optical Satellite Refrigeration Module" Aerospace 10, no. 1: 52. https://doi.org/10.3390/aerospace10010052

APA Style

Luo, H., Liu, Z., Fan, C., & Zhao, F. (2023). Analysis and Experimental Study of Excitation Force Transfer Path of Optical Satellite Refrigeration Module. Aerospace, 10(1), 52. https://doi.org/10.3390/aerospace10010052

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