CN114297892A - Method, device and program for detailed analysis of contribution of tire deformation mode to rolling resistance - Google Patents

Method, device and program for detailed analysis of contribution of tire deformation mode to rolling resistance Download PDF

Info

Publication number
CN114297892A
CN114297892A CN202111613566.8A CN202111613566A CN114297892A CN 114297892 A CN114297892 A CN 114297892A CN 202111613566 A CN202111613566 A CN 202111613566A CN 114297892 A CN114297892 A CN 114297892A
Authority
CN
China
Prior art keywords
tire
rolling resistance
analysis
contribution
deformation mode
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111613566.8A
Other languages
Chinese (zh)
Inventor
王剑波
崔志博
胡德斌
郭磊磊
李进
王毅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhongce Rubber Group Co Ltd
Original Assignee
Zhongce Rubber Group Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhongce Rubber Group Co Ltd filed Critical Zhongce Rubber Group Co Ltd
Priority to CN202111613566.8A priority Critical patent/CN114297892A/en
Publication of CN114297892A publication Critical patent/CN114297892A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Tires In General (AREA)

Abstract

The invention relates to the technical field of intelligent design of tires, in particular to a method, equipment and a computer program for refining and analyzing contribution of tire deformation modes to rolling resistance. The method comprises the following steps: 1) establishing a tire finite element model, and carrying out load analysis; 2) extracting real stress and strain values of units at the same position and different angles around the tire, and fitting sigma and epsilon values by using a 100-order trigonometric series; 3) calculating the energy loss density e of the rubber material unit in the corresponding deformation mode by using the fitted parameter data; 4) and drawing a tire section unit grid graph according to the node coordinates and the unit composition by using a program, and filling different colors into the units according to the energy loss density value of each unit in a single deformation mode. The invention realizes the accurate positioning of the part which has larger contribution to the rolling resistance, obtains the contribution of the deformation mode of each part of the tire to the rolling resistance, and refines the analysis of the tire rolling resistance.

Description

Method, device and program for detailed analysis of contribution of tire deformation mode to rolling resistance
Technical Field
The invention relates to the technical field of intelligent design of tires, in particular to a method, equipment and a computer program for refining and analyzing contribution of tire deformation modes to rolling resistance.
Background
With the implementation of the European tire labeling law, the industry development trend of green tires and the requirement of endurance mileage of electric vehicles, the low-rolling-resistance tire has long been a consensus pursuit in the tire industry, the rolling resistance is one of key indexes for measuring the performance of the tire, and the tire cannot be sold if the rolling resistance performance does not reach the standard. The factors influencing the rolling resistance of the tire mainly include materials, structures, road conditions, use conditions and the like, wherein the road conditions and the use conditions belong to uncontrollable factors, and for tire design engineers, the structures and the materials are two main entry points for regulating and controlling the rolling resistance of the tire. The method mainly comprises a test method and a finite element simulation method, the test method can accurately obtain the rolling resistance value of the tire, but the method needs to manufacture and test the tire, has high cost and long period, cannot obtain the internal deformation state of the tire, and can find the generation reason of high or low rolling resistance by carrying out multiple tests on tires with different structures. Compared with the finite element simulation method, the method has the advantages of low cost and short period, and can visually obtain the deformation and stress state of the tire, thereby being convenient for analyzing the influence of the structure and the material on the deformation and the rolling resistance of the tire. Currently, simulation analysis research on the rolling resistance of a tire mainly focuses on the proportion of each component of the tire to the overall rolling resistance of the tire and the rolling resistance prediction after the tire structure or material is adjusted, however, the guidance effect on a tire design engineer is limited, specific positions with large rolling resistance of the tire cannot be accurately positioned, for example, the tire tread structure covers the running surface of the whole tire, only the proportion of the rolling resistance of the whole tread can be analyzed, and the proportion (the tread center or the tread edge) in the tread is large cannot be judged.
Meanwhile, the stress state of a rubber material in a tire is very complex, six strain (stress) components are contained according to the traditional mechanical analysis, the current rolling resistance analysis method cannot judge the rolling resistance contribution of various deformation modes to each part of the tire, so that the material test cannot be guided (the loss tangent value of the rubber material is generally required to be tested by calculating the rolling resistance, the loss tangent value is related to the deformation mode (stretching or shearing)), a main deformation mode with large rolling resistance contribution of each part is required to be obtained, and the performance of the material is tested according to the deformation mode, so that the method has very important guidance for tire material design engineers.
Disclosure of Invention
The present invention is intended to solve the above problems of the prior art, that is, the present simulation analysis research on the rolling resistance of the tire mainly focuses on the ratio of each component of the tire to the overall rolling resistance of the tire, and the rolling resistance prediction after adjusting the tire structure or material, however, the guidance of the tire design engineer is limited, and the specific location where the rolling resistance of the tire is large cannot be precisely located, for example, the tire tread structure covers the running surface of the whole tire, and at present, only the rolling resistance ratio of the whole tread can be analyzed, and which portion (tread center or tread edge) in the tread contributes greatly cannot be judged. Meanwhile, the stress state of the rubber material in the tire is very complex, six strain (stress) components are contained according to the traditional mechanical analysis, and the current rolling resistance analysis method cannot judge the rolling resistance contribution of various deformation modes to each part of the tire, so that the material test cannot be guided (the loss tangent value of the rubber material is generally required to be tested when the rolling resistance is calculated, and the loss tangent value is related to the deformation mode (stretching or shearing)).
And further provides a method for analyzing the contribution of the tire deformation mode to the rolling resistance in a refining manner, so that the accurate positioning of the part with larger contribution to the rolling resistance is realized, the contribution of the deformation mode of each part of the tire to the rolling resistance is obtained, and the analysis of the rolling resistance of the tire is refined.
In order to achieve the purpose, the invention adopts the following technical scheme:
a method for refining analysis of the contribution of tire deformation modes to rolling resistance, the method comprising the steps of:
1) establishing a tire finite element model, and carrying out load analysis: carrying out grid division on a design drawing (material distribution diagram) of the tire, giving material properties, applying rated inflation pressure, carrying out inflation analysis, applying rated load to the tire on the basis of the inflation analysis, carrying out tire load analysis, and outputting stress and strain values of all rubber material units;
2) extracting real stress and strain values of units at the same position and different angles in a circle of the tire by using a program, respectively recording the real stress and strain values as a sigma value and an epsilon value, fitting the sigma and epsilon values by using a 100-order trigonometric series by taking an angle as an x value and taking x as a radian as a unit:
Figure RE-GDA0003510448080000021
Figure RE-GDA0003510448080000022
3) calculating the energy loss density e of the rubber material unit in the corresponding deformation mode by using the fitted parameter data, wherein the tan delta value is the loss tangent of the rubber material,
Figure RE-GDA0003510448080000023
4) and drawing a tire section unit grid graph according to the node coordinates and the unit composition by using a program, wherein the units are closed triangles or quadrangles, and filling the units with different colors according to the energy loss density values of the units under the single deformation mode, which are calculated in the third step, by using the program.
Preferably, in the step 1), the tire material distribution map is subjected to meshing division, the tire material distribution map is divided into triangular or quadrilateral units, the framework material is divided into 2-node one-dimensional units, and material attributes are given to each component material to establish a tire finite element model.
Preferably, the step 1) applies an air pressure of 0.25MPa to the inner surface of the tire inner liner, and performs inflation analysis using Abaqus software.
Preferably, in the step 1), the tire section is rotated by 360 degrees and divided in the circumferential direction in addition to the two-dimensional inflation analysis, one circle is divided into 110 sections, the tire rim is fixed, the road surface is moved in the rim direction by applying 5881N load to the road surface, the load analysis is performed, and the stress and strain values of all the rubber material units are output.
Preferably, the step 2) utilizes a Python writing program to extract the real stress and strain values of the unit at the same position and different angles in a circle of the tire.
Preferably, the stress and strain values of step 2) include 11 direction, 22 direction, 33 direction, 12 direction or 21 direction, 13 direction or 31 direction and 23 direction or 32 direction as shown in fig. 5.
Preferably, in the step 4), a cell grid map of the tire section is drawn according to the node coordinates and the cell composition by utilizing a matplotlib program in Python.
Further, the present invention also discloses a computer device, which comprises a memory, a processor and a computer program stored on the memory, wherein the processor executes the computer program to implement the steps of the method according to any one of the above technical solutions.
Further, the present invention also discloses a computer readable storage medium, on which a computer program or instructions are stored, which when executed by a processor implement the steps of the method according to any of the above technical solutions.
A computer program product comprising computer programs or instructions which, when executed by a processor, carry out the steps of the method according to any one of the preceding claims.
By adopting the technical scheme, the invention realizes the accurate positioning of the part which has larger contribution to the rolling resistance, obtains the contribution of the deformation mode of each part of the tire to the rolling resistance and refines the analysis of the rolling resistance of the tire.
Drawings
FIG. 1 is a 21550R15 tire material distribution diagram;
FIG. 2 is a 21550R15 tire section grid and material components;
FIG. 3 shows the inflated deformation results of 21550R15 tire;
FIG. 4 shows the results of 21550R15 tire three-dimensional circumferential meshing;
FIG. 5 shows the three-dimensional load results and direction indications for the 21550R15 tire;
FIG. 6 is the strain output after 21550R15 tire load analysis;
FIG. 7 shows the 11-way stress in a 21550R15 tire for cell number 3 in one revolution;
FIG. 8 shows the 11-way strain for cell No. 3 in the 21550R15 tire over one week;
FIG. 9 is a trigonometric series fit of 11 direction strain for unit No. 3 in a 21550R15 tire over one week;
FIG. 10 is a graph of the rolling resistance distribution of a 21550R15 tire resulting from a 11-way deformation;
FIG. 11 is a graph of the rolling resistance distribution of a 21550R15 tire resulting from a 22-way deformation;
FIG. 12 is a graph of the rolling resistance distribution of a 21550R15 tire resulting from a 33-way deformation;
FIG. 13 is a graph of the rolling resistance distribution of a 21550R15 tire resulting from a 12-way deformation;
FIG. 14 is a graph of the rolling resistance distribution of the 21550R15 tire resulting from a 13-way deformation;
fig. 15 shows the rolling resistance distribution of the 21550R15 tire due to the 23-direction deformation.
Detailed Description
The invention will be described in further detail below with reference to the accompanying drawings: the present embodiment is implemented on the premise of the technical solution of the present invention, and a detailed implementation is given, but the scope of the present invention is not limited to the following embodiments.
Take 21550R15 tire as an example:
firstly, a tire material distribution diagram (shown in figure 1) is subjected to meshing division to be divided into triangular or quadrilateral units (shown in figure 2), framework materials are divided into 2-node one-dimensional units, material attributes (shown in table 1) are given to each component material to establish a tire finite element model, 0.25MPa air pressure is applied to the inner surface of a tire inner liner, and inflation analysis is carried out by utilizing Abaqus software, wherein the result is shown in figure 3. In addition to the two-dimensional inflation analysis, the tire section was rotated 360 degrees and divided into sections in the circumferential direction, one circle was divided into 110 sections (fig. 4), the tire rim was fixed, a 5881N load was applied to the road surface, the road surface was moved in the rim direction, the load analysis was performed (fig. 5), and the stress and strain values of all the rubber material cells were output (fig. 6).
Secondly, extracting the real stress and strain values of the unit (including 11 direction, 22 direction, 33 direction, 12 direction (equal to 21 direction), 13 direction or 31 direction and 23 direction or 32 direction, as shown in fig. 5) of the tire at the same position and different angles in a circle by using a Python writing program, respectively recording the values as a sigma value (stress) and a epsilon value (strain), drawing a relation curve of x, sigma and epsilon (fig. 7 and 8) by taking the angle as a x value (x is in radian), and fitting the sigma and epsilon values by using a 100-order trigonometric series:
Figure RE-GDA0003510448080000041
Figure RE-GDA0003510448080000042
the first 10 th order parameters for a strain fit in the 11 direction for cell No. 3 are shown in table 2.
The third step: using the fitted 100 σ s obtained for each deformation mode cell datancParameter, 100 εnsParameter, 100 εncParameter and 100 εnsThe energy loss density e of the rubber material unit in the corresponding deformation mode is calculated by parameters, the tan delta value is the loss tangent of the rubber material of the corresponding unit (shown in the table 1),
Figure RE-GDA0003510448080000043
partial results of the calculations are shown in table 3.
TABLE 121550R 15 materials Properties for tire components
Figure RE-GDA0003510448080000044
Figure RE-GDA0003510448080000051
TABLE 2 partial trigonometric series fitting coefficients for 21550R15 tire No. 3 element strain in circumferential direction 11
Order of the scale εnc εns
n=1 0.048677128161277006 -0.01178
n=2 0.005762 0.034562
n=3 -0.02496 0.006944
n=4 -0.01306 -0.01652
n=5 0.010247 -0.0155
n=6 0.009846 0.008189
n=7 -0.00664 0.005613
n=8 -0.00468 -0.00443
n=9 0.00364 -0.00312
n=10 0.000865 0.003969
TABLE 3 rolling loss values for 6 deformation modes of 21550R15 tire section unit
Figure RE-GDA0003510448080000052
Figure RE-GDA0003510448080000061
The fourth step: and drawing a tire section cell grid graph according to the node coordinates and the cell composition by utilizing a matplotlib program in Python, wherein the cells are closed triangles or quadrangles, and filling the cells with different colors according to the energy loss density values of each cell calculated in the third step under a single deformation mode by utilizing the program, and the result is shown in fig. 11-15.
Through the detailed analysis of the tire rolling resistance realized by the above technology, the main cause of the rolling resistance generated at different parts of the tire can be intuitively obtained, for example, as can be seen from fig. 13, the rolling resistance of the bead and the sidewall caused by the 33 direction is very small, that is, the 33 direction deformation with too large amplitude basically does not occur, the direction of reducing the rolling resistance of the bead region is to reduce the shear deformation of the 13 direction (fig. 14), and by comparing fig. 14 with the other 5 results, the main cause of the rolling resistance generated at the sidewall region is the shear deformation of the 13 direction, so that the shear mode should be adopted when measuring the tan δ value of the sidewall rubber material. Through the above brief analysis, it can be known that the rolling resistance of the tire can be rapidly and inexpensively refined and analyzed by using the patent technology, and the product development efficiency and the product structural design direction can be greatly improved.
The above description is only a preferred embodiment of the present invention, and these embodiments are based on different implementations of the present invention, and the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (10)

1. A method for refining the analysis of the contribution of tire deformation modes to rolling resistance, characterized in that it comprises the following steps:
1) establishing a tire finite element model, and carrying out load analysis: carrying out grid division on a design drawing (material distribution diagram) of the tire, giving material properties, applying rated inflation pressure, carrying out inflation analysis, applying rated load to the tire on the basis of the inflation analysis, carrying out tire load analysis, and outputting stress and strain values of all rubber material units;
2) extracting real stress and strain values of units at the same position and different angles in a circle of the tire by using a program, respectively recording the real stress and strain values as a sigma value and an epsilon value, fitting the sigma and epsilon values by using a 100-order trigonometric series by taking an angle as an x value and taking x as a radian as a unit:
Figure DEST_PATH_IMAGE002
Figure DEST_PATH_IMAGE004
3) calculating the energy loss density e of the rubber material unit in the corresponding deformation mode by using the fitted parameter data, wherein the tan delta value is the loss tangent of the rubber material,
Figure DEST_PATH_IMAGE006
4) drawing a tire section unit grid graph according to the node coordinates and the unit composition by using a program, wherein the units are closed triangles or quadrangles, and filling the units with different colors according to the energy loss density values of the units under the single deformation mode, which are calculated in the step 3), by using the program.
2. The method for refining and analyzing the contribution of the tire deformation mode to the rolling resistance as claimed in claim 1, wherein the step 1) is to perform meshing on the tire material distribution diagram, divide the tire material distribution diagram into triangular or quadrilateral units, divide the framework material into 2-node one-dimensional units, and endow material properties to each component material to establish a tire finite element model.
3. The method for the refined analysis of the contribution of the tire deformation mode to the rolling resistance as claimed in claim 1, wherein the step 1) is to apply an air pressure of 0.25MPa to the inner surface of the tire liner and perform the inflation analysis by using the Abaqus software.
4. The method for analyzing the contribution of the tire deformation mode to the rolling resistance of claim 1, wherein the step 1) is to rotate the tire section by 360 degrees on the basis of the two-dimensional inflation analysis, divide the tire section in the circumferential direction, divide one circumference into 110 sections, fix the tire rim, apply 5881N load to the road surface, move the road surface in the direction of the rim, analyze the load, and output the stress and strain values of all rubber material units.
5. The method for the detailed analysis of the contribution of the tire deformation mode to the rolling resistance according to claim 1, wherein the step 2) utilizes a Python programming to extract the real stress and strain values of the cells at the same position and different angles around the tire.
6. A method for refining analysis of contribution of tire deformation mode to rolling resistance as claimed in claim 1, wherein step 2) stress and strain values comprise 11 direction, 22 direction, 33 direction, 12 direction or 21 direction, 13 direction or 31 direction and 23 direction or 32 direction as shown in fig. 5.
7. The method for the detailed analysis of the contribution of the tire deformation mode to the rolling resistance according to claim 1, wherein the step 4) is performed by utilizing a matplotlib program in Python to draw a cell grid diagram of the tire section according to the node coordinates and the cell composition.
8. A computer device comprising a memory, a processor and a computer program stored on the memory, characterized in that the processor executes the computer program to implement the steps of the method of any of claims 1-7.
9. A computer-readable storage medium, on which a computer program or instructions are stored, which, when executed by a processor, carry out the steps of the method of any one of claims 1 to 7.
10. A computer program product comprising a computer program or instructions for implementing the steps of the method of any one of claims 1 to 7 when executed by a processor.
CN202111613566.8A 2021-12-27 2021-12-27 Method, device and program for detailed analysis of contribution of tire deformation mode to rolling resistance Pending CN114297892A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111613566.8A CN114297892A (en) 2021-12-27 2021-12-27 Method, device and program for detailed analysis of contribution of tire deformation mode to rolling resistance

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111613566.8A CN114297892A (en) 2021-12-27 2021-12-27 Method, device and program for detailed analysis of contribution of tire deformation mode to rolling resistance

Publications (1)

Publication Number Publication Date
CN114297892A true CN114297892A (en) 2022-04-08

Family

ID=80970266

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111613566.8A Pending CN114297892A (en) 2021-12-27 2021-12-27 Method, device and program for detailed analysis of contribution of tire deformation mode to rolling resistance

Country Status (1)

Country Link
CN (1) CN114297892A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114889372A (en) * 2022-05-23 2022-08-12 中策橡胶集团股份有限公司 Method for judging pulling and pressing state of rubber material at tire bead part, application and program product

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114889372A (en) * 2022-05-23 2022-08-12 中策橡胶集团股份有限公司 Method for judging pulling and pressing state of rubber material at tire bead part, application and program product
CN114889372B (en) * 2022-05-23 2023-08-08 中策橡胶集团股份有限公司 Judgment method and application of tension and compression state of rubber material at tire bead part

Similar Documents

Publication Publication Date Title
Li et al. State of charge prediction algorithm of lithium-ion battery based on PSO-SVR cross validation
CN107024663B (en) Lithium battery screening technique based on charging curve feature KPCA cluster
US9020798B2 (en) Real-time capable battery cell simulation
CN111873804B (en) Electric automobile working condition data generation method and energy consumption evaluation working condition generation method
CN104677640B (en) A kind of fuel cell hybrid car economic testing method
CN105115692B (en) The pneumatic Forecasting Methodology that a kind of CFD numerical simulations are combined with wind tunnel test
CN108732508B (en) Real-time estimation method for lithium ion battery capacity
CN110568360B (en) Lithium battery aging diagnosis method based on fuzzy logic algorithm
CN114297892A (en) Method, device and program for detailed analysis of contribution of tire deformation mode to rolling resistance
Einhorn et al. Comparison of electrical battery models using a numerically optimized parameterization method
CN102829948B (en) Method for simulating large-proportion scaling of ultrasonic wind tunnel experiment of grid fin
CN112763180B (en) Method for quickly predicting aerodynamic noise of pantograph model of high-speed train in acoustic wind tunnel
CN106446343A (en) Method for automatically extracting parameterized profile line of straight-line blade of radial-flow impeller
CN112163372A (en) SOC estimation method of power battery
CN104390794A (en) Method of predicating flat-pavement mechanical properties of tyre based on test data of rotary drum test bed
CN111191186A (en) Multi-cell filtering method for positioning position of mobile robot in production workshop
CN114297891A (en) Method, device and computer program for calculating ratio of bending energy to compression energy of tread
CN114707374A (en) Method, application and program product for calculating influence of rubber material modulus on tire rolling resistance
CN104615835A (en) Engine intercooler analyzing method
CN115659625A (en) Method for predicting rolling resistance of periodic pattern tire, application and computer program product
CN111611646A (en) Rapid design method for aerodynamic configuration of aircraft forebody
CN116593915A (en) Method, system, storage medium and terminal for monitoring voltage in battery module
CN113125974B (en) Method and device for detecting lithium precipitation of battery
Zhou et al. Modelling and aerodynamic prediction of fish bone active camber morphing airfoil
CN107764493A (en) A kind of acquisition methods and acquisition device for vibrating long duration test PSD

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination