CN108166985A - 一种基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法 - Google Patents
一种基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法 Download PDFInfo
- Publication number
- CN108166985A CN108166985A CN201711277987.1A CN201711277987A CN108166985A CN 108166985 A CN108166985 A CN 108166985A CN 201711277987 A CN201711277987 A CN 201711277987A CN 108166985 A CN108166985 A CN 108166985A
- Authority
- CN
- China
- Prior art keywords
- cutter
- cpi
- formula
- drive shaft
- model
- 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 46
- 230000003416 augmentation Effects 0.000 claims abstract description 23
- 230000001133 acceleration Effects 0.000 claims abstract description 10
- 238000013178 mathematical model Methods 0.000 claims abstract description 7
- 244000309464 bull Species 0.000 claims description 37
- 238000013016 damping Methods 0.000 claims description 15
- 230000015572 biosynthetic process Effects 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000006467 substitution reaction Methods 0.000 claims description 6
- 230000005641 tunneling Effects 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 4
- 230000003321 amplification Effects 0.000 claims description 3
- 230000003190 augmentative effect Effects 0.000 claims description 3
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000005070 sampling Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000009412 basement excavation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 235000013399 edible fruits Nutrition 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21D—SHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
- E21D9/00—Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
- E21D9/06—Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining
- E21D9/08—Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining with additional boring or cutting means other than the conventional cutting edge of the shield
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21D—SHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
- E21D9/00—Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
- E21D9/06—Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining
- E21D9/093—Control of the driving shield, e.g. of the hydraulic advancing cylinders
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Environmental & Geological Engineering (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
本发明公开了一种基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法,首先建立包含齿隙非线性的刀盘驱动***数学模型;接着对刀盘驱动***数学模型进行线性化和简化获得状态空间模型;然后通过把转动角加速度引入到状态变量中,建立刀盘驱动***增广模型;并根据刀盘驱动***的动力学方程,设计驱动轴扭矩的估计算法;最后以刀盘转速跟踪误差和驱动轴扭矩分配不均匀程度为性能指标,设计模型预测控制器,实现刀盘驱动***的控制;并设计了状态观测器以实现对转动角加速度的估计。本发明通过对驱动轴扭矩的估计,并结合模型预测控制器,同时实现了刀盘转速控制与驱动轴扭矩均衡控制,并使***具有更好的抗干扰能力和鲁棒性。
Description
技术领域
本发明盾构机刀盘驱动***的控制领域,具体涉及一种基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法。
背景技术
盾构隧道掘进机是一种专门用于开挖地下隧道工程的掘进装备,盾构隧道掘进机技术集成了计算机、新材料、自动化、信息化、***科学、管理科学等多学科技术,具有掘进速度快、施工周期短、安全可靠性高、生态环境影响小等优点。盾构隧道掘进机刀盘驱动***包括刀盘、齿轮传动***、驱动电机及控制***,通常整个刀盘驱动***需要多台电机同时驱动刀盘旋转。
盾构掘进机的掘进过程中,为了保证施工的安全高效进行,需要控制刀盘转速稳定。但是盾构掘进机施工的地质条件具有很强的随机性,刀盘驱动***面临大范围的突变载荷,因此刀盘驱动***的控制器要具有良好的载荷适应性。
由于刀盘驱动***本质是一个多电机冗余驱动***,刀盘负载扭矩通过大齿圈与小齿轮的啮合分配到多根驱动轴上,在刀盘掘进面载荷剧烈波动的情况下,由于各个电机、齿轮传动***的差异,特别是小齿轮与大齿圈间啮合齿隙的差异,引起驱动轴扭矩分配不均衡,部分驱动轴承受的扭矩过大,将会引起断轴事故,影响隧道掘进施工的高效进行。因此需要协调控制各根驱动轴扭矩分配。
目前国内外刀盘驱动***的控制多采用转速或是转矩的PID并行控制、主从控制,以及采用转速差、电机扭矩差等方式的同步补偿控制,还没有基于驱动轴扭矩估计所设计的直接针对驱动轴扭矩均衡分配的控制算法。
发明内容
针对现有技术的不足,本发明提供一种基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法,该方法通过建立刀盘驱动控制***的增广模型,并引入驱动轴扭矩估计方法和状态观测器,使得驱动轴扭矩分配更加均衡,减少使用加速度传感器。具体技术方案如下:
一种基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法,其特征在于,所述的刀盘驱动***包括驱动电机、减速器、小齿轮、大齿圈,所述的控制方法包括以下步骤:
步骤一、建立包含齿隙的非线性的刀盘驱动***数学模型;
步骤二、因刀盘驱动***正常工况下需要保持刀盘的转动方向一致,因此将步骤一的非线性模型转化为线性化的模型,并对其进行简化,得到刀盘驱动***的状态空间模型;
步骤三、把转动角加速度引入到状态变量中,建立刀盘驱动***的增广模型;
步骤四、根据刀盘驱动***的动力学方程,设计驱动轴扭矩的估计算法;
步骤五、以刀盘转速跟踪误差和驱动轴扭矩分配不均匀程度为性能指标,设计模型预测控制器,实现刀盘驱动***的控制;
进一步地,步骤一中所述的包含齿隙非线性的刀盘驱动***数学模型如下;
(1)驱动电机的电磁力矩与输出力矩关系为:
式中,Tei为第i号驱动电机的电磁力矩,Jmi为第i号驱动电机转动惯量,bmi为第i号驱动电机转动阻尼系数,Mpmi为第i号驱动电机受到减速器的负载力矩,θmi为第i号驱动电机的转动角度;
减速器的减速和力矩放大作用表示为:
θmi=impθpi (2)
impMpmi=Mmpi (3)
式中,imp为减速器的传动比,θpi为经第i号减速器减速后的转动角度,也是对应小齿轮的转动角度,Mmpi为第i号减速器的输出力矩,也是对应的小齿轮的输入力矩。
(2)小齿轮的转动平衡方程为:
式中,Jpi为第i号小齿轮的转动惯量,bpi为第i号小齿轮的转动阻尼系数,Mcpi为第i号小齿轮受到大齿圈的阻力矩。
(3)刀盘的转动平衡方程为:
式中,Jc为刀盘以及大齿圈的转动惯量,bc为刀盘以及大齿圈的转动阻尼系数,TL为盾构掘进中刀盘旋转所受到的阻力矩,Mpci为大齿圈受到第i号小齿轮的啮合力矩,θc为刀盘的转动角度,n为刀盘驱动***的驱动电机个数。
根据齿轮传动的性质有:
ipcMcpi=Mpci (6)
式中,ipc为小齿轮到大齿圈的传动比;
(4)小齿轮与大齿圈间的啮合过程可描述为:
xcpi=rbpθpi-rbcθc-ecpi(t) (7)
式中,xcpi为第i个齿轮副沿着啮合线上动态相对位移,rbp为小齿轮的基圆半径,rbc为大齿圈的基圆半径,ecpi(t)为第i号小齿轮与大齿圈啮合时的当量累计啮合误差,kcpi为第i号小齿轮与大齿圈之间的啮合刚度,Ccpi为第i号小齿轮与大齿圈之间的啮合阻尼,bcpi为第i号小齿轮与大齿圈啮合时存在的齿隙。
进一步地,步骤二中刀盘驱动***的状态空间模型的建立具体为:
在盾构机掘进过程中,刀盘驱动***正常工况下需要保持刀盘的转动方向一致,即保持 xcpi>bcpi,或xcpi<bcpi,下述过程为xcpi>bcpi时的刀盘驱动***模型的简化过程,对于 xcpi<bcpi的情况,仅需将bcpi替换为-bcpi进行简化;
(1)当xcpi>bcpi时,公式(8)简化为:
bcpi是一个固定的值,令ei=ecpi+bcpi,则上式可表示为:
(2)定义状态变量和干扰变量则公式(11)可写成:
Mcpi=MXi·X+Mdi·d (12)
其中,
Mdi=[O1×i-rbpkcpiOn-1-CcpiO1×(n-i)]
(3)对公式(1)-(4)合并简化如下:
其中,表示第i号驱动电机以及其所连接的减速器、小齿轮构成的驱动子***的等效转动惯量,表示驱动子***的等效转动阻尼;
(4)将(12)代入(13)得到
将(12)代入(5)得到:
(5)定义控制变量u=[Te1,...Ten],则(14)、(15)可以写成:
其中,
(6)根据公式(16)、(17),同时定义刀盘转速和驱动轴扭矩y2=[Mcp1,...,Mcpn]T,则***的状态空间模型为
y1=C1X (19)
y2=MX·X+Md·d (20)
其中,
A0=[O(n+1)×(n+1)I(n+1)×(n+1)]
C1=[O1×(2n+1) 1]
进一步地,步骤三中把转动角加速度引入到状态变量中,建立刀盘驱动***的增广模型具体如下:
(1)定义构造增广状态变量Xaug=[XT ZT]T;
记则公式(18)可转换为:
其中:
(2)定义
则公式(21)可转换为如下的连续时间增广模型:
进一步地,步骤四中设计驱动轴扭矩的估计算法具体为;
由公式(13)可得到驱动轴扭矩的估计公式:
将公式(23)表示成矩阵的形式,
Mcpi=MXuaugiXaug·+Muaugi·u (24)
其中,MXaugi=[O1×(n+i)-bEiO1×n-JEiO1×(n+1-i)],Muaugi=[O1×(i-1)impO1×(n-i)];
结合y2的定义和公式(24)可得
y2=MXaugXaug·+Muaug·u (25)
其中,
进一步地,步骤五中以刀盘转速跟踪误差和驱动轴扭矩分配不均匀程度为性能指标,设计模型预测控制器,实现刀盘驱动***的控制具体为;
(1)以采样周期Ts对连续时间的增广模型公式(22)进行离散化,可得到离散时间的増广模型
y1(k)=Cd1augXaug(k) (27)
y2(k)=MdXaugXaug(k)·+Mduaug·u(k) (28)
其中
Cd1aug=C1aug,MdXaug=MXaug,Mduaug-Muaug
(2)定义平均驱动轴扭矩值期望刀盘转速则k时刻的优化性能指标为
(3)设置u(k)和Δu(k)添加约束条件,
umin≤u(k)≤umax (30)
Δumin≤Δu(k)≤Δumax (31)
(4)选择控制时域M和预测时域P后,可得到总的性能指标为:
(5)运用优化算法,在满足公式(30)、(31)的约束条件下求解使公式(32)的总的性能指标最小的优化变量{Δuk},得到k时刻最优的控制输入增量序列{Δu(k),Δu(k 1)...Δu(k+M-1)},选取第一项Δu(k)作为k时刻的***输入,即可实现对刀盘驱动***的控制。
进一步地,所述的增广模型中的Xaug通过状态观测器估计得到,具体如下:
(1)定义刀盘驱动***的可测输出:
(2)由所述的连续时间增广模型公式(22)可得到状态观测器的状态空间模型:
其中,为Xaug的估计值,为yθ的估计值,Haug=[I(2n+2)×(2n+2)O(n+1)×(n+1)],进一步得到全阶观测器:
其中,Laug为状态观测器比例增益矩阵,通过配置(Aaug-LaugHaug)的极点的方法获得Laug各个系数。
与现有技术相比,本发明的有益效果如下:
(1)本发明建立了刀盘驱动控制增广模型,引入一个积分环节,对于恒定负载可以消除稳态误差,***更具鲁棒性;
(2)本发明引入驱动轴扭矩估计方法,解决了驱动轴扭矩难以测量的问题,使控制器设计可以更加充分考虑刀盘驱动***的机电耦合特性,从而使驱动轴扭矩分配更加均衡;
(3)本发明引入状态观测器,可实现刀盘驱动***中转动角加速度估计,可减少使用加速度传感器;
(4)本发明提出将模型预测控制方法与驱动轴扭矩估计方法相结合的方法,在考虑约束条件下同时实现刀盘转速控制目标与驱动轴扭矩均衡控制目标。
附图说明
图1为本发明的控制方法的示意图;
图2为使用本发明的控制方法获得的转速跟踪和负载转矩随时间变化的曲线图;
图3为使用本发明的控制方法获得驱动轴的扭矩随时间变化的曲线图;
图4为使用PI方法控制获得的转速跟踪和负载转矩随时间变化的曲线图;
图5为使用PI方法控制获得的驱动轴的扭矩随时间变化的曲线图。
具体实施方式
如图1所示,一种基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法,其特征在于,所述的刀盘驱动***包括驱动电机、减速器、小齿轮、大齿圈,所述的控制方法包括以下步骤:
步骤一、建立包含齿隙的非线性的刀盘驱动***数学模型;
(1)驱动电机的电磁力矩与输出力矩关系为:
式中,Tei为第i号驱动电机的电磁力矩,Jmi为第i号驱动电机转动惯量,bmi为第i号驱动电机转动阻尼系数,Mpmi为第i号驱动电机受到减速器的负载力矩,θmi为第i号驱动电机的转动角度;
减速器的减速和力矩放大作用表示为:
θmi=impθpi (2)
impMpmi=Mmpi (3)
式中,imp为减速器的传动比,θpi为经第i号减速器减速后的转动角度,也是对应小齿轮的转动角度,Mmpi为第i号减速器的输出力矩,也是对应的小齿轮的输入力矩。
(2)小齿轮的转动平衡方程为:
式中,Jpi为第i号小齿轮的转动惯量,bpi为第i号小齿轮的转动阻尼系数,Mcpi为第i号小齿轮受到大齿圈的阻力矩。
(3)刀盘的转动平衡方程为:
式中,Jc为刀盘以及大齿圈的转动惯量,bc为刀盘以及大齿圈的转动阻尼系数,TL为盾构掘进中刀盘旋转所受到的阻力矩,Mpci为大齿圈受到第i号小齿轮的啮合力矩,θc为刀盘的转动角度,n为刀盘驱动***的驱动电机个数。
根据齿轮传动的性质有:
ipcMcpi=Mpci (6)
式中,ipc为小齿轮到大齿圈的传动比;
(4)小齿轮与大齿圈间的啮合过程可描述为:
xcpi=rbpθpi-rbcθc-ecpi(t) (7)
式中,xcpi为第i个齿轮副沿着啮合线上动态相对位移,rbp为小齿轮的基圆半径,rbc为大齿圈的基圆半径,ecpi(t)为第i号小齿轮与大齿圈啮合时的当量累计啮合误差,kcpi为第i号小齿轮与大齿圈之间的啮合刚度,Ccpi为第i号小齿轮与大齿圈之间的啮合阻尼,bcpi为第i号小齿轮与大齿圈啮合时存在的齿隙。
步骤二、因刀盘驱动***正常工况下需要保持刀盘的转动方向一致,因此将步骤一的非线性模型转化为线性化的模型,并对其进行简化,得到刀盘驱动***的状态空间模型;
在盾构机掘进过程中,刀盘驱动***正常工况下需要保持刀盘的转动方向一致,即保持xcpi>bcpi,或xcpi<bcpi,下述过程为xcpi>bcpi时的刀盘驱动***模型的简化过程,对于 xcpi<bcpi的情况,仅需将bcpi替换为-bcpi进行简化;
(1)当xcpi>bcpi时,公式(8)简化为:
bcpi是一个固定的值,令ei=ecpi+bcpi,则上式可表示为:
(2)定义状态变量和干扰变量则公式(11)可写成:
Mcpi=MXi·X+Mdi·d (12)
其中,
Mdi=[O1×i-rbpkcpi On-1-Ccpi O1×(n-i)]
(3)对公式(1)-(4)合并简化如下:
其中,表示第i号驱动电机以及其所连接的减速器、小齿轮构成的驱动子***的等效转动惯量,表示驱动子***的等效转动阻尼;
(4)将(12)代入(13)得到
将(12)代入(5)得到:
(5)定义控制变量u=[Te1,...Ten],则(14)、(15)可以写成:
其中,
(6)根据公式(16)、(17),同时定义刀盘转速和驱动轴扭矩y2=[Mcp1,...,Mcpn]T,则***的状态空间模型为
y1=C1X (19)
y2=MX·X+Md·d (20)
其中
A0=[O(n+1)×(n+1) I(n+1)×(n+1)]
C1=[O1×(2n+1)1]
步骤三、把转动角加速度引入到状态变量中,建立刀盘驱动***的增广模型;
(1)定义构造增广状态变量Xaug=[XT ZT]T;
记则公式(18)可转换为:
其中:
(2)定义
则公式(21)可转换为如下的连续时间增广模型:
步骤四、根据刀盘驱动***的动力学方程,设计驱动轴扭矩的估计算法;
由公式(13)可得到驱动轴扭矩的估计公式:
将公式(23)表示成矩阵的形式,
Mcpi=MXaugiXaug·+Muaugi·u (24)
其中,MXaugi=[O1×(n+i)-bEiO1×n-JEiO1×(n+1-i)],Muaugi=[O1×(i-1)impO1×(n-i)];
结合y2的定义和公式(24)可得
y2=MXaugXaug·+Muaug·u (25)
其中,步骤五、以刀盘转速跟踪误差和驱动轴扭矩分配不均匀程度为性能指标,设计模型预测控制器,实现刀盘驱动***的控制;
(1)以采样周期Ts对连续时间的增广模型公式(22)进行离散化,可得到离散时间的増广模型
y1(k)=Cd1augXaug(k) (27)
y2(k)=MdXaugXaug(k)·+Mduaug·u(k) (28)
其中
Cd1aug=C1aug,MdXaug=MXaug,Mduaug-Muaug
(2)定义平均驱动轴扭矩值期望刀盘转速则k时刻的优化性能指标为
(3)设置u(k)和Δu(k)添加约束条件,
umin≤u(k)≤umax (30)
Δumin≤Δu(k)≤Δumax (31)
(4)选择控制时域M和预测时域P后,可得到总的性能指标为:
(5)运用优化算法,在满足公式(30)、(31)的约束条件下求解使公式(32)的总的性能指标最小的优化变量{Δuk},得到k时刻最优的控制输入增量序列{Δu(k),Δu(k 1)...Δu(k+M-1)},选取第一项Δu(k)作为k时刻的***输入,即可实现对刀盘驱动***的控制。
优选的,所述的增广模型中的Xaug通过状态观测器估计得到,具体如下:
(1)定义刀盘驱动***的可测输出:
(2)由所述的连续时间增广模型公式(22)可得到状态观测器的状态空间模型:
其中,为Xaug的估计值,为yθ的估计值,Haug=[I(2n+2)×(2n+2) O(n+1)×(n+1)],进一步得到全阶观测器:
其中,Laug为状态观测器比例增益矩阵,通过配置(Aaug-LaugHaug)的极点的方法获得Laug各个系数。
为了验证本发明的有效性,在MATLAB R2017a软件中搭建仿真环境,进行仿真验证。
考虑一个由四个驱动电机共同驱动的盾构掘进机,其刀盘驱动***的各项仿真参数如表 1所示。
按照前面所述方法设计状态观测器、增广模型、模型预测控制器,进而得到仿真结果,如图2和3所示。
从图2中可以看出,在负载干扰下,刀盘转速可以维持在期望值上,刀盘驱动***的负载适应性能良好。从图3中可以看出,在所设计的控制器控制下各个驱动轴扭矩很快就收敛到一致,实现驱动轴扭矩的均衡分配。
为便于比较仿真性能,另设计一组基于PI控制器的转速并行控制的仿真作为比较。设计 PI控制器使闭环极点位于-10处,得到刀盘转速仿真结果图4和驱动轴扭矩仿真结果图5。对比图2与图4可以发现,本发明的控制器可以更快使刀盘转速收敛到期望值;对比图3与图5可以发现,本发明的控制器控制下刀盘驱动轴扭矩可以迅速实现均衡,而使用PI控制器则驱动轴扭矩分配会不均衡,无法收敛到同一值。因此,使用本发明的控制器可以较大提升刀盘驱动控制***的性能。
表1刀盘驱动***的各项仿真参数列表
综上所述,仿真结果验证了本发明所提出的控制方法可以在负载干扰、啮合误差、***个组成动力学性质不一致等情况下保证刀盘驱动***满足转速控制和驱动轴扭矩均衡分配的性能要求,实现良好的控制效果。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (7)
1.一种基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法,其特征在于,所述的刀盘驱动***包括驱动电机、减速器、小齿轮、大齿圈,所述的控制方法包括以下步骤:
步骤一、建立包含齿隙的非线性的刀盘驱动***数学模型;
步骤二、因刀盘驱动***正常工况下需要保持刀盘的转动方向一致,因此将步骤一的非线性模型转化为线性化的模型,并对其进行简化,得到刀盘驱动***的状态空间模型;
步骤三、把转动角加速度引入到状态变量中,建立刀盘驱动***的增广模型;
步骤四、根据刀盘驱动***的动力学方程,设计驱动轴扭矩的估计算法;
步骤五、以刀盘转速跟踪误差和驱动轴扭矩分配不均匀程度为性能指标,设计模型预测控制器,实现刀盘驱动***的控制。
2.根据权利要求1所述的基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法,其特征在于,步骤一中所述的包含齿隙非线性的刀盘驱动***数学模型如下;
(1)驱动电机的电磁力矩与输出力矩关系为:
式中,Tei为第i号驱动电机的电磁力矩,Jmi为第i号驱动电机转动惯量,bmi为第i号驱动电机转动阻尼系数,Mpmi为第i号驱动电机受到减速器的负载力矩,θmi为第i号驱动电机的转动角度;
减速器的减速和力矩放大作用表示为:
θmi=impθpi (2)
impMpmi=Mmpi (3)
式中,imp为减速器的传动比,θpi为经第i号减速器减速后的转动角度,也是对应小齿轮的转动角度,Mmpi为第i号减速器的输出力矩,也是对应的小齿轮的输入力矩。
(2)小齿轮的转动平衡方程为:
式中,Jpi为第i号小齿轮的转动惯量,bpi为第i号小齿轮的转动阻尼系数,Mcpi为第i号小齿轮受到大齿圈的阻力矩。
(3)刀盘的转动平衡方程为:
式中,Jc为刀盘以及大齿圈的转动惯量,bc为刀盘以及大齿圈的转动阻尼系数,TL为盾构掘进中刀盘旋转所受到的阻力矩,Mpci为大齿圈受到第i号小齿轮的啮合力矩,θc为刀盘的转动角度,n为刀盘驱动***的驱动电机个数。
根据齿轮传动的性质有:
ipcMcpi=Mpci (6)
式中,ipc为小齿轮到大齿圈的传动比;
(4)小齿轮与大齿圈间的啮合过程可描述为:
xcpi=rbpθpi-rbcθc-ecpi(t) (7)
式中,xcpi为第i个齿轮副沿着啮合线上动态相对位移,rbp为小齿轮的基圆半径,rbc为大齿圈的基圆半径,ecpi(t)为第i号小齿轮与大齿圈啮合时的当量累计啮合误差,kcpi为第i号小齿轮与大齿圈之间的啮合刚度,Ccpi为第i号小齿轮与大齿圈之间的啮合阻尼,bcpi为第i号小齿轮与大齿圈啮合时存在的齿隙。
3.根据权利要求2所述的基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法,其特征在于,步骤二中刀盘驱动***的状态空间模型的建立具体为:
在盾构机掘进过程中,刀盘驱动***正常工况下需要保持刀盘的转动方向一致,即保持xcpi>bcpi,或xcpi<bcpi,下述过程为xcpi>bcpi时的刀盘驱动***模型的简化过程,对于xcpi<bcpi的情况,仅需将bcpi替换为-bcpi进行简化;
(1)当xcpi>bcpi时,公式(8)简化为:
bcpi是一个固定的值,令ei=ecpi+bcpi,则上式可表示为:
(2)定义状态变量和干扰变量则公式(11)可写成:
Mcpi=MXi·X+Mdi·d (12)
其中,
Mdi=[O1×i-rbpkcpiOn-1-CcpiO1×(n-i)]
(3)对公式(1)-(4)合并简化如下:
其中,表示第i号驱动电机以及其所连接的减速器、小齿轮构成的驱动子***的等效转动惯量,表示驱动子***的等效转动阻尼;
(4)将(12)代入(13)得到
将(12)代入(5)得到:
(5)定义控制变量u=[Te1,...Ten],则(14)、(15)可以写成:
其中,
(6)根据公式(16)、(17),同时定义刀盘转速和驱动轴扭矩y2=[Mcp1,...,Mcpn]T,则***的状态空间模型为
y1=C1X (19)
y2=MX·X+Md·d (20)
其中,
A0=[O(n+1)×(n+1)I(n+1)×(n+1)]
C1=[O1×(2n+1)1]
。
4.根据权利要求3所述的基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法,其特征在于,步骤三中把转动角加速度引入到状态变量中,建立刀盘驱动***的增广模型具体如下:
(1)定义构造增广状态变量Xaug=[XT ZT]T;
记则公式(18)可转换为:
其中:
(2)定义
则公式(21)可转换为如下的连续时间增广模型:
。
5.根据权利要求4所述的一种基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法,其特征在于,步骤四中设计驱动轴扭矩的估计算法具体为;
由公式(13)可得到驱动轴扭矩的估计公式:
将公式(23)表示成矩阵的形式,
Mcpi=MXaugiXaug·+Muaugi·u (24)
其中,MXaugi=[O1×(n+i)-bEi O1×n-JEi O1×(n+1-i)],Muaugi=[O1×(i-1) imp O1×(n-i)];
结合y2的定义和公式(24)可得
y2=MXaugXaug·+Muaug·u (25)
其中,
6.根据权利要求5所述的基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法,其特征在于,步骤五中以刀盘转速跟踪误差和驱动轴扭矩分配不均匀程度为性能指标,设计模型预测控制器,实现刀盘驱动***的控制具体为;
(1)以采样周期Ts对连续时间的增广模型公式(22)进行离散化,可得到离散时间的增广模型
y1(k)=Cd1augXaug(k) (27)
y2(k)=MdXaugXaug(k)·+Mduaug·u(k) (28)
其中
Cd1aug=C1aug,MdXaug=MXaug,Mduaug=Muaug
(2)定义平均驱动轴扭矩值期望刀盘转速则k时刻的优化性能指标为
(3)设置u(k)和Δu(k)添加约束条件,
umin≤u(k)≤umax (30)
Δumin≤Δu(k)≤Δumax (31)
(4)选择控制时域M和预测时域P后,可得到总的性能指标为:
(5)运用优化算法,在满足公式(30)、(31)的约束条件下求解使公式(32)的总的性能指标最小的优化变量{Δuk},得到k时刻最优的控制输入增量序列{Δu(k),Δu(k+1)...Δu(k+M-1)},选取第一项Δu(k)作为k时刻的***输入,即可实现对刀盘驱动***的控制。
7.根据权利要求6所述的控制器设计方法,所述的增广模型中的Xaug通过状态观测器估计得到,具体如下:
(1)定义刀盘驱动***的可测输出:
(2)由所述的连续时间增广模型公式(22)可得到状态观测器的状态空间模型:
其中,为Xaug的估计值,为yθ的估计值,Haug=[I(2n+2)×(2n+2) O(n+1)×(n+1)],进一步得到全阶观测器:
其中,Laug为状态观测器比例增益矩阵,通过配置(Aaug-LaugHaug)的极点的方法获得Laug各个系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711277987.1A CN108166985B (zh) | 2017-12-06 | 2017-12-06 | 一种基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711277987.1A CN108166985B (zh) | 2017-12-06 | 2017-12-06 | 一种基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108166985A true CN108166985A (zh) | 2018-06-15 |
CN108166985B CN108166985B (zh) | 2019-06-11 |
Family
ID=62525392
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711277987.1A Active CN108166985B (zh) | 2017-12-06 | 2017-12-06 | 一种基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108166985B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114706322A (zh) * | 2022-03-30 | 2022-07-05 | 西南交通大学 | 一种盾构机姿态自动控制仿真*** |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6173041A (ja) * | 1984-09-18 | 1986-04-15 | Hitachi Constr Mach Co Ltd | トルク演算装置 |
JPH0860983A (ja) * | 1994-08-25 | 1996-03-05 | Hitachi Constr Mach Co Ltd | トンネル掘削機のカッタトルクの制御装置及び制御方法 |
CN103016014A (zh) * | 2012-07-31 | 2013-04-03 | 中国科学院沈阳自动化研究所 | 盾构掘进机刀盘驱动装置及其同步控制方法 |
CN203702159U (zh) * | 2013-12-27 | 2014-07-09 | 北京市三一重机有限公司 | 盾构机及其刀盘驱动*** |
CN103941586A (zh) * | 2014-02-08 | 2014-07-23 | 上海交通大学 | 一种掘进机主驱动***转速的预测函数控制方法 |
CN107165648A (zh) * | 2017-06-28 | 2017-09-15 | 中船重型装备有限公司 | 用于盾构机的刀盘设计方法 |
CN107246271A (zh) * | 2017-06-28 | 2017-10-13 | 浙江大学 | 一种盾构机刀盘驱动***的驱动轴扭矩协调分配控制方法 |
-
2017
- 2017-12-06 CN CN201711277987.1A patent/CN108166985B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS6173041A (ja) * | 1984-09-18 | 1986-04-15 | Hitachi Constr Mach Co Ltd | トルク演算装置 |
JPH0860983A (ja) * | 1994-08-25 | 1996-03-05 | Hitachi Constr Mach Co Ltd | トンネル掘削機のカッタトルクの制御装置及び制御方法 |
CN103016014A (zh) * | 2012-07-31 | 2013-04-03 | 中国科学院沈阳自动化研究所 | 盾构掘进机刀盘驱动装置及其同步控制方法 |
CN203702159U (zh) * | 2013-12-27 | 2014-07-09 | 北京市三一重机有限公司 | 盾构机及其刀盘驱动*** |
CN103941586A (zh) * | 2014-02-08 | 2014-07-23 | 上海交通大学 | 一种掘进机主驱动***转速的预测函数控制方法 |
CN107165648A (zh) * | 2017-06-28 | 2017-09-15 | 中船重型装备有限公司 | 用于盾构机的刀盘设计方法 |
CN107246271A (zh) * | 2017-06-28 | 2017-10-13 | 浙江大学 | 一种盾构机刀盘驱动***的驱动轴扭矩协调分配控制方法 |
Non-Patent Citations (1)
Title |
---|
刘然: "盾构机刀盘同步驱动***自适应均载控制研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114706322A (zh) * | 2022-03-30 | 2022-07-05 | 西南交通大学 | 一种盾构机姿态自动控制仿真*** |
CN114706322B (zh) * | 2022-03-30 | 2023-06-23 | 西南交通大学 | 一种盾构机姿态自动控制仿真*** |
Also Published As
Publication number | Publication date |
---|---|
CN108166985B (zh) | 2019-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Shi et al. | Energy saving of cutterhead hydraulic drive system of shield tunneling machine | |
Wang et al. | Energy-minimum optimization of the intelligent excavating process for large cable shovel through trajectory planning | |
CN103135512B (zh) | 盾构电液控制***综合试验平台 | |
CN102354146B (zh) | 一种运动控制***及其位置控制方法 | |
CN104454350A (zh) | 风力涡轮机及其降低转子不平衡载荷的控制方法 | |
Huo et al. | Electromechanical coupling dynamics of TBM main drive system | |
Huang et al. | Development of a dual displacement controlled circuit for hydraulic shovel swing motion | |
Kongshu et al. | Force transmission characteristics for the non-equidistant arrangement thrust systems of shield tunneling machines | |
CN103941586A (zh) | 一种掘进机主驱动***转速的预测函数控制方法 | |
CN107246271B (zh) | 一种盾构机刀盘驱动***的驱动轴扭矩协调分配控制方法 | |
Dolipski et al. | Investigating the simulated control of the rotational speed of roadheader cutting heads, relating to the reduction of energy consumption during the cutting process | |
CN108166985A (zh) | 一种基于驱动轴扭矩估计的盾构机刀盘驱动***的控制方法 | |
CN107092189B (zh) | 基于模型预测控制的多变量输入eha***的控制方法 | |
CN101532517A (zh) | 随动***负载的电液伺服模拟方法 | |
CN107608211A (zh) | 一种掘进机多电机刀盘***的分布式预测控制方法 | |
Liu et al. | Data-driven Optimal Control of Earth Pressure Balance for Shield Tunneling Machine. | |
Qi et al. | Research on new intelligent pump control based on sliding mode variable structure control | |
Liu et al. | Review on inlet/outlet oil coordinated control for electro-hydraulic power mechanism under sustained negative load | |
CN108303871B (zh) | 一种掘进机刀盘驱动***转速的分数阶控制方法 | |
CN203950135U (zh) | 一种模拟风力机的实验装置 | |
Chen | [Retracted] Hydraulic Lifting and Rotating System Lifting Machinery Transmission Control Design | |
Yue et al. | Load observer-based integral sliding mode trajectory tracking control of shield automatic tunnelling via hyperbolic tangent function | |
Liu et al. | Trajectory tracking of hard rock tunnel boring machine with cascade control structure | |
Zhang et al. | Dynamic Modeling and CAE Cosimulation Method for Heavy‐Duty Concrete Spreader | |
Kong et al. | Research on dynamic characteristics of implement mechanism of mining excavator |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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 |