发明名称 一种基于启发式算法的高速主轴回转精度评定方法
摘要 本发明提供了一种基于启发式算法的高速主轴回转精度评定方法,首先,根据交叉准则建立最小区域圆法的高速主轴回转精度评定的数学模型,构建目标函数;其次,随机初始化种群的速度、位置、温度等,并计算种群的个体适应度和全局最优值,进而计算当前温度下各个个体适应度值的适配值;再次,将适配值的求和函数作为子目标进行一维搜索,从所有个体适应度值中找到全局最优的替代值,更新种群中每个个体的速度和位置;最后,若未达到预设的最大迭代次数则进行退温操作,继续迭代,直到得出最小区域圆法评价函数的最优解。本发明将高速主轴回转精度的评定问题转换为非线性无约束优化问题,并利用启发式算法对其进行求解,具有较高精度及良好的鲁棒性。
申请公布号 CN105701285A 申请公布日期 2016.06.22
申请号 CN201610015948.3 申请日期 2016.01.12
申请人 西安电子科技大学 发明人 章云;魏乔;张大兴;袁帅
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 西安吉盛专利代理有限责任公司 61108 代理人 何锐
主权项 一种基于启发式算法的高速主轴回转精度评定方法,其特征在于,包括如下步骤:步骤一,根据交叉准则建立基于最小区域圆法的高速主轴回转精度评价模型,由此确定非线性优化的目标函数,如下所示:ψ=max{R<sub>1</sub>,R<sub>2</sub>...R<sub>m</sub>}‑min{R<sub>1</sub>,R<sub>2</sub>...R<sub>m</sub>}其中ψ为目标函数,<img file="FDA0000904689540000011.GIF" wi="634" he="107" />为被测点到圆心的距离,i=1,2...m,m为被测轮廓上的测点个数,(a<sub>i</sub>,b<sub>i</sub>)为最小区域圆法确定的两同心圆的圆心坐标,(xx<sub>i</sub>,yy<sub>i</sub>)为被测轮廓上任一测点;步骤二,将被测轮廓上的所有点看作一个种群,随机初始化种群的位置x<sub>i</sub>、速度v<sub>i</sub>(i=1,2...m);步骤三,根据目标函数计算每个个体的适应度值,将其存储在p<sub>i</sub>中,然后在所有个体适应度值中选取最优值作为全局最优的初始值,存储于p<sub>g</sub>中;步骤四,设置温度的初始值,采用公式t<sub>0</sub>=ψ(p<sub>g</sub>)/ln5来确定,进而根据下式来确定当前温度下的各p<sub>i</sub>的适配值:<maths num="0001"><math><![CDATA[<mrow><mi>T</mi><mi>f</mi><mi>i</mi><mi>t</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>=</mo><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><mi>&psi;</mi><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>&psi;</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>g</mi></msub><mo>)</mo></mrow><mo>)</mo><mo>/</mo><msub><mi>t</mi><mi>i</mi></msub></mrow></msup><mo>/</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msup><mi>e</mi><mrow><mo>-</mo><mrow><mo>(</mo><mi>&psi;</mi><mo>(</mo><msub><mi>p</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mi>&psi;</mi><mrow><mo>(</mo><msub><mi>p</mi><mi>g</mi></msub><mo>)</mo></mrow><mo>)</mo><mo>/</mo><msub><mi>t</mi><mi>i</mi></msub></mrow></msup></mrow>]]></math><img file="FDA0000904689540000012.GIF" wi="850" he="147" /></maths>其中,Tfit(p<sub>i</sub>)即为相应的p<sub>i</sub>所对应的适配值,t<sub>i</sub>为当前的温度;步骤五,将适配值的求和函数作为子目标函数,对其进行一维搜索,找到目标函数的最优值;步骤六,由步骤五的结果得到全局最优值p<sub>g</sub>的替代值<img file="FDA0000904689540000013.GIF" wi="307" he="87" />其中[α<sub>n</sub>,β<sub>n</sub>]为步骤五中最优值所在区间,取(α<sub>n</sub>+β<sub>n</sub>)/2为最优值的近似取值,由此更新粒子的速度和位置;步骤七,迭代次数增加i=i+1,判断是否达到最大迭代次数,即是否满足i<M<sub>iter</sub>,若满足则转步骤八,否则停止搜索,输出结果;步骤八,进行退温操作t<sub>i+1</sub>=λ·t<sub>i</sub>,转步骤三。
地址 710071 陕西省西安市太白南路2号西安电子科技大学