CN114861335B - Calibration method of automobile dynamics calculation model and related equipment - Google Patents

Calibration method of automobile dynamics calculation model and related equipment Download PDF

Info

Publication number
CN114861335B
CN114861335B CN202210808083.1A CN202210808083A CN114861335B CN 114861335 B CN114861335 B CN 114861335B CN 202210808083 A CN202210808083 A CN 202210808083A CN 114861335 B CN114861335 B CN 114861335B
Authority
CN
China
Prior art keywords
load
model
calibration
wheel center
calculation 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.)
Active
Application number
CN202210808083.1A
Other languages
Chinese (zh)
Other versions
CN114861335A (en
Inventor
丁鼎
韩广宇
张永仁
卢放
赵永钊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Lantu Automobile Technology Co Ltd
Original Assignee
Lantu Automobile Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Lantu Automobile Technology Co Ltd filed Critical Lantu Automobile Technology Co Ltd
Priority to CN202210808083.1A priority Critical patent/CN114861335B/en
Publication of CN114861335A publication Critical patent/CN114861335A/en
Application granted granted Critical
Publication of CN114861335B publication Critical patent/CN114861335B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

The application discloses a calibration method of an automobile dynamics calculation model and related equipment. The method comprises the following steps: acquiring a component actual measurement load and a component simulation load based on a target load applied by a wheel center; performing first calibration on the preset automobile dynamics calculation model according to the component actual measurement load and the component simulation load to obtain a first calibrated automobile dynamics calculation model, wherein the first calibration is used for calibrating the preset automobile dynamics calculation model except the tire model; acquiring wheel center load actual measurement data and wheel center load simulation data based on the test yard characteristic road surface; and carrying out second calibration on the tire model of the preset automobile dynamics calculation model after the first calibration is finished according to the wheel center load actual measurement data and the wheel center load simulation data to obtain a target automobile dynamics model. The method for calibrating the automobile dynamic model, provided by the embodiment of the application, enables the calibration work to be more efficient and accurate, has important engineering significance, and can improve the accuracy of calibrating the automobile dynamic model.

Description

Calibration method of automobile dynamics calculation model and related equipment
Technical Field
The present disclosure relates to the field of automobile simulation, and more particularly, to a calibration method for an automobile dynamics calculation model and related devices.
Background
The automobile can perform strength durability simulation analysis work on parts in a research and development stage to predict the strength durability of a design scheme, so that structural design improvement is performed in a data stage, and the test verification period and the test cost of the automobile can be shortened. In the part strength endurance simulation analysis, the strength endurance load of the part is the most important input, and the precision of the strength endurance load directly determines the precision of the part strength endurance simulation. And the precision of the strength and durability load is ensured by the precision of a dynamic load calculation model of the automobile multi-body dynamics. At present, a plurality of commercial software can realize the establishment of an automobile multi-body dynamic load calculation model, but how to set model parameters of the automobile multi-body dynamic load calculation model can ensure that the precision of the strength durable load meets the use requirement on engineering is always a difficult point of industry. Because the parameters of the calculation model are numerous, a plurality of parameters have great influence on the precision of the strength and durability load of the automobile, and if the precision of the strength and durability load cannot meet the requirement, the precision of the calculation model cannot be improved by depending on engineering experience or manually adjusting the parameters of the model so as to meet the requirement of engineering use.
Disclosure of Invention
In this summary, concepts in a simplified form are introduced that are further described in the detailed description. The summary of the invention is not intended to identify key features or essential features of the claimed subject matter, nor is it intended to be used as an aid in determining the scope of the claimed subject matter.
In order to improve the accuracy of the vehicle dynamics calculation model, in a first aspect, the present invention provides a calibration method of a vehicle dynamics calculation model, where the method includes:
acquiring a component actual measurement load and a component simulation load based on a target load applied by a wheel center;
performing first calibration on a preset automobile dynamics calculation model according to the component actual measurement load and the component simulation load to obtain a first calibrated automobile dynamics calculation model, wherein the first calibration is used for calibrating the preset automobile dynamics calculation model except the tire model;
acquiring wheel center load actual measurement data and wheel center load simulation data based on the test yard characteristic road surface;
and carrying out second calibration on the tire model of the preset automobile dynamics calculation model after the first calibration is finished according to the wheel center load actual measurement data and the wheel center load simulation data so as to obtain a target automobile dynamics model.
Optionally, the obtaining of the component actual measurement load and the component simulation load based on the target load applied by the wheel center includes:
acquiring the actual measurement load of the part measured by the automobile shaft coupling test bench test based on the target load applied to the wheel center;
and acquiring the part simulation load corresponding to the preset automobile dynamics calculation model under the condition that the target load is applied to the wheel center.
Optionally, the first calibration of the preset vehicle dynamics calculation model according to the actual measurement load of the component and the simulation load of the component to obtain a first calibrated vehicle dynamics calculation model includes:
acquiring a first load root mean square error of the actual measurement load of the component and the simulation load of the component;
and under the condition that the first load root mean square error is larger than a first preset root mean square error, carrying out first calibration on the preset automobile dynamics calculation model to obtain a first calibrated automobile dynamics calculation model.
Optionally, the first calibrating the preset vehicle dynamics calculation model to obtain a first calibrated vehicle dynamics calculation model includes:
under the condition that the first load root mean square error is larger than the first preset root mean square error, acquiring the model parameter sensitivity of the preset automobile dynamics calculation model;
and performing first calibration on the model parameters with the sensitivity greater than the first preset sensitivity by adopting a genetic algorithm to obtain a first calibrated automobile dynamics calculation model.
Optionally, the above-mentioned wheel center load measured data and wheel center load simulation data obtained based on the test yard characteristic road surface includes:
testing and obtaining the wheel center load actual measurement data on the characteristic road surface of the test yard based on the target vehicle;
and carrying out simulation calculation on the preset automobile dynamics calculation model completing the first calibration on the test yard characteristic road surface to obtain the wheel center load simulation data.
Optionally, the second calibrating the tire model of the preset vehicle dynamics calculation model completing the first calibrating according to the wheel center load actual measurement data and the wheel center load simulation data to obtain the target vehicle dynamics model includes:
acquiring a second load root mean square error of the wheel center load actual measurement data and the wheel center load simulation data;
and under the condition that the second load root mean square error is larger than a second preset root mean square error, performing second calibration on the tire model of the preset automobile dynamics calculation model subjected to the first calibration so as to obtain the second calibrated automobile dynamics calculation model.
Optionally, the second calibrating the tire model of the preset vehicle dynamics calculation model after the first calibrating is performed to obtain the second calibrated vehicle dynamics calculation model includes:
acquiring tire parameter sensitivity;
and performing second calibration on the tire parameters with the tire parameter sensitivity greater than a second preset sensitivity by adopting a genetic algorithm to obtain a target automobile dynamics calculation model.
In a second aspect, the present invention further provides a calibration apparatus for a vehicle dynamics calculation model, including:
a first acquisition unit configured to acquire a component actual measurement load and a component simulation load based on a target load applied by a wheel center;
the first calibration unit is used for carrying out first calibration on a preset automobile dynamics calculation model according to the component actual measurement load and the component simulation load, wherein the first calibration is used for calibrating the preset automobile dynamics calculation model except the tire model;
the second acquisition unit is used for acquiring wheel center load actual measurement data and wheel center load simulation data based on the test yard characteristic road surface;
and the second calibration unit is used for carrying out second calibration on the tire model of the preset automobile dynamics calculation model after the first calibration is finished according to the wheel center load actual measurement data and the wheel center load simulation data so as to obtain a target automobile dynamics model.
In a third aspect, an electronic device includes: a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor is configured to implement the steps of the method for calibrating a computational model of vehicle dynamics as described in any one of the first aspect above when the computer program stored in the memory is executed.
In a fourth aspect, the present invention further provides a computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the calibration method for the vehicle dynamics calculation model according to any one of the first aspect.
To sum up, the calibration method for the vehicle dynamics calculation model provided by the embodiment of the application comprises the following steps: acquiring a component actual measurement load and a component simulation load based on a target load applied by a wheel center; performing first calibration on a preset automobile dynamics calculation model according to the component actual measurement load and the component simulation load to obtain a first calibrated automobile dynamics calculation model, wherein the first calibration is used for calibrating the preset automobile dynamics calculation model except the tire model; acquiring wheel center load actual measurement data and wheel center load simulation data based on the test yard characteristic road surface; and carrying out second calibration on the tire model of the preset automobile dynamics calculation model after the first calibration is finished according to the wheel center load actual measurement data and the wheel center load simulation data so as to obtain a target automobile dynamics model. According to the automobile dynamics model calibration method provided by the embodiment of the application, calibration work is divided into two parts by means of the whole automobile shaft coupling test bed, and the preset automobile dynamics calculation model except the tire model and the model parameters of the tire are calibrated respectively, so that the calibration work is more efficient and accurate. According to the automobile multi-body dynamic load calculation model calibrated by the scheme, accurate strength durable load can be obtained, so that the accuracy of strength durable simulation input of parts is ensured, and the method has important engineering significance. The method for automatically calibrating the automobile dynamic model is provided, the workload of calibration work of the automobile multi-body dynamic load calculation model is greatly reduced, and the accuracy of calibration of the automobile dynamic model can be improved by the scheme.
Additional advantages, objects, and features of the invention will be set forth in part in the description which follows and in part will become apparent to those having ordinary skill in the art upon examination of the following or may be learned from practice of the invention.
Drawings
Various additional advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the specification. Also, like reference numerals are used to refer to like parts throughout the drawings. In the drawings:
fig. 1 is a schematic flowchart of a calibration method for an automotive dynamics calculation model according to an embodiment of the present application;
FIG. 2 is a schematic structural diagram of a calibration apparatus of an automotive dynamics calculation model according to an embodiment of the present disclosure;
fig. 3 is a schematic structural diagram of an electronic device of a calibration method for an automotive dynamics calculation model according to an embodiment of the present application.
Detailed Description
According to the automobile dynamics model calibration method provided by the embodiment of the application, calibration work is divided into two parts by means of the whole automobile shaft coupling test bed, and the preset automobile dynamics calculation model except the tire model and the model parameters of the tire are calibrated respectively, so that the calibration work is more efficient and accurate. The automobile multi-body dynamic load calculation model calibrated according to the scheme can obtain accurate strength durable load, so that the accuracy of strength durable simulation input of parts is ensured, and the method has important engineering significance. The method for realizing the automatic calibration of the automobile dynamic model is provided, the workload of the calibration work of the automobile multi-body dynamic load calculation model is greatly reduced, and the accuracy of the calibration of the automobile dynamic model can be improved by the scheme.
The terms "first," "second," "third," "fourth," and the like in the description and in the claims of the present application and in the drawings described above, if any, are used for distinguishing between similar elements and not necessarily for describing a particular sequential or chronological order. It will be appreciated that the data so used may be interchanged under appropriate circumstances such that the embodiments described herein may be practiced otherwise than as specifically illustrated or described herein. Furthermore, the terms "comprises," "comprising," and "having," and any variations thereof, are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of steps or elements is not necessarily limited to those steps or elements expressly listed, but may include other steps or elements not expressly listed or inherent to such process, method, article, or apparatus. The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all of the embodiments.
Referring to fig. 1, a schematic flow chart of a calibration method for an automotive dynamics calculation model provided in an embodiment of the present application may specifically include:
s110, acquiring a component actual measurement load and a component simulation load based on a target load applied by a wheel center;
illustratively, a measured measurement of the vehicle, i.e., a measured component load, is obtained by the entire axle coupling test stand based on a target load applied at the wheel center. And applying the same target load in a preset automobile dynamics calculation model, and acquiring the corresponding part simulation load. The target load comprises a force load signal, a torque load signal and a displacement load signal;
s120, performing first calibration on a preset automobile dynamics calculation model according to the component actual measurement load and the component simulation load, wherein the first calibration is used for calibrating the preset automobile dynamics calculation model except the tire model;
for example, a first calibration based on the component measured loads and the component simulated loads will be performed without consideration of a preset vehicle dynamics model of the vehicle tires. It will be appreciated that the first calibration is to calibrate other model parameters than the tire model parameters. It should be noted that the preset automotive multi-body dynamics calculation model may include a powertrain model, a front suspension model, a rear suspension model, a braking system model, a body system model, a steering system model, and an FTire tire model. The first calibration automobile dynamics calculation model preset automobile dynamics calculation model is a preset automobile multi-body dynamics model excluding a tire model. The first calibration of the model parameters of the preset automobile dynamics calculation model comprises the following steps: the numerical value of the hard point coordinate, the mass of the part, the barycenter coordinate of the part, the rotational inertia of the part, the spring stiffness, the damper damping, the stiffness and the damping of the elastic part, the clearance of the limiting block and the like.
S130, acquiring wheel center load actual measurement data and wheel center load simulation data based on the test yard characteristic road surface;
illustratively, the characteristic road surface of the test yard may include belgium, twisted, pebble, washboard, stone-strip, noise, resonance, etc., the wheel center load actual measurement data is obtained by a road test of the target vehicle on the characteristic road surface of the test yard, and the wheel center load simulation data is obtained by a simulation calculation of 3D digital data corresponding to the characteristic road surface of the test yard through a preset vehicle dynamics calculation model after a first calibration. The wheel center load measured data and wheel center load simulated data include each wheel center X, Y and forces and moments in the Z direction.
And S140, performing second calibration on the tire model of the preset automobile dynamics calculation model subjected to the first calibration according to the wheel center load actual measurement data and the wheel center load simulation data to obtain a target automobile dynamics model.
Illustratively, the tire model of the preset automobile dynamics calculation model which completes the first calibration is subjected to second calibration according to the wheel center load actual measurement data and the wheel center load simulation data, so that the error between the load simulation data of the tire model after calibration and the wheel center load actual measurement data is within an acceptable range, and the parameters of the tire model are optimized, thereby obtaining the target automobile dynamics model after the first calibration and the second calibration.
In summary, according to the method for calibrating the automobile dynamics model provided by the embodiment of the application, the calibration work is divided into two parts by virtue of the whole automobile shaft coupling test bed, and the calibration work is more efficient and accurate by respectively calibrating the preset automobile dynamics calculation model except the tire model and the model parameters of the tire. The automobile multi-body dynamic load calculation model calibrated according to the scheme can obtain accurate strength durable load, so that the accuracy of strength durable simulation input of parts is ensured, and the method has important engineering significance. The method for realizing the automatic calibration of the automobile dynamic model is provided, the workload of the calibration work of the automobile multi-body dynamic load calculation model is greatly reduced, and the accuracy of the calibration of the automobile dynamic model can be improved by the scheme.
In some examples, the obtaining of the component measured load and the component simulated load based on the target load applied by the wheel center includes:
acquiring the actual measurement load of the part measured by the automobile shaft coupling test bench test based on the target load applied to the wheel center;
and acquiring the part simulation load corresponding to the preset automobile dynamics calculation model under the condition that the target load is applied to the wheel center.
Illustratively, measured data of the vehicle, i.e., measured loads of the components, are measured by a vehicle axle coupling test stand. It should be noted that, the entire vehicle axle coupling test bed: the method applies a road simulation test technology, adopts a shaft coupling mode, and simultaneously provides 6 degrees of freedom for each angle of a passenger car: the input of vertical, horizontal, lateral, camber, steering and braking realizes the force and displacement simulation drive of the test yard road load in a laboratory, reproduces more than 90% of damage of the road, and performs relatively quick and comprehensive examination on the whole chassis and the structural member of the vehicle body. And corresponding sensors are arranged on the measured automobile to acquire wheel center load signals and component measured loads corresponding to the automobile. The wheel center load signal of the automobile is obtained by arranging a wheel center six-component sensor on the wheel center of the automobile, acquiring the six-component load of the wheel center of the automobile and defining the six-component load signal as the wheel center load signal of the automobile.
The component loads of the vehicle are defined as follows: an acceleration sensor is arranged on the wheel center of the automobile, an acceleration sensor is arranged at the installation point of the shock absorber of the automobile, and the like, so that the acceleration load signal of the part of the automobile can be obtained. Strain gauge sensors are arranged on parts such as a two-force rod piece of an automobile, a ball pin of a part, a guide rod of a shock absorber, a spring and the like, and a strain signal is calibrated into a force signal through a drawing press machine, so that a part force load signal of the automobile can be obtained. A stay wire displacement sensor and the like are arranged between a body and a chassis swing arm of the automobile, so that a component displacement load signal of the automobile can be obtained.
The preset automobile multi-body dynamics calculation model can comprise a power assembly model, a front suspension model, a rear suspension model, a braking system model, a vehicle body system model, a steering system model and an FTire tire model. The first calibrated automobile dynamics calculation model is a preset automobile multi-body dynamics model without tire model calibration.
Respectively applying the target load to the wheel center of the test vehicle in the automobile shaft coupling test bed and the wheel center of the preset automobile dynamics calculation model, which specifically includes:
and A, the target load comprises a force load signal, a torque load signal and a displacement load signal. Applying X-direction force load signal to wheel center of automobile
Figure 532467DEST_PATH_IMAGE001
: amplitude of
Figure 715187DEST_PATH_IMAGE002
In a frequency range of
Figure 822820DEST_PATH_IMAGE003
A white noise signal of (a); applying Y-direction force load signal to wheel center of automobile
Figure 280346DEST_PATH_IMAGE004
: amplitude of
Figure 124412DEST_PATH_IMAGE005
In a frequency range of
Figure 743612DEST_PATH_IMAGE006
The white noise signal of (a); applying Z-direction load transfer signal to wheel center of automobile
Figure 541804DEST_PATH_IMAGE007
: amplitude of
Figure 271863DEST_PATH_IMAGE008
In a frequency range of
Figure 471900DEST_PATH_IMAGE009
A white noise signal of (a); applying X-direction moment load signal to wheel center of automobile
Figure 996422DEST_PATH_IMAGE010
: amplitude of
Figure 547489DEST_PATH_IMAGE011
In a frequency range of
Figure 317124DEST_PATH_IMAGE012
The white noise signal of (a); applying Y-direction moment load signal to wheel center of automobile
Figure 371668DEST_PATH_IMAGE013
: amplitude of
Figure 535933DEST_PATH_IMAGE014
In a frequency range of
Figure 308717DEST_PATH_IMAGE015
A white noise signal of (a); applying Z-direction moment load signal to wheel center of automobile
Figure 380578DEST_PATH_IMAGE016
: amplitude of
Figure 555207DEST_PATH_IMAGE017
In a frequency range of
Figure 624795DEST_PATH_IMAGE018
White noise signal of (2). On the axle coupling test bed of the automobile, the load signal defined above is applied to the wheel center of the automobile, and the signal sensor arranged in the above step is combinedThe component load signal set of the automobile can be obtained through testing
Figure 884875DEST_PATH_IMAGE019
The expression is as follows:
Figure 291585DEST_PATH_IMAGE020
wherein:
Figure 842694DEST_PATH_IMAGE021
representing a set of component measured load signals;
Figure 145499DEST_PATH_IMAGE022
Figure 627296DEST_PATH_IMAGE023
Figure 40960DEST_PATH_IMAGE024
representing a component measured load signal obtained by testing;
Figure 924602DEST_PATH_IMAGE025
representing the total number of the measured load signals of the component obtained by the test.
B, applying the same target load in the step A to a preset automobile dynamics calculation model to calculate and obtain a part simulation load signal set of the automobile
Figure 929467DEST_PATH_IMAGE026
The expression is as follows:
Figure 101823DEST_PATH_IMAGE027
wherein:
Figure 53598DEST_PATH_IMAGE028
representing a part simulation load signal set obtained by calculation;
Figure 558791DEST_PATH_IMAGE029
Figure 203399DEST_PATH_IMAGE030
Figure 659788DEST_PATH_IMAGE031
representing the part simulation load signal of the automobile obtained by calculation;
Figure 680834DEST_PATH_IMAGE032
the total number of the part simulation load signals of the automobile obtained by the test is represented.
In summary, according to the method for calibrating the automobile dynamics model provided by the embodiment of the application, the target load is applied to the wheel center, the influence of the automobile tire on the calibration of the automobile dynamics is isolated, and the preset automobile dynamics calculation model outside the tire can be calibrated well.
In some examples, the first calibrating the preset vehicle dynamics calculation model according to the measured component load and the simulated component load to obtain a first calibrated vehicle dynamics calculation model includes:
acquiring a first load root mean square error of the actual measurement load of the component and the simulation load of the component;
and under the condition that the first load root mean square error is larger than a first preset root mean square error, carrying out first calibration on the preset automobile dynamics calculation model to obtain a first calibrated automobile dynamics calculation model.
Illustratively, the set of measured loads of the component is obtained from measurements
Figure 273489DEST_PATH_IMAGE033
And the simulated load set of the simulated component
Figure 88998DEST_PATH_IMAGE034
By calculating the first load mean square root error
Figure 235946DEST_PATH_IMAGE035
And has a first predetermined root mean square error
Figure 795103DEST_PATH_IMAGE036
Compare if, if
Figure 773424DEST_PATH_IMAGE037
If the calculation accuracy of the simulation load signal of the component meets the requirement of the preset accuracy, the preset automobile dynamics simulation model does not need to be calibrated firstly, and if the calculation accuracy of the simulation load signal of the component meets the requirement of the preset accuracy, the first calibration is not needed to be carried out on the preset automobile dynamics simulation model
Figure 992790DEST_PATH_IMAGE038
If the calculation accuracy of the simulated load signal of the component does not meet the requirement of the preset accuracy, the preset automobile dynamics simulation model needs to be calibrated firstly.
It should be noted that, in the following description,
Figure 423771DEST_PATH_IMAGE039
representing the simulated load signal curve of the part of the automobile obtained by solving calculation
Figure 724303DEST_PATH_IMAGE040
Obtaining a part test load signal curve of the automobile corresponding to the corresponding test
Figure 557130DEST_PATH_IMAGE041
The root mean square error function of;
Figure 448862DEST_PATH_IMAGE042
representing a first pre-set root mean square error,
in summary, the method for calibrating the automobile dynamics model provided by the embodiment of the application evaluates the actual measurement load and the simulation load of the part of the automobile through the root-mean-square error, and does not perform the first calibration under the condition of meeting the error requirement, so that the workload of the calibration work of the automobile multi-body dynamics dynamic load calculation model can be greatly reduced.
In some examples, the first calibrating the preset vehicle dynamics calculation model to obtain a first calibrated vehicle dynamics calculation model includes:
under the condition that the first load root mean square error is larger than the first preset root mean square error, acquiring the model parameter sensitivity of the preset automobile dynamics calculation model;
and performing first calibration on the model parameters with the sensitivity greater than the first preset sensitivity by adopting a genetic algorithm to obtain a first calibrated automobile dynamics calculation model.
Illustratively, if the first loading root mean square error is greater than the first predetermined root mean square error, i.e. the first loading root mean square error is greater than the first predetermined root mean square error
Figure 632719DEST_PATH_IMAGE043
At this time, it is stated that the calculation accuracy of the component simulation load signal does not meet the requirement of the preset accuracy, the parameter sensitivity in the preset automobile dynamics calculation model needs to be obtained, and the model parameter with the model parameter sensitivity greater than that of the first preset sensitivity model is subjected to first calibration by adopting a genetic algorithm to obtain a first calibrated automobile dynamics calculation model. Wherein the model parameters include: the coordinate value of the hard point, the mass of the part, the mass center coordinate of the part, the rotational inertia of the part, the rigidity of the spring, the damping of the shock absorber, the rigidity and the damping of the elastic part, the clearance of the limiting block and the like.
The specific operation of performing the first calibration may include:
s210, if
Figure 533679DEST_PATH_IMAGE044
If the load signal does not meet the requirement, the calculation accuracy of the load signal is judged to be not met; synchronously marking and reordering the part simulation load signals which do not meet the conditions, and defining the part simulation load signals as the load signals for improving the calculation precision by carrying out the first calibration of the parameters of the whole vehicle model subsequently, wherein the expression is as follows:
Figure 158695DEST_PATH_IMAGE045
Figure 486908DEST_PATH_IMAGE046
wherein:
Figure 626903DEST_PATH_IMAGE047
representing a part simulation load signal set of the automobile, which is obtained by calculation and needs to be subjected to precision improvement;
Figure 98598DEST_PATH_IMAGE048
to represent
Figure 640437DEST_PATH_IMAGE049
Correspondingly testing and obtaining an obtained actual measurement load signal set of the automobile component;
Figure 139552DEST_PATH_IMAGE050
Figure 970105DEST_PATH_IMAGE051
Figure 212867DEST_PATH_IMAGE052
to represent
Figure 609213DEST_PATH_IMAGE053
Calculating and obtaining a part simulation load signal of the automobile in the set;
Figure 544808DEST_PATH_IMAGE054
Figure 659395DEST_PATH_IMAGE055
Figure 643531DEST_PATH_IMAGE056
to represent
Figure 159963DEST_PATH_IMAGE057
And acquiring the obtained actual measurement load signal of the part of the automobile by the test in the set.
Figure 499416DEST_PATH_IMAGE058
Representing the total number of load signals.
S220, screening model parameters of a preset automobile multi-body dynamic load calculation model except the tire model: carrying out detailed review on the combination of the model parameters and the previous parameter test results; for some model parameters, because the reliability of the test result obtained by the test is high, the model parameters do not participate in calibration optimization in the following process, for example: model parameters such as component mass, spring stiffness and the like; the sensitivity can be determined for other model parameters
Figure 101298DEST_PATH_IMAGE059
The expression of (a) is as follows:
Figure 951442DEST_PATH_IMAGE060
wherein:
Figure 322381DEST_PATH_IMAGE061
expressing an absolute value function;
Figure 537462DEST_PATH_IMAGE062
a value representing the sensitivity of a model parameter to a load signal of a component of the vehicle;
Figure 361061DEST_PATH_IMAGE063
Figure 749317DEST_PATH_IMAGE064
Figure 974762DEST_PATH_IMAGE065
a larger value, a smaller value and a middle value representing the setting of the model parameters;
Figure 423061DEST_PATH_IMAGE066
Figure 501001DEST_PATH_IMAGE067
Figure 427368DEST_PATH_IMAGE068
the parameters of the representation model are respectively
Figure 710582DEST_PATH_IMAGE069
Figure 329782DEST_PATH_IMAGE070
And a component load signal of the automobile obtained by the calculation under P0.
S230, if the sensitivity of the model parameters to the load signals needing to be improved in calculation precision is less than or equal to a first preset sensitivity
Figure 127974DEST_PATH_IMAGE071
If the parameter is not sensitive to the load signal, the parameter does not participate in the subsequent calibration calculation.
S240, if the sensitivity of the model parameters to the load signals needing to be improved in calculation accuracy is larger than a first preset sensitivity
Figure 123612DEST_PATH_IMAGE071
It is said that the parameter is sensitive to the load signal and if the sensitivity of the model parameter to the load signal that has reached computational accuracy is less than a threshold constant
Figure 58070DEST_PATH_IMAGE072
Then the model parameters need to be calculated for the first calibration.
S250, first calibration: defining parameter sets of preset automobile multi-body dynamic load calculation models except tire models after model parameter screening is finished
Figure 848171DEST_PATH_IMAGE073
The expression is as follows:
Figure 336922DEST_PATH_IMAGE074
wherein:
Figure 605092DEST_PATH_IMAGE075
representing a parameter set of a preset automobile multi-body dynamic load calculation model except the tire model after model parameter screening is completed;
Figure 158171DEST_PATH_IMAGE076
Figure 119173DEST_PATH_IMAGE077
Figure 891957DEST_PATH_IMAGE078
parameters representing a computational model;
Figure 963819DEST_PATH_IMAGE079
representing the total number of computational model parameters.
And for parameter sets to
Figure 138448DEST_PATH_IMAGE080
Sets a value range, whose expression is as follows:
Figure 4773DEST_PATH_IMAGE081
wherein:
Figure 264853DEST_PATH_IMAGE082
a parameter representing a predetermined automotive multi-body dynamic load calculation model other than the tire model;
Figure 78088DEST_PATH_IMAGE083
Figure 107224DEST_PATH_IMAGE084
and expressing the minimum constant and the maximum constant of the parameters, and taking the value as the maximum value if the data is greater than the maximum value in subsequent iterative calculation, and taking the value as the minimum value if the data is less than the minimum value. Step S250 is embodiedStep S2501 to step S2508 may also be included.
S2501, setting the maximum iteration number as
Figure 645915DEST_PATH_IMAGE085
An iteration objective function vector F and an iteration termination error vector
Figure 658870DEST_PATH_IMAGE086
The expression formula is as follows:
Figure 72534DEST_PATH_IMAGE087
Figure 956176DEST_PATH_IMAGE088
Figure 429883DEST_PATH_IMAGE089
Figure 602238DEST_PATH_IMAGE090
Figure 554014DEST_PATH_IMAGE091
the iteration termination condition is as follows:
Figure 557742DEST_PATH_IMAGE092
Figure 936771DEST_PATH_IMAGE093
representing an iterative objective function vector;
Figure 891695DEST_PATH_IMAGE094
Figure 912741DEST_PATH_IMAGE095
Figure 770975DEST_PATH_IMAGE096
the method comprises the steps of calculating a part load signal curve of an automobile and obtaining a root mean square error value of the part load signal curve of the automobile through corresponding test;
Figure 789747DEST_PATH_IMAGE097
representing an iteration termination error vector;
Figure 733432DEST_PATH_IMAGE098
Figure 27010DEST_PATH_IMAGE099
Figure 739751DEST_PATH_IMAGE100
representing an iteration end error vector
Figure 726162DEST_PATH_IMAGE097
Each column of data value parameter constants.
S2502, model parameter set
Figure 422722DEST_PATH_IMAGE101
Within the element value range of (2), randomly producing K model parameter sets
Figure 21456DEST_PATH_IMAGE101
Generating collections
Figure 588704DEST_PATH_IMAGE102
Any individual in set SK
Figure 746016DEST_PATH_IMAGE103
Corresponding objective function vector
Figure 664293DEST_PATH_IMAGE104
The expression formula is as follows:
Figure 299674DEST_PATH_IMAGE105
Figure 987007DEST_PATH_IMAGE106
Figure 518483DEST_PATH_IMAGE107
wherein:
Figure 658477DEST_PATH_IMAGE108
representing a set of model parameters
Figure 363128DEST_PATH_IMAGE109
A set of (a);
Figure 403503DEST_PATH_IMAGE110
representing a set of model parameters
Figure 902617DEST_PATH_IMAGE109
The number of (2) is summed up;
Figure 529908DEST_PATH_IMAGE111
representing a model parameter set of a preset automobile multi-body dynamic load calculation model except a tire model;
Figure 507091DEST_PATH_IMAGE112
Figure 169016DEST_PATH_IMAGE113
Figure 104611DEST_PATH_IMAGE114
representing a set of model parameters
Figure 219198DEST_PATH_IMAGE111
Parameter data of (3);
Figure 72DEST_PATH_IMAGE115
representing model parameters of
Figure 17969DEST_PATH_IMAGE111
A time-corresponding objective function vector;
Figure 858886DEST_PATH_IMAGE116
Figure 460769DEST_PATH_IMAGE117
Figure 45334DEST_PATH_IMAGE118
representing model parameters
Figure 416272DEST_PATH_IMAGE111
The first, e, and F columns of values.
Will be assembled
Figure 428090DEST_PATH_IMAGE108
Each model parameter set in (1)
Figure 251690DEST_PATH_IMAGE109
The parameter values are converted from decimal values to binary values.
S2503, sorting the initially generated set SK aiming at the non-dominance solution, and putting the non-dominance solution into an external set SW, wherein the specific steps are as follows:
(1) For individuals within set SK
Figure 639946DEST_PATH_IMAGE119
Wherein
Figure 363926DEST_PATH_IMAGE120
(2) For all other individuals
Figure 546646DEST_PATH_IMAGE121
Figure 123121DEST_PATH_IMAGE122
And is
Figure 49488DEST_PATH_IMAGE123
Comparing individuals
Figure 129440DEST_PATH_IMAGE124
And individuals
Figure 748640DEST_PATH_IMAGE121
Dominant and non-dominant relationships between;
(3) If none exists
Figure 546832DEST_PATH_IMAGE121
Corresponding objective function vector
Figure 214573DEST_PATH_IMAGE125
Each column of values of (a) is less than unity
Figure 414611DEST_PATH_IMAGE124
Corresponding objective function vector
Figure 706177DEST_PATH_IMAGE126
Each column value of (1) is marked
Figure 726085DEST_PATH_IMAGE124
Is a non-dominant solution; if an individual is present
Figure 259835DEST_PATH_IMAGE121
Corresponding objective function vector
Figure 314379DEST_PATH_IMAGE127
Each column of values of (1) is less than unity
Figure 806540DEST_PATH_IMAGE124
Corresponding objective function vector
Figure 579324DEST_PATH_IMAGE128
Each column value of (1) is recorded synchronously
Figure 651185DEST_PATH_IMAGE124
Is administered once.
(4) Order to
Figure 58770DEST_PATH_IMAGE129
And (3) turning to the step (2) until all non-dominated solutions are found, and putting the non-dominated solutions into the outer set SW.
S2504, setting a selection algorithm for the set SK, specifically including the following steps:
defining any individual in set SK
Figure 190674DEST_PATH_IMAGE124
Fitness function of
Figure 185175DEST_PATH_IMAGE130
The expression formula is as follows:
Figure 60727DEST_PATH_IMAGE131
wherein:
Figure 89863DEST_PATH_IMAGE132
representing an individual
Figure 127089DEST_PATH_IMAGE124
A fitness function of;
Figure 608886DEST_PATH_IMAGE133
means when the individual is
Figure 22550DEST_PATH_IMAGE124
In the case of a non-dominant solution,
Figure 407657DEST_PATH_IMAGE134
(ii) a When the individual is
Figure 615785DEST_PATH_IMAGE124
To govern the solution, take
Figure 584878DEST_PATH_IMAGE133
= number of times dominated +1;
Figure 802233DEST_PATH_IMAGE135
represents a positive constant parameter; e denotes a natural constant.
Selection of individuals in set SK by roulette
Figure 540381DEST_PATH_IMAGE124
Individual, individual
Figure 184989DEST_PATH_IMAGE124
Probability of being selected
Figure 22406DEST_PATH_IMAGE136
The expression formula is as follows:
Figure 512293DEST_PATH_IMAGE137
wherein:
Figure 370527DEST_PATH_IMAGE138
representing an individual
Figure 186037DEST_PATH_IMAGE124
Probability of being selected;
Figure 864143DEST_PATH_IMAGE139
representing an individual
Figure 423300DEST_PATH_IMAGE124
A fitness function of;
Figure 434244DEST_PATH_IMAGE140
the fitness function representing all individuals in the set SK totals up.
S2505, setting a cross algorithm aiming at the set SK: aiming at the individuals who have finished the selection algorithm in the previous step, randomly selecting the individuals to participate in the cross operation,random probability set as
Figure 420654DEST_PATH_IMAGE141
The expression formula is as follows:
Figure 851636DEST_PATH_IMAGE142
wherein:
Figure 948905DEST_PATH_IMAGE143
Figure 516152DEST_PATH_IMAGE144
respectively represent the number of times of the iteration
Figure 407885DEST_PATH_IMAGE145
And t is the random probability of participating in the cross operation
Figure 60583DEST_PATH_IMAGE146
Figure 961543DEST_PATH_IMAGE147
(ii) a t represents the current iteration number; t represents the maximum number of iterations;
Figure 412991DEST_PATH_IMAGE148
represents a constant parameter; random probability of participating in cross operation in first iteration
Figure 741204DEST_PATH_IMAGE149
Is a constant set initially;
the calculation rule of the cross operation is as follows: individuals
Figure 881198DEST_PATH_IMAGE150
And individual
Figure 320270DEST_PATH_IMAGE151
Are all randomly selected individuals, and are aimed at
Figure 862109DEST_PATH_IMAGE150
Corresponding binary data and individuals
Figure 361224DEST_PATH_IMAGE151
The corresponding binary data are exchanged at random points to generate new individuals
Figure 722935DEST_PATH_IMAGE150
And new individuals
Figure 965698DEST_PATH_IMAGE151
S2506, setting a variation algorithm for the set SK, specifically including the following steps: aiming at the individuals who have finished the cross algorithm in the previous step, randomly selecting the individuals to participate in the cross operation, and setting the random probability as
Figure 129088DEST_PATH_IMAGE152
The expression formula is as follows:
Figure 595841DEST_PATH_IMAGE153
wherein:
Figure 444849DEST_PATH_IMAGE154
Figure 225723DEST_PATH_IMAGE155
respectively represent the number of times of the iteration
Figure 742155DEST_PATH_IMAGE156
T is a random probability of participating in mutation operation
Figure 583072DEST_PATH_IMAGE157
Figure 919375DEST_PATH_IMAGE158
(ii) a t represents the current iteration number; t represents the maximum number of iterations; random participation in variation operation in first iterationProbability of
Figure 2476DEST_PATH_IMAGE159
Is a constant set initially;
the calculation rule of the variation operation is as follows: individuals
Figure 373414DEST_PATH_IMAGE160
For randomly selected individuals, for the individuals
Figure 181970DEST_PATH_IMAGE160
The corresponding binary data is inverted at random point positions (namely 0 is changed to 1,1 is changed to 0), so that new individuals are generated
Figure 271149DEST_PATH_IMAGE160
S2507, after the set SK completes the selection algorithm, the cross algorithm and the mutation algorithm, the set SK is updated to be a new set
Figure 393826DEST_PATH_IMAGE161
(ii) a Set after updating
Figure 884850DEST_PATH_IMAGE161
And external set
Figure 67569DEST_PATH_IMAGE162
Taking a union set, then carrying out non-dominated sorting on the union set, and updating to generate a new external set
Figure 879930DEST_PATH_IMAGE162
S2508, if the new external set is generated by updating
Figure 806298DEST_PATH_IMAGE162
If the iteration termination condition is satisfied, the external set is output
Figure 886249DEST_PATH_IMAGE162
Namely, the calibrated parameters of the whole vehicle model except the tire model are obtained; such asGenerating a new external set of fruit updates
Figure 239870DEST_PATH_IMAGE162
If the iteration termination condition is not met, updating to generate a new external set
Figure 38062DEST_PATH_IMAGE162
As a new set SK, the newly produced set SK is calculated again according to the steps from S2504 to S2507 until the maximum iteration number reaches the set maximum iteration number
Figure 768120DEST_PATH_IMAGE163
Or an iteration termination condition is satisfied.
In summary, the method for calibrating the automotive dynamics model provided by the embodiment of the application evaluates the actual measurement load and the simulation load of the component of the vehicle through the root mean square error, evaluates the sensitivity of the parameter of the component model under the condition that the error requirement is not met, and calibrates the model parameter with higher sensitivity by adopting the genetic algorithm, so that the workload of the calibration work of the automotive multi-body dynamics dynamic load calculation model can be greatly reduced.
In some examples, the obtaining of the wheel center load actual measurement data and the wheel center load simulation data based on the test yard characteristic road surface includes:
testing and obtaining the wheel center load actual measurement data on the characteristic road surface of the test yard based on the target vehicle;
and carrying out simulation calculation on the preset automobile dynamics calculation model completing the first calibration on the test yard characteristic road surface to obtain the wheel center load simulation data.
Illustratively, the measured wheel center load data is obtained by testing a target vehicle on a test yard characteristic road surface, the test yard characteristic road side may include belgium, twisted, pebble, washboard, stone, noise, resonance, and the like, and the measured wheel center load data may include: the front wheel comprises a left front wheel center X-direction force/Y-direction force/Z-direction force/X-direction moment/Y-direction moment/Z-direction moment, a right front wheel center X-direction force/Y-direction force/Z-direction force/X-direction moment/Y-direction moment/Z-direction moment, a left rear wheel center X-direction force/Y-direction force/Z-direction moment/Y-direction moment/Z-direction moment, and a right rear wheel center X-direction force/Y-direction force/Z-direction force/X-direction moment/Y-direction moment/Z-direction moment.
The wheel center load simulation data is obtained by performing simulation calculation on 3D digital data corresponding to the test yard characteristic road surface by using a first calibrated preset automobile dynamics calculation model, and the wheel center load simulation data comprises the following steps: the front wheel comprises a left front wheel center X-direction force/Y-direction force/Z-direction force/X-direction moment/Y-direction moment/Z-direction moment, a right front wheel center X-direction force/Y-direction force/Z-direction force/X-direction moment/Y-direction moment/Z-direction moment, a left rear wheel center X-direction force/Y-direction force/Z-direction moment/Y-direction moment/Z-direction moment, and a right rear wheel center X-direction force/Y-direction force/Z-direction force/X-direction moment/Y-direction moment/Z-direction moment.
In summary, according to the method for calibrating the automobile dynamics model provided by the embodiment of the application, the second calibration is performed on the tire model of the preset automobile dynamics calculation model after the first calibration is completed, and the target automobile dynamics calculation model after the second calibration is completed is higher in precision.
In some examples, the second calibrating the tire model of the preset vehicle dynamics calculation model with the first calibration completed according to the wheel center load measured data and the wheel center load simulation data to obtain the target vehicle dynamics model includes:
acquiring a second load root mean square error of the wheel center load actual measurement data and the wheel center load simulation data;
and under the condition that the second load root mean square error is larger than a second preset root mean square error, performing second calibration on the tire model of the preset automobile dynamics calculation model after the first calibration is completed to obtain the second calibrated automobile dynamics calculation model.
Illustratively, the wheel center measured load signals of the automobile are respectively obtained through testing
Figure 968157DEST_PATH_IMAGE164
And calculating to obtain a wheel center simulation load signal set of the automobile
Figure 492680DEST_PATH_IMAGE165
For comparison, the expression formula is as follows:
Figure 73440DEST_PATH_IMAGE166
wherein:
Figure 607190DEST_PATH_IMAGE167
Figure 661734DEST_PATH_IMAGE168
respectively representing forces and moments in three directions of a wheel center X, Y and Z of each wheel, which are obtained by solving, calculating and simulating and obtained by test actual measurement;
Figure 91578DEST_PATH_IMAGE169
representing the curve of the actual measurement load signal of the wheel center of the automobile obtained by solving and calculating
Figure 864362DEST_PATH_IMAGE170
Obtaining the wheel center simulation load signal curve of the automobile by corresponding test
Figure 936223DEST_PATH_IMAGE171
A second loading root mean square error function of;
Figure 845273DEST_PATH_IMAGE172
representing a second predetermined root mean square error. If the root mean square error of the second load meets the conditions, the calculation accuracy of the load signal meets the requirements, the preset automobile dynamics simulation model does not need to be calibrated secondarily, and if the root mean square error of the second load meets the conditions, the preset automobile dynamics simulation model does not need to be calibrated secondarily
Figure 275380DEST_PATH_IMAGE173
If the calculation accuracy of the load signal does not meet the requirement of the preset accuracy, the preset automobile dynamics simulation model needs to be calibrated for the second time.
In summary, the method for calibrating the automobile dynamic model provided by the embodiment of the application judges the actual tire load and the simulated tire load of the vehicle through the root-mean-square error, and does not perform the second calibration under the condition of meeting the error requirement, so that the workload of the calibration work of the automobile multi-body dynamic load calculation model can be greatly reduced.
In some examples, the second calibrating the tire model of the preset vehicle dynamics calculation model with the first calibration completed to obtain the second calibrated vehicle dynamics calculation model includes:
acquiring tire parameter sensitivity;
and carrying out second calibration on the tire parameters with the tire parameter sensitivity greater than a second preset sensitivity by adopting a genetic algorithm to obtain a target automobile dynamics calculation model.
For example, if the tested load signal and the calculated load signal do not satisfy the above condition, the calculation accuracy of the load signal is not satisfied; synchronously marking and reordering the load signals which do not meet the conditions, and defining the load signals as the load signals for subsequently calibrating tire model parameters so as to improve the calculation accuracy, wherein the expression formula is as follows:
Figure 473143DEST_PATH_IMAGE174
Figure 145433DEST_PATH_IMAGE175
wherein:
Figure 908989DEST_PATH_IMAGE176
representing a wheel center load signal set of the automobile obtained by calculation needing precision improvement;
Figure 834964DEST_PATH_IMAGE177
to represent
Figure 316761DEST_PATH_IMAGE053
Correspondingly testing to obtain a wheel center load signal set of the automobile;
Figure 730424DEST_PATH_IMAGE178
Figure 348488DEST_PATH_IMAGE179
Figure 822194DEST_PATH_IMAGE180
to represent
Figure 588025DEST_PATH_IMAGE181
Calculating in the set to obtain a wheel center load signal of the automobile;
Figure 306844DEST_PATH_IMAGE182
Figure 779414DEST_PATH_IMAGE183
Figure 220760DEST_PATH_IMAGE184
to represent
Figure 677149DEST_PATH_IMAGE185
The tests in the set obtain the wheel load signals of the cars,
Figure 432615DEST_PATH_IMAGE186
representing the total number of load signals.
S310, screening tire model parameters: carrying out detailed review on the test result of the tire model parameters; for some model parameters, because the reliability of the test result obtained by the test is high, the model parameters do not participate in calibration optimization in the following process, for example: geometric dimension, mass, moment of inertia and other model parameters; the sensitivity can be calculated for other model parameters, and the expression formula is as follows:
Figure 290850DEST_PATH_IMAGE187
wherein:
Figure 840780DEST_PATH_IMAGE188
expressing an absolute value function;
Figure 17421DEST_PATH_IMAGE189
representing a value of sensitivity of a certain tire model parameter to a certain wheel center load signal of the automobile;
Figure 576578DEST_PATH_IMAGE190
Figure 289319DEST_PATH_IMAGE191
Figure 72468DEST_PATH_IMAGE192
a larger value, a smaller value, and a middle value representing the setting of a certain tire model parameter;
Figure 237870DEST_PATH_IMAGE193
Figure 335139DEST_PATH_IMAGE194
Figure 167966DEST_PATH_IMAGE195
representing a certain tire model parameter as
Figure 826742DEST_PATH_IMAGE190
Figure 479441DEST_PATH_IMAGE191
Figure 114821DEST_PATH_IMAGE192
And obtaining the wheel center load signal of the automobile through the following calculation.
S320, if the sensitivity of a certain tire model parameter to a load signal needing to be improved in calculation precision is smaller than a threshold constant
Figure 802154DEST_PATH_IMAGE196
If the parameter is not sensitive to the load signal, the parameter does not participate in the subsequent calibration meterAnd (4) calculating.
S330, if the sensitivity of the model parameters to the load signals needing to be improved in calculation precision is larger than a threshold constant
Figure 927105DEST_PATH_IMAGE196
Then the parameter is shown to be sensitive to the load signal, and if the sensitivity of the model parameter to the load signal that has achieved computational accuracy is less than a threshold constant
Figure 270362DEST_PATH_IMAGE197
Then the model parameters need to participate in the subsequent calibration calculation.
S340, after the model parameter screening is completed, defining parameter set combination of the tire model
Figure 506171DEST_PATH_IMAGE198
The expression formula is as follows:
Figure 546546DEST_PATH_IMAGE199
wherein the content of the first and second substances,
Figure 780081DEST_PATH_IMAGE198
a set of parameters representing a tire model;
Figure 407372DEST_PATH_IMAGE200
Figure 650134DEST_PATH_IMAGE201
Figure 843218DEST_PATH_IMAGE202
parameters representing a tire model;
Figure 247655DEST_PATH_IMAGE203
a number of parameters representing a tire model;
and for parameter sets to
Figure 362241DEST_PATH_IMAGE198
Sets a value range, which is expressed as follows:
Figure 706897DEST_PATH_IMAGE204
wherein:
Figure 957750DEST_PATH_IMAGE205
a parameter representing a tire model;
Figure 533088DEST_PATH_IMAGE206
Figure 134970DEST_PATH_IMAGE207
and expressing the minimum constant and the maximum constant of the parameter, and taking the value as the maximum value if the data is greater than the maximum value in subsequent iterative computation, and taking the value as the minimum value if the data is less than the minimum value.
S360: and carrying out second calibration by adopting a genetic algorithm to calibrate the model parameters of the tire.
Referring to step S250, calibrating the same algorithm in the preset vehicle dynamics calculation model by using a genetic algorithm may specifically include:
s3601, setting the maximum iteration times, the iteration target function vector and the iteration termination error vector of the subsequent iteration.
S3602, tire model parameter set
Figure 719536DEST_PATH_IMAGE208
Within the range of the element values of (A), randomly producing
Figure 887212DEST_PATH_IMAGE209
Set of tire model parameters
Figure 131986DEST_PATH_IMAGE208
Generating collections
Figure 221165DEST_PATH_IMAGE210
S3603 aiming at initially generated set
Figure 609421DEST_PATH_IMAGE210
Ordering against non-dominant solutions and putting non-dominant solutions into an outer set
Figure 569287DEST_PATH_IMAGE211
In (1).
S3604, setting a target set
Figure 814323DEST_PATH_IMAGE210
Selection algorithm, crossover algorithm, and mutation algorithm.
S3605, in the set
Figure 125219DEST_PATH_IMAGE210
After the selection algorithm, the cross algorithm and the variation algorithm are completed, the algorithm is updated into a new set
Figure 786007DEST_PATH_IMAGE210
(ii) a Will be provided with
Figure 633003DEST_PATH_IMAGE210
And external set
Figure 783361DEST_PATH_IMAGE211
Taking a union set, and then carrying out non-dominated sorting on the union set to generate a new external set
Figure 784815DEST_PATH_IMAGE211
S3606, if the generated external set
Figure 514874DEST_PATH_IMAGE211
If the iteration termination condition is satisfied, the external set is output
Figure 714911DEST_PATH_IMAGE211
The tire model parameters are calibrated; if the generated external set
Figure 505013DEST_PATH_IMAGE211
If the iteration termination condition is not met, updating to generate a new external set
Figure 790500DEST_PATH_IMAGE211
As a new set
Figure 58671DEST_PATH_IMAGE212
Set to be updated
Figure 369608DEST_PATH_IMAGE212
The operation is carried out again according to the steps S3603 to S3605 until the maximum iteration number reaches the set maximum iteration number
Figure 330611DEST_PATH_IMAGE213
Or an iteration termination condition is satisfied.
In summary, the method for calibrating the automobile dynamic model provided by the embodiment of the application evaluates the actual tire measurement load and the tire simulation load of the vehicle through the root mean square error, evaluates the sensitivity of the tire model parameters under the condition that the error requirements are not met, and calibrates the model parameters with higher sensitivity by adopting a genetic algorithm, so that the workload of the calibration work of the automobile multi-body dynamic load calculation model can be greatly reduced.
Referring to fig. 2, an embodiment of a calibration apparatus for an automotive dynamics calculation model in an embodiment of the present application may include:
a first acquisition unit 21 configured to acquire a component actual measurement load and a component simulation load based on a target load applied from a wheel center;
a first calibration unit 22, configured to perform a first calibration on a preset vehicle dynamics calculation model according to the component actual measurement load and the component simulation load to obtain a first calibrated vehicle dynamics calculation model, where the first calibration is to calibrate the preset vehicle dynamics calculation model except for the tire model;
the second obtaining unit 23 is configured to obtain wheel center load actual measurement data and wheel center load simulation data based on the test yard characteristic road surface;
and the second calibration unit 24 is configured to perform a second calibration on the tire model of the preset vehicle dynamics calculation model after the first calibration is completed according to the wheel center load actual measurement data and the wheel center load simulation data, so as to obtain a target vehicle dynamics model.
As shown in fig. 3, an electronic device 300 is further provided in the present embodiment, which includes a memory 310, a processor 320, and a computer program 311 stored in the memory 320 and operable on the processor, and when the computer program 311 is executed by the processor 320, the steps of any one of the above-mentioned calibration methods for a vehicle dynamics calculation model are implemented.
Since the electronic device described in this embodiment is a device used for implementing a calibration apparatus of a vehicle dynamics calculation model in this embodiment, based on the method described in this embodiment, a person skilled in the art can understand a specific implementation manner of the electronic device of this embodiment and various variations thereof, so that how to implement the method in this embodiment by the electronic device is not described in detail herein, and as long as the person skilled in the art implements the device used for implementing the method in this embodiment, the scope of protection intended by this application is included.
In a specific implementation, the computer program 311 may implement any of the embodiments corresponding to fig. 3 when executed by a processor.
It should be noted that, in the foregoing embodiments, the descriptions of the respective embodiments have respective emphasis, and reference may be made to relevant descriptions of other embodiments for parts that are not described in detail in a certain embodiment.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded computer, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Embodiments of the present application further provide a computer program product, where the computer program product includes computer software instructions, and when the computer software instructions are executed on a processing device, the processing device executes a flow of a calibration method of an automobile dynamics calculation model in the embodiment corresponding to fig. 1.
The computer program product includes one or more computer instructions. The procedures or functions according to the embodiments of the present application are all or partially generated when the computer program instructions are loaded and executed on a computer. The computer may be a general purpose computer, a special purpose computer, a network of computers, or other programmable device. The computer instructions may be stored on a computer readable storage medium or transmitted from one computer readable storage medium to another computer readable storage medium, for example, the computer instructions may be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via wired (e.g., coaxial cable, fiber optic, digital Subscriber Line (DSL)) or wireless (e.g., infrared, wireless, microwave, etc.) means. Computer-readable storage media can be any available media that a computer can store or a data storage device, such as a server, data center, etc., that includes one or more available media. The usable medium may be a magnetic medium (e.g., floppy disk, hard disk, magnetic tape), an optical medium (e.g., DVD), or a semiconductor medium (e.g., solid State Disk (SSD)), among others.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the above-described systems, apparatuses and units may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other manners. For example, the above-described apparatus embodiments are merely illustrative, and for example, a division of a unit is merely a logical division, and an actual implementation may have another division, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may be in an electrical, mechanical or other form.
The units described as separate parts may or may not be physically separate, and parts displayed as units may or may not be physical units, may be located in one position, or may be distributed on a plurality of network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of the embodiment.
In addition, functional units in the embodiments of the present application may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit. The integrated unit can be realized in a form of hardware, and can also be realized in a form of a software functional unit.
The integrated unit, if implemented in the form of a software functional unit and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present application may be substantially implemented or contributed to by the prior art, or all or part of the technical solution may be embodied in a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method of the embodiments of the present application. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, an optical disk, or other various media capable of storing program codes.
The above embodiments are only used to illustrate the technical solutions of the present application, and not to limit the same; although the present application has been described in detail with reference to the foregoing embodiments, it should be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; and such modifications or substitutions do not depart from the spirit and scope of the corresponding technical solutions in the embodiments of the present application.

Claims (10)

1. A calibration method of an automobile dynamics calculation model is characterized by comprising the following steps:
acquiring a component actual measurement load and a component simulation load based on a target load applied by a wheel center;
performing first calibration on a preset automobile dynamics calculation model according to the component actual measurement load and the component simulation load to obtain a first calibrated automobile dynamics calculation model, wherein the first calibration is used for calibrating the preset automobile dynamics calculation model except the tire model;
acquiring wheel center load actual measurement data and wheel center load simulation data based on the test yard characteristic road surface;
and carrying out second calibration on the tire model of the preset automobile dynamics calculation model after the first calibration is finished according to the wheel center load actual measurement data and the wheel center load simulation data so as to obtain a target automobile dynamics model.
2. The method of claim 1, wherein obtaining component measured loads and component simulated loads based on the target load applied by the wheel center comprises:
acquiring the actual measurement load of the part measured by the automobile shaft coupling test bench test based on the target load applied to the wheel center;
and acquiring the part simulation load corresponding to the preset automobile dynamics calculation model under the condition that the target load is applied to the wheel center.
3. The method of claim 1, wherein the first calibrating the pre-set vehicle dynamics calculation model based on the component measured loads and the component simulated loads to obtain a first calibrated vehicle dynamics calculation model comprises:
acquiring a first load root mean square error of the component actual measurement load and the component simulation load;
and under the condition that the first load root mean square error is larger than a first preset root mean square error, performing first calibration on the preset automobile dynamics calculation model to obtain a first calibrated automobile dynamics calculation model.
4. The method of claim 3, wherein the first calibrating the predetermined vehicle dynamics calculation model to obtain a first calibrated vehicle dynamics calculation model comprises:
obtaining model parameter sensitivity of the preset automobile dynamics calculation model;
and performing first calibration on the model parameters with the sensitivity greater than a first preset sensitivity by adopting a genetic algorithm to obtain a first calibrated automobile dynamics calculation model.
5. The method of claim 1, wherein the obtaining of wheel center load measured data and wheel center load simulation data based on the test yard characteristic road surface comprises:
testing and acquiring the wheel center load actual measurement data on the characteristic road surface of the test yard based on a target vehicle;
and carrying out simulation calculation on the preset automobile dynamics calculation model completing the first calibration on the test yard characteristic road surface to obtain the wheel center load simulation data.
6. The method according to claim 1, wherein the second calibrating the tire model of the pre-set vehicle dynamics calculation model completing the first calibration according to the measured wheel center load data and the wheel center load simulation data to obtain the target vehicle dynamics model comprises:
acquiring a second load root mean square error of the wheel center load actual measurement data and the wheel center load simulation data;
and under the condition that the second load root mean square error is larger than a second preset root mean square error, performing second calibration on the tire model of the preset automobile dynamics calculation model after the first calibration is completed to obtain a second calibrated automobile dynamics calculation model.
7. The method according to claim 6, wherein the second calibrating the tire model of the preset vehicle dynamics calculation model after the first calibrating to obtain the second calibrated vehicle dynamics calculation model comprises:
acquiring tire parameter sensitivity;
and performing second calibration on the tire parameters with the tire parameter sensitivity greater than a second preset sensitivity by adopting a genetic algorithm to obtain a target automobile dynamics calculation model.
8. A calibration device for an automotive dynamics calculation model is characterized by comprising:
a first acquisition unit configured to acquire a component actual measurement load and a component simulation load based on a target load applied by a wheel center;
the first calibration unit is used for carrying out first calibration on a preset automobile dynamics calculation model according to the component actual measurement load and the component simulation load so as to obtain a first calibrated automobile dynamics calculation model, wherein the first calibration is used for calibrating the preset automobile dynamics calculation model except the tire model;
the second acquisition unit is used for acquiring wheel center load actual measurement data and wheel center load simulation data based on the test yard characteristic road surface;
and the second calibration unit is used for carrying out second calibration on the tire model of the preset automobile dynamics calculation model after the first calibration is finished according to the wheel center load actual measurement data and the wheel center load simulation data so as to obtain a target automobile dynamics model.
9. An electronic device, comprising: memory, processor and computer program stored in the memory and executable on the processor, characterized in that the processor is adapted to carry out the steps of the method for calibrating a computational model of vehicle dynamics according to any one of claims 1 to 7 when executing the computer program stored in the memory.
10. A computer-readable storage medium having stored thereon a computer program, characterized in that: the computer program, when executed by a processor, implements a method for calibrating a computational model of vehicle dynamics according to any one of claims 1-7.
CN202210808083.1A 2022-07-11 2022-07-11 Calibration method of automobile dynamics calculation model and related equipment Active CN114861335B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210808083.1A CN114861335B (en) 2022-07-11 2022-07-11 Calibration method of automobile dynamics calculation model and related equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210808083.1A CN114861335B (en) 2022-07-11 2022-07-11 Calibration method of automobile dynamics calculation model and related equipment

Publications (2)

Publication Number Publication Date
CN114861335A CN114861335A (en) 2022-08-05
CN114861335B true CN114861335B (en) 2022-10-21

Family

ID=82625767

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210808083.1A Active CN114861335B (en) 2022-07-11 2022-07-11 Calibration method of automobile dynamics calculation model and related equipment

Country Status (1)

Country Link
CN (1) CN114861335B (en)

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19910967C1 (en) * 1999-03-12 2000-09-21 Avl Deutschland Gmbh Method for simulating the behavior of a vehicle on a road
EP2729927A4 (en) * 2011-10-06 2014-07-30 Cae Inc Method of developing a mathematical model of dynamics of a vehicle for use in a computer-controlled vehicle simulator
CN107991103A (en) * 2017-10-20 2018-05-04 开沃新能源汽车集团有限公司 A kind of batteries of electric automobile pack arrangement Prediction method for fatigue life based on true road spectrum
CN112329133B (en) * 2020-10-20 2023-03-24 东风汽车集团有限公司 Suspension dynamics model K & C performance calibration method
CN113010964B (en) * 2021-03-16 2024-04-16 慧勒智行汽车技术(昆山)有限公司 Virtual test field-based vehicle bench test load spectrum analysis method
CN113434953B (en) * 2021-06-07 2022-10-28 江铃汽车股份有限公司 Method for correcting whole vehicle attitude of multi-body dynamic model of vehicle
CN113656943B (en) * 2021-07-15 2023-10-31 桂林电子科技大学 Method for extracting fatigue load spectrum of whole chassis part of commercial vehicle
CN114462188A (en) * 2021-12-22 2022-05-10 北京新能源汽车股份有限公司 Road load testing method and device, electronic equipment and storage medium

Also Published As

Publication number Publication date
CN114861335A (en) 2022-08-05

Similar Documents

Publication Publication Date Title
CN113010964B (en) Virtual test field-based vehicle bench test load spectrum analysis method
CN106897527B (en) Method and device for analyzing endurance load of vehicle suspension rack
CN107292013B (en) Method and device for testing strength of suspension system
CN113570057B (en) Vehicle wheel center vertical displacement measuring method and device based on model training
CN114676598B (en) Acceleration method and device for whole vehicle road endurance test of vehicle body system
CN109791094A (en) Method and system for the identification of efficient load
CN112131672A (en) Durable load spectrum simulation method, device, storage medium and device
CN203881541U (en) Measuring device
CN114065373A (en) Automobile control arm rack endurance test method, device and equipment
CN103822789A (en) Method and system for measuring wheel center six-component force
CN115169167A (en) Method and system for optimizing and matching motion stroke parameters of automobile plate spring
CN114861335B (en) Calibration method of automobile dynamics calculation model and related equipment
JP4094885B2 (en) Vehicle simulation method
Nam et al. Durability prediction for automobile aluminum front subframe using nonlinear models in virtual test simulations
CN111090959B (en) Vehicle load spectrum acquisition method and system
CN113607392B (en) Spring arm endurance test method and device
CN114676648B (en) Vehicle load spectrum prediction method and device based on machine learning
CN113656994B (en) Suspension force acquisition method and device for automobile suspension system
Kanchwala et al. Model Building, Hardpoint Optimization & Experimental Correlation of a Single Seater EV-Toyota COMS
Lugo et al. Test-driven full vehicle modelling for ADAS algorithm development
CN111985044B (en) Analysis method and device for rigidity of transverse stabilizer bar
Hatekar et al. Silent block bush design and optimization for pick-up truck leaf spring
Lefèvre et al. Quasi-static and dynamic suspension measurements vs. multi-body and real‑time simulation results
Kersten et al. Modern chassis development as a result of skilfully combining testing and simulation
Ghasemiazar et al. Sensitivity analysis and optimization of an off road car vibration performance using DOE and RSM methods

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