Rescheduling Method for Starting Time Delay in Steelmaking and Continuous Casting Production Processes
-
摘要: 在炼钢-连铸生产过程中, 因铁水或废钢供应不及时经常发生钢水在转炉设备上开工延迟, 导致相邻炉次在同一设备上产生作业冲突或同一浇次内的相邻炉次在连铸机上断浇, 使原调度计划失效.由于炼钢-连铸生产存在多台转炉、多台精炼炉和多台连铸机, 以及多重精炼方式且存在着可重入生产情况, 加上要保证已开工和完工炉次按原调度计划执行, 使得已有炼钢-连铸生产重调度方法难以适用.本文建立了开工时间延迟下的炼钢-连铸生产重调度模型, 将模型分解和启发式方法相结合, 提出了由炉次加工设备指派和作业时间决策组成的启发式重调度方法.通过实际工业数据仿真实验和在某大型钢铁企业的炼钢-连铸生产重调度的工业应用表明所提出的重调度方法显著减少了炉次加工冗余等待时间.Abstract: In steelmaking and continuous casting production processes, the starting time delay often occurs, which may lead to casting break or processing conflict so that the initial scheduling plan becomes unrealizable. There are multiple converters, multiple refining machines and multiple continuous casters in the whole production process. It is more difficult to reschedule because of multi-refining number and reentrant line. Existing rescheduling methods are not suitable. A mixed programming model is presented which is decomposed into an equipment assignment submodel and an operation time decision submodel. A novel multi-stage heuristic method is proposed. Simulation experiment is studied for comparing the proposed method with other heuristic methods, which have different ordering rules and assignment rules. The results show that the proposed method is most effective. Finally, the proposed method is successfully applied to a large steel plant to show that the method can obviously shorten the redundant waiting time for all charges.
-
Key words:
- Steelmaking and continuous casting /
- rescheduling /
- heuristic method /
- reentrant
-
滑模变结构控制因其强鲁棒性、模型依赖性低及参数灵活等特点被广泛应用于航天器姿态控制等工程问题中[1-2].该方法存在的较大缺陷是控制过程中的抖振现象[3-4].鉴于此, 国内外学者采用边界层法[5]、与智能控制相结合的方法[6]及基于趋近律的控制方法[7]等来削弱抖振.基于趋近律的滑模变结构控制方法由于可改善系统趋近过程动态品质, 因此自提出以来国内外学者纷纷对其进行研究并取得一定成果. 1996年高为炳[8]提出了经典的指数趋近律, 并通过调整趋近律参数改善了趋近过程的动态品质且一定程度地削弱了抖振, 但仍具备改善空间. Yu等[9]将指数趋近律与单幂次趋近律线性组合, 既缩短了趋近时间又使趋近过程更加平滑, 但到达滑模面后系统仍存在抖振现象. Bandyopadhyay等[10]提出了一种改进的单幂次趋近律及一种积分滑模面, 结合非线性反馈思想, 将其应用于执行机构输出饱和情况下的控制问题中, 其抗干扰能力较强, 但趋近时间较长.
由于指数趋近律及单幂次滑模趋近律可调节的参数较少, 为使系统具备更优越的性能, 往往需要与神经网络、自适应及反馈控制等其他算法相结合, 尽管这些方法可以在一定程度上改善系统品质, 但仍然无法摆脱指数趋近律或单幂次趋近律固有的抖振及响应慢等问题.因此, 国内外学者致力于通过提高趋近律阶次来使系统具备强鲁棒性及响应快速性等性能[10-11].梅红等[11]首次提出了双幂次滑模趋近律并应用于机器人跟踪问题中, 使系统状态量快速收敛, 提高了趋近速率, 但并没有给出趋近速率的表达式或定性分析.张合新等[12]提出了一种改进的双幂次趋近律方法, 在双幂次趋近律后面加一指数项, 并将其运用于非线性单入单出系统仿真, 进一步加快了系统的趋近速率, 其将两个幂次项参数分别选为${\alpha > 1}$及${0 < \beta < 1}$, 即将系统状态分为大于1及小于1两部分分别进行速率调节, 但是, 当系统初始状态远大于1或远小于$-1$时, 根据指数函数性质, 趋近速率仍具备改善空间.
综上, 已有的趋近律在控制过程中仍然存在抖振[8, 13-14]及收敛速度有待加强[11-12, 15]的问题.本文提出一种多幂次滑模趋近律, 增加的第三个幂次项将系统的趋近过程在双幂次趋近律[12]的基础上进一步细分并根据细分后的趋近阶段进行有针对性的速率调节, 提高系统趋近过程全程速率.尤其当系统初始状态远大于1或远小于$-1$时, 第三个幂次项将使系统状态迅速收敛.此外, 在系统具备不确定性及受到外加干扰情况下, 给出了干扰稳定界.
1. 多幂次滑模趋近律
1.1 多幂次滑模趋近律设计
对于连续时间系统, 设计多幂次滑模趋近律如下:
$ \begin{align}\label{eq1} \dot{s}=-k_1|s|^\alpha {\rm sgn}(s)-k_2|s|^\beta {\rm sgn}(s)-k_3|s|^\gamma {\rm sgn}(s)-k_4s \end{align} $
(1) 其中, $k_1>0$, $k_2>0$, $k_3>0$, $k_4>0$, $\alpha >1$, $0 > \beta < 1$. $\gamma $的取值情况为
$ \begin{align}\label{eq2} \gamma =\begin{cases} \max\{\alpha, |s|\}, &|s|\geq1\\ \min\{\beta, |s|\}, &|s| < 1 \end{cases} \end{align} $
(2) 由指数函数性质可知, 当系统状态满足条件$|s| < 1$时, 趋近律(1)主要受$-k_2|s|^\beta {\rm sgn}(s) - k_3|s|^\gamma {\rm sgn}(s)$影响; 当系统状态满足条件$|s|\geq1$时, 趋近律(1)主要受$-k_1|s|^\alpha {\rm sgn}(s)$ $-$ $k_3|s|^\gamma {\rm sgn}(s)$影响.$\gamma $的取值(2)可保证系统在状态满足$|s|$ \onecolumn%这一行用于第一页末尾$ > $ $\alpha $或$|s| < \beta $时自适应地改变趋近律中的指数参数, 从而获得较快的收敛速率.
1.2 存在性及可达性证明
定理1.对于滑模趋近律(1), 系统状态$s$在其作用下可达到平衡点$s=0$.
证明.根据式(1)可得到关系式
$ s\dot{s}=-k_1|s|^{\alpha +1}-k_2|s|^{\beta +1}-k_3|s|^{\gamma +1}-k_4s^2\leq0 $
(3) 当且仅当$s=0$时, 有$s\dot{s}=0$.
根据连续系统滑模趋近律存在且可达性条件[16], 若满足$s\dot{s}$ $\leq$ $0$, 则所设计的滑模趋近律为存在且可达的, 即系统状态$s$在趋近律(1)作用下可达到平衡点$s=0$.
2. 多幂次滑模趋近律特性分析
2.1 稳态抖振分析
对于传统指数趋近律[4], 其表达式为
$ \dot{s}_\alpha =-k_\alpha s-\varepsilon _\alpha {\rm sgn}(s) $
(4) 其中, $k_\alpha >0$, $\epsilon _\alpha >0$.当$s=0^+$时, 式(4)可写为$\dot{s}_\alpha =-\epsilon _\alpha $, 即临近稳态时, 在正方向趋近稳态的系统状态以$\dot{s}_\alpha =-\epsilon _\alpha $的速率运动.当$s=0^-$时, 式(4)可写为$\dot{s}_\alpha =\epsilon _\alpha $, 即临近稳态时, 在负方向趋近稳态的系统状态以$\dot{s}_\alpha =\epsilon _\alpha $的速率运动.因此, 系统不会稳定在平衡点处, 而是在平衡点处进行幅值为$\varepsilon _\alpha $的抖振.
对于多幂次趋近律(1), 当$s=0^+$及$s=0^-$时, 式(1)均可写为$\dot{s}=0$, 即系统在临近稳态时不会产生抖振现象.
2.2 多幂次滑模趋近律趋近速率
定理2.设$s$初始状态为$s_0$, 则采用滑模趋近律(1)对应的系统到达滑模面的时间小于$[T_1+T_3+T_3+T_4]$, 其中, $T_1$, $T_2$, $T_3$, $T_4$为
$ \begin{align} T_1=&\ \frac{1}{(1-\alpha)k_4} \left[ {\ln}\left(\alpha ^{1-\alpha }+\dfrac{k_1}{k_4}{\rm sgn}(s_0)\right)-\right.\notag \\ & \left.{\ln}\left(\alpha ^{1-\alpha }\right)+\dfrac{k_1}{k_4}{\rm sgn}(s_0) \right]\notag\\[2mm] T_2=&\ \frac{1}{(1-\alpha )k_4} \left[ {\ln}\left(1+\dfrac{k_1+k_3}{k_4}{\rm sgn}(s_0)\right)-\right.\notag \\ & \left.{\ln}\left(\alpha ^{1-\alpha} +\dfrac{k_1+k_3}{k_4}{\rm sgn}(s_0)\right)\right]\notag\\[2mm] T_3=&\ \dfrac{1}{(\beta -1)k_4} \left[ {\ln}\left(\beta ^1-\beta \right)+\dfrac{k_2}{k_4}{\rm sgn}(s_0)-\right.\notag \\ & \left.{\ln}\left(1+\dfrac{k_2}{k_4}{\rm sgn}(s_0)\right)\right]\notag\\[2mm] T_4=&\ \frac{1}{(\beta -1)k_4} \left[ {\ln}\left(1+\dfrac{k_2+k_3}{k_4}{\rm sgn}(s_0)\right)-\right.\notag \\ & \left.{\ln}\left(\beta ^{1-\beta }+\dfrac{k_2+k_3}{k_4}{\rm sgn}(s_0)\right)\right] \end{align} $
证明.不妨设$s_0>\alpha >1$, 将系统的趋近过程分为四个阶段, 分别为
1)系统从$s_0$到达$s(t_1)=\alpha $
由于$\alpha >1$, $0 < \beta < 1$, $\gamma ={\max}\{{\alpha, |s|}\}$, 则$k_1|s|^\alpha {\rm sgn}(s)$ $+$ $k_3|s|^\gamma {\rm sgn}(s)>k_2|s|^\beta {\rm sgn}(s)$, 尤其当参数满足$k_1|s|^{\alpha -\beta }+k_3|s|^{\gamma -\beta }>10k_2$时, 有
$ k_1|s|^\alpha {\rm sgn}(s)+k_3|s|^\gamma {\rm sgn}(s)\gg k2|s|^\beta {\rm sgn}(s) $
成立, 此阶段趋近速率主要受$-k_1|s|^\alpha {\rm sgn}(s)-k_3|s|^\gamma {\rm sgn}(s)$影响, 则趋近律(1)可写为
$ \dot{s}=-k_1s^\alpha -k_3s^\gamma -k_4s $
(6) 下面将求解式(6), 含有两个幂次项的方程求解具有一定难度, 故将其考虑为$\dot{s}+k_4s=-k_1s^\alpha $及$\dot{s}+k_4s=-k_3s^\gamma $两个方程分别求解, 由于$\dot{s}+k_4s=-k_1s^\alpha $及$\dot{s}+k_4s=-k_3s^\gamma $的解分别表示系统仅受$-k_1|s|^\alpha {\rm sgn}(s)$及仅受$-k_3|s|^\gamma {\rm sgn}(s)$的影响情况下所需要的趋近时间, 因此, 趋近时间一定小于两方程解中的任一值.求解$\dot{s}+k_4s=-k_1s^\alpha $的具体过程如下:
$ s^{-\alpha }\frac{{\rm d}s}{{\rm d}t}+k_4s^{1-\alpha }=-k_1 $
(7) 设中间变量$y=s^{1-\alpha }$, 则式(6)可写为
$ \frac{{\rm d}y}{{\rm d}t}+(1-\alpha )k_4y=-(1-\alpha )k_1 $
(8) 对于一阶线性非齐次方程(8), 采用常数变易法来求解[17], 令$y=u(t){\rm e}^{-\int(1-\alpha )k_4 {\rm d}t}$, 其中, $u=u(t)$为$t$的连续函数, 代入式(8)中得到
$ frac{{\rm d}u}{{\rm d}t}=-(1-\alpha )k_1{\rm e}^{-\int(1-\alpha )k_4 {\rm d}t} $
(9) 对式(9)两边积分可得
$ u=-\int{(1-\alpha )k_1{\rm e}^{-\int(1-\alpha )k_4 {\rm d}t} {\rm d}t}+c_1 $
(10) 其中, $c$为常数.则可求得式(8)的解为
$ \begin{array}{*{35}{l}} y=\left[-\int{(1-\alpha ){{k}_{1}}{{\text{e}}^{-\int{(1-\alpha )}{{k}_{4}}\text{d}t}}\text{d}t}+{{c}_{1}} \right]{{\text{e}}^{-\int{(1-\alpha )}{{k}_{4}}\text{d}t}}= \\ \ \ \ \ \ {{c}_{1}}{{\text{e}}^{-\int{(1-\alpha )}{{k}_{4}}\text{d}t}}-\frac{{{k}_{1}}{{\text{e}}^{\int{(1-\alpha )}{{k}_{4}}\text{d}t}}}{{{k}_{4}}}{{\text{e}}^{-\int{(1-\alpha )}{{k}_{4}}\text{d}t}}= \\ \ \ \ \ \ {{c}_{1}}{{\text{e}}^{-(1-\alpha ){{k}_{4}}t}}-\frac{{{k}_{1}}}{{{k}_{4}}} \\ \end{array} $
(11) 由于$y=s^{1-\alpha}$, 则有
$ s^{1-\alpha}=c_1{\rm e}^{-(1-\alpha)k_4t}-\frac{k_1}{k_4} $
(12) 由于$t=0$时, $s=s_0$, 则可求得常数$c_1$为
$ c_1=s_0^{1-\alpha}\frac{k_1}{k_4} $
(13) 结合式(12)及式(13), 可求出式(7)的解为
$ t=\frac{1}{(1-\alpha )k_4}\left\{{\ln}\left(s^{1-\alpha}+\frac{k_1}{k_4}\right)-{\ln}\left(s_0^{1-\alpha}+\frac{k_1}{k_4}\right)\right\} $
(14) 则由$s_0$到达$s(t_1)$所需时间为
$ t_1=\frac{1}{(1-\alpha )k_4}\left\{{\ln}\left(\alpha ^{1-\alpha}+\frac{k_1}{k_4}\right)-{\ln}\left(s_0^{1-\alpha}+\frac{k_1}{k_4}\right)\right\} $
(15) 因此, 系统由$s_0$到达$s(t_1)$所需时间小于$t_1$.
2)系统从$s(t_1)=\alpha $到达$s(t_2)=1$
趋近律(1)可写为
$ \dot{s}=-k_1|s|^\alpha {\rm sgn}(s)-k_3|s|^\gamma {\rm sgn}(s)-k_4s $
(16) 这一阶段由于$|s| < \alpha $, 因此, 有$\gamma =\alpha $, 则式(16)可写为
$ \dot{s}+k_4s=-(k_1+k_3)|s|^\alpha {\rm sgn}(s) $
(17) 对应的趋近时间可由上一阶段推导过程得出, 所需时间为
$ t=\frac{1}{(1-\alpha )k_4}\left[{\ln}\left(s^{1-\alpha}+\frac{k_1+k_3}{k_4}\right)-{\ln}\left(\alpha ^{1-\alpha}+\frac{k_1}{k_4}\right)\right] $
(18) 可以得出, 系统由$s(t_1)=\alpha $到达$s(t_2)=1$所需时间为
$ t= \frac{1}{(1-\alpha )k_4}\left[{\ln}\left(1+\frac{k_1+k_3}{k_4}\right)-{\ln}\left(\alpha ^{1-\alpha}+\frac{k_1+k_3}{k_4}\right)\right] $
(19) 因此, 系统由$s(t_1)$到达$s(t_2)=1$所需时间小于$t_2$.
3)系统从$s(t_2)=1$到达$s(t_3)=\beta $
由于$\alpha >1$, $0 < \beta < 1$, $\gamma ={\min}\{\beta, |s|\}$, 则$k_3|s|^\gamma {\rm sgn}(s)$ $+$ $k_2|s|^\beta {\rm sgn}(s)>k_1|s|^\alpha {\rm sgn}(s)$, 尤其当参数满足$k_3|s|^{\gamma -\alpha }+k_2|s|^{\beta -\alpha }>10k_1$时, 有
$ k_3|s|^\gamma {\rm sgn}(s)+k_2|s|^\beta {\rm sgn}(s)\gg k1|s|^\alpha {\rm sgn}(s) $
成立, 此阶段趋近速率主要受$-k_3|s|^\gamma {\rm sgn}(s)-k_2|s|^\beta {\rm sgn}(s)$影响, 则趋近律(1)可写为
$ \dot{s}=-k_2|s|^\beta -k_3|s|^\gamma -k_4s $
(20) 下面求解式(20), 同样将其考虑为$\dot{s}+k_4s=-k_2s^\beta $及$\dot{s}$ $+$ $k_4s$ $=$ $-k_3s^\gamma $两个方程分别求解, 由于$\dot{s}+k_4s=-k_2s^\beta $及$\dot{s}$ $+$ $k_4s=-k_3s^\gamma $的解分别表示系统仅受$-k_2|s|^\beta {\rm sgn}(s)$及仅受$-k_3|s|^\gamma {\rm sgn}(s)$的影响情况下所需要的趋近时间, 因此, 趋近时间一定小于两方程解中的任一值.求解$\dot{s}+k_4s=-k_2s^\beta $的具体过程如下.
$ s^{1-\beta }=c_2{\rm e}^{-(1-\beta )k_4t}-\frac{k_2}{k_4} $
(21) 其中, $t=0$时, $s=1$, 则可求得常数$c_2$为
$ c_2=1+\frac{k_2}{k_4} $
(22) 结合式(21)及式(22), 可求得趋近时间为
$ t=\frac{1}{(\beta -1)k_4}\left[{\ln}\left(s^{1-\beta }+\frac{k_2}{k_4}\right)-{\ln}\left(1+\frac{k_2}{k_4}\right)\right] $
(23) 则系统由$s(t_2)=1$到达$s(t_3)=\beta $所需时间为
$ t_3=\frac{1}{(\beta -1)k_4}\left[{\ln}\left(\beta ^{1-\beta }+\frac{k_2}{k_4}\right)-{\ln}\left(1+\frac{k_2}{k_4}\right)\right] $
(24) 因此, 系统由$s(t_2)=1$到达$s(t_3)=\beta $所需时间小于$t_3$.
4)系统从$s(t_3)=\beta $到达$s(t_4)=0 $
趋近律(1)可写为
$ \dot{s}=-k_2|s|^\beta {\rm sgn}(s) -k_3|s|^\gamma {\rm sgn}(s)-k_4s $
(25) 这一阶段由于$|s|>\beta $, 因此, 有$\gamma =\beta $, 则式(16)可写为
$ \dot{s}+k_4s=-(k_2+k_3)|s|^\beta {\rm sgn}(s) $
(26) 对应的趋近时间可由上一阶段推导过程得出, 所需时间为
$ \begin{align}\label{eq26} t=&\ \frac{1}{(\beta -1)k_4}\left[{\ln}\left(s^{1-\beta }+\frac{k_2+k_3}{k_4}\right)-\right.\notag \\[0mm] & \left.{\ln}\left(\beta ^{1-\beta }+\frac{k_2+k_3}{k_4}\right)\right] \end{align} $
(27) 则系统由$s(t_3)=\beta $到达$s(t_4)=0 $所需时间为
$ \begin{align}\label{eq27} t_4=&\ \frac{1}{(\beta -1)k_4}\left[{\ln}\left(1+\frac{k_2+k_3}{k_4}\right)-\right.\notag \\[0mm] & \left.{\ln}\left(\beta ^{1-\beta }+\frac{k_2+k_3}{k_4}\right)\right] \end{align} $
(28) 因此, 系统由$s(t_3)=\beta $到达$s(t_4)=0 $所需时间小于$t_4$.
此外, 若$s_0 < -\alpha <-1$, 则系统由$s_0$到达$s(t_4^*)=0$所需时间小于$[t_1^*+t_2^*+t_3^*+t_4^*]$, 其中,
$ t_{1}^{*}=\frac{1}{(1-\alpha ){{k}_{4}}}\left[\ln \left( {{\alpha }^{1-\alpha }}-\frac{{{k}_{1}}}{{{k}_{4}}} \right)-\ln \left( s_{0}^{1-\alpha }+\frac{{{k}_{1}}}{{{k}_{4}}} \right) \right] $
(29) $ \begin{align} & t_{2}^{*}=\frac{1}{(1-\alpha ){{k}_{4}}}\left[\ln \left( 1-\frac{{{k}_{1}}+{{k}_{3}}}{{{k}_{4}}} \right)-\right. \\ & \left. \ \ \ \ \ \ \ln \left( {{\alpha }^{1-\alpha }}-\frac{{{k}_{1}}+{{k}_{3}}}{{{k}_{4}}} \right) \right] \\ \end{align} $
(30) $ t_{3}^{*}=\frac{1}{(\beta-1){{k}_{4}}}\left[\ln \left( {{\beta }^{1-\beta }}+\frac{{{k}_{2}}}{{{k}_{4}}} \right)-\ln \left( 1+\frac{{{k}_{2}}}{{{k}_{4}}} \right) \right] $
(31) $ \begin{align} & t_{4}^{*}=\frac{1}{(\beta-1){{k}_{4}}}\left[\ln \left( 1-\frac{{{k}_{2}}+{{k}_{3}}}{{{k}_{4}}} \right)-\right. \\ & \left. \ \ \ \ \ \ \ln \left( {{\beta }^{1-\beta }}-\frac{{{k}_{2}}+{{k}_{3}}}{{{k}_{4}}} \right) \right] \\ \end{align} $
(32) 因此, 将式(15)与式(29)、式(19)与式(30)、式(24)与式(31)、式(28)与式(32)整理合并, 写为如下等式:
$ \begin{align} & {{T}_{1}}=\frac{1}{(1-\alpha ){{k}_{4}}}\left[\ln \left( {{\alpha }^{1-\alpha }}+\frac{{{k}_{1}}}{{{k}_{4}}}\text{sgn}({{s}_{0}}) \right)-\right. \\ & \ \ \ \ \ \left. \ln \left( s_{0}^{1-\alpha }+\frac{{{k}_{1}}}{{{k}_{4}}}\text{sgn}({{s}_{0}}) \right) \right] \\ \end{align} $
(33) $ \begin{array}{l} {T_2} = \frac{1}{{(1- \alpha ){k_4}}}\left[{\ln \left( {1 + \frac{{{k_1} + {k_3}}}{{{k_4}}}{\rm{sgn}}({s_0})} \right)-} \right.\\ \left. {\;\;\;\;\;\;\ln \left( {{\alpha ^{1-\alpha }} + \frac{{{k_1} + {k_3}}}{{{k_4}}}{\rm{sgn}}({s_0})} \right)} \right] \end{array} $
(34) $ \begin{array}{l} {T_3} = \frac{1}{{(\beta- 1){k_4}}}\left[{\ln \left( {{\beta ^{1-\beta }} + \frac{{{k_2}}}{{{k_4}}}{\rm{sgn}}({s_0})} \right)-} \right.\\ \left. {\;\;\;\;\;\;\ln \left( {1 + \frac{{{k_2}}}{{{k_4}}}{\rm{sgn}}({s_0})} \right)} \right] \end{array} $
(35) $ \begin{array}{l} {T_4} = \frac{1}{{(\beta- 1){k_4}}}\left[{\ln \left( {1 + \frac{{{k_2} + {k_3}}}{{{k_4}}}{\rm{sgn}}({s_0})} \right)-} \right.\\ \;\;\;\;\;\;\left. {\ln \left( {{\beta ^{1-\beta }} + \frac{{{k_2} + {k_3}}}{{{k_4}}}{\rm{sgn}}({s_0})} \right)} \right] \end{array} $
(36) 综上所述, 对于滑模趋近律(1), 系统状态的收敛时间小于$[T_1+T_2+T_3+T_4]$. $\square$
上文给出了趋近速率的具体求解过程, 亦可利用Matlab的dsolve语句求解.
2.3 趋近速率分析
现阶段幂次趋近律中趋近速率较快的为张合新等提出的改进的双幂次趋近律[12], 如式(37)所示.此处将多幂次趋近律与其对比, 以分析趋近速率的快速性.
$ \begin{align}\label{eq36} \dot{s}_c=-k_{c1}|s|^{\alpha_c} {\rm sgn}(s)-k_{c2}|s|^{\beta_c} {\rm sgn}(s) \end{align} $
(37) 其中, $k_{c1}>0$, $k_{c2}>0$, $0 <\beta_c <1$, $\alpha_c>1$.
双幂次趋近律设计的主要思想为利用指数函数的性质来对趋近过程的不同阶段针对性地加快趋近速率.多幂次趋近律(1)在此基础上进一步加快了趋近速率.双幂次趋近律[12]在$|s|>1$和$|s| < 1$范围内用常值幂次项系数$\alpha_c$和$\beta_c$进行调节, 趋近速率随这两个系数的选取而变化.针对系统状态与平衡点的情况, 多幂次趋近律引入的第三个幂次项系数$\gamma $进一步提升了趋近速率.根据式(2)可知, 当系统状态满足$|s|$ $>$ $\alpha $时, $\gamma =|s|$, 即趋近律的指数项大于双幂次趋近律的指数项, 从而实现了在$s>\alpha $阶段加快趋近速率的目的.同样地, 在$|s| < \beta $阶段趋近速率也高于双幂次趋近律.
值得说明的是, 如果增大双幂次趋近律[12]的参数$k_{c1}$及$k_{c2}$也可以在一定程度上加快系统趋近速率, 但由于其在$|s|$ $>$ $\alpha$及$|s|<\beta $阶段并没有进行针对性的速率调节, 而是将其分别与$1<|s|\leq\alpha $及$\beta \leq|s|<1$阶段组成一个整体, 即$|s|>1$及$|s|<1$阶段, 进行统一的速率调节, 因此, 在改善趋近速率方面具备局限性.然而, 多幂次趋近律(1)将系统趋近过程分为四个阶段, 即$|s|<\beta$, $\beta \leq|s|<1$, $1<|s|\leq\alpha$及$|s|>\alpha$ (此处说明的是系统所能经历最多的阶段, 即初始状态满足$|s_0|>\alpha >1$, 若系统初始状态不满足该条件, 则所经历的阶段为上述四阶段的子集), 并针对四个阶段分别设计了趋近幂次项的系数, 保证了系统在各个趋近阶段均具备保持较快趋近速率的自适应能力, 并且幂次项系数$\gamma $在$|s|>\alpha $取为$|s|$, 导致控制量随系统收敛而逐渐减小, 在提高趋近速率的同时也满足了实际系统对控制量幅值的限制要求.
3. 多幂次滑模趋近律干扰稳定界分析
上文证明了多幂次滑模趋近律可使系统在有限时间内到达平衡点, 且给出了趋近时间, 但当系统存在不确定性及外加干扰时, 该趋近律仅可使系统收敛到平衡点的一个邻域内, 称为干扰稳定界.
引理1.令x ∈ D ⊂Rn, $\dot{x}=f(x)$, $f:\boldsymbol R^n$ $\rightarrow$ ${{\boldsymbol R}^n}$为定义在平衡点邻域$D$内的连续函数, 假设存在连续函数$V$满足以下条件, 则函数$\dot{x}=f(x)$关于平衡点有限时间收敛[18]:
1) $V$是正定的;
2) $\dot{V}$除平衡点外是负定的;
3)存在实数$\varepsilon >0, ~~0 < \lambda < 1$和一邻域$N\subset{D}$使得$\dot{V}$ $+$ $\varepsilon V^\lambda$ $\leq$ $0$.
定理3.考虑如下存在不确定性及外加干扰的系统
$ \begin{align}\label{eq37} \dot{s}= &-k_1|s|^\alpha {\rm sgn}(s)-k_2|s|^\beta {\rm sgn}(s) -\notag\\ &\ k_3|s|^\gamma {\rm sgn}(s)-k_4s+{\rm d}(t) \end{align} $
(38) 若满足$|d(t)|\leq\delta $, $\delta >0$为常数, 则系统(38)的状态在有限时间内收敛到以下区域
$ |s| \le \min \left\{ {{{\left( {\frac{\delta }{{{k_1}}}} \right)}^{\frac{1}{\alpha }}}, {{\left( {\frac{\delta }{{{k_2}}}} \right)}^{\frac{1}{\beta }}}, {{\left( {\frac{\delta }{{{k_3}}}} \right)}^{\frac{1}{\gamma }}}} \right\} $
(39) $ \begin{array}{l} |\dot s| \le \min \left\{ {\delta, {{\left( {{k_1}\frac{\delta }{{{k_2}}}} \right)}^{\frac{\alpha }{\beta }}}, {{\left( {{k_1}\frac{\delta }{{{k_3}}}} \right)}^{\frac{\alpha }{\gamma }}}} \right\} + \\ \;\;\;\;\;\;\min \left\{ {\delta, {{\left( {{k_2}\frac{\delta }{{{k_1}}}} \right)}^{\frac{\beta }{\alpha }}}, {{\left( {{k_2}\frac{\delta }{{{k_3}}}} \right)}^{\frac{\beta }{\gamma }}}} \right\} + \\ \;\;\;\;\;\;\min \left\{ {\delta, {{\left( {{k_3}\frac{\delta }{{{k_1}}}} \right)}^{\frac{\gamma }{\alpha }}}, {{\left( {{k_3}\frac{\delta }{{{k_2}}}} \right)}^{\frac{\gamma }{\beta }}}} \right\} + \\ \;\;\;\;\;\;{k_4}\min \left\{ {{{\left( {\frac{\delta }{{{k_1}}}} \right)}^{\frac{1}{\alpha }}}, {{\left( {\frac{\delta }{{{k_2}}}} \right)}^{\frac{1}{\beta }}}, {{\left( {\frac{\delta }{{{k_3}}}} \right)}^{\frac{1}{\gamma }}}} \right\} + \delta \end{array} $
(40) 证明.选取Lyapunov函数为式(40)形式, 则有
$ \begin{align}\label{eq40} V=\frac{1}{2}s^2 \end{align} $
(41) 对其求导,
$ \begin{align} &\dot{V}=s\dot{s}=\notag\\ &\ -k_1|s|^{\alpha +1}-k_2|s|^{\beta +1}-k_3|s|{\gamma+1}-k_4s^2+s\times {\rm d}(t)\leq\notag\\ &\ -k_1|s|^{\alpha +1}-k_2|s|^{\beta+1}-k_3|s|{\gamma +1}-k_4s^2+|s|\!\times\!|{\rm d}(t)|\!\leq\notag\\ &\ -k_1|s|^{\alpha+1}-k_2|s|^{\beta +1}-k_3|s|{\gamma +1}-k_4s^2+|s|\times\delta \end{align} $
上式可变形为
$ \begin{align}\label{eq42} \dot{V}+k_4s^2\leq-k_1|s|^{\alpha +1}-k_2|s|^{\beta +1}-k_3|s|{\gamma +1}+|s|\times\delta \end{align} $
(43) 若$k_1|s|^{\alpha +1}-k_2|s|^{\beta +1}-k_3|s|{\gamma +1}\geq\delta$成立, 则
$ \dot{V}+2k_4V\leq0 $
(44) 进一步地, 若满足$k_1|s|^\alpha\geq\delta $, $k_2|s|^\beta\geq\delta$或$k_2|s|^\gamma\geq\delta$其中之一, 亦有式(37)成立, 而对应的区域分别为
$ {|s| \le {{\left( {\frac{\delta }{{{k_1}}}} \right)}^{\frac{1}{\alpha }}}} $
(45) $ {|s| \le {{\left( {\frac{\delta }{{{k_2}}}} \right)}^{\frac{1}{\beta }}}} $
(46) $ {|s| \le {{\left( {\frac{\delta }{{{k_3}}}} \right)}^{\frac{1}{\gamma }}}} $
(47) 因此, 由引理1可得, 系统状态$S$在有限时间内收敛到区域(39)内.将式(39)代入式(38)中, 可得
$ \begin{array}{l} |\dot s| = {k_1}|s{|^\alpha }{\rm{sgn}}\left( s \right) + {k_2}|s{|^\beta }{\rm{sgn}}\left( s \right) + {k_3}|s{|^\gamma }{\rm{sgn}}\left( s \right){\rm{ + }}\\ \quad {k_4}s + d\left( t \right) \le {k_1}|s{|^\alpha } + {k_2}|s{|^\beta } + {k_3}|s{|^\gamma } + {k_4}|s| + |d\left( t \right)| \le \\ \quad {k_1}\min {\left\{ {{{\left( {\frac{\delta }{{{k_1}}}} \right)}^{\frac{1}{\alpha }}}, {{\left( {\frac{\delta }{{{k_2}}}} \right)}^{\frac{1}{\beta }}}, {{\left( {\frac{\delta }{{{k_3}}}} \right)}^{\frac{1}{\gamma }}}} \right\}^\alpha } + \\ \quad {k_2}\min {\left\{ {{{\left( {\frac{\delta }{{{k_1}}}} \right)}^{\frac{1}{\alpha }}}, {{\left( {\frac{\delta }{{{k_2}}}} \right)}^{\frac{1}{\beta }}}, {{\left( {\frac{\delta }{{{k_3}}}} \right)}^{\frac{1}{\gamma }}}} \right\}^\beta } + \\ \quad {k_3}\min {\left\{ {{{\left( {\frac{\delta }{{{k_1}}}} \right)}^{\frac{1}{\alpha }}}, {{\left( {\frac{\delta }{{{k_2}}}} \right)}^{\frac{1}{\beta }}}, {{\left( {\frac{\delta }{{{k_3}}}} \right)}^{\frac{1}{\gamma }}}} \right\}^\gamma } + \\ \quad {k_4}\min \left\{ {{{\left( {\frac{\delta }{{{k_1}}}} \right)}^{\frac{1}{\alpha }}}, {{\left( {\frac{\delta }{{{k_2}}}} \right)}^{\frac{1}{\beta }}}, {{\left( {\frac{\delta }{{{k_3}}}} \right)}^{\frac{1}{\gamma }}}} \right\} + \delta = \\ \min \left\{ {\delta, {{\left( {{k_1}\frac{\delta }{{{k_2}}}} \right)}^{\frac{\alpha }{\beta }}}, {{\left( {{k_1}\frac{\delta }{{{k_3}}}} \right)}^{\frac{\alpha }{\gamma }}}} \right\} + \\ \min \left\{ {\delta, {{\left( {{k_2}\frac{\delta }{{{k_1}}}} \right)}^{\frac{\beta }{\alpha }}}, {{\left( {{k_2}\frac{\delta }{{{k_3}}}} \right)}^{\frac{\beta }{\gamma }}}} \right\} + \\ \min \left\{ {\delta, {{\left( {{k_3}\frac{\delta }{{{k_1}}}} \right)}^{\frac{\gamma }{\alpha }}}, {{\left( {{k_3}\frac{\delta }{{{k_2}}}} \right)}^{\frac{\gamma }{\beta }}}} \right\} + \\ {k_4}\min \left\{ {{{\left( {\frac{\delta }{{{k_1}}}} \right)}^{\frac{1}{\alpha }}}, {{\left( {\frac{\delta }{{{k_2}}}} \right)}^{\frac{1}{\beta }}}, {{\left( {\frac{\delta }{{{k_3}}}} \right)}^{\frac{1}{\gamma }}}} \right\} + \delta \end{array} $
(48) $\square$
4. 数值仿真
由于航天器姿态机动控制问题中动态响应的快速性及平滑性是主要考虑的性能指标, 因此数值仿真以航天器姿态机动控制系统为例[19], 选取俯仰轴单轴机动情况进行仿真.将其考虑为非线性系统
$ \begin{align}\label{eq48} \ddot{\vartheta}=&\ \frac{\cos{\gamma}}{J_z}u+f(t)\\ f(t)=&\ \dot{\omega_y}\sin{\gamma} +\omega_y\dot{\gamma}\cos\gamma +\omega_z\dot{\gamma}\sin\gamma -\notag\\ &\ \frac{\cos\gamma}{J_z}(J_y-J_x)\omega_x\omega_y \end{align} $
(49) 其中, $f(t)$表示除机动轴外其余两轴带来的耦合影响, $\vartheta$, $\psi$, $\gamma$分别为飞行器俯仰角, 偏航角及滚转角, $\omega_x$, $\omega_y$, $\omega_z$分别为飞行器本体坐标系相对地面坐标系旋转角速度在本体坐标系$Ox$, $Oy$及$Oz$轴的分量, $J_x$, $J_y$, $J_z$分别为飞行器绕本体坐标系$Ox$, $Oy$及$Oz$轴的转动惯量, $u$为航天器本体坐标系下机动轴的外加力矩.将耦合项$f(t)$, 系统不确定性及外加干扰统一考虑为系统的干扰${ d}(t)$, 则航天器单轴姿态机动的模型为
$ \begin{align}\label{eq50} \ddot{\vartheta}=\frac{\cos\gamma}{J_z}u+{ d}(t) \end{align} $
(51) 数值仿真主要分为两部分: 1)在无干扰情况下, 验证本文多幂次趋近律的有效性及优越性; 2)在干扰作用下, 验证定理3给出干扰稳定界的正确性.在数值仿真中, 滑模面统一选择较为简单的形式[8]: $\sigma =\lambda e$, 其中, $\lambda >0$, $e=\vartheta_d-\vartheta$为航天器俯仰角跟踪误差, $\vartheta_d$为俯仰角控制指令.
4.1 各趋近律对比仿真
当${ d}(t)=0$时, 分别采用指数趋近律、单幂次趋近律、双幂次趋近律以及本文提出的多幂次趋近律设计控制律$u$, 并进行对比仿真, 结果如图 1、图 2及图 3所示.
a)指数趋近律
$ \dot{s}_a=-k_as-\varepsilon_a{\rm sgn}(s) $
(52) b)单幂次趋近律
$ \begin{align}\label{eq52} \dot{s}_b=-k_b|s|^{\alpha_b}-{\rm sgn}(s) \end{align} $
(53) c)双幂次趋近律
$ \begin{align}\label{eq53} \dot{s}_c=-k_{c1}|s|^{\alpha_c}{\rm sgn}(s)-k_{c2}|s|^{\beta_c}{\rm sgn}(s) \end{align} $
(54) d)多幂次趋近律
$ \begin{align}\label{eq54} \dot{s}=-k_1|s|^{\alpha}{\rm sgn}(s)-k_2|s|^{\beta}{\rm sgn}(s)-k_3|s|^{\gamma}{\rm sgn}(s)-k_4s \end{align} $
(55) 趋近律参数选择如表 1所示.
表 1 各趋近律仿真参数Table 1 Simulation parameters of every reaching lawa b c d ka=0.9 kb=0.7 kc1=0.7 k1=0.7 εa=0.5 αb=1.1 kc2=0.9 k2=0.9 αc=1.1 k3=1.1 βc=0.8 k4=1.4 α=1.7 β=0.5 控制目标为令俯仰角单轴大角度机动至90$^\circ$即$\pi /2 {\rm rad}$, 即指令为$\vartheta_d={\pi}/{2} {\rm rad}$, 滚转角及偏航角在此过程中保持为0 rad.航天器姿态跟踪系统参数为$J_z=324 {\rm kg\cdot m^2}$, 初始状态$s(0)=6$, $\dot{s}(0)=6$, $\vartheta_d(0)=0 {\rm rad}$, $\dot{\vartheta_d}(0)=0 {\rm rad/s}$.
由图 1可见, 指数趋近律下的姿态跟踪在到达期望角度后, 姿态角在期望值附近往返跳动, 形成抖振.单幂次趋近律及双幂次趋近律在趋近速率方面有所提升, 且双幂次趋近律下的跟踪速率较快.多幂次趋近律下的航天器单轴姿态跟踪速度较其他三种趋近律快, 且跟踪上期望信号后可稳定在期望角度, 不发生抖振.
由图 2可见, 指数趋近律(52)在趋近阶段动态响应平滑, 但从初值到达滑模面过程时间较长, 可通过不断修正参数$k_a$及$\varepsilon_a$改善趋近过程动态品质, 因常数项$-\varepsilon_a{\rm sgn}(s)$的存在系统不会稳定在平衡零点处, 而是在平衡零点附近产生幅值大小为$\varepsilon_a$的抖振, 而其余三个趋近律均不含此项, 故可稳定在平衡零点处, 不产生抖振现象, 图 3可验证此结论.由图 2和图 3可见, 单幂次趋近律(53)在系统状态远离滑模面时收敛速度较小, 双幂次趋近律(54)在一定程度上增大了收敛速度, 这是由于双幂次趋近律具备两个可调节的幂次项系数, 使系统在接近滑模面及远离滑模面时均具备较快的收敛速度.本文提出的多幂次趋近律具备最快的收敛速度, 且系统可收敛于平衡零点, 无抖振现象.这是由于多幂次趋近律增加了可调节的幂次项系数个数, 且增加了一阶项$k_4s$, 加快了系统在远离滑模面时的收敛速度.此外, 由图 2可见, 在航天器做单轴机动运动时, 基于多幂次滑模趋近律的控制器输出最大幅值较其他三种趋近律小, 对控制器的最大输出值要求较低.
4.2 干扰稳定界仿真
假设系统的不确定性及所受外界干扰可统一写为
$ d(t) = 0.5\cos t + \sin 2t $
(56) 由定理3可得$\delta\geq|{\rm d}(t)|=\sqrt{1.25}$, 取$\delta =\sqrt{1.25}$.趋近律参数选为$\alpha =2$, $\beta =0.5$, $\gamma =-2$, $k_1=0.7$, $k_2=0.57$, $k_3$ $=$ $0.52$, $k_4=0.81$.代入式(38)及式(39), 可得$|s|\leq 0.682$, $|\dot{s}|$ $\leq$ $2.559$.仿真结果如图 4及图 5所示.
由图 4和图 5可见, 在系统存在不确定性或受到外加干扰作用时, 系统稳态过程中无法收敛至平衡零点, 而是在平衡零点附近波动, 造成稳态误差, 且稳态误差最大幅值分别满足$|s|\leq0.682$, $|\dot{s}|\leq2.559$.与本文定理3给出的干扰误差界一致.因此, 在系统存在不确定性或外加干扰时, 在多幂次趋近律作用下系统的稳态误差一定在本文给出的干扰误差界内.
5. 结论
本文针对航天器姿态机动等对响应快速性具备较高要求的系统提出了一种多幂次趋近律, 理论证明了其存在性、可达性及稳定性, 并详细推导了趋近速率及干扰稳定界, 仿真验证了该趋近律的有效性及优越性, 且检验了干扰稳定界的推导结果.虽然在受到扰动时, 系统状态可收敛在干扰稳定界内, 但其仍然不能稳定在平衡零点, 可以考虑神经网络或干扰观测器来克服系统的受扰影响, 此外, 多幂次趋近律中参数之间满足某种关系时将会使系统的干扰稳定界缩小, 从而使得系统具备更好的抗干扰性, 以上内容有待进一步研究.
-
图 3 开工时间延迟扰动下的炼钢–连铸生产重调度策略
Fig. 3 Rescheduling strategy for the steelmaking and continuous casting
图 12 不同重调度方法的设备负荷均衡率比较
Fig. 12 Comparison of load balancing rate of equipment for methods
表 1 各种重调度方法
Table 1 Rescheduling methods
序号 重调度排序 指派规则优先级 重调度方法 1 Seq1 Assign1 H11 2 Seq1 Assign2 H12 3 Seq1 Assign3 H13 4 Seq1 Assign4 H14 5 Seq1 Assign5 H15 6 Seq1 Assign6 H16 7 Seq2 Assign1 H21 8 Seq2 Assign2 H22 9 Seq2 Assign3 H23 10 Seq2 Assign4 H24 11 Seq2 Assign5 H25 12 Seq2 Assign6 H26 13 Seq3 Assign1 H31 14 Seq3 Assign2 H32 15 Seq3 Assign3 H33 16 Seq3 Assign4 H34 17 Seq3 Assign5 H35 18 Seq3 Assign6 H36 
下载: 导出CSV
表 2 各种启发式方法的RPI值
Table 2 RPI values of rescheduling methods
问题 参数 H11 H12 H13 H14 H15 H16 H21 H22 H23 H24 H25 H26 H31 H32 H33 H34 H35 H36 P1 $C_{1}=1$ 2 3 322 327 33 41 0 0 132 134 2 3 2 2 132 140 7 10 P2 $C_{1}=2$ 4 13 102 010 101 702 186 270 6 5 88 560 88 515 6 13 10 15 88 386 88 847 54 51 P3 $C_{1}=3$ 30 71 877 070 875 109 422 522 19 20 819 627 819 112 26 37 30 27 820 876 819 229 96 98 P4 $C_{2}=1$ 12 65 744 525 735 497 800 818 0 0 694 356 693 395 12 28 12 21 695 642 693 220 188 190 P5 $C_{2}=2$ 45 77 885 908 887 138 244 542 0 0 817 767 817 774 3 17 24 26 818 366 818 702 44 43 P6 $C_{2}=3$ 1 1 242 279 5 5 8 12 230 249 2 2 18 6 125 247 2 3 P7 $C_{2}=4$ 1 1 1 111 1 106 13 16 21 16 1 015 1 020 23 25 9 10 1 016 1 024 16 15 P8 $C_{2}=5$ 0 1 553 542 6 7 15 15 497 497 16 18 8 9 508 501 11 14 P9 $C_{3}=1$ 29 58 694 030 693 295 303 417 19 19 654 301 654 034 11 22 28 28 654 200 653 089 99 97 P10 $C_{3}=2$ 6 25 272 029 270 402 242 313 6 6 244 777 244 773 10 17 12 15 245 962 245 842 48 53 P11 $C_{3}=3$ 1 5 13 345 13 440 96 104 0 0 9 241 9 254 13 15 1 1 9 232 9 248 10 9 P12 $C_{4}=1$ 2 1 544 354 541 868 464 477 9 10 504 493 504 403 12 15 6 5 505 085 504 079 115 113 P13 $C_{4}=2$ 33 84 386 237 386 413 117 257 9 7 358 574 357 953 13 21 23 30 359 090 358 857 27 25 P14 $C_{4}=3$ 1 2 48 812 48 857 60 99 7 9 45 253 45 404 9 18 13 9 45 219 45 280 15 22
下载: 导出CSV
-
[1] Ouelhadj D, Petrovic S. A survey of dynamic scheduling in manufacturing systems. Journal of Scheduling, 2009, 12(4):417-431 doi: 10.1007/s10951-008-0090-8 [2] Abumaizar R J, Svestka J A. Rescheduling job shops under random disruptions. International Journal of Production Research, 1997, 35(7):2065-2082 doi: 10.1080/002075497195074 [3] Sanmartí E, Espuña A, Puigianer L. Batch production and preventive maintenance scheduling under equipment failure uncertainty. Computer and Chemical Engineering, 1997, 21(10):1157-1168 doi: 10.1016/S0098-1354(96)00327-4 [4] Vieira G E, Herrmann J W, Lin E. Analytical models to predict the performance of a single-machine system under periodic and event-driven rescheduling strategies. International Journal of Production Research, 2000, 38(8):1899-1915 doi: 10.1080/002075400188654 [5] Roy R, Adesola B A, Thornton S. Development of a knowledge model for managing schedule disturbance in steel-making. International Journal of Production Research, 2004, 42(18):3975-3994 doi: 10.1080/00207540410001716453 [6] Dorn J. Case-based Reactive Scheduling. London, UK:Chapman & Hall, 1994. 32-50 [7] Dorn J. Cooperating scheduling systems. International Journal of Manufacturing Technology and Management, 2001, 3(6):570-585 doi: 10.1504/IJMTM.2001.001430 [8] Dorn J, Kerr R, Thalhammer G. Reactive scheduling:improving the robustness of schedules and restricting the effects of shop floor disturbances by fuzzy reasoning. International Journal of Human-Computer Studies, 1995, 42(6):687-704 doi: 10.1006/ijhc.1995.1031 [9] Ouelhadj D, Cowling P I, Petrovic S. Utility and stability measures for agent-based dynamic scheduling of steel continuous casting. In:Proceeding of the 2003 IEEE International Conference on Robotics & Automation. Taipei, China:IEEE, 2003. 175-180 [10] Ouelhadj D, Petrovic S, Cowling P I, Meisels A. Inter-agent cooperation and communication for agent-based robust dynamic scheduling in steel production. Advanced Engineering Informatics, 2004, 18(3):161-172 doi: 10.1016/j.aei.2004.10.003 [11] Guo D F, Li T K. Rescheduling algorithm for steelmaking-continuous casting. In:Proceeding of the 2nd IEEE Conference on Industrial Electronics and Applications. Harbin, China:IEEE, 2007. 1421-1425 [12] Li T K, Guo D F. Constraint-based approach for steelmaking-continuous casting rescheduling. In:Proceedings of the 20th International Conference on Industrial, Engineering and Other Applications of Applied Intelligent Systems. Kyoto, Japan:Springer, 2007. 1108-1117 [13] Zhang C S, Li T K, Wang B L, Zhang W X, Sun B. Dynamic modeling method for the scheduling problem in steelmaking-continuous casting with disturbance of product quality. Energy Procedia, 2011, 13:253-261 doi: 10.1016/S1876-6102(14)00454-8 [14] Worapradya K, Buranathiti T. Production rescheduling based on stability under uncertainty for continuous slab casting. In:Proceedings of ASIMMOD 2009. Bangkok, Thailand, 2009. 170-176 [15] Ozoe Y, Konishi M. Agent based scheduling of steel making processes. In:Proceedings of the 2009 International Conference on Networking, Sensing and Control. Okayama, Japan:IEEE, 2009. 278-281 [16] Chen K, Zheng Z, Liu Y, Gao X Q. Real-time scheduling method for steelmaking-continuous casting. In:Proceedings of the 2010 IEEE International Conference on Industrial Engineering and Engineering Management. Macao, China:IEEE, 2010. 2366-2370 [17] Mao K, Pan Q K, Pang X F, Chai T Y. An effective Lagrangian relaxation approach for rescheduling a steelmaking-continuous casting process. Control Engineering Practice, 2014, 30:66-77 [18] Mao K, Pan Q K, Pang X F, Chai T Y. A novel Lagrangian relaxation approach for a hybrid flowshop scheduling problem in the steelmaking-continuous casting process. European Journal of Operational Research, 2014, 236(1):51-60 doi: 10.1016/j.ejor.2013.11.010 [19] Tang L X, Zhao Y, Liu J Y. An improved differential evolution algorithm for practical dynamic scheduling in steelmaking-continuous casting production. IEEE Transactions on Evolutionary Computation, 2014, 18(2):209-25 doi: 10.1109/TEVC.2013.2250977 [20] Kumar P R. Re-entrant lines. Queueing Systems, 1993, 13(1-3):87-110 doi: 10.1007/BF01158930 [21] 俞胜平, 郑秉霖, 柴天佑.炼钢连铸混合智能调度方法及其应用.华东理工大学学报(自然科学版), 2006, 32(7):844-848 http://www.cnki.com.cn/Article/CJFDTOTAL-HLDX200607024.htmYu Sheng-Ping, Zheng Bing-Lin, Chai Tian-You. Hybrid intelligent scheduling method and its application for steelmaking and continuous casting. Journal of East China University of Science and Technology (Natural Science Edition), 2006, 32(7):844-848 http://www.cnki.com.cn/Article/CJFDTOTAL-HLDX200607024.htm [22] Pang X F, Yu S P, Zheng B L, Chai T Y. Complete modification rescheduling method and its application for steelmaking and continuous casting. In:Proceedings of the 17th IFAC World Congress. Seoul, Korea:IFAC, 2008. 12837-12842 -
下载:
计量
- 文章访问数: 2510
- HTML全文浏览量: 149
- PDF下载量: 1030
- 被引次数: 0
