[go: up one dir, main page]

CN116088312A - Multi-axis cooperative control method for dual-drive gantry platform of chip mounter - Google Patents

Multi-axis cooperative control method for dual-drive gantry platform of chip mounter Download PDF

Info

Publication number
CN116088312A
CN116088312A CN202310059796.7A CN202310059796A CN116088312A CN 116088312 A CN116088312 A CN 116088312A CN 202310059796 A CN202310059796 A CN 202310059796A CN 116088312 A CN116088312 A CN 116088312A
Authority
CN
China
Prior art keywords
axis
chip mounter
gantry platform
matrix
gantry
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.)
Pending
Application number
CN202310059796.7A
Other languages
Chinese (zh)
Inventor
于兴虎
史朋威
龚见素
赵宇哲
孙昊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Ningbo Yitang Intelligent Technology Co ltd
Original Assignee
Ningbo Yitang Intelligent Technology Co ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Ningbo Yitang Intelligent Technology Co ltd filed Critical Ningbo Yitang Intelligent Technology Co ltd
Priority to CN202310059796.7A priority Critical patent/CN116088312A/en
Publication of CN116088312A publication Critical patent/CN116088312A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02PCLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
    • Y02P90/00Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
    • Y02P90/02Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Feedback Control In General (AREA)

Abstract

贴片机双驱龙门平台的多轴协同控制方法,涉及贴片机龙门平台的运动控制技术领域。本发明是为了解决现有贴片机双驱龙门平台控制系统中,由于无法保证电机一致性、机械横梁同步运动与负载定位运动存在耦合动力学关系、以及外界干扰都会对系统的控制造成误差,导致系统不确定性风险高且控制难度大的问题。本发明所述的贴片机双驱龙门平台的多轴协同控制方法,分别设计期望模型补偿控制器、鲁棒反馈控制器和神经网络控制器,并将设计的控制器输出之和作为贴片机龙门平台的总控制信号,利用该总控制信号对贴片机龙门平台的X轴、Y1轴和Y2轴进行多轴协同控制。

Figure 202310059796

The invention discloses a multi-axis cooperative control method for a double-drive gantry platform of a placement machine, and relates to the technical field of motion control of the gantry platform of a placement machine. The present invention is to solve the problem that in the existing dual-drive gantry platform control system of the placement machine, due to the inability to ensure the consistency of the motor, the coupling dynamic relationship between the synchronous movement of the mechanical beam and the positioning movement of the load, and external interference will cause errors in the control of the system, It leads to problems with high risk of system uncertainty and difficult control. The multi-axis cooperative control method of the dual-drive gantry platform of the placement machine according to the present invention, respectively designs the expected model compensation controller, the robust feedback controller and the neural network controller, and uses the sum of the designed controller outputs as the placement The general control signal of the gantry platform of the placement machine is used to perform multi-axis coordinated control of the X-axis, Y1-axis and Y2-axis of the gantry platform of the placement machine.

Figure 202310059796

Description

贴片机双驱龙门平台的多轴协同控制方法Multi-axis collaborative control method for dual-drive gantry platform of placement machine

技术领域Technical Field

本发明属于贴片机龙门平台的运动控制技术领域,尤其涉及多轴协同控制。The invention belongs to the technical field of motion control of a gantry platform of a placement machine, and in particular relates to multi-axis coordinated control.

背景技术Background Art

双驱龙门平台是贴片机中典型的驱动结构,双驱龙门平台具有推力大,负载强,精度高等特点,在现代工业高端制造领域愈加受到关注和应用。在贴片机龙门平台系统中,机械横梁横跨在平台两端并由两直线电机共同驱动以实现Y轴方向进给运动,移动工作台在另一电机驱动下沿横梁所在X轴方向进行定位运动,通过两轴之间的密切配合与协调运动共同完成复杂的生产任务。需要注意的是,相比于工作台的单驱直线运动控制,XY多轴协同运动控制显得更为复杂和重要,其控制精度直接影响着贴片机龙门平台系统的运动性能和加工工艺。因此,研究贴片机双驱龙门平台多轴协同控制问题具有重要的实际意义。The dual-drive gantry platform is a typical drive structure in the placement machine. The dual-drive gantry platform has the characteristics of large thrust, strong load, and high precision. It has received more and more attention and application in the field of modern industrial high-end manufacturing. In the gantry platform system of the placement machine, the mechanical beam spans across the two ends of the platform and is driven by two linear motors to realize the Y-axis feeding movement. The mobile worktable is driven by another motor to perform positioning movement along the X-axis direction where the beam is located. The complex production tasks are completed through close cooperation and coordinated movement between the two axes. It should be noted that compared with the single-drive linear motion control of the worktable, the XY multi-axis coordinated motion control is more complex and important, and its control accuracy directly affects the motion performance and processing technology of the gantry platform system of the placement machine. Therefore, it is of great practical significance to study the multi-axis coordinated control problem of the dual-drive gantry platform of the placement machine.

在整个控制系统中,多个直线电机驱动方式能够使得系统获得较大推力和速度,但同时也给系统增加了不确定性风险和控制难度。这主要是由于两方面因素,其一,现实条件下很难保证三个直线电机具备完全相同的机械参数和机电特性;其二,机械横梁同步运动与负载定位运动相互影响且具有复杂的耦合动力学关系。另外,贴片机龙门平台系统中不可避免会存在参数不确定性及外界未知非线性扰动的问题,例如负载变化,推力波动,测量噪声等,这些都很容易造成较大的运动误差和系统震荡。严重时甚至会导致机械设备损坏,威胁操作人员的生命安全。In the entire control system, multiple linear motor drive modes can enable the system to obtain greater thrust and speed, but at the same time it also increases the uncertainty risk and control difficulty of the system. This is mainly due to two factors. First, it is difficult to ensure that the three linear motors have exactly the same mechanical parameters and electromechanical characteristics under realistic conditions; second, the synchronous motion of the mechanical beam and the load positioning motion affect each other and have a complex coupling dynamic relationship. In addition, there will inevitably be problems of parameter uncertainty and external unknown nonlinear disturbances in the SMT gantry platform system, such as load changes, thrust fluctuations, measurement noise, etc., which can easily cause large motion errors and system oscillations. In severe cases, it may even cause damage to mechanical equipment and threaten the life safety of operators.

发明内容Summary of the invention

本发明是为了解决现有贴片机双驱龙门平台控制系统中,由于无法保证电机一致性、机械横梁同步运动与负载定位运动存在耦合动力学关系、以及外界干扰都会对系统的控制造成误差,导致系统不确定性风险高且控制难度大的问题,现提供贴片机双驱龙门平台的多轴协同控制方法。The present invention aims to solve the problems in the existing control system of the dual-drive gantry platform of a placement machine, in which the motor consistency cannot be guaranteed, the synchronous movement of the mechanical crossbeam and the load positioning movement have a coupled dynamic relationship, and external interference will cause errors in the control of the system, resulting in high system uncertainty risk and great control difficulty. A multi-axis collaborative control method for the dual-drive gantry platform of a placement machine is now provided.

贴片机双驱龙门平台的多轴协同控制方法,分别设计期望模型补偿控制器、鲁棒反馈控制器和神经网络控制器,并将设计的控制器输出之和作为贴片机龙门平台的总控制信号,利用该总控制信号对贴片机龙门平台的X轴、Y1轴和Y2轴进行多轴协同控制;The multi-axis cooperative control method of the dual-drive gantry platform of the placement machine is designed. The expected model compensation controller, the robust feedback controller and the neural network controller are designed respectively, and the sum of the designed controller outputs is used as the total control signal of the gantry platform of the placement machine. The total control signal is used to perform multi-axis cooperative control on the X-axis, Y1-axis and Y2-axis of the gantry platform of the placement machine.

所述期望模型补偿控制器表达式为:The desired model compensation controller expression is:

Figure BDA0004061058310000011
Figure BDA0004061058310000011

其中,υa为期望模型补偿控制器的输出,Φd为含有期望信息的回归量矩阵,

Figure BDA0004061058310000012
Figure BDA0004061058310000013
的估计值,
Figure BDA0004061058310000028
为贴片机龙门平台系统耦合动力学模型中不确定参数构成的系数矩阵,Among them, υ a is the output of the expected model compensation controller, Φ d is the regression matrix containing the expected information,
Figure BDA0004061058310000012
for
Figure BDA0004061058310000013
The estimated value of
Figure BDA0004061058310000028
is the coefficient matrix of uncertain parameters in the coupled dynamics model of the SMT gantry platform system.

Figure BDA0004061058310000021
Figure BDA0004061058310000021

