Research on the Influence of Multiple Parameters on the Responses of a B-type Subway Train

To obtain improved comprehensive crashworthiness criteria for a B-type subway train, the influence laws of the vehicle design collision weight M and empty stroke D on the train’s collision responses were investigated, and multi-objective optimization and decision-making were performed to minimize TS (total compression displacement along the moving train) and TAMA (the overall mean acceleration along the moving train). Firstly, a one-dimensional train collision dynamics model was established and verified by comparing with the results of the finite element model. Secondly, based on the dynamics model, the influence laws of M and D on the collision responses, such as the energy-absorbing devices’ displacements and absorbed energy, vehicles’ velocity and acceleration, TS , TAMA and the coupling correlation effect were investigated. Then, surrogate models for TS and TAMA were developed using the optimal Latin hypercube method (OLHD) and response surface method (RSM), and multi-objective optimization was conducted using the particle swarm optimization algorithm method (MPOSO). Finally, the entropy method was used to obtain the weight coefficients for TS and TAMA , and multi-objective decision-making was performed. The results indicate that D and M significantly affect the compression displacements and energy absorption of the first three collision interfaces, but have limited impact on the last three collision interfaces. The velocity versus time curves of vehicle M1 and M2 are shifted and parallel with different D . However, the velocity versus time curves of all the vehicles are shifted but gradually divergent with different M . The maximum collision instantaneous accelerations of the vehicles are directly determined by M , but are only slightly affected by D . Under the coupling effect, all concerned collision responses are strongly correlated with M ; however, the responses are weakly correlated with D except for the compression displacement at the M2-M3 collision interface and the maximum collision instantaneous acceleration of vehicle M2. The comprehensive crashworthiness criteria of the B-type subway train were significantly improved after multi-objective optimization and decision-making. The research provides more theoretical and engineering application references for the subway train crashworthiness design.


Introduction
The crashworthiness of a train is the last barrier to protect the lives of the crew and passengers, and this assessment has become a mandatory requirement in Europe and North America [1][2][3]. Owing to the characteristics of long marshaling and heavy weight, the enormous kinetic energy transferred during a train collision causes plastic deformation of the vehicle structure and at the same time generates impact acceleration, resulting in injury or death to the occupants [4]. Maintaining as much structural integrity as possible after the first collision of a train or vehicle, reducing impact acceleration and mitigating secondary collisions of occupants and the vehicle's internal structure are important measures to ensure the safety of occupants and passengers [5,6].
The research methods used to assess a train's crashworthiness can be divided into theoretical analysis, experimental studies and numerical simulation. Lu [7] established train collision equations based on linear theory assumptions, and studied the influence laws of different vehicles marshaling on the first peak force of the train collision interfaces. According to the simulation results of five different projects, theoretical equations for energy absorption requirements at the front and middle ends of the train were presented. While the theoretical analysis method is usually suitable for simple and linear models, it is very difficult to solve large deformation and nonlinear problems. The experimental studies method can most realistically reflect the crashworthiness characteristics of structures and is favored by many researchers [8][9][10][11], but it is expensive and has poor repeatability. The experimental studies method is usually only employed for energy-absorbing devices or small structures, and is not a realistic option for studying the crashworthiness of the entire train. Numerical simulations can compensate for the shortcomings of theoretical analysis and experimental research, and with the maturity of commercial software, this method has been rapidly developed and applied in recent years. Numerical simulation methods include the finite element method and the multi-body dynamics method [4]. Xie et al. [12,13] investigated the collision responses of two colliding trains and vehicles impacting obstacles using the finite element method respectively. Compared with the finite element method, the multi-body dynamics method can significantly shorten the calculation time, and the calculation accuracy is within the acceptable range in engineering. The multi-body dynamics method has been widely used to study train collision laws and parameter optimization design [14][15][16][17][18].
A great deal of research has been carried out by domestic and international scholars around the crashworthiness of trains. The research conducted in different fields can be summarized as follows. At a global level, the collision response and energy management of the marshaling train have been studied, including the energyabsorption structure distribution and energy-absorbing order optimization [19]; furthermore, cab and interior structure optimization, and occupant secondary collision damage assessments have been conducted [20][21][22][23][24]. At a local level, based on the maximum energy absorption guideline, special energy-absorption devices have been designed to optimize and improve the energy-absorbing performance [25][26][27][28][29][30].
The vehicle design collision weight M and train collision speed V are two key factors that significantly affect the train collision behavior [7]. However, only few related studies in the literature have investigated their influence laws on the train collision response. Meanwhile, it has been observed in Ref. [31] that the empty stroke D (the running stroke of the train after the semi-automatic coupler shearing off but before the anti-climber devices contacting) has a clear impact on the compression stroke of the energy-absorbing devices at different collision interfaces. In this study, a typical six-car marshaling B-type subway train is considered as the research object, and a one-dimensional train collision dynamics model is established. Next, the train's comprehensive crashworthiness criteria of TS (total compression displacement along the moving train) and TAMA (the overall mean acceleration along the moving train) are defined. Then, the influence laws of D and M on the train collision response and the correlation of the coupling effect under the standard collision condition (V = 25 km/h) [1] are investigated. Based on the optimal Latin hypercube method (OLHD) and the response surface method (RSM), surrogate models for TS and TAMA are established, and multi-objective optimization is performed using the particle swarm optimization algorithm (MPOSO). Furthermore, the entropy method is used to obtain the weight coefficients for TS and TAMA and multi-objective decision-making is performed.

Typical B-type Subway Train and Energy-absorbing Devices
The layout of typical energy-absorbing devices of a B-type subway train consisting of six cars is shown in Figure 1. Considering three cars as a basic unit, the overall arrangement is symmetrical. The semi-automatic coupler and anti-climbing energy-absorption devices (anti-climbers) are located at the front end of the lead vehicle. At the collision interfaces between the head vehicle and second vehicle, a semi-permanent coupler A-B is arranged. The semi-permanent coupler C-D is located at the other collision interfaces. The anti-climbing energy-absorption devices are metal cutting-style tubes. The semi-automatic and semi-permanent couplers are combined structures constituting of a rubber buffer and expandable crush tubes. Figure 2 shows the structure and mechanical characteristics of the energy-absorption devices. All couplers have the same coupler buffer performance with a compression stroke of 73 mm and maximum resistance force of 800 kN. The semi-automatic coupler is equipped with a shear protection device, and the shearing force is 1200 kN.
After the semi-automatic coupler has been sheared off, the train runs with 22 mm of empty travel, after which

Model Building
A one-dimensional train collision dynamics model was formulated using the LS-DYNA software, as shown in Figure 3(a). M1 to M6 and S1 to S6 represent the first to the sixth vehicle of the moving and stationary trains respectively. The train collision dynamics model only considers the degree of freedom in the train's direction of movement.
The modeling method used for the one-dimensional train collision dynamics model is listed in detail in Table 1. The weight of the carbody and bogie is equivalent to the same mass block with #MAT 20 material, and only the degree of freedom in the moving direction of the train is retained. The weight of the vehicle can be changed by adjusting the material's density. The anticlimbing energy-absorption devices, semi-automatic couplers and semi-permanent couplers are simulated by discrete beams with #MAT S08 material.
In the actual collision of a rail vehicle, the contact between the wheels and the track is a complex contact process dominated by rolling friction and supplemented by the Coulomb sliding friction. The rolling friction force is much small, but due to the existence of Coulomb sliding friction caused by the deformation of the wheel or the track, there is still a certain Coulomb sliding friction force between the wheel and rail when the wheels roll on the track. To bring the simulation model closer to the actual vehicle situation, it is necessary to introduce an equivalent Coulomb sliding friction coefficient in the wheel-rail contact.
It is very difficult to accurately measure the Coulomb sliding friction coefficient and the force between the wheel and the rail. In addition, the Coulomb sliding friction force and the coefficient are influenced by the condition of the track and environmental factors, e.g. if the track is covered with rain and snow, the equivalent Coulomb sliding friction coefficient becomes smaller. The equivalent Coulomb sliding friction coefficient between wheel and rail is considered as 0.008 in Ref. [19].
In this study, the wheel-rail contact force is simulated by discrete beams with #MAT S03 material and the equivalent friction coefficient is also considered as 0.008 referring to Ref. [19].

Model Validation
The BS EN 15227:2020 (E) standard [1] clearly clarifies that evaluating the behavior of a complete train by physical testing is impractical, but it could be validated by dynamic simulation. Also, in Refs. [16] and [19], FE  The train collision finite element model was established using the LS-DYNA software, as showed in Figure 3(b), and the detailed modeling method used for the FE model is shown in Table 1. The wheel-rail contact equivalent friction coefficient is also considered as 0.008 [19].
Different from the wheel-rail contacting, the contacts between the two colliding trains and the adjacent vehicles in the marshalling train during the collision process are completely the Coulomb sliding contacts, and there must be Coulomb sliding friction forces. The materials used for the vehicles are stainless steel. Ref. [32] gives the Coulomb sliding friction coefficients of common used materials and the equivalent Coulomb sliding friction coefficient between the steel and the steel is 0.15 [32].
In order to better simulate the actual collision situation, the contacts between the impacting trains and adjacent vehicles are added in the three-dimensional train collision simulation model, and the equivalent Coulomb sliding friction coefficients are introduced to simulate the contact friction forces. Referring to Ref. [32], the equivalent Coulomb sliding friction coefficients are also taken as 0.15.
In Ref. [33], the accuracy of the one-dimensional train collision dynamics model has been verified by comparing the simulation results with the energy absorbing devices' test results and the three-dimensional train collision FE model's simulation results including the compression stroke of each collision interface, the train collision energy conversion and the vehicles' velocity versus time curves. In this study, to further verify the accuracy of the onedimensional train collision dynamics model, numerical calculations were performed on both the train collision models using the same boundary conditions, based on the two collision scenarios including different collision weight and collision speed, as shown in Table 2. The compression displacement of each collision interface, the energy collision conversion curves, and vehicles' velocity versus time curves were compared between the two models.
The comparison result of the train collision energy conversion between the collision models based on the collision scenario 1 is presented in Figure 4(a), based on the collision scenario 2 is presented in Figure 4(d). From Figure 4(a) and (d), it can be seen that in the train collision process, the kinetic energy decreases gradually and converts into the internal energy, the total energy keeps the same and finally the internal energy is about half of the initial kinetic energy. For both collision scenarios, the collision energy conversion curves between the train collision dynamics model and the finite element model are very consistent and the energy change trends and times are almost identical.
The comparison of the compression displacements of each collision interface for the two models in both two collision scenarios are shown in Figure 4(b) and (e). And it can be seen that the compression displacements of each collision interface of the two models have the same distribution, the S1-M1 collision interface between the two trains has the largest compression displacement, while the compression displacements of the other collision interfaces are symmetrically distributed around the S1-M1 collision interface and gradually decrease towards the two ends. The kinetic energy of the train collision is mainly absorbed by the energy absorbers located at the first three collision interfaces. The compression displacements of the two models are slightly different at the S1-M1, S2-S3 and M2-M3 collision interfaces, which may be due to the elastic energy absorption of car-bodies in the FE model, while the other collision interfaces are essentially the same.
The comparison of the vehicles' velocity versus time curves of the two models in both two collision scenarios are shown in Figure 4(c) and (f ). And it can be seen that similar to the energy conversion curves, the vehicles' velocity versus time curves of the two models also exhibit the same trend and balance time. After balance, the final vehicles' velocity is approximately half the initial velocity. Compared to the dynamics model, the final vehicles' velocity of the FE model is slightly lower, and the curves have higher fluctuates, which can be mainly attributed to the elastic energy absorption of the vehicles. In general, the two models show good agreement.
In summary, by comparing the dynamics model and the FE model from above two collision scenarios and three aspects, the one-dimensional train collision dynamics model always shows good consistency with the FE model. It can be proved that the calculation results of the onedimensional train collision dynamics model are accurate and reliable, and the proposed model can be used as a basic model for further research.

Vehicle Design Parameters and Train
Crashworthiness Criteria

Definition of Vehicle Design Parameters
It is assumed that the design collision weight (vehicle weight in work order plus half weight of the seated passengers) of each vehicle in the B-type train is the same, which is M. The vehicle design collision weight M, the train collision speed V, and the empty stroke D between the semi-automatic coupler and the anti-climber all have a significant influence on the train's collision response. The standard collision speed of a subway train is 25 km/h [1]. The vehicle design collision weight M and the empty stroke D are inherent design parameters of the vehicle and are not restricted by the standard. Thus, M and D The weight of the Tc car is 46.56 t, the weight of the Mp/M car is 50.24 t, and the total weight of the six-car marshaling train is 294.08 t.

Collision velocity 15 km/h
Wheel-rail friction The wheel/rail friction coefficient is set as 0.008 [19].

Definition of Train Crashworthiness Evaluation Criteria
In Ref. [19], two comprehensive crashworthiness evaluation criteria have been proposed for the entire train, namely TS (total crash displacement along the train) and TMA (total mean acceleration along the train). The smaller is the TS, the higher is the integrity of the occupant's survival space; the smaller is the TMA, the lower is the train collision acceleration.
In this study, due to the symmetry of the collision response between the moving and stationary trains, only the moving train was considered. Furthermore, the maximum collision instantaneous acceleration was considered as the vehicle collision acceleration response. In this study, based on Ref. [19], TS and TAMA were selected as the train's crashworthiness evaluation criteria, and the definitions of TS and TAMA are shown in Eqs. (1) and (2).
The total compression displacement along the moving train, TS, is defined as: where n is the number of collision interfaces of the moving train, and S i is the compression displacement of the energy-absorbing device at the collision interface j of the moving train.
The overall mean acceleration along the moving train, TAMA, can be defined as: where MA i is the maximum collision instantaneous acceleration of the vehicle i in the moving train, and m is the number of vehicles of the moving train.

Influence of the Empty Stroke D on Train Collision Response
The vehicle design collision weight M was fixed at 33 t, and the empty stroke D was taken nine points at equal intervals within the defined range. Under the standard subway train collision scenario [1], based on the numerical simulation results of the one-dimensional train collision dynamics model, the influence law of the empty stroke D on the train collision response was investigated.

Influence on the Response of the Energy-absorbing Device
Under different empty stroke D values, the compression displacement and energy absorption versus time curves of the energy-absorbing devices at each collision interface along the moving train are shown in Figure 5. And it can be seen that with increase in D, the compression displacement and energy absorption at the collision interface S1-M1 (anti-climbers) and M1-M2 generally show an upward trend, but are not strictly monotonous.
In contrast, the compression displacement and energy absorption at the M2-M3 collision interface gradually decrease. The compression displacement and energy absorption at the other three collision interfaces change slightly, and the curves are nearly horizontal. Table 3 lists the compression displacement and energy absorption of the energy-absorbing devices located at the first three collision interfaces under different D values. It can be concluded from Table 3 that the semi-automatic coupler is always fully compressed and sheared off at the initial velocity of 25 km. With increase in D from 20 to 50 mm, compared to the initial value, the relative changes in the compression displacement at the first three collision interfaces are 20.59%, 4.03% and − 25.41% respectively; furthermore, the relative changes in the energy absorption are 9.60%, 4.48% and − 27.25% respectively.
While the empty stroke D has a significant impact on the compression displacement and energy absorption of the first three collision interfaces, it has little effect on the latter three collision interfaces. In general, D has a positive effect on the S1-M1 and M1-M2 collision interfaces, while it has a negative effect on the M2-M3 collision interface. Among the first three collision interfaces, D has the greatest effect on the M2-M3 collision interface, followed by S1-M1, and has the least effect on the M1-M2 interface. The empty stroke D exhibits the same influence laws on the compression displacement and energy absorption at different train collision interfaces.

Influence on the Vehicle Response
Under different empty stroke D, the velocity versus time curves of each vehicle of the moving train are shown in Figure 6. It can be seen from the figure that with increase in D, the velocity versus time curves of vehicles M1 and M2 begin to shift at 0.2 s. The velocity profiles at different D are shifted to parallel each other, and the time at which the final velocity reaches balance is also shifted. In contrast, the velocity curves of vehicles M3 to M6 remain the same under different D. With increase in D from 20 to 50 mm, the velocity balance times of vehicles M1 and M2 increase from 0.396 to 0.433 s and 0.556 to 0.595 s respectively. Table 4 lists the maximum collision instantaneous acceleration of each vehicle of the moving train under different empty stroke D. It can be concluded that with increase in D, the maximum collision instantaneous acceleration of vehicles M1 to M3 remains the same; the maximum collision instantaneous acceleration of vehicles M4 to M6 shows a minor change compared to the initial value, and the relative changes are − 0.47%, − 1.16% and − 0.11% respectively.
The maximum collision instantaneous acceleration of a vehicle is not directly affected by the empty stroke D, which is completely determined by the compression force of the energy-absorbing devices located at both ends of the vehicle and the vehicle's weight. A slight change in the maximum collision instantaneous acceleration of vehicles M3 to M6 is caused by the different compression displacements of the energy-absorbing devices located at both ends of the vehicle, which is within the range of the coupler buffer, and different compression displacements produce different compression forces.

Influence on the Responses of TS and TAMA
The collision responses of TS and TAMA under different D along the moving train are listed in Table 5. It is known that both TS and TAMA show a downward trend with increase in D. With increase in D from 20 to 50 mm, compared to the initial value, the relative changes in TS and TAMA are − 2.11% and − 0.22% respectively.
The empty stroke D has a negative impact on TS and TAMA. Similarly, D has no direct effect on TAMA. The decrease in TAMA can also be attributed to the slight Table 3 Collision responses of the energy-absorbing devices under different empty stroke D

Influence of the Vehicle Design Collision Weight M on Train Collision Response
The empty stroke D was fixed at 35 mm, and the vehicle design collision weight M was taken nine points at equal intervals within the defined range. Under the standard subway train collision scenario [1], based on  the numerical simulation results of the one-dimensional train collision dynamics model, the influence laws of the vehicle design collision weight M on the train collision response were investigated.

Influence on the Response of the Energy-absorbing Device
Under different vehicle design collision weights M, the compression displacement and energy absorption versus time curves of the energy-absorbing devices at each collision interface along the moving train are shown in Figure 7. And it can be observed that with an increase in M, the compression displacement and energy absorption at the collision interfaces S1-M1, M1-M2, and M2-M3 increase in general but are not strictly monotonous. Similar to the situation under different D, the compression stroke and energy absorption at the other three collision interfaces change slightly, and the curves are almost horizontal. Table 6 lists the compression displacement and energy absorption of the energy-absorbing devices located at the first three collision interfaces under different M. It can be seen from Table 6 that, in a similar manner, the semi-automatic coupler is always completely compressed and sheared off at the initial velocity of 25 km/h. With increase in M from 30 to 36 t, compared with the initial value, the relative changes in the compression displacement at the first three collision interfaces are 42.17%, 19.43% and 24.32% respectively; the relative changes in the energy absorption are 18.78%, 22.49% and 26.57% respectively.
The vehicle design collision weight M has a significant impact on the compression displacement and energy absorption of the first three collision interfaces, but has little effect on the latter three collision interfaces. In general, M has a positive effect on all the first three collision interfaces, in terms of the compression displacement and energy absorption. Considering the compression displacement, the impact on the S1-M1 collision interface is the greatest, followed by M2-M3, with the smallest impact on M1-M2. With respect to energy absorption, M has the greatest effect on the M2-M3 collision interface, followed by M1-M2, with the least effect on S1-M1.     Table 7 lists the maximum collision instantaneous accelerations of each vehicle in the moving train under different M. It can be seen that with increase in M from 30 to 36 t, the maximum collision instantaneous accelerations of all vehicles decrease linearly, and compared with the initial value, the relative changes are all around − 16%.
The vehicle design collision weight M directly determines the vehicles' maximum collision instantaneous acceleration. An increase in M leads to an extension of the balance time of the vehicle velocity, and therefore directly leads to a decrease in the average acceleration of the vehicle in the train collision.

Influence on the Responses of TS and TAMA
The collision responses of TS and TAMA under different M along the moving train are listed in Table 8. It is known that with increase in M, TS gradually increases, while TAMA gradually decreases. With an increase in M from 30 to 36 t, TS increases from 1248.459 to 1459.425 mm, indicating a relative change of 16.90%; TAMA decreases from 1.4422 to 1.2122 g, indicating a relative change of − 15.95%. M has a positive effect on the TS, while it has a negative effect on the TAMA.

Coupling Influence and Correlation Analysis of the Design Parameters
To investigate the coupling influence and correlation of D and M on the train collision responses, the OLHD method was applied to the experimental design within the defined range. Both D and M were selected as integers. The number of experiments was 10. In a standard subway train collision scenario (25 km/h), considering the compression displacements of the first three collision interfaces, the maximum collision instantaneous acceleration of the first two cars, and the responses of TS and TAMA, the one-dimensional train collision dynamics model was used for numerical simulation to obtain the train collision response matrix, as shown in Table 9.
The correlation coefficient is a statistical indicator that reflects the closeness of the relationship between variables, and the correlation coefficient ranges between − 1 and 1. If the correlation coefficient is 1, this indicates that the two variables are completely linearly correlated; − 1 indicates that the two variables are completely negatively correlated; 0 indicates that the two variables are not correlated. The closer the data is to 0, the weaker the correlation.
The correlation coefficient r xy can be calculated as Eq. (3) r xy = s xy s x s y ,  where s xy represents the sample covariance; s x represents the sample standard deviation of x, and s y represents the sample standard deviation of y.
The sample covariance and sample standard deviation can be calculated as Eqs. (4)- (6).
where x i is the sample ith value of x; x is the sample average value of x; y i is the sample ith value of x; y is the sample average value of y, and n is the number of sample points.
The correlation coefficient matrix between the vehicle design parameters and train collision responses is shown in Table 10. It can be concluded that under the coupling

Formulation of the Surrogate Model
According to the train collision response matrix in Table 9, the surrogate models of TS and TAMA are  established based on the RSM, as shown in Eqs. (7) and (8).
To evaluate the accuracy of the obtained cubic polynomial surrogate models, three evaluation criteria are introduced, namely the average relative error (ARE), the root mean square error (RSME) and the fitting coefficient (R 2 ), as defined in Eqs. (9)-(11) [26,27,34,35].
where, y i is the response predicted value of the surrogate model corresponding to each input variable experiment; y i is the response value calculated from the train collision dynamics model corresponding to each input variable experiment, and y i is the average response value calculated from the train collision dynamics model corresponding to all input variable experiments.
The results of the accuracy evaluation of the surrogate models are shown in Table 11. The closer the fit coefficient R 2 is to 1, the more accurate the surrogate model will be.
The smaller the sum, the better the accuracy of the surrogate model.
The smaller the ARE and RSME , the higher accuracy of the surrogate models. It can be concluded from Table 11 that the accuracy of the cubic polynomial surrogate models established by the RSM is sufficiently high and can be used as basic models for further research.

Multi-objective Optimization
Based on the established surrogate models, under the standard subway train collision scenario (25 km/h), the particle swarm optimization algorithm (MOPSO) was applied to perform multi-objective optimization of TS and TAMA. The empty stroke D and vehicle design collision weight M can be arbitrarily selected within the defined range, and the initial collision velocity of the train is considered to be 25 km/h. From Table 7, it can be seen that the maximum collision instantaneous acceleration of each vehicle is below 5 g, which meets the requirement of the standard [1]; therefore, the acceleration is no longer restricted. According to Figure 2, the maximum compression stroke of the anti-climber in S1-M1 interface S anti is 265 mm, and the maximum compression stroke of the energyabsorbing devices at the other collision interfaces S i (i = 1-5) is 418 mm. Based on the above objectives and constraints, the multi-objective optimization problem of the comprehensive train crashworthiness evaluation criteria TS and TAMA can be defined as Eq. (12).
By optimizing and solving the multi-objective problem, the Pareto frontier optimal set of TS and TAMA can be obtained, as showed in Figure 9. And it can be seen that as TS decreasing, TAMA gradually increases. In contrast, as TAMA decreasing, TS gradually increases. The minimum TS can be derived from one end of the optimal set, while the minimum TAMA can be gained from the other end. TS and TAMA constitute a pair of relatively contradictory goals, which are consistent with the conclusion obtained in Ref. [19]. The values of TS and TAMA need to be determined according to the actual engineering requirements.

Multi-objective Decision Making
The entropy method is an objective weighting method. The entropy value can be used to measure the amount of data information contained in the evaluation criterion. The greater the data information contained in the (12)  evaluation criterion, the smaller the entropy value [36,37].
To minimize the influence on the evaluation results owing to the different dimensions and magnitude of the evaluation criteria, it is necessary to conduct dimensionless processing for each evaluation criterion [36]. In this study, a standardized approach [37] is employed, as defined in Eq. (13).
where x j (j = 1, 2, 3, · · · , m) is the average value of the observation of the item index j; s j (j = 1, 2, 3, · · · , m) is the mean squared error of the observation of the item index j; x ij and x * ij are the original and observation values obtained after standardized processing of the evaluation index j for the evaluated object i.
After standardized processing, the data may contain both positive and negative values, which are inconvenient for data processing. The T score function should be introduced to modify the standardized data [36,38], as shown in Eq. (14).
where x T ij is the transformed data. The weighting factor calculation for the entropy method mainly consists of the following steps [37,38]: (1) Calculate the feature weight (13) x * ij = x ij − x j s j , where P ij is the characteristic proportion of the evaluated object (scheme) i under the evaluation criterion j.
(2) Calculate the entropy value of the evaluation criterion j (3) Calculate the difference coefficient of the evaluation criterion j (4) Calculate the weight coefficient of the evaluation criterion j where the weight coefficient w j is the weight of each evaluation criterion after normalization. Using the entropy method, the weight coefficients of the train crashworthiness criteria are obtained, as shown in Table 12. The weight coefficients of TS and TAMA are 0.50003 and 0.49997 respectively.
The respective weight coefficients are assigned to TS and TAMA, and multi-objective optimization is performed by applying the MPOSO method to obtain the optimal solution, as shown in Table 13. It is known that when the empty stroke D is 50 mm, and the vehicle design collision weight M is 32.4 t, the optimal solution can be obtained. Here, TS is 1307.5 mm, and TAMA is 1.4 g.
Based on the same design parameters, the one-dimensional train collision dynamics model is used for the numerical calculations to obtain TS and TAMA. Compared with the predictions of the surrogate models, the (16) e j = − n i=1 p ij ln(p ij ).