CN105628416A - Upper stage single machine sine scan test condition refining design method - Google Patents
Upper stage single machine sine scan test condition refining design method Download PDFInfo
- Publication number
- CN105628416A CN105628416A CN201410637581.XA CN201410637581A CN105628416A CN 105628416 A CN105628416 A CN 105628416A CN 201410637581 A CN201410637581 A CN 201410637581A CN 105628416 A CN105628416 A CN 105628416A
- Authority
- CN
- China
- Prior art keywords
- test
- cla
- maximum
- amplitude
- frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Testing Of Engines (AREA)
- Monitoring And Testing Of Transmission In General (AREA)
Abstract
The invention discloses an upper stage single machine sine scan test condition refining design method. With combination of theory and experiment, upper stage/ rocket interface acceleration of satellite and rocket load coupled analytic results is used as forcing function, the ratio of acceleration response to excitation at an upper stage system-level sine scan test single machine is used as a transmission characteristic, and the forcing function is multiplied by the transmission characteristic to obtain the maximum single machine environmental envelope. By use of the method, environmental envelopes corresponding to different frequencies can be provided so that the test condition is finer and the severity of environmental condition design is avoided, and the test condition is accurate and is with certain margin. Furthermore, the system-level sine scan test data processing and test condition envelopes can both be achieved by digital batch processing and thus the efficiency and fine degree of mechanical environmental condition design can be effectively increased.
Description
Technical field
The present invention relates to Upper Stage unit dynamic environmental test condition design field, particularly relate to a kind of Upper Stage unit sine sweep test condition minute design method.
Background technology
Unit dynamic environmental test condition is the important evidence of Spacecraft guidance and control, production and verification experimental verification, dynamic environmental test condition formulate reasonability and fine degree, directly affect the weight of spacecraft unit, funds, Development Schedule and flight reliability.
Upper Stage, between rocket and spacecraft, is the uppermost one-level of carrier rocket. Owing to Upper Stage many stars allotter layout is complicated, task time between rocket and spacecraft, develop the division of labor and belong to the features such as carrier rocket, the unit dynamic environmental test condition formulating method of Upper Stage can not simply indiscriminately imitate the design pattern of carrier rocket or spacecraft.
Unit sine sweep test is the one of unit dynamic environmental test, unit sine sweep test condition according to mainly include close model test data analogize, load coupling analysis result, system-level low-frequency vibration result of the test etc., above method is used alone all certain limitation:
The border impedance of system-level low-frequency vibration test is bigger with Live Flying state difference, cause that unit low-frequency vibration experimental condition is harsher, Upper Stage unit sine sweep test condition is formulated according to system-level low-frequency vibration result of the test, unit qualification magnitude can be made to be likely to be breached 20g, and (g is acceleration of gravity, normal value 9.8m/s2), this will make unit more for meeting environmental suitability weightening finish;
Low-frequency vibration experimental condition in interface phase of adjustment that satellite and the rocket load coupling analysis result being primarily adapted for use in, for unit (on device) on arrow, owing to its locally fine point details is difficult to simulation accurately, generally not in application load coupling analysis result unit place response envelope as unit sine sweep test condition;
And adopting close model to analogize mode, it is not accurate enough that the limitation owing to deducing may result in band segment experimental condition.
Summary of the invention
Present invention solves the technical problem that it is formulate experimental condition according to the formulating method of unit sine sweep test condition in prior art to there is the problems such as experimental condition is harsher, band segment experimental condition is not accurate enough, so provide a kind of can be more fine formulate experimental condition and it can be avoided that the Upper Stage unit sine sweep test condition minute design method of severity of environmental condition design.
In order to solve above-mentioned technical problem, the technical solution used in the present invention is as follows:
Upper Stage unit sine sweep test condition minute design method, comprises the following steps:
S01: in the frequency range needed for the sine sweep test condition of the unit that sets objectives, provides the Frequency point of certain intervals;
Obtain Upper Stage/rocket interface X/Y/Z tri-directional acceleration shock response spectrum-frequency curve under different operating modes in satellite and the rocket load coupling analysis result, acceleration shock response spectrum-frequency curve under all operating modes on same direction is carried out maximum envelope, find the maximum in the acceleration shock response spectrum under operating modes on same direction, the same described Frequency point of correspondence, different respectively, obtain Upper Stage/rocket interface X/Y/Z tri-directional acceleration shock response spectrum maximum-frequency curve, be expressed as Xcla��YclaAnd Zcla;
Under acquisition appointment magnitude, in Stage Vehicle System level sine sweep test, X/Y/Z tri-tests the acceleration transmission ratio amplitude-frequency curve in tri-directions of X/Y/Z at target unit place under direction, is expressed as:
X-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
X-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
X-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
Described appointment magnitude selects between 0.05g��0.2g;
S02: calculate the X obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Xtest��Y-Xtest��Z-XtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three result of calculations-frequency curve Xcla-X-Xtest��Xcla-Y-XtestAnd Xcla-Z-Xtest;
Calculate the Y obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Ytest��Y-Ytest��Z-YtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three result of calculations-frequency curve Ycla-X-Ytest��Ycla-Y-YtestAnd Ycla-Z-Ytest;
Calculate the Z obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Ztest��Y-Ztest��Z-ZtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three result of calculations-frequency curve Zcla-X-Ztest��Zcla-Y-ZtestAnd Zcla-Z-Ztest;
Described specified multiple selects between 1.20 times��1.40 times, and with fixing in a unit sine sweep test condition design;
S03: X will be obtained in step S02cla-X-Xtest��Ycla-X-Ytest��Zcla-X-ZtestArticle three, curve carries out envelope on frequency domain, finds the maximum in described result of calculation on three curves of corresponding same described Frequency point respectively, obtains result of calculation maximum-frequency curve as the X frequency domain maximum envelope to test;
Step S02 will obtain Xcla-Y-Xtest��Ycla-Y-Ytest��Zcla-Y-ZtestArticle three, curve carries out envelope on frequency domain, finds the maximum in described result of calculation on three curves of corresponding same described Frequency point respectively, obtains the frequency domain maximum envelope that result of calculation maximum-frequency curve is tested as Y-direction;
Step S02 will obtain Xcla-Z-Xtest��Ycla-Z-Ytest��Zcla-Z-ZtestArticle three, curve carries out envelope on frequency domain, finds the maximum in described result of calculation on three curves of corresponding same described Frequency point respectively, obtains the frequency domain maximum envelope that result of calculation maximum-frequency curve is tested as Z-direction;
S04: be smoothed by the described frequency domain maximum envelope drawn, obtains the sine sweep test condition covering under this correspondence direction of target unit of whole described frequency domain maximum envelope.
The present invention also provides for another kind of Upper Stage unit sine sweep test condition minute design method, comprises the following steps:
S01: in the frequency range needed for the unit sine sweep test condition that sets objectives, provide the Frequency point of certain intervals;
Obtain Upper Stage/rocket interface X/Y/Z tri-directional acceleration shock response spectrum-frequency curve under different operating modes in satellite and the rocket load coupling analysis result, acceleration shock response spectrum-frequency curve under all operating modes on same direction is carried out maximum envelope, find the maximum in the acceleration shock response spectrum under operating modes on same direction, the same described Frequency point of correspondence, different respectively, obtain Upper Stage/rocket interface X/Y/Z tri-directional acceleration shock response spectrum maximum-frequency curve, be expressed as Xcla��YclaAnd Zcla;
Under acquisition appointment magnitude, in Stage Vehicle System level sine sweep test, X/Y/Z tri-tests the acceleration transmission ratio amplitude-frequency curve in tri-directions of X/Y/Z at target unit place under direction, is expressed as:
X-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
X-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
X-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
Described appointment magnitude selects between 0.05g��0.2g;
S02: calculate the X obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Xtest��Y-Xtest��Z-XtestOn acceleration transmission than the product of amplitude, obtain three products-frequency curve (Xcla-X-Xtest)�@��(Xcla-Y-Xtest) and (Xcla-Z-Xtest);
Calculate the Y obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Ytest��Y-Ytest��Z-YtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three products-frequency curve (Ycla-X-Yest)�@��(Ycla-Y-Ytest) and (Ycla-Z-Ytest);
Calculate the Z obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Ztest��Y-Ztest��Z-ZtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three products-frequency curve (Zcla-X-Ztest)�@��(Zcla-Y-Ztest) and (Zcla-Z-Ztest); Described specified multiple selects between 1.20 times��1.40 times, and with fixing in a unit sine sweep test condition design;
S03: (X will be obtained in step S02cla-X-Xtest)�@��(Ycla-X-Ytest)�@��(Zcla-X-Ztest) three curves carry out envelope on frequency domain, find the maximum in described product on three curves of corresponding same described Frequency point respectively, and calculate the result of calculation of the specified multiple of the maximum of described product, obtain result of calculation maximum-frequency curve as the X frequency domain maximum envelope to test;
Step S02 will obtain (Xcla-Y-Xtest)�@��(Ycla-Y-Ytest)�@(Zcla-Y-Ztest) three curves carry out envelope on frequency domain, find the maximum in described product on three curves of corresponding same described Frequency point respectively, and calculate the result of calculation of the specified multiple of the maximum of described product, obtain the frequency domain maximum envelope that result of calculation maximum-frequency curve is tested as Y-direction;
Step S02 will obtain (Xcla-Z-Xtest)�@��(Ycla-Z-Ytest)�@(Zcla-Z-Ztest) three curves carry out envelope on frequency domain, find the maximum in described product on three curves of corresponding same described Frequency point respectively, and calculate the result of calculation of the specified multiple of the maximum of described product, obtain the frequency domain maximum envelope that result of calculation maximum-frequency curve is tested as Z-direction;
S04: the described frequency domain maximum envelope drawn is smoothed, obtain covering test on this correspondence direction of target unit of whole described frequency domain maximum envelope sine sweep test condition.
Preferably, in step S01, described appointment magnitude is selected in 0.05g.
Preferably, in step S04, adopt limited step line segment that described frequency domain maximum envelope is smoothed, provide the sinusoidal experiments condition under this same test direction of target unit.
Preferably, described specified multiple is 1.25 times.
Preferably, in step S02, during calculating, after in Stage Vehicle System level sine sweep test, acceleration transmission in target unit place is multiplied by specified multiple than amplitude, it is interpolated by the described Frequency point provided in step S01, Upper Stage in satellite and the rocket load coupling analysis result/rocket interface acceleration shock response spectrum maximum is interpolated by the described Frequency point provided in S01, then is multiplied by the above-mentioned two interpolated data on same described Frequency point.
Preferably, in step S02, during calculating, acceleration transmission in target unit place in Stage Vehicle System level sine sweep test is interpolated by the described Frequency point provided in step S01 than amplitude, after Upper Stage in satellite and the rocket load coupling analysis result/rocket interface acceleration shock response spectrum maximum is multiplied by specified multiple, it is interpolated by the described Frequency point provided in S01, then the above-mentioned two interpolated data on same described Frequency point is multiplied.
Preferably, in step S01, being spaced between 0.01Hz��0.2Hz of the described Frequency point provided.
Preferably, in step S01, the described Frequency point provided be spaced apart 0.05Hz.
Beneficial effects of the present invention is as follows:
The Upper Stage unit sine sweep test condition minute design method of the present invention, compare existing Upper Stage dynamic environmental test condition, the Upper Stage sine surface sweeping experimental condition provided is more fine, meet experimental condition to its actual experimental environment spreadability while, it is more reasonable that test magnitude and surplus select, reduce stand-alone environment adaptability design difficulty, reduce unit and cross the risk of experiment. And system-level sinusoidal surface sweeping experimental data processing, the experimental condition envelope in method for designing provided by the invention all can adopt digitized batch processing, it is effectively increased efficiency and the fine degree of mechanical environment condition design.
Accompanying drawing explanation
Fig. 1 is the flow chart of the Upper Stage unit sine sweep test condition minute design method of the present invention;
Fig. 2 is certain unit Y-direction low frequency sinusoidal condition of scanning envelope comparison diagram of the present invention.
Detailed description of the invention
Below in conjunction with the drawings and specific embodiments, technical scheme and beneficial effect are further illustrated.
It needs to be determined that coordinate system when the method for the present invention uses, the selection of coordinate system does not affect the enforcement of the method for the present invention, concrete, it is possible to choose rocket coordinate system, other coordinate system can certainly be chosen, as long as coordinate system is unified in the design process of whole experimental condition.
The Upper Stage unit sine sweep test condition minute design method of the present invention, comprises the following steps:
S01: in the frequency range needed for the sine sweep test condition of the unit that sets objectives, provides the Frequency point of certain intervals;
Obtain Upper Stage/rocket interface X/Y/Z tri-directional acceleration shock response spectrum-frequency curve under different operating modes in satellite and the rocket load coupling analysis result, acceleration shock response spectrum-frequency curve under all operating modes on same direction is carried out maximum envelope, find the maximum in same Frequency point on same direction, corresponding, acceleration shock response spectrum under different operating modes respectively, obtain Upper Stage/rocket interface X/Y/Z tri-directional acceleration shock response spectrum maximum-frequency curve, be expressed as Xcla��YclaAnd ZclaSuch as, the Upper Stage that can obtain in satellite and the rocket load coupling analysis result under three operating modes/rocket interface X/Y/Z tri-directional acceleration shock response spectrum-frequency curve, for X-direction, corresponding each Frequency point has three acceleration shock response modal data, find the maximum in these three data, the maximum of corresponding each Frequency point is sequentially connected with and namely obtains X-direction acceleration shock response spectrum maximum-frequency curve Xcla, Y-direction, Z-direction are also such.
Under acquisition appointment magnitude, in Stage Vehicle System level sine sweep test, under X/Y/Z tri-test direction, amplitude-frequency curve is compared in the acceleration transmission in tri-directions of X/Y/Z at target unit place, namely obtain X to test the acceleration transmission in tri-directions of X/Y/Z at target unit place under direction and test, than amplitude-frequency curve, Y, the acceleration transmission that the acceleration transmission in tri-directions of X/Y/Z at target unit place under direction tests tri-directions of X/Y/Z at target unit place under direction than amplitude-frequency curve, Z and compare amplitude-frequency curve, concrete, represent each curve with following symbol:
X-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
X-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
X-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
Magnitude is specified to select between 0.05g��0.2g;
S02: the acceleration transmission under the X/Y/Z direction at the target unit place obtained with X-direction test respectively by the acceleration shock response spectrum maximum of X-direction respectively is multiplied than amplitude, the acceleration transmission that the acceleration shock response spectrum maximum of Y-direction is tested with Y-direction under the X/Y/Z direction at the target unit place that direction obtains respectively is multiplied than amplitude, acceleration transmission under the X/Y/Z direction at the target unit place that the acceleration shock response spectrum maximum of Z-direction obtains with Z-direction test respectively is multiplied than amplitude, corresponding each direction obtains three result of calculation-frequency curves, concrete:
Calculate the X obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Xtest��Y-Xtest��Z-XtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three result of calculations-frequency curve Xcla-X-Xtest��Xcla-Y-XtestAnd Xcla-Z-Xtest;
Calculate the Y obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Ytest��Y-Ytest��Z-YtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three result of calculations-frequency curve Ycla-X-Ytest��Ycla-Y-YtestAnd Ycla-Z-Ytest;
Calculate the Z obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Ztest��Y-Ztest��Z-ZtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three result of calculations-frequency curve Zcla-X-Ztest��Zcla-Y-ZtestAnd Zcla-Z-Ztest;
Specified multiple selects between 1.20 times��1.40 times, and with fixing in a unit sine sweep test condition design;
S03: X will be obtained in step S02cla-X-Xtest��Ycla-X-Ytest��Zcla-X-ZtestArticle three, curve carries out envelope on frequency domain, finds the maximum in described result of calculation on three curves of corresponding same described Frequency point respectively, obtains result of calculation maximum-frequency curve as the X frequency domain maximum envelope to test;
Step S02 will obtain Xcla-Y-Xtest��Ycla-Y-Ytest��Zcla-Y-ZtestArticle three, curve carries out envelope on frequency domain, finds the maximum in described result of calculation on three curves of corresponding same described Frequency point respectively, obtains the frequency domain maximum envelope that result of calculation maximum-frequency curve is tested as Y-direction;
Step S02 will obtain Xcla-Z-Xtest��Ycla-Z-Ytest��Zcla-Z-ZtestArticle three, curve carries out envelope on frequency domain, finds the maximum in described result of calculation on three curves of corresponding same described Frequency point respectively, obtains the frequency domain maximum envelope that result of calculation maximum-frequency curve is tested as Z-direction;
S04: be smoothed by the frequency domain maximum envelope drawn, obtains the sine sweep test condition covering under this correspondence direction of target unit of whole frequency domain maximum envelope.
The method for designing of the present invention, it is as forcing function using Upper Stage/rocket interface acceleration shock response spectrum in satellite and the rocket load coupling analysis result, using the transmission of Stage Vehicle System level sine sweep test unit place acceleration than amplitude as transmission characteristic, both are multiplied, and take the safety coefficient of specified multiple, such as 1.25 times, draw the maximum environment envelope of unit, owing to being multiplication relation, the opportunity that takes into of this specified multiple can be arranged according to convenience of calculation, thus step S02, step S03 can make following adjustment:
S02: calculate the X obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Xtest��Y-Xtest��Z-XtestOn acceleration transmission than the product of amplitude, obtain three products-frequency curve (Xcla-X-Xtest)�@��(Xcla-Y-Xtest) and (Xcla-Z-Xtest);
Calculate the Y obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Ytest��Y-Ytest��Z-YtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three products-frequency curve (Ycla-X-Yest)�@��(Ycla-Y-Ytest) and (Ycla-Z-Ytest);
Calculate the Z obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Ztest��Y-Ztest��Z-ZtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three products-frequency curve (Zcla-X-Ztest)�@��(Zcla-Y-Ztest) and (Zcla-Z-Ztest); Described specified multiple selects between 1.20 times��1.40 times, and with fixing in a unit sine sweep test condition design;
S03: (X will be obtained in step S02cla-X-Xtest)�@��(Ycla-X-Ytest)�@��(Zcla-X-Ztest) three curves carry out envelope on frequency domain, find the maximum in described product on three curves of corresponding same described Frequency point respectively, and calculate the result of calculation of the specified multiple of the maximum of described product, obtain result of calculation maximum-frequency curve as the X frequency domain maximum envelope to test;
Step S02 will obtain (Xcla-Y-Xtest)�@��(Ycla-Y-Ytest)�@(Zcla-Y-Ztest) three curves carry out envelope on frequency domain, find the maximum in described product on three curves of corresponding same described Frequency point respectively, and calculate the result of calculation of the specified multiple of the maximum of described product, obtain the frequency domain maximum envelope that result of calculation maximum-frequency curve is tested as Y-direction;
Step S02 will obtain (Xcla-Z-Xtest)�@��(Ycla-Z-Ytest)�@(Zcla-Z-Ztest) three curves carry out envelope on frequency domain, find the maximum in described product on three curves of corresponding same described Frequency point respectively, and calculate the result of calculation of the specified multiple of the maximum of described product, obtain the frequency domain maximum envelope that result of calculation maximum-frequency curve is tested as Z-direction.
Same reason, the calculating process in step S02 can adopt the following two kinds:
Can when calculating, after in Stage Vehicle System level sine sweep test, acceleration transmission in target unit place is multiplied by specified multiple than amplitude, it is interpolated by the Frequency point provided in step S01, Upper Stage in satellite and the rocket load coupling analysis result/rocket interface acceleration shock response spectrum maximum is interpolated by the Frequency point provided in S01, then is multiplied by the above-mentioned two interpolated data on same Frequency point;
Can also when calculating, acceleration transmission in target unit place in Stage Vehicle System level sine sweep test is interpolated by the Frequency point provided in step S01 than amplitude, after Upper Stage in satellite and the rocket load coupling analysis result/rocket interface acceleration shock response spectrum maximum is multiplied by specified multiple, it is interpolated by the Frequency point provided in S01, then the above-mentioned two interpolated data on same Frequency point is multiplied.
The selection of magnitude influences whether the height of result of calculation, and magnitude is turned down, it is impossible to complete excitation, magnitude is too big, and result of calculation is too big, and experimental condition can be made harsher, appointment magnitude is selected in 0.05g through research by inventor, it is possible to reduce the severity of unit sine surface sweeping experimental condition.
If the frequency domain maximum envelope obtained in step S03 can be made bad operation in test directly as experimental condition, therefore, it is smoothed by the present invention in step S04, specifically, as shown in Figure 2, can adopt limited step line segment that frequency domain maximum envelope is smoothed, provide the sinusoidal experiments condition under this same test direction of target unit. When arranging step line segment, it is possible to according to crest location and the height segmentation of frequency domain maximum envelope, each section according to the multiple of more than 1 times of peak value, for instance can be 1.25 times, it is determined that experimental condition.
The present inventor considers flight remote measurement and the relation of load coupling analysis, and the specified multiple in step S01 is set to 1.25 times, so that while avoiding the severity of environmental condition design, it is ensured that environmental condition design possesses certain surplus.
In the present invention, the interval between Frequency point in step S01 can determine according to the designing requirement of experimental condition, generally can be arranged between 0.01Hz��0.2Hz, it is such as 0.05Hz, if requiring higher can reduce interval further, certainly, it is desirable to relatively low words can also expand interval.
Fig. 2 gives the design result of the unit sine surface sweeping experimental condition designed for some unit distinct methods, wherein, the method that the curve being made up of step line segment is the present invention obtains, in remaining three curves, a bottom curve is to obtain in satellite and the rocket loading analysis, topmost one is adopt traditional obtaining according to system-level low-frequency vibration test data, middle one is that the method for the present invention is in the step S03 frequency domain maximum envelope obtained, as we can see from the figure, between result that the maximum envelope obtained according to the method for the present invention obtains at satellite and the rocket load coupling analysis and the result that tradition obtains according to the method for system-level low-frequency vibration experimental condition, thus being provided simultaneously with certain surplus what avoid Environment Design severity, reduce stand-alone environment adaptability design difficulty, reduce unit and cross the risk of experiment.
To sum up, the method for designing of the present invention, when level unit sine sweep test condition is formulated herein above, adopt theoretical and test combination method, using Upper Stage/rocket interface acceleration in satellite and the rocket load coupling analysis result as forcing function, Stage Vehicle System level sine sweep test unit place acceleration responsive with excitation ratio as transmission characteristic, both are multiplied, and take the safety coefficient of specified multiple, draw the maximum environment envelope of unit. The method can provide the environment envelope that different frequency place is corresponding, and experimental condition is more fine, thus avoiding the severity of environmental condition design, experimental condition relatively accurately and has certain surplus simultaneously. And in the present invention, system-level sine sweep test data process, experimental condition envelope all can adopt datumization batch processing, is effectively increased efficiency and the fine degree of mechanical environment condition design.
Claims (9)
1. Upper Stage unit sine sweep test condition minute design method, it is characterised in that comprise the following steps:
S01: in the frequency range needed for the sine sweep test condition of the unit that sets objectives, provides the Frequency point of certain intervals;
Obtain Upper Stage/rocket interface X/Y/Z tri-directional acceleration shock response spectrum-frequency curve under different operating modes in satellite and the rocket load coupling analysis result, acceleration shock response spectrum-frequency curve under all operating modes on same direction is carried out maximum envelope, find the maximum in the acceleration shock response spectrum under operating modes on same direction, the same described Frequency point of correspondence, different respectively, obtain Upper Stage/rocket interface X/Y/Z tri-directional acceleration shock response spectrum maximum-frequency curve, be expressed as Xcla��YclaAnd Zcla;
Under acquisition appointment magnitude, in Stage Vehicle System level sine sweep test, X/Y/Z tri-tests the acceleration transmission ratio amplitude-frequency curve in tri-directions of X/Y/Z at target unit place under direction, is expressed as:
X-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
X-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
X-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
Described appointment magnitude selects between 0.05g��0.2g;
S02: calculate the X obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Xtest��Y-Xtest��Z-XtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three result of calculations-frequency curve Xcla-X-Xtest��Xcla-Y-XtestAnd Xcla-Z-Xtest;
Calculate the Y obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Ytest��Y-Ytest��Z-YtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three result of calculations-frequency curve Ycla-X-Ytest��Ycla-Y-YtestAnd Ycla-Z-Ytest;
Calculate the Z obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Ztest��Y-Ztest��Z-ZtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three result of calculations-frequency curve Zcla-X-Ztest��Zcla-Y-ZtestAnd Zcla-Z-Ztest;
Described specified multiple selects between 1.20 times��1.40 times, and with fixing in a unit sine sweep test condition design;
S03: X will be obtained in step S02cla-X-Xtest��Ycla-X-Ytest��Zcla-X-ZtestArticle three, curve carries out envelope on frequency domain, finds the maximum in described result of calculation on three curves of corresponding same described Frequency point respectively, obtains result of calculation maximum-frequency curve as the X frequency domain maximum envelope to test;
Step S02 will obtain Xcla-Y-Xtest��Ycla-Y-Ytest��Zcla-Y-ZtestArticle three, curve carries out envelope on frequency domain, finds the maximum in described result of calculation on three curves of corresponding same described Frequency point respectively, obtains the frequency domain maximum envelope that result of calculation maximum-frequency curve is tested as Y-direction;
Step S02 will obtain Xcla-Z-Xtest��Ycla-Z-Ytest��Zcla-Z-ZtestArticle three, curve carries out envelope on frequency domain, finds the maximum in described result of calculation on three curves of corresponding same described Frequency point respectively, obtains the frequency domain maximum envelope that result of calculation maximum-frequency curve is tested as Z-direction;
S04: be smoothed by the described frequency domain maximum envelope drawn, obtains the sine sweep test condition covering under this correspondence direction of target unit of whole described frequency domain maximum envelope.
2. Upper Stage unit sine sweep test condition minute design method, it is characterised in that comprise the following steps:
S01: in the frequency range needed for the unit sine sweep test condition that sets objectives, provide the Frequency point of certain intervals;
Obtain Upper Stage/rocket interface X/Y/Z tri-directional acceleration shock response spectrum-frequency curve under different operating modes in satellite and the rocket load coupling analysis result, acceleration shock response spectrum-frequency curve under all operating modes on same direction is carried out maximum envelope, find the maximum in the acceleration shock response spectrum under operating modes on same direction, the same described Frequency point of correspondence, different respectively, obtain Upper Stage/rocket interface X/Y/Z tri-directional acceleration shock response spectrum maximum-frequency curve, be expressed as Xcla��YclaAnd Zcla;
Under acquisition appointment magnitude, in Stage Vehicle System level sine sweep test, X/Y/Z tri-tests the acceleration transmission ratio amplitude-frequency curve in tri-directions of X/Y/Z at target unit place under direction, is expressed as:
X-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Xtest, represent that under X test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
X-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Ytest, represent that under Y test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
X-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the X-direction at target unit place;
Y-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the Y-direction at target unit place;
Z-Ztest, represent that under Z test direction, amplitude-frequency curve is compared in the acceleration transmission of the Z-direction at target unit place;
Described appointment magnitude selects between 0.05g��0.2g;
S02: calculate the X obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Xtest��Y-Xtest��Z-XtestOn acceleration transmission than the product of amplitude, obtain three products-frequency curve (Xcla-X-Xtest)�@��(Xcla-Y-Xtest) and (Xcla-Z-Xtest);
Calculate the Y obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Ytest��Y-Ytest��Z-YtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three products-frequency curve (Ycla-X-Yest)�@��(Ycla-Y-Ytest) and (Ycla-Z-Ytest);
Calculate the Z obtaining on same described Frequency point respectivelyclaOn acceleration shock response spectrum maximum and X-Ztest��Y-Ztest��Z-ZtestOn acceleration transmission than the product of amplitude, and the specified multiple calculating this product obtains result of calculation, obtains three products-frequency curve (Zcla-X-Ztest)�@��(Zcla-Y-Ztest) and (Zcla-Z-Ztest); Described specified multiple selects between 1.20 times��1.40 times, and with fixing in a unit sine sweep test condition design;
S03: (X will be obtained in step S02cla-X-Xtest)�@��(Ycla-X-Ytest)�@��(Zcla-X-Ztest) three curves carry out envelope on frequency domain, find the maximum in described product on three curves of corresponding same described Frequency point respectively, and calculate the result of calculation of the specified multiple of the maximum of described product, obtain result of calculation maximum-frequency curve as the X frequency domain maximum envelope to test;
Step S02 will obtain (Xcla-Y-Xtest)�@��(Ycla-Y-Ytest)�@(Zcla-Y-Ztest) three curves carry out envelope on frequency domain, find the maximum in described product on three curves of corresponding same described Frequency point respectively, and calculate the result of calculation of the specified multiple of the maximum of described product, obtain the frequency domain maximum envelope that result of calculation maximum-frequency curve is tested as Y-direction;
Step S02 will obtain (Xcla-Z-Xtest)�@��(Ycla-Z-Ytest)�@(Zcla-Z-Ztest) three curves carry out envelope on frequency domain, find the maximum in described product on three curves of corresponding same described Frequency point respectively, and calculate the result of calculation of the specified multiple of the maximum of described product, obtain the frequency domain maximum envelope that result of calculation maximum-frequency curve is tested as Z-direction;
S04: the described frequency domain maximum envelope drawn is smoothed, obtain covering test on this correspondence direction of target unit of whole described frequency domain maximum envelope sine sweep test condition.
3. Upper Stage unit sine sweep test condition minute design method as claimed in claim 1 or 2, it is characterised in that: in step S01, described appointment magnitude is selected in 0.05g.
4. Upper Stage unit sine sweep test condition minute design method as claimed in claim 3, it is characterized in that: in step S04, adopt limited step line segment that described frequency domain maximum envelope is smoothed, provide the sinusoidal experiments condition under this same test direction of target unit.
5. Upper Stage unit sine sweep test condition minute design method as claimed in claim 4, it is characterised in that: described specified multiple is 1.25 times.
6. Upper Stage unit sine sweep test condition minute design method as claimed in claim 5, it is characterized in that: in step S02, during calculating, after in Stage Vehicle System level sine sweep test, acceleration transmission in target unit place is multiplied by specified multiple than amplitude, it is interpolated by the described Frequency point provided in step S01, Upper Stage in satellite and the rocket load coupling analysis result/rocket interface acceleration shock response spectrum maximum is interpolated by the described Frequency point provided in S01, then is multiplied by the above-mentioned two interpolated data on same described Frequency point.
7. Upper Stage unit sine sweep test condition minute design method as claimed in claim 5, it is characterized in that: in step S02, during calculating, acceleration transmission in target unit place in Stage Vehicle System level sine sweep test is interpolated by the described Frequency point provided in step S01 than amplitude, after Upper Stage in satellite and the rocket load coupling analysis result/rocket interface acceleration shock response spectrum maximum is multiplied by specified multiple, it is interpolated by the described Frequency point provided in S01, then the above-mentioned two interpolated data on same described Frequency point is multiplied.
8. Upper Stage unit sine sweep test condition minute design method as claimed in claim 7, it is characterised in that: in step S01, being spaced between 0.01Hz��0.2Hz of the described Frequency point provided.
9. as claimed in claim 8 Upper Stage unit sine sweep test condition minute design method, it is characterised in that: in step S01, the described Frequency point provided be spaced apart 0.05Hz.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410637581.XA CN105628416B (en) | 2014-11-06 | 2014-11-06 | Upper Stage unit sine sweep test condition minute design method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410637581.XA CN105628416B (en) | 2014-11-06 | 2014-11-06 | Upper Stage unit sine sweep test condition minute design method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105628416A true CN105628416A (en) | 2016-06-01 |
| CN105628416B CN105628416B (en) | 2017-12-22 |
Family
ID=56043556
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410637581.XA Active CN105628416B (en) | 2014-11-06 | 2014-11-06 | Upper Stage unit sine sweep test condition minute design method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105628416B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109946030A (en) * | 2019-03-11 | 2019-06-28 | 北京星际荣耀空间科技有限公司 | A kind of carrier rocket vibration data detection method, device, system and storage medium |
| CN111707528A (en) * | 2020-07-16 | 2020-09-25 | 上海宇航系统工程研究所 | Dynamic environment time-frequency conversion method based on cumulative damage equivalence |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070012821A1 (en) * | 2004-08-11 | 2007-01-18 | Buehler David B | Launch vehicle crew escape system |
| CN101354251A (en) * | 2008-09-12 | 2009-01-28 | 航天东方红卫星有限公司 | Method for determining deep space detector equivalent transfer orbit |
| RU2422335C1 (en) * | 2009-12-21 | 2011-06-27 | Открытое акционерное общество "Ракетно-космическая корпорация "Энергия" имени С.П. Королева" | Spacehead |
| US20120227374A1 (en) * | 2011-03-09 | 2012-09-13 | United Launch Alliance, Llc | Integrated vehicle fluids |
| CN103885456A (en) * | 2014-03-31 | 2014-06-25 | 航天东方红卫星有限公司 | Method for distributing and transferring multiple satellites launched by one rocket to different orbit planes |
| CN105659841B (en) * | 2012-03-30 | 2014-07-09 | 北京空间飞行器总体设计部 | A kind of aeronautical satellite thermal vacuum test method |
-
2014
- 2014-11-06 CN CN201410637581.XA patent/CN105628416B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20070012821A1 (en) * | 2004-08-11 | 2007-01-18 | Buehler David B | Launch vehicle crew escape system |
| CN101354251A (en) * | 2008-09-12 | 2009-01-28 | 航天东方红卫星有限公司 | Method for determining deep space detector equivalent transfer orbit |
| RU2422335C1 (en) * | 2009-12-21 | 2011-06-27 | Открытое акционерное общество "Ракетно-космическая корпорация "Энергия" имени С.П. Королева" | Spacehead |
| US20120227374A1 (en) * | 2011-03-09 | 2012-09-13 | United Launch Alliance, Llc | Integrated vehicle fluids |
| CN105659841B (en) * | 2012-03-30 | 2014-07-09 | 北京空间飞行器总体设计部 | A kind of aeronautical satellite thermal vacuum test method |
| CN103885456A (en) * | 2014-03-31 | 2014-06-25 | 航天东方红卫星有限公司 | Method for distributing and transferring multiple satellites launched by one rocket to different orbit planes |
Non-Patent Citations (2)
| Title |
|---|
| 王建民 等: "上面级结构模态试验与分析技术", 《航天器环境工程》 * |
| 胡伟峰 等: "先进上面级动特性分析及结构动力学改进设计", 《导弹与航天运载技术》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109946030A (en) * | 2019-03-11 | 2019-06-28 | 北京星际荣耀空间科技有限公司 | A kind of carrier rocket vibration data detection method, device, system and storage medium |
| CN111707528A (en) * | 2020-07-16 | 2020-09-25 | 上海宇航系统工程研究所 | Dynamic environment time-frequency conversion method based on cumulative damage equivalence |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105628416B (en) | 2017-12-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fujii et al. | Modelling the Milky Way as a dry Galaxy | |
| CN105912772B (en) | A kind of model modification method integrating frequency, the vibration shape and frequency response function | |
| CN106055486A (en) | Automatic operation maintenance method and platform of distributed file system | |
| CN103616149B (en) | Whole star micro-vibration ground experiment free boundary analog supports frequency test method | |
| CN105628416A (en) | Upper stage single machine sine scan test condition refining design method | |
| CN102810127A (en) | Virtual vibration test system for spacecraft | |
| CN106777478A (en) | The method of vehicle cantilever Strength Analysis for Structural Parts | |
| CN102567578A (en) | Spacecraft vibration test fixture evaluation system | |
| CN105510891B (en) | For the in-orbit modal identification method of satellites with synthetic aperture radar | |
| CN108090263B (en) | Verification method for pneumatic damping of deicing vibration scale test of power transmission line | |
| CN112836351B (en) | Design method and system for sine vibration test of commercial microsatellite satellite-mounted product | |
| CN110455477A (en) | A kind of acquisition methods of solid-rocket cargo tank structure oscillating load spectrum | |
| CN106168924A (en) | Software testing method and device and electronic equipment | |
| CN107727350B (en) | Micro-nano satellite vector vibration test method | |
| CN105953816A (en) | Distributed atmosphere data system | |
| CN102270249A (en) | Method for identifying characteristic frequency of parts | |
| CN108981864A (en) | A kind of liquid level calibrating installation and calibration method | |
| CN105548921B (en) | A kind of method and device of the determining bearing containing magnetic | |
| CN103759951A (en) | Tester capable of simultaneously testing brake performance of multiple agricultural vehicles and test method | |
| CN109944581A (en) | The acquisition methods and device of oil well set damage | |
| CN1229647C (en) | Automatic positioning system/global positioning system receiver comprehensive testing system | |
| CN208238889U (en) | A kind of tank level calibrating installation | |
| Steffen et al. | The sensitivity of GNSS measurements in Fennoscandia to distinct three-dimensional upper-mantle structures | |
| CN103185563A (en) | Angle calculating device | |
| CN105570377A (en) | Micro-vibration passive vibration isolation device for control moment gyro |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |