CN105750570B - 一种电主轴铣削颤振时延主动控制方法及其系统 - Google Patents
一种电主轴铣削颤振时延主动控制方法及其系统 Download PDFInfo
- Publication number
- CN105750570B CN105750570B CN201610228617.8A CN201610228617A CN105750570B CN 105750570 B CN105750570 B CN 105750570B CN 201610228617 A CN201610228617 A CN 201610228617A CN 105750570 B CN105750570 B CN 105750570B
- Authority
- CN
- China
- Prior art keywords
- mtd
- mrow
- msub
- mtr
- mfrac
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q1/00—Members which are comprised in the general build-up of a form of machine, particularly relatively large fixed members
- B23Q1/70—Stationary or movable members for carrying working-spindles for attachment of tools or work
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Mechanical Engineering (AREA)
- External Artificial Organs (AREA)
- Feedback Control In General (AREA)
Abstract
本公开涉及一种电主轴铣削颤振时延主动控制方法及其系统,应用于包含作动装置的电主轴中,其中所述方法包括步骤:将采集电主轴系统相互正交方向上的振动信号、测定的电主轴系统参数以及切削力系数输入到包含状态时延的连续铣削颤振主动控制状态方程;通过将所述连续铣削颤振主动控制状态方程离散化后,使用离散最优控制算法计算电主轴系统的最优主动控制力,进而获得驱动压电作动器对电主轴系统施加主动控制力的驱动电压。采用本公开方法或系统可有效抑制铣削加工时的颤振,扩大铣削稳定域,改善加工效果并提高了加工效率。
Description
技术领域
本公开属于机械切削加工控制技术领域,更具体地,涉及一种铣削加工颤振时延主动控制方法及其系统。
背景技术
铣削加工广泛应用于航空、航天、模具等行业。在铣削加工中,颤振是制约加工表面质量和生产率的关键因素。颤振是机床-刀具-工件构成的闭环切削系统的动态不稳定现象,根源于切削力的异常变化,造成刀具与工件之间十分强烈而持续的相对振动。颤振是切削加工中的自激振动现象,其中再生型颤振(因铣削厚度变化所产生的动态铣削力激起的自激颤振)是引起振动的主要因素。切削颤振会在零件表面留下振纹,从而降低零件加工质量;由于颤振的发生,在加工过程中不得不采取保守的切削用量从而降低机床加工效率;颤振会导致加工刀具的提前报废,同时产生大量的噪声污染环境。
针对铣削加工过程中的颤振现象,振动主动控制技术被应用于对铣削颤振的控制上。振动主动控制是指在振动控制过程中,根据所检测到的振动信号,应用一定的控制策略,经过实时计算,进而驱动作动器对控制目标施加一定的影响,达到抑制或消除振动的目的。通过对铣削系统进行数学建模,实验测定系统各动力学参数,再设计控制算法是颤振主动控制的核心;但是,铣削动力学系统比较复杂,铣削过程的复杂性不仅在于系统中时延的存在,还在于铣削过程中刀齿旋转切削运动的周期性,这导致铣削动力学方程中的系数呈周期性变化,从而加大了控制算法的设计难度。因此,铣削颤振主动控制方法的设计具有很大的挑战性。
发明内容
针对上述问题,本公开提供了一种电主轴铣削颤振时延主动控制方法及其系统,用于稳定包含时延和非线性动态切削力的铣削动力学系统,从而达到颤振控制,扩大铣削稳定域和改善加工效果的目的。
一种电主轴铣削颤振时延主动控制方法,所述方法包括下述步骤:
S100、建立连续铣削颤振主动控制状态方程;
S200、将所述连续铣削颤振主动控制状态方程转换为离散控制系统状态方程;
S300、利用离散最优控制算法计算所述离散控制系统状态方程,获得最优主动控制力;
S400、利用压电作动器驱动电压-输出力模型,根据计算得到的最优主动控制力计算压电作动器所需的驱动电压,进而驱动压电作动器对电主轴系统施加主动控制力。
进一步地,所述步骤S100之前还包括采集刀柄末端两个正交方向的振动信号、测定电主轴系统参数和切削力系数。
可选地,所述连续铣削颤振主动控制状态方程的建立包括下述步骤:
S101、建立包含状态时延的铣削颤振主动控制力学模型:
式中:
其中:M∈R2×2为电主轴系统质量矩阵,C∈R2×2为电主轴系统阻尼矩阵,K∈R2×2为电主轴系统刚度矩阵,X(t)为刀具正交的x,y方向的位移向量,为刀具正交的x,y方向的速度向量,为刀具正交的x,y方向的加速度向量,b为轴向切削深度,τ为切削系统时延,H(t)∈R2×2为切削力变化矩阵,其周期与切削系统时延相同为τ,U(t)为主动控制力,N为刀齿数,ω为电主轴转速;
S102、利用傅里叶零阶展开对切削力变化矩阵H(t)在一个周期内平均化:
其中:
K11,K12,K21,K22为平均化后的切削力常数;mx,mv为电主轴系统质量矩阵在正交的x,y方向上的分量;cx,cv为电主轴系统阻尼矩阵在正交的x,y方向上的分量;kx,kv为电主轴系统刚度矩阵在正交的x,y方向上的分量;
S103、利用步骤S102中的方程建立连续铣削颤振主动控制状态方程:
其中:
优选地,所述步骤S400中计算获得的驱动电压经功率放大器放大后驱动压电作动器工作。
优选地,所述压电作动器驱动电压-输出力模型为:
其中:
K′ex为外部结构的动刚度;Ka为压电作动器等效刚度;V为施加在压电叠层作动器两端的电压;n为压电片的层数;d33为压电常数。
基于所述方法,实现了一种电主轴铣削颤振时延主动控制系统,所述系统包括实时控制器;所述实时控制器包括下述单元:
参数接收单元:将采集到的刀柄末端两个正交方向的振动信号、以及测定的电主轴系统参数和切削力系数输入到包含状态时延的连续铣削颤振主动控制状态方程;
转换单元:基于连续铣削颤振主动控制状态方程的参数获取离散控制系统方程的参数;
最优主动控制力计算单元:利用离散最优控制算法计算使所述离散控制系统稳定的最优主动控制力;
驱动电压计算单元:利用压电作动器驱动电压-输出力模型,根据计算得到的最优主动控制力计算压电作动器所需的驱动电压。
可选地,所述振动信号使用传感器采集;所述传感器包括位移传感器和加速度传感器,并将采集到的振动信号输出给信号调理模块以进行放大滤波处理。
优选地,所述连续铣削颤振主动控制状态方程的建立包括下述步骤:
S101、建立包含状态时延的铣削颤振主动控制力学模型:
式中:
其中:M∈R2×2为电主轴系统质量矩阵,C∈R2×2为电主轴系统阻尼矩阵,K∈R2×2为电主轴系统刚度矩阵,X(t)为刀具正交的x,y方向的位移向量,为刀具正交的x,y方向的速度向量,为刀具正交的x,y方向的加速度向量,b为轴向切削深度,τ为切削系统时延,H(t)∈R2×2为切削力变化矩阵,其周期与切削系统时延相同为τ,U(t)为主动控制力,N为刀齿数,ω为电主轴转速;
S102、利用傅里叶零阶展开对切削力变化矩阵H(t)在一个周期内平均化:
其中:
K11,K12,K21,K22为平均化后的切削力常数;mx,mv为电主轴系统质量矩阵在正交的x,y方向上的分量;cx,cv为电主轴系统阻尼矩阵在正交的x,y方向上的分量;kx,kv为电主轴系统刚度矩阵在正交的x,y方向上的分量;
S103、利用步骤S102中的方程建立连续铣削颤振主动控制状态方程:
其中:
进一步地,所述驱动电压计算单元中计算获得输出给功率放大器放大电压,进而驱动压电作动器工作。
可选地,所述实时控制器包括DSP、FPGA或者安装了实时系统的计算机。
在本公开技术方案中,可以根据从刀具传递到刀柄末端的两个方向(正交的x方向和y方向)的振动信号、电主轴系统动力学模型和时延最优控制算法实时计算出最优主动控制力,并根据最优主动控制力和作动器驱动电压-输出力模型实时计算出作动器所需的驱动电压,进而驱动作动器对电主轴系统施加主动控制力,从而达到提高加工效果,扩大铣削稳定域和抑制切削颤振的目的。所述系统采用实时控制器对方法中涉及的计算进行实时计算,可以实现实时控制驱动作动器对控制目标施加实时影响,达到实时抑制或消除振动的目的。
附图说明
图1为本公开一个实施例中的电主轴铣削颤振时延主动控制方法工作流程图;
图2为本公开一个实施例中的电主轴铣削颤振主动控制系统构成示意图;
图中:1-刀具,2-位移/加速度传感器,3-信号调理单元,4-实时控制单元,5-功率放大器,6-作动器;
图3为本公开一个实施例中不施加主动控制力时铣削稳定域图,横轴表示转速,单位为转/分钟;纵轴为轴向切削深度,单位为毫米;
图4为本公开一个实施例中不施加主动控制力时主轴系统响应,横轴为时间,单位为转/分钟;纵轴为振动幅值,单位为毫米;
图5为本公开一个实施例中在时延主动控制下的铣削稳定域图;
图6为本公开一个实施例中的时延主动控制下的主轴系统响应。
具体实施方式
为了抑制电主轴在高速铣削过程中的颤振,提高铣削加工精度和效率,在一个实施例中,如图1所示,本公开揭示了一种电主轴铣削颤振时延主动控制方法,所述方法包括下述步骤:
S100、建立连续铣削颤振主动控制状态方程;
S200、将所述连续铣削颤振主动控制状态方程转换为离散控制系统状态方程;
S300、利用离散最优控制算法计算所述离散控制系统状态方程,获得最优主动控制力;
S400、根据压电作动器驱动电压-输出力模型,利用计算得到的最优主动控制力来计算压电作动器所需的驱动电压,进而驱动压电作动器对电主轴系统施加主动控制力。
在本公开中,最优主动控制力是指:在铣削加工过程中,为了达到颤振控制的目的,能够驱动作动装置对控制目标施加的力。
优选地,在所述步骤S100之前还包括:
S000、采集刀柄末端两个正交方向的振动信号、测定电主轴系统参数和切削力系数。所述电主轴系统参数和切削力系数用于输入到包含状态时延的连续铣削颤振主动控制状态方程。其中,所述振动信号包括位移和加速度。根据所述铣削颤振主动控制状态方程的需要,需要实验测定的电主轴系统参数包括质量、阻尼、刚度等。当采集的振动信号为模拟信号时,需将采集的振动信号的模拟信号转变为数字信号;进一步地,步骤S400还包括:将计算出的驱动电压信号的数字信号转变为模拟信号。
优选地,所述振动信号包括刀柄末端x方向的位移信号、速度信号,y方向的位移信号、速度信号。
在一个实施例中,提供了连续铣削颤振主动控制状态方程的建立过程,包括下述步骤:
S101、建立包含状态时延的铣削颤振主动控制力学模型:
式中:
其中:M∈R2×2为电主轴系统质量矩阵,C∈R2×2为电主轴系统阻尼矩阵,K∈R2×2为电主轴系统刚度矩阵,X(t)为刀具正交的x,y方向的位移向量,为刀具正交的x,y方向的速度向量,为刀具正交的x,y方向的加速度向量,b为轴向切削深度,τ为切削系统时延,H(t)∈R2×2为切削力变化矩阵,其周期与切削系统时延相同为τ,U(t)为主动控制力,N为刀齿数,ω为电主轴转速;
S102、利用傅里叶零阶展开对切削力变化矩阵H(t)在一个周期内平均化:
其中:
K11,K12,K21,K22为平均化后的切削力常数;mx,mv为电主轴系统质量矩阵在正交的x,y方向上的分量;cx,cv为电主轴系统阻尼矩阵在正交的x,y方向上的分量;kx,kv为电主轴系统刚度矩阵在正交的x,y方向上的分量;
S103、利用步骤S102中的方程建立连续铣削颤振主动控制状态方程:
其中:
在将包含状态时延的连续铣削颤振主动控制状态方程进行离散化时,为了简化问题,优选地,假设时延为采样周期整数倍,即τ=lT连续控制对象和离散的控制器之间采用零阶保持器连结。
离散化后的状态方程为:
Z(k+1)=FZ(k)+G1Z(k-l)+G2U(k)
其中:
取增广状态为:
得标准离散状态方程:
其中:
对于离散控制系统,通过离散最优控制算法,在获得所述离散控制系统的权矩阵后,优选地,使用matlab中的dlqr函数求解获得所述离散控制系统的最优反馈增益矩阵。
具体地,对上述获得的离散控制系统,定义该离散控制系统的性能目标函数:
其中权矩阵Q、R分别为半正定矩阵和正定矩阵:
在优化权矩阵Q、R后,可使用matlab中dlqr函数直接求解得到最优反馈增益矩阵L1,...Ll+1,则对于t=kT时刻电主轴系统的最优控制力为:
U(k)=-L1Z(k)-L2Z(k-l)-…-Ll+1Z(k-1)。
在一个实施例中,开环仿真结果如图3、图4所示,在不施加主动控制力时铣削稳定域较小,主轴系统响应发散,系统不稳定而发生颤振;闭环仿真结果如图5、图6所示,在时延主动控制下铣削稳定域被扩大,主轴系统响应稳定,系统处于稳态切削,从而达到颤振控制和改善加工效果的目的。
仿真参数:
刀具齿数N | 4 |
质量mx=my | 0.014kg |
阻尼比ζ | 0.049 |
固有频率ωn | 2π×778rad/s |
切向切削力系数kt | 6×108kg/(ms2) |
法向切削力系数kn | 2×108kg/(ms2) |
径向切深a | 2.5mm |
轴向切深b | 1.2mm |
刀具半径r | 2.5mm |
主轴转速ω | 18000r/min |
采样周期T | 50us |
优选地,利用作动器驱动电压-输出力模型计算作动器的驱动电压,并输出以驱使压电作动器工作。所述作动器驱动电压-输出力模型的表达式为:
其中:
K′ex为外部结构的动刚度;Ka为压电作动器等效刚度;V为施加在压电叠层作动器两端的电压;n为压电片的层数;d33为压电常数。
优选地,对步骤S400中计算获得的驱动电压,经功率放大器放大后驱动压电作动器工作。
基于所述方法,在一个实施例中,揭示了一种电主轴铣削颤振时延主动控制系统,以抑制电主轴在高速铣削过程中的颤振,提高铣削加工精度和效率。所述系统包括实时控制器;所述实时控制器包括下述单元:
参数接收单元:将采集到的刀柄末端两个正交方向的振动信号、以及测定的电主轴系统参数和切削力系数输入到包含状态时延的连续铣削颤振主动控制状态方程;
转换单元:基于连续铣削颤振主动控制状态方程的参数获取离散控制系统方程的参数;
最优主动控制力计算单元:利用离散最优控制算法计算使所述离散控制系统稳定的最优主动控制力;
驱动电压计算单元:利用利用压电作动器驱动电压-输出力模型根据计算得到的最优主动控制力计算压电作动器所需的驱动电压。
在一个实施例中,进一步给出了实时控制器中连续铣削颤振主动控制系统的建立过程,其步骤包括:
S101、建立包括状态时延的铣削颤振主动控制力学模型:
式中:
其中:M∈R2×2为电主轴系统质量矩阵,C∈R2×2为电主轴系统阻尼矩阵,K∈R2×2为电主轴系统刚度矩阵,X(t)为刀具正交的x,y方向的位移向量,为刀具正交的x,y方向的速度向量,为刀具正交的x,y方向的加速度向量,τ为切削系统时延,H(t)∈R2×2为切削力变化矩阵,其周期与切削系统时延相同为τ,U(t)为主动控制力,N为刀齿数,ω为电主轴转速;
S102、利用傅里叶零阶展开对切削力变化矩阵H(t)在一个周期内平均化:
其中:
K11,K12,K21,K22为平均化后的切削力常数;mx,mv为电主轴系统质量矩阵在正交的x,y方向上的分量;cx,cv为电主轴系统阻尼矩阵在正交的x,y方向上的分量;kx,kv为电主轴系统刚度矩阵在正交的x,y方向上的分量;
S103、利用步骤S102中的方程建立连续铣削颤振主动控制状态方程:
其中:
在一个实施例中,提供了如图2所示的电主轴铣削颤振主动控制系统。其中,主动作动装置位于电主轴前端轴承或刀柄末端。当刀具1在高速铣削发生颤振时,振动传递到前端轴承上,主动控制装置将最优主动控制力施加在刀柄末端或前端轴承上,使得加工从颤振区域进入稳定区域。根据检测到的振动信号,利用实时控制器进行实时计算,计算后的控制信号经过功率放大器后驱动压电作动器工作。如图2所示,电主轴铣削颤振主动控制系统包括:
传感器2,用于采集传递到刀柄末端两个方向的振动信号(正交的x方向和y方向),传感器的安装方向与作动器的施力方向相同,x、y方向各布置一个位移传感器和加速度传感器;
信号调理模块3,为传感器提供正常工作所需的激励电源,同时对传感器输出的信号进行放大、滤波等处理;
实时控制器4:用于根据采集到的振动信号,利用内嵌的时延最优控制算法实时计算出主动控制力和驱动电压;所述实时控制器优选是可编程的嵌入式硬件,比如DSP、FPGA或者安装了实时系统的计算机。采用实时控制器对方法中涉及的计算进行实时计算,可以实现实时控制驱动作动器对控制目标施加实时影响,达到实时抑制或消除振动的目的。
功率放大器5:放大驱动电压,驱动作动器工作;
作动器6:对电主轴颤振控制系统施加主动控制力,例如压电作动器。
以上对本公开进行了详细介绍,本文中应用了具体个例对本公开的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本公开的方法及其核心思想;同时,对于本领域技术人员,依据本公开的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本公开的限制。
Claims (8)
1.一种电主轴铣削颤振时延主动控制方法,其特征在于,所述方法包括下述步骤:
S100、建立连续铣削颤振主动控制状态方程;
S200、将所述连续铣削颤振主动控制状态方程转换为离散控制系统状态方程;
S300、利用离散最优控制算法计算所述离散控制系统状态方程,获得最优主动控制力;
S400、根据压电作动器驱动电压-输出力模型,利用计算得到的最优主动控制力计算压电作动器所需的驱动电压,进而驱动压电作动器对电主轴系统施加主动控制力;
所述连续铣削颤振主动控制状态方程的建立包括下述步骤:
S101、建立包含状态时延的铣削颤振主动控制力学模型:
<mrow>
<mi>M</mi>
<mover>
<mi>X</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>C</mi>
<mover>
<mi>X</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>K</mi>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>b</mi>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&lsqb;</mo>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>+</mo>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
式中:
其中:M∈R2×2为电主轴系统质量矩阵,C∈R2×2为电主轴系统阻尼矩阵,K∈R2×2为电主轴系统刚度矩阵,X(t)为刀具正交的x,y方向的位移向量,为刀具正交的x,y方向的速度向量,为刀具正交的x,y方向的加速度向量,b为轴向切削深度,τ为切削系统时延,H(t)∈R2×2为切削力变化矩阵,其周期与切削系统时延相同为τ,U(t)为主动控制力,N为刀齿数,ω为电主轴转速;
S102、利用傅里叶零阶展开对切削力变化矩阵H(t)在一个周期内平均化:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>c</mi>
<mi>x</mi>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msub>
<mi>c</mi>
<mi>y</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>k</mi>
<mi>x</mi>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msub>
<mi>k</mi>
<mi>y</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>x</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>y</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>K</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>K</mi>
<mn>12</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>K</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>K</mi>
<mn>22</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>x</mi>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
<mo>-</mo>
<mi>x</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>y</mi>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
<mo>-</mo>
<mi>y</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中:
K11,K12,K21,K22为平均化后的切削力常数;mx,my为电主轴系统质量矩阵在正交的x,y方向上的分量;cx,cy为电主轴系统阻尼矩阵在正交的x,y方向上的分量;kx,ky为电主轴系统刚度矩阵在正交的x,y方向上的分量;
S103、利用步骤S102中的方程建立连续铣削颤振主动控制状态方程:
<mrow>
<mover>
<mi>Z</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>B</mi>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中:
<mrow>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>11</mn>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>x</mi>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>c</mi>
<mi>x</mi>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>12</mn>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>21</mn>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>22</mn>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>y</mi>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>c</mi>
<mi>y</mi>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<msub>
<mi>K</mi>
<mn>11</mn>
</msub>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<msub>
<mi>K</mi>
<mn>12</mn>
</msub>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<msub>
<mi>K</mi>
<mn>21</mn>
</msub>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<msub>
<mi>K</mi>
<mn>22</mn>
</msub>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
1
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
2.根据权利要求1所述的方法,其特征在于,在所述步骤S100之前还包括:
S000、采集刀柄末端两个正交方向的振动信号、测定电主轴系统参数和切削力系数。
3.根据权利要求1所述的方法,其特征在于,步骤S400中计算获得的驱动电压经功率放大器放大后驱动压电作动器工作。
4.根据权利要求1所述的方法,其特征在于,所述压电作动器驱动电压-输出力模型为:
<mrow>
<mi>V</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>K</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<msubsup>
<mi>K</mi>
<mrow>
<mi>e</mi>
<mi>x</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
</mrow>
<mrow>
<msub>
<mi>nd</mi>
<mn>33</mn>
</msub>
<msub>
<mi>K</mi>
<mi>a</mi>
</msub>
<msubsup>
<mi>K</mi>
<mrow>
<mi>e</mi>
<mi>x</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
</mrow>
</mfrac>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
</mrow>
其中:
K′ex为外部结构的动刚度;Ka为压电作动器等效刚度;V为施加在压电叠层作动器两端的电压;n为压电片的层数;d33为压电常数;U为第k个周期的电主轴系统的最优控制力。
5.一种电主轴铣削颤振时延主动控制系统,其特征在于,所述系统包括实时控制器;所述实时控制器包括下述单元:
参数接收单元:将采集到的刀柄末端两个正交方向的振动信号、以及测定的电主轴系统参数和切削力系数输入到包含状态时延的连续铣削颤振主动控制状态方程;
转换单元:基于所述连续铣削颤振主动控制方程转换为离散控制系统状态方程;
最优主动控制力计算单元:利用离散最优控制算法计算使所述离散控制系统稳定的最优主动控制力;
驱动电压计算单元:利用压电作动器驱动电压-输出力模型,根据计算得到的最优主动控制力计算压电作动器所需的驱动电压;
所述连续铣削颤振主动控制状态方程的建立包括下述步骤:
S101、建立包含状态时延的铣削颤振主动控制力学模型:
<mrow>
<mi>M</mi>
<mover>
<mi>X</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>C</mi>
<mover>
<mi>X</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>K</mi>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>b</mi>
<mi>H</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&lsqb;</mo>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>+</mo>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
式中:
其中:M∈R2×2为电主轴系统质量矩阵,C∈R2×2为电主轴系统阻尼矩阵,K∈R2×2为电主轴系统刚度矩阵,X(t)为刀具正交的x,y方向的位移向量,为刀具正交的x,y方向的速度向量,为刀具正交的x,y方向的加速度向量,b为轴向切削深度,τ为切削系统时延,H(t)∈R2×2为切削力变化矩阵,其周期与切削系统时延相同为τ,U(t)为主动控制力,N为刀齿数,ω为电主轴转速;
S102、利用傅里叶零阶展开对切削力变化矩阵H(t)在一个周期内平均化:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mover>
<mi>x</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mover>
<mi>y</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>c</mi>
<mi>x</mi>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msub>
<mi>c</mi>
<mi>y</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>k</mi>
<mi>x</mi>
</msub>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<msub>
<mi>k</mi>
<mi>y</mi>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>x</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>y</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>K</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>K</mi>
<mn>12</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>K</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>K</mi>
<mn>22</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>x</mi>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
<mo>-</mo>
<mi>x</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>y</mi>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
<mo>-</mo>
<mi>y</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>+</mo>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中:
K11,K12,K21,K22为平均化后的切削力常数;mx,my为电主轴系统质量矩阵在正交的x,y方向上的分量;cx,cy为电主轴系统阻尼矩阵在正交的x,y方向上的分量;kx,ky为电主轴系统刚度矩阵在正交的x,y方向上的分量;
S103、利用步骤S102中的方程建立连续铣削颤振主动控制状态方程:
<mrow>
<mover>
<mi>Z</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>B</mi>
<mi>U</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
其中:
<mrow>
<mi>Z</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow>
<mover>
<mi>y</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
<mo>;</mo>
</mrow>
<mrow>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>11</mn>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>x</mi>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>c</mi>
<mi>x</mi>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>12</mn>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>21</mn>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>22</mn>
</msub>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>y</mi>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>c</mi>
<mi>y</mi>
</msub>
</mrow>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
<msub>
<mi>A</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<msub>
<mi>K</mi>
<mn>11</mn>
</msub>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<msub>
<mi>K</mi>
<mn>12</mn>
</msub>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<msub>
<mi>K</mi>
<mn>21</mn>
</msub>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<msub>
<mi>K</mi>
<mn>22</mn>
</msub>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mi>x</mi>
</msub>
</mfrac>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<msub>
<mi>m</mi>
<mi>y</mi>
</msub>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>.</mo>
</mrow>
6.根据权利要求5所述的系统,其特征在于:
所述振动信号使用传感器采集;
所述传感器包括位移传感器和加速度传感器,并将采集到的振动信号输出给信号调理模块以进行放大滤波处理。
7.根据权利要求5所述的系统,其特征在于,所述驱动电压计算单元将计算获得的驱动电压输出给功率放大器放大电压,进而驱动压电作动器工作。
8.根据权利要求5所述的系统,其特征在于,所述实时控制器包括DSP、FPGA或者安装了实时系统的计算机。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610228617.8A CN105750570B (zh) | 2016-04-13 | 2016-04-13 | 一种电主轴铣削颤振时延主动控制方法及其系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610228617.8A CN105750570B (zh) | 2016-04-13 | 2016-04-13 | 一种电主轴铣削颤振时延主动控制方法及其系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105750570A CN105750570A (zh) | 2016-07-13 |
CN105750570B true CN105750570B (zh) | 2017-10-20 |
Family
ID=56334924
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610228617.8A Active CN105750570B (zh) | 2016-04-13 | 2016-04-13 | 一种电主轴铣削颤振时延主动控制方法及其系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105750570B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106475581B (zh) * | 2016-09-28 | 2018-06-26 | 西安交通大学 | 一种主轴轴承外圈径向液压加载抑制振动的装置和方法 |
CN106903547B (zh) * | 2017-04-20 | 2023-03-28 | 中南大学 | 适用于高速铣削加工的电主轴主动振动控制装置及系统 |
CN107423502B (zh) * | 2017-07-20 | 2019-07-12 | 西安交通大学 | 基于铣刀螺旋角和弯曲效应的铣削优化方法及优化系统 |
CN108490789A (zh) * | 2018-05-08 | 2018-09-04 | 西安交通大学 | 一种电主轴铣削颤振鲁棒主动控制方法及其主动控制系统 |
CN109396948A (zh) * | 2018-11-16 | 2019-03-01 | 华中科技大学 | 一种提高五轴铣削系统加工稳定性的主动控制方法及系统 |
JP6901051B1 (ja) * | 2019-08-09 | 2021-07-14 | 住友電気工業株式会社 | 転削工具、モジュールおよび切削システム |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0315469A2 (en) * | 1987-11-06 | 1989-05-10 | Seiko Seiki Kabushiki Kaisha | Machining apparatus and method |
CN103252680A (zh) * | 2013-04-24 | 2013-08-21 | 华中科技大学 | 一种铣削加工颤振主动控制系统及其方法 |
CN104076734A (zh) * | 2014-06-26 | 2014-10-01 | 天津职业技术师范大学 | 一种铣削加工颤振在线寻优方法 |
CN104298799A (zh) * | 2014-05-07 | 2015-01-21 | 上海交通大学 | 基于加工振动的立铣刀关键几何参数设计方法 |
CN104647132A (zh) * | 2014-12-22 | 2015-05-27 | 华中科技大学 | 一种基于磁悬浮轴承电主轴的铣削颤振主动控制方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8805910B2 (en) * | 2011-03-04 | 2014-08-12 | Kylowave Inc. | System and method for generating discrete-time model (DTM) of continuous-time model (CTM) for a dynamical system |
-
2016
- 2016-04-13 CN CN201610228617.8A patent/CN105750570B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0315469A2 (en) * | 1987-11-06 | 1989-05-10 | Seiko Seiki Kabushiki Kaisha | Machining apparatus and method |
CN103252680A (zh) * | 2013-04-24 | 2013-08-21 | 华中科技大学 | 一种铣削加工颤振主动控制系统及其方法 |
CN104298799A (zh) * | 2014-05-07 | 2015-01-21 | 上海交通大学 | 基于加工振动的立铣刀关键几何参数设计方法 |
CN104076734A (zh) * | 2014-06-26 | 2014-10-01 | 天津职业技术师范大学 | 一种铣削加工颤振在线寻优方法 |
CN104647132A (zh) * | 2014-12-22 | 2015-05-27 | 华中科技大学 | 一种基于磁悬浮轴承电主轴的铣削颤振主动控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105750570A (zh) | 2016-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105750570B (zh) | 一种电主轴铣削颤振时延主动控制方法及其系统 | |
CN104647132B (zh) | 一种基于磁悬浮轴承电主轴的铣削颤振主动控制方法 | |
CN104517033B (zh) | 一种面向能量效率的数控加工工艺参数多目标优化方法 | |
CN108490789A (zh) | 一种电主轴铣削颤振鲁棒主动控制方法及其主动控制系统 | |
CN107457609A (zh) | 基于刚度变化的铣削颤振抑制方法及铣削颤振优化系统 | |
CN103268430B (zh) | 基于机床刀具动刚度测量的铣削工艺参数优化方法 | |
CN103984285A (zh) | 一种多约束五轴加工进给率定制方法 | |
CN113946922B (zh) | 一种五轴联动铣削过程的动力学集成建模和加工精度预测方法 | |
Al-Zaharnah | Suppressing vibrations of machining processes in both feed and radial directions using an optimal control strategy: the case of interrupted cutting | |
Chen et al. | Cutting mechanism investigation in vibration-assisted machining | |
CN114112291B (zh) | 一种模型适应型振动主动抑制方法 | |
CN113156885B (zh) | 基于自适应模糊推理的电主轴铣削颤振调控方法及装置 | |
Zhou et al. | Energy consumption prediction model of plane grinder processing system based on BP neural network | |
Omar et al. | Control of interval parameter overhead cranes via Monte Carlo simulation | |
CN116107267B (zh) | 一种数控机床控制参数优化方法及装置 | |
Ben Jdidia et al. | Optimizing cutting conditions in single pass face milling for minimum cutting energy, time, cost, and surface roughness | |
Kleinwort et al. | Simulation-based dimensioning of the required actuator force for active vibration control | |
Long-Xiang et al. | Active control experiment study of a flexible beam with multiple time delays | |
Xu et al. | Research of gyroscopic effects on the stability of high speed milling | |
CN111015346B (zh) | 一种刀具切削颤振控制方法、装置以及设备 | |
Hino et al. | Prediction of chatter in high-speed milling by means of fuzzy neural networks | |
CN103869751A (zh) | X-c直驱磨削廓形误差非线性耦合控制系统及方法 | |
Tolochko et al. | Suppression of structural vibrations of an overhead crane using Takagi-Sugeno-Kang fuzzy controller | |
Xu et al. | Singularity analysis of Jeffcott rotor-magnetic bearing with time delays | |
Huijuan et al. | The study of adaptive control based on fuzzy PID for laser cutting |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |