CN105512424A - Method for obtaining engineering machinery tyre vertical characteristic parameters based on pulse test - Google Patents
Method for obtaining engineering machinery tyre vertical characteristic parameters based on pulse test Download PDFInfo
- Publication number
- CN105512424A CN105512424A CN201510969939.3A CN201510969939A CN105512424A CN 105512424 A CN105512424 A CN 105512424A CN 201510969939 A CN201510969939 A CN 201510969939A CN 105512424 A CN105512424 A CN 105512424A
- Authority
- CN
- China
- Prior art keywords
- tire
- vertical
- finite element
- pulse
- element model
- 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
Classifications
-
- 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]
-
- 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)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Tires In General (AREA)
Abstract
The invention provides a method for obtaining engineering machinery tyre vertical characteristic parameters based on a pulse test. The method includes the steps that by building a tyre-road surface coupling three-dimensional finite element model and pulse working condition virtual dynamic simulation, the tyre center point vertical acceleration corresponding to multiple sets of estimation material parameters is obtained; secondly, data are obtained through a real vehicle road pulse test, and the material parameters of all components of a tyre are distinguished out through a neural network; finally, a vertical rigidity damping parameter of an oversized tyre is obtained through a vertical static load test and a free vibration attenuation test based on a virtual test rig. According to the method, a simulation method is combined with the real vehicle test, a tyre acceleration feature curve is obtained, accordingly the tyre vertical rigidity damping parameter is obtained, restraints of the tyre size and load are avoided, a bench test is not needed, and test cost is low; the obtained tyre material parameters are more close to reality compared with component material parameters measured by a separated material performance test; the method is suitable for obtaining mechanical feature parameters of the oversized tyre.
Description
Technical field
The present invention relates to a kind of method obtaining off-the-road tyre Vertical Characteristic parameter based on pulse testing, refer to a kind of method obtaining off-the-road tyre vertical stiffness damping parameter based on pulse testing especially; Belong to field of automobile design.
Background technology
The sole support that tire contacts with road surface as vehicle, not only carry tare, and as the first order vibration insulating system of vehicle, buffering and decay are in the process of moving from the shock and vibration on road surface, and its mechanical characteristic is the basis of Vehicle Dynamics Analysis and research.Over size tyre is large because of its bearing capacity, can effective Lifting Convey efficiency, is widely used in the large engineering vehicles such as large-tonnage quarry tipper, wheel loader.Its structural constituent is similar to tyres for passenger cars, is compositely suppressed to form by multiple material, mainly comprises casingply, belt, wire loop and rubber layer.Its Main Mechanical has: determine to a great extent and affect the driveability of vehicle and the vertical sliding characteristic of tire of braking safety; Determine the handling maneuver of vehicle and the tire cornering characteristics of directional stability; NVH (noise, vibration, the ride comfort) characteristic of vehicle is had to the tire vertical stiffness damping characteristic of material impact; Affect the tire drag characteristic etc. of VE Vehicle Economy.Thus, build reasonable and accurate tire mechanics model, for raising Full Vehicle Dynamics simulation accuracy, carry out the dynamics of vehicle exploitation relevant to tire, significant.Domestic and international experts and scholars have carried out a large amount of research work in this respect, construct various informative tire theoretical model or experience semiempirical model, and these models can describe many-sided mechanical characteristic of tire to a certain extent simultaneously and be widely used.Theoretical model is on the basis that physical tire model carries out Rational Simplification, builds mathematical model to be described mechanics of tire characteristic.Representative has: Fiala tire model, Pacejka string model and Ftire tire model etc.Theoretical model highly versatile, solving precision is high, accurately can reflect that each parameter of tire model is on the impact of mechanics of tire performance, but shortcoming to be model expression loaded down with trivial details, on the low side for counting yield during vehicle dynamics simulation.Empirical model is on the basis of a large amount of tire mechanical property testing data, builds the dynamics that Related Mathematical Models expresses tire, as magic formula (Magicformula) tire model, MF-Tyre model and Swift tire model etc.The advantage of empirical model is that structure is comparatively simple, and computing is very fast, but shortcoming is larger to the dependence of test findings.Semiempirical model is then a kind of model between theoretical model and empirical model, as unification (UniTire) tire model.In addition, also have and directly obtain tire vertical stiffness damping parameter by bench test (being generally German SCHENCK company to manufacture).But the common problem being described tire dynamics characteristic by these models is, the physical significance of model parameter is clear and definite not, can not directly follow the design feature of tire itself and material parameter to connect; They all depend on a large amount of test figures targetedly simultaneously, and these testing equipments generally can only carry out dynamic test for the passenger vehicle tyre of diameter of tyres within 1.5 meters, and the test load that testing equipment can provide is no more than 1 ton.
Over size tyre build is huge, and single tire bearing capacity is close to 100 tons, and maximum gauge, close to 4 meters, is generally inner tube free tire.At present except the external giant tire companies of a few family such as French Michelin and Japanese BRIDGESTONE tyre company can for except relevant special vehicle provides the testing equipment of over size tyre, China especially carries out the mechanical property testing of over size tyre for test condition almost not relevant car load producer, and tire production producer is unwilling the test figure of being openly correlated with.
Although along with the rise of Virtual Prototype Technique, CAE technology is introduced into the research of mechanics of tire characteristic, adopts finite element method to carry out to tire the study hotspot that mechanical characteristic analysis has become tire art in recent years.But these analyze the structural strength analysis also mainly concentrating on small light passenger tire itself, seldom have and relate to the research that the tire stiffness damping characteristic be concerned about in Automobile Design process obtains problem.Due to complicacy and the singularity of tire construction, material character, load and distortion, the mechanics property analysis of over size tyre is very difficult work always, existing analytical approach seems helpless, does not almost have pertinent literature to carry out open report to the theory needed for its mechanical characteristic analysis or test method.The domestic large engineering vehicle being representative with large-tonnage electric power wheel self-discharging vehicle is in the stage of being domesticized by the external monopolization scale that enters at present, in product design process, as the parts with ground mere contact, the mechanical characteristic of over size tyre directly affects the dynamics of car load, is basis and the key of automobile research design.But reality is that developer often has to do a large amount of simplification to it in the design process because the mechanics parameters that do not have tire production producer to provide, or directly avoids actual tire model in analytical calculation, only by virtue of experience value.Result is that car load producer is difficult to confirm the computational accuracy in development and Design process, even causes very large error, seriously constrains the fiduciary level of product development, and the process that similar large engineering vehicle domesticizes.Therefore, for the specific design object such as such as vehicle NVH performance, seek over size tyre and to be correlated with the new method that mechanics parameters obtains, Promoting Enterprise grasps the core technology in product development process, promotes localization process, just seems particularly important.
Summary of the invention
For solving the problem, the present invention proposes a kind of method obtaining off-the-road tyre Vertical Characteristic parameter based on pulse testing.
The method obtaining off-the-road tyre Vertical Characteristic parameter based on pulse testing that the present invention proposes; First set up tire-road coupling three-dimensional finite element model, and carry out the virtual dynamics simulation of pulse operation, tire centerline point vertical acceleration corresponding under obtaining many group estimation material parameters; Secondly, set up the BP neural network parameter discrimination method based on single hidden layer, and utilize l-G simulation test data to train it; Then, bring by tire centerline point vertical acceleration real vehicle road pulse testing obtained the neural network trained into, pick out the material parameter of each component of tire; Finally, by based on the vertical static loading test of virtual test platform and the vertical stiffness damping parameter of free vibration die-away test acquisition over size tyre.
The present invention is based on the method that pulse testing obtains off-the-road tyre Vertical Characteristic parameter, comprise the steps:
Step one, set up tire-road coupling three-dimensional finite element model
First, according to the tire size Plotting data tire CAD cross section profile figure of design, by forming tire material different in kind, be wire loop, casing ply, belt, rubber layer four parts by tire CAD cross section profile diagram root, in described four parts, be the length of side by 1-5mm, in each part, divide quadrilateral mesh respectively; By tire CAD cross section profile figure around tire centerline axle rotating 360 degrees, obtain the tire three-dimensional finite element model by material character subregion; Contact is set up by adhesion mode between heterogeneity material area;
Secondly, connect leverage with real vehicle and boundary condition is constrained to tire, tire three-dimensional finite element model is retrained, obtains the tire three-dimensional finite element model connecting leverage constraint;
Finally, run the constraint of time channel in the face of tire as boundary condition using real vehicle, the tire three-dimensional finite element model connecting leverage constraint is retrained, obtain tire-road coupling three-dimensional finite element model;
Step 2, pulse operation Virtual Experiment Simulation
Triangular pulse block is set in the tire-road coupling three-dimensional finite element model that step one obtains, obtains the tire-road coupling three-dimensional finite element model comprising pulse packet; Triangular pulse tile height size h is:
D/15≤h≤d/12, d is diameter of tyres; Width >=the tyre width of triangular pulse block.
According to the wire loop, casing ply, belt, the Poisson ratio of rubber layer, the empirical data of elastic modulus that form tire, determine the Poisson ratio of wire loop, casing ply, belt, rubber layer, the base value of elastic modulus respectively;
The Poisson ratio of wire loop, casing ply, belt, rubber layer, the base value of elastic modulus ± 10-20% within the scope of, optimization Latin Hypercube Sampling method is adopted to choose at least 40 group wire loops, casing ply, belt, the Poisson ratio of rubber layer, the numerical value of elastic modulus, as the input parameter of pulse operation Virtual Experiment Simulation; The input parameter of every set of pulses operating mode Virtual Experiment Simulation is substituted in the tire-road coupling three-dimensional finite element model comprising triangular pulse block, carry out pulse operation Virtual Experiment Simulation with travelling speed V, obtain at least 40 tires central point vertical acceleration curves; Travelling speed V span is 20-60Km/h;
Step 3, neural metwork training
The type selecting neural network is BP network, and comprise an input layer, input layer number is n, a hidden layer, and node in hidden layer is 10, output layer, and output layer nodes is 8; N node of input layer represents n accekeration respectively, and 8 nodes of output layer represent Poisson ratio and the elastic modulus of wire loop, casing ply, belt, rubber layer in tire finite element model respectively; Input layer is tanh (tansig) function to the transport function of hidden layer, and hidden layer is non-linear (sigmoid) transport function to the transport function of output layer;
At least 35 curves in the tire centerline point vertical acceleration curve selecting step 2 to obtain, n decile is carried out to the accelerating curve section of the impulse response section of each selected curve, obtain tire centerline point n vertical acceleration value at least 35 group pulse response sections, to often organize the input amendment of n vertical acceleration value as BP network, after the BP network operation, obtain network output valve;
Repetition training is carried out to BP network;
During error≤0.5-1% between the Poisson ratio of the wire loop of the pulse Virtual Experiment Simulation tire finite element model corresponding with corresponding vertical acceleration curve when the output valve of BP network, casing ply, belt, rubber layer and the input value of elastic modulus, complete neural metwork training; Obtain the neural network trained;
Step 4, acquisition tire material parameter
By travelling speed V, pulse testing is carried out to actual vehicle, obtains tire centerline vertical acceleration curve; Then, the accelerating curve section of impulse response section in the tire centerline vertical acceleration curve of acquisition is carried out n decile, obtain tire centerline point n vertical acceleration value in impulse response section, using the input parameter of n vertical acceleration value as BP network, in the BP neural network that substitution step 3 trains, the output valve of network is wire loop in actual tire, casing ply, belt, the Poisson ratio of rubber layer, elastic mould value;
Step 5, set up virtual simulation test platform
First, wire loop in the actual tire that step 4 is obtained, casing ply, belt, the Poisson ratio of rubber layer, elastic mould value substitutes into the tire three-dimensional finite element model by material character subregion that step one obtains, obtain the tire three-dimensional finite element model imparting tire authentic material parameter, then, in the tire three-dimensional finite element model imparting tire authentic material parameter, fixed constraint is applied to tire centerline, simultaneously, set up the horizontal supporting plate with the vertical lower surface rigid contact of tire, obtain the tire three-dimensional finite element model virtual simulation test platform applying fixed constraint,
Fixed constraint is applied to tire centerline, refers to that the degree of freedom controlling tire centerline 6 directions is 0;
In step 6, the tire three-dimensional finite element model virtual simulation test platform of applying fixed constraint that obtains in step 5, by horizontal supporting plate, tire lower surface is applied to the static pressure of vertical direction, measure squeegee action, obtain with the static pressure of vertical direction as ordinate, squeegee action amount is the relation curve of horizontal ordinate, and the curve that in curve, the tangent slope at each point place is formed is the vertical non-linear rigidity curve of tire;
Step 7, in the tire three-dimensional finite element model virtual simulation test platform of the applying fixed constraint obtained in step 5, the vertical freedom degree of release tire centerline, after the vertical range be raised to by the tire three-dimensional finite element model imparting tire authentic material parameter apart from horizontal supporting plate is at least the 1/12mm of diameter of tyres d, freely fall on horizontal supporting plate, measure the vertical free vibration response of tire three-dimensional finite element model central point displacement, obtain the free vibration attenuation curve of tire centerline point displacement, the vertical damping ratio of tire is calculated by free vibration attenuation curve.
The present invention is based on pulse testing obtains in the method for off-the-road tyre Vertical Characteristic parameter, described diameter of tyres d >=1500mm.
The present invention is based on pulse testing obtains in the method for off-the-road tyre Vertical Characteristic parameter, the constraint of real vehicle structure to tire comprise real vehicle suspension leverage, main shaft, steering linkage to the constraint of tire degree of freedom and real vehicle deadweight the load restraint to tire.
The present invention is based on pulse testing obtains in the method for off-the-road tyre Vertical Characteristic parameter, and the constraint of tire faced by real vehicle operation time channel refers to the contact of road surface and tire, by rigid contact process.
The present invention is based on the method that pulse testing obtains off-the-road tyre Vertical Characteristic parameter, in step 3, impulse response section refers in the accelerating curve that pulse operation Virtual Experiment Simulation obtains, and contacts to the accelerating curve section corresponding to the time period departed from triangular pulse BOB(beginning of block) with tire.
The present invention is based on the method that pulse testing obtains off-the-road tyre Vertical Characteristic parameter, in step 3, adopt remaining data of not participating in the accelerating curve of training in tire centerline point vertical acceleration curve to verify the neural network trained.
The present invention is based on the method that pulse testing obtains off-the-road tyre Vertical Characteristic parameter, in step 3, before the accelerating curve section of paired pulses response section carries out decile, the accelerating curve section of paired pulses response section carries out acceleration signal denoising.
The present invention is based on the method that pulse testing obtains off-the-road tyre Vertical Characteristic parameter, the accelerating curve section of paired pulses response section carries out denoising, the wavelet toolbox in matlab software is used to carry out wavelet analysis to described segment of curve, select db13 small echo, Decomposition order is 5 layers, carries out soft-threshold denoising.
Compared to prior art, advantage of the present invention is as follows:
1, the test that the present invention carries out is real train test, is not therefore subject to the constraint of tire construction size and institute's bearing load size, as long as this tire has practical application;
2, the present invention utilizes experiment to obtain tyre accelerations family curve, can obtain tire parameter, not need to carry out special on-the-spot bench test to tire, also do not need to carry out independent test to the component materials of tire, effectively can reduce experimentation cost;
3, the present invention, considers the interaction between tire material each component in finite element model, and the component materials parameter that the material property testing that the tire material parameter therefore obtained by discrimination method is comparatively separated records will be more realistic;
4, the present invention, on the basis of the material parameter obtained, can to obtain under different charge pressure tire vertical stiffness damping parameter comparatively accurately further, and can further investigate dynamic perfromances such as tire natural frequency, the vibration shapes, the parameter physical significance obtained definitely.
Accompanying drawing explanation
Accompanying drawing 1 is the inventive method process flow diagram;
Accompanying drawing 2 is tire cross-section profile and material distribution drawing;
Accompanying drawing 3 is tire cross-section two dimensional finite element grid chart;
Accompanying drawing 4 is the tire three-dimensional finite element model pressing material character subregion;
Accompanying drawing 5 is tire-road coupling three-dimensional finite element model;
The 1 tires central point vertical acceleration curve that accompanying drawing 6 obtains for pulse operation Virtual Experiment Simulation;
Accompanying drawing 7 is BP schematic network structure;
Accompanying drawing 8 is the accelerating curve after impulse response section in accompanying drawing 6 takes soft-threshold denoising;
Accompanying drawing 9 is actual vehicle pulse testing, the tire centerline impulse response section vertical acceleration curve obtained;
Accompanying drawing 10 is squeegee action and static pressure relation curve;
Accompanying drawing 11 is tire free vibration attenuation curve;
Embodiment
Embodiment 1
See accompanying drawing 1-11;
Obtain the method for off-the-road tyre Vertical Characteristic parameter based on pulse testing, its process flow diagram, as Fig. 1, comprises the steps:
Step one, set up tire-road coupling three-dimensional finite element model
Described diameter of tyres 3866mm.
First, according to tire size Plotting data tire CAD cross section profile figure, see accompanying drawing 2, by forming tire material different in kind, being wire loop, casing ply, belt, rubber layer four parts by tire CAD cross section profile diagram root, in described four parts, is the length of side by 5-10mm, in each part, divide quadrilateral mesh respectively, see accompanying drawing 3; By tire CAD cross section profile figure around tire centerline axle rotating 360 degrees, obtain the tire three-dimensional finite element model by material character subregion, see accompanying drawing 4; Contact is set up by adhesion mode, tire pressure 1Mpa between heterogeneity material area;
Secondly, connect leverage with real vehicle and boundary condition is constrained to tire, tire three-dimensional finite element model is retrained, obtains the tire three-dimensional finite element model connecting leverage constraint;
Real vehicle connect the constraint of leverage to tire comprise real vehicle suspension leverage, main shaft, steering linkage to the constraint of tire degree of freedom and real vehicle deadweight the load restraint to tire;
Finally, run the constraint of time channel in the face of tire as boundary condition using real vehicle, the tire three-dimensional finite element model connecting leverage constraint is retrained, obtain tire-road coupling three-dimensional finite element model; The constraint of tire faced by real vehicle operation time channel refers to the contact of road surface and tire, by rigid contact process;
Step 2, pulse operation Virtual Experiment Simulation
Triangular pulse block is set in the tire-road coupling three-dimensional finite element model that step one obtains, obtains the tire-road coupling three-dimensional finite element model comprising pulse packet, see accompanying drawing 5; Triangular pulse tile height size h is:
258≤h≤322mm; Width=the tyre width of triangular pulse block.
According to the wire loop, casing ply, belt, the Poisson ratio of rubber layer, the empirical data of elastic modulus that form tire, determine the Poisson ratio of wire loop, casing ply, belt, rubber layer, the base value of elastic modulus respectively, in table 1;
The Poisson ratio of wire loop, casing ply, belt, rubber layer, the base value of elastic modulus ± 20% scope in, optimization Latin Hypercube Sampling method is adopted to choose at least 40 group wire loops, casing ply, belt, the Poisson ratio of rubber layer, the numerical value of elastic modulus, in table 1, as the input parameter of pulse operation Virtual Experiment Simulation; The input parameter of every set of pulses operating mode Virtual Experiment Simulation is substituted in the tire-road coupling three-dimensional finite element model comprising triangular pulse block, pulse operation Virtual Experiment Simulation is carried out with travelling speed 30km/h, obtain 40 tires central point vertical acceleration curves, accompanying drawing 6 is an accelerating curve wherein; Travelling speed V span is 30Km/h;
Table 1
Step 3, neural metwork training
The type selecting neural network is BP network, comprise an input layer, input layer number is 103, a hidden layer, node in hidden layer is 10, an output layer, output layer nodes is 8, n node of input layer represents n accekeration respectively, 8 nodes of output layer represent wire loop in tire finite element model respectively, casing ply, belt, the Poisson ratio of rubber layer and elastic modulus input layer are tanh (tansig) function to the transport function of hidden layer, hidden layer is non-linear (sigmoid) transport function to the transport function of output layer, see accompanying drawing 7,
At least 35 curves in the tire centerline point vertical acceleration curve selecting step 2 to obtain, 102 deciles are carried out to the accelerating curve section of the impulse response section of each selected curve, before carrying out decile, the accelerating curve section of paired pulses response section carries out denoising, uses the wavelet toolbox in matlab software to carry out wavelet analysis to described segment of curve, select db13 small echo during denoising, Decomposition order is 5 layers, take soft-threshold denoising, obtain the accelerating curve after after denoising, see accompanying drawing 8.Carry out decile after denoising and obtain tire centerline point 103 vertical acceleration values in impulse response section, using the input parameter of 103 vertical acceleration values as BP network, after the BP network operation, obtain network output valve, repetition training is carried out to BP network; Impulse response section refers in the accelerating curve that pulse operation Virtual Experiment Simulation obtains, and contacts to the accelerating curve section corresponding to the time period departed from tire with triangular pulse BOB(beginning of block);
In the pulse Virtual Experiment Simulation tire finite element model corresponding with corresponding vertical acceleration curve when the output valve of BP network the Poisson ratio of wire loop, casing ply, belt, rubber layer and elastic modulus input value between error≤0.5-1% time, complete neural metwork training; Obtain the neural network trained;
Remaining data of not participating in the accelerating curve of training in tire centerline point vertical acceleration curve are adopted to verify the neural network trained;
Step 4, acquisition tire material parameter
By travelling speed 30km/h, pulse testing is carried out to actual vehicle, obtains tire centerline vertical acceleration impulse response section curve, see accompanying drawing 9; Then, the accelerating curve section of impulse response section in the tire centerline vertical acceleration curve of acquisition is carried out 102 deciles, obtain tire centerline point 103 vertical acceleration values in impulse response section, using the input parameter of 103 vertical acceleration values as BP network, in the neural network that substitution step 3 trains, obtain wire loop in actual tire, casing ply, belt, the Poisson ratio of rubber layer, elastic modulus, in table 2;
Table 2
| Title | Elastic modulus/MPa | Poisson ratio |
| Wire loop | 1.703×10 5 | 0.3565 |
| Casing ply | 1159 | 0.4134 |
| Belt | 1.120×10 4 | 0.3608 |
| Rubber layer | 11.66 | 0.4975 |
Step 5, set up virtual simulation test platform
In actual tire step 4 obtained, the Poisson ratio of wire loop, casing ply, belt, rubber layer, elastic modulus substitute into the tire three-dimensional finite element model by material character subregion that step one obtains, obtain the tire three-dimensional finite element model imparting tire authentic material parameter, in the tire three-dimensional finite element model imparting tire authentic material parameter, fixed constraint is applied to tire centerline, simultaneously, set up the horizontal supporting plate with the vertical lower surface rigid contact of tire, obtain the tire three-dimensional finite element model virtual simulation test platform applying fixed constraint;
In step 6, the tire three-dimensional finite element model virtual simulation test platform of applying fixed constraint that obtains in step 5, by horizontal supporting plate, tire lower surface is applied to the static pressure of vertical direction, measure squeegee action, obtain with the static pressure of vertical direction as ordinate, squeegee action amount is the relation curve of horizontal ordinate, see accompanying drawing 10, the curve that in curve, the tangent slope at each point place is formed is the vertical non-linear rigidity curve of tire;
Step 7, in the tire three-dimensional finite element model virtual simulation test platform of the applying fixed constraint obtained in step 5, the vertical freedom degree of release tire centerline, after the vertical range be raised to by the tire three-dimensional finite element model imparting tire authentic material parameter apart from horizontal supporting plate is at least the 1/12mm of diameter of tyres d, freely fall on horizontal supporting plate, measure the vertical free vibration response of tire three-dimensional finite element model central point displacement, obtain the free vibration attenuation curve of tire centerline point displacement, see accompanying drawing 11, the vertical damping ratio of tire is calculated by free vibration attenuation curve,
Tire three-dimensional finite element model central point displacement measurement is monitored automatically by simulation software.
Claims (10)
1. obtain the method for off-the-road tyre Vertical Characteristic parameter based on pulse testing, comprise the steps:
Step one, set up tire-road coupling three-dimensional finite element model
First, according to the tire size Plotting data tire CAD cross section profile figure of design, by forming tire material different in kind, be wire loop, casing ply, belt, rubber layer four parts by tire CAD cross section profile diagram root, in described four parts, grid division respectively; By tire CAD cross section profile figure around tire centerline axle rotating 360 degrees, obtain the tire three-dimensional finite element model by material character subregion; Contact is set up by adhesion mode, tire pressure 1Mpa between heterogeneity material area;
Secondly, connect leverage with real vehicle and boundary condition is constrained to tire, tire three-dimensional finite element model is retrained, obtains the tire three-dimensional finite element model connecting leverage constraint;
Finally, run the constraint of time channel in the face of tire as boundary condition using real vehicle, the tire three-dimensional finite element model connecting leverage constraint is retrained, obtain tire-road coupling three-dimensional finite element model;
Step 2, pulse operation Virtual Experiment Simulation
Triangular pulse block is set in the tire-road coupling three-dimensional finite element model that step one obtains, obtains the tire-road coupling three-dimensional finite element model comprising pulse packet;
According to the wire loop, casing ply, belt, the Poisson ratio of rubber layer, the empirical data of elastic modulus that form tire, determine the Poisson ratio of wire loop, casing ply, belt, rubber layer, the base value of elastic modulus respectively
The Poisson ratio of wire loop, casing ply, belt, rubber layer, the base value of elastic modulus ± 10-20% within the scope of, optimization Latin Hypercube Sampling method is adopted to choose at least 40 group wire loops, casing ply, belt, the Poisson ratio of rubber layer, the numerical value of elastic modulus, as the input parameter of pulse operation Virtual Experiment Simulation; The input parameter of every set of pulses operating mode Virtual Experiment Simulation is substituted in the tire-road coupling three-dimensional finite element model comprising triangular pulse block, carry out pulse operation Virtual Experiment Simulation with travelling speed V, obtain at least 40 tires central point vertical acceleration curves; Travelling speed V span is 20-60Km/h;
Step 3, neural metwork training
The type selecting neural network is BP network;
At least 35 curves in the tire centerline point vertical acceleration curve selecting step 2 to obtain, n decile is carried out to the accelerating curve section of the impulse response section of each selected curve, obtain tire centerline point n vertical acceleration value at least 35 group pulse response sections, to often organize the input amendment of n vertical acceleration value as BP network, after the BP network operation, obtain network output valve;
Repetition training is carried out to BP network;
During error≤0.5-1% between the Poisson ratio of the wire loop of the pulse Virtual Experiment Simulation tire finite element model corresponding with corresponding vertical acceleration curve when the output valve of BP network, casing ply, belt, rubber layer and the input value of elastic modulus, complete neural metwork training; Obtain the neural network trained;
Step 4, acquisition tire material parameter
By travelling speed V, pulse testing is carried out to actual vehicle, obtains tire centerline vertical acceleration curve; Then, the accelerating curve section of impulse response section in the tire centerline vertical acceleration curve of acquisition is carried out n decile, obtain tire centerline point n vertical acceleration value in impulse response section, using the input parameter of n vertical acceleration value as BP network, in the BP neural network that substitution step 3 trains, the output valve of network is wire loop in actual tire, casing ply, belt, the Poisson ratio of rubber layer, elastic mould value;
Step 5, set up virtual simulation test platform
First, wire loop in the actual tire that step 4 is obtained, casing ply, belt, the Poisson ratio of rubber layer, elastic mould value substitutes into the tire three-dimensional finite element model by material character subregion that step one obtains, obtain the tire three-dimensional finite element model imparting tire authentic material parameter, then, in the tire three-dimensional finite element model imparting tire authentic material parameter, fixed constraint is applied to tire centerline, simultaneously, set up the horizontal supporting plate with the vertical lower surface rigid contact of tire, obtain the tire three-dimensional finite element model virtual simulation test platform applying fixed constraint,
In step 6, the tire three-dimensional finite element model virtual simulation test platform of applying fixed constraint that obtains in step 5, by horizontal supporting plate, tire lower surface is applied to the static pressure of vertical direction, measure squeegee action, obtain with the static pressure of vertical direction as ordinate, squeegee action amount is the relation curve of horizontal ordinate, and the curve that in curve, the tangent slope at each point place is formed is the vertical non-linear rigidity curve of tire;
Step 7, in the tire three-dimensional finite element model virtual simulation test platform of the applying fixed constraint obtained in step 5, the vertical freedom degree of release tire centerline, after the vertical range be raised to by the tire three-dimensional finite element model imparting tire authentic material parameter apart from horizontal supporting plate is at least the 1/12mm of diameter of tyres d, freely fall on horizontal supporting plate, measure the vertical free vibration response of tire three-dimensional finite element model central point displacement, obtain the free vibration attenuation curve of tire centerline point displacement, the vertical damping ratio of tire is calculated by free vibration attenuation curve.
2. the method obtaining off-the-road tyre Vertical Characteristic parameter based on pulse testing according to claim 1, is characterized in that: in step one, and grid is that the length of side divides by 5-10mm.
3. the method obtaining off-the-road tyre Vertical Characteristic parameter based on pulse testing according to claim 1, is characterized in that: described diameter of tyres d >=1500mm.
4. the method obtaining off-the-road tyre Vertical Characteristic parameter based on pulse testing according to claim 1, it is characterized in that: in step 2, triangular pulse tile height size h is:
D/15≤h≤d/12, d is diameter of tyres; Width >=the tyre width of triangular pulse block.
5. the method obtaining off-the-road tyre Vertical Characteristic parameter based on pulse testing according to claim 1, is characterized in that: the constraint of real vehicle structure to tire comprise real vehicle suspension leverage, main shaft, steering linkage to the constraint of tire degree of freedom and real vehicle deadweight the load restraint to tire.
6. the method obtaining off-the-road tyre Vertical Characteristic parameter based on pulse testing according to claim 1, is characterized in that: the constraint of tire faced by real vehicle operation time channel refers to the contact of road surface and tire, by rigid contact process.
7. the method obtaining off-the-road tyre Vertical Characteristic parameter based on pulse testing according to claim 1, it is characterized in that: in step 3, impulse response section refers in the tire centerline point vertical acceleration curve that pulse operation Virtual Experiment Simulation obtains, contact as starting point with tire with triangular pulse BOB(beginning of block), tire and triangular pulse block depart from the accelerating curve section corresponding to the time period of terminal.
8. the method obtaining off-the-road tyre Vertical Characteristic parameter based on pulse testing according to claim 1, it is characterized in that: in step 3, paired pulses response section accelerating curve section carry out decile before, paired pulses response section accelerating curve section carry out acceleration signal denoising.
9. the method obtaining off-the-road tyre Vertical Characteristic parameter based on pulse testing according to claim 8, it is characterized in that: the accelerating curve section of paired pulses response section carries out denoising, the wavelet toolbox in matlab software is used to carry out wavelet analysis to described segment of curve, select db13 small echo, Decomposition order is 5 layers, carries out soft-threshold denoising.
10. the method obtaining off-the-road tyre Vertical Characteristic parameter based on pulse testing according to claim 1, is characterized in that: in step 5, applies fixed constraint to tire centerline, refers to that the degree of freedom controlling tire centerline 6 directions is 0.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510969939.3A CN105512424B (en) | 2015-12-22 | 2015-12-22 | The method that off-the-road tyre Vertical Characteristic parameter is obtained based on pulse testing |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510969939.3A CN105512424B (en) | 2015-12-22 | 2015-12-22 | The method that off-the-road tyre Vertical Characteristic parameter is obtained based on pulse testing |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105512424A true CN105512424A (en) | 2016-04-20 |
| CN105512424B CN105512424B (en) | 2018-01-26 |
Family
ID=55720404
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510969939.3A Expired - Fee Related CN105512424B (en) | 2015-12-22 | 2015-12-22 | The method that off-the-road tyre Vertical Characteristic parameter is obtained based on pulse testing |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105512424B (en) |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107066715A (en) * | 2017-03-31 | 2017-08-18 | 柳州顺联科技有限公司 | Tire Optimization Design in engineering machinery |
| CN108416178A (en) * | 2018-05-02 | 2018-08-17 | 南京航空航天大学 | A kind of mechanical elastic vehicle wheel Design of Structural parameters method |
| CN109670266A (en) * | 2018-12-29 | 2019-04-23 | 青岛慧拓智能机器有限公司 | A method of for creating mine vehicle kinetic model |
| CN109982872A (en) * | 2016-10-05 | 2019-07-05 | 索莱拉控股股份有限公司 | Vehicle tyre monitors system and method |
| CN110751633A (en) * | 2019-10-11 | 2020-02-04 | 上海眼控科技股份有限公司 | Multi-axis cart braking detection method, device and system based on deep learning |
| CN111624004A (en) * | 2020-06-11 | 2020-09-04 | 哈尔滨工业大学 | Rapid prediction method for radial tire braking distance |
| CN112287492A (en) * | 2020-12-28 | 2021-01-29 | 四川大学 | Dynamic characteristic identification method for five-axis linkage machining center swinging turntable |
| CN114279659A (en) * | 2020-09-27 | 2022-04-05 | 中车株洲电力机车研究所有限公司 | Hinge system test method, system, medium and equipment |
| CN116522510A (en) * | 2023-06-28 | 2023-08-01 | 中信戴卡股份有限公司 | Inversion method and system for tire material parameters for wheel performance simulation |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090012689A1 (en) * | 2007-07-04 | 2009-01-08 | Michelin Recherche Et Technique S.A. | Method of estimating an available grip margin of a tire when rolling |
| CN101840449A (en) * | 2010-04-13 | 2010-09-22 | 北京农业信息技术研究中心 | Tyre stress simulation method and system thereof |
| CN103473395A (en) * | 2013-08-21 | 2013-12-25 | 奇瑞汽车股份有限公司 | Finite element analysis method for simulating ride comfort of entire vehicle under pulse input |
| CN104390794A (en) * | 2014-12-02 | 2015-03-04 | 吉林大学 | Method of predicating flat-pavement mechanical properties of tyre based on test data of rotary drum test bed |
| CN105138796A (en) * | 2015-09-16 | 2015-12-09 | 重庆长安汽车股份有限公司 | Modal tire modeling method for whole-vehicle vibration noise simulation |
-
2015
- 2015-12-22 CN CN201510969939.3A patent/CN105512424B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090012689A1 (en) * | 2007-07-04 | 2009-01-08 | Michelin Recherche Et Technique S.A. | Method of estimating an available grip margin of a tire when rolling |
| CN101840449A (en) * | 2010-04-13 | 2010-09-22 | 北京农业信息技术研究中心 | Tyre stress simulation method and system thereof |
| CN103473395A (en) * | 2013-08-21 | 2013-12-25 | 奇瑞汽车股份有限公司 | Finite element analysis method for simulating ride comfort of entire vehicle under pulse input |
| CN104390794A (en) * | 2014-12-02 | 2015-03-04 | 吉林大学 | Method of predicating flat-pavement mechanical properties of tyre based on test data of rotary drum test bed |
| CN105138796A (en) * | 2015-09-16 | 2015-12-09 | 重庆长安汽车股份有限公司 | Modal tire modeling method for whole-vehicle vibration noise simulation |
Non-Patent Citations (5)
| Title |
|---|
| GUARNERI P，ROCCA G，GOBBI M．: "A neural-network-based model for the dynamic simulation of the tire／suspension system while traversing road irregularities", 《NEURAL NETWORKS，IEEE TRANSACTIONS》 * |
| 于德介，雷慧: "一种基于神经网络的结构参数识别方法", 《湖南大学学报(自然科学版)》 * |
| 束晓琳，李红，郫孔辉: "基于改进的粒子群优化算法的轮胎参数", 《科技导报》 * |
| 王和毅, 谷正气: "汽车轮胎模型研究现状及其发展分析", 《橡胶工业》 * |
| 谷正气等: "基于ＧＡ－ＢＰ网络的矿山路面不平度辨识", 《中国机械工程》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109982872A (en) * | 2016-10-05 | 2019-07-05 | 索莱拉控股股份有限公司 | Vehicle tyre monitors system and method |
| CN107066715A (en) * | 2017-03-31 | 2017-08-18 | 柳州顺联科技有限公司 | Tire Optimization Design in engineering machinery |
| CN108416178A (en) * | 2018-05-02 | 2018-08-17 | 南京航空航天大学 | A kind of mechanical elastic vehicle wheel Design of Structural parameters method |
| CN109670266A (en) * | 2018-12-29 | 2019-04-23 | 青岛慧拓智能机器有限公司 | A method of for creating mine vehicle kinetic model |
| CN109670266B (en) * | 2018-12-29 | 2023-04-18 | 青岛慧拓智能机器有限公司 | Method for creating mining vehicle dynamics model |
| CN110751633A (en) * | 2019-10-11 | 2020-02-04 | 上海眼控科技股份有限公司 | Multi-axis cart braking detection method, device and system based on deep learning |
| CN111624004A (en) * | 2020-06-11 | 2020-09-04 | 哈尔滨工业大学 | Rapid prediction method for radial tire braking distance |
| CN114279659A (en) * | 2020-09-27 | 2022-04-05 | 中车株洲电力机车研究所有限公司 | Hinge system test method, system, medium and equipment |
| CN114279659B (en) * | 2020-09-27 | 2023-06-23 | 中车株洲电力机车研究所有限公司 | Hinge system test method, system, medium and equipment |
| CN112287492A (en) * | 2020-12-28 | 2021-01-29 | 四川大学 | Dynamic characteristic identification method for five-axis linkage machining center swinging turntable |
| CN116522510A (en) * | 2023-06-28 | 2023-08-01 | 中信戴卡股份有限公司 | Inversion method and system for tire material parameters for wheel performance simulation |
| CN116522510B (en) * | 2023-06-28 | 2023-09-01 | 中信戴卡股份有限公司 | Inversion method and system for tire material parameters for wheel performance simulation |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105512424B (en) | 2018-01-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105512424A (en) | Method for obtaining engineering machinery tyre vertical characteristic parameters based on pulse test | |
| CN110298125B (en) | Fatigue analysis method based on virtual test field technology | |
| CN103942392B (en) | Automotive chassis technical parameter robust design method based on full life circle | |
| CN109426685B (en) | Component fatigue analysis method | |
| CN111159938B (en) | Underwater vibration noise prediction method for river-crossing tunnel of road | |
| CN106441531A (en) | Dynamic weighing method and system on condition of uniform motion of vehicle | |
| CN111737816A (en) | Lightweight design method of non-inflatable explosion-proof wheel | |
| CN112347541B (en) | Method for rapidly testing bearing capacity of asphalt pavement under normal traffic of vehicle | |
| CN108228994B (en) | Stress calculation method for vehicle and working device under random excitation of off-road | |
| CN103759954A (en) | Method and device for precisely testing rolling resistance of tyres | |
| CN104458287A (en) | Tire lateral relaxation length measurement method | |
| CN103164577B (en) | Determine the method for harbour bridge automobile load rated figure | |
| Lu et al. | Analysis of asphalt pavement mechanical behaviour by using a tire-pavement coupling model | |
| CN114417630A (en) | Suspension design, simulation and analysis integrated method based on MWorks | |
| CN105973189A (en) | Method for obtaining surface roughness load spectrum of soft road surface | |
| Marjani | Development of FEA wide-base truck tire and soil interaction models | |
| Su | A road load data processing technique for durability optimization of automotive products | |
| CN109918843B (en) | Method for obtaining vibration compaction value of non-adhesive roadbed soil based on discrete element method | |
| Ragheb et al. | Development of a combat vehicle FEA tire model for off-road applications | |
| CN104166778A (en) | Method for optimizing rigidity and damping of contact part of automobile seat guide rail | |
| CN110955933B (en) | Mechanical structure fuzzy fatigue reliability calculation method based on response surface method | |
| Tong et al. | Study on the simulation of radial tire wear Characteristics | |
| Delvecchio et al. | Estimation of the tire vertical deflection in a motorcycle suspension via Kalman-filtering techniques | |
| CN104330269A (en) | Method for testing damping ratio of tyre-rim combination body | |
| Anghelache et al. | Truck tyre-road contact stress measurement and modelling |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180126 Termination date: 20211222 |