其中,Mhk为归一化后贴片机中移动工作台的质量,Bxk为归一化后贴片机龙门平台X轴处的粘滞摩擦系数,Axk为归一化后贴片机龙门平台X轴处的库仑摩擦系数,Myk为归一化后贴片机龙门平台中横梁的总质量,Byk为归一化后贴片机龙门平台Y1轴和Y2轴导轨处的粘滞摩擦系数之和,Ayk为归一化后贴片机龙门平台Y1轴和Y2轴导轨处的库仑摩擦系数之和,Jk为归一化后的转动惯量,Kαk为归一化后的等效旋转刚度,Bsk为归一化后贴片机龙门平台Y2轴和Y1轴导轨处的粘滞摩擦系数之差,Ask为归一化后贴片机龙门平台Y2轴和Y1轴导轨处的库仑摩擦系数之差,

Figure BDA0004061058310000022
Figure BDA0004061058310000023
分别为归一化后贴片机龙门平台X轴和Y轴动力学中非线性扰动的常值,
Figure BDA0004061058310000024
为归一化后贴片机龙门平台旋转动力学中非线性扰动的常值;Among them, M hk is the mass of the moving workbench in the SMT machine after normalization, B xk is the viscous friction coefficient at the X-axis of the SMT gantry platform after normalization, A xk is the Coulomb friction coefficient at the X-axis of the SMT gantry platform after normalization, Myk is the total mass of the beam in the SMT gantry platform after normalization, Byk is the sum of the viscous friction coefficients at the Y1-axis and Y2-axis guide rails of the SMT gantry platform after normalization, A yk is the sum of the Coulomb friction coefficients at the Y1-axis and Y2-axis guide rails of the SMT gantry platform after normalization, J k is the normalized moment of inertia, K αk is the normalized equivalent rotational stiffness, B sk is the difference in viscous friction coefficients at the Y2-axis and Y1-axis guide rails of the SMT gantry platform after normalization, A sk is the difference in Coulomb friction coefficients at the Y2-axis and Y1-axis guide rails of the SMT gantry platform after normalization,
Figure BDA0004061058310000022
and
Figure BDA0004061058310000023
are the normalized constants of the nonlinear disturbances in the X-axis and Y-axis dynamics of the SMT gantry platform,
Figure BDA0004061058310000024
is the normalized constant value of the nonlinear disturbance in the rotational dynamics of the gantry platform of the placement machine;

所述鲁棒反馈控制器表达式为:The robust feedback controller expression is:

υe=-MqZ-1G-Kds-Keε,υ e =-M q Z -1 GK d sK e ε,

其中,υe为鲁棒反馈控制器的输出,Mq为归一化后贴片机龙门平台的惯性矩阵,s为滑模变量,Z=diag[z1,z2,z3],

Figure BDA0004061058310000025
Where, υ e is the output of the robust feedback controller, M q is the normalized inertia matrix of the gantry platform of the placement machine, s is the sliding mode variable, Z = diag [z 1 , z 2 , z 3 ],
Figure BDA0004061058310000025

Si为误差转换函数且有: Si is the error transfer function and has:

Figure BDA0004061058310000026
Figure BDA0004061058310000026

ε=[ε123]为无约束的转换误差,中间变量ri=ei(t)/ρi(t),ei(t)为t时刻的误差矢量,ei(0)为误差矢量初始值,

Figure BDA0004061058310000027
为超调量调节系数,ε=[ε 123 ] is the unconstrained conversion error, the intermediate variable riei (t)/ ρi (t), ei (t) is the error vector at time t, ei (0) is the initial value of the error vector,
Figure BDA0004061058310000027
is the overshoot adjustment coefficient,

ρi为性能函数且表达式如下:ρ i is the performance function and is expressed as follows:

Figure BDA0004061058310000029
Figure BDA0004061058310000029

ρi0为稳态误差初始值,ρi∞为稳态误差最大值,

Figure BDA0004061058310000037
为误差收敛速率,ρ i0 is the initial value of the steady-state error, ρ i∞ is the maximum value of the steady-state error,
Figure BDA0004061058310000037
is the error convergence rate,

Figure BDA0004061058310000031
λ为控制增益矩阵,Kd为比例反馈矩阵,Ke为鲁棒增益矩阵;
Figure BDA0004061058310000031
λ is the control gain matrix, Kd is the proportional feedback matrix, and Ke is the robust gain matrix;

所述神经网络控制器表达式为:The neural network controller expression is:

Figure BDA0004061058310000032
Figure BDA0004061058310000032

其中,υnn为神经网络控制器的输出,

Figure BDA0004061058310000033
为神经网络权重矢量,χ为神经网络输入矢量,χ=[χ12,...,χm],
Figure BDA0004061058310000034
m为神经网络的神经元总层数,j=1,2,...,m,R(χ)为神经网络多层神经元的总输出,η为鲁棒项系数,sign(·)表示符号函数;Among them, υ nn is the output of the neural network controller,
Figure BDA0004061058310000033
is the neural network weight vector, χ is the neural network input vector, χ=[χ 12 ,...,χ m ],
Figure BDA0004061058310000034
m is the total number of neuron layers in the neural network, j = 1, 2, ..., m, R(χ) is the total output of the multi-layer neurons in the neural network, η is the robust term coefficient, and sign(·) represents the sign function;

所述贴片机龙门平台的总控制信号v=[v1,v2,v3]:The total control signal v of the gantry platform of the placement machine is [v 1 ,v 2 ,v 3 ]:

v=υaennv=υ aenn .

进一步的,根据下式利用贴片机龙门平台的总控制信号实现对贴片机龙门平台X轴、Y1轴和Y2轴的协同控制:Furthermore, the total control signal of the SMT gantry platform is used to realize the coordinated control of the X-axis, Y1-axis and Y2-axis of the SMT gantry platform according to the following formula:

Figure BDA0004061058310000035
Figure BDA0004061058310000035

其中,Sm为贴片机龙门平台横梁两侧电机之间的距离,α为贴片机龙门平台横梁的实际旋转角度,Kxk为归一化后的X轴电机推力常数,Kmy=K2/K1,K1和K2分别为贴片机龙门平台Y1轴和Y2轴电机推力常数,ux、u1和u2分别为X轴、Y1轴和Y2轴的控制输入。Among them, Sm is the distance between the motors on both sides of the crossbeam of the gantry platform of the placement machine, α is the actual rotation angle of the crossbeam of the gantry platform of the placement machine, Kxk is the normalized thrust constant of the X-axis motor, Kmy = K2 / K1 , K1 and K2 are the thrust constants of the Y1-axis and Y2-axis motors of the gantry platform of the placement machine respectively, and ux , u1 and u2 are the control inputs of the X-axis, Y1-axis and Y2-axis respectively.

进一步的,上述贴片机龙门平台系统耦合动力学模型为:Furthermore, the coupled dynamics model of the above-mentioned SMT machine gantry platform system is:

Figure BDA0004061058310000036
Figure BDA0004061058310000036

其中,Cq、Bq、Kq和Aq分别为归一化后的科氏力矩阵、粘滞摩擦矩阵、等效刚度矩阵和库伦摩擦矩阵,q为贴片机龙门平台系统状态矢量矩阵且q=[x(t),y(t),α(t)],x(t)和y(t)分别为t时刻贴片机龙门平台X轴和Y轴的实际位移,α(t)为t时刻贴片机龙门平台系统横梁的实际旋转角度,

Figure BDA0004061058310000041
Figure BDA0004061058310000042
分别为贴片机龙门平台系统轴运动的速度和加速度矢量矩阵,非线性函数
Figure BDA0004061058310000043
Figure BDA0004061058310000044
Figure BDA0004061058310000045
分别为归一化后非线性扰动的常值矩阵和时变矩阵。Among them, Cq , Bq , Kq and Aq are the normalized Coriolis force matrix, viscous friction matrix, equivalent stiffness matrix and Coulomb friction matrix respectively, q is the state vector matrix of the SMT gantry platform system and q=[x(t),y(t),α(t)], x(t) and y(t) are the actual displacements of the X-axis and Y-axis of the SMT gantry platform at time t respectively, α(t) is the actual rotation angle of the beam of the SMT gantry platform system at time t,
Figure BDA0004061058310000041
and
Figure BDA0004061058310000042
They are the velocity and acceleration vector matrices of the axis motion of the SMT gantry platform system, and the nonlinear function
Figure BDA0004061058310000043
Figure BDA0004061058310000044
and
Figure BDA0004061058310000045
are the constant matrix and time-varying matrix of the normalized nonlinear perturbation respectively.

进一步的,上述贴片机龙门平台系统的误差矢量ei(t)的约束条件为:Furthermore, the constraint condition of the error vector e i (t) of the gantry platform system of the above-mentioned placement machine is:

Figure BDA0004061058310000046
Figure BDA0004061058310000046

其中,ei(t)=q-qd=[x(t)-xd(t),y(t)-yd(t),α(t)],Among them, e i (t) = qq d = [x (t)-x d (t), y (t) - y d (t), α (t)],

qd=[xd(t),yd(t),αd(t)]为贴片机龙门平台系统的期望轨迹信号矩阵,xd(t)和yd(t)分别为t时刻贴片机龙门平台X轴和Y轴的期望位移,αd(t)为贴片机龙门平台横梁的期望旋转角度且αd(t)=0。 qd =[ xd (t), yd (t), αd (t)] is the expected trajectory signal matrix of the gantry system of the placement machine, xd (t) and yd (t) are the expected displacements of the X-axis and Y-axis of the gantry at time t, respectively, αd (t) is the expected rotation angle of the crossbeam of the gantry and αd (t)=0.

进一步的,上述

Figure BDA0004061058310000047
Figure BDA0004061058310000048
Furthermore, the above
Figure BDA0004061058310000047
Figure BDA0004061058310000048

Figure BDA0004061058310000049
Figure BDA0004061058310000049

其中,h为贴片机中移动工作台质心的偏移距离,x(t)为贴片机龙门平台系统横梁X轴的实际位移,Sg为贴片机龙门平台Y1轴和Y2轴导轨之间的距离。Among them, h is the offset distance of the center of mass of the mobile worktable in the placement machine, x(t) is the actual displacement of the X-axis of the beam of the gantry platform system of the placement machine, and Sg is the distance between the Y1-axis and Y2-axis guide rails of the gantry platform of the placement machine.

进一步的,贴片机龙门平台系统的回归量

Figure BDA00040610583100000410
包括残差矩阵部分和含有期望信息部分,表达式如下:Furthermore, the regression quantity of the SMT gantry system
Figure BDA00040610583100000410
Including the residual matrix part and the part containing expected information, the expression is as follows:

Figure BDA00040610583100000411
Figure BDA00040610583100000411

其中,

Figure BDA0004061058310000051
为仅含有期望信息的回归量矩阵,
Figure BDA0004061058310000052
为回归量残差矩阵。in,
Figure BDA0004061058310000051
is the regressor matrix containing only the expected information,
Figure BDA0004061058310000052
is the residual matrix of the regressors.

进一步的,上述滑模变量s的表达式为:Furthermore, the expression of the sliding mode variable s is:

Figure BDA0004061058310000053
Figure BDA0004061058310000053

其中,

Figure BDA0004061058310000054
in,
Figure BDA0004061058310000054

进一步的,上述神经网络中第j层神经元的输出R(χj)为:Furthermore, the output R(χ j ) of the j-th layer neuron in the above neural network is:

Figure BDA0004061058310000055
Figure BDA0004061058310000055

其中,cj和bj分比为神经网络核函数的中心坐标矢量和宽度。Among them, cj and bj are the central coordinate vector and width of the neural network kernel function respectively.

进一步的,上述神经网络权重矢量

Figure BDA0004061058310000056
的自适应率
Figure BDA0004061058310000057
为:Furthermore, the above neural network weight vector
Figure BDA0004061058310000056
Adaptive rate
Figure BDA0004061058310000057
for:

Figure BDA0004061058310000058
Figure BDA0004061058310000058

其中,φ学习率更新矩阵。Where φ is the learning rate update matrix.

本发明的有益效果:Beneficial effects of the present invention:

本发明提出了一种适用于贴片机龙门平台的精密多轴协同控制方法,解决了复杂生产环境下,抗干扰能力差,控制性能不能预先设定等难题。本发明所述的贴片机双驱龙门平台的多轴协同控制方法,是基于预设性能技术与自适应神经网络控制方法共同实现的。其中,提前设计预设性能函数以规划系统的期望性能,使用误差变换技术将约束问题进行简化处理;自适应环节中参数估计算法能够实现系统参数快速收敛,然后通过所设计的自适律对模型进行精确补偿;鲁棒反馈模块是基于实际误差信号所设计的控制器,主要用于镇定名义系统;神经网络补偿器用来逼近未建模动态与外部环境中的干扰,进一步提高系统的稳定性。最后经过在龙门运动实验台上进行对比实验,验证了所提的同步控制方法的有效型与优越性。The present invention proposes a precise multi-axis collaborative control method suitable for the gantry platform of a placement machine, which solves the problems of poor anti-interference ability and inability to pre-set control performance in a complex production environment. The multi-axis collaborative control method of the dual-drive gantry platform of the placement machine described in the present invention is implemented based on the preset performance technology and the adaptive neural network control method. Among them, the preset performance function is designed in advance to plan the expected performance of the system, and the constraint problem is simplified using the error transformation technology; the parameter estimation algorithm in the adaptive link can achieve rapid convergence of the system parameters, and then the model is accurately compensated by the designed adaptive law; the robust feedback module is a controller designed based on the actual error signal, which is mainly used to stabilize the nominal system; the neural network compensator is used to approximate the interference in the unmodeled dynamics and the external environment, and further improve the stability of the system. Finally, a comparative experiment was carried out on the gantry motion test bench to verify the effectiveness and superiority of the proposed synchronous control method.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

图1为本发明所述一种预设性能保证的贴片机龙门平台多轴协同控制方法流程图;FIG1 is a flow chart of a multi-axis collaborative control method for a gantry platform of a placement machine with preset performance guarantee according to the present invention;

图2为本发明涉及的贴片机龙门平台运动示意图;FIG2 is a schematic diagram of the movement of the gantry platform of the chip placement machine according to the present invention;

图3为实施例中XY轴跟随期望轨迹时的运动位置变化曲线图;FIG3 is a graph showing a motion position change curve when the XY axis follows a desired trajectory in an embodiment;

图4为实施例中X轴跟随期望轨迹时的跟踪误差变化曲线图;FIG4 is a graph showing a tracking error variation when the X-axis follows a desired trajectory in an embodiment;

图5为实施例中Y轴跟随期望轨迹时的跟踪误差变化曲线图;FIG5 is a graph showing a tracking error variation when the Y-axis follows a desired trajectory in an embodiment;

图6为实施例中Y轴跟随期望轨迹时的旋转角度变化曲线图。FIG. 6 is a curve diagram showing the rotation angle variation when the Y-axis follows the desired trajectory in the embodiment.

具体实施方式DETAILED DESCRIPTION

下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其它实施例,都属于本发明保护的范围。需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。The following will be combined with the accompanying drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work belong to the scope of protection of the present invention. It should be noted that the embodiments of the present invention and the features in the embodiments can be combined with each other without conflict.

现有对龙门控制开展研究大多以单轴定位和双驱同步跟踪为主,并且控制算法在进一步提高系统的动态响应和稳态精度缺乏灵活性和前瞻性。现有的研究在提升系统的动态响应和稳态精度方面存在不足、控制效果差等问题,本实施方式提供一种预设性能保证的精密多轴协同控制方法,能够有效地处理贴片机龙门系统中存在的各种不确定性和未知干扰,同时提高各个轴的瞬态性能和稳态性能,保证系统平稳安全运行。参照图1和图2具体说明本实施方式。Most of the existing research on gantry control focuses on single-axis positioning and dual-drive synchronous tracking, and the control algorithm lacks flexibility and foresight in further improving the dynamic response and steady-state accuracy of the system. Existing research has deficiencies in improving the dynamic response and steady-state accuracy of the system, and poor control effects. This embodiment provides a precise multi-axis collaborative control method with preset performance guarantees, which can effectively handle various uncertainties and unknown interferences in the gantry system of the placement machine, while improving the transient performance and steady-state performance of each axis, ensuring the smooth and safe operation of the system. Refer to Figures 1 and 2 to explain this embodiment in detail.

步骤一、根据贴片机龙门平台几何关系和受力分析,建立耦合动力学模型。Step 1: Establish a coupled dynamics model based on the geometric relationship and force analysis of the placement machine gantry platform.

如图2所示,定义广义坐标系OXY和旋转坐标系

Figure BDA0004061058310000061
各参数坐标的定义为:As shown in Figure 2, the generalized coordinate system OXY and the rotation coordinate system are defined.
Figure BDA0004061058310000061
The definition of each parameter coordinate is:

定义横梁自身质心Gy的广义位置坐标为(0,y),移动平台(包括负载)质心Gx的旋转位置坐标为(x,-h)。根据上述坐标进一步定义整个横梁(包括移动工作台)质心位置为G,横梁在Y1轴侧的光栅反馈坐标为(x1,y1),横梁在Y2轴侧的光栅反馈坐标为(x2,y2)。横梁左右两个电机之间的距离为Sm,左右两个导轨之间的距离为Sg,左右两个光栅之间的距离为Se;整个横梁质心到Y1电机距离为Sm1,到Y1导轨的距离为Sg1,到Y1光栅的距离为Se1;整个横梁质心到Y2电机距离为Sm2,到Y2导轨的距离为Sg2,到Y2光栅的距离为Se2Define the generalized position coordinates of the center of mass G y of the beam itself as (0, y), and the rotational position coordinates of the center of mass G x of the mobile platform (including the load) as (x, -h). Based on the above coordinates, further define the center of mass position of the entire beam (including the mobile workbench) as G, the grating feedback coordinates of the beam on the Y1 axis side as (x 1 , y 1 ), and the grating feedback coordinates of the beam on the Y2 axis side as (x 2 , y 2 ). The distance between the left and right motors of the beam is S m , the distance between the left and right guide rails is S g , and the distance between the left and right gratings is Se ; the distance from the center of mass of the entire beam to the Y1 motor is S m1 , the distance to the Y1 guide rail is S g1 , and the distance to the Y1 grating is Se1 ; the distance from the center of mass of the entire beam to the Y2 motor is S m2 , the distance to the Y2 guide rail is S g2 , and the distance to the Y2 grating is Se2 .

贴片机龙门平台正常工作时,期望工作状态下横梁与左右侧的导轨是垂直的,即横梁与X轴无夹角,但实际工作移动平台的高速运动对横梁形造成一定的冲击影响,横梁两端一旦不同步即会出现旋转形成夹角,在较大内力约束下容易引起系统的震荡失稳。因此,为了全面考虑负载及横梁的因素,本实施方式增加了对横梁旋转内因和轴间耦合的分析。When the SMT gantry platform is working normally, it is expected that the beam is perpendicular to the left and right guide rails, that is, there is no angle between the beam and the X-axis. However, the high-speed movement of the actual working mobile platform has a certain impact on the beam shape. Once the two ends of the beam are not synchronized, they will rotate to form an angle, which is easy to cause the system to oscillate and become unstable under the constraint of large internal forces. Therefore, in order to fully consider the factors of load and beam, this implementation adds the analysis of the internal cause of beam rotation and the coupling between axes.

贴片机龙门平台横梁的实际旋转角度为α,考虑到实际情况α≈0,则有:The actual rotation angle of the beam of the gantry platform of the placement machine is α. Considering the actual situation that α≈0, we have:

Figure BDA0004061058310000071
Figure BDA0004061058310000071

定义整个贴片机龙门系统状态矢量矩阵q=[x(t),y(t),α(t)]。x(t)和y(t)分别为t时刻贴片机龙门平台X轴和Y轴的实际位移,α(t)为t时刻贴片机龙门平台系统横梁的实际旋转角度,可由Y1轴和Y2轴上光栅编码器反馈获得。考虑X、Y1和Y2轴的运动学关系,并结合Lagrange(拉格朗日)方程,进行动力学建模为:Define the state vector matrix q = [x(t), y(t), α(t)] of the entire SMT gantry system. x(t) and y(t) are the actual displacements of the X-axis and Y-axis of the SMT gantry platform at time t, respectively. α(t) is the actual rotation angle of the beam of the SMT gantry platform system at time t, which can be obtained by the feedback of the grating encoders on the Y1 and Y2 axes. Considering the kinematic relationship between the X, Y1 and Y2 axes, and combining the Lagrange equation, the dynamic modeling is as follows:

Figure BDA0004061058310000072
Figure BDA0004061058310000072

其中,

Figure BDA0004061058310000073
Figure BDA0004061058310000074
分别为贴片机龙门平台系统轴运动的速度和加速度矢量矩阵,
Figure BDA0004061058310000075
为非线性函数通常由
Figure BDA0004061058310000076
来替代。v=Tqur表示所设计的控制输入矩阵,r=x,1,2,ux、u1、u2作为三个电机输入电压是本实施方式最终要调制输出的参数。
Figure BDA0004061058310000077
Figure BDA0004061058310000078
分别为归一化后非线性扰动的常值矩阵和时变矩阵,且
Figure BDA0004061058310000079
in,
Figure BDA0004061058310000073
and
Figure BDA0004061058310000074
They are the velocity and acceleration vector matrices of the axis motion of the SMT gantry platform system,
Figure BDA0004061058310000075
As a nonlinear function, it is usually
Figure BDA0004061058310000076
v = T q u r represents the designed control input matrix, r = x, 1, 2, u x , u 1 , u 2 as the three motor input voltages are the parameters to be modulated and output in the present embodiment.
Figure BDA0004061058310000077
and
Figure BDA0004061058310000078
are the constant matrix and time-varying matrix of the normalized nonlinear perturbation, respectively, and
Figure BDA0004061058310000079

归一化后的惯性矩阵Mq、科氏力矩阵Cq、粘滞摩擦矩阵Bq、库仑摩擦矩阵Aq、等效刚度矩阵Kq和推力系数矩阵Tq分别为:The normalized inertia matrix M q , Coriolis force matrix C q , viscous friction matrix B q , Coulomb friction matrix A q , equivalent stiffness matrix K q and thrust coefficient matrix T q are:

Figure BDA00040610583100000710
Figure BDA00040610583100000710

Figure BDA00040610583100000711
Figure BDA00040610583100000711

Figure BDA00040610583100000712
Figure BDA00040610583100000712

Figure BDA00040610583100000713
Figure BDA00040610583100000713

定义参数“F”的归一化处理公式表示为Fk=F/k1,则其中,Mhk为归一化后贴片机中移动工作台的质量,Bxk为归一化后贴片机龙门平台X轴处的粘滞摩擦系数,Axk为归一化后贴片机龙门平台X轴处的库仑摩擦系数,Myk为归一化后贴片机龙门平台中横梁的总质量,Byk为归一化后贴片机龙门平台Y1轴和Y2轴导轨处的粘滞摩擦系数之和,Ayk为归一化后贴片机龙门平台Y1轴和Y2轴导轨处的库仑摩擦系数之和,Jk为归一化后的转动惯量,Kαk为归一化后的等效旋转刚度,Bsk为归一化后贴片机龙门平台Y2轴和Y1轴导轨处的粘滞摩擦系数之差,Ask为归一化后贴片机龙门平台Y2轴和Y1轴导轨处的库仑摩擦系数之差,The normalized processing formula for defining the parameter "F" is expressed as F k =F/k 1 , where M hk is the mass of the moving workbench in the SMT machine after normalization, B xk is the viscous friction coefficient at the X-axis of the SMT machine gantry platform after normalization, A xk is the Coulomb friction coefficient at the X-axis of the SMT machine gantry platform after normalization, Myk is the total mass of the beam in the SMT machine gantry platform after normalization, Byk is the sum of the viscous friction coefficients at the Y1-axis and Y2-axis guide rails of the SMT machine gantry platform after normalization, A yk is the sum of the Coulomb friction coefficients at the Y1-axis and Y2-axis guide rails of the SMT machine gantry platform after normalization, J k is the normalized moment of inertia, K αk is the normalized equivalent rotational stiffness, B sk is the difference in the viscous friction coefficients at the Y2-axis and Y1-axis guide rails of the SMT machine gantry platform after normalization, and A sk is the normalized difference in Coulomb friction coefficient between the Y2-axis and Y1-axis guide rails of the SMT gantry platform.

Kmy=K2/K1K my =K 2 /K 1 ,

K1和K2分别代表Y1轴和Y2轴上电机的推力常数,Kx为X轴上电机的推力常数,可通过查找手册获得。另外,By=B1+B2,Bs=B2-B1,Ay=A1+A2,As=A2-A1,Bθ和Aθ,θ=1,2,分别代表Y1和Y2轴导轨处的粘滞摩擦系数和库仑摩擦系数。x为移动平台的位置反馈,由X轴光栅编码器获得。h表示移动平台质心偏移距离。 K1 and K2 represent the thrust constants of the motors on the Y1 and Y2 axes, respectively. Kx represents the thrust constant of the motor on the X axis, which can be obtained by looking up the manual. In addition, By = B1 + B2 , Bs = B2 - B1 , Ay = A1 + A2 , As = A2 - A1 , Bθ and , θ = 1, 2, represent the viscous friction coefficient and Coulomb friction coefficient at the Y1 and Y2 axis guide rails, respectively. x is the position feedback of the mobile platform, obtained by the X-axis grating encoder. h represents the displacement distance of the center of mass of the mobile platform.

步骤二、依据期望进行性能函数设计,并对复杂约束进行误差变换。Step 2: Design the performance function based on expectations and perform error transformation on complex constraints.

设定期望轨迹信号矩阵qd为:Set the expected trajectory signal matrix qd to:

qd=[xd(t),yd(t),αd(t)],q d =[x d (t), y d (t), α d (t)],

式中,期望横梁旋转角度αd(t)=0。Wherein, the desired beam rotation angle α d (t)=0.

定义系统误差为:Define the systematic error as:

ei(t)=q-qd=[x(t)-xd(t),y(t)-yd(t),α(t)],e i (t) = qq d = [x (t)-x d (t), y (t) - y d (t), α (t)],

式中,x(t)和y(t)分别为t时刻贴片机龙门平台X轴和Y轴的实际位移,α(t)为t时刻贴片机龙门平台系统横梁的实际旋转角度,这些参数均通过光栅反馈后计算获得。Where x(t) and y(t) are the actual displacements of the X-axis and Y-axis of the SMT gantry platform at time t, respectively, and α(t) is the actual rotation angle of the beam of the SMT gantry platform system at time t. These parameters are calculated after grating feedback.

引入预设性能函数:Introduce preset performance functions:

Figure BDA0004061058310000081
Figure BDA0004061058310000081

其中,ρi0为稳态误差初始值,ρi∞为稳态误差最大值,

Figure BDA0004061058310000082
为误差收敛速率,均是待设计的常值。Among them, ρ i0 is the initial value of the steady-state error, ρ i∞ is the maximum value of the steady-state error,
Figure BDA0004061058310000082
is the error convergence rate, all of which are constants to be designed.

设定系统误差满足约束:Set the system error to satisfy the constraints:

Figure BDA0004061058310000091
Figure BDA0004061058310000091

式中,ei(0)为误差矢量初始值,

Figure BDA0004061058310000092
为超调量调节系数并决定着系统的超调量包络。Where e i (0) is the initial value of the error vector,
Figure BDA0004061058310000092
It is the overshoot adjustment coefficient and determines the overshoot envelope of the system.

引入误差转换函数简化控制器设计:Introducing the error transfer function simplifies the controller design:

Figure BDA00040610583100000915
Figure BDA00040610583100000915

优选的,双曲正切函数为误差转换函数Si,并进行求其反函数为:Preferably, the hyperbolic tangent function is the error conversion function S i , and its inverse function is obtained as follows:

Figure BDA0004061058310000093
Figure BDA0004061058310000093

式中,ri=ei(t)/ρi(t)。In the formula, r i =e i (t)/ρ i (t).

对上式进行微分有:Differentiating the above formula gives:

Figure BDA0004061058310000094
Figure BDA0004061058310000094

其中,

Figure BDA0004061058310000095
in,
Figure BDA0004061058310000095

此时带约束的系统误差ei(t)经过变换,成为无约束的转换误差εi(t),以供后面控制器设计需要。At this time, the constrained system error e i (t) is transformed into an unconstrained conversion error ε i (t) for the subsequent controller design needs.

步骤三、结合给定的期望信号,对动力学模型进行自适应控制器设计。Step 3: Design an adaptive controller for the dynamic model based on the given expected signal.

对步骤一系统动力学模型中不确定矢量进行参数化处理,所构成的系数矩阵

Figure BDA0004061058310000096
为:The coefficient matrix formed by parameterizing the uncertainty vector in the system dynamics model in step 1 is
Figure BDA0004061058310000096
for:

Figure BDA0004061058310000097
Figure BDA0004061058310000097

对于实际的物理系统,参数矢量和集总不确定性满足以下界限:For real physical systems, the parameter vector and lumped uncertainty satisfy the following bounds:

Figure BDA0004061058310000098
Figure BDA0004061058310000098

Figure BDA0004061058310000099
Figure BDA0004061058310000099

式中,

Figure BDA00040610583100000910
Figure BDA00040610583100000911
分别为
Figure BDA00040610583100000912
的上下边界,D为集总不确定性的上界,令
Figure BDA00040610583100000913
为参数自适应估计值,
Figure BDA00040610583100000914
为参数估计误差,则:In the formula,
Figure BDA00040610583100000910
and
Figure BDA00040610583100000911
They are
Figure BDA00040610583100000912
The upper and lower boundaries of , D is the upper bound of the lumped uncertainty, let
Figure BDA00040610583100000913
is the adaptive estimate of the parameter,
Figure BDA00040610583100000914
is the parameter estimation error, then:

Figure BDA0004061058310000101
Figure BDA0004061058310000101

优选设计有界不连续投影自适应律:Optimally design bounded discontinuous projection adaptive law:

Figure BDA0004061058310000102
Figure BDA0004061058310000102

其中,Γ为待设计的自适应学习率矩阵,τ是自适应函数,

Figure BDA0004061058310000103
为常见的标准投影自适应率。
Figure BDA0004061058310000104
是参数
Figure BDA0004061058310000105
的估计值,
Figure BDA0004061058310000106
为设计的自适应率上界,
Figure BDA0004061058310000107
Figure BDA0004061058310000108
的一阶导数。Among them, Γ is the adaptive learning rate matrix to be designed, τ is the adaptive function,
Figure BDA0004061058310000103
is the common standard projection adaptation rate.
Figure BDA0004061058310000104
is a parameter
Figure BDA0004061058310000105
The estimated value of
Figure BDA0004061058310000106
is the upper bound of the designed adaptation rate,
Figure BDA0004061058310000107
for
Figure BDA0004061058310000108
The first derivative of .

设计带有遗忘因子和协方差重置参数估计算法为:Design parameter estimation algorithm with forgetting factor and covariance reset as:

Figure BDA0004061058310000109
Figure BDA0004061058310000109

其中,

Figure BDA00040610583100001010
是滤波后的线性回归矩阵,
Figure BDA00040610583100001011
表示估计误差;μ代表归一化因子且δ为遗忘因子;
Figure BDA00040610583100001012
是设计的自适应更新率上界,tr(·)则表示求迹运算。in,
Figure BDA00040610583100001010
is the filtered linear regression matrix,
Figure BDA00040610583100001011
represents the estimation error; μ represents the normalization factor and δ is the forgetting factor;
Figure BDA00040610583100001012
is the upper bound of the designed adaptive update rate, and tr(·) represents the trace operation.

定义贴片机龙门平台系统的回归量矩阵为:The regression matrix of the placement machine gantry platform system is defined as:

Figure BDA00040610583100001013
Figure BDA00040610583100001013

其中,

Figure BDA00040610583100001014
为回归量残差矩阵,
Figure BDA00040610583100001015
Figure BDA00040610583100001016
分别为期望速度和期望加速度;in,
Figure BDA00040610583100001014
is the residual matrix of the regression variables,
Figure BDA00040610583100001015
and
Figure BDA00040610583100001016
are the expected speed and expected acceleration respectively;

Figure BDA00040610583100001017
是仅含期望信息的回归量矩阵。
Figure BDA00040610583100001017
is the regressor matrix containing only the expected information.

设计期望模型补偿控制器为:Design the desired model compensation controller as:

Figure BDA00040610583100001018
Figure BDA00040610583100001018

实现对动力学模型参数补偿。Realize compensation of dynamic model parameters.

本步骤是对步骤一中系统动力学模型的13个参数进行在线估计,实现对模型的补偿。This step is to perform online estimation of the 13 parameters of the system dynamics model in step 1 to achieve compensation for the model.

步骤四、通过获取到的位置/速度信息,设计鲁棒反馈控制器;Step 4: Design a robust feedback controller based on the acquired position/speed information;

定义滑模函数:Define the sliding mode function:

Figure BDA00040610583100001019
Figure BDA00040610583100001019

其中,λ=diag[λ123]为控制增益矩阵,λi为正的常值;系统的转换误差ε随着类滑模变量s减小而减小。Wherein, λ=diag[λ 123 ] is the control gain matrix, λ i is a positive constant; the conversion error ε of the system decreases as the quasi-sliding mode variable s decreases.

定义半正定能量函数:Define the semi-positive definite energy function:

Figure BDA0004061058310000111
Figure BDA0004061058310000111

其中,Z=diag[z1,z2,z3],Ke为鲁棒增益矩阵。Wherein, Z = diag[z 1 ,z 2 ,z 3 ], Ke is the robust gain matrix.

鲁棒反馈控制器υe为:The robust feedback controller υ e is:

υe=-MqZ-1G-Kds-Keε,υ e =-M q Z -1 GK d sK e ε,

其中,

Figure BDA0004061058310000112
Kd为比例反馈矩阵,Ke为鲁棒增益矩阵,且同时满足:in,
Figure BDA0004061058310000112
Kd is the proportional feedback matrix, Ke is the robust gain matrix, and they both satisfy:

Figure BDA0004061058310000113
Figure BDA0004061058310000113

Figure BDA0004061058310000114
Figure BDA0004061058310000114

其中,σmin(·)和σmax(·)为求矩阵·的最小与最大特征根运算。Among them, σ min (·) and σ max (·) are operations for finding the minimum and maximum eigenvalues of the matrix ·.

考虑到

Figure BDA0004061058310000115
Figure BDA0004061058310000116
有界,得到以下合理假设:Considering
Figure BDA0004061058310000115
and
Figure BDA0004061058310000116
is bounded, and the following reasonable assumptions are obtained:

Figure BDA0004061058310000117
Figure BDA0004061058310000117

其中,ξ为一任意小正数,

Figure BDA0004061058310000118
为待设计的理想神经网络最优补偿。Among them, ξ is an arbitrarily small positive number,
Figure BDA0004061058310000118
The optimal compensation for the ideal neural network to be designed.

对V(t)求导并简化可得:Taking the derivative of V(t) and simplifying it, we get:

Figure BDA0004061058310000119
Figure BDA0004061058310000119

其中,κ=[||ε||,||s||]T,系数β=2σmin(H)/{max[σmax(MqZ-1),σ(Ke)]}。由李雅普诺夫定理可知,所有信号都是有界的。Among them, κ=[||ε||,||s||] T , coefficient β=2σ min (H)/{max[σ max (M q Z -1 ),σ(K e )]}. From Lyapunov's theorem, we know that all signals are bounded.

本步骤是设计鲁棒反馈控制器的过程,鲁棒反馈控制器是基于实际误差信号所设计的控制器,主要用于镇定名义系统。This step is the process of designing a robust feedback controller. The robust feedback controller is a controller designed based on the actual error signal and is mainly used to stabilize the nominal system.

步骤五、对未建模动态和补偿残差进行神经网络逼近。Step 5: Perform neural network approximation on the unmodeled dynamics and compensation residuals.

RBF(Radial Basis Function,径向基函数)神经网络包含输入层,隐含层和输出层,其中隐含层的激活函数使用高斯基函数,第j层神经元的输出表示为:The RBF (Radial Basis Function) neural network consists of an input layer, a hidden layer, and an output layer. The activation function of the hidden layer uses the Gaussian basis function, and the output of the j-th layer neuron is expressed as:

Figure BDA00040610583100001110
Figure BDA00040610583100001110

其中,χ为神经网络输入矢量,χ=[χ12,...,χm],

Figure BDA0004061058310000121
m为神经网络的神经元总层数,cj和bj分别为核函数的中心坐标矢量与宽度。Where χ is the neural network input vector, χ=[χ 12 ,...,χ m ],
Figure BDA0004061058310000121
m is the total number of neuron layers in the neural network, c j and b j are the center coordinate vector and width of the kernel function respectively.

使用RBF神经网络逼近残差及未知误差,其理想神经网络的总输出为:Using RBF neural network to approximate residual and unknown errors, the total output of the ideal neural network is:

F=W*TR(χ)+∈,F=W * TR(χ)+∈,

其中,W*为最优权重矢量,∈表示逼近输出误差,R(χ)为输入矢量χ对应的多层神经元输出。Among them, W * is the optimal weight vector, ∈ represents the approximate output error, and R(χ) is the multi-layer neuron output corresponding to the input vector χ.

实际设计过程中,RBF神经网络输出为:In the actual design process, the output of the RBF neural network is:

Figure BDA0004061058310000122
Figure BDA0004061058310000122

定义

Figure BDA0004061058310000123
为权重误差矢量,其中
Figure BDA0004061058310000124
为所设计的权重矢量。definition
Figure BDA0004061058310000123
is the weight error vector, where
Figure BDA0004061058310000124
is the designed weight vector.

在实际系统中,合理假设逼近误差∈是有上界的,即:In practical systems, it is reasonable to assume that the approximation error ∈ is bounded, namely:

∈≤|∈M|,∈≤|∈ M |,

其中,∈M为待设计的逼近误差上界值。Among them, ∈ M is the upper bound of the approximation error to be designed.

优选设计神经网络控制器为:The optimal design of the neural network controller is:

Figure BDA0004061058310000125
Figure BDA0004061058310000125

其中,η为鲁棒项系数,sign(·)表示符号函数系数η≥|∈M+ξ|。Where η is the robust term coefficient and sign(·) represents the sign function coefficient η ≥ |∈ M +ξ|.

进一步设计RBF神经网络的权重自适应律:Further design the weight adaptation law of RBF neural network:

Figure BDA0004061058310000126
Figure BDA0004061058310000126

其中,φ为待设计的学习率更新矩阵。Among them, φ is the learning rate update matrix to be designed.

重新定义李雅普诺夫能量函数Vq(t),即:Redefine the Lyapunov energy function V q (t), that is:

Figure BDA0004061058310000127
Figure BDA0004061058310000127

其中,tr(·)表示为对矩阵·进行求迹运算。Here, tr(·) represents the trace operation of the matrix ·.

对式

Figure BDA0004061058310000128
进行求导并简化可得:Pair
Figure BDA0004061058310000128
Taking the derivative and simplifying we get:

Figure BDA0004061058310000129
Figure BDA0004061058310000129

因此闭环系统渐进稳定,即贴片机龙门平台运动轨迹会渐进收敛到期望参考信号。Therefore, the closed-loop system is asymptotically stable, that is, the motion trajectory of the gantry platform of the placement machine will gradually converge to the expected reference signal.

本步骤对未建模动态和补偿残差进行神经网络逼近,进一步提高系统的稳定性。This step performs neural network approximation on the unmodeled dynamics and compensation residuals to further improve the stability of the system.

步骤六、优化并综合所设计的控制器,对系统进行实时控制。Step 6: Optimize and integrate the designed controller to control the system in real time.

综合步骤三至五中设计的控制器单元,所述贴片机龙门平台的协同控制器可表示为:Combining the controller units designed in steps 3 to 5, the collaborative controller of the placement machine gantry platform can be expressed as:

v=υaennv=υ aenn .

根据式上式所述总控制器输出矩阵v及推力系数,获得直线电机控制输入为:According to the total controller output matrix v and thrust coefficient described in the above formula, the linear motor control input is obtained as:

Figure BDA0004061058310000131
Figure BDA0004061058310000131

Sm为贴片机龙门平台横梁两侧电机之间的距离。Kxk与Kmy可通过离线辨识获得。将设计的控制系统应用到具体的双驱龙门实验平台上,设计控制器参数以及优化调节各个参数自适应律系数,使得龙门平台的多轴协同跟踪误差与同步控制精度满足预设性能要求。 Sm is the distance between the motors on both sides of the beam of the gantry platform of the placement machine. Kxk and Kmy can be obtained through offline identification. The designed control system is applied to the specific dual-drive gantry experimental platform, the controller parameters are designed and the adaptive law coefficients of each parameter are optimized and adjusted, so that the multi-axis collaborative tracking error and synchronous control accuracy of the gantry platform meet the preset performance requirements.

下面给出一个实施例:An embodiment is given below:

设定Sm=0.975m,Sg=1.095m,Se=1.185m,期望输入xd(t)=0.08*sin(πt),yd(t)=0.04*[1-cos(πt)],αd(t)=0。使用最小二乘法进行参数辨识,得到名义值及参数上下限设定分别为:Set S m = 0.975m, S g = 1.095m, Se = 1.185m, and expect input x d (t) = 0.08*sin(πt), y d (t) = 0.04*[1-cos(πt)], α d (t) = 0. Use the least squares method to identify the parameters, and get the nominal value and the upper and lower limits of the parameters as follows:

Figure BDA0004061058310000132
Figure BDA0004061058310000132

Figure BDA0004061058310000133
Figure BDA0004061058310000133

Figure BDA0004061058310000134
Figure BDA0004061058310000134

辨识得到Kxk≈1.08,Kmy≈1.02,选择光滑函数

Figure BDA0004061058310000135
设计控制器增益矩阵λ=diag[120,100,80]。The identification results are K xk ≈1.08, K my ≈1.02, and the smooth function is selected
Figure BDA0004061058310000135
Design the controller gain matrix λ = diag[120,100,80].

