CN106096193A - Car load idling, acceleration simulating engine excitation solving system and method - Google Patents

Abstract

The invention discloses a kind of car load idling, accelerate simulating engine excitation solving system and method, module, computing module, display module and result output module are set including solving；Solve and module is set for arranging electromotor basic parameter, Engine Block Test cylinder pressure and analytical model；Computing module includes time-domain calculation and frequency-domain calculations；Time-domain calculation is based on electromotor basic parameter and Engine Block Test cylinder pressure, and combine inertia force computation model, moment of inertia computation model gentle moment of volume computation model calculates time domain excitation；Frequency-domain calculations is to combine time-domain analysis result and analytical model, by fast fourier transform, time domain excitation is converted into frequency domain and encourages, and then extracts the excitation amplitude that Relative degree frequency is corresponding；And carry out image output by display model, and utilize result output module output drive amplitude and phase place.The present invention can solve excitation for specific idling, accelerating mode, and computational efficiency is high, and demand parameter model is simpler.

Description

Car load idling, acceleration simulating engine excitation solving system and method

Technical field

The invention belongs to vehicle complete vehicle NVH technology, be specifically related to a kind of car load idling, accelerate simulating engine excitation and solve System and method.

Background technology

Car load NVH performance study is one of current automotive community priority research areas, and Engine Excitation and road surface are swashed Encourage and be the main driving source of car load vibration and noise.When running car is in flat road surface or idle running, only electromotor Itself being driving source, the car load NVH problem therefore caused by Engine Excitation always is domestic and international main engine plants and colleges and universities scholar Primary study direction.In electromotor works, accurately say it is that piston is doing reciprocal up and down motion with in combustion process, create attached Afterburning and moment of torsion, they excite automobile chassis vibrate by power assembly suspension system, thus the vibration in generation car load car with Noise.Therefore, can predict and analyze, in automobile research early stage, the car load vibration noise caused by Engine Excitation to seem particularly Important.

At present, domestic main engine plants are mainly by suspending the loading actual measurement vibrational excitation of active/passive side at car load CAE model, Thus analyze car internal vibration noise response.But, the method is the harshest to the demand of excitation, tends not to accomplish at car load item Mesh is carried out early stage and is just predicted by CAE simulation analysis means.Additionally, due to excitation is derived from the active/passive side of electromotor, Can not analyze suspension rigidity and the installation site impact on car load NVH performance, thus cause can not problem analysis generation comprehensively Root.So analyzing car load NVH performance from " source path response ", it is necessary to be analyzed realizing road at car load CAE environment Footpath optimizes, therefore needs an Engine Excitation to solve software platform system and carry out excitation calculating.

It is less that present engine excitation solves software, wherein fairly perfect with AVL Excite function, applies wider.Although AVL Excite can solve accurate Engine Excitation, and it is main or is used for analyzing the vibration NVH of engine crankshaft itself Performance, is rarely used in car load NVH Engine Excitation and calculates.Additionally, AVL Excite need accurate bent axle FEM (finite element) model and Detailed electromotor basic parameter, needs to consume substantial amounts of analysis time；Car load idling, acceleration are analyzed needs and are entered excitation itself The quick order analysis of row, for analyzing the order reason producing car load NVH problem, and AVL Excite is not met by such Demand.And AVL Excite is as a business software platform, is mainly used in engine body vibration analysis, so to whole Engine Excitation solution efficiency needed for car NVH performance evaluation is not high enough.

Therefore, it is necessary to develop a kind of new car load idling, accelerate simulating engine excitation solving system and method.

Summary of the invention

It is an object of the invention to provide a kind of car load idling, accelerate simulating engine excitation solving system and method, with gram Take that existing business software exists efficiency is low, demand parameter model complicated, can not solve sharp for specific idling, accelerating mode The shortcoming encouraged.

Car load idling of the present invention, acceleration simulating engine excitation solving system, including:

Solving and arrange module, for user setup engine parameter and analytical model, this engine parameter includes electromotor Basic parameter and Engine Block Test cylinder pressure；

Computing module, including to the time-domain calculation of Engine Excitation and frequency-domain calculations；Described time-domain calculation is based on described Electromotor basic parameter and Engine Block Test cylinder pressure, and combine inertia force computation model, moment of inertia computation model and gas force Square computation model calculates time domain excitation, i.e. inertia force in time domain, the excitation of moment of inertia gentle moment of volume；Described frequency-domain calculations It is to combine time-domain analysis result and analytical model, by fast fourier transform, time domain excitation is converted into frequency domain and encourages, then Extract the excitation amplitude that Relative degree frequency is corresponding；

Display module, the form graphically changed shows all kinds of excitation results that described computing module is calculated, including time Territory excitation, frequency domain excitation and acceleration excitation order analysis；

Result output module, encourages output display by the frequency domain of electromotor in table form, including excitation amplitude and Phase place.

Described inertia force computation model is:

∑F_ι2=-M_recλrω²[cos2 θ+cos2 (θ+180 °)+cos2 (θ+180 °)+cos2 θ]=-4M_recλrω²cos2θ；

Wherein, F_t2It is two grades of reciprocal inertia forces, M_recRepresenting piston mass and 1/3 connecting rod quality sum, λ represents that crank is long The half of degree and length of connecting rod ratio, ω represents engine speed, and θ is current time crank rotation angle, and r represents crank length Half.

Described moment of inertia computation model is:

∑T_i=2M_recr²ω²sin2θ；

Wherein, T_iRepresent engine inertia moment, M_recRepresenting piston mass and 1/3 connecting rod quality sum, r represents crank The half of length, ω represents engine speed, and θ is current time crank rotation angle.

Described gas moment computation model is:

$T_g=F_{g}r(sin\mathrm{\θ}+cos\mathrm{\θ}\frac{\mathrm{\λ}sin\mathrm{\θ}}{\sqrt{1-{\mathrm{\λ}}^2{\sin }^2\mathrm{\θ}}});$

$F_g=P_g\frac{{\mathrm{\πD}}^2}{4};$

Wherein, T_gRepresent engine air moment of volume, F_gRepresenting the gas pressure to piston, r represents the half of crank length, θ For current time crank rotation angle, λ represents half and length of connecting rod ratio, the P of crank length_gRepresent Engine Block Test cylinder Pressure, D represents piston diameter.

Described solving arranges module, display module and result output module and is respectively provided with user interface.

A kind of car load idling of the present invention, acceleration simulating engine excitation method for solving, comprise the following steps:

(1) setting is solved:

Arranging engine parameter and analytical model, this engine parameter includes electromotor basic parameter and Engine Block Test cylinder Pressure, described electromotor basic parameter includes idling speed, piston diameter, crank length, length of connecting rod, piston mass, connecting rod matter Amount, described analytical model includes idling analysis and accelerates to analyze；

(2) solve:

Read electromotor basic parameter and Engine Block Test cylinder pressure, calculate mould in conjunction with inertia force computation model, moment of inertia Type gentle moment of volume computation model, calculates time domain excitation；Discriminatory analysis pattern, if idling analysis, then by quick Fourier Time domain excitation is converted into frequency domain excitation by leaf transformation, and provides each order of electromotor excitation result；If accelerating to analyze, first Transfer the time domain excitation of each for accelerator rotating speed to frequency domain by fast fourier transform to encourage, obtain each rank of accelerator Secondary excitation result, obtains producing the problem order of acceleration risk, then exports acceleration corresponding under this order by order analysis Frequency domain encourages；

(3) display output:

The form graphically changed shows all kinds of excitation results that described computing module is calculated, including time domain excitation, frequently Territory excitation and acceleration excitation order analysis；

(4) result output:

In table form the frequency domain of electromotor is encouraged output display, including amplitude and the phase place of excitation.

Described inertia force computation model is:

∑F_ι2=-M_recλrω²[cos2 θ+cos2 (θ+180 °)+cos2 (θ+180 °)+cos2 θ]=-4M_recλrω²cos2θ；

Wherein, F_t2It is two grades of reciprocal inertia forces, M_recRepresenting piston mass and 1/3 connecting rod quality sum, λ represents that crank is long The half of degree and length of connecting rod ratio, ω represents engine speed, and θ is current time crank rotation angle, and r represents crank length Half.

