CN105157942A - Implementation method of vibration table integrated with shock excitation module and virtual prototype - Google Patents

Implementation method of vibration table integrated with shock excitation module and virtual prototype Download PDF

Info

Publication number
CN105157942A
CN105157942A CN201510528489.4A CN201510528489A CN105157942A CN 105157942 A CN105157942 A CN 105157942A CN 201510528489 A CN201510528489 A CN 201510528489A CN 105157942 A CN105157942 A CN 105157942A
Authority
CN
China
Prior art keywords
road
frequency
vibration table
suspension
virtual prototype
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510528489.4A
Other languages
Chinese (zh)
Other versions
CN105157942B (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.)
LIUZHOU XINGTA INDUSTRY AND TRADE Co Ltd
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201510528489.4A priority Critical patent/CN105157942B/en
Publication of CN105157942A publication Critical patent/CN105157942A/en
Application granted granted Critical
Publication of CN105157942B publication Critical patent/CN105157942B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Vehicle Body Suspensions (AREA)

Abstract

The invention discloses an implementation method of a vibration table integrated with a shock excitation module and a virtual prototype. Through the MATLAB algorithm-based programming process, a transfer function/user-defined function extracted via the interface of the virtual prototype is adopted as excitation transformation characteristics (provided with a dedicated database). After that, the transfer function/user-defined function of a road excitation spectrum extracted via the virtual prototype and conveyed from a vehicle wheel to a connection point of a vehicle body and a vehicle suspension is adopted as an excitation conversion function. In this way, the road excitation is converted into the excitation at the connection point of the vehicle body and the vehicle suspension. Therefore, the excitation test for vehicles on a vibration table is enabled without any extra suspension system.

Description

The implementation method of the vibration table of a kind of integrated exciting module and virtual prototype
Technical field
The present invention relates to car consumption credit field, be specifically related to the implementation method of the vibration table of a kind of integrated exciting module and virtual prototype.
Background technology
In Automobile Reliability Test, vehicle road test can reflect test findings really, but takes time and effort, and economy is not high.In order to the update accelerating automobile product is simultaneously cost-saving, increasing Automobile Enterprises and automotive research institute/institute select to carry out indoor bench test.Indoor bench test mainly utilizes road simulator stand to reproduce equivalent road surface thus carries out the fail-test of car load or parts, and it can lifting test efficiency, and better repeatable.
But existing vibration table only has hardware facility, there is no corresponding software support, and original vibration table function is fairly simple, discrete exciting module integration is not got up, can not the Vibration Condition of simulated automotive on road preferably, cause the result of testing acquisition on this testing table to there is very large error, do not meet actual conditions.The more important thing is, current vibration table exciting object is the tie point of suspension and vehicle body or the lower swing arm of suspension/subframe, and its pumping signal not acts on the spectrum of road surface roughness information on tire, so need to carry out conversion processing to it.
Existing vibration table function is fairly simple, general only have sinusoidal excitation and the arbitrary excitation simulating signal as road spectrum, discrete exciting module integration is not got up, that is be separate between each exciting module, do not consider the coherence between four-wheel and late effect, therefore existing testing table can not maximize the vibration shape situation of reduction vehicle complete vehicle actual road test, cause the result of testing acquisition on this testing table and actual value to there is very large error, do not meet actual conditions.
Existing similar testing table can not form data channel with extracting directly vehicle suspension transport function with virtual prototype, and directly can not be imported by open interface by user, cause the equivalent spectrum of road surface roughness directly acting on suspension and vehicle body tie point place in bench test to encourage and cannot realize
Such as Dong Yuan is bright, Ma Fangwu, the patent of invention " a kind of auto test bed (CN103063447A) " of the people such as Yang Guobin, its hardware device such as frock and fixture etc. mainly for car test bench are set forth, and software are not embedded vibration table; Dong Yuan is bright simultaneously, Ma Fangwu, another patent of invention " road simulator stand (CN103512757A) " of the people such as Yang Guobin, its using car load as experimental subjects, and road spectrum incentive action is simulated vehicle road test in wheel, although the method can carry out bench test preferably, it does not consider that the excitation phase between coherent relationships between wheel and front and back wheel is poor, and wheel is introduced Bench Test System and add experiment work amount.
Summary of the invention
For solving the problem, the invention provides the implementation method of the vibration table of a kind of integrated exciting module and virtual prototype, by MATLAB arithmetic programming, devise and a kind of extract transport function or User-Defined Functions as excitation conversion characteristic (being provided with private database) by virtual prototype interface, thus reach the transport function (or by User-Defined Functions) of suspension and vehicle body tie point as encouraging transfer functions by extracting road spectrum excitation from virtual prototype by wheel, so just road excitation can be converted into the excitation at suspension and vehicle body tie point place, thus on vibration table, vehicle body is encouraged, and do not need additional suspension system.
For achieving the above object, the technical scheme that the present invention takes is:
An implementation method for the vibration table of integrated exciting module and virtual prototype, comprises the steps:
The generation of S1, spectrum of road surface roughness, under the prerequisite considering four-wheel coherence, generates spectrum of road surface roughness Model in Time Domain;
S11, definition algorithm parameter;
S12, based on the near front wheel, generate standard spectrum of road surface roughness according to four-wheel coherence;
S13, generation Stochastic Road Irregularities: utilize MATLAB to programme, white Gaussian noise signal is generated satisfactory random signal by wave filter, carry out the arbitrary excitation input of simulated roadway with this.
S2, by MATLAB arithmetic programming, designing tyre, to the transport function of suspension, processing testing the data of suspension transport function recorded, drawing the displacement signal after transport function respectively, and being delivered on suspension by displacement signal;
S3, according to the displacement of test determination and the corresponding relation of voltage, the displacement signal of spectrum of road surface roughness is converted into voltage signal.
Wherein, described algorithm parameter comprises
Spatial frequency is worth nl up and down, nh: the upper and lower bound of effective band will within the spectral limit of state specified standards road nl=0.011m -1nh=2.83m -1;
Reference frequency: national regulations value n0=0.1m -1;
Road roughness coefficient Gq (n0): corresponding with the standard spectrum of road surface roughness grade of national Specification, selects the road roughness coefficient that different brackets is corresponding as required;
Vehicle velocity V: speed of a motor vehicle 10-50m/s commonly used by automobile;
Wheelbase L1, wheelspan B: determine according to test vehicle;
The upper lower limit value fl of frequency, fh: according to formula f=n*v, as v=10m/s, the frequency range calculated is: 0.11 ~ 28.3Hz, as v=50m/s, the frequency range calculated is: 0.55 ~ 141.5Hz, only needs the natural frequency of covering car based on signal, and the frequency range of the number of winning the confidence is: 0.55 ~ 28.3Hz;
Sample frequency fs: according to sampling thheorem, sample frequency should be greater than the twice f of original signal maximum frequency s> 2f h=2*28.3=56.6Hz, considers again resolution and the algorithm speed of equipment, gets fs=512Hz;
Sampling number N.
Wherein, described transport function is tried to achieve according to formula H=|mag|*exp (i*pha), and in formula, mag is amplitude, and pha is phase place.
Wherein, the data that amplitude (mag) and phase place (pha) are extremely put to suspension by the wheel that test records calculate gained.
Wherein, spectrum of road surface roughness is spatial frequency spectrum, i.e. the relation of road roughness and road surface spatial frequency, can be converted to the relation of road roughness with road surface length, considers the operation of automobile, can transfer the relation between road roughness and time further to.The smooth-going performance evaluation of vehicle is mainly carried out in time domain, therefore the input on road surface also should represent according to pavement grade time domain data corresponding in standard, and road surface Model in Time Domain is the mathematical model of the relation by certain method for numerical simulation, arbitrary excitation being expressed as road roughness and time
The present invention has following beneficial effect:
For vehicle pavement simulation excitation (four-way) testing table, generate road roughness simulation excitation signal based on GB/T7031-2005 and be converted into voltage signal; Road roughness simulation excitation signal is converted to the excitation of vehicle body tie point via suspension transport function; Suspension transport function is extracted by the virtual prototype of corresponding vehicle, and design specialized data-interface sets up the data channel between testing table and virtual prototype for this reason; Suspension transport function is also imported by User Defined by open interface; Set up vehicle suspension transfer function data storehouse simultaneously, and by the transport function extracting or import call in order to testing table stored in this database; When the coherent relationships considered between four-wheel and the excitation between suspension and wheel transform, make automobile can complete experimental test on vibration table, and suspension system comparatively serious for nonlinear characteristic need not be introduced testing table, ensure that higher test accuracy and test reproduction degree.
Accompanying drawing explanation
Fig. 1 is the time domain unevenness of B level spectrum of road surface roughness in the embodiment of the present invention.
Fig. 2 is B level Road Surface Power Spectrum Density in the embodiment of the present invention.
Fig. 3 is the drive singal figure that in the embodiment of the present invention, four-wheel is relevant.
Embodiment
In order to make objects and advantages of the present invention clearly understand, below in conjunction with embodiment, the present invention is further elaborated.Should be appreciated that specific embodiment described herein only in order to explain the present invention, be not intended to limit the present invention.
As shown in Figure 1, embodiments provide the implementation method of the vibration table of a kind of integrated exciting module and virtual prototype, comprise the steps:
The generation of S1, spectrum of road surface roughness, under the prerequisite considering four-wheel coherence, generates spectrum of road surface roughness Model in Time Domain;
S11, definition algorithm parameter;
S12, based on the near front wheel, generate standard spectrum of road surface roughness according to four-wheel coherence;
S13, generation Stochastic Road Irregularities: utilize MATLAB to programme, white Gaussian noise signal is generated satisfactory random signal by wave filter, carry out the arbitrary excitation input of simulated roadway with this.
S2, by MATLAB arithmetic programming, designing tyre, to the transport function of suspension, processing testing the data of suspension transport function recorded, drawing the displacement signal after transport function respectively, and being delivered on suspension by displacement signal;
S3, according to the displacement of test determination and the corresponding relation of voltage, the displacement signal of spectrum of road surface roughness is converted into voltage signal.
Described algorithm parameter comprises
Spatial frequency is worth nl up and down, nh: the upper and lower bound of effective band will within the spectral limit of state specified standards road nl=0.011m -1nh=2.83m -1;
Reference frequency: national regulations value n0=0.1m -1;
Road roughness coefficient Gq (n0): corresponding with the standard spectrum of road surface roughness grade of national Specification, selects the road roughness coefficient that different brackets is corresponding as required;
Vehicle velocity V: speed of a motor vehicle 10-50m/s commonly used by automobile;
Wheelbase L1, wheelspan B: determine according to test vehicle;
The upper lower limit value fl of frequency, fh: according to formula f=n*v, as v=10m/s, the frequency range calculated is: 0.11 ~ 28.3Hz, as v=50m/s, the frequency range calculated is: 0.55 ~ 141.5Hz, only needs the natural frequency of covering car based on signal, and the frequency range of the number of winning the confidence is: 0.55 ~ 28.3Hz;
Sample frequency fs: according to sampling thheorem, sample frequency should be greater than the twice f of original signal maximum frequency s> 2f h=2*28.3=56.6Hz, considers again resolution and the algorithm speed of equipment, gets fs=512Hz; Sampling number N.
Described transport function is tried to achieve according to formula H=|mag|*exp (i*pha), and in formula, mag is amplitude, and pha is phase place.The data that amplitude (mag) and phase place (pha) are extremely put to suspension by the wheel that test records calculate gained.
Embodiment
B grade standard spectrum of road surface roughness is example:
1, input parameter: emulation speed of a motor vehicle u=50km/h; Wheelspan B=1.6m; Front and back wheel wheelbase is 2.5m; Sampling T.T. is 10s, adjustment parameter P=100
2, the relevant standard spectrum of road surface roughness (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 into voltage signal and drive singal (as Fig. 3):
4, this signal is loaded into the Road Simulation Test that automobile suspended point can realize automobile by testing table
The above is only the preferred embodiment of the present invention; it should be pointed out that for those skilled in the art, under the premise without departing from the principles of the invention; can also make some improvements and modifications, these improvements and modifications also should be considered as protection scope of the present invention.

Claims (4)

1. an implementation method for the vibration table of integrated exciting module and virtual prototype, is characterized in that, comprise the steps:
The generation of S1, spectrum of road surface roughness, under the prerequisite considering four-wheel coherence, generates spectrum of road surface roughness Model in Time Domain;
S11, definition algorithm parameter;
S12, based on the near front wheel, generate standard spectrum of road surface roughness according to four-wheel coherence;
S13, generation Stochastic Road Irregularities: utilize MATLAB to programme, white Gaussian noise signal is generated satisfactory random signal by wave filter, carry out the arbitrary excitation input of simulated roadway with this.
S2, by MATLAB arithmetic programming, designing tyre, to the transport function of suspension, processing testing the data of suspension transport function recorded, drawing the displacement signal after transport function respectively, and being delivered on suspension by displacement signal;
S3, according to the displacement of test determination and the corresponding relation of voltage, the displacement signal of spectrum of road surface roughness is converted into voltage signal.
2. the implementation method of the vibration table of a kind of integrated exciting module according to claim 1 and virtual prototype, it is characterized in that, described algorithm parameter comprises
Spatial frequency is worth n1 up and down, nh: the upper and lower bound of effective band will within the spectral limit of state specified standards road n1=0.011m -1nh=2.83m -1;
Reference frequency: national regulations value n0=0.1m -1;
Road roughness coefficient Gq (n0): corresponding with the standard spectrum of road surface roughness grade of national Specification, selects the road roughness coefficient that different brackets is corresponding as required;
Vehicle velocity V: speed of a motor vehicle 10-50m/s commonly used by automobile;
Wheelbase L1, wheelspan B: determine according to test vehicle;
The upper lower limit value f1 of frequency, fh: according to formula f=n*v, as v=10m/s, the frequency range calculated is: 0.11 ~ 28.3Hz, as v=50m/s, the frequency range calculated is: 0.55 ~ 141.5Hz, only needs the natural frequency of covering car based on signal, and the frequency range of the number of winning the confidence is: 0.55 ~ 28.3Hz;
Sample frequency fs: according to sampling thheorem, sample frequency should be greater than the twice f of original signal maximum frequency s> 2f h=2*28.3=56.6Hz, considers again resolution and the algorithm speed of equipment, gets fs=512Hz; Sampling number N.
3. the implementation method of the vibration table of a kind of integrated exciting module according to claim 1 and virtual prototype, it is characterized in that, described transport function is tried to achieve, in formula according to formula H=|mag|*exp (i*pha), mag is amplitude, and pha is phase place.
4. the implementation method of the vibration table of a kind of integrated exciting module according to claim 3 and virtual prototype, it is characterized in that, the data that amplitude (mag) and phase place (pha) are extremely put to suspension by the wheel that test records calculate gained.
CN201510528489.4A 2015-08-21 2015-08-21 A kind of implementation method of the vibration mechine of integrated exciting module and virtual prototype Active CN105157942B (en)

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 true CN105157942A (en) 2015-12-16
CN105157942B 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)

Cited By (4)

* Cited by examiner, † Cited by third party
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
CN112528419A (en) * 2020-12-24 2021-03-19 上海理工大学 Target structural member and external excitation load database construction method based on correlation
CN113239558A (en) * 2021-05-21 2021-08-10 中国工程物理研究院总体工程研究所 Mechanism and data combined driving transportation vibration modeling method
CN114199593A (en) * 2021-12-16 2022-03-18 东风汽车集团股份有限公司 Device and method for measuring vibration sensitivity of tire rotation excitation to attachment point

Citations (10)

* Cited by examiner, † Cited by third party
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

Patent Citations (10)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
田力军: "道路模拟试验台及其控制策略研究", 《流体传动与控制》 *

Cited By (6)

* Cited by examiner, † Cited by third party
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
CN112528419A (en) * 2020-12-24 2021-03-19 上海理工大学 Target structural member and external excitation load database construction method based on correlation
CN112528419B (en) * 2020-12-24 2022-06-17 上海理工大学 Target structural member and external excitation load database construction method based on correlation
CN113239558A (en) * 2021-05-21 2021-08-10 中国工程物理研究院总体工程研究所 Mechanism and data combined driving transportation vibration modeling method
CN114199593A (en) * 2021-12-16 2022-03-18 东风汽车集团股份有限公司 Device and method for measuring vibration sensitivity of tire rotation excitation to attachment point
CN114199593B (en) * 2021-12-16 2023-09-05 东风汽车集团股份有限公司 Device and method for measuring vibration sensitivity from rotation excitation of tire to attachment point

Also Published As

Publication number Publication date
CN105157942B (en) 2018-10-16

Similar Documents

Publication Publication Date Title
CN105157942A (en) Implementation method of vibration table integrated with shock excitation module and virtual prototype
CN103318180B (en) A kind of vehicle road out-of-flat automatic recognition system and method
Yao et al. Comparison of electric vehicle’s energy consumption factors for different road types
Kamble et al. Development of real-world driving cycle: Case study of Pune, India
Fotouhi et al. Tehran driving cycle development using the k-means clustering method
Günther et al. Big data driven dynamic driving cycle development for busses in urban public transportation
CN104634443B (en) The processing method of HVDC transmission line corona audible noise under high ambient noise
CN112487622B (en) Method and device for locating and sizing electric vehicle charging pile and terminal equipment
Fiori et al. Microscopic series plug-in hybrid electric vehicle energy consumption model: Model development and validation
CN102890866B (en) Traffic flow speed estimation method based on multi-core support vector regression machine
CN102074043B (en) Three-dimensional road surface elevation data generation method for vehicle dynamics research
CN107832910B (en) Method for evaluating influence of road traffic characteristics on carbon monoxide emission concentration
CN110399880A (en) Recognition methods, device and the equipment of a kind of characters on license plate and license plate classification
WO2018232991A1 (en) Selected parking space based parking charging management method and system
CN108593314B (en) Parameter calculation method for test bed of vehicle suspension device
Duarte et al. A methodology to estimate real-world vehicle fuel use and emissions based on certification cycle data
Schipper et al. Dazzled by diesel? The impact on carbon dioxide emissions of the shift to diesels in Europe through 2009
CN104239658A (en) Inverse solution method for nonlinear stiffness characteristic parameters and curve of suspension of air spring seat
CN107123173A (en) For the method and apparatus for the transport condition for monitoring unmanned vehicle
CN102890750A (en) Data analysis method for transportation safety recorder
CN109342073A (en) Acquisition methods, device and the realization device of road excitation load
Alvarez et al. Effect of hybrid system battery performance on determining CO2 emissions of hybrid electric vehicles in real-world conditions
CN203283230U (en) Automatic road unevenness identification system for vehicle
CN115496275A (en) Urban traffic carbon emission prediction method, device, electronic equipment and storage medium
CN102521430A (en) Joint simulation method of electric automobile

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