Intelligent Time-Domain Parameters Matching for Shock Response Spectrum and Its Experimental Validation in Active Vibration Control Systems

Shock response spectrum (SRS) test is a kind of vibration test that uses the principle of equivalent damage to simulate complex shock and vibration environment for evaluating product fragility. However, for a certain SRS, neither an analytical nor a unique inverse time-domain shock waveform (TSW) exists, making the SRS test a very challenging problem. Synthesis of a TSW using a group of wavelets with diﬀerent frequencies, amplitudes, and phases is considered to be a very promising method. However, it is challenging to ﬁnd an optimal group of wavelets since there are hundreds of optimization variables and objective functions. In this paper, a novel intelligent parameter matching method for TSW synthesis is proposed for the SRS test. A variable-weight method has been introduced to combine the eﬀects of hundreds objective functions so that the optimization problem can be solved by using the genetic algorithm. The frequency response function of the shaking system has been taken into consideration in the calculation of the objective function, so the synthesized TSW can be applied directly to the SRS test. In the optimization process, normalized control factors have been formulated for the optimization variables so that they can be tuned in a reasonable small range to avoid irrationality and accelerating convergence in searching process. The eﬀectiveness of the proposed method has been veriﬁed experimentally in two active vibration control systems, in which one contains a 500N minishaker and the other contains a 1-ton shaking table. It can be seen from the experimental results that the proposed method can accurately synthesize the input TSW for the shaking system, where the output TSW and SRS can accurately meet the time-and frequency-domain speciﬁcation.


Introduction
Shock load is a kind of load generally all kinds of products will undergo in the process of manufacture, transportation, and use.Shock often causes natural frequency response or damage to equipment, which greatly reduces product reliability.e shock test is performed to determine the environmental adaptability, structural reliability, and shock resistance of the product under polytypic mechanical shock load.e control strategy for the shock test can be divided into classic waveform control and shock response spectrum (SRS) control.e former one mainly uses drop hammer or pendulum shock machine [1] to produce simple shock waveforms, such as half-sine wave, rectangular wave, trapezoidal wave, and sawtooth wave [2,3].It pays much attention to the time-domain shock waveform (TSW) itself, which is not able to effectively simulate the real environmental shock [4].Moreover, this method has larger lowfrequency energy, which is easier to damage the experimental equipment.e SRS control focuses on the effect of the shock to the device under test (DUT) and adopts the principle of equivalent damage to generate complex TSWs, which can truly and effectively simulate the real environmental shock.
SRS is widely applied for evaluating the product fragility and the vibration environment for the protective packaging system [5], electronic equipment [3,6], weapons [7], spacecraft and launch vehicle [4], high-speed craft, bridges, etc.For example, Gannon [8] developed a repeatable laboratorial method using SRS for testing shock mitigating seats on high-speed watercraft.George et al. [9] developed an SRS test system of a launch vehicle using LabVIEW.Wang et al. [10] conducted pyroshock response analysis of separation nuts in the aerospace launch system based on the finite element (FE) model and SRS.Haciefendioglu et al. [11] conducted the dynamic assessment of partially damaged historic masonry bridges under blast-induced ground motion using the multipoint SRS method and FE model.Stewart et al. [12] analyzed the response of the hydraulic blast simulator used for generating high-g shock based on the SRS method.Botta et al. [13] investigated the SRS in plates under impulse loads, considering the Mindlin-Reissner model.Li et al. [14] studied the damage boundary of a beam under shock based on the structural dynamics and the SRS analysis.
SRS can also be applied for vibration/noise suppression and even be used for fault diagnosis.For example, Okuyama et al. [15] investigated a design method based on SRS analysis for reducing acoustic noise of a seeking control system in a hard disk driver.Atsumi [16] evaluated the settling vibration and acoustic noise caused by the acceleration feedforward input of the hard disk driver using SRS analysis.Shin et al. [17] investigated two simple methods to suppress the residual vibrations of a translating or rotating flexible cantilever beam, where both methods are based on the SRS of the system.Hussain et al. [18,19] proposed a gearbox fault detection method using the real coded genetic algorithm (GA), adaptive wavelet filtering, and novel SRS feature extraction.
With the development of computer, digital signal processing, and active vibration control technology [20][21][22][23][24][25][26][27][28][29], the SRS test with electrodynamic shaker has become the mainstream.Compared with the traditional shock test equipment, the electrodynamic shaker has the advantages of a simple setup, better controllability and repeatability, convenient adjustment, and high accuracy, which helps in easy realization of the control of SRS [4].In the SRS test, accurate generation of the TSW is the key problem, where there is lot of research devoted to it.For example, Monti et al. [30] investigated a GA-based dynamic load synthesis method for SRS and verified the FE model of space equipment.James et al. [31] proposed a stochastic SRS decomposition method based on probabilistic definitions of temporal peak acceleration, spectral energy, and phase lag distributions of mechanical impact pyrotechnic shock test data.Brake [32] investigated the inverse of SRS by using the Morlet wavelet as the basis function and using GA to find the optimal parameters.Soize [33] solved the inverse problem of SRS from the perspective of maximum entropy principle, in which a random generator of accelerograms was developed.Wattiaux et al. [34] investigated the prediction of the vibration levels generated by pyrotechnic shocks using an approach based on SRS specification and equivalent mechanical shock.
It can be seen from the literature review that the timedomain parameters matching, i.e., finding the inverse of SRS, is the critical and difficult problem since a specific SRS specification may not have unique inverse TSW.Moreover, the dynamic of the shaking system may greatly affect the control of SRS.Synthesis of a TSW using a group of wavelets with different frequencies, amplitudes, and phases is considered to be the very promising method.However, it is challenging to find an optimal group of wavelets since there are hundreds of optimization variables and objective functions.In this paper, the intelligent parameters matching method of TSW considering the SRS specification and frequency response function (FRF) of the shaking system is proposed.e multiple objective problem, considering the restriction of shock duration, shaking forces, and SRS specification, is simplified through a variable-weight method so that the optimization problem can be solved based on GA.
e TSW matching corresponding to the SRS specification is successfully obtained and has been experimentally validated in an active vibration control (AVC) system under different working conditions. is paper is structured as follows: Section 2 introduces the principle of SRS and SRS specification.Section 3 presents the parameters matching process of TSW with regard to SRS.Section 4 presents the parameters optimization using the proposed objective function and GA.Section 5 conducts the experimental validation of SRS in an AVC system with minishaker.Section 6 conducts the experimental validation of SRS in an AVC system with 1-ton shaking table.Section 6 briefly concludes the whole paper.

Principle of SRS
Shock response spectrum (SRS) is a graphical representation of a shock, or any other transient acceleration input, in terms of how a single-degree-of-freedom (SDOF) system (like a mass on a spring) would respond to that input.In real application, a continuous system can be approximately represented as a series of the discrete linear SDOF system [35].Applying the same shock excitation to these SDOF systems, the relationship between the peak responses and the nature frequency of these SDOF systems is SRS, as shown in Figure 1. e horizontal axis in Figure 1 shows the natural frequency of the hypothetical SDOF systems, and the vertical axis shows the peak acceleration of these SDOF systems with same shock input.In Figure 1, there are M SDOF systems fixed on the same foundation, where m i , k i , c i , and f i are the mass, stiffness, damping, and natural frequency of the i th (i � 1, 2, . .., M) SDOF system, € x s (t) is the foundation acceleration excitation, € x i (t) is the acceleration response of the i th SDOF system, and A i is the maximum amplitude of the acceleration response of the i th SDOF system.e natural frequency f i and maximum amplitude A i will determine one point in SRS.Connecting all the M points (f i , A i ) i�1,2,...,M forms the integrated SRS.
It can be seen that SRS is not simply a spectrum of shock response of the continuous system.It describes maximum responses of each frequency component of the system.us, it is able to reflect the harmful effect that shock brings about to the system and offers a criterion to the damage of the shock test.SRS specification is the standard to judge whether the product has passed certain strength examination in the shock test [30].Figure 2 shows a typical SRS specification.It includes a standard curve and its upper and lower tolerance limit envelope, which is composed of low-frequency rising 2 Shock and Vibration section and high-frequency horizontal section.In Figure 2, f max and f min are the maximum and minimum values of analysis frequency, f c is the frequency of inflection point, and Φ 1 and Φ 2 are the slope of the curve.
In the SRS test, the analysis frequency f i is selected in a logarithmic interval.
e intervals between analysis frequencies is generally selected to be 1/6 octave or smaller (1/ 12, 1/24, etc.), according to the requirement of accuracy.For simplicity, the damping ratio of each SDOF system is specified to be equal according to the damping characteristics of the DUT.e response of the SDOF system is the basis for the calculation of SRS. e dynamic model of the SDOF system is illustrated in Figure 3, where € x s (t) is the acceleration of the foundation (equivalent to the shock signal) and € x i (t) is the acceleration response of the SDOF system.e dynamic equation of the SDOF system can be written as Because the actual excitation is not simply harmonic or impulse, the traditional Duhamel integral method (DIM) is of large computational complexity.In this paper, an improved recursive digital filtering method (IRDFM), which is firstly proposed by Smallwood [36], has been applied.We set the discrete sampling value of acceleration excitation € x s (t) as X(k • Δt) and that of acceleration response € x i (t) as Y(k • Δt), respectively, where Δt is the sampling interval.e recursive relation can be written as (2) e coefficients of the above equation are defined as Shock excitation SDOF systems Shock response of each system Shock response spectrum (SRS) Amplitude Frequency (Hz) Acceleration (g)

Upper and lower tolerance limit envelopes
Frequency (Hz)

Shock and Vibration
where is the undamped natural frequency, ζ � c/2mω n is the damping ratio, and

􏽱
is the damped natural frequency.e above method will be better than the DIM and finite difference method (FDM) in speed and accuracy.
e error of sinusoidal response and SRS between theoretical solution and the three numerical solutions are compared, as shown in Figure 4.
It can be seen that the error of the IRDFM is obviously smaller than that of the other two methods, and the error of the FDM is the largest.At the same time, the time consumption of the three numerical methods are compared, where the code is running on the PC with Intel Core i7-4790 CUP @ 3.60 GHz and 16G RAM, number of analysis frequency is 288, sampling frequency is 48000 Hz, and sampling number is 7680.
e result shows that DIM cost 1.346 s, FDM cost 0.121 s, and IRDFM cost 0.098 s.In summary, the IRDFM is superior in both accuracy and efficiency.

Parameters Matching of TSW
In the SRS test, the first step is to determine test specification and then solve an inverse of SRS, i.e., TSW, to driving the electrodynamic shaker and making the SRS of DUT matching the predefined test specification.e SRS corresponding to a certain TSW is unique, while the TSW corresponding to a certain SRS is not unique.In order to achieve the test goal, an applicable TSW needs to be constructed according to a determined SRS specification, which is known as time-domain parameters matching.ree principles should be followed in the time-domain parameters matching.ey are (1) the duration of the TSW should be short and not exceed the prescribed limit; (2) e maximum acceleration of the TSW should be small and not exceed the limit of the shaker; (3) e SRS corresponding to the TSW should conform to the test specification and not exceed the prescribed tolerant limit envelope.
e commonly used time-domain parameters matching methods include attenuated sine wave method, attenuated cosine wave method, and wavelet synthesis method.Because the wavelet synthesis method is more flexible in adjusting the parameters and the time of the shock waveform is easier to control, it is more preferable in application.e wavelet synthesis method was first proposed by the U.S. military.A series of windowed sinusoidal waves are superimposed to form the TSW, which is expressed as where € x s (t) is the synthesized acceleration TSW and u i (t) is the windowed sinusoidal wave (called wavelet), which can be expressed as follows:  Shock and Vibration where A i is the amplitude of the wavelet, f i is the frequency of the wavelet, τ i is the delay of the wavelet, N i is the number of half sinusoids in the wavelet, and T i is the duration of the wavelet.In order to satisfy the above three principles, the parameters of each wavelet must be selected reasonably.e methods to determine the parameters of wavelet are listed as follows.
3.1.Wavelet Frequency f i .e wavelet frequency f i is determined according to analysis bandwidth.It is sampled at an equal logarithmic interval within specified analysis bandwidth of the test specification.Each wavelet frequency is corresponding to an analysis point in the SRS.If the interval of the wavelet frequency is 1/6 dB/Oct, the number of total analysis points is where f min and f max are the lower and upper frequency limits of the SRS specification and ⌈ • ⌉ is the upward integer symbol.en, analysis frequency can be represented as

Number of Half Sinusoids N i and Duration of Wavelet
T i .e number of half sinusoids N i is the sum of peaks and troughs of each wavelet.e relationship between N i and the duration of wavelet T i is Generally, the duration of TSW in the SRS test is not greater than 20 ms, i.e., T max ≤ 20 ms.Since the duration of the wavelet cannot exceed that of the TSW, it can be evaluated according to the following formula, i.e., In addition, in order to satisfy the symmetry of wavelet, N i is required to take an odd number which is greater than or equal to 3.After the value of N i is determined, the value of T i can be determined according to equation (8).

Wavelet Delay τ i .
e wavelet delay τ i determines the phase of each wavelet in superposition, which is a key parameter affecting the maximum value of the synthesized TSW.For simplicity, each wavelet is generally placed in the middle of time range of the TSW, i.e., And the wavelet delay should satisfy 3.4.Wavelet Amplitude A i .e wavelet amplitude is generally evaluated according to the following empirical formula, i.e., where A SRS,i and A FRF,i are the values of the SRS specification and frequency response function (FRF) of the shaker, respectively, corresponding to the analysis frequency.
It can be seen that the synthesized TSW is affected by five parameters.However, because the analysis frequency f i is determined according to the SRS specification, and there is a clear relationship between the duration of the wavelet T i and the number of half-sinusoidal waves N i , there are mainly three independent variable parameters, i.e., A i , N i , and τ i , that affect the shock waveform.After these parameters are determined, a TSW can be synthesized.However, the SRS corresponding to this waveform often does not conform to the three principles above.Generally, there are two main problems: one is that the calculated value of the SRS does not meet the requirements of the SRS specification; the other one is that the maximum value of the TSW exceeds the limit of the shaker.It can be seen from the analysis that the maximum value of the TSW and the corresponding SRS are affected by the parameters of A i , N i , and τ i at the same time.Modifying a certain parameter will affect both targets simultaneously, so it is difficult to match proper parameters of TSW corresponding to a certain SRS specification.

Parameters Optimization Based on GA
e genetic algorithm (GA) is applied in this paper for parameters optimization of TSW. is algorithm is a kind of the random search optimization algorithm that simulates the evolution rule in nature.It can be seen from the above analysis of wavelet synthesis that there is a great degree of subjectivity and randomness in parameter determination and it is difficult to find the appropriate parameters artificially.If the advantages of automatic searching by using GA are combined with wavelet synthesis, the shortcomings in parameter selection can be overcome.As presented before, the optimization problem is a multiple objective problem, where there are two main restrictions.One is that the maximum value of TSW is restricted by the limit of the shaker.Actually, we want the amplitude of TSW as small as possible, to reduce the energy consumption.e other restriction is that the SRS of the shock waveform must agree with the test specification.Although we defined the tolerant limit to the SRS specification, we still want the error as smaller as possible.at is to say, the error between actual SRS and the specification at every analysis frequency should meet the requirement.So the purpose of GA is to obtain a Shock and Vibration suboptimal solution of time-domain parameters by considering the above restrictions.e optimization process of this problem using GA is as follows.

Chromosome Coding (Variable Regularization).
According to the above analysis, the parameters affecting the optimization objectives are A i , N i , and τ i , which can be used as optimization variables.Since the empirical formula for calculating the three factors has been given in the wavelet synthesis method, three control factors, i.e., α A i , α N i , and α τ i , are introduced for optimization variables.So the optimization variables can be redefined as where in equation ( 14), a odd denotes the largest odd number under a.In this way, α A i , α N i , and α τ i can be used as optimization factors to control the three parameters A i , N i , and τ i so that the three parameters can be adjusted in a relatively reasonable small range to avoid the irrationality in searching process.At the same time, in order to facilitate operations on chromosomes, the three control factors are normalized, i.e., e control factors are also called gene.en, the chromosome could consist of genes, which is expressed as e chromosome is a (1 × 3M)-dimensional vector containing the normalized control factors (gene) of all variables at all analysis frequencies; M is the number of analysis frequencies.In the optimization process, the chromosome is also called individuals.

Optimization Objective.
e parameters matching of TSW is not simply a single objective optimization problem.
ere are multiple optimization objectives.e first objective is that the maximum amplitude of the TSW must be within the range of the shaker, and the smaller the better.
us, the first objective can be defined as where g max is the maximum amplitude of the shaker and a max is the maximum amplitude of the TSW.e second objective is that the SRS corresponding to the TSW must conform to the SRS specification.at is to say, the error between actual SRS and the specification at every analysis frequency must be within the tolerance limit, and the smaller the better.Since there are M analysis frequency, the second objective can be subdivided into M subobjective, i.e., where A SRS,i is the amplitude of the SRS specification at the i th analysis frequency, A FRF,i is the amplitude of the shocking system FRF at the i th analysis frequency, and A S,i is the amplitude actual output SRS at the i th analysis frequency.e parameters should be optimized so that all the objectives s i , (i � 0, 1, 2, . . ., M) can be minimized.In the SRS test, there are generally hundreds of analysis frequencies, so there are also hundreds of objective functions, making the parameters matching a very challenging problem.

Fitness Function.
Since there are (M + 1) objective functions, a reasonable fitness function f(s 0 , s 1 , s 2 , . . ., s M ) should be designed so that under the condition minf(s 0 , s 1 , s 2 , . . ., s M ), the minimum value of s i , (i � 0, 1, 2, . . ., M) can be guaranteed.In this paper, the variable-weight method is applied to combine multiple objectives.e fitness function can be obtained by combining all the objectives, i.e., where f(s 0 , s 1 , s 2 , . . ., s M ) is the fitness function and c i (i � 0, 1, 2, . . ., M) are the weight coefficients of each objective function.If the weight coefficients are constant, the fitness function may not be able to fully reflect the situation of the objectives, i.e., M + 1 objectives may not simultaneously be minimized when the combined function are minimized.erefore, it is necessary to use an adaptive weighting method to amplify the biased objective, i.e., to adjust the weight adaptively according to the deviation value of each objective.
e greater the deviation, the greater the weight will be.e exponential weighted fitness function is proposed by making use of the explosive growth characteristic of the exponential function to give greater weight to the objective function with large deviation so as to increase the weight of the objective function in the fitness function.e weight coefficients can be expressed as is weight coefficient can guarantee that the fitness function is small only when all objective functions are small.
In order to verify the performance of the above exponential weighted fitness function, the following three weighted fitness functions are compared: f(s 0 , s 1 ) �  i�0,1 exp(s i ) • s i (exponential weight), f(s 0 , s 1 ) �  i�0,1 s 2 i • s i (square weight), and f(s 0 , s 1 ) �  i�0,1 s i (linear weight).In order to quantitatively 6 Shock and Vibration compare the differences among the three fitness functions, the contour lines of the two target values of the three at the fitness function of 0.5, 1, 2, and 3 are also illustrated.e comparison results are shown in Figure 5.
It can be seen from Figure 5 that the value of the exponential weighted fitness function is much greater than that of the other two fitness functions, which amplifies the deviation of each objective function.It can be predicted that the exponential weighted fitness function is more effective than the other two when the deviation of the arbitrary objective function is larger.It can be seen from the contour map that, when the fitness function value is 2, the maximum deviation of the single target in the exponential weight is about 0.8, the maximum deviation of the single target in the square weight is about 1.2, and the maximum deviation of the single target in the linear weight is 2.
erefore, when the fitness function has the same value, the exponential weighted fitness function can guarantee that all objective functions are near their minimum value.