Described moment of inertia computation model is:

∑T_i=2M_recr²ω²sin2θ；

Wherein, T_iRepresent engine inertia moment, M_recRepresenting piston mass and 1/3 connecting rod quality sum, r represents crank The half of length, ω represents engine speed, and θ is current time crank rotation angle.

Described gas moment computation model is:

$T_g=F_{g}r(sin\mathrm{\θ}+cos\mathrm{\θ}\frac{\mathrm{\λ}sin\mathrm{\θ}}{\sqrt{1-{\mathrm{\λ}}^2{\sin }^2\mathrm{\θ}}});$

$F_g=P_g\frac{{\mathrm{\πD}}^2}{4};$

Wherein, T_gRepresent engine air moment of volume, F_gRepresenting the gas pressure to piston, r represents the half of crank length, θ For current time crank rotation angle, λ represents half and length of connecting rod ratio, the P of crank length_gRepresent Engine Block Test cylinder Pressure, D represents piston diameter.

The invention have the advantages that the present invention does not relies on existing business software, by integrated engine inertia force, Moment of inertia gentle moment of volume solving model, develops four module, it may be assumed that solve and arrange module, computing module, display module and knot Really output module, carries out rapid solving and analysis to car load idling, accelerating mode Engine Excitation, contributes to car load NVH performance Analyze.Additionally, the present invention has taken into full account that the Engine Excitation under different operating mode produces principle, idling is encouraged each order Excitation superposition excitation solves；Acceleration excitation completes formulate order excitation solve, and realize encouraging order analysis.For the ease of whole Car NVH Performance Analysis, present invention achieves output module and directly exports calculating excitation result, generate Nastran analyzing card Sheet, decreases analysis project's teacher operation element.

Accompanying drawing explanation

Fig. 1 is the system architecture diagram of the present invention；

Fig. 2 is to solve in the present invention to arrange module GUI function distribution schematic diagram；

Fig. 3 is the execution flow process of computing module in the present invention；

Fig. 4 is the execution logic of algoritic module in the present invention；

Fig. 5 is time domain excitation GUI display layout in the present invention；

Fig. 6 is that frequency domain of the present invention encourages GUI display layout；

Fig. 7 is result output module GUI layout in the present invention；

Fig. 8 is instantiation time domain excitation outcome graphic display in the present invention；

Fig. 9 is example frequency domain excitation outcome graphic display in the present invention；

Figure 10 is instantiation Engine Excitation result output in the present invention.

Detailed description of the invention

The invention will be further described below in conjunction with the accompanying drawings.

Car load idling as shown in Figure 1, acceleration simulating engine excitation solving system, arrange module, calculating including solving Module, display module and result output module.Formally begin for carried electromotor and calculate its output drive.Algorithm performs Data stream as shown in Figure 4, for time-domain calculation, read electromotor basic parameter and Engine Block Test cylinder pressure, in conjunction with inertia force Computation model, one cycle intrinsic motivation excitation curve of moment of inertia computation model gentle moment of volume analysis of calculation models, and by Display module exports；For frequency-domain calculations, in conjunction with Engine Excitation result in time domain and analytical model, become by fast Flourier Change, time domain excitation is transformed into frequency domain, then extract the excitation amplitude that Relative degree frequency is corresponding, and carried out by display model Image exports, and utilizes result output module output drive amplitude and phase place.Heretofore described solving arranges module, display Module and result output module are respectively provided with user interface, it is simple to car load NVH simulation analysis engineer's intuitive operation, problem analysis and Extract result；And algoritic module mainly comprises time domain and frequency-domain calculations solves, wherein result in time domain is intuitively analyzed for user and is sent out Excitation on motive, and frequency-domain calculations is analyzed software (such as: Nastran, OptStruct) for structure simulation and is analyzed car load on frequency domain Car internal vibration noise problem.

There is in the user interface of module " engine parameter " button as in figure 2 it is shown, solve to arrange and " solve calculating mould Formula " button, user, by clicking on " engine parameter " button, enters engine parameter and arranges interface, start problem analysis Basic parameter and the Engine Block Test cylinder pressure of machine are configured, such as: idling speed is 750rpm, piston diameter 76mm, crank Length 82.6mm, length of connecting rod 133.2mm, piston mass 242g, connecting rod quality 450g.User " solves calculating mould by clicking on Formula " button enter analytical model interface is set, analytical model include idling analysis and accelerate analyze.

As shown in Figure 3 and Figure 4, computing module includes Engine Excitation (reciprocal inertia force, moment of inertia and gas force Square) time-domain calculation and frequency-domain calculations.This computing module support increase income algorithmic code update, revised excitation can be embedded and solve Method, improves car load NVH simulation analysis precision.Time-domain calculation is to survey based on described electromotor basic parameter, engine test Examination data, and combine embedded inertia force computation model, moment of inertia computation model gentle moment of volume computation model is when calculating Territory is encouraged, i.e. inertia force in time domain, the excitation of moment of inertia gentle moment of volume.Frequency-domain calculations is to combine time-domain analysis result and divide Analysis pattern, it is judged that analytical model, if idling analysis, is then converted into frequency domain by fast fourier transform by time domain excitation and swashs Encourage, and each order of electromotor excitation result, if accelerating to analyze, first pass through fast fourier transform by each for accelerator turn The time domain excitation of speed is converted into frequency domain excitation, thus obtains each order excitation result of accelerator, is obtained by order analysis To producing the problem order accelerating risk, then export acceleration frequency domain excitation corresponding under this order.

As shown in Figure 5 and Figure 6, all kinds of excitations knot that the form display computing module that display module is graphically changed is calculated Really, including time domain excitation, frequency domain excitation and acceleration excitation order analysis.As shown in Figure 8, main in the display of time domain excitation result GUI Gas moment to be included, inertia force, inertia force and phase angle.As it is shown in figure 9, described frequency domain excitation result GUI shows main bag Include gas moment order analysis, moment of inertia order analysis, inertia force order analysis and excitation spectrum figure.Display module main Effect is to provide patterned display interface and is easy to engineer and analyzes Engine Excitation in time domain and frequency domain, it is determined that produces and encourages The frequency range of peak value, it was predicted that the risk of car load NVH problem；Shown by order analysis, the main excitation of directviewing description accelerator Order.

As it is shown in fig. 7, the frequency domain excitation that result output module is mainly the electromotor the most just solved is defeated Go out display, including amplitude and the phase place of excitation.Excitation relates to: gas moment, inertia force and moment of inertia, and the frequency of correspondence Value.This module can show the amplitude of all kinds of excitation of electromotor under each frequency, and the phase place of correspondence, as shown in Figure 10, just Excitation peak point is accurately judged in engineer.Output result can directly be processed by engineer, derives.

A kind of car load idling of the present invention, acceleration simulating engine excitation method for solving, comprise the following steps:

(1) setting is solved:

Arranging engine parameter and analytical model, engine parameter includes electromotor basic parameter and Engine Block Test cylinder Pressure, described electromotor basic parameter includes idling speed, piston diameter, crank length, length of connecting rod, piston mass, connecting rod matter Amount, described analytical model includes idling analysis and accelerates to analyze.

(2) solve:

Read electromotor basic parameter and Engine Block Test cylinder pressure, calculate mould in conjunction with inertia force computation model, moment of inertia Type gentle moment of volume computation model, calculates time domain excitation；Discriminatory analysis pattern, if idling analysis, then by quick Fourier Time domain excitation is converted into frequency domain excitation by leaf transformation, and provides each order of electromotor excitation result；If accelerating to analyze, first Transfer the time domain excitation of each for accelerator rotating speed to frequency domain by fast fourier transform to encourage, obtain each rank of accelerator Secondary excitation result, obtains producing the problem order of acceleration risk, then exports acceleration corresponding under this order by order analysis Frequency domain encourages.Wherein:

Described inertia force computation model is:

∑F_ι2=-M_recλrω²[cos2 θ+cos2 (θ+180 °)+cos2 (θ+180 °)+cos2 θ]=-4M_recλrω²cos2θ (1)；

Wherein, F_t2It is two grades of reciprocal inertia forces, M_recRepresenting piston mass and 1/3 connecting rod quality sum, λ represents that crank is long The half of degree and length of connecting rod ratio, ω represents engine speed, and θ is current time crank rotation angle, and r represents crank length Half.

Described moment of inertia computation model is:

∑T_i=2M_recr²ω²sin2θ (2)；

Wherein, T_iRepresent engine inertia moment, M_recRepresenting piston mass and 1/3 connecting rod quality sum, r represents crank The half of length, ω represents engine speed, and θ is current time crank rotation angle.

Described gas moment computation model is:

$T_g=F_{g}r(sin\mathrm{\θ}+cos\mathrm{\θ}\frac{\mathrm{\λ}sin\mathrm{\θ}}{\sqrt{1-{\mathrm{\λ}}^2{\sin }^2\mathrm{\θ}}})---\left(3\right);$

$F_g=P_g\frac{{\mathrm{\πD}}^2}{4}---\left(4\right);$

Wherein, T_gRepresent engine air moment of volume, F_gRepresenting the gas pressure to piston, r represents the half of crank length, θ For current time crank rotation angle, λ represents half and length of connecting rod ratio, the P of crank length_gRepresent Engine Block Test cylinder Pressure, D represents piston diameter.

(3) display output:

The form graphically changed shows all kinds of excitation results that described computing module is calculated, including time domain excitation, frequently Territory excitation and acceleration excitation order analysis.

(4) result output:

By output module, the frequency domain of electromotor is encouraged output display, including amplitude and the phase place of excitation.

In the present invention, the Gas torque in Fig. 8 to Figure 10 represents that gas moment, Inertia Force represent inertia force, Inertia torque represents that moment of inertia, Frequence represent that frequency, angle/degree represent angle/degree, Excite Value represents excitation value.

Claims (9)

1. a car load idling, acceleration simulating engine excitation solving system, it is characterised in that including:

Solving and arrange module, for user setup engine parameter and analytical model, this engine parameter includes that electromotor is basic Parameter and Engine Block Test cylinder pressure；

Computing module, including to the time-domain calculation of Engine Excitation and frequency-domain calculations；Described time-domain calculation is to start based on described Machine basic parameter and Engine Block Test cylinder pressure, and combine inertia force computation model, moment of inertia computation model gentle moment of volume meter Calculate model and calculate time domain excitation, be i.e. inertia force in time domain, the excitation of moment of inertia gentle moment of volume；Described frequency-domain calculations is knot Close time-domain analysis result and analytical model, by fast fourier transform, time domain excitation is converted into frequency domain and encourages, then extract The excitation amplitude that Relative degree frequency is corresponding；

Display module, the form graphically changed shows all kinds of excitation results that described computing module is calculated, and swashs including time domain Encourage, frequency domain encourages and accelerates excitation order analysis；

Result output module, encourages output display by the frequency domain of electromotor in table form, including amplitude and the phase place of excitation.

Car load idling the most according to claim 1, acceleration simulating engine excitation solving system, it is characterised in that: described Inertia force computation model is:

∑F_ι2=-M_recλrω²[cos2 θ+cos2 (θ+180 °)+cos2 (θ+180 °)+cos2 θ]=-4M_recλrω²cos2θ；

Wherein, F_t2It is two grades of reciprocal inertia forces, M_recRepresenting piston mass and 1/3 connecting rod quality sum, λ represents crank length Half and length of connecting rod ratio, ω represents engine speed, and θ is current time crank rotation angle, and r represents the one of crank length Half.

Car load idling the most according to claim 1, acceleration simulating engine excitation solving system, it is characterised in that: described Moment of inertia computation model is:

∑T_i=2M_recr²ω²sin2θ；

Wherein, T_iRepresent engine inertia moment, M_recRepresenting piston mass and 1/3 connecting rod quality sum, r represents crank length Half, ω represents engine speed, and θ is current time crank rotation angle.

Car load idling the most according to claim 1, acceleration simulating engine excitation solving system, it is characterised in that: described Gas moment computation model is:

$T_g=F_{g}r\left(sin\mathrm{\θ}+cos\mathrm{\θ}\frac{\mathrm{\λ}sin\mathrm{\θ}}{\sqrt{1-{\mathrm{\λ}}^2{\sin }^2\mathrm{\θ}}}\right);$

$F_g=P_g\frac{{\mathrm{\πD}}^2}{4};$

Wherein, T_gRepresent engine air moment of volume, F_gRepresenting the gas pressure to piston, r represents the half of crank length, and θ is for working as Front moment crank rotation angle, λ represents half and length of connecting rod ratio, the P of crank length_gRepresent Engine Block Test cylinder pressure, D table Show piston diameter.

5. according to the arbitrary described car load idling of Claims 1-4, accelerating simulating engine excitation solving system, its feature exists In: described in solve and module, display module and result output module are set are respectively provided with user interface.

6. a car load idling, acceleration simulating engine excitation method for solving, it is characterised in that comprise the following steps:

(1) setting is solved:

Arranging engine parameter and analytical model, this engine parameter includes electromotor basic parameter and Engine Block Test cylinder pressure, Described electromotor basic parameter includes idling speed, piston diameter, crank length, length of connecting rod, piston mass, connecting rod quality, Described analytical model includes idling analysis and accelerates to analyze；

(2) solve:

Read electromotor basic parameter and Engine Block Test cylinder pressure, in conjunction with inertia force computation model, moment of inertia computation model and Gas moment computation model, calculates time domain excitation；Discriminatory analysis pattern, if idling analysis, is then become by fast Flourier Time domain excitation of changing commanders is converted into frequency domain excitation, and provides each order of electromotor excitation result；If accelerating to analyze, first pass through The time domain excitation of each for accelerator rotating speed is transferred to frequency domain excitation by fast fourier transform, and each order obtaining accelerator swashs Encourage result, obtain producing the problem order of acceleration risk by order analysis, then export acceleration frequency domain corresponding under this order Excitation；

(3) display output:

The form graphically changed shows all kinds of excitation results that described computing module is calculated, and swashs including time domain excitation, frequency domain Encourage and accelerate to encourage order analysis；

(4) result output:

In table form the frequency domain of electromotor is encouraged output display, including amplitude and the phase place of excitation.

Car load idling the most according to claim 6, acceleration simulating engine excitation method for solving, it is characterised in that:

Described inertia force computation model is:

∑F_ι2=-M_recλrω²[cos2 θ+cos2 (θ+180 °)+cos2 (θ+180 °)+cos2 θ]=-4M_recλrω²cos2θ；

Wherein, F_t2It is two grades of reciprocal inertia forces, M_recRepresenting piston mass and 1/3 connecting rod quality sum, λ represents crank length Half and length of connecting rod ratio, ω represents engine speed, and θ is current time crank rotation angle, and r represents the one of crank length Half.

8. according to the car load idling described in claim 6 or 7, accelerate simulating engine excitation method for solving, it is characterised in that: institute Stating moment of inertia computation model is:

∑T_i=2M_recr²ω²sin2θ；

Wherein, T_iRepresent engine inertia moment, M_recRepresenting piston mass and 1/3 connecting rod quality sum, r represents crank length Half, ω represents engine speed, and θ is current time crank rotation angle.

9. according to the car load idling described in claim 6 or 7, accelerate simulating engine excitation method for solving, it is characterised in that: institute Stating gas moment computation model is:

$T_g=F_{g}r\left(sin\mathrm{\θ}+cos\mathrm{\θ}\frac{\mathrm{\λ}sin\mathrm{\θ}}{\sqrt{1-{\mathrm{\λ}}^2{\sin }^2\mathrm{\θ}}}\right);$

$F_g=P_g\frac{{\mathrm{\πD}}^2}{4};$

Wherein, T_gRepresent engine air moment of volume, F_gRepresenting the gas pressure to piston, r represents the half of crank length, and θ is for working as Front moment crank rotation angle, λ represents half and length of connecting rod ratio, the P of crank length_gRepresent Engine Block Test cylinder pressure, D table Show piston diameter.