Determination of Optimal Drop Height in Free-Fall Shock Test Using Regression Analysis and BackPropagation Neural Network

The primary purpose of this study is to provide methods that can be used to determine the most suitable drop height for shock testing military equipment, in an efficient and cost ineffective manner. Shock testing is widely employed to assess the performance of electronic systems, including military devices and civilian systems. Determining the height of the drop for a test item is an important step prior to performing the test. Dropping a test item from excessive height leads high G-peak values to damage the test equipment. On the other hand, dropping an item from a low height may not generate the required G-peak value and duration. Therefore, prior to performing shock tests, an optimal drop height must be established to ensure that the resulting G-peak value and duration timematch the required test values.The traditional trial-and-error methods are time-consuming and cost-ineffective, often requiring many attempts. To improve the conventional approaches, this study proposes the application of regression analysis and back-propagation neural network for determining the most suitable drop height for free-fall shock tests. A new method is suggested for determining the drop test height.The results of the model are verified, using the results of a series of experiments that are conducted to verify the accuracy of the suggested approaches.The results of the study indicate that both approaches are equally effective in providing an effective guideline for performing drop tests from heights that would result in the peak Gs and duration needed for testing electronic devices.


Introduction
Shock testing has an increasingly important role in ensuring that military electronic devices meet the requirements for effective and reliable operation in harsh environments.Shock specifications of military electronic devices are generally expressed in terms of a simple acceleration pulse, such as a half-sine wave or a sawtooth wave in milliseconds, to simulate the shocks that the devices may experience in military environments.A shock test may include an assessment of the overall system integrity for safety purposes during handling, transportation, or use.The provision and application of shock simulation methods present many problems for laboratorybased teams [1].Shock loading problems cover a very wide range for physical parameters and are associated with many complex types of damage and malfunction.Shock tests will also be conducted to assess the life time from physical and functional performance of a device.Guo and Zhang [2] utilized the vibration transmissibility of a honeycomb structure which is applied regression analysis to design a shock absorber.Various methods have also designed shock waveforms and the conical lead pellet; however, these studies made the determination of the drop height in a shock machine disappear [3][4][5].In the conventional trial-and-error method, a test item is placed on the shock machine table and then dropped from different heights to obtain various G-peak values and durations until the required shock specifications are met.The required shock specifications are obtained from  measuring the shock to a device.A device must be tested to verify its environmental adaptability to serve its owner.The test item was placed on the impact table of a free-fall shock test machine, as shown in Figure 1 and Figure 2 shows a typical shock waveform in the time domain.The dotted line represents the actual waveform resulting from the free-fall shock test item impacting the elastomers with an imperfect rebound [6]; the solid line indicates the specifications defined based on the peak Gs and duration of a half-sine waveform.The duration defines that a test item makes contact with the top programmer until the time it breaks contact.The shock event is often simulated using a dummy test item to adjust the height of the shock machine and determine the type of programmers needed via repetitive trials before a formal shock test is performed with the free-fall shock machine in the laboratory.This method is time-consuming and costineffective as many attempts are needed and elastomers are consumed before the optimal drop height is determined.Thus far, a simple, convenient method is not suggested for estimating the optimal drop height.Regression analysis and a back-propagation neural network (BPNN) are reliable estimating approaches for determining the optimal drop height to assure the shock test quality in this crucial problem.

Shock Motion
The purpose of the shock test is to obtain a waveform to match the solid line in Figure 2. The solid line in Figure 2 represents the acceleration history measured by an accelerometer fixed to the table.The relationship between the shock specifications (peak Gs and duration), drop height, and the programmer are elucidated in a shock test.Shock motion analysis is applied to use the dynamical model depicted in Figure 3.
In Figure 3,  is the total mass, including that of the impact table, test item, and fixture; () is the impact table displacement; and  is the drop height from the top of the programmer to the bottom of the impact table.
The motion of the test item can be divided into two stages: before impact and after impact.According to Newton's second law the equation for motion in the preimpact stage can be expressed as where  is gravity.The downward velocity is defined as negative and is derived as In the postimpact stage, as the table impacts the programmer, the programmer will rebound with a reactive force  to the impact table.Therefore, the equation for the motion of the test item during the impact can be written as A detailed analysis of the reactive force is complex as the reactive force varies depending on the mechanical properties of the programmer material.In many cases, the reactive force is nonlinear result in computing the exact displacement, (), which is difficult.A shock test is designed so as to verify the specifications, and an appropriate programmer is used to generate a waveform in which the peak Gs and duration mimic the solid line and match the shape in Figure 2. When the table impacts the programmer, its acceleration is given as (4) in which the initial condition of zero acceleration (gravity is neglected) is assumed.Consider the following: The initial velocity of ( 4) is found based on (2); that is, where   is the impact velocity of the test item and programmer at the instant of impact.Equation ( 4) can be integrated and the result is substituted into (5) to yield Let   be the velocity at time  =  as the test item breaks contact with the programmer.Let  =  be substituted into (6), and   can be derived as in The area under the curve of acceleration is equal to the change in velocity between   and   , as in The function of the programmer is to act like a spring.Let the elastic coefficient of the programmer be given by .The formula for the conservation of energy is given below: where  is the velocity of the test item and  is the displacement of the test item;  = 0 implies that the test item is in contact with the programmer.The force is defined such that the upward direction is positive.Equation (9) can be differentiated with respect to time to obtain Finally, (10) can be solved to yield Let  = √/; thus, the solution to (11) is From ( 12) and ( 5), the acceleration can be rewritten as where  and  are the magnitude and phase of the acceleration, respectively.Consider ) . ( When √2 is much larger than ,  and  can simplified as follows: From ( 13)-( 15), the acceleration can be represented as in where The peak Gs is determined from  and , and duration is determined from .

Regression Analysis and Back-Propagation
Neural Network (BPNN) 3.1.Regression Analysis.Regression analysis is used to fit the curve of the relationship between the input and output database.In this study, polynomial analysis and the ration linear polynomial regression analysis were performed on the data sets.The polynomial is denoted as follows: where  is the input data;  is the output data; and  0 ,  1 ,  2 , . . .,  −1 ,   are the parameters to be determined.All data can be reformulated into an equation as follows: where  = 1, . . .,  for different cases.
The square error can be written as where  = 1, . . .,  for different cases and  is the polynomial order.
To minimize the total square error , one must set the partial differential of  to  0 ,  1 ,  2 , . . .,   at zero: The optimal solutions of  0 ,  1 ,  2 , . . .,  −1 ,   can be derived when the data sets are individually input into (20) to obtain the following matrix equation for the different cases: or where  is the polynomial order.
It is clear that (23) has  unknown variables over all cases.Substituting (21) into (20) yields the square error as follows: Another type is the rational linear polynomial.The rational linear polynomial is often used to obtain minimum errors in many cases.This mathematical expression can be written as where  is the case number;   ,   , and   are the parameters;   is the output value; and   is the input value of the case .

BPNN.
where   is the desired output and   is the actual output.The partial derivative ((  )/(  )) determines the weight adjustment characterizing the gradient descendent of the algorithm.This partial derivation measures the weight   contribution to the BPNN error function for pattern .If this derivation is positive, the error amount is considered to be increasing and the weight   should be reduced in order to decrease the difference between the actual output and the desired output.In the case of the partial derivation being negative, the weight   contributes to the generation of an output smaller than the desired output (and, therefore, to reduce this difference, the weight should be increased).Equation (27) defines the way in which the connection weights in the network are adjusted: where   is the weight of the connection between neurons  and , and  is the learning rate.
The derivation was shown in (27) to provide a generic rule used to adjust the weight of the connections.Equation (28) illustrates this generic rule of weight adjustment below: where   = Σ    if neuron  is a hidden layer unit.The choice of an appropriate learning parameter  will significantly influence the convergence rate of the algorithm.If this parameter is too small, it may take several steps to reach an acceptable convergence.On the other hand, a high learning rate can lead to oscillations that could impede convergence.A possible solution is to use a momentum term  that would quantify the importance of the weight variation at a step.Equation (29) illustrates the weight adjustment rule with a momentum term: The momentum term can make the learning algorithm more stable and accelerate the convergence in flat regions of the error functions.However, determining the optimal value of the momentum parameter  introduces the same difficulties as encountered in the case of the learning rate parameter.
This study applies BPNNs to determine the drop height and duration for the shock test with programmers.
The structure of the BPNN (Figure 4) is as follows.
(1) The BPNN has an input layer expressing the input parameters in the network.The number of neurons depends on the problem's complexity.The input layer is the first (bottom) layer in the structure.(2) Above the bottom layer is a middle hidden layer expressing the interaction between input parameters and processed neurons.The number of neurons cannot be precisely determined.The number of neurons is typically determined based on when the optimal result is obtained.
(3) The third (top) layer of the BPNN is an output layer which denotes the network output.The number of neurons is also determined based on the problem's complexity.
(4) The BPNN includes a transfer function.The sigmoid function is chosen as the nonlinear transfer function; it is expressed as follows: The following procedure is used to analyze the BPNN: (1) determine the number of neurons at each layer; (2) set the initial weights and bias values in the network randomly; (3) insert input and output vectors into the network for training the weights; (4) estimate the output values of the hidden and output layers; (5) calculate the difference in output values between the hidden layer and output layer; (6) establish the adjustment coefficients for weights and bias values; (7) update bias and weight values; (8) repeat steps (3)-( 7) until the network converges.

Proposed Scheme
The proposed prediction scheme is composed of data collected under different conditions.The peak Gs and duration of the shock test are determined for different programmers.
From previous (and new) experiments, the data sets are arranged so as to fit the optimal curve by regression analysis; the BPNN is also trained to develop a knowledge database until the learning structure is robust.The data sets for four cases of different programmers are gathered to train and recall the target using either a BPNN or regression analysis.Finally, the degree of accuracy of targets is compared for either BPNN or regression analysis in various situations.

Data Collected and Error Evaluating
Processes.Before collecting data sets, all programmers and free-fall shock machines are calibrated for testing.Figure 5 shows the complete experimental setup of the free-fall shock machine.Figures 6 and 7 show photographs of different programmers.In the experiment, the same test item, type, and number of programmers were used to obtain data sets at seven heights (10, 20, 40, 60, 80, 100, and 120 cm); each test was run twice at each height to obtain peak Gs and durations that were close to the specifications.The experiment was run 28 times to ensure that the 28 data sets collected were within tolerances.Bad data sets were eliminated as necessary during data processing.Tables 1 and 2 present all 28 data sets which indicate the peak Gs and durations . Figure 8 shows the four representative shock-test results.The data sets were used to verify the effectiveness of the proposed method for determining the optimal drop height and duration.Tables 3 and 4 also showed the resultant mean square errors (MSEs).The MSE is defined in (31) and the maximum error (ME) is also defined in (32): where   and   are the target height and duration, respectively;   and   are, respectively, the height and duration estimated by either regression analysis or BPNN;  is the test number; and  is the total number of data sets in each case.

Regression Analysis Application.
The 24 data sets were selected for drop heights of 10, 20, 40, 60, 80, 100, and 120 cm (Table 1), and the remaining four data sets were used for verification of the results via regression analysis (Table 2).The relationship between  and  in cases 1 and 3 used a linear polynomial for estimating duration as (23); however, the relationship between  and  in cases 2 and 4 used rational linear polynomial for estimating duration as (25).The optimal parameters for ( 25 In Figures 11 and 12, the solid line and black dotted line are the fitting curves for cases 2 and 4, respectively, as also determined by regression analysis.The circles in these figures are the fitting data sets; the crosses represent the test data sets.Tables 3 and 4 show the MSEs of height and duration values.
According to the estimated results, the average MSE heights of the regressive patterns and test patterns are 0.02218 cm and 0.9368 cm, respectively.Meanwhile, the average MSE duration for the regressive patterns and test patterns is 0.00656 ms and 0.024 ms, respectively.The maximum average error of height and average error of duration in the regressive and test patterns is individually less than 1 cm and 1 ms, respectively.Regression analysis can be used to determine the equations for drop height and duration; however, only three variables can be used to fit the curves.

BPNN Application.
In the same arrangement used for the training network in Table 1, the remaining four data sets are used for verification of the results obtained via the BPNN test (Table 2).Each selected data set includes two elements, the peak Gs and the coefficient of elasticity, .The coefficient of elasticity, , and the weight of a test item for the four cases are set as constants.The data sets are normalized with respect to        peak Gs, height, and duration to 0.2−0.8.This study uses two neurons in the input layer, 12 neurons in the hidden layer, and one neuron in the output layer.After the BPNN was trained 163 times, the network converged (Figure 13).The sudden drop means a local minimum error.Using early stopping can overcome this phenomenon and the results were not affected through real experiments verified.Tables 3 and 4 5 and 6).The four data sets were further tested to estimate the height and duration by BPNN (Figures 9−12), and Tables 3 and 4  0.6252 cm and 0.0168 ms, respectively.Thus, the BPNN can be applied successfully to determine the optimal drop height and duration for the free-fall shock test.The BPNN values have small errors as they obtained a good fitting curve in this study.Tables 3−6 summarize the conditions in the different cases and computational results.these curves indicate that some difference exists between the two methods.The computational results obtained by regression analysis and BPNN are very close to the actual drop heights and durations.A further comparison of the two approaches (Tables 3−6) indicates that the BPNN performs better than regression analysis in terms of average MSE and average maximum error.Moreover, the BPNN outperforms regression analysis based on its ability to accommodate more than three variables.Additionally, the order of precision obtained when using regression analysis is not particularly high.It is noteworthy that the test value for duration in case 2 obtained by the BPNN and regression analysis is 0.0016 ms and 0.301 ms, respectively (Table 6).Case 2 demonstrates that the test point on the curve obtained by BPNN estimated with a higher degree of accuracy than obtained by regression analysis (Figure 12) (the application of regression analysis is limited to only three variables).A summary of this study's conclusions and a performance comparison are given in Table 7.

Figure 5 :
Figure 5: The setup of the free-fall shock machine.

Figure 6 :
Figure 6: The elastomer of felt on the machine.

Figure 7 :
Figure 7: The plastic elastomer on the machine.
) are  2 = 0.0895,  2 = 1095, and  2 = 102.7 when  = 2 and  4 = 0.925,  4 = 90.82,and  4 = −43.33 when  = 4. Cases 1 and 3 and cases 2 and 4 are coupled when drawing the fitting curves.Figures 9 and 10 depict the relation between drop height  and .Figures 11 and 12 depict the relation between  and .The solid line and black dotted line are the fitting curves for cases 1 and 3 in Figures 9 and 10, respectively, as determined by regression analysis.

Figure 9 :
Figure 9: Fitting curves of  versus  in cases 1 and 3.

Figure 10 :
Figure 10: Fitting curves of  versus  in cases 2 and 4.

4. 4 .
Comparison of the Two Approaches.The fitting curves obtained by regression analysis do not completely match the fitting curves obtained by the BPNN.The curves of the BPNN are spiral on regressive curves (Figures 10−12);

Figure 13 :
Figure 13: Training of epochs by ANN.

Table 1 :
Training data sets by free-fall shock machine.

Table 2 :
Test data sets by free-fall shock machine.

Table 3 :
MSE of heights by BPNN and regression analysis.

Table 4 :
MSE of time duration by BPNN and regression analysis.

Table 5 :
Maximum errors (ME) of height by BPNN and regression analysis.

Table 6 :
Maximum errors (ME) of time duration by BPNN and regression analysis.