CN110532716A - A kind of Numerical Simulation of Working Process in Gasoline Engine system - Google Patents
A kind of Numerical Simulation of Working Process in Gasoline Engine system Download PDFInfo
- Publication number
- CN110532716A CN110532716A CN201910830136.8A CN201910830136A CN110532716A CN 110532716 A CN110532716 A CN 110532716A CN 201910830136 A CN201910830136 A CN 201910830136A CN 110532716 A CN110532716 A CN 110532716A
- Authority
- CN
- China
- Prior art keywords
- main unit
- exhaust
- gasoline engine
- cylinder
- compression
- 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.)
- Pending
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 title claims abstract description 16
- 230000006835 compression Effects 0.000 claims abstract description 62
- 238000007906 compression Methods 0.000 claims abstract description 62
- 238000002485 combustion reaction Methods 0.000 claims description 13
- 239000000446 fuel Substances 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 7
- 238000006073 displacement reaction Methods 0.000 claims description 2
- 241000208340 Araliaceae Species 0.000 claims 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 claims 1
- 235000003140 Panax quinquefolius Nutrition 0.000 claims 1
- 235000008434 ginseng Nutrition 0.000 claims 1
- 238000009413 insulation Methods 0.000 claims 1
- 238000011161 development Methods 0.000 abstract description 5
- 238000009795 derivation Methods 0.000 description 3
- 230000000977 initiatory effect Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000002459 sustained effect Effects 0.000 description 1
- 238000010977 unit operation Methods 0.000 description 1
Landscapes
- Output Control And Ontrol Of Special Type Engine (AREA)
- Combined Controls Of Internal Combustion Engines (AREA)
Abstract
The invention discloses a kind of Numerical Simulation of Working Process in Gasoline Engine systems, more particularly to simulation system field, including system main program, the system main program includes compression main unit, burn main unit, expand main unit, it is vented main unit, it is folded to open main unit and air inlet main unit, the compression main unit, burn main unit, expand main unit, it is vented main unit, it is folded to open main unit and air inlet main unit is sequentially connected, the system main program is connected with call subroutine, the system main program is also configured with worksheet, the worksheet includes hve_xy table and hvs_xy table.System of the invention is not limited to the limitation of molding software, and suitable for different models, versatility is high, can greatly improve working efficiency, saves the development time.
Description
Technical field
The present invention relates to simulation system fields, it is more particularly related to a kind of Gasoline Engine Working Process Numerical-Mode
Quasi- system.
Background technique
Current China automobile has entered the period that maintains sustained and rapid growth, and the huge market demand is the hair of auto Parts Industry
Exhibition provides huge development space.The development and automobile production and consumption of auto parts industry are closely bound up.It is run according to the world
Standard, automobile industry vehicle and components macro ratio are about 1:1.7, and in prediction on such basis, Chinese auto Parts Industry is still
Have greatly improved space.The purposes of this software is mainly used for engine development initial stage progress Engine Process Simulation, fastly
Speed obtains the performance parameter of engine, preferred to carry out scheme.Its main feature is that arithmetic speed is fast, operator is required it is low, no
As with business software, it is desirable that operator skillfully grasps the application method of software, and modeling and operation time are long, provides to operation
Source requires high, low efficiency, and the engine development period is long.Meanwhile software can be used for the engine of various types, versatility is good.
Summary of the invention
In order to overcome the drawbacks described above of the prior art, the embodiment of the present invention provides a kind of Gasoline Engine Working Process Numerical-Mode
Quasi- system passes through setting compression main unit, burning main unit, expansion main unit, row using the platform based on MATLAB
Gas main unit, it is folded open main unit and air inlet main unit, can the course of work effectively under the various operating conditions of analog gasoline machine, packet
Six integrated stages such as compression, burning, expansion, exhaust, inlet and exhaust valve overlapping and air inlet are included, correlation performance parameters, system are obtained
It is not limited to the limitation of molding software, suitable for different models, versatility is high, can greatly improve working efficiency, saves exploitation
Time.
To achieve the above object, the invention provides the following technical scheme: a kind of Numerical Simulation of Working Process in Gasoline Engine system,
The system is based on MATLAB platform, including system main program, and the system main program includes compression main unit, burning main
Unit, exhaust main unit, is folded and opens main unit and air inlet main unit expansion main unit, the compression main unit, combustion
Burn main unit, expansion main unit, exhaust main unit, it is folded open main unit and air inlet main unit is sequentially connected, it is described
System main program is connected with call subroutine, and the system main program is also configured with worksheet, and the worksheet includes hve_xy
Table and hvs_xy table;
It is provided with gasoline engine structural parameters, each stroke phase angle, original state parameter in the compression main unit and tests
Parameter is demonstrate,proved, inlet valve lift parameter is provided in the hvs_xy table, is provided with exhaust valve lift parameter in the hve_xy table.
In a preferred embodiment, the compression main unit, burning main unit, expansion main unit, exhaust
Main unit, the folded source code file for opening main unit, air inlet main unit and worksheet are located under same English path, described
The system file name of source code file is made of English and number.
In a preferred embodiment, system specifically comprises the following steps: when being simulated
Step 1, when operation, system main program is based on MATLAB platform, MATLAB platform is opened, by MATLAB platform
Operating path is specified on the root of source code file, and modifies inlet valve lift parameter and exhaust valve lift parameter;
Step 2 opens compression main unit, system parameter is arranged, specifically:
%% structural parameters
S=0.082;% piston stroke
Lambda=0.3066566941;% crank to connecting rod length ratio
Epsilon=10.1;% compression ratio
D=0.074;% cylinder diameter
N=3601;% revolving speed
L0=14.8;% chemically correct fuel
Vs=S*pi*D^2/4;% cylinder capacity
Cylinder_number=4;% cylinder number
Ze=2;% exhaust valve number
Zs=2;% inlet valve number
Beta_e=45;% exhaust valve cone angle
Beta_s=45;% inlet valve cone angle
Dve=0.0198;% exhaust valve seat throat diameter
Dvs=0.0243;% suction valve cone throat diameter
Dme=0.0229;% exhaust valve average diameter
Dms=0.0275;% inlet valve average diameter
Doe=0.0055;% exhaust valve rod diameter
Dos=0.0055;% inlet valve shank diameter
Alpha_phi_initial=0.9;% compression stroke starts excess air coefficient in cylinder
Ta=374;% compression stroke starts Temperature of Working in cylinder
Pa=0.132850435031838;% compression stroke starts initial pressure in cylinder
Angle_1=250;
Va=(pi*D^2/4) * [S/ (epsilon-1)+(S/2) * ((1+1/Lambda)-cos (pi*angle_1/
180)-(1/Lambda)*sqrt(1-Lambda^2*(sin(pi*angle_1/180)).^2))];% compression stroke starts cylinder
Interior initialization volume
ML=pa*1e6*va/ (9.81* (29.2647-0.0402/alpha_phi_initial) * ta);% ideal gas
Body state equation pv=mRgT obtains working medium quality in the cylinder that compression stroke starts
Each stroke angles of %%
Angle_1=250;
Angle_2=331;
Angle_3=405;
Angle_4=487.6;
Angle_5=708;
Angle_6=741.6;
Angle_7=970;
Phi_z=angle_3-angle_2;% combustion continuation angle;
%% initiation parameter
Ta=374;% compression stroke starts Temperature of Working in cylinder
Alpha_phi_initial=0.9;% compression stroke starts excess air coefficient in cylinder
Pa=0.132850435031838;% compression stroke starts initial pressure in cylinder
A1=pi*D^2/4;A2=a1;% cylinder cap and piston area
Tw1=580;Tw2=550;Tw3=480;% compression stroke starts cylinder cap, piston, cylinder sleeve initial temperature
Gf=mL/ (1+L0*alpha_phi_initial);% compression stroke starts initial fuel quality in cylinder
Pr=0.1058904;Exhaust pipe initial pressure after % exhaust valve
Ps=0.096264;Air inlet pipe initial pressure after % inlet valve
Ts=320;Air inlet pipe initial pressure after % inlet valve
Tr=1150;Exhaust pipe initial temperature after % exhaust valve
M=1.717;% combustion quality index;
It modifies to the parameter for including among the above;
Step 3 runs MATLAB platform, according to some number of prompt in command Window input 1~4, with determination
Instantaneous average heat transfer coefficient calculation formula, then some number in 1~3 is inputted, with the instantaneous specific heat at constant volume calculation formula of determination;
Step 4, operation burning main unit, display in command Window " * * * burn * * * " indicate operation completion, according to
Secondary opening follow-up expansion main unit, is folded and opens main unit and air inlet main unit exhaust main unit, and operation button is clicked,
" * * * has expanded * * * ", " * * * is pure to be vented * * * ", " * * * valve is folded to end * * * " and " * * * can be shown after the completion of operation respectively
The pure complete * * * " of air inlet, when all system main program operations finish, command Window can show that gasoline engine performance emulates key data
As a result.
In a preferred embodiment, in the step 1, gasoline engine compression top center correspond to time point be 360 °
CA, with 720 °C of A for a working cycles, including four air inlet, compression, expansion and exhaust strokes.
In a preferred embodiment, after the completion of to gasoline engine simulation trial, MATLAB platform is according to simulation
Emulation data draw different data images, including Temperature of Working curve, pressure curve, figure-of-merit curve, excess air system in cylinder
Number curve, cylinder inner volume, working medium specific enthalpy at inlet and exhaust valve, combustion heat release rate, instantaneous average heat transfer coefficient, shows function at piston displacement
Figure, equivalent amount of fuel, the rate of pressure rise, temperature rise rate, adiabatic exponent, gas constant, intake and exhaust rate of discharge, inlet and exhaust valve liter
Journey, air inlet and exhaust valve discharge coefficient, the effective flow area of inlet and exhaust valve, density, specific internal energy, specific enthalpy etc. in instantaneous cylinder.
Technical effect and advantage of the invention:
The present invention is based on MATLAB platform, by setting compression main unit, burning main unit, expansion main unit,
Exhaust main unit, it is folded open main unit and air inlet main unit, can the course of work effectively under the various operating conditions of analog gasoline machine,
Including six integrated stages such as compression, burning, expansion, exhaust, inlet and exhaust valve overlapping and air inlets, correlation performance parameters are obtained, are
System is not limited to the limitation of molding software, and suitable for different models, versatility is high, can greatly improve working efficiency, saving is opened
Send out the time.
Specific embodiment
Below in conjunction with the embodiment of the present invention, technical scheme in the embodiment of the invention is clearly and completely described,
Obviously, described embodiments are only a part of the embodiments of the present invention, instead of all the embodiments.Based in the present invention
Embodiment, every other embodiment obtained by those of ordinary skill in the art without making creative efforts, all
Belong to the scope of protection of the invention.
Embodiment 1
A kind of Numerical Simulation of Working Process in Gasoline Engine system, the system are based on MATLAB platform, including system main program,
The system main program is including compressing main unit, burning main unit, expansion main unit, exhaust main unit, folding and open
Main unit and air inlet main unit, the compression main unit, burning main unit, expansion main unit, exhaust main are mono-
Member, it is folded open main unit and air inlet main unit is sequentially connected, the system main program is connected with call subroutine, the system
Main program is also configured with worksheet, and the worksheet includes hve_xy table and hvs_xy table;
It is provided with gasoline engine structural parameters, each stroke phase angle, original state parameter in the compression main unit and tests
Parameter is demonstrate,proved, inlet valve lift parameter is provided in the hvs_xy table, is provided with exhaust valve lift parameter in the hve_xy table;
The compression main unit, burning main unit, expansion main unit, exhaust main unit, it is folded open main unit,
The source code file of air inlet main unit and worksheet is located under same English path, the system file name of the source code file
It is made of English and number.
Embodiment 2
On the basis of embodiment 1, system of the invention specifically comprises the following steps: when being simulated
Step 1, when operation, system main program is based on MATLAB platform, MATLAB platform is opened, by MATLAB platform
Operating path is specified on the root of source code file, and modifies inlet valve lift parameter and exhaust valve lift parameter;
Step 2 opens compression main unit, system parameter is arranged, specifically:
%% structural parameters
S=0.082;% piston stroke
Lambda=0.3066566941;% crank to connecting rod length ratio
Epsilon=10.1;% compression ratio
D=0.074;% cylinder diameter
N=3601;% revolving speed
L0=14.8;% chemically correct fuel
Vs=S*pi*D^2/4;% cylinder capacity
Cylinder_number=4;% cylinder number
Ze=2;% exhaust valve number
Zs=2;% inlet valve number
Beta_e=45;% exhaust valve cone angle
Beta_s=45;% inlet valve cone angle
Dve=0.0198;% exhaust valve seat throat diameter
Dvs=0.0243;% suction valve cone throat diameter
Dme=0.0229;% exhaust valve average diameter
Dms=0.0275;% inlet valve average diameter
Doe=0.0055;% exhaust valve rod diameter
Dos=0.0055;% inlet valve shank diameter
Alpha_phi_initial=0.9;% compression stroke starts excess air coefficient in cylinder
Ta=374;% compression stroke starts Temperature of Working in cylinder
Pa=0.132850435031838;% compression stroke starts initial pressure in cylinder
Angle_1=250;
Va=(pi*D^2/4) * [S/ (epsilon-1)+(S/2) * ((1+1/Lambda)-cos (pi*angle_1/
180)-(1/Lambda)*sqrt(1-Lambda^2*(sin(pi*angle_1/180)).^2))];% compression stroke starts cylinder
Interior initialization volume
ML=pa*1e6*va/ (9.81* (29.2647-0.0402/alpha_phi_initial) * ta);% ideal gas
Body state equation pv=mRgT obtains working medium quality in the cylinder that compression stroke starts
Each stroke angles of %%
Angle_1=250;
Angle_2=331;
Angle_3=405;
Angle_4=487.6;
Angle_5=708;
Angle_6=741.6;
Angle_7=970;
Phi_z=angle_3-angle_2;% combustion continuation angle;
%% initiation parameter
Ta=374;% compression stroke starts Temperature of Working in cylinder
Alpha_phi_initial=0.9;% compression stroke starts excess air coefficient in cylinder
Pa=0.132850435031838;% compression stroke starts initial pressure in cylinder
A1=pi*D^2/4;A2=a1;% cylinder cap and piston area
Tw1=580;Tw2=550;Tw3=480;% compression stroke starts cylinder cap, piston, cylinder sleeve initial temperature
Gf=mL/ (1+L0*alpha_phi_initial);% compression stroke starts initial fuel quality in cylinder
Pr=0.1058904;Exhaust pipe initial pressure after % exhaust valve
Ps=0.096264;Air inlet pipe initial pressure after % inlet valve
Ts=320;Air inlet pipe initial pressure after % inlet valve
Tr=1150;Exhaust pipe initial temperature after % exhaust valve
M=1.717;% combustion quality index;
It modifies to the parameter for including among the above;
Step 3 runs MATLAB platform, according to some number of prompt in command Window input 1~4, with determination
Instantaneous average heat transfer coefficient calculation formula, then some number in 1~3 is inputted, with the instantaneous specific heat at constant volume calculation formula of determination;
Step 4, operation burning main unit, display in command Window " * * * burn * * * " indicate operation completion, according to
Secondary opening follow-up expansion main unit, is folded and opens main unit and air inlet main unit exhaust main unit, and operation button is clicked,
" * * * has expanded * * * ", " * * * is pure to be vented * * * ", " * * * valve is folded to end * * * " and " * * * can be shown after the completion of operation respectively
The pure complete * * * " of air inlet, when all system main program operations finish, command Window can show that gasoline engine performance emulates key data
As a result;
In above-mentioned steps, it is that a job follows with 720 °C of A that gasoline engine compression top center, which corresponds to time point as 360 °C of A,
Ring, including four air inlet, compression, expansion and exhaust strokes;
Further, after the completion of to gasoline engine simulation trial, MATLAB platform is drawn not according to the emulation data of simulation
Same data image, including Temperature of Working curve, pressure curve, figure-of-merit curve, excess air coefficient curve, piston position in cylinder
Working medium specific enthalpy at shifting, cylinder inner volume, inlet and exhaust valve, combustion heat release rate, instantaneous average heat transfer coefficient, indicator card, equivalent amount of fuel,
The rate of pressure rise, temperature rise rate, adiabatic exponent, gas constant, intake and exhaust rate of discharge, inlet and exhaust valve lift, air inlet and exhaust valve stream
Density, specific internal energy, specific enthalpy etc. in the effective flow area of coefficient of discharge, inlet and exhaust valve, instantaneous cylinder.
Embodiment 3
System of the invention is in the process of running, specific as follows to the modification of each parameter and the relevant source code of unit operation:
%% structural parameters
S=0.082;% piston stroke
Lambda=0.3066566941;% crank to connecting rod length ratio
Epsilon=10.1;% compression ratio
D=0.074;% cylinder diameter
N=3601;% revolving speed
L0=14.8;% chemically correct fuel
Vs=S*pi*D^2/4;% cylinder capacity
Cylinder_number=4;% cylinder number
Ze=2;% exhaust valve number
Zs=2;% inlet valve number
Beta_e=45;% exhaust valve cone angle
Beta_s=45;% inlet valve cone angle
Dve=0.0198;% exhaust valve seat throat diameter
Dvs=0.0243;% suction valve cone throat diameter
Dme=0.0229;% exhaust valve average diameter
Dms=0.0275;% inlet valve average diameter
Doe=0.0055;% exhaust valve rod diameter
Dos=0.0055;% inlet valve shank diameter
Alpha_phi_initial=0.9;% compression stroke starts excess air coefficient in cylinder
Ta=374;% compression stroke starts Temperature of Working in cylinder
Pa=0.132850435031838;% compression stroke starts initial pressure in cylinder
Angle_1=250;
Va=(pi*D^2/4) * [S/ (epsilon-1)+(S/2) * ((1+1/Lambda)-cos (pi*angle_1/
180)-(1/Lambda)*sqrt(1-Lambda^2*(sin(pi*angle_1/180)).^2))];% compression stroke starts cylinder
Interior initialization volume
ML=pa*1e6*va/ (9.81* (29.2647-0.0402/alpha_phi_initial) * ta);% ideal gas
Body state equation pv=mRgT obtains working medium quality in the cylinder that compression stroke starts
Each stroke angles of %%
Angle_1=250;
Angle_2=331;
Angle_3=405;
Angle_4=487.6;
Angle_5=708;
Angle_6=741.6;
Angle_7=970;
Phi_z=angle_3-angle_2;% combustion continuation angle
%% initiation parameter
Ta=374;% compression stroke starts Temperature of Working in cylinder
Alpha_phi_initial=0.9;% compression stroke starts excess air coefficient in cylinder
Pa=0.132850435031838;% compression stroke starts initial pressure in cylinder
A1=pi*D^2/4;A2=a1;% cylinder cap and piston area
Tw1=580;Tw2=550;Tw3=480;% compression stroke starts cylinder cap, piston, cylinder sleeve initial temperature
Gf=mL/ (1+L0*alpha_phi_initial);% compression stroke starts initial fuel quality in cylinder
Pr=0.1058904;Exhaust pipe initial pressure after % exhaust valve
Ps=0.096264;Air inlet pipe initial pressure after % inlet valve
Ts=320;Air inlet pipe initial pressure after % inlet valve
Tr=1150;Exhaust pipe initial temperature after % exhaust valve
M=1.717;% combustion quality index
%% certificate parameter
Origin_p_e=45.3;% complete machine effective power
Origin_pmax=6.469;% highest cylinder pressure, initial data
%***************************** is counter to push away eta_m (mechanical efficiency) * * * * * * * * * * * * * * * * * * * * * * * * * *
*************
Origin_pm_i=1.205;%IMEP (mean indicated pressure (MIP) is obtained by image, initial data
Origin_pm_e=origin_p_e*30*4/ (0.25*pi*D^2*S*1000*n*4);% initial data
Origin_eta_m=origin_pm_e/origin_pm_i;% mechanical efficiency
%***************************** is counter to push away eta_m (mechanical efficiency) * * * * * * * * * * * * * * * * * * *
********************
%*********************************************** inlet and outlet cam lift is located in advance
Manage * * * * * * * * * * * * * * * * * * * *
Hve_y=xlsread (' hve_xy.xlsx');% angle -- mm, exhaust valve lift illustrate that xlsread can also be with
Read the excel file that suffix is xlsx.
Hvs_y=xlsread (' hvs_xy.xlsx');% angle -- mm, inlet valve lift
%*********************************************** inlet and outlet cam lift is located in advance
Manage * * * * * * * * * * * * * * * * * * * *
%% other parameters
Hu=44100;% gasoline engine low heat value
Va=(pi*D^2/4) * [S/ (epsilon-1)+(S/2) * ((1+1/Lambda)-cos (pi*angle_1/
180)-(1/Lambda)*sqrt(1-Lambda^2*(sin(pi*angle_1/180)).^2))];% compression stroke starts cylinder
Interior initialization volume
Choice1=input (selection of formula ' instantaneous average heat transfer coefficient ');
The choice2=input selection of r, u, cv, cp, h calculation (' ');
Step=.1;% time step, 0.1CA.
Miu_e_coff_a=0.98;Miu_e_coff_b=3.3;Miu_e_coff_c=2;Miu_s_coff_a=
0.98;Miu_s_coff_b=3.3;Miu_s_coff_c=2;% discharge coefficient relevant parameter;
The %% condition of convergence
While round (choice1)~=choice1 | | choice1≤0 | | choice1 >=5
Choice1=input (' the positive integer of 1-4, re-enter choice1')
end
While round (choice2)~=choice2 | | choice2≤0 | | choice2 >=4
Choice2=input (' the positive integer of 1-3, re-enter choice2')
end
%% function call
Options=odeset (' abs', 1e-10, ' rel', 1e-10);
[phi1, t]=ode45 (' dt_dphi1', [angle_1:step:angle_2], ta, options);%% should
Functional simulation gasoline engine compression stroke is called, the ODE of dT/dphil is obtained.The ODE is solved by ode45
Obtain the t under different crank angle phi1.
%******************************************************* * * * * * * * * * * * *
***************
[f1,t1,p1,m_phi1,alpha_phi1,x_phi1,v_phi1,alpha_g1,k_phi1,r_phi1,dv_
Dphi1, dqw_dphi1, u_phi1, h_phi1, cv_phi1]=dt_dphi1 (phi1, t);% exports intermediate variable, obtains
Function dt_dphi1 all output variables
Dp1_dphi1=((m_phi1*r_phi1) ./v_phi1.^2) .* (v_phi1.*f1-t1.*dv_dphi1) *
1e-006;% is to compression stroke in-cylinder pressure derivation
Ax=4*min (v_phi1)/D;% minimum cylinder sleeve area
Disp (' * * * * * * * * * * * * * * * * * * * * * * * * compressed * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * ')
global angle_2angle_3step
Options=odeset (' abs', 1e-10, ' rel', 1e-10);
[phi2, t]=ode45 (' dt_dphi2', [angle_2:step:angle_3], t1 (end), options);%
The primary condition of burning is Temperature of Working t1 (end) in the cylinder for compressing and ending
[f2,t2,p2,m_phi2,alpha_phi2,x_phi2,v_phi2,dv_dphi2,alpha_g2,xx_phi2,
Dxx_dphi2, k_phi2, r_phi2, u_phi2, h_phi2, cv_phi2, dqw_dphi2]=dt_dphi2 (phi2, t);%
Export intermediate variable
Dp2_dphi2=((m_phi2*r_phi2) ./v_phi2.^2) .* (v_phi2.*f2-t2.*dv_dphi2) *
1e-006;% is to combustion stroke in-cylinder pressure derivation
%******************************************************* * * * * * * * * * * * *
******
Index_tmax=find (t2==max (t2));The corresponding crank angle of working medium maximum temperature in % cylinder
Tmax_phi2=phi2 (index_tmax);Working medium maximum temperature in % cylinder
Index_pmax=find (p2==max (p2));The corresponding crank angle of working medium maximum pressure in % cylinder
Pmax_phi2=phi2 (index_pmax);Working medium maximum pressure in % cylinder
max(dp2_dphi2);The maximum rate of pressure rise in % cylinder
Disp (' * * * * * * * * * * * * * * * * * * * * * * * * burnt * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * ')
global angle_3angle_4step
Options=odeset (' abs', 1e-10, ' rel', 1e-10);
[phi3, t]=ode45 (' dt_dphi3', [angle_3:step:angle_4], t2 (end), options);%
The primary condition of expansion stroke is Temperature of Working t1 (end) in the cylinder for burning and ending
[f3,t3,p3,m_phi3,alpha_phi3,x_phi3,v_phi3,dv_dphi3,alpha_g3,k_phi3,r_
Phi3, u_phi3, h_phi3, cv_phi3, dqw_dphi3]=dt_dphi3 (phi3, t);
Dp3_dphi3=((m_phi3*r_phi3) ./v_phi3.^2) .* (v_phi3.*f3-t3.*dv_dphi3) *
1e-006;% is to expansion stroke in-cylinder pressure derivation
Disp (' * * * * * * * * * * * * * * * * * * * * * * * * expanded * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * ')
global gf angle_4angle_5pr step
% [filename pathname]=uigetfile (' * .xlsx';' * .xls'}, ' select file ');
%str=[pathname filename];
%xls=xlsread (str);
%hv_x=xls (1 :)
%hv_y=xls (2 :)
Options=odeset (' abs', 1e-10, ' rel', 1e-10);
[phi4, t]=ode45 (' dt_dphi4', [angle_4:step:angle_5], [0;gf;t3(end)],
options);The primary condition of the pure exhaust stroke of % includes temperature in exhaust quality, internal combustion oil of cylinder quality and the cylinder of expansion end of a period
Degree finally obtains the corresponding exhaust quality of each crank angle of pure exhaust stroke, internal combustion oil of cylinder quality and cylinder temperature.
Last: the foregoing is only a preferred embodiment of the present invention, is not intended to restrict the invention, all in the present invention
Spirit and principle within, any modification, equivalent replacement, improvement and so on, should be included in protection scope of the present invention it
It is interior.
Claims (5)
1. a kind of Numerical Simulation of Working Process in Gasoline Engine system, which is characterized in that the system is based on MATLAB platform, including is
Unite main program, the system main program include compression main unit, burning main unit, expansion main unit, exhaust main it is mono-
Member folds and opens main unit and air inlet main unit, the compression main unit, burning main unit, expansion main unit, exhaust
Main unit, it is folded open main unit and air inlet main unit is sequentially connected, the system main program is connected with call subroutine, institute
The system main program of stating is also configured with worksheet, and the worksheet includes hve_xy table and hvs_xy table;
Gasoline engine structural parameters, each stroke phase angle, original state parameter and verifying ginseng are provided in the compression main unit
It counts, is provided with inlet valve lift parameter in the hvs_xy table, is provided with exhaust valve lift parameter in the hve_xy table.
2. a kind of Numerical Simulation of Working Process in Gasoline Engine system according to claim 1, which is characterized in that the compression
Main unit, expansion main unit, exhaust main unit, is folded and opens main unit, air inlet main unit and work burning main unit
The source code file for making table is located under same English path, and the system file name of the source code file is by English and digital group
At.
3. a kind of Numerical Simulation of Working Process in Gasoline Engine system according to claim 2, it is characterised in that: system is carrying out
When simulation, specifically comprise the following steps:
Step 1, when operation, system main program is based on MATLAB platform, MATLAB platform is opened, by the work of MATLAB platform
Path is specified on the root of source code file, and modifies inlet valve lift parameter and exhaust valve lift parameter;
Step 2 opens compression main unit, system parameter is arranged, specifically:
It modifies to the parameter for including among the above;
Step 3 runs MATLAB platform, instantaneous to determine according to some number of prompt in command Window input 1~4
Average heat transfer coefficient calculation formula, then some number in 1~3 is inputted, with the instantaneous specific heat at constant volume calculation formula of determination;
Step 4, operation burning main unit, display in command Window " * * * burn * * * " indicate operation completion, successively beat
Follow-up expansion main unit is opened, exhaust main unit, folds and opens main unit and air inlet main unit, clicks operation button, operation
Can show respectively after the completion " * * * has expanded * * * ", " * * * is pure to be vented * * * ", " * * * valve folded ends * * * " and " * * * it is pure into
The complete * * * " of gas, when all system main program operations finish, command Window can show that gasoline engine performance emulates key data result.
4. a kind of Numerical Simulation of Working Process in Gasoline Engine system according to claim 3, it is characterised in that: the step 1
In, it is 360 °C of A that gasoline engine compression top center, which corresponds to time point, with 720 °C of A for a working cycles, including it is air inlet, compression, swollen
Swollen and four strokes of exhaust.
5. a kind of Numerical Simulation of Working Process in Gasoline Engine system according to claim 3, it is characterised in that: to gasoline engine
After the completion of simulation trial, MATLAB platform draws different data images, including working medium temperature in cylinder according to the emulation data of simulation
It writes music line, pressure curve, figure-of-merit curve, excess air coefficient curve, piston displacement, cylinder inner volume, working medium ratio at inlet and exhaust valve
Enthalpy, combustion heat release rate, instantaneous average heat transfer coefficient, indicator card, equivalent amount of fuel, the rate of pressure rise, temperature rise rate, insulation refer to
Number, gas constant, intake and exhaust rate of discharge, inlet and exhaust valve lift, air inlet and exhaust valve discharge coefficient, the effective flow area of inlet and exhaust valve,
Density, specific internal energy, specific enthalpy etc. in instantaneous cylinder.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910830136.8A CN110532716A (en) | 2019-09-04 | 2019-09-04 | A kind of Numerical Simulation of Working Process in Gasoline Engine system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910830136.8A CN110532716A (en) | 2019-09-04 | 2019-09-04 | A kind of Numerical Simulation of Working Process in Gasoline Engine system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN110532716A true CN110532716A (en) | 2019-12-03 |
Family
ID=68666550
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910830136.8A Pending CN110532716A (en) | 2019-09-04 | 2019-09-04 | A kind of Numerical Simulation of Working Process in Gasoline Engine system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110532716A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112861316A (en) * | 2021-01-11 | 2021-05-28 | 西安交通大学 | Method for characterizing cold emission process by using supercritical carbon dioxide as working medium and computer program product |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104102781A (en) * | 2014-07-18 | 2014-10-15 | 上汽通用五菱汽车股份有限公司 | Multi-objective optimization based engine performance optimizing method |
CN104794258A (en) * | 2015-03-03 | 2015-07-22 | 联合汽车电子有限公司 | Automobile hardware-in-loop simulation system |
US20150347648A1 (en) * | 2014-05-30 | 2015-12-03 | Avl List Gmbh | Method for Simulation of an Internal Combustion Engine |
-
2019
- 2019-09-04 CN CN201910830136.8A patent/CN110532716A/en active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150347648A1 (en) * | 2014-05-30 | 2015-12-03 | Avl List Gmbh | Method for Simulation of an Internal Combustion Engine |
CN104102781A (en) * | 2014-07-18 | 2014-10-15 | 上汽通用五菱汽车股份有限公司 | Multi-objective optimization based engine performance optimizing method |
CN104794258A (en) * | 2015-03-03 | 2015-07-22 | 联合汽车电子有限公司 | Automobile hardware-in-loop simulation system |
Non-Patent Citations (1)
Title |
---|
袁佳斌: "基于MATLAB的内燃机工作过程模拟软件开发", 《中国优秀硕士学位论文全文数据库(工程科技Ⅱ辑)》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112861316A (en) * | 2021-01-11 | 2021-05-28 | 西安交通大学 | Method for characterizing cold emission process by using supercritical carbon dioxide as working medium and computer program product |
CN112861316B (en) * | 2021-01-11 | 2022-12-09 | 西安交通大学 | Method for characterizing cold emission process by using supercritical carbon dioxide as working medium and computer program product |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN100535407C (en) | Adjustable high boost system with structure of series-parallel connection of turbochargers in different size | |
US5103645A (en) | Internal combustion engine and method | |
CN102278193B (en) | Auxiliary pressurization system of internal combustion engine | |
CN110543715B (en) | Two-stage supercharging system matching method based on simulation calculation | |
CN103352754A (en) | Multi-mode large and small turbocharger diesel engine sequential supercharging structure | |
CN101949324A (en) | Turbosupercharging through-flow matching method for internal combustion engine | |
CN102678283A (en) | Engine exhaust gas bypass device | |
CN106777553A (en) | A kind of turbocharger matching selection method based on vehicle performance | |
CN110532716A (en) | A kind of Numerical Simulation of Working Process in Gasoline Engine system | |
CN102678282A (en) | Regulating system with air intake and discharge combined mechanism | |
CN102678272A (en) | Regulating device for air intake and discharge flow of supercharger | |
CN2419376Y (en) | Natural gas engine | |
CN102678269A (en) | Supercharged engine air inlet pipe deflating system | |
CN110725745A (en) | Two-stroke engine | |
CN200996322Y (en) | Air supplier for improving tail-gas turbine booster engine response | |
CN206309471U (en) | A kind of two-stage supercharging device of engine | |
CN203201653U (en) | Automobile engine boost device and automobile engine with same | |
CN104632362B (en) | Exhaust gas expanding mechanism for two-stroke pure oxygen engine | |
CN2622403Y (en) | Gas turbine | |
CN2511795Y (en) | Two stroke IC engine with scavenging cylinder | |
CN2811572Y (en) | Intake valve and exhaust valve co-used engine | |
CN114738109B (en) | Matching method for key parameters of two-stage turbocharged internal combustion engine | |
CN104612811B (en) | Mechanical double plate rotary apparatus | |
CN206917735U (en) | A kind of double cylinder combined type two-stage engines | |
RU2196901C2 (en) | Combination internal combustion engine |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20191203 |