A high-load impact-resistant all-steel radial tire structure design method and system

Abstract

The present application relates to the computer aided design technical field of all-steel radial tire, and particularly relates to a high-load-bearing impact-resistant all-steel radial tire structure design method and system, comprising collecting data, constructing a parameterized geometric model, demarcating space, belt end point and tire side coordinate interval, and configuring transition reinforcing layer and auxiliary support rubber block digital model; performing discretization processing to generate finite element grid model data; loading overload mechanics and road impact parameters for nonlinear dynamics solving, outputting transient displacement response and stress tensor data; extracting belt end point stress distribution and tire side strain energy density sequence, and identifying stress concentration node coordinates; iteratively updating structure size parameters through optimization of objective function, outputting final size parameter matrix and generating target all-steel radial tire digital model file. The present application effectively improves the stress concentration problem under overload and impact working conditions, and improves the impact resistance of all-steel radial tire.

Description

Technical Field

[0001] This invention relates to the field of computer-aided design technology for all-steel radial tires, and in particular to a structural design method and system for high load-bearing and impact-resistant all-steel radial tires. Background Technology

[0002] Heavy-duty truck tires, made of steel radial steel, bear enormous loads and face complex road impacts. Researchers typically employ computer-aided engineering and finite element simulation (FEM) technology for structural design. FEM first establishes a digital geometric model and meshes the nodes. Then, under specific stress boundary conditions, it calculates the stress and strain distribution at each node, allowing designers to assess impact resistance based on stress distribution contour maps. Traditional analysis methods require repeated modifications to structural dimensions and multiple runs of simulation programs. This computational process provides a theoretical simulation basis for the initial tire development.

[0003] When all-steel radial tires are subjected to road impacts under overload conditions, they undergo extreme nonlinear physical deformation, directly causing severe stress concentration in the sidewall and belt layer. Existing finite element models struggle to accurately calculate the transient mechanical changes during complex deformation processes. The industry has yet to resolve the technical challenge of stress concentration in the sidewall and belt layer caused by overload, making it difficult for computer finite element simulations to accurately optimize impact-resistant structures. For example, when a tire rolls over the edge of a hard pothole, the belt layer endpoints and sidewall region undergo instantaneous and severe flexural deformation, causing severe distortion of the corresponding analysis mesh nodes within the finite element model due to excessive geometric deformation. This mesh distortion causes the stiffness matrix of the simulation system to exhibit a singular state. This abnormal matrix value ultimately leads to the solutionr's calculation interruption or output deviating from the actual stress distribution data, preventing computer-aided design systems from automatically completing iterative optimization of impact-resistant structures based on distorted data results. Summary of the Invention

[0004] To address the shortcomings of existing technologies, this invention provides a design method and system for a high-load-bearing and impact-resistant all-steel radial tire structure. This invention solves the technical problem that stress concentration in the tire sidewall and belt layer caused by overload makes it difficult to accurately optimize the impact-resistant structure through computer finite element simulation.

[0005] To solve the above-mentioned technical problems, the specific contents of the present invention are as follows: In a first aspect, the present invention provides a high-load-bearing and impact-resistant all-steel radial tire structure design method, comprising: Collect geometric contour coordinate data, material constitutive parameters, surcharge mechanical parameters, and road impact parameters; A parametric geometric model is constructed based on the geometric contour coordinate data and the material constitutive parameters. Within the parametric geometric model, spatial coordinate intervals, belt layer endpoint coordinate intervals, and tire sidewall coordinate intervals are defined. Within the spatial coordinate intervals, a digital model of a transition reinforcement layer with a first size parameter and a digital model of an auxiliary support block with a second size parameter are configured. The parametric geometric model configured with the digital model of the transition reinforcement layer and the digital model of the auxiliary support block are discretized to output finite element mesh model data with node coordinate sequences and material property arrays. The overload mechanical parameters and the road impact parameters are loaded into the finite element mesh model data. Nonlinear dynamics calculations are performed on the finite element mesh model data after loading the overload mechanical parameters and the road impact parameters, and transient displacement response data and stress tensor data are output. Based on the transient displacement response data and the stress tensor data, extract the first transient stress distribution data corresponding to the coordinate interval of the belt layer endpoint and the second strain energy density sequence corresponding to the coordinate interval of the tire sidewall, and extract the coordinates of the stress concentration nodes based on the first transient stress distribution data and the second strain energy density sequence. The first transient stress distribution data, the second strain energy density sequence, and the coordinates of the stress concentration nodes are input into a preset optimization objective function. The first size parameter and the second size parameter are iteratively updated using the preset optimization objective function, and the final size parameter matrix is ​​output. The target all-steel radial tire digital model file is generated based on the final size parameter matrix.

[0006] Furthermore, the high-load-bearing and impact-resistant all-steel radial tire structure design method of the present invention includes defining spatial coordinate intervals, belt layer end point coordinate intervals, and sidewall coordinate intervals within the parametric geometric model; configuring a transition reinforcement layer digital model with a first size parameter and an auxiliary support block digital model with a second size parameter within the spatial coordinate intervals; discretizing the parametric geometric model configured with the transition reinforcement layer digital model and the auxiliary support block digital model; and outputting finite element mesh model data with node coordinate sequences and material property arrays, including: Extract the boundary coordinates of the carcass ply and the belt layer from the geometric contour coordinate data. Define a first spatial coordinate interval between the outer periphery of the carcass ply boundary coordinates and the inner periphery of the belt layer boundary coordinates. Set the digital model of the transition reinforcement layer loaded with the first size parameters within the first spatial coordinate interval. The first size parameters include an initial cord angle variable and an initial width variable. Extract the inner liner boundary coordinates from the geometric contour coordinate data, define a second spatial coordinate interval between the outer perimeter of the inner liner boundary coordinates and the inner perimeter of the carcass ply boundary coordinates, set the auxiliary support block digital model loaded with the second size parameter in the second spatial coordinate interval, the second size parameter includes an initial hardness variable and an initial height variable, and the spatial coordinate interval includes the first spatial coordinate interval and the second spatial coordinate interval; Extract the belt layer endpoint coordinates and sidewall coordinates from the geometric contour coordinate data, and define the belt layer endpoint coordinate interval and the sidewall coordinate interval according to the belt layer endpoint coordinates and the sidewall coordinates respectively; The parametric geometric model loaded with the digital model of the transition reinforcement layer and the digital model of the auxiliary support block is discretized by a finite element mesh generation program, and the output is the finite element mesh model data containing the node coordinate sequence and the material property array.

[0007] Furthermore, the high-load-bearing and impact-resistant all-steel radial tire structure design method of the present invention, wherein the nonlinear dynamic solution calculation is performed on the finite element mesh model data after loading the overload mechanical parameters and the road impact parameters, and the transient displacement response data and stress tensor data are output, includes: Identify the nodes corresponding to the digital model of the transition reinforcement layer and the digital model of the auxiliary support block in the finite element mesh model data, extract the coordinates corresponding to the identified nodes as a sequence of force-bearing node coordinates, and continuously extract the mesh distortion data corresponding to the sequence of force-bearing node coordinates during the nonlinear dynamic solution calculation. Obtain the input dynamic mesh reconstruction rule parameters and mesh distortion threshold; When the mesh distortion data exceeds the mesh distortion threshold, the dynamic mesh reconstruction rule parameters are triggered, and the updated sequence of force-bearing node coordinates is recalculated and output according to the nonlinear geometric large deformation rule. The updated sequence of force-bearing node coordinates is used to assemble the global stiffness matrix, which is then output as an updated global stiffness matrix. The updated global stiffness matrix is ​​then substituted into the dynamic equations to continue nonlinear dynamic calculations, thereby obtaining the transient displacement response data and the stress tensor data.

[0008] Furthermore, in the high-load-bearing and impact-resistant all-steel radial tire structure design method of the present invention, the step of loading the overload mechanical parameters and the road impact parameters into the finite element mesh model data includes: The internal constant pressure load and radial overload force vector are obtained by analyzing the overload mechanical parameters. The inner cavity boundary nodes and wheel axle center coordinate points are extracted from the finite element mesh model data. The internal constant pressure load is applied to the inner cavity boundary nodes, and the radial overload force vector is applied to the wheel axle center coordinate points. The road impact parameters are analyzed to obtain a rigid road digital model containing the geometric coordinates of the protruding obstacles on the surface. The tread nodes around the finite element mesh model data are extracted, and the contact state of the rigid road digital model in contact with the tread nodes is set as the contact constraint condition. The internal constant pressure load, the radial overload force vector, and the contact constraint conditions are combined as mechanical and impact loading conditions and loaded into the finite element mesh model data.

[0009] Furthermore, the high load-bearing and impact-resistant all-steel radial tire structure design method of the present invention, wherein the step of extracting the first transient stress distribution data corresponding to the coordinate interval of the belt layer endpoint and the second strain energy density sequence corresponding to the coordinate interval of the tire sidewall based on the transient displacement response data and the stress tensor data, and extracting the coordinates of stress concentration nodes based on the first transient stress distribution data and the second strain energy density sequence, includes: Numerical low-pass filtering is performed on the transient displacement response data and the stress tensor data to filter out high-frequency oscillation noise data, and the filtered tensor data array is output. The stress values ​​located in the coordinate interval of the endpoint of the belt layer are extracted from the filtered tensor data array and spliced ​​together to form the first transient stress distribution data. Extract strain energy density values ​​located in the tire side coordinate interval from the filtered tensor data array and splice them into the second strain energy density sequence. Obtain a preset reference stress safety threshold and a preset reference strain energy safety threshold. Compare the stress values ​​in the first transient stress distribution data with the reference stress safety threshold. Compare the strain energy density values ​​in the second strain energy density sequence with the reference strain energy safety threshold. Extract and combine the corresponding node coordinates of values ​​greater than the reference stress safety threshold and values ​​greater than the reference strain energy safety threshold to form the stress concentration node coordinates.

[0010] Furthermore, the high load-bearing and impact-resistant all-steel radial tire structure design method of the present invention, wherein the first transient stress distribution data, the second strain energy density sequence, and the stress concentration node coordinates are input into a preset optimization objective function, and the first dimensional parameter and the second dimensional parameter are iteratively updated using the preset optimization objective function, includes: The preset optimization objective function includes the sub-objective of minimizing the maximum stress value at the end of the belt layer and the sub-objective of minimizing the peak strain energy density value in the tire sidewall coordinate range. Extract the current cord angle and width of the transition reinforcement layer digital model as the first size parameter of the current iteration step, and extract the current height dimension and rubber hardness of the auxiliary support block digital model as the second size parameter of the current iteration step; The first size parameter and the second size parameter of the current iteration step are input into the preset optimization objective function for numerical evaluation. When the numerical evaluation result does not reach the preset numerical convergence standard, the first size parameter and the second size parameter are cross-mutated according to the genetic optimization algorithm to generate a new generation of variable data array. The corresponding node coordinates in the finite element mesh model data are updated according to the new generation variable data array. The nonlinear dynamic solution calculation is retried using the finite element mesh model data with updated coordinates, and the updated transient stress distribution data and updated strain energy density sequence of the next iteration are output.

[0011] Furthermore, the high load-bearing and impact-resistant all-steel radial tire structure design method of the present invention, wherein the step of performing cross-mutation calculations on the first size parameter and the second size parameter according to a genetic optimization algorithm to generate a new generation of variable data array includes: Within the genetic optimization algorithm, weighting coefficients are set to control the evolution of the first size parameter and the second size parameter; The first numerical adjustment amount is output according to the weighting coefficient. The first numerical adjustment amount is used to change the cord angle and width corresponding to the first size parameter. The stress transfer matrix corresponding to the coordinate interval of the belt layer endpoint is calculated according to the changed cord angle and width. The stress tensor remainder term overflowing from the stress transfer matrix is ​​imported into the spatial coordinate interval. The second numerical adjustment amount is output according to the weight coefficient. The height dimension and rubber hardness corresponding to the second size parameter are changed by the second numerical adjustment amount. The equivalent stiffness matrix of the auxiliary support block digital model resisting the stress tensor remainder is calculated according to the changed height dimension and rubber hardness. The equivalent stiffness matrix is ​​superimposed on the coordinates of the stress concentration node to calculate the new generation variable data array.

[0012] Furthermore, the high load-bearing and impact-resistant all-steel radial tire structure design method of the present invention, wherein the step of iteratively updating the first dimension parameter and the second dimension parameter using the preset optimization objective function to output the final dimension parameter matrix includes: Collect the numerical sequence of residual force vector, the numerical sequence of maximum stress at the end of the belt layer, and the numerical sequence of peak strain energy density in the tire sidewall coordinate range from the output of the preset optimization objective function in multiple consecutive iteration steps. Obtain a preset Newton iteration convergence threshold, a preset reference stress safety threshold, and a preset reference strain energy safety threshold. Compare the residual force vector numerical sequence with the preset Newton iteration convergence threshold, compare the maximum stress numerical sequence at the end of the belt layer with the preset reference stress safety threshold, and compare the peak strain energy density numerical sequence in the tire sidewall coordinate interval with the preset reference strain energy safety threshold. When the residual force vector numerical sequence is determined to be less than or equal to the preset Newton iteration convergence threshold, and the maximum stress numerical sequence at the end of the belt layer is less than or equal to the preset reference stress safety threshold, and the peak strain energy density numerical sequence in the tire sidewall coordinate interval is less than or equal to the preset reference strain energy safety threshold, an iteration termination command is issued, and the variable data array corresponding to the current iteration step is locked as the final size parameter matrix.

[0013] Furthermore, the high load-bearing and impact-resistant all-steel radial tire structure design method of the present invention, wherein generating a target all-steel radial tire digital model file based on the final size parameter matrix includes: The final dimension parameter matrix is ​​decomposed using a matrix analytical algorithm to extract target cord angle data, target width dimension data, target height dimension data, and target rubber hardness data. The target cord angle data, target width dimension data, target height dimension data, and target rubber hardness data are collectively referred to as the design variable optimization values. Read the preset blank drawing template data, and write the optimized design variable values ​​into the attribute slots of the blank drawing template data according to the standard computer-aided manufacturing digital drawing protocol specification; Execute the assembly and packaging command, compile and convert the blank drawing template data after writing the optimized values ​​of the design variables, and generate the target all-steel radial tire digital model file; The target all-steel radial tire digital model file is transferred to a terminal non-volatile storage medium for storage.

[0014] Secondly, the present invention provides a high-load-bearing and impact-resistant all-steel radial tire structure design system, applied to the aforementioned high-load-bearing and impact-resistant all-steel radial tire structure design method, comprising: The data acquisition module is used to collect geometric contour coordinate data, material constitutive parameters, overload mechanical parameters, and road impact parameters. The model building module is used to construct a parametric geometric model based on the geometric contour coordinate data and the material constitutive parameters. Within the parametric geometric model, spatial coordinate intervals, belt layer endpoint coordinate intervals, and tire sidewall coordinate intervals are defined. Within the spatial coordinate intervals, a digital model of the transition reinforcement layer with a first size parameter and a digital model of the auxiliary support block with a second size parameter are configured. The parametric geometric model configured with the digital model of the transition reinforcement layer and the digital model of the auxiliary support block are discretized, and finite element mesh model data with node coordinate sequence and material property array are output. The dynamics solution module is used to load the overload mechanical parameters and the road impact parameters into the finite element mesh model data, perform nonlinear dynamics solution calculations on the finite element mesh model data after loading the overload mechanical parameters and the road impact parameters, and output transient displacement response data and stress tensor data. The data extraction module is used to extract the first transient stress distribution data corresponding to the coordinate interval of the end point of the belt layer and the second strain energy density sequence corresponding to the coordinate interval of the tire sidewall based on the transient displacement response data and the stress tensor data, and to extract the coordinates of the stress concentration nodes based on the first transient stress distribution data and the second strain energy density sequence. The iterative optimization module is used to input the first transient stress distribution data, the second strain energy density sequence, and the coordinates of the stress concentration nodes into a preset optimization objective function, iteratively update the first size parameters and the second size parameters using the preset optimization objective function, output the final size parameter matrix, and generate a target all-steel radial tire digital model file based on the final size parameter matrix.

[0015] Beneficial effects of this invention: When performing nonlinear dynamic calculations on the finite element mesh model data after applying overload mechanical parameters and road impact parameters, the mesh distortion data corresponding to the stress node coordinate sequence is continuously extracted. Because heavy-duty vehicle tires undergo extreme nonlinear physical deformation under overload conditions due to road impact, dynamic mesh reconstruction rule parameters are triggered when the mesh distortion data exceeds a set mesh distortion threshold. Subsequently, the stress node coordinate sequence is recalculated according to the nonlinear geometric large deformation rule, and the updated stress node coordinate sequence is used to reassemble the global stiffness matrix. The updated global stiffness matrix is ​​then substituted into the dynamic equations to continue the nonlinear dynamic calculation. This processing logic coordinates the dynamic mesh reconstruction rule with the real-time update mechanism of the global stiffness matrix, effectively overcoming the defect in existing technologies where mesh distortion causes the simulation system's stiffness matrix to exhibit a singular state, leading to the interruption of the equation solver calculation. This ensures that the nonlinear dynamic calculation continues to converge, ultimately continuously outputting transient displacement response data and stress tensor data that conform to the objective force distribution law.

[0016] After obtaining objective and accurate transient displacement response data and stress tensor data, numerical low-pass filtering is performed on the transient displacement response data and stress tensor data to filter out high-frequency oscillation noise data and output a smooth filtered tensor data array. From the filtered tensor data array, the first transient stress distribution data located in the coordinate interval of the belt layer endpoint and the second strain energy density sequence located in the tire sidewall coordinate interval are extracted. The stress values ​​in the first transient stress distribution data are compared with a preset benchmark stress safety threshold, and the strain energy density values ​​in the second strain energy density sequence are also compared with the preset benchmark strain energy safety threshold. The corresponding node coordinates with values ​​greater than the aforementioned safety thresholds are then extracted and combined into stress concentration node coordinates. The transient mechanical changes during physical deformation are objectively transformed into quantifiable numerical comparison logic. The aforementioned comparison and extraction method overcomes the shortcomings of existing finite element models in accurately calculating transient mechanical changes during complex deformation processes, enabling precise calibration of the locations of extreme nonlinear physical deformations inside the tire. This provides accurate three-dimensional spatial coordinate data guidance for subsequent targeted optimization of impact-resistant structures.

[0017] The first transient stress distribution data, the second strain energy density sequence, and the coordinates of stress concentration nodes are input into a preset optimization objective function. When the numerical evaluation result of the preset optimization objective function fails to meet the preset numerical convergence criterion, a genetic optimization algorithm is used to perform cross-mutation calculations on the first and second size parameters. The change in the first size parameter reconstructs the stress transfer matrix corresponding to the coordinate interval of the belt layer endpoint, objectively altering the physical transfer path of external road impact energy. The change in the second size parameter calculates the equivalent stiffness matrix of the auxiliary support block digital model resisting the stress tensor remainder. By superimposing the equivalent stiffness matrix onto the coordinates of the stress concentration nodes, the local deformation amplitude is effectively attenuated at the underlying mechanical logic. The aforementioned iterative method deeply coordinates the stiffness orientation adjustment of the transition reinforcement layer digital model with the physical support strength of the auxiliary support block digital model, effectively overcoming the shortcomings of traditional analysis schemes that require repeated modification of structural dimensions and multiple runs of simulation calculation programs. The numerical calculation can automatically output the final size parameter matrix that meets the mechanical safety boundary conditions. The target all-steel radial tire digital model file generated according to the final size parameter matrix possesses excellent high load-bearing capacity and road impact resistance physical performance. Attached Figure Description

[0018] To more clearly illustrate the technical solution of the present invention, the drawings used in the embodiments will be briefly introduced below. Obviously, those skilled in the art can obtain other drawings based on the drawings without creative effort.

[0019] Figure 1 This is a flowchart illustrating a high-load-bearing and impact-resistant all-steel radial tire structure design method provided by the present invention. Detailed Implementation

[0020] To make the technical solution of the present invention clearer, the present invention will be clearly and completely described below with reference to specific embodiments and corresponding drawings. Obviously, the described embodiments are only a part of the embodiments of the present invention, and not all of them. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0021] Firstly, please refer to Figure 1 This invention provides a high-load-bearing and impact-resistant all-steel radial tire structure design method, comprising: Step 1: Collect geometric contour coordinate data, material constitutive parameters, overload mechanical parameters, and road impact parameters; Step 2: Construct a parametric geometric model based on the geometric contour coordinate data and the material constitutive parameters. Define a spatial coordinate interval, a belt layer endpoint coordinate interval, and a tire sidewall coordinate interval within the parametric geometric model. Configure a digital model of the transition reinforcement layer with a first size parameter and a digital model of the auxiliary support block with a second size parameter within the spatial coordinate interval. Discretize the parametric geometric model configured with the digital model of the transition reinforcement layer and the digital model of the auxiliary support block to output finite element mesh model data with node coordinate sequence and material property array. Step 3: Load the overload mechanical parameters and the road impact parameters into the finite element mesh model data, perform nonlinear dynamic solution calculation on the finite element mesh model data after loading the overload mechanical parameters and the road impact parameters, and output transient displacement response data and stress tensor data; Step 4: Based on the transient displacement response data and the stress tensor data, extract the first transient stress distribution data corresponding to the coordinate interval of the belt layer endpoint and the second strain energy density sequence corresponding to the coordinate interval of the tire sidewall, and extract the coordinates of the stress concentration nodes based on the first transient stress distribution data and the second strain energy density sequence. Step 5: Input the first transient stress distribution data, the second strain energy density sequence, and the stress concentration node coordinates into a preset optimization objective function. Use the preset optimization objective function to iteratively update the first size parameter and the second size parameter, output the final size parameter matrix, and generate the target all-steel radial tire digital model file based on the final size parameter matrix.

[0022] In the technical application scenario of heavy-duty vehicle all-steel radial tires facing overload conditions and complex road surface impacts, the research and development process relies on computer-aided engineering and finite element simulation technology to perform impact-resistant structural design. The system's underlying layer collects geometric contour coordinate data, material constitutive parameters, overload mechanical parameters, and road impact parameters. The model building component constructs a parametric geometric model based on the geometric contour coordinate data and material constitutive parameters. Within the parametric geometric model, spatial coordinate intervals, belt layer endpoint coordinate intervals, and sidewall coordinate intervals are defined. Within the spatial coordinate intervals, a digital model of the transition reinforcement layer with a first-dimensional parameter and a digital model of the auxiliary support block with a second-dimensional parameter are configured. The parametric geometric model containing the transition reinforcement layer and auxiliary support block digital models is discretized to generate and output finite element mesh model data with node coordinate sequences and material property arrays.

[0023] The system loads overload mechanical parameters and road impact parameters onto the finite element mesh model data. It then performs nonlinear dynamic calculations on the finite element mesh model data after loading these external parameters, objectively outputting transient displacement response data and stress tensor data. Based on the transient displacement response data and stress tensor data, it extracts the first transient stress distribution data corresponding to the coordinate interval of the belt layer endpoints and the second strain energy density sequence corresponding to the coordinate interval of the tire sidewall. Based on the first transient stress distribution data and the second strain energy density sequence, it extracts the coordinates of stress concentration nodes. The first transient stress distribution data, the second strain energy density sequence, and the coordinates of stress concentration nodes are uniformly input into a preset optimization objective function. The preset optimization objective function is then used to iteratively update the first and second size parameters. The final size parameter matrix is ​​calculated and output, and finally, a digital model file of the target all-steel radial tire is generated based on the final size parameter matrix.

[0024] Regarding the steps of constructing the parametric geometric model and applying mechanical boundary conditions, the underlying spatial topology division and force constraint mechanism are specifically manifested as follows: The boundary coordinates of the carcass ply and belt layer are extracted from the geometric contour coordinate data. A first spatial coordinate interval is defined between the outer perimeter of the carcass ply boundary coordinates and the inner perimeter of the belt layer boundary coordinates. A digital model of the transition reinforcement layer, loaded with first-dimensional parameters, is set within the first spatial coordinate interval; the first-dimensional parameters include initial cord angle variables and initial width variables defining the mesh topology orientation. Simultaneously, the boundary coordinates of the inner liner layer are extracted from the geometric contour coordinate data. A second spatial coordinate interval is defined between the outer perimeter of the inner liner layer boundary coordinates and the inner perimeter of the carcass ply boundary coordinates. A digital model of the auxiliary support block, loaded with second-dimensional parameters, is set within the second spatial coordinate interval; the second-dimensional parameters include initial hardness variables and initial height variables. The first and second spatial coordinate intervals are collectively referred to as the spatial coordinate interval.

[0025] Subsequently, the endpoint coordinates of the belt layer and the sidewall coordinates were extracted from the geometric contour coordinate data. Based on these coordinates, the endpoint coordinate intervals of the belt layer and the sidewall coordinate intervals were precisely defined. A finite element mesh generation program was then performed on the parametric geometric model containing the digital models of the transition reinforcement layer and the auxiliary support blocks to discretize it, outputting finite element mesh model data containing node coordinate sequences and material property arrays. During the model loading stage, the overload mechanical parameters were analyzed to obtain the internal constant pressure load and the radial overload force vector. The inner cavity boundary nodes and the wheel axle center coordinate points were extracted from the finite element mesh model data. The internal constant pressure load was applied to the inner cavity boundary nodes, and the radial overload force vector was applied to the wheel axle center coordinate points. The road impact parameters were analyzed to obtain a rigid road surface digital model containing the geometric coordinates of protruding obstacles. The tread nodes on the periphery of the finite element mesh model data were extracted, and the physical collision state of the rigid road surface digital model contacting the tread nodes was set as the mechanical contact constraint condition. Finally, the internal constant pressure load, the radial overload force vector, and the contact constraint conditions were used together as the mechanical and impact loading conditions input into the simulation calculation layer.

[0026] In the aforementioned nonlinear dynamic solution calculation of the finite element mesh model data, the computational processing logic is as follows: A node set filtering algorithm is used to identify the nodes corresponding to the digital models of the transition reinforcement layer and the auxiliary support blocks in the finite element mesh model data. The three-dimensional coordinates corresponding to the identified nodes are extracted as a sequence of force-bearing node coordinates. During the nonlinear dynamic solution calculation, the computational core continuously extracts the mesh distortion data corresponding to the force-bearing node coordinate sequence; the system acquires the input dynamic mesh reconstruction rule parameters and the set mesh distortion threshold.

[0027] When the mesh distortion data exceeds the mesh distortion threshold, it objectively indicates that the physical shape of the local mesh element has reversed and collapsed. The system immediately triggers the dynamic mesh reconstruction rule parameters, recalculates and outputs the updated stress node coordinate sequence according to the nonlinear geometric large deformation rule. The updated stress node coordinate sequence is used to reassemble the global stiffness matrix, and the updated global stiffness matrix is ​​output. The updated global stiffness matrix is ​​then seamlessly substituted into the dynamic equation to continue the nonlinear dynamic solution calculation until the integration time step is completed, and the transient displacement response data and stress tensor data containing global deformation information are obtained.

[0028] Regarding the steps for extracting the coordinates of stress concentration nodes, the underlying data signal processing and physical threshold determination logic is as follows: Numerical low-pass filtering is performed on the acquired transient displacement response data and stress tensor data. High-frequency oscillation noise data is filtered out using an algorithm, and a smoothed tensor data array is output. Stress values ​​located within the coordinate intervals of the belt layer endpoints are extracted from the filtered tensor data array and concatenated according to time node order to form the first transient stress distribution data. Simultaneously, strain energy density values ​​located within the tire sidewall coordinate interval are extracted and concatenated to form the second strain energy density sequence.

[0029] The system is pre-configured with a reference stress safety threshold of 150 MPa and a reference strain energy safety threshold of 5.0 MJ / m³. The calculation module objectively compares the stress values ​​in the first transient stress distribution data with the 150 MPa reference stress safety threshold, and compares the strain energy density values ​​in the second strain energy density sequence with the 5.0 MJ / m³ reference strain energy safety threshold. Through numerical comparison logic, the system accurately extracts the corresponding node coordinates for reference stress safety thresholds greater than 150 MPa and reference strain energy safety thresholds greater than 5.0 MJ / m³, and combines the extracted over-limit node coordinates into stress concentration node coordinates in the three-dimensional topological space.

[0030] In the aforementioned step of iteratively updating dimensional parameters using a preset optimization objective function, the specific data processing path of the genetic optimization algorithm is as follows: The system inputs the sub-objectives of minimizing the maximum stress at the belt layer endpoints and minimizing the peak strain energy density in the sidewall coordinate interval into the preset optimization objective function. The system extracts the current cord angle and width of the transition reinforcement layer digital model as the first dimensional parameter for the current iteration step, and simultaneously extracts the current height and rubber hardness of the auxiliary support block digital model as the second dimensional parameter for the current iteration step. The first and second dimensional parameters of the current iteration step are then input into the preset optimization objective function for multi-dimensional numerical evaluation.

[0031] When the numerical evaluation results fail to meet the preset numerical convergence criteria, weighting coefficients are set within the genetic optimization algorithm to control the co-evolution of the first and second size parameters. Based on these weighting coefficients, a first numerical adjustment is output, which changes the cord angle and width corresponding to the first size parameter. The stress transfer matrix corresponding to the endpoint coordinate interval of the belt layer is recalculated based on the changed cord angle and width, and the stress tensor remainder that overflows from the stress transfer matrix and cannot be absorbed is imported into the spatial coordinate interval. Simultaneously, a second numerical adjustment is output based on the weighting coefficients, which changes the height dimension and rubber hardness corresponding to the second size parameter. The equivalent stiffness matrix of the auxiliary support block digital model resisting the stress tensor remainder is calculated based on the changed height dimension and rubber hardness. The equivalent stiffness matrix is ​​directly superimposed at the stress concentration node coordinates at the mathematical level to perform cross-mutation calculations, thereby generating a new generation of variable data arrays containing multiple combinations of structural parameter adjustments. The system automatically updates the corresponding node coordinates in the finite element mesh model data based on the new generation variable data array, and uses the updated finite element mesh model data to re-trigger the nonlinear dynamics solution closed loop, outputting the updated transient stress distribution data and updated strain energy density sequence for the next iteration step.

[0032] Regarding the steps of outputting the final dimensional parameter matrix and generating the target all-steel radial tire digital model file, the underlying multi-dimensional constraint convergence verification and drawing compilation logic specifically involves: collecting the residual force vector numerical sequence, the maximum stress numerical sequence at the belt layer endpoints, and the peak strain energy density numerical sequence in the sidewall coordinate interval from the preset optimization objective function output in multiple consecutive iteration steps. The system's underlying layer retrieves a Newton iteration convergence threshold set to 0.001, a reference stress safety threshold set to 150 MPa, and a reference strain energy safety threshold set to 5.0 MJ / m³. The data comparison component compares the residual force vector numerical sequence with the Newton iteration convergence threshold of 0.001, simultaneously compares the maximum stress numerical sequence at the belt layer endpoints with the reference stress safety threshold of 150 MPa, and compares the peak strain energy density numerical sequence in the sidewall coordinate interval with the reference strain energy safety threshold of 5.0 MJ / m³.

[0033] When the residual force vector numerical sequence is less than or equal to the Newton iteration convergence threshold of 0.001, the maximum stress numerical sequence at the belt layer endpoint is less than or equal to the reference stress safety threshold of 150 MPa, and the peak strain energy density numerical sequence in the tire sidewall coordinate interval is less than or equal to the reference strain energy safety threshold of 5.0 MJ / m³, it objectively indicates that the current structural spatial topology has completely absorbed and resisted external impact loads. The system then issues an iteration termination command, locking the variable data array corresponding to the current iteration step in memory as the final dimension parameter matrix. The built-in matrix analysis algorithm is used to decompose the final dimension parameter matrix, extracting the target cord angle data, target width dimension data, target height dimension data, and target rubber hardness data; the target cord angle data, target width dimension data, target height dimension data, and target rubber hardness data are merged and collectively referred to as the design variable optimization value. The preset blank drawing template data is read through the system data bus, and according to the standard computer-aided manufacturing digital drawing protocol specification, the design variable optimization value is input as a deterministic parameter and written into the attribute slot of the blank drawing template data. The drawing output calculation module is called to execute the assembly and packaging instructions. The blank drawing template data after writing the design variable optimization values ​​is pushed into the compilation and conversion engine for low-level format reconstruction, generating a target all-steel radial tire digital model file that can directly guide industrial manufacturing. Finally, the target all-steel radial tire digital model file is losslessly transmitted to the terminal non-volatile storage medium for long-term physical storage.

[0034] The dynamics solution module calls the microprocessor unit in the hardware architecture to perform nonlinear dynamics solution equations on the finite element mesh model data after loading overload mechanical parameters and road impact parameters:

[0035] In solving equations for nonlinear dynamics, This represents the global quality matrix of the finite element mesh model data. This represents the global damping matrix of the finite element mesh model data; This represents the global stiffness matrix of the finite element mesh model data. This represents the transient acceleration response matrix at a given time point. This represents the transient velocity response matrix at a given time point. This represents the output transient displacement response data; This represents the time-varying column vector of external loads generated by the transformation of overload mechanical parameters; This represents the time-varying impact force column vector generated from the transformation of road impact parameters. The iterative optimization module calls the computation scheduling core, substituting the first transient stress distribution data, the second strain energy density sequence, and the coordinates of stress concentration nodes into the preset optimization objective function:

[0036] In the preset optimization objective function, This represents the numerical evaluation result output by the preset optimization objective function; Indicates the first dimension parameter; Indicates the second dimension parameter; This represents the stress weighting coefficient of the belt layer built into the memory; This represents the weighting coefficient of the tire side strain energy built into the memory; This represents the maximum stress value at the end of the belt layer extracted from the first transient stress distribution data; This indicates the preset reference stress safety threshold; This represents the peak strain energy density value extracted from the second strain energy density sequence within the tire sidewall coordinate range; This indicates the preset reference strain energy safety threshold.

[0037] Secondly, the present invention provides a high-load-bearing and impact-resistant all-steel radial tire structure design system, applied to the aforementioned high-load-bearing and impact-resistant all-steel radial tire structure design method, comprising: The data acquisition module is used to collect geometric contour coordinate data, material constitutive parameters, overload mechanical parameters, and road impact parameters. The model building module is used to construct a parametric geometric model based on the geometric contour coordinate data and the material constitutive parameters. Within the parametric geometric model, spatial coordinate intervals, belt layer endpoint coordinate intervals, and tire sidewall coordinate intervals are defined. Within the spatial coordinate intervals, a digital model of the transition reinforcement layer with a first size parameter and a digital model of the auxiliary support block with a second size parameter are configured. The parametric geometric model configured with the digital model of the transition reinforcement layer and the digital model of the auxiliary support block are discretized, and finite element mesh model data with node coordinate sequence and material property array are output. The dynamics solution module is used to load the overload mechanical parameters and the road impact parameters into the finite element mesh model data, perform nonlinear dynamics solution calculations on the finite element mesh model data after loading the overload mechanical parameters and the road impact parameters, and output transient displacement response data and stress tensor data. The data extraction module is used to extract the first transient stress distribution data corresponding to the coordinate interval of the end point of the belt layer and the second strain energy density sequence corresponding to the coordinate interval of the tire sidewall based on the transient displacement response data and the stress tensor data, and to extract the coordinates of the stress concentration nodes based on the first transient stress distribution data and the second strain energy density sequence. The iterative optimization module is used to input the first transient stress distribution data, the second strain energy density sequence, and the coordinates of the stress concentration nodes into a preset optimization objective function, iteratively update the first size parameters and the second size parameters using the preset optimization objective function, output the final size parameter matrix, and generate a target all-steel radial tire digital model file based on the final size parameter matrix.

Claims

1. A design method for a high-load-bearing and impact-resistant all-steel radial tire structure, characterized in that, include: Collect geometric contour coordinate data, material constitutive parameters, surcharge mechanical parameters, and road impact parameters; A parametric geometric model is constructed based on the geometric contour coordinate data and the material constitutive parameters. Within the parametric geometric model, spatial coordinate intervals, belt layer endpoint coordinate intervals, and tire sidewall coordinate intervals are defined. Within the spatial coordinate intervals, a digital model of a transition reinforcement layer with a first size parameter and a digital model of an auxiliary support block with a second size parameter are configured. The parametric geometric model configured with the digital model of the transition reinforcement layer and the digital model of the auxiliary support block are discretized to output finite element mesh model data with node coordinate sequences and material property arrays. The overload mechanical parameters and the road impact parameters are loaded into the finite element mesh model data. Nonlinear dynamics calculations are performed on the finite element mesh model data after loading the overload mechanical parameters and the road impact parameters, and transient displacement response data and stress tensor data are output. Based on the transient displacement response data and the stress tensor data, extract the first transient stress distribution data corresponding to the coordinate interval of the belt layer endpoint and the second strain energy density sequence corresponding to the coordinate interval of the tire sidewall, and extract the coordinates of the stress concentration nodes based on the first transient stress distribution data and the second strain energy density sequence. The first transient stress distribution data, the second strain energy density sequence, and the coordinates of the stress concentration nodes are input into a preset optimization objective function. The first size parameter and the second size parameter are iteratively updated using the preset optimization objective function, and the final size parameter matrix is ​​output. The target all-steel radial tire digital model file is generated based on the final size parameter matrix.

2. The high load-bearing and impact-resistant all-steel radial tire structure design method according to claim 1, characterized in that, The process involves defining spatial coordinate intervals, belt layer endpoint coordinate intervals, and sidewall coordinate intervals within the parametric geometric model. Within these spatial coordinate intervals, a digital model of a transition reinforcement layer with a first size parameter and a digital model of an auxiliary support block with a second size parameter are configured. The parametric geometric model containing these two digital models is then discretized to output finite element mesh model data with a sequence of node coordinates and an array of material properties, including: Extract the boundary coordinates of the carcass ply and the belt layer from the geometric contour coordinate data. Define a first spatial coordinate interval between the outer periphery of the carcass ply boundary coordinates and the inner periphery of the belt layer boundary coordinates. Set the digital model of the transition reinforcement layer loaded with the first size parameters within the first spatial coordinate interval. The first size parameters include an initial cord angle variable and an initial width variable. Extract the inner liner boundary coordinates from the geometric contour coordinate data, define a second spatial coordinate interval between the outer perimeter of the inner liner boundary coordinates and the inner perimeter of the carcass ply boundary coordinates, set the auxiliary support block digital model loaded with the second size parameter in the second spatial coordinate interval, the second size parameter includes an initial hardness variable and an initial height variable, and the spatial coordinate interval includes the first spatial coordinate interval and the second spatial coordinate interval; Extract the belt layer endpoint coordinates and sidewall coordinates from the geometric contour coordinate data, and define the belt layer endpoint coordinate interval and the sidewall coordinate interval according to the belt layer endpoint coordinates and the sidewall coordinates respectively; The parametric geometric model loaded with the digital model of the transition reinforcement layer and the digital model of the auxiliary support block is discretized by a finite element mesh generation program, and the output is the finite element mesh model data containing the node coordinate sequence and the material property array.

3. The high load-bearing and impact-resistant all-steel radial tire structure design method according to claim 2, characterized in that, The nonlinear dynamic solution calculation is performed on the finite element mesh model data after loading the overload mechanical parameters and the road impact parameters, and the transient displacement response data and stress tensor data are output, including: Identify the nodes corresponding to the digital model of the transition reinforcement layer and the digital model of the auxiliary support block in the finite element mesh model data, extract the coordinates corresponding to the identified nodes as a sequence of force-bearing node coordinates, and continuously extract the mesh distortion data corresponding to the sequence of force-bearing node coordinates during the nonlinear dynamic solution calculation. Obtain the input dynamic mesh reconstruction rule parameters and mesh distortion threshold; When the mesh distortion data exceeds the mesh distortion threshold, the dynamic mesh reconstruction rule parameters are triggered, and the updated sequence of force-bearing node coordinates is recalculated and output according to the nonlinear geometric large deformation rule. The updated sequence of force-bearing node coordinates is used to assemble the global stiffness matrix, which is then output as an updated global stiffness matrix. The updated global stiffness matrix is ​​then substituted into the dynamic equations to continue nonlinear dynamic calculations, thereby obtaining the transient displacement response data and the stress tensor data.

4. The high load-bearing and impact-resistant all-steel radial tire structure design method according to claim 3, characterized in that, The loading of the overload mechanical parameters and the road impact parameters into the finite element mesh model data includes: The internal constant pressure load and radial overload force vector are obtained by analyzing the overload mechanical parameters. The inner cavity boundary nodes and wheel axle center coordinate points are extracted from the finite element mesh model data. The internal constant pressure load is applied to the inner cavity boundary nodes, and the radial overload force vector is applied to the wheel axle center coordinate points. The road impact parameters are analyzed to obtain a rigid road digital model containing the geometric coordinates of the protruding obstacles on the surface. The tread nodes around the finite element mesh model data are extracted, and the contact state of the rigid road digital model in contact with the tread nodes is set as the contact constraint condition. The internal constant pressure load, the radial overload force vector, and the contact constraint conditions are combined as mechanical and impact loading conditions and loaded into the finite element mesh model data.

5. The high load-bearing and impact-resistant all-steel radial tire structure design method according to claim 4, characterized in that, The step of extracting first transient stress distribution data corresponding to the coordinate interval of the belt layer endpoint and second strain energy density sequence corresponding to the coordinate interval of the tire sidewall based on the transient displacement response data and the stress tensor data, and extracting stress concentration node coordinates based on the first transient stress distribution data and the second strain energy density sequence, includes: Numerical low-pass filtering is performed on the transient displacement response data and the stress tensor data to filter out high-frequency oscillation noise data, and the filtered tensor data array is output. The stress values ​​located in the coordinate interval of the endpoint of the belt layer are extracted from the filtered tensor data array and spliced ​​together to form the first transient stress distribution data. Extract strain energy density values ​​located in the tire side coordinate interval from the filtered tensor data array and splice them into the second strain energy density sequence. Obtain a preset reference stress safety threshold and a preset reference strain energy safety threshold. Compare the stress values ​​in the first transient stress distribution data with the reference stress safety threshold. Compare the strain energy density values ​​in the second strain energy density sequence with the reference strain energy safety threshold. Extract and combine the corresponding node coordinates of values ​​greater than the reference stress safety threshold and values ​​greater than the reference strain energy safety threshold to form the stress concentration node coordinates.

6. The high load-bearing and impact-resistant all-steel radial tire structure design method according to claim 5, characterized in that, The step of inputting the first transient stress distribution data, the second strain energy density sequence, and the stress concentration node coordinates into a preset optimization objective function, and iteratively updating the first size parameter and the second size parameter using the preset optimization objective function, includes: The preset optimization objective function includes the sub-objective of minimizing the maximum stress value at the end of the belt layer and the sub-objective of minimizing the peak strain energy density value in the tire sidewall coordinate range. Extract the current cord angle and width of the transition reinforcement layer digital model as the first size parameter of the current iteration step, and extract the current height dimension and rubber hardness of the auxiliary support block digital model as the second size parameter of the current iteration step; The first size parameter and the second size parameter of the current iteration step are input into the preset optimization objective function for numerical evaluation. When the numerical evaluation result does not reach the preset numerical convergence standard, the first size parameter and the second size parameter are cross-mutated according to the genetic optimization algorithm to generate a new generation of variable data array. The corresponding node coordinates in the finite element mesh model data are updated according to the new generation variable data array. The nonlinear dynamic solution calculation is retried using the finite element mesh model data with updated coordinates, and the updated transient stress distribution data and updated strain energy density sequence of the next iteration are output.

7. The high load-bearing and impact-resistant all-steel radial tire structure design method according to claim 6, characterized in that, The step of performing crossover and mutation calculations on the first size parameter and the second size parameter according to the genetic optimization algorithm to generate a new generation of variable data array includes: Within the genetic optimization algorithm, weighting coefficients are set to control the evolution of the first size parameter and the second size parameter; The first numerical adjustment amount is output according to the weighting coefficient. The first numerical adjustment amount is used to change the cord angle and width corresponding to the first size parameter. The stress transfer matrix corresponding to the coordinate interval of the belt layer endpoint is calculated according to the changed cord angle and width. The stress tensor remainder term overflowing from the stress transfer matrix is ​​imported into the spatial coordinate interval. The second numerical adjustment amount is output according to the weight coefficient. The height dimension and rubber hardness corresponding to the second size parameter are changed by the second numerical adjustment amount. The equivalent stiffness matrix of the auxiliary support block digital model resisting the stress tensor remainder is calculated according to the changed height dimension and rubber hardness. The equivalent stiffness matrix is ​​superimposed on the coordinates of the stress concentration node to calculate the new generation variable data array.

8. The high load-bearing and impact-resistant all-steel radial tire structure design method according to claim 7, characterized in that, The step of iteratively updating the first size parameter and the second size parameter using the preset optimization objective function to output the final size parameter matrix includes: Collect the numerical sequence of residual force vector, the numerical sequence of maximum stress at the end of the belt layer, and the numerical sequence of peak strain energy density in the tire sidewall coordinate range from the output of the preset optimization objective function in multiple consecutive iteration steps. Obtain a preset Newton iteration convergence threshold, a preset reference stress safety threshold, and a preset reference strain energy safety threshold. Compare the residual force vector numerical sequence with the preset Newton iteration convergence threshold, compare the maximum stress numerical sequence at the end of the belt layer with the preset reference stress safety threshold, and compare the peak strain energy density numerical sequence in the tire sidewall coordinate interval with the preset reference strain energy safety threshold. When the residual force vector numerical sequence is determined to be less than or equal to the preset Newton iteration convergence threshold, and the maximum stress numerical sequence at the end of the belt layer is less than or equal to the preset reference stress safety threshold, and the peak strain energy density numerical sequence in the tire sidewall coordinate interval is less than or equal to the preset reference strain energy safety threshold, an iteration termination command is issued, and the variable data array corresponding to the current iteration step is locked as the final size parameter matrix.

9. The high load-bearing and impact-resistant all-steel radial tire structure design method according to claim 8, characterized in that, The step of generating the target all-steel radial tire digital model file based on the final size parameter matrix includes: The final dimension parameter matrix is ​​decomposed using a matrix analytical algorithm to extract target cord angle data, target width dimension data, target height dimension data, and target rubber hardness data. The target cord angle data, target width dimension data, target height dimension data, and target rubber hardness data are collectively referred to as the design variable optimization values. Read the preset blank drawing template data, and write the optimized design variable values ​​into the attribute slots of the blank drawing template data according to the standard computer-aided manufacturing digital drawing protocol specification; Execute the assembly and packaging command, compile and convert the blank drawing template data after writing the optimized values ​​of the design variables, and generate the target all-steel radial tire digital model file; The target all-steel radial tire digital model file is transferred to a terminal non-volatile storage medium for storage.

10. A high-load-bearing and impact-resistant all-steel radial tire structure design system, applied to the high-load-bearing and impact-resistant all-steel radial tire structure design method as described in any one of claims 1 to 9, characterized in that, include: The data acquisition module is used to collect geometric contour coordinate data, material constitutive parameters, overload mechanical parameters, and road impact parameters. The model building module is used to construct a parametric geometric model based on the geometric contour coordinate data and the material constitutive parameters. Within the parametric geometric model, spatial coordinate intervals, belt layer endpoint coordinate intervals, and tire sidewall coordinate intervals are defined. Within the spatial coordinate intervals, a digital model of the transition reinforcement layer with a first size parameter and a digital model of the auxiliary support block with a second size parameter are configured. The parametric geometric model configured with the digital model of the transition reinforcement layer and the digital model of the auxiliary support block are discretized, and finite element mesh model data with node coordinate sequence and material property array are output. The dynamics solution module is used to load the overload mechanical parameters and the road impact parameters into the finite element mesh model data, perform nonlinear dynamics solution calculations on the finite element mesh model data after loading the overload mechanical parameters and the road impact parameters, and output transient displacement response data and stress tensor data. The data extraction module is used to extract the first transient stress distribution data corresponding to the coordinate interval of the end point of the belt layer and the second strain energy density sequence corresponding to the coordinate interval of the tire sidewall based on the transient displacement response data and the stress tensor data, and to extract the coordinates of the stress concentration nodes based on the first transient stress distribution data and the second strain energy density sequence. The iterative optimization module is used to input the first transient stress distribution data, the second strain energy density sequence, and the coordinates of the stress concentration nodes into a preset optimization objective function, iteratively update the first size parameters and the second size parameters using the preset optimization objective function, output the final size parameter matrix, and generate a target all-steel radial tire digital model file based on the final size parameter matrix.

Images