CN116432320B - Method for rapidly calculating rolling resistance of tire - Google Patents
Method for rapidly calculating rolling resistance of tire Download PDFInfo
- Publication number
- CN116432320B CN116432320B CN202310472577.1A CN202310472577A CN116432320B CN 116432320 B CN116432320 B CN 116432320B CN 202310472577 A CN202310472577 A CN 202310472577A CN 116432320 B CN116432320 B CN 116432320B
- Authority
- CN
- China
- Prior art keywords
- tire
- rubber material
- rolling resistance
- calculating
- energy density
- 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.)
- Active
Links
- 238000005096 rolling process Methods 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 28
- 239000000463 material Substances 0.000 claims abstract description 48
- 238000004458 analytical method Methods 0.000 claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims abstract description 11
- 238000004088 simulation Methods 0.000 claims abstract description 10
- 238000010586 diagram Methods 0.000 claims abstract description 5
- 238000009826 distribution Methods 0.000 claims abstract description 5
- 238000009825 accumulation Methods 0.000 claims abstract description 4
- 230000006872 improvement Effects 0.000 abstract description 4
- 238000005457 optimization Methods 0.000 abstract description 3
- 230000008569 process Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- SYHGEUNFJIGTRX-UHFFFAOYSA-N methylenedioxypyrovalerone Chemical compound C=1C=C2OCOC2=CC=1C(=O)C(CCC)N1CCCC1 SYHGEUNFJIGTRX-UHFFFAOYSA-N 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Tires In General (AREA)
Abstract
The invention discloses a method for rapidly calculating rolling resistance of a tire, which comprises the following steps: step one, carrying out load simulation analysis on a tire: meshing the material distribution diagram of the tire; giving material properties to the materials of all the parts, and establishing a tire finite element model; applying rated inflation pressure to the inner surface of the inner lining layer of the tire to perform inflation analysis; applying rated load to the tire on the basis of inflation analysis, performing tire load analysis, and outputting strain energy density values of all rubber material units; step two, extracting strain energy density values of rubber material units at the same position and at different angles around the tire, and calculating strain energy density accumulation w of each rubber material unit i The method comprises the steps of carrying out a first treatment on the surface of the Step three, calculating the energy loss h of each rubber material unit i The method comprises the steps of carrying out a first treatment on the surface of the And step four, calculating the rolling resistance value R of the tire. The method can realize the rapid calculation of the rolling resistance of the tire and provide support for the optimization or improvement of the rolling resistance of the tire.
Description
Technical Field
The invention belongs to the technical field of radial tires, and relates to a tire rolling resistance calculation method.
Background
Rolling resistance is one of the important basic indexes for measuring the performance of the tire, and with the improvement of the requirements of the electric automobile on the endurance capacity, the requirements on the rolling resistance of the tire are higher and higher. How to quickly optimize or improve the rolling resistance of a tire is an important technical problem faced by tire design engineers. The basis of the rapid optimization or improvement of the tire rolling resistance is to realize rapid prediction of the tire rolling resistance, and the measurement and evaluation of the tire rolling resistance mainly comprises a test and a finite element simulation method, wherein the test method can accurately obtain the tire rolling resistance value, but the test is needed to manufacture the tire, and the cost is high, the period is long, and compared with the finite element simulation method, the cost is low, and the period is short. At present, simulation prediction of tire rolling resistance is mainly based on a Fourier transform method, simulation result data is required to be subjected to Fourier transform and then processed, for a common car tire, 10 minutes are generally required, and for a car tire, up to 30 minutes are required, so that the defect of long data processing period exists, particularly when a program is automatically optimized, hundreds of schemes are generally required to be calculated, only tens of hours are required in the Fourier transform process, and the efficiency is low. According to the calculation principle of the tire rolling resistance, the tire rolling resistance is generated due to energy loss in the deformation process of the rubber material, and the deformation energy loss and the whole stored energy have a proportional relation, so that the rolling resistance can be calculated by using the energy of the rubber material for extracting the simulation calculation result, the rolling resistance of the tire can be calculated within a few seconds, and the tire rolling resistance has enough precision compared with a Fourier transformation method.
Disclosure of Invention
In order to solve the problems of the simulation prediction based on the Fourier transform method, the invention provides a method for rapidly calculating the rolling resistance of a tire. The method can realize the rapid calculation of the rolling resistance of the tire and provide support for the optimization or improvement of the rolling resistance of the tire.
The invention aims at realizing the following technical scheme:
a method of rapidly calculating tire rolling resistance comprising the steps of:
step one, carrying out load simulation analysis on a tire:
step one, carrying out grid division on a design drawing (material distribution diagram) of the tire;
step one, endowing material properties to the materials of all the parts, and establishing a tire finite element model;
step three, applying rated inflation pressure to the inner surface of the tire lining layer to perform inflation analysis;
step four, applying rated load to the tire on the basis of inflation analysis, performing tire load analysis, and outputting strain energy density values of all rubber material units;
step two, reading strain energy density values of rubber material units at the same position and at different angles around the tire, and calculating strain energy density accumulation w of each rubber material unit i ：
Wherein i is the number of rubber material units, k is the total number of tire sections, n is the number of tire sections, abs () is the absolute value, e i,n A strain energy density value of the rubber material unit with the number i at the nth section;
step three, calculating the energy loss h of each rubber material unit i For the strain of the simple harmonic period and the non-simple harmonic period, the ratio of the loss energy to the elastic energy when the rubber material is deformed is (2 pi tan delta) 1/1.75 The energy loss of the rubber material unit can thus be obtained as:
wherein, tan delta i The loss tangent value of the material to which the ith rubber material unit belongs;
step four, calculating a rolling resistance value R of the tire:
wherein V is i The volume of the ith rubber material unit around the tire is r, which is the rolling radius of the tire.
Compared with the prior art, the invention has the following advantages:
(1) The rolling resistance of the tire can be obtained without complex Fourier transformation;
(2) The method is suitable for calculating rolling resistance of all tires, including tires with light, tires with longitudinal grooves and tires with complex patterns;
(3) Compared with the Fourier transform method, the rolling resistance calculation time is shortened by at least 2 times, and the precision can be kept unchanged or slightly improved.
Drawings
FIG. 1 is a tire material distribution diagram;
FIG. 2 is a grid division;
FIG. 3 is an inflation analysis;
FIG. 4 is a section division;
FIG. 5 is a load analysis;
FIG. 6 is a graph of strain energy density values for rubber material units;
FIG. 7 is a schematic view of a method and an angle for dividing a tire circumferential section;
FIG. 8 is a graph of cell strain energy density values at different angles at the same location;
FIG. 9 is the energy loss of a portion of a cell;
FIG. 10 is a graph comparing the method of the present invention with the Fourier transform method and the measured results.
Detailed Description
The following description of the present invention is provided with reference to the accompanying drawings, but is not limited to the following description, and any modifications or equivalent substitutions of the present invention should be included in the scope of the present invention without departing from the spirit and scope of the present invention.
Taking the simulation of 21550R15 tires under Abaqus software as an example:
the first step: the tire material distribution diagram (fig. 1) is subjected to grid division into triangular or quadrilateral units (fig. 2), the framework material is divided into 2-node one-dimensional units, material properties (table 1) are given to each component material, a tire finite element model is built, air pressure of 0.25MPa is applied to the inner surface of the tire inner liner, and inflation analysis is performed by using Abaqus software, and the result is shown in fig. 3. On the basis of the two-dimensional inflation analysis, the tire cross section was rotated 360 degrees, and divided into 110 cross sections in the circumferential direction (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 strain energy density values (SENER values) of all the rubber material units were output (fig. 6).
TABLE 1
And a second step of: strain energy density values (fig. 8) of units (fig. 7) of the same position and different angles around the tire were extracted by using Python program, and strain energy density accumulation w of each rubber material unit was calculated i ：
Where i is a unit number, n is a tire cross-section number, abs () is an absolute value, e i,n The strain energy density value at the nth section for the cell numbered i. Taking the example of unit 1, the tire model has 110 sections, k=110, e 1,1 ＝0.0018996，e 1,2 ＝0.0020277，e 1,3 ＝0.002156，e 1,4 ＝0.002305，e 1,5 ＝0.0024541，e 1,6 ＝0.0026104，e 1,7 ＝0.0027667，e 1,8 ＝0.0029258，e 1,9 ＝0.0030821，e 1,10 ＝0.0032433，e 1,11 = 0.0034032, etc.), then:
and a third step of: calculating the energy loss h of each unit i For the strain of the simple harmonic period and the non-simple harmonic period, the ratio of the loss energy to the elastic energy when the rubber material is deformed is (2 pi tan delta) 1/1.75 The energy loss of the unit can thus be obtained as:
wherein, tan delta i Is the loss tangent of the material to which the ith cell belongs. Taking the example of unit number 1, which belongs to the carcass part, the loss tangent tan delta of the material 1 ＝0.1，w 1 ＝0.02004052:
Fourth step: the rolling resistance value R of the tire was calculated, and the tire model used in this example had 1380 rubber material units in total, and the energy loss of a part of the units was shown in table 2:
wherein V is i The volume of the ith unit around the tire is extracted from a model calculation file, and r is the rolling radius of the tire and r is the inflation radius of the tire model 0 And the load deflection d, in this embodiment the inflation radius r 0 R= 287.44mm was calculated as 315.877mm and d= 28.437mm for load sinking. The rolling resistance of the tire can thus be calculated as:
R＝(h 1 V 1 +h 2 V 2 +h 3 V 3 +h 4 V 4 +h 5 V 5 +h 6 V 6 +…)/(2π×287.44)＝27.918N
TABLE 2
| Unit numbering | Loss of unit energy |
| 1 | 0.009405785 |
| 3 | 0.013866373 |
| 4 | 0.011053977 |
| 6 | 0.013533912 |
| 7 | 0.007356863 |
| 9 | 0.010121829 |
| 10 | 0.010825596 |
| 12 | 0.006243324 |
| 13 | 0.01072208 |
| 15 | 0.008950428 |
| 16 | 0.013483618 |
| 18 | 0.008960802 |
| 19 | 0.014359517 |
| 20 | 0.00584299 |
| 22 | 0.017246872 |
| 25 | 0.012827823 |
| 26 | 0.010716181 |
| 27 | 0.013931255 |
| 28 | 0.009329001 |
| 29 | 0.013060762 |
| 32 | 0.009725073 |
| 33 | 0.00923388 |
| 34 | 0.00844858 |
The tire rolling resistance is calculated rapidly through the method, the whole calculation process only needs 2 seconds, and the rolling resistance calculation of the unified tire model based on the Fourier transform method needs 3 minutes, so that the efficiency is improved by 90 times. In order to verify the validity of the calculation result, the method proposed by the invention is fully verified, as shown in fig. 10. As can be seen from fig. 10, the rolling resistance of the tire predicted by the method of the present invention is consistent with the trend of the measured result, and the error is smaller in value than that of the fourier transform method, which proves the effectiveness of the method of the present invention.
Claims (1)
1. A method for rapidly calculating the rolling resistance of a tyre, characterized in that it comprises the steps of:
step one, carrying out load simulation analysis on a tire:
step one, carrying out grid division on a material distribution diagram of the tire;
step one, endowing material properties to the materials of all the parts, and establishing a tire finite element model;
step three, applying rated inflation pressure to the inner surface of the tire lining layer to perform inflation analysis;
step four, applying rated load to the tire on the basis of inflation analysis, performing tire load analysis, and outputting strain energy density values of all rubber material units;
step two, reading strain energy density values of rubber material units at the same position and at different angles around the tire, and calculating strain energy density accumulation w of each rubber material unit i The strain energy density of each rubber material unit is accumulated w i The calculation formula of (2) is as follows:
wherein i is the number of rubber material units, k is the total number of tire sections, n is the number of tire sections, abs () is the absolute value, e i,n A strain energy density value of the rubber material unit with the number i at the nth section;
step three, calculating the energy loss h of each rubber material unit i ：
Wherein, tan delta i The loss tangent value of the material to which the ith rubber material unit belongs;
step four, calculating a rolling resistance value R of the tire, wherein the calculation formula of the rolling resistance value R of the tire is as follows:
wherein V is i The volume of the ith rubber material unit around the tire is r, which is the rolling radius of the tire.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310472577.1A CN116432320B (en) | 2023-04-27 | 2023-04-27 | Method for rapidly calculating rolling resistance of tire |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310472577.1A CN116432320B (en) | 2023-04-27 | 2023-04-27 | Method for rapidly calculating rolling resistance of tire |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116432320A CN116432320A (en) | 2023-07-14 |
| CN116432320B true CN116432320B (en) | 2023-11-14 |
Family
ID=87094365
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310472577.1A Active CN116432320B (en) | 2023-04-27 | 2023-04-27 | Method for rapidly calculating rolling resistance of tire |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116432320B (en) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010249527A (en) * | 2009-04-10 | 2010-11-04 | Yokohama Rubber Co Ltd:The | Method for evaluating rolling resistance of tire, system for evaluating tire using the same, and program for evaluating rolling resistance of tire |
| CN114330059A (en) * | 2021-12-27 | 2022-04-12 | 中策橡胶集团股份有限公司 | Method, apparatus and computer program for calculating the contribution of tire tread compression to rolling resistance |
| CN114692297A (en) * | 2022-03-14 | 2022-07-01 | 中策橡胶集团股份有限公司 | Method for analyzing influence degree of rubber matrix on tire rolling resistance, application and program product |
| CN114707374A (en) * | 2022-03-14 | 2022-07-05 | 中策橡胶集团股份有限公司 | Method, application and program product for calculating influence of rubber material modulus on tire rolling resistance |
| CN115659625A (en) * | 2022-10-20 | 2023-01-31 | 中策橡胶集团股份有限公司 | Method for predicting rolling resistance of periodic pattern tire, application and computer program product |
-
2023
- 2023-04-27 CN CN202310472577.1A patent/CN116432320B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2010249527A (en) * | 2009-04-10 | 2010-11-04 | Yokohama Rubber Co Ltd:The | Method for evaluating rolling resistance of tire, system for evaluating tire using the same, and program for evaluating rolling resistance of tire |
| CN114330059A (en) * | 2021-12-27 | 2022-04-12 | 中策橡胶集团股份有限公司 | Method, apparatus and computer program for calculating the contribution of tire tread compression to rolling resistance |
| CN114692297A (en) * | 2022-03-14 | 2022-07-01 | 中策橡胶集团股份有限公司 | Method for analyzing influence degree of rubber matrix on tire rolling resistance, application and program product |
| CN114707374A (en) * | 2022-03-14 | 2022-07-05 | 中策橡胶集团股份有限公司 | Method, application and program product for calculating influence of rubber material modulus on tire rolling resistance |
| CN115659625A (en) * | 2022-10-20 | 2023-01-31 | 中策橡胶集团股份有限公司 | Method for predicting rolling resistance of periodic pattern tire, application and computer program product |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116432320A (en) | 2023-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103245516B (en) | For the system characterizing the machine of tyre evenness and the method using this sign | |
| CN109426685B (en) | Component fatigue analysis method | |
| CN116432320B (en) | Method for rapidly calculating rolling resistance of tire | |
| Cho et al. | Finite element estimation of hysteretic loss and rolling resistance of 3-D patterned tire | |
| CN114707374A (en) | Method, application and program product for calculating influence of rubber material modulus on tire rolling resistance | |
| CN111191397B (en) | Rapid prediction method for static radial stiffness of radial tire | |
| Yang et al. | Materials testing for finite element tire model | |
| KR101758048B1 (en) | Tire uniformity improvement using estimates based on convolution/deconvolution with measured lateral force variation | |
| CN115659625A (en) | Method for predicting rolling resistance of periodic pattern tire, application and computer program product | |
| CN116451349A (en) | Prediction method for vertical rigidity of spoke type non-pneumatic tire and application thereof | |
| WO2014051613A1 (en) | Improvement of the uniformity of a tire using tooling signature analysis | |
| Golbakhshi et al. | Evaluating the effect of dissipated viscous energy of a rolling tire on stress, strain and deformation fields using an efficient 2D FE analysis | |
| CN112734203B (en) | Vehicle damage calculation method, device and system based on road surface and storage medium | |
| CN115391916A (en) | Wheel double-shaft fatigue simulation analysis method, device, equipment and medium | |
| CN114235448A (en) | Rail vehicle bogie wheel fatigue damage assessment method and system | |
| JP4581671B2 (en) | Method and apparatus for predicting cross-sectional shape during tire manufacture, computer program and information storage medium therefor | |
| EP3065941B1 (en) | Tire uniformity improvement through modified sampling of uniformity parameters | |
| JP2007131209A (en) | Model for numerical analysis of tire, its forming method, method for analyzing rolling resistance of tire, and model for numerical analysis of rolling resistance of tire | |
| Heinrietz et al. | Identification of parametric tire models for the fatigue evaluation of suspension components | |
| Vinchu et al. | Weight Optimization of Vehicle Rim with Fatigue Tests | |
| Wang et al. | Nummerical Simulation-Based Analysis of the Effect of Curing Bladders on Shaping of Tires | |
| EP3065942A1 (en) | Improvement of the uniformity of a tire using estimation of transient effects | |
| CN114201902A (en) | Method, device and program for predicting contribution of uniform distribution of ground contact pressure to tire rolling resistance | |
| Jackowski et al. | Comparative analysis of small size non-pneumatic tires and pneumatic tires-radial stiffness and hysteresis, selected parameters of the contact patch | |
| EP2979074B1 (en) | Tire uniformity estimation using estimates based on convolution/deconvolution |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |