高压油管广泛应用于发动机制造,一种高效合理的高压油管压力控制方案,对推动发动机制造业和经济社会发展具有十分重要的意义。通过对高压燃油系统工作原理的研究,提出了一种新的高压油管压力控制方案。首先,根据高压油管进油和喷油质量守恒关系,分别推导出了高压油管在恒压和不恒压条件下的压强与时间的微分方程模型。其次,推导出了高压油管油泵凸轮柱塞运动的角速度计算公式。进一步,考虑到高压油管恒压条件下多个喷油嘴工作的情况,分别讨论了多个喷油嘴不同工作情况下的控制方案。最后,通过对数据进行拟合,检验了模型的准确性。数值结果表明,在恒定压力下单向阀每次开启时长为0.64 ms,高压泵凸轮柱塞角速度为2.7543 rad/ms。多个喷油嘴同时和不同时工作时,高压泵凸轮柱塞角速度大小分别为5.5086 rad/ms、1.37715 rad/ms。此外,本文根据结果制定了一种高压油管压力控制方案,为发动机高压油管的发展提供了理论依据。 High-pressure tubing is widely used in engine manufacturing. An efficient and reasonable high- pressure tubing pressure control scheme is of great significance for promoting engine manufacturing and economic and social development. Through the research on the working principle of high-pressure fuel system, a new high-pressure fuel pipe pressure control scheme is proposed. Firstly, according to the mass conservation relationship between high-pressure oil pipe feeding and spraying, the differential equation model of pressure and time of high-pressure oil pipe under constant pressure and non-constant pressure conditions is derived respectively. Secondly, the formula for calculating the angular velocity of the cam plunger movement of the high-pressure tubing oil pump is derived. Further, considering the operation of multiple injectors under the constant pressure condition of high pressure tubing, the control schemes of different injectors under different working conditions are discussed respectively. Finally, the accuracy of the model was tested by fitting the data. The numerical results show that the constant value of the check valve is 0.64 ms at a constant pressure and the angular velocity of the high pressure pump cam plunger is 2.7543 rad/ms. When multiple injectors are working at the same time and at different times, the angular velocity of the high pressure pump cam plunger is 5.5086 rad/ms and 1.37715 rad/ms, respectively. In addition, based on the results, a high pressure fuel pipe pressure control scheme is developed, which provides a theoretical basis for the development of high-pressure fuel pipe.
陈兴志,黄飞翔,王代文,乐文涛,刘乃瑶
西南科技大学理学院,四川 绵阳
收稿日期: 2019年10月17日;录用日期:2019年11月4日;发布日期:2019年11月11日
高压油管广泛应用于发动机制造,一种高效合理的高压油管压力控制方案,对推动发动机制造业和经济社会发展具有十分重要的意义。通过对高压燃油系统工作原理的研究,提出了一种新的高压油管压力控制方案。首先,根据高压油管进油和喷油质量守恒关系,分别推导出了高压油管在恒压和不恒压条件下的压强与时间的微分方程模型。其次,推导出了高压油管油泵凸轮柱塞运动的角速度计算公式。进一步,考虑到高压油管恒压条件下多个喷油嘴工作的情况,分别讨论了多个喷油嘴不同工作情况下的控制方案。最后,通过对数据进行拟合,检验了模型的准确性。数值结果表明,在恒定压力下单向阀每次开启时长为0.64 ms,高压泵凸轮柱塞角速度为2.7543 rad/ms。多个喷油嘴同时和不同时工作时,高压泵凸轮柱塞角速度大小分别为5.5086 rad/ms、1.37715 rad/ms。此外,本文根据结果制定了一种高压油管压力控制方案,为发动机高压油管的发展提供了理论依据。
关键词 :质量守恒,非线性拟合,模拟仿真,高压油管压力
Copyright © 2019 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
内燃机对世界工业革命和社会经济发展起到了巨大的推动作用 [
在以往的研究中,许多学者分别对高压油管进行了不同的研究。在国内,1990年,赵学亭等人 [
本文通过对高压油管压力控制的研究,提出了一种新的基于质量守恒方程的高压油管压力控制模型。赵学亭等人 [
(一) 高压油管恒压压力控制模型
燃油的压力变化量和密度变化量成正比。比例系数为 E ( P ) ρ ,与弹性模量E有关 [
力变化量和密度变化量的分析,建立弹性模量随压强变化模型如下:
E ( p ) = k e λ p (1)
由此关系建立压力变化量 Δ P 与密度变化量 Δ ρ 关系模型:
d P = k d ρ = E ( P ) ρ d ρ (2)
考虑高压油管工作过程中存在的质量守恒:
m i n + m o l d = m o u t + m n o w (3)
又 m o l d = m n o w , m i n = m o u t 。建立进出气体质量相等的微分方程模型如下( ρ 160 = 0.8712 为160 MPa下的密度):
d ρ = ρ 160 C A 2 ( P A − P ) ρ 160 ∫ v ¯ d t d t (4)
建立了压力与时间关系微分方程模型得到:
d P = k ρ 160 C A 2 ( P A − P ) ρ 160 ∫ v ¯ d t ⋅ exp ( − λ P − e − λ P k λ − D k ) d t (5)
(二) 高压油管不恒压压力控制模型
随着高压油管工作的进行,( m i n + m o l d )一直在不断增加,( m o u t + m n o w )也在同步增加,记( m o u t + m n o w )为m。建立高压油管变动压力的微分模型如下( m o l d 是初始状态下管内的气体质量):
{ d m = ρ 160 C A 2 [ P A + ln ( − λ k ln m V + v ¯ t − λ D ) λ ] d t m 0 = ρ 100 V (6)
由 m i n = m o u t ,建立出高压油管变动压力的微分模型:
d P = k ρ 0 C A 2 × ( P A − P ) ρ 0 ∫ v ¯ d t ⋅ exp ( − λ − e − λ P k λ − D k ) d t (7)
(三) 高压油管恒压单喷油嘴柱塞运动模型
考虑凸轮驱动柱塞由低压状态旋转至高压状态这个转一圈的过程,低压状态下,高压油泵的体积为 V b ,高压油泵在低压状态下气体的质量 m b 。高压状态下,残余体积为v mm3,残留气体的质量 Δ m 。进入高压油管中气体的质量为:
m i n = m b − Δ m (8)
高压油管内气体的密度:
ρ g = m i n + m o l d V o l d (9)
引用高压油管工作过程中存在的质量等式:
M i n + m o l d = m o u t + m n o w (10)
由 m o l d = m n o w 即 M i n = m o u t 。从而得出质量守恒模型:
{ M i n = n m o u t m i n = d ρ d v (11)
建立针阀运动高度模型如下:
h ( t ) = a t 5 + b t 4 + c t 3 + d t 2 + e t + ε (12)
建立喷油体积和时间的关系模型如下:
V = 1 3 ∗ π ∗ ( a t 5 + b t 4 + c t 3 + d t 2 + e t + ε ) ∗ ( R 2 + r 2 + R ∗ r ) (13)
图1. 喷油嘴的动态工作过程及体积示意图
将喷油嘴的形状抽象补充成一个完整的圆锥,如图1所示,令喷孔半径为r,则喷孔到圆锥顶点的距离为:
h 1 = r tan 9 ∘ (14)
设升值为h,则阴影部分圆台的下底面半径为:
R = ( h 1 + h ) tan 9 ∘ (15)
然后分别计算每个区间段的气体体积 V t ,以及针阀的体积 V z ,H为针阀高度,便可以得出如下函数:
V = V t + d V z (16)
考虑一个周期 t = [ 0 , 100 ] 这个时间段,将其分为四个区间,分别为: [ 0 , 0.44 ] , [ 0.45 , 2 ] , [ 2.01 , 2.45 ] , [ 2.46 , 100 ] 将每一个时间段的总体积V表示出来:
{ d V z = π R 2 d H V t = 1 3 π h ( R 2 + r 2 + R r ) (17)
所以,在一个周期内,可以得出其转动次数n,有了转动次数,就可以求出柱塞转一圈所需的时间,进而算出角速度,则有:
ω = 2 π ∗ M i n T ∗ ρ d v (18)
(四) 高压油管恒压多喷油嘴柱塞运动模型
当多个喷油嘴工作时,同样的时间内,释压量翻倍,质量变为原来的两倍,必须加快柱塞的转动次数方能达到平衡。由气体质量的守恒得:
M i n + m o l d = m o u t + m n o w (19)
由 m o l d = m n o w , m o u t = ρ d v 应有:
ω = 2 π ∗ M i n T ∗ ρ d v (20)
(一) 高压油管压力控制模型检验
本文通过微分方程的数值解法和利用MATLAB软件模拟仿真分别对本文模型进行有效性和精确性的检验。根据2019高教社杯全国大学生数学建模竞赛A题的数据 [
图2. 弹性模量随压强变化的曲线图
通过图像可以得知,趋势是一个递增指数型函数,进一步验证了本文所建立的弹性模量与压力模型 E ( p ) = k e λ p 的模型的有效性和准确性。做出拟合图像和原图对比,得到的图像效果如图3所示:
图3. 函数拟合与原图像的对比图
从图中的拟合效果来看,拟合函数图像和原始数据图像几乎重合,所以 E ( p ) 函数的表达式可以确定为一个指数型函数,对其进行非线性函数的拟合,得到模型中的 k = 1476.3877 , λ = 0.00401 ,所以
E ( p ) = 1476.3877 e 0.00401 p
通过对高压油管恒定压力的微分模型的求解,利用计算机模拟仿真画出模型所得到的压强P和时间t的图像,如图4所示:
由图可知,本文选择压强为105 MPa的上限波动最大值,图中的时间是1 s内的开启时间,换算成100 ms的周期,此时的时间为0.64 ms。即每100 ms内开启0.64 ms,压强的增加值为0.5 MPa。
(二) 高压油管恒压喷油嘴柱塞运动模型检验
对所给出的针阀运动曲线数据进行模拟仿真,得到针阀升值随时间的变化曲线如图5所示:
图4. 压强随时间变化图像
图5. 针阀升值随时间的变化曲线
由图5可见,针阀开启和关闭两个过程是相互对称的,为一个分段函数,通过对一个周期时间的划分,经过MATLAB数据拟合仿真得到针阀开启时升程高度变化曲线如下图所示:
由图6可知,原始数据曲线和针阀运动升程高度模型拟合曲线完全重合。针阀运动升程高度模型拟合的效果非常好,进一步验证了针阀运动升程高度模型,对其进行非线性函数的拟合,得到模型中的 a = − 772.3439 、 b = 570.2069 、 c = − 96.5000 、 d = 8.8246 、 e = − 0.2908 、 ε = 0.0017 。所以针阀运动升程高度模型为:
h ( t ) = − 772.3439 t 5 + 570.2069 t 4 − 96.5 t 3 + 8.8246 t 2 − 0.2908 t + 0.0017
图6. 针阀开启时升程高度变化拟合曲线图
综上所述,通过利用现有数据对模型进行检验,进一步说明了本文所建立基于质量守恒的高压油管压力微分模型具有较高准确性、有效性。
(一) 高压油管恒压压力控制模型求解结果
针对2019高教社杯全国大学生数学建模竞赛A题的数据 [
图7. 压强随时间变化图像
由图可知,本文选择压强为105 MPa的上限波动最大值,图中的时间是1 s内的开启时间,换算成100 ms的周期,此时的时间为0.64 ms。即每100 ms内开启0.64 ms,压强的增加值为0.5 MPa。
(二) 高压油管不恒压压力控制模型求解结果
通过对高压油管不恒定压力的微分模型的求解,分别得到2 s、5 s、10 s的压强与时间的对应关系。分别如下图8、图9、图10所示:
图8. 0~2秒压强随时间的变化曲线
图9. 0~5秒压强随时间的变化曲线
图10. 0~10秒压强随时间的变化曲线
三个时间段下m-t,ρ-t,P-t关系的图像汇总如图11所示:
图11. 三个变量随时间t的分段变化曲线
根据三种时间段下,压强和时间的关系曲线,分别找出数据集中所有的开启时间段以及关闭时间段并求和取平均。得到的高压油管油泵和喷油嘴的开关策略为如下表1所示:
由图像和表中数据可知,在0~2 s的调整期内,喷油嘴在一个喷油周期内开启时长为148.0 ms,关闭时间为63.5 ms。在0~5 s的调整期内,喷油嘴在一个喷油周期内开启时长为146.5 ms,关闭时间为62.0 ms。在0~10 s的调整期内,喷油嘴在一个喷油周期内开启时长为148.0 ms,关闭时间为 61.0 ms。
调整时段 | 开启的时间 | 关闭的时间 |
---|---|---|
0~2 s | 148.0 ms | 63.5 ms |
0~5 s | 146.5 ms | 62.0 ms |
0~10 s | 148.0 ms | 61.0 ms |
表1. 高压油管油泵和喷油嘴的开关策略表
(三) 高压油管恒压单喷油嘴柱塞运动模型求解结果
通过对高压油管恒定压力下单喷油嘴柱塞运动模型的求解,计算出公式中需要的常数量: R = 1.25 mm ,
r = 0.7 mm , h 1 = r tan 9 ∘ = 4.4192 mm 常数量代入模型中,得道各个m的质量:
{ m 1 = 132.3578 1 , t = [ 0 , 0.44 ] m 2 = 1293.4526 , t = [ 0.45 , 2 ] m 3 = 12.5873 , t = [ 2.01 , 2.45 ] m 4 = 0 , t = [ 2.46 , 100 ]
则有: n = m 1 + m 2 + m 3 + m 4 m i n = 1.074 , t z = T n = 2.2812 ms , ω = 2 π t z = 2.7543 rad / ms 。
当高压油管单喷油嘴柱塞在恒定压力下运动时,在一个喷油周期内,高压油管喷油嘴喷油一次,相当于凸轮柱塞转动1.074次,高压油泵凸轮柱塞每转动一圈用时2.2812 ms。高压油泵凸轮柱塞转动的角速度为2.7543 rad/ms。
(四) 高压油管恒压多喷油嘴柱塞运动模型求解结果
通过对高压油管恒定压力下多喷油嘴柱塞运动模型的求解,计算出公式中需要的常数量: R = 1.25 mm ,
r = 0.7 mm , h 1 = r tan 9 ∘ = 4.4192 mm 代入上面公式中,计算出各个m的质量:
{ m 1 = 132.3578 1 , t = [ 0 , 0.44 ] m 2 = 1293.4526 , t = [ 0.45 , 2 ] m 3 = 12.5873 , t = [ 2.01 , 2.45 ] m 4 = 0 , t = [ 2.46 , 100 ]
分为两种情况进行讨论:
(1) 多个喷油嘴同时开关,则:
n = 2 ( m 1 + m 2 + m 3 + m 4 ) m i n = 2.148
t z = T n = 1.1406 ms
ω = 2 π t z = 5.5086 rad / ms
(2) 多个喷油嘴不同时工作,则:
n = m 1 + m 2 + m 3 + m 4 m i n = 1.074
t z = T 2 n = 1.1406 ms
ω = 2 π t z = 1.37715 rad / ms
当高压油管多个喷油嘴柱塞在恒定压力下运动,多个喷油嘴同时开关时,在一个喷油周期内,高压油管喷油嘴喷油一次,相当于凸轮柱塞就转动2.148次,高压油泵凸轮柱塞每转动一圈用时1.1406 ms。高压油泵凸轮柱塞转动的角速度为5.5086 rad/ms。多个喷油嘴不同时工作时,在一个喷油周期内,高压油管喷油嘴喷油一次,高压油泵凸轮柱塞就转动1.074次,高压油泵凸轮柱塞每转动一圈用时1.1406 ms。高压油泵凸轮柱塞转动的角速度为1.37715 rad/ms。
本文通过所建立的基于质量守恒方程的高压油管压力控制模型对2019高教社杯全国大学生数学建模竞赛A题数据 [
1) 本文通过建立的基于质量守恒方程的高压油管压力控制模型进而推导得到压强与时间、压强与密度的数学模型,在质量守恒的条件下,通过计算机模拟仿真可以直观的得到压强与时间、压强与密度的关系和变化规律。
2) 燃油进入和喷入高压油管是燃油机的工作原理,研究燃油机的高压油管压力控制意义深刻,本文通过对燃油机高压油管压力控制研究后,解决了高压油管压力控制方面的问题,得出了高压油管的压力控制相关方案。
3) 在分析高压油管工作原理和燃油压力关系的基础上,本研究运用微分方程建立了基于质量守恒方程的高压油管压力控制模型对压力变化进行模拟仿真,并获得了较为理想的仿真结果结果。在今后可将此模型用于压力变化的问题之中。
西南科技大学大学生创新基金项目(项目编号:CX19-061)。
陈兴志,黄飞翔,王代文,乐文涛,刘乃瑶. 基于质量守恒方程的高压油管压力控制研究Research on Pressure Control of High Pressure Tubing Based on Mass Conservation Equation[J]. 应用数学进展, 2019, 08(11): 1732-1744. https://doi.org/10.12677/AAM.2019.811203