CN112287465A - Method for realizing real-time simulation of in-plane tire flexible ring model - Google Patents
Method for realizing real-time simulation of in-plane tire flexible ring model Download PDFInfo
- Publication number
- CN112287465A CN112287465A CN202011319509.4A CN202011319509A CN112287465A CN 112287465 A CN112287465 A CN 112287465A CN 202011319509 A CN202011319509 A CN 202011319509A CN 112287465 A CN112287465 A CN 112287465A
- Authority
- CN
- China
- Prior art keywords
- model
- tire
- bending
- flexible ring
- spring
- 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.)
- Granted
Links
Images
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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Tires In General (AREA)
Abstract
A method for realizing real-time simulation of an in-plane tire flexible ring model comprises the steps of establishing the tire flexible ring model, simplifying a tire belt part into mass points connected by a series of spring units, and simplifying compressed air in a tire side part and a tire into a spring damping unit; a linear expression mode is adopted during modeling of the bending characteristic of the beam-containing part, so that the non-linear degree of the model is reduced, and the calculation efficiency of the model is improved; in the model solving process, a Newton iteration method with second-order convergence speed and a Newmark method with unconditional convergence characteristics are combined, so that the model solving speed can meet the requirement of real-time simulation.
Description
Technical Field
The invention relates to a method for realizing real-time simulation of an in-plane tire flexible ring model, and belongs to the technical field of tire flexible ring model modeling.
Background
The tire is the only part of the vehicle contacting with the road surface, transmits the acting force of the road surface to the vehicle, and the performance index of the whole vehicle is closely related to the mechanical property of the tire, so that the vehicle dynamics simulation analysis needs an accurate and efficient tire model.
At present, tire models are classified into the following categories according to different applicable frequency ranges: the method is suitable for a tire steady-state model with the tire motion input frequency lower than 8 Hz; the method is suitable for a tire rigid ring model with the tire motion input frequency lower than 80 Hz; the method is suitable for the tire flexible ring model with the tire motion input frequency higher than 80 Hz; the method is theoretically suitable for the tire finite element model with the tire motion input frequency of any value. The tire steady-state model and the rigid ring model are limited in tire motion input frequency range, the tire finite element model is huge in calculation amount and cannot be used in the whole vehicle dynamics simulation, and therefore the establishment of the tire flexible ring model with high calculation efficiency and high calculation precision is always a research hotspot of the tire model
The modeling method of the tire flexible ring model mainly comprises the following steps: the first is to discretize the tire belt into a series of mass points or concentrated masses connected by spring damping units, the sidewall portions and the compressed air being represented by spring damping units; the second is to represent the tire belt portion by a beam unit, and the sidewalls and compressed air by a spring damping system; and thirdly, simplifying the finite element model of the tire into a layer of composite shell elements. The beam unit and the shell unit are added into the model, so that the calculation consumption is increased, and real-time simulation is difficult to perform, so that the first tire flexible ring model modeling method is more common. In order to enable the tire model to have real-time simulation capability, a reasonable model structure and a rapid solving algorithm are of great importance.
In many documents relating to tire flex-ring models, modeling relating to belt tensile mechanical characteristics and modeling relating to compressed air in the sidewalls and tires are substantially linear, but there is a large room for improvement in modeling of belt bending characteristics and in model solution algorithms. At present, inverse trigonometric functions are introduced into the modeling of the belt bending characteristics, so that the nonlinearity degree of the model is increased, and the calculation efficiency of the model is reduced; in order to avoid the computational consumption caused by modeling of the belt bending characteristics, even without modeling the belt bending characteristics, some researchers do not accurately express the ground contact characteristics of the tire, which is obviously inaccurate. In the aspect of solving the algorithm, some documents adopt ODE solvers carried by commercial software MATLAB, and some documents do not refer to the specific solving algorithm of the model.
As described above, in order to obtain a tire flexible ring model which has high calculation efficiency, high calculation accuracy, and can be simulated in real time, a belt bending characteristic modeling method having a high degree of linearization and a calculation algorithm having high calculation efficiency are required.
Disclosure of Invention
The invention provides a method for realizing real-time simulation of an in-plane tire flexible ring model, and particularly introduces a belt bending characteristic modeling method with high linearization degree in the model and a solving algorithm with high calculation efficiency. The established tire flexible ring model not only has higher precision, but also has real-time simulation capability, and can be used for the complete vehicle dynamics simulation.
A method for realizing real-time simulation of an in-plane tire flexible ring model comprises the following steps:
building a flexible ring model of the tire, simplifying a belt part of the tire into mass points connected by a series of spring units, and simplifying compressed air in a sidewall part and the tire into spring damping units. A linear expression mode is adopted during modeling of the bending characteristic of the beam part, the non-linear degree of the model is reduced, and the calculation efficiency of the model is improved. In the model solving process, a Newton iteration method with second-order convergence speed and a Newmark method with unconditional convergence characteristics are combined, so that the model solving speed can meet the requirement of real-time simulation. The method specifically comprises the following steps:
firstly, dispersing a belt beam part into a plurality of mass points which are uniformly distributed, wherein the belt beam mass points are connected through an extension spring and a bending spring; the mass point and the rim are connected through a series of spring damping units to simulate compressed air in a tire sidewall part and a tire; the tread is discretized into a series of elastic elements for calculating the contact force of the tire with the road surface.
And secondly, the acting force between the mass points, which is received due to the deformation of the bending spring, is subjected to linearization treatment, so that the acting force of the bending spring, which is received by the mass points, is only linearly related to the displacement of the mass points, an inverse trigonometric function with high nonlinear degree is not introduced any more, and the nonlinear degree of belt bending characteristic modeling is reduced.
And thirdly, combining and using a Newton iteration method with second-order convergence speed and a Newmark method with unconditional convergence characteristics in the model resolving process. And converting the model kinetic equation set into a linear equation set by a Newmark method, and then solving the linear equation set by a Newton iteration method.
The invention has the beneficial effects that:
1. the invention has the advantages of high applicable frequency range, accurate calculation and capability of carrying out real-time operation.
2. The invention solves the problem of high nonlinearity of the traditional tire flexible ring model with the restraint bending characteristic modeling, reduces the complexity of the model, and plays an important role in accurate and real-time operation of the model.
3. The invention can carry out the solving algorithm of the real-time operation of the model; the time for solving is saved, and the real-time performance is realized.
Drawings
FIG. 1 is a schematic structural diagram of an in-plane tire flexible ring model built by the present invention.
FIG. 2 is a schematic view of a tread model in the tire model established by the present invention
FIG. 3 is a graphical illustration of bending force linearization.
Detailed Description
A method for realizing real-time simulation of an in-plane tire flexible ring model comprises the following steps:
first, as shown in FIG. 1, the belt portion is discretized into a plurality of uniformly distributed mass points miMass point miConnected by springs along the circumferential direction and the radial direction of the tire, expresses the tensile property of the belt, and has three adjacent mass points miThe bending characteristic of the belt is expressed by connecting the two through a bending spring, and the mass point miThe rim accommodates compressed air in the tire and the sidewall of the tire via the spring damping units in the circumferential direction and the radial direction of the tire. As shown in fig. 2, the tread is discretized to be evenly distributed at the mass point miBy the deformation of the spring, the tyre and the road surface are calculatedThe contact force of (2).
Two, to mass point miThe mass point m is linearized by the acting force of the bending spring deformationiThe applied force of the bending spring is only linearly related to the displacement of the mass point, an inverse trigonometric function with high nonlinear degree is not introduced any more, and the nonlinear degree of belt bending characteristic modeling is reduced.
The structure shown in figure 3, two rigid rods are connected by a bending spring, and when three points ABC respectively generate displacement uk-1ukuk+1In time, according to the calculation method of the prior literature, the point C is subjected to the acting force F of a bending spring due to the change of & lt ABCk+1The calculation formula is as follows:
in the above formula: k is a radical ofbkIs the bending spring rate; thetakAnd thetak-1The included angles between the BC rod and the AB rod and the Z axis are respectively formed; lkThe length of the projection of the BC pole on the Z axis; wherein with respect to thetakThe solution of (2) involves the application of an inverse trigonometric function, increasing the degree of model non-linearity.
The potential energy of the bending spring is derived by applying the clip theorem to obtain the acting force of the bending spring on the point C due to the change of the & lt ABC. Firstly, calculating the displacement u generated when three points of ABC are respectively displacedk-1ukuk+1The calculation formula of the potential energy V of the time bending spring is as follows:
in the above formula: k is a radical ofbkIs the bending spring rate; u. ofk-1、ukAnd uk+1Displacement of points A, B and C, respectively; lk-1And lkThe lengths of the projection of the AB rod and the BC rod on the Z axis respectively
The action force of a bending spring on the C point obtained by applying the clip theorem due to the change of & lt ABC is as follows:
method for calculating bending force of the invention requires thetakSatisfies thetak≈tan(θk) In the process of model simulation, when the number of the mass points reaches a certain value, the requirement of theta is completely metk≈tan(θk)。
And thirdly, combining and using a Newton iteration method with second-order convergence speed and a Newmark method with unconditional convergence characteristics in the model resolving process. And converting the model kinetic equation set into a linear equation set by a Newmark method, and then solving the linear equation set by a Newton iteration method.
The basic assumptions of the Newmark method are as follows:
in the above formulaRespectively representing the displacement state, the speed state and the acceleration state of each mass point in the t moment model;respectively representing the displacement state, the speed state and the acceleration state of each mass point in the t + delta t moment model; the time for simulating is not long at; γ and β are parameters of the Newmark method, typically β is 0.25 and γ is 0.5;
expressing the speed state and the acceleration state in the t + delta t moment model as a function of the displacement state in the t + delta t moment model through the known states of all the mass points in the t moment model, so that a motion equation set of the t + delta t moment model is converted into an algebraic equation set, and then the algebraic equation set is solved by a Newton iteration method;
the basic format of the Newton iteration method is as follows
In the above formula xkAnd xk+1The values of the variables in the k-th iteration step and the k + 1-th iteration step equation set, f' (x)k) Jacobian matrix f (x) of the kth iteration step equation setk) Is the value of the kth set of iteration step equations.
Due to the local convergence characteristic of the Newton iteration method, the calculation efficiency of the model, the contact change condition of the model, the rolling speed of the model and the like need to be comprehensively considered when the time step is selected
By the bending force calculation method introduced by the invention, the non-linear degree of the model is reduced, so that the solving algorithm of the invention can be smoothly implemented; the currently common bending force calculation method cannot use the solving algorithm mentioned in the invention due to the highly nonlinear characteristic, and generally adopts Runge-Kutta solving algorithm. Solving 1s simulation, the same computer equipment (CPU: Intel i 59600 kf) and the same compiling software (VS2019), the currently common bending force calculation method and Runge-Kutta solving algorithm need about 13s of calculation time, the bending force calculation method and the solving algorithm introduced in the invention need about 0.65s of calculation time, the solving time is saved, and the method has real-time performance.
Claims (4)
1. A method for realizing real-time simulation of an in-plane tire flexible ring model is characterized by comprising the following steps: the method comprises the following steps:
building a flexible ring model of the tire, simplifying a belted part of the tire into mass points connected by a series of spring units, and simplifying compressed air in a sidewall part and the tire into a spring damping unit; a linear expression mode is adopted during modeling of the bending characteristic of the beam-containing part, so that the non-linear degree of the model is reduced, and the calculation efficiency of the model is improved; in the model solving process, a Newton iteration method with second-order convergence speed and a Newmark method with unconditional convergence characteristics are combined, so that the model solving speed can meet the requirement of real-time simulation.
2. The method for realizing the in-plane tire flexible ring model real-time simulation of the claim 1, is characterized in that: the method comprises the following steps:
firstly, dispersing a belt beam part into a plurality of mass points which are uniformly distributed, wherein the belt beam mass points are connected through an extension spring and a bending spring; the mass point and the rim are connected through a series of spring damping units to simulate compressed air in a tire sidewall part and a tire; discretizing the tread into a series of elastic elements for calculating the contact force of the tyre with the road surface;
secondly, the acting force between the mass points, which is received due to the deformation of the bending spring, is subjected to linearization treatment, so that the acting force of the bending spring, which is received by the mass points, is only linearly related to the displacement of the mass points, an inverse trigonometric function with high nonlinear degree is not introduced any more, and the nonlinear degree of belt bending characteristic modeling is reduced;
combining and using a Newton iteration method with second-order convergence speed and a Newmark method with unconditional convergence characteristics in the model resolving process; and converting the model kinetic equation set into a linear equation set by a Newmark method, and then solving the linear equation set by a Newton iteration method.
3. The method for realizing the in-plane tire flexible ring model real-time simulation of the claim 2 is characterized in that: obtaining the acting force of the bending spring on the point C due to the change of the & lt ABC by applying the card theorem to derive the potential energy of the bending spring; firstly, calculating the displacement u generated when three points of ABC are respectively displacedk-1ukuk+1The calculation formula of the potential energy V of the time bending spring is as follows:
in the above formula: k is a radical ofbkIs the bending spring rate; u. ofk-1、ukAnd uk+1Displacement of points A, B and C, respectively; lk-1And lkThe lengths of the projection of the AB rod and the BC rod on the Z axis respectively
The action force of a bending spring on the C point obtained by applying the clip theorem due to the change of & lt ABC is as follows:
in the above formula: k is a radical ofbkIs the bending spring rate; u. ofk-1、ukAnd uk+1Displacement of points A, B and C, respectively; lk-1And lkThe lengths of the projection of the AB rod and the BC rod on the Z axis respectively
Method of calculating bending force requires θkSatisfies thetak≈tan(θk) In the process of model simulation, when the number of the mass points reaches a certain value, the requirement of theta is completely metk≈tan(θk)。
4. The method for realizing the in-plane tire flexible ring model real-time simulation of the claim 2 is characterized in that: the basic assumptions of the Newmark method are as follows:
in the above formula, U (t)Respectively representing the displacement state, the speed state and the acceleration state of each mass point in the t moment model; u (t + delta t)Displacement of each mass point in the model at time t + Δ tA state, a velocity state, and an acceleration state; the time for simulating is not long at; γ and β are parameters of the Newmark method, typically β is 0.25 and γ is 0.5;
expressing the speed state and the acceleration state in the t + delta t moment model as a function of the displacement state in the t + delta t moment model through the known states of all the mass points in the t moment model, so that a motion equation set of the t + delta t moment model is converted into an algebraic equation set, and then the algebraic equation set is solved by a Newton iteration method;
the basic format of the Newton iteration method is as follows
In the above formula xkAnd xk+1The values of the variables in the k-th iteration step and the k + 1-th iteration step equation set, f' (x)k) Jacobian matrix f (x) of the kth iteration step equation setk) Is the value of the kth set of iteration step equations.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011319509.4A CN112287465B (en) | 2020-11-23 | 2020-11-23 | Method for realizing real-time simulation of in-plane tire flexible ring model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011319509.4A CN112287465B (en) | 2020-11-23 | 2020-11-23 | Method for realizing real-time simulation of in-plane tire flexible ring model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112287465A true CN112287465A (en) | 2021-01-29 |
CN112287465B CN112287465B (en) | 2022-03-18 |
Family
ID=74425183
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011319509.4A Active CN112287465B (en) | 2020-11-23 | 2020-11-23 | Method for realizing real-time simulation of in-plane tire flexible ring model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112287465B (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113239450A (en) * | 2021-03-11 | 2021-08-10 | 华南理工大学 | Interconnected air suspension multi-objective optimization method based on interval analysis |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002092364A2 (en) * | 2001-05-15 | 2002-11-21 | Wilson Kitchener C | Vehicle and vehicle tire monitoring system, apparatus and method |
CN101001762A (en) * | 2004-07-16 | 2007-07-18 | 倍耐力轮胎股份公司 | High-performance tyre for a motor vehicle |
WO2012074527A1 (en) * | 2010-12-02 | 2012-06-07 | Michelin Recherche Et Technique, S.A. | Method for prediction and control of tire uniformity parameters from crown thickness variation |
EP2604990A1 (en) * | 2011-12-13 | 2013-06-19 | Snap-on Equipment Srl a unico socio | Method and system for determining and improving running characteristics of a pneumatic tyre of a vehicle wheel |
CN106515348A (en) * | 2016-12-23 | 2017-03-22 | 长春孔辉汽车科技股份有限公司 | Intelligent accelerated speed damping semi-active control method for vehicle suspension system |
CN111767653A (en) * | 2020-07-01 | 2020-10-13 | 中国人民解放军火箭军工程大学 | Heavy-duty tire modeling method based on analytic sidewall stiffness |
-
2020
- 2020-11-23 CN CN202011319509.4A patent/CN112287465B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002092364A2 (en) * | 2001-05-15 | 2002-11-21 | Wilson Kitchener C | Vehicle and vehicle tire monitoring system, apparatus and method |
CN101001762A (en) * | 2004-07-16 | 2007-07-18 | 倍耐力轮胎股份公司 | High-performance tyre for a motor vehicle |
WO2012074527A1 (en) * | 2010-12-02 | 2012-06-07 | Michelin Recherche Et Technique, S.A. | Method for prediction and control of tire uniformity parameters from crown thickness variation |
EP2604990A1 (en) * | 2011-12-13 | 2013-06-19 | Snap-on Equipment Srl a unico socio | Method and system for determining and improving running characteristics of a pneumatic tyre of a vehicle wheel |
CN106515348A (en) * | 2016-12-23 | 2017-03-22 | 长春孔辉汽车科技股份有限公司 | Intelligent accelerated speed damping semi-active control method for vehicle suspension system |
CN111767653A (en) * | 2020-07-01 | 2020-10-13 | 中国人民解放军火箭军工程大学 | Heavy-duty tire modeling method based on analytic sidewall stiffness |
Non-Patent Citations (9)
Title |
---|
DANG LU 等: "Tire side slip rolling Simulation with flexible carcass", 《2011 3RD INTERNATIONAL CONFERENCE ON ADVANCED COMPUTER CONTROL》 * |
YI XIONG 等: "Case Studies in Mechanical Systems and Signal Processing", 《JOURNALS & BOOKS》 * |
ZHIMING LIU: "Tire Longitudinal Mechanics Properties under Various BraKing Conditions", 《APPLIED MECHANICS AND MATERIALS》 * |
刘志浩 等: "基于弹性基础柔性梁的重载轮胎面内胎体与胎侧耦合建模及参数辨识", 《振动与冲击》 * |
刘志浩 等: "重载轮胎面内欧拉梁模型及耦合振动特性研究", 《振动工程学报》 * |
史石荣: "重载车辆条件下的车路耦合动力分析模型和车路相互作用分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
郭孔辉: "适用于高频激励的面内轮胎动态模型", 《吉林大学学报(工学版)》 * |
黄世庆: "适用于短波长不平路面的轮胎柔性环模型研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
黄新艺: "混凝土连续曲线梁桥在车辆荷载作用下的动力响应研究", 《中国优秀博士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113239450A (en) * | 2021-03-11 | 2021-08-10 | 华南理工大学 | Interconnected air suspension multi-objective optimization method based on interval analysis |
CN113239450B (en) * | 2021-03-11 | 2022-12-16 | 华南理工大学 | Interconnected air suspension multi-objective optimization method based on interval analysis |
Also Published As
Publication number | Publication date |
---|---|
CN112287465B (en) | 2022-03-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Nariman-Zadeh et al. | Evolutionary design of generalized polynomial neural networks for modelling and prediction of explosive forming process | |
Oden et al. | Numerical analysis of certain contact problems in elasticity with non-classical friction laws | |
CN108446445B (en) | Composite material wing optimization design method based on aerodynamic reduced order model | |
CN107844627A (en) | It is a kind of only to export Time variable structure modal parameter Bayesian Estimation method | |
CN106096257A (en) | A kind of non-linear cable elements analyzes method and system | |
CN110162826B (en) | Analysis method for thermoelasticity dynamic response of thin-wall structure | |
Dang et al. | Accelerated loosely-coupled CFD/CSD method for nonlinear static aeroelasticity analysis | |
CN107766670B (en) | Periodical chiral cellular structural material equivalent elastic modulus prediction technique | |
CN112287465B (en) | Method for realizing real-time simulation of in-plane tire flexible ring model | |
CN103154703B (en) | Method for predicting elastic response performance of rubber product, method for design, and device for predicting elastic response performance | |
CN113723027A (en) | Static aeroelasticity calculation method for elastic aircraft | |
CN104281730A (en) | Great-rotating-deformation plate shell structure dynamic response finite element analysis method | |
CN111125963A (en) | Numerical simulation system and method based on Lagrange integral point finite element | |
CN110717216A (en) | Method for forecasting rolling response of helicopter with flexible air bag under irregular wave | |
CN114997027A (en) | Method for intelligently solving random signals of axle system | |
CN109271655B (en) | Material scale effect analysis method based on asymmetric finite element algorithm | |
CN106055733B (en) | The kinetic parameter of multifunction structure determines method | |
JP2009045627A (en) | Method and program for optimizing press forming condition | |
CN111177960B (en) | Thin shell collision contact calculation method based on ANCF | |
CN117131703A (en) | Airplane sliding power response acquisition method based on airplane-runway interaction | |
CN103154702B (en) | Method for predicting elastic response performance of rubber product, method for design, and device for predicting elastic response performance | |
CN111563297A (en) | Supersonic aircraft thermal environment calculation method based on BP network | |
CN105117519A (en) | Electrostatically driven stepped micro-cantilever structure evaluation method and system | |
CN115455754A (en) | Mine hydraulic support design method based on digital twinning | |
CN114297890A (en) | Tire modal simulation method, device and program based on rubber material equivalent density calculation |
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 |