A tire vulcanization temperature field calculation method, device, equipment and medium

Abstract

The application discloses a tire vulcanization temperature field calculation method, device, equipment and medium, relates to the tire vulcanization temperature field calculation technical field, and includes: the target tire model is divided to obtain the initial grid; the initial grid information, the grid division period and the calculation time step are determined, and are imported to the heat conduction solving program to calculate to obtain each finite element node temperature of the first time step; the error of each finite element node temperature is calculated to obtain each unit error, and the target unit error is determined; whether the target unit error meets the preset condition is judged, if it meets, the above steps are repeatedly executed to obtain the finite element node temperature of all time steps; whether the difference between the finite element node temperatures corresponding to adjacent time steps of the current period is within the preset temperature difference range is judged, if it is within the range, the calculation is ended. Can increase the calculation efficiency, reduce the tire vulcanization temperature field calculation error, satisfy the business demand diversity.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of tire vulcanization temperature field calculation, and in particular relates to a tire vulcanization temperature field calculation method, device, equipment and medium. BACKGROUND

[0002] Vulcanization is the last process in the tire production process, which is a process of changing the molecular structure of the rubber by heating. The tire is a thick product of poor thermal conductor, and the internal temperature field cannot be uniform during the heating vulcanization process, which can easily cause uneven vulcanization of each part of the tire, resulting in over-vulcanization and under-vulcanization, reducing the mechanical properties of the tire, and having a great influence on the quality of the tire. If the distribution of the internal temperature field of the tire during the entire vulcanization process can be obtained, the production process can be accurately adjusted to ensure product quality. In recent years, compared with other calculation methods, the finite element method has been proved to have higher precision and more extensive use in the field of tire vulcanization temperature field calculation with complex shape and multi-layer structure. And as a key step of finite element calculation preprocessing, the division of finite element mesh directly affects the precision of the calculation result and the size of the calculation scale. Generally speaking, the denser the mesh, the higher the calculation precision, but the time cost will also increase. Therefore, when dividing the mesh, the pros and cons generated should be considered comprehensively. At present, the main method of dividing the mesh in this field is: first, the density of the mesh of each part of the tire is determined artificially, and then the software is used to automatically divide. The defects of this method are: first, when determining the density of the mesh of each part of the tire, there is strong subjectivity, and there is no accurate mesh division method, but only the mesh of each part of the tire is encrypted or sparse processed according to experience; second, only one set of mesh is used in the entire vulcanization process, which cannot adapt to the changing temperature field, and there are problems of over-dense or over-sparse in some parts, thereby causing the decrease of calculation efficiency and the deviation of calculation result.

[0003] From the above, in the process of tire vulcanization temperature field calculation, how to divide the mesh, increase the efficiency of tire vulcanization temperature field calculation, reduce the error of tire vulcanization temperature field calculation, and meet the diversity of business needs is a problem to be solved in the field. SUMMARY

[0004] Therefore, the purpose of the present application is to provide a tire vulcanization temperature field calculation method, device, equipment and medium, which can effectively increase the efficiency of tire vulcanization temperature field calculation, reduce the error of tire vulcanization temperature field calculation, and meet the diversity of business needs. The specific scheme is as follows:

[0005] In a first aspect, the present application discloses a tire vulcanization temperature field calculation method, comprising:

[0006] acquire a target tire model, perform finite element meshing on the target tire model to obtain an initial mesh;

[0007] determine initial mesh information, a preset meshing period, and a calculation time step, and import the initial mesh information, the meshing period, and the calculation time step into a preset heat conduction solving program to perform node temperature calculation to obtain each finite element node temperature of a first time step of a current period;

[0008] perform error calculation on each of the finite element node temperatures to obtain each element error, and determine a target element error from all of the element errors;

[0009] determine whether the target element error meets a preset condition, and if the target element error meets the preset condition, repeat the above steps to obtain each finite element node temperature of all time steps of the current period;

[0010] determine whether a difference between the finite element node temperatures corresponding to adjacent time steps of the current period is within a preset temperature difference range, and if the difference between the finite element node temperatures corresponding to adjacent time steps of the current period is within the preset temperature difference range, end the tire vulcanization temperature field calculation.

[0011] Optionally, the acquiring of the target tire model and the performing of the finite element meshing on the target tire model to obtain the initial mesh comprises:

[0012] filtering, according to a tire type, a target tire model from all tire models;

[0013] performing finite element meshing on the target tire model by using a preset finite element simulation software to obtain the initial mesh.

[0014] Optionally, the determining of the initial mesh information, the preset meshing period, and the calculation time step comprises:

[0015] determining the preset meshing period and the calculation time step, and extracting the initial mesh to obtain initial mesh information including index and coordinate information of initial mesh nodes, global number information of each element node, and global number information of inner and outer boundary points;

[0016] saving the initial mesh information including the index and coordinate information of the initial mesh nodes, the global number information of each element node, and the global number information of the inner and outer boundary points to a local device in the form of an array.

[0017] Optionally, the initial mesh information, the mesh division period and the calculation time step are introduced into a preset heat conduction solving program for node temperature calculation to obtain each finite element node temperature of a first time step of a current period.

[0018] The heat conduction solving program is determined based on a weighted residual method.

[0019] The initial mesh information, the preset mesh division period and the time step are introduced into the heat conduction solving program for node temperature calculation to obtain each finite element node temperature of a first time step of a current period.

[0020] Optionally, the error of each finite element node temperature is calculated to obtain each element error, and a target element error is determined from all the element errors, including:

[0021] The error of each finite element node temperature is calculated by using a preset error calculation method to obtain each element error.

[0022] The largest element error is selected from all the element errors as the target element error.

[0023] Optionally, after determining whether the target element error meets a preset condition, the method further includes:

[0024] If the target element error does not meet the preset condition, the initial mesh element corresponding to the target element error is marked and refined by using a preset marking strategy program to form a new mesh for re-calculation until the target element error meets the preset condition.

[0025] Optionally, after the tire vulcanization temperature field calculation is completed, the method further includes:

[0026] All finite element node temperatures are obtained, and a tire size is determined based on the target tire model.

[0027] A cross-sectional temperature field image of the target tire model is drawn based on all the finite element node temperatures and the tire size.

[0028] In a second aspect, the application discloses a tire vulcanization temperature field calculation device, including:

[0029] A mesh division module is configured to obtain a target tire model, perform finite element mesh division on the target tire model, and obtain an initial mesh.

[0030] a grid temperature calculation module, configured to determine initial grid information, a preset grid division period and a calculation time step, and to input the initial grid information, the grid division period and the calculation time step into a preset heat conduction solving program to calculate node temperatures to obtain each finite element node temperature of a first time step of a current period;

[0031] an error calculation module, configured to calculate errors of each of the finite element node temperatures to obtain each element error, and to determine a target element error from all the element errors;

[0032] a first judging module, configured to judge whether the target element error meets a preset condition, and if the target element error meets the preset condition, to repeat the above steps to obtain each finite element node temperature of all time steps of the current period;

[0033] a second judging module, configured to judge whether a difference between the finite element node temperatures corresponding to adjacent time steps of the current period is within a preset temperature difference range, and if the difference between the finite element node temperatures corresponding to adjacent time steps of the current period is within the preset temperature difference range, to end the tire vulcanization temperature field calculation.

[0034] In a third aspect, the present application discloses an electronic device, comprising:

[0035] a memory, configured to save a computer program;

[0036] a processor, configured to execute the computer program to implement the tire vulcanization temperature field calculation method.

[0037] In a fourth aspect, the present application discloses a computer storage medium, configured to save a computer program; wherein the computer program is executed by a processor to implement the steps of the tire vulcanization temperature field calculation method disclosed above.

[0038] It can be seen that the application provides a tire vulcanization temperature field calculation method, which comprises the following steps: obtaining a target tire model, performing finite element grid division on the target tire model to obtain an initial grid; determining initial grid information, a preset grid division period and a calculation time step, and importing the initial grid information, the grid division period and the calculation time step into a preset heat conduction solving program for node temperature calculation to obtain the finite element node temperature of the first time step of the current period; performing error calculation on each finite element node temperature to obtain each unit error, and determining a target unit error from all the unit errors; judging whether the target unit error meets a preset condition, if the target unit error meets the preset condition, repeating the above steps to obtain the finite element node temperature of all time steps of the current period; judging whether the difference between the finite element node temperatures corresponding to adjacent time steps of the current period is within a preset temperature difference range, if the difference between the finite element node temperatures corresponding to adjacent time steps of the current period is within the preset temperature difference range, ending the tire vulcanization temperature field calculation. In the process of tire vulcanization temperature field calculation, the grid can be automatically adjusted according to the solution change and the need according to the set calculation period, the grid spacing in the large gradient area of the temperature field can be effectively reduced, thus the error can be reduced and the precision can be improved; the grid spacing in the area with small gradient is large, thus the boundary value information can be transmitted faster, the convergence can be accelerated, and the calculation efficiency can be improved. BRIEF DESCRIPTION OF DRAWINGS

[0039] In order to more clearly illustrate the technical solutions in the embodiments of the present application or the prior art, the following will briefly introduce the drawings needed to be used in the embodiments or the prior art description. Obviously, the drawings in the following description are only the embodiments of the present application, and for those skilled in the art, other drawings can also be obtained without creative labor based on the provided drawings.

[0040] Figure 1 A flow chart of a tire vulcanization temperature field calculation method disclosed by the present application;

[0041] Figure 2 A flow chart of a tire vulcanization temperature field calculation method disclosed by the present application;

[0042] Figure 3 A detailed method schematic diagram of a tire vulcanization temperature field calculation method disclosed by the present application;

[0043] Figure 4 A comparison diagram of a tire vulcanization temperature field calculation method disclosed by the present application;

[0044] Figure 5 A tire vulcanization temperature field cross-section temperature schematic diagram disclosed by the present application;

[0045] Figure 6 A tire vulcanization temperature field calculation method disclosed by the application is a specific flowchart;

[0046] Figure 7 A tire vulcanization temperature field calculation device structure schematic diagram disclosed by the application;

[0047] Figure 8 An electronic device structure diagram provided by the application. DETAILED DESCRIPTION

[0048] The technical solutions in the embodiments of the application will be clearly and completely described in connection with the drawings in the embodiments of the application. Obviously, the described embodiments are only part of the embodiments of the application, rather than all the embodiments of the application. Based on the embodiments in the application, all other embodiments obtained by those skilled in the art without creative labor fall within the protection scope of the application.

[0049] At present, vulcanization is the last process in the tire production process, which is a process of changing the molecular structure of the rubber by heating. The tire is a thick product of poor thermal conductor, and the internal temperature field cannot be uniform during the heating vulcanization process, which is easy to cause uneven vulcanization of each part of the tire, resulting in over-vulcanization and under-vulcanization, reducing the mechanical properties of the tire, and having a great influence on the quality of the tire. If the distribution of the internal temperature field of the tire during the entire vulcanization process can be obtained, the production process can be accurately adjusted, so as to ensure the product quality. In recent years, compared with other calculation methods, the finite element method for solving the heat conduction equation has been proved to have higher precision and more extensive use in the field of tire vulcanization temperature field calculation with complex shape and multi-layer structure. And as a key step of finite element calculation preprocessing, the division of finite element mesh directly affects the precision of calculation results and the size of calculation scale. Generally speaking, the denser the mesh, the higher the calculation precision, but the time cost will also increase. Therefore, when dividing the mesh, the pros and cons generated should be considered comprehensively. At present, the main method of dividing the mesh in this field is: first, the density of the mesh of each part of the tire is determined artificially, and then the software is used to automatically divide. The defects of this method are: first, when determining the density of the mesh of each part of the tire, there is strong subjectivity, there is no accurate mesh division method, only the mesh of each part of the tire is encrypted or sparse processed according to experience; second, only one set of mesh is used in the entire vulcanization process, which cannot adapt to the changing temperature field, and there are problems of over-dense or over-sparse in some parts, thereby leading to the decrease of calculation efficiency and the deviation of calculation results. As can be seen from the above, in the process of tire vulcanization temperature field calculation, how to divide the mesh, increase the efficiency of tire vulcanization temperature field calculation, reduce the error of tire vulcanization temperature field calculation, and meet the diversity of business needs is a problem to be solved in the field.

[0050] Referring to Figure 1 As shown in the drawings, the embodiment of the present application discloses a tire vulcanization temperature field calculation method, which can specifically include:

[0051] Step S11: Obtain a target tire model, and perform finite element mesh division on the target tire model to obtain an initial mesh.

[0052] In this embodiment, the target tire model is selected from all tire models according to the tire type; and the finite element mesh division is performed on the target tire model by using a preset finite element simulation software to obtain the initial mesh. In this embodiment, the tire model is a three-dimensional model; and the unit mesh is a three-dimensional tetrahedral unit mesh.

[0053] In this embodiment, the obtained three-dimensional tire model is imported into the ABAQUS finite element simulation software, the approximate global size of the global seed is set to 50 by using the ready-to-use automatic mesh division module of the software, the unit shape is set to tetrahedron, a set of meshes with 1920 nodes and 6200 tetrahedral units are divided, and the meshes are used as the initial meshes for calculation.

[0054] Step S12: Determine the initial mesh information, a preset mesh division period and a calculation time step, and import the initial mesh information, the mesh division period and the calculation time step into a preset heat conduction solving program to perform node temperature calculation, so as to obtain the finite element node temperature of the first time step in the current period.

[0055] In this embodiment, after obtaining each unit mesh, the preset mesh division period and the calculation time step are determined, and the initial mesh is extracted to obtain the initial mesh information including the index and coordinate information of the initial mesh nodes, the global number information of each unit node and the global number information of the inner and outer boundary points; and the initial mesh information including the index and coordinate information of the initial mesh nodes, the global number information of each unit node and the global number information of the inner and outer boundary points is saved to the local in the form of an array.

[0056] In this embodiment, the initial mesh information is specifically stored in the form of an array: the array Vertices contains the index and coordinate information of the mesh nodes, the array Elements contains the global number information of each unit, and the array Boundarynodes contains the global number information of the inner and outer boundary points. The theoretical basis of the program is that the control equation of the heat transfer problem is converted into the finite element analysis formula of the transient heat transfer problem based on the weighted residual method.

[0057] In the embodiment, the specific initial grid information includes: index and coordinates of the unit grid nodes, global number of each unit grid, and global number of the inner and outer boundary points. The calculation period refers to the period of grid re-partitioning. The theoretical basis of the program is that the control equation of the heat transfer problem is converted into a finite element analysis equation of the transient heat transfer problem based on the weighted residual method. The finite element analysis equation is:

[0058]

[0059] wherein N is the number of finite element nodes in each unit, N is 4 because it is a three-dimensional tetrahedral linear unit, λ is the thermal conductivity of the tire material, and the value is 0.233 W / m / ℃, ρ is the material density, the unit is kg / m3, c T is the specific heat capacity of the material, the unit is J / kg / ℃, ρc T is 1750000, φ i and φ j are shape functions of the unit.

[0060] Next, the interpolation of and is performed on the time t, and if the entire time t is divided into n time steps, n linear equations can be obtained: m = 1, 2,..., n;

[0061] wherein,

[0062]

[0063] In the embodiment, the main steps of the heat conduction solving program are as follows: firstly, the unit base functions of the three-dimensional tetrahedral linear unit are created; then the heat transfer matrix of each unit is calculated according to each unit base function; then the total heat transfer matrix A is synthesized; then the initial condition and the boundary condition are set, and the temperature field at the first time step is calculated, that is,

[0064] Specifically, the initial temperature is set to room temperature 25℃, the outer boundary condition: the model is in contact with the hot plate, the heat resistance is ignored, and is simplified as the first type of boundary condition, and the temperature is set to 141℃; the inner boundary condition: the capsule and the overheated water exchange heat through convection, since the thermal conductivity of rubber is very small, the main thermal resistance of heat transfer is the thermal resistance of rubber, and it is also simplified as the first type of boundary condition, and the temperature is set to 165℃. A calculation time step of 0.04s is set here, since the temperature field is constantly changing in the tire vulcanization process, and considering the time cost, a separate set of calculation grids cannot be divided for each time step, so the period of grid re-partitioning is set to 10s;

[0065] Step S13: error calculation is performed on each of the finite element node temperatures to obtain a unit error, and a target unit error is determined from all of the unit errors.

[0066] Step S14: it is determined whether the target unit error meets a preset condition, and if the target unit error meets the preset condition, the above steps are repeatedly performed to obtain each finite element node temperature of all time steps in a current period.

[0067] In this embodiment, after it is determined whether the target unit error meets the preset condition, the method further comprises:

[0068] If the target unit error does not meet the preset condition, the initial mesh unit corresponding to the target unit error is marked and refined by using a preset marking strategy program, a new mesh is formed to perform calculation again until the target unit error meets the preset condition. In this embodiment, the marking strategy mainly includes the following three kinds:

[0069] M r = {e | e ∈ M, η(e) > θ · max e∈M η(e)} ;

[0070] M r is a minimum subset to meet ; M r = {e | e ∈ M, η(e) > θ · tol} ;

[0071] Wherein, M r is a refined unit set, θ is a refinement ratio, and tol is a self-set limit condition. In this embodiment, the refinement method is a binary refinement method of a specified refinement edge.

[0072] Step S15: it is determined whether a difference between the finite element node temperatures corresponding to adjacent time steps in the current period is within a preset temperature difference range, and if the difference between the finite element node temperatures corresponding to adjacent time steps in the current period is within the preset temperature difference range, the calculation of the tire vulcanization temperature field is ended.

[0073] In the embodiment, a target tire model is acquired, finite element meshing is performed on the target tire model to obtain an initial mesh, initial mesh information, a preset meshing period and a calculation time step are determined, and the initial mesh information, the meshing period and the calculation time step are imported into a preset heat conduction solving program for node temperature calculation to obtain the finite element node temperature of the first time step of the current period. Errors of the finite element node temperatures are calculated to obtain unit errors, and a target unit error is determined from all the unit errors. It is judged whether the target unit error meets a preset condition. If the target unit error meets the preset condition, the above steps are repeatedly performed to obtain the finite element node temperature of all time steps of the current period. It is judged whether the difference between the finite element node temperatures corresponding to adjacent time steps of the current period is within a preset temperature difference range. If the difference between the finite element node temperatures corresponding to adjacent time steps of the current period is within the preset temperature difference range, the tire curing temperature field calculation is ended. In the process of tire curing temperature field calculation, the grid can be automatically adjusted according to the solution change and the need according to the preset meshing period, the grid spacing in the large gradient area of the temperature field can be effectively reduced, and thus the error can be reduced and the precision can be improved. The grid spacing in the area with small gradient is large, and thus the boundary value information can be transmitted faster, the convergence can be accelerated, and the calculation efficiency can be improved.

