CN113094816B - Method for constructing comprehensive working condition vibration spectrum and long-life test spectrum of armored vehicle - Google Patents
Method for constructing comprehensive working condition vibration spectrum and long-life test spectrum of armored vehicle Download PDFInfo
- Publication number
- CN113094816B CN113094816B CN202110405031.5A CN202110405031A CN113094816B CN 113094816 B CN113094816 B CN 113094816B CN 202110405031 A CN202110405031 A CN 202110405031A CN 113094816 B CN113094816 B CN 113094816B
- Authority
- CN
- China
- Prior art keywords
- spectrum
- working condition
- coefficient
- reserve
- comprehensive
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computer Hardware Design (AREA)
- Data Mining & Analysis (AREA)
- Software Systems (AREA)
- Mathematical Physics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Probability & Statistics with Applications (AREA)
- Operations Research (AREA)
- Evolutionary Biology (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The invention discloses a method for constructing a comprehensive working condition vibration spectrum and a long service life test spectrum of an armored vehicle, which relates to the field of data processing methods, wherein the method for constructing the comprehensive working condition vibration spectrum comprises the steps of obtaining a first relation and a second relation among a single working condition actual measurement spectrum, a first storage coefficient and the comprehensive working condition actual measurement spectrum; obtaining the relation between the first reserve coefficient and the measured spectrum of the comprehensive working condition according to the first relation and the second relation; training a first RBF neural network to obtain a current first reserve coefficient; obtaining a first relation and a second relation among the single working condition standard spectrum, the second reserve coefficient and the comprehensive working condition standard spectrum; obtaining the relation between the second reserve coefficient and the comprehensive working condition standard spectrum according to the first relation and the second relation; training a second RBF neural network to obtain a current second reserve coefficient; and obtaining a reserve coefficient according to the current first and second reserve coefficients, thereby obtaining an actual measured spectrum and a standard spectrum of the actual comprehensive working condition. The invention improves the accuracy and can simulate the real vibration environment more effectively.
Description
Technical Field
The invention relates to the technical field of armored vehicle vibration spectrum data processing methods, in particular to a method for constructing an armored vehicle comprehensive working condition vibration spectrum and a long service life test spectrum.
Background
The armored vehicle belongs to military vehicles, can run on various road surfaces, is mainly used for transportation, investigation, maneuverability attack and the like, and therefore the vibration of the armored vehicle not only affects the fatigue degree and the working efficiency of fighters, but also affects the attack precision of weapons arranged on the vehicle. If the whole vehicle vibrates excessively, the service life of the device can be shortened, and even the device is damaged.
In order to analyze the vibration performance of armored vehicles, a test method is mainly adopted at present to test the structural strength, durability and the like of armored vehicle frames and the like. For this reason, long-life simulation tests are often used for quality inspection to determine the vibration performance and safety and reliability of the entire armored vehicle.
However, the random vibration environment measurement data induction method specified in the GJB/Z126-99 vibration and impact environment test data induction method can only perform data induction on random vibration under single working condition to obtain a random vibration standard spectrum. And the armored vehicle can have different road surfaces and different vehicle speeds and other comprehensive working conditions when actually running, so the method aiming at the single working condition can not be directly used for the comprehensive working conditions, and can not be directly used for the long-service-life simulation test.
Disclosure of Invention
Therefore, in order to overcome the defects, the embodiment of the invention provides a method for constructing a comprehensive working condition vibration spectrum and a long-life test spectrum of an armored vehicle, which can realize reconstruction of a comprehensive working condition vibration actual measurement spectrum and a standard spectrum, and construct the long-life test spectrum suitable for the comprehensive working condition based on the comprehensive working condition vibration standard spectrum.
Therefore, the method for constructing the comprehensive working condition vibration spectrum of the armored vehicle comprises the following steps:
obtaining the road surface occupation ratio b of the armored vehicle under the ith road condition under the test comprehensive working conditioniAnd the vehicle speed ratio a of the jth vehicle speed used for driving on the ith road conditionijWherein, in the step (A),1, p is the total number of road surface types, and q is the total number of vehicle speed types;
respectively obtaining the single-working-condition actual measurement spectrum of the mth sample under the road condition i and the vehicle speed jFrom the measured spectrum of the simplex conditionFirst reserve factor cij (1)Inducing to obtain comprehensive working condition actual measurement spectrum FGS of all vehicle speeds and road conditions of the measuring points, and obtaining single working condition actual measurement spectrumFirst reserve factor cij (1)And a first relation between the measured spectrum FGS and the comprehensive working condition;
from the measured spectrum of the simplex conditionAnd a first reserve coefficient cij (1)Calculating to obtain comprehensive working condition actual measurement spectrum FGS of all vehicle speeds and road conditions of the measuring points, and obtaining single working condition actual measurement spectrumFirst reserve factor cij (1)And a second relation between the measured spectrum FGS and the comprehensive working condition;
single working condition actual measurement spectrumFirst reserve factor cij (1)Substituting a second relation between the measured spectrum and the comprehensive working condition actual measurement spectrum FGS into the single working condition actual measurement spectrumFirst reserve factor cij (1)Obtaining a first storage coefficient c according to a first relation between the measured spectrum and the comprehensive working condition FGSij (1)And the relation between the measured spectrum FGS and the comprehensive working condition;
the measured spectrum FGS of the comprehensive working condition is used as the input of the first RBF neural network and is based on the first reserve coefficient cij (1)A first reserve coefficient c obtained by the relation between the measured spectrum FGS and the comprehensive working conditionij (1)As expected output of the first RBF neural network, training the first RBF neural network by utilizing input and expected output, and inputting the actually measured spectrum of the current comprehensive working condition into the trained first RBF neural network to obtain a current first reserve coefficient;
respectively obtaining the standard spectrum of the single working condition of the mth sample under the conditions of the road condition i and the vehicle speed jAccording to a single working condition standard spectrumSecond reserve factor cij (2)Obtaining the comprehensive working condition standard spectrum FGG of all vehicle speeds and road conditions of the measuring points by induction, and obtaining the single working condition standard spectrumSecond reserve factor cij (2)And a first relation between the combined working condition normalized spectrum FGG;
according to a single working condition standard spectrumAnd a second reserve coefficient cij (2)Calculating to obtain comprehensive condition standard spectrum FGG of all vehicle speeds and road conditions of the measuring points, and obtaining single condition standard spectrumSecond reserve factor cij (2)And a second relationship between the combined operating mode normalized spectrum FGG;
standard spectrum of single working conditionSecond reserve factor cij (2)And the second relation generation between the comprehensive working condition specification spectrum FGGEntry condition standard spectrumSecond reserve factor cij (2)And obtaining a second reserve coefficient c according to a first relation between the comprehensive working condition normalized spectrum FGGij (2)And the relation between the comprehensive working condition standard spectrum FGG;
taking the comprehensive working condition standard spectrum FGG as the input of a second RBF neural network, and according to a second reserve coefficient cij (2)Second reserve coefficient c obtained by relation with comprehensive working condition standard spectrum FGGij (2)As expected output of the second RBF neural network, training the second RBF neural network by utilizing input and expected output, and inputting the current comprehensive working condition standard spectrum into the trained second RBF neural network to obtain a current second reserve coefficient;
calculating to obtain a reserve coefficient c according to the current first reserve coefficient and the current second reserve coefficientij;
According to the reserve coefficient cijActual measurement spectrum G of single working condition under road condition i and vehicle speed jshice(i, j) and the simplex condition normalized spectrum Gguifan(i, j), respectively calculating and obtaining comprehensive working condition actual measurement spectrums FGS and comprehensive working condition standard spectrums FGG of all vehicle speeds and road conditions of the measured points.
Preferably, the measured spectrum is based on a single-working conditionAnd a first reserve coefficient cij (1)The step of inducing to obtain the comprehensive working condition actual measurement spectrum FGS of all the vehicle speeds and road conditions of the measuring points comprises the following steps:
from the measured spectrum of the simplex conditionAnd a first reserve coefficient cij (1)Calculating to obtain the mean value and the variance under the road condition i and the vehicle speed j, wherein the calculation formula is as follows:
wherein N isi,jThe number of samples under a single working condition;
calculating to obtain weighted mean value and variance according to the mean value and variance under the road condition i and the vehicle speed j, wherein the calculation formula is as follows:
and calculating to obtain the comprehensive working condition actual measurement spectrum FGS of all vehicle speeds and road conditions of the measuring points according to the weighted mean value and the variance, wherein the calculation formula is as follows:
wherein, F11The first tolerance upper limit coefficient is calculated by the following formula:
wherein Z isβTo satisfy Prob [ Z.ltoreq.Zβ]A normally distributed quantile of beta,is a degree of freedom of (N-1) ×2Distribution of alpha quantites, tN-1;(1-α)1-alpha quantile distributed at the center t with the degree of freedom (N-1), N is the total number of samples,
preferably, the measured spectrum is based on a single-working conditionAnd a first reserve coefficient cij (1)The calculation formula of the comprehensive working condition actual measurement spectrum FGS for calculating and obtaining all vehicle speeds and road conditions of the measuring points is as follows:
wherein N isi,jIs the number of samples in the simplex case.
Preferably, the normalized spectrum according to single-working-conditionAnd a second reserve coefficient cij (2)The step of summarizing and obtaining the comprehensive working condition standard spectrum FGG of all vehicle speeds and road conditions of the measured points comprises the following steps:
according to a single working condition standard spectrumAnd a second reserve coefficient cij (2)Calculating to obtain the mean value and the variance of each spectral line under the road condition i and the vehicle speed j, wherein the calculation formula is as follows:
wherein N isi,jThe number of samples under a single working condition;
calculating to obtain a weighted mean value and a weighted variance according to the mean value and the variance of each spectral line under the road condition i and the vehicle speed j, wherein the calculation formula is as follows:
obtaining a statistic F by weighted mean and variance calculationnAnd tnThe calculation formula is as follows:
according to statistic FnAnd tnFrequency division is carried out under the confidence coefficient (1-alpha) to obtain H1A frequency band, wherein the spectral line number of two end points of the h-th frequency band is kh1、kh2,h=1,2,…,H1The number of spectral lines is Nh=kh2-kh1+1;
And calculating to obtain a weighted mean value and a weighted variance in the h frequency band according to the weighted mean value and the weighted variance, wherein the calculation formula is as follows:
calculating according to the weighted mean and the variance in the h frequency band to obtain the straight spectrum tolerance upper limit estimation in the h frequency band, wherein the calculation formula is as follows:
wherein, F12The second tolerance upper limit coefficient with the confidence coefficient of (1-alpha) and the quantile point of beta is calculated by the following formula:
wherein Z isβTo satisfy Prob [ Z.ltoreq.Zβ]A normally distributed quantile of beta,with a degree of freedom of (N.N)hX of-1)2The distribution of the alpha quantile points is,with a degree of freedom of (N.N)h-1) central t distribution 1-a quantiles;
and taking the upper limit estimation of the tolerance of the straight spectrum in the h-th frequency band as the amplitude of the frequency band, and connecting each adjacent straight spectrum by using a straight line under a log-log coordinate to obtain a comprehensive working condition standard spectrum FGG of all vehicle speeds and road conditions of the measured point.
Preferably, the normalized spectrum according to single-working-conditionAnd a second reserve coefficient cij (2)The calculation formula of the comprehensive working condition standard spectrum FGG for calculating and obtaining all vehicle speeds and road conditions of the measured points is as follows:
wherein N isi,jIs the number of samples in the simplex case.
Preferably, the reserve coefficient cijThe rms values of the current first reserve coefficient and the current second reserve coefficient.
Preferably, said reserve coefficient c is used for determining a reserve valueijActual measurement spectrum G of single working condition under road condition i and vehicle speed jshice(i, j) and the simplex condition normalized spectrum Gguifan(i, j), respectively calculating and obtaining the calculation formulas of the comprehensive working condition actual measurement spectrum FGS and the comprehensive working condition standard spectrum FGG of all vehicle speeds and road conditions of the measured points, wherein the calculation formulas are respectively as follows:
the method for constructing the long-life test spectrum of the armored vehicle comprises the following steps:
obtaining the time T of the jth vehicle speed used for driving under the ith road conditionijI is 1,2, …, p, j is 1,2, …, q, p is the total number of road surface types, q is the total number of vehicle speed types;
time T of jth vehicle speed used for driving according to ith road conditionijAnd calculating to obtain a segment acceleration proportional coefficient beta under the road condition i and the vehicle speed jij;
Constructing an armored vehicle according to the aboveComprehensive working condition standard spectrum obtained by method of vehicle comprehensive working condition vibration spectrumAnd a segment acceleration proportionality coefficient betaijAnd calculating to obtain a long-life test spectrum LF.
Preferably, the piecewise acceleration proportionality coefficient betaijThe calculation formula of (2) is as follows:
wherein, TsFor the actual run time, m is an index.
Preferably, the long-life test spectrum LF is calculated by the formula:
the technical scheme of the embodiment of the invention has the following advantages:
1. according to the method for constructing the armored vehicle comprehensive working condition vibration spectrum, the measured spectrum and the standard spectrum of the comprehensive working condition are respectively deduced by using the induction method and the formula calculation method, the relationship between the first storage coefficient and the measured spectrum of the comprehensive working condition and the relationship between the second storage coefficient and the standard spectrum of the comprehensive working condition are respectively obtained by using the conclusions of the two methods, so that the RBF neural network is respectively trained based on the relationships to obtain the first storage coefficient estimation value and the second storage coefficient estimation value, and then the storage coefficient adopted when the vibration spectrum is actually constructed is calculated according to the estimation values, so that the construction precision is improved.
2. According to the method for constructing the long-life test spectrum of the armored vehicle, provided by the embodiment of the invention, the acceleration proportional coefficient is segmented according to the road condition and the vehicle speed, so that the method is more consistent with the actual vibration data measurement test condition, the accuracy is improved, the real vibration environment can be more effectively simulated, and the test result is more real and reliable.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a flow chart of a specific example of a method for constructing an armored vehicle combination behavior vibration spectrum according to embodiment 1 of the present invention;
FIG. 2 is a flow chart of another specific example of a method for constructing an armored vehicle combination vibration spectrum according to embodiment 1 of the present invention;
FIG. 3 is a flow chart of yet another specific example of a method for constructing a comprehensive working condition vibration spectrum of an armored vehicle in embodiment 1 of the present invention;
fig. 4 is a flowchart of a specific example of the method of constructing the long life test spectrum of the armored vehicle in embodiment 2 of the present invention.
Detailed Description
The technical solutions of the present invention will be described clearly and completely with reference to the accompanying drawings, and it is to be understood that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
In describing the present invention, it is to be understood that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. The terms "comprises" and/or "comprising," when used in this specification, are intended to specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof. The term "and/or" includes any and all combinations of one or more of the associated listed items.
Furthermore, certain drawings in this specification are flow charts illustrating methods. It will be understood that each block of the flowchart illustrations, and combinations of blocks in the flowchart illustrations, can be implemented by computer program instructions. These computer program instructions may be loaded onto a computer or other programmable apparatus to produce a machine, such that the instructions which execute on the computer or other programmable apparatus create means for implementing the functions specified in the flowchart block or blocks. These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart block or blocks. The computer program instructions may also be loaded onto a computer or other programmable apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart block or blocks.
Accordingly, blocks of the flowchart illustrations support combinations of means for performing the specified functions and combinations of steps for performing the specified functions. It will also be understood that each block of the flowchart illustrations, and combinations of blocks in the flowchart illustrations, can be implemented by special purpose hardware-based computer systems which perform the specified functions or steps, or combinations of special purpose hardware and computer instructions.
In addition, the technical features involved in the different embodiments of the present invention described below may be combined with each other as long as they do not conflict with each other.
Example 1
The embodiment provides a method for constructing a vibration spectrum of comprehensive working conditions of an armored vehicle, wherein the vibration spectrum comprises an actually measured spectrum and a standard spectrum, and as shown in figure 1, the method comprises the following steps:
102-1, actually measuring the spectrum according to the simplex conditionAnd a first reserve coefficient cij (1)Calculating to obtain the mean value and the variance under the road condition i and the vehicle speed j, wherein the calculation formula is as follows:
wherein N isi,jThe number of samples under a single working condition;
102-2, calculating to obtain a weighted mean value and a weighted variance according to the mean value and the variance under the road condition i and the vehicle speed j, wherein the calculation formula is as follows:
102-3, calculating and obtaining comprehensive working condition actual measurement spectrums FGS of all vehicle speeds and road conditions of the measured points according to the weighted mean value and the variance, wherein the calculation formula is as follows:
wherein, F11The first tolerance upper limit coefficient is calculated by the following formula:
wherein Z isβTo satisfy Prob [ Z ≦ Zβ]A normally distributed quantile of beta,is a degree of freedom of (N-1) ×2Distribution of alpha quantites, tN-1;(1-α)1-alpha quantile distributed at the center t with the degree of freedom (N-1), N is the total number of samples,obtaining single working condition actual measurement spectrumFirst reserve factor cij (1)And a first relation between the measured spectrum FGS of the comprehensive working condition.
103, actually measuring the spectrum according to the single working conditionAnd a first reserve coefficient cij (1)Calculating to obtain comprehensive working condition actual measurement spectrum FGS of all vehicle speeds and road conditions of the measuring points, and obtaining single working condition actual measurement spectrumFirst reserve factor cij (1)And the second relation between the measured spectrum FGS under the comprehensive working condition, the calculation formula is as follows:
104, single working condition actual measurement spectrumFirst reserve coefficient cij (1)Substituting a second relation between the measured spectrum and the comprehensive working condition actual measurement spectrum FGS into the single working condition actual measurement spectrumFirst reserve factor cij (1)And obtaining a first reserve coefficient c by a first relation between the measured spectrum FGS and the comprehensive working conditionij (1)And the relation between the measured spectrum FGS under the comprehensive working condition.
105, taking the measured spectrum FGS of the comprehensive working condition as the input of the first RBF neural network, and according to the first reserve coefficient cij (1)A first reserve coefficient c obtained by the relation between the measured spectrum FGS and the comprehensive working conditionij (1)And as the expected output of the first RBF neural network, training the first RBF neural network by utilizing the input and the expected output, and inputting the actually measured spectrum of the current comprehensive working condition into the trained first RBF neural network to obtain a current first reserve coefficient.
106, respectively obtaining the standard spectrum of the single working condition of the mth sample under the road condition i and the vehicle speed jAccording to the sheetOperating mode standard spectrumAnd a second reserve coefficient cij (2)Obtaining the comprehensive working condition standard spectrum FGG of all vehicle speeds and road conditions of the measuring points by induction, and obtaining the single working condition standard spectrumSecond reserve factor cij (2)And the behaviour-by-wire normalization spectrum FGG, as shown in fig. 3, the generalization procedure is as follows:
step 106-1, normalizing the spectrum according to the single working conditionAnd a second reserve coefficient cij (2)Calculating to obtain the mean value and the variance of each spectral line under the road condition i and the vehicle speed j, wherein the calculation formula is as follows:
wherein N isi,jThe number of samples under a single working condition;
106-2, calculating to obtain a weighted mean value and a weighted variance according to the mean value and the variance of each spectral line under the road condition i and the vehicle speed j, wherein the calculation formula is as follows:
step 106-3, calculating to obtain statistic F according to weighted mean and variancenAnd tnThe calculation formula is as follows:
step 106-4, according to the statistic FnAnd tnFrequency division is carried out under the confidence coefficient (1-alpha) to obtain H1A frequency band, wherein the spectral line number of two end points of the h-th frequency band is kh1、kh2,h=1,2,…,H1The number of spectral lines is Nh=kh2-kh1+1;
Preferably, step 106-4 comprises:
judging whether the frequency bands are the same under a given confidence coefficient (1-alpha), wherein the judgment conditions are as follows:
F(N-1,N-1);α/2≤Fn(k,k+1)≤F(N-1,N-1);(1-α/2)
|tn(k,k+1)|≤t2(N-1);(1-α/2)
wherein, t2(N-1);(1-α/2)The central t with the degree of freedom of 2(N-1) is distributed with 1-alpha/2 quantile points, F(N-1,N-1);α/2F with degree of freedom of (N-1) is distributed at alpha/2 quantile, if Fn(k, k +1) obeys F distribution with degree of freedom (N-1), tnAnd (k, k +1) obeys the central t distribution with the degree of freedom of 2(N-1), and when the judgment condition is met, the Power Spectral Densities (PSDs) of adjacent spectral lines k and k +1 of the characteristic sample belong to the same frequency band, otherwise, the power spectral densities do not belong to the same frequency band.
Merging adjacent spectral lines belonging to the same frequency band in the same frequency band to form H1A frequency band, wherein the spectral line number of two end points of the h-th frequency band is kh1、kh2,h=1,2,…,H1The number of spectral lines is Nh=kh2-kh1+1。
106-5, calculating to obtain a weighted mean value and a weighted variance in the h frequency band according to the weighted mean value and the weighted variance, wherein the calculation formula is as follows:
106-6, calculating according to the weighted mean and the variance in the h frequency band to obtain the straight spectral tolerance upper limit estimation in the h frequency band, wherein the calculation formula is as follows:
wherein, F12The second tolerance upper limit coefficient with the confidence coefficient of (1-alpha) and the quantile point of beta is calculated by the following formula:
wherein Z isβTo satisfy Prob [ Z.ltoreq.Zβ]A normally distributed quantile of beta,with a degree of freedom of (N.N)hX of-1)2The distribution of the alpha quantile points is as follows,with a degree of freedom of (N.N)h-1) central t distribution 1-a quantiles;
106-7, taking the upper limit estimation of the tolerance of the flat spectrum in the h-th frequency band as the amplitude of the frequency band, connecting each adjacent flat spectrum by a straight line under a double logarithmic coordinate to obtain a comprehensive working condition standard spectrum FGG for measuring all vehicle speeds and road conditions of the points, thereby obtaining a single working condition standard spectrum FGGSecond reserve factor cij (2)And a first relation between the behaviour-synthesis specification spectrum FGG.
108, normalizing the spectrum under the single working conditionSecond reserve factor cij (2)Substituting a second relation between the comprehensive working condition standard spectrum FGG and the single working condition standard spectrumSecond reserve factor cij (2)And obtaining a second reserve coefficient c according to a first relation between the comprehensive working condition normalized spectrum FGGij (2)And the behaviour-wide specification spectrum FGG.
according to the method for constructing the armored vehicle comprehensive working condition vibration spectrum, the measured spectrum and the standard spectrum of the comprehensive working condition are respectively deduced by using an induction method and a formula calculation method, the relationship between the first storage coefficient and the measured spectrum of the comprehensive working condition and the relationship between the second storage coefficient and the standard spectrum of the comprehensive working condition are respectively obtained by using the conclusions of the two methods, so that the RBF neural network is respectively trained based on the relationships to obtain the first storage coefficient estimation value and the second storage coefficient estimation value, and then the storage coefficient adopted when the vibration spectrum is actually constructed is calculated according to the estimation values, so that the construction precision is improved.
Preferably, the process of training the first/second RBF neural network using the inputs and the desired outputs is:
definition of X ═ (X)1,x2,…,xn)TFor the network input vector, Y ═ Y1,y2,…,ys)TIs output by the network, [ phi ]i() is the radial basis function of the ith hidden layer node. The distribution function of the RBF neural network is:
where m is the number of hidden layer neuron nodes, i.e. the number of radial basis function centers, and the coefficient wiIs a connection weight;
wherein φ (#) is a radial basis function, | | x-ciI is the Euclidean norm, ciIs the ith center, xi, of the RBFiIs the ith of RBFRadius, the available network output is:
thus, the matrix expression for an RBF network can be expressed as:
D=HW+E,
wherein the desired output vector is D ═ D (D)1,d2,…,dp)TThe error vector between the desired output and the network output is E ═ E (E)1,e2,…,ep)TThe weight vector is W ═ W1,w2,…,wm)TThe regression matrix is H ═ H1,h2,…,hm)T;
Taking into account the influence of all training samples, ci、ξiAnd wiThe adjustment amounts of (a) and (b) are:
wherein phi isi(xj) For the ith implicit node pair xjInput of η1、η2、η3Respectively corresponding learning rates, ci(t) and ci(t +1) c at the t-th and t + 1-th iterations, respectivelyi,ξi(t) and xii(t +1) is ξ for the t-th and t + 1-th iterations, respectivelyi,wi(t) and wi(t +1) w at the t-th and t + 1-th iterations, respectivelyi(ii) a Obtaining a mean square error according to the cost function E so as to finish the training condition; when actually outputtingAnd the mean square error of the expected output is less than the set threshold, the network is considered to be trained.
Example 2
The embodiment provides a method for constructing a long life test spectrum of an armored vehicle, as shown in FIG. 4, comprising the following steps:
wherein, TsSetting the actual running time according to the actual service life requirement of the armored vehicle, wherein m is an index and is generally 3-9;
according to the method for constructing the long-life test spectrum of the armored vehicle, the acceleration proportional coefficient is segmented according to the road condition and the vehicle speed, so that the method is more consistent with the actual vibration data measurement test condition, the accuracy is improved, the real vibration environment can be more effectively simulated, and the test result is more real and reliable.
It should be understood that the above examples are only for clarity of illustration and are not intended to limit the embodiments. Other variations and modifications will be apparent to persons skilled in the art in light of the above description. And are neither required nor exhaustive of all embodiments. And obvious variations or modifications derived therefrom are intended to be within the scope of the invention.
Claims (10)
1. A method for constructing a comprehensive working condition vibration spectrum of an armored vehicle is characterized by comprising the following steps:
obtaining the road surface occupation ratio b of the armored vehicle under the ith road condition under the test comprehensive working conditioniAnd the vehicle speed ratio a of the jth vehicle speed used for driving on the ith road conditionijWherein, in the step (A), p is the total number of road surface types, and q is the total number of vehicle speed types;
respectively obtaining the single-working-condition actual measurement spectrum of the mth sample under the road condition i and the vehicle speed jFrom the measured spectrum of the simplex conditionAnd a first reserve coefficient cij (1)Obtaining the comprehensive working condition actual measurement spectrum FGS of all vehicle speeds and road conditions of the measuring points by induction, and obtaining the single working condition actual measurement spectrumFirst reserve factor cij (1)And a first relation between the measured spectrum FGS and the comprehensive working condition;
from the measured spectrum of the simplex conditionAnd a first reserve coefficient cij (1)ComputingObtaining the comprehensive working condition actual measurement spectrum FGS of all vehicle speeds and road conditions of the measuring points, and obtaining the single working condition actual measurement spectrumFirst reserve coefficient cij (1)And a second relation between the measured spectrum FGS and the comprehensive working condition;
single working condition actual measurement spectrumFirst reserve factor cij (1)Substituting a second relation between the measured spectrum and the comprehensive working condition actual measurement spectrum FGS into the single working condition actual measurement spectrumFirst reserve factor cij (1)Obtaining a first storage coefficient c according to a first relation between the measured spectrum and the comprehensive working condition FGSij (1)And the relation between the measured spectrum FGS of the comprehensive working condition;
the measured spectrum FGS of the comprehensive working condition is used as the input of the first RBF neural network and is based on the first reserve coefficient cij (1)A first reserve coefficient c obtained by the relation between the measured spectrum FGS and the comprehensive working conditionij (1)As expected output of the first RBF neural network, training the first RBF neural network by utilizing input and expected output, and inputting the actually measured spectrum of the current comprehensive working condition into the trained first RBF neural network to obtain a current first reserve coefficient;
respectively obtaining the standard spectrum of the single working condition of the mth sample under the conditions of the road condition i and the vehicle speed jAccording to a single working condition standard spectrumAnd a second reserve coefficient cij (2)Obtaining the comprehensive working condition standard spectrum FGG of all vehicle speeds and road conditions of the measuring points by induction, and obtaining the single working condition standard spectrumSecond reserve factor cij (2)And a first relation between the comprehensive working condition specification spectrum FGG;
according to a single working condition standard spectrumAnd a second reserve coefficient cij (2)Calculating to obtain comprehensive condition standard spectrum FGG of all vehicle speeds and road conditions of the measured points, and obtaining single-condition standard spectrumSecond reserve factor cij (2)And a second relation between the comprehensive working condition specification spectrum FGG;
standard spectrum of single working conditionSecond reserve factor cij (2)Substituting a second relation between the comprehensive working condition standard spectrum FGG and the single working condition standard spectrumSecond reserve factor cij (2)And obtaining a second reserve coefficient c according to a first relation between the comprehensive working condition normalized spectrum FGGij (2)And the relation between the comprehensive working condition standard spectrum FGG;
taking the comprehensive working condition standard spectrum FGG as the input of a second RBF neural network, and according to a second reserve coefficient cij (2)Second reserve coefficient c obtained by relation with comprehensive working condition standard spectrum FGGij (2)As expected output of the second RBF neural network, training the second RBF neural network by utilizing input and expected output, and inputting the current comprehensive working condition standard spectrum into the trained second RBF neural network to obtain a current second reserve coefficient;
calculating to obtain a reserve coefficient c according to the current first reserve coefficient and the current second reserve coefficientij;
According to the reserve coefficient cijActual measurement spectrum G of single working condition under road condition i and vehicle speed jshice(i, j) and the Standard spectra of Single working cases Gguifan(i, j), respectively calculating and obtaining comprehensive working condition actual measurement spectrums FGS and comprehensive working condition standard spectrums FGG of all vehicle speeds and road conditions of the measured points.
2. The method of claim 1, wherein the measured spectra are based on single-site conditionsAnd a first reserve coefficient cij (1)The step of inducing to obtain the comprehensive working condition actual measurement spectrum FGS of all the vehicle speeds and road conditions of the measuring points comprises the following steps:
from the measured spectrum of the simplex conditionAnd a first reserve coefficient cij (1)Calculating to obtain the mean value and the variance under the road condition i and the vehicle speed j, wherein the calculation formula is as follows:
wherein N isi,jThe number of samples under a single working condition;
calculating to obtain weighted mean value and variance according to the mean value and variance under the road condition i and the vehicle speed j, wherein the calculation formula is as follows:
and calculating according to the weighted mean value and the variance to obtain a comprehensive working condition actual measurement spectrum FGS of all the vehicle speeds and road conditions of the measuring points, wherein the calculation formula is as follows:
wherein, F11The first tolerance upper limit coefficient is calculated by the following formula:
3. the method of claim 1, wherein the measured spectra are based on single-site conditionsAnd a first reserve coefficient cij (1)The calculation formula of the comprehensive working condition actual measurement spectrum FGS for calculating and obtaining all vehicle speeds and road conditions of the measuring points is as follows:
wherein N isi,jIs the number of samples in the simplex case.
4. The method of claim 1, wherein the spectrum is normalized according to a single regimeAnd a second reserve coefficient cij (2)Comprehensive worker for obtaining all vehicle speeds and road conditions of measuring points through inductionThe step of the condition-normalized spectrum FGG comprises:
according to a single working condition standard spectrumAnd a second reserve coefficient cij (2)Calculating to obtain the mean value and the variance of each spectral line under the road condition i and the vehicle speed j, wherein the calculation formula is as follows:
wherein N isi,jThe number of samples under a single working condition;
calculating to obtain a weighted mean value and a weighted variance according to the mean value and the variance of each spectral line under the road condition i and the vehicle speed j, wherein the calculation formula is as follows:
obtaining a statistic F by weighted mean and variance calculationnAnd tnThe calculation formula is as follows:
according to statistic FnAnd tnFrequency division is carried out under the confidence coefficient (1-alpha) to obtain H1A frequency band, wherein the spectral line number of two end points of the h-th frequency band is kh1、kh2,h=1,2,…,H1The number of spectral lines is Nh=kh2-kh1+1;
And calculating to obtain a weighted mean value and a weighted variance in the h frequency band according to the weighted mean value and the weighted variance, wherein the calculation formula is as follows:
calculating according to the weighted mean and the variance in the h frequency band to obtain the straight spectrum tolerance upper limit estimation in the h frequency band, wherein the calculation formula is as follows:
wherein, F12The second tolerance upper limit coefficient with the confidence coefficient of (1-alpha) and the quantile point of beta is calculated by the following formula:
wherein Z isβTo satisfy Prob [ Z.ltoreq.Zβ]A normally distributed quantile of beta,with a degree of freedom of (N.N)hX of-1)2The distribution of the alpha quantile points is,with a degree of freedom of (N.N)h-1) central t distribution 1-a quantiles;
and taking the upper limit estimation of the tolerance of the straight spectrum in the h-th frequency band as the amplitude of the frequency band, and connecting each adjacent straight spectrum by using a straight line under a log-log coordinate to obtain a comprehensive working condition standard spectrum FGG of all vehicle speeds and road conditions of the measured point.
5. The method of claim 1, wherein the spectrum is normalized according to a single regimeAnd a second reserve coefficient cij (2)The calculation formula of the comprehensive working condition standard spectrum FGG for calculating and obtaining all vehicle speeds and road conditions of the measured points is as follows:
wherein, Ni,jIs the number of samples in the simplex case.
6. Method according to claim 1, characterized in that the reserve coefficient cijThe root mean square values of the current first reserve coefficient and the current second reserve coefficient.
7. The method of claim 1, wherein said determining is based on said reserve coefficient cijActual measurement spectrum G of single working condition under road condition i and vehicle speed jshice(i, j) and the simplex condition normalized spectrum Gguifan(i, j), respectively calculating and obtaining the calculation formulas of the comprehensive working condition actual measurement spectrum FGS and the comprehensive working condition standard spectrum FGG of all vehicle speeds and road conditions of the measured points, wherein the calculation formulas are respectively as follows:
8. a method for constructing a long life test spectrum of an armored vehicle is characterized by comprising the following steps:
obtaining the time T of the jth vehicle speed used for driving under the ith road conditionijI is 1,2, …, p, j is 1,2, …, q, p is the total number of road surface types, q is the total number of vehicle speed types;
time T of jth vehicle speed used for driving according to ith road conditionijAnd calculating to obtain a segment acceleration proportional coefficient beta under the road condition i and the vehicle speed jij;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110405031.5A CN113094816B (en) | 2021-04-15 | 2021-04-15 | Method for constructing comprehensive working condition vibration spectrum and long-life test spectrum of armored vehicle |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110405031.5A CN113094816B (en) | 2021-04-15 | 2021-04-15 | Method for constructing comprehensive working condition vibration spectrum and long-life test spectrum of armored vehicle |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113094816A CN113094816A (en) | 2021-07-09 |
CN113094816B true CN113094816B (en) | 2022-07-05 |
Family
ID=76677816
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110405031.5A Active CN113094816B (en) | 2021-04-15 | 2021-04-15 | Method for constructing comprehensive working condition vibration spectrum and long-life test spectrum of armored vehicle |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113094816B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115950647B (en) * | 2022-12-19 | 2023-10-27 | 中国人民解放军陆军装甲兵学院 | Monitoring, diagnosing and analyzing method for scrapping of electromechanical equipment |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110929334A (en) * | 2019-10-30 | 2020-03-27 | 中国北方车辆研究所 | Method for calculating transmission efficiency of track armored vehicle under gear shifting dynamic working condition |
CN111597673A (en) * | 2019-02-21 | 2020-08-28 | 株洲中车时代电气股份有限公司 | Random vibration fatigue acceleration test method and system |
CN112595411A (en) * | 2020-12-15 | 2021-04-02 | 华域麦格纳电驱动系统有限公司 | Vibration spectrum data acquisition and processing method and device of electric drive system |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8195358B2 (en) * | 2008-09-11 | 2012-06-05 | Deere & Company | Multi-vehicle high integrity perception |
US8818567B2 (en) * | 2008-09-11 | 2014-08-26 | Deere & Company | High integrity perception for machine localization and safeguarding |
-
2021
- 2021-04-15 CN CN202110405031.5A patent/CN113094816B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111597673A (en) * | 2019-02-21 | 2020-08-28 | 株洲中车时代电气股份有限公司 | Random vibration fatigue acceleration test method and system |
CN110929334A (en) * | 2019-10-30 | 2020-03-27 | 中国北方车辆研究所 | Method for calculating transmission efficiency of track armored vehicle under gear shifting dynamic working condition |
CN112595411A (en) * | 2020-12-15 | 2021-04-02 | 华域麦格纳电驱动系统有限公司 | Vibration spectrum data acquisition and processing method and device of electric drive system |
Non-Patent Citations (2)
Title |
---|
基于BP神经网络的装甲车辆发动机使用状态评价;张会奇等;《兵工自动化》;20140115(第01期);第32-34页 * |
装甲车辆发动机载荷谱仿真与试验研究;骆清国等;《兵工学报》;20121215(第12期);第1442-1447页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113094816A (en) | 2021-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110417011B (en) | Online dynamic security assessment method based on mutual information and iterative random forest | |
US20080228338A1 (en) | Automated engine data diagnostic analysis | |
CN111582718A (en) | Cable channel fire risk assessment method and device based on network analytic hierarchy process | |
CN110263429B (en) | Method for adjusting shape and surface precision of cable net reflecting surface antenna based on machine learning | |
EP4075281A1 (en) | Ann-based program test method and test system, and application | |
CN112883497B (en) | Space valve reliability assessment method based on multi-source information fusion | |
CN110516821B (en) | Ship maintenance resource allocation optimization method and device | |
CN113094816B (en) | Method for constructing comprehensive working condition vibration spectrum and long-life test spectrum of armored vehicle | |
CN114583767B (en) | Data-driven wind power plant frequency modulation response characteristic modeling method and system | |
CN110850358A (en) | Electric energy meter comprehensive verification method and system based on stepwise regression algorithm | |
CN109389517B (en) | Analysis method and device for quantifying line loss influence factors | |
CN114064459A (en) | Software defect prediction method based on generation countermeasure network and ensemble learning | |
CN113255880A (en) | Method and system for judging electricity stealing data based on improved neural network model | |
CN117076887A (en) | Pump station unit running state prediction and health assessment method and system | |
CN113779893B (en) | Reliability evaluation method for guaranteed simulation output result based on echo state network | |
CN113063589B (en) | Gear microscopic error vibration prediction method based on neural network | |
CN115755831A (en) | Fault diagnosis method and system based on fuzzy hierarchical analysis | |
CN108280299A (en) | truss structure sensitivity analysis method | |
CN108985564B (en) | QFD (quad Flat No-lead) and binary semantic-based FMEA (failure mode and effects analysis) method | |
CN110852536A (en) | Storage tank maintenance decision determining method and device | |
CN107886248A (en) | A kind of equipment state big data computational methods and equipment based on Golstein BP algorithms | |
CN110321239B (en) | System reliability evaluation method based on system survival characteristics | |
CN112653129B (en) | Transient power angle stability margin estimation method, device and system | |
Küçüksakarya et al. | PANEL DATA ANALYSIS OF RELATIONSHIP BETWEEN ECONOMIC GROWTH, FOREIGN DIRECT INVESTMENT, EXCHANGE RATE AND TRADE OPENNESS IN NEWLY INDUSTRIALIZED COUNTRIES | |
Li et al. | Simulation of Reliability Prediction Based on Multiple Factors for Spinning Machine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |