一种基于动态测试数据的进给系统闭环时域仿真方法技术方案

本发明专利技术提供一种基于动态测试数据的进给系统闭环时域仿真方法，包括如下步骤：1.基于伺服电机自激励进行动态测试与频响估计，对激励与响应信号进行频响估计得到1入3出频响阵；2.基于进给系统控制结构，得到指令位移输入到工作台面位移响应输出的闭环频率特性；3.利用傅里叶逆变换求取进给系统的单位脉冲响应；将单位脉冲响应与指定位移指令进行时域卷积，求得闭环时域响应输出作为进给系统闭环时域仿真输出。本发明专利技术克服了相关方法的建模复杂性与建模精度或辨识精度的问题，具有简单直接、精度高的特点。提供了一种进给系统高阶动态特性的实验建模与时域仿真方法，可用于单轴运动精度的预测与评价，也可用于控制器的设计与参数优化设置。

【技术实现步骤摘要】
一种基于动态测试数据的进给系统闭环时域仿真方法
本专利技术属于进给系统动态特性建模与时域仿真技术，涉及动态测试技术、频响估计方法与时域仿真方法，具体为一种基于动态测试数据的进给系统闭环时域仿真方法。
技术介绍
数控机床是重要的机械制造装备，进给轴的性能是机床实现高速高精加工的前提与基础。在机床设计与调试阶段，往往需要通过对进给轴的指令输入与实际响应进行时域考察，从而分析其稳态误差、动态误差，进行性能评价，用于进一步的优化设计与控制器参数优化设置。进行时域分析最直接的方法是利用系统的内部信号采集功能，即利用伺服或数控系统提供的扭矩指令、位移指令和位移反馈等信号进行分析。这种方法的一个局限是无法采集闭环外的工作台处或主轴端部的运动信号，而这些点的误差响应正是直接影响加工质量的；另一个局限是该法无法应用于设计阶段，必须在机械、控制完全安装好后才能施行，此时进行的分析得到的结果很难再用于机械的结构优化与控制器的选取。另一个最常用的方法是采用多刚体动力学或有限元方法对机械部分进行建模，再考虑控制作用构成机电集成模型，往往获得传递函数形式或状态空间方程形式的模型，再进行时域仿真分析。这种方法的局限是存在模型复杂性与准确性的矛盾，采用低阶模型计算效率高、便于分析，但误差大，特别是在动态误差的预测上误差更大；采用高阶或有限元模型则带来计算效率低、模型参数辨识困难等问题。第三种方法是将动态特性测试得到的脉冲响应或频率响应辨识为传递函数或状态空间模型，再进行时域分析。这种方法同样面临待辨识模型阶次的选择、辨识效率与精度的权衡等问题。
技术实现思路
针对现有技术中存在的问题，本专利技术提供一种基于动态测试数据的进给系统闭环时域仿真方法，直接应用动态测试数据进行时域仿真，具有计算过程简洁和精度高的特点，同时避免了直接建模或模型辨识所带来的效率与精度问题。本专利技术是通过以下技术方案来实现：一种基于动态测试数据的进给系统闭环时域仿真方法，包括如下步骤，步骤1，基于伺服电机自激励进行动态测试与频响估计，对激励与响应信号进行频响估计得到1入3出频响阵；步骤2，根据被控对象模型和步骤1中得到的1入3出频响阵，基于进给系统的控制结构，得到指令位移输入到工作台面位移响应输出的闭环频率特性；步骤3，根据步骤2中的闭环频率特性，利用傅里叶逆变换求取进给系统的单位脉冲响应；将单位脉冲响应与指定位移指令进行时域卷积，求得闭环时域响应输出作为进给系统闭环时域仿真输出。优选的，步骤1中，使用进给系统的伺服电机产生一定幅值与频宽的扫频激励扭矩，分别利用编码器、长光栅与加速度传感器采集电机转子、滑台底部与工作台面处的动态响应，对激励与响应信号进行频响估计求取电机-编码器频响、电机-长光栅频响和电机-工作台频响，并构成1入3出频响阵。进一步，采用的电机扫频激励扭矩如下，T1(t)＝Asin(ωt(t)t)ωt(t)＝a+bt,t∈[0,t0]式中，A为激励幅值，ωt为时变频率，a，b为时变频率的线性变化参数，t0为激励时长。进一步，测得编码器处速度响应ω1(t)，长光栅处位移响应x2(t)，工作台面处加速度响应a3(t)；使用功率谱估计计算三个位移频响分别为，式中，h11，h21，h31依次为电机扭矩输入下的编码器处、长光栅处和工作台面处的位移频响，j为虚数单位，ω为角频率，GT1T1,GT1ω1,GT1x2,GT1a3分别为输入扭矩的自能量谱与输入扭矩分别于三个输出的互能量谱；上述三个位移频响构成进给系统机械部分的1入3出位移频响阵模型如下，进一步，激励幅值的选择原则是能够克服静摩擦，激励起工作台的宏观移动；激励频率范围能够覆盖最大加速度参数下运动指令加减速过程的频宽，并包含最大主轴转速下的切削力干扰频率。进一步，步骤2中，考虑伺服电机与控制时延下的被控对象模型为：式中，kI为电流放大系数，kT为扭矩常数，e-jωτ为控制时延，isc为丝杠螺母副的传动比。进一步，基于包括位置环控制器与速度环控制器的控制结构，得到指令位移输入到工作台面位移响应输出的闭环频率特性如下，式中，r(jω)，x3(jω)分别为指令位移与工作台处位移响应的傅里叶变换，Gcp，Gcv分别为位置环控制器与速度环控制器的传递函数。进一步，步骤3中，根据步骤2中的闭环频率特性，利用傅里叶逆变换求取系统的单位脉冲响应如下，进一步，将单位脉冲响应gx3_r(t)与指定位移指令r(t)进行时域卷积，求得闭环时域响应输出为，相应地，有误差响应为，e(t)＝x3(t)-r(t)。与现有技术相比，本专利技术具有以下有益的技术效果：本专利技术直接应用动态测试数据进行闭环时域仿真，无需进行基于机理的动态特性建模，也不必进行模型形式转换或辨识，克服了相关方法的建模复杂性与建模精度或辨识精度的问题，具有简单直接、精度高的特点。提供了一种进给系统高阶动态特性的实验建模与时域仿真方法，可用于单轴运动精度的预测与评价，也可用于控制器的设计与参数优化设置。附图说明图1为基于动态测试数据的进给系统时域仿真方法流程图。图2为动态测试与频响估计原理框图。图3为实施例测试得到的扭矩-编码器频响h11。图4为实施例测试得到的扭矩-工作台频响h31。图5为一般的进给系统级联控制结构框图。图6为计算得到不同控制器下的闭环频率特性Gx3-r。图7为计算得到不同控制器下的闭环时域仿真结果。图中：1为电机与编码器，2为长光栅与读数头，3为置于工作台上的加速度计，4为考虑了电机与控制时延的被控对象H’，5为速度环控制器Gcv，6为位置环控制器Gcp。具体实施方式下面结合具体的实施例对本专利技术做进一步的详细说明，所述是对本专利技术的解释而不是限定。本专利技术一种基于动态测试数据的进给系统闭环时域仿真方法，如图1所示，其包括三个步骤：基于伺服电机自激励的动态测试与频响估计；基于控制结构的闭环频率特性计算；基于傅里叶逆变换与时域卷积的闭环时域仿真计算输出。具体的如下所述。步骤1，基于伺服电机自激励的动态测试与频响估计，如图2所示，使用进给系统的伺服电机产生一定幅值与频宽的扫频激励扭矩，分别利用编码器、长光栅与加速度传感器采集电机转子、滑台底部与工作台面处的动态响应，对激励与响应信号进行频响估计求取电机-编码器频响、电机-长光栅频响和电机-工作台频响，并构成1入3出频响阵。采用的电机扫频激励扭矩：T1(t)＝Asin(ωt(t)t)ωt(t)＝a+bt,t∈[0,t0]式中，A为激励幅值，ωt为时变频率，a，b为时变频率的线性变化参数，t0为激励时长。激励幅值的选择原则是能够克服静摩擦，激励起工作台的宏观移动，本例中A取为30％的电机额定扭矩；激励频率范围应覆盖最大加速度参数下运动指令加减速过程的频宽，并包含最大主轴转速下的切削力干扰频率；最大主轴转速下的切削力干扰频率为主轴转速与刀齿数之乘积；本例中取为0～1KHz，此时相关参数为a＝0rad/s,b＝314000rad/s2,t0＝0.2s。在上述激励下，测得编码器处速度响应ω1(t)，长光栅处位移响应x2(t)，工作台面处加速度响应a3(t)。使用功率谱估计计算三个位移频响：式中，h11，h21，h31依次为电机扭矩输入下的编码器处、长光栅处、工作台面处的位移频响，j为虚数单位，ω为角频率，GT1T1,GT1本文档来自技高网...

【技术保护点】
一种基于动态测试数据的进给系统闭环时域仿真方法，其特征在于，包括如下步骤，步骤1，基于伺服电机自激励进行动态测试与频响估计，对激励与响应信号进行频响估计得到1入3出频响阵；步骤2，根据被控对象模型和步骤1中得到的1入3出频响阵，基于进给系统的控制结构，得到指令位移输入到工作台面位移响应输出的闭环频率特性；步骤3，根据步骤2中的闭环频率特性，利用傅里叶逆变换求取进给系统的单位脉冲响应；将单位脉冲响应与指定位移指令进行时域卷积，求得闭环时域响应输出作为进给系统闭环时域仿真输出。

【技术特征摘要】
1.一种基于动态测试数据的进给系统闭环时域仿真方法，其特征在于，包括如下步骤，步骤1，基于伺服电机自激励进行动态测试与频响估计，对激励与响应信号进行频响估计得到1入3出频响阵；步骤2，根据被控对象模型和步骤1中得到的1入3出频响阵，基于进给系统的控制结构，得到指令位移输入到工作台面位移响应输出的闭环频率特性；步骤3，根据步骤2中的闭环频率特性，利用傅里叶逆变换求取进给系统的单位脉冲响应；将单位脉冲响应与指定位移指令进行时域卷积，求得闭环时域响应输出作为进给系统闭环时域仿真输出。2.根据权利要求1所述的一种基于动态测试数据的进给系统闭环时域仿真方法，其特征在于，步骤1中，使用进给系统的伺服电机产生一定幅值与频宽的扫频激励扭矩，分别利用编码器、长光栅与加速度传感器采集电机转子、滑台底部与工作台面处的动态响应，对激励与响应信号进行频响估计求取电机-编码器频响、电机-长光栅频响和电机-工作台频响，并构成1入3出频响阵。3.根据权利要求2所述的一种基于动态测试数据的进给系统闭环时域仿真方法，其特征在于，采用的电机扫频激励扭矩如下，T1(t)＝Asin(ωt(t)t)ωt(t)＝a+bt,t∈[0,t0]式中，A为激励幅值，ωt为时变频率，a，b为时变频率的线性变化参数，t0为激励时长。4.根据权利要求3所述的一种基于动态测试数据的进给系统闭环时域仿真方法，其特征在于，测得编码器处速度响应ω1(t)，长光栅处位移响应x2(t)，工作台面处加速度响应a3(t)；使用功率谱估计计算三个位移频响分别为，式中，h11，h21，h31依次为电机扭矩输入下的编码器处、长光栅处和工作台面处的位移频响，j为虚数单位，ω为角频率，GT1T1,GT1ω1,GT1x2,GT1a3分别为输入扭矩的自能量谱与输入扭矩分别于三个输出的互能量谱；上述三个位移频响构成进给系统机械部分的1入3出位移频响阵模型如下，

【专利技术属性】
技术研发人员：刘辉，赵万华，杨晓君，吕盾，张俊，张会杰，
申请(专利权)人：西安交通大学，
类型：发明
国别省市：陕西,61