发明名称 可调金属切削系统的鲁棒自适应动态面控制方法
摘要 一种可调金属切削系统的鲁棒自适应动态面控制方法,其特点是,包括可调金属切削系统数学模型、磁滞模型、神经网络系统描述和自适应动态面控制器的设计等内容,通过引入的误差转换函数可任意指定控制精度,实现了对切刀切割深度的精确控制;系统模型未完全可知的情况下,通过RBF神经网络实现未知项的逼近;通过估计未知参数向量的范数来代替估计未知参数地向量,极大地减轻了系统的计算负担;采用自适应动态面技术结合误差转换函数与RBF的策略,放宽了对时滞的假设同时保证了跟踪误差和过渡过程能够在任意指定的指标内,消除了反推法的“微分爆炸”问题,保证系统信号半全局一致最终有界,具有方法科学合理,控制精度高,适用等优点。
申请公布号 CN103792849B 申请公布日期 2017.02.08
申请号 CN201410083566.5 申请日期 2014.03.09
申请人 东北电力大学 发明人 张秀宇;王建国;刘平;孔庆新;燕鹏;胡石磊
分类号 G05B13/04(2006.01)I 主分类号 G05B13/04(2006.01)I
代理机构 吉林市达利专利事务所 22102 代理人 陈传林
主权项 可调金属切削系统的鲁棒自适应动态面控制方法,其特征是,它包括以下内容:1)可调金属切削系统数学模型可调金属切削系统的数学模型包括切割系统装置本身模型、用来描述磁滞作动器的磁滞模型及未知延时环节模型,建立可调金属切削系统:<img file="FDA0001071809270000011.GIF" wi="1342" he="70" />其中x表示切割深度的波动部分,也称作模具厚度的偏移量;m,c,k分别代表切割器械的质量、阻尼系数和弹性硬度;k<sub>a</sub>为压电作动器的弹性当量,u为压电作动器的输出,F<sub>Δ</sub>为切割器械的切力变化F<sub>Δ</sub>=h×(x‑μx(t‑τ))+f(x‑μx(t‑τ))其中<img file="FDA0001071809270000012.GIF" wi="195" he="62" /><img file="FDA0001071809270000013.GIF" wi="56" he="39" />h和μ为正常数,τ为未知时间延迟;基于上述对金属切削系统的描述,令x=x<sub>1</sub>,<img file="FDA0001071809270000014.GIF" wi="167" he="63" />则式(1)表述成反馈形式:<img file="FDA0001071809270000015.GIF" wi="157" he="62" /><img file="FDA0001071809270000016.GIF" wi="1373" he="235" />y=x<sub>1</sub>其中<img file="FDA0001071809270000017.GIF" wi="169" he="121" />对于可调金属切削系统(1),作如下假设:假设1:未知延时项<img file="FDA0001071809270000018.GIF" wi="155" he="63" />i=1,…,n,满足如下不等式:<img file="FDA0001071809270000019.GIF" wi="1382" he="135" />其中<img file="FDA00010718092700000110.GIF" wi="110" he="71" />为未知连续函数;假设2:参考信号y<sub>r</sub>光滑、有界,且对于t≥0,<img file="FDA00010718092700000111.GIF" wi="150" he="70" />属于某一紧集;2)磁滞模型采用Prandtl‑Ishlinskii模型来表示磁滞非线性,Prandtl‑Ishlinskii模型是Preisach模型的一个子集,由stop算子和play算子组成,该模型的优势主要是能准确刻画磁滞现象,同时便于实现控制器的设计和实时控制,由于滞环的非对称性,play算子使得磁滞输出w随输入u增加或减小,设C<sub>m</sub>[0,t<sub>E</sub>]为一个单调分段连续的泛函空间,对于任意输入u,定义f<sub>r</sub>:R→R, f<sub>r</sub>(u,w)=max(u‑r,min(u+r,w)),其中r表示滞环阈值且满足r≥0,则play算子F<sub>r</sub>[·]定义为如下形式:F<sub>r</sub>[u](0)=f<sub>r</sub>(u(0),0)F<sub>r</sub>[u](t)=f<sub>r</sub>(u(t),F<sub>r</sub>[u](t<sub>i</sub>))     (4)for t<sub>i</sub>≤t≤t<sub>i+1</sub>and 0≤i≤N‑1其中,0=t<sub>0</sub>&lt;t<sub>1</sub>&lt;…&lt;t<sub>N</sub>=t<sub>E</sub>为[0,t<sub>E</sub>]的一个分割,使得输入u在每一个[t<sub>i</sub>,t<sub>i+1</sub>]内具有单调性,即递增或者递减,PI模型表述为如下形式:<img file="FDA0001071809270000021.GIF" wi="1397" he="103" />其中p(r)是密度函数,满足p(r)≥0且<img file="FDA0001071809270000022.GIF" wi="347" he="70" />密度函数p(r)通过实验辨识得到,由于p(r)随着r的增大而趋于零,为简便,选择D=∞作为式(5)的积分上界,λ为由p(r)确定的未知常数,λ=1.5,<img file="FDA0001071809270000023.GIF" wi="358" he="71" />r∈[0,10]、输入u(t)=7sin(3t)/(1+t),t∈[0,2π]且F<sub>r</sub>[u](0)=0,得到PI磁滞模型;磁滞输出表示为:<img file="FDA0001071809270000024.GIF" wi="1345" he="103" />将(6)式代入(2)得:<img file="FDA0001071809270000025.GIF" wi="158" he="63" /><img file="FDA0001071809270000026.GIF" wi="1366" he="235" />y=x<sub>1</sub>其中:β=bλ    (8)p<sub>g</sub>(r)=bp(r)     (9)且p<sub>g</sub>(r)为未知磁滞密度函数,p<sub>g</sub>(r)和β将在接下来被在线估计,针对上述带有磁滞驱动的非线性系统,控制目的是采用反推法的改进‑动态面方法使得闭环系统保证全局稳定,同时,通过调整控制参数实现系统输出x<sub>1</sub>=y能够跟踪参考信号y<sub>r</sub>;3)神经网络系统描述一般地,神经元网络是一个多输入单输出的系统,其数学表达式为:Y=θ<sup>T</sup>ψ(ξ)    (10)其中,ξ∈R<sup>n</sup>是神经元网络的输入;Y∈R是神经元网络的输出,θ∈R<sup>N</sup>是N维可调参数 向量,这里N是神经元节点数;ψ(ξ):R<sup>n</sup>→R<sup>N</sup>是非线性向量函数,且ψ(ξ)=[ψ<sub>1</sub>(ξ),…,ψ<sub>N</sub>(ξ)]<sup>T</sup>,这里,ψ<sub>k</sub>(ξ)是高斯函数,引理1:对紧集<img file="FDA00010718092700000316.GIF" wi="185" he="66" />上的任意连续函数<img file="FDA00010718092700000317.GIF" wi="323" he="72" />通过选择适当的参数σ<sub>k</sub>和ζ<sub>k</sub>,对于足够大的正整数N,存在形如式(10)的神经网络系统,使得,<img file="FDA0001071809270000031.GIF" wi="478" he="74" /><img file="FDA0001071809270000032.GIF" wi="1270" he="71" />其中,θ<sup>*</sup>是权参数向量θ的最优值,定义δ(ξ)为逼近误差,<img file="FDA0001071809270000033.GIF" wi="1365" he="111" /><img file="FDA00010718092700000318.GIF" wi="1273" he="70" />4)自适应动态面控制器的设计定义系统跟踪误差e:=x<sub>1</sub>‑y<sub>r</sub>   (14)其中y<sub>r</sub>为理想跟踪轨迹,性能函数和误差转换函数将进行如下定义:性能函数<img file="FDA0001071809270000034.GIF" wi="459" he="70" />被定义为一个光滑递减的正函数,使得对于所有t≥0<img file="FDA0001071809270000035.GIF" wi="1403" he="167" />其中0&lt;σ&lt;1且<img file="FDA0001071809270000036.GIF" wi="411" he="68" /><img file="FDA0001071809270000037.GIF" wi="78" he="54" />为系统稳定时跟踪误差所允许的最大值,为将式(15)转换成一等价函数,引入误差转换函数如下:<img file="FDA0001071809270000038.GIF" wi="1253" he="70" />其中,S<sub>1</sub>是由误差转换函数进行转换后的转换误差,Φ(S<sub>1</sub>)是某一光滑、严格单调递增的函数,其反函数具有如下性质:<img file="FDA0001071809270000039.GIF" wi="1366" he="166" />且<img file="FDA00010718092700000310.GIF" wi="1502" he="191" />由式(18)可知,如果<img file="FDA00010718092700000311.GIF" wi="194" he="63" />那么式(17)成立,同时再考虑<img file="FDA00010718092700000312.GIF" wi="187" he="59" />和式(16),当e(0)&gt;0时,<img file="FDA00010718092700000313.GIF" wi="670" he="63" />成立,或者当e(0)&lt;0时,<img file="FDA00010718092700000314.GIF" wi="667" he="62" />成立,从而公式(15)成立,因此,经分析可知,若欲实现给定的性能指标,只需证明<img file="FDA00010718092700000315.GIF" wi="170" he="63" />即可,由Φ(S<sub>1</sub>)的严格单调递增的性质,可得:<img file="FDA0001071809270000041.GIF" wi="1222" he="134" />需要注意的是,将e(0)=0时的情况包含到e(0)&gt;0或e(0)&lt;0中处理,同时,σ不可取零,因为这会使得S<sub>1</sub>(0)无界,第1步:令第1个面误差S<sub>1</sub>为式(19)所定义,考虑式(1)和式(14),有<img file="FDA0001071809270000042.GIF" wi="1333" he="119" /><img file="FDA0001071809270000043.GIF" wi="1270" he="138" />根据Φ和<img file="FDA0001071809270000044.GIF" wi="51" he="39" />的定义,可知Ψ&gt;0,考虑如下二次型方程<img file="FDA0001071809270000045.GIF" wi="1182" he="118" />因此,可得<img file="FDA0001071809270000046.GIF" wi="1518" he="116" />其中x<sub>2d</sub>为待设计的虚拟控制信号,设计虚拟控制律:<img file="FDA0001071809270000047.GIF" wi="1469" he="118" />令x<sub>2d</sub>经一阶低通滤波器获得的新变量z<sub>2</sub>如下:<img file="FDA0001071809270000048.GIF" wi="1350" he="71" />其中,τ<sub>2</sub>为滤波器时间常数;第2步:定义第2个面误差为:S<sub>2</sub>=x<sub>2</sub>‑z<sub>2</sub>    (26)令<img file="FDA0001071809270000049.GIF" wi="1699" he="122" />考虑二次型函数:<img file="FDA00010718092700000410.GIF" wi="1390" he="143" />其中,定义:<img file="FDA00010718092700000411.GIF" wi="1070" he="87" />γ<sub>ζ</sub>,γ<sub>pr</sub>和<img file="FDA00010718092700000412.GIF" wi="59" he="61" />为正设计参 数;<img file="FDA0001071809270000051.GIF" wi="302" he="79" />分别为<img file="FDA0001071809270000052.GIF" wi="439" he="126" /><img file="FDA0001071809270000053.GIF" wi="109" he="77" />的估计值;对于未知延时项<img file="FDA0001071809270000054.GIF" wi="162" he="70" />根据假设1,下述不等式成立:<img file="FDA0001071809270000055.GIF" wi="1446" he="118" />因此,由式(27),可得:<img file="FDA0001071809270000056.GIF" wi="1662" he="263" />其中<img file="FDA0001071809270000057.GIF" wi="157" he="55" />为假设1中的未知连续函数;<img file="FDA0001071809270000058.GIF" wi="150" he="70" />定义为如下形式:<img file="FDA0001071809270000059.GIF" wi="1382" he="135" />由于式(29)中有未知项,因此,需用神经网络作为逼近器,在紧集<img file="FDA00010718092700000510.GIF" wi="59" he="75" />上来逼近未知项:<img file="FDA00010718092700000511.GIF" wi="1518" he="135" />其中<img file="FDA00010718092700000512.GIF" wi="526" he="84" />则有:<img file="FDA00010718092700000513.GIF" wi="1309" he="151" />其中α<sub>2</sub>为正设计参数,根据式(28)‑(32)有:<img file="FDA00010718092700000514.GIF" wi="1933" he="299" />设最终的控制律v为:<img file="FDA00010718092700000515.GIF" wi="1099" he="71" />其中<img file="FDA00010718092700000516.GIF" wi="1381" he="231" /><img file="FDA00010718092700000517.GIF" wi="140" he="78" />和<img file="FDA00010718092700000518.GIF" wi="139" he="70" />的含义由式(29))的调参律设计为如下形式<img file="FDA0001071809270000061.GIF" wi="1222" he="95" /><img file="FDA0001071809270000062.GIF" wi="1294" he="134" /><img file="FDA0001071809270000063.GIF" wi="1398" he="119" />其中σ<sub>ζ</sub>,σ<sub>pr</sub>和<img file="FDA0001071809270000064.GIF" wi="59" he="61" />为正设计参数。
地址 132012 吉林省吉林市船营区长春路169号