数字科氏质量流量计正负阶跃交替激励启振方法

1、数字科氏质量流量计正负阶跃交替激励启振方法第31卷 第1期 仪 器 仪 表 学 报 Vol.31 No. 1 2010年1月 Chinese Journal of Scientific Instrument Jan. 2010 数字科氏质量流量计正负阶跃 *交替激励启振方法 11,2111李 苗，徐科军，侯其立，朱永强，李 叶 (1 合肥工业大学电气与自动化工程学院 合肥 230009； 2 工业自动化安徽省工程技术研究中心 合肥 230009) 摘 要：数字驱动由于其控制算法的灵活性可以有效地提高科氏质量流量计的工作性能。针对数字驱动的启振问题，本文基于模型分析和实验仿真，提出了正负阶跃交替

2、激励启振方法。该方法通过跟踪检测传感器自激输出信号相位进行正负阶跃交替激励，使流量管的振动幅值持续、可靠地快速增大，进而结合非线性幅值控制算法和频率估计，可使流量管启振时间大幅度地缩短，克服了传统启振方法的局限性。针对高准型CNG050型科氏质量流量传感器，研制了基于DSP的二次仪表，并进行实验研究，验证了该启振方法的可靠性和优越性。 关键词：科氏流量计；数字驱动；正负阶跃信号；启振；两相流；批料流 中图分类号：TB931 文献标识码：A 国家标准学科分类代码：510.8040 Startup method of digital Coriolis mass flowmeter using al

3、ternating exciting of positive-negative step signal 11,2111Li Miao，Xu Kejun，Hou Qili，Zhu Yongqiang，Li Ye (1. School of Electrical and Automation Engineering, Hefei University of Technology, Hefei 230009, China； 2. Engineering Technology Research Center of Industrial Automation, Anhui Province, Hefei

4、 230009, China) Abstract：Digital driving can improve the performance of Coriolis mass flowmeter effectively due to the flexibility of its control algorithm. To solve the startup problem, this paper presents a startup method using alternating exciting of positive-negative step signal based on model a

5、nalysis and experimental simulation. The method tracks the sensor output signal phase and alternately stimulates the flow tube with positive-negative step signal, which can increase the vibration amplitude rapidly and reliably. Combined with non-linear amplitude control and frequency estimation, the

6、 startup time of the flow tube can be shortened significantly, and the limitations of traditional startup methods are overcome. For Micromotion-CNG050 type Coriolis mass flow sensors, a secondary instrument was developed based on DSP, and experimental studies were made to verify the reliability and

7、superiority of the method. Key words：Coriolis flowmeter; digital driving; positive-negative step signal; startup; two-phase flow; batch flow 1- 2个数量级以上，而实际测量要求流量管以固有频率稳1 引 言 幅振动，因而，在批料流或两相流发生时，驱动系统如何进行频率、相位和幅值的快速跟踪控制至关重要。流随着科技和工业生产的发展，特别是过程工业的发量管的启振过程包含了这些参量的跟踪控制，是一个由展，需要采用科氏流量计测量批料流和气液两相流等复静止到稳定工作的

8、过程，能完全反映出驱动系统的跟踪杂流体。但是，由于市场现有型号科式流量管空满管固控制性能，因而研究流量管的快速启振具有重要意义。有频率大多相差近10 Hz，两相流发生时阻尼比变化有两收稿日期：2009-07 Received Date：2009-07 *基金项目：“863”计划项目（2009AA04Z128）、合肥工业大学专利转化基金（2008HGCZ0523）资助 第1期 李 苗 等：数字科氏质量流量计正负阶跃交替激励启振方法 173 3-4现有的模拟驱动启振时间长，且无法维持流量管在复式中：=2f, f为驱动频率，A为正弦驱动的幅值，为ddd杂流体工况下稳幅振动，甚至可能造成流量管停振。为

9、方便计算，设初始相位=90。 05-9此，人们研究数字驱动方法。数字驱动即在控制回路对式(2)进行拉普拉斯变换，再与式(1)相乘，通过反中引入数字系统环节，通过数字合成的方法输出信号驱拉普拉斯变换可得流量管的正弦响应函数为： 2Ak(bcost?asint)动流量管，控制方法灵活。目前，国内的研究重点还是11d1d y(t)=+22 (3) b+a11放在提高科氏流量计的测量精度方面，很少研究数字驱?t2Ak(acost+bsint)e22m2m动问题；国外虽然研究了数字驱动的相关内容，如非线2式中：=1?为流量管固有频率，=为与10-11mnn性幅值控制，数字科氏流量计驱动系统的FPGA实阻

10、尼比及自然频率有关的指数衰减因子。 12现方法，但对数字驱动的启振过程及实现细节未做详可见，传感器输出含有两个频率分量：稳定的驱动13-15细披露。为此，本文采用数字驱动方法研究启振过信号频率分量和衰减的固有频率分量。在开始时，由于程并通过比较提出一种基于相位跟踪的正负阶跃启振方无法得到准确的固有频率，激振信号频率只能给定一个法，结合频率估计和非线性幅值控制，使流量管快速平大概值，因此，传感器输出信号中必然含有两个频率分稳启振。 量。这时固有频率的频率估计难度较大，对算法要求较高。为了准确估计出固有频率，当传感器输出信号的幅2 启振方法 值等于或者大于a值时，使流量管进入零驱动模式，则2.1

11、数字驱动 传感器输出信号中只含衰减的固有频率分量，又由于阻科氏流量计的振动体系为无限自由度受迫振动体尼比很小，从而衰减因子很小，流量管就会在较长时间系，有无限多个主振型。激振系统一般都采用第一主振内以固有频率作近似稳幅振荡。 型。因而，这里采用一维有阻尼受迫振动系统来描述流由于频率估计算法需要传感器信号的幅值达到一定量管模型，即其传递函数可表示为： 的数值，当初始正弦自激信号频率与固有频率相差较远KSs时，可能需要较长时间才能达到或始终达不到阈值，从G(S)= (1) 22S+2S+nn而导致启振很慢甚至失败，因此用正弦波作为自激信号式中：K为与流量管特性有关的参数，、分别为阻尼sn具有一定的

12、局限性。用三角波和方波作为自激信号也存比和自然振荡角频率。 在同样的问题。虽然这些波形是由多种频率成份组成的，数字驱动原理框图如图1所示。驱动起始阶段，由但是，均为基频的奇数倍。若其基频与固有频率相差甚驱动模块产生系统自激信号激振流量管；当传感器输出远或与之不成奇数倍关系，则流量管同样无法可靠地启信号幅值u达到频率估计要求幅值a后，停止激励，进振。随机波作为自激信号时，其工作原理与模拟驱动相入零驱动模式；由频率估计算法获得准确固有频率，合同。因此，启振时间较长，不利于提升数字驱动的启振成正弦驱动信号；结合非线性幅值控制算法，使流量管性能。流量管固有频率虽可通过实验建模获知，但很难迅速启振至期望

13、幅值，进入正常工作阶段。 用于批料流或两相流情况。 2.3 正负阶跃交替激励启振 从分析流量管的阶跃响应入手，提出正负阶跃交替 激励启振法。当输入信号为正阶跃信号时，其拉普拉斯图1 数字驱动原理框图 变换为： Fig.1 Schematic of digital driving A0X(S)= (4) S由图1可知，自激信号的选择以及频率估计算法和式中：A为阶跃幅值，将该式与式(1)相乘，再作拉普拉0幅值控制都将影响着流量管的快速平稳启振。又由仿真斯反变换，即可得流量管的阶跃响应函数为： 研究及实验可知，自激驱动使流量管输出增大至要求幅AK?t0sy(t)=sin(t)e (5) m值a的时间

14、是制约启动速度的关键，因而需研究如何产m生自激信号。 可见，阶跃信号能激励出只含固有频率的输出信号，2.2 正弦波启振 幅值相同的正负阶跃对流量管的作用等值反相。据此，当初始激励为正弦波时，即： 搭建Simulink仿真图，如图2所示。 x(t)=Asin(t+) (2) d0 174 仪 器 仪 表 学 报 第31卷 输出信号的频率，它由两个格型滤波器级联而成。整个格型IIR陷波器的传递函数为： ?1?21+2kz+z0H(z)= (6) ?1?21+k(1+)z+z0格型自适应算法计算量小，收敛速度快，并且收敛后数值稳定。该算法要求被测信号频率成份单一，而实际中不能完全滤除噪声，故信号幅值

15、必须达到一定要求才能满足其估计精度。 图2 流量管正、负阶跃响应 Fig.2 Positive-negative step response of a flow tube 由图2可知，单位阶跃激励，流量管的输出幅值在 图4 格型IIR陷波器 0.07 mV左右，而频率估计算法要求幅值在3 mV以上，且Fig.4 Lattice IIR notch filter 实际工作环境中存在噪声，因而很难估算出准确的固有频率。但是，若在合适的位置，用正负阶跃交替激励，则流2.5 非线性幅值控制方法 量管输出将不断加强。通过MATLAB仿真可知，当流量管非线性幅值控制框图如图5所示。在外循环中，首先输出信号

16、进入(90+k360)(+90+k360)范围内时施加通过计算，得到传感器输出信号的幅值A(t)；对给定幅值负阶跃，可使流量管输出信号得到加强；当信号进入A和A(t)取自然对数，分别得到a和a(t)，求其差值；然00(+90+k360)(+270+k360)范围内时施加正阶跃同样可后，通过控制器C(s)求出驱动增益K(t)。取自然对数单元0使信号得到加强；在信号到达(0+k360)时施加负阶跃、可以使求的差值比单纯相减来得大，这样就加快了控制的到达(180+k360)时施加正阶跃可使信号得到最大幅度加速度。为了得到零稳态误差，控制器C(s)采用PI(比例积强。考虑实际测量中的可操作性，由幅值检

17、测来实现相位分)控制器。该控制器优于模拟驱动中简单的比例控制。 判定，即通过过零检测实现(0+k360)和(180+k360)相位点检测。在实际中，由于环境噪声的存在，过零点很难找准，而前面所述的相位段的确定则要容易的多。因此，可在零点附近设置一个滞环带b。滞环带b的设定可避免在零点由于噪声叠加而导致反复误动作的情况。b值大小可根据实际噪声大小设定，略大于噪声即可。当输出信号大于b时施加负阶跃，而当其小于b时施加正阶跃。 根据上述原理，用MATLAB仿真流量管启振情况， 并在输出信号上叠加幅值大小为1 mV的随机噪声，仿图5 非线性幅值控制框图 真结果如图3所示。在仿真中，改变流量管的阻尼比、

18、Fig.5 Schematic of Non-linear control 自然频率等参数时，正负阶跃启振法均能快速激励出符3 系统实现 合频率估计算法要求的输出信号。 针对高准CNG050型科氏流量计一次仪表，选用片上资源丰富、带浮点协处理器、信号处理能力强的TI公司DSP芯片TMS320F28335作为核心处理器，研制了二次仪表，总体硬件框图如图6所示。由驱动模块、信号采集模块、人机接口模块、信号输出模块及通讯模块等组成。其中，驱动实现主要考虑两种波形合成方法，即 直接数字合成和递归算法合成。由于递归算法合成输出图3 正负阶跃启振法仿真图 每更新一点都需处理器计算完成，因而处理工作量比较F

19、ig.3 Simulation diagram of positive-negative step startup 大。考虑到二次仪表的负荷能力，系统中采用直接数字2.4 频率估计算法 合成的方法，由DSP控制驱动模块中的DDS器件实现。本文采用如图4所示的格型IIR陷波器估计传感器 第1期 李 苗 等：数字科氏质量流量计正负阶跃交替激励启振方法 175 DDS器件只需频率和相角参数即能合成所需波形，结合启振过程软件操作流程为：1)初始阶段跟踪传感器MDAC进行幅值控制，进而实现驱动信号的更新。 信号相位用正负阶跃激励流量管；由于实际噪声在0.1 mV左右，因此，设定滞环带b值为0.15 mV

20、；2)频率估计算法要求幅值为3 mV以上，且相对稳定，所以，当检测幅值达到5 mV以上时，进入零驱动模式；3)由频率估计算法计算固有频率；4)根据固有频率，结合非线性幅值控制算法更新正弦驱动，完成快速启振。 4 实验结果 根据上述实验系统，通过示波器采得数字驱动正负图6 基于DDS的全数字驱动框图 阶跃启振实验曲线波形如图8所示（为方便图示，驱动Fig.6 Schematic of all-digital driving based on DDS 信号及传感器信号曲线做了偏置和增益处理）。图中，测试启振时间电平信号的高电平时间为启振时间。由图可基于该系统，我们进行了正负阶跃启振实验。实验知，数

21、字驱动基于其合成信号的灵活性，初始阶段就可中，以启振开始到输出信号幅值达到期望幅值80%时间根据流量管信号特点产生正负阶跃交替激励信号，经过段作为启振时间，采用TektronicDPO4054型号示波器同约0.07 s的激励时间即可获得符合频率估计算法要求的时采集流量管传感器输出信号、驱动信号及测试启振时输出信号；获得固有频率后结合非线性幅值控制算法，间电平信号。实验步骤为：1)将传感器两端接上钢丝软可使流量管在4.2 s内成功启振至期望幅值，并稳幅振荡。管水平放置，充满水后关闭上游阀门，保证传感器内满为验证数字驱动正负阶跃激励启振方法的优越性，我们管无流量；2)将传感器与二次仪表的驱动输出及

22、采集输同时研制了模拟驱动系统，并基于同样的实验条件进行入相连；3)将示波器接上前述三路信号输出；4)示波器上了启振实验，实验结果曲线如图9所示。由图可知，因电准备好后，给系统上电，并同时由示波器采集三路信其有限的驱动增益和简单的增益控制算法致使流量管启号；5)启振完成后由示波器保存所采信号，并给系统断振缓慢，模拟驱动启振时间长达22 s，比数字驱动时间电。系统软件程序流程图如图7所示。 长达5倍左右。 图8 正负阶跃启振实验曲线图 Fig.8 Experimental curve of positive-negative step startup 图9 模拟驱动实验结果曲线图 图7 启振操作流

23、程图 Fig.9 Experimental curve of analog driving Fig.7 Flow chart of startup 176 仪 器 仪 表 学 报 第31卷 可见，数字驱动基于其控制算法的灵活性，采用正 XU K J, XUN W F . Research on analog driving methods 负阶跃交替激励启振法可使流量管启振时间大幅缩短，of coriolis mass flow meterJ. Acta Metrologica Sinica. 从而提升了跟踪控制性能和流量管测量效率。快速启振2005, 26(2): 149-154. 方法结合

24、先进的幅值控制等算法，可使流量管在各种工5 张瀚, 徐科军. 新型数字式科氏质量流量变送器J. 况下稳定工作，进而为扩大科式流量计的应用范围奠定自动化仪表. 2005, 26(1): 25-28. 了基础。 ZHANG H, XU K J . A new digital coriolis mass flow transmitterJ. Process Automation Instrumentation. 2005, 5 结 论 26(1): 25-28. 6 HENRY M P, CLARKE D W, VIGNOS J H. Digital 通过上述研究和实验，得出以下结论： Flowme

25、terP. U. S Patent: No.2004/6754594B2. 1)正负阶跃启振方法能快速可靠地启振流量管，实7 HENRY M. Self-validating digital Coriolis mass flow 现过程简单。其中，滞环带的设置避免了在零点由于噪meterJ. Computing & Control Engineering Journal. 声叠加而导致反复误动作的情况，而零驱动模式的引入2001, 10: 219-226. 保证了频率估计算法的精度； 8 HENRY M P, CLARKE D W, VIGNOS J H, et al. Digital

26、 2)实验结果表明，正负阶跃启振的数字驱动能使流FlowmeterP. U. S Patent: No.2002/003 8186A1. 量管在4.2 s内完成启振，与模拟驱动相比，大幅提升了9 徐科军, 张瀚. 一种科氏流量计的数字信号处理与驱流量管启振的可靠性和快速性； 动方法研究J. 计量学报. 2004, 25(4): 339-343, 379. 3)若用正弦波启振，则要求驱动信号频率与流量管XU K J , ZHANG H. A study on the digital signal proc-固有频率之差在一定范围内，以保证能激励出有效信号。essing and driving m

27、ethod for coriolis mass flowmetersJ. 即使通过实验建模获知固有频率也很难用于批料流或两Acta Metrologica Sinica. 2004, 25(4): 339- 343, 379. 相流情况； 10 CLARKE D W. Non-linear control of the oscillation am-4)理论分析表明，流量管固有频率变化时，频率估plitude of a coriolis mass-flow meterJ. European Jour-计算法能快速跟踪计算并更新驱动输出，维持流量管稳nal of Control. 1998,4:

28、196-207. 定工作。即使瞬间停振，也能重新快速启振至稳定状态。11 李祥刚，徐科军.科氏质量流量管非线性幅值控制方法从而大幅提升了驱动系统的跟踪控制性能，为科式流量研究J.电子测量与仪器学报，2009，6：82-86. 计用于批料流和气液两相流的测量奠定了基础。 LI X G，XU K J.Research on non-linear amplitude con-本文研制的系统已经申请国家发明专利。 trol method of coriolis mass flow-tubeJ.Journal of Elec-参考文献 tronic Measurement and Instrument,

29、2009,6:82-86. 12 MAYELA Z, HENRY M P. An FPGA implementation of 1 MANUS H, TOMBS M, DUTA M, et al. Two-phase flow metering of heavy oil using a coriolis mass flow meter: Aa digital coriolis mass flow metering drive systemJ. IEEE Trans Ind. Electron., 2008, 55(7): 2820-2831. case studyJ. Flow Measure

30、ment and Instrumentation. 13 HENRY M P, ZAMORA M E. Startup and Operational 2005, 17: 399-413. 2 HENRY Mtechniques for a digital flowmeterP. U. S Patent: P, CLARKE D W, ARCHER N, et al. A No.2005/6950760B2. self-validating digital Coriolis mass-flow meter: an over-6view14 HENRY M P, ZAMORA M E. Star

31、tup and operational J. Control Engineering Practice. 2000, 8(5): 487-506. 3 徐科军, 于翠欣, 苏建徽, 等. 科里奥利质量流量计激振techniques for a digital flowmeterP. U. S Patent： No. 电路的研制2006/7146280B2. J. 合肥工业大学学报. 2002, 23(1): 37-40. XU K J, YU C X, SU J H,15 马金尼斯R L. 科里奥利流量计中用于驱动控制的初 et al. Development of excita-始化算法P.