CN114609975A - 一种基于复合自抗扰控制的fts控制方法 - Google Patents
一种基于复合自抗扰控制的fts控制方法 Download PDFInfo
- Publication number
- CN114609975A CN114609975A CN202210321697.7A CN202210321697A CN114609975A CN 114609975 A CN114609975 A CN 114609975A CN 202210321697 A CN202210321697 A CN 202210321697A CN 114609975 A CN114609975 A CN 114609975A
- Authority
- CN
- China
- Prior art keywords
- fts
- time
- disturbance
- neural network
- estimation
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 42
- 239000002131 composite material Substances 0.000 title claims abstract description 33
- 238000005520 cutting process Methods 0.000 claims abstract description 58
- 238000013528 artificial neural network Methods 0.000 claims abstract description 46
- 239000000919 ceramic Substances 0.000 claims description 22
- 230000008859 change Effects 0.000 claims description 13
- 239000002184 metal Substances 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 9
- 238000013016 damping Methods 0.000 claims description 7
- 238000005070 sampling Methods 0.000 claims description 4
- 230000005284 excitation Effects 0.000 claims description 3
- 238000012886 linear function Methods 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 210000002569 neuron Anatomy 0.000 claims description 3
- 238000012549 training Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000003044 adaptive effect Effects 0.000 abstract description 3
- 238000003754 machining Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 9
- 238000007514 turning Methods 0.000 description 6
- 238000013178 mathematical model Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012545 processing Methods 0.000 description 2
- 230000001172 regenerating effect Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000011217 control strategy Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 230000010365 information processing Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000005459 micromachining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/414—Structure of the control system, e.g. common controller or multiprocessor systems, interface to servo, programmable interface controller
- G05B19/4142—Structure of the control system, e.g. common controller or multiprocessor systems, interface to servo, programmable interface controller characterised by the use of a microprocessor
-
- 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
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/34—Director, elements to supervisory
- G05B2219/34013—Servocontroller
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Computer Hardware Design (AREA)
- Operations Research (AREA)
- Manufacturing & Machinery (AREA)
- Algebra (AREA)
- Automation & Control Theory (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于复合自抗扰控制的FTS控制方法,属于精密加工领域。所述方法在构建FTS***模型时,引入了完全未知的迟滞非线性函数、状态时滞和时变时滞的双时滞***以及内部未建模动态非线性,针对所构建的FTS***模型设计的控制器将迟滞时滞及***内部的动态非线性性视为对FTS***的总扰动,由线性观测器实时估计,并利用自适应神经网络逼近未知的非线性误差函数,实现迟滞时滞补偿功能。从而较好描述精密切削过程中的颤振现象,解决了现有FTS模型无法解释和预测切削颤振导致的幅值失真现象、稳定振幅与切削参数的关系、跳变现象对切削过程稳定性的影响及外界干扰对切削精度影响等问题。
Description
技术领域
本发明涉及一种基于复合自抗扰控制的FTS控制方法,属于精密加工领域。
背景技术
在追求产品小型化、集成化的今天,微机械零件被广泛应用于信息处理、空间技术和光纤通讯等诸多领域。微机械零件的加工过程通常是在车床上安装有快速响应能力的刀具微进给伺服***完成的,刀具微进给伺服***又称快速刀具伺服(Fast Tool Servo,FTS)***,车削过程中,FTS驱动刀具完成高频高精度跟踪运动,是加工微纳米级零件的关键部件。
现有的FTS包括压电陶瓷型FTS、磁致伸缩型FTS、洛伦兹力FTS以及麦克斯韦力FTS;其中压电陶瓷型FTS具有响应速度快,加速度高、频响范围可达数千赫兹等优点,因此在非轴对称车削中得到了广泛应用并取得了良好效果,但是由于精密车削过程往往因为自激振动而产生切削颤振,不仅会影响加工质量,限制切削效率,还会产生噪声污染,增大***机械磨损,大大降低了机床和刀具的使用寿命。由于切削颤振机理的复杂性以及切削过程众多干扰因素影响,很难对其进行有效的预测和消除。因此,研究和分析FTS***切削过程中切削力与切削厚度的关系,建立相关数学模型,进而设计控制器来抑制切削颤振具有重要的理论意义和工程实用价值。
近年来,在快速刀具伺服***建模与控制领域中,非线性时滞FTS模型渐渐取代线性时滞FTS模型成为主流,进而对非线性时滞FTS模型反演求其逆模型进行控制,但这往往需要获得被控对象的精确数学模型。而现有的非线性时滞FTS模型不能解释和预测切削颤振导致的幅值失真现象、稳定振幅与切削参数的关系、跳变现象对切削过程稳定性的影响及外界干扰对切削精度影响等,因此无法进一步提供更为精确的数学模型,也就导致切削精度无法进一步提高。
自抗扰控制是一种新的无模型抗扰控制,然而传统自抗扰控制参数整定较为复杂,并且一组参数只适用于单一***。线性自抗扰控制解决了传统自抗扰控制参数整定复杂的难题,然而,在实际车削过程中,控制器带宽会受到如生产成本、噪声灵敏度等诸多因素的限制。因此,当观测器带宽无法达到实际需要数值时,切削精度无法得到保证。
发明内容
为了进一步提高切削精度,本发明提供了一种基于复合自抗扰控制的FTS控制方法,所述方法应用于压电陶瓷型FTS中,所述压电陶瓷型FTS包括主轴、工作间、切刀以及压电作动器,其中切刀通过弹簧与压电作动器相连;所述方法包括:
步骤1:基于主轴旋转下的状态时滞β1、压电作动器迟滞非线性以及内部未建模动态非线性的结合,以及切刀横向给进产生的时变时滞β2建立FTS***模型;所述内部未建模动态非线性指压电陶瓷型FTS中除时变时滞、状态时滞以及迟滞非线性之外的其他未能明确的动态非线性总和;
步骤2:针对步骤1建立的FTS***模型,将***中主轴旋转下的状态时滞β1、压电作动器迟滞非线性以及内部未建模动态非线性的结合,以及切刀横向给进产生的时变时滞β2视为总扰动,采用线性自抗干扰控制器对总扰动进行初步估计;
步骤3:针对步骤2采用线性自抗干扰控制器对总扰动进行初步估计所产生的估计误差,采用BP神经网络控制器逼近对总扰动的估计误差,从而确定对总扰动的精确估计;
步骤4:根据步骤2和3最终得到的总扰动的精确估计对***进行补偿,实现对压电陶瓷型FTS的精确跟踪。
可选的,所述步骤1中针对压电作动器迟滞非线性,采用Backlash模型描述。
可选的,所述步骤1包括:
压电陶瓷型FTS***的动力学过程描述为:
其中x表示工件切削厚度的波动部分,m、c和k分别是金属切削机的等效质量、阻尼系数和弹簧刚度,ka是压电作动器的弹簧刚度,u是应用于压电作动器的控制输入信号,F1表示切削力的变化量,表示为:
F1=kmv(t) (2)
其中km是常数,v(t)是切削厚度随时间的变化率,表示为:
v(t)=x-αx(t-β) (3)
其中β是连续切削之间的时间间隔,α是重叠因子;
切削力变化与切削厚度变化相关的非线性函数表示为:
F1=H(v)=H(x-αx(t-β)) (4)
其中H(v)由线性函数f(v(t))和有界非线性函数h·v(t)构成,分别表示时滞特性和迟滞特性,即:
H(v)=h·v(t)+f(v(t)) (5)
其中,h表示压电作动器迟滞非线性;
因此,将式(3)代入式(1)可得:
其中,状态时滞β1为主轴旋转产生的时间延迟,当主轴转速Ω变化时,Ω和β1满足如下关系:
时变时滞β2为0到1之间均匀分布的随机数,采样时间为0.1s;
其中Γ表示FTS***内部未建模动态非线性;
变换后的FTS***表示为:
式(9)中,y表示FTS***输出的跟踪信号。
可选的,所述步骤2包括:
针对式(9)描述的FTS***,将***所有不确定性视为总扰动x3,即令:
采用线性自抗干扰控制器对总扰动x3进行初步估计,包括:
将式(9)定义为一个扩展状态空间表达式:
其中
将x3当作一个状态变量,通过LADRC线性扩张状态观测器实时估计x3;
LADRC线性扩张状态观测器为:
其中η=[η1 η2 η3]T,θ=[θ1 θ2 θ3]T;
对公式(13)做拉氏变换,得到特征方程:
λ0(s)=s3+θ1s2+θ2s+θ3=(s+ω0)3
式(9)描述的FTS***被简化为
其中,η3表示采用线性自抗干扰控制器对总扰动x3进行初步估计得到的初步估计值。
可选的,所述步骤3包括:
定义跟踪误差:
e=x-xd
其中,xd表示期望信号;
定义滤波跟踪误差:
γ=[ΛT 1]e (15)
对式(15)求导得:
其中,
采用伪控制反馈线性化技术,定义输入信号u为:
其中φ为BP神经网络的伪控制输入;
利用FTS***的历史输入和输出数据作为数据集训练BP神经网络,以确定BP神经网络的伪控制输入φ,从而确定BP神经网络的权值自适应律;
利用确定了权值自适应律的BP神经网络实现对于对总扰动的估计误差逼近,从而确定对总扰动的精确估计。
可选的,所述BP神经网络中隐层节点数为30;
权值自适应律为:
其中P=PT>0和Q=QT>0是任意常数矩阵;
xbp为神经网络的输入向量,σ(·)是隐层神经元的sigmod激励函数。
本申请还提供一种压电陶瓷型FTS,包括主轴、工作间、切刀以及压电作动器,其中切刀通过弹簧与压电作动器相连;所述压电陶瓷型FTS采用上述基于复合自抗扰控制的FTS控制方法进行控制。
本发明有益效果是:
(1)本申请在构建FTS***模型时,引入了完全未知的迟滞非线性函数、状态时滞和时变时滞的双时滞***以及内部未建模动态非线性,针对所构建的FTS***模型设计的控制器将迟滞时滞及***内部的动态非线性性视为对FTS***的总扰动,由线性观测器实时估计,并利用自适应神经网络逼近未知的非线性误差函数,实现迟滞时滞补偿功能。从而较好描述精密切削过程中的颤振现象,解决了现有FTS模型无法解释和预测切削颤振导致的幅值失真现象、稳定振幅与切削参数的关系、跳变现象对切削过程稳定性的影响及外界干扰对切削精度影响等问题。
(2)本申请采用两种控制方案结合形成的复合控制策略,无需精确模型就能估计和补偿***的未知干扰,不仅可以实现高精度控制,而且具有较强的抗干扰能力,选择带宽作为复合控制器性能的度量,减少了需要整定的控制器参数,鲁棒性更好,具有较高的工程应用价值。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是金属切削过程示意图。
图2是本发明一个实施例中提及的Backlash迟滞模型示意图。
图3是本发明提供的快速道具伺服***总扰动图示意图。
图4是本发明提供的复合自抗扰控制框图。
图5是本发明一个实施例提供的不同带宽下线性自抗扰跟踪曲线图。
图6是本发明一个实施例提供的不同带宽下线性自抗扰跟踪曲线误差图。
图7是本发明一个实施例提供的采用线性自抗扰控制方法和本申请提出的复合控制方法在不同带宽下跟踪曲线图。
图8是本发明一个实施例提供的采用线性自抗扰控制方法和本申请提出的复合控制方法在不同带宽下跟踪曲线误差图。
图9是本发明一个实施例提供的现有控制方法中只采用基于神经网络、只采用线性自抗扰和采用本申请提供的复合控制的时变时滞输出信号对比图。
图10是本发明一个实施例提供的现有控制方法中只采用基于神经网络、只采用线性自抗扰和采用本申请提供的复合控制的迟滞环曲线对比图。
图11是本发明提供的现有控制方法中只采用基于神经网络、只采用线性自抗扰和采用本申请提供的复合控制的跟踪曲线对比图。
图12是本发明提供的现有控制方法中只采用基于神经网络、只采用线性自抗扰和采用本申请提供的复合控制的跟踪误差对比图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
本申请中涉及的名词解释如下:
死区宽度,在控制***中,某些执行机构如果动作频繁,会导致小幅震荡,造成严重的机械磨损。从控制要求来说,很多***又允许被控量在一定范围内存在误差。而允许被控量的误差大小,被称为死区宽度。
时变时滞:切刀横向给进产生的时滞。
状态时滞:主轴旋转产生的时间延迟。
迟滞非线性:压电陶瓷等智能材料固有特性,输入输出表现为迟滞环。
未建模动态非线性:压电陶瓷型FTS中除时变时滞、状态时滞以及迟滞非线性之外的其他未能明确的动态非线性总和。
图1为压电陶瓷型FTS的金属切削过程示意图,压电陶瓷型FTS包括主轴、工作间、切刀,压电作动器等部件,切刀通过弹簧与压电作动器相连;车削过程中,主轴通过减震弹簧带动工件轴向旋转,切刀在压电作动器带动下,以与主轴转速相关联的频率沿工件径向或轴向做往复进给运动,从而加工出工件的端面。
实施例一:
本实施例提供一种基于复合自抗扰控制的FTS控制方法,
步骤1:基于主轴旋转下的状态时滞β1、压电作动器迟滞非线性以及内部未建模动态非线性的结合,以及切刀横向给进产生的时变时滞β2建立FTS***模型;所述内部未建模动态非线性指压电陶瓷型FTS中除时变时滞、状态时滞以及迟滞非线性之外的其他未能明确的动态非线性总和;
步骤2:针对步骤1建立的FTS***模型,将***中主轴旋转下的状态时滞β1、压电作动器迟滞非线性以及内部未建模动态非线性的结合,以及切刀横向给进产生的时变时滞β2视为总扰动,采用线性自抗干扰控制器对总扰动进行初步估计;
步骤3:针对步骤2采用线性自抗干扰控制器对总扰动进行初步估计所产生的估计误差,采用BP神经网络控制器逼近对总扰动的估计误差,从而确定对总扰动的精确估计;
步骤4:根据步骤2和3最终得到的总扰动的精确估计对***进行补偿,实现对压电陶瓷型FTS的精确跟踪。
实施例二:
本实施例提供一种基于复合自抗扰控制的FTS控制方法,所述方法包括:
Step1:考虑主轴旋转下的状态时滞β1、压电作动器迟滞非线性以及内部未建模动态非线性结合的情况,同时又引入了再生效应下的时变时滞β2建立FTS***模型;
与现有的FTS***模型只考虑状态时滞β1和时变时滞β2不同,本申请在对FTS***建模时还考虑了压电作动器迟滞非线性与内部未建模动态非线性;且本申请在建模时,以一定数量的Backlash模型来仿真迟滞非线性h。
建模过程如下:
首先,FTS***的动力学过程描述为
其中x表示工件切削厚度的波动部分,m、c和k分别是金属切削机的等效质量、阻尼系数和弹簧刚度,ka是压电作动器的弹簧刚度,u是应用于压电作动器的控制输入,F1表示刀具切削力与切削厚度的关系,表示为
F1=kmv(t) (2)
其中km是常数,v(t)是切削厚度随时间的变化率,表示为:
v(t)=x-αx(t-β) (3)
其中β是连续切削之间的时间间隔,α是重叠因子。切削力取决于当前切削厚度,而当前车削是位于前次切削的表面进行,切削力随切屑厚度波动从而产生再生颤振。机床的切削力变化与切削厚度变化相关的非线性函数可以表示为:
F1=H(v)=H(x-αx(t-β)) (4)
其中H(v)由线性函数和有界非线性函数构成,分别表示迟滞特性和时滞特性,即:
H(v)=h·v(t)+f(v(t)) (5)
因此,将式(3)代入式(1)可得:
其中,状态时滞β1为主轴旋转产生的时间延迟,当主轴转速Ω变化时,Ω和β1满足如下关系:
时变时滞β2为0到1之间均匀分布的随机数,采样时间为0.1s。h表示未知的迟滞非线性。本申请采用一定数量的Backlash模型来仿真迟滞非线性h,其数学模型表示为如果并且τ=mu-md+,或者并且τ=mu-md_,则其他情况则为0。
Backlash非线性模型是一个一阶速度驱动的动态***,其中u和是输入,τ是***的状态,d是模型的死区宽度。当u的运动方向改变,τ的运动方向不会立即改变,而是会滞后一段时间,由此产生迟滞现象。Backlash迟滞模型如附图2所示。
其中Γ表示FTS***内部未建模动态非线性,变换后的金属切削***可表示为
式(9)所示的本发明构建的FTS***模型主要考虑了主轴旋转下的状态时滞、压电作动器迟滞非线性以及内部未建模动态非线性结合的情况,同时又引入了再生效应下的时变时滞。该模型描述了切削力变化与切屑厚度变化之间的非线性关系,迟滞和双时滞的结合对控制器的设计提出了更高的要求。
Step2:针对步骤1建立的FTS***模型,将***中主轴旋转下的状态时滞β1、压电作动器迟滞非线性以及内部未建模动态非线性的结合,以及切刀横向给进产生的时变时滞β2视为总扰动,设计复合自抗扰器中的线性自抗干扰控制器LADRC,采用线性自抗干扰控制器对总扰动进行初步估计;
线性自抗干扰控制器LADRC将被控对象视为级联积分模型,将所有其他对象信息和外部扰动视为广义扰动,利用线性扩展状态观测器估计广义扰动,并将其纳入线性状态反馈控制律,以快速抑制扰动,快速刀具伺服***总扰动图如附图3所示。
对于式(9)描述的FTS***,将***所有不确定性视为总扰动x3,即令:
LADRC的核心是将***的不确定性当做总扰动来估计和补偿,如果能够实时估计x3,那么就可以主动提取扰动信息并抵消扰动。
式(9)可以被定义为一个扩展状态空间表达式:
其中
将x3当作一个状态变量,通过LADRC线性扩张状态观测器实时估计x3;
LADRC线性扩张状态观测器为:
其中η=[η1 η2 η3]T,θ=[θ1 θ2 θ3]T;
对公式(13)做拉氏变换,得到特征方程:
λ0(s)=s3+θ1s2+θ2s+θ3=(s+ω0)3
式(9)描述的金属切削***可以被简化为
现有技术中,通常将观测器带宽假设为足够大,因此将式(14)简化为:
但实际应用中,观测器带宽有限,因此,本申请后续采用BP神经网络控制器逼近对总扰动的估计误差,即x3-η3,进而进行补充,以实现高精度跟踪控制。
Step3针对步骤2采用线性自抗干扰控制器对总扰动进行初步估计所产生的估计误差,采用BP神经网络控制器逼近对总扰动的估计误差,从而确定对总扰动的精确估计;
定义跟踪误差:
e=x-xd
其中,xd表示期望信号;
定义滤波跟踪误差:
γ=[ΛT 1]e (15)
通过选择适当的系数向量Λ=[f1,f1]T,使得当γ→0时,e→0。
对式(15)求导可得:
其中,
采用伪控制反馈线性化技术,定义输入信号:
逼近误差表示为:
设计伪控制输入为:
其中δ是自适应鲁棒项。将式(20)和(22)代入式(16)可得
自适应BP神经网络能够最大限度地逼近复杂非线性关系并自动修正神经网络参数。其处理结构是并行的,不仅运算速度快,而且容错性强。
利用FTS***的历史输入和输出数据作为数据集训练BP神经网络,以确定BP神经网络的伪控制输入φ,从而确定BP神经网络的权值自适应律;
利用确定了权值自适应律的BP神经网络来逼近估计误差:
其中G=[tij]T和K=[wjb]T为BP神经网络的两个可调权值,
为神经网络的输入向量,σ(·)是隐层神经元的sigmod激励函数,逼近误差ι(xbp)是有界的,|ι(xbp)|≤ιN,ιN>0是一个未知的常数。神经网络的输出可以表示为:
其中P=PT>0和Q=QT>0是任意常数矩阵。基于神经网络的复合自抗扰控制框图如附图4所示。
Step4 FTS***颤振抑制控制
用50个backlash模型构成Backlash-like迟滞模型,死区宽度分布为1/50-50/50。状态时滞Ω设定每两秒转速变化为[1 305 610 560 760 560 800 1120],时变时滞设定为0到1之间的随机数,采样时间为0.1s。神经网络的输入向量为隐层节点数为30,初始权值设置为零,神经网络的权值通过自适应律动态调整,P=8IG,Q=5Ik,IG,Ik均为单位矩阵。使用的输入信号为xd(t)=0.1π(sin 2t-0.1cos t)。
为了证明本申请提出的复合控制方法的控制效果,本申请通过仿真实验对现有技术中采用线性自抗扰控制和本申请所提供的复合控制的跟踪特性进行了对比分析:
仿真实验中,观测器带宽ω0分别设定为40、70和100。传统线性自抗扰采用PD控制,因此还需设置两个可调参数kp和kd,为三个不同带宽的线性自抗扰控制器设置三种依带宽调整的PD参数,即kp40=64,kd40=16,kp70=196,kd70=28,kp100=400,kd100=40。
下面结合仿真图进行说明:
附图5显示了不同带宽下线性自抗扰跟踪曲线图,附图6显示了不同带宽下线性自抗扰跟踪曲线误差图。根据图5和图6可以明显看出现有技术中的线性自抗扰控制器的跟踪精度受带宽影响较大,带宽ω0=100时的跟踪精度远大于ω0=40时的跟踪精度,即带宽越大,跟踪精度越高。
附图7为线性自抗扰控制和本申请提出的复合控制在不同带宽下跟踪曲线图,附图8显示线性自抗扰控制和复合控制在不同带宽下跟踪曲线误差图。从图8中可以看出,本申请所提出的复合自抗扰控制方案受带宽影响较小,带宽ω0=100时的跟踪精度和ω0=70时的跟踪精度几乎相同,这是因为当观测器带宽受限时,本申请所提出的复合自抗扰控制方案中的神经网络能够及时对估计误差进行补偿,实现高精度跟踪控制,从而不受带宽的影响。
附图9分别为现有技术中分别采用线性自抗扰、神经网络和本申请所提出的复合控制的时变时滞输出信号图,附图10分别为线性自抗扰、神经网络和复合控制形成的迟滞环曲线,其输入数据均由***内部自身数据产生。从图中可以看到本申请所构建的新型FTS模型有效描述了非对称迟滞非线性特性。
附图11为现有技术中分别采用线性自抗扰、神经网络和本申请所提出的复合控制方案的跟踪曲线图,附图12为线性自抗扰、神经网络和复合控制方案的跟踪误差比较图。从图12中可以看到虽然线性自抗扰控制、神经网络控制和复合控制方案对于迟滞时滞快速刀具伺服***均能实现有效跟踪控制,但是本申请所提出的复合控制的初始误差最小且相较线性自抗扰控制与神经网络控制具有最小的跟踪误差,最高的跟踪精度,进而能够使得所加工出的工具质量和效率得到提升。
本发明实施例中的部分步骤,可以利用软件实现,相应的软件程序可以存储在可读取的存储介质中,如光盘或硬盘等。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于复合自抗扰控制的FTS控制方法,所述方法应用于压电陶瓷型FTS中,所述压电陶瓷型FTS包括主轴、工作间、切刀以及压电作动器,其中切刀通过弹簧与压电作动器相连;其特征在于,所述方法包括:
步骤1:基于主轴旋转下的状态时滞β1、压电作动器迟滞非线性以及内部未建模动态非线性的结合,以及切刀横向给进产生的时变时滞β2建立FTS***模型;所述内部未建模动态非线性指压电陶瓷型FTS中除时变时滞、状态时滞以及迟滞非线性之外的其他未能明确的动态非线性总和;
步骤2:针对步骤1建立的FTS***模型,将***中主轴旋转下的状态时滞β1、压电作动器迟滞非线性以及内部未建模动态非线性的结合,以及切刀横向给进产生的时变时滞β2视为总扰动,采用线性自抗干扰控制器对总扰动进行初步估计;
步骤3:针对步骤2采用线性自抗干扰控制器对总扰动进行初步估计所产生的估计误差,采用BP神经网络控制器逼近对总扰动的估计误差,从而确定对总扰动的精确估计;
步骤4:根据步骤2和3最终得到的总扰动的精确估计对***进行补偿,实现对压电陶瓷型FTS的精确跟踪。
2.根据权利要求1所述的方法,其特征在于,所述步骤1中针对压电作动器迟滞非线性,采用Backlash模型描述。
3.根据权利要求1所述的方法,其特征在于,所述步骤1包括:
压电陶瓷型FTS***的动力学过程描述为:
其中x表示工件切削厚度的波动部分,m、c和k分别是金属切削机的等效质量、阻尼系数和弹簧刚度,ka是压电作动器的弹簧刚度,u是应用于压电作动器的控制输入信号,F1表示切削力的变化量,表示为:
F1=kmv(t) (2)
其中km是常数,v(t)是切削厚度随时间的变化率,表示为:
v(t)=x-αx(t-β) (3)
其中β是连续切削之间的时间间隔,α是重叠因子;
切削力变化与切削厚度变化相关的非线性函数表示为:
F1=H(v)=H(x-αx(t-β)) (4)
其中H(v)由线性函数f(v(t))和有界非线性函数h·v(t)构成,分别表示时滞特性和迟滞特性,即:
H(v)=h·v(t)+f(v(t)) (5)
其中,h表示压电作动器迟滞非线性;
因此,将式(3)代入式(1)可得:
其中,状态时滞β1为主轴旋转产生的时间延迟,当主轴转速Ω变化时,Ω和β1满足如下关系:
时变时滞β2为0到1之间均匀分布的随机数,采样时间为0.1s;
其中Γ表示FTS***内部未建模动态非线性;
变换后的FTS***表示为:
式(9)中,y表示FTS***输出的跟踪信号。
4.根据权利要求3所述的方法,其特征在于,所述步骤2包括:
针对式(9)描述的FTS***,将***所有不确定性视为总扰动x3,即令:
采用线性自抗干扰控制器对总扰动x3进行初步估计,包括:
将式(9)定义为一个扩展状态空间表达式:
其中
将x3当作一个状态变量,通过LADRC线性扩张状态观测器实时估计x3;
LADRC线性扩张状态观测器为:
其中η=[η1 η2 η3]T,θ=[θ1 θ2 θ3]T;
对公式(13)做拉氏变换,得到特征方程:
λ0(s)=s3+θ1s2+θ2s+θ3=(s+ω0)3
式(9)描述的FTS***被简化为
其中,η3表示采用线性自抗干扰控制器对总扰动x3进行初步估计得到的初步估计值。
8.一种压电陶瓷型FTS,包括主轴、工作间、切刀以及压电作动器,其中切刀通过弹簧与压电作动器相连;其特征在于,所述压电陶瓷型FTS采用权利要求1-7任一所述的基于复合自抗扰控制的FTS控制方法进行控制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210321697.7A CN114609975B (zh) | 2022-03-25 | 2022-03-25 | 一种基于复合自抗扰控制的fts控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210321697.7A CN114609975B (zh) | 2022-03-25 | 2022-03-25 | 一种基于复合自抗扰控制的fts控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114609975A true CN114609975A (zh) | 2022-06-10 |
CN114609975B CN114609975B (zh) | 2024-04-30 |
Family
ID=81867467
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210321697.7A Active CN114609975B (zh) | 2022-03-25 | 2022-03-25 | 一种基于复合自抗扰控制的fts控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114609975B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107132759A (zh) * | 2017-02-24 | 2017-09-05 | 沈阳工业大学 | 一种直线电机驱动fts基于eso滑模改进重复控制方法 |
CN110045677A (zh) * | 2019-02-27 | 2019-07-23 | 沈阳工业大学 | 一种麦克斯韦快刀伺服轨迹跟踪控制方法 |
CN112379595A (zh) * | 2020-11-14 | 2021-02-19 | 西北工业大学 | 非线性严格反馈切换***复合干扰学习控制方法 |
-
2022
- 2022-03-25 CN CN202210321697.7A patent/CN114609975B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107132759A (zh) * | 2017-02-24 | 2017-09-05 | 沈阳工业大学 | 一种直线电机驱动fts基于eso滑模改进重复控制方法 |
CN110045677A (zh) * | 2019-02-27 | 2019-07-23 | 沈阳工业大学 | 一种麦克斯韦快刀伺服轨迹跟踪控制方法 |
CN112379595A (zh) * | 2020-11-14 | 2021-02-19 | 西北工业大学 | 非线性严格反馈切换***复合干扰学习控制方法 |
Non-Patent Citations (4)
Title |
---|
WU, D等: "Design of a normal stress electromagnetic fast linear actuator", 《IEEE TRANSACTIONS ON MAGNETICS》, 30 April 2010 (2010-04-30), pages 1007 - 1014, XP011297658, DOI: 10.1109/TMAG.2009.2036606 * |
吴丹等: "快速刀具伺服***自抗扰控制的研究与实践", 《控制理论与应用》, vol. 30, no. 12, 31 December 2013 (2013-12-31), pages 1534 - 1542 * |
桂思齐: "快速刀具伺服的控制***搭建与实验", 《中国优秀硕士学位论文全文数据库 工程科技I辑》, 15 January 2021 (2021-01-15) * |
王丽君等: "时滞***的自抗扰控制综述", 《控制理论与应用》, vol. 30, no. 12, 31 December 2013 (2013-12-31), pages 1521 - 1533 * |
Also Published As
Publication number | Publication date |
---|---|
CN114609975B (zh) | 2024-04-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Atkeson et al. | Robot trajectory learning through practice | |
CN105159069B (zh) | 一种压电陶瓷作动器的位移控制方法 | |
Hajikolaei et al. | Spindle speed variation and adaptive force regulation to suppress regenerative chatter in the turning process | |
CN107577146B (zh) | 基于摩擦整体逼近的伺服***的神经网络自适应控制方法 | |
Kleinwort et al. | Comparison of different control strategies for active damping of heavy duty milling operations | |
Heng et al. | Design of super twisting algorithm for chattering suppression in machine tools | |
CN113641181B (zh) | 基于在线性能评估飞行器增益自适应姿态控制方法和*** | |
CN114278695B (zh) | 一种基于磁流变阻尼器的薄壁件加工振动半主动控制方法 | |
CN114114928A (zh) | 一种压电微定位平台的固定时间自适应事件触发控制方法 | |
Wang et al. | Study on application of model reference adaptive control in fast steering mirror system | |
Qian et al. | Youla parameterized adaptive vibration control against deterministic and band-limited random signals | |
CN107037726B (zh) | 一种非同元分数阶***滑膜干扰观测器设计方法 | |
Finney et al. | A pole-assignment controller for an electrohydraulic cylinder drive | |
CN114609975B (zh) | 一种基于复合自抗扰控制的fts控制方法 | |
CN109702554B (zh) | 一种自适应振动塑形铣削颤振抑制方法 | |
Li et al. | Multi-modal vibration control for all-clamped plate subjected to periodic disturbances by ESO-based frequency-shaped LQR | |
Fischer et al. | Simulation-based stability analysis of a thin-walled cylinder during turning with improvements using an adaptronic turning chisel | |
Moradi et al. | Sliding mode control of machining chatter in the presence of tool wear and parametric uncertainties | |
CN113110021B (zh) | 一种用于伺服***辨识与控制器设计的方法 | |
CN112114517A (zh) | 3d打印温度控制算法 | |
Wu et al. | Adaptive CNC machining process optimization of near-net-shaped blade based on machining error data flow control | |
Ke et al. | Tangential velocity tracking-based task coordinate frame approach for contouring control of biaxial motion systems | |
Eze et al. | Servo position control in hard disk drive of a computer using MRAC integrating PID algorithm | |
Taktak-Meziou et al. | A prediction-based optimal gain selection in RISE feedback control for hard disk drive | |
CN112147898B (zh) | 一种仅依赖控制方向信息的刚体***抗扰控制方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |