-
本文中所采用的测试系统主要是由主控制器、大功率半导体激光器、光学整形模块、椭球面反射镜、热电偶、红外辐射测温模块、信号调理电路、数据采集模块、模/数转换器(analog-to-digital converter,ADC)和反馈控制模块等组成。系统框图如图 1所示。
在此实验系统中,半导体激光器作为被测热电偶的阶跃激励源,其内部的激光二极管可以被高频电流所调制,改变输出功率,在很短的时间内实现热电偶测温端的温度阶跃。所使用的半导体激光器的工作模式有固定功率输出和变功率输出两种模式,在实验中采用变功率输出模式。被测热电偶与红外探测模块分别置于椭球面反射镜的两个共轭焦点处,在热电偶受到激光激励后,椭球面反射镜可以将其所产生的红外辐射汇聚到红外探测模块,由于红外辐射探测器的优点在于其反应速度远远快于热电偶,因此可以将其作为热电偶表面温度的真值,通过ADC输入到控制模块,然后控制模块通过PID控制器或者Tornambe控制器反馈控制激光器的输出功率,以此在热电偶表面达到期望温度。最后信号调理电路和数据采集模块采集热电偶在工作中的时间-电压响应曲线,通过分析热电偶的时间-电压响应曲线,得到热电偶在不同控制器下的时间常数。
热电偶一般被作为1阶系统,在时间常数的测试实验中,不考虑热电偶传热能量和辐射能量的散失,可得如下公式:
$ Q_{1}=\frac{\rho c V \mathrm{d} T_{1}}{\mathrm{~d} t} $
(1) $ Q_{2}=h S\left(T_{2}-T_{1}\right) $
(2) 式中,Q1为热电偶测量端储热量,Q2为热电偶对流传热量,ρ为热电偶测量端密度,c为热电偶测量端比热容,V为热电偶测量端体积,T1为热电偶测量端温度,T2为介质温度,h为对流换热系数,S为测量端换热面积。
设τ=ρcV/(hS),根据(1)式、(2)式整理可得1阶微分方程:
$ \tau\left(\frac{\mathrm{d} T_{1}}{\mathrm{~d} t}\right)+T_{1}=T_{2} $
(3) 1阶微分方程(3)式的通解为${T_1} = C{{\rm{e}}^{ - \frac{1}{\tau }}} + {T_2}$,C为常数。假设t=0时,T1=T1, 0,则可得:
$ \frac{T_{1}-T_{1, 0}}{T_{2}-T_{1, 0}}=1-\mathrm{e}^{-\frac{1}{\tau}} $
(4) 根据(4)式可以得到,当t=τ时,热电偶表面温度达到阶跃温度的63.2%,τ即为热电偶的时间常数。
热电偶时间常数测试的理想输出曲线如图 2所示。横坐标为时间t,纵坐标为温度T,T0为热电偶起始时刻的温度,T3为稳定温度,t1为起始时刻,t2为热电偶表面温度达到阶跃温度的63.2%的时刻。
-
Tornambe控制算法是一种基于主动补偿的控制方法,其基本原理是利用扩张状态观测器,实现对非线性系统各种不确定性的主动估计和动态补偿,进而实现对系统的有效控制。它与其它控制算法的最大区别是不依赖于精确的系统数学模型,核心在于构造观测器来对扰动进行估计和补偿,最终实现对动态系统的反馈控制,因此,对这一控制算法的研究一直受到人们的广泛关注[19-21]。
单变量非线性数学模型表示如下:
$ \left\{ {\begin{array}{*{20}{l}} {\dot x = f(x) + g(x)u,\left( {x \in {{\bf{R}}^n},u \in {\bf{R}}} \right)}\\ {y = h(x),(y \in {\bf{R}})} \end{array}} \right. $
(5) 式中,R是实数集,x表示系统状态,f(x)和g(x)为系统已知的有界非线性函数,u为系统的控制输入,h(x)为系统输出。对于(5)式所示系统,假设相对阶r已知且输出y及y(i)(1≤i≤(r-1))可测,则有微分变换:
$ \left\{\begin{array}{l} z_{i}=\varphi_{i}(x)=L_{f}^{i-1} h(x), (i=1, 2, \cdots, r) \\ w_{i}=\varphi_{i}(x), (i=1, 2, \cdots, n-r) \end{array}\right. $
(6) 式中,φi(x)为中间函数,Lf为所代表的非线性系统,f为光滑函数, 则可以得到下式:
$ \left\{\begin{array}{l} \dot{z}_{1}=z_{2} \\ \dot{z}_{2}=z_{3} \\ \vdots \\ \dot{z}_{r-1}=z_{r} \\ \dot{z}_{r}=a(\boldsymbol{z}, \boldsymbol{w})+b(\boldsymbol{z}, \boldsymbol{w}) u \\ \dot{\boldsymbol{w}}=c(\boldsymbol{z}, \boldsymbol{w}) \\ y=z_{1} \end{array}\right. $
(7) 式中,z=(z1, z2, …, zr)T;w=(w1, w2, …, wn-r)T;a(z, w),b(z, w)和c(z, w)可由f(x),g(x),h(x)和微分变换函数得到。
假设$\dot{\boldsymbol{w}}=c(\boldsymbol{z}, \boldsymbol{w})$是动态稳定的,则有:
$ y^{(r)}+h_{r-1} y^{(r-1)}+\cdots+h_{1} y^{(1)}+h_{0} y=r $
(8) 式中,h0,h1,…,hr-1为系统的动力学参量,在保证系统稳定性的前提要求下,hi(i=0, 1, 2, …, r-1)的选取必须保证系统的闭环极点位于s的开左半平面,即特征方程sr+hr-1sr-1+…+h2s2+h1s+h0=0的根位于s的开左半平面。
在满足上述条件下,将扩张状态定义为:
$ d(\boldsymbol{z}, \boldsymbol{w}, u)=a(\boldsymbol{z}, \boldsymbol{w})+[b(\boldsymbol{z}, \boldsymbol{w})-1] u $
(9) 则可以用${\mathit{\boldsymbol{\dot z}}_r} = d(\mathit{\boldsymbol{z}}, \mathit{\boldsymbol{w}}, u) + u$来作为系统的动态性能方程。如果a(z, w)和b(z, w)均可测,为了实现系统的线性化,可以采用精确反馈线性化方法将d(z, w, u)抵消掉,得到:
$ u=-h_{0} z_{1}-h_{1} z_{2}-\cdots-h_{r-1} z_{r}-d(\boldsymbol{z}, \boldsymbol{w}, u) $
(10) 但是在实际系统中,不可避免地会存在大量的内外扰动,状态变量a(z, w)和b(z, w)经常会因为无法得到精确值而降低系统的控制性能。为了消除这种不确定性因素,利用包含积分环节的$\hat d$代替(9)式中的d(z, w, u),这样就可以利用扩张状态观测器来对系统的总体不确定扩张状态d(z, w, u)给予估计和补偿,得到如下的Tornambe控制器:
$ \left\{\begin{array}{l} u=-h_{0} z_{1}-h_{1} z_{2}-\cdots-h_{r-1} z_{r}-\hat{d} \\ \hat{d}=\xi+\sum\limits_{i=0}^{r-1} k_{i} z_{i+1} \\ \dot{\xi}=-k_{r-1} \xi-k_{r-1} \sum\limits_{i=0}^{r-1} k_{i} z_{i+1}-\sum\limits_{i=0}^{i-2} k_{i} z_{i+2}-k_{r-1} u \end{array}\right. $
(11) 式中,k0, k1, …, kr-2为任意常数,kr-1=sgn(b(z, w))μ,μ为一正值常数,μ的选取决定了系统的稳定性;ki(i=0, 1, …, r-2)为任意常数,主要影响系统的控制性能。如果存在常数μ*>0,当μ>μ*时,方程组(5)式和(11)式所示闭环系统是渐进稳定的。
取r=2,得到2阶Tornambe控制器的控制方程为:
$ \left\{\begin{array}{l} u=-h_{0}\left(y-y_{r}\right)-h_{1}\left(y-y_{r}\right)^{(1)}-\hat{d} \\ \hat{d}=\xi+k_{0}\left(y-y_{r}\right)+k_{1}\left(y-y_{r}\right)^{(1)} \\ \dot{\xi}=-k_{1} \xi-k_{1} k_{0}\left(y-y_{r}\right)- \\ \qquad k_{1}{}^{2}\left(y-y_{r}\right)^{(1)}-k_{0}\left(y-y_{r}\right)^{(1)}-k_{0} u \end{array}\right. $
(12) 式中,y为系统输出,yr为系统给定,$\hat d$为扩张状态观测函数,用来实现对系统内部扰动的实时估计和补偿,ξ为中间变量,h0和h1为可调参量,决定系统的灵敏度,k0和k1为可调参量,决定系统的抗干扰能力。
-
在温度25℃、湿度55%、标准大气压环境下,以上面所提到的时间常数测试系统为实验基础,所采用的半导体激光器的光谱曲线图如图 5所示,工作参量如表 1所示。分别采用Tornambe控制器和PID控制器对CO1-K型热电偶进行测试,数据采集卡采集得到热电偶的时间-电压输出电压响应曲线如图 6所示。
Table 1. Operating parameters of semiconductor laser
wavelength bandwidth instantaneous power modulation of time load current 980nm 3.9nm 594.2W 10μs 10A~200A 在LabVIEW中可以很直观地得到在这两种不同的控制方法下,热电偶的电压信号上升时间、超调量和时间常数,分析结果如表 2所示。
Table 2. CO1-K thermocouple measurement results
controller rise time/ ms overshoot/ % time constant/ ms stable voltage/ V PID 830.6 0 456.2 6.93 Tornambe 540.3 6.47 284.6 6.95 通过上述实验数据对比可以发现:在同一硬件系统下,CO1-K型热电偶在传统PID控制下所测得的时间常数为456.2ms,在Tornambe控制下的时间常数测试结果为284.6ms; 从实验结果中可以发现,相较于传统的PID控制算法,Tornambe控制下的系统达到稳态所需要的时间更短,阶跃温升源的上升时间更短,因此在实验中,热电偶受到激励后的响应时间更短,速度更快。在控制器的参量设置上,PID控制器中参量设置合理,所以超调量为0,Tornambe控制器中h参量设置较大,使得系统可以更快达到预期设定目标。虽然在Tornambe控制下出现了6.47%的超调量,但是在允许范围内,对实验结果没有太大影响。
基于半导体激光器的Tornambe控制算法研究
Research on Tornambe control algorithm based on semiconductor laser
-
摘要: 为了给热电偶时间常数测试中提供更加准确的阶跃温升信号,优化控制效果,提高热电偶时间常数测量的准确性,采用Tornambe控制器来反馈控制半导体激光器的输出功率,在MATLAB下的Simulink模块中实现整个系统的搭建,并对系统进行同一输入信号下的仿真研究,通过实验来比较2阶Tornambe控制器与比例-积分-微分(PID)控制在系统运行过程中的控制效果。结果表明,采用PID控制器测得CO1-K型热电偶的时间常数为456.2ms,而采用Tornambe控制器测得的热电偶时间常数为284.6ms。2阶Tornambe控制器能够有效缩短热电偶达到平衡温度的时间,且在控制器结构参量整定上也更加简便,具有较强的实用价值。
-
关键词:
- 测量与计量 /
- 时间常数测试 /
- Tornambe控制 /
- 热电偶 /
- 比例-积分-微分
Abstract: In order to provide more accurate step temperature rise signal in thermocouple time constant testing, optimize the control effect, and improve the accuracy of thermocouple time constant measurement, Tornambe controller was used to feedback control the output power of semiconductor laser. The system was constructed in the Matlab Simulink module, and the simulation was carried out under the same input signal. The control effects of the second-order Tornambe controller and the proportional-integral-differential (PID) were compared in the running process of system. The results show that the time constant of CO1-K thermocouple measured by PID controller is 456.2ms, while that measured by Tornambe controller is 284.6ms. The second-order Tornambe controller can effectively shorten the time of the thermocouple temperature balance. Moreover, it is more convenient to set the parameters of the controller structure, which has strong practical value. -
Table 1. Operating parameters of semiconductor laser
wavelength bandwidth instantaneous power modulation of time load current 980nm 3.9nm 594.2W 10μs 10A~200A Table 2. CO1-K thermocouple measurement results
controller rise time/ ms overshoot/ % time constant/ ms stable voltage/ V PID 830.6 0 456.2 6.93 Tornambe 540.3 6.47 284.6 6.95 -