1. Introduction
At present, with the continuous development of computer hardware facilities, the calculation and analysis speed of computers has been rapidly increased. Because of the difficulty and high cost of testing heat conduction experiments in the actual process [1], computer simulation is usually used to simulate the experiments before the actual experiments, and the improvement of computer hardware performance nowadays makes the simulation more accurate and detailed [2].
This experiment is an experimental simulation of data from a thermal analysis experiment conducted on a model of the head of a projectile. Considering the complexity of the actual heat transfer situation, ANSYS/LS-DYNA [3] is used as the simulation software.
2. Simulation Experiment Preparation
2.1. Geometric Modeling
Before the simulation analysis, the geometric model of the device needs to be established, including the geometric models of the radome, heat shield, thermal protection layer, antenna, metal connection ring and other components, as shown in Figure 1. Since the structure of the device has axisymmetric characteristics, only its axisymmetric cross section is needed to build the geometric model [4], and the three-dimensional model is simplified to a two-dimensional plane model in order to reduce the calculation. The software used to build the geometric model is Solidworks [5], which has a good fit with ANSYS/LS-DYNA.
2.2. Delineating the Finite Element Mesh
Import the geometric model from the previous step into the meshing software. The finite element meshing is performed for the geometric characteristics of the simulation object. For the above device, its structure and thermal boundary conditions have axisymmetric characteristics, so the heat conduction process in the device can be considered as an axisymmetric problem [6], and the model cross-section is meshed using planar quadrilateral cells [7], as shown in Figure 2. This can greatly reduce the simulation calculation volume and speed up the simulation analysis.
Figure 1. Example of a figure caption (figure caption).
(a) (b)
Figure 2. A device finite element model meshing. (a) Overall meshing of the finite element model; (b) local meshing of the finite element model.
2.3. Simulation Model of Heat Transfer Process inside a Device
The delineated finite element mesh is imported into LS-PrePost, the pre-processing software of ANSYS/LS-DYNA, and keywords are added to establish a complete heat transfer simulation model. The model of a device after importing LS-PrePost is shown in Figure 3.
In order to reduce the computational effort and speed up the computer simulation, the cell type used for this heat transfer analysis is a planar quadrilateral cell.
Firstly, material models were added to the air domain (density of 1.29 kg/m3) in LS-PrePost as listed in Table 1; secondly, material models were added to the radome; then, material models were added to the heat shield; next, material models [8] were added to the thermal shield, Yin steel ring, and antenna (Table 2).
The initial temperature of the model is 313 K (39.85˚C). Add thermal boundary conditions to the model [9], as shown in Figure 4.
Table 1. Air thermal physical model parameters.
Table 2. Thermophysical parameters of model materials.
Figure 3. Model of a device after importing LS-PrePost.
(a)(b)
Figure 4. Model thermal boundary condition setting. (a) Model thermal boundary conditions; (b) thermal boundary conditions temperature change curve.
Set other parameters required for thermal simulation analysis in LS-PrePost, such as thermal conductivity between components, simulation step size, solver parameters, etc.
3. Simulation Experiment Results
By conducting thermal analysis simulation experiments on the model, the internal temperature distribution states of the model at different moments are obtained [10]. To ensure the uniformity and continuity of the overall simulation state, the temperature distribution clouds at five time points of 50 s, 100 s, 150 s, 200 s, and 250 s are intercepted by using the equal difference time division. As shown in Figure 5.
After analyzing the temperature distribution at each time point of the model, the temperature of the inner surface of the top of the heat shield [11] and the maximum temperature of the antenna were calculated to obtain the following Figure 6 and Figure 7.
Figure 6. Temperature variation of the inner surface of the top of the heat shield.
Figure 7. Variation of maximum antenna temperature.
4. Experimental Conclusion
1) After the thermal analysis of the warhead model, we found that the temperature difference between the exterior and interior of the heat shield is very obvious, before 100 s, the heat shield and the internal temperature remain basically unchanged, and after 100 s, the temperature gradually increases and the rate of temperature rise remains basically unchanged, with a high degree of linearity of temperature change.
2) The maximum temperature of the antenna is maintained at about 313 K before 50 s, and it starts to warm up after 50 s with a high degree of linearity of the temperature rise curve.
3) During the experiment, the temperature transfer between the model shell to the internal heat shield showed a clear stratification, the closer the location to the surface of the projectile the faster the rate of heat transfer, and conversely, the closer the location to the heat shield the slower the heat transfer rate of the part.