The design of unmanned aircraft vehicles during operation is exposed to significant temperature effects, which causes changes in the dielectric characteristics of radio technical elements and can lead to failures of on-board radio electronic equipment, as a result of which the task of determining the amount of temperature change in the compartments where the radio electronic equipment is located is performed. At the same time, some external elements of the aircraft vehicle can have complex structures with shapes close to spherical and different thicknesses in each cross-sectional structure.
 There is a need to develop a method for calculating the non-stationary temperature field with the purpose of finding out the real conditions in which all the elements of the aircraft structure work, which would help to determine how the temperature of the structure changes during any flight. At the same time, it is necessary to take into account the three dimensions of temperature gradients on spherical parts of structures of various shapes, namely nose part, rounded wing and others. In them, as a rule, each cross-section of the wall has its own thickness and, according to each cross-section, the temperature distribution along the thickness of the wall will correspond, that is, it is necessary to ensure the redistribution of temperature along the structural element, the calculation of which month. In addition, it is necessary to note that the wall of the aircraft can be multi-layered and the properties of its material can change along with the change in temperature. Considering the fact that it is impossible to take into account all the listed features by analytical methods, when choosing a calculation method, it will be advisable to stop at numerical methods.
 One of the most well-known numerical methods is the finite difference method with an implicit calculation scheme, which has an indisputable advantage for modeling - absolute stability with sufficient changes in time and coordinate steps. We will use the level of thermal conductivity in spherical coordinates in order to increase the accuracy of the calculation, for which solution, after finding out the boundary conditions, we will perform a one-dimensional calculation of the temperature by thickness in each mesh section of the structure, and also perform a calculation of the two-dimensional "heat flow" along the solid between the points, that are located on the same thickness, but have a temperature difference. Thus, when calculating the temperature field, the sphericity and multi-lay redness of the structural elements is taken into account, and the property of the material of the aircraft structure changes depending on the temperature.