发明名称 基于区间二型模糊逻辑的微细电火花两阶模糊控制方法
摘要 一种基于区间二型模糊逻辑的微细电火花两阶模糊控制方法,属于微细电火花加工技术领域,涉及一种对微细电火花加工放电状态的检测和控制方法。其特征是首先将采集的电压和电流经过第一阶区间二型模糊系统,得到采样点放电状态矢量;对一个周期内的放电状态矢量进行统计,得到该周期内各放电状态的比率;将其作为输入量再经过第二阶区间二型模糊系统,最终输出微细电极的进给速度。本发明的效果和益处是所提出的控制方法有效克服了平均电压控制方法过于简单化、调节量单一化等缺点及传统一型模糊逻辑在处理不确定性方面的局限性,保证了控制系统和加工过程的稳定性和准确性,显著提高了加工效率,是一种非常适用于微细电火花加工的控制方法。
申请公布号 CN102069245A 申请公布日期 2011.05.25
申请号 CN201010544725.9 申请日期 2010.11.08
申请人 大连理工大学 发明人 贾振元;张玲瑄;刘巍;李安超;魏力
分类号 B23H7/20(2006.01)I;B23H9/14(2006.01)I;G05B13/04(2006.01)I 主分类号 B23H7/20(2006.01)I
代理机构 大连理工大学专利中心 21200 代理人 侯明远
主权项 1.一种基于区间二型模糊逻辑的微细电火花两阶模糊控制方法,其特征在于,将传统一型模糊系统扩展为二型模糊系统,并建立两阶模糊控制器,第一阶模糊控制系统用来准确判别采样点的放电状态,第二阶模糊控制系统用来输出电机的进给速度;首先将加工过程中实时采集到的间隙电压和间隙电流作为第一阶模糊控制系统的输入,经过区间二型模糊逻辑的模糊化和规则表进行推理,得到包含四种放电状态概率的采样点放电状态矢量,通过对一个分析周期内的所有放电状态矢量进行统计,得到该分析周期内各放电状态的比率,并将其作为第二阶模糊系统的输入,再经过完整的区间二型模糊系统,最终输出电机的运动速度,将该速度乘以速度系数,即可得到微细电极的进给速度;基于区间二型模糊逻辑的微细电火花加工两阶模糊控制方法具体步骤如下:(1)第一阶控制系统输入输出隶属度函数的建立将加工过程中实时采集到的间隙电压和间隙电流,除以数据采集卡的增益η,得到值在[0,5]之间的标准化电压U<sub>i</sub>和电流I<sub>i</sub>,作为第一阶模糊控制系统的输入;结合加工经验和神经网络优化方法,初步建立一型模糊输入输出隶属度函数,并制定相应模糊规则,然后将输入隶属度函数进一步扩展为区间二型模糊隶属度函数:针对每一个隶属度函数μ<sup>i</sup>,选取容量为N的样本,样本数据记为z<sub>i</sub>=[U<sub>i</sub>,I<sub>i</sub>,S<sub>Pi</sub>],其中S<sub>Pi</sub>为经过模糊推理得到的放电状态值;将样本分成k组,并保证每组中的样本数量大于1且||z<sub>i</sub>-z<sub>j</sub>||≤δ,δ为根据实际样本情况选取的一个极小值;找出每组样本中隶属度最大和最小的两个数据,并将其隶属度值记为max<sup>l</sup>(μ)和min<sup>l</sup>(μ),l=1,2,...,k;最后以原隶属度函数μ<sup>i</sup>为中心,向两侧分别进行扩展,扩展量为:<maths num="0001"><![CDATA[<math><mrow><mo>|</mo><mi>&Delta;</mi><msup><mi>&mu;</mi><mi>i</mi></msup><mo>|</mo><mo>=</mo><mfrac><mn>1</mn><mi>k</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><mrow><mo>(</mo><msup><mi>max</mi><mi>l</mi></msup><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>-</mo><msup><mi>min</mi><mi>l</mi></msup><mrow><mo>(</mo><mi>&mu;</mi><mo>)</mo></mrow><mo>)</mo></mrow><mo>/</mo><mn>2</mn></mrow></math>]]></maths>扩展后,隶属度函数形成了上下两个边界,其中上边界的函数表达式为<img file="FSA00000346399900022.GIF" wi="292" he="61" />下边界的函数表达式为<u>μ</u><sup>i</sup>=μ<sup>i</sup>-|Δμ<sup>i</sup>|,并且当<img file="FSA00000346399900023.GIF" wi="319" he="61" />时,取<u>μ</u><sup>i</sup>=0,当<u>μ</u><sup>i</sup>=μ<sup>i</sup>-|Δμ<sup>i</sup>|<0时,取<img file="FSA00000346399900024.GIF" wi="138" he="61" />(2)采样点放电状态辨识方法使用输入变量隶属度的最小化运算,对第一阶模糊推理的前提隶属度进行计算;对于模糊规则l,其前提隶属度F<sup>l</sup>为一个区间<img file="FSA00000346399900025.GIF" wi="254" he="70" /><u>f</u><sup><u>l</u></sup>和<img file="FSA00000346399900026.GIF" wi="48" he="60" />为其上下边界,因此一阶模糊系统的输出是一个区间<img file="FSA00000346399900027.GIF" wi="188" he="73" />设输出论域Y<sub>i</sub>内共有n条隶属度函数,其重心的最大及最小值分别为<u>b</u><sup><u>l</u></sup>和<img file="FSA00000346399900028.GIF" wi="58" he="62" />则<u>S</u><sub><u>p</u></sub>和<img file="FSA00000346399900029.GIF" wi="46" he="64" />可通过下面两式计算<maths num="0002"><![CDATA[<math><mrow><munder><msub><mi>S</mi><mi>p</mi></msub><mo>&OverBar;</mo></munder><mo>=</mo><mfrac><mrow><munder><msup><mi>b</mi><mn>1</mn></msup><mo>&OverBar;</mo></munder><munder><mo>&Integral;</mo><msub><mi>Y</mi><mn>1</mn></msub></munder><mover><msup><mi>f</mi><mn>1</mn></msup><mo>&OverBar;</mo></mover><msub><mi>dy</mi><mn>1</mn></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>2</mn></mrow><mi>n</mi></munderover><munder><msup><mi>b</mi><mi>l</mi></msup><mo>&OverBar;</mo></munder><munder><mo>&Integral;</mo><msub><mi>Y</mi><mi>l</mi></msub></munder><munder><msup><mi>f</mi><mi>l</mi></msup><mo>&OverBar;</mo></munder><msub><mi>dy</mi><mi>l</mi></msub></mrow><mrow><munder><mo>&Integral;</mo><msub><mi>Y</mi><mn>1</mn></msub></munder><mover><msup><mi>f</mi><mn>1</mn></msup><mo>&OverBar;</mo></mover><msub><mi>dy</mi><mn>1</mn></msub><mo>+</mo><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>2</mn></mrow><mi>n</mi></munderover><munder><mo>&Integral;</mo><msub><mi>Y</mi><mi>l</mi></msub></munder><munder><msup><mi>f</mi><mi>l</mi></msup><mo>&OverBar;</mo></munder><msub><mi>dy</mi><mi>l</mi></msub></mrow></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0003"><![CDATA[<math><mrow><mover><msub><mi>S</mi><mi>p</mi></msub><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><mover><msup><mi>b</mi><mi>l</mi></msup><mo>&OverBar;</mo></mover><munder><mo>&Integral;</mo><msub><mi>Y</mi><mi>l</mi></msub></munder><mover><mi>f</mi><mo>&OverBar;</mo></mover><msub><mi>dy</mi><mi>l</mi></msub><mo>+</mo><mover><msup><mi>b</mi><mi>n</mi></msup><mo>&OverBar;</mo></mover><munder><mo>&Integral;</mo><msub><mi>Y</mi><mi>n</mi></msub></munder><mover><mi>f</mi><mo>&OverBar;</mo></mover><msub><mi>dy</mi><mi>n</mi></msub></mrow><mrow><munderover><mi>&Sigma;</mi><mrow><mi>l</mi><mo>=</mo><mn>1</mn></mrow><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></munderover><munder><mo>&Integral;</mo><msub><mi>Y</mi><mi>l</mi></msub></munder><mover><mi>f</mi><mo>&OverBar;</mo></mover><msub><mi>dy</mi><mi>l</mi></msub><mo>+</mo><munder><mo>&Integral;</mo><msub><mi>Y</mi><mi>n</mi></msub></munder><mover><mi>f</mi><mo>&OverBar;</mo></mover><msub><mi>dy</mi><mi>n</mi></msub></mrow></mfrac><mo>,</mo></mrow></math>]]></maths>如果输出隶属度函数仍保持一型,则<img file="FSA000003463999000212.GIF" wi="207" he="63" />利用加权求和方法可得到最终输出的放电状态值为<img file="FSA000003463999000213.GIF" wi="385" he="73" />0<λ<1;(3)采样点放电状态矢量映射及统计方法使用四维单位矢量p来描述采样点的四种放电状态,将采样点输出论域按放电状态取值范围划分成4个类,记为<img file="FSA000003463999000214.GIF" wi="260" he="69" />其中i=1,2,3,4;保留完整的采样点放电状态信息,设计映射规则为:如果S<sub>p</sub>∈M<sub>i</sub>∩M<sub>i+1</sub>,则<img file="FSA000003463999000215.GIF" wi="250" he="141" />且<img file="FSA000003463999000216.GIF" wi="296" he="130" />其余元素为0;如果s<sub>p</sub>仅包含于M<sub>i</sub>,则p<sub>i</sub>=1,其余元素为0;由此完成采样点放电状态矢量映射过程;假设一个分析周期内采样点数量为N,则该分析周期内的短路率ζ<sub>SH</sub>,火花率ζ<sub>SP</sub>,脉间率ζ<sub>OFF</sub>和开路率ζ<sub>O</sub>可按下式得到<maths num="0004"><![CDATA[<math><mrow><msub><mi>&zeta;</mi><mi>SH</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>p</mi><mn>1</mn><mi>j</mi></msubsup></mrow><mi>N</mi></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0005"><![CDATA[<math><mrow><msub><mi>&zeta;</mi><mi>SP</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>p</mi><mn>2</mn><mi>j</mi></msubsup></mrow><mi>N</mi></mfrac><mo>,</mo></mrow></math>]]></maths><maths num="0006"><![CDATA[<math><mrow><msub><mi>&zeta;</mi><mi>OFF</mi></msub><mo>=</mo><mfrac><mrow><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>p</mi><mn>3</mn><mi>j</mi></msubsup></mrow><mi>N</mi></mfrac><mo>,</mo></mrow></math>]]></maths>ζ<sub>O</sub>=1-ζ<sub>D</sub>-ζ<sub>S</sub>-ζ<sub>OFF</sub>,式中ζ<sub>SH</sub>,ζ<sub>SP</sub>,ζ<sub>OFF</sub>,ζ<sub>O</sub>∈[0,1];忽略脉间状态,其余三种放电状态可按下式统计得到<img file="FSA00000346399900034.GIF" wi="410" he="109" /><img file="FSA00000346399900035.GIF" wi="406" he="108" />ζ<sub>O</sub>′=1-ζ′<sub>SH</sub>-ζ′<sub>SP</sub>;(4)第二阶控制系统输入输出隶属度函数的建立第二阶控制系统的输入为前述(3)中所得的短路率ζ<sub>SH</sub>′和火花率ζ<sub>SP</sub>′,输出为电极的运动命令速度v,将其乘以速度系数即可得到微细电极的进给速度;第二阶区间二型模糊控制系统的设计方法同第一阶,在输出速度前首先得到输出速度v的一型模糊区间<img file="FSA00000346399900036.GIF" wi="109" he="58" /><u>v</u>计算方法如下:①首先将各输出论域上的每条隶属度函数重心的左端点<u>v</u><sup><u>l</u></sup>按由大到小的顺序排序,即<u>v</u><sup><u>1</u></sup>≤<u>v</u><sup><u>2</u></sup>≤...≤<u>v</u><sup><u>n</u></sup>,令<img file="FSA00000346399900037.GIF" wi="308" he="70" />且v<sup>l</sup>=<u>v</u><sup><u>l</u></sup>,然后利用<img file="FSA00000346399900038.GIF" wi="205" he="221" />计算v且令<u>v</u>=v;②找到一个值k<sub>1</sub>,其中1≤k<sub>1</sub>≤n-1,使得<img file="FSA00000346399900039.GIF" wi="257" he="53" />③对于l≤k<sub>1</sub>,令<img file="FSA000003463999000310.GIF" wi="166" he="61" />对于l>k<sub>1</sub>,令f<sup>l</sup>=<u>f</u><sup><u>l</u></sup>,然后再利用公式<img file="FSA000003463999000311.GIF" wi="206" he="224" />计算v;④如果<u>v</u>=v,则停止计算,否则令<u>v</u>=v返回②重新进行计算;<img file="FSA000003463999000312.GIF" wi="23" he="46" />的计算方法同<u>v</u>,仅需要更改各变量的上下标;由此得到输出速度的降型集为<img file="FSA000003463999000313.GIF" wi="926" he="257" />最后取平均值得到输出电机命令速度为<img file="FSA000003463999000314.GIF" wi="246" he="96" />
地址 116024 辽宁省大连市甘井子区凌工路2号