Satellite three-dimensional vibration response data processing method and system

Abstract

The application provides a satellite three-dimensional vibration response data processing method and system, comprising the following steps: S1, arranging three-dimensional acceleration sensors at preset positions on the satellite; S2, performing a single-direction sinusoidal sweep vibration test on the satellite, and collecting vibration response frequency spectrums of the measuring points in three orthogonal directions by using the three-dimensional acceleration sensors; and S3, performing synthesis processing on the obtained three-dimensional vibration response frequency spectrums to obtain a three-dimensional synthesis vibration response amplitude spectrum, which is used for processing three-dimensional vibration response data of components on the satellite. Compared with the traditional three-dimensional independent vibration response amplitude spectrum or the amplitude spectrum of only the main vibration direction, the three-dimensional synthesis vibration response amplitude spectrum obtained by the application can more truly and comprehensively reflect the actual mechanical environment of the products on the satellite in the whole-satellite sinusoidal sweep vibration test, improve the comprehensiveness of the whole-satellite vibration data analysis, and improve the accuracy of the coverage of the whole-satellite mechanical environment by the sinusoidal sweep vibration test conditions of the products on the satellite and the accuracy of the margin evaluation.

Description

Technical Field

[0001] This invention relates to the field of satellite sinusoidal sweep frequency vibration test data processing, specifically to a satellite three-dimensional vibration response data processing method and system. Background Technology

[0002] The sinusoidal sweep vibration test (hereinafter referred to as the vibration test) is a mechanical environment simulation test that must be carried out during the satellite development process as specified in GJB 1027 "Test Requirements for Launch Vehicles, Upper Stages and Spacecraft". Its purpose is to verify the satellite's ability to withstand low-frequency vibration environment during launch and to operate normally. At the same time, the on-board vibration response data obtained from the test can be used to evaluate the coverage and margin of the sinusoidal sweep vibration test conditions of the on-board products on the overall mechanical environment of the satellite.

[0003] In the analysis of whole-satellite vibration test data, there are currently two main methods for processing vibration response data.

[0004] The first method involves processing the vibration response data in the three orthogonal directions separately to obtain three-dimensional independent vibration response amplitude spectra. Because this method yields three-dimensional independent amplitude spectra, it cannot directly reflect the overall mechanical environment and vibration response at the measurement point. In practice, the vibration responses in the three directions are typically treated as occurring independently. For example, when evaluating the coverage and margin of vibration test conditions for onboard products, a coverage analysis of the vibration response amplitude spectra in each direction is generally performed directly. This is equivalent to treating the vibration responses in each direction as occurring at different times, which is inconsistent with reality. Furthermore, the vibration response in each direction is only the component of the actual total vibration response at that location in the three orthogonal directions; the vibration response in each direction is necessarily smaller than the total vibration response. These two reasons result in an inaccurate and unrealistic analysis of the actual mechanical environment experienced by onboard products during the overall satellite vibration test, leading to inaccurate evaluations of the coverage and margin of the vibration test conditions for the overall satellite mechanical environment.

[0005] Another method is to process only the vibration response data in the principal vibration direction to obtain the amplitude spectrum of the vibration response in that direction. The vibration response in the principal vibration direction (i.e., the same direction as the test direction) is generally larger than in the other two directions, so this method directly ignores the vibration responses in the other two directions. This method is essentially a simplification of the first method, and therefore suffers from the same drawbacks. Furthermore, for some large-mass products, due to complex local modes, phenomena such as torsional vibration and oscillation often occur, causing the vibration response in some frequency bands outside the principal vibration direction to be greater than that in the principal vibration direction. In such cases, the error in evaluating the overall mechanical environment and vibration response using this method will be larger than that of the first method.

[0006] To address the shortcomings of existing technologies, this invention provides a satellite triaxial vibration response data processing method that obtains the triaxial composite vibration response amplitude spectrum by acquiring and synthesizing triaxial vibration response spectra. Currently, no similar descriptions or reports have been found, nor have similar domestic or international materials been collected.

[0007] A search for "vibration data processing method" yielded Chinese patent application CN201510280536.8, "A Method, Apparatus, and System for Processing Vibration Data of Equipment." This patent describes a method for processing vibration data of equipment. The method includes: receiving vibration data of an equipment collected in a first sampling period, the vibration data being associated with the equipment's rotational speed and rotational speed fluctuation in that first sampling period; determining whether the equipment's rotational speed and rotational speed fluctuation associated with the vibration data satisfy a first rule; if the rotational speed and rotational speed fluctuation associated with the vibration data satisfy the first rule, then caching the vibration data; repeating the above receiving, determining, and caching steps in subsequent sampling periods until vibration data for another sampling period satisfying the first rule is cached; and selecting one vibration data from the vibration data in the two sequentially cached sampling periods according to a second rule associated with at least rotational speed and rotational speed fluctuation, and retaining the selected vibration data. This is completely different from the content of this patent. Another Chinese patent, CN201210341427.9, "A Data Processing Method in a Distributed Optical Fiber Vibration System," describes a data processing method for a distributed optical fiber vibration sensor, which can effectively solve the problem of spectrum tilt and achieve threshold-triggered alarm for disturbance information. First, optical path noise is eliminated by normalizing the two signals. Second, wavelet multi-layer decomposition and reconstruction technology is used to reduce the noise of the defect spectrum. Then, curve fitting is used to correct the tilted defect spectrum to achieve a flat defect spectrum, which improves the accuracy of the defect point and the positioning accuracy of the disturbance measurement. This is completely different from the content of this patent invention. Chinese patent application number CN201410039208.4, entitled "A Data Processing System and Method," applies to a cloud environment for information processing / exchange in a distributed information system. When using the data processing system of this invention to perform a data processing method, firstly, a basic information directory is set, where a basic information directory with conditions / rules is defined according to a specific target. Next, reciprocal interaction is performed, where users repeatedly interact with the basic information directory. Through this reciprocal interaction, condition / rule iteration is carried out, generating a derived information directory and accompanying changes in group composition. Finally, a specific group composed of users is formed for the specific target. Here, specific groups composed of users are formed for each of the basic information directory and / or derived information directory related to the specific target, and users agree on transaction conditions for the specific target. This is completely different from the content of this patent invention.

[0008] Patent document CN110929342A (application number: CN201911228109.X) discloses a method for processing helicopter vibration response data. It obtains data blocks that meet analysis requirements by preprocessing measured vibration data and dividing it into time periods; it obtains the main influencing frequencies based on the specific location of the components; it uses FFT transformation to obtain the vibration amplitude at each main influencing frequency; it combines the helicopter flight spectrum and uses statistical methods to obtain the distribution of amplitudes at each main influencing frequency; it uses parameter estimation and hypothesis testing to obtain the probability distribution and obtain vibration response amplitudes that meet certain confidence and reliability requirements. However, this patent mainly addresses the problem of compiling helicopter vibration response spectra, which is completely different from the content of this invention. Summary of the Invention

[0009] To address the shortcomings of existing technologies, the purpose of this invention is to provide a method and system for processing satellite three-dimensional vibration response data.

[0010] A satellite triaxial vibration response data processing method provided by the present invention includes:

[0011] Step S1: Deploy a triaxial accelerometer at a predetermined location on the satellite;

[0012] Step S2: Conduct a single-direction sinusoidal frequency sweep vibration test on the satellite, and use a triaxial accelerometer to collect the vibration response spectrum of the measuring point in three orthogonal directions;

[0013] Step S3: The obtained triaxial vibration response spectrum is synthesized to obtain the triaxial synthesized vibration response amplitude spectrum, which is used to process the triaxial vibration response data of the on-board components.

[0014] Preferably, in step S1:

[0015] The deployed triaxial accelerometers can simultaneously acquire the vibration response of the measured part in three orthogonal directions.

[0016] Preferably, in step S2:

[0017] The spectrum consists of the real part and the imaginary part.

[0018] Preferably, in step S2:

[0019] The vibration response spectra of the measurement point in the three orthogonal directions (x, y, z) are as follows:

[0020] X(f) = a x (f)+ib x (f)

[0021] Y(f) = a y (f)+ib y (f)

[0022] Z(f) = a z (f)+ib z (f)

[0023] In the formula, f is the sinusoidal scanning frequency, i is the imaginary number, and a x (f), a y (f), a z (f) represents the real part of the spectrum, b x (f), b y (f), b z (f) Imaginary part spectrum;

[0024] The amplitude spectra of the vibration response in the x, y, and z directions are as follows:

[0025]

[0026]

[0027]

[0028] The phase spectra of the vibration response in the x, y, and z directions are as follows:

[0029]

[0030]

[0031]

[0032] Preferably, in step S3:

[0033] The triaxial vibration response spectrum obtained in step 2 is synthesized to obtain the triaxial synthesized vibration response amplitude spectrum.

[0034]

[0035] In the formula, A m , Let A be the amplitude spectrum and phase spectrum of the combined vibrational response in the x and y directions, respectively; α be the angle between the combined vibration in the x and y directions and the x-axis; β be the angle between the combined vibration in the x, y, and z directions and the combined vibration in the x and y directions; and A be the amplitude spectrum and phase spectrum of the combined vibrational response in the x and y directions, respectively. z The amplitude spectrum of the vibration response in the z-direction. The phase spectrum of the vibration response in the z-direction;

[0036]

[0037]

[0038]

[0039]

[0040] α is the angle between the combined vibration in the x and y directions and the x-axis, A x For the x-axis vibration response amplitude spectrum, A y The amplitude spectrum of the vibration response in the y-direction; The phase spectrum of the vibration response in the x-direction. The phase spectrum of the vibrational response in the y-direction; β is the angle between the combined vibrations in the x, y, and z directions and the combined vibrations in the x and y directions, A m , The amplitude and phase spectra of the synthesized vibrational responses in the x and y directions are respectively, A z The amplitude spectrum of the vibration response in the z-direction. The phase spectrum of the vibration response in the z-direction;

[0041] First calculate α; then substitute α into the calculation to get A. m , A m , Substitute into the calculation to obtain β; then use A m , Substituting β into the calculation yields the final triaxial composite vibration response amplitude spectrum A. n .

[0042] A satellite triaxial vibration response data processing system provided by the present invention includes:

[0043] Module M1: A three-dimensional accelerometer sensor is placed at a predetermined location on the satellite;

[0044] Module M2: Conducts a single-direction sinusoidal sweep vibration test on the satellite, and uses a triaxial accelerometer to collect the vibration response spectrum of the measuring point in three orthogonal directions;

[0045] Module M3: Synthesizes the obtained triaxial vibration response spectrum to obtain the triaxial synthetic vibration response amplitude spectrum, which is used to process the triaxial vibration response data of on-board components.

[0046] Preferably, in module M1:

[0047] The deployed triaxial accelerometers can simultaneously acquire the vibration response of the measured part in three orthogonal directions.

[0048] Preferably, in module M2:

[0049] The spectrum consists of the real part and the imaginary part.

[0050] Preferably, in module M2:

[0051] The vibration response spectra of the measurement point in the three orthogonal directions (x, y, z) are as follows:

[0052] X(f) = a x (f)+ib x (f)

[0053] Y(f) = a y (f)+ib y (f)

[0054] Z(f) = a z (f)+ib z (f)

[0055] In the formula, f is the sinusoidal scanning frequency, i is the imaginary number, and a x (f), a y (f), a z (f) represents the real part of the spectrum, b x (f), b y (f), b z (f) Imaginary part spectrum;

[0056] The amplitude spectra of the vibration response in the x, y, and z directions are as follows:

[0057]

[0058]

[0059]

[0060] The phase spectra of the vibration response in the x, y, and z directions are as follows:

[0061]

[0062]

[0063]

[0064] Preferably, in module M3:

[0065] The triaxial vibration response spectrum obtained in step 2 is synthesized to obtain the triaxial synthesized vibration response amplitude spectrum.

[0066]

[0067] In the formula, A m , Let A be the amplitude spectrum and phase spectrum of the combined vibrational response in the x and y directions, respectively; α be the angle between the combined vibration in the x and y directions and the x-axis; β be the angle between the combined vibration in the x, y, and z directions and the combined vibration in the x and y directions; and A be the amplitude spectrum and phase spectrum of the combined vibrational response in the x and y directions, respectively. z The amplitude spectrum of the vibration response in the z-direction. The phase spectrum of the vibration response in the z-direction;

[0068]

[0069]

[0070]

[0071]

[0072] α is the angle between the combined vibration in the x and y directions and the x-axis, A x For the x-axis vibration response amplitude spectrum, A y The amplitude spectrum of the vibration response in the y-direction; The phase spectrum of the vibration response in the x-direction. The phase spectrum of the vibrational response in the y-direction; β is the angle between the combined vibrations in the x, y, and z directions and the combined vibrations in the x and y directions, A m , The amplitude and phase spectra of the synthesized vibrational responses in the x and y directions are respectively, A z The amplitude spectrum of the vibration response in the z-direction. The phase spectrum of the vibration response in the z-direction;

[0073] First calculate α; then substitute α into the calculation to get A. m , A m , Substitute into the calculation to obtain β; then use A m , Substituting β into the calculation yields the final triaxial composite vibration response amplitude spectrum A. n .

[0074] Compared with the prior art, the present invention has the following beneficial effects:

[0075] 1. Compared with the traditional three-dimensional independent vibration response amplitude spectrum or the amplitude spectrum of vibration response only in the main vibration direction, the three-dimensional composite vibration response amplitude spectrum obtained by this invention can more realistically and comprehensively reflect the actual mechanical environment experienced by the on-board products in the whole-satellite sinusoidal sweep frequency vibration test.

[0076] 2. This invention improves the comprehensiveness of the analysis of whole-satellite vibration data, and enhances the coverage of the sinusoidal sweep frequency vibration test conditions of on-board products on the whole-satellite mechanical environment and the accuracy of margin evaluation. Attached Figure Description

[0077] Other features, objects, and advantages of the present invention will become more apparent from the following detailed description of non-limiting embodiments with reference to the accompanying drawings:

[0078] Figure 1 This is a flowchart of the steps of the present invention;

[0079] Figure 2 This is a schematic diagram of the time-domain synthesis of vibration response in the x, y, and z directions;

[0080] Figure 3 This is a schematic diagram of the synthesis of time-domain vibration response in the x and y directions;

[0081] Figure 4 This is a schematic diagram of the time-domain synthesis of the vibration response in the x′ (x and y directions) and z-direction directions. Detailed Implementation

[0082] The present invention will now be described in detail with reference to specific embodiments. These embodiments will help those skilled in the art to further understand the present invention, but do not limit the invention in any way. It should be noted that those skilled in the art can make several changes and improvements without departing from the concept of the present invention. These all fall within the scope of protection of the present invention.

[0083] Example 1:

[0084] According to the present invention, a satellite triaxial vibration response data processing method is provided, such as... Figures 1-4 As shown, it includes:

[0085] Step S1: Deploy a triaxial accelerometer at a predetermined location on the satellite;

[0086] Step S2: Conduct a single-direction sinusoidal frequency sweep vibration test on the satellite, and use a triaxial accelerometer to collect the vibration response spectrum of the measuring point in three orthogonal directions;

[0087] Step S3: The obtained triaxial vibration response spectrum is synthesized to obtain the triaxial synthesized vibration response amplitude spectrum, which is used to process the triaxial vibration response data of the on-board components.

[0088] Specifically, in step S1:

[0089] The deployed triaxial accelerometers can simultaneously acquire the vibration response of the measured part in three orthogonal directions.

[0090] Specifically, in step S2:

[0091] The spectrum consists of the real part and the imaginary part.

[0092] Specifically, in step S2:

[0093] The vibration response spectra of the measurement point in the three orthogonal directions (x, y, z) are as follows:

[0094] X(f) = a x (f)+ib x (f)

[0095] Y(f) = a y (f)+ib y (f)

[0096] Z(f) = a z (f)+ib z (f)

[0097] In the formula, f is the sinusoidal scanning frequency, i is the imaginary number, and a x (f), a y (f), a z (f) represents the real part of the spectrum, b x (f), b y (f), b z (f) Imaginary part spectrum;

[0098] The amplitude spectra of the vibration response in the x, y, and z directions are as follows:

[0099]

[0100]

[0101]

[0102] The phase spectra of the vibration response in the x, y, and z directions are as follows:

[0103]

[0104]

[0105]

[0106] Specifically, in step S3:

[0107] The triaxial vibration response spectrum obtained in step 2 is synthesized to obtain the triaxial synthesized vibration response amplitude spectrum.

[0108]

[0109] In the formula, A m , Let A be the amplitude spectrum and phase spectrum of the combined vibrational response in the x and y directions, respectively; α be the angle between the combined vibration in the x and y directions and the x-axis; β be the angle between the combined vibration in the x, y, and z directions and the combined vibration in the x and y directions; and A be the amplitude spectrum and phase spectrum of the combined vibrational response in the x and y directions, respectively. z The amplitude spectrum of the vibration response in the z-direction. The phase spectrum of the vibration response in the z-direction;

[0110]

[0111]

[0112]

[0113]

[0114] α is the angle between the combined vibration in the x and y directions and the x-axis, A x For the x-axis vibration response amplitude spectrum, A y The amplitude spectrum of the vibration response in the y-direction; The phase spectrum of the vibration response in the x-direction. The phase spectrum of the vibrational response in the y-direction; β is the angle between the combined vibrations in the x, y, and z directions and the combined vibrations in the x and y directions, A m , The amplitude and phase spectra of the synthesized vibrational responses in the x and y directions are respectively, A z The amplitude spectrum of the vibration response in the z-direction. The phase spectrum of the vibration response in the z-direction;

[0115] First calculate α; then substitute α into the calculation to get A. m , A m , Substitute into the calculation to obtain β; then use A m , Substituting β into the calculation yields the final triaxial composite vibration response amplitude spectrum A. n .

[0116] This invention can also be used for processing triaxial vibration response data of large components on spacecraft.

[0117] Example 2:

[0118] Example 2 is a preferred embodiment of Example 1, and is used to illustrate the present invention in more detail.

[0119] Those skilled in the art can understand the satellite triaxial vibration response data processing method provided by the present invention as a specific implementation of the satellite triaxial vibration response data processing system, that is, the satellite triaxial vibration response data processing system can be implemented by executing the steps of the satellite triaxial vibration response data processing method.

[0120] A satellite triaxial vibration response data processing system provided by the present invention includes:

[0121] Module M1: A three-dimensional accelerometer sensor is placed at a predetermined location on the satellite;

[0122] Module M2: Conducts a single-direction sinusoidal sweep vibration test on the satellite, and uses a triaxial accelerometer to collect the vibration response spectrum of the measuring point in three orthogonal directions;

[0123] Module M3: Synthesizes the obtained triaxial vibration response spectrum to obtain the triaxial synthetic vibration response amplitude spectrum, which is used to process the triaxial vibration response data of on-board components.

[0124] Specifically, in module M1:

[0125] The deployed triaxial accelerometers can simultaneously acquire the vibration response of the measured part in three orthogonal directions.

[0126] Specifically, in module M2:

[0127] The spectrum consists of the real part and the imaginary part.

[0128] Specifically, in module M2:

[0129] The vibration response spectra of the measurement point in the three orthogonal directions (x, y, z) are as follows:

[0130] X(f) = a x (f)+ib x (f)

[0131] Y(f) = a y (f)+ib y (f)

[0132] Z(f) = a z (f)+ib z (f)

[0133] In the formula, f is the sinusoidal scanning frequency, i is the imaginary number, and a x (f), a y (f), a z (f) represents the real part of the spectrum, b x (f), b y (f), b z (f) Imaginary part spectrum;

[0134] The amplitude spectra of the vibration response in the x, y, and z directions are as follows:

[0135]

[0136]

[0137]

[0138] The phase spectra of the vibration response in the x, y, and z directions are as follows:

[0139]

[0140]

[0141]

[0142] Specifically, in module M3:

[0143] The triaxial vibration response spectrum obtained in step 2 is synthesized to obtain the triaxial synthesized vibration response amplitude spectrum.

[0144]

[0145] In the formula, A m , Let A be the amplitude spectrum and phase spectrum of the combined vibrational response in the x and y directions, respectively; α be the angle between the combined vibration in the x and y directions and the x-axis; β be the angle between the combined vibration in the x, y, and z directions and the combined vibration in the x and y directions; and A be the amplitude spectrum and phase spectrum of the combined vibrational response in the x and y directions, respectively. z The amplitude spectrum of the vibration response in the z-direction. The phase spectrum of the vibration response in the z-direction;

[0146]

[0147]

[0148]

[0149]

[0150] α is the angle between the combined vibration in the x and y directions and the x-axis, A x For the x-axis vibration response amplitude spectrum, A y The amplitude spectrum of the vibration response in the y-direction; The phase spectrum of the vibration response in the x-direction. The phase spectrum of the vibrational response in the y-direction; β is the angle between the combined vibrations in the x, y, and z directions and the combined vibrations in the x and y directions, A m , The amplitude and phase spectra of the synthesized vibrational responses in the x and y directions are respectively, A z The amplitude spectrum of the vibration response in the z-direction. The phase spectrum of the vibration response in the z-direction;

[0151] First calculate α; then substitute α into the calculation to get A. m , A m , Substitute into the calculation to obtain β; then use A m , Substituting β into the calculation yields the final triaxial composite vibration response amplitude spectrum A. n .

[0152] Example 3:

[0153] Example 3 is a preferred example of Example 1, and is used to illustrate the present invention in more detail.

[0154] To address the shortcomings of existing technologies, the purpose of this invention is to provide a method for processing satellite triaxial vibration response data. This method, by acquiring and synthesizing triaxial vibration response spectra, obtains a triaxial synthesized vibration response amplitude spectrum. This provides a more realistic and comprehensive reflection of the actual mechanical environment experienced by onboard products during whole-satellite sinusoidal sweep vibration tests, improving the comprehensiveness of whole-satellite vibration data analysis and enhancing the coverage and accuracy of margin evaluation of the sinusoidal sweep vibration test conditions on the whole-satellite mechanical environment. Furthermore, this invention can also be used for processing triaxial vibration response data of large onboard components.

[0155] In this embodiment, the satellite triaxial vibration response data processing method provided by the present invention includes the following steps:

[0156] Step 1: Deploy an accelerometer on the satellite at the area of ​​interest, capable of simultaneously acquiring the vibration response of the measured area in three orthogonal directions.

[0157] Step 2: Conduct a single-direction sinusoidal frequency sweep vibration test on the satellite. During the test, use a set of triaxial accelerometers to collect the vibration response spectrum of the measuring points in three orthogonal directions. This spectrum should include both the real and imaginary parts.

[0158] The vibration response spectra of the measurement point in the three orthogonal directions (x, y, z) are recorded as follows:

[0159] X(f) = a x (f)+ib x (f)

[0160] Y(f) = a y (f)+ib y (f)

[0161] Z(f) = a z (f)+ib z (f)

[0162] In the formula, f is the sinusoidal scanning frequency, i is the imaginary number, and a x (f), a y (f), a z (f) represents the real part of the spectrum, b x (f), b y (f), b z (f) Imaginary part spectrum. Then the amplitude spectra of the vibration response in the x, y, and z directions are respectively...

[0163]

[0164]

[0165]

[0166] The phase spectra of the vibration response in the x, y, and z directions are respectively

[0167]

[0168]

[0169]

[0170] Step 3: Perform synthesis processing on the triaxial vibration response spectrum obtained in Step 2 to obtain the triaxial synthesized vibration response amplitude spectrum.

[0171] Since direct synthesis in the frequency domain is not possible, it is necessary to convert to the time domain for relevant derivation. Based on the triaxial vibration response spectrum, the time-domain signals of the triaxial vibration response at a sinusoidal scanning frequency f are as follows:

[0172]

[0173]

[0174]

[0175] A schematic diagram of the time-domain synthesis of vibration response in the x, y, and z directions is shown below. Figure 2 In the figure, the combined vibration of x(t) and y(t) is m(t), with an angle α between m(t) and x(t); the combined vibration of m(t) and z(t) is n(t), with an angle β between n(t) and m(t). This combination process consists of two steps:

[0176] First, the vibration response spectra in the x and y directions are synthesized. A schematic diagram of the time-domain synthesis of the x and y direction vibration responses is shown below. Figure 3 In the diagram, rotating the OX axis counterclockwise by an angle α yields the OX′ axis. The sum of the components of x(t) and y(t) along the x′ direction is m(t).

[0177]

[0178] In the formula, A m , Let m(t) be the amplitude spectrum and phase spectrum, respectively.

[0179]

[0180]

[0181] When m(t) is the combined vibration of x(t) and y(t), A mTo achieve the maximum value (because the synthesized vibration is always greater than its components), that is, we have...

[0182]

[0183] In the formula,

[0184]

[0185] From this we can obtain

[0186]

[0187] In summary, we can calculate α; substituting α into the calculation yields A. m , At this time, A m , That is, the amplitude spectrum and phase spectrum of the combined vibration response in the x and y directions.

[0188] Then, similarly, the vibration response spectra in the x′ and z directions are synthesized, where the x′ direction represents the direction of the synthesized vibration m(t) in the x and y directions. A schematic diagram of the time-domain synthesis of the vibration response in the x′ and z directions is shown below. Figure 4 In the diagram, rotating the OX′ axis counterclockwise by an angle β, we obtain the OX″ coordinate system. The sum of the components of m(t) and z(t) in the x″ direction is n(t).

[0189]

[0190] In the formula, A n , Let n(t) be the amplitude spectrum and phase spectrum, respectively.

[0191]

[0192]

[0193] When n(t) is the combined vibration of m(t) and z(t), A n To achieve the maximum value (because the synthesized vibration is always greater than its components), that is, we have...

[0194]

[0195] In the formula,

[0196]

[0197] From this we can obtain

[0198]

[0199] In conclusion, A m , Substitute the values ​​to calculate β; substitute β into the calculation to calculate A. n , At this time, A n , That is, the amplitude spectrum and phase spectrum of the synthesized vibration response in the x, y, and z directions.

[0200] Those skilled in the art will understand that, in addition to implementing the system, apparatus, and their modules provided by this invention in purely computer-readable program code, the same program can be implemented in the form of logic gates, switches, application-specific integrated circuits, programmable logic controllers, and embedded microcontrollers by logically programming the method steps. Therefore, the system, apparatus, and their modules provided by this invention can be considered a hardware component, and the modules included therein for implementing various programs can also be considered structures within the hardware component; alternatively, modules for implementing various functions can be considered both software programs implementing the method and structures within the hardware component.

[0201] Specific embodiments of the present invention have been described above. It should be understood that the present invention is not limited to the specific embodiments described above, and those skilled in the art can make various changes or modifications within the scope of the claims, which do not affect the essence of the present invention. Unless otherwise specified, the embodiments and features described in this application can be arbitrarily combined with each other.

Claims

1. A method for processing satellite three-dimensional vibration response data, characterized in that, include: Step S1: Deploy a triaxial accelerometer at a predetermined location on the satellite; Step S2: Conduct a single-direction sinusoidal frequency sweep vibration test on the satellite, and use a triaxial accelerometer to collect the vibration response spectrum of the measuring point in three orthogonal directions; Step S3: The obtained triaxial vibration response spectrum is synthesized to obtain the triaxial composite vibration response amplitude spectrum, which is used to process the triaxial vibration response data of the on-board components; In step S3: The triaxial vibration response spectrum obtained in step 2 is synthesized to obtain the triaxial synthesized vibration response amplitude spectrum. （7） In the formula, , They are respectively , The amplitude and phase spectra of the synthesized vibrational response, for , Towards the synthesis vibration and The included angle of the axis, for , , Towards the synthesis vibration and , The angle between the synthesized vibrations, for To the vibration response amplitude spectrum, for To the phase spectrum of vibration response; （8） （9） （10） （11） for , Towards the synthesis vibration and The included angle of the axis, for To the vibration response amplitude spectrum, for To the vibration response amplitude spectrum; for To the phase spectrum of vibration response, for To the phase spectrum of vibration response; for , , Towards the synthesis vibration and , The angle between the synthesized vibrations, , They are respectively , The amplitude and phase spectra of the synthesized vibrational response, for To the vibration response amplitude spectrum, for To the phase spectrum of vibration response; First calculate Then Substitute into the calculation to obtain , ;Will , Substituting into the calculation yields ;Will , , Substituting the values ​​into the calculation, we obtain the final triaxial composite vibration response amplitude spectrum. .

2. The satellite triaxial vibration response data processing method according to claim 1, characterized in that, In step S1: The deployed triaxial accelerometers can simultaneously acquire the vibration response of the measured part in three orthogonal directions.

3. The satellite triaxial vibration response data processing method according to claim 1, characterized in that, In step S2: The spectrum consists of the real part and the imaginary part.

4. The satellite triaxial vibration response data processing method according to claim 3, characterized in that, In step S2: Record the measurement point at , , The vibration response spectra in the three orthogonal directions are as follows: In the formula, The frequency is a sinusoidal scanning frequency. It is an imaginary number. , , For the real part of the spectrum, , , Imaginary part spectrum; , , The vibration response amplitude spectra are as follows: （1） （2） （3） , , The phase spectra of the vibration response are as follows: （4） （5） （6）。 5. A satellite three-dimensional vibration response data processing system, characterized in that, include: Module M1: A three-dimensional accelerometer sensor is placed at a predetermined location on the satellite; Module M2: Conducts a single-direction sinusoidal sweep vibration test on the satellite, and uses a triaxial accelerometer to collect the vibration response spectrum of the measuring point in three orthogonal directions; Module M3: Synthesizes the obtained triaxial vibration response spectrum to obtain the triaxial composite vibration response amplitude spectrum, which is used to process the triaxial vibration response data of on-board components; In module M3: The triaxial vibration response spectrum obtained in step 2 is synthesized to obtain the triaxial synthesized vibration response amplitude spectrum. （7） In the formula, , They are respectively , The amplitude and phase spectra of the synthesized vibrational response, for , Towards the synthesis vibration and The included angle of the axis, for , , Towards the synthesis vibration and , The angle between the synthesized vibrations, for To the vibration response amplitude spectrum, for To the phase spectrum of vibration response; （8） （9） （10） （11） for , Towards the synthesis vibration and The included angle of the axis, for To the vibration response amplitude spectrum, for To the vibration response amplitude spectrum; for To the phase spectrum of vibration response, for To the phase spectrum of vibration response; for , , Towards the synthesis vibration and , The angle between the synthesized vibrations, , They are respectively , The amplitude and phase spectra of the synthesized vibrational response, for To the vibration response amplitude spectrum, for To the phase spectrum of vibration response; First calculate Then Substitute into the calculation to obtain , ;Will , Substituting into the calculation yields ;Will , , Substituting the values ​​into the calculation, we obtain the final triaxial composite vibration response amplitude spectrum. .

6. The satellite triaxial vibration response data processing system according to claim 5, characterized in that, In module M1: The deployed triaxial accelerometers can simultaneously acquire the vibration response of the measured part in three orthogonal directions.

7. The satellite triaxial vibration response data processing system according to claim 5, characterized in that, In module M2: The spectrum consists of the real part and the imaginary part.

8. The satellite triaxial vibration response data processing system according to claim 7, characterized in that, In module M2: Record the measurement point at , , The vibration response spectra in the three orthogonal directions are as follows: In the formula, The frequency is a sinusoidal scanning frequency. It is an imaginary number. , , For the real part of the spectrum, , , Imaginary part spectrum; , , The vibration response amplitude spectra are as follows: （1） （2） （3） , , The phase spectra of the vibration response are as follows: （4） （5） （6）。

Images