Proceedings ECT 2009
Proceedings ECT 2009
Proceedings ECT 2009
PROCEEDINGS
OF THE 4 INTERNATIONAL CONFERENCE ON
ELECTRICAL AND CONTROL TECHNOLOGIES
TH
2009
KAUNAS
ISSN 1822-5934
PROCEEDINGS
OF THE 4 INTERNATIONAL CONFERENCE ON
ELECTRICAL AND CONTROL TECHNOLOGIES
TH
Sponsored by:
Lithuanian State Science and Studies Foundation
ABB UAB
Conference organizers
Kaunas University of Technology, Lithuania
IFAC Committee of National Lithuanian Organisation
Lithuanian Electricity Association, Full Member of EURELECTRIC
Conference Chairman
Prof. J. Daunoras, Dean of Kaunas University of Technology,
Electrical and Control Engineering faculty, Lithuania
Honorary Chairman
Prof. R. J. Kays, Vice Rector of Kaunas University of Technology, Lithuania
Conference Program Advisory Committee:
Prof. R. Deksnys, Kaunas University of Technology, Lithuania
Prof. L. Fara, Polytechnic University of Bucharest, Romania
Prof. J. Gerhards, Riga Technical University, Latvia
Prof. E. Heinemann, Schmalkalden University, Germany
Dr. B. Jaekel, Siemens AG, Germany
Prof. T. Jokinen, Helsinki Technical University, Finland
Prof. V. Kaminskas, Vytautas Magnus University, Lithuania
Prof. S. Kauinis , Kaunas University of Technology, Lithuania
Prof. J. Laugis, Tallinn University of Technology, Estonia
Assoc. Prof. V. Maerauskas, Kaunas University of Technology, Lithuania
Prof. L. A. Markeviius, Kaunas University of Technology, Lithuania
Prof. V. Mikinis, Lithuanian Energy Institute, Lithuania
Prof. A. Morkvnas, Kaunas University of Technology, Lithuania
Prof. A. Narglas, Kaunas University of Technology, Lithuania
Prof. A. Nemura, Lithuanian Energy Institute, Lithuania
V. Pakeviius, President of Lithuanian Electricity Association, Lithuania
Prof. L. Ribickis, Riga Technical Univesity, Latvia
Assoc. Prof. J. Sa da Costa, Technical University of Lisbon, Portugal
Prof. A. Sauhats, Riga Technical University, Latvia
Prof. R. Simutis, Kaunas University of Technology, Lithuania
Prof. A. W. Sowa, Bialystok Technical University, Poland
Prof. R. Strzelecki, Zielona Gora University, Poland
Prof. M. Valdma, Tallinn University of Technology, Estonia
Conference Editorial Committee:
General Editor: Prof. A. Navickas, Kaunas University of Technology, Lithuania
Prof. A. Sauhats, Riga Technical University, Latvia
Prof. A. Virbalis, Kaunas University of Technology, Lithuania
Assoc. Prof. V. Aubalis, Kaunas University of Technology, Lithuania
Assoc. Prof. V. Galvanauskas, Kaunas University of Technology, Lithuania
Assoc. Prof. A. Jonaitis, Kaunas University of Technology, Lithuania
Organizing committee:
Assoc. Prof. M. Aubalis, Kaunas University of Technology, Lithuania
Assoc. Prof. V. Galvanauskas, Kaunas University of Technology, Lithuania
Assoc. Prof. A. Jonaitis, Kaunas University of Technology, Lithuania
CONTENTS
A new approach based on support vector machine to analyse and control the REM-sleep-stage
Maria Krautwald, Rudolf Baumgart-Schmitt, Daniel Trommer, Christian Walther
11
Adaption of fuzzy rules by different multi criteria optimization procedures to control stages of anaesthesia
Christian Walther, Rudolf Baumgart-Schmitt, Maria Krautwald, Daniel Trommer
15
19
23
27
31
36
40
44
48
54
58
62
66
70
75
79
A new wireless solution to transmit and evaluate biomedical signals under low signal noise ratio to control
the state of relaxation
Daniel Trommer, Markus Greulich, Rudolf Baumgart-Schmitt, Maria Krautwald
Christian Walther, Andreas Wenzel
83
87
Kazimierz Puaski Technical University of Radom, Faculty of Transport and Electrical Engineering, Poland
93
99
Kaunas University of Technology, Lithuania; Halmstad University, Sweden; Kaunas University of Medicine, Lithuania
103
107
111
115
Mathematical model for adaptive control of dissolved oxygen concentration in biological wastewater
treatment process
Kstutis Jonelis, Kstutis Brazauskas, Donatas Leviauskas
121
Adaptive control system for dissolved oxygen concentration in Kaunas wastewater treatment plant
Jolanta Repyt, Ieva Cebatorien, Andr Wenzel
125
131
Different simulation models in exemplary employment for pump aggregate pipeline systems
Anna Koziorowska, Jacek Bartman
University of Rzeszw, Institute of Technology, Poland
135
139
143
146
152
156
162
166
Analysis of Latvian power system development in the Kurzeme region using multicriterion optimization
Anatoly Mahnitko, Alexander Gavrilov
170
176
182
186
192
196
200
205
211
215
219
The edge effects influence on the braking force of linear induction motor
Bronius Karalinas, Edvardas Matkeviius
223
227
231
235
239
243
247
A comparative study between small traditional induction motors and small conducting sleeve induction
motors using 3D finite element transient analysis
Seyed Mohsen Hosseini, Javad S. Moghani, Akram Hosseini
251
Graphical method of estimating optimum values for the main dimensions of a D.C. Machine
C. Easwarlal, V. Palanisamy, M.Y. Sanavullah
256
The research of causes and consequences of parameters changes of grounding electrical installations
Marina Kizhlo, Arvids Kanbergs
Riga Technical University, Institute of Power Engineering, Latvia
261
A simple transient-based wide-area selectivity technique for earth-faults in isolated and compensated
neutral medium voltage networks
Mohamed F. Abdel-Fattah, Matti Lehtonen
266
Sensors for on-line partial discharge detection in covered-conductor overhead distribution networks
G. Murtaza Hashmi, Matti Lehtonen, Mikael Nordman, Muzamir Esa
272
Power Systems and High Voltage Engineering, Helsinki University of Technology (TKK), Espoo, Finland
Ensto Utility Networks, Porvoo, Finland
278
Influence of the distribution of the axial, radial and circumferential component of the magnetic induction
vector within the space between switches on rupturing capacity of vacuum switches
Bogdan Kwiatkowski, Robert Pkala
282
Rzeszow University , Technology and IT Institute; Rzeszow University, Technology and IT Institute, Poland
The research of the influence of the boundary conditions of the atmospheric overvoltages on the insulation
of high voltage equipment
Saulius Gudius, Linas Andronis Markeviius, Alfonsas Morkvnas
286
290
The models of sources of electromagnetic interferences over UHV line for transient processes investigation
Linas A. Markeviius, Vytautas ioinys
294
298
302
308
312
318
323
327
Potential of non-traditional energy sources for cold recovery in absorption refrigerating machines
Michail G. Berengarten, Semyon I. Vainshtein, Aleksandr S. Pushnov, Vytautas Adomaviius
333
337
Application of renewable energy sources for feeding of water recycling systems technological equipment
Michail G. Berengarten, Semyon I. Vainshtein, Aleksandr S. Pushnov, Vytautas Adomaviius
343
Wind power plant with redox flow battery inverter system converter energy conversion processes
research
Povilas Balinas, Povilas Norkeviius
347
Hybrid water heating system based on solar collectors and heat exchanger fed from the centralised
heat supply grid
Vytautas Adomaviius, Gytis Petrauskas
352
Optimisation of hybrid water heating system based on solar collectors and heat exchanger in order to maximize
356
solar heat energys share
Vytautas Adomaviius, Gytis Petrauskas
The Centre of Renewable Energy Technologies at KTU, Lithuania
Kaunas University of Technology, Lithuania
Control of energy fluxes in the grid-tied domestic power system based on wind energy
eslovas Ramonas, Vytautas Adomaviius, Vytautas Kepalas
362
Research of the power conversion processes in the system of power supply from a number of wind
turbines over the one grid-tied inverter
eslovas Ramonas, Vytautas Adomaviius, Vytautas Kepalas
368
374
The Centre for Renewable Energy Technologies at Kaunas University of Technology, Lithuania
Lithuanian-German joint venture UAB Tometa, Germany
Joint Stock Company UAB Saules energija, Lithuania
Kaunas Business Machines Company, Lithuania
10
380
1. Introduction
2. Methods
11
Figure 1 shows a two dimensional room with a trainingdataset {xi,yi} illustrate by a set of white circles (class
1) and a set of black circles (class 2).
w
b
d H , x i yi
, xi
,
w
w
d+
d-
H+
w, w .
(6)
d 1 w, x b
H-
class 1 (y = +1)
support vector
d 1 w, x b
M d d
(1)
H 2 x | w, x b 1.
(4)
(9)
w, x i b 1 if yi 1
(10)
w, x i b 1 if yi 1 .
(11)
w
2
.
, x x
w
w
and
(2)
H 1 x | w, x b 1
w
b
1 , (8)
, x
w
w
w
(7)
H x | w, x b 0,
b
w
1
, x
w
w
w
and
class 2 (y = -1)
and
(5)
n
1 2
w Ci ,
2
i 1
12
(12)
Ld ( ) i
w, x i b 1 i if yi 1
(13)
w, x i b 1 i if yi 1 .
(14)
i 1
1 n
i j yi y j x i , x j ,
2 i , j 1
(19)
and
i 0
and
y
i 1
0.
(20, 21)
class 1 (y = +1)
(22)
support vector
class 2 (y = -1)
L ( x, ) f ( x ) h ( x ) ,
(15)
db
(17)
K ( x i , x j ) exp x i x j
i 1
and
dL p
dw
0 w i yi x i .
(23)
0 i yi 0
dL p
class 2 (y = -1)
1 2 n
w i yi w, x i b 1 , (16)
2
i 1
w ,b
class 1 (y = +1)
class 2 (y = -1)
L p ( w, b, )
class 1 (y = +1)
(18)
(24)
0.
i 1
13
4. Conclusion
The results of the SVM have to observe critically
because of a possible preference of one class.
The SVM is suitable for a two class classification. Even
for classification of biological signals it can be used but
it doesnt achieve excellent percent values.
The suitability of the SVM concerning two-class
classification of biological signals was also shown
before with anaesthetic classification between awake not awake and stage A4 - not A4. The results were over
80% by using different datasets [4].
In multiclass-classification a SVM show similar results
compared to Fuzzy-Logic and neuronal networks [5].
Every learnt SVM only need one prediction algorithm
and appending upon Training-data different Input data.
That is the reason why a SVM could be very efficient
on a microcontroller platform.
5. References
1. Rasch B., Pommer J., Diekelmann S. & Born J.
Pharmacological REM sleep suppression
paradoxically improves rather than impairs skill
memory Nature Neuroscience advance online
publication Published online: 5 October 2008
doi:10.1038/nn.2206
2. Suykens J.A.K., Van Gestel T., De Brabanter J., De
Moor B. & Vandewalle J. Least squares support
vector machines World Scientific Publishing Co.
Pte. Ltd: 2002 ISBN 9812381511
3. Hsu C.-W., Chang C.-C. and Lin C.-J. A Practical
Guid to Support Vector Classification. Dec 20 2007.
www.csie.ntu.edu.tw/~cjlin
4. Walther C., Baumgart-Schmitt R. and Backhaus K.
Support Vector Machines and Optimized Neural
Networks Adaptive Tools for Monitoring and
Controlling the Depth of Anaesthesia. Proceedings
of International Conference Electrical and Control
Technologies 2008, Kaunas, 2008. p. 11-14.
5. Baumgart-Schmitt R., Walther C. and Backhaus K.
Design and Assessment of Fuzzy Rules by Multi
Criteria Optimization to Classifiy Anaesthetic
Stages. 4^th European Congress of the International
Federation for Medical and Biological Engineering,
Antwerp, Belgium, 2008.
3. Results
Because of the individualism of every human it is
difficult to classify EEG-signals on the basis of training
data from another patient.
In consideration of this the result of 83.07% correct
classification is satisfying.
The results have to observe critically. Often one class
was preferred. And even if the SVM can classify to
90% it is disuse because it only classified Non-REM
which was present to 90%. There are two possibilities
to avoid this situation. First is to test the SVM with
Non-REM and REM in a data record ratio 1:1. Another
is to adjust prediction of the SVM that both classes are
predicted percent similar.
The results were analyzed by calculating how much
percent of each class was classified correctly. Like it
shown in Table 1 class Non-REM were classified by
83.84 % and REM by 79.26 %. These values are nearly
similar. Then the SVM is suitable for online
classification.
Table 1. Results of classification by the Support Vector
Machine
testing-data
number of epochs
number of epochs in %
right classify epochs
right classify epochs in %
Non-REM
10793
83.23
9049
83.84
REM
2175
16.77
1724
79.26
14
15
Data Set
Gender
Age
Epochs
0066
0007
0006
0010
0011
0009
0072
0012
0021
0024
F
F
F
F
F
F
F
F
F
F
52
62
63
45
66
36
75
37
42
41
1180
1438
449
1286
947
355
1828
614
1068
995
16
Data Set
Gender
Age
Number of
Epochs
0121
0118
0117
0116
0125
0127
0126
0134
0135
0128
F
F
M
M
F
F
F
F
M
M
73
77
47
65
84
24
45
45
67
55
1745
385
546
371
399
440
283
299
467
1389
(1)
Criteria
q1
q2
q3
q4
Table
Weight
mode Aw
Weight
mode Bw
0.4
0.3
0.3
0.0
4,6,7,8
0.2
0.4
0.4
0.0
W (E)
A1
(E)
A2
(E)
A3
(E)
A4
(E)
3. Results
To get the final results, which are represented in the
tables 4, 5, 6 and 7 a sub set of fuzzy models was
extracted from the Pareto-optimal front for every
MOEA-Run. At the end of a multi-criteria optimization
process a suitable number of solutions have to be
selected from the Pareto set. In our case a general
selection criterion was defined by the sum of the
weighted criteria q1 to q4 for every individual of the
Pareto-optimal front. If this selection criterion has a
high value the evaluated individual implies a high
performance within the search space. The used weight
modes Aw and Bw with different weight vectors are
presented in table 3.
After the selection of the individuals, these different
optimized fuzzy models were used to classify the
validation data (table 2). To explain the problem of
finding the best solutions using the two different weight
modes Aw and Bw the validation results of the respective
individuals are presented in table 4 and 5. Both sets of
individuals were trained in the optimization process
using the NSGAII algorithm. The fuzzy models used to
calculate the solutions presented in table 5 show a
slightly better performance in comparison to the
solutions of table 4. The performance is measured by
the degrees of concordance between the automatic
classification in each column and the given ratings of
the anaesthesiologists in each row. Beside the diagonal
fields there are also some grey highlighted rectangles in
the tables 4, 5, 6 and 7. This means that confusions
indicated by these values inside the rectangles are
negligible. The fuzzy models which have performed the
results of table 5 were only found by using the weight
mode Bw.
W
19.8
3.2
A1
47.8
43.6
A2
9.9
16.5
A3
2.7
21.4
A4
19.8
15.3
0.1
4.9
32.7
60.4
1.7
0.0
0.3
14.3
82.4
3.0
0.0
0.8
5.7
31.3
62.1
mode Bw
NSGAII
W (E)
A1
(E)
A2
(E)
A3
(E)
A4
(E)
W
21.6
5.2
A1
44.1
42.3
A2
8.1
17.4
A3
1.8
13.3
A4
24.3
21.8
0.1
5.4
35.4
57.7
1.5
0.0
0.3
16.7
78.5
1.5
0.0
1.4
4.3
27.4
66.9
W (E)
A1
(E)
A2
(E)
A3
(E)
A4
(E)
17
W
27.9
5.24
A1
43.2
44.8
A2
4.5
18.2
A3
7.2
16.1
A4
17.1
15.7
0.6
4.9
27.6
63.3
3.5
0.1
0.9
8.4
86.6
4.0
0.3
2.0
7.6
30.8
59.4
W (E)
A1
(E)
A2
(E)
A3
(E)
A4
(E)
5. Conclusions
W
63.1
10.5
A1
10.8
41.5
A2
0.0
17.3
A3
2.7
5.7
A4
23.4
25.0
0.0
4.4
43.1
46.5
6.0
0.0
1.8
20.6
70.9
6.8
0.1
6.0
11.9
17.2
64.8
6. References
1. Walther C., Baumgart-Schmitt R. and Backhaus K.
Support Vector Machines and Optimized Neural
Networks Adaptive Tools for Monitoring and
Controlling the Depth of Anaesthesia. Proceedings
of International Conference Electrical and Control
Technologies 2008, Kaunas, 2008. p. 11-14.
2. Baumgart-Schmitt R., Walther C., Backhaus K.,
Reichenbach R., Sturm K.-P. and Jger U. Robust
Nonlinear Adaptive Network Classification of
Anaesthesia. 2008 IAPR Workshop on Cognitive
Information Processing, Santorini, Greece, 2008.
3. Baumgart-Schmitt R., Walther C. and Backhaus K.
Design and Assessment of Fuzzy Rules by Multi
Criteria Optimization to Classifiy Anaesthetic
Stages. 4th European Congress of the International
Federation for Medical and Biological Engineering,
Antwerp, Belgium, 2008.
4. Hoffmann
T.
Erprobung
von
Methoden
mehrkriterieller Optimierung zur nichtlinearen
Klassifikation am Beispiel von simulierten
Evolutionsstrategien. Diploma thesis. University of
Applied Sciences Schmalkalden, Germany, 2009.
5. Kim M.-S., Kim C.-H. and Lee J.-J.
Structure/Parameter Optimization of Fuzzy Models
by Evolutionary Algorithm. Proceedings of the
IASTED International Conference on Intelligent
Systems and Control, Salzburg, Austria, 2003.
6. Laabidi K., Bouani F. and Ksouri M. Multi-criteria
optimization in nonlinear predictive control.
Mathematics and Computers in Simulation, Volume
76 , Issue 5-6, 2008. p. 363-374.
7. Deb K., Pratap A., Agarwal S. and Meyarivan T. A
Fast and Elitist Multiobjective Genetic Algorithm:
NSGA-II. IEEE Transactions on Evolutionary
Computation, Volume 6, Issue 2, 2002. p. 182-197.
8. Zitzler E., Laumanns, M. and Thiele L. SPEA2:
Improving the Strength Pareto Evolutionary
Algorithm. TIK-Report 103, Swiss Federal Institute
of Technology Zurich, Switzerland, 2001.
9. Zitzler E. and Knzli S. Indicator-Based Selection in
Multiobjective Search. Proceedings of the 8th
International Conference on Parallel Problem
Solving from Nature (PPSN VIII), Birmingham,
UK, 2004. p. 832-842.
W,A1,A2,A3 (E)
A4 (E)
W,A1,A2,A3
89.4
31.6
A4
10.6
68.4
4. Discussion
Despite using different multi criteria procedures the
results differ only slightly in classification of the
validation data.
Differences in the degrees of concordance obviously
only depend on the chosen individuals and the used
selection criteria, which were defined by a weight
mode. Mainly the three multi-objective evolutionary
algorithms differ in sorting the individuals and
calculating the density value for each fuzzy model. This
causes a different speed of the optimization process
itself for each algorithm. In every case a population of
optimized individuals is generated as a result.
Measurements of the optimization speed of such
algorithms (NSGAII, SPEAII and IBEA) were not
analyzed within this work. To estimate the efficiency of
the multi criteria optimization procedures defined test
functions were employed by coworker Hoffmann [4].
The tables 4, 5, 6 and 7 contain higher degrees of
confusions of the stage W and A1 to the stage A4.
These values are a result of the higher dimensionality of
the fuzzy models in the training procedure in
comparison to [3]. It is possible to distinguish correctly
between Wake, A1 and A4, which was shown in [1] and
[3]. This is underlined by the results shown in Table 8.
18
2. Theoretical background
General fuzzy inference is used in systems whose rule
structures are predetermined by the user and
membership functions are fixed and arbitrarily chosen.
ANFIS allows generating fuzzy inference system
automatically, without having predetermined structure
based on characteristics of variables in the system.
Membership function parameters are set during learning
process of the inference system. Adaptive network
based fuzzy inference systems (ANFIS) is a FIS
implemented in the framework of an adaptive fuzzy
neural network. Neuro-adaptive learning techniques
provide a method for the fuzzy modelling procedure to
learn information about a data set, in order to compute
the membership function parameters that best allow the
associated fuzzy inference system to track the given
input/output data. ANFIS constructs a fuzzy inference
19
RGB
images
of an eye
ANFIS learning
Segment
images?
No
Yes
Image
segmentation
(1)
z = ax + by + c
Selection of the
requisite segment
ANFIS
learning
data
selection
Noise reduction
Data for
analysis
Finish
4. Results of experiments
0
1
A and B
w z
i i
FinalOutput
i 1
N
(2)
i 1
3. Segmentation algorithm
An algorithm for image segmentation using different
trained fuzzy set structures is presented in Fig. 2. It
detects iris of a human eye. The color composition of
20
Fig. 3. The data base of the images with human eye positioned
in five different positions
Learnin
g epochs
Number
of MFs
MF shape
Training
time, s
Processin
g speed,
fps
1
2
3
20
Triangular
Trapezoidal
Gaussian
curve
Triangular
Trapezoidal
Gaussian
curve
Triangular
Trapezoidal
Gaussian
curve
Triangular
Trapezoidal
Gaussian
curve
Triangular
Trapezoidal
Gaussian
curve
Triangular
Trapezoidal
Gaussian
curve
1.49
1.44
8.14
8.27
1.53
10.78
10.83
7.28
5.09
5.26
10.96
3.33
3.51
4.65
8.17
8.26
3.72
26.71
26.96
7.23
5.10
5.32
27.66
6.79
6.71
4.74
8.23
8.35
7.49
54.05
54.14
7.30
5.02
5.19
55.50
4.60
4
5
6
7
8
9
50
10
11
12
13
14
15
16
17
18
21
100
2. Theoretical background
Background subtraction is widely used to detect moving
object from static cameras. It is usually regard as one of
the most important step in applications such as traffic
monitoring, human motion capture, recognition and
video surveillance, etc. In order to get a good
performance of the whole system, the background
subtraction method should not be so time and space
consuming, and the accuracy is also required. Gaussian
mixture model is a robust background subtraction
method and is widely used ever since it is proposed.
23
1
(2)v/ 2 detKi
(1)
1
exp( (c i )T Ki1(c i ))
2
Mahanalobis distance between new pictures pixel and ith region is calculated as follows in (2) expression:
M (c, i) (c i ) T K i1 (c i )
Fig. 2. Cursor control systems web-cam position: webcams capturing region horizontal to xy plane; 1 vertically
down directed web-cam; 2 users hand. Hand movement
directions are horizontal to plane
(2)
0 otherwise
(3)
Input
layer
Hidden
layer
Output
layer
Output
B
Fig. 1. Neural networks structure
24
1
2
100
50
20
40
60
80
100
120
140
160
180
200
Y coord. Recognised and real hand cordinates when picture size: 120x160
100
80
60
40
20
2
0
20
40
60
80
100
120
140
160
180
200
-5
20
40
60
80
100
120
140
160
180
200
180
200
20
40
60
80
100
120
140
160
25
6. Acknowledgment
0.2
0.1
0
160x120
320x240
640x480
7. References
1. Zhen Tang, Zhenjiang Miao. Fast Background
Subtraction and Shadow Elimination Using
Improved Gaussian Mixture Model. Haptic, Audio
and Visual Environments and Games. IEEE
International Workshop on 12-14 Oct. 2007. p. 38
41.
2. Shi Peijun, Jones E.G., Zhu Qiuming. Median
model for background subtraction in intelligent
transportation system. Proceedings of the SPIE,
Volume 5298, 2004. p. 168-176.
3. Alberto Broggi, Pietro Cerri, Paolo Medici, Pier
Paolo Porta. Real Time Road Signs Recognition.
Proceedings of the 2007 IEEE Intelligent Vehicles
Symposium Istanbul, Turkey, June 13-15, 2007. p.
981 986.
4. Shuangdong Zhu, Lanlan Liu. Traffic Sign
Recognition based on Color Standardization.
Proceedings of the 2006 IEEE International
Conference on Information Acquisition August 20 23, 2006. p. 951 955.
5. Yang Kun, Zhu Hong, Pan Ying-jie. Human Face
Detection Based on SOFM Neural Network.
Proceedings of the 2006 IEEE International
Conference on Information Acquisition August 20 23, 2006. p. 1253 1257.
6. Paliy I., Kurylyak Y., Kapura V., Sachenko A,
Lamovsky D, Sadykhov R. Combined Approach to
Face Detection for Video-Surveillance. IEEE
International Workshop on Intelligent Data
Acquisition and Advanced Computing Systems:
Technology and Applications 6-8 September, 2007.
p. 594 598.
7. Andy Chiem, Adel Al-Jumaily, Rami N. Khushaba.
A Novel Hybrid System for Skin Lesion Detection.
1-4244-1502-0/07, 2007. p 567 572.
160x120
320x240
640x480
t, s
0.005
0.184
0.665
2(x)
2(y)
1.241
6.797
59.911
0.667
2.393
16.185
X step
9
5
2
Y step
8
4
2
26
27
(1)
-1
Reinforcement learning can be classified as a nonsuper- wised method of machine learning. In our case
the agent (robot) tries to minimize a long term negative
reward (return) in goal directed moving through an
unknown environment. The on-line learning takes place
during the complete actions. The state specific decisions
are the results of (1).
3. Cellular Automata
Cellular automata were successfully applied in image
processing [10]. The optimization procedure bases on a
model which is created by sensor data and modified by
cellular automata. We solved three different image
processing problems by a cellular automat. The first
problem is connected with recognized gaps between the
hindrances which are too small for passing by the
selected robots. The rules of the defined automat are
responsible to close the gates before the optimization
procedure starts. The next problem arise if the selected
resolution of the two dimensional model of the world is
too small compared to space which is necessary for one
robot. The optimal path calculated by DDP can lead
beside the hindrances too narrow for passing by. The
cellular automat is used to make the walls virtually
6. Simulation
To test the reinforcement learning and the
implementation in C++ the world for path finding was
simulated by software. The simulated robots were
placed on one square only according the pre selected
spatial resolution. Therefore there was no need to apply
the rules of the cellular automat. By means of a drawing
program the model can be created on the selected
spatial resolution without any limitation. Up to 6 robots
can be started. Because of the high speed of this
simulation the different models of the world can be
tested and evaluated very efficiently. Fig. 2 presents the
output of the program in simulation mode. The path of
28
29
30
2. Robot E-puck
E-puck [4] is small and elegant robot. It is comfortable to
use in research because it is small, light and fast. Its view
from the top is showed in Figure 1.
The diameter of it is 70mm. The structure of this robot can
be defined by three parts the body, motors and printed
circuit. The body ensures other parts, including the
battery. The battery is connected with printed circuit by
two contacts. The battery is based on LiION technology,
has 5Wh capacity and is sufficient for about 2-3 hours of
intensive use.
The motor uses steppers with gear reduction. There are 20
steps per revolution, the gear reduction is 50:1. The Robot
31
Sensor values
3500
3000
2500
2000
1500
1000
500
0
0.5
1.5
2.5
3.5
4.5
5.5
Distance
(1)
32
global
, local [90;180]
.
, local (180;90]
d right dleft
cos( local ),
x x0
2
2
d d
local _ new local left right ,
d
(2)
where D d baseline
200
global angle
100
50
0
-50
-100
-150
-100
100
200
local angle
Analogously:
, global [90;180]
90 global
.
local
270
global
global ( 180;90]
y global x local
(3)
(9)
k k 0 90 local ,
(10)
(5)
local)
x x0 D sin(local) D sin(
d
d
d
y y D cos( ) D cos( right left )
0
local
local
d
; dbaseline is distance
(4)
local)
x x0 D sin(local) D sin(
d
d right dleft
y y D cos( ) D cos(
local)
0
local
d right d left
(8)
150
-200
-200
d right d left
(7)
(6)
33
4. Results
We get better results but there also is some noise and the
slide of the wheels. Perhaps robot must turn around slower
because it has problems then.
34
not enter the territory but wants get the map of that
territory
6. Acknowledgements
35
38
Fig. 8. Path of the third size robot in the maze. Robot finds
first too narrow opening and searches for new shortest path
Fig. 11. Path of the third size robot in the maze. Search of the
shortest path was successful. Circles in maze shows places,
where robot does not pass obstacle
7. Conclusions
Suggested method allows using of vector marks to find
nearly shortest path for a robot with finite dimensions as
well as in the case of a robot represented by a material
point. That makes the finding the path quite simple and
quick.
Too narrow openings can be detected and closed by
artificial obstacles to avoid in future. The tree of vector
marks is modified cutting not useful and adding new
necessary branches, when artificial obstacles added.
The algorithm finds the shortest possible way with
respect to too narrow openings in quite complicated
environments (mazes), as shown in examples.
Fig. 9. Path of the third size robot in maze. Robot finds
another (too narrow opening) obstacle and searches for new
shortest path
8. References
1. Baranauskas V., arkauskas K., Bartkeviius S.
Target Retrieval in Known Enviroment. Electronics
and Electrical Engineering, 2008, No. 8(88). p.9194.
2. Baranauskas V., Bartkeviius S., arkauskas K.
Modified Algorithm of Vector Marks Formation.
Proceedings of the Conference on Electrical and
Control Technologies, 2008, Kaunas, 2008. p. 1518.
3. Pimenta L.C.A., Fonseca A.R., Pereira G.A.S.,
Mesquita R.C., Silva E.J., Caminhas W.M., Campos
M.F.M. Robot Navigation based on Electrostatic
Field Computation. IEEE Transactions on
magnetics, Vol. 42, No.4, April 2006. p. 1459-1462.
Fig. 10. Path of the third size robot in maze. Robot finds new
possible shortest path, but another to narrow corridor is found
and robot must found new shortest path, according to vector
marks
39
40
41
cos 1
sin
1
T0,1
0
0
0
sin 1
cos 1
cos 2
sin
2
T1,2
0
sin 2
cos 2
0
0
cos 3
sin
3
T2,3
0
sin 3
cos 3
0
0
0
0
,
l1
1
l2 cos 2
l2 sin 2
,
l3 cos 3
l3 sin 3
.
(1)
(2)
(3)
multiplied:
T T0,1 T1,2 T2,3 .
(4)
0]
0]
-l1]
1]
T2 =
[ cos(t2), -sin(t2),
[ sin(t2), cos(t2),
[
0,
0,
[
0,
0,
Kinematic
pair type
Rotational
Link
No.
1
1, 2
Rotational
2, 3
Rotational
Parameters
d
a
0 90
l1
l2
l3
0,
0,
1,
0,
l2*cos(t2)]
l2*sin(t2)]
0]
1]
l3*cos(t3)]
l3*sin(t3)]
0]
1]
T=T1*T2*T3
Fig. 4. MatLab task for matrix multiplication
42
is
homogenous
C1C2 C3 C1S2S3
S C C S S S
1
2
3
T 1 2 3
S2 C3 C2S3
C1C2S3 C1S2 C3
S1C2S3 S1S2 C3
S1
C1
S2S3 C2 C3
(5)
x cos 1 cos 2l3 cos 3 cos 1 sin 2l3 sin 3 cos 1l2 cos 2 ;
y sin 1 cos 2l3 cos 3 sin 1 sin 2l3 sin 3 sin 1l2 cos 2 ;
z sin 2l3 cos 3 cos 2l3 sin 3 l2 sin 2 l1 .
7. Step forming using derived homogenous
transformation matrix
8. Conclusions
43
2. Algorithm overview
The main idea of the algorithm is that the decisions are
immediately evaluated for their success and a special
subsystem tries to adjust them so that they would
perform better in the next generation.
The robot obtains proximity sensor information and
passes it to an artificial neural network, saving a copy of
sensor values in an array. After this, both the sensor
information and the decision of the network are
forwarded to a primary IF-THEN test. The exact
content of this test is unimportant for our purposes,
44
4. Experiment setup
For our experiment the e-puck robot was connected to
MATLAB computing environment, which was then
used to implement the algorithm and obtain various
debugging information in an accessible graphical
representation. The connection was achieved via
Bluetooth technology, with all data sent to and read
from an emulated COM serial port.
The robot was put into a box approximately 60
centimetres long and 40 centimetres wide, with several
cardboard obstacles inside (Fig. 3). Both the inside
walls of the box and obstacles were covered with sheets
of white paper to enhance their reflectivity to infrared
waves. Obstacles were not fixed to the box in any way
to make changing the configuration of the test maze in
the middle of an experiment possible.
Fluorescent lights were turned off in the surrounding
area to reduce the noise picked by infrared sensors.
45
S
2
S
3
Lef
t
S
4
S
5
S
6
S
7
S
8
Righ
t
Bad
1st
Generation
20th
46
6. Future work
8. Acknowledgements
7. Conclusions
3.
With well formed primary and secondary IF-THEN
blocks, this algorithm is capable of converging to a
good solution very quickly. In our tests it was revealed
that the robot was capable of navigating the maze
sufficiently fast with a minor percentage (below 10%)
of overridden network decisions in very few iterations
of the algorithm without any further input from a human
expert. However, after a while the quality of the
artificial network started degrading considerably, with
increasing twitchiness of control by neural network,
which led to IF-THEN overrides and thus decreased
average speed. A further work on avoiding such
degradation or finding a best neural network and
stopping the teaching process is required.
4.
5.
6.
47
control,
Lithuanian
voice
1. Introduction
Speech based interfaces for the control of various home
appliances and as a basic modality for the control of
telephony based services has enormous potential and
advantages. In recent years we saw significant number
of new such type services and prototypes presented and
developed in various countries and in various fields.
There were two important factors that influenced
development and caused growing number of
applications: appearance of well established and easily
available telephony application program interface
(TAPI) and tools for development of such type
interfaces and even more importantly appearance of
speech application programming interfaces (SAPI) with
appropriate speech engines. Such speech recognition
engines or engines integrated into Microsoft Speech
Server (OCS-SS) enables to achieve highly accurate
speaker independent recognition of speech commands
necessary for the control of various speech based
telephony services.
Speech recognition engines have one very important
factor: they are not language independent and depend on
the language spoken by the users of the system [1].
Microsoft and other major speech technologies
providers develop and market speech engines for
popular languages having many native speakers
(English, Spanish, French, etc.). Many other less widely
used languages remains out of the scope of interest for
the major speech recognition solution providers.
In such situation businesses and state institutions, in
countries were such less popular languages are used as a
48
49
4. Experiments
The prototype of Lithuanian speech recognition system
for telephony services aims to recognize a set of
Lithuanian names and family names. The examples of
such services could be the phone directory of small or
medium enterprise or information service aimed to the
companys personnel. Since the variety of names and
particularly family names is very big weve limited the
number to 100 names and surnames.
4.1. Speech corpora
For the Lithuanian first and family names speech
database 40 different speakers uttered each name and
surname once. The total number of combinations first
name + family name was equal to 100 in this case (4000
utterances total). Two speech databases were applied in
the experiments: the initial corpora and the corrected
corpora. Corrected database was freed from various
inadequacies and mistakes that were present in the initial
database. Between inadequacies and mistakes were the
pronunciation errors. Most of such pronunciation errors
were situations when speaker used other phoneme than
the phoneme present in the original family name, but
still getting grammatically correct and often really
existing name (in example, speaker said Pojarskaite
instead of Bojarskaite or Laimis instead of Laimius).
Other errors were related with such problems as a
stammer near family name or some technical spoilage
such as a truncated part of the word (most often at the
end of the name the error of automatic segmentation
system).
4.2. Lithuanian
experiments
proper
names
recognition
50
Speakers
20 males
20 females
all 40
speakers
all 40
speakers
C, %
89.2
91.2
90.1
90.8
Recognition results
U, % S, %
6.8
3.8
4.5
3.9
5.7
3.8
5.4
O, %
0.2
0.4
0.4
3.6
Butkus Valdas
Daugla Andrius
Domarkien Inga
Diaugys Audrius
Giedraitis Laimius
Gudas Audrius
Ignataviius Gintaras
Ivanceviit Erika
0.2
Recognized
as
Gedutis
Eugenijus
Butkus
Valdas
Dutkeviit Ieva
Should be
Giedra
Nerijus
Gudas
Audrius
Ivanceviit
Erika
Corpora error
Adomonyt Asta
Aliaustien Jadvyga
Reading in syllables
Andiulyt Loreta
Averjovien Violeta
Bartu.kait Aura
Bilat Audeniu
Binnkeviius Ramnas
Pojarskait Egl
Background noise, pant
Buinskas Arnas
Buloctien
Danguol,
background noise
background noise
Daudla Andrius
Domartien Inga
Diaugys Andrius
Giedraitis Laimis
background noise
First letter cut
...aneviit Erika
Aleksandraviit Ramut
Anereviit Giedr
Andziulyt Loreta
Bartukait Aura
Biaigo Sandra
Bojarskait Egl
Bulaxien Danguol
Domarkien Inga
Drungilien Gintar
Recognized Should be
as
Gedutis
Giedra
Eugenijus
Nerijus
Butkus
Gudas
Valdas
Audrius
Balcius
Ernestas
Gailiunas
Rytis
Dolgij Andrej
Grigonis
Audrius
Baublys Algis
Giedra Nerijus
10
51
SNR 20 lg
smax
,
nmax
(1)
smax max(abs ( si )) ,
i1, N
(3)
Clipped signal
equations (4):
s,
sic c
sc ,
if
si sc
if
si sc
Table 5. The four most frequently recognized transcription for the two speaker uttering 2 Lithuanian names
Name
Beliukeviit
speaker
ALGRUD
transcr.
occur.
111
23
99
20
15
13
3
10
Varanauskas
speaker
RYTMAS
transcr.
occur.
505
18
121
12
504
11
507
11
ALGRUD
transcr.
occur.
67
24
130
11
4
8
70
8
52
RYTMAS
transcr.
occur.
19
10
166
8
144
8
6
8
(4)
Corpora
unmodified
SNR 40 dB
SNR 30 dB
SNR 20 dB
Clipping 0,3
Clipping 0,1
C, %
96.2
95.8
91.1
43.0
92.5
81.0
Recognition results
U, %
S, %
O, %
1.3
2.3
0.2
1.8
2.3
0.1
1.9
6.9
0.1
12.0
30.4
14.6
5
2.5
0.0
13.2
4.8
1.0
6. References
1. Fung, P., Schultz, T. Multilingual spoken language
processing. In: IEEE Signal Processing Magazine,
Vol. 25, Issue 3, pp. 89-97 (2008).
2. Zgank A., et al. The COST278 MASPER initiative
crosslingual speech recognition with large telephone
databases. In: Proc. of 4th International Conference
on Language Resources and Evaluation LREC04,
pp. 2107 2110, Lisbon (2004).
3. Kasparaitis P. Lithuanian Speech Recognition Using
the English Recognizer. In: INFORMATICA 2008,
Vol. 19, No 4, pp. 505516, Institute of Mathematics
and Informatics, Vilnius (2008).
4. Beyerlein, P., et.al. Towards Language Independent
Acoustic Modeling, In: Proc. of IEEE ICASSP00,
Vol. 2, pp. II1029-II1032, Istanbul (2000).
5. Lindberg B., et al. Noise Robust Multilingual
Reference recognizer Based on Speech Dat. In: Proc.
Of ICSLP 2000, vol. 3, pp. 370373, Beijing (2000).
6. Villaseor-Pineda L., et. al. Toward Acoustic
Models for Languages with Limited Linguistic
Resources, In: CICLing-2005, Lecture Notes in
Computer Science, Vol. 3406/2005, pp. 433-436,
Springer-Verlag (2005).
7. Schultz T. and Kirchhoff, K. (Eds.), Multilingual
Speech Processing, Academic Press, 2006.
8. Maskeliunas R. Modeling Aspects of Multimodal
Lithuanian Human - Machine Interface. A. Esposito
et al. (Eds.): Multimodal Signals, LNAI 5398, pp.
7582, Springer-Verlag Berlin, Heidelberg (2009).
5. Conclusions
Speech recognition engines are not language
independent and depend on the language spoken by the
users of the system. Two major ways for Lithuanian
speech recognition solution are: a) developing own
speech recognition engine from scratch; b) adapting the
foreign language based engine for the recognition of the
native language.
The application of the English based engine for the 100
Lithuanian names recognition allowed achieving
recognition accuracy of more than 90 %. These results
show that the implementation of longer commands and
transcription generation method proposed by [3] were
confirmed.
The 5 most confusing names produced almost 40 % of
all substitution errors and slightly more than 35 % of all
insertion errors. Most of those names dont have
particularly difficult phonetic structure (the family name
Biaigo may be treated as the more difficult case). The
bigger number of errors obtained by the name Gailiunas
Rytis may be explained by the presence of the name
Gailiunas Vytautas in the same list. But most of those
errors cant be explained straightforwardly. For
example, the name Gudas often was confused with the
name Butkus. Such situation shows, that theres still
plenty of room for the optimization and further
investigation and the better recognition accuracy is
obtainable.
The most often recognized (as the most similar)
transcriptions arent the same for the different speakers.
A big number of transcriptions were recognized as the
most similar ones for each of the speaker (27
transcription variations of the name Beliukeviciute 27
for the first speaker and 21 for the second speaker).
Each speaker tended to activate more different
transcriptions than the second.
53
Client
Client
Client
1. Introduction
Client
Client
Ethernet
SERVER
MES
Systems
SCADA/HMI
Systems
PLC/IPC
PLC
Fielbus
ERP
Systems
Fielbus
Industrial
devices
54
DA
DA
DX
SERVER
SERVER
DA
DA
CLIENT
OPC DA
OPC
SERVER
Communication
protocol
Computer
network
The system shown in Fig. 4 is employed to test loadcarrying capacity of a computer network using OPC.
PC1 and PC2 computers are connected to Access Point
AP TP-LINK, TL-WR543G. Access Point acts as a link
in the process of information sharing among computers
in a LAN following TCP/IP protocol.
OPC server: "Matrikon Simulation Server" process
data simulator - is installed on PC1. Readout of one type
of process datum 4-bite integer is assumed, as shown
in Fig. 5. Software of OPC client, created by Delphi
compiler of Borland, educational version, is installed on
PC2. OPC "Sentrol" components are used to create the
application, accelerating its construction at the stage of
design and sending orders to OPC server in the
Value
OPC Item
Timestamp
Quality
55
AP
4000
[kbps]
3500
PC 1
PC 2
3000
2500
2000
1500
1000
500
0
0
2000
4000
6000
8000
10000
12000
Item s/s
Dow nload DL
Upload UL
OPC DL+UL
N 24 13 n ,
(1)
where:
N the size of 4-bytes communication protocol data,
n quantity of process data,
24 protocol header including the table name = 8 +
16 = 24 bytes,
13 value + quality + timestamp = 4 + 1 + 8 = 13
bytes.
CIM (Crystal Internet Meter) is used to measure loadcarrying capacity of the network. This software enables
to measure quantities of UL (Upload) and DL
(Download). The measurement is possible in KB/s
(kilobyte per second) and kbps (kilobit per second). A
process of LAN load-carrying capacity is shown in
Fig. 6.
The measurements are carried out for several different
loads of process data read out. The measurement results
are presented in Fig. 7.
56
[kbps]
3500
3000
2500
2000
1500
1000
500
0
0
2000
4000
6000
8000
10000
12000
Item s/s
OPC DL+UL
Measured DataSet
7. References
Counted DataSet
1.
Fig. 8. Comparative characteristic curves of LAN loadcarrying capacity for OPC and information transmission by
means of DataSet technology
2.
3.
6. Conclusion
4.
5.
6.
7.
8.
57
Inlet water
TP
NP2
MV1
Outlet water
1. Introduction
NP1
PIC
TT
TIC
B1
B2
Supplement water
Gas
Air
NS1
NS2
58
T2 T2
T1 T1
Time
T2
T1
Time
PV
M1
Reservation
Digital
Analogue
SP
Controller
M2
M3
Faults
Fig. 4. Pressure control algorithm block diagram
59
Q1
Q2
Q3
SP
F2
A
2
3
F3 Fn A1A2 An A3 B1B2 Bn
A NN AA NNN A AA NN
IIS
1
F1
A
FB
Typical FB
Fig. 5. Automated control algorithm functional scheme
Continuing the development of an automated subprogram exported KKS signal list must be examined.
For this purpose a prototype program using visual Basic
60
File.csv
A=Filter(File, KKS)
A(n)
F
n > -1
FB(i)=C
FB(i)=A(n)
6. Conclusions
File.exp
7. References
1. Gimbutis G. Thermo-technology. (In Lithuanian),
Vilnius, 1993. 333 p.
2. Ezerskis D. Development and Application of
Programming Methodology for Automation Control
Systems, Based on Design Patterns. Dissertation,
Kaunas, 2006. 128 p.
3. LRS. Requirements for boiler-house equipment
installation. (In Lithuanian), Vilnius, 2006.
http://www3.lrs.lt/.
4. IEC. IEC 61131-3 Programmable controllers Part
3: Programming languages, Second edition 2003.
http://www.plcopen.org/.
5. VGB. Identification System for Power Stations
(KKS). http://www.vgb.org/.
61
62
tin ten
0.75
K=
t3
,
te
(6)
(2)
t t
2n tin te ,
(3)
t B = 3n
t
t
2
in en
(4)
Q = c g (t3 t2 ) = sp Fsp 3 2 ti ,
2
63
60
12
10
50
8
T, C
40
6
4
30
2
20
Tfeed, C
Tret, C
Tair, C
-2
10
-4
:3
0
:2
0
:0
5
:5
5
:4
5
:4
0
03
01
23
21
18
:2
5
:1
0
:0
0
:5
0
:3
5
16
14
12
10
08
05
03
01
:4
0
-6
:3
0
time
64
65
W2
-
W1
-
A1
E1
A2
E2
A1
A2
E1
E2
D1
D2
D3
1. Introduction
D4
66
2.
3.
4.
5.
6.
PFDavg T / 2 ,
GPS
PIC4
(1)
Steam
PIS5
LIC1
water
1/ MTBF ,
(2)
Natural
gas
SIL
SIL1
SIL2
SIL3
SIL4
V1
PIC3
PS7
Burner
No.2
Burner
No.1
air
No.
1
2
3
4
5
6
PS6
PIC2
PFDavg
10-2<PFDavg<10-1
10-3<PFDavg<10-2
10-4<PFDavg<10-3
10-5<PFDavg<10-4
67
K 2.1
K 52.1
(4)
K 46.1
E 63,08
E 62,15
E 62,14
E 62,13
SR 3
322
320
318
1.1
323
K 51.1
314
313
312
316
T T p n1 n 24 365 10 87600 h
E 62,12
E 62,11
E 62,10
E 62,09
E 62,08
ZP1
07 KT 97
8 DI
DG 1.1
DG 2.1
325
K2
(11)
68
, (17)
Device
Differential pressure relay
Pressure relay
Flame scanner
Natural gas shut-off valve
Boiler fail safe PLC
5. Conclusions
The risk analysis of middle power steam boiler shows
that protections system of such boiler should fulfil SIL3
according to EN 61508. Analysis of a typical existing
protection system shows demand for improving existing
protection systems.
Relay based systems should be used only for boilers
those require SIL2 protection equipment. In this case
only certified SIL sensors and actuators could be used.
Using standard PLC in a start-up sequence of a boiler
for flame sensor bypass improves overall system PFD
and makes system unsafe.
Only fail safe PLCs should be used for protection
systems and boiler start-up system.
8.710-9
6.110-9
1.2710-8
4.2810-8
7.6910-9
6. References
1. American Institute of Chemicals Engineers
Guidelines for process equipment reliability data
with data tables, New-York, 1989. 312 p.
2. Boercsoek J. Functional Safety. Basic Principles of
Safety-related Systems Heidelberg, 2007. 538 p.
3. Ezerskis D. Development and Application of
Programming Methodology for Automation Control
Systems, Based on Design Patterns. Dissertation,
Kaunas, 2006. 128 p.
4. Exida Safety engineering Exida, 2005. 275 p.
5. IEC
61508-1.
Functional
safety
of
electrical/electronic/programmable electronic safetyrelated systems. General requirements. Geneva,
1998. 115 p.
6. IEC
61508-5.
Functional
safety
of
electrical/electronic/programmable electronic safetyrelated systems. Examples of methods for the
determination of safety integrity levels. Geneva,
1998. 57 p.
7. OREDA Ofshore reliability data book, 4th edition,
Trodheim, 2002. 835 p.6.
8. http://www.exida.com/applications/sael/index2.asp fail safe equipment database
(16)
69
70
X R cos( ),
Y R sin( ),
Z H.
(1)
(2)
(3)
71
V s t ,
H V H s t.
(4)
(5)
72
Rt
Rabs Rt _ max
Rabs _ max
l ,
(6)
73
7. Conclusions
In the paper, stepper motor positioning and position
tracking problems, typical in the programming of 3D
scanning systems are presented and the solutions to the
problems are suggested.
It was experimentally shown that, motors, used in the
3D scanner should have encoders in order to increase
scanning speed and quality. It was also demonstrated,
that even if it is possible to precisely control the rotation
of the stepper motor, it is impossible to precisely
calculate its shaft rotation position evaluating the time in
PC-based systems without the encoders.
In case stepper motors are not equipped with encoders,
the modelling errors can be reduced by installing
external proximity switches that indicate some positions
of the motor shaft. An intermediate solution can also be
used where the motor control based positioning is
applied. This method is slow; nevertheless rate of the
positioning errors can be reduced to the desirable level
at the expense of the scanning speed. It should be
noticed that not all PLC systems supports this method
also called referencing or homing.
Time controlled positioning should be avoided in nonreal time systems but can be used in PLC systems were
equal cycle execution periods are ensured.
8. References
1. TwinCat Information System. Documentation for
KL2531 and KL2541, One Channel Stepper Motor
Terminals, Version 1.5, 2006.
2. Davidsonas A., Liutkeviius R. A PLC Application
for Data Scanning and Visualization. Proceedings of
International Conference Electrical and Control
Technologies 2008, Kaunas, 2008. p. 144-149.
74
d cos
,
2l
(1)
luminance contrast
A target can only be seen if it differs in luminance from
background. On this understanding luminance contrast
which defines visibility may be calculated. There are
several ways of defining luminance contrast, though
here is traditional one:
C = (Lt Lb) / Lb ,
(2)
(3)
77
l w
3.44 W
(4)
6. Conclusion
1. Optimal value of visual angle in slide is 2 for
subscripts and superscripts and 3 for main text. If
its aesthetic and there is free space in slide larger
fonts can be used.
2. In case of good lighting (300500 lx) contrast
between text and background must be no less than
0.5, i. e. screen illuminance created by projector
must be no less than twice bigger than illuminance
created by surroundings.
3. If visual angle is greater than 3 and contrast is
higher than 0.5, identifying speed with confidence
probability 0.99 exceeds practicable avarage reading
rate of english text.
7. References
1. Murdoch J. B. Illumination Engineering from
Edisons Lamp to the Laser. New York, 1985.
p. 542.
2. The IESNA lighting handbook. New York, 2000.
p. 800.
3. Masiokas S., Kriuglait M., Otas K. Colour
luminance contrast of digital projection image.
Electronics and Electrical Engineering, no. 7(79).
Kaunas, 2007. p. 19-22.
4. Shaikevic
A. S.
About classifying visual
perfomances. Lighting engineering, no 5. Moscow,
1958. p. 13-20 (in Russian).
5. Shaikevic A. S. Nomogram for assessing identifying
speed and confidence probability. Lighting
engineering, no 1. Moscow, 1964. p. 19-22
(in Russian).
78
2. Decision Trees
Goal of the classification is the assignment of a target
variable attribute to an unknown object. The
classification is done by a classifier which learns the
classification based on classified training data. To build
a classifier can be several methods like neuronal
networks [9], support vector machines [4], fuzzy
systems [8,6] or trees [2,3,5,7] used. The classifier in
this work based on trees.
Trees parting a feature space into a subset of rectangles.
An attribute of a target variable, a label, is assigned to
every rectangle. There exist two different types of tree
based models; regression trees and classification trees.
Regression trees work with continuous target variables
and try to minimize the distance between the label and
the target variable values in the feature space subset.
Classification trees work with discrete target variables.
They try to achieve that a subset only contains attribute
for the target variable which are identical with the
subset label. [5]
This work uses classification trees and as learning
method Quinlans C4.5 algorithm.
The C4.5 algorithm is an improved version of the ID34
algorithm.[7] The learning algorithm searches for a
hierarchical split structure to divide the data set T into
subsets Ti with the goal that a subset at the lowest split
layer only contains objects with the same class Ck.
1. Introduction
Decision trees are widely used for classification
problems. Their advantages compared with other
classification techniques are:
Decision trees can be effectively applied to any data
structure, in particular to discrete, continuous or
mixed data. [3]
Decision trees are very robust against outliers. [3]
Prediction rules are easy to interpret. [3]
Decision trees perform a stepwise variable selection
and complexity reduction. [3]
A decision tree is generated out of evaluated feature
data. The most common methods to create a decision
tree are C4.5 [7], CART1 [3] and CHAID2 [2]. The
generated tree is often very complex and overfits the
classification problem. Hence to improve the
classification ability is a tree pruning necessary. Pruned
trees are the base for a real word classification
application. Therefore is a decision tree implementation
on a target environment required, for example on a
microcontroller or a PLC3.
A manual conversion of the develop decision tree into a
target source code is often not expedient. Automatic
1
2
3
79
Iterative Dichotomiser 3
(1)
info information content
pm probability distribution of the classes
(3)
gain information gain of the split
(5)
gain ratio normalized information gain
80
Fuzzy
classifier
Bayes classifier
direct
crestline
regression regression
Tree
2.0
2.0
2.0
0.0
2.0
2.0
2.0
2.0
4.5
38.5
3.5
4.0
81
7. References
1. Norm. IEC 61131-3:2003. Speicherprogrammierbare
Steuerungen Teil 3: Programmiersprachen
2. Bankhofer U. Datenanalyse und Statistik: Eine
Einfhrung fr konomen im Bachelor. Gabler.
Wiesbaden, 2008.
3. Breiman L. Classification and regression trees.
Chapman & Hall. New York, 1993.
4. Cristianini N. An introduction to support vector
machines and other kernel-based learning methods.
Cambridge Univ. Press. Cambridge, 2002.
5. Hastie T. The elements of statistical learning: Data
mining, inference, and prediction. Springer. New
York, 2001.
6. Koch M. Klassifikatorkonzepte zur Steuerung
dynamischer Prozesse. Ilmenau, 1995.
7. Quinlan J. C4.5 programs for machine learning.
Kaufmann. San Mateo, 1993.
8. Russo M., Jain C. Fuzzy learning and applications.
CRC Press. Boca Raton, 2001.
9. Seiffert U. Self-organizing neural networks: Recent
advances
and
applications.
Physica-Verl.
Heidelberg, 2002.
10. Seitz M. Speicherprogrammierbare Steuerungen:
Von den Grundlagen der Prozessautomatisierung bis
zur vertikalen Integration. Fachbuchverl. Leipzig im
Hanser-Verl. Mnchen, 2003.
82
83
relPalpha
Palpha
Ptotal
100% ,
2. Methods
2.1. Integration of the SD Card device
For integrating a SD card slot into a microcontroller
driven device it is necessary to implement an algorithm
to manage the file system. After this the SD card can be
easily red by a personal computer (PC) as well as in
other handheld devices with SD card support.
For implementing the file allocation table (FAT)
management on the microcontroller device we use the
Portable FAT Library for MCU Applications [3].
There are two different ways to communicate with a SD
card. On the one hand in the fast 4- bit communication
mode and on the other hand in the simple serial
peripheral interfaces (SPI) mode. We decided to use the
SPI mode. This modus turned out to be as fast enough
to transmit the data during one sampling period. This
mode is efficient enough to implement it on a
microcontroller with small program memory.
The device with integrated SD- card data storing is
shown in Fig.2.
(1)
f sample 2 f max ,
(2)
84
85
relP0.1Hz
P0.1Hz
100% ,
Ptotal
4. References
1. Baumann, C. RFT: a simplified fast real time sliding
DFT algorithm, Convict Episcopal de Luxembourghttp://www.restena.lu/convict/Jeunes/ultimate_stuff
,RFT.pdf, 2005.
2. Baumgart-Schmitt, R. systemanalytische Anstze
in der kognitiven Psychophysiologie unter
besonderer Bercksichtigung stationrer visueller
Strukturbildungsprozesse, Humboldt Universitt
Berlin, 1986.
3. Sham I., Hue W., Rizun P. Portable FAT Library for
MCU
Applicationshttp://www.circuitcellar.com/
magazine #176, March 2005.
4. Peng, C.-K. Heart rate dynamics during three forms
of meditation, International Journal of Cardiology,
February 2003.
5. Oida E., Moritani T., Yamori Y. Tone-entropy
analysis on cardiac recovery after dynamic exercise
Proceeding of 4th eMBEC 2008
6. Kamiya N., Yokoyama K., Niwa S. Evaluation of
Parasympathetic Nerve Activity Level Monitoring
System, IEIC Technical Report, VOL.102, NO.727,
2003.
(3)
PI (n)
(4)
for
1 n N 1,
(5)
A( s )
4 2 f 22 s 2
(1)
2 2 f 1
2 f 2
s
s 4 2 f 12 s 2
s 4 2 f 22
Q1
Q1
87
Filter work
Input
signal
Time [s]
Amplitude
Time [s]
Amplification [dB]
> - 12 dB/octave
>12 dB/
/octave
Tolerance [dB]
0.4
+ 1 dB
0.25
10
+ 0.5 dB
0.4
[Hz]
+ 1 dB
10
100
[Hz]
88
x ( k )
u1( k 1)
1
y ( k ) u1( k )
b2
11u1( k ) 12 u 2( k ) 1 x ( k )
(3)
u 2( k 1) 21u1( k ) 22 u 2( k ) 2 x ( k )
Fig. 4. Block diagram of system correction, awe transducer
input acceleration, awz standard acceleration, ap transducer
output acceleration, af acceleration after filtration, ak
acceleration after correction
y(k-1)
+-
u^1(k-1)
delay
11
u^1(k)
(2)
12
u^2(k-1)
22
+
+
+
+++
delay
G ( s)
^
x(k-1)
1/b2
u^2(k)
Wr
89
af
Acceleration [m/s2]
awz
ak
Time [s]
Ap
Af
Ak
3.15
16.66
5.89
3.31
Acceleration [m/s2]
ak
Time [s]
b)
af
ak
awz
Time [s]
90
Velocity [m/s]
Standard accelerometer
Measuring accelerometer
Before correction
*1
*2
After correction
Standard accelerometer
Measuring accelerometer
Before correction
After correction
*1
*2
Time [s]
(ym(t)- yrz(t))
Before correction
After correction
Time [s]
91
Before correction
3.
4.
After correction
5.
6.
Time [s]
7.
After
correction
15
128.89 129.28
0.17 0.06
1102
263
14.5
217.7
325.2
0.01
8.
9.
2.054 2.196
7.143
6.674
10.
0.0026
0.0035
4. Conclusion
11.
12.
13.
5. References
1.
2.
92
F ( x) ai f i ( x) ,
y G (x) ,
(2)
i 1
2
a arg min y i F (x i ) .
a
(3)
i 1
(1)
F ( x) ai x i .
(4)
i 0
r
f i ( x) x jij ,
(5)
j 1
94
j 1
RMSE
y i F (x i ) , (6)
100% i
2
SD
y i y
2
RRMSE 100%
3. Applications
This section describes a number of regression modelling
applications where the ABFC methods were used as
well as demonstrates results of empirical comparisons
of the ABFC methods to other well known and widely
used regression modelling methods.
The methods compared are the following: 1) full
polynomials (FP) of different degree; 2) sparse
polynomials of different degree built using one of the
most popular subset selection methods SFS utilizing
AICC as the evaluation criterion; 3) F-ABFC; 4) EFABFC; 5) Locally Weighted Polynomials (LWP) of
different degree utilizing cross-validated Gaussian
95
processes
in
Switched
4. Conclusions
In this paper we shortly reviewed the ABFC approach,
briefly characterized two of its special cases F-ABFC
and EF-ABFC, and demonstrated the efficiency of the
methods in a number of applications in the field
electrical and control technologies while comparing
them to other widely used methods.
The ABFC methods, in contrast to the methods of
subset selection, do not require the user to predefine the
maximal degree or the dictionary of the basis functions.
Instead they automatically adapt to the particular data at
hand.
Overall, in the performed empirical experiments it is
concluded that in the considered problems the adaptive
model building methods of ABFC can outperform the
other polynomial regression modelling methods and
compare rather well to the other state-of-the-art
methods MARS, RBF, and Kriging. The ABFC
methods can be viewed as competitive tools for
regression modelling.
Directions of future research include employment of
ABFC or derivations of it in different real-world
applications, further theoretical studies concerning the
efficiency of ABFC, as well as adaptation of ABFC for
other types of tasks, e.g., classification.
F-ABFC and EF-ABFC together with many other
regression modelling methods are implemented in
VariReg software tool freely available for noncommercial research and educational purposes at
http://www.cs.rtu.lv/jekabsons/.
Method
FP, p = 1
FP, p = 2
FP, p = 3
FP, p = 4
FP, p = 5
SFS, p = 1
SFS, p = 2
SFS, p = 3
SFS, p = 4
SFS, p = 5
LWP, p = 1
LWP, p = 2
LWP, p = 3
LWP, p = 4
LWP, p = 5
F-ABFC
EF-ABFC
MARS
MARS, CV
RBF
Kriging
P1
91.88
57.79
31.48
27.65
18.06
91.88
56.83
35.65
29.66
27.12
28.70
20.35
23.01
16.17
16.38
17.60
12.77
14.69
15.60
11.74
14.45
P2
P3
29.94 79.29
35.37 52.30
- 44.44
- 61.85
- 382.65
29.94 79.29
28.58 54.31
49.67 49.36
53.91 41.06
72.73 39.59
25.75 48.84
35.68 45.90
- 48.81
- 127.84
- 305.94
41.15 41.66
26.25 42.04
35.54 38.63
34.88 38.24
27.19 39.67
28.02 34.02
P4
12.33
13.04
12.02
16.04
12.32
15.98
9.36
17.77
15.13
16.60
9.38
P5
36.19
29.83
13.20
12.76
7.80
36.19
28.79
13.21
9.49
10.09
12.89
15.51
9.43
10.05
12.41
8.00
14.58
36.31
28.72
12.60
4.66
5. References
1. Breiman L. Heuristics of instability and stabilization
in model selection. Annals of Statistics, vol. 24,
1996. p. 2350-2383.
2. Cherkassky V., Mulier F.M. Learning from Data:
Concepts, Theory, and Methods, 2nd ed. WileyIEEE Press, 2007. 538 p.
3. Cleveland W.S., Devlin S.J. Locally weighted
regression: An approach to regression analysis by
local fitting, Journal of the American Statistical
Association, vol. 83, 1988. p. 596-610.
4. Friedman J.H. Multivariate Adaptive Regression
Splines (with discussion), The Annals of Statistics,
vol. 19, no. 1, 1991. p. 141.
5. Gutmann H.-M. A Radial Basis Function Method
for Global Optimization, Journal of Global
Optimization, vol. 19, 2001. p. 201-227.
6. Hastie T., Tibshirani R., Friedman J. The Elements
of Statistical Learning. Springer, 2003. 552 p.
97
98
cation accuracy of over 95% was achieved when categorizing the images into one healthy and two pathological (nodular and diffuse) classes. Fig. 1 presents
characteristic examples from the three decision classes
considered, namely, nodular, diffuse, and healthy. As
can be seen from Fig. 1, examples representing the
different classes are rather different. However, it is
worth noting that due to the large variety of appearance of vocal fold mass lesions, the classification task
can sometimes be difficult to solve even for a trained
physician [10].
Introduction
Videostroboscopy
99
vocal fold vibrations exhibit only one single fundamental frequency. Multiple tones (fundamental frequencies) may be recorded in the case of some diseases, such as polyps, nodules, and cysts [14]. In such
cases, a clear view of the vibrating vocal folds can
not be obtained with the single-flash-timing laryngeal
videostroboscopy. A multiple-flash-timing laryngeal
videostroboscopy technique was proposed by Deguchi
et al. [14], to deal with such cases. Multiple light
emitting diodes are used as illumination sources.
In [15] image sequences recorded with the stroboscopy
system have been used to measure the glottic angle and the angular velocities of vocal fold abduction and adduction. The authors point out that semiautomated edge tracking would be an important improvement of the technique.
It is worth mentioning that not only edge tracking
but also other tasks usually carried out when analyzing video data need automated or semi-automated
procedures. Decision making is one of such tasks. In
clinical practice, decision making is quite often based
on subjective evaluation of video data. Quantitative
measures of motion, colour distribution and geometry
of vocal folds can provide objective information and
be useful in medical treatment planning and greatly
facilitate tracing progress over time.
The long-term goal of this work is a decision support
system for diagnostics of laryngeal diseases. A voice
signal, sequences of colour vocal folds images obtained
from videostroboscopy, and questionnaire data [16]
are the information sources to be used in the analysis. This paper is concerned with automated characterization of image sequences obtained from laryngeal
videostroboscopy into a healthy class and two classes
of disorders.
2.
e i)
i (
(2)
i=1
e i ), (
e j )i of the
of the covariance matrix Kij = h(
centered data points, where the expansion coefficients
i of the eigenvector are found from the eigenvalue
problem
= K,
(3)
The data
The medical task considered in this paper concerns laryngeal colour image sequences, based automated categorization of laryngeal disorders into three decision
classes: healthy and two pathological classes, namely
diffuse and nodular mass lesions of vocal folds [7].
Mass lesions of vocal folds could be categorized into
six classes namely, polypus, papillomata, carcinoma,
cysts, keratosis, and nodules. This categorization is
based on clinical signs and histological structure of
the mass lesions of vocal folds. We distinguished
two groups of mass lesions of vocal folds i.e. nodular
lesions (localized thickenings)nodules, polyps, and
cystsand diffuse lesionspapillomata, hyperplastic
laryngitis with keratosis, and carcinoma.
The data have been recorded at the Department
of Otolaryngology, Kaunas University of Medicine,
Lithuania. The image sequences were acquired during
routine videostroboscopy. The duration of one image
sequence was equal to 4s. The resolution of 720 576
pixels was used to record the image sequences.
3.
M
X
M
X
e i ), ()i.
e
i h(
(4)
i=1
The classifier
Depending on the definition of the optimization problem, several forms of support vector machine (SVM)
can be distinguished, for example, 1-norm or 2-norm
SVM. Assuming that (x) is the non-linear mapping
of the data point x into the new space, the 1-norm
soft margin SVM can be constructed by solving the
following minimization problem:
M
X
1
i
(5)
min wT w +
w,b, 2
i=1
subject to
Features
yi (hw, (xi )i + b) 1 i ,
i 0,
i = 1, ..., M ,
Various types of features characterizing colour, texture, and geometry of biological structures seen in
100
(6)
f (x) = H
i yj k(x, xi ) + b ,
(7)
between the spaces and derive a positive definite kernel based on that concept. The "QR" factorization of
the matrices (A, B) and the kernel Gram-Schmidt orthogonalization process are used to derive the kernel.
Applying the "QR" factorization the matrices (A, B)
can be written as A = QA RA and B = QB RB , where
Q is an orthonormal basis and R is an upper-diagonal
matrix of size M M of the Gram-Schmidt coefficients representing the columns of the original matrix
in the new basis. The principal angles cos(i ) are
given by the singular values i of the matrix QTA QB ,
cos(i ) = i , i = 1, ..., M . It was shown that
i=1
(A, B) = det(QTA QB )2 =
1X 2
1 T
w w+
e
2
2 i=1 i
5.
i = 1, ..., M.
5.1.
Experimental setup
(9)
1
where Zij = yi yj (xi , xj ), I is the identity matrix,
1 = [11 , ..., 1M ], y = [y1 , ..., yM ], and = [1 , ..., M ].
Since an SVM is a binary classifier while the task is
to distinguish between three classes, the one-againstone scheme is used to make the categorization in this
work.
4.2.
Experimental investigations
(8)
subject to
yi (hw, (xi )i + b) = 1 ei ,
(11)
w,b,e
cos(i )2
i=1
min
M
Y
5.2.
Results
Kernel function
For (xi , xj ), one usually uses the linear: xTi xj , Gaussian: exp{||xi xj ||2 /} or polynomial: (xTi xj +1)d
kernel. The kernel is defined over a pair of vectors.
In this work, classification is based on a set of vectors
rather than on a single vector. For example, a sequence of images is recorded from a patient. Each
image is represented by a feature vector. Feature
vectors are then collected into a matrix (each vector constitutes a matrix column) and used to make
a decision. Therefore, a kernel function utilized by
the LS-SVM classifier is defined over a pair of matrices (A, B) rather than over a pair of vectors. A
positive definite kernel of such type has been recently
proposed by Wolf and Shashua [22]. The authors use
the principal angles between the two column spaces
defined by the matrices (A, B) to assess the matching
The graph presented in Fig. 2 plots the correct classification rate of the test set data as a function of the
percentage of the data variance accounted for by the
number of the kernel principal directions used. As
can be seen from Fig. 2, the percentage of the data
variance accounted for by the optimal number of the
components is close to 90. The graph presented in
Fig. 3 relates classification accuracy the regularization constant , and the number of the kernel principal components used.
101
90
7.
Accuracy, %
80
70
8.
60
50
9.
40
30
10
12
14
16
Number of kPC
10.
Fig. 2.
The correct classification rate of the test set
data as a function of the percentage of the data variance
accounted for by the number of the kernel principal
components used
11.
Accuracy, %
100
80
12.
60
40
13.
20
15
20
10
15
10
0
14.
Number of kPC
6.
15.
Conclusions
16.
17.
References
1.
2.
3.
4.
5.
6.
18.
19.
20.
21.
22.
23.
102
x Ax Bu, x(0) x 0
(1)
y Cx
x A (t ) x B (t )u K (t )( y y )
y [0 1]x
(2)
Plant
A , B
Estimation
Matlab/Simulink
State
Observer
2. State observer
103
0 2
2
A
, B 0, C 0 1
0.154 0.154
(3)
C1
R2
1 n
e(i)
n i 1
R4
y
R1
R3
C2
0 1.9681
4.1032
A
, B 0.0498
0.0291 0.0194
(4)
are chosen:
104
C2
R2
u
R1
C3
R4
R3
R6
R5
a ii , i 1,2
-20; -1
3.0254 10 4
-200; -10
3.5039 10 4
-0.2; -0.01
1.6916 10 4
-0.2; -0.2
0
0
0.3
1
A 0
0.2
0 , B 0, C 0 0 1
0
0
0.5
0
AME
2.538 10
-20; -20
4.6 10 4
3.9136 10 4
0
0.2
0.19547
A
0
0.2
0
0.145
0
0.898
36.204 0
A
3
0
bi , i 1,2
5; 2
3.3043 10 5
0.1; 5
2.8067 10 5
0.2; 0.2
3.0379 10 5
20; 20
4.0738 10 5
AME
105
4.252 10 4
a ii , i 1,2
-10; -5
1.9856 10 4
1.3555 10 4
-1; -0.5
9.6398 10 5
2.2967 10 4
1.33 10 4
1.0642 10 4
a ii , i 1,2,3
-3; -0.2; -4
AME
AME
6. Conclusions
10 0 1
A
, B , C 0 1
5
0
0
31.343 8.9445
7. References
106
107
e r .
e*
rr` iqs
e r ` e* .
Lr p ids
`*r ,
*
2
idse* may be
`*r
The desired torque of
`
r m
`
`
r
r
rL
idse* .
r L p
Teme*
(1)
(2)
i as*
i qse *
Teme*
(3)
*
rm
maintaining
flux,
`dre Lmidse*
idr`e is
i bs*
`e
dr
i qse * , i dse *
i cs*
e*
ds
2*
2*
Fig. 2. Indirect field-oriented control of a current regulated PWM inverter induction motor drive
4. Simulation of transients
Simulation results witch were got from Fig. 3 are shown
in Fig 4. The greatest current value at starting reaches
110 A, but after motor starting, steady state currents
amplitudes are approximately 25 A at motor starting
without load. The torque change low is from 0.5 s
suddenly grows tile nominal, from 1 s its goes down tile
haft of nominal torque of machine, next from 1.25 s its
grows tile nominal and after 1.5 s its drop tile zero.
Response to that function can be seen in the second
graph. During the load torque changing, motor windings
currents change also. They repeat load changing
behaviour.
(4)
108
Voltage, V
Tem, Nm
100
50
0
-50
IPsist, Wb
150
100
50
Currents, A
0
100
50
0
-50
0
0,2
0,4
0,6
0,8
1
1,2
Time, s
1,4
1,6
1,8
109
6. Conclusions
Simulink model of current space vector control is
created.
Acquired simulation results allow making these
conclusions:
1. Developed motor torque shows that during the startup a slight overshoot exists. After the system start-up,
the applied torque is changed by proper low for
investigation of changes of motor.
2. The current composes constant part conditioned by
motor resistances and variable part conditioned by
motor load.
3. The speed error at motor starting due to static friction
is insignificant.
4. PWM space vector control way allows using by order
slower microprocessors practically without loss
converter output signal quality.
7. References
1. Bierke S. Enchanced control of an alternating current
motor using fuzzy logic and a TMS320 digital signal
processor, Texas Instruments, 1996. 524 p.
2. Valentine R., Motor control electronics handbook
[CD-ROM] 1998.
3. Kliuchev I.V. Theory of electrical drives. (In Russian)
Moscow: Energoatomizdat, 1985. 253 p.
4. Holts J. Pulsewidth modulation - a survey, IEEE,
1992, 13 p.
5. Caron J.P., Hautier J.P. Modelling and Control of
Induction Machine. Technip Edition, France, 1995, 312
p.
6. Field Orientated Control of 3-Phase AC-Motors.
Application report, Texas Instruments, USA, 1998, 473
p.
7. Hazzab A., Bousserhane I.K., Kamli M., Rahli M.
New Adaptive fuzzy PI-Sliding Mode Controller for
Induction Machine Speed Control. Third IEEE
International Conference on Systems, Signals & Devices
SSD'05, Tunisia, 2005. 125 p.
8. Hazzab A., Bousserhane I.K., Kamli M., Rahli M.
Design of fuzzy sliding mode controller by genetic
algorithms for induction machine speed control. Third
IEEE International Conference on Systems, Signals &
Devices SSD'05, Tunisia, 2005. 56 p.
9. Lorenz R.D., Lawson D.B. Simplified A Approach to
Continuous On-Line Tuning of Field-Oriented Induction
Machine Drives. IEEE Trans. On Industry application,
Vol.26, Issue 3, May/June 1990. 425 p.
10. Geleeviius V., Kriinas K., Kubilius V. Control
systems of electrical drives. (In Lithuanian) Vilnius:
Mokslas, 190. 360 p.
110
M di1 L di2 i R ,
2
2 2
12 dt
dt
(1)
111
d 2 i1 R1 L2 R2 L1 di1
R1 R2
i1
dt 2
L1 L2 M 122 dt L1 L2 M 122
(2)
L2
du1
R2
u1 .
2
L1 L2 M 12 dt
L1 L2 M 122
d 2i1
dt
R1 L2 R2 L1 di1
R1 R2
i 0.
2 dt
2 1
L1 L2 M12
L1 L2 M12
(3)
1,2
(4)
(5)
(6)
A U L
max 2 .
B U max R2
R1 L2 L1 R2
R1 R2 L1 L2 M
2
2
12
c R1 R2 / L1 L2 M .
2
12
R1 R2 2 L1 L2 M 122
R1 L2 L1 R2
b R1 L2 R2 L1 / L1 L2 M 122 ,
(7)
(9)
R1 L2 R2 L1
2
L L M 12
1 1 2
2
2 RL R L
4 R1 R2
1 2
2 1
2
2
L1 L2 M 12
L1 L2 M 12
1 C1 A
.
2 C 2 B
/ L L
d L2 / L1 L2 M 122 ,
e R2
2
12
(8)
112
U m 0
Um
0
0
0
0
0
0
0
1
0 d A c 2 Bb
0 e Ab B c 2
1 C1
A
2 C 2
B
. (11)
0.15
0.1
0.05
Current [A]
-0.05
-0.1
-0.15
-0.2
-0.25
-0.3
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
0.09
0.1
Time [s]
(12)
0.06
0.05
0.04
300 points
0.03
0.02
0.01
-0.01
1000 points
-0.02
-0.03
-0.04
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
0.09
0.1
Time [s]
Parameter
1 [1/s]
2 [1/s]
A [A]
B [A]
C1 [A]
C2 [A]
1000 pt.
-18.1
-290
0.133
0.022
0.028
0.1
300 pt.
-18.95
-45.8
0.151
0.01
0.025
0.045
Real value
-15.2
-31.7
0.116
0.02
0.028
0.05
113
0.022 0.02
0.02
100 % 14.6 %,
0.116
6. Conclusions
100 % 10 %.
18.96 15.2
15.2
45.8 31.7
31.7
100% 24.7 %,
7. References
100 % 44.5%.
0.1
0.05
Current [A]
-0.05
-0.1
-0.15
-0.2
-0.25
-0.3
0.01
0.02
0.03
0.04
0.05
0.06
0.07
0.08
0.09
0.1
Time [s]
114
1. Introduction
In the context of this paper modelling of ammonia
production process are representations of chemical
process incorporating mass balance systems for key
components, and mechanistic sub-models for
modelling of absolute reaction rates.
Due to complex reaction systems in this process, like in
all biotechnological and chemical processes [1, 2, 3], it
is not possible to describe the process in details using
only mechanistic models. Thus it is necessary to use
methods that are able to fill some gaps in special
chemical or engineering knowledge. Further
expansions of the data and knowledge driven
approaches used in chemistry engineering are mass
balances systems. In the recently published literature,
some interesting examples show, that hybrid
combination of artificial neural networks, mechanistic
kinetics and mass balance equations can lead to
considerable advantages [1, 2, 3, 4]. Though, in
chemistry one must work with many data of off-line or
quasi-off-line types, and for this reason the above
mentioned data must be generated using interpolation
techniques. Regrettably, it may be challenged not only
by some artificial disturbances, but also by accuracy of
such interpolation methods. Moreover, in chemical
technology one is most often interested in reaction
rates, because they define the behaviour of the
conversions in chemical processes. Rates have to be
determined from the measured or estimated
concentrations by differentiation. In this respect noisy
(1)
CH4+O2 CO2+2H2+Q.
(2)
115
(3)
d (Ciout )
F
F
qi Ci Ciin in Ciout in ,
dt
V
V
(4)
2.4 Methanation
The small amounts of CO and CO2, remaining in the
synthesis gas, are poisonous for the ammonia synthesis
catalyst and must be removed by conversion to CH4 in
the methanator. The composition of the gas leaving the
methanator follows the chemical equilibriums (6) and
(7):
CO+ 3H2 CH4+H2O + 206.3 [kJ/mol],
(6)
(7)
(8)
k1
aA bB cC ,
116
(9)
k2
aA cC dD bB,
(10)
where a, b, c, and d are empirical coefficients of
chemical equilibriums; A, B, C, and D components of
chemical equilibriums; k1 and k2 constants of reaction
rates.
The absolute reaction rate for components in (9) and
(10) are:
rA ak1 Aa B b ak2 Aa C c ,
(11)
(12)
(13)
rD dk2 AaC c ,
(14)
(15)
E
b
) K exp(a ),
k k0 exp(
RT
T
nRT
,
p
(16)
1(in)
492.00
3115.6
17242
4.99e-5
0.0521
0.0917
2(in)
830.00
3066.5
60676
0.0281
0.027
0.196
0.0321
0.2042
-
CO2 conversion
Shift conversion
3(in)
350.45
2970.4
82376
0.0241
0.0416
0.1845
8.264e-4
0.1804
-
4(in)
199.63
2868.4
84032
0.0715
0.0133
0.279
1.067e-3
0.2055
-
5(in)
65.01
2711.5
60199
0.1684
1.93e-3
0.5819
2.31e-3
0.0108
-
6(in)
50.00
2611.5
77.92
0.2273
-
Methanation
7(in)
50.00
2611.5
151.25
0.2273
-
8(in)
295.00
2518.3
49120
1.59e-4
1.27e-3
0.3901
1.43e-3
2.97e-3
-
d (C H 2 )
F
rH 2 (C H 2 in 2 CH 2 ) 2 ,
dt
V1
d (CCH 4 )
F
rCH 4 (CCH 4in 2 CCH 4 ) 2 ,
dt
V1
(19)
d (CH 2O )
F
F
rH 2O (CH 2Oin1 CH 2O ) 1 (CH 2Oin 2 CH 2O ) 2 ,
dt
V1
V1
(20)
d (CO 2 )
F
rO 2 (CO 2in1 CO 2 ) 1 ,
dt
V1
(17)
117
(18)
(21)
rates
of
the
modelled
(22)
(23)
(24)
rH 2O k2CH 2CO0.52 ,
rO 2 k1CCH 4CO 2 0.5k2CH 2CO0.52 ,
d (CCO 2 )
F
rCO 2 (CCO 2in 3 CCO 2 ) 3 ,
dt
V2
(25)
(26)
(27)
d (CCO )
F
rCO (CCOin 3 CCO ) 3 ,
dt
V2
(28)
F
d (C H 2 )
rH 2 (C H 2in 3 C H 2 ) 3 ,
dt
V2
(29)
d (CH 2O )
F
rH 2O (CH 2Oin 3 CH 2O ) 3 ,
dt
V2
(30)
rCO 2 kCCO C H 2O ,
(31)
rCO kCcoCH 2O ,
(32)
rH 2 kCcoC H 2O ,
(33)
rH 2O kCCOCH 2O ,
(34)
118
rCO 2 k1CCO 2C H 2O ,
(45)
rH 2O k1CCO 2CH 2O ,
(46)
d (CCO 2 )
F
rCO 2 (CCO 2in 4 CCO 2 ) 4 ,
dt
V3
(35)
d (CCO )
F
rCO (CCOin 4 CCO ) 4 ,
dt
V3
(36)
d (C H 2 )
F
rH 2 (C H 2in 4 C H 2 ) 4 ,
dt
V3
d (CCO 2 )
F
rCO 2 (CCO 2 in 8 CCO 2 ) 8 ,
dt
V5
(47)
(37)
(48)
(38)
d (CCO )
F
rCO (CCOin8 CCO ) 8 ,
dt
V5
d (CH 2O )
F
rH 2O (CH 2Oin 4 CH 2O ) 4 ,
dt
V3
(39)
rCO kCcoCH 2O ,
(40)
rH 2 kCcoC H 2O ,
(41)
rH 2O kCCOCH 2O ,
(42)
F
F
(CCO 2 in 6 CCO 2 ) 6 (CCO 2in 7 CCO 2 ) 7 ,
V4
V4
d (CH 2O )
F
rH 2O (C H 2Oin5 C H 2O ) 5 ,
dt
V4
(49)
d (CCH 4 )
F
rCH 4 (CCH 4in8 CCH 4 ) 8 ,
dt
V5
(50)
d (CH 2O )
F
rH 2O (CH 2Oin8 CH 2O ) 8 ,
dt
V5
(51)
F
d (C H 2 )
rH 2 (C H 2in 8 C H 2 ) 8 ,
dt
V5
(43)
rCO 2 k2Cco 2C H4 2 ,
(52)
rCO k1CcoCH3 2 ,
(53)
(54)
(55)
(56)
(44)
119
5. Modelling results
The most important concentration components of
ammonia production process are shown in Fig. 3.
Concentrations of the components were chosen for
illustration based on the importance to the ammonia
production process and available measurement points
in the plant. The model was created taking into account
the data provided by ammonia technological rules of
procedure (see Table 1) and the modelling results were
compared with the reference values. The results, i. e.
8
6
4
CH4 conc.,[mol/l]
-3
x 10
time,[s]
CH4 concentrations after methanator
3
2
-4
0.16
0.14
0.12
0.1
3
4
5
time,[s]
CO concentration after low temp. CO converter
3
time,[s]
0.011
0.01
0.009
5
modelled
0
-4
9.7
2
x 10
9.75
0.012
CO2 conc.,[mol/l]
9.8
time,[s]
x 10
CO conc.,[mol/l]
x 10
CO conc.n,[mol/l]
CH4 conc.,[mol/l]
-4
10
3
4
time,[s]
CO2 concentration after CO2 convertor
3
2
1
0
measured
tech.reg.
time,[s]
2.
3.
4.
5.
7. References
1. Galvanauskas V., Georgieva P. and Feyo de
Azevedo S. Dynamic Optimisation of Industrial
Sugal Crystallization Process based on a Hybrid
(Mechanistic+ANN) Model. Proceedings of IEEE
World Congress on Computational Intelligence,
6.
7.
120
0.8
DOC, mg/l
0.7
Keywords:
mathematical
model,
identification,
dissolved oxygen concentration, biological wastewater
treatment process.
1.
DOC
Set point
0.9
0.6
0.5
0.4
0.3
0.2
0.1
0
Introduction
10
20
30
40
50
60
70
Time, h
2.
Technological process
121
ksed, %
100
50
0
Fig. 2. Technological scheme of wastewater treatment process
The relationship
polynomial:
is
approximated
by
ksed = 0.01034g2-2.034g+10.
quadratic
(1)
Fsed
k sed
Fin .
100
S bio
(2)
(3)
S sed S in
k dec
S in .
100
k
k
(6)
Fair 4 Fair
(4)
4
,
1 2 3 4
(7)
122
3.
X
Fbio ,
dt
Yx / s K s
V
max s
OUR
Yx / o K s
X,
k L a N Q ,
(8)
(9)
(10)
h-1
max
(11)
mg/l
mg/mg
mg/mg
mg/l
-
K
Yx/s
Yx/o
X0
k
Value
0.16
32.28
5.09
0.138
0.10
35.36
2.80
0.052
140
-0.0006
1020
Air flow
1000
980
m3/h
960
940
900
where
880
860
interval of observation.
For integration of the model equations (8), (9) along the
observation time interval, the discrete observation data
of technological parameters Fin , Fret , Q and c were
840
0.05
0.1
0.15
0.2
Time, h
0.25
0.3
0.35
0.4
Inlet flow
Returned flow
Total flow
900
800
700
sbio is estimated by
Total flow
m3/h
920
600
500
400
X X 0 Fret k ,
(13)
Inlet flow
300
200
100
Returned flow
0
0.05
0.1
0.15
0.2
Time, h
0.25
0.3
0.35
123
0.4
176
175
174
mg/l
173
172
6.
171
1.
170
169
0
0.05
0.1
0.15
0.2
Time, h
0.25
0.3
0.35
0.4
DOC, mg/l
0.25
0.2
0.15
0.1
0.05
0.05
0.1
0.15
0.2
Time, h
0.25
0.3
0.35
0.4
DOC, mg/l
0.16
Observed DOC
Model predicted DOC
0.14
0.12
0.1
0.08
0.05
0.1
0.15
0.2
0.25
Time, h
Fig. 7b. Model validation
0.3
0.35
0.4
DOC, mg/l
0.25
Observed DOC
Model predicted DOC
0.2
0.15
0.1
0.05
0.05
0.1
0.15
0.2
Time, h
0.25
0.3
0.35
0.4
References
Conclusions
124
2. Object Description
Two new settlers, methanol dosing station, blower room,
4 bioreactors, and distribution chamber with the pump
station were built in Kaunas waste water treatment plant
in 2008 [8].
Wastewater from the primary settlers is transported to
the distribution chamber. Here, wastewater is mixed
with the circulated sludge and dosed methanol. Sewage
and sludge mixture circulation is divided into 4 channels
for sewage and sludge mixture circulation of the
biological treatment tanks (bioreactors).
Kaunas waste water treatment plant has 4 biological
treatment tanks, where the mixture stirred 24 mixers
ensuring adequate mixing of sewage. Each mixer is
equipped in the bottom of the tanks. Effluent at the
origin of organic material will be removed biologically
active sediment phase, seeks to phosphorus, nitrogen
and carbon degradation.
125
W0 ( s ) =
0.0018
e 12600 s , mg/l/%
3.5962 103 s + 1
(1)
0.1
0.09
0.08
0.07
0.06
0.05
0.04
0.03
Experimental data
Modeling data
Filter data
0.02
0.01
0
10
12
14
Time, h
W0 ( s ) =
0.0029
e 12000 s , mg/l/%
4.9796 103 s + 1
(2)
0.2
Experimental data
Modeling data
Filter data
0.18
0.16
0.14
0.12
0.1
0.08
0.06
0.04
0.02
0
10
12
14
Time, h
0.35
Experimental data
Modeling data
Filter data
0.3
0.25
0.2
0.15
0.1
0.05
10
12
Time, h
126
14
0.965 T pr 0.855
( )
K pr pr
T pr
I=
0.0039
W0 ( s ) =
e 12240 s , mg/l/%
6.9547 103 s + 1
(3)
= 0.01179
pr
0.796 0.147
T pr
0.965 T pr 0.855
D=
( )
0.308 T pr
K pr pr
(7)
pr
T pr
0.929
(8)
= 5.527 105
The dissolved oxygen control system was realized in the
software SIMULINK (MATLAB).
The control systems with different regulators (PI and
PID) reaction to process interference (Fig. 5) and set
point change (Fig. 6) were presented and analyzed.
I=
pr
1.03 0.165
T pr
(4)
= 0.01179
(5)
0.24
0.22
I model
II model
III model
0.24
I model
II model
III model
0.25
0.18
0.16
0.24
PI
0.23
0.22
0.21
0.2
20
40
60
80
I model
II 120
model
III model
100
Time, h
0.19
0.18
0.17
0.16
20
40
60
80
100
120
Time, h
0.22
0.2
I model
II model
III model
0.18
0.16
20
40
60
80
100
120
Time, h
0.18
0.16
0.005
0.2
0.15
0.2
PID
0.22
0.26
0.26
0.24
0.586 T pr 0.916
( )
K pr pr
T pr
0.26
0.586 T pr 0.916
P=
( )
= 89.53
K pr pr
0.01
0.015
0.02
0.025
0.03
0.035
Time, h
P=
0.965 T pr 0.855
( )
= 152.6
K pr pr
(6)
127
Tss , g/l
7.99
8.44
8.72
0.0018
W0 ( s ) =
e 12600 s
3.5962 103 s + 1
0,0029
W0 ( s ) =
e 12000 s
4,9796 103 s + 1
0.0039
W0 ( s ) =
e 12240 s
6.9547 103 s + 1
(9)
Tss , g/l
106.00
88.66
89.53
0.0138
0.0111
0.0095
7.99
Tss , g/l
7.99
8.44
8.72
8.44
188.1
154.2
152.6
0.01547
0.01346
0.01179
D
6.501105
5.455105
5.527105
8.72
3. Adaptive Dissolved Oxygen Concentration
Control System in the software SIMULINK
(MATLAB)
Ko, mg/l/%
ko
7.99
Tss
Tssfcn
Differential
Equation
Editor
to
Object
Param eters
Process
Interference
Out1
Object
Set Point
0.15
P
K
ko
fcn
to
Transport
Delay
Scope
u
I
D
D
PID
Regulator
Parameters
128
0.19
0.18
0.22
0.21
0.2
0.24
Adaptive PI regulator
PI regulator
0.22
0.2
0.18
Adaptive PI regulator
PI regulator
0.16
Adaptive PI regulator
PI regulator
0.19
20
40
60
80
100
120
Time, h
0.18
0.17
0.16
0.15
0.17
20
40
60
80
100
120
0.21
0.2
0.19
0.18
0.22
0.21
0.2
0.19
0.18
0.24
0.22
0.2
T =8.72 g/l
0.18
0.16
20
Adaptive0 PID regulator
PID regulator
0.17
40
60
80
100
120
Time, h
Tss=8.44 g/l
0.16
0.15
20
40
0.17
60
80
100
120
Time, h
Tss=7.99 g/l
0.16
0.15
Time, h
0.23
0.26
20
40
60
80
100
120
Time, h
0.16
20
40
60
80
100
120
Time, h
0.26
0.25
0.24
0.21
0.2
0.19
0.18
0.17
0.16
0.15
0.25
0.24
Adaptive PI regulator
PI regulator
0.23
0.25
0.22
0.21
0.2
0.19
0.18
20
0.17
40
0.16
0.15
0.22
0.23
0.24
0.22
0.2
0.18
0.16
0.24
0
0.23
0.22
Adaptive PI regulator
PI regulator
80
100
120regulator
Adaptive
PI
0.2
Time, h
PI regulator
0.21
60
40
60
80
100
0.23
0.22
0.25
0.24
0.21
0.2
0.19
20
0.18
0.17
0.16
0.15
0.19
0.18
20
0.17
0.25
120
0.15
0.23
0.2
0.24
0.26
0.21
0.24
0.22
40
0.21
60
80
100
120
Time, h
0.2
0.19
0.18
20
40
0.17
0.16
0.15
Time, h
0.23
20
40
60
80
100
120
Time, h
0.16
0.15
20
40
60
80
100
120
Time, h
129
0.2
0.19
0.18
0.21
0.24
0.23
0.22
0.21
0.2
Tss=8.72 g/l
0.24
0.22
Adaptive PI regulator
Adaptive PI regulator
Adaptive PIDAdaptive
regulatorPID regulator
0.2
0.18
0.16
Adaptive
PI regulator
Adaptive PID regulator
0.19
0.18
20
40
60
80
100
120
Time, h
0.17
Tss=8.44 g/l
0.16
0.15
0.17
20
40
60
80
100
120
Time, h
5. Acknowledgment
Tss=7.99 g/l
0.16
0.15
0.26
20
40
60
80
100
120
Time, h
0.24
0.2
0.18
0.16
0.25
0.24
0.23
0.22
0.21
0.2
0.19
20
0.18
0.17
0.16
0.15
0.22
6. References
0.26
0.25
Adaptive PI regulator
Adaptive PID regulator
0.24
0.23
0.22
60 Adaptive
80 PI regulator
100
120
PID regulator
Time,Adaptive
h
0.21
40
0.2
0.19
0.18
20
0.17
40
60
Adaptive PI regulator
Adaptive
PID regulator
80
100
120
Time, h
0.16
0.15
20
40
60
80
100
120
Time, h
4. Conclusions
130
H ti H sti H i ,
(1)
H sti H ni H sti ,
(2)
The equations (1) and (2) define the required head at the
station's suction branch. Depending on the value of head
at the suction branch, the head of the pump set Hdi is
expressed:
H di H ti H s ,
(3)
where: Hs inlet water head.
131
H
n=nN
Htrma
Htri
Htmax
Hstma
Htrmin
Hti
Htri -
(Qi, Hti)
Hsti
Hstmi
Qmin
Qi
Qmax
zi fi si pi .
132
(5)
Q1
(7)
C fi si .
Qm
.
Q zi
(8)
n zi
nN
i 1
zi
2 a 2 Q b2 Q
2
for
Q 0.5Q N
for
Q 0.5Q N .
(14)
Q
1 zi
Q1
0, 09
(15)
1 a1Q 2 b1Q c1
0, 09
pi 1
H zad
A 2 4 H 0 B
Q zi2
H zad
2 B
Q zi2
c f Q z1 , Q z 2 ...Q zm .
(16)
(9)
Hti [m]
Qm [m3/h]
Hti [m]
Qm [m3/h]
Hti [m]
Qm [m3/h]
(11)
(12)
70
85
70
80
70
75
70
70
70
65
70
60
70
50
70
40
70
30
65
90
65
85
65
80
65
75
65
70
65
65
65
60
65
50
65
40
65
30
60
90
60
85
60
80
60
75
60
70
60
65
60
60
60
50
60
40
60
30
70
90
(13)
133
n zi n z1 n z 2 n z 3 .
(17)
(19)
5. Conclusion
cmax
6. References
cmax
134
Experience data
Real
system
Software
Reality
modelling
Simulation
Model
characteristic package;
neuron models;
fuzzy models;
etc.
Real object
Definition of essentials
and inessential for
analyze object features
Separation
compact
subobjects
Ideological schema
Model
description
Analyses of sub-objects
determination of type
of model
Simulation model
Fig. 2. Periods of development of simulation model
135
they describe
concrete
object
(phenomenon).
Fuzzy models have character descriptivelyexpert.
They do not require acquaintance of parameter of object
directly but experimental acquaintance of its behaviour
is indispensable.
It is possible to divide process of modelling into several
periods (Fig. 2).
Firstly, the ideological scheme of modelled object
should be constructed, provided that all the particles
have been previously distinguished. It will be necessary
to carry analysis of the problem what kind of simulation
model will be proper for description of chosen element.
In final period it is necessary to describe simulation
models of individual elements and define manner of
exchange of data among them.
a)
b)
1 1 0 s A (t)
u AB (t)
u (t) U 0 1 1 s (t) ,
d
B
BC
1 0 1 s C (t)
u CA (t)
(1)
where U d 2U ac .
These switching functions were used to control
transistors in simulation model of PWM rectifier.
136
5. Pump characteristic
dI
1
A( U RI ) ,
dt K
(2)
where:
K = r + s + m;
m inverse of magnetising inductance;
s, r inverse of stator and rotor inductances;
s r
0
0
s r m
s r ;
0
s r m
0
A
s r
0
r s m
0
s r
r
s
m
U u sa
u sb
0 0 T ;
usa = Um sin(2ft + );
usb = Um sin(2f t + + 120o);
I i sa i sb ira i rb T
1
1
3 2rsa rsb 3 rsc rsb
1
1
2rsb rsc
R rsc rsa
3
3
0
0
0
0
0 0 ;
rr 0
0 rr
0
m
r
;
i sb irb irb
r
m
0
0
0
0
2
.
3
3
2
3
3
7. Algorithm of calculations
Te T ,
dt
J
Te
(3)
(4)
where:
T() mechanical moment;
p number of pole pairs;
J moment of inertia;
Te electromagnetic moment.
137
8. Simulation results
9. Conclusions
10. References
1.
2.
3.
Fig. 8. Results of simulation test motor speed
4.
5.
6.
7.
8.
138
139
Input variables
Output variables
140
3. Simulator
A new GrowSim simulation tool for drum granulationdrying process is under development. The simulator is
based on the extended process model and is intended to
be used for advanced process control and knowledge
acquisition. Using powerful MATLAB software the tool
simulates the physical process combining mass and heat
transfer, including major units such as a drum-granulator
dryer, classifier, crusher, transportation system. The
graphical user interface (Fig. 3) has been built to mimic
the process control environment available to the process
operator in the real plant, with important additional
information provided.
The simulation environment is composed of sections
where the operator can change process parameters,
observe the current or past states, get some advice on
how the process should be controlled by the skilled
operator. The simulator allows to introduce common
4. Results
The most common description of particles in means of
their size is given by particle size distribution (PSD).
Strong evidence of normal granule size distribution can
be observed in Fig. 6, where the measured initial nuclei
and granules out of granulator distribution data is fitted
to the cumulative normal distribution curve with
parameters defined.
Some simulation scenarios have been carried out to
check the model and simulator using the same data as in
Fig. 4. The results obtained are in fair accordance.
141
1
Measured distribution (initial seeds)
Measured distribution (granules from granulator)
Fitted Normal distribution (=2.95 mm; =0.79 mm)
Fitted Normal distribution (=2.75 mm; =0.67 mm)
0.9
0.8
0.7
F(x)
0.6
0.5
0.4
0.3
0.2
6. Conclusions
0.1
0
0
Granule Diameter, mm
Percentage
Simulated
Measured
0.8
0.6
0.4
0.2
0
0
0.4
Percentage
Particle Size, mm
Histogram
0.3
0.2
7. Acknowledgment
0.1
0
0
0.5
1.5
2.5
3.5
4.5
Particle Size, mm
8. References
1
0.8
0.6
0.4
0.2
0
0
Particle Size, mm
Histogram
0.4
0.3
0.2
0.1
0
0
Particle Size, mm
Granule Mass
0.9
d0.5-2, %
0.8
Normalized Parameters
Flow,out
d2-4, %
d4-8, %
0.7
0.6
0.5
0.4
0.3
0.2
0.1
0
0
10
15
20
25
30
35
40
45
50
Time, min
142
Introduction
2.
143
Relative moisture
Target moisture
value
5
4
3
3
Mixing phase
1
0
0
20
40
60
80
t,s
100
120
140
160
180
120
90
60
30
Material
number
1
2
3
4
5
Added water, l
Mixing time, s
150
6 7 8 9 10 11 12 13 14 15 16
Experiment number
Transmission
line under test
Signal
reading
circuit
100
75
50
25
0
Sand
0-2 mm, %
ANN
144
As it is shown in
Fig. 5, the possibility to identify the type of material
was found. Most of the samples were recognized. Only
a few samples were assigned to wrong material type. It
shows the possibility to use this method in real time
moisture control systems.
5.
Learning
result
-1
-0.5
0.5
4.
6. References
1.
2.
Results
3.
4.
M5
Material
M4
M3
M2
M1
14
21
28
Conclusions
35
Sample Nr
145
146
dimensional process.
4. Fuzzy control
The results of static tuned PID control of two
dimensional process showed that this control method
was not suitable for the plant. Consequently, we needed
another control method which could control the plant in
suitable manner. One of possible control techniques to
air pressure plant is a fuzzy control.
147
1
N
i, j
j t N 1
(2)
(1)
148
i 1
( y)
( y)
(3)
i 1
149
5. Experimental results
The experiment was performed with different sets of
reference signals, but the main control challenges
existed when the reference signals were changed
significantly, especially when they were reduced.
Therefore the experiment was concentrated to those
problems. In the section about the structure of control
system was mentioned that the air channels dependence
of each other was decided by the number of opened air
feeding vents. The experiment was performed with such
a fixed number of opened air feeding vents that ensured
the necessary amount of air to channels fans since this
configuration of opened air feeding vents guaranteed a
certain amount of interactions between air channels and
the corks were able to reach the reference signal values
physically.
In Fig. 8 and Fig. 9 are demonstrated the results of PID
and fuzzy controls, respectively, when the reference
signals for both air channels are 0, then changes to 75
and later to 60, the abbreviation of such a reference
signal configuration is marked as 0->75->60.
1
N
( y(t ) x(t ))
(4)
t 1
PID <1>
A
PID <2>
0 -> 75
50
28.771
50
4.48
26
50
1.095
28.263
75 -> 60
31
50
1.13
4.11
30
0.549
38
0.44
Fuzzy <1>
A
Fuzzy <2>
0 -> 75
50
21.535
50
3.025
32
50
1.084
20.938
75 -> 60
43
50
0.991
2.731
33
1.589
43
0.782
150
PID <1>
A
PID <2>
0 -> 75
50
23.975
50
8.056
16
50
0.299
23.401
75 -> 40
50
6.435
21
1.041
28
1.652
Fuzzy <1>
A
Fuzzy <2>
0 -> 75
50
21.395
50
9.527
35
50
0.72
19.91
75 -> 40
32
50
0.954
10.044
28
0.438
35
0.636
6. Conclusions
The analysis of experimental results proves that fuzzy
control is more effective control method to air pressure
plant than PID. The implementation of humans
heuristic information about how to control a system,
enables to estimate the critical system control problems
(the dynamic and nonlinearities of system, interactions
between the channels) better than PID which directly
guides to better system control. The crucial factors of
successful plants control are the estimation of corks
moving directions and proportions. The synthesis of
adaptive PID control could be one of possible
replacements to static tuned PID in order to control the
plant with PID-based control.
7. References
1. lstrm K.J., Hgglund T.H. PID controllers:
Theory, Design and Tuning. North Carolina, 1995.
343 p.
2. Passino K., Yurkovich S. Fuzzy Control. California,
1998. p. 9-10.
3. Kovaic Z., Bogdan S. Fuzzy Controller Design
Theory and Applications. Portland, 2005. 416 p.
4. Abonyi J. Fuzzy Model Identification for Control.
Boston, 2003. 273 p..
5. Albertos P., Sala A. Multivariable Control Systems.
London, 2004. p. 35-38.
151
dCH
(1 x 2 x)
dt
FpH (CH0 CH ) Fs CH
,
V
V
(1)
where:
CH and C H0 concentration of hydrogen-ions in the
cultivation medium and in the feeding solution,
respectively, [mol/l]; the latter concentration can differ
from the theoretically calculated and is subjected to
model based identification;
x biomass concentration in the cultivation medium,
[g/l];
specific biomass growth rate, [1/h];
FpH flow of the alkali solution for pH control, [l/h];
Fs flow of the feeding solution, [l/h];
V cultivation medium volume, [l];
1 , 2 model parameters to be identified.
152
0,
(4)
(1 2 ) xV *
.
(CH0 CH )
(5)
dt
FpH
level corresponds to pH 7.
The biomass growth in the fed-batch process can be
modelled by means of the differential equation
F FpH
dx
x s
x.
dt
V
(2)
(3)
FpH
V
(CH )
CH0 CH
V
FpH .
(6)
CH0 CH
V d ( CH )
C
)
FpH .
H
FpH
dt
FpH
(7)
W ( s)
C H ( s )
FpH ( s )
K0
,
T0 s 1
(8)
K0
CH0 CH
FpH
T0
(CH0 CH ) 2
(1 2 ) xV *
(CH0 CH )
V
.
FpH (1 2 ) x
(9)
(10)
Ti
1
,
2 OUR
(11)
PID discr
f(u)
Error
Fcn1
pH SP
kr
Ti
kr
Un
Td
Td
Saturation Fb
Tdiscr
Tdiscr
mu
Tsp
Ti
fcn
Fb
Fb2 offline
w
Ti const / adapt
Fb1 online
Fb offline
x
1
1.e-7
20/3600s+1
pH sensor pH sensor T ransfer Fcn
Delay
C(H+)(0)
t
C(H+)
mu
fcn
-log10(u(1))
Fs
dx/dt
Saturation C(H+)
C(H+)->pH
pH
Fcn3
-6.e-5
f(u)
Fs1
OH-
pH err %
Fs
6.e-5
Switch
Base-Acid
Process x, C(H+), V
Fcn4
f(u)
H+
C(H+) err %
10
Clock
mu set
t
Ti adapt
fcn
mu
rH
OUR fcn
Fcn2
mu
f(u)
T ransport
Delay OUR
OUR1
H+ production rate
OUR
due
to
exponentially
growing
biomass
concentration that has acted as a disturbance. The
simulation results (see Fig. 3) show that the
tracking error of the adaptive PI controller was
significantly lower than the one in the case of a PI
controller with constant tuning parameters.
At the end of the process the adaptive controller
still was able to keep pH value of 6.999 as
compared to 6.96 in the case of the control system
based on the standard PI controller.
5. Conclusions
The proposed adaptive control algorithm proved to
be a robust and stable solution ensuring high
quality pH control in the analyzed biotechnological
fed-batch process. It can become a real alternative
for standard PID control. Further practical
investigations should be carried out in order to
identify model parameters from the experimental
data and to practically implement the proposed
adaptive approach into a real control system.
6. References
1. Sonnleitner B.
Instrumentation
of
Biotechnological Processes. Advances in
Biochemical Engineering/ Biotechnology, 2000,
66. p. 1-64.
2. Carr-Brion K.G. (Ed.) Measurement and
Control in Bioprocessing. Elsevier Science
Publishers, London, 1991. 262 p.
3. McMillan G.K., Cameron R.A. Advanced pH
Measurement and Control. ISA, 2005. 294 p.
4. Ljung L. System Identification: Theory for the
User, 2nd Edition. Prentice Hall, New Jersey,
1999. 672 p.
5. Levine W.S. (Ed.) The Control Handbook.
IEEE/CRC Press, 1996. 1566 p.
155
Abstract: Application of software agents and multiagents systems for decision making are the new
approach to the capital markets. The flow of
information is increasing rapidly and is extremely
important to evaluate the news on time and make the
right decisions. The papers propose the basic structure
of the multi-agent system that will function in real time
dynamic environment in capital markets. This IT system
will be able to suggest investment solutions taking into
consideration variable characteristics of individual
investors. IT system will acquire information from and
monitor changes to stock reporting databases, interpret
stock information, and predict the near future of an
investment. Broad spectrum of securities valuation
methods will be employed and profile of the user will
be incorporated in the decision making process. A
multi-agent system approach is ideal in cases that
require the coordination of several component tasks
across distributed and flexible networks of information
sources and expertise. The paper will present the basic
structure of Securities valuation IT system using multiagents approach.
Keywords: software agents, multi-agents systems,
securities valuation.
1. Introduction
Deep knowledge of fundamental, technical or other
securities valuation methodology is not enough for
successful trading in financial markets. Since trading
moved to electronic space, new challenges have
emerged. First, financial markets have become global.
On one hand, the investing in any region of the world
has become very simply. On the other hand, there is a
need to know the specifics of the region adapt to them
and learn to work with. Second, the financial markets
have become dynamic. The abundance of information
could quickly change the direction of the markets.
Therefore, there is a need to continuously track and
interpret the information and make decisions based on
156
2. Agents properties
The software agents in our IT system will possess the
following basic properties [7, 8]:
Autonomy: the ability of the agent to act
without any external interference. The agent
must be able to make its own choices and make
its own decisions.
Communication skills: the agent must be able
to establish relation with environment. The
most important is the interaction between
similar agents and the user.
Reactivity: surrounding environment is
typically dynamic and requires the adaptation
ability from the agent.
Mental notions: during his life cycle, the agent
is learning from different situations. Thus,
through the memories and through the
interaction with others agents, the agent is able
to gain its own knowledge.
Persistence: the mobile agents life is usually
longer than the length of his basic tasks. Thus,
it is created to satisfy one of the following
criteria:
1. To finish after having done the basic task;
2. To finish after having exhausted the
assigned internal resources.
Social ability: the agent exists among the other
agents, thus it must be able to exchange the
information in order to achieve its assigned
goals or help the other agents.
Mobility: a mobile agent has the inner attitude
to vary its communication partners. It is able to
interact with both a similar agent and a user.
Vitality: in order to achieve its goal, the agent
must cope with anomalous situation that could
create a state of instability that could damage it
irreparably. An agent who has a strong vitality
is able to solve the most adverse situations it
meets.
Pro-activity: a mobile agent is able to start new
interactions with other agents and coordinate
various agents activities stimulating their
certain responses.
Truthfulness: the agents use the information
received during the exchange and not viewed
157
3. Types of agents
The securities evaluation IT system will consist of 3
types of agents:
o Interface agent,
o Information agent,
o Co-operating agent (due the functional clarity it
will be named the Task agent in the further
article.)
The system could have more agents, too [10, 11].
The Interface agent is the mediator between the user
(investor) and the securities evaluation IT system [17,
18]. Interface agent receives the information from the
user and accepts its target, then transmits it to other
system elements. On the contrary, it presents the
systems generated answers to the user in an
understandable form.
The main functions of the Interface agent would be the
following: to collect and analyze the information about
the user and share it with other IT system agents, to
submit the system generated responses to the user, to
request for additional information and if it is necessary
to inform the user. Therefore, the IT systems user
interacts with interface agent only, no other agents are
involved.
The main function of the Information agent is to provide
the information to the Task agents [10, 15]. For
example, in order to learn the fundamental value of
securities, the information agent needs to gather the
financial data about the issuer using the Internet data
search. Then the Information agent transmits the
information to the Task agent, which is responsible for
fundamental analysis. The Information agents could
vary and perform different information gathering
strategies, use different sources of information.
Information agents distinguish their selves with proactivity feature. They are able to generate events in the
surrounding environment and to start new interactions
with other agents instead of waiting for assignment of
the data search.
Task agents perform a lot of work in shaping the
recommendation to the user. Each Task agent has
certain knowledge about certain security evaluation
approach. Depending on users request, which is
expressed through the interface agent, IT system may
take different methods and their combinations. When
performing the tasks, the agents are cooperating with
others thus communication skills and social ability is
the main property. Task agents are more complex than
the interface and information agents.
The security evaluation IT system obtains more
possibilities when agents are able to use the mental
notions. It is likely to have a more efficient IT system
when it is able to accumulate the history of
recommendations and to select the best decisions [18].
158
159
160
161
2. Data analysis
1. Introduction
Automatic taller machines (ATM) are equipment for
money supply for consumers. More than one thousand
ATMs are in Lithuania and over one billion in whole of
the world. So it is not a difficult to imagine what
amount of cash flow is needed for loading and ensuring
of ATMs correct action. At present the analysis and
forecast of money demand is made by human expert. He
is analyzing money flow with his own empirical model
witch estimate input parameters like money demand few
days ago, how much days till holidays and others.
Depending on strategies the ATM loading is executed
daily or weekly by fixed size of money load.
The main objective is to make load of ATM exactly
what it needs to be because money has a cost depending
on interest rate.
Studies of ATM turnover simulation and prediction
made in recent years were stimulated by private
companies and were unpublished. The main attention in
current publications in ATM field is concentrated to
transaction security. Recently investigation in ATM
turnover modelling field was performed in Kaunas
University of Technology [1, 2]. The main goal of the
research was to optimize ATM cash upload. It was
formed and identified an artificial neural network model
for each ATM and forecasted demand of cash for
specified period of time. These articles focus on an
optimization routine and wasnt briefly discussed
significant variables, ANN model structure and
parameters identification. This article tries to eliminate
these deficiencies.
Fig. 1. Real and mean cash flow of the all ATMs dependence
on week day number
162
Fig. 2. Real and mean cash flow of the all ATMs dependence
on month
TFout
2
1.
1 e2in
(1)
Inputs
[D M MD M7D]
[D M MD M5D]
[D M MD C1]
[D M MD C1 C2]
[D M MD]
[D M MD H]
[D M MD M5D C1]
[D M MD M5D C1 C2]
[D M MD C1 C2 C3]
3. Model identification
Artificial neural networks were selected as proper
technique for nonlinear modelling. Neural networks
were successfully used in many areas such as pattern
recognition, classification, nonlinear modelling and
others [3-7]. It is also convenient to use ANN when
necessary to adapt the model settings or select an
appropriate model structure. Requirements for model:
Mean
prediction
error of 5
learning
restarts,
MSE
2.1671*105
2.7676*105
6.7330*105
3.5555*105
2.3275*105
2.6056*105
2.3444*105
2.6756*105
2.4678*105
Best of 5
learning
restarts,
MSE
1.8137*105
1.8214*105
1.8473*105
1.9098*105
1.7813*105
2.1436*105
1.8689*105
1.9264*105
1.9393*105
163
Nave prediction
hypothesis
9.5504*104
Prediction
error
ANN model
prediction
5.8679*104
x 10
2.2
1.8
1.6
S u m a b s o lu te e rro r
1.4
1.2
0.8
0.6
0.4
0.2
630
640
650
660
670
680
690
700
710
720
730
ATM number
y t 1 yt ,
(2)
164
8. References
1. Simutis R., Dilijonas D., Bastina L. Cash demand
forecasting for ATM using neural networks and
support vector regression algorithms. In: Selected
papers of EurOPT-2008: the 20th International
Conference "Euro Mini Conference on Continuous
Optimization and Knowledge-Based Technologies",
Vilnius, Technika, 2008. p. 416421.
2. Simutis R., Dilijonas D., Bastina L., Friman J.,
Drobinov P. Optimization of cash management for
ATM network. Information technology and control,
2007, vol. 36, p. 117-121.
3. Gupta M.M., Jin L., Homma N. Static and dynamic
neural networks: from fundamentals to advanced
theory. John Wiley & Sons, 2003. 700 p.
4. Gnoth S., Jenzsch M., Simutis R., Lbbert A.,
Control of cultivation processes for recombinant
protein production: a review. Bioprocess and
Biosystems Engineering, 2008, vol. 31, no. 1. p. 2139.
5. Simutis R. Exploratory analysis of biochemical
processes using hybrid modeling methods,
Proceedings of the First European Symposium
Principles of Data Mining and Knowledge
Discovery Trondheim, Norway, June 24-27, 1997.
p. 200-210.
6. Wei H., Shouyang W., Lean Y., Yukun B., Lin W. A
new computational method of input selection for
stock market forecasting with neural networks.
Proceedings of International Conference ICCS 2006,
Part IV, 2006. p. 308-315.
7. Huang W., Lai K.K., Nakamori Y., Wang S.Y.
Forecasting foreign exchange rates with artificial
neural networks: a review. International Journal of
Information Technology & Decision Making, 3(1),
2004. p.145-165.
7. Conclusions
Computational method for ATM cash flow forecasting
using artificial neural network was proposed. Influence
of various input sets to model accuracy was
investigated. The best input set was used for model
identification. ANN model prediction efficiency was
estimated by comparison with Nave prediction
hypothesis. The error of ANN model was about 39%
smaller. The result is hopeful, but further experimental
investigation about practical implementation and more
accurate cash turnover prediction is necessary.
Acknowledgments
This research work has been supported by means of EU
structural funds (Project No. BPD2004-ERPF-3.17-0606/0045).
165
166
2008
1917
2010
2144
2015
2471
4599
3466
3466
1183
2425
1722
344
161
7
148
2425
1722
344
161
7
148
9
2425
1722
344
161
7
148
9
34
34
34
73
73
73
887
887
837
29
52
2384
287
1183
102
52
29
52
1498
287
289
110
52
760
29
52
1519
287
289
131
52
760
760
1968
1947
2214
-176
-525
Fuel
Lt/t
Lt/ toe
Price of
import,
production
1130
1183
Price
at
Lithuanian
power
plants
1191
1247
Lt/t
Lt/ toe
Lt/th.m3
Lt/ toe
Lt/m3
Lt/toe
Lt/t
Lt/ toe
Lt/t
Lt/ toe
2105
2151
1000
1250
140
837
1100
1152
1200
1256
2145
2099
1080
1350
Unit
1161
1216
1261
1320
168
1400 ktoe of natural gas and 320 ktoe of fuel oil. Since
2012 a great role will fall on new combined cycle gas
turbines at Lithuanian and Kaunas power plants.
Therefore natural gas demand at Lithuanian TPP will
decrease till 760-860 ktoe in 2015, from which 430 ktoe
will fall on the new unit of the Lithuanian TPP that will
produce 2.8 TWh of electricity.
If natural gas supply through pipeline Minsk-Vilnius
will be lost, lack of natural gas could be displaced by
supply from Latvian underground gas storage. In
extreme case when electricity import is not possible and
gas supply Minsk-Vilnius is interrupted for two winter
months Lithuanian power plants from Latvian gas
storage will require 300 ktoe in 2010 and 380 ktoe in
2015. One third of this gas (130 ktoe in 2010 and 112
ktoe in 2015) will be consumed by units of Lithuanian
TPP, that have no flue gas desulphurisation equipments.
In this a case electricity to Lithuanian consumers could
be provided without any limits and there were no
environmental violations on SO2 emissions. If natural
gas supply from Latvian gas storage to Lithuanian
power plants werent higher than 120 million m3 (about
100 ktoe), then fuel oil and other liquid fuel would be
fired and this would lead to environmental violations on
SO2 emissions or to limitation of electricity
consumption.
In 2010 nuclear fuel in Lithuania will be changed to
expensive fossil fuel. Mainly because of increased
expenditures for fuel procurement, the costs of power
plants operation will greatly increased 2.6-3.1 times
and will make about 3.0-3.7 billion Lt. This would have
an impact on electricity price, which will increase.
Average annual electricity production price in
Lithuanian power system in 2010-2011 will increase till
31 Lct/kWh (or 89.8 Euro/MWh). When new combined
cycle gas turbine power units start operation, the
average annual electricity production price would
decrease till 26-27 Lct/kWh (or 75.3-78.2 Euro/MWh).
However this price would be 3 times higher than the
price of Ignalina NPP. This electricity production price
will be much higher than electricity price in
neighbouring countries. For example, it is forecasted
that electricity production price in Scandinavian market
will range from 30-45 Euro/MWh in 2010-2015, e.g. it
would be 2 times lower than in Lithuania. In the
beginning of February 2009 the average daily electricity
production price in Scandinavian market was
approximately 40.8 Euro MWh, in Poland and Germany
correspondingly 44 Euro MWh and 50.7 Euro/MWh.
Contrary to these prices the production price of
electricity in Europe-Ural price zone in Russia was
equal to only 12.6 Euro/MWh in January 2009.
The presented information about electricity prices in
neighbouring markets illustrates the importance of
connection of Lithuanian power system with Poland and
Sweden systems. After connection of these systems,
Lithuania will have ability to choose electricity
suppliers and therefore will avoid price intrusion, which
is substantial if there is the only link with one country.
It is worth to note that electricity price to final
consumers is not equal to the price of electricity
production. Electricity transmission, distribution, supply
169
170
K j = K0 l j ,
(1)
i
K j .
100
(2)
sea
CL
and
p = 2% .
Then constant operational costs are, /year:
CEk 1 =
100
( pa + p ) .
(3)
j =1
STOCKHOLM
Pk j =
Pk 0 max + Pk 0 min
2
lj ,
(4)
171
P
3 U n cos
demand current, A.
(5)
P3 f = I 2 R0 l (3 + 11 + 12 + 1V ) + 3 Piz , (14)
(6)
(7)
)]
R = R0 1 + 20 p 20 ,
PKL = Pvad + P3 f .
(8)
1000
conductor resistance at 20C;
where R0 =
F
conductor
specific
conductivity
Cu = 32 m mm 2 ;
20 = 0.00393 1 C temperature coefficient for
copper;
p = 90C maximum allowed temperature for
xs2 =
8 f ks
10 7 ,
R
CEm L j = PL j ( + ) .
(9)
(10)
CE j = CEkj + CEmL j .
;
2
x
1.18
d dz
d dz
yt =
0.312
+
4
4
s
192 + 0.8 xt s
8 f kt
xt2 =
10 7 ,
R
(18)
f = 50 Hz frequency.
Closeness effect factor:
2
(17)
4
t
(16)
(15)
C j = CK j + CE j .
(19)
(11)
C j = C j +
(12)
i
KP j ,
100
(20)
(13)
172
C 1 ,
Kd
CEmL1,
CE1,
C1,
1
1.015
1.03
1.045
1.06
962715
990036
1017765
1045901
1074443
2628715
2656036
2683765
2711901
2740443
6793715
6821036
6848765
6876901
6905443
6943715
6971036
6998765
7026901
7055443
Kd
Kd
1
1.015
1.03
1.045
1.06
1
1.015
1.03
1.045
1.06
C 2 ,
CEmL2,
CE2,
C2,
121351
124795
128290
131836
135434
331351
334795
338290
341836
345434
856351
859795
863290
866836
870434
3356351
3359795
3363290
3366836
3370434
AR 1 ,
AR 2 ,
AR 3 ,
MWh
1500000
1522500
1545338
1568518
1592045
MWh
1160.25
1177.65
1195.32
1213.25
1231.45
MWh
146.25
148.89
151.27
153.70
156.16
Kd
1
1.015
1.03
1.045
1.06
CEmL3,
356249
360516
364845
369239
373696
S j = a l j ,
C 3 ,
CE3,
15356249
15360516
15364845
15369239
15373696
45356249
45360516
45364845
45369239
45373696
Y j = S j cY ,
CKY j =
CY j = S j y ,
d j =
(22)
Cec j = C KY j + CY j .
.
(28)
(29)
(23)
Assuming that annual cost increases by 5% due to
inflation there were gotten annual damage expenses for
strategies 1 , 2 and 3 , accordingly, Cec 1 , Cec 2
AR j = d j Pmax Tmax .
(27)
d j Td , L
i
Yj .
100
(21)
(26)
d j = d ,sp l j .
(25)
(24)
173
and 3
Year
Kd
Kd
0.15309
0.07400
1.015
0.15368
0.07407
1.03
0.15428
0.07414
Cec 1 ,
Cec 2 ,
Cec 3 ,
1552500
180000
1578375
183000
1.045
0.15488
0.07421
0.15550
0.07428
12316500
1.015
12519722
R1 C , S
1.03
12712681
1604250
186000
1.06
1.045
12911798
1630125
189000
0.12605
0.27925
1.06
13106809
1656000
192000
1.015
0.12643
0.28009
R2
~
AR , S
1.03
0.12655
0.28037
1.045
0.12668
0.28065
1.06
0.12681
0.28093
0.12605
0.01461
R3 Cec , S
S
Kp
0
1
0.28
1.015
0.51
1.03
0.15
1.045
0.06
1.06
0.12607
0.01462
0.12619
0.01463
1.045
0.12625
0.01464
1.06
0.12635
0.01465
0.15428
0.07414
0.51
AR
0.51
0.12655
0.28037
Cec
0.51
0.12619
0.01463
aij is
1.03
1.015
normalized
values
using
the
expression
R (K , S ) = K (S )i K (S )max in the table 7, where K
corresponding partial criterion denomination [8].
Indirect target possessive functions (table 8) were gotten
using expressions:
C = max[min ( R~1 , S~ )] ;
Ar = max[min ( R~2 , S~ )] ;
174
3.
7. Conclusions
4.
5.
8. References
1.
2.
//
:
,
,
.
III
-
13-16 2008. 2, III
-
, Krievija, Ekaterinburga, 13. 16.
oktobris,
2008.
274277.
lpp.
6.
7.
8.
9.
175
..
//
. 2003. 3 14 22 .
.. . ,
, 1975. 320 .
176
Actions
Calculates transmission capacity limits,
losses and forecasts consumption
Receives information on the volumes of
fixed transactions. Receives regulation
power bids
Monitors national balance, regulates if
needed and registers regulating energy
transactions
Collects measurement data and balance
information
Calculates the imbalance for each hour
177
178
Date
May 18, 2000
July 19, 2000
December 18,
2001
January 1, 2002
April 18, 2003
June 1, 2004
179
180
7.
8.
9.
10.
11.
12.
6. References
1.
2.
3.
4.
5.
6.
13.
14.
15.
16.
181
Baltrel
annual
report
2007.
http://www.baltrel.org/annual_report_2007.pdf
Nord Pool Spot Baltic project No.97/2008.
http://www.nordpoolspot.com/en/Market_Informat
ion/Exchange-information/No972008-Nord-PoolSpot-Baltic-project-ends-31-December-2008/
No.55/2008 Introduction of implicit auction via
the Estlink cable by 1 July 2009 at the earliest .
http://www.nordpoolspot.com/en/Market_Informat
ion/Exchange-information/No522008Introduction-of-implicit-auction-via-the-Estlinkcable-by-1-July-2009-at-the-earliest/
No.27/2007 Joint Baltic-Nordic project launched.
http://www.nordpoolspot.com/en/Market_Informat
ion/Press-releases-list/No272007-Joint-BalticNordic-project-launched/
The Nordic Power Exchange and the Nordic
Model for a Liberalised Power Market
http://www.nordpoolspot.com/upload/Nordic%20p
ower%20market/Nordic%20power%20market.pdf
Sweden grid company official internet site.
www.svk.se
Electrical Power System Code. The Act of the
Government of Estonian Republic No.184,June.
http://www.mkm.ee/failid/The_Grid_Code.doc
Latvian grid code.
Legal acts about electricity market in Lithuanian.
http://www.ukmin.lt/lt/veiklos_kryptys/energetika
/elektra/
Official internet site of Estonian Power system.
http://www.pohivork.ee/index.php?id=311&L=1
Balancing principles for Latvia and Estonia, Arnis
Staltmanis Nowember 2006. Vilnius, Mini forum.
http://www.energyegulators.eu/portal/page/portal/EER_HOME/EER
_INITIATIVES/ERI/Baltic/Meetings1/SG_meetin
gs/3rd_Baltic%20SG/DD/A._Staltmanis_Balanci
ng.ppt#330,1,Slide 1
3.2. Latvia
Latvia already has in place legislation for
implementation of GO for RES-E and CHP-E:
GO for RES-E adopted on 24 July 2007 by
Regulation of Cabinet of Ministers No. 503
Regulations on Electricity Generation from
Renewable Energy Sources,
GO for CHP-E adopted on 6 November 2006 by
Regulation of Cabinet of Ministers No. 503
Regulations Regarding Electricity Production in
Cogeneration.
According to these regulations Issuing Body of GO is
Ministry of Economy. Rules of GO issuing and the
status with regard to the registry are still unclear.
Therefore GO was not issued yet.
Legislation regarding disclosure of generation attributes
is not adopted in Latvia.
Latvia is net electricity importer. It produces only
around 2/3 of electricity consumed. Major part of
electricity import comes from Lithuania and Estonia,
some from Russia and Scandinavia (Finland). Latvian
TSO, Latvenergo, provides information about local
electricity generation and amounts of imported
electricity. Information about CO2 or nuclear waste is
not available. Currently Latvia doesnt have plans on
introduction of disclosure.
3.3. Estonia
Estonia already has in place legislation
implementation of GO for RES-E and CHP-E:
3.1. Lithuania
Lithuania already has in place legislation
implementation of GO for RES-E and CHP-E:
for
184
for
6. Conclusions
1. Implementation of the proposed E-TRACK standard
could be beneficial for all three Baltic States
accommodating integration into EU electricity and
its attributes markets.
2. Taking into account current level of available
information E-TRACK standard could be
implemented at the relatively low cost.
3. Based on current electricity market structure of
Baltic States it is preferable to have one domain for
all three countries.
4. Despite current cautions of all stakeholders groups
implementation of E-TRACK standard could have
positive impact on RES-E electricity and overall
market liberalization.
7. References
1. Konstantinaviciute I., Tarvydas D. Tracking
Systems for Electricity. Energy 2007, Nr. 3. P. 33
39.
2. A European Standard for the Tracking of Electricity.
Final report from the E-TRACK Project. - August
2007. 114 p.
3. The E-TRACK Standard August 2007, 22 p.
www.e-track-project.org
5. Position of stakeholders
From provided discussions with major stakeholders and
from questionnaires it is evident that in Baltic States
185
186
Nuclide
Thermal neutron
(0,025 eV) , b
U235
584,4
U238
1,17710-7
Pu238
17.89
Pu239
747,4
Pu240
0,0588
Pu241
1012
Pu242
2,55710-3
Am241
600,4
Am242
2100
1 barn ( 1b 10 24 cm 2 )
210
205
200
66
195
190
Am242
Am241
Pu242
Pu241
Pu240
Pu239
U238
U235
185
3.2
3.0
2.8
2.6
2.4
2.2
Pu242(F)
Pu241(F)
Pu241(T)
Pu240(F)
Pu239(F)
Pu239(T)
U238(F)
U235(F)
2.0
U235(T)
Number of neutrons
Fission neutron
, b
1,235
0,3284
1.994
1,800
1,357
1,648
1,127
0,2296
1,756
187
Reactor
type
PWR
CANDU
BWR
HTGR
FBR
Fuel
Moderator
Coolant
UO2
enriched U235
Natural UO2
(0.7 % U235
UO2
enriched U235
UO2
( 8-19 %)
UO2/PuO2
( 16-20 %)
Light
water
Heavy
water
Light
water
graphite
Light
water
Heavy
water
Light
water
Gas
none
Liquid
metal
Number of Units
300
250
200
150
100
50
0
PWR
BWR
CANDU
HTGR
FBR
60
0.6
Composition, %
Number of Units
50
40
30
20
10
0.5
0.4
0.3
0.2
0.1
0
Pu238
Pu239
Pu240
Pu241
Pu242
0
PWR
BWR
PHWR
FBR
HTR
235
Composition, %
1.E+00
1.E- 01
1.E-02
1.E-03
1.E-04
1.E-05
1.E-06
1.E-07
1.E-08
Pu236
Pu238
Pu239
Pu240
Pu241
Pu242
Composition, %
100
10
0.1
U235
Pu238
Pu239
Pu240
189
Pu241
Pu242
Am241
18000
Contribution, %
16000
14000
12000
10000
8000
PWR
6000
RBMK
4000
2000
0
0
10
20
30
40
50
60
Burnup, MWd/kg
Coolant
Fuel
GCFR
LCR
MSR
fast
fast
Epithermal
fast
helium
Pb-Bi
Fluoride
salts
sodium
Thermal
or fast
thermal
water
U238
U238
UF in
salt
U238,
MOX
UO2
helium
UO2
SCR
SWCR
VHTR
Power,
MWe
288
50-1200
1000
1501500
1500
20
18
16
14
12
10
8
6
4
2
0
EV
To A
s
G hib
E/
a
H
ita
ch
R
i
us
si
Ko an MH
I
In
re
an du
s
tr
I
C
hi ndu y
ne
st
y1
se
0
in
In
du
di
s
an
t
i n ry
du
st
ry
AE
C
L
Neutron
AR
Type
Number of contracts
250
190
1
2
3
4
5
6
7
ABWR
APR-1400
AP-600
AP-1000
ESBWR
APWR
US-APWR/EUAPWR
ATMEA
EPR
SWR 1000
Advanced BWR
AES-92
(VVER-1000)
AES-2006
(VVER-1200)
VVER-1500
CANDU-9
Enhanced
CANDU-6
ACR-700
ACR-1000
PBMR
GT-MHR
8
9
10
11
12
13
14
15
16
17
18
19
20
Power,
MWel
1400
600
1000
1520
1538
1700
Available
for
construction
Yes
From 2010
Yes
Yes
From 2010-2011
Yes
From 2012
1100
1600-1750
1200-1290
1500
1000
Not yet
Yes
Yes
Not yet
Yes
1150-1200
Yes
1500
925-1300
750
Not yet
Yes
Yes
700
1080-1200
165
285
Not yet
Not yet
After 2013
Not yet
48
Burnup, MWd/kg
No
46
44
42
40
38
36
3.2
3.4
3.6
3.8
Enrichm ent, %
11. References
1. W.M. Stacey. Nuclear reactor physics. John Wiley
and Sons, Inc. 2007. 735 p.
2. A.R. Foster, R.L. Wright, Basic nuclear engineering.
Boston. 1993. 611 p.
3. D. Bodansky. Nuclear energy. Principles, practices
and prospects. Springer, 2004. 701 p.
4. H-D Berger, High Fuel Burn-ups and related safety
consideration, Presentation at 2006 FJOHSS, 23
August 1 September, 2006.
5. WNA information on Uranium Markets,
www.world-nuclear.org
Abstract: Nuclear energy is a key component of lowcarbon and secures energy supply policy in the context
of the increasing demand of energy, depletion of fossil
resources, increasing energy prices, growing energy
dependence and greenhouse emission.
This article presents the overview of the electricity
generation mix in the Baltic Sea region. It is based on
the analysis of the nuclear energy contribution to the
total electricity generation in different countries.
3.0
Generation, TWh/a
2.5
100
80
60
40
20
2001-2007
1996-2000
1991-1995
1986-1990
1981-1985
2006
2002
1998
1994
1990
1986
1982
1978
120
1976-1980
1974
1970
140
1971-1975
1.0
0.0
1965-1970
1.5
0.5
1. Introduction
Number of Reactors
2.0
Despite the fact that after this period there was a sharp
decrease in the number of nuclear power reactors put
into service, nuclears share of global electricity
generation all time was increasing (Fig.2).
192
60
Estonia
Lithuania
50
Latvia
Total
40
30
20
10
0
1990
1995
2000
2005
8
70
60
50
Nuclear share, %
80
40
30
20
10
Czech Rep.
Bulgaria
Hungary
Switzerland
Slovenia
Sweden
Belgium
Slovakia
Lithuania
France
-2
-4
Estonia
Latvia
Lithuania
Total
-6
1992
1996
2000
2004
193
2008
120
Production
100
2.5
Electricity, MWh
Power, TWh
3.0
2.0
1.5
1.0
0.5
Export
80
60
40
0.0
LT-KGD
BY-LT
LT-LV
LV-RU
EE-LV
RU-EE
EE-RU
20
0
0
10
20
30
40
50
Weeks
60
Hydro
Nuclear
Coal
Gas
50
Share, %
40
30
20
10
0
1990
1994
1998
2002
2006
6
Exchange, TWh
Import
Export
4
6
3
2
1
0
Denmark
Finland
Norway
Poland
Germany
194
Electricity generation in Sweden is dominated by CO2free hydro and nuclear power (Fig.10).
70
Nuclear
Hydro
Coal
Gas
60
Schare, %
50
40
30
20
10
0
1990
1994
1998
2002
2006
NC,
MWe
1167
1240
1370
771
1329
1275
1360
1284
1288
878
1400
1346
785
1305
890
1392
1345
REC,
TWh
13.7
23.0
128.4
11.0
141.8
69.6
113.8
79.7
89.1
24.2
139.1
88.6
10.2
152.2
28.7
113.5
43.6
60
50
Share, %
1975
1977
1986
1977
1988
1982
1985
1984
1985
1979
1988
1984
1976
1989
1980
1985
1979
GC,
MWe
1225
1300
1440
806
1400
1345
1430
1344
1344
912
1475
1402
840
1400
926
1458
1410
Nuclear
Coal
Year
40
30
7. References
20
10
0
1990
1994
1998
2002
2006
195
Production, GWh
40
35
30
25
20
15
10
5
0
Heat
Hot water
10 11 12
Months
196
(gCeq/kWh)
Total
228-366
206-357
149-246
106-188
6.3-64.6
8.4-16.6
8.2-76.4
2.5-13.1
2.5-5.7
Share, %
30
LT
25
EU-27
20
15
10
5
0
EI
197
30
EE
LV
LT
20
45
15
40
10
Production, GWh
CO2 Emission, Mt
25
5
0
1990
1994
1998
2002
2006
EE
35
LV
LT
30
25
20
15
10
5
0
1996
1998
2000
2002
2004
2006
45
CO2 emission, Mt
40
35
30
25
6
20
15
10
5
0
DK
EE
LV
LI
FI
SE
198
199
200
WE
RE
W W
E
Q
WQ
RQ
W W
E
Q
(1)
I ,
WE
RE
W W
Q
E
WQ
RQ
W W
Q
E
S i S i I , R R1 , R2 ,...Rn (3)
(2)
i, s i S i Ri ( s i , s i ) Ri ( s i/ , s i )
(4)
s i S i Ri ( s i , s / i ) Ri ( s i/ , s / i )
(5)
201
c i R( S
i ) R( S )
(6)
iS N
S ! (n 1 S )!
n!
R( S
i ) R( S ) (7)
1 2 R( S
i ) R( S ) / 2
(8)
DISTRICT Z
Electricity
Electricity
DISTRIBUTIVE
GRID
Name
1
2
Investments
Summary costs of fuel (per
hour)
Price of electricity exported
Price of electricity imported
Amount of heat supplied to
the grid
Amount of electricity
supplied to the grid
3
4
5
No.
Amount,
Measuring units
6834 thous. EUR
980622 EUR
105 EUR
126 EUR
183,6 thous. MWh
16,6 thous. MWh
Heat
CHP
(9)
Electricity
...
...
Heat
CHP
203
3.
126,00
100
100
EUR/ MWh
120
EUR/ MWh
120
80
60
40
38,50
Price of heat
80
Price of electricity
60
40
20
115,00
37,60
5. Acknowledgements
20
0
Variant B1
Variant B2
6. References
Fig. 4 schematically shows the comparison of heat and
electricity tariffs in Variants B1 and B2. The mentioned
prices, for which the consumers of district Z will buy
the thermal and electric energy, differ markedly from
one another. Hence, the application of the method based
on the game approach results in electricity prices that
are sufficiently close to the market prices as well as to
relatively low heat prices. It appears that such price
ratios stimulate the extensive application of
cogeneration processes in energy production.
Fig. 5 represents the prices of heat and electricity,
excluding the State support. In this case, the heat prices
appear to be higher than the average heat prices in Riga.
Hence, the considered Variant turns out to be noncompetitive.
126,00
100
100
80
60
55,40
EUR/ MWh
120
EUR/ MWh
120
115,00
Price of heat
80
60
40
40
20
20
54,50
Price of electricity
0
Variant B1`
Variant B2`
4. Conclusions
1.
2.
204
205
R MC T ATC ;
where: R congestion revenues; MC difference of
marginal electricity prices between interconnected
parties; T duration of congestions; ATC available
transfer capacity of interconnection.
Economic principles of cross-border trade between
regions A and B are shown in Fig. 1. Without
interconnection between those regions in region A
electricity price would be MC A1 and in B region MC
B1, and MC A1 > MC B1.
If available transfer capacity exists between those
regions, marginal price of A region would decrease till
MC A2 and in B region would increase till MC B2. In
case that there are no congestions between regions,
electricity price would be equal MC AB in both regions.
Price
Net import curve region A (deficit)
MC A1
MC A2
Welfare
Welfaregain
gain
MC AB
Congestion rent
Welfare loss
MC B2
Welfare gain
Net export curve region B (surplus)
MC B1
ATC
Volume
206
Price,
EUR/MWh
Market area B
Import
Quantity, MWh
X MW
Export
Quantity, MWh
Export from B/
Import to A
Market area A
Market area B
Market price P2
0 MWh
0 MWh
Available transfer
capacity[ATC]
Price,
EUR/MWh
Import
Quantity, MWh
Market price P2
Effect of
interconnection
Import to B/
Export from A
Market price P1
Market price P1
Market area A
Export
Quantity, MWh
207
Price,
EUR/MWh
GENERATION
CONSUMPTION
A area market
price P1
P2
CONSUMPTION
GENERATION
Benefit for
consumers
P1
Benefit for
generators
B area market
price P2
Decrease of generation
price = decrease of
customers costs
Increase of generation
price = increase of
customers costs
SELL BIDS
(GENERATION)
Export f rom B
Import to A
Quantity, MWh
Import to A
Export f rom B
Quantity, MWh
4.3. Third
capacity
scenario
limited
interconnection
With the condition that available transfer capacity ATC
is equal to actually needed capacity X. Available
transfer capacity ATC is lower than the quantity of
proposals to sell electricity using interconnection
capacity, what is corresponding to demand curve in the
higher price area.
Having price difference between external and internal
market, congestion rent is being cumulated by electricity
exchange, which could be distributed among
transmission system operators at the end of the year. The
level of congestion rent is indicating how often
congestions occur and is it necessary to increase
interconnection capacity. On the other side, these
revenues could guarantee the return of investments of
new interconnection, therefore further analysis of
interconnection return of investments is related with
possibility to forecast electricity price in two analyzed
markets and to distribute congestion rent for the
investment.
Market area A
TSO revenues
Market area B
Import
Quantity, MWh
Market price P2
Effect of
interconnection
0 MWh
X MW
Export
Quantity, MWh
Eksportas i B
Importas A
Import to B/
Export from A
Available transfer
capacity[ATC]
208
Table 1. Forecast of Nord Pool electricity physical exchange price in 20102019, EUR/MWh
Year
Maximal Nordic electricity
market price
Minimal Nordic electricity
market price
2010
2011
2012
2013
2014
2015
2016
2017
2018
2019
53.0
55.1
54.0
56.2
55.0
57.2
56.0
58.2
57.0
59.3
37.9
38.1
41.0
45.0
46.0
47.8
49.8
51.7
53.8
56.0
2010
2011
2012
2013
2014
2015
2016
2017
2018
2019
98.0
101.9
106.0
110.2
114.6
119.2
124.0
129.0
134.1
139.5
72.1
75.0
78.0
81.1
84.3
87.7
91.2
94.9
98.7
102.6
209
6. Conclusions
1. Success, economic and social benefit of cross-border
interconnection projects depend on the electricity
market conditions and system control possibilities,
because these subjects enable objectively evaluate
pay-back time of the investments, its security, social
and economic benefit in the long run.
2. Marginal electricity price in the isolated electricity
market according to the first analyzed scenario,
without physical connection would be developing
independently from the price of neighboring
countries. Social and economic benefit of common
electricity market is not achieved and benefit is
obtained only through the deployment of internal
market resources and internal competition.
3. Simulation of interconnection according to the
second scenario with sufficient interconnection
capacity shows that generators in the exporting
country could increase their revenues, but costs for
consumers would be higher. And on the contrary,
consumers costs would be lower in the importing
country, because of the lower electricity price. There
is no congestion rent for interconnection operator,
because there is no price difference between the
markets due to sufficient interconnection capacity.
4. Simulation of interconnection according to the
second scenario with limited interconnection
capacity shows that different price areas would
develop. Having price difference between external
and Baltic market, congestion rent is being
cumulated by electricity exchange, which could be
distributed among transmission system operators at
7. References
1. Analysis of Cross-Border Congestion Management
Methods for the EU Internal Electricity Market.
Study commissioned by the European Commission.
CONSENTEC, Final report, June 2004.
2. Market Operator (Regulatory Oversight), Final
Version. Energy Regulators Regional Association
(ERRA), Prepared by KEMA International B.V.,
December 2007.
3. Towards a common coordinated regional congestion
management method in Europe. Study commissioned
by European Commission, CONSENTEC, October
2007.
4. Market information. www.nordpoolspot.com, 2008.
210
211
xR xm cos ;
x R ym cos y R xm sin
x R2 y R2
where:
xm, ym speed of rotor mass centre in coordinates x, y,
x m xm ;
dS
dt
dR
yR
cR
dt
dR
dt
u p0R S rS iS ;
(4)
p0 ( R ) R rRiR
where:
ik, iik (i=S, R; k= , ) full linkage magnetic fluxes
and currents,
u, u stator voltages,
rS, rR
resistance,
p0
number of pole pairs (indexes S and R
indicate belonging to stator and rotor, and and belong
to separate electric contours).
Stator voltages we find as:
u U m sin(0t p0 R ) (5)
0 xR xm
(3)
p0 (R ) R rRiR ;
u U m cos(0t p0 R );
ym
u p0R S rS iS ;
dt
dS
y m ym ;
cm
(1)
y R ym sin
imk m mk k , m S , R; k ,
R2
R1
(6)
S aS bS R aR bR ;
S mS nS R mR nR
where:
212
(7)
1 TN M TS 2
cos p0
M b
d ;
( )
0
2
N b
0
sin 2 p0
d ;
( )
1 TM N TS 2
sin p0 cos p0
S b
d ;
( )
0
2
T S R
1 TM 1 TM TS 2 ;
p0
3w2
d xm 1
Fx c x xR x xm sin ;
dt
m
d ym 1
Fy gm c y y R y ym cos ;
dt
m
dym
d
d 1
ME MM m xm sin m
mg cos ,
dt
dt
dt
(13)
(8)
Electromagnetic moment
ME
3
p0 iS iS
2
(9)
where
2
2
;
F
c
sind ,
cos
0
2
)2
0
0 ( )
2
F c0
There
3w
iS iR
w
p0
(11)
cos p0
3w
iS iR
w
p0
sin p0 ,
(12)
213
Efficiency [%]
100
S um -Squared E rror
10
50
1,00E-01
100
200
300
400
E poc h
500
600
700
1,00E-04
1,00E-05
1,00E-06
1,00E-07
5. Conclusions
800
8
Learning Rate
1,00E-03
-5
-3
1,00E-02
x 10
30
60
25
70
10
20
80
15
90
10
10
6
4
6. References
2
0
0
100
200
300
400
500
E poc h
600
700
800
900
214
215
i1d = i1q = 0.
d 1d
U 1d = R1i1d 0 el 1q +
d
d 1q
U 1q = R1i1q + 0 el 1q +
d
U 2d
U 2q
,
d 2 d
= R2 i2 d ( 0 el ) 2 q +
d
d 2 q
= R2 i2 q + ( 0 el ) 2 q +
d
TM
d
= [ M em M l ] ,
d
1d = X 1i1d + X ad i2 d
1q = X 1i1q + X ad i2 q
2d
2q
,
= X 2i2 d + X ad i1d
= X 2i2 q + X ad i1q
(4)
1
X2
X 1 ad
X2
X 1i2 d
X 1 0el i2 q ,
2
X
X 1 ad
X2
(1)
di2 q
1
=
X 1 0el i2 d
2
dt
X 1 X ad
X2
(2)
R2
X2
1
X1
2
X ad
X2
(6)
X 1i2 q ,
1
d
M l.
=
dt
TM
(3)
(5)
(7)
216
R2 X ad
i2 d X ad i2 q ,
X2
R X
= X ad i2 d 2 ad i2 q ,
X2
U 1d res =
U 1q res
(8)
(9)
217
5. Conclusions
6. References
218
Fig. 1. Two-core axial inductor machine: 1 body, 2 armature cores with windings, 3 rotor cores, 4 rotor shell,
5 excitation winding, 6 shaft, ( is a magnetic flux)
x y x y x y
(1)
( AAB ACD ) l A l ,
2
(2)
220
where
In order to obtain optimal geometry of teeth area of twocore inductor generator having axial excitation and
toothed windings it has to be used objective function,
representing correlation between generators nominal
power and its mass (3):
(3)
Ik
2
E0
xS 2
(5)
mE02
,
2 xS
(6)
max min
,
2
a
a
120G
k w2 nN Z1Z 2
max min 2
max(k ) min( k )
(11)
(7)
IN
(10)
Fk
(4)
2Z1
,
ma
(9)
but the number of coils in the parallel branch of twocore generator (10) is
nk
mU N I N
,
C
G
UN
Z 2 nN
,
60
(8)
221
5. Conclusion
1. .., .., ..
. Latvian
Journal of Physics and Technical Science, Riga,
1991, Nr.1, pp.87-94.
2. .. .
- , 1984, .247.
3. Orlova S., Pugachov V., Levin N. Magnetic field in
the tooth zone of an axial inductor machine (AIM).
5th annual Conference of Young Scientists on Energy
Issues 29 May, 2008, Kaunas, Lithuania, CD
proceedings. II-20 II-27
4. QuickField User Guide, version 5.5. WEB
http://www.quickfield.com
222
1. Introduction
In new technologies there are widely used both rotary
motors of common structure and linear electric motors.
In cases when special technologies require linear or
sliding motion there are applied LIM with the secondary
element shorter than an inductor. Such an executing
device with the help of a sliding magnetic field not only
expels the secondary element conductive to electricity
from the active zone of the inductor but has the capacity
to perform the functions of an electromagnetic brake, a
damper, a spring, an absorber and executes the other
functions of an energy converter [1, 2]. The motors
mentioned above could be successfully installed in
compressors, in repulsion and shok hammering
mechanisms, in aggressive and hazardous technologies
such as for extraction of conductive metallic parts when
utilizing and decontaminating mercury, luminescence,
daylight lamps and other gas discharge tubes.
223
KF =
1+
0, 63
2c ( 2b 2c )
(1)
KF =
1 Re th(c1 1 is 0 ) sv 0
1 sv 0 + s 2 02 Re th(c1 1 is 0 )
(2)
H = H1 + H 2 ;
is a
(b c)
.
(3)
Felm =
(b c)
1.3 .
0.41 .
0
2
Re H 2 &j ( x, t )dV ;
(7)
(V )
(4)
(6)
Felm = F1 + F2 + F3 ;
(5)
(8)
224
F3 = F0 Re
1
p 3
(2 K + 1)
0 v0 n sin 2
K =0
( )
( 0v )
+ nK2
2
1 exp p 0 v + ( 0 v ) + nK
+
2
K 3 ( 0 v0 n K 3 )(1 i K 3 )
1 exp p 0 v + ( 0 v ) + nK2
+
2
K 4 ( 0 v0 n K 4 )(1 i k 4 )
(9)
1.5
F3 / F0
1.1
0.9
active zone; = b
= 0v
K3
v+
K 4 = 0
1.3
( 0v )
( 0v )
0.7
0 = 50
0.5
0.3
0.1 0.2
0
0.4 0.6 0.8
-0.1
+ nk2
0 = 10
-0.3
(10)
2.2
0 = 1
-0.5
+ nK2
c /
nK =
2K + 1
;
c
K = 0; 1; 2; 3; and etc.
1.5
1.3
F3 / F0
0.9
0 = 50
0.7
0.5
0.3
sh ( c 1 i 0 v )
hs = Re 1
c
1 i 0 v
ch c ( 1) 1 i 0 v
ch ( c 1 i 0 v )
0 = 10
1.1
0.1 1
0
-0.1
-0.3
-0.5
(11)
0 = 1
p
225
5. Conclusions
1.75
F3 / F0 1.5
1.25
1
0.75
0.5
0.25
0
1 1.2
1.4
1.6
1.8
2.2
2.4
2.6
2.8
-0.25
-0.5
4
2
c/
0
-2
-4
0.2
0.4
0.6
0.8
1.2
1.4
1.6
1.8
6. References
-6
-8
-10
226
d 1
M M st1 M12 = J1 dt ;
M + M = J d 2 .
12
2
st 2
dt
(1)
(2)
1. Introduction
Electromechanical system as object of investigation
comprises
electrical
and
mechanical
parts.
Electromechanical power converter and its control
system depends to electrical part as well as all moving
masses coupled between themselves form mechanical
part. Electromechanical system includes various
mechanical chains, with infinite or finite elasticity and
clearance. Systems with infinite stiffness and without
clearance are one-mass systems and are quite well
analyzed [1, 2, 3]. Systems with capable to deform
chains and clearance are more complex. They are
described by high order nonlinear differential equations,
and without essential simplifying of problem they
cannot be solved in analytical way. In these cases
computer models of solved program must be developed,
using specialized software, and system responses
simulated.
Some problems of two-mass system were considered in
[4, 5, 6, 7, 8, 9, 10].
The paper presents general model of two-mass system,
developed in Simulink. Results simulations at different
stiffness of two-mass system and different clearance are
presented and considered. Dependences of maximal
speed values against elasticity factor and clearance are
formed and presented.
c12 =
(3)
(4)
M st 2 + M12 = J 2 s2 ;
sM12 = c12 (1 2 ) .
(5)
227
Parameter
Motor power, kW
Number of pole pairs
Phase voltage, V
Power factor
Rated torque, Nm
Rated current, A
Inertia, kgm2
M 12
b)
Value
1,1
2
230
0,81
7
3,56
0,00262
J=
M M st1 M 12 = J1s1 ;
M st 2 + M12 = J 2 s2 ;
M 12 = c12 1 2
, when
2
M 12 = 0, when 1 2
.
2
(7)
(6)
;
1 2 >
2
m r2
,
2
250
200
300
150
100
60
400
200
30
motor 600
150
10
100
50
0.01
0.02
0.03
0.04
0.05
Time, s
0.06
0.07
228
0.08
0.09
0.1
Simulation results are compared with results of onemass system operating at no-load for exploring
dependences of system parameters influence. Clearance
is determined by time duration after which the first mass
(motor) starts to move the second mass.
Stiffness value was chosen freely for getting
dependences of speed maximum value against stiffness.
Comparison of simulation results indicate motor speed
lag in two-mass system, depending on stiffness, with
speed of motor, operating in one-mass system. That
depends on the load torque, developed by second mass.
After the shaft dislodges the second mass, the motor
speed reaches the greater value than that at no load due
to inertia of the second mass accelerating the motor.
Oscillations of speed in steady-state also are greater
because elasticity of the shaft. The shaft under rotation
of motor is being twisted and deformed. In this way it
accumulates energy and delays the influence of the
X: 0.02182
Y: 194.4
200
Rotation speed, rad/s
X: 0.02638
X: 0.02873
Y: 216.3
Y: 211.8
X: 0.03237
Y: 197.7
X: 0.02458
Y: 214.2
220
X: 0.02262
Y: 204.3
X: 0.01542
Y: 171.6
180
X: 0.01435
Y: 150.4
160
140
X: 0.01382
Y: 122.5
120
X: 0.03736
Y: 177.3
X: 0.02087
Y: 183.1
X: 0.01488
Y: 163.5
X: 0.03948
Y: 163
X: 0.04065
Y: 158.6
X: 0.01435
Y: 135.9
X: 0.01329
Y: 111
100
0.005
0.01
0.015
0.02
0.025
0.03
0.035
Time, s
0.04
0.045
0.05
0.055
dz=-5-5
dz=-2-2
dz=-0.5-0.5
dz=0
250
200
dz=-10-10
dz=-20-20
dz=-15-15
dz=-25-25
dz=-30-30
dz=-35-35
dz=-40-40
motor
150
100
50
dz=-45-45
dz=-50-50
0
0.05
0.1
0.15
0.2
Time, s
229
0.25
5.
250
1.
200
150
100
synchronous
1 peak
2 peak
50
0
0
10
30
300
250
200
150
100
synchronous
50
peak
0
0 0,5 2
10 15 20 25 30 35 40 45 50
Clearence
1.
2.
3.
4.
230
Ui vzWz d .
(1)
1. Introduction
The quantitative characteristics of fluid flow: local and
mean velocity, the mass and volume flow and quantity
are important for many technological processes and
commercial accounts. Its measurements compose about
40% of all technological measurements which are
performed in industry. The conventional problems are
the accounting of drinking and industrial water and heat
measurement. The actual problems with arising of
ecology requirements are the accounting of sewages, the
batching of different special fluids and other. The
significant part of meters using for this purposes are the
electromagnetic fluid flow meters EMFM. The EMFM
are convenient for use in different informational and
control systems and in the devices, intended for different
combined measurements, par example, in heat meters.
The modern EMFM are intellectual devices, in which the
microcomputers are used for measurement cycle
organization, for subsidiary parameters calculation, for
parasitic signals suppressing, for measurement error
diminishing.
The mean advantages of EMFM are absence of
resistance to fluid flow, practically non-inert
measurement, linear transfer characteristic and small
measurement uncertainty. There is what for commercial
U KQ, K K1 ,
(2)
231
Ie
Ue
VCC
CC1
Ur1
Ur2
CC2
R0
U0
Ue
EC1
R0
2
EC2
Ie
Ur
U0
rKie,
(3)
ie I 0 (1 e
);
(4)
232
T0
Tm
N0KADCUI0KADCKIAI0R0T0/i,
Tp
(5)
Ur
U Ir
ADC
Tm
K IS K I 0
K IS K I 0
dt
K IS K I 0
Tm
NrKADCUI0KADCKIAKI0/i.
Ue
(10)
K R0T0
There are in this circuit: IA input amplifier, I integrator, ADC converter analog-to-digital, P processor, C commutator, M memory unity.
At first the voltage U0R0I0 is connected via commutator
C to amplifier IA input. This voltage acts in resistor R0
(see Fig. 2 or Fig. 3), when the transient process is finite
and excitation current ie is equal to I0 (interval T0 in the
time diagram, presented in Fig. 4).
Therefore the voltage U0 is connected to IA input during
the interval T0. This voltage amplifies into amplifier IA
and integrates during all interval T0. After the end of
interval T0 the integrator output voltage is equal:
UI0KIAI0R0T0/i,
(8)
(9)
C
IA
K I 0 t
di
d r
e.
K e
e
dt
dt
e
Ur
U0
(7)
Nr
.
N0
(11)
(6)
233
2.
3.
a
When the two control coils are mounted we can not only
observe the variation of transfer coefficient, but we can
understand the reason of this variation. It is usefully to
mount the coils CC1 and CC2, as it is showed in Fig. 2:
one under the central part of magnetic core, other
under the edge of excitation coil. When the signals of
both coils decrease, therefore, we can suppose, that
permeability of magnetic circuit decreases. When the
signal at least of one coil increases and the signal of
central coil increases more than the signal of peripheral
coil, therefore, the magnetic admixtures appear in the
flow (see Fig. 6, b). When, contrary, the signal of
peripheral coil increases more than the signal of central
1.
2.
3.
234
1. Introduction
235
H max
2
Z
H R2
2
Z
2
H R2 4 H Z2 H R2 sin 2 Z R , (1)
In
r2
2
x2
1,5
(2)
B, mT
b)
2
rx cos
d
1.5
2
2
In 0 r x 2r cos
H
,
4 2
r2
d
1.5
2
2
0 r x 2rx cos
75
70
61,97 mT
65
60
55
50
45
40
35
30
25
20
15
10
1,214 mT
0,505 mT
5
00
1600 A
1,25 mT
0,5 mT
31,354 mT
1,214 mT
-X -10 -9 -8 -7 -6 -5 -4 -3 -2 -1 0
0,505 mT
9 10 X
Distance, m
(3)
100 ,
X r 00
(4)
a)
b)
a)
b)
237
H ZPK
I R0
N 1 2
I R0EK
K 12
j 0 0
HRPK
R0EK R cos
2
2
2
z g K jk R R0EK RR0EK cos
j0 0
(5)
d,
3
4 i0 0
2
2
2
z in R R0 2RR0 cos
. K12
z hp cos
I R0 N12
0EK
(6)
z gK hk cos
3
2
2
2
z gK jk R R0EK 2RR0EK cos
(7)
6. Conclusions
I R
R R cos
0
d
3
4 i0 0
2
2
2
z in R R0 2RR0 cos
d.
7. References
.//
3-
- 2008, 2008. .
81-90.
5. J. Schaffer. Commutating Current-Limiters an effective
alternative for high current protection.// NETA World, vol.
18, Winter 1996-97. - No. 4. - P. 7-18.
6. .., .., ..,
..
. //
IX-
" 2030.
-.", 2007.
7. ABB AG. Calor Emag Medium Voltage Products. Islimiter.// Leaflet no. DEABB 2243 08 E Leaflet no.
DEABB 2243 08 E. - Germany, 2008. - P. 24.
238
239
4. Switching regulator
240
5. A charging stage
At this stage, the inductor [L], main assignment is to
limit the current slew rate through the power switch. It
creates a fast energy storage which is expressed as
follows:
E 0.5 LI L2 .
di VL
dt L
(1)
(2)
VL L
di
.
dt
VIN
tOFF
VL
tON
(3)
t
VOUT VIN 1 ON
tOFF
6. Discharging
(4)
241
t ON
,
t OFF t ON
(5)
VOUT
VIN
.
1 D
(6)
9. References
1. Marcel P.J., Gaudreau P.E., Casey J.A. , Hawkey
T.J., Kempkes M.A., Mulvaney J.M., Ver Planck P.
Solid State High Voltage DC Power distribution and
control. In proceeding of the Particle Accelerator
Conference, New York, 1999.
2. Adoko S. Effective Paradigm for direct current (DC)
distribution using light emitting diodes (LEDs) in
architectural lighting system. The 3rd Regional
World Renewable Energy Congress (AISEC-10),
Kuwait City, Kuwait, 2008.
3. Adoko S., Irowa O. Integrated power supplys
design and contruction practice for solar lighting
applications. Energy and Environment (EE09)
WSEAS Coference, Cambridge, UK, 2009.
4. Theraja B.L., Theraja A.K Textbook of Electrical
Technology. Ninja Construction & development Co.
(p) Ltd Ram Nagan, New Delhi, 1997.
5. Bose B.K. Adjustable speed AC Drives A
technology status review. IEEE Trans Ind. App.,
1982, vol. 70, no. 2. p. 116-135.
6. Herrit B., Morogavel R. PWM DC Motor Control
Using Timer A of the MSP 430. Application
Report SLAA 120 Texas Instruments December
2000.
7. Schaefer C. Field Weakening of Brushless
Permanent Magnet Motors with Rectangular Current.
Proc. Env. Conf. Power Application (EPE),
Florence, 1991, vol. 3, p. 429-434.
8. Karunadasa J.P., Renfrew A.C. A flexible fast Digital
Controller for a Brushless DC motor. 4th Int. Conf.
Power Electronic and Variable speed Drives, 1990.
p. 429-434.
9. Matsoi N. Sensorless PM Brushless DC Motor
Drives. IEEE Trans Ind. Electronics, April 1996,
vol. 43, no. 2. p 300-308.
8. Conclusion
DC motors speed control technique has been hitherto
very difficult and energy wasteful method to implement.
The modern switching technique as adopted in this work
242
243
T,K
306
2 min
302
3 min
300
298
- corona
T,K
310
1 min
308
2 min
306
3 min
4 min
304
5 min
302
6 min
8 min
304
300
298
operating
cycle
5 min
306
302
4 min
0
1
2
3
4
5
6 n
Fig. 2. Variation of temperature at 6 equally distant points in
overall length of the condenser tube influenced by negative
corona field for pause of the compressor operation
control test
308
1 min
304
T,K
310
308
pause
t,min
300
0
1
2
3
4
5
6
Fig. 3. Variation of temperature at 6 equally distant points in
overall length of the condenser tube during the control test for
pause of the compressor operation
0
2
4
6
8
10
12
Fig. 1. Variation of temperature at the beginning of condenser
tube
T,K
308
3,5 min
2,7 min
1,8 min
306
304
302
0,9 min
300
298
296
0
n
0
1
2
3
4
5
6
Fig. 4. Variation of temperature at 6 equally distant points in
overall length of the condenser tube influenced by negative
corona field for operating cycle of the compressor
244
T,K
4 min
310
3 min
308
T,K
289
287
2 min
306
304
302
300
285
50Hz corona
283
+ corona
281
1 min
control test
279
- corona
277
275
t,min
273
0
1
2
3
4
5
6
Fig. 5. Variation of temperature at 6 equally distant points in
overall length of the condenser tube during the control test for
operating cycle of the compressor
T,K
275
271
273
10
20
30
40
50
60
t,h
50Hz corona
271
control test
269
- corona
267
+ corona
265
263
Fig. 7. Variation of frozen water temperature
6. References
T,K
283
+ corona
50Hz corona
278
control test
273
- corona
40
80
120
160
200
t,min
240
268
263
Fig. 8. Variation of melting ice temperature
246
2. Experimental arrangement
We use the direct corona electrode system a set of
parallel corona wires above the plane electrode (Fig.1)
Plane electrode is represented by the surface of drying
clay panel. Test drying clay sample is a plate of 1 cm
thickness and of 10x10 cm2 area. Radius of corona
wires is ro= 0.05 mm, the step of the grid d = 15 mm,
and the height h is variable in the interval 10-14 mm.
1. Introduction
r0
x
Fig. 1. Electrode system of the experimental device
We use the calculated data for volt-ampere characteristics of corona in the discussed electrode system at
several values of the height h to determine the proper
values of test voltage and current corresponding the
source voltage 10 kV and the current 1200 A. The
used formulas are the same as in [9].
Initial field strength of the discharge depends upon the
radius of the wire r0 and the relative air density :
247
E0 = 3030000 (1 + 0,0298
(ro )),
(1)
= pT0 / p0T ,
(2)
U 0 = r0 E0 2 ln 2 0 .
d
d
(3)
1 2 ln 2 0 .
=
d
U 0 d
(4)
r
d
(5)
10
6.36
11
6.69
12
7.02
13
7.35
14
7.68
Fig. 3. Computed volt-ampere characteristics for positive
corona
248
Fig. 4. Comparison of computed and experimental voltampere characteristics for negative corona
Fig. 6. Drying process in negative corona field
Fig. 5. Comparison of computed and experimental voltampere characteristics for positive corona
3. Results of an experiment
Drying test is performed at 19 C ambient air temperature, 45 % relative humidity and 101-102 kPa of an
atmospheric pressure. Drying clay sample is arranged
horizontally on the grounded aluminium sheet. Initial
weight of wet clay sample is 0.2410.001 kg.
Periodicity of sample weighing is 1 h. Each drying
session is performed twice: under free convection
conditions and with a sample subjected by direct current
corona field. Total discharge current is adjusted to 200
A at the voltage 7.5 kV (approximately, the value of
the voltage is slightly increased to compensate the
increasing of resistance of drying clay sample). Final
weight of dry clay sample is 0.2060.001 kg. Total
weight of water in wet clay sample is 0.035 kg. Curves
of drying process with negative and positive corona
field in comparison with natural drying conditions in
still air are presented in Fig. 6 and Fig. 7. Letter m
denotes the mass of drying clay sample, and t is the time
measured from the beginning of drying. Drying
experiment is finished at 10 h from the beginning when
the sample drying in corona field has no moisture.
mnd, %
0.6
2.8
8.6
14.3
22.8
31.4
m+cd, %
31.0
54.0
88.6
94.3
97.1
100.0
m-cd, %
28.6
51.4
82.8
91.4
97.1
100.0
r+, times
51.7
19.3
10.3
6.6
4.3
3.2
r-, times
47.7
18.3
9.6
6.4
4.3
3.2
249
5. References
4. Conclusions
250
1. Introduction
Conducting Sleeve Induction Motor (CSIM) is a special
induction motor with a solid rotor which is made of two
parts: ferromagnetic core and hollow conducting
cylinder. The material of hollow conducting cylinder
could be copper or aluminum. This motor has particular
applications such as helical motion induction machines
[1], [2], [3], [4]. In helical motion induction machines,
induced currents on the rotor must move in helical
direction, so squirrel cage structures can not be used.
Therefore using solid structure for the rotor is
unavoidable. One of important application of solid rotor
induction machine is high speed motors. Due to solid
(1)
(2)
(3)
(4)
251
Fig. 1. Three and two dimensional views of a three phase CSIM structure
permeability.
The following equation directly results from the above
equations:
H +
B
=0
t
(5)
2 (mm)
4 (mm)
15 (mm)
1 (mm)
0.5 (mm)
2000
250
5.8e+7 (S/m)
6
20
3 (A/mm2)
2
50
3
2
1.5
Torque (N.m)
2
1.5
0.5
0
0
0.02
0.04
0.06
0.08
0.1
0.12
0.1
0.12
-0.5
Torque (N.m)
-1
1
-1.5
Time (second)
0.5
0
0
0.01
0.02
0.03
0.04
0.05
0.06
-0.5
4500
Time (Second)
4000
3500
Speed (rpm)
3000
2500
2000
1500
1000
3.49E+03
500
2.99E+03
0
0
0.02
0.04
Speed (rpm)
2.49E+03
0.06
0.08
Time (Second)
1.99E+03
1.49E+03
9.90E+02
1.5
4.90E+02
1.00E-02
2.00E-02
3.00E-02
4.00E-02
5.00E-02
Currents (Ampere)
1
-1.00E+01
0.00E+00
6.00E-02
Time (Second)
0.5
Phase A
0
Phase B
0
0.02
0.04
0.06
0.08
0.1
0.12
Phase C
-0.5
-1
-1.5
Time (Second)
Currents (Ampere)
1
0.5
Phase A
0
-0.5
Phase B
0
0.01
0.02
0.03
0.04
0.05
0.06
Phase C
-1
-1.5
-2
Time (Second)
253
3500
3000
Speed (rpm)
2500
2000
1500
1000
500
0
0
0.01
0.02
0.03
0.04
0.05
0.06
Time (Second)
Currents (Ampere)
1
0.5
Phase A
Phase B
0
0
0.01
0.02
0.03
0.04
0.05
0.06
Phase B
-0.5
-1
-1.5
Time (Second)
2.2
1.7
Torque (N.m)
1.2
2.00E+00
1.50E+00
2D
0.7
3D
0.2
Torque (N.m)
1.00E+00
0.00E+00
-0.3
5.00E+02
1.00E+03
1.50E+03
2.00E+03
2.50E+03
3.00E+03
5.00E-01
-0.8
Speed (rom)
0.00E+00
0
0.01
0.02
0.03
0.04
0.05
0.06
-5.00E-01
-1.00E+00
Time (Second)
6. Conclusions
254
7. References
255
I.
II.
III.
Cost
Durability
Compliance with performance criteria as laid
down in specifications.
256
(2)
On substitution, Co 62.1 .
From equation (1),
D2 L = Pa /Co n .
Substituting the values,
D2 L = 0.0037037 m3
(3)
D2 L 0.0037 m3 .
(4)
or
2. Symbols
Pa =
D =
L =
n =
Bav=
ac =
P =
=
=
=
Co =
b =
bp =
=
Lf =
X =
p =
(1)
i.
ii.
iii.
(5)
257
D +L = X;
where X is a constant, which must be minimum.
D +L = X
(objective function)
D2 L = 0.0037 m3 (constraint equation).
(7)
(6)
L,
m
0.01
0.02
0.05
0.07
0.09
0.1
0.12
0.14
0.16
0.2
0.25
0.3
D2,
m2
0.37
0.185
0.074
0.053
0.041
0.037
0.0308
0.026
0.023
0.0185
0.0148
0.0123
D,
m
0.608
0.43
0.272
0.23
0.202
0.192
0.1756
0.162
0.152
0.136
0.121
0.111
258
Table 2. D + L =X (= 0.2)
D,m
0.2
L,m
0.2
Table 3. D + L = X (= 0.3)
D,m
0.3
L,m
0.3
D,m
0.4
L,m
0.4
Table 4. D + L = X (= 0.4)
Design
variables
D, m
L, m
(D+L), m
L f
L
L f
=0
=0
Methods
Graphical
0.192
0.100
0.292
Lagrangian
0.1948
0.0974
0.2922
L f
Conventional
0.195
0.100
0.295
8. Conclusions
A d.c.machine design problem is solved by
conventional method, which employs assumptions,
approximations, etc., Further it requires an equation in
terms of pole proportions for the separation of D and L,
which results in different values for D and L depending
upon the separating condition equation. None of these
values would be optimum. Hence the conventional
design problem is converted as an optimization
problem.
The procedure for estimating optimum
values for the main dimensions are explained with the
help of graphical method. The values are verified with
the values obtained by popularly available optimizing
method. The values are agreeable. The speciality of
this method is that optimum values for both the free
design variables are obtained simultaneously. Such
method can be used in intermediate steps in the process
of design of an electric machine, provided the constraint
equation and objective function contain only two same
variables. The example may be the estimation of
window dimensions of a transformer. The performance
of the d.c. machine would improve because the
optimum values are estimated at the first step itself
which will be used in further calculations. Hence such
graphical method may be adopted for obtaining
optimum values for the two variables at the time of
design of an electric machine at appropriate places.
(8)
(9)
(10)
= 0
7. Comparison of results
The values for armature diameter and core length
obtained by all methods are given in Table 5 for
comparison.
259
9. References
1.
2.
3.
4.
5.
6.
10. Biographies
C. Easwarlal was born in 1947.
He studied his B.E (Electrical and
Electronics Engineering) degree
course and
M.Sc (Engg)
(Electrical Machines) degree course
at PSG College of Technology,
Coimbatore - 641 004 and received
the degrees in the years 1971 and
1973 from University of Madras. He has worked as an
Electrical Engineer in a company for two years. Then he
took up the teaching assignment. He has put in 28
Years of experience. Currently he is working towards
his Ph.D degree in Anna University, Chennai, India. He
is a professor of Electrical and Electronics Engineering,
Sona College of Technology, Salem 636 005, India.
His area of interests is design of Electrical Machines,
optimization in design, power system analysis.
260
2. Soil resistivity
The property of resistivity can be defined for any
material. As applied to soil, resistivity is an indication
of a given soils ability to carry electric current. The
flow of electricity in the soil is largely electrolytic,
determined by the transport of ions dissolved in
moisture. The soil resistivity is the resistance measured
between two opposing surfaces of a 1m3 cube of
homogeneous soil material, usually measured in m.
The soil resistivity has a direct effect on the resistance
of the grounding system. The typical features of soils
with certain resistivity values are shown in Table 1.
These changes are caused by the influence of
admixtures and by different structure of the mineral
grains, on which the measurements were conducted.
Microscopic cracks and oxidations of surface within the
limits of individual grains produce significant changes
in the values of the measured resistances. The values of
the resistivity of the ground water are given in Table 2.
261
Resistivity
(m)
100
1001000
1000
Features
Lowland at the mouths of
rivers or by the sea. Usually
abundant in water.
Midland plains where
ground water is not so
difficult to obtain.
Hilly zones, piedmont
districts and high lands,
where drainage is good.
Resistivity (m)
200000
50000
200
70
20 ~ 70
2
0.3
0.2 ~ 0.25
0.15
20C
10C
0C
0C (ice)
-5C
-15C
2000 ~ 5000
104 ~ 107
1000 ~ 5000
Temperature
Resistivity (m)
10 ~ 150
10 ~ 200
50 ~ 100
100 ~ 1000
200 ~ 2000
Rate
1.0
1.4
1.8
4.2
11.0
45.9
262
Rl
0,366 lg
2l
d
(1)
Characteristic of
the climatic zone
I
Average longstanding lowest
-20 -15
temperature of air
(January), C
Average longstanding highest
16 18
temperature of air
(July), C
Average annual
amount of
40
precipitation, cm
Duration of the
freezing of water, 190 170
days
1 2 ... n
-14 -10
-10 0 0 +5
18 22
22 24 24 - 26
50
50
30 - 50
150
100
(2)
Climatic zones
eq avg k ,
IV
Climatic zones
avg
Climatic zones
II
III
(3)
(4)
However, the obtained equivalent soil resistivity of twolayered structure of the ground cannot be used directly
for calculating the resistance of the grounding system.
In fact, the grounding system consists of different type
of the grounding conductors, which occupy different
position relative to the equivalent layers of the ground
of the calculated two-layered structure of the ground.
Therefore for the basic types of grounding conductors
and grounding system tables and drawings have been
developed. That makes it possible to calculate the
equivalent soil resistivity respecting the soil resistivity
of two-layered structure of the ground and relative
placement of the grounding conductors in it.
VES 2Ra ,
where
a l1 2
l 2 3
(5)
6. Conclusions
VES
R(
l12 4 l 223
2l 223
).
(6)
7. References
1.
2.
3.
264
4.
5.
6.
7.
8.
265
1. Introduction
As a result of climate change, windiness and thunder
storms in particular are considerably increasing. This in
turn leads to a higher fault probability. Together with
the climate change, the modern societys growing
dependence on uninterrupted power supply will greatly
influence the development of distribution networks and
organizations. Therefore, the requirement of reliable
and efficient protection system is increased. Earth faults
detection in isolated and compensated neutral medium
voltage distribution networks presents a challenge to
protection system due to small levels of fault currents
that can not be detected by conventional relays [1].
Although several techniques have been proposed, and
much progress has been made, a reliable and flexible
solution has still not been found. Transient-based
protective schemes utilize the transient component in
the fault signals, which provides fast information about
266
120
ir (t ) ia (t ) ib (t ) ic (t ) .
(1)
1
N
100
80
60
40
0.5 ms
1.0 ms
1.5 ms
2.0 ms
2.5 ms
20
500
1000
1500
(i
k 1
r ,k
)2 .
(2)
3500
4000
4500
5000
120
F=100 Hz
RMS=100%
100
80
60
40
0.5 ms
1.0 ms
1.5 ms
2.0 ms
2.5 ms
20
50
100
150
200
250
300
Frequency (Hz)
350
400
450
500
267
F5
Background Network
4 Feeders (198 km)
F3
Background
Network
198 km
Ir,back
F4
66/20 kV
/Y
F2
66/20 kV
/Y
Fault
Indicator
CB F1
Communication Link
7 km
7 km
Fault
Point
Ir
4 km 3
5 km
5 km
Wireless Communication
Remote Terminal Unit (RTU)
5 km
5 km
5 km
Ce
Faulted Feeder, 53 km
66 kV
Supply
B 10 km
A
66 kV
Supply
Fault
Point
I r ,back
r
back feeders
(3)
Ir
I r ,back
100 .
(4)
1
100 .
2n 3
(5)
268
1
100 .
m (n 1)
100
60
40
20
0
-20
-40
-60
(6)
ir(t) (F1)
ir(t) (F2)
ir(t) (F3)
ir(t) (F4)
ir(t) (F5)
80
10
15
Time (ms)
20
25
30
50
45
Ir (F1)
Ir (F2)
Ir (F3)
Ir (F4)
Ir (F5)
40
35
30
25
20
15
10
5
0
10
15
Time (ms)
20
25
30
160
140
W
W
W
W
W
120
(F1)
(F2)
(F3)
(F4)
(F5)
100
80
60
40
20
0
269
10
15
Time (ms)
20
25
30
160
Transient Fault Weights for Sections of F1(%)
140
6. Conclusions
120
100
W
W
W
W
W
W
W
W
W
80
60
40
20
0
10
15
Time (ms)
20
25
(AB)
(B1)
(BC)
(BD)
(C1)
(C2)
(C3)
(D1)
(D2)
30
7. Appendix
Background Network;
4 Feeders (2,3,4 & 5)
Total Length of 198 km
66/20 kV
Feeder 3 (50 km)
B1=10 km
BC=5 km
C1=5 km
C2=7 km
C3=4 km
AB=5 km
BD=7 km
D1=5 km
D2=5 km
D
140
SLG Fault
(Phase A)
120
Fault Resistance
251 km MV Distribution Network
100
80
60
Min.
Min.
Max.
Max.
W-faulty (Iso.)
W-faulty (Comp.)
W-healthy (Iso.)
W-healthy (Comp.)
Setting Weight = 60 %
1.1 m
40
20
0
0
10
8.1 m
1
10
10
10
Fault resistance (ohm)
10
10
270
1.1 m
1.1 m Raven
d= 10.11 mm
r= 0.536 /km
Soil resistivity: 250 m
8. References
9. Biographies
Mohamed F. Abdel-Fattah was born in Qualiobia, Egypt on
June 11, 1972. He is a lecturer and researcher at Zagazig
University, Egypt. He received his
B.Sc. with distinction and first class
honors, M.Sc. and Ph.D. degrees from
Zagazig University (Faculty of
Engineering, Department of Electrical
Power and Machines Engineering) in
1995, 2000 and 2006, respectively.
He is a member of IEEE. Currently, he
is a Post-Doctoral Fellow at Helsinki
University of Technology (TKK),
Department of Electrical Engineering, Espoo, Finland.
His main research interest is protection of electrical power
systems including fault diagnosis and location of high-voltage
transmission lines, transient/ultra-high-speed and wide-are
protective schemes and compensation of the high resistance
effect in earth faults. (e-mail: mohamed.abdel-fattah@tkk.fi ,
mohamed.f.abdelfattah@ieee.org ).
Matti Lehtonen (1959) was with VTT Energy, Espoo,
Finland from 1987 to 2003, and since 1999 has been a
professor at the Helsinki University of
Technology (TKK), where he is now
head of Power Systems and High
Voltage Engineering.
Matti Lehtonen received both his
Masters and Licentiate degrees in
Electrical Engineering from Helsinki
University of Technology in 1984 and
1989, respectively, and the Doctor of
Technology degree from Tampere
University of Technology in 1992.
His main activities include power system planning and
asset management, power system protection including earth
fault problems, harmonic related issues and applications of
information technology in distribution systems (e-mail:
matti.lehtonen@tkk.fi ).
271
SENSORS FOR ON-LINE PARTIAL DISCHARGE DETECTION IN COVEREDCONDUCTOR OVERHEAD DISTRIBUTION NETWORKS
G. Murtaza HASHMI*, Matti LEHTONEN*, Mikael NORDMAN**, Muzamir ESA*
*Power Systems and High Voltage Engineering, Helsinki University of Technology (TKK), Espoo, Finland
**Ensto Utility Networks, Porvoo, Finland
272
273
v rc (t ) = M
di (t )
dt
(1)
Ll
+
vrc(t)=Mdi(t)/dt
Cl
Zout vout(t)
Vout
Z out
(2)
= 2
V rc
s Ll Z out C l + s ( Ll + Rl Z out C l ) + ( Rl + Z out )
d2
d1
drc
i(t)
Conductor
Rogowski coil
turns
Internal return
wire loop
v out (t ) v rc (t ) M
di (t )
dt
(3)
Zout
vout(t)
Signal to the measuring instrument
i (t ) =
274
1
M
out (t ) dt
(4)
FDATool
vout (t)
Voltage signal obtained
from Rogowski coil
1
M
out (t ) dt
|u|
1
s
Abs
Integrator-2
Division
PD measured in
coulombs
q(t)
1
s
Integrator-1
200e-9
4. PD measurement methodology
q (t ) = i (t )dt =
(5)
7. References
6. Conclusions
277
Dielectric
fault
of
insulation
Problems caused by a
fault
Winding inter-turn short
circuit
Winding shorted to core
Leads shorted to tank
OLTC failure due to
short circuit
Gas generation
Electrical
connection
fault
Winding
failure
Tank
leakage,
mechanical
fault
Bushing
fault
OLTC
failure
Cooling
system
failure
Overheating
between
defective contacts at
terminal bushings
Winding open circuit
Winding deformation or
displacement
Oil leaks
High
impurities
or
moisture content
Cracked bushing
OLTC not functioning
Overheating of winding
and insulation
Paper degradation
279
Diagnostic test
Turns-ratio
measurement
Insulation
resistance
measurement
Insulation and
winding resistance
Dissolved
gas
analysis (DGA)
Thermography
Winding
resistance
measurement
Frequency
response analysis
Visual inspection
Oil quality
analysis, Tan delta
test
Tan delta,
inspection
Tap changer
functional tests
Thermography,
fans functional
test
Furfural analysis
Degree of
polymerization
( t )
n(t )
,
N t
(1)
(2)
Number of
outage
n(t)
3
3
3
5
2
4
2
3
6
3
(t),
per 1
year
0.03
0.03
0.03
0.05
0.02
0.04
0.02
0.03
0.06
0.03
p
0.97
0.97
0.97
0.95
0.98
0.96
0.98
0.97
0.94
0.97
Basic level
Level 1
Level 2
Level 3
Diagnostic method
Design analysis:
manufacturer;
winding
configuration;
materials etc
Operating envirnoment:
protection scheme;
ambient temperatures;
lightning level ect.
Maintenance:
loading and
overloads;
fault frequency;
previous repairs etc.
Dissolved gas analysis
Oil quality analysis
Infrared thermography
External inspections:
cooling system;
bushings;
tank;
lightning arresters ect.
Turns-ratio measurement
Winding resistance
measurement
Tan delta measurement
Insulation resistance
measurement
Internal inspection:
coils and clamp;
LTC switch;
oil level ect.
Furfural analysis
Degree of polymerization
Frequency response analysis
Partial discharge analysis
Use
Databas
e
Databas
e
Databas
e
6. References
1.
On-line
On-line
On-line
On-line
2.
3.
Off-line
Off-line
Off-line
Off-line
4.
Off-line
On-line
Off-line
Off-line
On-line
5.
6.
..
()
. // .: - ,
2002.
Transformer
fire.
The
internet:
http://www.electrical-contractor.net/forums/
ubbthreads.php/ topics/114034/1
M.Kosaki,
M.Ieda,
K.Itoh,
Y.Tsutsumi,
Application experience of expert system to
insulation diagnosis of power apparatuses in Japan
// CIGRE Session paper, 1990, Ref.No: 15/33-04.
Young Zaidey bin Yang Ghazali, Mohd Aizim bin
Talib, Hannah binti Ahmad Rosli. Condition
Assessment of Power Transformers in TNB
Distribution System and Determination of
Transformer Condition Index // 17th Conference on
the Electric Power Supply Industry, October 27-31,
2008.
Facilities Instructions, Standards, and Techniques
Vol. 3-31 // Hydroelectric Research and Technical
Services Group, US Dept. of the Interior, Bureau
of Reclamation, June 2003.
Franchek M.A., Woodcock D.J. Life cycle
considerations of loading transformers above
nameplate rating // Proceedings of the Sixty-Fifth
Annual International Conference of Doble Clients,
1998, Sec 8-10.1
cewk
i
Pytka stykowa
rdo napicia
Cewka pomiarowa
Cewka
pomiarowa
rejestrator
1. Introduction
282
Measurement coils
Vn S n d
potentiometer
(1)
Bars
t
Bn
dt
200
abs(u offset ) * f
i 1
the
t
offset
4* z *S
2 * 50 *103
max(abs ( I ))
(2)
min (u ) max(u )
1 200
1 200
(3)
283
c)
Styk bipolarny (skadowa osiowa)
600
500
400
B [mT]
r=12 mm
r=20 mm
300
r=30 mm
200
100
0
0
20
40
60
80
100
120
140
160
180
kt [stopnie]
d)
Styk bipolarny (skadowa osiowa)
450
400
350
B [mT]
300
r=12 mm
250
r=20 mm
200
r=30 mm
150
100
50
0
20
40
60
80
100
120
140
160
180
kt [stopnie]
1000
B [mT]
800
600
osiowa
400
promieniowa
200
obwodowa
0
12 14 16 18 20 22 24 26 28 30
r [mm]
b)
600
500
B [mT]
400
r=12 mm
r=20 mm
300
r=30 mm
200
B [mT]
100
0
0
20
40
60
80
100
120
140
160
180
kt [stopnie]
800
700
600
500
400
300
200
100
0
b)
osiowa
promieniowa
obwodowa
12 14 16 18 20
22 24 26 28 30
r [mm]
c)
400
350
300
r=12 mm
r=20 mm
200
r=30 mm
150
B [mT]
B [mT]
250
100
50
0
0
20
40
60
80
100
120
140
160
180
kt [stopnie]
800
700
600
500
400
300
200
100
0
osiowa
promieniowa
obwodowa
12 14 16 18 20
22 24 26 28 30
r [mm]
284
d)
6. Conclusions
The presented construction of the vacuum switches
construction does not exhaust all possibilities within this
range. The research has been conducted in the world all
the time on new interesting constructions that allow for
the increase of the rupturing capacity and rated voltage
often with the simultaneous decrease of the vacuum
chamber size. There are works on new constructions and
research methods of the switch systems done because
vacuum switches are ecological, have huge mechanical
and breaking capacity and allow for the switching off
the huge short-circuit currents.
The measurements results done on a chosen switch
system in the vacuum chambers of the switches shown
that there is a huge influence of the switches setting
towards themselves on magnetic field components
distribution. The measurement results of the magnetic
field distribution between switches used in the vacuum
chambers of the switches determine the original
approach in switch construction research. They also
allow for the evaluation of the magnetic field strength
value between switches and for the choice of the optimal
field distribution on the switch surface for achieving the
most profitable solution.
B [mT]
800
600
osiowa
400
promieniowa
200
obwodowa
0
12 14 16 18 20 22 24 26 28 30
r [mm]
B [mT]
500
ustawienie 1
400
ustawienie 2
300
ustawienie 3
200
ustawienie 4
7. References
100
0
12
14
16
18
20
22
24
26
28
30
r [mm]
B [mT]
250
ustawienie 1
200
ustawienie 2
150
ustawienie 3
100
ustawienie 4
50
0
12
14
16
18
20
22
24
26
28
30
r [mm]
0
B [mT]
1000
ustawienie 1
800
ustawienie 2
600
ustawienie 3
400
ustawienie 4
200
0
12
14
16
18
20
22
24
26
28
30
r [mm]
285
1. Introduction
The strategy of the protection of the high voltage
equipment insulation against overvoltages has been
changed when overvoltages arresters started being used
instead of dischargers. The valve dischargers can protect
the insulation of the equipment mounted in the
substation against violation in the zone of 100 120 m.
for its feature during its switching on, when the
discharge gap is broken down, to generate the returning
electromagnetic wave of the opposite direction which
dampers the overvoltages wave that comes from the
line.
During the process of the replacement of valve
discharges by overvoltages arresters in the electric
energy system, there might be a situation when a valve
discharger is left for the protection of the equipment
insulation, and the overvoltages arrester is connected to
the transformer. In such a combination of the protection
devices, a dangerous level of overvoltages can occur in
the case when the overvoltages arrester lowers the level
of overvoltages and the valve discharger cannot switch
on. In such a situation the insulation of the equipment,
mounted in the substation near the incoming lines, can
be damaged.
The aim of the work is to analyse the conditions of the
possible violation of the insulation, to evaluate the
change of the characteristics of the network elements
L5
42
6
Ci2
Ci1
43 L4
2
7 C
i3
Cr1
17
10
Ci6
45 L2
40
4
Cr2
12 Ci8
11
Ci5
Ci4
16
L3
44
40
40
40
1
24
30
86
Cui1
19
31
27
21
36
CT1
37
CT3
26
32
Cui2
86
38
CT2
Cui3
Cr4
21
Ci10
20
28
22
21
29
33
39
15
14
Ci9
25
23
40
Cr3 5
13
Ci7
18
46 L1
Cui4
34
40
Cu1
35
41
Cu2
CT4
286
C Cik
1
N,
2 0
u 3 x1 TU1 u 3 x1 ,
(1)
, q, 0;
= , b, c.
Current coordinates are calculated as follows:
L Lik 0 N F M ,
2
i 3 x1 TI1 i 3 x1 ,
here: TU1 ir TI1 inverted modal and TI . matrix values.
Phase coefficients and decrement of dissipation are
calculated as follows:
diag G TU1 LC TU ,
Surge resistance:
(2)
F = (Fik) ;
Fik
H1 1 1 1 12 H12 1 2 1 ;
2
21
1 22
arg k z2 k 02 ;
ik arctg
bi bk
hi hk
k0 0 0 ;
k z2 k 02
j 0
k1 j
(4)
rik D ik
(3)
1
;
1
287
288
2.
3.
4.
5. References
1.
2.
4. Conclusions
1.
3.
4.
289
and magnetic H M M 2 .
2
(1)
(2)
(3)
(4)
(5)
(6)
conversion
~
r
C R (T
E 2 ) w.
c
(7)
1
( D S ) am
~
r 2
E )
w R (T
c
.
e
(8)
R ( D S ) 1am
z 1 P1 , [h]
(9)
Activation temperature
~
Ta w / R P2 ; [K]
(10)
r
P3 ; [Km3J-1].
c
Second m aterial
(11)
800
700
T BE P e P T P E .
2
Resurce R, h
R R eTa
600
(12)
500
R given
400
R calculated
300
200
100
0
1
Test No.
First material
700
600
Resource R, h
500
400
R given
R calculated
300
200
100
0
1
Test No.
292
5.
4. Conclusions
1.
2.
1.
2.
3.
4.
, ..; , ..
. : , 1979.
304 .
, . .
.
:
, 1988. 128 c.
Kreuger, F. H. Partial discharges detection high
voltage equipment. London: Butterworth & Co.
Ltd. 1989, 192 p.
Aubin, J.; and Gervais, P. On-line monitoring of
key fault gases in power transformers and shunt
reactors. Sixty-second Annual International
Conference of Doble Clients, 1995. P. 122 135.
293
ts
(1)
there u t - the gap voltage; g - the gap length in m; Bconstant 0.42 for positive voltage and 0.5 for negative.
The leader propagation time tl could be calculated from
[4]:
u t
dl
k l u t
E str ,
g
l
dt
Introduction
0.5
,
u (t )
B
g
(2)
(3)
10 kV / cm .
for the pulse voltage could seek E str
The assumption is made for leader propagation time
calculation. The idea is that the leader stops develop,
than the voltage gradient inside the unabridged part of
gap drops below value E str . If the applied pulse voltages
rise time f is short, than gap flashover time become
294
2 x 10
a)
u(t)
u(t)[50%], s=100cm
-5
t,s
-5
t,s
10
1.5
b)
l(t)
0.5
0 -6
10
10
S=100 cm
S=150 cm
S=200 cm
2.2
2
x 10
1.8
Um, V
I max
1.6
1.4
1.2
Teup
Te
0.8
0.6
0
0.2
0.4
0.6
0.8
1.2
1.4
-5
t, s x 10
Ted
1
Tf
0
0
0.5
1.5
2.5
3.5
4
x 10
-5
295
layers on the moment of channel neutralization; Te electron spreading time in streamers, formed in the
middle of channel; Te d - electron spreading time in
(4)
(5)
I max
2 0 ve E str
E str
H l 10,8 H l ,
I max
4 0 ve E str
E str
H l 21,7 H l .
0,828
0,0834 I max
, I max 20 kA;
,
Tf
0,624
0,154 I max
, I max 20 kA.
4.
a max
a max
I max
Tf
I max
4 0 ve v E str
E str
6,5 109 A / s .(7)
Tf
z
a Req a Req 0 / 1
z
lightning
z
,
T f Req T f Req 0 1
z lightning
u u0 z I ,
there
(9)
I max
t / T f 2 e t / ,
1 t / T f 2
z lightning z line / 2
- wave
z zlightning / 2 .
(14)
(10)
Then the wave propagates approximately 100 m
distance, insulation breakdown after t 3,3 ns as it was
described in part 2 starts. After insulator breakdown,
tower and ground impedances should be involved in to
the equivalent lightning channel surge impedance:
surge
impedance.
The Heidler function [4] describes negative polarity
lightning current:
i t
2 z line / 2
(13)
(8)
2 0 ve v E str
E str
3,25 109 A / s (6)
z
I Req I max / 1
zlightning
(12)
(11)
296
(15)
5.
z 1 z lightning z x / z x z lightning ,
2.
(17)
3.
1.
(19)
2.
3.
4.
5.
6.
297
Conclusions
(18)
z 1 z lightning 2 / 1 z x1 zlightning ,
Acknowledgment
References
, . .; , . .; , .
.
.
, 75, . 9, 2005, p. 4451.
, . .
: .
, 76, . 12,
2006, p. 101-105.
, . . .
, , 1975.
Dudurych, I. M.; Gallagher, T. J. EMTP analysis of
the lightning performance of a HV transmission
line. IEE Proc. Generation Transmmission
Distribution. Vol. 150, No. 4, July 2003, p. 501506.
, . . .
, , 1973.
Shindo, T. Suzuki, T. A New Calculation Method
of Brekdown Voltage Time Characteristics of
Long Air Gaps. IEEE Trans. on Power Aparatus
and Systems. Vol. PAS-104, No. 6, June 1985,
p.1556-1563.
(1)
298
(3)
1 km
L1
UFLS
1 73 k m
L8
Nr.1
77MW
Nr.2
200MW
L2
P2-6 =1 5 0MW
144MW
380MW
45MW
P6 -5 =29 7MW
P2-6 =18 8M W
1 51 k m
62km
L6
22 k m
L5
UFLS
UFLS
70 km
P4 -5 =154 MW
144MW
144MW
L3
P3-6 =26 9M W
1 50 k m
P3-2 =1 8 7M W
1 km
Nr.3
787MW
P3 -4 =3 06MW
16 7 k m
L4
UFLS
25MW
144MW
299
Load maximum
MW
1200
1200
1580
51.6
Special setting of
UFLS, fact =
49.2Hz, tact = 0.5s
UFLS-1, fact =
48.8 47.2Hz, tact
= 0.5s
MW
36
36
47
50.8
MW
504
504
664
42
42
42
Frequency (Hz)
Parameters of
UFLS
50
49.2
48.4
47.6
46.8
39
40
41
42
43
44
45
46
47
48
49
50
Time (s)
200MW trip
130MW trip
300MW trip
40
40.5
(4)
41
41.5
42
42.5
43
Time (s)
200MW trip
130MW trip
300MW trip
300
(6)
n =P /PUFLS_act,
360
300
240
180
120
60
0
-60
39.5
40.5
41
41.5
42
42.5
43
Time (s)
200MW trip
130MW trip
300MW trip
(8)
4. Conclusions
Under-frequency load shedding operation during
emergency situation was analyzed using current
operating principle of UFLS automation in Baltic power
systems, as well as new suggested by authors active
power deficiency application principle.
Results showed that UFLS operating under new
principle more efficiently prevents frequency drop,
faster and more accurately restores frequency to rated
values, as well as prevents from frequency
overcorrection.
Proper selection of load to be shed can improve the
operation of UFLS automation preventing deep
frequency drop and avoiding over-frequency situation
after control.
51
50.5
Frequency (Hz)
40
50
49.5
49
48.5
48
47.5
39
40
41
42
43
44
45
5. References
46
Time (s)
200MW trip
300MW trip
130MW trip
301
3. Intelligent networks
303
304
305
306
9. Conclusions
1. Automation, as it is in process, will transform
distribution networks to digital (intelligent) distribution
networks where computer (microprocessor) network
will overlie power network.
2. Intelligent network vision is gradually turning to
emerging network technologies that will help convert
data to intelligence, intelligence to knowledge and
knowledge to intelligent business decisions.
3. The major difficulty related to future intelligent
substations and networks is referred to as interoperability of proprietary devices (different manufacturers
keep to different communication protocols).
4. Sometimes views of distribution system operators on
future intelligent networks are overcome by
apprehensions of complexity in introduction and
operation of intelligent technologies. Operators desire to
implement plug-and-play type technologies and are
apprehensive of superabundant information flows.
5. Distribution automation, particularly on feeders, has
proved its high effectiveness in urban areas in terms of
reliability performance (SAIDI minutes).
6. Automated meter reading systems are progressing to
advanced metering infrastructures (AMI) featuring twoway communications for switching control applications.
7. With the communication technology advance leading
to low-cost products, the automation of rural distribution
networks, including Distribution SCADA or microSCADA, is entering into cost-effective applications.
10. References
1.
307
generators,
distributed
1. Introduction
Power of generating units in electricity supply system
varies from few kilowatts to hundreds megawatts.
Those units can be a part of small, large and
interconnected electric power systems. Generating unit
connection to the grid can be different. Unit can be
connected straight to the load, both to the load and to
the distribution network, only to the distribution
network and to the transmission network. If unit is
connected to the electricity supply system, consumer
can be near, not far and far away. This variety is the
problem for exact definition identification for
distributed generators. Electricity supply structures and
energetic institutions give very different definitions for
distributed generators. The reason of different
distributed generators definitions is a diverse
understanding of traditional big power units. Mostly, if
unit or units are not recognized as a traditional, they are
called distributed generators. Sustaining definitions of
electricity supply structures and energetic institutions, it
is not possible to detect exact limits for the definition of
distributed generators. Definitions are selected by
subjective understanding. Objective definition must be
chosen by arguments and proofs.
308
T4
l(2n-6)
l6
l4
I2
I3
3
T(n-2)
I4
4
T(n-1)
In
n-2
l(2n-4)
T3
I1
l2
FL
T2
l(2n-3)
l(2n-5)
l5
I(n-1)
V1
l3
I(n-2)
l1
n-1
312
from all line terminals: I1, I2, , In, together with the
locally measured three-phase voltage phasor (V1) are
taken as the fault locator input signals (Fig. 1). In
natural way, these measurements are considered further
as synchronized.
The developed fault location algorithm consists of
(2n-3) subroutines designated for locating faults on
particular line sections (as depicted in Fig.1) where n
is the number of terminals. Since the position of a fault
is a random factor, the faulted line section is not known
in advance. Therefore, the valid subroutine (for the
faulted section) will be chosen using a special selection
procedure.
d1
I 1i
n-1
Fault type
a1
a2
a-g
a0
b-g
0.5 j0.5 3
0.5 + j0.5 3
c-g
0.5 + j0.5 3
0.5 j0.5 3
a-b, a-b-g
a-b-c, a-b-cg
1.5 + j0.5 3
1.5 j0.5 3
b-c, b-c-g
c-a, c-a-g
j 3
1.5 + j0.5 3
j 3
1.5 j0.5 3
0
0
Z iL(2n4)
I Ap = a1 I 11 + a 2 I 12 + a 0
(1 d1 ) Z iL1
n
Z iL(2n3)
where
d1
V 1p = a1V 11 + a 2 V 12 + a 0 V 10
T(n-1)
Z 1L1
T2
Z iL2
R1F
d1 Z iL1
where
a F1 , a F2 , a F0 share coefficients (Table2).
The i-th sequence component of the total fault current is
determined as a sum of the i-th sequence components of
currents from all line terminals (1, 2, , n-1, n):
(5)
I F1 = I 1i + I 2i + I 3i + ... + I (n -1)i + I ni
(2)
(3)
where
a1 , a 2 , a 0 - weighting coefficients (Table I);
313
a F1
a F2
a F0
a-g
b-g
1.5 + j1.5 3
c-g
1.5 j1.5 3
a-b
1.5 j0.5 3
b-c
j 3
c-a
1.5 j0.5 3
a-b-g
3 j 3
j 3
b-c-g
j2 3
j 3
c-a-g
3 j 3
j 3
V T(n -1)2 , V T(n -1)0 in the end tap point T(n-1) are
where
transf.
transf.
a-b-c
1.5 + j0.5 3
1.5 j0.5 3 *)
0
a-b-c-g
*) a F2 0 , however, negative sequence component is not
transf.
calculated:
1
Z iL1 T2
(1 d(2n 3) ) ZiL(2n3)
T(n-1)
Z iL2
d(2n3)
Z iL(2n4)
n-1
n2
1 I
transf.
j=
i =1, 2 , 3
ji
(10)
where:
Z iL(2n -3) impedance of line section L(2n-3) for the
d(2n 3) Z iL(2n3)
d1 =
I ni
(7)
V Tki
k 1
1 I
j=
i =1, 2 , 3
ji
(11)
where:
transf.
V T(k -1)i - calculated voltage at point (k-1),
symmetrical components.
transf.
transf.
transf.
The values of current I T(n -1)n1 , I T(n -1)n2 , I T(n -1)n0 , flowing
from tap point T(n-1) to station n in line section L(2n3) are calculated, as follows:
I T(n -1)ni =
transf.
n 1
1, 2 , 0
i=
ji
(12)
j =1
The fault loop equation for faults on last line section has
the following form:
V T(n -1)np (1 d (2n -3) ) Z 1L(2n -3) I T(n -1)np R(2n -3)F I F = 0 (13)
314
where
V T(n -1)np = a1V
transf.
T(n -1)1
+ a 2V
transf.
T(n -1)2
Z 0L(2n -3)
+ a 0V
transf.
T(n -1)0
transf.
Z 1L(2n -3)
(14)
transf.
where
Z 1L(2n -3) impedance of section L(2n-3) for the positive
Z 0 L(2n -3)
Symmetrical
seq.,
impedance of section L(2n-3) for the zero seq.
transf.
components
of
transf.
I TkT(k +1)1 ,
currents
transf.
k 1
transf.
1, 2 , 0
i=
ji
(20)
j =1
obtained.
d (2k1)
1
Z iL1
transf.
I TkT(k+
1)i
(1 d (2k k) ) Z iL(2k1)
Z iL(2k2)
Z iL2
transf.
Z iL(2n3) n
transf.
V T(k+
1)i
Z iL(2n4)
Z iL(2k)
V Tk0
Z iL1
transf.
I Tkki
Z iL(2k)
transf.
V T(k+
1)i
Z iL(2n3) n
(21)
where
V TkT(k +1)p = a1V Tk1 + a 2V Tk2 + a 0V Tk0
transf.
Z iL(2k)
Z iL(2n4)
transf.
k+1
n-1
transf.
transf.
i=
j =1, j k
(16)
ji
transf.
Z 0L(2k-2)
Z 1L(2k-2)
(18)
transf.
I Tkk0
(19)
transf.
Z 1L(k +1)
sequence component,
k number of the tap point.
After substituting (22)(23) into (21) and the respective
mathematical manipulations one obtains the searched
distance to fault d (2k 1) and fault resistance R(2k 1)F .
transf.
1, 2 , 3
Z 0L(k +1)
(22)
transf.
transf.
I Tkki =
transf.
d ( 2 k 2) Z iL(2k2)
d (2k2)
n-1
k+1
(1 d (2k 2) ) Z iL(2k 2)
Z iL2
315
(25)
4. ATP-EMTP evaluation
actual
5. Conclusion
316
4000
1.5
0.5
0.5
2000
2000
4000
1.5
0
20
40
60
80
100
120
20
40
Time (ms)
600
b
200
100
100
200
300
20
40
60
Time (ms)
80
100
120
200
200
400
20
80
100
120
3000
400
600
0
300
60
Time (ms)
40
60
80
100
120
2000
a
1000
1000
2000
3000
20
40
Time (ms)
60
Time (ms)
80
100
120
Fig. 6. Fault location on four-terminal network the example: input signals of the fault locator.
11
0.74
SUB_4 (10.41 )
0.72
Fault resistance ()
10.5
10
SUB_2 (9.83 )
9.5
8.5
0.66
0
10
20
30
40
50
60
10
20
30
40
50
60
Fig. 7. Fault location on four-terminal network the example: results of the subroutines SUB_2 (valid), SUB_4 (invalid).
6. References
1.
2.
3.
4.
5.
6.
317
318
Id
Im
Vdc
inet
abc
dq
Vq
Ws L Id +
Ws L Id -
PI
PI
Iqms
Idms -
PI
Vd
+
0
Vdc
d id
v dr = R id + L d t s L iq + Vdn
v = R i + L d iq L i + V
q
s
d
qn
qr
dt
(2)
The network represents the simple model with a threephase voltage source. The transmission line between the
network and the transformer is modelled with its
resistance, inductance and capacitance by equivalent
circuit p.
The transformer is considered to be linear (the magnetic
saturation has not been taken into account).
4. Grid side converter control
The role of the grid side converter is to keep the DClink voltage constant regardless to the magnitude and
direction of the rotor power. A vector control approach
is used, allowing independent control of the active and
reactive powers flowing between the supply and the
GSC.
Fig.2 shows a simplified representation of the grid
connected converter. AC-side series inductances
accounts for transformer leakage reactance, and series
resistances represents inverter and transformer
conduction losses
The model of the system is given by the following
equation:
P = vq iq and Q = v q id
(3)
v dc idc = v q iq and C
vdcC
dv dc
= idc im
dt
dvdc
= P Pm
dt
(4)
(5)
Va Van
1 0 0 ia
ia
d
R
1
ib = 0 1 0 ib + Vb Vbn (1)
dt
L
L
Vc Vcn
0 0 1 ic
ic
319
d dqs
m s . qds
Vdqs =R s Idqs +
dt
V =R I + d dqr m .
r
qdr
dqr s dqr
dt
(6)
(7)
1
0
g(x,t)= 0
0
1
0
0
1
0
Lr
0
0
0 (III.31) (13)
1
Lr
f(x,t)=
M
Rs.ads +Rs.cLrIdr + Vqs +sqs
Lss
Lss
1
M
M
Rr.bLrIdr + Vqs +Rr.cds +r LrIqr + Vds
Lss
Lss
Lr
1 R .bLI + M V +R.c LI + M V
Lr r r qr Lss ds r ds r r dr Lss qs
Vds = 0
Vqs s ds
(8)
dr =L r Idr + L . Vqs
s
s
=L I + M V
qr
ds
qr r
Ls .s
(9)
(14)
The general model of the machine is given by the
following equations:
Vds
V
qs =
Vdr
Vqr
d
M
Rs.ads Rs.cLrIdr + Vqs sqs + ds
Lss
dt
dqs
Lss
dt
dIdr
M
M
Rr.bLrIdr +L Vqs Rr.cds r LrIqr +L Vds +Lr dt
s s
s s
dIqr
M
M
Rb
.r LrIqr + Vds Rr.cqs +r LrIdr + Vqs +Lr
Lss
Lss
dt
d = ( I drref I dr )
q = ( I qrref I qr )
Where
dq = 0
f ( x, t ) + g ( x, t ) Vdq = 0
(17)
(16)
M
M
Rr.bLrIdr + Vqs +Rr.cds +r LrIqr + Vds
Lss
Lss
Ueqdq =
R.bLI + M V +R.c LI + M V
r r qr L ds r qs r r dr L qs
s s
s s
M
c=
LsLr
X =f(x,t)+g(x,t)U dq
Vdr and Vqr are the two control vectors of, they
(10)
1
1
, b=
,
With: a=
Ls
Lr
(15)
(11)
(12)
320
d-
Abbreviations
-a-
: Mechanical speed
turbine : is the turbine angular frequency
Vds , Vqs , Vdr , Vdr Tow-phase stator and rotor currents
voltage
ds , qs , dr , qr : Tow-phase stator and rotor flux
-b-
321
9. References
1.
2.
3.
4.
5.
6.
7.
322
Source
Frequen
cy range
Power range
efficiency
Advantages
Disadvantages
Comments
Cycloconvertors
16Hz
25Hz
Several
dozen kW
hundreds
of kW
=70. .80%
Simple desing of
the power part and
control system
Many valves
Power mains
50(60)
Hz
Several
dozen kW
hundreds
of kW
MW
Low frequency,
low power
density, nonsymmetrical
mains loading
Static magnetic
frequency
multipliers
150Hz
450Hz
Several
dozen kW
hundreds
of kW
=60...80%
Nature of the
operation, simple
design, low
investment, high
reliability, simple
operation
Constant
frequency, large
dimensions, heavy
weight
Electric
mechanical
generators
960 Hz
2000 Hz
8000 Hz
10000 Hz
Hundreds
of kW
2MW
= 60...80%
High reliability
High investment,
troublesome
operation, large
dimensions, heavy
weight, necessary
to build
foundations, long
start-up, water
cooling, noise
Semiconductors
statistical
converters
300 Hz
500 Hz
Severel kW
Hundreds
kW for MW
= 90...95%
Completely static
power and
frequency, high
reliability,
small overall
dimensions and
weight,
immediately ready
to work
Requirements
qualified persons
for operation,
Gentle
Basic sources of
power supply to
induction heaters
Lamp
generators
300 kHz
MHz
Tens of kW
hundreds kW
= 50... 70%
Completely static
operation, high
reliability,
immediately ready
to work
Low power,
constant
frequency, lamp
life 4000...6000h
Hz to as
323
a)
b)
capacitor voltage
T1 -
< <
>
< <
direct voltage
. This method is the most common, yet
it requires application of a regular-value source of direct
voltage.
b) the frequency method involves varying the frequency
f of transistor switching in a range comprising the
operation ranges, usually range 2 and 3, where
, for
instance, remains constant. It is not employed in
inverters working at high frequencies since the
switching losses rise during resonant tuning.
c) the method of voltage pulse band width regulation at
It involves an appropriate phase shift of
the signals controlling the transistors T1 and T4 in
relation to the signals for T2 and T3 at a constant
frequency of their turn-on.
324
a)
b)
a)
b)
a)
b)
<
inverter
5. Conclusion
This paper was intended to present methods of
supplying power to induction heating systems in the
range from medium to very high (MHz) frequencies. In
practice, their designs vary. The differences are
primarily determined by the intended use and operating
conditions of the inverters, and by the type of power
325
326
327
U0
U"0
U t U U 0 U e
or
U 0 U 0 e
U t U 0 I 0 xd U 0 I 0 xd U e
U 0 I 0 xd U 0 I 0 xd e
Td//0
I 0 xd xd e
At the subsequent voltage curve point (when the shortterm voltage component extinct) derived tangent to the
line that corresponds settled voltage U a direct axis
transient time constant T'd0 can be determined.
Extrapolating the curve of the exponential transient
voltage U' to the stoppage time the value of U'0 is
determined. According to it, the direct axis transient
inductive resistance x'd value is determined as follows:
xd
(1)
xd
(2)
(3)
where U0, I0, U are initial voltage and current and the
steady state voltage of the disconnected generator; U'0,
U0 initial values of transient and subtransient
voltages, T'd0, Td0 direct axis open circuit transient
and subtransient time constants, xd, x'd, xd direct axis
synchronous,
transient and subtransient inductive
reactances.
The variation in voltage when the generator that is
loaded with capacitive reactive load was switched off is
shown in fig. 1.
The initial value of the subtransient voltage U0 can be
expressed from the first voltage jump and a subtransient
inductive resistance value can be found:
U 0 U 0 I 0 xd ;
xd
U 0 U 0
.
I0
(6)
U0 U
.
I0
(7)
t
Td//0
U 0 U 0
.
I0
t
Td/0
or
U t U 0 I 0 xd I 0 xd xd e
T'd0
T"d0
t
Td/0
I 0 xd
I0 x'd
U'0
t
/
Td 0
I0 x"d
Parameters
xq
x'q
xq
T'q0
Tq0
(4)
Non-salient pole
generator
0.9 xd
1.5 x'd
xd
0.3 T'd0
Td0
Salient pole
generator
(0.6-0.7) xd
xd
0
Td0
(5)
328
234.5 B
Td0 Td0 B
;
234.5 N
rf rfB
Lq ( s ) Lq
(8)
234.5 N
.
234.5 B
q ( s ) Lq ( s ) iq ( s ) ,
(11)
Ld ( s)
Lq ( s)
Equivalent
Scheme of
d axis
id
b)
1 sT1d
Ld ( s ) Ld
(1 sTd ) (1 sTd )
(1 sTd0 ) (1 sTd0 )
(15)
Td Td
Td0 Td0
Tq Tq
Tq0 Tq0
Ld Ld ( j) Ld
(1 sTd0 ) (1 sTd0 )
(14)
(16)
(17)
(18)
G ( s) G0
Z q ( s) ra
Lq Lq ( j) Lq
Equivalent
Scheme of
q axis
iq
Z d ( s) ra
,
s
Ld Ld ( j) Ld
efd
a)
(1 sTq0 ) (1 sTq0 )
(9)
d ( s ) G ( s ) efd ( s ) Ld ( s ) id ( s ) ,
(1 sTq ) (1 sTq )
Lq Lq ( j) Lq
Td0
Tq
Tq0
(19)
(20)
Rotor poles of the hydro units that are made of the steel
sheets and the free currents closes through damping
windings of rotor transverse axis. Hydro generators that
is usually modelled with 2.1 model, now are designed
with one contour in the transverse axis. In case, when
damping windings time constants are smaller than the
excitation windings time constants, it is considered that
(12)
Td
(13)
329
Tq
Tq0
I fi Eqi D e
(21)
ln
C
(22)
I fM EqM
EqK EqM
I fK EqK
e
C EqK Eqt
(25)
(26)
Eq
1.2
1.0
Eq,t
I
I
If1.0(0)f1.0 f1.1
I fK EqK
If,,t
(24)
D
C E E
e qi qt ;
Eqi
P
.
d
t 0
dt
(23)
TJ
C ( E qi E qt )
If1.2 If
330
EqN I fN
I fN
.
I f, te1
(27)
xd
xd
(30)
(31)
EQN
(32)
xd
U N* cos N 2 U N* sin N I N* k qd xd 2
I fN* 1 k qd I dN* xd
0.
1 .1
.
K scr
(39)
Ratio
x"d / xd
Turbogenerator
0.10
Hidrogenerator
0.2
x'd / xd
0.16
0.3
TJ
(35)
J 2N 0.00274 GD2 n 2
,s
SN
SN
(40)
(38)
(34)
xq
2
xd 5 I fN
3I fN 2 I fN 2 .
2
2
xd I fN
* cos N sin N .
2
kqd
2
kqd sinN IfN* 1 kqd IdN*
IfN
* 1
.
2 2
2
2
2
(28)
2
2
xd I fN
* cos N ra * sin N .
2
1 kqd I dN
*
2 I dN2 *
2
xd 2kqd sinN IfN* 1 kqd IdN* IfN
* 1 0 . (36)
2
xd2 k qd
1 k qd
331
TJ TJPN cos N , s.
(41)
6. References
H TJ / 2 , s.
1.
(42)
5. Conclusions
2.
3.
4.
332
333
Gas, nm3/h
Efficiency, %
Diesel fuel,
kg/hour
Fuel consumption
width
21
MW
22
MW
23
MW
35
MW
415
MW
Overall
dimensions, m
length
Boiler plant
5,3
9,8
198
200
220
922
6,9
9,8
396
400
440
922
6,9
9,8
594
600
660
922
10,0
13,8
1485
1500
1650
922
24,3
14,3
5940
6000
6600
922
Q0
0,75 ,
Qh
(1)
Q0
January
February
March
April
May
June
July
August
September
October
November
December
Annually
(2)
24 Q g T k l k md
Northern Europe
0.8
1.5
2.6
3.4
4.2
5.0
4.4
4.0
3.3
2.1
1.2
0.8
2.8
Central Europe
1.7
3.2
3.6
4.7
5.3
5.9
6.0
5.3
4.4
3.3
2.1
1.7
3.9
(3)
334
P2
Solar
collector
Tank
P1
Evaporator
Absorber
Pump of cooling
agents recirculation
P2
Tank
P1
Solar
collector
Evaporator
Heater
Heat
storage
drums
system
Absorber
Pump of cooling
agents recirculation
Ground level
Fig. 2. Single-stage absorption refrigerating unit using hot water from the solar collector with heat storage drums system in the seal
under ground; 1 polyethylene tubes in the seal
335
T21
T3
T11
Biogas
T1
T2
B1
4
5
K1
Fig. 4. Scheme of installation, which uses biogas and heat of sewage for the heat supply system: 1- methane tank; 2 solar
collectors; 3 block of biogas cleaning; 4 internal combustion engine; 5- heat pump plant; 6 condenser; 7 evaporator; T1
water from cooling system of internal combustion engine; T11 water in heating system, t = 80-85 0C; T21 water from the heating
system, t = 55-60 0C; T3 hot water in hot water supply system, t = 50-55 0C; B1 cool water from the running water supply
system; K1 sewage
3.
4. Conclusions
1.
2.
3.
4.
5.
6.
7.
5. References
1.
2.
..
//
. --,
3(19), 2005. p.p. 20-24.
www.wartsila.com (called on 02.16.2009).
8.
336
WRTSIL.W-P 09.2004/Bock/s
Office Waasa Graphics.
Lundh M. Dalenbck J.-O. Solar heated
residential area Anneberg initial evaluation //
LTMA mokslo darbai, 2. Klaipda,
Klaipdos universitetas, 2005. p.p. 142-153.
Kytra S. Atsinaujinantys energijos altiniai.
Kaunas, KTU, 2007. 140 p.
..
. 2267060.
., . c
.
.
2- (2-3
2000 . .-), 2. ., 2000. p.p.
405-410.
Rabinovych M.D. Renewable energy sources
for sustainable development of historical cities.
Environment Protection Engineering, vol.32,
1, 2006. Wroclaw, 2006. p.p.81-86.
Keywords: cooling towers with natural draught, nontraditional energy sources, wind turbines, electric
energy.
1. Introduction
Ecological, economical and energy crisis currently
taking place in the world makes us to look for new
possibilities and new approaches in order to find out
new energy sources for heat and power generation. First
of all, we look for usage of renewable energy sources,
which are universally available, free from payment and
ecologically clean.
Some possible ways to use the cover of acting or
suspended cooling tower for the installation of small
scale wind turbines and (or) photovoltaic modules are
described in this paper. The suggestions discussed
bellow are based on the analysis of peculiarities in
aerodynamic and heat situation in the cooling towers of
various types.
General view and scheme of counter-flow cooling
towers with natural draught are shown at Fig. 1 and Fig.
2.
Reference
[7]
[4]
[6]
Gage
diameter, m.
125
Height of
offtake
windows, m.
10
Height of
tower, m.
150
120
[13]
104,5 64,4
80
82
97
108
90
102
110
130
On the
upper level
of nozzles
At the
outgoing
opening
Height of offtake
windows, m.
Diameter of tower,
m.
-576-65
40
25,2
17,6
3,3
Steel
-725-66
39
32,5
21,8
3,55
//
-900-66
43,4
31,2
21,8
3,3
//
-1520-70
53,5
47,6
25,6
3,3
Ferroconcrete
-1200-70
48,3
40
26
3,3
Steel
-2100-70
64,8
52,2
33
3,8
//
-2600-70
71
58,2
37.9
4,3
//
-3200-73
81
64,4
40,4
//
-4000-71
91,5
72
41,6
5,5
//
No. of
project
Height of the
cooling tower, m.
Material of
tower frame
338
339
Q1 Q1T Q1C ;
(1)
Q1T K g F 10 3 K g F t H t B 10 3 , (2)
1
B
(3)
Q1C K g FC tC 10 3 ,
(5)
tB
150 m
the
75 m
after
W0
calculated
Q1
was
tH
W/(m2K)
125 m
tC 90 0 200 = 70 0;
340
Altitude, m
m 0,6 V 0,77 .
(8)
where
200
2,6
(7)
100
2,1
(6)
z
V1 V2 1 ,
z2
40
1,55
V2 V1
10
1,0
(9)
150 m
1
k k V 2
Pv 1 2
K,
2
2
(10)
Wind
341
2.
3.
1
4.
Wind
5.
6.
Parameters
Wind speed, m/s
Capacity, W/m2
Wind speed, m/s
Capacity, W/m2
Wind speed, m/s
Capacity, W/m2
7.
Values of parameters
depending on the localitys
surface roughness
0
1
2
3
7,6
5,3
4,7
3,7
456
181
119
58
8,9
7,3
6,6
5,8
702
400
309
204
9,5
8,5
7,8
6,9
890
595
462
319
8.
9.
4. Conclusions
1.
2.
3.
10.
11.
12.
13.
14.
5. References
1.
15.
.., ..,
.. // .:
, 1976. 210 p.
342
..
//
. .. , , 236, 2000.
p.p. 242-247.
Israel Ziga-Gonzles. Modeling heat and
mass transfer in cooling towers with natural
convection. Summary of dissertation. Czech
Technical University in Prague, faculty of
mechanical engineering. Prague. 2005. 23 p.
Hoffmann J.E., Krger D.G. Analysis of heat
mass and momentum transfer in the rain zone
of a natural draft counter-flow cooling tower.//
7th COOLING TOWER AND SPRAYING
POND SYMPOSIUM. Leningrad, USSR, May
29-June 2, 1990. p.p. A10-(18).
Grange J.L. A computer code for particular
operating conditions of wet cooling tower// 7th
COOLING TOWER AND SPRAYING POND
SYMPOSIUM. Leningrad, USSR, May 29June 2, 1990. p.p. 8-(111).
.., ..
. . //
.. . .:
. 1998. 278 p.
.., ..
2005.
Kugeleviius J.A., Kuprys A., Kugeleviius J.
Aplinkosaugini ir energetini energijos
gamybos rodikli vertinimai // Journal of
Environmental Engineering and Landscape
Management, 2005, Vol. XIII, 4. p.p.
192a-199a.
Katinas V., Tumosa A. The possible usage of
wind energy in Lithuania. Energy Agency,
Directorate of the Energy Saving Programme,
Kaunas, 1996. 27 p.
..
. 3-
., . .:
, 1989. 223 p.
http://www.chv-praha.cz
(called
on
02.02.2009)
http://www.uniserv.com.pl
(called
on
02.16.2009)
.., ..
: // .:
, 2007. 52 p.
Kytra S. Atsinaujinantys energijos altiniai.
Kaunas, 2007. 140 p.
..,
..
: ,
, . .: ,
2005. 264 p.
1. Introduction
Industry
2000
2003
2004
2005
133,5
135,6
135,0
135,5
6155
5842
6391
6369
2000
135
0,05
2003
1130
163
2004
786
8,4
100
Temperature of
recycled water 0
Cooling
tower in
Cooling
tower
out
4045
25-28
40-45
25-30
Supply of
recycled
water in
major
enterprises,
thousands
m3/hour
100
2005
1090
0,5
343
Electric drive
capacity,
kW
11
30
75
200
Rotation
frequency,
min-1
365
178
178
110
250
91
250
800
3
3
4
127
65
1410
950
730
1
38 m
22 m
2
20 m
40 m
Fig.3. .Scheme of positioning of solar collectors on the wall of
CK-400 single fan cooling tower: 1 the cooling tower wall; 2
PV modules
344
2
55
1230
1230
33
33
Fig. 5 Scheme of solar module operating activity with symmetric parabola-cylinder concentrator in stationary position under
various Sun height: a) midday, June 1st , b) midday March 16 or 8.30 and 16.30, Moscow latitude; 1 symmetric concentrating
reflector; 2 double-sided photoelectric module; 3 solar rays path.
(1)
145
plain
angle
h0
37
plain
angle
1230
h = h0 + (90 - ),
(2)
345
2,5
I,
1,5
hm
145
0,5
h0
37
0
0
10
1230
U,
4. Conclusions
1.
Fig. 7. Change of the concentrators middle plain position in
relation to horizon plain under experimental modeling of
various sun heights
2.
1.
2.
3.
346
..
// .2004,
10, C. 31-33.
Pushnov A., Berengarten M., Ryabushenko A.
Influence of Aerodynamic Conditions I
Cooling Tower or Losses of Water in System
of Defensive Water in System of Defensive
Water Supply. May 22-23, 2008. The 7-th
International
Conference
Environmental
Engineering, VGTU Press Technika. Vilnius, 2008. - P.290-296.
.., ..
: /
.. . - .:
, 1998. - 376.
Grid
Wind
CS
Load
ISC
AC/DC
WT
SG
AC/DC
1. Introduction
RFB
Fig. 1. Wind power plant with ISC and RFB electrical energy
conversion structural scheme
R (Load)
ISC
AC/DC
Ired
RFB
Grid
Ug
0,20
Uo
0,15
Io
0,10
Io=f(Uo)
0,05
0,00
c
1
Breaker1
Out1
20
40
60
80
100
Conn 2
Clock1
120
Uo (V)
Out2
W ind power
plant
Super
capacitor
Switch
Conn 1
Conn 3
Step2
WPP internal
resistance
Conn 4
ISC
1,00
0,80
Clock2
Switch1
Breaker3
Breaker4
0,60
c
1
c
1
Step1
U
c
Step3
I
0,40
Step4
Breaker2
0,20
Grid1
R
(Load)
Out1
Conn 2
Out2
RFB
Conn 1
Conn 3
Conn 4
AC/DC
0,00
0
Grid
20
40
60
80
100
RFB
internal
resistance
120
Uo (V)
348
Uo=60V
0,8
0,4
0,2
0,0
0
200
300
400
500
1,4
1,2
Uo=24V
Ug (V)
500
0,6
Uo=5V
0,4
100
R ()
1,6
1,0
0,8
Uo=5V
0,6
Uo=60V
IR (A)
1,8
Uo=24V
1,0
400
0,2
300
0,0
0
100
200
300
400
500
R ()
200
100
0
20
40
60
80
100
120
Uo (V)
349
350
5. Conclusions
1. Research results confirmed voltage system
(U=const., I=var.) electrical energy conversion
possibility to current system (U=var., I=const.)
electrical energy. Converter operates at inverter
system converter mode.
2. Wind energy generated voltage value variation has
low effect on ISC efficiency factor at wind power
plant with RFB.
3. Stochastic character wind power can be steadily
transmitted to grid after ISC and grid frequency f
and phase synchronization when grid voltage Ug
is equal to standard grid voltage RMS values
(Ug=230 V; 400 V). Electrical energy steadily is
transmitted to grid when Uo variation range is -5
100%.
4. ISC can operate autonomously. In this case RFB
have to be connected to ISC for energy storage.
RFB set autonomous load AC voltage effective
value.
5. Energy conversion processes results shows that
RFB connected to ISC input could determine
7. References
1.
2.
3.
4.
5.
351
1. Introduction
180
160
140
120
100
80
60
40
20
0
EH,
kWh/m 2
10 11 12
Months
352
P2
E, W/m 2
T2
T3
DHW HEAT
EXCHANGER
CHSS
outlet
DHW users
P2
Preheated water
T1
SC
CU
CV
P1
P1
CHSS
inlet
Fig. 2. Structure of hybrid water heating system based on solar collectors and centralized heat supply grid
353
where
dimensions,
c = 4190 the specific heat of water, J/kg .K,
E, kWh/m 2
100
80
E S C = E1h S SC ,
60
40
where
20
2
9 10 11 12
Months
array, m2.
0
1
(2)
S SC
k d c V (T1 T2 )
E1h
(3)
(1)
354
6. Conclusions
2.
3.
4.
T
Ton = Tmax .1 1
T1N
where
E
,
E ST
(4)
5.
6.
the temperatures, C,
T1 the instantaneous temperature in the tank,
C,
T1N the rated temperature in the tank, C,
7.
Ton =
E (60 T1 )
.
2000
1.
2.
(5)
3.
4.
355
2. Object of research
The suggested experimental hybrid water heating
system is intended to install in condominium (Fig. 1),
which is located in Kaunas, Lithuania. The main
technical characteristics of this condominium are given
in Table 1.
Presently domestic hot water (DHW) is being prepared
in the condominium only by means of heat exchanger
running on the centralized heat supply grid. The main
disadvantages of the existing DHW system are
considerable price for water heating and blackouts of
hot water supply in summertime when the system can
356
450
V,m3
400
381
350
302
300
290
290
255
273
227
250
200
150
85
100
50
0
1
10
11
12
13
10 11 12 Months
Nr.
1
2
3
4
5
6
7
Parameters of house
Units
Values
Number of apartments
75
Number of stairwells
7
Number of stories
6
Number of dwellers
270
Total heated area
m2
6 750
Total heated space
m3
16 500
Length x width x height
m
60 x 20 x14,60
Annual consumption of
m3
3 428
DHW
Average monthly
m3
285,67
consumption of DHW
Average daily
m3
9,39
consumption of DHW
Annual consumption of
MWh
175
heat for DHW
Average amount of heat
energy for production of
kWh/m3
51
1 m3 DHW, heat losses
and 75 serpentine pipes
Tariff for water heating
Lt/m3
12,19
Rate of Litas against Euro
is constant and equal to 3,4528 LTL for 1
o
o
Nr.
1
2
3
4
5
6
7
Parameters of heat
Units
Values
exchanger P 0,15 - 200
Producer Joint Stock Company
UAB Pergals koncernas
Rated thermal capacity
kW
200
Rated water flow
m3/min
0,15
Operating pressure
MPa
0,6
Testing pressure
MPa
1,25
Maximal operating
C
120
temperature
Length x width x height
cm
37 x 28 x 93
Weight
kg
280
357
Fig. 3. Block diagram of control process for the hybrid water heating system
358
359
Switching delays;
Pump blocking protection;
Operating hours / limiting value;
Command execution check;
Malfunction catch / malfunction handling.
360
3.
4.
5.
6.
Fig. 9. View of trend curve on CU screen for actual difference
of temperatures
7.
5. Discussion
8.
9.
7. References
1.
2.
3.
6. Conclusions
1.
2.
4.
5.
361
2. Research object
Small scale efficient wind turbines [3] of sufficient
capacity can be used for the heating purposes as well
[4]. Simplified scheme of wind turbines power
conversion system used for all domestic energy needs is
shown in Fig. 1. The wind turbine (WT) has a
permanent magnet synchronous generator SG. As it is
shown in Fig. 1, the synchronous generator SG is
connected to the electric grid over the rectifier, the
shorting circuit, the inverter, the transformer and the
grids automatic switch 1Q,F. Shorting circuit consist of
shorting transistor Vt, reactor Ld and disjunctive diode
Vat. Wind energy converted by the WT is used for the
domestic electrical appliances, for charging the battery,
for the domestic hot water (DHW) preparation and
362
Fig. 1. Simplified scheme of wind turbines power conversion system used for the domestic energy needs
U 2 p U1p
UB
3 ks
(1)
IR
Pi
,
kis k m g
(2)
Pi Ptm
U d2
,
Rh
(3)
363
(4)
364
Ld
s2
id
rd
L1
s3
iL
r1
Rh
eg0
ih
PMSG
s4
s1
+
b
ii Lf
rf
iv
rk
rb
Vat1
ik
ib
iv
Ck Eb
Vat2
Lt
ii
rt
ei
SIC
dii 1 u e r i ;
dt Li ab i i i
du
1
c U c0
ik ;
dt
C
k
u e u u r i i ;
g0
g
K
1
h
d
d
d ih id
u g 2 r1 ih id 2 L1
;
dt
m xK
ih id ;
u K
2
U c 0 k s Eim ;
u u r i ;
c
k k
ab
1
ib r u ab Eb ;
b
ik id ib ii .
(5)
365
d g
1
Mt M g M f ;
J
dt
Ptm
M t ;
g
M g k m iL 2 ;
3 3
M k .
F
g
f
6
5
Pt 0.0461 v w
2.1317 v w
36.5307 v w4
3
271.5777 v w
731.0533 v w2 531.6232 v w
(7)
8.5174.
(6)
5. Results of research
As it was shown in Fig. 2, according to the proposed
algorithm the system can supply the wind energy
converted into the power to the water heating unit or
into the electric grid, if the heat storage tank is fully
charged. If generator produces more energy than heating
units require, the surplus power can be supplied to the
electric grid as well.
Quality of the load control in the proposed system was
researched by means of mathematical simulation.
Results of the research are presented below in Fig. 5 and
Fig.6.
366
2.
3.
4.
5.
1.
2.
3.
4.
5.
6.
7.
367
368
Fig. 1. Simplified scheme of power supply from park of small-scale WTs into EG over the one inverter
369
Ld1
Ld2
rd1
Ldn
rd2
Lg1
S1
rf
Lg2
S2
Vat
Lgn
ign
rg1
rg2
rgn
ug2
Sn
+
ugn
SG2
Lf
ii
ig2
SG1
Si
rdn
ig1
ug1
b
SGn
rk
Vat
ik
iv
Ck
Lt
ii
rt
ei
SIU
When the wind turbines generators ug1, rg1 ugn, rgn are
switched into the scheme of short circuit by means of
the switches S1 Sn (in case of AC the rectifier is
necessary), the power is being stored in the reactors Ld1
Ldn. When the switches S1 Sn are switched out, the
power flows into the storing capacity Ck. If the Si is
switched on, the inverter I supplies the power into
electric grid 7. Pulse width modulation realised by
means of the switch Si is used in order to adjust voltages
of the wind generators and the electric grid. The grid
voltage is used for the modulation. In this way
operation of converter is being synchronized with the
electric grid. Carrier frequency of the modulated signal
usually makes up from few to 20 kHz. Operation of the
switches Sn can be synchronized with the electric grid or
not. In our case switches S1 Sn operate independently
from the electric grid voltage.
Mathematical description of the systems power circuits
was made up with reference to the equivalent scheme of
the proposed system (Fig.2).
Every generator (Fig. 2) can operate in the two modes:
shorting mode when the switches S1 Sn are switched
on, and power inverting mode when the switches S1
Sn are switched out. When the generators are operating
in the shorting mode (S1 Sn are switched on), their
power circuits are not related electrically (independent).
Then currents of the generatorspower circuits can be
deduced from the following independent equations:
dign
dt
1
u gn rn ign .
Ln
dig1 1
dt L u g1 uab r1 ig1 ;
1
dig 2 1
u g 2 uab r2 ig 2 ;
L2
dt
1
dign
u gn uab rn ign ;
dt
Ln
di
1
i uab ei ri ii ;
dt Li
U k E ;
s
im
c0
du
1
c
U c0
ik ;
dt
Ck
u u r i ;
ab n c k k
i i i .
k i 1 gi i
(1)
370
(2)
d gn
dt
1
Jn
Ptn
Tgn ign ,
i
gn n
Pt 16.517 10 5 v w5 0.01881 v w4
0.6797 v w3 9.1408 v w2 27.8 v w 4.
(3)
Pt 16.517 10 5 v w5 0.01881 v w4
0.6797 v w3 9.1408 v w2 27.8 v w 4.
(5)
(4)
371
372
3.
4.
5.
6.
7.
1.
2.
3.
4.
5.
6.
5. Conclusions
1.
2.
7.
373
374
K ew
C
; /W
P10
K wm
K te
(1)
P10
; W/kg
m
K wm P10 2
; W2/m2
K ew S C
(2)
16,54
Kew, E/W
12,37
1,21
1,29
Proven 15
CB-3,1/200
1,34
13,27
1,47
1,45
Superwind
Fig. 1. The 5 lowest and the 3 highest prices of small HAWTs for 1 watt
Kwk, W/kg
45,00
40,00
35,00
30,00
25,00
20,00
15,00
10,00
5,00
0,00
42,80
40,6
38,3
36,8
35,4
8,2
(4)
18
16
14
12
10
8
6
4
2
0
(3)
K wk
P10
; W/m2
S
Hannevind
30
Airdolphin
Aerocraft
120
7,4
5,0
Fig.2. The 5 highest and the 3 lowest capacities of small HAWT per 1 kg of their weight (without tower)
375
350,00
Kwm, W/m2
308,64
300,00
256,64
255,56
250,00
244,09
240,33
200,00
150,00
100,00
87,50
85,71
72,43
Superwind
350
GW-11
Aimpair 300
50,00
0,00
Keith_2K
ARI-450
Aerocraft
1002
Airdolphin
Proven 15
Fig. 3. The 5 highest and the 3 lowest capacities of small HAWT per 1 m2 of their swept area
250
Kte, W 2/Em2
198,62
200
173,35
150,08
150
138,83
136,59
100
50
9,29
6,92
5,29
Ampair 300
Superwind
350
0
Proven 15
ICO-GE 750
ARI-450
Fig. 4. The 5 highest and the 3 lowest technical-econimical efficiencies of the 45 researched small HAWT
Name of
HAWT
Proven 15
ICO-GE 750
ARI-450
WindMax-H5
WindMax-H8
376
Name of parameter
Rated
Rated
Swept
Price,
wind
power,
area,
speed,
2
W
m
m/s
15 000
9,8
63,58
18 511
750
12
2,54
789
450
11,8
1,54
495
485
13
1,54
450
775
12,5
2,54
750
125,00
Kew, E/W
120
97,22
100
89,87
80
60
40
20
2,63
3,00
3,09
3,67
3,78
0
Vertikalrotor Windspire FSW Gyro-5 FSW Gyro FSW Gyro10 kW
Extreme
WS-4B
WS-2B
WS-0.3C
Fig. 5. The 5 lowest and the 3 highest prices of small VAWTs for 1 watt
Kwk, W/kg
18
16,67
16
14
12
10
8
7,60
6,62
4,23
3,82
3,20
2,86
0,93
2
0
WindGen Vertikalrotor
1000
0.3
Turby
WS-0.3C
Fig. 6. The 5 highest and the 3 lowest capacities of small VAWT per 1 kg of their weight (without tower)
377
Kwm, W/m 2
450
384,62
400
350
293,71
300
250
181,52
200
169,81
171,43
121,13
150
72,00
72,00
100
50
W
S2B
W
S4B
1
FS
W
G
yr
W
in
ds
pi
re
Tu
rb
y
kW
o5
G
yr
FS
W
FS
W
G
yr
W
in
dG
oEx
tre
en
m
e
10
00
Fig.7. The 5 highest and the 3 lowest capacities of small VAWT per 1 m2 of their swept area
Kte, W2/Em2
90,00
80,00
70,00
60,00
50,00
40,00
30,00
20,00
10,00
0,00
79,60
43,95
W
S0.
3C
W
S4B
0,80
0,74
FS
W
yr
o
10
kW
en
10
00
W
in
dG
0.
3
er
tik
al
ro
to
r
FS
W
e5
yr
o5
1,07
W
S2B
49,37
yr
oEx
tre
m
G
FS
W
57,79
58,74
Fig. 8. The 5 highest and the 3 lowest technical-econimical efficiencies of the 14 researched small VAWT
Name of
VAWT
FSW GyroExtreme 5
FSW Gyro-5
Verticalrotor 0.3
WindGen1000
FSW Gyro-10
Name of parameter
Rated
Rated
Swept
Price,
wind
power,
area,
speed,
W
m2
m/s
500
12,5
1,43
1 550
500
300
1 000
10 000
10,72
12,5
12,5
11,62
3,03
2,5
1,56
37,2
1 700
998
4 675
22 000
378
7. Conclusions
1.
2.
6. Comparison of HAWTs against VAWTs
HAWT-5b
VAWT-5b
4.
38,8
5.
20
10
3.
3,22
7,79
1,35
Kew
Kwk
8. References
300
250
200
150
100
50
0
261,05 240,22
1.
2.
HAWT-5b
3.
VAWT-5b
159,49
4.
57,89
Kwm
Kte
5.
6.
7.
8.
379
http://www.allsmallwindturbines.com (called
on 03.02.2009).
Small Wind Turbine Global Market Study
2008. Published by the American Wind Energy
Association. AWEA, June 2008. 26 p.
Christensen K. Catalogue of Small Windmills
from 5 W to 50 kW, 2006. Nordic Folkecenter
for Renewable Energy. 105 p.
Adomaviius V. Sizing of wind-based energy
system for homesteads power and heat.
Proceedings of International conference
Electrical and Control Technologies 2007.
KTU, Kaunas, 2007. pp. 157-162.
Adomaviius V. The concept of wind energy
consumption for heat and power in the
dwelling houses. Proceedings of the 12th
International conference Technical and
Technological Progress in Agriculture.
Raudondvaris: IAE at LUA, 20-21 September,
2007. pp. 226-230.
Martinez F., Herrero L.C., Santiago P. G. and
Gonzales J. M. Analysis of the efficiency
improvement in small wind turbines when
speed
is
controlled.
1-4244-07559/07/$20.002007 IEEE.
Murray D. Big Opportunities for the Small
Wind Industry. Go Green Energy, LLC. 2008.
5 p. (http://www.gogreenenergyllc.com
called on 02.16.2009).
Migliore P., van Dam J. and Huskey A.
Acoustic tests of small wind turbines. AIAA2004-1185. pp. 14.
380
381
382
40
P [MW]
30
20
10
0
0
10
12
14
Taime [s]
P [MW]
50
40
30
20
10
0
Taime [s]
Generator power 40MW and
area av erage wine speed
Generator power 25MW and
area av erage wine speed
Generator power 17MW and
area av erage wine speed
Generator power 40MW and
10
12
14
f [HZ]
50
48
46
0
Taime [s]
Generator power 40MW and
area av erage wine speed
Generator power 25MW and
area av erage wine speed
Generator power 17MW and
area av erage wine speed
Generator power 40MW and
10
12
14
383
140
120
P [MW]
100
80
60
P [MW]
40
20
0
0
Taime [s]
Voltage in Maeikiai power plant point
10
12
14
Taime [s]
140
120
P [MW]
100
80
60
40
20
0
0
Taime [s]
Voltage in Maeikiai power plant point
Voltage in Klaipda port point
Voltage in Bting oil terminal point
P [MW]
100
80
60
40
20
0
0
Taime [s]
Voltage in Maeikiai power plant point
Voltage in Klaipda port point
Voltage in Bting oil terminal point
Fig. 15. Wind power park generated power during the period
384
Fig. 16. The concept of power grid with wind power plants,
fuel cells and energy storage system
2.
2.
3.
4.
385
Deveikis, Tomas,
Nevardauskas,
Enrikas
Stanioniene
Renata.
Investigation
of
electromagnetic compatibility of wind power station
and electrical network user // EMD 2008 :
proceedings of the XVIII International Conference
on Electromagnetic Disturbances, September 25-26
2008, Vilnius, Lithuania. - ISSN 1822-3249. Vilnius. - 2008, p. 155-160
Assessment of Islanded Operation of Distribution
Networks and Measures for Protection. ETSU
Report ETSU K/EL/00235/REP. 2001.
www.enercon.de
www.vestas.com