Method and equipment for determining vibration contribution degree of internal vibration source of structure

By monitoring the vibration signals of vibration sources and precision instruments, identifying the main frequencies and removing noise interference, and establishing a vibration transfer matrix, the problem of difficulty in judging the contribution of vibration sources in industrial buildings is solved, and the normal working requirements of precision instruments are met.

CN120764280APending Publication Date: 2025-10-10HARBIN INSTITUTE OF TECHNOLOGY (SHENZHEN) (INSTITUTE OF SCIENCE AND TECHNOLOGY INNOVATION HARBIN INSTITUTE OF TECHNOLOGY SHENZHEN)
View PDF 0 Cites 0 Cited by

Patent Information

Application Number
CN202510937459.2
Authority / Receiving Office
CN · China
Patent Type
Applications(China)
Current Assignee / Owner
Filing Date
2025-07-08
Publication Date
2025-10-10

AI Technical Summary

Technical Problem

Existing technologies make it difficult to accurately determine the vibration contribution of various vibration sources in industrial buildings to precision instruments, resulting in the inability to effectively control the vibration response of precision instruments, affecting their normal operation.

Method used

By monitoring the vibration input and output signals of the vibration source and precision instruments, the fast Fourier transform is used to identify the main frequency components and remove noise interference, establish a vibration transfer matrix, decompose the vibration contribution and calculate the contribution of each vibration source.

Benefits of technology

The accurate calculation of the vibration contribution of the vibration source at the output end of the precision instrument is achieved, ensuring the normal operation of the precision instrument in a high vibration environment.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN120764280A_ABST
    Figure CN120764280A_ABST
Patent Text Reader

Abstract

The invention relates to the technical field of civil engineering, in particular to a method and equipment for determining the vibration contribution degree of a vibration source in a structure, and the method comprises the steps: firstly obtaining an input load of a vibration source arranged under a vibration source arrangement point location in a target structure and a fast Fourier transform result of the input load; the method comprises the following steps: determining the main load frequency and spectrum of a vibration source, then obtaining the output time-history response of a precise instrument under a precise instrument arrangement point location in a target structure and the fast Fourier transform result of the output time-history response, determining the main vibration frequency and spectrum of the output end of the precise instrument, and finally determining the main vibration frequency and spectrum of the output end of the precise instrument according to the main load frequency and spectrum of the vibration source. Analyzing based on finite element simulation, establishing a dynamic flexibility function of the input end of the vibration source and the output end of the precise instrument, establishing a transfer matrix according to the dynamic flexibility function, and decomposing the frequency spectrum of the output end of the precise instrument according to the transfer matrix to obtain the frequency spectrum of the output end of the precise instrument. And determining the vibration contribution degree of the internal vibration source of the target structure.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] The present invention relates to the technical field of civil engineering, and in particular to a method and device for determining the vibration contribution of a vibration source inside a structure. Background Art

[0002] Industrial relocation is a new high-quality development model emerging to address the accelerated industrial transformation and upgrading and increasingly limited land resources. It also involves a stacked factory model that integrates upstream and downstream industrial chains, rather than the single-sector, single-story factory building. In this stacked factory model, upstream and downstream industrial chains are relocated together within the building to carry out production and processing activities. This results in kN-level excitation forces from downstream equipment in industries such as metal smelting, automotive manufacturing, and textile manufacturing being transmitted to upstream equipment in industries such as ultra-high-definition video display, precision instruments, and high-end intelligent equipment manufacturing, where vibration tolerances reach µm levels. Industrial relocation has the following characteristics: First, the downstream industry chain's vibration sources are diverse and numerous, resulting in high vibration input. Second, compared to single-story factories, upstream and downstream equipment are arranged in layers, allowing source vibration to be transmitted vertically along the structure. Third, precision instruments have higher vibration tolerance levels (VCs). Failure to clearly define the contribution of each source's vibration response to the precision instrument's vibration can cause vibration at the instrument to exceed the tolerance level, thereby impacting production and processing activities.

[0003] To determine the extent of impact of various vibration sources on precision instruments when their vibration tolerances are exceeded, existing methods involve setting different start-up and shutdown conditions for the power equipment's vibration sources and analyzing the vibration response of the precision instruments under these conditions to determine the extent of impact. While this method can identify the cause of abnormal vibration in precision instruments, it is time-consuming and labor-intensive, and it cannot quantitatively analyze the contribution of each vibration source to the vibration of the precision instrument. Summary of the Invention

[0004] In view of this, an object of the present invention is to provide a method and apparatus for determining the vibration contribution of internal vibration sources of a structure, so as to overcome the problems existing in the current prior art.

[0005] To achieve the above objectives, the present invention adopts the following technical solutions: In one aspect, the present application provides a method for determining the vibration contribution of a vibration source inside a structure, comprising: Step S1, obtaining the input load of the vibration source arranged at the vibration source arrangement point inside the target structure and the fast Fourier transform result of the input load, and determining the main load frequency and spectrum of the vibration source; Step S2, obtaining the output time history response of the precision instrument at the precision instrument arrangement point inside the target structure and the fast Fourier transform result of the output time history response, and determining the main vibration frequency and spectrum of the output end of the precision instrument; Step S3, performing analysis based on finite element simulation based on the main load frequency and spectrum of the vibration source and the main vibration frequency and spectrum of the output end of the precision instrument, establishing a dynamic flexibility function of the input end of the vibration source and the output end of the precision instrument, and establishing a transfer matrix based on the dynamic flexibility function; Step S4: Decomposing the frequency spectrum of the output end of the precision instrument according to the transfer matrix to determine the vibration contribution of the internal vibration source of the target structure.

[0006] Furthermore, in the above method, step S1 includes: Obtaining a load input time history of the vibration source through a vibration pickup installed on a vibration source arranged below a vibration source arrangement point inside the target structure; Decomposing the load input time history of the vibration source to determine the simple harmonic excitation signal and noise signal of the vibration source; Denoising the load input time history of the decomposed vibration source; Generating a load input time history vector of the internal vibration source of the target structure based on a denoising result of the load input time history of the vibration source; Performing Fourier transform on the load input time history vector to obtain an input load spectrum vector containing a main load frequency of an internal vibration source of the target structure; The main load frequency and spectrum of the vibration source are determined according to the input load spectrum vector.

[0007] Furthermore, in the above method, step S2 includes: Acquiring a time history response of an output end of the precision instrument through a vibration pickup installed at a precision instrument arrangement point inside the target structure, and generating a response vector of the output end of the precision instrument; Performing the denoising process on the response vector of the output end of the precision instrument, and performing Fourier transform on the denoising process result to obtain the frequency spectrum vector of the output end of the precision instrument; According to the frequency spectrum vector of the output end of the precision instrument, the main vibration frequency and frequency spectrum of the output end of the precision instrument are determined.

[0008] Furthermore, in the above method, step S3 includes: According to the main load frequency and spectrum of the vibration source, as well as the main vibration frequency and spectrum of the output end of the vibration source, analysis is performed based on finite element simulation to determine the response vectors of the output end of the precision instrument inside the target structure under the excitation of the vibration source load inside the target structure, one by one, and each response vector is recorded as a matrix element of a row of the transfer matrix; All the generated response vectors are arranged in the order of the internal vibration sources of the target structure to generate the transfer matrix.

[0009] Furthermore, in the above method, step S4 includes: Decomposing the frequency spectrum vector at the output end of the precision instrument according to the linear system characteristics and the transfer matrix, and determining an expression of any element in the frequency spectrum vector at the output end of the precision instrument; Determining a vibration contribution calculation formula of any vibration source inside the target structure based on an expression of any element in the frequency spectrum vector at the output end of the precision instrument; The vibration contribution of the internal vibration source of the target structure is determined according to the vibration contribution calculation formula.

[0010] Furthermore, the above method, the denoising process, includes: Determining a signal to be denoised, and performing Fourier transform on the signal to be denoised; Sort the corresponding main frequencies from large to small according to the transformation amplitude; Performing inverse Fourier transform on the simple harmonic signals corresponding to the main frequencies according to the order of the main frequencies to obtain a reconstructed signal; Subtracting the reconstructed signal from the vibration measured signal of the signal to be denoised to obtain a residual signal; If the autocorrelation function of the residual signal has a non-zero point equal to 0, the residual signal is filtered out.

[0011] On the other hand, the present application provides a device for determining the vibration contribution of a vibration source inside a structure, including a processor and a memory, wherein the processor is connected to the memory: The processor is configured to call and execute the program stored in the memory; The memory is used to store the program, and the program is used to at least execute any of the above methods for determining the vibration contribution of internal vibration sources of a structure.

[0012] The beneficial effects of the present invention are: This application monitors the vibration input load input time history of the internal vibration source of the target structure and the vibration output time history response of the precision instrument end, and identifies the main frequency components of the load through FFT transformation and removes the noise interference components; constructs the vibration transfer relationship between the vibration load input end and the vibration response output end by establishing a vibration transfer matrix; decomposes the vibration output spectrum into a linear superposition of the contributions of each vibration load based on the vibration transfer matrix and the vibration load input spectrum, and determines the main vibration source of the vibration output end based on the contribution of each vibration source component to the output end spectrum, thereby realizing the accurate calculation of the vibration contribution of the vibration source at the output end of the precision instrument and determining the main vibration source, thereby ensuring the normal operation requirements of the precision instrument. BRIEF DESCRIPTION OF THE DRAWINGS

[0013] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0014] Figure 1 This is a flow chart provided by an embodiment of a method for determining the vibration contribution of a vibration source inside a structure according to the present invention; Figure 2 This is a denoising process flow chart provided by an embodiment of a method for determining the vibration contribution of a vibration source inside a structure according to the present invention; Figure 3 It is a structural schematic diagram provided by an embodiment of a device for determining the vibration contribution of a vibration source inside a structure according to the present invention. DETAILED DESCRIPTION

[0015] To make the objectives, technical solutions, and advantages of the present invention more apparent, the technical solutions of the present invention will be described in detail below. Obviously, the embodiments described are only some of the embodiments of the present invention, not all of them. Based on the embodiments of the present invention, all other implementations obtained by those of ordinary skill in the art without inventive effort are within the scope of protection of the present invention.

[0016] Figure 1 This is a flow chart of an embodiment of a method for determining the vibration contribution of a structure's internal vibration source. Figure 1 , this embodiment may include the following steps: S1. Obtain the input load of the vibration source arranged at the vibration source arrangement point inside the target structure and the fast Fourier transform result of the input load to determine the main load frequency and spectrum of the vibration source.

[0017] S2. Obtain the output time history response of the precision instrument at the precision instrument arrangement point inside the target structure and the fast Fourier transform result of the output time history response, and determine the main vibration frequency and spectrum of the output end of the precision instrument.

[0018] S3. Based on the main load frequency and spectrum of the vibration source and the main vibration frequency and spectrum of the output end of the precision instrument, an analysis is performed based on finite element simulation to establish the dynamic flexibility function of the input end of the vibration source and the output end of the precision instrument, and a transfer matrix is ​​established based on the dynamic flexibility function.

[0019] S4. Decompose the frequency spectrum at the output of the precision instrument according to the transfer matrix to determine the vibration contribution of the internal vibration source of the target structure.

[0020] It can be understood that the application identifies the main load frequency components and removes noise interference components by monitoring the target structure internal vibration source vibration input load input time history and the precision instrument end vibration output time history response and through FFT transformation; the vibration transmission relationship between the vibration load input end and the vibration response output end is constructed by establishing the vibration transmission matrix; the vibration output frequency spectrum is decomposed into the linear superposition of each vibration load contribution according to the vibration transmission matrix and the vibration load input frequency spectrum, and the main vibration source of the vibration output end is determined according to the contribution degree of each vibration source component to the output end frequency spectrum, so that the vibration source vibration contribution degree of the precision instrument output end is accurately calculated and the main vibration source is determined, and the normal working requirement of the precision instrument is ensured.

[0021] It should be noted that the target structure can be an industrial building.

[0022] Preferably, step S1 comprises: obtaining the load input time history of the vibration source through the vibration pickup instrument installed on the vibration source arranged at the point position inside the target structure; decomposing the load input time history of the vibration source to determine the simple harmonic excitation signal and the noise signal of the vibration source; de-noising the load input time history of the decomposed vibration source; generating the load input time history vector of the vibration source inside the target structure according to the de-noising result of the load input time history of the vibration source; performing Fourier transform on the load input time history vector to obtain the input load frequency spectrum vector of the vibration source inside the target structure containing the main load frequency; determining the main load frequency and the frequency spectrum of the vibration source according to the input load frequency spectrum vector.

[0023] It can be understood that the internal power equipment of the industrial building is generally large in size and has fixed load action form, so it is relatively easy to identify the load, and the vertical vibration acceleration of the equipment surface can be monitored through the vibration pickup instrument installed thereon and the acceleration signal can be converted into the load signal through Newton's second law.

[0024] According to the vibration source load input time history of each internal power equipment measured at different arrangement points, the vector is written as: , wherein represents the load input time history of the i-th vibration source, and n represents the number of vibration sources.

[0025] For rotating equipment such as fans, motors, pumps and the like, the vibration signal can be decomposed into a series of simple harmonic excitation signals and noise signals: , wherein is the noise signal, ​Indicates the number of simple harmonic components.

[0026] In order to reduce the impact of noise on monitoring accuracy, the randomness of the noise signal will cause its autocorrelation function to be almost 0 at non-zero points. Perform noise identification and denoising to obtain the denoised input time history vectors of each vibration source load:

[0027] right Perform FFT transformation on the vector to obtain the input load spectrum vector of the internal vibration source containing the main load frequency: ,in Indicates the The spectrum of the vibration source load.

[0028] Preferably, step S2 includes: The time history response of the precision instrument output end is obtained by using a vibration pickup installed at the precision instrument arrangement point inside the target structure, and a response vector of the precision instrument output end is generated; De-noising the response vector at the output of the precision instrument, and performing Fourier transform on the de-noising result to obtain the frequency spectrum vector at the output of the precision instrument; According to the frequency spectrum vector of the output end of the precision instrument, the main vibration frequency and frequency spectrum of the output end of the precision instrument are determined.

[0029] It can be understood that by arranging vibration pickups at each precision instrument arrangement point to monitor the output end response, the precision instrument output end response vector is obtained: ,in Indicates the The output time response, The number of output points.

[0030] Similarly, after denoising and FFT transformation, the spectrum vectors of each output end can be obtained: ,in Indicates the The output spectrum.

[0031] Preferably, step S3 includes: Based on the main load frequency and spectrum of the vibration source, as well as the main vibration frequency and spectrum of the vibration source output end, finite element simulation is used to analyze and determine the response vectors of the output end of the precision instrument inside the target structure under the excitation of the vibration source load inside the target structure. Each response vector is recorded as a matrix element in a row of the transfer matrix. All the generated response vectors are arranged according to the order of the internal vibration sources of the target structure to generate a transfer matrix.

[0032] It can be understood that by normalizing the output-end response spectrum by dividing it by the input-end load spectrum, the precision instrument-end response spectrum (dynamic compliance function) under unit force input of the power equipment is obtained:

[0033] The first vibration source excitation is obtained by actual measurement The output terminal response is transformed by FFT and recorded as the first row of the transfer matrix:

[0034] Get the first Under the excitation of a vibration source The output terminal response is recorded as the transfer matrix after FFT transformation. OK:

[0035] in Indicates the Under the excitation of the vibration source unit, the The response spectrum of each output terminal is shown in Figure 1.

[0036] Establish an industrial building finite element model, extract the dynamic flexibility curves of each vibration source load input point with an amplitude of 1 and each output point in the sweep frequency range including the main frequency, and write , which is the transfer matrix.

[0037] Preferably, step S4 includes: Decomposing the spectrum vector of the precision instrument output end according to the linear system characteristics and the transfer matrix, and determining the expression of any element in the spectrum vector of the precision instrument output end; According to the expression of any element in the frequency spectrum vector at the output end of the precision instrument, the vibration contribution calculation formula of any vibration source inside the target structure is determined; According to the vibration contribution calculation formula, the vibration contribution of the internal vibration source of the target structure is determined.

[0038] It is understandable that the structure of the industrial building is in the linear elastic stage during its service life, and the vibration input of the vibration source meets the requirements of superposition and homogeneity for the output response of the precision instrument. According to the linear system characteristics, it can be obtained:

[0039] According to this relationship, we can get Any element in a vector Both The result of linear superposition of multiple vibration sources.

[0040] Will Decompose into The expression for the superposition of the individual vibration sources is:

[0041] Rule No. Source pair The spectrum contribution of each output endpoint bit is: .

[0042] Preferably, the denoising process includes: Determine the signal to be denoised, and perform Fourier transform on the signal to be denoised; Sort the corresponding main frequencies from large to small according to the transformation amplitude; According to the order of the main frequencies, the simple harmonic signals corresponding to the main frequencies are subjected to inverse Fourier transform to obtain the reconstructed signals; The reconstructed signal is subtracted from the vibration measured signal of the signal to be denoised to obtain a residual signal; If the autocorrelation function of the residual signal is not 0 at any point, the residual signal is filtered out.

[0043] It is understandable that if Figure 2 As shown in the figure, the original signal is first subjected to FFT transformation, and the corresponding main frequencies are sorted from large to small according to the transformation amplitude. The simple harmonic signal with a larger amplitude is selected and subjected to inverse Fourier transform to obtain the reconstructed signal. The reconstructed signal is subtracted from the measured vibration signal to obtain the residual signal. If the autocorrelation function of the residual signal is 0 or almost 0 at a non-zero point, then the residual signal is the noise component that needs to be filtered out.

[0044] The present invention also provides a device for determining the vibration contribution of a vibration source inside a structure, which is used to implement the above method embodiment. Figure 3 This is a schematic diagram of a structure provided by an embodiment of a device for determining the vibration contribution of a structure internal vibration source according to the present invention. Figure 3 As shown, the apparatus for determining the vibration contribution of internal vibration sources in a structure in this embodiment includes a processor 21 and a memory 22, wherein the processor 21 is connected to the memory 22. The processor 21 is configured to call and execute a program stored in the memory 22; the memory 22 is configured to store the program, which is configured to at least execute the method for determining the vibration contribution of internal vibration sources in the structure in the above embodiment.

[0045] The specific implementation scheme of the device for determining the vibration contribution of a structure's internal vibration source provided in the embodiments of the present application can refer to the implementation scheme of the method for determining the vibration contribution of a structure's internal vibration source in any of the above embodiments, and will not be repeated here.

[0046] It can be understood that the same or similar parts of the above embodiments can be referenced to each other, and the contents not described in detail in some embodiments can refer to the same or similar contents in other embodiments.

[0047] It should be noted that, in the description of the present invention, the terms "first", "second", etc. are used for descriptive purposes only and should not be understood as indicating or implying relative importance. In addition, in the description of the present invention, unless otherwise specified, the meaning of "plurality" is at least two.

[0048] Any process or method description in a flowchart or otherwise described herein may be understood to represent a module, segment or portion of code comprising one or more executable instructions for implementing the steps of a specific logical function or process, and the scope of the preferred embodiments of the present invention includes alternative implementations in which functions may be performed out of the order shown or discussed, including performing functions in a substantially simultaneous manner or in the reverse order depending on the functions involved, which should be understood by those skilled in the art to which the embodiments of the present invention pertain.

[0049] It should be understood that various components of the present invention may be implemented using hardware, software, firmware, or a combination thereof. In the above-described embodiments, multiple steps or methods may be implemented using software or firmware stored in a memory and executed by a suitable instruction execution system. For example, if implemented using hardware, as in another embodiment, any one of the following technologies known in the art or a combination thereof may be used: a discrete logic circuit having logic gate circuits for implementing logic functions on data signals, an application-specific integrated circuit having suitable combinational logic gate circuits, a programmable gate array (PGA), a field-programmable gate array (FPGA), etc.

[0050] Those skilled in the art will understand that all or part of the steps in the method of the above embodiment can be completed by instructing related hardware through a program, and the program can be stored in a computer-readable storage medium. When the program is executed, it includes one or a combination of the steps of the method embodiment.

[0051] In addition, the functional units in the various embodiments of the present invention may be integrated into a single processing module, or each unit may exist physically separately, or two or more units may be integrated into a single module. The aforementioned integrated modules may be implemented in the form of hardware or in the form of software functional modules. If the integrated modules are implemented in the form of software functional modules and sold or used as independent products, they may also be stored in a computer-readable storage medium.

[0052] The storage medium mentioned above can be a read-only memory, a magnetic disk or an optical disk, etc.

[0053] Throughout this specification, reference to terms such as "one embodiment," "some embodiments," "examples," "specific examples," or "some examples" means that a specific feature, structure, material, or characteristic described in conjunction with that embodiment or example is included in at least one embodiment or example of the present invention. In this specification, schematic representations of the above terms do not necessarily refer to the same embodiment or example. Furthermore, the specific features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples.

[0054] Although the embodiments of the present invention have been shown and described above, it will be understood that the above embodiments are illustrative and are not to be construed as limitations on the present invention. A person skilled in the art may change, modify, replace and modify the above embodiments within the scope of the present invention.

Claims

1. A method for determining the vibration contribution of a vibration source inside a structure, characterized in that: include: Step S1, obtaining the input load of the vibration source arranged at the vibration source arrangement point inside the target structure and the fast Fourier transform result of the input load, and determining the main load frequency and spectrum of the vibration source; Step S2, obtaining the output time history response of the precision instrument at the precision instrument arrangement point inside the target structure and the fast Fourier transform result of the output time history response, and determining the main vibration frequency and spectrum of the output end of the precision instrument; Step S3, performing analysis based on finite element simulation based on the main load frequency and spectrum of the vibration source and the main vibration frequency and spectrum of the output end of the precision instrument, establishing a dynamic flexibility function of the input end of the vibration source and the output end of the precision instrument, and establishing a transfer matrix based on the dynamic flexibility function; Step S4: Decomposing the frequency spectrum of the output end of the precision instrument according to the transfer matrix to determine the vibration contribution of the internal vibration source of the target structure.

2. The method according to claim 1, characterized in that The step S1 includes: Obtaining a load input time history of the vibration source through a vibration pickup installed on a vibration source arranged below a vibration source arrangement point inside the target structure; Decomposing the load input time history of the vibration source to determine the simple harmonic excitation signal and noise signal of the vibration source; Denoising the load input time history of the decomposed vibration source; Generating a load input time history vector of the internal vibration source of the target structure based on a denoising result of the load input time history of the vibration source; Performing Fourier transform on the load input time history vector to obtain an input load spectrum vector containing a main load frequency of an internal vibration source of the target structure; The main load frequency and spectrum of the vibration source are determined according to the input load spectrum vector.

3. The method according to claim 2, characterized in that The step S2 includes: Acquiring a time history response of an output end of the precision instrument through a vibration pickup installed at a precision instrument arrangement point inside the target structure, and generating a response vector of the output end of the precision instrument; Performing the denoising process on the response vector of the output end of the precision instrument, and performing Fourier transform on the denoising process result to obtain the frequency spectrum vector of the output end of the precision instrument; According to the frequency spectrum vector of the output end of the precision instrument, the main vibration frequency and frequency spectrum of the output end of the precision instrument are determined.

4. The method according to claim 3, characterized in that The step S3 comprises: According to the main load frequency and spectrum of the vibration source, as well as the main vibration frequency and spectrum of the output end of the vibration source, analysis is performed based on finite element simulation to determine the response vectors of the output end of the precision instrument inside the target structure under the excitation of the vibration source load inside the target structure, one by one, and each response vector is recorded as a matrix element of a row of the transfer matrix; All the generated response vectors are arranged in the order of the internal vibration sources of the target structure to generate the transfer matrix.

5. The method according to claim 4, characterized in that The step S4 comprises: Decomposing the frequency spectrum vector at the output end of the precision instrument according to the linear system characteristics and the transfer matrix, and determining an expression of any element in the frequency spectrum vector at the output end of the precision instrument; Determining a vibration contribution calculation formula of any vibration source inside the target structure based on an expression of any element in the frequency spectrum vector at the output end of the precision instrument; The vibration contribution of the internal vibration source of the target structure is determined according to the vibration contribution calculation formula.

6. The method according to claim 5, characterized in that The denoising process includes: Determining a signal to be denoised, and performing Fourier transform on the signal to be denoised; Sort the corresponding main frequencies from large to small according to the transformation amplitude; Performing inverse Fourier transform on the simple harmonic signals corresponding to the main frequencies according to the order of the main frequencies to obtain a reconstructed signal; Subtracting the reconstructed signal from the vibration measured signal of the signal to be denoised to obtain a residual signal; If the autocorrelation function of the residual signal has a non-zero point equal to 0, the residual signal is filtered out.

7. A device for determining the vibration contribution of a vibration source inside a structure, characterized in that: The device comprises a processor and a memory, wherein the processor is connected to the memory: The processor is configured to call and execute the program stored in the memory; The memory is used to store the program, and the program is used at least to execute the method for determining the vibration contribution of the internal vibration source of the structure as described in any one of claims 1-6.