Generate the Initial Population.
e initial population is generated by the chaotic search algorithm.First, a chromosome (individual) is generated randomly and then the initial population is obtained by the chaotic iterative algorithm.In this paper, the following logistic mapping algorithm is adopted for chaos calculation of gene, i.e., where μ is the control parameter and λ is the random variable (gene), where 0 < λ < 1 and λ ≠ 0.25, 0.5, 0.75.Suppose the length of the chromosome is 3M, firstly select 3M different values in the interval (0, 1) as the gene of the first chromosome of the first generation.en iteratively generate new genes through equation ( 22) to form new chromosome, until the intimal population of first generation are generated.

Selection.
In GA, selection is the process of preserving the dominant individuals (chromosomes) and eliminating the inferior individuals.In this paper, controlling population and roulette strategy are adopted for selection.In this process, first calculate the fitness of each chromosome, and then pick up the part of the advantage individuals as control population, which do not take part in crossing and variation.e rest of the individuals will be involved in crossing and variation normally.Finally, compare the results of the two parts and eliminate disadvantage individuals by "roulette wheel" way, getting new species for the next generation  Shock and Vibration of evolution.By controlling the population, the dominant individuals in the population can be protected from dying out due to crossing and variation, and the convergence speed of the algorithm can be improved.

Crossing and Variation.
Crossing and variation is an important part to generate population diversity, which can avoid the nonconvergence caused by population monotony in the process of evolution.e crossing and variation probabilities are important parameters, which directly affect the diversity of population.e adaptive probability is applied, i.e., the crossing and variation probabilities are kept updating with evolution.
e specific formula is as follows: where p c and p m are, respectively, the crossing and variation probabilities, p c0 and p m0 are, respectively, the initial crossing and variation probabilities, and n is the current generation.
Besides the probability, the way of crossing and variation is also an important factor affecting the diversity of population.e commonly used methods of arithmetic crossing and multipoint variation are applied.Arithmetical crossing is to randomly select two parent individuals, which are divided into two parts according to the location of the intersection, namely, . Four children can be produced by crossing, i.e., where α and β are the random numbers between [0, 1]. e fitness function of two parents and four children is calculated after the crossing, and the best two are selected to the next generation.Multipoint variation is to randomly select m variation locations on a stochastic chromosome and reproduces random numbers to replace them.e fitness of the two individuals before and after the variation is also compared to choose the better one to the next generation.

Simulation and Optimization Result.
e optimization process of TSW using GA is shown in Figure 6.Step 1 is coding the chromosome, i.e., regulate the optimization variable (also called gene) so that they are suitable for the optimization process of GA.Step 2 is to generate initial population, i.e., generate an adequate numbers of the chromosome (also called individuals).Step 3 is calculating the fitness function, considering the identified FRF and the SRS specification.Step 4 is to judge whether the fitness function of the optimal individual is smaller enough to meet the stopping condition.If YES, decode the chromosome to calculate the TSW.IF NO, continue the evolution process.Step 5 is the genetic operations, i.e., selection, crossing, and variation; Step 6 is to judge whether the optimization process is too long.If NO, continue the evolution process.If YES, break the evolution process.
In order to verify the effect of parameters matching of TSW using GA, a simulation was conducted.During the simulation, the main parameters of the GA are set as follows.Population size (i.e., the number of the chromosome in one generation) is 500.Control population size (i.e., the number of the dominating chromosome to be revered) is 20.Initial crossing probability is 0.7.Initial variation probability is 0.2.e maximum number of generation is 100.e trend of the fitness function is shown in Figure 7(a), and the synthesized TSW is shown in Figure 7(b).
It can be seen from Figure 7(a) that the GA converges very fast in the process of evolution and meets the requirement at about 80 th generation.It can be seen from Figure 7(b) that the time-domain parameters matching using GA can meet the time-domain requirement, i.e., the duration and the maximum value of TSW do not exceed the predefined limit.Although we have considered the dynamic of the shaking system in the optimization process, due to the environmental disturbance, the frequency-domain requirement, i.e., the requirement to match the SRS specification, should be validated in AVC system experimentally.

