CN105157942B - A kind of implementation method of the vibration mechine of integrated exciting module and virtual prototype - Google Patents
A kind of implementation method of the vibration mechine of integrated exciting module and virtual prototype Download PDFInfo
- Publication number
- CN105157942B CN105157942B CN201510528489.4A CN201510528489A CN105157942B CN 105157942 B CN105157942 B CN 105157942B CN 201510528489 A CN201510528489 A CN 201510528489A CN 105157942 B CN105157942 B CN 105157942B
- Authority
- CN
- China
- Prior art keywords
- road
- suspension
- frequency
- virtual prototype
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Abstract
The invention discloses the implementation methods of a kind of integrated exciting module and the vibration mechine of virtual prototype, pass through MATLAB arithmetic programmings, it devises one kind and excitation conversion characteristic (being equipped with private database) is used as by virtual prototype interface extraction transmission function or User-Defined Functions, to encourage the transmission function (or by User-Defined Functions) for reaching suspension and vehicle body tie point by wheel to be used as excitation conversion function by extracting road spectrum from virtual prototype, can thus road excitation be converted to suspension and the excitation at vehicle body tie point, to on vibration mechine to vehicle body into row energization, and suspension system need not be added.
Description
Technical field
The present invention relates to car consumption credit fields, and in particular to a kind of vibration examination of integrated exciting module and virtual prototype
Test the implementation method of platform.
Background technology
In Automobile Reliability Test, vehicle road test can really reflect test result, but take time and effort, economical
Property is not high.In order to which the update for accelerating automobile product is cost-effective simultaneously, more and more Automobile Enterprises and automotive research institute/
Institute's selection carries out indoor bench test.Indoor bench test mainly using road simulator stand reproduce equivalent road surface into
The reliability test of row vehicle or parts, it can be with lifting test efficiency, and repeatability is preferably.
But existing vibration mechine only has hardware facility, without corresponding software support, and original vibration tries
It tests that platform function is fairly simple, does not integrate discrete exciting module, it cannot preferable simulated automotive shaking on road
Emotionally condition causes the result for testing acquisition on the testing stand there are prodigious error, does not meet actual conditions.It is prior
It is that current vibration mechine exciting object is suspension and the tie point of vehicle body or the lower swing arm of suspension/subframe, excitation
Signal is not the road surface spectrum information acted on tire, so needing to carry out conversion processing to it.
Existing vibration mechine function is fairly simple, general only to have sinusoidal excitation with arbitrary excitation as road spectrum
Analog signal does not integrate discrete exciting module, that is to say, that is independent from each other, does not have between each exciting module
In view of between four-wheel coherence with delay effect, therefore existing testing stand cannot maximize reduction vehicle complete vehicle road road test
The vibration shape situation tested causes the result for testing acquisition on the testing stand, there are prodigious error, not meet reality with actual value
Situation.
Existing similar testing stand cannot with virtual prototype formed data channel directly to extract vehicle suspension transmission function,
And it can not be introduced directly by open interface by user, cause to directly act on suspension and vehicle body tie point in bench test
The equivalent spectrum of road surface roughness excitation at place cannot achieve
Such as Dong Yuan is bright, Ma Fangwu, a kind of patent of invention " auto test bed of Yang Guobin et al.
(CN103063447A) " it, is illustrated mainly for the hardware device of car test bench such as tooling and fixture etc., will not
Software is embedded in vibration mechine；Dong Yuanming simultaneously, another patent of invention " Road Simulation Test of Ma Fangwu, Yang Guobin et al.
Road spectrum incentive action using vehicle as experimental subjects, and is simulated vehicle road by platform (CN103512757A) " in wheel
Experiment, although this method can preferably carry out bench test, it does not account for coherent relationships between wheel and preceding
Excitation phase between trailing wheel is poor, and wheel introducing Bench Test System is increased test work load.
Invention content
To solve the above problems, the present invention provides the realities of a kind of integrated exciting module and the vibration mechine of virtual prototype
Existing method is devised a kind of by virtual prototype interface extraction transmission function or User Defined by MATLAB arithmetic programmings
Function is as excitation conversion characteristic (being equipped with private database), to by extracting road spectrum excitation from virtual prototype by wheel
Suspension is reached with the transmission function (or by User-Defined Functions) at vehicle body tie point as excitation conversion function, thus
The excitation that road excitation can be converted at suspension and vehicle body tie point, on vibration mechine to vehicle body into row energization,
And suspension system need not be added.
To achieve the above object, the technical solution that the present invention takes is：
A kind of implementation method of the vibration mechine of integrated exciting module and virtual prototype, includes the following steps：
The generation of S1, spectrum of road surface roughness generate spectrum of road surface roughness Model in Time Domain under the premise of considering four-wheel coherence；
S11, algorithm parameter is defined；
S12, based on the near front wheel, according to four-wheel coherence generate standard spectrum of road surface roughness；
S13, Stochastic Road Irregularities are generated：It is programmed using MATLAB, white Gaussian noise signal is met by filter generation
It is required that random signal, with this come simulated roadway arbitrary excitation input.
S2, pass through MATLAB arithmetic programmings, the transmission function of designing tyre to suspension, the suspension transmission letter that experiment is measured
Several data are handled, and obtain the displacement signal after transmission function respectively, and displacement signal is transmitted on suspension；
S3, according to experiment measure displacement and voltage correspondence, by the displacement signal of spectrum of road surface roughness be converted into voltage letter
Number.
Wherein, the algorithm parameter includes
Value n1, nh above and below spatial frequency：The upper and lower bound of effective band will state specified standards road spectral limit it
Interior n1=0.011m-1, nh=2.83m-1；
Reference frequency：National regulations value n0=0.1m-1；
Road roughness coefficient Gq (n0)：It is corresponding with the standard spectrum of road surface roughness grade of national Specification, it selects as needed
Select the corresponding road roughness coefficient of different brackets；
Vehicle velocity V：Automobile often uses speed 10-50m/s；
Wheelbase L1, wheelspan B：It is determined according to experiment vehicle；
The upper lower limit value fl, fh of frequency：According to formula f=n*v, as v=10m/s, calculated frequency range is：
0.11~28.3Hz, as v=50m/s, calculated frequency range is：0.55~141.5Hz only needs to cover based on signal
The frequency range of the intrinsic frequency of vehicle, the number of winning the confidence is：0.55~28.3Hz；
Sample frequency fs：According to sampling thheorem, sample frequency should be greater than twice of f of original signal maximum frequencys＞ 2fh=2*
28.3=56.6Hz, it is also contemplated that the resolution ratio and algorithm speed of equipment, take fs=512Hz；
Sampling number N.
Wherein, the transmission function is according to formula H=| mag | * exp (i*pha) are acquired, and in formula, mag is amplitude, pha
For phase.
Wherein, amplitude (mag) and phase (pha) calculate institute by testing data of the wheel measured to suspension up to point
.
Wherein, spectrum of road surface roughness is spatial frequency spectrum, i.e. the relationship of road roughness and road surface spatial frequency, can be converted to road
Face unevenness considers the operation of automobile with the relationship of road surface length, can further switch between road roughness and time
Relationship.The smooth-going performance evaluation of vehicle mainly carries out in the time domain, therefore the input on road surface also should be according to corresponding in standard
Pavement grade indicates that road surface Model in Time Domain is to be expressed as arbitrary excitation by certain method for numerical simulation with time domain data
The mathematical model of the relationship of road roughness and time
The invention has the advantages that：
For vehicle pavement simulation excitation (four-way) testing stand, it is imitative that road roughness is generated based on GB/T 7031-2005
True pumping signal is simultaneously converted into voltage signal；Road roughness simulation excitation signal is converted into vehicle body via suspension transmission function
Tie point encourages；Suspension transmission function is extracted in the virtual prototype by corresponding to vehicle, and design specialized data-interface establishes examination thus
Test the data channel between platform and virtual prototype；Suspension transmission function can also be imported by open interface by User Defined；Together
Shi Jianli vehicle suspension transfer function datas library, and the transmission function extracted or imported is stored in the database in case testing stand
It calls；Between considering four-wheel coherent relationships and suspension and wheel between excitation conversion in the case of, so that automobile is existed
Experimental test is completed on vibration mechine, and without the more serious suspension system of nonlinear characteristic is introduced testing stand, is ensured
Higher test accuracy and experiment reproduction degree.
Description of the drawings
Fig. 1 is the time domain unevenness of B grades of spectrum of road surface roughness in the embodiment of the present invention.
Fig. 2 is B grades of Road Surface Power Spectrum Densities in the embodiment of the present invention.
Fig. 3 is the drive signal figure that four-wheel is relevant in the embodiment of the present invention.
Specific implementation mode
In order to make objects and advantages of the present invention be more clearly understood, the present invention is carried out with reference to embodiments further
It is described in detail.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, it is not used to limit this hair
It is bright.
As shown in Figure 1, an embodiment of the present invention provides the vibration mechines of a kind of integrated exciting module and virtual prototype
Implementation method includes the following steps：
The generation of S1, spectrum of road surface roughness generate spectrum of road surface roughness Model in Time Domain under the premise of considering four-wheel coherence；
S11, algorithm parameter is defined；
S12, based on the near front wheel, according to four-wheel coherence generate standard spectrum of road surface roughness；
S13, Stochastic Road Irregularities are generated：It is programmed using MATLAB, white Gaussian noise signal is met by filter generation
It is required that random signal, with this come simulated roadway arbitrary excitation input.
S2, pass through MATLAB arithmetic programmings, the transmission function of designing tyre to suspension, the suspension transmission letter that experiment is measured
Several data are handled, and obtain the displacement signal after transmission function respectively, and displacement signal is transmitted on suspension；
S3, according to experiment measure displacement and voltage correspondence, by the displacement signal of spectrum of road surface roughness be converted into voltage letter
Number.
The algorithm parameter includes
Value n1, nh above and below spatial frequency：The upper and lower bound of effective band will state specified standards road spectral limit it
Interior n1=0.011m-1, nh=2.83m-1；
Reference frequency：National regulations value n0=0.1m-1；
Road roughness coefficient Gq (n0)：It is corresponding with the standard spectrum of road surface roughness grade of national Specification, it selects as needed
Select the corresponding road roughness coefficient of different brackets；
Vehicle velocity V：Automobile often uses speed 10-50m/s；
Wheelbase L1, wheelspan B：It is determined according to experiment vehicle；
The upper lower limit value fl, fh of frequency：According to formula f=n*v, as v=10m/s, calculated frequency range is：
0.11~28.3Hz, as v=50m/s, calculated frequency range is：0.55~141.5Hz only needs to cover based on signal
The frequency range of the intrinsic frequency of vehicle, the number of winning the confidence is：0.55~28.3Hz；
Sample frequency fs：According to sampling thheorem, sample frequency should be greater than twice of f of original signal maximum frequencys＞ 2fh=2*
28.3=56.6Hz, it is also contemplated that the resolution ratio and algorithm speed of equipment, take fs=512Hz；Sampling number N.
The transmission function is according to formula H=| mag | * exp (i*pha) are acquired, and in formula, mag is amplitude, and pha is phase
Position.Amplitude (mag) and phase (pha) are as testing obtained by data calculating of the wheel measured to suspension up to point.
Embodiment
For B grades of standard spectrum of road surface roughness：
1, input parameter：Emulate speed u=50km/h；Wheelspan B=1.6m；Front and back wheel base is 2.5m；Sample total time
For 10s, adjusting parameter P=100
2, the relevant standard spectrum of road surface roughness (such as Fig. 1,2) of four-wheel is generated by programming：
3, the spectrum of road surface roughness displacement signal of above-mentioned generation is converted to voltage signal i.e. drive signal (such as Fig. 3)：
4, the signal is loaded into the Road Simulation Test that automobile can be realized on automobile suspended point by testing stand
The above is only a preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art
For member, without departing from the principle of the present invention, it can also make several improvements and retouch, these improvements and modifications are also answered
It is considered as protection scope of the present invention.
Claims (4)
1. a kind of implementation method of the vibration mechine of integrated exciting module and virtual prototype, which is characterized in that including walking as follows
Suddenly：
The generation of S1, spectrum of road surface roughness generate spectrum of road surface roughness Model in Time Domain under the premise of considering four-wheel coherence；
S11, algorithm parameter is defined；
S12, based on the near front wheel, according to four-wheel coherence generate standard spectrum of road surface roughness；
S13, Stochastic Road Irregularities are generated：It is programmed using MATLAB, white Gaussian noise signal is met the requirements by filter generation
Random signal, with this come simulated roadway arbitrary excitation input；
S2, by MATLAB arithmetic programmings, the transmission function of designing tyre to suspension, the suspension transmission function that experiment is measured
Data are handled, and obtain the displacement signal after transmission function respectively, and displacement signal is transmitted on suspension；
S3, according to experiment measure displacement and voltage correspondence, convert the displacement signal of spectrum of road surface roughness to voltage signal.
2. the implementation method of the vibration mechine of a kind of integrated exciting module according to claim 1 and virtual prototype,
It is characterized in that, the algorithm parameter includes
Value n1, nh above and below spatial frequency：The upper and lower bound of effective band will within the spectral limit of state specified standards road n1
=0.011m-1, nh=2.83m-1；
Reference frequency：National regulations value n0=0.1m-1；
Road roughness coefficient Gq (n0)：It is corresponding with the standard spectrum of road surface roughness grade of national Specification, it selects as needed not
The corresponding road roughness coefficient of ad eundem；
Vehicle velocity V：Automobile often uses speed 10-50m/s；
Wheelbase L1, wheelspan B：It is determined according to experiment vehicle；
The upper lower limit value fl, fh of frequency：According to formula f=n*v, as v=10m/s, calculated frequency range is：0.11~
28.3Hz, as v=50m/s, calculated frequency range is：0.55~141.5Hz only needs to cover consolidating for vehicle based on signal
There are the frequency, the frequency range for the number of winning the confidence to be：0.55~28.3Hz；
Sample frequency fs：According to sampling thheorem, sample frequency should be greater than twice of f of original signal maximum frequencys＞ 2fh=2*28.3
=56.6Hz, it is also contemplated that the resolution ratio and algorithm speed of equipment, take fs=512Hz；
Sampling number N.
3. the implementation method of the vibration mechine of a kind of integrated exciting module according to claim 1 and virtual prototype,
It is characterized in that, the transmission function is according to formula H=| mag | * exp (i*pha) are acquired, and in formula, mag is amplitude, and pha is phase
Position.
4. the implementation method of the vibration mechine of a kind of integrated exciting module according to claim 3 and virtual prototype,
It is characterized in that, amplitude (mag) and phase (pha) are as testing obtained by data calculating of the wheel measured to suspension up to point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510528489.4A CN105157942B (en) | 2015-08-21 | 2015-08-21 | A kind of implementation method of the vibration mechine of integrated exciting module and virtual prototype |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510528489.4A CN105157942B (en) | 2015-08-21 | 2015-08-21 | A kind of implementation method of the vibration mechine of integrated exciting module and virtual prototype |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105157942A CN105157942A (en) | 2015-12-16 |
| CN105157942B true CN105157942B (en) | 2018-10-16 |
Family
ID=54798878
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510528489.4A Active CN105157942B (en) | 2015-08-21 | 2015-08-21 | A kind of implementation method of the vibration mechine of integrated exciting module and virtual prototype |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105157942B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109520756A (en) * | 2018-11-20 | 2019-03-26 | 中车齐齐哈尔车辆有限公司 | The test method and device of goods train vehicle safety |
| CN112528419B (en) * | 2020-12-24 | 2022-06-17 | 上海理工大学 | Target structural member and external excitation load database construction method based on correlation |
| CN113239558B (en) * | 2021-05-21 | 2022-12-20 | 中国工程物理研究院总体工程研究所 | Mechanism and data combined driving transportation vibration modeling method |
| CN114199593B (en) * | 2021-12-16 | 2023-09-05 | 东风汽车集团股份有限公司 | Device and method for measuring vibration sensitivity from rotation excitation of tire to attachment point |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61164133A (en) * | 1985-01-16 | 1986-07-24 | Honda Motor Co Ltd | Vibration testing method |
| JP2007218790A (en) * | 2006-02-17 | 2007-08-30 | Railway Technical Res Inst | Rolling-stock-car vibrating system |
| JP2009014370A (en) * | 2007-07-02 | 2009-01-22 | Imv Corp | Vibration testing device for carrying out multiple loop controls |
| CN101510230A (en) * | 2009-03-11 | 2009-08-19 | 同济大学 | Vehicle road load emulation method |
| CN101770533A (en) * | 2008-12-29 | 2010-07-07 | 王萌 | Vehicle-ground interaction dynamic simulation system |
| JP2011169866A (en) * | 2010-02-22 | 2011-09-01 | Shimadzu Corp | Vibration testing device |
| CN102519692A (en) * | 2011-11-28 | 2012-06-27 | 重庆长安汽车股份有限公司 | Rigid-body mode integration test method for automobile power assembly and suspension |
| CN102944430A (en) * | 2012-11-28 | 2013-02-27 | 重庆理工大学 | Motorcycle frame actual running load spectrum fatigue reproduction test method |
| JP2013167557A (en) * | 2012-02-16 | 2013-08-29 | Kayaba System Machinery Kk | Vibration testing machine |
| CN104849039A (en) * | 2014-12-18 | 2015-08-19 | 北汽福田汽车股份有限公司 | Automobile gearbox suspension test method and system |
-
2015
- 2015-08-21 CN CN201510528489.4A patent/CN105157942B/en active Active
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPS61164133A (en) * | 1985-01-16 | 1986-07-24 | Honda Motor Co Ltd | Vibration testing method |
| JP2007218790A (en) * | 2006-02-17 | 2007-08-30 | Railway Technical Res Inst | Rolling-stock-car vibrating system |
| JP2009014370A (en) * | 2007-07-02 | 2009-01-22 | Imv Corp | Vibration testing device for carrying out multiple loop controls |
| CN101770533A (en) * | 2008-12-29 | 2010-07-07 | 王萌 | Vehicle-ground interaction dynamic simulation system |
| CN101510230A (en) * | 2009-03-11 | 2009-08-19 | 同济大学 | Vehicle road load emulation method |
| JP2011169866A (en) * | 2010-02-22 | 2011-09-01 | Shimadzu Corp | Vibration testing device |
| CN102519692A (en) * | 2011-11-28 | 2012-06-27 | 重庆长安汽车股份有限公司 | Rigid-body mode integration test method for automobile power assembly and suspension |
| JP2013167557A (en) * | 2012-02-16 | 2013-08-29 | Kayaba System Machinery Kk | Vibration testing machine |
| CN102944430A (en) * | 2012-11-28 | 2013-02-27 | 重庆理工大学 | Motorcycle frame actual running load spectrum fatigue reproduction test method |
| CN104849039A (en) * | 2014-12-18 | 2015-08-19 | 北汽福田汽车股份有限公司 | Automobile gearbox suspension test method and system |
Non-Patent Citations (1)
| Title |
|---|
| 道路模拟试验台及其控制策略研究;田力军;《流体传动与控制》;20110930(第5(2011)期);第16-19页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105157942A (en) | 2015-12-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105157942B (en) | A kind of implementation method of the vibration mechine of integrated exciting module and virtual prototype | |
| CN101393079B (en) | Automobile vehicle body structure fatigue life predicting system | |
| CN103318180B (en) | A kind of vehicle road out-of-flat automatic recognition system and method | |
| WO2021119947A1 (en) | Method for quick detection of damage to bridge, and related device | |
| CN105005694B (en) | A kind of bridge fatigue life frequency-domain analysis method based on dynamic weighing system | |
| CN113010964B (en) | Virtual test field-based vehicle bench test load spectrum analysis method | |
| CN103411743B (en) | Movable equipment is utilized to carry out the method for testing of portable impact shock to bridge | |
| CN106323451B (en) | The method and device thereof of displacement signal are obtained by acceleration signals | |
| CN101625285A (en) | Virtual car crash detecting system and construction method thereof | |
| CN110231181A (en) | A kind of vehicle physical method for parameter estimation based on vibration-testing information | |
| CN110119577A (en) | A method of extracting pure electric vehicle truck batteries bag bracket structure load spectrum | |
| CN104239658A (en) | Inverse solution method for nonlinear stiffness characteristic parameters and curve of suspension of air spring seat | |
| ATE403894T1 (en) | METHOD FOR OPTIMIZING VEHICLES AND MOTORS FOR DRIVING SUCH VEHICLES | |
| CN110411757A (en) | Spindle nose dynamic load calculation method, device, equipment and storage medium | |
| CN108593314B (en) | Parameter calculation method for test bed of vehicle suspension device | |
| CN102890750B (en) | Data analysis method for transportation safety recorder | |
| Burdzik et al. | Adapted linear decimation procedures for TFR analysis of non-stationary vibration signals of vehicle suspensions | |
| CN203283230U (en) | Automatic road unevenness identification system for vehicle | |
| CN104132792A (en) | Device and method for using laser displacement signals for testing bridge floor flexibility | |
| CN104455157A (en) | Obtaining method of car seat suspension hydraulic buffer nonlinear speed characteristic parameter | |
| CN105973189A (en) | Method for obtaining surface roughness load spectrum of soft road surface | |
| CN106771746B (en) | A kind of electric car dynamic operation condition electromagnetic disturbance fast appraisement method | |
| CN107292034A (en) | The in-car vibration noise simulating analysis of road excitation based on Hyperworks and system | |
| CN103425827B (en) | A kind of simulating analysis of train wheel flat | |
| CN109342078A (en) | The economic testing method, apparatus and test device of electric four-wheel drive vehicle |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20190930 Address after: 545005 the Guangxi Zhuang Autonomous Region District, Liuzhou, Yufeng Bay Road, No. 32-1 Patentee after: LIUZHOU XINGTA INDUSTRY AND TRADE CO., LTD. Address before: 545005 No. 32-1, community Bay Road, Liuzhou, the Guangxi Zhuang Autonomous Region Patentee before: Miao Zhi Gang |
|
| TR01 | Transfer of patent right |