CN106017118B - 一种带有输入时滞的电熔镁炉电极调节优化控制方法 - Google Patents
一种带有输入时滞的电熔镁炉电极调节优化控制方法 Download PDFInfo
- Publication number
- CN106017118B CN106017118B CN201610596894.4A CN201610596894A CN106017118B CN 106017118 B CN106017118 B CN 106017118B CN 201610596894 A CN201610596894 A CN 201610596894A CN 106017118 B CN106017118 B CN 106017118B
- Authority
- CN
- China
- Prior art keywords
- input delay
- electrode
- furnace
- phase electric
- mrow
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F27—FURNACES; KILNS; OVENS; RETORTS
- F27D—DETAILS OR ACCESSORIES OF FURNACES, KILNS, OVENS, OR RETORTS, IN SO FAR AS THEY ARE OF KINDS OCCURRING IN MORE THAN ONE KIND OF FURNACE
- F27D19/00—Arrangements of controlling devices
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F27—FURNACES; KILNS; OVENS; RETORTS
- F27D—DETAILS OR ACCESSORIES OF FURNACES, KILNS, OVENS, OR RETORTS, IN SO FAR AS THEY ARE OF KINDS OCCURRING IN MORE THAN ONE KIND OF FURNACE
- F27D19/00—Arrangements of controlling devices
- F27D2019/0028—Regulation
Abstract
本发明提供一种带有输入时滞的电熔镁炉电极调节优化控制方法,该方法首先建立三相电熔镁炉的带有输入时滞的电极调节系统线性模型和无输入时滞的电极调节系统线性模型;求解三相电熔镁炉的无输入时滞的电极调节系统线性模型,得到无输入时滞的电极调节系统的黎卡提矩阵以及电熔镁炉电极调节最优控制率;确定三相电熔镁炉的带有输入时滞的电熔镁炉电极调节最优控制器;在三相电熔镁炉熔化阶段工作过程中,将三相电熔镁炉设定熔化电流值与实际熔化电流值的偏差作为带有输入时滞的电熔镁炉电极调节最优控制器的输入,得到三相电机转速,控制电熔镁炉电机转子工作。该方法有效的改善了电熔镁炉电流回路的控制效果。
Description
技术领域
本发明属于电熔镁炉过程控制技术领域,具体涉及一种带有输入时滞的电熔镁炉电极调节优化控制方法。
背景技术
电熔镁砂又称电熔MgO,是一种重要的镁的氧化物,通常以菱镁矿石为主要生产原料。采用三相交流电熔镁炉进行熔炼,通过调整三相电极与熔池液面之间的位置来控制三相电极电流,使之产生电弧,通过电弧放热使炉内原料受热熔化形成熔液,熔液再经过冷却结晶后生成成品。
时滞对电熔镁砂品位的影响尤为值得关注。在三相电极熔化阶段,通过电机来调节电极位置进而改变电弧长度,使得实际熔化电流能够快速有效的跟踪理想熔化电流,但是由于电机齿轮与齿轮之间存在空隙,电机转子本身由于风阻、轴承摩擦等原因,总会使电机不能快速的调节电极的升降,即存在一定的时滞,这对电熔镁砂的品味高低有着很大的影响,实际电熔镁炉工业生产过程中,采用常用的PID控制算法,这是因为PID控制算法结构简单,参数调整相对来说非常容易,容易被工厂工人接受,但是PID控制算法忽略了时滞的因素,这在一定程度上大大降低了电熔镁砂品位。
发明内容
针对现有技术的不足,本发明提出一种带有输入时滞的电熔镁炉电极调节优化控制方法。
本发明的技术方案是:
一种带有输入时滞的电熔镁炉电极调节优化控制方法,包括以下步骤:
步骤1:根据三相电熔镁炉的系统状态变量、带有输入时滞的电机转速以及实际熔化电流的关系,建立三相电熔镁炉的带有输入时滞的电极调节系统线性模型和无输入时滞的电极调节系统线性模型;
步骤2:采用极小值原理求解三相电熔镁炉的无输入时滞的电极调节系统线性模型,得到无输入时滞的电极调节系统的黎卡提矩阵,从而确定无输入时滞的电熔镁炉电极调节最优控制率;
步骤3:根据最优控制理论,确定三相电熔镁炉的无输入时滞的电熔镁炉电极调节最优控制器;
步骤4:根据三相电熔镁炉的带有输入时滞的电极调节系统线性模型的二次性能指标和其对应的无输入时滞的电极调节系统线性模型的二次性能指标的数值相等,得到三相电熔镁炉的带有输入时滞的电熔镁炉电极调节最优控制器;
步骤5:在三相电熔镁炉熔化阶段工作过程中,将三相电熔镁炉设定熔化电流值与实际熔化电流值的偏差作为带有输入时滞的电熔镁炉电极调节最优控制器的输入,得到三相电机转速,控制电熔镁炉电机转子工作。
可选地,所述三相电熔镁炉的带有输入时滞的电极调节系统线性模型如下所示:
其中,x(t)为带有输入时滞的三相电熔镁炉的系统状态变量,A为三相电机传动系数矩阵,B为带有输入时滞的执行机构传动系数矩阵,u(t-τ)为带有输入时滞的三相电机转速,y(t)为三相电熔镁炉实际熔化电流值,C为电弧系数矩阵,t为时间,τ为时滞;
所述三相电熔镁炉的无输入时滞的电极调节系统线性模型如下所示:
其中,X(t)为无输入时滞的三相电熔镁炉的系统状态变量,B0=e-AτB为无输入时滞的执行机构传动系数矩阵,u(t)为无输入时滞的三相电机转速。
可选地,所述无输入时滞的电熔镁炉电极调节最优控制器如下式所示:
其中,为无输入时滞的三相电机最优转速,K0为无输入时滞的电熔镁炉电极调节最优控制率,X(t)为无输入时滞的三相电熔镁炉的系统状态变量,t为时间。
可选地,所述三相电熔镁炉的带有输入时滞的电熔镁炉电极调节最优控制器如下所示:
u*(t-τ)=-Kdx(t-τ);
其中,u*(t-τ)为带有输入时滞的三相电机最优转速,t为时间,τ为时滞,x(t-τ)为带有输入时滞的三相电熔镁炉的系统状态变量,K0为无输入时滞的电熔镁炉电极调节最优控制率,和为系统转移状态矩阵,A为三相电机传动系数矩阵,B为带有输入时滞的执行机构传动系数矩阵,K0为无输入时滞的电熔镁炉电极调节最优控制率,tf为终点时刻。
本发明的有益效果:
本发明提出一种带有输入时滞的电熔镁炉电极调节优化控制方法,该方法有效的改善了电熔镁炉电流回路的控制效果,进而提高了电熔镁砂的品位;由于电流控制更加稳定,在一定程度上避免了由于电流剧烈波动而引起的生产故障,从而确保了生产安全。
附图说明
图1为本发明实施方式中带有输入时滞的电熔镁炉电极调节优化控制方法的流程图;
图2为本发明实施方式中采用带有输入时滞的电熔镁炉电极调节最优控制器控制三相电熔镁炉熔化阶段工作过程的示意图;
图3为本发明实施方式中三相电机转速仿真图;
图4为本发明实施方式中三相电极实际电流值仿真图。
具体实施方式
下面结合附图对本发明具体实施方式加以详细说明。
本实施方式中,三相电熔镁炉熔化过程中,电极电压固定,三相电极电流与对应的电弧电阻在电熔镁炉工作点附近呈线性特征,电弧电阻的大小与三相电极的升降快慢成比例关系,电极的升降又与电机的转速成一个比例关系。因此,三相电极电流与三相电机的转速在工作点附近是一个线性的关系。
一种带有输入时滞的电熔镁炉电极调节优化控制方法,如图1所示,包括以下步骤:
步骤1:根据三相电熔镁炉的系统状态变量、带有输入时滞的电机转速以及实际熔化电流的关系,建立三相电熔镁炉的带有输入时滞的电极调节系统线性模型和无输入时滞的电极调节系统线性模型。
本实施方式中,在三相电熔镁炉熔化阶段,以三相电机转速为输入,三相电流为输出的电熔镁炉电极调节系统可以视为一个线性模型,该线性模型可采用常规的带有输入时滞的线性状态空间模型表示,得到建立的三相电熔镁炉的带有输入时滞的电极调节系统线性模型如式(1)所示:
其中,x=[x1 x2 x3 x4]T为系统状态变量,分别表示电熔镁炉的炉内燃料厚度、原料电阻率、熔池电阻率以及电极端部球面半径。x(t)为带有输入时滞的三相电熔镁炉的系统状态变量,A为三相电机传动系数矩阵,B为带有输入时滞的执行机构传动系数矩阵,u(t-τ)为带有输入时滞的三相电机转速,y(t)为三相电熔镁炉实际熔化电流值,C为电弧系数矩阵,t为时间,τ=0.5s为时滞。
本实施方式中,采用系统辨识方法求取三相电机传动系数矩阵A、带有输入时滞的执行机构传动系数矩阵B和电弧系数矩阵C。三相电机传动系数矩阵A如式(2)所示:
带有输入时滞的执行机构传动系数矩阵B如式(3)所示:
B=[1 0 0 0]T (3)
电弧系数矩阵C如式(4)所示:
C=[0 0 0 4.0482e+06] (4)
根据最优控制理论,可得对应公式如式(5)所示:
其中,X(t)为无输入时滞的三相电熔镁炉的系统状态变量,A(t′-t+τ)为坐标平移后的转移矩阵,u(t′)为三相电机转速,t′为时间。
将公式(5)左右两边同时对时间t求导,得到三相电熔镁炉的无输入时滞的电极调节系统线性模型如式(6)所示:
其中,B0=e-AτB为无输入时滞的执行机构传动系数矩阵,u(t)为无输入时滞的三相电机转速。
步骤2:采用极小值原理求解三相电熔镁炉的无输入时滞的电极调节系统线性模型,得到无输入时滞的电极调节系统的黎卡提矩阵,从而确定无输入时滞的电熔镁炉电极调节最优控制率。
本实施方式中,采用极小值原理求解三相电熔镁炉的无输入时滞的电极调节系统线性模型,得到无输入时滞的电极调节系统的黎卡提矩阵P0如式(7)所示:
确定无输入时滞的电熔镁炉电极调节最优控制率K0如式(8)所示:
K0=-R-1P0(t)X(t)=[0.0007 0.0056 0.0395 0] (8)
步骤3:根据最优控制理论,确定三相电熔镁炉的无输入时滞的电熔镁炉电极调节最优控制器。
本实施方式中,对于三相电熔镁炉的无输入时滞的电极调节系统线性模型,根据最优控制理论,确定三相电熔镁炉的无输入时滞的电熔镁炉电极调节最优控制器如式(9)所示:
其中,为无输入时滞的三相电机最优转速。
步骤4:根据三相电熔镁炉的带有输入时滞的电极调节系统线性模型的二次性能指标和其对应的无输入时滞的电极调节系统线性模型的二次性能指标的数值相等,得到三相电熔镁炉的带有输入时滞的电熔镁炉电极调节最优控制器。
本实施方式中,建立三相电熔镁炉的带有输入时滞的电极调节系统线性模型状态空间表达式和无输入时滞的电极调节系统线性模型状态空间表达式的对应关系如下:
对于三相电熔镁炉的带有输入时滞的电极调节系统线性模型,在给定初始状态x(0)=[1 1 1 1]T和外输入作用转速下,可以求出其状态方程的解,如式(10)所示:
将公式(10)中的t都用t+τ来代替,得到三相电熔镁炉的带有输入时滞的电极调节系统线性模型状态空间表达式和无输入时滞的电极调节系统线性模型状态空间表达式的等价关系式如式(11)所示:
通过最优控制理论,定义三相电熔镁炉的带有输入时滞的电极调节系统线性模型的二次性能指标J1如式(12)所示:
其中,tf为终点时刻。
通过最优控制理论,定义三相电熔镁炉的无输入时滞的电极调节系统线性模型的二次性能指标J2如式(13)所示:
其中,
将三相电熔镁炉的带有输入时滞的电极调节系统线性模型状态空间表达式和无输入时滞的电极调节系统线性模型状态空间表达式的等价关系式公式(11)代入三相电熔镁炉的带有输入时滞的电极调节系统线性模型的二次性能指标J1公式(12)中,得到如式(14)所示:
由公式(14)可知,根据三相电熔镁炉的带有输入时滞的电极调节系统线性模型的二次性能指标和其对应的无输入时滞的电极调节系统线性模型的二次性能指标的数值相等。
得到三相电熔镁炉的带有输入时滞的电熔镁炉电极调节最优控制器如式(15)所示:
u*(t-τ)=-Kdx(t-τ) (15)
其中,u*(t-τ)为带有输入时滞的三相电机最优转速,kd1=7.1136e-04*e-7.1136e-04t+(7.1136e-04-710.8)*0.5+1.688,kd2=7.1136e-04*e-0.0056t-63374.995+0.071,kd3=7.1136e-04*e-0.0395t-48484.98+1.0451, 为系统转移状态矩阵。
步骤5:在三相电熔镁炉熔化阶段工作过程中,将三相电熔镁炉设定熔化电流值与实际熔化电流值的偏差作为带有输入时滞的电熔镁炉电极调节最优控制器的输入,得到三相电机转速,控制电熔镁炉电机转子工作。
采用带有输入时滞的电熔镁炉电极调节最优控制器控制三相电熔镁炉熔化阶段工作过程如图2所示。
三相电机转速仿真结果如图3所示,三相电极实际电流值仿真结果如图4所示,由图3和图4可知,本发明方法电流控制更加稳定,在一定程度上避免了由于电流剧烈波动而引起的生产故障,从而确保了生产安全。
Claims (4)
1.一种带有输入时滞的电熔镁炉电极调节优化控制方法,其特征在于,包括以下步骤:
步骤1:根据三相电熔镁炉的系统状态变量、带有输入时滞的电机转速以及实际熔化电流的关系,建立三相电熔镁炉的带有输入时滞的电极调节系统线性模型和无输入时滞的电极调节系统线性模型;
步骤2:采用极小值原理求解三相电熔镁炉的无输入时滞的电极调节系统线性模型,得到无输入时滞的电极调节系统的黎卡提矩阵,从而确定无输入时滞的电熔镁炉电极调节最优控制率;
步骤3:根据最优控制理论,确定三相电熔镁炉的无输入时滞的电熔镁炉电极调节最优控制器;
步骤4:根据三相电熔镁炉的带有输入时滞的电极调节系统线性模型的二次性能指标和其对应的无输入时滞的电极调节系统线性模型的二次性能指标的数值相等,得到三相电熔镁炉的带有输入时滞的电熔镁炉电极调节最优控制器;
步骤5:在三相电熔镁炉熔化阶段工作过程中,将三相电熔镁炉设定熔化电流值与实际熔化电流值的偏差作为带有输入时滞的电熔镁炉电极调节最优控制器的输入,得到三相电机转速,控制电熔镁炉电机转子工作。
2.根据权利要求1所述的带有输入时滞的电熔镁炉电极调节优化控制方法,其特征在于,所述三相电熔镁炉的带有输入时滞的电极调节系统线性模型如下所示:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>B</mi>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>-</mo>
<mi>&tau;</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>C</mi>
<mi>x</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>;</mo>
</mrow>
其中,x(t)为带有输入时滞的三相电熔镁炉的系统状态变量,A为三相电机传动系数矩阵,B为带有输入时滞的执行机构传动系数矩阵,u(t-τ)为带有输入时滞的三相电机转速,y(t)为三相电熔镁炉实际熔化电流值,C为电弧系数矩阵,t为时间,τ为时滞;
所述三相电熔镁炉的无输入时滞的电极调节系统线性模型如下所示:
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mi>X</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>=</mo>
<mi>A</mi>
<mi>X</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>0</mn>
</msub>
<mi>u</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>y</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>=</mo>
<mi>C</mi>
<mi>X</mi>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>;</mo>
</mrow>
其中,X(t)为无输入时滞的三相电熔镁炉的系统状态变量,B0=e-AτB为无输入时滞的执行机构传动系数矩阵,u(t)为无输入时滞的三相电机转速。
3.根据权利要求1所述的带有输入时滞的电熔镁炉电极调节优化控制方法,其特征在于,所述无输入时滞的电熔镁炉电极调节最优控制器如下式所示:
<mrow>
<msubsup>
<mi>u</mi>
<mn>0</mn>
<mo>*</mo>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>K</mi>
<mn>0</mn>
</msub>
<mi>X</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>;</mo>
</mrow>
其中,为无输入时滞的三相电机最优转速,K0为无输入时滞的电熔镁炉电极调节最优控制率,X(t)为无输入时滞的三相电熔镁炉的系统状态变量,t为时间。
4.根据权利要求1所述的带有输入时滞的电熔镁炉电极调节优化控制方法,其特征在于,所述三相电熔镁炉的带有输入时滞的电熔镁炉电极调节最优控制器如下所示:
u*(t-τ)=-Kdx(t-τ);
其中,u*(t-τ)为带有输入时滞的三相电机最优转速,t为时间,τ为时滞,x(t-τ)为带有输入时滞的三相电熔镁炉的系统状态变量,K0为无输入时滞的电熔镁炉电极调节最优控制率,和为系统转移状态矩阵,A为三相电机传动系数矩阵,B为带有输入时滞的执行机构传动系数矩阵,K0为无输入时滞的电熔镁炉电极调节最优控制率,tf为终点时刻。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610596894.4A CN106017118B (zh) | 2016-07-27 | 2016-07-27 | 一种带有输入时滞的电熔镁炉电极调节优化控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610596894.4A CN106017118B (zh) | 2016-07-27 | 2016-07-27 | 一种带有输入时滞的电熔镁炉电极调节优化控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106017118A CN106017118A (zh) | 2016-10-12 |
CN106017118B true CN106017118B (zh) | 2017-12-29 |
Family
ID=57115446
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610596894.4A Active CN106017118B (zh) | 2016-07-27 | 2016-07-27 | 一种带有输入时滞的电熔镁炉电极调节优化控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106017118B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107526293B (zh) * | 2017-09-28 | 2020-08-28 | 东北大学 | 一种基于补偿信号的电熔镁炉电极电流切换pid控制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4815175B2 (ja) * | 2005-09-22 | 2011-11-16 | 東芝三菱電機産業システム株式会社 | アーク炉用電極昇降装置 |
CN103409639B (zh) * | 2013-07-31 | 2015-03-18 | 沈阳化工大学 | 电渣炉电极熔速最优控制诊断系统 |
-
2016
- 2016-07-27 CN CN201610596894.4A patent/CN106017118B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106017118A (zh) | 2016-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103741214B (zh) | 一种多晶硅铸锭工艺 | |
CN106643194B (zh) | 玄武岩电助熔窑炉温度控制系统与方法 | |
CN104164574A (zh) | 一种电渣重熔制备大型发电机护环用空心钢锭的方法 | |
CN202030861U (zh) | 一种多晶硅晶体生长炉加热装置 | |
CN102278893B (zh) | 一种用于高温炉温度控制系统的温度控制方法 | |
CN106017118B (zh) | 一种带有输入时滞的电熔镁炉电极调节优化控制方法 | |
CN105414509A (zh) | 一种电缆用无氧铜杆及其上引法制备工艺 | |
CN204058565U (zh) | 用于真空自耗电弧炉水套式结晶器的冷却装置 | |
CN103741213A (zh) | 一种多晶硅铸锭用熔料工艺 | |
CN103909238A (zh) | 一种自动化多坩埚组合式非晶制带机 | |
CN102409402A (zh) | 650kg多晶硅铸锭工艺 | |
CN202785966U (zh) | 一种控制平板玻璃成型温度场的装置 | |
CN105541108B (zh) | 基于3d打印技术的硫系玻璃元件制备方法 | |
CN202141955U (zh) | 基于plc的高温炉温度控制系统 | |
CN203397236U (zh) | 全电熔窑炉温度自动控制装置 | |
CN107526293B (zh) | 一种基于补偿信号的电熔镁炉电极电流切换pid控制方法 | |
CN203534209U (zh) | 极心圆可调的矿热炉电极装置 | |
CN103898603A (zh) | 一种双电源多晶硅铸锭工艺 | |
CN204494798U (zh) | 一种盐浴炉及其电控系统 | |
CN103590109A (zh) | 直拉单晶炉磁场装置及使用该磁场装置的拉晶方法 | |
CN202975874U (zh) | 应用于蓄电池铸焊模具温度plc及温控表自动控制装置 | |
CN102503357A (zh) | 高纯а-氧化铝的一种预熔化结晶方法 | |
CN109014107A (zh) | 大规格高碳铬轴承钢缩孔缺陷的控制方法 | |
CN208791814U (zh) | 一种氧化锆晶体生产装置 | |
CN202975876U (zh) | 应用于蓄电池铸焊模具温度plc全自动控制装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |