Method for the simulation of the physical behavior of a tire rolling on the ground

• US 8,560,289 B2
• Filed: 07/26/2007
• Issued: 10/15/2013
• Est. Priority Date: 09/01/2006
• Status: Active Grant

• Current Assignee: Michelin Recherche et Technique SA (Compagnie Generale DES Etablissements Michelin SCA)
• Original Assignee: Michelin Recherche et Technique SA (Compagnie Generale DES Etablissements Michelin SCA)
• Inventors: Fevrier, Pierre, Martin, Herv
• Primary Examiner(s): Rivas, Omar Fernandez
• Assistant Examiner(s): LUU, CUONG V
• Application Number:
  US12/439,683
• Publication Number:
  US 20090292515A1
• Time in Patent Office: 2,273 Days
• Field of Search: 702/42, 702/34, 703/2, 703/6, 703/8, 731/46, 73/8, 701/1, 701/71
• US Class Current: 703/8
• CPC Class Codes: B60C 99/006 Computer aided tyre design or simulation
• CPC Class Codes: G01M 17/02 Tyres
• CPC Class Codes: G06F 2119/08 Thermal analysis or thermal optimisation
• CPC Class Codes: G06F 30/15 Vehicle, aircraft or watercraft design
• CPC Class Codes: G06F 30/20 Design optimisation, verification or simulation optimisation, verification or simulation of circuit designs G06F30/30
• CPC Class Codes: G06F 30/23 using finite element methods [FEM] or finite difference methods [FDM]

Abstract

A method for simulating the behavior of a tire mounted on a vehicle in running conditions on the ground, wherein a mechanical model is provided for essentially computing the longitudinal (Fx) and transverse (Fy) stresses transmitted by the tire between the ground and the vehicle in accordance with dynamic parameters related to the physical conditions of the tire running and use and in accordance with physical tire-specific parameters. The mechanical model is set and solved in an iterative manner, under the assumption that the tire in contact with the surface of the ground has an adherence contact area and a sliding contact area and under the assumption that there is a unique x-coordinate point b that is indicative of the transition between the two contact areas.

Claims

1. A method of simulating the physical behavior of a tire equipping a vehicle rolling on the ground, wherein the tread of the tire has a contact area which includes an adherent contact area and a sliding contact area, wherein the method comprises:

a computer exercising the steps of;

calculating longitudinal forces and lateral forces transmitted by the tire, between the ground and the vehicle, based on dynamic parameters related to physical rolling comprising a slip angle and a longitudinal slip rate of the tire and operating conditions of the tire and based on specific physical parameters of the tire;

establishing, over the course of a preliminary phase, by applying physical laws which are at least one of known and established by specific experiments, and as a first model, a model of the longitudinal forces, the lateral forces, a self-aligning torque related to the intensity of these forces and to the distribution of same within the contact area, and of an equilibrium of the basic shear and sliding forces of the tire at a presumed unique point of passage between the adherent and sliding contact regions, this first model forming a system of respective equations expressed in relation to the dynamic parameters, the specific parameters, and the abscissa of the point of passage;

assigning values to the dynamic parameters and to the specific parameters, with a view to a digital application;

solving, by successive approximations, over the course of a digital application iterative phase including a succession of computation cycles, and at least from previously known or estimated values of the abscissa of the point of passage, of the lateral forces, and of the self-aligning torque of new values for the abscissa of the point of passage, of the lateral forces, and of the self-aligning torque, which solve the system of equations of the first model for the values assigned to the dynamic parameters and to the specific parameters, the result of which is that the longitudinal forces, the lateral forces and the self-aligning torque can be computed in real time as the vehicle is traveling;

taking account, in the first model established over the course of the preliminary phase, of an influence of temperature on the values of at least some of the specific parameters;

establishing, over the course of this preliminary phase, by applying physical laws which are at least one of known and established by specific experiments, and as a second model, a local temperature rise model expressing variations in a contact temperature of the tire tread with the ground from a beginning to an end of the contact area, via contact of the tire tread with the ground and sliding of the tire over the ground;

establishing, over the course of this preliminary phase, by applying physical laws which are at least one of known and established by specific experiments, and as a third model, a global temperature rise and thermal flux model, this third model forming a system of equations expressing variations in a peripheral tire tread temperature and an internal tire temperature, over one wheel revolution period, in relation to previously known or estimated values of the peripheral and internal temperatures, thermal conductivity of the tire tread, and thermodynamic component phenomena including one or more of internal deformations of the tire, heat exchanges between the tire and an environment thereof, and the sliding of the tire tread over the ground;

accepting, as additional dynamic parameters, a temperature of the ground, air temperature of the tire, and an initial temperature of the tire;

computing, at each computing cycle for each iterative phase, new values for the abscissa of the point of passage, for the lateral forces and for the self-aligning torque, by using the first model enhanced by the influence of temperature; and

upon interruption of each iterative phase, updating the values for the peripheral and internal temperatures to take into account changes undergone by these temperatures since an end of a preceding iterative phase;

wherein the longitudinal forces are obtained as a function of at least the slip angle and the longitudinal slip rate.

2. The simulation method of claim 1, further comprising an operation implemented after the end of the iterative phase and including updating at least the dynamic parameters, in order to take account of the changes undergone by these parameters during execution of the iterative phase, and to initiate a new iterative phase.

3. The simulation method of claim 2, wherein each new computing cycle for each iterative phase includes the steps of

computing a new temporary value for the abscissa of the point of passage from the equilibrium of basic forces equations and from the previously known or estimated values for the lateral forces and the self-aligning torque;

computing, from the new temporary value of the abscissa of the point of passage and from the equations linking the lateral forces and the self-aligning torque to the dynamic parameters, to the specific parameters, and to the abscissa of the point of passage, new values for the lateral forces and the self-aligning torque, which can be used for a possible subsequent computing cycle;

conditionally interrupting said iterative phase, at least when the deviation between the new temporary value of the abscissa of the point of passage and the previously known or estimated value of this abscissa is lower than the predetermined accuracy limit; and

upon interruption of said iterative phase, and as values for this phase, assigning, to the lateral forces and self-aligning torque, the new values for the lateral forces and the self-aligning torque, which were obtained during the last computing cycle.

4. The simulation method of claim 1, wherein the first model is enhanced by at least taking account of the influence of change in the contact temperature on the value of the friction coefficient of the tire rubber on the ground.

5. The simulation method of claim 4, wherein the first model is enhanced by taking account of the influence of the internal tire temperature on the value of the shear modulus of the rubber compound of the tire tread.

6. The simulation method of claim 1, wherein the second computing model uses the peripheral temperature of the tire to compute the contact temperature at the entry to the contact surface.

7. The simulation method of claim 1, wherein each iterative phase is preceded by a preparatory phase, over the course of which quota quantities are computed, including the dimensions of the contact area, by means of a sub-set of said equations and based on values assigned to the dynamic parameters and to the specific parameters, each quota quantity being used in said iterative phase with a value assigned thereto during the preparatory phase.

8. Application of the method of claim 1 to the real-time simulation of the dynamic behavior of a vehicle equipped with a chassis and with several tires rolling on the ground, an application in which at least the first model is used for each tire and is associated with a single dynamic chassis model, in which, for each tire, the chassis model provides at least said first model with the values of at least some of the dynamic parameters, and in which, for each tire, the chassis model uses the values for the longitudinal forces, the lateral forces and the self-aligning torque, which were obtained by implementing at least the first model.