3.1. The Numerical Analysis of the Safety Structure
The calculation consists of two stages. In the first stage, the kinematical calculations for establishing the tilting angle when the structure becomes unstable, the time until it is falling onto the concrete surface, the angle of the impact of the frame structure against the concrete surface, and the angular velocity of the structure around the rollover axis at the moment of the impact were carried out. In the second stage, calculations of the stressed state of the deformed structure at the moment of the impact upon marginal conditions of the task using data from the results of the kinematical analysis were carried out. In this way, dynamic analysis of the safety frame was performed using realistic model boundary conditions.
The kinematical calculations were carried out by applying the annex Motion Dassault systems of SolidWorks (SolidWorks 2019 version). SolidWorks Motion is the light version of MSC ADAMS, as the world’s most famous and widely used Multibody Dynamics (MBD) software. Multibody Dynamics simulations are a powerful method to study both the kinematic and the dynamic behavior of complex systems. MBD simulations are suited to study the dynamic behavior of interconnected rigid and/or flexible bodies undergoing large translational or rotational displacements. The motion of those bodies is calculated based on applied loads and boundary conditions defined. Typically, these simulations have short calculation times, making them the preferred tool to conduct parameter studies or optimizations very efficiently.
The analysis of the impact was performed by applying the Ansys Autodyn system (Ansys Autodyn 2019R1 version) from ANSYS Inc. Ansys Autodyn simulates the response of materials to short duration severe loadings from impact, high pressure, or explosions. It is best suited for simulating large material deformation or failure. Complex physical phenomena, such as the interaction of liquids, solids, and gases, can be modeled within Autodyn. Integrated within Ansys Workbench with its own native user interface, this program has for decades led the industry in ease of use, enabling consumers to get accurate results with the least amount of time and effort.
A spatial geometrical model of the safety structure for the verification is presented in
Figure 3 below. In the calculations, several computational models were formed. The model for kinematic analysis was a spatial tubular structure of a bus safety frame, on the basis of which a natural safety frame was also produced. For dynamic (impact) analysis, a shell model was developed, which allows the significant reduction of the number of finite elements, thus shortening the computation time.
The natural bus safety frame fragment was welded using square and rectangular cross-section profiles. The thicknesses of these profiles were 2 and 3 mm. In addition, stiffening edges and V-shaped hinges were used for welding the frame. The thicknesses of these elements were 1.5, 3, and 4 mm. The computational model for kinematic analysis was created using a direct spatial model of the frame, on the basis of which the frame was designed for a natural experiment, that is, the geometry of the model corresponds exactly to the real frame model. All frame elements were fixed to each other. Distance sensors were modeled at three corners of the model frame to accurately determine the moment of contact. The frame rollover computational model is presented in
Figure 4 below. Here, 1—the concrete surface (base); 2—the tip-up platform; 3—the tip-up spatial model of the safety structure with an auxiliary stand; the width of the stand conforms to the width of the rear axle of the coach and the lower part conforms to the point of a wheel’s support on the base; 4—an additional mass of 610 kg; and 5—sensors for contact detection.
In
Figure 5, the initial position of a rollover is shown. It was established by calculations on varying the angle until the structure of the frame loses its stability and starts turning over. To determine this initial angle, the position of the center of gravity of the frame structure was calculated and the frame was tilted until this center of gravity was extended beyond the axis of the tip-up platform hinge (about which tip-up was made).
In course of the kinematical analysis, it was found that the initial angle when the structure of the frame (loaded with an additional mass of 610 kg) loses its stability and starts turning equals 40 degrees. At the said moment, the velocity and the angular velocity of the top point of the frame equals 0. An additional mass of 610 kg was loaded to bring the frame as close as possible to real operating conditions (simulating seats and passengers). The double profile was mounted at a height to match the center of gravity of the passenger bus.
In the course of the analysis of the frame rollover, it was found that the top point of the frame structure reached the base within 1.64 s. The velocity of the angle at the moment of the impact against the base was about 5.97 m/s. In
Figure 6, the alteration of the angular velocity of the safety frame around the rollover axis before the impact is shown. In the course of an analysis, it was found that the angular velocity of the frame at the moment of the impact was 129 degrees per second (2.24 rad/s). The obtained values were used in the dynamic (impact) analysis of the frame as the boundary conditions of the model.
In the second stage, an analysis of the falling of the safety frame was carried out, striving to establish the permanent deformations on a simulation of a structure rollover test according to UNECE Regulation No.66/02.
As already mentioned, a surface geometric model of the frame based on a spatial model had been developed for the analysis of finite elements. In addition, in the development of this model, the frame tubular elements were correctly interconnected, eliminating the gaps between the frame elements for welding, which was necessary in the spatial model. The elements of the surface model were of different thicknesses: 1.5, 2, 3, and 4 mm. All elements were interconnected using bonded type component contact. In addition, a three-dimensional base was modeled at an angle of 18.9 degrees to the vertical plane of the safety frame. Thus, the computational model for impact analysis was constructed using mixed (“shell” and “solid” type) finite elements (base and H-beam “solid” type finite elements, the rest of the structure was “shell” type finite elements). “No penetration” type contact with a coefficient of friction of 0.6 was used for the contact between the frame and the base.
The analysis of the safety frame with an additional mass of 610 kg was performed using a true (computational) nonlinear steel S355 material model with kinematic hardening, as shown in
Figure 7 below. The true stress–strain curve for simulation of S355 material was made according Al-Thairy et al. [
21]. The mechanical properties of concrete were used to model the base material (Young’s modulus E = 23.0 × 10
9 Pa, density 2400 kg/m
3 and Poisson’s ratio 0.24).
The chosen geometry and the marginal conditions for the model conform to the requirements of the rollover test to the maximum possible extent. The angle of the model’s frame with a plane conforms to the angle where the structure becomes unstable (found in the course of the analysis of the falling). The initial angular velocity of the safety frame was chosen to be 2.24 rad/s, as established in the course of the above-mentioned analysis of the falling.
In
Figure 8, the presented computational model is divided into shell type finite elements. In order to select the right grid size of the finite element, test calculations were performed with a partial frame model and the optimal size of the finite element was sought. The results of the calculations obtained by shredding a finite element grid were more accurate, but only up to a certain level. With a large increase in the number of finite elements, the computational error increases rapidly due to the approximate calculation method. In addition, as the size of the finite element decreases, the integration step decreases proportionally, resulting in a significant (non-linear) increase in the time to solve the problem. It consists of 37,249 nodal points and 31,782 finite elements.
In
Figure 9, the summarized frame displacement fields 0.24 s after the beginning of the impact are shown; in
Figure 10, fields of equivalent plastic deformations of the structure’s frame at the same time moment are shown. The summarized frame displacement fields make it possible to judge the suitability of the frame in the sense that the structural elements of the frame do not fall into the safety space, as shown in
Figure 2. Even if the elements of the frame structure cross this safety space as little as possible, the safety frame structure would be considered unsuitable. Equivalent plastic deformations indicate how irreversibly the frame structure will be deformed. If the equivalent deformations exceed the critical values (for steel S355 ε > 0.20), the disintegration of the structure would start. The formulation of the problem would make it possible to see such a process of structural disintegration.
Figure 10 shows the equivalent plastic deformation fields of the structure’s frame, 0.24 s after the beginning of the impact. Furthermore, we can see that the largest deformations occur in the upper corners of the frame structure, as well as on the sides in the zone of maximum overall dimension. It is these local deformation zones that will determine the plastic deformations of the frame. In order to minimize the plastic deformations of the frame in these places, the structure should be stiffened.
3.2. A Natural Experiment of the Safety Frame Falling
A natural experiment of the safety frame falling was carried out, striving to establish the permanent deformations of the frame on the rollover test of the structure according to UNECE Regulation No. 66/02.
For the experiment, a platform that conforms to provisions of UNECE Regulation No. 66/02 was used. In
Figure 11, a safety frame on a tip-up platform is shown. The initial position of the frame (the point where the experiment starts from) was chosen in a such way that the position of the frame is close to the one losing stability.
The rollover of the safety frame was shot by a Phantom v711 high-speed camera.
The safety frame rollover experiment was fixed with the frequency of 1100 shots per second (1 shot: 909.09 μs).
In an analysis of the filmed materials, it was found that the frame was falling from the initial position until its corner touched the ground for about 1.66 s, and after four rebounds with decreasing amplitudes, became stable. In
Figure 12 below, the filmed image of the safety frame 0.24 s after the beginning of the impact is shown. The number of frame rebounds allows us to judge the stiffness of the structure and the efficiency of shock absorption. In addition, this parameter is very valuable in verifying the computational model because the actual structure is welded (the frame has residual stresses) while the computational model is idealized.
In the course of the analysis of the data provided in the said materials, it was found that the maximum deformation of the frame took place 0.105 s after the beginning of the impact and it rebounded from the base 0.157 s after the beginning of the impact. The frame achieved the maximum amplitude of the first rebound 0.356 s after the beginning of the impact and reached the ground once more after 0.674 s.
After the rollover test, the frame was measured, and it was found that the maximum displacements of it were 132.46 mm. In
Figure 13a, the silhouettes of the deformed frame and the frame before deformation obtained in course of the said measurement are provided. In
Figure 13b, the deformed frame with safety space provided is shown. It is clear that the safety clearance was not compromised, and the frame was stiff enough to pass the test according to UNECE Regulation No. 66/02.
A numerical analysis of the results obtained from the safety structure rollover test and the data provided in the filmed materials of the above-mentioned natural experiment affirm that a good coincidence of the results was obtained, so the computational model should be considered correct. For example, in the natural test (the results were obtained from the filmed materials), the maximum deformation of the frame took place 0.103 s after the beginning of the impact and in the numerical experiment, the time was equal to 0.096 s; this is less than 7% difference. In the natural test, the frame rebounded from the base 0.157 s after the beginning of the impact and in the numerical experiment, the time was equal to 0.155 s; this is less than 1.2% difference. A good coincidence of the results of the calculations with the results of the natural experiment was confirmed by the permanent deformed contour of the safety frame.
Figure 14 shows juxtaposed deformed forms of the safety frames, which were obtained during a natural experiment and numerical analysis. The maximum difference in results is less than 2%.
The contour of the safety frame prior to and after the deformation appeared during the rollover test (contour before deformation—black line) and the deformed state after the rollover (numerical analysis—orange line, experimental test—black dotted line). The results of the numerical analysis of the safety structure rollover test and experimental test shows a good coincidence of the results which are very similar.