燃气轮机燃油调节剂的计算机仿真研究.docxVIP

燃气轮机燃油调节剂的计算机仿真研究 1 燃油控制系统研究的意义 近年来，随着科学技术的进步，电动汽车的燃油控制系统也变得简单和复杂，从机械液压式燃油控制发展到完全授权数字电子控制（fadc）。目前, 欧美等技术先进国家新研制的发动机几乎都采用全权限数字电子控制系统, 但我国由于种种原因, 目前现役的燃气轮机仍然大量使用机械液压式燃油控制系统。而燃油控制系统中最重要的组成部分就是燃油调节器, 因而对其开展研究对燃油控制系统的改造升级来说具有重大的工程实际意义。 本文以某舰用燃气轮机燃油调节器为研究对象, 通过分析对象的原理和特点, 采用机理分析的建模方法建立了调节器的精确数学模型, 并以matlab/simulink为仿真平台对其进行了仿真研究, 考虑了随动活塞的动态平衡过程, 更加贴切实际的工作过程。 2 定压差流量控制系统模型 某舰用燃气轮机的燃油调节采用的是定压差流量调节, 通过保持油门前后的油压差恒定来实现对燃气轮机的燃油调节。这种调节的优点在于能够使燃气轮机在各种工况下都能够保持较好的自稳定性。 定压差流量调节器主要包括油门开关、定压差活门、计量器随动活塞三个组件组成, 如图1所示。 本文所建立的模型是基于以下的假设: 1) 忽略油液粘度的变化;2) 不计分油活门阻力和管路损失;3) 忽略活门的间隙泄漏。 2.1 弹簧下活门运动方程 定压差活门的工作原理如图2所示, 当油门前后压力超过设定值, 伺服油压进入计量器的活塞腔, 计量器活塞随之运动, 使压差回复到设定值。设初始平衡位置时弹簧的伸缩量为c0, 则定压差活门的运动微分方程为 p1Ac=mcd2ydt2+p2Ac+kc(-c0+y)+mcg(1)p1Ac=mcd2ydt2+p2Ac+kc(?c0+y)+mcg(1) 式中:p1为油门前的压力;p2为油门后的压力;Ac为定压差活门的截面积;mc为定压差活门的质量;kc为定压差活门弹簧的弹性系数;y为定压差活门的开度。 2.2 td2ytdr2+btdyt仍是空白 设计量器的初始位置为yj0, 此时弹簧的预紧压缩量为cj0, 则可以得出滑阀的运动微分方程: pft(Ah-Af)=mtd2ytdt2+Btdytdt+pttAh+kj(y+cj0)(2)pft(Ah?Af)=mtd2ytdt2+Btdytdt+pttAh+kj(y+cj0)(2) 式中:pft为滑阀腔的压力;ptt为弹簧腔的压力;Ah为计量器活塞的截面积;Af为计量器滑阀的截面积;mt为活塞与滑阀的总质量;Bt为活塞运动的粘性阻尼系数;kj为计量器弹簧的弹性系数。 2.3 弹簧腔内的容积变化量q 当定压差活门向上运动时, 经节流口1流入阀杆腔内的燃油应与计量器阀杆腔的容积变化量相等;同理, 弹簧腔中的燃油经节流口2排出的燃油也应与弹簧腔的容积变化量相等, 则: Q1=u1A1√2g(pd-pft)γ=dyjdt(Ah-Af)(3)Q2=u2A2√2g(ptt-p0)γ=dyjdtAh(4) 式中:Q1、Q2分别为分油窗口1和2的燃油体积流量;u1、u2分别为分油窗口1和2的流量系数。 当定压差活门向下运动时, 计量器阀杆腔中的燃油经节流口排出的燃油应等于阀杆腔的容积变化量;经过节流口3流进计量器弹簧腔的燃油流量应等于弹簧腔的容积变化量;同理, 弹簧腔中的燃油经节流口4排出的燃油也应与弹簧腔的容积变化量相等, 则: Q3=u1A1√2g(pft-p0)γ=dyjdt(Ah-Af)(5)Q4=u2A2√2g(pd-ptt)γ=dyjdtAh(6) 式中:Q3、Q4分别为分油窗口3和4的燃油体积流量;u3、u4分别为分油窗口3和4的流量系数。 由于在初始平衡位置时, 节流口1、2、3和4的初始开度相等且为零, 所以它们开启时的通流面积应相等, 则: A1=A2=A3=A4=r2θ1-r(r-|y|) sinθ1(7) 式中:θ1=arcosr-yr;r为圆孔半径。 2.4 压力位移特性的数学模型 压力位移特性是值计量器的活塞随定压差活门压力变化的特性。定压差活门从中立位置向上运动时, 则由式 (3) 与式 (4) 相除可得: pft=pd-(Ah-Af)2A2h(ptt-p0)(8) 把式 (8) 代入到式 (2) 可得: -[(Ah-Af)3A2h+Ah]ptt+pd(Ah-Af)+(Ah-Af)3A2hp0=mtd2yjdt2+Btdyjdt+kj(yj+c0)(9) 把式 (7) 代入到式 (4) 可得: u2[r2θ1-r(r-y)sinθ1]√2g(ptt-p0)γ=dyjdtAh(10) 由此可得定压差活门向上运动时的压力位移特性的数学模型: {pft=pd-(Ah-Af)2A2h(ptt-p0)-[(Ah-Af)3A2h+Ah]ptt+pd(Ah