A Study on the Electrode Protection System of an 100t Ultra High
Power Electric Furnace
UnChol Choe, YongChol Sin, CholGuk Yun

Kim Il Sung University,
Pyongyang,
Democratic People’s Republic of Korea
Abstract
In this paper a method to design the electrode protection system using arc current and voltage and
pressure inside the hydraulic cylinder of electrode lift device of an 100t ultra high power furnace
(UHP furnace) and confident factor(CF) and Lyapnov function, and then we validated our method
through the field experimentation.
Key words: UHP furnace, electrode protection, electrode lift control, electrode broking, confident
factor, Lyapnov function

1. Introduction
One of main instruments making steel, an UHP furnace has the low productive cost and very high
productivity.
The electrode protection system to prevent the non-normal conditions such as the electrode
broking which occurs during steel-melting is very important for the safe operation and saving
electrode.
The previous works considered about control of melting in UHP are classified into approaches
using the mathematical modeling of plant and approaches using soft-computing.
In the work [4] they proposed the electrode lift control method of the mode settling the arc
current in which the on-off control is combined with the PID control, in the work [5] they proposed
the parameter identification method and the advanced electrode lift PID control algorithm using
least square method and pole-placement adaptive algorithm.
In the work [6] they proposed a reference model adaptive control using partly linear model
instead of non-linear model and in the works[1,2,8,12] they framed the intelligent control system
using soft-computing technology and then considered the electrode controlling issues.
In the previous works above, the core control value is arc current or voltage or impedance or arc
length.
The necessary problem for electrode control and establishing the protective system of electrode
is a model representing dynamics according to control operation [7~10].
The design theory and method about hybrid intelligent system aimed to multi-functions allow to
stand the defaults of individual computing and to quantify the complicated events exactly [13,14].
In the work[15], they introduced the concept of confident factor into uncertainty inference and
assessed the inference results with the value of confident factor .
Lyapnov function used in somedeciplines such as automatic control theory reflects the dynamic
1

changes quantitatively to assessed stability of a dynamic system.
They didn't study to use confident factor and Lyapnov function together in predictive control
or predictive control used soft computing.
That’s why in the classical predictive control they don’t need uncertainty inference because they
use its mathematical model and the research on model for inference dynamic change viewed from
change of energy is insufficient in the predictive control using soft computing.
In this paper we use arc current and the 2nd voltage and confident factor and Lyapnov function
being used in inference to assess the change of melting process of 100t UHP and proposed the
protective system for preventing electrode broking and validated through the introducing
experimentation.

2. The principle of electrode protection during steel-melting
The principle diagram of the electrode lift control system in UHP is as figure 1.

Figure 1. The principle diagram of the electrode lift control system in UHP

The electrode lift automatic control system is used to control the electrode lift automatically to
settle arc length constantly in order to settle arc current and arc voltage and so on about their set
points.
In this paper, arc current and the 2nd phase voltage and deviation of arc current and pressure
inside hydraulic cylinder for electrode lift are selected as factors representing dynamic change of
UHP.
And we spill the electrode lift control range into arc-generating range and normal range about
set current and over-current range where arc current is over set current to be able to interrupt the
normal operation and the dangerous current range.
From factors representing dynamic change of UHP and the electrode lift control ranges, the
protection system is framed as follow.
① A protection loop to prevent accident when arc current is dangerous.
② A protection loop to prevent accident when arc current is more than over current.
③ A loop to identify non-metal by pressure change inside the hydraulic cylinder for electrode lift.
④ A Loop to identify charge collapse when voltage is low.
2

⑤ Predicting charge collapse using intelligent method.

3. A protection loop by means of arc current and pressure inside
hydraulic device for electrode lift.
When arc current is over the dangerous current, maxim load current of furnace transformer, this
causes an electrical accident in electrode lift device and power system.
And the frequent operation at lower current than dangerous current causes the serious oscillating
of electrode that an electrode accident may occur.
The 1st and 2nd protection loops are realized by setting of the dangerous current and over current
and rules setting the out of electrode lift controller.
The setting of over current and dangerous current and control out for electrode lift is done
considering the maxim load current of furnace transformer whose capacity is 55,000KVA.
The set point of arc current in electrode lift control is from 35 000A to 45 000A.
I s = 35,000 A ~ 45,000 A .
When arc current is more than dangerous current, we must lift the electrode fastest,
I a ≥ I danger = ( I s + 15,000 A) ⇒ to lift the electrode fastest.

where I a is arc current and I danger is the dangerous current.
When the arc current is over current, we lift electrode at half of maximum velocity,
I a ≥ I over = ( I s + 7,500 A) ⇒ to lift electrode at half of maximum velocity,
where I over is bias of over-current.
If arc current suspends about over current more than 1s, we lift electrode fastest.
The 3rd protection loop prevents the electrode broking when electrode meets non-metal.
When electrode presses the non-metal during electrode control, electrode is forced physically and
at this time hydraulic pressure inside hydraulic cylinder lowers.
Therefore if the pressure inside electrode lift cylinder is lower than 0.9 times of normal pressure
for 1.5s, we agree that electrode meets non-metal.

4. Identification of charge collapse and electrode protection by
prediction.
The 4th and 5th protection loops prevent the electrode broking by physical and electrical impact
when charge collapses into the hole been formed during melting stage.

1) Identification of charge collapse by low-voltage and electrode lift control.
Phase voltage E 2 is the vector sum of voltage at reactance U X and voltage at resistance U R .
E2 = U X + U R
(1)
U X = Ia X
(2)
3

(

U R = I a Ra + R f

)

(3)

Here Ra , R f are arc resistances and system resistances respectively.
If measured the 2nd phase voltage be close to U X , Ra decreases and it means that electrode met
charge, namely the charge collapses.
We agree that the charge collapse when the following condition folds.
E 2 ≤ 1.2 ⋅ I a X ⇒ charge collapses
Then we must lift electrode fastest for 2s.

2) Prediction of charge collapse using intelligent method and electrode
protection.
Fuzzy rule R j whose conclusion part is a constant is of the form,
r j : if ( x1 is A j1 ) and L and ( x k is A jk )and L and ( x n is A jn )
then ( y j1 is a j1 ) and L and ( y jk is a jk ) and L and ( y jp is a jp )

(5)

where xk (k = 1, n) is the input variable of primary part of j ( j = 1, m) th fuzzy rule , A jk is the

membership function for xk (k = 1, n) in j ( j = 1, m) th fuzzy rule, y jk (k = 1, p ) is the output
variable of the primary part in j ( j = 1, m) th fuzzy rule, a jk is the constant for y jk (k = 1, p ) in
j ( j = 1, m) th fuzzy rule.

If we define the fuzzy membership functions as Gaussian function, it is possible to be correspond
the j ( j = 1, m) th rule to a RBF neural network.
In this paper we represent the dynamic change as the quantified assessment rule using the
confident factor and the change of Lyapnov function and consider a problem making a predictive
model with RBF NN network reflecting assessment rule.
(1). Making assessment rules using confident factor and Lyapnov function.
We define the language set assessing the dynamic change of plant at a moment as {to change
well, to change badly, to not change}.
In this paper we introduced the confident factor being used in inference uncertainty and
Lyapnov function representing varying energy of dynamic characteristics.

We select error of arc current from set current, e ∈ R 3×1 and error deviation e& ∈ R 3×1 at 2
moments previous and now.
We selected the deviation of Lyapnov function to be predicted at t moment, V&t and the confident
4

factor about dynamic change, CFt .
At this time we made rule for assessment about dynamic change of plant as follows,
rj : if x11 = μ11 and L x1n = μ1n and x21 = μ21 L and x2n = μ2n
then y j1 = a1 and y j 2 = a2 and y j 3 = a3

(6)

where r j is the j th rule, x ki is a variable, μ ki is the property function value according x ki , a k
is the property function value about the conclusion y jk , k = 1,3 to represent the predictive state about
dynamic change in rules.
The range of index is (k = 1,2, i = 1, n) and at the predict moment a1 denotes bigness of
orientation change of dynamics, a 2 denotes the confident factor of orientation change, a3
denotes the change orientation at next moment.
We selected Lyapnov function as follows,
V = X 1T PX 1
P>0

(7)

where X 1 = ( x11 ,L x1n ) T = (e(t − 1),Le(t − n))T , X 1 ∈ R n×1 .
We denote the variable representing the Lyapnov function at last moment as y 0 = a0 and
predicted about change of dynamic characteristics as follows.
Table2. qualitative results of change orientation of dynamic characteristics by qualitative
inference
Sign combination of

a 0 , a1

change orientation of dynamic characteristics

a0 ≤ 0, a1 < 0

To change to fine orientation

a0 ≥ 0, a1 > 0

To change to bad orientation

a 0 = 0, a1 = 0

No change

a 0 ≤ 0, a1 > 0

To oscillate and change to bad orientation

a 0 ≥ 0, a1 < 0

To oscillate to change to good orientation

The 2nd variable of conclusion part in the rule allows to predict the change of dynamic
characteristics with confident factor CF .
Therefore we make the 6th assessment rule as follows,
R1 : if x11 = μ11 L and x1n = μ1n and x21 = μ 21 L

L and x2 n = μ 2 n then y j1 = a1 and y j 2 = a2 and y j 3 = a3

5

(8)

R2 : if x11 = μ11

L and x1n = μ1n and x21 = μ 21 L

L and x2 n = μ 2 n then y j1 = a1 and y j 2 = a2 and y j 3 = a3
R3 : if x11 = μ11 L and x1n = μ1n and x21 = μ 21 L
L and x2 n = μ 2 n then y j1 = a1 and y j 2 = a2 and y j 3 = a3

(9)
(10)

where Ri , i = 1,3 denotes the set of rules formed with several rules.
V&
We represent a0 =
and select the property functions of conditional variables of the jth rules
V&
max

in the rule set Ri , i = 1,3 as follows,

μ1i = exp[−

( x1i − θ j1 ) 2

μ 2i = exp[−

2σ 2j1

( x2i − θ j 2 ) 2
2σ 2j 2

(11)

]

(12)

]

where θ ji denotes the mean and σ ji denotes the band of Gaussian function.
The property function value of each variable in conclusion part of rule is as follows.
a1 =

V&t
V&

(13)

max

a2 = CF =

μ~1 =

1
⋅ ( a1 + μ~1 + μ~2 )
3
1
n

n

∑

μ1i , μ~2 =

i =1

(14)
1
n

n

∑μ

2i

i =1

The value of a 3 representing the change orientation of dynamic characteristics is as follows.
if a0 ≤ 0, a1 < 0 , a 3 =1. if a0 ≥ 0, a1 > 0 , a 3 =-1
if a 0 = 0, a1 = 0 , a 3 =0. if a 0 ≤ 0, a1 > 0 , a 3 =-0.5
if a 0 ≥ 0, a1 < 0 , a 3 =0.5
Equations above allow to predict the change of dynamic characteristics with the valued value in
“how much” according to CF.
Because we must to analyze lots of data, we manipulate the issue of making the predictive model
by RBF NN network reflected the assessment rules.
(2) Construction of RBF NN network reflected assessment rules.
x ki

(k = 1,2, i = 1, n) ,the variables of primary part are selected as input variables of NN

network.
And a1 and a 2 , a 3 , the property function values of the out of conclusion part y i , i = 1,3

are

selected as outs of NN network.
The fuzzy rules are expressed as equations (8)~(10) and the membership functions of primary
6

part are expressed as equations (11)~(12) and then the firing strength of j , ( j = 1, M ) th rule is as
follows,
2

w j = ∏ w jk

(15)

k =1

n

w j1 = ∏ exp[−

( x1i − θ j1 ) 2
2σ 2j1

i =1

n

], w j 2 = ∏ exp[−

( x2i − θ j 2 ) 2

i =1

2σ 2j 2

]

where the inference result at the j th rule of rule set with M rules is as follows,
M

∑w
yk =

j

⋅y jk

j =1

=

M

∑w

M

∑ w~

j

y jk , (k = 1,3)

(16)

j =1

j

j =1

~ = wj .
where w
j
M
∑ wj
j =1

The structure of RBF NN network combined with the valued rules is as follows.

Figure2. The structure of RBF NN network combined with the valued rules

The error of out of the NN network from target about Nth training sample is as follows,
En =
M

where y k =

∑w
j =1

j

⋅y jk

M

∑w
j =1

1
2

3

∑ (T

k

− yk ) 2

(17)

k =1

M

~ y and M is number of cells in the hidden layer.
=∑w
j jk
j =1

j

We update the parameters in the form,
U new
= U old
+ ηΔU j
j
j

(18)

where η is the training rate, ΔU i is the modification of parameter.
The modification of parameter are as follows.
7

~ ⋅ (a − y ) ⋅ ( x j − θ ji )
Δθ ji = −(Tk − y k ) ⋅ w
j
k
k
2

(19)

~ ⋅ (a − y ) ⋅ ( x j − θ ji )
Δσ ji = −(Tk − y k ) ⋅ w
j
k
k
3

(20)

~ ⋅
Δa jk = −(Tk − y k ) ⋅ w
j

(21)

σ ji

σ ji

So we use the calculation results of the predictive model with RBF NN network to protect
electrode.

5. Introducing experiments
1) Simulation results by MATLAB and analysis
We show a curve predicting the charge collapse by CF which is out of the NN network as
figure 3.

Figure 3. a curve predicting charge collapse

In this figure, the points whose CF equals 1 predict the charge collapse and the circles are the real
charge collapse. There are charge collapses after 13 minutes and big charge collapse after 13
minutes 40 seconds and the accuracy of the predicting collapse is about 90%.

2) Introductions and effectiveness
We use the fuzzy self-tuning PI control in the electrode lift automatic control.
From the simulation above we made the identification loop of charge collapse and rules setting
control out.
If the change rate of variance of arc current increase in the same sign, a1 is big and the CF of
orientation change a 2 is over 0.6,we see the change orientation at the next moment a3 = -1 and
alert the charge collapse and lift the electrode which changes fastest at 30% of maxim lift velocity.
[Step 1]: to save the derivations of arc currents of 5 moments previous from now.
[Step 2]: if each derivation above is bigger than 300A, go to Step 3, else go to Step 5.
[Step 3]: summing the arc current derivations at 2 moments previous from now and name it Sum1,
and summing the arc current derivations from 3 moments previous to 5 moments previous and
name it sum2.
[Step 4]: if Sum1-Sum2>6000, alert the charge collapse and lift an electrode which changes
8

fastest at 30% of the maximum velocity.
[Step 5]: end
In 100t UHP furnace we examine 218 charges, among them we examined 98 charges with the
Fuzzy self-tuning PI control using our electrode protection system and we examined 120 charges
with the old method.
As result we shared 2.7kg/t of electrode per 1t of steel.

Conclusion
In this paper we studied to use arc current and the 2nd phase voltage and the pressure inside
hydraulic system for electrode lift and the confident factor and Lyapnov function assessing the
change of dynamic characteristics of 100t UHP and proposed the electrode protection system in
order to prevent the electrode broking.
We also validated our method through the field experimentation in 100t UHP.
References

[1] Thomas Brehm, Kuldip S. Rattan, “Hybrid Fuzzy Logic PID Controller”, IEEE Con
trol Systems Magazine, 807─814, 2005
[2] Baogang Hu, George K. I, “New Methodology for Analytical and Optimal Design of
Fuzzy PID Controllers”, IEEE Transaction on Fuzzy Systems, VOL. 7, NO. 5, October 20
02, 521─540
[3] Y. Oysal, “Fuzzy PID Controller Design For Load Frequency Control Using Gain Sc
aling Techique”, IEEE Trans, Systems Man Cybern, vol 20, 419─435, 2001
[4] 杨文华, 交流电弧炉电极调节控制系统. 重型机械,1994.4:31─49

[5]. 杨文华. 电

弧炉电极调节参考模型自造应系统的研究, 重型机械,1994.4:31─49

[6] 刘小河，刘康平，电弧炉电极调节系统自制的仿真研究，电工技术学报，1999
[7] T. Zheng and E.B. Makram, “An adaptive arc furnace model”, IEEE Transactions on
Power Delivery, Vol. 15, No. 3, July 2000, pp 931-939.
[8] A. M. O. Haruni, M. Negnevitsky, M. E. Haque, Kashem M. Muttaqi,
“Implementation of Artificial Intelligence Technique to Model Arc Furnace
Responses”, 2008 Australasian Universities Power Engineering Conference (AUPEC'08), p

9

p 205[9] Vito Logar, Dejan Dovzˇan and Igor Sˇkrjanc, “Mathematical modeling and experim
ental validation of an electric arc furnace”, ISIJ International, Vol. 51 (2011), No. 3, pp. 38
2–391
[10] Labar Hocine, Dgeghader Yacine, Kelaiaia Mounia Samira, and Bounaya Kamel,
“Closely Parametrical Model for an Electrical Arc Furnace”, World Academy of Science, En
gineering and Technology 16 2008, pp96-100. [11] G. Bisio et.al, “Heat transfer, energy sav
ing and pollution control in UHP electric-arc furnaces”, Energy 25 (2000), pp 1047–1066.
[12] G. Shabib, K. Hassan, “Fuzzy Logic Control of Three Phase Submerged Arc Ferro
silicon Furnace”, Proceedings of the 14th International Middle East Power Systems Conferen
ce (MEPCON’10), Cairo University, Egypt, December 19-21, 2010, Paper ID 156.
[13] Orest Lozynskyi, Roman Paranchuk, Yaroslav Paranchuk, “Multicriterion intelligent
control system and optimal stabilization of arc steel-melting furnace electrical regimes coor
dinates”, Computational Problems of Electrical Engineering, Vol. 1, No. 1, 2011, pp 35-44.
[14] Z. Zhang, C. Zhang, “Agent-Based Hybrid Intelligent Systems(An Agent-Based Framework for
Complex Problem Solving)”, Springer, 2005, pp 15-33
[15] Adrian A. Hopgood, “Intelligent Systems for Engineers and Scientists”, CRC Press, 2001, pp 20
~42

10

