供油系统对其柴油机的性能具有重要的影响。根据高压油管的相关工作原理,本文给出了针对燃油进入和喷出过程对应的动力系统数学模型,并利用MATLAB对高压油管的运作过程进行了数值模拟。相比现存的模型,本文的模型具有更好的可操作性和有效性。 The fuel supply system plays an extremely important part in the performance of diesel engine. One mathematical model of dynamical system for the process of fuel entering and ejecting is given in this paper based on the operating principle of high-pressure tubing. Simulation of the operation process of high-pressure oil tube is added to illustrate our results. Compared with the models in the literature, ours is more feasible and useful.
顾喆奕,郑子豪,朱一曾雄,张慧增,李安水*
杭州师范大学数学系,浙江 杭州
收稿日期:2019年11月9日;录用日期:2019年11月28日;发布日期:2019年12月6日
供油系统对其柴油机的性能具有重要的影响。根据高压油管的相关工作原理,本文给出了针对燃油进入和喷出过程对应的动力系统数学模型,并利用MATLAB对高压油管的运作过程进行了数值模拟。相比现存的模型,本文的模型具有更好的可操作性和有效性。
关键词 :流通模型,高压油泵运作模型,曲线拟合,参数优化
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/
在燃油喷射系统中,低压燃油通过高压油泵的压缩进入高压油管,再从喷油嘴处喷出。由于燃油的流入和喷出是间歇性的,这会导致高压油管内压力产生波动,从而使得每次喷出的燃油量都不同,进而影响到发动机的工作效率。
对此我们要对以下问题进行探讨研究,以实现对高压油管的压力控制。
1) 该问题讨论的高压油管其内腔长500 mm,直径10 mm。高压油管左侧有一供油入口A,小孔直径1.4 mm,供油持续的时间长短由单向阀控制,且单向阀每打开一次后需要关闭10 ms。在高压油管右侧有一喷油器每次喷油2.4 ms,工作周期为100 ms,从题目中可以得到喷油器工作时喷油的速率变化图。高压油泵在入口A提供160 MPa的恒定压力,且高压油管内的初始压力为100 MPa。在第一问中,我们需要找到一个合适的单向阀每次的开启时长t,使得高压油管内的压力尽可能稳定在100 MPa左右。而在第二问中,我们需要找到合适的单向阀开启时长使得管内压力分别在经过约2 s、5 s和10 s后从100 MPa增加到150 MPa,并稳定下来。
2) 在系统的运行过程中,高压油管处的燃油来自高压油泵的柱塞腔出口,柱塞腔的内直径5 mm。通过喷油嘴的针阀控制喷油。每次从油泵流入油管的燃油流量受凸轮的控制,题目所给的附件1 [
3) 在第三个问题中,首先,在题目二的基础上增加一个喷油嘴,对此情况进行分析讨论,找到最合适的供油与喷油方案。其次,在上述基础上,往高压油管左侧增加一个单向减压阀,其直径为1.4 mm,该阀可以在打开后将高压油管的燃油释放到低压管道中。我们将对此状况进行建模分析讨论,给出合适的供油与减压阀控制方案。
本论文所要解决的具体问题,可参见全国大学生数学建模竞赛网站2019年试题 [
针对问题一,由题意,我们根据燃油进出高压油管系统的结构特征以及其工作原理,结合题目中注1所给提示,利用密度公式、流体在管道中流动的特性、流量相关公式以及注2所给公式,对包含燃油的压力变化量与密度变化量的公式进行推导整理,最终构建高压油管燃油流通计算模型。当进入高压油管中的燃油流量与喷出的燃油流量相同时,管内的压力能保持稳定状态,由此我们可以借助建立的模型对满足上述条件下的单向阀开启时长进行求解。第二小问,利用MATLAB对高压油管内压力的变化过程进行计算,得到在经过约2 s、5 s和10 s后,管内压力从100 MPa增加到150 MPa,并稳定下来的单向阀开启时长。
针对问题二,利用高压油泵的工作原理,建立类似问题一中模型的高压油泵模型。结合模型一、二,为了使得高压油管内的压力尽量稳定在初始压力,需要让高压油泵向油管输送的燃油量与喷油嘴处喷出的量尽可能相等。喷出量我们可以通过针阀的周期性运动计算得到,输送量与凸轮的角速度有着直接关系,由此我们便可确定满足题目条件的凸轮角速度。
针对问题三,首先增加一个喷油嘴,即增加了燃油的流出流量,这使得原先的供油可能不足以满足需求,结合问题二的分析研究,我们可以通过调整凸轮的角速度增加燃油的输送量。其次再安装一个单向减压阀,它可以使燃油回流到外部低压油路处,从而使得油管内压力减小,即导致燃油的流出量的增加。我们设计了两个方案,一是通过高压油管内的压力,二是设置开启的时机和时长来控制单向减压阀。
为了模型建立的方便,我们提出了以下这几条基本假设:
1) 假设在一次喷射过程中燃油的温度不变;
2) 不考虑高压系统中各个零件的形变状况的影响;
3) 忽略压力传导所消耗的时间;
4) 不考虑加工问题导致的泄露;
5) 假设题目所给的数据真实可靠。
我们给文中所用符号及其说明,见表1。
符号 | 意义 | 符号 | 意义 |
---|---|---|---|
ρ | 燃油密度 | E | 弹性模量 |
R | 高压油管内半径 | L | 高压油管内腔长度 |
S | 截面面积 | x | 流体流动位移 |
v | 流体流速 | Q 1 | 流入高压油管的燃油流量 |
Q 2 | 从高压油管流出的燃油流量 | d A | 供油入口A处小孔直径 |
n | 单向阀每秒工作次数 | t | 单向阀每次开启时长 |
喷油器每秒工作次数 | P | 燃油压力 | |
h | 柱塞高度 | d z | 柱塞腔内的直径 |
ω | 凸轮的角速度 | α | 喷油嘴密封座圆锥半角 |
表1. 相关符号
根据题意,我们希望建立关于燃油进入和喷出间歇性工作对高压油管内压力造成影响的模型,用以描述高压油管内燃油的流动。通过查阅大量资料 [
燃油的压力变化量与密度变化量成正比,其比例系数为 E / ρ ,其中 ρ 为燃油的密度,密度与压力有关,E为弹性模量,它与压力的关系可以通过题目所给的附件3 [
根据相关原理,我们可以得到如下公式:
d P = E P ⋅ d ρ (1)
通过密度公式 ρ = m / v ,又由于高压油管在燃油运输过程中是刚性的,我们将密度公式代入公式(1)后得到:
d P = E ρ ⋅ d m V (2)
其中V由高压油管的内半径R与管道内腔长度L决定,表达式为 V = π R 2 L 。
由于流体在管道中流动满足如下公式:
m = S ρ x (3)
其中S代表截面面积, ρ 为流体密度,x为流体在管道中流动位移,我们可将公式(2)改写为:
d P = E V ⋅ S d x (4)
又因为 v = d x / d t ,我们可以得到:
d P = E V ⋅ S v ⋅ d t (5)
通过流量 = 截面积 × 流体速率,即 Q = S v ,将此代入公式(5)并整理可得到:
d P d t = E V ⋅ Q (6)
其中Q代表流入高压油管的燃油流量 Q 1 与流出的燃油流量之差
1) 高压油泵向高压油管流入的燃油流量 Q 1 :
Q 1 = C A 1 2 Δ P ρ = C A 1 2 | P 2 − P 1 | ρ (7)
式中,流量系数 C = 0.85 ,流通面积 A 1 = π ∗ d A 2 / 4 ( d A 为供油入口A处小孔直径), | P 2 − P 1 | 为高压油泵与油管交界处前后燃油压力差, ρ 为高压侧成燃油密度。
2) 高压油管从喷油嘴流出的燃油流量
Q 2 = C A 2 2 | P 2 − P 1 | ρ (8)
式中流量系数 C = 0.85 ,流通面积 A 2 = π ∗ d B 2 / 4 (
综上所述,我们可以得到高压油管燃油流通模型:
d P d t = E V ⋅ ( Q 1 − Q 2 ) = E V ⋅ ( C A 1 2 | P 2 − P 1 | ρ − C A 2 2 | P 3 − P 2 | ρ ) (9)
通过查阅相关文献 [
V E ⋅ d P d t = C A 1 2 | P 2 − P 1 | ρ − C A 2 2 | P 3 − P 2 | ρ (10)
通过公式(9)我们可以发现,在高压油管内压力保持不变的情况下,即
由题目中给出的喷油速率-时间图像,我们可以计算得到高压油管从喷油嘴流出的燃油流量 Q 2 为:
Q 2 = ( 2 + 2.4 ) × 20 2 × n 2 (11)
其中 n 2 为喷油器每秒工作次数。由此计算得到 Q 2 = 440 mm 3 / s ,它表示一秒内喷油嘴喷出的燃油流量。
而高压油泵向高压油管流入的燃油流量 Q 1 ,根据公式(7)并将具体数值代入,我们可以得到:
正常情况下,燃油密度 ρ 跟随压力变化的关系式 [
ρ = ρ 0 ( 1 + 0.6 × 10 − 9 P 1 + 1.7 × 10 − 9 P ) (13)
其中 ρ 0 表示在压力为100 MPa时燃油的密度, ρ 0 = 0.850 mg / mm 3 。将 P = 160 MPa 代入,即可得到此时 ρ = 0.914 mg / mm 3 。将此结果代入公式(12),即可得到流入的燃油流量 Q 1 = 14.9916 mm 3 / ms ,它表示在单向阀开启时间下,每毫秒(ms)喷出的燃油流量,由于需要 Q 1 = Q 2 ,且单位相同,我们可以得出 Q 1 ∗ n t = Q 2 ,表示每秒钟喷出的燃油流量与进入的流量相同,n表示一秒内开启单向阀次数,t表示每次开启单向阀时间。再与 1000 ms = n t + 10 n 进行联立,最终计算得到单向阀每次开启的时长约为0.302 ms。
问题一的第二小问是通过调整单向阀每次开启的时长,使得高压油管内的压力在分别经过约2 s、5 s和10 s后稳定在150 MPa。我们采用MATLAB迭代的方式,模拟出高压油管中的压力从100 MPa到150 MPa直至稳定的过程。根据问题要求,分别在2 s、5 s、10 s内完成该过程,并求出不同时间要求下单向阀开启的时间。下面我们以2 s为例,对该问题进行求解。
在第一小问中我们已经求得当单向阀开启时间约为0.302 ms时,此时高压油管内的压力几乎处于稳定状态,即管内压力稳定在100 MPa左右。因此我们将单向阀开启时间约束在0.3~1.2 ms范围内,并以0.001 ms为步长,观察不同的时长对管内压力的影响,求得既满足在2 s内完成高压油管内压力从100 MPa增加到150 MPa并稳定在150 MPa,又满足单向阀开启时间最短的阀开启时间t。其次,我们以0.1 ms为步长设置0~12,000 ms的循环来表示时间的推移,记为T。在MATLAB中设置矩阵 P ( i , k ) 用于存储不同阀值不同时刻下燃油的压力。根据前文公式(6)得到燃油压力的计算公式为:
P = E V ⋅ ( Q 1 − Q 2 ) + P 前 (14)
其中E为弹性模量,V为高压油管体积, Q 1 为流入高压油管的燃油流量, Q 2 为流出的燃油流量, P 前 为前一时刻燃油压力。
1) 弹性模量E
我们通过对题目所给的附件3 [
E = 1532 + 1 × 10 − 4 P 3 − 1.082 × 10 − 3 P 2 + 5.474 P (15)
下图图1为我们利用题目所给的附件3 [
图1. 高压油管体积V
V = L π R 2 (16)
2) 流入高压油管的燃油流量 Q 1 :
Q 1 = C A 1 2 | P 2 − P 1 | ρ = 0.85 × π 1.4 2 4 2 ⋅ | 160 − P 前 | ρ 前 (17)
其中 P 前 为前一时刻燃油压力,
3) 高压油管流出的燃油流量 Q 2 :
计算 Q 2 时我们需要判断当前时间喷油器是否开启或是处于开启时的什么状态。喷油器每秒工作10次,每次工作时喷油时间为2.4 ms,由此将100 ms记为一个工作周期,根据题目中所给的喷油速率示意图,我们通过对T除以100的余数(记为rem)跟图示时间进行比较,取该时刻下与前一时刻时的燃油流量的均值,即 ( Q 前 + Q 后 ) / 2 作为最终所要求的 Q 2 ,整理化简后我们得到如下公式:
Q 2 = { 0 , r e m > 2.4 or r e m = 0 − 100 × r e m + 245 , 2.2 < r e m ≤ 2.4 20 , 0.2 < r e m ≤ 2.2 100 × r e m − 5 , 0 < r e m ≤ 0.2 (18)
综上所述,利用公式(14)~(18)可求得不同阀值不同时刻下燃油的压力P在矩阵
同理在经过约5 s和10 s的计算方法与以上相同,完整的MATLAB程序代码见附录2,由此我们得到分别经过约2 s、5 s和10 s后高压油管内压力从100 MPa增加到150 MPa并稳定时,单向阀开启时长的设置,结果如下表2所示:
时间(s) | 2 | 5 | 10 |
---|---|---|---|
单向阀开启时长(ms) | 0.883 | 0.750 | 0.699 |
表2. 单向阀开启时长
图2. 单向阀不同开启时长压力随时间变化示意图
上图2是我们利用MATLAB迭代对高压油管内压力进行仿真模拟,得到满足问题一第二小问要求的不同单向阀开启时长下,高压油管内压力随时间变化的示意图,其中横轴表示时间(ms),纵轴表示压力(MPa)。橙色线表示经过约2 s时间达到要求的单向阀开启时长,黄色线表示经过约5 s,紫色线表示经过约10 s的情形。从图中我们可以发现刚开始时,高压油管内的压力上升幅度很大,很快就到达150 MPa附近,之后逐渐向150 MPa逼近,最后围绕150 MPa上下波动,即高压油管压力值在合理的范围内波动,表明燃油喷射的稳定性。
1) 凸轮
通过查阅题目所给的附件1 [
图3. 凸轮的边缘曲线图
该凸轮对柱塞位置h的影响可以记为 h ( θ ) ,又因为 θ = ω t ,所以h是一个关于 ω 和t的函数,并呈现出周期性的性质,其周期 T = 2 π 。
根据题目中的描述以及查阅关于高温油泵的资料 [
图4. 高压油泵简单示意图
由上图4可以分析出高压油泵的工作原理:随着凸轮的旋转,柱塞上下起伏,从而对上方的柱塞腔内燃油进行压缩,增大压力直至高压油泵内燃油压力大于高压油管内的燃油压力,打开单向阀,将高压油泵中的燃油输入到高压油管中。
2) 针阀
由题意,我们绘制出针阀的运作图像:
图5. 针阀的运作图
上图5中,D代表针阀的直径,d代表最下端喷口的直径, α 为密封座圆锥的半角,h为针阀距最下端喷口的距离,在图中用棕色线条表示。其中燃油流通面积S可以表示为 S = S R − S r , S R 表示大圆锥的侧面积, S r 表示小圆锥的侧面积。
圆锥侧面积计算公式为 S = π R L ,R为圆锥半径,L为圆锥底面周长。
{ L R = D 2 ⋅ 1 cos α + h ⋅ sin α L r = D 2 ⋅ 1 cos α R R = cos α ⋅ L R R r = D 2 (19)
通过上述公式,我们可以表示出燃油的流通面积S:
S = π R R L R − π R r L r = π cos α ⋅ ( D 2 ⋅ 1 cos α + h sin α ) 2 − π D 2 4 ⋅ 1 cos α (20)
假设柱塞与凸轮的强度和刚性足够大且没有变形,在运动过程中没有冲击震动的现象,根据高压油泵的原理,我们可以建立了类似模型一的高压油泵模型:
d P c d t = E V c ( Q c 1 − Q c 2 ) (21)
其中
1) 压入高压油泵的燃油流量 Q c 1 :
通过题目我们可以发现当柱塞处于下止位点时,柱塞腔内充满了低压燃油,且压力为0.5 MPa,随着凸轮的运作,柱塞向上运动时,低压燃油被迅速压缩,柱塞腔内压力增大,在这个过程中进入的燃油量可以用公式 Q = S v 表示,其中S代表柱塞的底面积,利用公式 S = π ⋅ d z 2 / 4 得到,v可以通过 d h / d t 计算得到,即
Q c 1 = π ⋅ d z 2 4 ⋅ d h d t (22)
其中
2) 流入高压油管的燃油流量 Q c 2 :
当柱塞向上挤压柱塞腔内的燃油时,柱塞腔容积减小,燃油的压力增加。当高压油泵内的燃油压力 P c 大于高压油管内的燃油压力
Q 2 = C A 2 Δ P ρ = C A 2 | P c − P m | ρ c (23)
其中A为小孔的截面面积,
将模型一、二结合起来,以 Q c 2 = Q m 1 作为两个模型间的桥梁,其中
{ d P c d t = E V c ( Q c 1 − Q c 2 ) d P m d t = E V m ( Q m 1 − Q m 2 ) Q c 1 = π ⋅ d z 2 4 ⋅ d h ( ω , t ) d t Q c 2 = C A c 2 | P c − P m | ρ c Q m 1 = Q c 2 Q m 2 = C A m 2 | P m − P o u t | ρ m (24)
其中
根据题意我们知道凸轮控制柱塞的上下运动。柱塞向上运动时柱塞腔容积减少,高压油泵内压强增大,当高压油泵内压强高于高压油管内压强时(高压油管内压力初始值100 MPa),开启单向阀,向高压油管中输入燃油;柱塞向下运动时柱塞腔容积增大,从外部吸入0.5 MPa的低压燃油,当柱塞达到下止点时柱塞腔内充满0.5 MPa的低压燃油。当凸轮角速度增大时,单位时间内凸轮带动柱塞上下运动的频率上升,由于单位时间内单向阀开启的时间与柱塞上下运动的频率呈正相关,这使得单位时间内从高压油泵流到油管中得燃油流量增加。又由于喷油嘴B处的开启与关闭呈现周期性变化,因此单位时间内喷油嘴B处的燃油流量基本处于定值。
根据公式 P = E / V ⋅ Δ Q ⋅ T ,存在一个凸轮角速度 ω 0 ,使得单位时间内,开启单向阀时,高压油泵向高压油管输入的燃油流量与喷油嘴B处喷出的燃油流量几乎相等,使得高压油管内压强基本稳定在100 MPa。
由此对于问题二的求解,我们以上面的模型为基础,在MATLAB中进行迭代,计算得到能够使管道压力 P m 稳定在100 MPa周围的角速度 ω 0 = 0.048 rad / ms 。下图为不同角速度下,油管内压力随时间的变化图。
图6. 不同凸轮角速度压力随时间的变化
上图6中,横轴为时间推移,纵轴为高压油管内压力变化。根据上图我们可以发现:
当凸轮角速度小于 ω 0 时,高压油泵向高压油管输入的燃油流量小于喷油嘴B处放出的燃油流量,高压油管内压力处于逐渐降低的状态。
当凸齿轮角速度大于 ω 0 时,高压油泵向高压油管输入的燃油流量大于喷油嘴B处放出的燃油流量,高压油管内压强处于逐渐升高的状态。
在问题三中,我们对模型一进行改进,使其拥有两个喷射口。
其中 Q m 1 代表从高压油泵模型中流入的燃油流量, Q m 2 代表从喷射器1流出的燃油流量,
1) 双喷油嘴同时运作
我们先使两个喷油嘴在同一时刻运作,得到的结果如下:
图7. 两个喷油嘴共同工作时油管内压力随时间变化图像
上图7横轴表示时间从0~10 s,纵轴表示高压油管内的压力。从图中我们发现由于增加了一个喷油嘴,喷油流量超过了供油流量,导致高压油管内的压力随时间逐渐下降,为了保持压力稳定,我们需要对高压油泵的供油流量进行调整。
通过问题二的求解过程,我们知道随着凸轮角速度
图8. ω = 0.096 rad / ms 时管内压力随时间的变化图像
从图8中可以发现当 ω 为0.096 rad/ms时,油管内压力围绕100 MPa上下波动。
接下来为了衡量管内压力的波动程度,我们计算了方差
σ 2 = ( P m i − P s ) 2 , i = 0 , 1 , ⋯ , 100000 (26)
当
2) 两个喷油嘴邻接运作
第一个喷油嘴运行结束后,第二个喷油嘴紧接着开始运作的状况下,管内压力变化如图所示:
图9. 两个喷油嘴邻接运作油管内压力随时间变化图像
从图9中我们可以发现,此运作状态下油管内压力也围绕100 MPa上下波动,这种状况下的 ω 等于0.192 rad/ms,其方差为1.6890 × 107。由此可以发现邻接工作时的方差大于同时运作时的方差,即邻接工作时管内压力的波动程度更大,说明这种方式劣于同时运作。
两个喷油嘴均匀运作
由题目所给的附件2 [
图10. 两个喷油嘴均匀运作油管内压力随时间变化图像
从图10中我们可以发现,此运作状态下油管内压力也围绕100 MPa上下波动,这种状况下的 ω 等于0.192 rad/ms,其方差为4.9560 × 106。由此可以发现均匀工作时的方差小于同时运作时的方差,即均匀工作时管内压力的波动程度更小,说明这种方式优于同时运作。
综上所述,通过测试不同的喷油嘴工作方案,我们发现当两个喷油嘴均匀工作时的方差最小,且这种方案下的凸轮角速度为0.192 rad/ms。
通过查阅文献 [
d P m d t = E V m ( Q m 1 − Q m 2 − Q m 3 − Q m 4 ) (27)
其中新加入的变量 Q m 4 可以表示为如下形式:
其中 δ = { 0 , 单 向 阀 关 闭 1 , 单 向 阀 开 启 。
基于上述模型改进,我们可以找到减压阀的两种控制方案:
方案一:由时间t决定是否开启单向减压阀;
方案二:由高压油管内的压力 P m 决定是否开启单向减压阀。
接下来我们将分别对这两种状况利用MATLAB进行迭代计算,通过调整对应的约束条件观察两种控制方案是否能够有效地使管内压力保持稳定,并计算出满足这个条件的情况下的供油状况。
1) 方案一:由时间t决定是否开启单向减压阀
在方案一中我们控制单向减压阀的条件为开启时长与开启时间间隔。利用MATLAB计算在不同开启时长和开启时间间隔情况下凸轮角速度 ω 和方差
图11. 方案一时管内压力随时间变化图像
从图11中我们可以观察到方案一时压力波动十分剧烈,所以我们认为方案一的效果并不好。
2) 方案二:由高压油管内的压力 P m 决定是否开启单向减压阀
该方案中,我们将设置通过不断改变单向减压阀开关条件,来观察高压油管内压力 P m 的图像变化。当管内压力与低压燃油的压力差,即与单向减压阀出口内外压力差大于某个阈值时开启该减压阀。
Q m 5 = δ C A 5 2 | P m − P 低 压 燃 油 | ρ m , δ = { 0 , P m − P 低 压 燃 油 < P s 1 , P m − P 低 压 燃 油 ≥ P s (29)
其中 P s 为我们设定的压力差阈值。利用MATLAB计算得到能使方差最小的结果如图所示:
图12. 方案二时管内压力随时间变化图像
根据图12我们可以发现管内的压力上涨被减压阀严格控制,在该方案下的凸轮角速度为0.63 rad/ms,其方差为3.8973 × 106。
综上所述,针对该问我们最终确定的供油方案为凸轮角速度 ω = 0.63 rad / ms ,单向减压阀控制方案为当管内与低压油管压力差大于101 MPa时,开启单向减压阀。
1) 本文在模型的建立过程中,用到了密度方程、燃油流动的方程以及系统零件的运动方程等,为仿真提供了理论依据。
2) 本文通过对系统结构和工作原理的分析,清楚系统的整个物理过程,建立了系统的详细数学模型,在此基础上,利用MATLAB编写计算仿真的模型。
1) 本文在利用MATLAB模拟仿真时,对高压油泵、高压油管、喷油器等仿真模型都做了一些简化,会导致仿真结果有所误差,同时模拟仿真时内部也会出现系统的误差。
2) 本文在研究过程中,忽略了燃油发热等造成的系统内部温度变化对整个系统的影响,对模型的计算结果可能会造成一定的影响。
此文中在模型建立过程中并没有考虑温度的变化造成的影响,可以加入燃油的热效应对系统的影响,对模型的个别参数进行一定的修正。本文研究的整个系统其实是一项涉及多领域的复杂而系统的工作,在建模中应尽量全面的考虑影响该系统的因素,对整个数学模型进行细化研究,提高系统计算的准确性。
浙江省自然科学基金(LQ18A010007),国家自然科学基金(11901145)。
顾喆奕,郑子豪,朱一曾雄,张慧增,李安水. 高压油管的压力控制问题研究On the Problem of Pressure Control of High-pressure Oil Tube[J]. 应用数学进展, 2019, 08(12): 1921-1936. https://doi.org/10.12677/AAM.2019.812221