In the article discussed the mathematical formulation and numerical algorithm for solving the problem of calculating the temperature field in the process vulcanizing of the product, whose the thermal characteristics are depended on the temperature. As a mathematical model considered the system of differential equations of heat conduction, taking into account the change in the coefficients of thermal conductivity and heat density in multilayer product of the temperature. The system of equations is solved for a given initial distribution of temperature and for a given (time-dependent) temperatures on the border of the product to the press-mold and to the diaphragm. On the border of the contacts of adjacent layers are given the condition of continuity of temperature and heat flux. Change of the thermal conductivity from the time is approximated by linear functions. The activation energy of the vulcanization process is determined on the basis of experimental data obtained in the control test samples using a reometer. Considering the function representing the corresponding integrals of the thermal conductivity, the original system of differential equations is transformed to an equivalent system of differential equations convenient for constructing numerical algorithms for solving the problem. The resulting system of partial differential equations derived using the method of finite-difference approximation is replaced by a system of algebraic equations. Solution of the system of algebraic equations is carried out under the scheme explicit difference approximation. In the article calculated the temperature field for the tire at given initial and boundary conditions. Stability and accuracy of the numerical algorithm for solving the problem is demonstrated by the calculations performed with different sampling step along the time and space coordinates. Assessment of the degree of completion of the process is carried out by calculated equivalent time for temperature values, accepted as equivalent. The developed algorithm is an important part of the algorithm for solving the problem of determining the optimal mode of curing, providing the required quality products at the lowest cost, which is very important in conditions of continuous growth in energy prices.

vulcanization, temperature field, initial and boundary conditions, the pneumatic tire, stability of the solution

1. Lukomskaya A.I., Minaev N.T., Kepersha L.M., Milkova E. M. Otsenka stepeni vulkanizatsii rezin v izdeliyakh [Assessment of the degree of vulcanization of rubber in products]. Moscow, TsNIITEneftekhim, 1972. 43 p. (In Russ.).

2. Denisov A. P., Gromov Yu. Yu., Pokornyi Yu. V. Investigation of the vulcanization process at the local tire repair VGU, MGU. Voronezh. Informatsionnye tekhnologii i sistemy. [Voronezh State University, Moscow State University. Voronezh. Information technologies and systems], 1999, pp. 97. (In Russ.). 3

3. Karmanova O. V., Tikhomirov S. G., Pyatakov Yu.V., Kasperovich A. V. et al. Modelling of the kinetics of non-isometric massive rubber vulcanization. Proceedings BSTU. [Trudy BGTU], 2014, no. 4, pp. 100–104. (In Russ.).

4. Ishchenko V. A., Shaptala M. V. Features of the operating conditions of vulcanization of pneumatic tires with the three-dimensional structure. Sistemnye tekhnologii: regional'nyi mezhvuz. sb. nauch. trudov. [System Technology: Regional Interuniversity collection of scientific papers], 2008, issue 2 (55), pp. 147 – 158. (In Russ.).

5. Vlasko A. V., Sakharov M. E., Poritskaya Z. Effect of non-isothermal curing and mechanical properties of rubber and rubber-cord samples. Kauchuk i rezina. [Rubber and rubber], 1991, no. 6, pp. 6-8. (In Russ.).

6. Samarskii A.A., Tikhonov A.N. Uravneniya matematicheskoi fiziki [Equations of mathematical physics]. Moscow, Izd-vo MGU, 1999. 798 p. (In Russ.).

7.