[0074] Referring to Figure 2 The embodiment of the application discloses a tire curing temperature field calculation method, which can specifically include the following steps:

[0075] Step S21: acquiring a target tire model, performing finite element meshing on the target tire model to obtain an initial mesh.

[0076] Step S22: determining initial mesh information, a preset meshing period and a calculation time step, and determining a heat conduction solving program based on a weighted residual method, importing the initial mesh information, the meshing period and the time step into the heat conduction solving program for node temperature calculation to obtain the finite element node temperature of the first time step of the current period.

[0077] Step S23: using a preset error calculation method to calculate errors of all the finite element node temperatures to obtain unit errors, and screening the unit error with the largest value from all the unit errors as a target unit error.

[0078] In the embodiment, the unit error is calculated in the unit of each finite element. Specifically, in the embodiment, there is no heat source in the tire, so the right end term f h is 0, and the measurement method of the unit error η e is:

[0079]

[0080] Where h e Let h be the diameter of element e, F(e) be the set of faces of e, and h be the diameter of element e. f Let n be the diameter of surface f. f Let f be the unit normal vector of surface f.

[0081] Will Substituting into the above formula for calculating unit error, the error η of each unit can be calculated. e And find the maximum value among them, denoted as Max(η). e ).

[0082] Step S24: Determine whether the target element error meets the preset conditions. If the target element error meets the preset conditions, repeat the above steps to obtain the temperature of each finite element node in all time steps of the current cycle.

[0083] In this embodiment, the maximum unit error Max(η) is determined. e If the constraint condition is not met, the elements are marked and refined according to the marking strategy to form a new mesh for recalculation until the maximum element error meets the constraint condition. If the constraint condition is met, the mesh information is used as the unified mesh for finite element calculations within this mesh generation cycle and is used in the remaining time steps of this cycle. Specifically, the constraint condition here is Max(η) e )<=1E-12; the marking strategy is M r ={e∣e∈M,η(e)>θ·max e∈M η(e)}, θ is set to 0.01, i.e., M r This is the set of cells that need to be refined. The bisection refinement method: When refining a cell, one of its edges is designated as the refinement edge. The midpoint of the refinement edge is connected to the two opposite vertices, thus dividing the cell in two, such as... Figure 3 As shown. The final calculated mesh information for the first mesh generation cycle is: 5721 vertices and 22363 elements. A comparison with the initial mesh is shown below. Figure 4 As shown.

[0084] Step S25: Determine whether the temperature difference between the finite element nodes corresponding to adjacent time steps in the current cycle is within the preset temperature difference range. If the temperature difference between the finite element nodes corresponding to adjacent time steps in the current cycle is within the preset temperature difference range, then end the tire vulcanization temperature field calculation.

[0085] The above steps are repeated in the next cycle to readjust the initial mesh and perform calculation, and the process is repeated to complete the temperature field calculation of the tire throughout the vulcanization process. Then, it is determined whether the difference between the finite element node temperatures corresponding to adjacent time steps in the current cycle is within the preset temperature difference range. If the difference between the finite element node temperatures corresponding to adjacent time steps in the current cycle is within the preset temperature difference range, the tire vulcanization temperature field calculation is ended. After calculation, when the time reaches 2000s, i.e., 200 mesh redrawing cycles, the temperature field basically reaches a steady state. The temperature field of the tire cross section when reaching the steady state is shown in Figure 5 .

[0086] In this embodiment, after the tire vulcanization temperature field calculation is ended, all finite element node temperatures are obtained, and the tire size is determined based on the target tire model. The cross-sectional temperature field image of the target tire model is drawn based on all the finite element node temperatures and the tire size, as shown in Figure 5 . The abscissa and ordinate are the tire size of the target tire model, and the contour graph in the figure represents the temperature cross section of the tire.

[0087] In this embodiment, the specific process is shown in Figure 6 . First, a target tire model is obtained, and finite element mesh partitioning is performed on the target tire model to obtain an initial mesh. Then, the initial mesh information, mesh partitioning cycle, and calculation time step are determined, and the initial mesh information, mesh partitioning cycle, and calculation time step are imported into a preset heat conduction solving program for node temperature calculation to obtain the finite element node temperature of the first time step in the current cycle. Then, the error of each finite element node temperature is calculated for each finite element to obtain the unit error, and the maximum value of the unit error is determined as the target unit error. It is determined whether the target unit error meets the preset condition. If the target unit error meets the preset condition, the mesh information is used as the unified mesh for finite element calculation in the current cycle, and the remaining time steps in the current cycle are calculated. Then, the above steps are repeated to complete the temperature field calculation of the tire throughout the vulcanization process. If the target unit error does not meet the preset condition, the unit mesh information is marked and refined according to the marking strategy to obtain a new mesh, and the step of determining the mesh information, mesh partitioning cycle, and calculation time step is jumped to.

[0088] In this embodiment, during the vulcanization temperature field calculation process, the mesh can be automatically adjusted according to the solution change and needs, which can effectively reduce the grid spacing in the large gradient area of the temperature field, thereby reducing the error and improving the precision. The grid spacing in the area with small gradient is larger, so the boundary value information can be transmitted faster, the convergence is accelerated, and the calculation efficiency is improved.

[0089] In the embodiment, a target tire model is acquired, finite element meshing is performed on the target tire model to obtain an initial mesh, initial mesh information, a preset meshing period and a calculation time step are determined, and the initial mesh information, the meshing period and the calculation time step are input into a preset heat conduction solving program for node temperature calculation to obtain the finite element node temperature of the first time step of the current period, errors of the finite element node temperatures are calculated to obtain unit errors, and a target unit error is determined from all the unit errors, whether the target unit error meets a preset condition is judged, if the target unit error meets the preset condition, the above steps are repeatedly executed to obtain the finite element node temperature of all time steps of the current period, whether the difference between the finite element node temperatures corresponding to adjacent time steps of the current period is within a preset temperature difference range is judged, and if the difference between the finite element node temperatures corresponding to adjacent time steps of the current period is within the preset temperature difference range, the tire vulcanization temperature field calculation is ended. In the process of tire vulcanization temperature field calculation, the mesh can be automatically adjusted according to the solution change and the need according to the preset calculation period, the grid spacing in the large gradient area of the temperature field can be effectively reduced, thus the error can be reduced and the precision can be improved, the grid spacing in the area with small gradient is large, thus the boundary value information can be transmitted faster, the convergence can be accelerated, and the calculation efficiency can be improved.

[0090] Referring to Figure 7 The embodiment of the application discloses a tire vulcanization temperature field calculation device, which can specifically include:

[0091] The meshing module 11 is configured to acquire a target tire model, perform finite element meshing on the target tire model to obtain an initial mesh, determine initial mesh information, a preset meshing period and a calculation time step, and input the initial mesh information, the meshing period and the calculation time step into a preset heat conduction solving program for node temperature calculation to obtain the finite element node temperature of the first time step of the current period.

[0092] The mesh temperature calculation module 12 is configured to determine initial mesh information, a preset meshing period and a calculation time step, and input the initial mesh information, the meshing period and the calculation time step into a preset heat conduction solving program for node temperature calculation to obtain the finite element node temperature of the first time step of the current period.

[0093] The error calculation module 13 is configured to calculate errors of the finite element node temperatures to obtain unit errors, and determine a target unit error from all the unit errors.

[0094] The first judgment module 14 is configured to judge whether the target unit error meets a preset condition, and if the target unit error meets the preset condition, repeatedly execute the above steps to obtain the finite element node temperature of all time steps of the current period.

[0095] The second judging module 15 is configured to judge whether the difference between the finite element node temperatures corresponding to adjacent time steps in the current period is within the preset temperature difference range. If the difference between the finite element node temperatures corresponding to adjacent time steps in the current period is within the preset temperature difference range, the calculation of the tire vulcanization temperature field is ended.

[0096] In the embodiment, a target tire model is obtained, finite element meshing is performed on the target tire model to obtain an initial mesh, initial mesh information, a preset meshing period and a calculation time step are determined, and the initial mesh information, the meshing period and the calculation time step are input into a preset heat conduction solving program to calculate node temperatures to obtain the finite element node temperatures of the first time step in the current period. Errors of the finite element node temperatures are calculated to obtain unit errors, and a target unit error is determined from all the unit errors. It is judged whether the target unit error meets a preset condition. If the target unit error meets the preset condition, the above steps are repeatedly performed to obtain the finite element node temperatures of all time steps in the current period. It is judged whether the difference between the finite element node temperatures corresponding to adjacent time steps in the current period is within the preset temperature difference range. If the difference between the finite element node temperatures corresponding to adjacent time steps in the current period is within the preset temperature difference range, the calculation of the tire vulcanization temperature field is ended. In the process of calculating the tire vulcanization temperature field, the mesh can be automatically adjusted according to the set calculation period according to the solution change and the need, the grid spacing in the large gradient area of the temperature field can be effectively reduced, and thus the error can be reduced and the precision can be improved. The grid spacing in the area with small gradient is large, and thus the boundary value information can be transmitted faster, the convergence can be accelerated, and the calculation efficiency can be improved.

[0097] In some specific embodiments, the meshing module 11 can specifically include:

[0098] The target tire model screening module is configured to screen a target tire model from all tire models according to the tire type.

[0099] The meshing module is configured to perform finite element meshing on the target tire model by using a preset finite element simulation software to obtain an initial mesh.

[0100] In some specific embodiments, the mesh temperature calculation module 12 can specifically include:

[0101] The mesh extraction module is configured to determine a preset meshing period and a calculation time step, and extract the initial mesh to obtain initial mesh information including index, coordinate information of initial mesh nodes, global number information of each unit node and global number information of inner and outer boundary points.

[0102] An information storage module is configured to store initial grid information including indexes of initial grid nodes, coordinate information, global number information of each element node, and global number information of inner and outer boundary points in an array to the local device.

[0103] In some embodiments, the grid temperature calculation module 12 can specifically include:

[0104] A program determination module is configured to determine a heat conduction solving program based on a weighted residual method.

[0105] A temperature calculation module is configured to input the initial grid information, a preset grid division period, and a time step into the heat conduction solving program to calculate node temperatures, so as to obtain finite element node temperatures of a first time step in a current period.

[0106] In some embodiments, the error calculation module 13 can specifically include:

[0107] An error calculation module is configured to calculate errors of all the finite element node temperatures by using a preset error calculation method, so as to obtain element errors.

[0108] A target element error determination module is configured to select a maximum element error from all the element errors as a target element error.

[0109] In some embodiments, the first judgment module 14 can specifically include:

[0110] A new element grid determination module is configured to, if the target element error does not meet a preset condition, mark and refine the initial grid element corresponding to the target element error by using a preset marking strategy program, form a new grid, and perform calculation again until the target element error meets the preset condition.

[0111] In some embodiments, the second judgment module 15 can specifically include:

[0112] A tire size determination module is configured to obtain all the finite element node temperatures, and determine a tire size based on the target tire model.

[0113] A graph drawing module is configured to draw a cross-sectional temperature field image of the target tire model based on all the finite element node temperatures and the tire size.

[0114] Figure 8A structural schematic diagram of an electronic device is provided in the embodiments of the present application. The electronic device 20 can specifically include at least one processor 21, at least one memory 22, a power supply 23, a communication interface 24, an input / output interface 25, and a communication bus 26. The memory 22 is configured to store a computer program, and the processor 21 is configured to load and execute the computer program to implement the related steps in the tire vulcanization temperature field calculation method performed by the electronic device disclosed in any of the foregoing embodiments.

[0115] In the embodiments, the power supply 23 is configured to provide working voltage for each hardware device on the electronic device 20, the communication interface 24 is capable of creating a data transmission channel between the electronic device 20 and external devices, and the communication protocol followed by the communication interface 24 can be any communication protocol applicable to the technical solutions of the present application, which is not specifically limited herein. The input / output interface 25 is configured to obtain external input data or output data to the outside, and the specific interface type can be selected according to the specific application requirements, which is not specifically limited herein.

[0116] In addition, the memory 22 as a carrier for resource storage can be a read-only memory, a random access memory, a magnetic disk, or an optical disk, and the resources stored thereon include an operating system 221, a computer program 222, and data 223, etc., and the storage mode can be temporary storage or permanent storage.

[0117] The operating system 221 is configured to manage and control each hardware device on the electronic device 20 and the computer program 222, so as to implement the operation and processing of the processor 21 on the data 223 in the memory 22, and the operating system 221 can be Windows, Unix, Linux, etc. The computer program 222 can further include computer programs for completing other specific work in addition to the computer programs for completing the tire vulcanization temperature field calculation method performed by the electronic device 20 disclosed in any of the foregoing embodiments. The data 223 can include data transmitted by external devices and received by the tire vulcanization temperature field calculation device, and can also include data collected by the input / output interface 25, etc.

[0118] The steps of the methods or algorithms described in conjunction with the embodiments disclosed herein can be implemented directly by hardware, a software module executed by a processor, or a combination thereof. The software module can be placed in a random access memory (RAM), a memory, a read-only memory (ROM), an electrically programmable ROM, an electrically erasable programmable ROM, a register, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the technical field.

[0119] Furthermore, this application also discloses a computer-readable storage medium storing a computer program. When the computer program is loaded and executed by a processor, it implements the tire vulcanization temperature field calculation method steps disclosed in any of the foregoing embodiments.

[0120] Finally, it should be noted that in this document, relational terms such as "first" and "second" are used only to distinguish one entity or operation from another, and do not necessarily require or imply any such actual relationship or order between these entities or operations. Furthermore, the terms "comprising," "including," or any other variations thereof are intended to cover non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitations, an element defined by the phrase "comprising one..." does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0121] The above provides a detailed description of the tire vulcanization temperature field calculation method, apparatus, equipment, and storage medium provided by the present invention. Specific examples have been used to illustrate the principles and implementation methods of the present invention. The descriptions of the above embodiments are only for the purpose of helping to understand the method and core ideas of the present invention. At the same time, for those skilled in the art, there will be changes in the specific implementation methods and application scope based on the ideas of the present invention. Therefore, the content of this specification should not be construed as a limitation of the present invention.

Claims

