含分数阶项的双圆磨削再生颤振数值分析方法技术

发明专利技术的含分数阶项的双圆磨削再生颤振数值分析方法，首先做出加入分数阶项的受力模型，通过磨削力的计算表达式得到磨削模型，求解不含分数阶模型稳定性并绘制稳定性边界，在matlab程序中运用Runge‑Kutta法迭代编制数值模拟程序，进行数值模拟仿真，采用一组合理参数为例，依据仿真结果，具体分析加入分数阶项对磨削系统的影响；在加入分数阶项后的磨削系统中研究了不同分数阶阶数及对应的分数阶系数对系统响应及稳定性的影响，合理选取分数阶参数能够有效提高系统对再生颤振扰动的稳定性，研究结果能够为分数阶阻尼器在轧辊磨削系统中的应用和优化提供基础，也可以用于轧辊磨削系统参数的优化，还能对实际加工过程提供理论指导。

【技术实现步骤摘要】
含分数阶项的双圆磨削再生颤振数值分析方法
本专利技术属于精密磨削
，涉及含分数阶项的双圆磨削再生颤振数值分析方法。
技术介绍
分数阶微分广泛应用在信号处理、图像处理、流体力学、粘弹性阻尼器等领域，在机械动力学系统中申永军引入了分数阶微分的概念，主要研究了引入分数阶微分项对duffing系统动力学特性的影响，由于分数阶微分能够表征对历史数据的全局依赖性，而且分数阶微分项在不同分数阶系数和分数阶阶次的影响下，能够通过等效线性阻尼和等效线性刚度分别影响系统的振幅和频率特性，相对于分数阶项的特点，与轧辊磨削相对比，轧辊磨削过程颤振公认的其中一种模型为时滞引起的再生颤振，由于时滞的影响，轧辊磨削过程历史数据全局依赖性大，更符合实际，用于颤振预测能够提高准确度，而且轧辊磨削需要对系统的振幅和频率进行调整抑制颤振，基于实际应用的角度，现阶段已经研究出了分数阶阻尼器，分数阶阻尼器在机械系统中的安装应用还不广泛，其可控的分数阶阶数和分数阶系数为控制系统参数扩展了思路，用以得到更好的磨削效果，提高磨削表面质量，对轧辊磨削颤振抑制方面提供理论基础和工程指导方面具有重大意义。解析解的意义不如数值解大，但它为数值解的研究提供了理论基础，实际中数值解应用较为广泛，现有技术中还缺少分数阶项对机械系统尤其是高精度要求的磨削系统的分析，没有对分数阶项与时滞引起的再生颤振耦合进行分析的方法，轧辊磨削是典型的双圆磨削过程，在磨削领域研究进展较滞后，无法从解析或数值的角度对分数阶项和双圆磨削再生颤振影响分析，没有相关理论进行指导，缺乏分数阶阻尼器在轧辊磨削系统中的应用和优化基础，更无法进一步对含分数阶项的轧辊磨削系统参数优化设计。
技术实现思路
为解决上述技术问题，本专利技术的目的是提供含分数阶项的双圆磨削再生颤振数值分析方法，用于轧辊磨削系统参数的优化，并对实际加工过程提供理论指导。含分数阶项的双圆磨削再生颤振数值分析方法，包括如下步骤：步骤1：构建轧辊磨床的砂轮和轧辊之间的加入分数阶项的磨削受力模型；步骤2：根据磨削受力模型推导出磨削力的计算表达式，再建立含分数阶项的双时滞轧辊磨削模型；步骤3：根据含分数阶项的双时滞轧辊磨削模型，建立不含分数阶项的模型，求解该模型稳定性，并绘制稳定性边界；步骤4：按照含分数阶项的双时滞轧辊磨削模型，依托经典四阶Runge-Kutta法，构造嵌套的迭代关系，基于matlab编制数值仿真程序，根据步骤3中解析解求得的稳定性边界进行验证；步骤5：改变输入参数，进行对比研究，并根据研究结果分析含分数阶项的轧辊磨削系统的动力学响应特性，并研究各参数对含分数阶项的轧辊磨削系统的影响，用于优化轧辊磨床的设计参数以及分数阶微分的控制参数。步骤1具体为：简化轧辊磨削过程的物理模型，由于轧辊质量较大以及轧辊的两端通过辊颈和托瓦配合定位，根据轧辊磨削特点，简化不必要的细节，做出假设条件，得到磨削受力模型，该磨削受力模型以轧辊辊颈为基准，砂轮架的位移量作为振动自由端，含时滞磨削力作为激振力；其中，假设条件1为：轧辊的轴线位移量为零；假设条件2为：轧辊不发生弹性变形。步骤2包括：步骤2.1：根据受力模型计算切削厚度的表达式，切削厚度计算方法为：δr＝μ1x3(t-τr)-x3(t)(1)δg＝x2(t-τg)-x2(t)(2)其中，δr为轧辊的磨削厚度，μ1为磨削交叠系数，x3(t)为轧辊表面磨削厚度变化量，τr为轧辊表面造成再生效应的时滞量；δg为砂轮切削厚度，x2(t)为砂轮表面的厚度变化，τg为砂轮的时滞量；步骤2.2：将磨削力与切削厚度关联，进而求取磨削力的计算表达式，根据弹性变形条件，恒定转速磨削条件下，磨削力的计算方法为：Pn＝k0Pt(4)其中，Pt为磨削力的切向分量，uch为磨削比能，Vr为轧辊表面线速度，b为砂轮宽度，Vg为砂轮表面线速度；Pn为磨削力的水平分量，k0为磨削力和正压力的转换系数；步骤2.3：根据公式(3)(4)得到磨削力与再生扰动之间的表达式：引入磨削系数G：G＝轧辊磨损量体积/砂轮磨损量体积＝bδrVr/bδgVg(6)x1(t)＝x2(t)+x3(t)(7)得到砂轮切屑厚度δg＝Vrδr/VgG，由于磨削系数一般取G＝1000，砂轮的线速度高于轧辊的线速度，得到砂轮的磨削磨损量x2(t)小于轧辊磨削磨损量x3(t)的千分之一，x1(t)代表砂轮架的位移量，据此得到磨削力与再生扰动之间的简化表达式：步骤2.4：根据磨削力与再生扰动之间的简化表达式获得含分数阶的双圆磨削过程的再生颤振模型为：其中，m代表砂轮架质量，c代表砂轮架阻尼，k代表砂轮架刚度，KDα[x1(t)]代表分数阶项，K为分数阶系数，α为分数阶阶数。步骤3包括：步骤3.1：建立不含分数阶项的模型，表达式如下：根据此模型定义状态向量然后将不含分数阶项的模型改写为状态向量的矩阵形式：步骤3.2：建立特征方程，通过mathmatica进行解析求解，根据特征方程的解，绘制稳定性边界。步骤4包括：步骤4.1：依托经典四阶Runge-Kutta法，构造嵌套的迭代关系，数值仿真迭代算法为：其中，z为时滞项，代表了旋转一圈时间内位移量为零矩阵，之后时滞项按照x的规律运动，l,n,i,j分别代表迭代次数，Dα为分数阶项的近似表达式，为分数阶二项式系数，其中ti＝ih为时间采样点，h为时间步长，y和分别为辅助构造函数，分别定义迭代初值y(t0)和k1代表起始点斜率，k2代表k1决定的中点斜率，k3代表k2决定的中点斜率，k4代表终点斜率，分配权重迭代计算，得到数值仿真结果；步骤4.2：依据步骤3求得的稳定性边界，选取稳定性边界附近的三个点，分别为：位于稳定性边界临近点、位于稳定性边界内的点和位于稳定性边界外的点，K取值为0，步骤4.1退化为只含时滞项的数值解，其他参数不变，确定转速Vg和砂轮宽度b的取值，代入步骤4.1进行数值求解；步骤4.3：位于稳定性边界临近点表现出了临界稳定特性，位于稳定性边界内的点表现出了收敛稳定特性，位于稳定性边界外的点表现出了发散失稳特性，解析解和数值解互相验证了求解过程的正确性。步骤5中，选取阶数为0时，由于分数阶项体现出很强的弹簧特性，分数阶系数的数量级和砂轮架刚度保持一致，在有效研究区间内，稳定性随着分数阶项系数的增大先减小再增大，证明了稳定性与分数阶项之间不存在线性相关关系，应分段考虑分数阶系数的选取。步骤5中，选取阶数为0.25时，分数阶项兼具弹簧特性和阻尼特性，但更偏向于弹簧特性，分数阶系数的数量级与砂轮架刚度保持一致，在有效研究区间内，稳定性随着分数阶项的增大先增大后减小，与阶数为零的时候正好相反，证明在此段内具有很强的非线性相关性。步骤5中，选取阶数为0.75时，分数阶项兼具弹簧特性和阻尼特性，但更偏向于阻尼特性，分数阶系数的数量级与砂轮架阻尼保持一致，在有效研究区间内，稳定性随着分数阶项的增大一直增加，只是变化的梯度较小，证明在此段内具有线性相关性。步骤5中，选取阶数为1时，分数阶项表现出很强的阻尼特性，分数阶系数的数量级与砂轮架阻尼保持一致，在有效研究区间内，稳定性随着分数阶项的增大一直增加，变化的梯度相对于阶数为0.75时明显增加，证明在此段内具有线性相关性。本专利技术的含分数阶项的双圆磨削再生颤本文档来自技高网...

【技术保护点】
1.含分数阶项的双圆磨削再生颤振数值分析方法，其特征在于，包括如下步骤：步骤1：构建轧辊磨床的砂轮和轧辊之间的加入分数阶项的磨削受力模型；步骤2：根据磨削受力模型推导出磨削力的计算表达式，再建立含分数阶项的双时滞轧辊磨削模型；步骤3：根据含分数阶项的双时滞轧辊磨削模型，建立不含分数阶项的模型，求解该模型稳定性，并绘制稳定性边界；步骤4：按照含分数阶项的双时滞轧辊磨削模型，依托经典四阶Runge‑Kutta法，构造嵌套的迭代关系，基于matlab编制数值仿真程序，根据步骤3中解析解求得的稳定性边界进行验证；步骤5：改变输入参数，进行对比研究，并根据研究结果分析含分数阶项的轧辊磨削系统的动力学响应特性，并研究各参数对含分数阶项的轧辊磨削系统的影响，用于优化轧辊磨床的设计参数以及分数阶微分的控制参数。

【技术特征摘要】
1.含分数阶项的双圆磨削再生颤振数值分析方法，其特征在于，包括如下步骤：步骤1：构建轧辊磨床的砂轮和轧辊之间的加入分数阶项的磨削受力模型；步骤2：根据磨削受力模型推导出磨削力的计算表达式，再建立含分数阶项的双时滞轧辊磨削模型；步骤3：根据含分数阶项的双时滞轧辊磨削模型，建立不含分数阶项的模型，求解该模型稳定性，并绘制稳定性边界；步骤4：按照含分数阶项的双时滞轧辊磨削模型，依托经典四阶Runge-Kutta法，构造嵌套的迭代关系，基于matlab编制数值仿真程序，根据步骤3中解析解求得的稳定性边界进行验证；步骤5：改变输入参数，进行对比研究，并根据研究结果分析含分数阶项的轧辊磨削系统的动力学响应特性，并研究各参数对含分数阶项的轧辊磨削系统的影响，用于优化轧辊磨床的设计参数以及分数阶微分的控制参数。2.如权利要求1所述的含分数阶项的双圆磨削再生颤振数值分析方法，其特征在于，步骤1具体为：简化轧辊磨削过程的物理模型，由于轧辊质量较大以及轧辊的两端通过辊颈和托瓦配合定位，根据轧辊磨削特点，简化不必要的细节，做出假设条件，得到磨削受力模型，该磨削受力模型以轧辊辊颈为基准，砂轮架的位移量作为振动自由端，含时滞磨削力作为激振力；其中，假设条件1为：轧辊的轴线位移量为零；假设条件2为：轧辊不发生弹性变形。3.如权利要求1所述的含分数阶项的双圆磨削再生颤振数值分析方法，其特征在于，步骤2包括：步骤2.1：根据受力模型计算切削厚度的表达式，切削厚度计算方法为：δr＝μ1x3(t-τr)-x3(t)(1)δg＝x2(t-τg)-x2(t)(2)其中，δr为轧辊的磨削厚度，μ1为磨削交叠系数，x3(t)为轧辊表面磨削厚度变化量，τr为轧辊表面造成再生效应的时滞量；δg为砂轮切削厚度，x2(t)为砂轮表面的厚度变化，τg为砂轮的时滞量；步骤2.2：将磨削力与切削厚度关联，进而求取磨削力的计算表达式，根据弹性变形条件，恒定转速磨削条件下，磨削力的计算方法为：Pn＝k0Pt(4)其中，Pt为磨削力的切向分量，uch为磨削比能，Vr为轧辊表面线速度，b为砂轮宽度，Vg为砂轮表面线速度；Pn为磨削力的水平分量，k0为磨削力和正压力的转换系数；步骤2.3：根据公式(3)(4)得到磨削力与再生扰动之间的表达式：引入磨削系数G：G＝轧辊磨损量体积/砂轮磨损量体积＝bδrVr/bδgVg(6)x1(t)＝x2(t)+x3(t)(7)得到砂轮切屑厚度δg＝Vrδr/VgG，由于磨削系数一般取G＝1000，砂轮的线速度高于轧辊的线速度，得到砂轮的磨削磨损量x2(t)小于轧辊磨削磨损量x3(t)的千分之一，x1(t)代表砂轮架的位移量，据此得到磨削力与再生扰动之间的简化表达式：步骤2.4：根据磨削力与再生扰动之间的简化表达式获得含分数阶的双圆磨削过程的再生颤振模型为：其中，m代表砂轮架质量，c代表砂轮架阻尼，k代表砂轮架刚度，KDα[x1(t)]代表分数阶项，K为分数阶系数，α为分数阶阶数。4.如权利要求1所述的含分数阶项...

【专利技术属性】
技术研发人员：付柯，李鹤，邓鹏飞，
申请(专利权)人：东北大学，
类型：发明
国别省市：辽宁,21