Simulation method, system, medium and program for lateral inclination, lateral deviation and longitudinal slip composite slip characteristics of tire

Abstract

The invention relates to the technical field of tire simulation design, in particular to a method, a system, a medium and a program for simulating lateral inclination, lateral deviation and longitudinal slip composite slip characteristics of a tire. According to the method, through the steps of finite element pretreatment, loading analysis modeling, accelerated rolling analysis, steady-state angular velocity extraction and roll-to-roll longitudinal slip characteristic simulation analysis, three-dimensional loading analysis, friction attribute setting and rolling analysis are carried out by utilizing an explicit solver and a self-compiled program, the angular velocity of steady-state rolling is extracted, and the angular velocity of roll-to-roll is obtained. And simulation and data analysis are carried out on tire characteristics under different slip rates. The system comprises a pre-processing module, a modeling module, a rolling analysis module, a characteristic simulation module and a solving and post-processing module, and can realize efficient simulation analysis of the tire roll-to-side deviation longitudinal slip characteristics. Through the method and the system, the dynamic performance of the tire under different working conditions can be accurately simulated, and reliable data support is provided for tire design and optimization.

Description

Technical Field

[0001] The present invention relates to the technical field of tire simulation design, and particularly to a simulation method, system, medium and program for the combined slip characteristics of tire roll, sideslip and longitudinal slip. Background Art

[0002] Tire virtual sampling is an indispensable and important part of vehicle virtual development. The tire dynamics model is the carrier of tire virtual sampling and an important link in realizing vehicle virtual development. The combined slip characteristics of tire roll, sideslip and longitudinal slip are the mechanical characteristics exhibited by the tire under the conditions of roll, sideslip and longitudinal slip, and are an important component in establishing the tire dynamics model. Currently, the data of the combined slip characteristics of tire roll, sideslip and longitudinal slip can be obtained through actual measurement and simulation. However, actual measurement requires the manufacture of physical tires and professional six-component force test equipment, with a long test cycle and high cost. With the development of simulation technology and the improvement of computer hardware level, obtaining the data of the combined slip characteristics of tire roll, sideslip and longitudinal slip through simulation has become the main means of tire virtual sampling.

[0003] In the prior art, the patent application for invention (Publication No.: CN114462282A, Publication Date: 20220510) discloses a finite element simulation method for the steady-state characteristics of tire roll, sideslip and longitudinal slip. This method uses the standard solver for three-dimensional loading and solving the steady-state angular velocity. Especially in the case of a more complex model, problems of non-convergence are likely to occur during the simulation process. Moreover, during dynamic simulation, it is necessary to extract the tire sinkage and update the road surface information again. In addition, the roll, sideslip and longitudinal slip solution adopts the scanning method, and the time required to calculate a complete working condition is long. If the finite element simulation model is large, the calculation efficiency is low and the simulation cost is high. As the calculation time increases, the probability of mesh distortion also increases. The patent application for invention of the applicant (Publication No.: CN117973143A, Publication Date: 20240503) discloses a finite element simulation method for tire sideslip characteristics. This method splits the sideslip simulation model file into multiple files for submission and calculation, avoiding the problem that the mesh is prone to distortion during long-term calculation, and having relatively high simulation accuracy and simulation efficiency. The patent application for invention of the applicant (Publication No.: CN118278248A, Publication Date: 20240702) discloses a finite element simulation method for tire longitudinal slip characteristics. This method splits the model file into multiple files for submission and calculation to perform longitudinal slip characteristic simulation analysis. During the calculation process, the mesh is not easily distorted, the simulation accuracy is high, and the simulation efficiency is also taken into account. However, the above two latter methods have not considered the conditions of simultaneous roll, sideslip and longitudinal slip. Summary of the Invention

[0004] To solve the above technical problems, the present invention provides a finite element simulation method for the combined slip characteristics of tire roll, sideslip, and longitudinal slip. Through steps such as finite element preprocessing, two-dimensional inflation and three-dimensional loading modeling analysis, rolling simulation analysis, combined slip characteristics simulation modeling analysis of tire roll, sideslip, and longitudinal slip, and simulation result extraction, the combined slip characteristics of the tire roll, sideslip, and longitudinal slip are simulated. In the three-dimensional loading, rolling, and combined slip finite element simulation of tire roll, sideslip, and longitudinal slip, explicit explicit analysis is adopted to avoid non-convergence problems. Since restart analysis is adopted in the analysis of tire roll, sideslip, and longitudinal slip, there is no need to update the road surface information, which simplifies the simulation method. The friction property between the tire and the road surface adopts a friction model obtained by fitting the friction test of the tread rubber material, and different slip rate values can be flexibly selected independently, so as to achieve the purpose of simulating the combined slip characteristics of tire roll, sideslip, and longitudinal slip with high precision and high efficiency.

[0005] To achieve the above object, the present invention adopts the following technical solutions: A finite element simulation method for the combined slip characteristics of tire roll, sideslip, and longitudinal slip, the method comprising the following steps: 1) Finite element preprocessing: 1.1) Draw the tire material distribution map; 1.2) Conduct fine mesh division to distinguish each component of the tire; 2) Loading analysis modeling: 2.1) Conduct two-dimensional static inflation simulation modeling and analysis; 2.2) Generate a three-dimensional model through the SMG command, set the tire central axis as the Y axis, the road surface loading along the Z axis direction, rotate the road surface by γ angle around the X axis, which is the roll angle of the tire, and rotate by α angle around the rotation axis perpendicular to the road surface within the plane of the road surface, which is the sideslip angle of the tire, establish a local coordinate system, the direction along the movement of the road surface is the X' axis, and the Z' axis is perpendicular to the road surface; 2.3) Import information including each component, constraint, surface, and road surface, set the mass and moment of inertia of the wheel and the road surface, and set the friction properties between the rim and the tire, and between the tire and the road surface; 2.4) Conduct three-dimensional loading analysis; 3) Accelerated rolling analysis step: Based on the loading analysis, establish an accelerated rolling analysis step, release the constraint in the 5 direction of the rim reference point, that is, the rotation along the Y axis, fix other directions, release the constraints in the X' axis and Z' axis directions of the road surface reference point, fix other directions, apply a moving speed to the horizontal direction of the road surface, that is, the X' axis direction, and set the output of historical variables; conduct restart analysis on the accelerated rolling analysis, establish a uniform rolling analysis step, apply a set moving speed to the horizontal direction of the road surface, that is, the X' axis direction, and set the output of historical variables; 4) Extract the angular velocity of the steady-state rolling analysis step, and average the angular velocity output by the historical variables to obtain the steady-state angular velocity; 5) Simulation modeling and analysis of roll, sideslip, and longitudinal slip characteristics: 5.1) Based on the steady-state rolling analysis, perform a restart analysis. Keep the road surface moving uniformly along the X' axis, apply an angular velocity to the 5 direction of the rim reference point, that is, the tire rotation direction around the Y axis. Set the amplitude function of the angular velocity to gradually change from the steady-state angular velocity to the given angular velocity from 0 s to t1 s, and keep the angular velocity unchanged from t1 s to t2 s. Set the historical variable output, and set the output time interval to t s, then output 1 / t data per second; 5.2) Select a series of angular velocities and generate a series of roll, sideslip, and longitudinal slip analysis input files using a program; 5.3) Submit for batch processing and solve, and automatically extract and post-process the results using a self-program.

[0006] Preferably, in step 1), perform a fine grid division on the tire material distribution map to distinguish each part of the tire.

[0007] Preferably, in step 2.1), all skeleton materials adopt the rebar model, and a lift equation is added to the belt. By designing the belt drum diameter and the design outer diameter, the elongation rate and the belt radius are calculated.

[0008] Preferably, in step 2.3), add a friction model between the road surface and the tire. The friction model includes slip velocity and contact pressure parameters, and set the friction coefficient range between the tire and the rim to be between 0.1 and 0.5.

[0009] Preferably, in step 2.4), the three-dimensional loading uses an explicit solver for simulation analysis.

[0010] Preferably, in step 3), the calculation time range for the accelerated rolling analysis is set between 0.05 s and 0.2 s, and the calculation time range for the uniform rolling analysis is set between 0.1 s and 0.3 s. The uniform rolling analysis needs to output the angular velocity variable of the rim reference point.

[0011] Preferably, in step 5.1), the calculation time t1 range for the roll, sideslip, and longitudinal slip analysis step is set between 0.05 s and 0.15 s, and t2 is between 0.05 s and 0.5 s. In the historical variable output, set 4 variables, namely RF1, RF2, RM3, and RF3 of the rim reference point. Among them, RF3 is used to reference whether the load is stable during the calculation process, and the maximum value of t is set to 0.0025 s, that is, at least 400 points are output per second.

[0012] Preferably, in step 5.2), among the selected slip ratios, the maximum and minimum slip ratios, i.e., the critical values, shall be included, and the others shall be selected according to the actual situation. Each slip ratio corresponds to a calculation file.

[0013] Preferably, in step 5.3), data extraction is performed on each calculation result file. The time history of the roll - camber - longitudinal slip analysis step is extracted as the last tL s of data. Each variable has tL / t data points. Averaging these data points gives the variable value corresponding to the slip ratio. The variable values include RF1, RF2, RM3, RF3, etc. Taking the slip ratio as the abscissa and the variable value RF1 (longitudinal force) as the ordinate to plot points, a scatter plot of the longitudinal force varying with the slip ratio can be obtained, or connecting with a straight line or curve in ascending order of the slip ratio to obtain a curve graph of the longitudinal force varying with the slip ratio. Similarly, scatter plots or curve graphs of variable values such as RF2 and RM3 varying with the slip ratio can be plotted.

[0014] Furthermore, the present invention also discloses a simulation method for the combined slip characteristics of tire roll - camber - longitudinal slip, which is applied in the simulation of the combined slip characteristics of tire roll - camber - longitudinal slip.

[0015] Furthermore, the present invention also discloses a system for finite - element simulation analysis of the combined slip characteristics of tire roll - camber - longitudinal slip. This system implements the above - mentioned method, and the system includes: a) A pre - processing module for drawing a tire material distribution map and performing fine mesh division; b) A modeling module for performing two - dimensional static inflation simulation modeling, three - dimensional loading analysis, and friction property setting; c) A rolling analysis module for performing accelerated rolling analysis, uniform rolling analysis, and angular velocity extraction; d) A characteristic simulation module for performing roll - camber - longitudinal slip characteristic simulation modeling and analysis based on the steady - state angular velocity and generating the corresponding input file; e) A solution and post - processing module for batch - processing solution, automatically extracting analysis result data for post - processing, and generating a scatter plot or curve graph between the variable value and the slip ratio.

[0016] Furthermore, the present invention also discloses a computer - readable storage medium, on which a computer program or instruction is stored. When the computer program or instruction is executed by a processor, the above - mentioned method is implemented.

[0017] Furthermore, the present invention also discloses a computer program product, including a computer program or instruction. When the computer program or instruction is executed by a processor, the above - mentioned method is implemented.

[0018] Due to the adoption of the above technical solution, the present invention takes into account all components of the tire, and precise measurement and fitting are carried out for each rubber and skeleton material. By using the friction model obtained through actual measurement and identification, the simulation accuracy is guaranteed. At the same time, by selecting multiple different slip rates for simulation, the number of slip rate analyses and the analysis time of analysis steps can be controlled, avoiding the problem of mesh distortion that may be caused by too long simulation time, and also taking into account the simulation efficiency, resulting in a significant improvement in simulation efficiency and a substantial reduction in simulation cost.

[0019] The beneficial effects of the present invention are as follows: 1. The present invention adopts the explicit display analysis method in three-dimensional loading and rolling analysis, avoiding the problem of non-convergence of the model during the simulation process. Since the restart analysis is adopted in the roll, sideslip and longitudinal slip analysis, there is no need to update the road surface information, simplifying the simulation method; 2. The present invention takes into account all components of the tire and the friction model between the road surface and the tire, with high simulation accuracy; 3. Except for several key slip rate values, other slip rate values can be freely and flexibly selected for simulation, with high simulation efficiency; 4. Since individual slip rate values are used for calculation respectively, the calculation time is short, thus avoiding the problem of mesh distortion that may be caused by too long simulation time. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] Figure 1 is the flowchart of the method of the present invention.

[0021] Figure 2 is the two-dimensional finite element simulation model of the 205 / 55R16 tire specification.

[0022] Figure 3 is the three-dimensional finite element simulation model of the 205 / 55R16 tire specification.

[0023] Figure 4 is the contact pressure distribution at a certain moment when the slip rate is -20% under the conditions of 230 kPa air pressure, 6500 N load, roll angle of -5° (IA - 5), and sideslip angle of -2° (SA - 2) for the 205 / 55R16 tire specification.

[0024] Figure 5 is the scatter plot of the simulation results with the slip rate as the abscissa and the longitudinal force as the ordinate under the conditions of 230 kPa air pressure, 6500 N load, IA - 5, and SA - 2 for the 205 / 55R16 tire specification.

[0025] Figure 6It is a scatter plot of simulation results with slip ratio as the abscissa and lateral force as the ordinate for the 205 / 55R16 tire specification under the conditions of 230 kPa air pressure, 6500 N load, and IA-5, SA-2.

[0026] Figure 7 It is a scatter plot of simulation results with slip ratio as the abscissa and aligning torque as the ordinate for the 205 / 55R16 tire specification under the conditions of 230 kPa air pressure, 6500 N load, and IA-5, SA-2. Specific implementation mode

[0027] The method of the present invention is used for the simulation analysis of the combined slip characteristics of roll, sideslip, and longitudinal slip of various different types of tires.

[0028] The present invention is further described below through an example. This embodiment is implemented on the premise of the technical solution of the present invention, but the protection scope of the present invention is not limited to the following embodiments.

[0029] Taking the 205 / 55R16 tire specification as an example, a finite element simulation analysis model of the combined slip characteristics of roll, sideslip, and longitudinal slip with a roll angle of -5° and a sideslip angle of -2° under 230 kPa air pressure and 6500 N load is established.

[0030] Use AutoCAD to draw the material distribution diagram, export it as a dxf file, and then import it into the Hypermesh software for mesh division. The minimum element size is preferably controlled above 0.5 mm. After division, export the inp file, and then import it into the ABAQUS / CAE software for two-dimensional inflation modeling analysis. It has longitudinal grooves, an inflation pressure of 230 kPa, the rubber material uses the yeoh model, the skeleton material uses the rebar model, and each rubber and skeleton component is distinguished and material parameters are measured. Through the belt drum diameter, designed outer diameter, and the dimensions of the material distribution diagram, the elongation rates of belt layer 1 and belt layer 2 can be calculated to be 4.92% and 4.90% respectively, and the belt radii are 291.75 and 292.75 respectively. The rim uses an analytical rigid body model. The two-dimensional inflation model is as Figure 2 shown. This model has a total of 1480 nodes and 1351 elements. Use the ABAQUS / Standard solver to perform inflation simulation under the condition of 230 kPa air pressure.

[0031] Generate a three-dimensional model through the SYMMETRIC MODEL GENERATION function, equally divide 60 cross-sections along the circumferential direction, rotate the road surface 5° around the X-axis, and then rotate 2° in the plane where the road surface is located. Establish a local coordinate system, and input the road surface information into the three-dimensional model inp file. The finite element three-dimensional model is as Figure 3 shown.

[0032] Import each component of the tire, import the constraints for the embedded skeletal materials, the road surface reference point constraints, and the rigid body constraints for the rim reference point. Input the geometric information of the rim and the road surface. Import the surfaces of the tire in contact with the rim, the inner cavity, and the rim. Set the weight of the wheel to 9.25 kg, and the moments of inertia in the X, Y, and Z directions to 110, 1200, and 110 respectively. Establish a friction model formula based on relevant literature and actual measurements, considering the influence of contact pressure and slip velocity on friction performance. Set the friction coefficient between the tire and the rim to 0.3.

[0033] Use the ABAQUS / explicit solver to perform a loading analysis under the working conditions of a load of 6500 N and 230 kPa.

[0034] Establish an accelerated rolling analysis step. Set the air pressure to 230 kPa and the load to 6500 N. Apply a moving speed in the horizontal direction of the road surface, i.e., the X' axis direction, in the local coordinate system. The speed gradually increases from 0 to 60 km / h, which is 16667 mm / s. Set the calculation time to 0.1 s. Establish a uniform rolling analysis step. Apply a constant moving speed of 60 km / h, which is 16667 mm / s, in the X' axis direction. Set the calculation time to 0.1 s. Set the historical variable output to the VR2 variable of the rim reference point. Set the time interval to 0.00125 s, i.e., output 800 points per second. After completion, submit the ABAQUS / Explicit calculation. A total of 80 VR2 variable data points are calculated, and their average value is taken to obtain a steady-state angular velocity of 54.3 rad / s.

[0035] Establish a roll, slip, and longitudinal slip analysis step. This step performs a restart analysis based on the uniform rolling analysis, maintaining the uniform motion of the road surface with a constant speed. Apply an angular velocity to the 5 direction of the rim reference point, i.e., the tire rotation direction around the Y axis. Set the amplitude function of the angular velocity to gradually change from 54.3 rad / s to the given angular velocity W from 0 s to 0.1 s, with the corresponding slip ratio SR. From 0.1 s to 0.25 s, the angular velocity W, i.e., the slip ratio SR, remains unchanged. Set the historical variable output to the four variables RF1, RF2, RM3, and RF3 of the rim reference point. Among them, RF3 is used to reference whether the load is stable during the calculation process. Set the time interval to 0.00125 s, i.e., output 800 points per second. SR takes -30%, -25%, -20%, -15%, -10%, -7.5%, -5%, -2.5%, 0%, 2.5%, 5%, 7.5%, 10%, 15%, 20%, 25%, 30%, a total of 17 points. Each point corresponds to an input file, that is, 17 input files need to be generated, which can be quickly generated through a program.

[0036] 17 input files are submitted to ABAQUS / Explicit for calculation in batches. Each input file takes 0.25 seconds to calculate, plus the rolling analysis step calculation time, a total of 4.45 seconds. If the slip rate is calculated by scanning, it will take about 15 seconds, and the total calculation time is less than one-third of the scanning method. Therefore, the calculation efficiency is significantly improved, and the mesh distortion problem that may be caused by too long simulation time is avoided. Figure 4 is the contact pressure distribution at a certain moment when the slip rate is -20%.

[0037] Data is extracted from each abaqus calculation result odb file, and the data of the time history of the lateral longitudinal slip analysis step from 0.15s to 0.25s are extracted. Each variable has a total of 80 data points. The corresponding variable value under the slip rate is obtained by averaging these 80 data points. Plot points with the slip rate as the horizontal coordinate and the longitudinal force (RF1) as the vertical coordinate to obtain a scatter plot with the slip rate as the horizontal coordinate and the longitudinal force as the vertical coordinate. Similarly, scatter plots of lateral force (RF2) and return moment (RM3) can be drawn. In addition, the vertical force RF3 results under each slip rate can also be extracted to measure whether the load control is stable. Figure 5 , Figure 6 and Figure 7 They are scatter plots of simulation results with slip rate as the horizontal axis and RF1, RF2 and RM3 as the vertical axis.

[0038] The above is a description of the embodiments of the present invention. Through the above description of the disclosed embodiments, professionals and technicians in the field can implement or use the present invention. Various modifications to these embodiments will be apparent to professionals and technicians in the field. The general principles defined herein can be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to these embodiments shown in this article, but will conform to the widest range consistent with the principles and novelties disclosed herein.

Claims

1. A finite element simulation method for tire roll, sideways and longitudinal slip composite slip characteristics, characterized in that: The method comprises the following steps: 1) Finite element pre-processing: 1.1) Draw a tire material distribution map; 1.2) Perform fine mesh division to distinguish the various components of the tire; 2) Loading analysis modeling: 2.1) Conduct two-dimensional static inflation simulation modeling and analysis; 2.2) Generate a 3D model using the SMG command, set the tire center axis as the Y axis, set the road surface load along the Z axis, rotate the road surface around the X axis by an angle of γ, which is the tire roll angle, and rotate the road surface around the rotation axis perpendicular to the road surface in the plane where the road surface is located by an angle of α, which is the tire slip angle. Establish a local coordinate system, with the direction of movement along the road surface as the X' axis, and the Z' axis perpendicular to the road surface; 2.3) Import information including components, constraints, surfaces, and road surfaces, set the mass and moment of inertia of the wheels and road surfaces, and set the friction properties between the rim and tire, and between the tire and road surfaces; 2.4) Conduct three-dimensional loading analysis; 3) Accelerate the rolling analysis steps: On the basis of loading analysis, an accelerated rolling analysis step is established, the constraint of the rim reference point 5 direction, i.e., rotation along the Y axis, is released, and other directions are fixed; the constraints of the road surface reference point in the X' and Z' axis directions are released, and other directions are fixed; a moving speed is applied to the horizontal direction of the road surface, i.e., the X' axis direction, and the historical variable output is set; a restart analysis is performed on the accelerated rolling analysis, a uniform rolling analysis step is established, a set moving speed is applied to the horizontal direction of the road surface, i.e., the X' axis direction, and the historical variable output is set; 4) Extract the angular velocity of the steady-state rolling analysis step and average the angular velocity output by the historical variables to obtain the steady-state angular velocity; 5) Simulation modeling and analysis of roll, yaw and longitudinal slip characteristics: 5.1) Perform restart analysis based on steady-state rolling analysis, keep the road surface moving at a constant speed along the X' axis, apply angular velocity to the 5 direction of the rim reference point, that is, the tire rotation direction around the Y axis, set the angular velocity amplitude function to change from the steady-state angular velocity to the given angular velocity from 0s to t1s, keep the angular velocity unchanged from t1s to t2s, set the historical variable output, set the output time interval to ts, and output 1 / t data per second; 5.2) Select a series of angular velocities and use the program to generate a series of roll, yaw and pitch analysis input files; 5.3) Submit for batch solution, and use self-programming to automatically extract data and post-process the results.

2. The finite element simulation method for tire roll, sideways and longitudinal slip composite slip characteristics according to claim 1, characterized in that: In step 1), the tire material distribution map is finely meshed to distinguish various tire components.

3. The finite element simulation method for tire roll, sideways and longitudinal slip composite slip characteristics according to claim 1, characterized in that: In step 2.1), the rebar model is used for all skeleton materials, in which the belt is added with the lift equation, and the elongation and belt radius are calculated by designing the belt drum diameter and the design outer diameter.

4. The finite element simulation method for tire roll, sideways and longitudinal slip composite slip characteristics according to claim 1, characterized in that: In step 2.3), add a friction model between the road surface and the tire. The friction model includes slip velocity and contact pressure parameters, and sets the friction coefficient between the tire and the rim to a range between 0.1 and 0.

5.

5. The finite element simulation method for tire roll, sideways and longitudinal slip composite slip characteristics according to claim 1, characterized in that: In step 2.4), the three-dimensional loading is simulated and analyzed using the explicit solver.

6. The finite element simulation method for tire roll, sideways and longitudinal slip composite slip characteristics according to claim 1, characterized in that: In step 3), the calculation time range of the accelerated rolling analysis is set to between 0.05s and 0.2s, and the calculation time range of the uniform rolling analysis is set to between 0.1s and 0.3s. The uniform rolling analysis requires the output of the angular velocity variable of the rim reference point.

7. The finite element simulation method for tire roll, sideways and longitudinal slip composite slip characteristics according to claim 1, characterized in that: In step 5.1), the calculation time of the roll, yaw and longitudinal slip analysis step is set to t1 between 0.05s and 0.15s, and t2 between 0.05s and 0.5s. In the historical variable output, 4 variables of rim reference points, RF1, RF2, RM3, and RF3, are set. RF3 is used to refer to whether the load is stable during the calculation process. The maximum setting of t is 0.0025s, that is, at least 400 points are output per second. And / or, in step 5.2), the selected slip rates must include the maximum and minimum slip rates, i.e., the critical values, and the rest are selected according to the actual situation, and each slip rate corresponds to a calculation file; And / or, in step 5.3), extract data from each calculation result file, extract the time history of the roll, yaw and longitudinal slip analysis step as the last t L s data, each variable has t L / t data points, and the corresponding variable value under the slip rate is obtained by averaging these data points. The variable values ​​include RF1, RF2, RM3, and RF3. Points are plotted with the slip rate as the horizontal coordinate and the variable value longitudinal force RF1 as the vertical coordinate to obtain a scatter plot of the longitudinal force changing with the slip rate, or they are connected with straight lines or curves in order from small to large slip rates to obtain a curve plot of the longitudinal force changing with the slip rate. Similarly, a scatter plot or curve plot of the RF2 and RM3 variable values ​​changing with the slip rate can be drawn.

8. A system for finite element simulation analysis of tire roll, sideways and longitudinal slip characteristics, characterized in that: The system implements the method described in any one of claims 1 to 7, and the system includes: a) a pre-processing module, which is used to draw a tire material distribution map and perform fine grid division; b) a modeling module, which is used to perform two-dimensional static inflation simulation modeling, three-dimensional loading analysis and friction property setting; c) a rolling analysis module, which is used to perform accelerated rolling analysis, uniform rolling analysis and angular velocity extraction; d) a characteristic simulation module, which is used to perform roll, yaw and longitudinal slip characteristic simulation modeling and analysis based on steady-state angular velocity, and generate corresponding input files; e) a solution and post-processing module, which is used to batch solve and automatically extract analysis result data for post-processing, and generate a scatter plot or curve graph between variable values ​​and slip rate.

9. A computer-readable storage medium having a computer program or instruction stored thereon, characterized in that: When the computer program or instruction is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

10. A computer program product comprising a computer program or instructions, characterized in that When the computer program or instruction is executed by a processor, the method according to any one of claims 1 to 7 is implemented.

Images