1. A method for calculating the vulcanization temperature field of a tire, characterized in that, include: Obtain the target tire model, and perform finite element mesh generation on the target tire model to obtain the initial mesh; The initial mesh information, the preset mesh generation period, and the calculation time step are determined, and the initial mesh information, the preset mesh generation period, and the calculation time step are imported into the preset heat conduction solution program to calculate the nodal temperature, so as to obtain the temperature of each finite element nodal in the first time step of the current period. Error calculations are performed on the temperature of each finite element node to obtain the error of each element, and the target element error is determined from all the element errors. Determine whether the target element error meets the preset conditions. If the target element error meets the preset conditions, repeat the above steps to obtain the temperature of each finite element node in all time steps of the current cycle. Determine whether the temperature difference between the finite element nodes corresponding to adjacent time steps in the current cycle is within the preset temperature difference range. If the temperature difference between the finite element nodes corresponding to adjacent time steps in the current cycle is within the preset temperature difference range, then end the tire vulcanization temperature field calculation.

2. The method for calculating the tire vulcanization temperature field according to claim 1, characterized in that, The process of obtaining the target tire model and performing finite element mesh generation on the target tire model to obtain an initial mesh includes: Select the target tire model from all tire models based on tire type; The target tire model is meshed using a pre-set finite element simulation software to obtain an initial mesh.

3. The method for calculating the tire vulcanization temperature field according to claim 1, characterized in that, The determination of the initial grid information, the preset grid division period, and the calculation time step includes: The preset mesh generation period and calculation time step are determined, and the initial mesh is extracted to obtain initial mesh information containing the index and coordinate information of the initial mesh nodes, the global number information of each unit node, and the global number information of the inner and outer boundary points; The initial grid information, including the index and coordinate information of the initial grid nodes, the global number information of each cell node, and the global number information of the inner and outer boundary points, is saved locally as an array.

4. The method for calculating the tire vulcanization temperature field according to claim 1, characterized in that, The step of importing the initial mesh information, the mesh generation period, and the calculation time step into a preset heat conduction solution program to calculate the nodal temperature, so as to obtain the finite element nodal temperature of the first time step of the current period, includes: The heat conduction solution procedure was determined based on the weighted residual method. The initial mesh information, the preset mesh generation period, and the time step are imported into the heat conduction solution program to calculate the nodal temperature, so as to obtain the temperature of each finite element nodal in the first time step of the current period.

5. The method for calculating the tire vulcanization temperature field according to claim 1, characterized in that, The step of calculating the error of the temperature at each finite element node to obtain the error of each element, and determining the target element error from all the element errors, includes: The temperature of all the finite element nodes is calculated using a preset error calculation method to obtain the error of each element. The unit error with the largest value is selected from all the unit errors as the target unit error.

6. The method for calculating the tire vulcanization temperature field according to claim 1, characterized in that, After determining whether the target unit error meets the preset conditions, the method further includes: If the target cell error does not meet the preset conditions, the initial grid cell corresponding to the target cell error is marked and refined using a preset marking strategy program to form a new grid for recalculation until the target cell error meets the preset conditions.

7. The method for calculating the tire vulcanization temperature field according to any one of claims 1 to 6, characterized in that, After completing the tire vulcanization temperature field calculation, the following is also included: Obtain all finite element node temperatures and determine the tire size based on the target tire model; Based on all the finite element node temperatures and the tire dimensions, a cross-sectional temperature field image of the target tire model is plotted.

8. A tire vulcanization temperature field calculation device, characterized in that, include: The mesh generation module is used to acquire the target tire model and perform finite element mesh generation on the target tire model to obtain the initial mesh. The mesh temperature calculation module is used to determine the initial mesh information, the preset mesh generation period, and the calculation time step, and import the initial mesh information, the mesh generation period, and the calculation time step into the preset heat conduction solution program to calculate the node temperature, so as to obtain the temperature of each finite element node in the first time step of the current period. The error calculation module is used to calculate the temperature error of each finite element node to obtain the error of each element, and to determine the target element error from all the element errors. The first judgment module is used to determine whether the target element error meets the preset conditions. If the target element error meets the preset conditions, the above steps are repeated to obtain the temperature of each finite element node in all time steps of the current cycle. The second judgment module is used to determine whether the temperature difference between the finite element nodes corresponding to adjacent time steps in the current cycle is within the preset temperature difference range. If the temperature difference between the finite element nodes corresponding to adjacent time steps in the current cycle is within the preset temperature difference range, the tire vulcanization temperature field calculation ends.

9. An electronic device, characterized in that, include: Memory, used to store computer programs; A processor for executing the computer program to implement the tire vulcanization temperature field calculation method as described in any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, Used to store computer programs; wherein, when the computer programs are executed by a processor, they implement the tire vulcanization temperature field calculation method as described in any one of claims 1 to 7.

Images