预设性能函数参数ρ0=[0.0003,0.0005,0.0002],ρ=[0.0001,0.0002,0.0001],

Figure BDA0004061058310000138
The preset performance function parameters are ρ 0 = [0.0003, 0.0005, 0.0002], ρ = [0.0001, 0.0002, 0.0001],
Figure BDA0004061058310000138

设定

Figure BDA0004061058310000136
自适应部分参数自适应学习率为Γ(0)=[100,100,...100,100]13×13。set up
Figure BDA0004061058310000136
The adaptive learning rate of the adaptive part parameters is Γ(0)=[100,100,...100,100] 13×13 .

自适应率上界

Figure BDA0004061058310000137
遗忘因子和标准化因子分别为δ=0.02,μ=0.1。Adaptation rate upper bound
Figure BDA0004061058310000137
The forgetting factor and normalization factor are δ=0.02 and μ=0.1 respectively.

设计鲁棒反馈部分参数为:Kd=diag[30,7,2],Ke=diag[100,200,300]。The parameters of the robust feedback part are designed as: K d = diag[30,7,2], Ke = diag[100,200,300].

神经网络部分高斯基函数的中心值与宽度分别为:The center value and width of the Gaussian basis function of the neural network are:

cj=[-0.3,-0.2,-0.1,0,0.1,0.2,0,3]2×7,bj=[10,10,10,10,10,10,10]Tc j =[-0.3,-0.2,-0.1,0,0.1,0.2,0,3] 2×7 , b j =[10,10,10,10,10,10,10] T .

学习率设定为φ=diag[5,5,5,5,5,5,5],鲁棒增益为η=0.05。The learning rate is set to φ = diag[5,5,5,5,5,5,5] and the robust gain is η = 0.05.

可选择的,将移动平台添加5kg额外负载,且初始启动位置位于广义坐标零点。Optionally, an additional load of 5 kg is added to the mobile platform, and the initial starting position is located at the generalized coordinate zero point.

图3所示XY轴跟随期望轨迹时的运动位置变化曲线。FIG3 shows the motion position change curve of the XY axis when following the desired trajectory.

图4和5所示分别为X和Y轴跟随期望轨迹时的跟踪误差变化曲线。与传统非线性控制算法(PID)相比,本实施例所提的自适应预设性能协同控制方法(APPC)在跟踪期望轨迹方面性能更优,所产生的误差更小,即控制精度更高。Figures 4 and 5 show the tracking error curves of the X and Y axes when they follow the desired trajectory. Compared with the traditional nonlinear control algorithm (PID), the adaptive preset performance cooperative control method (APPC) proposed in this embodiment has better performance in tracking the desired trajectory and generates smaller errors, that is, higher control accuracy.

图6所示为横梁在跟随期望轨迹时的旋转角度变化曲线。从实验结果可以看出,本实施例的同步控制方法抗干扰能力更强,横梁旋转角度波动更小,即使在有动态负载和横梁质心偏移影响下仍能保持较高的同步控制精度,确保贴片机龙门系统平稳安全运行。Figure 6 shows the curve of the rotation angle change of the beam when following the desired trajectory. From the experimental results, it can be seen that the synchronous control method of this embodiment has stronger anti-interference ability, smaller fluctuation of the beam rotation angle, and can maintain a high synchronization control accuracy even under the influence of dynamic load and beam center of mass offset, ensuring the smooth and safe operation of the placement machine gantry system.

虽然在本文中参照了特定的实施方式来描述本发明,但是应该理解的是,这些实施例仅仅是本发明的原理和应用的示例。因此应该理解的是,可以对示例性的实施例进行许多修改,并且可以设计出其他的布置,只要不偏离所附权利要求所限定的本发明的精神和范围。应该理解的是,可以通过不同于原始权利要求所描述的方式来结合不同的从属权利要求和本文中所述的特征。还可以理解的是,结合单独实施例所描述的特征可以使用在其它所述实施例中。Although the present invention is described herein with reference to specific embodiments, it should be understood that these embodiments are merely examples of the principles and applications of the present invention. It should therefore be understood that many modifications may be made to the exemplary embodiments and that other arrangements may be devised without departing from the spirit and scope of the present invention as defined by the appended claims. It should be understood that the various dependent claims and features described herein may be combined in a manner different from that described in the original claims. It will also be understood that the features described in conjunction with a single embodiment may be used in other described embodiments.

Claims (9)

1. The multi-axis cooperative control method of the dual-drive gantry platform of the chip mounter is characterized in that a desired model compensation controller, a robust feedback controller and a neural network controller are respectively designed, the sum output by the designed controllers is used as a total control signal of the gantry platform of the chip mounter, and the total control signal is utilized to carry out multi-axis cooperative control on an X axis, a Y1 axis and a Y2 axis of the gantry platform of the chip mounter;
the desired model compensation controller expression is:
Figure FDA0004061058300000011
wherein ,υa Compensating the output of the controller for the desired model, Φ d For a regression quantity matrix containing the desired information,
Figure FDA0004061058300000012
is an estimated value of theta and,θ is a coefficient matrix formed by uncertain parameters in a gantry platform system coupling dynamics model of the chip mounter,
Figure FDA0004061058300000013
wherein ,Mhk For normalizing the quality of the movable workbench in the chip mounter xk For normalizing the viscous friction coefficient of the X axis of the gantry platform of the post-chip mounter, A xk For normalizing the Coulomb friction coefficient, M, at the X axis of the gantry platform of the chip mounter yk For normalizing the total mass of the cross beam in the gantry platform of the rear chip mounter, B yk For normalizing the sum of viscous friction coefficients at Y1 axis and Y2 axis guide rails of the gantry platform of the post-chip mounter, A yk For normalizing the sum of Coulomb friction coefficients at Y1 axis and Y2 axis guide rails of a gantry platform of the post-chip mounter, J k To normalize the moment of inertia, K αk For normalized equivalent rotational stiffness, B sk For normalizing the difference between the viscous friction coefficients at the Y2 axis and Y1 axis guide rails of the gantry platform of the post-chip mounter, A sk For normalizing the difference of coulomb friction coefficients at the Y2 axis and Y1 axis guide rails of the gantry platform of the chip mounter,
Figure FDA0004061058300000014
and
Figure FDA0004061058300000015
Respectively taking the non-linear disturbance constant values in X-axis and Y-axis dynamics of a gantry platform of the normalized chip mounter as +.>
Figure FDA0004061058300000016
The method is a constant value of nonlinear disturbance in the rotation dynamics of a gantry platform of the normalized chip mounter;
the robust feedback controller expression is:
υ e =-M q Z -1 G-K d s-K e ε,
wherein ,υe For robust feedback controllersOutput, M q For the inertia matrix of the gantry platform of the normalized chip mounter, s is a sliding mode variable, and Z=diag [ Z ] 1 ,z 2 ,z 3 ],
Figure FDA0004061058300000017
i=1,2,3,
S i Is an error transfer function and has:
Figure FDA0004061058300000021
ε=[ε 123 ]as unconstrained conversion errors, the intermediate variable r i =e i (t)ρ i (t),e i (t) is an error vector at time t, e i (0) As an initial value of the error vector,
Figure FDA0004061058300000027
for the overshoot adjustment factor,
ρ i is a performance function and expressed as follows:
ρ i (t)=(ρ i0i∞ )e -lti∞
ρ i0 for the steady state error initial value ρ i∞ Is the steady state error maximum, l is the error convergence rate,
Figure FDA0004061058300000022
lambda is the control gain matrix, K d For proportional feedback matrix, K e Is a robust gain matrix;
the neural network controller expression is:
Figure FDA0004061058300000023
wherein ,υnn As an output of the neural network controller,
Figure FDA0004061058300000024
is a neural network weight vector, χ is a neural network input vector, χ= [ χ ] 12 ,...,χ m ],
Figure FDA0004061058300000025
m is the total number of layers of neurons of the neural network, j=1, 2, m, R (χ) is the total output of the multi-layer neurons of the neural network, η is the robust term coefficient, sign (·) represents a sign function;
total control signal v= [ v ] of gantry platform of chip mounter 1 ,v 2 ,v 3 ]:
v=υ aenn
2. The multi-axis cooperative control method of a dual-drive gantry platform of a chip mounter according to claim 1, wherein cooperative control of an X axis, a Y1 axis and a Y2 axis of the gantry platform of the chip mounter is realized by using a total control signal of the gantry platform of the chip mounter according to the following formula:
Figure FDA0004061058300000026
wherein ,Sm For the distance between motors at two sides of a gantry platform beam of the chip mounter, alpha is the actual rotation angle of the gantry platform beam of the chip mounter, K xk K is the normalized thrust constant of the X-axis motor my =K 2 K 1 ,K 1 and K2 Thrust constants, u, of motors of Y1 axis and Y2 axis of gantry platform of chip mounter respectively x 、u 1 and u2 Control inputs for the X-axis, Y1-axis and Y2-axis, respectively.
3. The multi-axis cooperative control method of the dual-drive gantry platform of the chip mounter according to claim 2, wherein the coupling dynamics model of the gantry platform system of the chip mounter is as follows:
Figure FDA0004061058300000031
wherein ,Cq 、B q 、K q and Aq Respectively normalized Coriolis force matrix, viscous friction matrix, equivalent stiffness matrix and coulomb friction matrix, wherein q is a gantry platform system state vector matrix of the chip mounter and q= [ x (t), y (t), and alpha (t)]X (t) and Y (t) are respectively the actual displacement of the X axis and the Y axis of the gantry platform of the chip mounter at the moment t, alpha (t) is the actual rotation angle of the cross beam of the gantry platform system of the chip mounter at the moment t,
Figure FDA0004061058300000032
and
Figure FDA0004061058300000033
The motion speed and acceleration vector matrix of the gantry platform system shaft of the chip mounter are respectively a nonlinear function +.>
Figure FDA0004061058300000034
Figure FDA0004061058300000035
and
Figure FDA0004061058300000036
Respectively a constant matrix and a time-varying matrix of the nonlinear disturbance after normalization.
4. The multi-axis cooperative control method of a dual-drive gantry platform of a chip mounter according to claim 3, wherein an error vector e of the gantry platform system of the chip mounter i The constraint of (t) is:
Figure FDA0004061058300000037
wherein ,ei (t)=q-q d =[x(t)-x d (t),y(t)-y d (t),α(t)],
q d =[x d (t),y d (t),α d (t)]Expected track signal matrix, x of gantry platform system of chip mounter d(t) and yd (t) the expected displacement of X axis and Y axis of gantry platform of chip mounter at t moment respectively, alpha d (t) is the expected rotation angle of the gantry platform beam of the chip mounter and alpha d (t)=0。
5. The multi-axis cooperative control method of the dual-drive gantry platform of the chip mounter according to claim 3, wherein,
Figure FDA0004061058300000038
Figure FDA0004061058300000041
Figure FDA0004061058300000042
wherein h is the offset distance of the mass center of the movable workbench in the chip mounter, X (t) is the actual displacement of the X axis of the cross beam of the gantry platform system of the chip mounter, and S g Is the distance between the Y1 axis guide rail and the Y2 axis guide rail of the gantry platform of the chip mounter.
6. The multi-axis cooperative control method of a dual-drive gantry platform of a chip mounter according to claim 1, wherein the regression of the gantry platform system of the chip mounter
Figure FDA0004061058300000043
The method comprises a residual matrix part and a part containing expected information, wherein the expression is as follows:
Figure FDA0004061058300000044
wherein ,
Figure FDA0004061058300000045
for regression matrix containing only the desired information, +.>
Figure FDA0004061058300000046
And the regression quantity residual error matrix.
7. The multi-axis cooperative control method of a dual-drive gantry platform of a chip mounter according to claim 1, wherein the expression of the sliding mode variable s is:
Figure FDA00040610583000000412
wherein ,
Figure FDA0004061058300000047
8. the multi-axis cooperative control method of a dual-drive gantry platform of a chip mounter according to claim 1, wherein the output R (χ j ) The method comprises the following steps:
Figure FDA0004061058300000048
wherein ,cj and bj The percentages are the central coordinate vector and width of the neural network kernel function.
9. The multi-axis cooperative control method of a dual-drive gantry platform of a chip mounter according to claim 8, wherein the neural network weight vector is
Figure FDA0004061058300000049
Adaptive rate of->
Figure FDA00040610583000000410
The method comprises the following steps:
Figure FDA00040610583000000411
wherein the phi-learning rate updates the matrix.
CN202310059796.7A 2023-01-17 2023-01-17 Multi-axis cooperative control method for dual-drive gantry platform of chip mounter Pending CN116088312A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310059796.7A CN116088312A (en) 2023-01-17 2023-01-17 Multi-axis cooperative control method for dual-drive gantry platform of chip mounter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310059796.7A CN116088312A (en) 2023-01-17 2023-01-17 Multi-axis cooperative control method for dual-drive gantry platform of chip mounter

Publications (1)

Publication Number Publication Date
CN116088312A true CN116088312A (en) 2023-05-09

Family

ID=86207941

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310059796.7A Pending CN116088312A (en) 2023-01-17 2023-01-17 Multi-axis cooperative control method for dual-drive gantry platform of chip mounter

Country Status (1)

Country Link
CN (1) CN116088312A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN119002383A (en) * 2024-10-25 2024-11-22 江苏集萃苏科思科技有限公司 Control method and control device of motion platform system and motion platform system
CN119730230A (en) * 2024-12-26 2025-03-28 北京博瑞先进科技有限公司 Gantry system of high-precision chip mounter and multi-axis coordination control method

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6496765B1 (en) * 2000-06-28 2002-12-17 Sandia Corporation Control system and method for payload control in mobile platform cranes

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6496765B1 (en) * 2000-06-28 2002-12-17 Sandia Corporation Control system and method for payload control in mobile platform cranes

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
PENGWEI SHI 等: "RBF Neural Network-Based Adaptive Robust Synchronization Control of Dual Drive Gantry Stage With Rotational Coupling Dynamics", IEEE TRANSACTIONS ON AUTOMATION SCIENCE AND ENGINEERING, vol. 20, no. 2, 30 May 2022 (2022-05-30), pages 1059 - 1068, XP011938284, DOI: 10.1109/TASE.2022.3177540 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN119002383A (en) * 2024-10-25 2024-11-22 江苏集萃苏科思科技有限公司 Control method and control device of motion platform system and motion platform system
CN119730230A (en) * 2024-12-26 2025-03-28 北京博瑞先进科技有限公司 Gantry system of high-precision chip mounter and multi-axis coordination control method

Similar Documents

Publication Publication Date Title
Shi et al. RBF neural network-based adaptive robust synchronization control of dual drive gantry stage with rotational coupling dynamics
CN116088312A (en) Multi-axis cooperative control method for dual-drive gantry platform of chip mounter
Xie et al. An acceleration-level data-driven repetitive motion planning scheme for kinematic control of robots with unknown structure
Li et al. An input dead zones considered adaptive fuzzy control approach for double pendulum cranes with variable rope lengths
Lin et al. Decoupled torque control of series elastic actuator with adaptive robust compensation of time-varying load-side dynamics
CN110007594B (en) Self-adaptive robust sliding mode control method of series-parallel robot for automobile electrophoretic coating and conveying
Cheng et al. Controller synthesis of tracking and synchronization for multiaxis motion system
Su et al. A neural-network-based controller for a single-link flexible manipulator using the inverse dynamics approach
Wen et al. Elman fuzzy adaptive control for obstacle avoidance of mobile robots using hybrid force/position incorporation
CN115157238B (en) A dynamic modeling and trajectory tracking method for multi-degree-of-freedom robots
CN108942924A (en) Model uncertainty mechanical arm motion control method based on multilayer neural network
CN105773623A (en) SCARA robot trajectory tracking control method based on prediction indirect iterative learning
Wang et al. Global iterative sliding mode control of an industrial biaxial gantry system for contouring motion tasks
CN101661294B (en) Self-adaptive fuzzy control method of strip centering system
Liu et al. Robust control of dual-linear-motor-driven gantry stage for coordinated contouring tasks based on feed velocity
CN114397820B (en) Self-adaptive control method of piezoelectric micro-motion platform based on Hopfield neural network estimator
CN109885077A (en) Attitude control method and controller of a quadrotor aircraft
CN110196555A (en) Neural network profile errors Coupling Control Unit and control method based on friciton compensation
CN117681186A (en) A robust control method for high-precision collaborative manipulator based on inequality constraints
Yao et al. A high-order fully actuated system approach to control of overhead cranes
Yang et al. Concurrent learning-based adaptive control of underactuated robotic systems with guaranteed transient performance for both actuated and unactuated motions
Wan et al. Modeling and synchronized control of a dual-drive “checkerboard” gantry with composite adaptive feedforward and RISE feedback
CN113703451B (en) Adaptive fault-tolerant control method for multi-mobile robot formation with preset performance
Li et al. RBF network dynamic sliding mode robust control for overhead cranes with uncertain parameters
Yang et al. Adaptive pid-like fault-tolerant control for uncertain h-type gantry stage with actuator failure

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