Procedures of SRS Test.
In the SRS test, the FRF between the input and output of the shaking system must be made clear to ensure that the waveform actually output to the DUT meets the test specifications.In application, it is difficult to model the shaking system in theory, and thus an identified model is applied.e optimization algorithm will optimize offline a group of parameters according to the SRS specification and identified FRF, then the TSW is synthesized to drive the shaker, an output TSW is measured from the shaker table, and then the measured TSW is transformed into SRS to see if it is in agreement with the SRS specification.If YES, the test is finished.If NO, identify the FRF again and optimize another group of parameters, as shown in Figure 9.

Case Study.
According to the characteristics of the shaker, the target SRS specification for the test is defined as 8 Shock and Vibration f min � 100 Hz, f c � 500 Hz, f max � 3000 Hz, Φ 1 � 6 dB/Oct, and Φ 1 � 25 g (250 m/s 2 ).e upper and lower tolerance envelop is with ±3 dB deviation of the specification line, the maximum TSW duration T max � 20 ms, frequency interval is 1/ 6 dB/Oct, damping ratio ζ � 0.05, and maximum acceleration of the shaker g max � 10 g(100 m/s 2 ).Before the test, the TSW conforming to the SRS test specification was simulated and calculated considering the FRF of the shaker based on GA. ree cases are investigated, which are no load, light load, and heavy load.In Case A, there is no DUT, but an accelerometer is placed on the table, as shown in Figure 10(a).e mass of movable parts of the shaker is 0.254 kg, including the weight of the accelerometer 0.004 kg. e output force under the acceleration of 10 g (100 m/s 2 ) is 21.5 N. e measured FRF is shown in Figure 10(b), and it can be seen that there is a first-

Shock and Vibration
order natural frequency about 46 Hz. e actual output TSW and its SRS of the shaking table are shown in Figures 10(c) and 10(d).It can be seen that the TSW meets the requirement of duration and maximum values.e SRS meets the requirement of test specification, i.e., the error at every analysis frequency is within the tolerance envelope.
In Case B, the shaking table is fixedly connected with a light load (0.057 kg) and an accelerator (0.004 kg), as shown in Figure 11(a).e mass of movable parts of the shaker is 0.311 kg, including the weight of DUT and accelerometer.
e output force of the shaking table under the acceleration of 10 g is 30.1 N. It can be seen that since the movable part is heavier, the output force must be larger to generate same acceleration.e measured FRF is shown in Figure 11(b), and it can be seen that there is a first-order natural frequency still around 46 Hz, which means that the light load will not strongly affect the natural frequency of the shaking system.e actual output TSW and its SRS of the shaking table are shown in Figures 11(c) and 11(d).It can be seen that the TSW meets the requirement of duration and maximum values.e SRS meets the requirement of test specification.
In Case C, the shaking table is connected with a heavy load (0.211 kg) and an accelerometer (0.004 kg), as shown in Figure 12(a).
e mass of movable parts of the shaker is 0.465 kg, including the weight of the DUT and accelerometer.
e output force of the shaking table under the acceleration of 10 g is 45.2 N. e measured FRF is shown in Figure 12(b), it can be seen that there is a first-order natural frequency about 45 Hz, which means that the heavy load will have a little bit influence on the natural frequency of the shaking system.e actual output TSW and its SRS of the shaking table are shown in Figures 12(c) and 12(d).It can be seen that the TSW meets the requirement of duration and maximum values.e SRS meets the requirement of test specification.
Furthermore, the shaking table will not work due to overcurrent protection if continues adding the load.It can be seen that under the allowable load of the shaking table, the time-domain matching result of the SRS perfectly agrees with the expected requirements.It also verifies that the timedomain matching strategy for the SRS proposed in this paper is effective.
e FRF of the system in this experiment is identified by the frequency sweep technique, and the result is shown in Figure 14.We can see that there is a resonance point at about 2400 Hz, and an antiresonance point at about 4100 Hz. e vibration controller plays the role of making the actual output of the shaker consistent with the target waveform according to the FRF.us, the vibration controller will exert larger input to the shaker at its antiresonance point.
is process is similar to Figure 9. Due to the limitation of the thrust and frequency range, the shaker cannot bear the SRS test with too high frequency or too large acceleration.Next, an experiment is carried out within the range of the shaker capacity.

Experimental Results.
According to the characteristics of the shaker, the target SRS specification for the test is defined as f min � 500 Hz, f c � 2000 Hz, f max � 6000 Hz, Φ 1 � 10 dB/Oct, and Φ 1 � 500 g (5000 m/s 2 ).e upper and lower tolerance envelop is with ±3 dB deviation of the specification line, the maximum TSW duration T max � 20 ms, frequency interval is 1/6 dB/Oct, damping ratio ζ � 0.05, and maximum acceleration of the shaker g max � 100 g (1000 m/s 2 ).
Before the test, the TSW conforming to the SRS specification was simulated and calculated based on GA. e parameters of GA are as follows.Population size is 500.Control population size is 30.Initial crossing probability is 0.7.Initial variation probability is 0.3.e maximum number of generation is 150.e simulated TSW and its SRS are shown in Figure 15.We can see that, in Figure 15(a), the maximum value of the simulated TSW is within the capacity of the shaker, and in Figure 15(b), the corresponding SRS is also within the tolerance of the standard SRS.
Next, the vibration controller to drive the shaker is used so that the actual output of the shaker confirms with the simulated TSW. Figure 16 shows the input and the actual output of the shaker.Due to the capacity limitation, there exists some error between the desired simulated TSW and the actual output of the shaker.
e final results are shown in Figure 17, where Figure 17(a) shows the actual output of the shaker and Figure 17(b) shows the corresponding SRS.We can find that there are some deviations between the actual SRS and the standard spectrum, especially at high frequencies.We guess that this may be due to the poor performance at high frequencies of the shaker.Shock and Vibration

Conclusion
In this paper, we proposed a novel intelligent time-domain shock waveform (TSW) parameters matching method for the shock response spectrum (SRS) test.e basic idea is to synthesis a TSW using a group of wavelets with different frequencies, amplitudes, and phases.However, it is challenging to find an optimal group of wavelets since there are hundreds of optimization variables and objective functions.We introduce a variable-weight method to the wavelet     Shock and Vibration synthesis to combine the effects of hundreds of objective functions so that the optimization problem can be solved by using the genetic algorithm.We also take carefully the consideration of the frequency response function of the shaking system into the calculation of the objective function, making the synthesized TSW being able to be applied to the SRS test directly.To accelerate the convergence of the optimization, we introduce normalized control factors to the optimization variables so that they can be tuned in a reasonable small range.Finally, we carry out experimental validation of the proposed algorithm in two active vibration control systems, i.e., system with 500 N minishaker and system with 1-ton shaking table.It can be seen from the experimental results that, under the allowable load of the shaking system, the proposed method can accurately synthesize the input TSW, where the output TSW and SRS of the shaker can accurately meet the predefined time-and frequency-domain specification.

Figure 3 :
Figure 3: e dynamic model of the SDOF system.

Figure 4 :
Figure 4: ree numerical calculation methods to calculate the error of impact response.

Figure 5 :
Figure 5: Performance curves of three weighted fitness functions.

Figure 16 :
Figure 16: (a) Input and (b) output of the shaker.