TJ 13 2019 3
TJ 13 2019 3
Editorial Office:
Sveučilište Sjever – Tehnički glasnik
Sveučilišni centar Varaždin
104. brigade 3, 42000 Varaždin, Hrvatska
Tel. ++385 42 493 328, Fax.++385 42 493 333
E-mail: tehnickiglasnik@unin.hr
https://tehnickiglasnik.unin.hr
https://www.unin.hr/djelatnost/izdavastvo/tehnicki-glasnik/
https://hrcak.srce.hr/tehnickiglasnik
Council of Journal:
Marin MILKOVIĆ, Chairman; Anica HUNJET, Member; Goran KOZINA, Member; Mario TOMIŠA, Member;
Vlado TROPŠA, Member; Damir VUSIĆ, Member; Milan KLJAJIN, Member; Anatolii KOVROV, Member
Editorial Board:
Chairman Damir VUSIĆ (1), Milan KLJAJIN (2)/(1), Marin MILKOVIĆ (1), Krešimir BUNTAK (1), Anica HUNJET (1), Živko KONDIĆ (1), Goran KOZINA (1), Ljudevit KRPAN (1), Krunoslav
HAJDEK (1), Marko STOJIĆ (1), Božo SOLDO (1), Mario TOMIŠA (1), Vlado TROPŠA (1), Vinko VIŠNJIĆ (1), Duško PAVLETIĆ (5), Branimir PAVKOVIĆ (5), Mile MATIJEVIĆ (3), Damir
MODRIĆ (3), Nikola MRVAC (3), Klaudio PAP (3), Ivana ŽILJAK STANIMIROVIĆ (3), Krešimir GRILEC (6), Biserka RUNJE (6), Predrag ĆOSIĆ (6), Sara HAVRLIŠAN (2), Dražan KOZAK
(2), Roberto LUJIĆ (2), Leon MAGLIĆ (2), Ivan SAMARDŽIĆ (2), Antun STOIĆ (2), Katica ŠIMUNOVIĆ (2), Goran ŠIMUNOVIĆ (2), Ladislav LAZIĆ (7), Ante ČIKIĆ (1)/(2),
Darko DUKIĆ (9), Gordana DUKIĆ (10), Srđan MEDIĆ (11), Sanja KALAMBURA (12), Marko DUNĐER (13), Zlata DOLAČEK-ALDUK (4), Dina STOBER (4)
Editor-in-Chief:
Milan KLJAJIN
Technical Editor:
Goran KOZINA
Graphics Editor:
Snježana IVANČIĆ VALENKO, Anja ZORKO
IT support:
Tomislav HORVAT
Print:
Centar za digitalno nakladništvo, Sveučilište Sjever
The journal is free of charge and four issues per year are published.
Circulation: 100 copies
Registration of journal:
The journal "Tehnički glasnik" is listed in the HGK Register on the issuance and distribution of printed editions on the 18th October 2007 under number 825.
Preparation ended:
September 2019
Legend:
(1) University North, (2) Mechanical Engineering Faculty in Slavonski Brod, (3) Faculty of Graphic Arts Zagreb, (4) Faculty of Civil Engineering Osijek, (5) Faculty of Engineering Rijeka, (6) Faculty of
Mechanical Engineering and Naval Architecture Zagreb, (7) Faculty of Metallurgy Sisak, (8) Tomas Bata University in Zlín, (9) Department of Physics of the University of Josip Juraj Strossmayer in
Osijek, (10) Faculty of Humanities and Social Sciences Osijek, (11) Karlovac University of Applied Sciences, (12) University of Applied Sciences Velika Gorica, (13) Department of Polytechnics - Faculty
of Humanities and Social Sciences Rijeka, (14) Faculty of Electrical Engineering and Computer Science - University of Maribor, (15) Faculty of Civil Engineering - University of Maribor, (16) University
College of Teacher Education of Christian Churches Vienna/Krems, (17) Mechanical Engineering Faculty Tuzla, (18) Mechanical Engineering Faculty Sarajevo, (19) University of Travnik - Faculty of
Technical Studies, (20) Higher Education Technical School of Professional Studies in Novi Sad, (21) University of Novi Sad - Faculty of Technical Sciences, (22) Faculty of Mechanical Engineering -
University of Montenegro, (23) Brno University of Technology, (24) Odessa State Academy of Civil Engineering and Architecture, (25) Faculty of Civil Engineering - University of Mostar, (26) Faculty of
Manufacturing Technologies with the seat in Prešov - Technical University in Košice
CONTENT I
Erol CAN
POWER REGULATION BY COUPLE HALF WAVE LPWM RECTIFIER AT THREE-PHASE LOADS 184
Željko BARAČ, Ivan PLAŠČAK, Tomislav JURIĆ, Pavo BALIČEVIĆ, Vinko DUVNJAK, Mladen JURIŠIĆ, Goran HEFFER, Monika MARKOVIĆ
INFLUENCE OF AN UNEVEN SURFACE ON THE VIBRATION OCCURRENCE AFFECTING THE TRACTOR OPERATOR 192
Gordana JAUKOVIĆ
ADVANTAGES AND DISADVANTAGES OF RAMAN SPECTROSCOPY IN TESTING PAPER BANKNOTES 226
Iryna SOLONENKO
THE USE OF CEMENT CONCRETE PAVEMENTS FOR ROADS, DEPENDING ON CLIMATIC CONDITIONS 235
Abstract: In this paper, a 23 factorial design analysis was used to study the parameters affecting the mechanical characteristics of polypropylene/linear low-density
polyethylene/nano-titanium dioxide (PP/LLDPE/TiO2) nanocomposites, and to optimize these factors in order to predict the maximum ultimate tensile strength (UTS),
elastic modulus (EM), and yield strength (YS) simultaneously. To do this, two levels of nano-titanium dioxide (TiO2), linear low-density polyethylene (LLDPE) and styrene-
ethylene-butylene-styrene (SEBS) as the coupling agent were selected and eight experiments were conducted for every response. The most effective factors influencing
the UTS, EM, and YS were found, and acceptable prediction regression models were taken. One noted that nanoparticles increased the elastic modulus. The attendance
of high levels of LLDPE and SEBS resulted in a decrease in YS and UTS. Moreover, the optimum values of variables were determined by using the contour plot.
Here, a 23 factorial design was used to supply a relation −1 and +1 range in such a way that the researcher could
for the YS, UTS, and EM of PP/LLDPE/TiO2 choose an experimental design from a list of designs.
nanocomposites as a mathematical function of parameters The 2k factorial design is one of the most widely applied
(SEBS, TiO2, and LLDPE). Moreover, the effect of every designs to investigate the effects of various parameters on a
agent on mechanical characteristics was studied and the particular response, where k is the number of parameters and
optimal range of each parameter was found in order to the base 2 shows the level of treatment for each discussed
achieve the best tensile properties. parameter [22]. The performed design is shown in Tab. 1
which briefly discusses the parameters and the change of
2 THE EXPERIMENTS their levels. Furthermore, different modes of combining
2.1 Materials materials via the software were determined. As presented in
Tab. 2, eight tests had to be prepared for each response
Polypropylene (PP-Z30S, MFR-25, 230 °C, 2.16 kg) and having three replicates.
linear low-density polyethylene (LLDPE-0209, MFR-0.9,
190 °C, 2.16 kg, and density of 0.920 gr ml-1) were procured Table 1 Level of factors applied to study the effect of LLDPE, TiO2 and SEBS
from the Arak Petrochemical Company, Iran. The nano-TiO2 Factors LLDPE TiO2 SEBS
−1 (low) 40 0 0
rutile structure with the mean size of 30 nm and density of Level
+1(high) 60 2 3
4.23 g/cm was obtained from the Iranian Nanomaterials
Pioneers (INP), Iran. Moreover, the KRATON polymer type Table 2 Results of running the software
G, namely styrene-ethylene/butylene-styrene (SEBS), was Full Factorial Design
used as a coupling agent. The compound was prepared by Sample No. 1 2 3 4 5 6 7 8
using a co-rotating screw extruder (ZSK 25 P8.2E WLE) PP (Wt. %) 80 60 78 58 77 57 75 55
LLDPE (Wt. %) 20 40 20 40 20 40 20 40
with a 170 to 190 °C temperature range. Next, they were TiO2 (Wt. %) 0 0 2 2 0 0 2 2
made as granules. Granules were injected with the help of an SEBS (Wt. %) 0 0 0 0 3 3 3 3
injection molding machine (IMAN MACHINE 125g) with Random order of samples 3 2 8 1 6 4 7 5
the temperature profiles of 190-200-210 °C, and the samples
were hence prepared. As three parameters at two levels were assumed, the
experimental design was named a 23 full factorial design that
2.2 Mechanical Testing needed eight test runs for every response. The average of
results is presented in Tab. 3.
The elastic modulus, yield, and ultimate tensile strength
were specified by a Zwick/Roell–Z100 machine (Germany) Table 3 The experimental results for YS, UTS and EM
Sample No. 1 2 3 4 5 6 7 8
due to the ASTM D638 standard with the strain rate of 50
YS (MPa) 19.63 12.9 15.33 11.56 21.86 13.33 16.16 11.9
mm/min at room temperature. Fig. 1 shows a sample before UTS (MPa) 23.63 19.5 21.86 18.3 23.93 18.03 20.33 19
and after the tensile tests. EM (MPa) 223 189 226.6 207.3 212 184.6 222.6 199.3
Table 4 ANOVA results for YS Relying on ANOVA for YS, UTS, and EM, a fitted
Source F P regression model with statistical importance was found as
Main Effects 249.31 0.000
LLDPE 555.84 0.000
follows:
TiO2 175.14 0.000
SEBS 16.95 0.001 Yield Strength = 15313 − 2.888LLDPE − 1.621TiO2 +
2-Way Interactions 21.94 0.000 0.504SEBS + 0.929LLDPE × TiO2 − 0.312LLDPE × SEBS
LLDPE×TiO2 57.56 0.000 (1)
LLDPE×SEBS 6.51 0.021
TiO2×SEBS 1.76 0.203
R-sq: 98.07%, R-sq(Pred): 95.67%, R-sq(Adj): 97.23%
The effects of the LLDPE×TiO2 interaction on YS are Ultimate Tensile Strength = 20.575 − 1.867LLDPE −
greater than the importance effect of a single-factor (i.e. 0.7TiO2 − 0.25SEBS + 0.642LLDPE × TiO2 (2)
factor SEBS). The important interaction of LLDPE and TiO2 R-sq: 98.28% R-sq(Pred): 96.13%, R-sq(Adj): 97.53%
shows that these parameters are related, i.e. if the level of one
parameter varies, the effect of the other one varies, too. Elastic Modulus = 208.08 − 13LLDPE + 5.92TiO2 −
The ANOVA results of UTS are presented in Tab. 5. It 3.42SEBS + 2.33LLDPE × TiO2 (3)
is clear that LLDPE (P = 0.000), TiO2 (P = 0.000), SEBS (P R-sq: 97.86%, R-sq(Pred): 95.18%, R-sq(Adj): 96.92%
= 0.003), and LLDPE×TiO2 (P = 0.000) whose p-values were
below or equal to 0.05 for UTS should be statistically From Eq. (2), one sees that all factors had a negative
important. Moreover, from the F-values in Tab. 5, the main effect on UTS. Thus, a lower factor setting (−1) would
relative importance of each factor and its interactions would result in a higher response. In the situation of Eq. (1), an
be: increase in LLDPE×TiO2 and SEBS from small to high levels
resulted in 2.77% and 6.76% increases in YS, whereas an
LLDPE>TiO2> LLDPE×TiO2> SEBS. increase in LLDPE and TiO2 resulted in a decrease in YS by
31.7% and 19.14%. Thus, LLDPE had a maximum effect on
Table 5 ANOVA results for UTS YS with a 31.7% contribution. In the case of Eq. (2), by
Source F P adding of LLDPE, TiO2, and SEBS decreased UTS by
Main Effects 261.86 0.000 16.63%, 6.58%, and 2.40%, respectively. Therefore, LLDPE
LLDPE 678.05 0.000 had a maximum effect on UTS with a 16.63% contribution.
TiO2 95.35 0.000
In Eq. (3), an increase in TiO2 and LLDPE × TiO2 from low
SEBS 12.16 0.003
2-Way Interactions 27.04 0.000 to high levels resulted in 2.57% and 5.85% increases in EM,
LLDPE×TiO2 80.12 0.000 whereas an increase in LLDPE and SEBS decreased EM by
LLDPE×SEBS 0.66 0.428 11.76% and 3.2%, respectively. Hence, LLDPE had a
TiO2×SEBS 0.34 0.569 maximum effect on EM with a 11.76% contribution.
A verified model must predict the response with good
Table 6 ANOVA results for EM accuracy with respect to the experimental data. Model
Source F P
Main Effects 235.29 0.000
adequacy is checked by R-Sq, R-Sq (adj), and R-Sq (pred).
LLDPE 553.29 0.000 A R-Sq value near 100% means a reliable fit to the
TiO2 114.57 0.000 experimental data [22]. Based on the ANOVA results, the R-
SEBS 38.20 0.000 squared of the regression equations was 98.07 % for YS,
2-Way Interactions 6.25 0.005 98.28% for UTS, and 97.86% for EM, which means that the
LLDPE×TiO2 17.82 0.001
LLDPE×SEBS 0.36 0.555
model is verified. The adjusted R-square was 97.23% for YS,
TiO2×SEBS 0.57 0.462 97.53% for UTS, and 96.92% for EM, which accounts for the
amount of predictors in the model. The prediction R-squared
The EM ANOVA results are presented in Tab. 6. It is statistic was calculated to be 95.67% for YS, 96.13% for
observed that LLDPE (P = 0.000), TiO2 (P = 0.000), SEBS UTS, and 95.18% for EM. Because the predicted R-square
(P = 0.000), and LLDPE×TiO2 (P = 0.001), whose p-values values were near the R-square and the adjusted R-square
were below or equal to 0.05 for EM would be statistically values for every response, none of the models appeared to be
important. Moreover, from the F-values in Tab. 6, the overfitting and none had an adequate predictive ability [22].
relative importance of each factor and its interactions would
be: 3.2 Main Effects and Interaction Plot for Yield Strength
LLDPE>TiO2> SEBS> LLDPE×TiO2. The main effects plot in Fig. 2 indicates that YS
decreases as LLDPE and TiO2 contents increase. Therefore,
The multiple regression analysis was performed on the the maximum YS of PP/LLDPE/TiO2 nanocomposites could
experimentally collected data for the YS, UTS, and EM of be found at a lower LLDPE and TiO2. Moreover, Fig. 2
the PP/LLDPE/TiO2 nanocomposites. Here, the analysis was shows that by increasing the amount of the SEBS factor
done by the Minitab® 16 software which applies the ordinary increases the yield strength. The relative strength of the effect
least squares technique to find the regression function. of different parameters can also be seen. The main effects
plot for YS (Fig. 2) showed that LLDPE was the most 3.3 Main Effects and Interaction Plot for Ultimate Tensile
important factor. Strength
Fig. 3 shows the interaction plot between the three
discussed factors, namely LLDPE, TiO2, and SEBS for YS, The main effects plot in Fig. 4 indicates that increasing
respectively. The plots, known as interaction plots, are the amount of all three major factors, especially
employed to explain important interactions between process polyethylene, the ultimate tensile strength is reduced.
parameters. The interaction plot summarizes the interaction LLDPE decreased UTS for LLDPE, which was much
between the maximum and minimum amounts of each factor. smoother than PP. Therefore, the maximum UTS of
From these plots, one sees that the initial interaction PP/LLDPE/TiO2 nanocomposites would be reached at a less
happened between LLDPE and TiO2 for YS, demonstrated amount of LLDPE, TiO2, and SEBS. The relative strength of
by non-parallel lines. the effect of different parameters could also be seen. The
main effects plot for UTS (Fig. 4) showed that LLDPE was
the most important factor.
According to Fig. 3, at a low level of LLDPE (20 wt. %) Figure 5 Interaction plot for UTS
and TiO2 (0 wt. %), the interaction was very important, but
adding both parameters produced a low interaction. Fig. 5 demonstrates the interaction plot between the three
However, the important interaction of LLDPE versus TiO2 discussed factors, namely LLDPE, TiO2, and SEBS, for UTS,
for YS showed that a lower TiO2 (0 wt. %) would result in an respectively. From the plot, one sees that the initial
improvement in the YS of PP/LLDPE/TiO2 nanocomposites interaction happened between LLDPE and TiO2 for UTS,
when factor LLDPE was under the low level (20 wt. %), showed by non-parallel lines.
while the influence of TiO2 was reduced at the high level of Fig. 5 shows that, at the low levels of LLDPE (20 wt. %)
LLDPE (40 wt. %). These data for YS suggested that the and TiO2 (0 wt. %), the interaction was very important, but
ideal TiO2 of PP/LLDPE/TiO2 nanocomposites differs from by adding both parameters, it produced a low interaction.
the LLDPE. The models with low TiO2 show a large YS However, the important interaction of LLDPE in respect to
when LLDPE is low, and the models with great TiO2 indicate TiO2 for UTS showed that lower TiO2 (0 wt. %) would result
a low YS when LLDPE is high. in an increase in the UTS of PP/LLDPE/TiO2
nanocomposites when LLDPE was under the low level (20
wt. %), while the effect of TiO2 was reduced at the high level LLDPE was the most important factor. The elastic modulus
of LLDPE (40 wt. %). These data for UTS suggested that the of polymer nanocomposites largely depends on the good
ideal TiO2 of PP/LLDPE/TiO2 nanocomposites differs with dispersion of nanoparticles in the matrix. Fig. 7 shows field
LLDPE. The models with low TiO2 show large UTS when emission scanning electron microscopy (FESEM) images
LLDPE is low, and the models with high TiO2 demonstrate taken from the samples’ fractured surface. It can be observed
small UTS when LLDPE is high. from Fig. 7 that nanoparticles are well dispersed in the
matrix. As a result, titanium oxide nanoparticles lead to an
3.4 Main Effects and Interaction Plot for Elastic Modulus increase in the elastic modulus.
Fig. 8 illustrates the interaction plot between the three
The main effects plot in Fig. 6 indicates that the elastic investigated parameters, namely LLDPE, TiO2, and SEBS,
modulus decreased as LLDPE and SEBS varied from a low for EM, respectively. One could see that the initial interaction
to a high level. Consequently, the maximum elastic modulus happened between LLDPE and TiO2 for EM, indicated by
of PP/LLDPE/TiO2 nanocomposites could be obtained at non-parallel lines (Fig. 8). Fig. 8 shows that, at the low levels
lower LLDPE and SEBS. of LLDPE (20 wt. %) and TiO2 (0 wt. %), the interaction was
highly significant, and that increasing both parameters
produced a good interaction and increased the EM. However,
the important interaction of LLDPE versus TiO2 for EM
showed that a high TiO2 (2 wt. %) would result in an increase
in the EM of PP/LLDPE/TiO2 nanocomposites when LLDPE
was under the high level (40 wt. %), while the effect of TiO2
was reduced at the low level of LLDPE (20 wt. %). These
data for EM suggested that the ideal TiO2 of PP/LLDPE/TiO2
nanocomposites varied with LLDPE. The models with high
TiO2 show large EM when LLDPE is high, and the models
with low TiO2 indicate small EM when LLDPE is low.
combination, the ultimate tensile strength of more than 32 been fixed in the middle level. It is clear that, to achieve the
MPa can be achieved. Based on Fig. 9, the presence of high best elastic modulus, low amounts of LLDPE and SEBS and
values of LLDPE (more than 35 wt.%) led to a significant high amounts of titanium dioxide nanoparticles should be
reduction in the ultimate tensile strength. used. By using this combination, the elastic modulus of more
Fig. 10 depicts the contour plot of the elastic modulus as than 220 MPa can be achieved. Based on Fig. 9, the presence
a function of TiO2*LLDPE, SEBS*LLDPE, and of high values of TiO2 led to a significant increase in the
TiO2*SEBS. In any of these three modes, the third factor has elastic modulus.
6.76% increase in YS, whereas an increase in LLDPE led [11] Daneshpayeh, S., Ashenai Ghasemi, F., Ghasemi, I., & Ayaz,
to a decrease in YS by 31.7%. An increase in LLDPE M. (2015). Predicting of mechanical properties of
and SEBS decreased the UTS to 16.63% and 2.40%, PP/LLDPE/TiO2 nanocomposites by response surface
methodology. Compos. Part. B, 84, 109-120.
respectively. Moreover, an increase in TiO2 and
https://doi.org/10.1016/j.compositesb.2015.08.075
LLDPE×TiO2 (from a low to a high level) resulted in [12] Garcia, M., Vilet, G. V., & Jain, S. (2004). Polypropylene/ SiO2
2.57% and 5.85% increases in EM, respectively. nanocomposites with improved mechanical properties. Rev
- The optimized ranges of variables on the tensile Adv Mater Science, 6(2), 169-175.
properties were found by using the contour plot. The [13] Selvin, T. P. & Kuruvilla, J. (2004) Mechanical properties of
results show that the most improved tensile properties titanium dioxide filled Polystyrene microcomposites. Mater
were obtained in the low level of LLDPE and SEBS and let, 58(3-4), 281-289.
the middle level of titanium dioxide nanoparticles. https://doi.org/10.1016/S0167-577X(03)00470-1
[14] Wacharawichananat, S., Thongyai, S., & Tipsri, T. (2009).
Effect of mixing conditions and particle sizes of titanium
5 REFERENCES
dioxide on mechanical and morphological properties of
Polypropylene/Titanium dioxide composites. Iran Polymer
[1] Zapata, P. A., Rabagliati, F. M., Lieberwirth, I., Catalina, F., &
journal, 18(8), 607-616.
Corrales, T. (2014). Study of the photodegradation of
[15] Altan, N. & Yildirim, M. (2010). Mechanical and
nanocomposites containing TiO2 nanoparticles dispersed in
Morgholigical properties of Polypropylene and High density
polyethylene and in poly (ethylene-cooctadecene). Polym
polyethylene matrix composites reinforced with surface
Degrad Stabil, 109, 106-114.
modified nano sized TiO2 particles. World Acad Sci Eng
https://doi.org/10.1016/j.polymdegradstab.2014.06.020 Technol, 4(10), 252-257.
[2] Chen, J. H., Zhong, J. C., Cai, Y. H., Su, W. B., & Yang, Y. B.
[16] Premphet, K. & Horanont, P. (2000). Phase structure of ternary
(2007). Morphology and thermal properties in the binaryblends
polypropylene/elastomer/filler composites effect of elastomer
of poly (propylene-co-ethylene) copolymer andisotactic
polarity. Polym, 41(26), 9283-9290.
polypropylene with polyethylene. Polymer, 48(10), 2946-
https://doi.org/10.1016/S0032-3861(00)00303-7
2957. https://doi.org/10.1016/j.polymer.2007.03.037
[17] Abughalia, M., Hassan, A., & Yussuf, A. (2011) Mechanical
[3] Owpradit, W. & Jongsomjit, B. (2008). A comparative study
and Thermal properties of Calcium Carbonate filled
on synthesis of LLDPE/TiO2 nanocomposites using different
PP/LLDPE Composite. J Appl Polym Sci, 121(4), 2413-2421.
TiO2 by in situ polymerization with zirconocene/dMMAO
https://doi.org/10.1002/app.33570
catalyst. Mater Chem Phys, 112(3), 954-961.
[18] Gang, L., Feng, T. Y., Yuan, Y. F., Cheng, Z. L., Xing, Z. Z.,
https://doi.org/10.1016/j.matchemphys.2008.07.050 & Ji, X. Q. (2005). Effect of nanoscale SiO2 and TiO2 as the
[4] Chaichana, E., Jongsomjit, B., & Praserthdam, P. (2007).
fillers on the mechanical properties and aging behavior of
Effect of nano-SiO2 particle size on the formation of
Linear Low Density Polyethylene/Low Density Polyethylene
LLDPE/SiO2 nanocomposite synthesized via the in situ
blends. J Polym Env, 13(4), 339-348.
polymerization with metallocene catalyst. Chem Eng Sci,
https://doi.org/10.1007/s10924-005-5528-x
62(3), 899-905. https://doi.org/10.1016/j.ces.2006.10.005
[19] Altan, M. & Yildirim, H. (2012). Mechanical and Antibacterial
[5] Zapata, P. A., Palza, H., Cruz, L. S., Lieberwirth, I., Catalina,
Properties of Injection Molded Polypropylene/ TiO2 Nano-
F., Corrales, T., & Rabagliati, M. (2013). Polyethylene and
Composites: Effects of Surface Modification. J Mater Sci
poly (ethylene-co-1-octadecene) composites with TiO2 based
Technol, 28(8), 686-692.
nanoparticles by metallocenic in situ polymerization. Polymer,
https://doi.org/10.1016/S1005-0302(12)60116-9
54(11), 2690-2698. [20] Guo, H. & Li, X. (2012). Preparation and Properties of
https://doi.org/10.1016/j.polymer.2013.03.048 Transparent SEBS/Titania Nanocomposite Films via
[6] Palza, H., Reznik, B., Kappes, M., Hennrich, F., Naue, I. F. C., Functionalization of SEBS and Sol-Gel Process. Polym
& Wilhelm, M. (2010). Characterization of melt flow Compos, 20(1/2), 155-160.
instabilities in polyethylene/carbon nanotube composites.
https://doi.org/10.1177/0967391112020001-230
Polymer, 51(16), 3753-3761. [21] Xue, B., Li, F., Xing, Y., Sun, M., Liu, D., & Jiang, Y. (2011).
https://doi.org/10.1016/j.polymer.2010.06.016 Preparation of Cu/OMMT/LLDPE nanocomposites and
[7] Palza, H., Vergara, R., & Zapata, P. A. (2011). Composites of synergistic effect study of two different nano materials in
polypropylene melt blended with synthesized silica polymer matrix. Polymer Bulletin, 67(8), 1463-1481.
nanoparticles. Compos Sci Technology, 71(4), 535-540.
https://doi.org/10.1007/s00289-011-0466-3
https://doi.org/10.1016/j.compscitech.2011.01.002 [22] Montgomery, D. C. (2011). Design and Analysis of
[8] Tong, Y., Li, Y., Xie, F., & Ding, M. (2000). Preparation and
Experiments. John Wiley & Sons, New York.
characteristics of polyimide–TiO2 nanocomposite film. Polym
[23] Box, G. P. G., Hunter, J. S., & Hunter, W. G., (2005). Statistics
Int, 49(11), 1543-1547.
for Experimenters: Design, Innovation and Discovery. John
https://doi.org/10.1002/1097-0126(200011)49:11<1543::AID- Wiley & Sons, USA.
PI535>3.0.CO;2-B
[9] Giang, V., Thai, H., Huynh, D., Trung, H., Lam, D. T., & Tuan,
M. V. (2013). Effect of titanium dioxide on the properties of
polyethylene/TiO2 nanocomposites. Compos B-Eng, 45(1),
1192-1198. https://doi.org/10.1016/j.compositesb.2012.09.058
[10] Norio, N. & Toyoharu, H. (2008). Preparation of TiO2
nanoparticles surface-modified by both carboxylic acid and
amine: dispersibility and stabilization in organic solvents.
Colloids Surf A, 317, 543-550.
https://doi.org/10.1016/j.colsurfa.2007.11.036
Authors’ contacts:
Sajjad DANEAHPAYEH
Faculty of Mechanical Engineering,
Shahid Rajaee Teacher Training University,
Shabanloo St., Lavizan, Tehran, Iran
Tel: +9821-22970060
sajad_danesh66@yahoo.com
Ismail GHASEMI
Department of Polymer Processing,
Iran Polymer and Petrochemical Institute (IPPI),
15 km Tehran-Karaj Highway, Pajuhesh Science and Technology Park,
Pajuhesh Boulevard, Tehran 14977-13115, Iran
Tel: +9821-44787050/
I.Ghasemi@ippi.ac.ir
Abstract: Ceramic materials are widely used in impact safekeeping systems. Ceramic is a heterogeneous material; its characteristics depend considerably both on specifications
of its ingredients and the material structure completely. The finite element method (FEM) can be a useful tool for strength computation of these materials. In this paper, the
mechanical properties of the ceramic composites are investigated, and the mechanical performance modeling of fiber-fortified ceramic matrix composites (CMC) is expressed by
the instance of aluminum oxide fibers in a matrix composite based on alumina. The starting point of the modeling is an infrastructure (primary cell) that contains a micromechanical
size, the statistical analysis characteristics of the matrix, fiber-matrix interface, fiber, and their reciprocal influences. The numeral assessment of the model is done using the FEM.
The numerical results of composite elastic modulus were computed based on the amount of the added fibers and the porosity was evaluated for empirical data of samples with a
similar composition. Various scanning electron microscope (SEM) images were used for each sample to specify the porosity. Also, the unit cell method presumed that the porous
ceramic substance is manufactured from an array of fundamental units, each with the same composition, material characteristic, and cell geometry. The results showed that when
the material consists of different pores and fibers, the amount of Young’s modulus reduces with the increment of porosity. The linear correlation model of elasticity versus porosity
value from experimental data was derived by MATLAB curve fitting. The experimental data from the mechanical test and numerical values were in good agreement.
Keywords: ceramic composite; elasticity; fiber; Finite Element Method (FEM); mechanical performance; thermal stress
and realized in making a specific material. Many modern possible: (1) a parallel arrangement of the fiber longitudinal
technologies need materials, such as metals, ceramics, metal axis in a single direction; and (2) an accidental adjustment.
alloys, and polymeric materials with unique combinations of Mechanical reactions of this unique composite depend on
properties that cannot be carried out. The required composite multiple parameters, such as strain-stress behaviors of the
materials for underwater, aerospace, and transportation matrix phases, fiber, the phase volume fractions, and the
applications have unique properties. Most of the composites direction a load or the stress is enforced to it [10]. In recent
typically have one construction in common, a binder or years, ceramic material has become increasingly substantial,
matrix mixed with a fortified material. A composite is a especially for unique uses that need great strength. Ceramics
material developed from two or more distinct phases, matrix show superior behavior at high temperatures in general. The
and dispersed. If the fibers are directionally oriented and difficulty of using ceramics is their fragile harm behavior. A
continuous, the material is called an advanced composite. single deficiency or defect can cause fragile damage in the
The composite material properties depend on the stiffness of whole structure [11].
bonds and the physical-mechanical features of their Fiber-reinforced or -modified composites have great
components. A specific advantage of composites is the mechanical performance, such as high unique strength and
capability of their components to be utilized fully. Such hardness. Particularly, laminated composite structures have
materials may take specified valuable properties that are not widely been used, where the in-plane properties are
found in the components. To obtain the optimal characteristic substantial [12].
in composites, their components are selected to have strongly In past decades, the finite element method as an efficient
various but perfect properties. The basis or matrix of tool has become a broadly accepted computational technique
composites may consist of metals or alloys (metallic not only in computational mechanics but also in certain
composites), carbon, polymers, and ceramics (non-metallic engineering majors. The approach has been extended to solve
composites) [8]. Technologically, the most substantial easy linear equations in mechanics of solid bodies along with
composites are generally those in which the dispersed phase severe non-linear problems in almost all fields of applied
is in the form of a fiber. Indeed, the powerful fibers physics or mathematics; and more recently, in many other
embedded in a softer matrix produce products with great branches of general science.
strength-to-weight proportions. The matrix materials Some packages that can be used are FLUENT,
transmit the load to the fibers that absorb the stress. ABAQUS, ADINA, ANSYS, NASTRAN, and LUSAS.
Also, the linear behavior of ordinary composites and its Finite element method gives a proximate distribution of the
frame and structure can be described analytically by Classical needed size in the considered field that is not possible to be
Laminate Theory (CLT). By investigating elastic mechanical specified analytically. A proximate solution is derived by
characteristics (Young’s moduli parameter, shear modulus of dividing the total area (bodies) to the minor elements. Then,
rigidity, and Poisson’s ratios), as well as influence of the related physical rules are applied to each element and
moisture and thermal impact, CLT can effectively compute total elements of the set of solutions to solve the existing
the stress-strain relations in composite become laminated up problems [13]. One of the novel numerical simulation
to fracture initiation. Through thickness, both the vertical and methods is finite element analysis (FEA) that can evaluate
shear forces are ignored in CLT, thus supposing composite the fracture probability caused by the characteristics and
plates to be extremely wide and long that causes neglecting performance of flaw distribution. It is considered an effective
edge effects. Finite element modeling (FEM) dominates tool to facilitate and increase the use of ceramics in
these issues and can perform the stress investigation of components and members; the FEA can also forecast the
structures more carefully. Computational or numerical distribution of ceramic rigidity. Specifically, the needed data
simulation results are being generally used in the design, on the microstructure distribution (i.e., relative density,
production, and investigation of several composites and porosity, and grain size) are used as the input parameters
structures. The important parameter in increment of [14]. Also, the finite element method is done to model the
industrial competitiveness is saving or reducing the design piezo-laminated plate panel as an appropriate composite
age. The capability to finish designs by a proper software on structure and show the static and dynamic analysis of
the computer effectively reduces the reliance on long-time structural vibrations subjected to an external force. The
methods and costly physical experiments. Composites can be piezo-laminated structure modeling is based on the
nearly optimized for suitable and special uses. Most assumptions of first-order shear deformation theory (FSDT)
commercial FE software possess user-friendly composite [15]. In general, the purpose of the numerical model is to
mathematical modeling capabilities, which can accurately investigate the effect of substantial parameters on the
obtain elastic action through anisotropic or orthotropic behavior of the fiber composite. This grid frame or
constructional modeling. Unlike CLT, FE models allow substructure (elementary cell) takes fiber, the fiber-matrix
accurate and regular stress analyses through the interface, the single components matrix, and their particular
accomplishment of an appropriate computerized sample harm behavior into an account. It is selected in such a way
model, leading to the anticipation of stress concentrations [9]. that a macrostructure property for the complete structure can
The arrangement or orientation of the fibers, the fiber be made by a proper number of substructures. The model has
density, and their distribution all have a significant influence been performed in a finite element method code to research
on the resistance and other performances of fiber-modified the effect of separate component factors on the behavior of
composites. With regards to direction, two extremes are the total.
In a finite element method modeling, several analyses and elasticity. This study focuses on the mechanical
can be done; for instance, dynamic and static analysis, which performance of ceramic composites as well as the
later could be modeled and involved forced vibrations. In comparison of experimental data and simulation results.
general, a static linear elastic analysis is the easiest and
default one but with limited use due to its closed scope [16]. 1.1 Methodology
Modeling of the accidental fiber distribution of a fiber-
resistant composite via simulation software is of high The present study deals with technical performances of
importance to investigate the approaching failure behavior of the elastic theory originated from FEA of the total
the material on the limitation of micro-scale. Fiber-reinforced representative volume elements in fiber-modified
composites are typically recognized as ordinal materials with composites. The fibers are put orderly in the square array that
three structural levels: macro-scale, mesoscale, and micro- is known as the uni-directional composite. An example of a
scale. The micro-scale defines the adjustment of composite uni-directional fiber composite is shown in Fig. 1.
fibers in the fiber bundle. The mesoscale typically relates to
the fabric/lamina shape or geometry, and the macro-scale
composite refers to engineering the structural action of the
stuff. In a multi-scale modeling and simulation of composite
stuff, micro-scale concepts (including both numerical-
analysis and mathematical methods) are usually applied to
anticipate the effective rigidity and strength performance of
axial symmetry of a composite material or transversely
isotropic fundamental properties of composites, which are
used as theoretical tools for technical structure design [17,
18]. Recently, a finite element analysis (FEA) modeling has
been performed for the long fiber- modified ceramic matrix
woven composites (LFMCWCs) sample. For instance, this
method has been validated to investigate the grinding process Figure 1 The schematic image of the unit cells.
of a 2.5D woven quartz such as fiber-reinforced silicon
dioxide ceramic matrix (SiO2/SiO2) composite. Concerning It is supposed that the fiber and matrix material is usually
the grinding tests, their results were directly confirmed for linearly elastic. Typically, a unit cell is employed for the
the precision of the FEA method [19]. analysis. The measure of the fiber volume in proportion to
The strength of these studies is their uniqueness, using the whole volume of the composite is served from the cross-
various materials and composites, as well as FEM use as a sectional areas of the fiber associated with the total cross-
useful tool to analyze and predict the mechanical sectional area of each unit cell. This fraction of composite is
performance of composites, which is a cost-effective and considered as a substantial parameter in a composite that is
time-saving approach. However, the past research named the fiber volume fraction (VF).
weaknesses are testing only a few limited materials and some
simulations without comparing with experimental data, in 2 EXPERIMENTAL METHOD
particular, microstructure SEM images that are very vital for
investigating mechanical performance. Some samples were made using the layered claystone,
In the present research, only one composite has been clay, and bauxite in the amounts defined by the Jahad
modeled with a total detailed account of its physical Daneshgahi Corporation according to scientific principles
morphology along with its thermo-physical performances. [21, 24]. Samples with the addition of 1% wt. content of low
The composite morphology at micro-scale has been made as alumina-based fibers (with the aspect ratio of l/d ≈ 17), were
close to reality as possible with measurements from SEM added to the composition in order to evaluate the results of
micrographs for the geometric details, such as the fiber radius numerical modeling to the empirical data. The material was
and length. A commercial FE code of ABAQUS has been blended in the ball mill with water to enable better
applied to create the computational model. The real distribution of components, proper fibers in particular. The
microstructure with fiber distribution has been specified via obtained composition was dried within 24 hours at 100 °C
an SEM image, for example. The science of image analysis and then pressed between pressures of 35000 kPa to 55000
has then been done on the main image to determine the fibers kPa to create two various porosities. Other samples were
based on various color threshold modules, histograms and thermally operated at about 1200 °C and then were tested via
algorithms. From this image, the data such as the distribution the compaction and the Brazilian tests [20] to get the values
of distances between neighboring fibers and the fiber radius of elasticity modulus (E) and strength (F). The test was done
extracted based on image pixels in a black and white image via the modern-hydraulic testing system (Instron testing
via MATLAB image processing toolbox. The obtained system, Model 1362) with a load cell of 5000 N and 100 kN
statistical parameters were utilized and employed mesh grids through a data logger. The calculated values of modulus were
to generate a finite element analysis that was used to served to evaluate the numerically obtained values [21]. The
anticipate the effective mechanical property and compare it whole samples were experimented using the scanning
with the experimental data, such as volume fraction, porosity,
electron microscope (SEM), and captured images were of the fibers was chosen while making the material and fixed
employed to determine their porosity [22]. on 1% volume. The fiber dimensions were obtained by
Porosity values were measured based on the various evaluating short fibers via the image processing program.
threshold level of darkness of the obtained images via the According to Fig. 1, the average fiber length was 0.1344 mm,
image processing method. Different SEM images of each and the average fiber radius was 0.003 mm. Thus, the
sample were used to determine the porosity of various formula for calculating the volume of the single cylindrical
samples. This method is appropriate since the pore, absorbs fiber was computed as follows:
more light, and it takes place as part of the image darker than
of the solid material. Given the porosity of the material and π r 2l =
Vf = π (0.003) 2 × 0.1344 =×
3.8 10−6 mm3 (1)
the proper method for determining the number and share of
pores, given in Fig. 2, the microstructure images were
The material volume was computed from total fiber
captured using an SEM system, (the Zeiss DSM-960A
constituents using the simple equation based on the cell
model, operated at the range of 21 kV).
volume of 0.0000038 mm3 and edge length of 0.068 mm.
Since the fiber length is 0.1345 mm and also the unit cell can
be big enough, the fiber can fit in each cell. The presumption
that eight fibers belong to one elementary cell provides a
cube with the dimensions of 0.14 mm × 0.14 mm, and the
depth of 0.14 mm; also, the fiber radius is r = 0.003425 mm,
and the length is l = 0.137388 mm. They are the same for
entire fiber volume fractions. The calculated volume
fractions of elementary cells for various fibers have an error
range of 0.04-0.1 in proportion to various fiber volumes in
ceramic.
All fibers were randomly put into a matrix, which
included 4 and 8 fibers. The schematic representation of fiber
has been shown in Fig. 3.
Fig. 5 shows the unit cell: a) the fibers with pores (porous been reported in Tab. 1 according to the density (without
fibers); and b) the fibers without pores. It can be realized that pores) from the Tab. 1.
the existing contact between matrix and fibers is the surface
to surface contact type with a friction factor of 0.1. In the Table 1 Physical properties of alumina (Al2O3).
software, it is determined that the surface to surface contact 𝜌𝜌 Cp Λ E ν Α
(g/cm3) (J/g·K) (W/m·K) (GPa) (1/K*106)
is as of the finite sliding and the friction formulation as of the 3.989 0.78 29 390 0.22 9.1
penalty.
Fig. 6 shows the matrix with 8 fibers. This model of
3.1 Calculation of Mechanical Properties by Analytical elementary mesh grid cell is originated from the matrix in
Solution which it does not have any pores.
The mechanical performance of the lamina is computed 3.2 Comparison of Experimental and Calculated Parameters
using the following equations of Elasticity Theory, as well as of Elastic Modulus
Halphin-Tsai’s formulae; Linear Young’s Modulus equation
in the fiber for each direction is as follows: To evaluate the achieved numerical simulation results,
the derived data were tested and analyzed via a mechanical
σ experiment. The obtained values of elastic modulus for
E= (2)
ε various combinations are divided into two main parts: one
with only pores, and one with fibers and pores
σ is the stress, and ε is the strain in each direction (x, y or z). simultaneously, which are shown in Fig. 7. The experimental
data were also added to the data diagram and fitted with the
model extended in this study. When the material consists of
ε2
ν 12 = − (3) pores and fibers, the quantity of Young’s modulus reduces
ε1 with the increment of porosity. The value of Young’s
modulus in the improved material with fiber has greater
ε1, ε2 are the strains in various directions (x, y, z). amount as compared to the material without fiber including
the same porosity. Poisson’s ratio, generally, decreased with
an increasing fiber volume fraction according to Eq. (4).
ν vf Vf + vmVm
= (4)
Figure 5 Elementary cells: a) with fibers and pores; and b) with fibers without
pores.
Figure 6 The impact of the load on the elementary cell. Fig. 8 shows the linear regression model of the elastic
modulus obtained by the experimental data.
Experimental data for pure Al2O3 used to simulate the This linear correlation was calculated by MATLAB
alumina ceramic are summarized in Tab. 1. It should be noted software of version R2014b as follows:
that the physical properties of the material (Density ρ,
Specific heat Cp, Thermal conductivity λ, Young’s modulus ε+
−0.08871 10.27
E= (5)
E, Poisson’s ratio ν, and thermal expansion factor α) have
4 CONCLUSION
Acknowledgments
Figure 9 Schematic flowchart describing the experimental procedure of samples. This research has been guided and supported by a
professional mechanical engineer, Seyed Hassan Alhoseini.
In the present research, the finite element analysis has
been done to anticipate the engineering and mechanical 5 REFERENCES
properties of the fibers in the field of a fiber-modified
ceramic composite. The simulation results were validated [1] James Paul, K. & Sd. Abdul Kalam. (2016). Mechanical
Properties of Okra Fiber Reinforced Composites Using Fem,
International Journal of Advanced Research in Mechanical [17] Wenzhi Wang, Yonghui Dai, Chao Zhang, Xiaosheng Gao &
Engineering & Technology (IJARMET), 2(2). Meiying Zhao. (2016). Micromechanical Modeling of Fiber-
[2] Shchurova, E. I. (2016). Modeling of the Ceramics Structure Reinforced Composites with Statistically Equivalent Random
for the Finite Element Analysis. Procedia Engineering, 150, Fiber Distribution. Journal of Materials (MDPI).
179-184. https://doi.org/10.1016/j.proeng.2016.06.744 https://doi.org/10.3390/ma9080624
[3] Sadik, C., El Amrani, I. E., & Albizane, A. (2014). Recent [18] Samadi, F. & Missagh Shamshiri, G. S. (2019). Numerical
advances in silica-alumina refractory: A review. Journal of Modeling for Simulation of Contaminant Migration of
Asian Ceramic Societies, 2(2), 83-96. Leachate in Soil Media. arXiv preprint arXiv:1905.02280.
https://doi.org/10.1016/j.jascer.2014.07.006 [19] Yuguo Wang, Haoji Wang, Jinhua Wei, Bin Lin, Jingyu Xu, &
[4] Choi, A. H., Heness, G., & Ben-Nissan, B. (2014). Using finite Sheng Fang. (2019). Finite element analysis of grinding
element analysis (FEA) to understand the mechanical process of long fiber reinforced ceramic matrix woven
properties of ceramic matrix composites. In Advances in composites: Modeling, experimental verification and material
Ceramic Matrix Composites, 286-311. removal mechanism. Ceramics International, 45(13), 5920-
https://doi.org/10.1533/9780857098825.2.286 15927. https://doi.org/10.1016/j.ceramint.2019.05.100
[5] Choi, A. H., Heness, G., & Ben-Nissan, B. (2018). Using finite [20] Andreev, G. E. (1991). A review of the Brazilian test for rock
element analysis (FEA) to understand the mechanical tensile strength determination. Part I: calculation formula.
properties of ceramic matrix composites. In Advances in Mining Science and Technology, 13(3), 445-456.
Ceramic Matrix Composites (pp. 375-400). https://doi.org/10.1016/0167-9031(91)91006-4
https://doi.org/10.1016/B978-0-08-102166-8.00016-5 [21] Dimitrijević, M. M., Veljović, D., Posarac-Marković, M.,
[6] Boccaccini, A. R. & Fan, Z. (1997). A new approach for the Jančić-Heinemann, R., Volkov-Husović, T., & Zrilić, M.
Young's modulus-porosity correlation of ceramic materials. (2012). Mechanical properties correlation to processing
Ceramics International, 23(3), 239-245. parameters for advanced alumina-based refractories. Science of
https://doi.org/10.1016/S0272-8842(96)00033-8 Sintering, 44(1), 25-33. https://doi.org/10.2298/SOS1201025D
[7] Farooqi, J. K. & Sheikh, M. A. (2006). Finite element [22] See http://thermalceramics.com/products/categories/bulk
modelling of thermal transport in ceramic matrix composites, [23] Yang, Z. (2010). Effect of Interphase on the Prediction of
Computational Materials Science 37, 361-373. Mechanical Properties for Unidirectional Composites. Ph.D.
https://doi.org/10.1016/j.commatsci.2005.11.001 Thesis, Harbin Institute of Technology, Harbin, China.
[8] Callister, W. D. & Rethwisch, D. G. (2011). Materials science [24] Bunshah, R. F. (2001). Handbook of hard coatings, Deposition
and engineering (Vol. 5, pp. 344-348). NY: John Wiley & Technologies, Properties. New York: Noyes Publication.
Sons. [25] Dimitrijevic, M., Heinemann, R. J., Husovic, T. V., Posarac,
[9] Reiner, J. & Vaziri, R. (2018). Structural Analysis of M., & Majstorovic, J. (2011). Morphological analysis of
Composites with Finite Element Codes: An Overview of surface degradation of advanced alumina based refactories
Commonly Used Computational Methods, Comprehensive subjected to thermal shock. Procedia Engineering, 10, 2153-
Composite Materials II, Vol. 8. 2157. https://doi.org/10.1016/j.proeng.2011.04.356
https://doi.org/10.1016/B978-0-12-803581-8.10050-5
[10] Ismar, H., Schröter, F., & Streicher, F. (2000). Modeling and
numerical simulation of the mechanical behavior of woven Authors’ contacts:
SiC/SiC regarding a three-dimensional unit cell.
Computational Materials Science, 19(1-4), 320-328. Majid AHMADI
https://doi.org/10.1016/S0927-0256(00)00170-1 Department of Mechanical Engineering,
[11] Kakani, S. L. & Kakani, A. (2010). Material Science, New Age Science and Research Branch,
Islamic Azad University (SRBIAU),
International Pvt Ltd Publishers.
Daneshgah Blvd, Simon Bulivar Blvd, Hesarak area, Tehran, Iran
[12] Tao Zeng, Lin-zhi Wu, & Li-cheng Guo. (2004). Mechanical
analysis of 3D braided composites: a finite element model. Seyed Hadi SEYEDIN
Composite Structures, 64, 399-404. Corresponding author
https://doi.org/10.1016/j.compstruct.2003.09.041 Department of Chemical Engineering,
[13] Dimitrijević, M. M., Tomić, N., Međo, B., Jančić-Heinemann, Science and Research Branch,
R., Rakin, M., & Volkov-Husović, T. (2014). Modeling of the Islamic Azad University (SRBIAU),
mechanical behavior of fiber-reinforced ceramic composites Daneshgah Blvd, Simon Bulivar Blvd, Hesarak area, Tehran, Iran
sh.seyedin@iran.ir
using (FEM). Science of Sintering, 46(3), 385-390.
https://doi.org/10.2298/SOS1403385D Seyed Vahid SEYEDIN
[14] Ozaki, S., Aoki, Y., Osada, T., Takeo, K., & Nakao, W. (2018). Department of Software Engineering,
Finite element analysis of fracture statistics of ceramics: Tehran South Branch, Islamic Azad University,
Effects of grain size and pore size distributions. Journal of the Deh haghi (Ahang) Ave, Abuzar Blvd, Afsariyeh Highway, Tehran, Iran
American Ceramic Society, 101(7), 3191-3204.
https://doi.org/10.1111/jace.15468
[15] Sharma, A. (2018). The finite element based study on active
attenuation of structural vibration using porous ferroelectric
ceramics. In AIP Conference Proceedings (Vol. 1975, No. 1, p.
030007). AIP Publishing. https://doi.org/10.1063/1.5042177
[16] Pedro Colmar Gonçalves da Silva Vellasco, Luciano Rodrigues
Ornelas de Lima, Luís Alberto Proença Simões da Silva.
(2017). Finite Element Modelling, Modelling Steel and
Composite Structures, Chapter 3, 209-211.
https://doi.org/10.1016/B978-0-12-813526-6.00003-9
Abstract: During the last two decades, technological development has led to an extreme increase in transport and connected industries. This has significantly increased the
production of automobile tires, which have their life span after which they go to waste. It is therefore necessary to extricate new products from recycled raw materials. For a product
to be created, it is necessary, especially in civil engineering, to examine and determine all the properties of this material (pressure, traction, stress resistance, etc.). The results
obtained by the research in this paper are mathematically processed by applying a logarithmic model. The aim of this process is to predict stress deformation in terms of tensile
force. The model obtained is significant with accuracy of 87.21% and has a very high accuracy of the deformation estimation in relation to the applied stress. Independent variables
were granulometric composition, binder (glue) and specific mass.
The results were processed using a logarithmic model. The amount of binder (polyurethane adhesive) ranges
Such a model is applied when we have multiple independent from 29.6 to 67 kg/m3 of the sample. A single manufacturer
variables each of which is subject to logarithm. Then, control may use a different type of binder within the material
of the model is done in terms of its correctness, significance, structure. The tested samples are based on polyurethane-
etc. Such actions are needed to prove its realistic application. based glue, while other binders such as epoxy can be used,
The proposed model predicts with certain accuracy the etc., which is especially important in civil engineering,
deformation of the real sample in relation to the assumed because of the economical acceptance, e.g. production of
stress provided that the known values for the specific mass construction panels for improvement of soundproofing
and binder amount are known. The minimum accuracy must properties on light metal structures and construction panels
be 80% to make the model applicable and acceptable. In such for making roadside sound barriers.
conditions we get a mathematical expression so that the
mentioned sizes can be calculated. 2.3 Results of Samples Testing
2.2 Main Characteristics of the Test Samples Twenty-seven samples were tested with different
properties. The changing variables are the granulometric
The samples were manufactured according to the valid composition, the specific weight, and the amount of the
norm [9] and as such were tested on the axial device. Load binder. In Fig. 3 we can see a graphic representation of the
appliance is manual, and monitored and read through stress ratio and the associated deformation of a particular
measuring devices. The following figure shows the test sample. In the same figure it can be noticed that there are
specimen according to HRN EN ISO 527-4: 2008 (for testing significant differences in values and that the influence of
of the tensile strength of polymeric materials) [9]. The individual variables in the sample on the results of the test is
detailed dimensions of the sample are shown in Fig. 2. visible.
Sample thickness
10 mm
15 mm
20 mm
with an accuracy of 87.21% we can predict the deformation such as specific mass and amount of binder, must remain
results in relation to the predicted stress if we apply this invariable, otherwise the model is not valid and cannot be
model. It should be noted that other independent variables, applied.
3.1 Testing Models of Expectation or Significance fulfilled, and we conclude that the model is correct and
reliable.
Testing the significance of the model is carried out to
determine the level of significance of the test. The test Table 5 The results of logarithm test results base-heteroscedasticity
Spearman Rank Order Correlations (Spreadsheet 14) MD pairwise
provides the information with which reliability is performed.
deleted Marked correlations are significant at p < .05000
If the data is such that the probability value of p < 0.05 we Valid - N Spearman - R t(N−2) p-value
say that the results are acceptable. Residuals &
254 0.032002 0.508281 0.611701
log(δ)
Table 3 The test results of the logarithm results-significance base Residuals &
254 0.024431 0.387948 0.698382
Analysis of Variance; DV: log deformation (Database) log(ρ)
Sums of - Mean - Residuals &
df F p-value 254 0.002687 0.042648 0.966016
Squares Squares log(pu)
Regress. 26.65913 3 8.886378 568.3461 0
Residual 3.90888 250 0.015636 4 MATHEMATICAL MODEL RESULTS
Total 30.56801
According to Tab. 3, the value of factor p is less than After making the logarithmic mathematical model and
0.05, and we conclude that the model is correct and reliable. performing the tests to verify the correctness of the applied
model, we conclude that it can be applied with a very high
3.2 Model Multicollinearity Test accuracy to predict deformation with known stress. The
condition is that the independent variables of the specific
Multicollinearity testing is performed to determine the mass and amount of polyurethane adhesive are known and
correlation of independent variables, or to show us what is unchangeable. This is ensured by the technological process
the interdependence of the same. of production for a given product. Based on the above-
mentioned requirements, the model defined the mathematical
Table 4 The results of the logarithm test results base-multicollinearity expression that we can use to make such calculations. By
Variables currently in the Equation; DV: log deformation studying the results of the multiple logarithmic model, we
(Crnoja database 25.05.2019.) came to the following mathematical expression:
Tolerance Logarithmic form of equation:
log(δ) 0.836293
log(ρ) (t/m3) 0.856332
log(pu) g/m2 0.97453 logε =
1.49 + 1.2996logδ − 1.289logρ − 0.134logpu (1)
3.3 Testing the Model for the Presence of Heteroscedasticity Where: ε – deformation, δ – stress (MPa), ρ – specific mass
(t/m3), pu – glue g/m2 (kg).
The presence of heteroscedasticity does not lead to the After defining the mathematical expression to determine
estimates obtained by the least squares method being biased, the value of the deformation value in relation to the given
but causes the ratings to have no minimal variation, i.e. are load, it is necessary to note the significant influence of certain
not effective. Hence, if there is heteroscedasticity, no areas in the test zone curve. These are the minimum and
prediction based on the original model’s ratings will be maximum values that have a significant impact on the actual
effective. The variation of the forecast, apart from the model. The actual model represents the test case from which
variation of residuals, includes the variation of parameter the model base was created. This means that the impact of
estimates. Therefore, it must be ensured that the model has the accuracy error we have noted initially depends to a great
no presence of heteroscedasticity [10]. extent on these values. In Tab. 6, we see the above values that
The condition that there is no heteroscedasticity is that have a significant impact on the mathematical model
the factor p > 0.05. In Tab. 5, we see that the condition is predicting the deformation size.
Table 6 Overview of the influential minimum and maximum from the model base 6 REFERENCES
Standard Residual: log log deformation (Database) Outliers
[1] European Tyre & Rubber Manufacturers’ Association. (2015).
Observed -
Predicted -
Standard -
Standard -
Std.Err. -
Residual
Residual
Pred.Val
Pred. v.
Number of End of life tyre report, http://www.etrma.org/uploads/
Value
Value
observation Modules/Documentsmanager/elt-report-v9a---final.pdf
[2] Marsili, L., Coppola, D., Bianchi, N., Maltese, S., Bianchi, M.,
& Fossi, M. C. (2014). Release of Polycyclic Aromatic
50 . . . |* . . . 0.59 0.30 0.28 -2.83 2.26 0.03
56 . . . | . * . . 1.51 1.09 0.43 -0.42 3.42 0.02
Hydrocarbons and Heavy Metals from Rubber Crump in
66 . . . * . . . 0.51 0.76 -0.25 -1.41 -2.02 0.02 Synthetic Turf Fields: Preliminary Assessment for Athletes.
89 . . . |* . . . 0.69 0.42 0.27 -2.48 2.15 0.02 Journal of Environmental and Analytical Toxicology, 5(2),
124 . .* . | . . . 0.37 0.86 -0.49 -1.12 -3.90 0.01 1000265. https://doi.org/10.4172/2161-0525.1000265
125 . . *. | . . . 0.74 1.15 -0.40 -0.23 -3.22 0.01 [3] Swedish Chemicals Inspecorate (2006). Synthetic turf from a
126 . . . * | . . . 0.91 1.25 -0.34 0.08 -2.73 0.01 chemical perspective – a status report 3/06 KEMI-
127 . . . * | . . . 1.06 1.38 -0.31 0.47 -2.50 0.01 Kemkalieinspektionen. Sundbyberg, Order. No. 510 834 July
128 . . . *| . . . 1.20 1.48 -0.28 0.79 -2.25 0.01 2006. e-mail: kemi@cm.se
129 . . . * . . . 1.28 1.54 -0.26 0.97 -2.07 0.01 [4] World Health Organization International Agency for Research
Minimum . .* . | . . . 0.37 0.30 -0.49 -2.83 -3.90 0.01
on Cancer - IARC (2010). Monographs on the Evaluation of
Maximum . . . | . * . . 1.51 1.54 0.43 0.97 3.42 0.03
Mean . . . |* . . . 0.89 1.02 -0.14 -0.62 -1.08 0.02
Carcinogenic Risks to Humans 1983. VOLUME 92 Some Non-
Median . . . *| . . . 0.83 1.12 -0.27 -0.33 -2.16 0.01 heterocyclic Polycyclic Aromatic Hydrocarbons and Some
Related Exposures, Lyon, France, 868.
[5] Simon, R. (2010). Review of the Impacts of Crumb Rubber in
Tab. 6 shows significant deviations of model values Artificial Turf Applications. UC Berkeley: Laboratory for
compared to actual measured values. As regards minimum Manufacturing and Sustainability. Retrieved from
and maximum, explanation of the occurrence is relatively https://escholarship.org/uc/item/9zp430wp
easy. We can notice that due to the nature of the test the initial [6] Birkholz, D. A., Belton, K. L., & Guidotti, T. L. (2012).
and final measurement are in the zone that is most difficult to Toxicological Evaluation for the Hazard Assessment of Tire
measure; the biggest error is in the initial deformation Crumb for Use in Public Playgrounds. Journal of the Air &
estimation and in the phase of plastic deformation or in the Waste Management Association, 53(7), 903-907.
area of material breakdown. https://doi.org/10.1080/10473289.2003.10466221
[7] Watterson, A. (2017). Artificial Turf: Contested Terrains for
The measurements were performed by a simple manual
Precautionary Public Health with Particular Reference to
stress method and the problem of recording the initial Europe? International Journal of Environmental Research and
deformation and deformation in the zone of plastic Public Health. https://doi.org/10.3390/ijerph14091050
deformation and breakdown of the material is extremely [8] http://gumiimpex.hr
difficult and subject to the subjective impression on the spot [9] HRN EN ISO 527-4:2008
at the time of measurement. [10] Nešić, N. (2014). Otkrivanje posledice prisustva
heteroskedastičnosti. Univerzitet u Novom Sadu,
5 CONCLUSION Prirodnomatematički fakultet, Departman za matematiku i
informatiku. Master rad, 90. (in Serbian)
The model is significant and shows with very high
precision the behavior of the material when it comes to the
Authors’ contacts:
stresses in real terms and the actual material that can be
obtained as a raw material on a certain geographical micro Anđelko CRNOJA, PhD student
location. Odessa State Academy of Civil Engineering and Architecture,
The accuracy of the applied multiple logarithmic model Didrihsona 4, 65029 Odessa, Ukraine
Tel: 00385992755466
is 87.21%. Therefore, we can conclude that the model is E-mail: acrnoja@hotmail.com
significant and applicable and has a very high accuracy of the
deformation estimation in relation to the applied stress. Željko KOS, PhD
The research will continue in the direction of examining Corresponding author
University North,
other properties such as the elasticity module etc. from which 31b, Jurja Krizanica, 42000 Varazdin, Croatia
information would be obtained that would be useful for Tel: 0038598757989
certain computer modelling when creating new products in a E-mail: zkos@unin.hr
wide range of human activities including civil engineering
Oleg POPOV, PhD, Assistant Professor
and construction industry. Department of Technology of Building Production,
The results obtained by applying a multiple logarithmic Odessa State Academy of Civil Engineering and Architecture,
model can be applied in defining and solving certain Didrihsona 4, 65029 Odessa, Ukraine
problems that require a response to the question of how Tel: 00380503339753
E-mail: oleg.a.popov@gmail.com
deformation is and is allowed if there is a certain stress effect.
Erol CAN
Abstract: The paper presents an application of power control by couple half wave LPWM rectifiers processing at the energy distribution line. The switching method with the
components of semiconductor for alternating energy control in the energy distribution line is used on alternating energy sources. In three-phase power line and load, power control
and regulation are done with the proposed modulation correction index (MCI). Therefore, the simulation model of 3-phase energy distribution is established after the relationships
between the energy sources and the circuit elements are determined by the circuit analysis method. The controls of the current and the voltage are done by trying at the different
loads in the line at the simulation. According to the modulation correction index (MCI), currents of unbalanced loads are balanced by the proposed method on the power line.
Then, the balance currents and powers that MCI provides are calculated. Finally, a power line of 15 km at length that is created by connecting the model of RLC is tested. Obtained
results show the effectiveness of the proposed method.
Keywords: balance currents; circuit analysis; modulation correction index; LPWM half rectifiers; simulation of power control
1 INTRODUCTION with making the circuit analysis of the power line. After the
design of the power line model, three-phase line operation is
As well as the production of energy, the distribution and performed in MATLAB Simulink for different load cases in
use of existing energy is of great importance. Therefore, the application phase.
many studies have been conducted on energy production and Then, obtained results of the simulation are verified in
distribution [1-3]. Direct current (DC) energy in power lines experimentation. The current and voltages of different loads
and different loads is controlled by converting it into without changing the sinus structure of current and voltage
changing energy with pulse width modulation (PWM) vary depending on the modulation index on the line. In the
inverters [4-9]. Therefore, the DC power source is derived simulation conducted for unbalanced loads, the impedance of
from either renewable energy sources or alternating current the third phase is half of the other phase impedances while
(AC) sources. The DC energy from the AC source is the two phases are equal to each other. Imbalance current of
converted back into AC current to control the power line. The the line with the low impedance value is balanced by
quality of the converted energy from the DC current regulating the modulation index value. These modulation
decreases since the current and voltage contain a high amount index values are named as modulation correction index
of harmonic distortion [10-12]. Multi-level inverters are used (MCI) for balancing the power and current values in the first
to eliminate or reduce these unwanted effects. This increases time. The imbalances due to different load conditions in the
the cost of the system and the complexity of the system while line are solved with different Modulation Correction Index
it makes it difficult to find the faults that can be on the device. values. These values are given in tables. When the results are
In addition, for some loads on the line controlled by the analyzed, the value of the voltage controlled on the load has
multi-level inverter, the quality of the current can be drawn lower distortion than that of the previous line studies with
to the desired level while the voltage of the generated energy inverters [15, 16]. In addition, these results can be obtained
cannot be obtained at the desired level [13, 14]. Thus, the in a less complicated and inexpensive way using the
paper presents a model of power control by inverse half proposed model. In addition, according to the results, the
rectifier processing method that uses the modulation power distribution line is simulated for efficient power
correction index at an energy distribution line. On the line, distribution and control without losing a lot of time and
there are two AC-source half-wave PWM rectifiers for each money. Then, the results are validated in experimentation. By
phase control current and voltage. Two GTO and two IGBT taking this model into consideration, applications can be
switches created are used on these half-wave rectifiers. This performed in less time and with less error rate. As the missing
provides a mixed switches structure at the PWM rectifiers and difficult sides of the power line in the simulation can be
while it is intended for demonstrating the availability of determined, the accidents and errors that may occur are
different accessible switch preferences when performing reduced and the loss of life and property is eliminated at the
energy control on the line. In the proposed control method, experimentation. According to the results obtained, it is
the effective values of the alternating current and voltage of observed that energy distribution and control are done
the load are controlled by PWM on the load with dividing as effectively.
the pulse width. Therefore, the modulation index values of
PWM control the effective values of the current and voltage
on the load without disturbing the sinus structures. In this
study, the structure of the proposed line model is described
Z A = Z 6 + R1 (4) 1
Z W= jω LL + RL + (13)
jωCC
The I1 is current of the first phase in ampere. I1 can be
expressed as in Eq. (5). Impedance (ZL) of the power distribution line can be
found as in Eq. (14).
DU m ⋅ sin(ωt )
I1 = (5)
ZA Z=
L ZA + ZW (14)
a)
b)
Figure 1 a) the power line of circuit model, b) LPWM
and distortion value of voltage are given in Fig. 3, while the index value. In Fig. 3b, the THD value of the voltage on the
current is as in Fig. 4. load is 4.49%. 104.7V is created with 0.95 of modulation
index on load. While the THD values of the voltage supplied
150 to the load with PWM inverters can exceed 100%, a level of
UZ5 UZ4 UZ6 4.59% can be achieved for the proposed method. In the
100 simulation, the current obtained is given at the 0.8 of
modulation index in Fig. 4.
For the balanced three-phase load that is ZA, ZB, ZC in
50
Fig. 4, a three-phase alternating current with 60° of the phase
difference is provided on load. 10.7 A is created with 0.95 of
0 modulation index on load. If loads in Fig. 2 are changed as
V
-150 20
0.05 0.06 0.07 0.08 0.09 0.1
t(s) 15
a) IZ5 IZ4 IZ6
10
Selected signal: 6 cycles. FFT window (in red): 2 cycles
100 5
50 0
0 A -5
V
-50 -10
-100 -15
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
Time (s) -20
Fundamental (60Hz) = 104.7 , THD= 4.59%
Mag (% of Fundamental)
15
5 IZ5 IZ4 IZ6
10
0 5
A
0
A
-5
-5
-10 -10
-15
-15
0.05 0.06 0.07 0.08 0.09 0.1
-20
t(s)
Figure 4 3-phase alternating current at 60 Hz of frequency for balanced load 0.05 0.06 0.07 0.08 0.09 0.1
t(s)
For the balanced three-phase loads that are Z4, Z5, Z6 in Figure 6 The balanced loads with the effective value of their current
Fig. 3a, a three-phase alternating voltage with 60 degree of
the phase difference is provided on load. These voltages are MCI is 0.5 according to Eq. 11. Therefore, when the
formed on the load according to arranging the modulation modulation index for IZ6 is made to 0.45, the modulation
index for IZ4 and IZ5 is 0.9. Then, IZ4 = IZ5 = IZ6 = 10.8A line imbalances. When the three-phase line is supplied with
for effective value of the current [7]. 330 volts; balancing currents and MCIs are given in Tab. 2
When the three-phase line is supplied with 220 volts; for different line imbalances. Z is 10 ohms as impedances.
balancing currents and MCIs are given in Tab. 1 for different
Table 1 Correction currents and MCIs of 3-phase lines fed with AC sources of 220V
D ZA ZB ZC MCI Imax – Ief for IZA and IZC Imax – Ief for IZB Ief f or IZA and IZC with MCI
0.8 Z 0.9Z Z 0.89 22A-17.6A 24.4A-19.6A 19.6A
0.8 Z 0.8Z Z 0.99 22A-17.6A 27.5A-22A 22A
0.7 Z 0.9Z Z 0.776 22A-15.4A 24.4A-17.08A 17.08A
0.7 Z 0.8Z Z 0.875 22A-15.4A 27.5A-19.25A 19.25A
0.7 Z 0.7Z Z 0.99 22A-15.4A 31.42A-22A 22A
0.7 Z 0.6Z Z X 22A-15.4A 36.6A-25.66A X
0.6 Z 0.9Z Z 0.666 22A-13.2A 24.4A-14.66A 14.66A
0.6 Z 0.8Z Z 0.756 22A-13.2A 27.5A-16.65A 16.65A
0.6 Z 0.7Z Z 0.857 22A-13.2A 31.42A-18.85A 18.85A
0.6 Z 0.6Z Z 0.99 22A-13.2A 36.66A-22A 22A
0.6 Z 0.5Z Z X 22A-13.2A 44A-26.4A X
0.6 Z 0.4Z Z X 22A-13.2A 55-33A X
Table 2 Correction currents and MCIs of 3-phase lines fed with AC sources of 330V
D ZA ZB ZC MCI Imax – Ief for IZA and IZC Imax – Ief for IZB Ief f or IZA and IZC with MCI
0.8 Z 0.9Z Z 0.89 33A-26.4A 36.6A-19.6A 29.4A
0.8 Z 0.8Z Z 0.99 33A-26.4A 41.25A-33A 33A
0.7 Z 0.9Z Z 0.776 33A-23.1A 36.6A-25.62A 25.62A
0.7 Z 0.8Z Z 0.875 33A-23.1A 41.25A-28.87A 28.87A
0.7 Z 0.7Z Z 0.99 33A-23.1A 47.3A-33A 33A
0.7 Z 0.6Z Z X 33A-23.1A 54.9A-38.49A X
0.6 Z 0.9Z Z 0.666 33A-19.8A 36.6A-21.99A 21.99A
0.6 Z 0.8Z Z 0.756 33A-19.8A 41.62A-24.97A 24.97A
0.6 Z 0.7Z Z 0.857 33A-19.8A 47.13A-28.27A 28.27A
0.6 Z 0.6Z Z 0.99 33A-19.8A 54.99A-33A 33A
0.6 Z 0.5Z Z X 33A-19.8A 66A-39.6A X
0.6 Z 0.4Z Z X 33A-19.8A 82.5A-44A X
Table 3 Balancing currents and balancing powers for three-phase system with 220 balanced in three cases. Although a current of 25.6 A is
volt input
required for the compensation of the line, the maximum
Impedance (Ohm) MI MCI Pb(VA) Ib(A)
0.9Z 0.8 0.89 440 2 current that the line can achieve is 22A at 0.99 of modulation
0.8Z 0.8 0.99 968 4.4 index. Current imbalance of the 6 ohm of the impedance
0.9Z 0.7 0.776 309 1.68 controlled by 0.7 modulation index cannot be compensated.
0.8Z 0.7 0.875 847 3.85 Thus, MCI is shown with X. Although a current of 26.4 A is
0.7Z 0.7 0.99 1452 6.6 required for the compensation of the line, the maximum
0.6Z 0.7 X X X
0.9Z 0.6 0.666 321,2 1.46
current that the line can achieve is 22A at 0.99 of modulation
0.8Z 0.6 0.756 759 3.45 index. Current imbalance of the 5 ohm impedance controlled
0.7Z 0.6 0.857 1243 5.65 by 0.6 modulation index cannot be compensated. So, MCI is
0.6Z 0.6 0.99 1936 8.8 shown with X for it.
0.5Z 0.6 X X X Although a current of 33 A is required for the
0.4Z 0.6 X X X
compensation of the line, the maximum current that the line
can achieve is 22 A at 0.99 of modulation index. Current
In Tab. 1, the imbalances generated in the impedances imbalance of the 4 ohm impedance controlled by 0.6
ranging from 10 ohm to 4 ohm on the line are controlled by modulation index cannot be compensated. Thus, MCI is
the modulation index that is from 0.8 to 0.6. While the shown with X for it. In Tab. 2, the imbalances generated in
effective value of the line current is balanced twice in the 0.8 the impedances ranging from 10 ohm to 4 ohm are controlled
modulation index, the effective value of the line current can on the line from 0.8 modulation index to 0.6 modulation
be set at three times for the 0.7 of modulation index value. In index. The line can be balanced in nine of twelve trials, while
the case of 0.6 modulation index, the effective value of the the line cannot be balanced in three cases. Tab. 3 shows the
line current can be compensated four times. For the 0.9Z balancing current and balancing power generated by
value of ZB at the 0.8 of the modulation index, the MCI value arranging of the MCI to compensate for the three-phase
is 0.89. For the 0.8Z value of ZB at the 0.8 of the modulation power line, which has 220V of the input voltage source. Tab.
index, the MCI value is 0.99. At the 0.7 of the modulation 4 has the balancing current and balancing power generated
index, the MCI value is 0.776 for the 0.9Z value of ZB while by arranging of the MCI to compensate for the three-phase
the MCI value is 0.875 for the 0.8Z value of ZB. The line can power line, which has 330V of the input voltage source. The
be balanced in nine of twelve trials, while the line cannot be balance current (Ib) is the difference between the high phase
current and the low phase current values. This value can be impedance value (Z) is 10 ohm, the balancing current and
expressed as follows. power cannot be provided with 0.6 of modulation index to
the system supplied with 220 volts for 0.5Z and 0.4Z
I b I ZA − I ZB
= (14) impedance values. When the impedance value (Z) is 10 ohm,
the balancing current and power cannot be provided with 0.7
I b I ZA − I ZB
= (15) of modulation index to the system supplied with 220 V for
0.6Z of the impedance values.
The balance power (Pb) can be found as in Eq. (16) while Table 4 Balancing currents and balancing powers for three-phase system with 330
Ui is input voltage for the line. volt input
Impedance (Ohm) MI MCI Pb(VA) Ib(A)
P= Ui ⋅ Ib (16) 0.9Z 0.8 0.89 880 4
b 0.8Z 0.8 0.99 1936 8.8
0.9Z 0.7 0.776 618 336
For the impedance of 0.9Z at Tab. 2, a compensation 0.8Z 0.7 0.875 847 7.7
power of 440 VA and a balancing current of 2 A are provided 0.7Z 0.7 0.99 2904 13.2
with a 0.89 MCI value when a compensation power of 0.6Z 0.7 X X X
0.9Z 0.6 0.666 642 2.92
968VA and a balancing current of 4.4 A are provided with a 0.8Z 0.6 0.756 1518 6.9
0.99 of MCI. The first greatest balancing power and current 0.7Z 0.6 0.857 2486 11.3
are 1936 VA and 8A that are provided with 0.99 of MCI at 0.6Z 0.6 0.99 3872 17.6
the 0.6Z of impedance, while the second greatest balancing 0.5Z 0.6 X X X
power and current are 1452 VA and 6.6 A that are provided 0.4Z 0.6 X X X
with 0.99 of MCI at 0.7Z of the impedance. When the
For the impedance of 0.9Z at Tab. 4, a compensation with 0.99 of MCI at the 0.6Z of impedance, while the second
power of 880 VA and a balancing current of 4 A are greatest balancing power and current are 2904 VA and 13.2
provided with a 0.89 MCI value when a compensation A that are provided with 0.99 of MCI at 0.7Z of the
power of 968VA and a balancing current of 8.8 A are impedance. As the unbalanced phase impedance decreases,
provided with a 0.99 of MCI. The first greatest balancing the applied MCI and the balancing current increase, which
power and current are 3872 VA and 8 A that are provided leads to an increase in the balancing power. When the
impedance value (Z) is 10 ohm, the balancing current and Parallel RL loads connected in series at the beginning and
power cannot be provided with 0.6 of modulation index to end of the 15 km long line are energized with alternating
the system supplied with 330 V for 0.5Z and 0.4Z source, a three-phase alternating current that occurs on the
impedance values. In line with the model in Figure 1, a line load is presented in Fig. 9.
length of 15 km in Fig. 7 can be obtained if a serial RLC
circuit is added to the series. Z1 = Z2 = Z3 = Z4 = Z5 = Z6 15
A
with 0.9 of modulation index. The three-phase alternating -5
voltage that occurs at the end of the power line is as in Fig.
8a, while the harmonic distortion of the voltage on the load -10
is as in Fig. 8b.
-15
0.05 0.06 0.07 0.08 0.09 0.1
150 t(s)
UZ5 UZ4 UZ6 Figure 9 A three-phase alternating current of the 15 km long line
100
On the power line, a 12-volt alternating current at the
50 maximum value is formed on loads. The effective value of
the voltage is 10.2 A because the modulation index is 0.9.
According to the results, the power distribution line is
0
V
5 REFERENCES [15] Raj, N., Jagadanand, G., & George, S. (2016). A Modified
Charge Balancing Scheme for Cascaded H-Bridge Multilevel
[1] Krauter, S. (2018). Simple and effective methods to match Inverter. Journal of Power Electronics, 16(6), 2067-2075.
photovoltaic power generation to the grid load profile for a PV https://doi.org/10.6113/JPE.2016.16.6.2067
based energy system. Solar Energy, 159, 768-776. [16] Zheng, C. F., Zhang, B., Qiu, D. Y., Zhang, X. H., & Xiao, L.
https://doi.org/10.1016/j.solener.2017.11.039 M. (2015). Wavelet PWM Technique for Single-Phase Three-
[2] Niu, B., Hwangbo, H., Zeng, L., & Ding, Y. (2018). Level Inverters. Journal of Power Electronics, 15(6), 1517-
Evaluation of alternative power production efficiency metrics 1523. https://doi.org/10.6113/JPE.2015.15.6.1517
for offshore wind turbines and farms. Renewable Energy, 128,
81-90. https://doi.org/10.1016/j.renene.2018.05.050
[3] Exel, L. & Frey, G. (2018). Modeling and simulation of local Author’s contacts:
flexibilities and their effect to the entire power
Erol CAN, PhD, Assoc. Prof.
system. Computer Science-Research and Development, 33(1- Aviation Electric Electronics, School of Civil Aviation,
2), 6. https://doi.org/10.1007/s00450-017-0346-7 Erzincan Binali Yıldırım University,
[4] Matsumori, H., Shimizu, T., Takano, K., & Ishii, H. (2018). Fatih Avenue, 24100 Erzincan, Turkey
Three‐Phase AC Filter Inductor Design for Three‐Phase cn_e@hotmail.com
PWM Inverter for Conversion Efficiency Improvement at
Low Load. Electrical Engineering in Japan, 203(1), 37-49.
https://doi.org/10.1002/eej.23053
[5] Ananda, A. S. (2018). Performance Analysis of Series-Passive
Filter in 5-Phase PWM Inverter Drive and Harmonic Study
Using Simulink/Matlab. In Advances in Power Systems and
Energy Management (pp. 139-146). Springer, Singapore.
https://doi.org/10.1007/978-981-10-4394-9_14
[6] Kumar, M. (2018). Time-Domain Characterization of
Digitized PWM Inverter with Dead-Time Effect. IEEE
Transactions on Circuits and Systems I: Regular Papers.
https://doi.org/10.1109/TCSI.2018.2827950
[7] Mukherjee, S., Giri, S. K., Kundu, S., & Banerjee, S. (2018).
A Generalized Discontinuous PWM Scheme for Three-Level
NPC Traction Inverter with Minimum Switching Loss for
Electric Vehicles. IEEE Transactions on Industry
Applications. https://doi.org/10.1109/TIA.2018.2866565
[8] Takahashi, H., Obara, H., & Fujimoto, Y. (2018, March).
Dead time compensation for three-level flying capacitor
inverter with phase shift PWM. In Advanced Motion Control
(AMC), 2018 IEEE 15th International Workshop on (pp. 229-
233). IEEE. https://doi.org/10.1109/AMC.2019.8371093
[9] Sivakumar, P., & Arutchelvi, M. S. (2018). Modified
composite power control strategy for grid connected wind‐PV
systems with unbalanced nonlinear current. International
Transactions on Electrical Energy Systems, e2587.
https://doi.org/10.1002/etep.2587
[10] Can, E., & Sayan, H. H. (2016). SSPWM three phase inverter
design and experimented on unbalanced loads. Tehnički
vjesnik, 23(5), 1239-1244.
https://doi.org/10.17559/TV-20150730222021
[11] Can, E. & Sayan, H. H. (2017). The increasing harmonic
effects of SSPWM multilevel inverter controlling load
currents investigated on modulation index. Tehnički
vjesnik, 24(2), 397-404.
https://doi.org/10.17559/TV-20151020134629
[12] Can, E. (2018). The modeling and analysis of a power
transmission line supplied by a solar power plant. Tehnički
glasnik, 12(3), 124-130.
https://doi.org/10.31803/tg-20180521111744
[13] Singh, S. N. & Singh, G. K. (2017). Modelling, design and
stability analysis of an improved SEPIC converter for
renewable energy systems. International Journal of
Electronics, 1-19.
[14] Pan, L., et al. (2017). One-cycle control for three-phase single-
capacitor Z-source inverter with unity power factor.
International Journal of Electronics, 104(4), 635-658.
https://doi.org/10.1080/00207217.2016.1242161
Željko BARAČ, Ivan PLAŠČAK, Tomislav JURIĆ, Pavo BALIČEVIĆ, Vinko DUVNJAK, Mladen JURIŠIĆ, Goran HEFFER,
Monika MARKOVIĆ
Abstract: The paper presents a study of the influence of an uneven surface and number of tractor hours on vibration that is transmitted to the operator's hand-arm system during
operation. The study was conducted on the asphalt, gravel and grass traversal during 2015 and 2016. The results indicate that the movement of a tractor on different agrotechnical
surfaces generates vibrations of varying intensity that affect the hand-arm system of the operator. The highest vibration values were recorded on the asphalt, while the least
vibrations were measured on the grass. The measured maximum vibration values are lower than the permissible 5 m·s−2 in the direction of all three axes, making it possible to
assume that they are not harmful to the health of the operator.
vibration transferred to the driver (e.g. low-pressure tyres, air 3 RESULTS AND DISCUSSION
sprung seat and cab, and front axle suspension).
Barač et al. [16] measured the vibration that affects the Compared to the results obtained in 2015, the results for
operator's hand-arm system when working with a sprayer and 2016 show a slight vibration deviation in the direction of all
mulcher. The lowest vibration values were measured on a three axes measured on all three types of surfaces (gravel,
sprayer in the direction of the x, y and z axes. Barač et al. [17] asphalt, and grass surface) (Tab. 1).
researched the impact of various agrotechnical surfaces
Table 1 Mean values (m·s−2)
(asphalt, gravel and grass surface) on the level of vibration
x axis y axis z axis
that affect the hand-arm system. The results show that the Surfaces Year
lowest vibration values were obtained for the grass surface in 2015 2016 2015 2016 2015 2016
the direction of all axes, whereas the highest values were Asphalt 0.079 0.081 0.071 0.073 0.069 0.071
obtained on the asphalt and gravel surfaces in the direction of Gravel 0.060 0.063 0.071 0.074 0.037 0.070
x and y axes, and the y axis, respectively. Grass 0.020 0.023 0.069 0.073 0.040 0.044
The aim of the research is to establish the level of
The highest mean vibration values were measured in the
mechanical vibration that affects the operator’s hand-arm
direction of x and z axes on the asphalt surface in 2015 and
system on various agrotechnical surfaces in relation to an
2016. The highest mean vibration values in 2015 were
increased amount of the tractor’s work hours. It is expected
obtained in the direction of the y axis on the asphalt and
that the increase in work hours will result in a heightened
gravel road, while the highest mean values were measured on
level of vibration that affects the hand-arm system of the
the gravel road in 2016. The lowest vibration values in the
driver of an agricultural tractor.
direction of the x axis were measured as follows: in 2015 and
2016 when the tractor moved along a grass surface; in 2015
2 MATERIALS AND METHODS
in the direction of the y axis along a grass surface; in 2016 on
the asphalt and grass surface in the direction of the z axis; in
The research was performed on a Landini Powerfarm
2015 on the gravel road; in 2016 on the grass surface.
100 tractor through a two-year period. Over the first and
Comparable results were obtained by the authors [23],
second research year (2015 and 2016), the tractor reached 5
who measured the vibrations produced by framesaws and
800 and 6 800 hours of work, respectively. The
bandsaws. The authors established that the obtained vibration
measurements were conducted on the agricultural fields and
level does not carry any health risk to the operator positioned
access roads of the School of Agriculture and Veterinary
on the framesaw stand even after an eight-hour exposure
Medicine in Osijek to establish the values of vibration
period. However, a one-hour exposure results in decreased
influencing the hand-arm system of the operator (in the
working conditions, which can afflict workers' wellbeing and
direction of the x, y and z axes). The tractor moved along
their efficiency. Deboli et al. [24] measured the vibrations on
gravel, asphalt, and grass surfaces. The measurement of
four tractors equipped with several types of tyres on three
vibration was repeated three times, which produced the
different surfaces, i.e. gravel road, asphalt, and a combination
mean value of vibration used for further research. Each
between asphalt and gravel road. In the last case, two tyres
measurement lasted for 30 minutes.
drove along the gravel road, and the other two along the
Vibration level measurement was performed in
asphalt road. Irrespective of the tyre type, the lowest
accordance with the standards [18, 19, 20], which define the
vibration was obtained when the tractor moved along the
limits of vibration exposure, their effects on the operator's
asphalt surface. The vibration level in this case did not reach
health, and include the guidelines for measurements at the
the allowed limit values, which is in accordance with the
workplace. The measurements were performed with a
results of our research. The results of the descriptive statistics
measuring device MMF VM30. A sensor was mounted on
related to the mean values of the obtained vibration in the
the steering wheel with the coordinates set in the direction
direction of the x and z axes on the asphalt surface (Tab. 2)
[20]: x axis: longitudinally (along the axis of motion –
show that the standard error was higher in the second year of
forward (positive) / backward (negative); y axis: laterally
measurement (2) compared to the first year (1), whereas it
(at the right angle to the motion direction); and z axis:
remained the same in the direction of the y axis for both
vertically, upward (positive) / downward, vertically to the
years. Furthermore, the variance analysis of the mean values
surface (negative). The weather conditions in 2015 were as
for the obtained vibration impacting the hand-arm system of
follows: air temperature ranged between 29° C and 31 °C;
the driver in the direction of all three axes (x, y and z) on the
relative humidity ranged from 62 % to 64 %; the impact of
asphalt surface shows no statistically significant difference
wind was minimal. The weather conditions in 2016 were
between the measurements in both the first and second year.
similar.
Tab. 3 shows the descriptive statistics related to the mean
The directives [21] and OG [22] contain both the limit
values of the obtained vibration in the direction of x and y
values and critical values for the exposure of the hand and
axes on the gravel surface. It is discernible that a higher
arm and they are as follows:
standard error occurs in the second (2) year of measurement
• The limit value for daily exposure based on the reference
compared with the first year (1). Variance analysis of the
time of eight hours is 5 m·s−2;
mean values of vibration in the direction of the x, y and z axes
• The critical value for daily exposure based on the
on the gravel surface shows no statistically significant
reference time of eight hours is 2.5 m·s−2.
difference concerning the mean value of the obtained determined in the first (1) measurement year. The variance
vibration in both the first and second measurement year. analysis of the mean values of vibration in the direction of
The highest standard error related to the obtained the x, y and z axes on grass impacting the hand-arm system
vibration was measured in the direction of the x and z axes on did not reveal any statistically significant differences
the grass surface in the second (2) measurement year (Tab. between the first (1) and second (2) measurement year.
4). A higher standard error in the direction of the y axis was
Table 2 Descriptive statistics and variance analysis for occurred vibrations on asphalt surface
N x σ C.V. (%) st. error F Sig. (95 %)
x axis
2015 3 0.07967 0.004509 5.65 0.002603
2016 3 0.08100 0.004583 5.65 0.002646 0.129 0.738
Total 6 0.08033 0.004131 5.14 0.001687
y axis
2015 3 0.07133 0.006110 8.56 0.003528
2016 3 0.07333 0.006110 8.33 0.003528 0.161 0.709
Total 6 0.07233 0.005574 7.70 0.002275
z axis
2015 3 0.06900 0.005568 8.06 0.003215
2016 3 0.07133 0.006028 8.45 0.003480 0.243 0.648
Total 6 0.07017 0.005345 7.61 0.002182
Table 3 Descriptive statistics and variance analysis for occurred vibrations on gravel surface
N x σ C.V. (%) st. error F Sig. (95 %)
x axis
2015 3 0.06033 0.002517 4.17 0.001453
2016 3 0.06300 0.003000 4.76 0.001732 1.391 0.304
Total 6 0.06167 0.002875 4.66 0.001174
y axis
2015 3 0.07100 0.004583 6.45 0.002646
2016 3 0.07400 0.005568 7.52 0.003215 0.519 0.511
Total 6 0.07250 0.004848 6.68 0.001979
z axis
2015 3 0.03667 0.032146 87.66 0.018559
2016 3 0.07000 0.020000 28.57 0.011547 2.326 0.202
Total 6 0.05333 0.030111 56.46 0.012293
Table 4 Descriptive statistics and variance analysis for occurred vibrations on grass surface
N x σ C.V. (%) st. error F Sig. (95 %)
x axis
2015 3 0.02067 0.003055 14.77 0.001764
2016 3 0.02333 0.003512 15.05 0.002028 0.985 0.377
Total 6 0.02200 0.003286 14.93 0.001342
y axis
2015 3 0.06900 0.006557 9.50 0.003786
2016 3 0.07333 0.005508 7.51 0.003180 0.768 0.430
Total 6 0.07117 0.005913 8.30 0.002414
z axis
2015 3 0.04000 0.002000 5 0.001155
2016 3 0.04400 0.004583 10.41 0.002646 1.920 0.238
Total 6 0.04200 0.003847 9.15 0.001571
4 CONCLUSION asphalt and gravel surfaces in 2015, and for the gravel surface
in 2016;
The following conclusions can be drawn based on our - The lowest vibration levels were obtained as follows: for
two-year research: the grass surface in the direction of the x axis in 2015 and
- The level of vibration that impacts the hand-arm system 2016; in the direction of the y axis in 2015; in the direction
for all agrotechnical surfaces in the direction of the x, y and z of the z axis in 2016. The lowest level for the asphalt surface
axes was higher in the second (2) measurement year was obtained in the direction of the y axis in 2016 (the same
compared to the first (1) measurement year, which is in result was obtained for the grass surface). The lowest level
accordance with the specified hypothesis; for the gravel surface was obtained in the direction of the z
- The highest levels of vibration were obtained for the axis in 2015;
asphalt surface in the direction of the x and z axes (in 2015 - The research shows that different agrotechnical surfaces
and 2016), as well as in the direction of the y axis for the (i.e. asphalt, gravel and grass surfaces) produce vibration of
various intensity levels that are transmitted to the hand-arm
Authors’ contacts:
THE EFFECTS OF DIELECTRIC VALUES, BREAST AND TUMOR SIZE ON THE DETECTION OF
BREAST TUMOR
Abstract: Although breast cancer is the second main cause of female deaths after lung cancer, early diagnosis plays a crucial role to diminish the death rate. Many techniques
have been improved to detect the cancerous cells. At different microwave frequencies, the malignant cells indicate different electrical characteristics as compared to the normal
cells. According to these frequencies, the breast tissue is more permeable than other tissues such as the brain and muscle. Due to this property of the breast tissue, microwaves
can be used for the detection of breast cancer. In this study, the breast prototype was modelled using the CST STUDIO SUITE electromagnetic simulation software with respect
to different breast size, tumor size and dielectric values tested at a range of the 0-3.0 GHz frequency. The objective of this paper is to investigate the effects of each factor and the
interactions of factors on detecting cancer cells using the factorial analysis. The results indicate that the factors such as fat and skin permittivity, tumor and breast sizes are more
effective in the detection of breast tumor. Although the effect of fibro permittivity is not significant alone, there are considerable interaction effects of a large breast size and small
tumor size through low-to-high values of fibro permittivity. Furthermore, the combinations of a breast radius smaller than almost 8.5 cm with a high level tumor radius and breast
radius larger than 8.5 cm with a low level tumor radius are desirable for lessening the return loss value.
Keywords: dielectric properties; factorial design; factor effects; microwave breast cancer detection
malignant tissue is increased because of the displacement and Design of Experiment (DOE) techniques have been used
rotations of this mobile charge by the microwave field. extensively in many different scientific fields. These
The relaxation times which are a third factor in malignant techniques are beneficial for getting accurate statistical
tissues are larger than those in the normal tissues resulting inferences for revealing the effects of factors on a process
from the increased motional freedom of water [13]. The last with less experimenting rather than considering all the
one is sodium concentration and water content [13]. The combinations of factors. One of the DOE technique
sodium concentration in cancer cells is higher than that in applications is found in [16] for analytical chemistry. In
normal cells. The high sodium concentration alone yields another study [17], a considerable decrease of patient
higher permittivity and conductivity values for malignant dissatisfaction was given for the emergency room. [18]
tissue. This higher sodium concentration also affects the cell applied the factorial design technique to see the effects of
membrane potential and causes the malignant tissue to retain prophylactic antiemetic interventions on the risk of
more fluid. This retained fluid is in the form of bound water, postoperative nausea and vomiting of the surgery patients.
which has higher permittivity and conductivity values. Another case study on the effect of treatment factors on the
Malignant tissue therefore has a higher water content ratio quality of life of patients with breast cancer is found in [19].
than that of normal tissue, which coincides with higher values In this study, 2k factorial design is used to reveal the effects
of permittivity and conductivity than normal breast tissue at of factors and their interactions on the success of detecting
the same microwave frequency. breast tumor.
The understanding and development of medical The outline of the paper starts with the introduction
microwave techniques needs knowledge of the microwave section. It is followed by illustrating the factors considered in
dielectric properties of human tissues. There have been many breast cancer detection and applying simulation and the
data on the dielectric properties of female human general factorial design. In the third section, significant
normal/malignant/benign and surrounding breast tissue. effects of factors and their interaction graphs are given.
Dielectric properties of breast carcinoma and the Finally, further suggestions on the improvement of breast
surrounding tissues were measured [14]. The experiments cancer detection are explained in the conclusion part.
were performed at the frequencies from 100 kHz to 100 MHz
at 37 °C using an automatic network analyzer and an end-of- Table 1 Computer fitted dielectric parameters of breast carcınoma and the
surrounding tissues in the frequency range from 100 kHz to 100 MHz [14]
the-line capacitive sensor. Using a computer program, the
Specimen Sample σDC
cole-cole parameters were found by curve fitting. Tissues Location No.
εδ ε∞ τ (μs) α
(ms/cm)
Diagnosis
were divided into three main categories. These were the 1 6090 50 0.633 0.326 2.20
central part of the tumor, the tumor surrounding tissue, and 2 8260 50 0.595 0.311 3.03 Infiltrating
the peripheral tissue. Tumor specimens were taken from 3 3720 50 0.267 0.283 6.47 ductal
A
4 5170 50 0.288 0.239 6.91 carcinoma
seven different patients. From each specimen, depending on 5 3440 50 0.238 0.238 7.12 metastases
its size, 2-6 samples in the form of thin disks were excised. 6 5710 50 0.263 0.277 6.60
The number of samples was 28. The duration of all Infiltrating
1 3330 50 1.30 0.380 3.3
measurements was 4 h after the surgery. As seen from Table B
2 5510 30 0.33 0.250 5.70
ductal
1, measurements showed significant differences in dielectric carcinoma
1 793 20 0.579 0.232 0.82 Infiltrating
properties between the samples taken from different C 2 482 30 0.603 0.281 0.38 ductal
locations. Surgical pathology reported the following 3 3310 30 0.557 0.310 3.05 carcinoma
diagnoses: A-infiltrating ductal carcinoma without 1 2740 30 0.269 0.358 3.89
Infiltrating
metastases, B-infiltrating ductal carcinoma, C-infiltrating 2 3285 30 0.3202 0.275 5.75
D lobular
3 2320 30 0.259 0.278 6.49
ductal carcinoma with focal papillary component, D- 4 1253 319 0.670 0.400 3.52
carcinoma
infiltrating carcinoma consistent with lobular carcinoma, E- 1 3125 167 0.525 0.223 1.02
Infiltrating
infiltrating lobular carcinoma, F-infiltrating ductal E
2 7060 30 1.605 0.362 4.99
lobular
carcinoma, G-infiltrating lobular carcinoma. 3 13100 30 1.396 0.349 6.59
carcinoma
4 3080 50 0.45 0.298 3.90
Using a resonant cavity technique, complex permittivity 1 8830 30 1.408 0.360 2.78 Infiltrating
of in vitro diseased and non-diseased human female breast F 2 10170 30 1.210 0.340 2.73 ductal
tissues were measured at 3.2 GHz [15]. Both dielectric 3 180 10 0.684 0.149 0.120 carcinoma
properties and water contents of these tissues were found. 1 288 10 0.983 0.10 0.161
Experimental data were compared with the models predicted 2 298 10 1.008 0.075 0.130
Infiltrating
3 675 10 1.391 0.163 0.444
from mixture equations. G
4 7930 50 1.216 0.367 5.33
lobular
Hitherto, breast cancer detection studies were carried out carcinoma
5 5490 50 0.768 0.379 4.15
using dielectric properties; however the effect of dielectric 6 8660 50 1.415 0.331 5.16
properties has not been investigated. A 3D-breast was
modeled by using the CST STUDIO SUITE electromagnetic 2 SIMULATION MODEL OF BREAST TUMOR
simulation software and analyzed by using different
dielectric properties and different breast and tumor sizes A hemisphere breast model was designed as showed in
based on experimental design analyses results in the Fig. 1 according to the dielectric properties of tissues which
frequency range from 0 GHz to 3 GHz. are permittivity (ε) and conductivity (σ). Microwave imaging
was carried out by transmitting a sequence of
electromagnetic waves through the breast model and factors A, C, J and K are less than 0.05 and they can be
measuring the scattered field on the breast. concluded as having significant effects with the confidence
of 0.95. Thus, Skin ε, Fat ε, BreastRadius and TumorRadius
are obtained as main concerns in the detection of breast
cancer.
BreastRadius and the small size of TumorRadius through the less than for the small size. However, this interpretation
low-to-high values of Fibro ε. becomes reverse after the 8.5 cm size BreastRadius.
Finally, in the last graph, while BreastRadius is smaller However, the values of S11 for the large size TumorRadius
than 8.5 cm, S11 values for the large size TumorRadius are are greater than the small size.
By using the numerical optimization tool of the Design The desired value of S11 is -21.8547. It is seen that the S11
Expert, the minimum value of S11 is found. In order to have value decreases at the given values of factors in Fig. 5, and
a minimum S11 value, the optimal design points are given in the increasing values of breast radius size and decreasing
Fig. 5. values of tumor size considerably help the detection of breast
The desirability graph based on two most effective tumor.
factors on breast tumor detection is demonstrated in Fig. 5.
Author contact:
Abstract: This paper focuses on the structure and performance of the pressure swirl nozzle and the study of liquid atomization. In this study, the atomizer has been designed and
some experiments have been performed on it. Since image processing is an efficient method for measuring the size of the droplet and since it considerably reduces the total
measuring time and eliminates the subjective observer’s error in sizing and counting spray drops, a digital camera has been used for capturing images and image processing has
been done by the MATLAB software. The results show that by increasing the atomization air pressure, the spray angle increases and the droplet’s size decreases. It is concluded
that the spray angle is a function of the atomization air pressure and orifice diameter. Moreover, when the distance from the spray centre line increases, the droplet’s average
velocity decreases.
Keywords: atomization; droplet size; image processing; pressure swirl nozzle; spray
darkness or by a black sheet at the back of spray. Moreover, where it gains a swirl motion, under which it leaves the exit
an investigation of spray characteristics such as the liquid orifice as a conical liquid sheet which consequently breaks
breakup length and spray cone angle of a charge-injected up into small droplets due to aerodynamic forces. The
electrostatic pressure-swirl nozzle has been conducted. The centrifugal motion of the swirling liquid makes a low-
liquid breakup length decreased, while the spray angle pressure area in the swirl chamber centre and generates an air
increased with an increase in the applied voltage and core along the centreline. The flow pattern inside the
injection pressure. An empirical equation to predict the atomizer is rather complex; it is a two-phase flow with
breakup length for the electrostatic pressure-swirl nozzle has secondary flow effects. There is a strong link between the
been suggested. The experimental result was within the range resulting spray characteristics and internal flow conditions;
of the predicted equation [6, 8].The above-mentioned two however, not all aspects of the internal flow are well realized
techniques are deducing the particle size information from [16]. Moreover, the pressure swirl and air-blast atomizers are
the optical signals scattered from an individual particle or a the two types of atomizers used in gas-turbine engines [17].
group of particles. In measuring volume, it is assumed that In recent years, some researchers studied the characteristics
all particles are spheres. Thus, basically, only the spherical of atomization for pressure-swirl atomizers of two various
particles should be processed with other techniques for geometries. The important effect of atomizer construction on
accurate measurement. On the other hand, in principle, the atomization process has been observed. Furthermore, a
various non-spherical particles can be processed through the guided experimental study reported that they aimed to
image processing technique because it is based on direct intensify the process of atomization and improve the
visualization and wide application is possible. Moreover, the properties of the spray [18].
measurement accuracy of the image processing techniques is
relatively insensitive to the optical properties of the particles 2 ATOMIZER DESIGN AND CONSTRUCTION
in comparison with the other techniques and optical
alignment is much easier [9]. In some cases, the simulation The two phases of atomizing including the following
program can be used in combination with calculation models stages and different sections of the atomizer have been shown
to predict the drop size depending on the hollow cone in Fig. 1.
nozzle's geometry and the volume flow for the atomization (1) Primary mixing of liquid and gas in the distributer
of Newtonian fluids [10]. In a new paper, the viscous flow (2) Secondary mixing in the mixing chamber
through the swirl chamber of a pressure-swirl atomizer has (3) Severe mixing in the spiral tube
been studied theoretically [11]. Additionally, the thermal (4) Liquid film formation in the cone cap
effects on the spray cone formation of a pressure swirl nozzle (5) Final atomization in the output of the nozzle tip.
in spray cooling have been investigated experimentally by
particle image velocimetry (PIV) [12]. Section (1) consists of two concentric tubes. Along the
The main purpose of the image processing is to develop internal tube, there are some holes to mix the liquid feed and
a simple appropriate algorithm to size and count the number atomizer gas. The liquid feed (water or vacuum gasoil) and
of spherical spray drops in the image frames, and emphasis the atomizer gas (air or steam) enter into the internal tube and
was placed on eliminating undesirable objects from the annular space, respectively, and mix together through the
image. Such objects included odd-shaped foreign materials holes existing on the internal tube. Then, the two-phase
and drops in contact or overlapping and parts of single drops mixture enters the section (3) which consists of a cylindrical
cut off by the image frame boundaries [6]. Therefore, this spiral surrounded by a connecting tube. This spiral produces
paper investigates the pressure swirl nozzle property and the a homogeneous mixture by circulatory movement and hard
application of a new image processing method for mixing at the opening of the orifice. Thus, the homogeneous
determining the droplet size distribution and droplet velocity mixture is atomized to a hollow conical spray after passing
distribution. The pressure of the liquid and gas mixture is the circular orifice. The orifice size and design have a main
measured in a mixing chamber by a pressure gauge. role in the droplet size and droplet velocity at the nozzle
outlet. The parameters for the design of the pressure swirl
1.1 Usage and Performance of Pressure Swirl Nozzle atomizer and operational variables are summarized in Tab. 1.
Pressure-swirl nozzles are used in a broad range of Table 1 Summary of the design parameters for the pressure swirl nozzle.
industrial uses, e.g. washing, cooling system, combustion Design data Spray Angle 60°
Length of the distribution tube 0.14 m
chamber, heavy fuel injection, painting and food processing, Number of holes in the distribution tube 15
etc. Their spray characteristics are approximately linked to Holes diameter 0.00075 m
the internal flow which predetermines the values of the liquid Length of the connecting tube 0.8 m
sheet formed at the outlet and discharge orifice. Pressure- Length of spiral 0.8 m
swirl atomizers are easy to make and prepare good Orifice shape conical
Cone angle 60o
atomization. They are frequently used in different Orifice hole diameter 0.001 m
applications where a large surface area of droplets is needed Orifice hole depth 0.005 m
or a surface must be coated by a liquid, e.g. combustion, fire Pressure of mixture 2 bar
suspension or air conditioning. In principle, the pumped Maximum liquid flow rate 6.5 L/h
liquid is injected via tangential ports into a swirl chamber Maximum gas flow rate 205 L/h
The schematic drawing of the pressure swirl nozzle and After the investigation and modification of images, they are
its design characteristics are shown in Fig. 1. Here the processed by a R2014a version of the MAT LAB software.
pressure swirl atomizer is designed for gas oil having the
maximum flow rate of 6.5 L/s, the average injection pressure
of 2.5 bar and a half spray cone angle of 30°, and the material
is made from stainless steel for the prevention of corrosion.
Figure 2 The schematic representation of the pressure swirl nozzle and the
performance analysis system
3 THE EXPERIMENTAL SET-UP The schematic image of spray droplets after photography
is shown in Fig. 3
The experimental set-up is shown in Fig. 2 .The air with
high pressure exits from the compressor, passes through the
tube and enters into the liquid pressure vessel. In Fig. 2, the
schematic images of various equipment used in the pressure
swirl nozzle’s performance analysis are shown.
When the sprayed stream exits the nozzle orifice, the
images of produced droplets are captured by a digital camera.
The initial image includes RGB image and grey levels Thus, the length and width of the centre of each drop are
(between 0 and 1) must be converted to the BW image. Each derived from Eq. (3).
image must be bmp (2D format) before processing.
Therefore, there is a matrix from the image including 0 and n n
1. Non-isolated droplets (i.e. droplets in overlap condition, or=
=
X =
∑ xi
i 1=
, Y
∑ y
i 1 i
(3)
droplets in contact) are all eliminated and the holes (empty n n
spaces) of drops are filled by the same pixels.
For segregation or better separation of pixels and proper Where xi and yi are the length and width of the droplet
image investigation, the Sobel function in MATLAB is used. pixels, respectively. n is number of pixels that include the
In Fig. 7, the Sobel modification method for accurate analysis droplet.
is shown. After this step, the image processing operation
starts and the information is derived. 4.3 Determination of the Spray Angle
a
θ = 2 arctan (4)
Figure 7 Image for displaying the Sobel method performance b
The Sobel method is a linear algorithm of edge detecting, The spray angle, θ, depends on the air pressure and air
which determines the boundaries of the light and dark points flow rate. In dilute solutions, the air flow rate influence is
of an image. The Sobel command is used in MATLAB as a more dominant than that of the liquid. Because the shear
function for image processing. The change in pixel gradients forces from the air flow are exerted on the liquid sheet in the
of an image is applied by the Sobel method for a complete nozzle head, where the mixture of air and liquid is spread out
separation of BW pixels before processing. In Fig. 8, from the nozzle orifice, it atomizes into small droplets.
different stages for preparing image (before processing) are
shown. The main image changed into BW and the Sobel
method is applied to separate the boundary between the dark
and light spots in an image for better processing and more
accurate results.
D10 D=
=
∑ 0Di (5)
mean
Ni
the droplet mean diameter has some fluctuations in the results show that as the spray cone section is set farther from
droplet size. the nozzle orifice, the droplet mean velocity decreases. The
liquid flow in the orifice centre is fast and near the orifice
Table 2 The mean diameter of the droplets in different distances from the nozzle tip walls, and it decreases due to the wall effects. As a result of
Distance (cm) Mean Diameter (μm) the mean droplet velocity, this velocity decreases. In the
8 25.02
14 24.5
spray cone centre line, the velocity reaches the maximum.
20 26.09 Since the spray jet is in the straight trajectory and away from
the spray cone centre, the spray stream angle causes the
5.2 The Results of the Mean Droplet’s Velocity droplets’ movement in more distances than the one of a
straight line and the mean droplet velocity decreases.
Due to the existence of air-cored vortex in the pressure
Table 3 The mean velocity of the droplets in different distances from the nozzle tip
swirl atomizer, the droplet number concentration is high at
Distance (cm) Mean Velocity (m/s)
the region near the axis of the spray. Fig. 13 shows the radial 8 14.9
distribution of the droplet mean velocity. At the centre of the 14 6.9
spray cone, the droplet velocity is high and reduces towards 20 5.05
the edge of the spray. A relatively uniform velocity
distribution of droplets in the spray is observed. At the axial 5.3 Calculation Number Density as Radial Distribution
distance of 8 cm, the calculated mean velocity of the droplets
at the centre is larger than the measured results at the edge of Fig. 14 shows the predicted number density of the
the spray. This is due to the over-estimation of the mean droplets at each cross-section. The number density
diameter which is depicted in Fig. 13. At the axial distance distributions are defined by dividing the number of the
of 20 cm, the calculated mean velocity of the droplets at the droplets in each interval by the total number of droplets
centre of the spray is smaller than the measured results. This accommodated on the circle radius of the cone section. The
might be because of the uncertainty associated with the number density is the percentage of the droplets’ number to
model implemented within the experimental results which the whole number, in different intervals. It is considerable
may have out layer points. The mean velocity of droplets at that the number of drops is more at the area near the spray
each interval in the spray cross-section is calculated by the centre, rather than other points, and is an indicator of the swirl
Eq. (6). motion of air-liquid mixture, which is counted as a feature of
the nozzle. Moreover, at intervals of 14 and 20 cm, the
N distance of droplets increases from the spray cone centre and
VAvg =
∑ i =1vi (6) the number of drops reduces at that area, which is why the
N number density reduces as well.
5.4 Air Pressure Changes and their Effect on the Pressure increasing the distance from the nozzle's head. The average
Swirl Nozzle velocity of drops is reduced gradually by increasing the air
pressure. By means of decreasing the droplet velocity, the
Fig. 15 shows the changes in the spray angle versus the number of collision is reduced and the average number of
air pressure for the pressure swirl nozzle system. The average drops from the centre to the edges of the spray cone increases
value of the spray angle for this nozzle is about 41 degrees, gently. Furthermore, near the edge of the spray cone, the
for Lacava et al., it is 27 degrees and for Digvijay et al., it is drops’ size varies as a result of the existing produced vortex
32 degrees. The percentages of the relative difference of the in the surrounding air of the jet spray, resulting in the drop
spray angle in this research in comparison with these two size distribution. The advantage of the image processing
works (Lacava and Digvijay) are 33% and 21%, respectively method is that it is cheap and available. The accuracy and
[15, 16]. quality of the images is related to the environmental
The reason behind this difference is the different conditions such as light and digital camera. The results
structure of nozzles, equipment type and experiment showed that by increasing the atomization air pressure, the
operational conditions. spray angle increases and it is predicted that the droplet size
decreases because of an increase in droplet velocity and more
collision of the droplets with each other and the orifice wall
at the outlet of the nozzle. It is found that the spray angle is a
function of air pressure and the orifice diameter of the
nozzle's head. Moreover, by increasing the distance from the
spray centre, the droplets’ average velocity decreases
because the outlet droplets’ contact the edges of the nozzle
orifice and air friction factor decrease their velocity.
Acknowledgement
Abstract: High speed steel (HSS) is a very important industrial tool material and has been constantly improved for different wear resistance applications and cutting tools, i.e.
drills, milling cutters, hobs and for the cutting tools in which the economical cutting speed is too low for choosing the carbide tools. The properties of HSS depend significantly on
the parameters of the conducted heat treatment. In this paper, the influence of deep cryogenic treatment in combination with nitriding of metallurgical powder metallurgy HSS on
the wear resistance was measured. Additionally, the cutting performance in a single point cutting tool machinability test at the configuration of the dry low-speed turning of steel
was investigated. The results showed that deep cryogenic treatment itself, and in combination with nitriding, resulted in the reduction of the wear rate. The results of the single
point cutting tool machinability test showed that deep cryogenic treated and nitrided HSS inserts performed worse than the classically heat-treated inserts and deep cryogenic
treated HSS inserts exhibited approximately the same flank wear as the nitrided ones.
Keywords: deep cryogenic treatment; cutting tools; high speed steel; nitriding; wear resistance
properties. The results showed that the most significant DCT, it could eliminate the advantageous influence of deep-
parameters in terms of mechanical properties are cryogenic treatment.
austenitizing and tempering temperatures. A low It is clear that an optimized combination of cryogenic
austenitizing temperature combined with a high tempering treatment temperature and time, cooling rate, and subsequent
temperature improves fracture toughness, bending strength tempering would make the best of sub-zero treatment while
and elongation, while the high austenitizing and low balancing the cost and performance of the material.
tempering temperatures improve hardness and wear The purpose of this paper is to study the influence of
behaviour. deep cryogenic treatment, in combination with nitriding, of
Authors in [14] analysed the influence of DCT with powder metallurgy high speed steel on the wear resistance
different austenitizing and tempering temperatures on and cutting performance in a single point cutting tool
fracture toughness, hardness and wear resistance of PM machinability test at the configuration of the dry low-speed
lower C, high W and Co cold work tool steel, high C and V turning of steel.
PM cold work tool steel and PM HSS. They reported that in
the case of low carbon cold-work tool steel, the DCT results 2 MATERIALS AND METHODS
greatly improved in fracture toughness, while maintaining
high hardness, which can be linked to the formation of finer The material used in this study was a high speed steel
needles such as martensite and martensitic transformation produced by powder metallurgy, grade PM S390 MC. The
accompanied by plastic deformation of the primary chemical composition of steel is presented in Tab. 1.
martensite. On the other hand, for high C cold-work tool
steel, the DCT has a negative effect, while for high-speed Table 1 Chemical composition of the PM S390 MC high speed steel
steel, the DCT has virtually no impact on its properties. They %C % Si % Mn % Cr % Mo % V % W % Co % Fe
1.64 0.60 0.30 4.80 2.00 4.80 10.40 8.00 rest
also reported an alteration in fracture toughness/Rockwell
hardness ratio which affects wear resistance of cold-work
For wear resistance-testing, disc shaped samples were
tool steel. It was concluded, however, that hardness is still a
prepared with the dimensions of ∅22×4 mm. For the purpose
major parameter affecting abrasive wear resistance.
of a single point cutting tool machinability test or turning test,
The application of DCT to HSS cutting tools has been
20 cutting tool inserts were prepared according to the ISO
enhanced in the recent years due to its effect in the prolonging
13399, CNMA 120408, Fig. 1, with the dimensions shown in
of tool life, dimension stability and stress relief. Authors in
Tab. 2.
[15] reported that the M2 HSS twist drills tool life increased
for approximately 280 % when drilling steel, taking the total
failure of tools as a criterion for determining the tool service
life. Authors in [16] reported a 77 % and 126 % improvement
tool life of DCT and tempered M2 tool steel drills.
Authors in [17] investigated the combined impact of the
cryogenic temperature with tempering on the mechanical
properties and wear resistance of the W9 HSS. They
Figure 1 Insert CNMA 120408 [18]
concluded that all varieties of DCT influenced the wear
resistance and hardness of W9 HSS. The performance of the Table 2 CNMA 120408 insert dimensions
W9 HSS taps during the tapping test was not notably Grade L, mm IC, mm S, mm RE, mm
influenced by the DCT temperature. It was shown that tool CNMA 120408 12.9 12.7 4.76 0.8
life was the longest when treated in the temperature range of
−196 to −160 °C with subsequent tempering, and it was The heat treatment was performed in an Ipsen vacuum
concluded that for treatment effectiveness and cost furnace with the preheating parameters 650 °C / 30 min, 850
efficiency, the deep cryogenic treatment at −160 °C could be °C / 20 min, 1050 °C / 15 min before the heating to the
the optimal choice. austenitization temperature. The heating rate to the
Authors [5] investigated the impact of different time and austenitization temperature of 1130 °C was 5 °C/min and the
temperature DCT parameters and the combination of DCT austenitization lasted 6 min, after what samples were
and plasma nitriding on the wear performance of powder– quenched in the nitrogen flow, p = 1050 mbar. The heat
metallurgy HSS considering the wear resistance and treatment parameters are presented in Tab. 3.
resistance to galling under dry sliding conditions. They
concluded that DCT improves the microstructure of PM HSS Table 3 Heat treatment parameters
Austenitization, Deep cryogenic Tempering, Nitriding,
generating the finer needle-like martensitic structure, which Batch
°C/min treatment, °C/h °C/h °C/h
results in higher hardness and better wear resistance 520 / 520 /
HT 1130 / 6 min -
properties, especially considering friction, and galling 490 / 2 h
resistance against stainless steel which considerably HTN 1130 / 6 min -
520 / 520 /
490 / 2.5 h
improves with longer cryogenic treatment times. They also 490 / 2 h
DCT 1130 / 6 min −196 °C / 24 h 520 / 2 h
concluded that plasma nitriding improves tribological DCTN 1130 / 6 min −196 °C / 24 h 520 / 2 h 490 / 2.5 h
properties of PM HSS and reduces the effect of austenitizing
temperature, but the results indicated that if combined with
After quenching, part of the test samples were soaked in treatment resulted in a 23.6 % reduction of the wear rate.
liquid nitrogen at a controlled speed and left for 24 h, Nitriding resulted in an approximately 18.4 % reduction of
followed by single tempering. Other samples were tempered the wear rate and in combination with DCT 25.6 %
three times. Part of the vacuum and deep cryogenic treated respectively.
samples were plasma nitrided (95% H2 and 5% N2). Nitriding
in the duration of 2.5 hours was carried out at a temperature
of 490 °C in a Metalplas Ionon BDAG GRUPPE BALCKE-
DÜRR AG apparatus. After heat treatment, thorough
analysis of the microstructure and mechanical properties was
performed on samples and the results were published in [19].
Wear resistance of specimens was determined by
simulating two-body abrasive wear by alternating the sliding
motion by using ball-on-flat contact with the Al2O3 ball as a
counter-body with a hardness of 1200 HV and a diameter of
∅10 mm. Wear tests were performed under dry sliding
conditions at an average sliding speed of 0.024 m/s with a
frequency of 5 Hz and amplitude of 2.4 mm, maximum
contact pressure of 1300 MPa (FN = 10 N) and total sliding
distance of 30 m, Fig. 2. Test results were evaluated in terms Figure 3 Flank wear of the tool, VB
of the samples’ wear rate (mm3/Nm) and average coefficient
of friction. Testing was performed at room temperature (22 The results show that nitriding after deep cryogenic
°C) and relative humidity of 40 %. treatment did not significantly influence the reduction in
wear rate as it influenced the wear rate after heat treatment.
Regarding the coefficient of friction, DCT samples had a bit
higher coefficient of friction, approximately 4 % higher than
classically heat-treated samples. Through nitrogen diffusion
layer formation [19, 20], plasma nitriding reduced friction
coefficient of the analysed steel in classic heat-treated
samples (HTN) and DCTN samples for 8 % and 11 %,
respectively.
Figure 2 Schematic diagram of ball-on-flat wear test
For every heat treatment batch, five HSS tool inserts Figure 4 Wear rate and Friction coefficient
were used with two cutting edges. Tool wear was measured
after every two minutes of turning by measuring the flank Fig. 5 presents the values of the average flank wear of
wear, VB, by using the Tescan Vega TS5136LS scanning the cutting tool inserts measured every two minutes of
electron microscope. The criterion for the end of tool life was turning. The results show that after two minutes, heat treated
the catastrophic failure of five 2-minute cycles of turning, i.e. (HT) inserts had the highest value of average flank wear, but
a total of 10 minutes. VB was the width of the flank wear after four minutes of turning, the cryogenically treated and
track at the edge of the flank, Fig. 3. The results represent the nitrided (DCTN) tool inserts had higher flank wear values
mean value of ten measurements. than the heat treated (HT) ones. The heat-treated inserts (HT)
and cryogenically treated and nitrided (DCTN) inserts had
3 RESULTS AND DISCUSSION total tool failure after eight minutes of turning with flank
wear values above 1.2 mm. The DCT inserts and heat-treated
The results of the abrasion wear testing show that DCT and nitrided (HTN) inserts exhibited similar wear behaviour
and nitriding resulted in lower wear rates of the high-speed with slightly lower flank wear values for the heat-treated and
steel specimens tested against alumina ball counter-body. nitrided (HTN) tool inserts. When compared with the
The results are presented in Figure 4. Deep cryogenic
laboratory wear test results, it can be observed that the DCT pronounced at a higher rate and they contributed to the
samples had a lower wear rate than the HTN samples and that intensity of the cutting tool insert flank face wear.
the DCTN samples exhibited the lowest wear rate of all Fig. 6 shows wear of the cutting inserts flank face and
samples, and in the machinability test, the DCTN tool inserts the formation of a build-up edge of the flank face in all four
had the poorest wear behaviour. Similar results were reported cutting inserts can be observed. The build-up edge is formed
in [5], where it was concluded by the authors that plasma as a result of indirect adhesion wear [21].
nitriding enhances the tribological properties of powder Adhesion wear is one of the tool wear mechanisms that
metallurgy HSS, but if combined with the DCT, it might can be present in a wider range of cutting temperatures. It can
reduce the advantageous effect of the deep-cryogenic be produced by two different ways. Direct adhesion wear is
treatment and the reasons for that are not clear at present. caused by the integration of tool particles to the chips, and
indirect or secondary adhesion wear is caused by an
integration of a fragment of the workpiece material to the
tool. This affects tool wear in the way that tool geometry is
altered by material integration and when these fragments are
removed, they can extract tool particles causing tool wear
[21]. As observed in Fig. 6 d, the highest build-up edge after
four minutes of turning is present at the DCTN insert even
though DCTN samples had the lowest coefficient of friction
in the laboratory wear test. This is also in accordance with
[5], where it was concluded after the laboratory wear tests
performed on samples, that plasma nitriding in combination
with DCT might eliminate the advantageous effect of DCT,
leading to reduced galling resistance. Galling resistance
reduction also had a large influence on the results in the
Figure 5 Average flank wear of differently heat treated HSS cutting tool inserts machinability test performed within this research since
galling is a very intense form of adhesion wear. Adhesion
wear leads to the formation of a build-up edge of the tool
flank face intensifying the wear of the cutting tool.
4 CONCLUSION
the analysed steel in the classically heat-treated samples and the properties of tool steel. J. Mater. Process. Technol., 118,
deep cryogenic treated samples for 8 % and 11 %, 350-355. https://doi.org/10.1016/S0924-0136(01)00973-6
respectively. Regarding the behaviour of HSS in the [11] Mohan Lal, D., Renganarayanan, S., Kalanidhi, A. (2001).
Cryogenic treatment to augment wear resistance of tool and die
machinability test, the heat-treated (HT) inserts and deep
steels. Cryogenics, 41, 149-155.
cryogenic treated and nitrided (DCTN) inserts exhibited a https://doi.org/10.1016/S0011-2275(01)00065-0
total tool failure after eight minutes of turning with flank [12] Baldissera, P. & Delprete, C. (2008). Deep Cryogenic
wear values above 1.2 mm. Treatment: A Bibliographic Review. The Open Mechanical
DCT inserts and heat-treated and nitrided (HTN) inserts Engineering Journal, 2, 1-11.
exhibited similar wear behaviour through the machinability https://doi.org/10.2174/1874155X00802010001
test with slightly lower flank wear values for the heat-treated [13] Oppenkowski, A., Weber, S., & Theisen, W. (2010).
and nitrided (HTN) tool inserts. When compared to the Evaluation of factors influencing deep cryogenic treatment that
laboratory wear test results, it can be observed that the DCT affect the properties of tool steels. J. Mater. Process. Technol.,
samples had a lower wear rate than the HTN samples and that 210, 1949-1955. https://doi.org/10.1016/j.jmatprotec.2010.07.007
[14] Podgornik, B., Paulin, I., Zajec, B., Jacobson, S., & Leskovšek,
the DCTN samples exhibited the lowest wear rate of all V. (2016). Deep cryogenic treatment of tool steels. Journal of
samples and in the machinability test, the DCTN tool inserts Materials Processing Technology, 229, 398-406.
had the poorest wear behaviour. https://doi.org/10.1016/j.jmatprotec.2015.09.045
Based on the results of the machinability test, it can be [15] Silva, F. J., Franco, S. D., Machado, A. R., Ezugwu, E. O.,
concluded that when high-speed steel is subjected to Souza, A. M. J. (2006). Performance of cryogenically treated
conditions in the real cutting process, classical heat treatment HSS tools. Wear, 261, 674-685.
and nitriding have the same effect as DCT and in https://doi.org/10.1016/j.wear.2006.01.017
combination with DCT, plasma nitriding eliminates the [16] Firouzdor, V., Nejati, E., & Khomamizadeh, F. (2008). Effect
of deep cryogenic treatment on wear resistance and tool life of
advantageous effect of the deep-cryogenic treatment.
M2 HSS drill. Journal of Materials Processing Technology,
206, 467-472. https://doi.org/10.1016/j.jmatprotec.2007.12.072
5 REFERENCES [17] Yan, X. G. & Li, D. Y. (2013). Effects of the sub-zero treatment
condition on microstructure, mechanical behavior and wear
[1] Holmberg, K. & Matthews, A. (1994). Coatings Tribology, resistance of W9Mo3Cr4V high speed steel. Wear, 302(1-2),
Elsevier Science B. V., Amsterdam. 854-862. https://doi.org/10.1016/j.wear.2012.12.037
[2] Soković, M., Barišić, B., & Sladić, S. (2009). Model of quality [18] www.iscar.com
management of hard coatings on ceramic cutting tools. Journal [19] Šolić, S., Cajner, F., & Leskovšek, V. (2012). Effect of Deep
of Materials Processing Technology, 209, 4207-4216. Cryogenic Treatment on Mechanical and Tribological
https://doi.org/10.1016/j.jmatprotec.2008.11.026 Properties of PM S390 MC High-Speed Steel. Materials
[3] Šolić, S. (2011). Influence of the microstructure on the Testing, 54(10), 688-693. https://doi.org/10.3139/120.110380
tribological properties of cutting tools, doctoral thesis, Zagreb, [20] Kato, H. & Eyre, T. S. (1994). Sliding wear characteristics of
UniZg, FAMENA. (in Croatian) nitrided steels. Surf. Eng., 10, 65-74.
[4] Šolić, S., Schauperl, Z., Godec, M. & Tropša, V. (2017). https://doi.org/10.1179/sur.1994.10.1.65
Microstructural changes in heat treatment of PM high-speed [21] Gómez-Parra, A., Álvarez-Alcón, M., Salguero, J., Batista, M.,
steels. Tehnički glasnik, 11(4), 166-170. & Marcos, M. (2013). Analysis of the evolution of the Built-
[5] Podgornik, B., Majdic, F., Leskovšek, V., & Vizintin, J. (2012). Up Edge and Built-Up Layer formation mechanisms in the dry
Improving tribological properties of tool steels through turning of aeronautical aluminium alloys. Wear, 302(1-2),
combination of deep-cryogenic treatment and plasma nitriding. 1209-1218. https://doi.org/10.1016/j.wear.2012.12.001
Wear, 288, 88-93. https://doi.org/10.1016/j.wear.2011.04.001
[6] Meng, F., Tagashira, K., Azuma, R., & Sohma, H. (1994). Role
of eta-carbide precipitations in the wear resistance Authors’ contacts:
improvements of Fe–12Cr–Mo–V–1.4C tool steel by
cryogenic treatment. ISIJ International, 34, 205-210. Dr. sc. Sanja ŠOLĆ, Assistant Professor
https://doi.org/10.2355/isijinternational.34.205 University North, Department of Mechanical Engineering
[7] Meng, F., Tagashira, K., Azuma, R., & Sohma, H. (1994). Role J. Križanića 31b, 42 000 Varaždin, Croatia
of eta-carbide precipitation’s in the wear resistance E-mail: sanja.solic@unin.hr
improvements Fe-12Cr-Mo-V-1,4C tool steel by cryogenic
Dr. sc. Zdravko SCHAUPERL, Full Professor
treatment. ISIJ International, 34(2), 205-210.
University of Zagreb, Faculty for Mechanical Engineering and Naval Architecture,
https://doi.org/10.2355/isijinternational.34.205 Ivana Lučića 5, 10000 Zagreb, Croatia
[8] Leskovšek, V., Kalin, M., & Vižintin, J. (2006). Influence of E-mail: zdravko.schauperl@fsb.hr
deep-cryogenic treatment on wear resistance of vacuum heat-
treated HSS. Vacuum, 80, 507-518. Dr. sc. Vlado TROPŠA, Associate Professor
https://doi.org/10.1016/j.vacuum.2005.08.023 University North, Department for Mechanical Engineering
[9] Leskovšek, V. & Ule, B. (1998). Improved vacuum heat- J. Križanića 31b, 42000 Varaždin, Croatia
treatment for fine-blanking tools from high speed steel M2. E-mail: vlado.tropsa@unin.hr
Journal of Material Processing Technology, 82, 89-94.
https://doi.org/10.1016/S0924-0136(98)00023-5
[10] Molinari, A., Pellizzari, M., Gialanella, S., Straffelini, G., &
Stiasny, K. H. (2001). Effect of deep cryogenic treatment on
Abstract: The food industry has recognized the important role that multi-layered, flexible packaging materials play and it uses them extensively within the packaging processes.
Trends within food processing and packaging continuously encourage manufacturers to develop new technologies which extend a product’s shelf-life. The samples that were the
subject of research described in this paper are retort stand-up pouches made of multi-layered flexible material, the layers of which have been bonded with a laminating process
that aims to both extend a product’s shelf-life and stability, permitting the final processing of the finished product under aggressive regimes. It is intended that packaging, when
subjected to graphic reproduction and high pressure processing, shall remain unchanged in all aspects, in particular the visual aspect that is the subject of the research. The
samples were printed on polyester film (PET) using a rotogravure printing technique and laminated on aluminum (Al), oriented polyamide (OPA) and finally on inner polypropylene
(PP) layers. The sample was monitored for the Vegeta Blue spot colorimetric difference in the CIE L*a*b* space throughout the individual graphic reproduction phases and the
final high pressure processing (HPP) of the finished product. The goal is to determine the size of the spot color colorimetric difference (∆𝐸𝐸 ∗ ) in relation to the defined formula and
the established standard. On the basis of measured results, guidelines and recommendations for the correction and quality assessment of spot color reproduction throughout the
entire process were defined, with the goal to minimize any deviation to the lowest possible level.
Keywords: CIE L*a*b*; colorimetric difference; food packaging; high pressure processing; spot color
treatment, and maintains its original freshness throughout the division as the stiffness of a material depends on the type of
shelf life [5]. polymer. Therefore, for example, PVC (polyvinylchloride)
or cellulose nitrate will already be rigid at 70 – 80 μ, whereas
LD PE (low density polyethylene) will not be rigid even at
Polyester (PET) 300 μ [10].
Printing Ink
Laminate adhesive Polymers are very rarely used in their original form, and
Aluminium (Al)
Laminate adhesive
in order to achieve good characteristics to be a packaging
Oriented polyamide (OPA) material, different additives are added during their
Laminate adhesive
Polypropylene (PP) production to improve their properties. In addition to the
combination of different polymers, multi-layered flexible
packaging is combined with other materials (e.g. paper and
Figure 1 Multi-layered Stand-up Pouch Material aluminum) that aim to increase endurance, packaging
resistance and barrier properties.
2 FLEXIBLE MULTI-LAYERED FOOD PACKAGING Currently, a new EU Directive is being developed, and it
is known as the SUP Directive (Single-Use Plastics). This
Food packaging experienced its greatest expansion and Directive shall also partly cover the packaging materials
flourishing in the second half of the 20th century. This is primarily intended for snack products, beverage industry and
partly due to the development and discovery of new polymer in general all food intended for direct consumption without
materials and partly due to increased food production as a further processing [11]. A great deal of pressure has been
result of population growth, urbanization, lifestyle changes, applied to manufacturers of plastics and polymer materials
female employment and increasingly pronounced hygiene relating to the adequate disposal and recycling of materials
habits [6]. with a view to reducing marine pollution in seas, marinas and
Packaging material is any material from which the environment in general. Regardless of its negative
packaging is made, such as: glass, plastic, paper, cardboard, reputation and poor image, plastic has a number of benefits
wood, metal multi-layered mixed materials and other over other packaging materials and it is impossible to
materials [7]. imagine life without it.
Multi-layered mixed flexible material plays an important Only 1.5% of all oil and gas consumed in Europe is used
role within the packaging material segment. Multi-layered as a raw material for the production of plastic packaging,
packaging is predicted to grow at an annual rate of 4.1%, while 90% is used for heating, transportation and energy
meaning an increase from 27.4 million metric tons in 2017 to production. If food is packaged with other packaging
33.5 million metric tons in 2022. The financial market share materials, the associated energy consumption would be
of flexible packaging material is growing annually by 4.3% doubled and greenhouse gas emissions would be almost
and by 2022 this is expected to amount to 283 billion USD tripled. This would also be accompanied by a 360% increase
[8]. in packaging weight. It is the lightest packaging material.
Historically, polymer materials have only been used as a More than 50% of all European products are packed in plastic
packaging material for a short period of time and, from their that only accounts for 17% of the total packaging weight of
initial mass use in the 1950s onto the present day, their all products on the market. Furthermore, this weight has been
properties have constantly improved, changed and been reduced by 28% over the previous 10 years. Lightweight
enhanced. Packaging materials are almost exclusively made packaging means lighter loads and fewer trucks required to
from plastomers. A considerably smaller amount of duromers deliver the same amount of products, thus reducing energy
and elastomers are used as an auxiliary material for sealing consumption for transportation, CO2 emissions and
lids and fasteners [9]. transportation costs [12].
The polymers most commonly used in plastic packaging From the perspective of companies and consumers
include: (users), packaging must achieve several goals:
• Polyethylene (PE), • Identify a brand,
• Polypropylene (PP), • Contain the product’s description and valid information,
• Styrene homo- and copolymers [polystyrene, (PS), • Facilitate smooth transportation and product protection,
styrene-acrylonitrile (SAN) and acrylonitrile butadiene • Facilitate storage for the consumer and
styrene (ABS), • Enable easy consumption of the product [13].
• Polyvinylchloride (PVC)
• Polyamide (PA) All of this is possible with multi-layer flexible packaging
• Polyethylene terephtalate (PET) and as certain polymer materials provide an excellent printing
• Polycarbonate (PC) [9]. substrate, look extremely attractive on the shelf and offer a
range of options for additional effects, such as various
Usually, these are divided into flexible and rigid lacquers, partial matte or glossy finishing, a metallic effect if
packaging materials. Everything that has a thickness of up to a printed surface is laminated on an aluminum surface
250 μ (0.010 inch) is considered a flexible material or film, without white primer, etc. Additionally, handling, production
and everything thicker is considered a rigid polymer. and transportation of this type of packaging material is
However, in practice, this is not the most effective method of simple and allows for storage in ambient temperatures for an
extended shelf life, provided of course that the packaging complex and requires the connectivity of designers,
remains hermetically sealed and secured. A relatively new marketing experts, development engineers and specialist
form of packaging on the Croatian market are stand-up expertise from suppliers of certain packaging materials.
pouches made of flexible multi-layered materials, suitable for Primary packaging is in most instances characterized by
sterilization and high pressure processing and practical for applied printing. In order to achieve excellent printing
use and application. Stand-up pouches can take various results, graphical preparation, customized to specific printing
forms, depending on the molding tool on the forming techniques and characteristics, is required. The printing
machine and the staples used, and they can include possible phase is followed by lamination processes, die-cutting
variants such as a sealing strip or laser perforations or forming and sealing the retort stand-up pouch in this case. In
indentations for easier opening. Depending on the addition to the processes listed above, the packaging that is
combination of packaging materials, these may be the subject of this research is also exposed to the processing
transparent or contain an aluminum layer, but the latter are of a finished product, which in this instance means exposure
not suitable for heating in a microwave oven. to high pressure processing during which unwanted changes
The retort stand-up pouches are the actual subject matter to the packaging material can occur, and that manifest in the
of this research. It defines the colorimetric difference of the form of de-lamination of the packaging material, separation
Vegeta Blue spot color in the CIE L*a*b* space applied to the of individual layers or changes in a polymer, which in turn
packaging design to detect the phase of the process that can affect colorimetric difference and subsequently
causes the greatest colorimetric difference in relation to the consumer perception.
defined standard and color formula. When shopping,
customers are inclined to select certain products according to 3.1 Design, Graphic Preparation and Printed Form
their individual habits and memorized characteristics. The Production
main feature of Vegeta branded products is their recognizable
blue color. Vegeta Blue is a color that should be kept within For the purposes of the research presented in this paper,
a strictly defined tolerance range (∆𝐸𝐸 ∗ ≤ 1), regardless of the a design has been applied to the already existing technical
type of packaging material, printing technique or finishing drawing of the retort stand-up pouch produced by the
process). This is not a simple task at all, as it requires Company Aluflexpack Novi from Umag. Packaging design
selecting an adequate supplier that understands the specifics usually begins with a good, clear and detailed brief for the
of the printing technique, its alignment within the graphic agency and the designer. Since it was necessary to apply the
preparation and design process and permanent quality control Vegeta Blue spot color on as large area of the pouch that was
and measurement of colorimetric values within the CIE the subject of this research as possible, the design was
L*a*b* space throughout the entire graphic and technological completed in the Podravka Design Studio and it included the
process. tailoring of an already existing design to new dimensions.
The Brief is usually drafted by the marketing team with
3 RESEARCH METHODOLOGY expert support from development specialists, who mostly
reduce the ideas of designers and marketing into a practicable
This chapter describes the methodology used in this and feasible framework.
research and the impact of the individual phases of the The stand-up pouch constructed from a multi-layered
graphic and technological process on the Vegeta Blue spot packaging material onto which the design was applied has a
colorimetric difference. Standard values for the Vegeta Blue total volume of 500 ml and measures 140×240 mm. The
spot color were prescribed in Podravka Packaging format of the packaging that is the subject of research is part
Development Dpt., coordinated with the supplier. The results of the standard Aluflexpack Novi product range and
of the measurement of the CIE L*a*b* color values and their formatted components of this shape are currently included in
colorimetric difference per processing stages were compared production. In addition to the design, a color measuring atlas
and colorimetric differences were observed compared to the has been applied on the back of the bag, which will be used
defined standard and the set Vegeta Blue spot color formula. in future studies relating to process of colorimetric
The measurement was taken after printing, lamination and difference.
high pressure processing. It indicated the stage at which the The design and graphical preparations were completed
greatest colorimetric difference outside the tolerance area in Adobe Illustrator CC. Color separation was completed in
occurred and whether it could have an impact on consumer Esko ArtPro, Version 16.0.0 in Company Anilox. Separation
perception. for printing was performed by using seven colors (cyan,
The graphical and technological process is divided into magenta, yellow, black, Vegeta Blue, reflex blue and white).
several stages. The first of these is the graphical design stage, A test print was performed on an Epson Stylus Pro WT7900
which in some instances implies solely the application of ink jet printer, which uses a Micro Piezo TFP head for
visual elements to the already existing form and packaging printing. The printer color management software is GMG
material, and in certain situations this may mean designing a RIP and GMG FlexoProof, Ver.5.8.1.312. A test print using
new form mainly in rigid plastic, glass, or metal packaging, digital technology (GMG proof) was performed for quality
which also requires the development of a new label, a shrink control of all design elements. Podravka, which in this
sleeve foil or another final package decoration. In such instance commissioned the design, received the final color
situations, this design process is more time-consuming and separation and digital print, for quality control purposes. A
signed digital print serves as a guarantee for the packaging robotic arms and are fully automated, replacing human labor
supplier that all is in order and that the next stage of the and they have led to significantly shorter delivery times,
process, namely the engraving of printed forms, can be often no longer than a few days [16].
embarked upon.
Separation was engraved from the application Collage,
Printed image
which was also used for installation and to define the
following parameters: Impression cylinder
• Lineature, (70 L/cm – 101 L/cm, depends on color),
• Individual separation angles, (33° – 60°, depends on the
color),
• Tool angle, (110° – 130°, depends on the color), Substrate
• Reproduction curve. Plate
Doctor blade
Preliminary prepared cylinders were then engraved on Plate
electronic engraving machines. The engraving was cylinder
performed on an Ohio Gravure Tec. Inc. USA, Gravostar Printing ink
Spectrum Engraving Machine. The process was tailored to
the parameters of the selected substrate and machine on Figure 2 Rotogravure Printing Schematic [17]
which the printing was performed. In order for the print
cylinders to be sufficiently durable to withstand printing on Aside from the high prices of the printing forms,
printing machines, they were finally processed with hard rotogravure printing is also known for evaporation ink
chrome plating, where a protective layer of chromium, finishing. Most printing that utilizes this technique involves
several micro millimeters thick, was applied. Final quality inks based on organic solvents. The vapors from organic
control was performed by a visual inspection of the finished solvents are harmful to human health and highly flammable.
printed forms and the pattern was printed on a test print For these reasons, the plants and facilities where rotogravure
machine utilizing the original inks and consumer substrate printing machines are used require special construction and
(PET 12 μ) at the Aluflexpack Novi Print Shop by simulating safety measures [14]. Companies operating within the EU are
the parameters of the printing machine. A test print was subject to stringent legislation and must be equipped with
performed on a manually processed printing press. The test special types of air purification and solvent recovery devices.
print, when obtained as described above, showed if there The Aluflexpack Novi Print Shop in Umag printed
were any visible defects and if any corrections on print forms design samples on the Schiavi Pulsar Printing Machine, 10
were necessary. The test print with engraved print forms was group, 1.240 mm wide with a maximum speed of 350 m/min,
performed on Rotogravure Proofing Press, a test print utilizing a rotogravure printing technique on transparent, 12
machine produced by JM Heaford of the UK, on a HA μ thick polyester (PET). Aluflexpack Novi utilizes a cutting
3500×1800 (Fixed Drum / Hydraulic Rotation) Model. edge solvent recovery plant and has been awarded the BRC
(the British Retail Consortium) and ISO (the International
3.2 Printing and Finishing Organisation for Standardisation) production certification.
Siegwerk printing inks (thermal-resistant ink),
The samples that are the subject of this paper have been specifically intended for printing packaging that can later be
printed by using the rotogravure printing technique. processed by various methods applied to finished products,
Rotogravure printing is a direct, intaglio printing technique. were used. Resistances of this ink series are archived by
The print face is hollowed out in its printed form and the using completely resistant pigments according to ISO 2836.
blank spaces are foremost [14]. The printing ink is applied on The printed samples were measured for the values of Vegeta
the entire print form during printing and later, a doctor blade Blue and were compared with the set standard values and pre-
removes the printing ink from the blank spaces (Fig. 2). By defined color formula. After the printing phase was
pressing the printing form onto the substrate, the printing ink completed, the laminating process was performed on the
is transferred to the substrate. The printing surface of the PET Rotomec Rotoconvert Machine with three unwinders (a
foil should be treated to be suitable for printing, which has a triplex in a single pass) 1.230 mm wide, with a maximum
significant influence on surface tension and enables the ink speed of 300 m/min. A Henkel agent and a catalyst were used
to connect with the printing substrate [15]. for lamination. A flexible retort stand-up pouch consists of
This printing technique allows the reproduction of very four different materials: PET, Al, OPA and PP. After
small details, small font-sized text, excellent photographs lamination, each material has to pass though the cross-linking
and product views. The usual perception of rotogravure of materials and solvent evaporation phase which aim to
printing within the printing industry is that it is used in the eliminate the problem of residual odor within the packaging
printing of top class brand packaging which must be high- material and any de-lamination and migration problems. The
quality. Quality remains the main feature of the printing Totani BH-60DLLS machine was used to form, seal on three
technology, but in recent years, the producers of printing sides (laterally and on the bottom) and die-cut the stand-up
forms, i.e. engraving companies, have sought to reduce costs pouches, operating at a maximum speed of 1.000
and automate the process. New, fast engraving heads use pouches/min.
3.3 High Pressure Processing (HPP) defined at what stage in the process and under which
processing regimes, the greatest colorimetric difference
The processes that may cause unnecessary changes in the occurred and what the most significant differences within the
packaging’s properties and appearance do not always end CIE L*a*b* color space were. The average colorimetric
upon the completion of the graphical and technological difference for each phase was defined in the statistical
process. The packaging, as an integral part of the product, is analysis of the measured values and it was expressed by
subject to various food processing regimes, the aggressive applying the method of the median mean. Parallel curves
parameters which can affect the packaging material’s were also shown after each measurement stage
properties. demonstrating the colorimetric difference values compared
The application of food high pressure processing in the to the set color (Figs. 4 and 5).
food industry has developed as an alternative to conventional The measurement was performed on three standard light
thermal processes, such as pasteurization and sterilization, sources: D65./10°, A.../10° and F11./10°. Aluflexpack
aiming to achieve a microbiologically safe product that does Novi’s Print Shop had previously recorded the CIE L*a*b*
not undergo any sensory, physical or nutritional value Vegeta Blue spot color values that are the subject of the
changes. research and the defined formula. For the purposes of this
High Pressure Processing is a cold pasteurization process research, the color formula was to a lesser extent modified in
in which the final packaged product is exposed in a given its composition and it did not use the original values.
chamber to high isostatic pressure (300–600 MPa / 43.500– The following standard measuring values were set for
87.000 psi), which is applied by water (Fig. 3) [18]. Vegeta Blue under three different light sources:
The presented Figs. 4 and 5 show the spectral curves of packaging manufacturers and their specific features (printing
the measured Vegeta Blue spot color regarding the performed inks, machine type, substrate, lamination agents, varnishing
stage of the process. The right part of the picture shows the L etc.), the challenge of maintaining color uniformity within
diagram with two wavelengths of color in two phases and the the CIE L*a*b* space is sometimes insurmountable. The
ab diagram shows the color difference between the two colorimetric values on a multi-layer flexible packaging
phases. The colorimetric difference after printing is within containing aluminum and on the packaging containing only
the tolerance limits and amounts to D65./10° light source transparent polymer materials are not identical. To achieve
CMC2:1 = 0.23. After lamination, a greater colorimetric the identical print, which would mean that the spot color is
difference occurred, the color surpassed the determined within the tolerance area, it is necessary to redefine the spot
tolerance range and the colorimetric difference under the color input formula for every substrate. This, of course,
same conditions equaled to D65./10° CMC2:1 = 1.56 (Fig. makes the process more complicated and more expensive for
4). After high pressure processing, an even higher print shops, but for the manufacturers of brands that are
colorimetric difference of Vegeta Blue occurred and the recognizable by their specific, protected color, the goals is to
average value on the measured samples equaled to D65./10° maintain identical colorimetric values. Specific features of
CMC2:1= 1.94. With lamination, Vegeta Blue became darker individual packaging materials require that the maximum
and crossed into the area of yellow-red. This clearly shows allowed tolerance area and the standard colorimetric values
that HPP also influences the Vegeta Blue spot color for the spot color be defined. During the measuring and
deviations. evaluation of samples, it was noticed that a rather small
deviation of the Vegeta Blue spot color, which occurred after
lamination (ΔE = 1.56), was more visible than greater
Printing
deviations of some other colors. After high pressure
Lamination processing, according to the measured values, Vegeta Blue
1
was darker (Fig. 5).
Օ D65./10Օ 1.56
X A.../10Օ 1.71
ΔF11./10Օ 1,43 5 CONCLUSION
Increasingly aware and informed consumers now have involvement of several different teams and experts is also
numerous requirements which they expect from packaging, necessary in order to prevent any problems on production
and these requirements have to be met. Packaging no longer lines or later during the handling and processing of a product.
has only a protective or storage function. It is a ‘silent seller’ The fact that industrial production is not the same as the small
on the shelf, contributing to distinction, increasing manufacturing plants of family proprietorships should not be
recognition and is interactive, easy to handle and provides a forgotten. Industry seeks to eliminate manual work,
range of additional benefits such as laser perforation, re- machines have high capacities and under ideal conditions, the
sealable options, microwave heating or direct consumption. output control of the packaging manufacturer becomes the
It must comply with the guidelines on sustainable input control of the ordering party, but all this presumes
development, environmental protection, etc. At first sight, error-free packaging items. The shelves in shops today are
what the consumer ranks the highest is the color perception overcrowded, full of special offer packaging and novelties. It
and an emotional link with the brand. In order to achieve is impossible to determine which of the packaging properties
attractiveness in a shop and draw the initial attention of the is the most important for the consumer, whether the
consumer, which only lasts a few seconds, a consumer makes his/her decision consciously or based on
multidisciplinary approach and a highly complex emotions linked to a familiar brand. The task for every
implementation process are necessary. manufacturer of recognized trademarked products is to
The goal of the research presented in this paper was to achieve uniformity of spot color on different packaging
determine the value of the Vegeta Blue spot colorimetric materials, with different printing techniques and with
difference that appears on packaging during printing, different suppliers. The parameters which influence changes
processing and high pressure processing of the completed in spot colors are the subject of this paper and they are really
product. The Vegeta Blue spot color was measured after these numerous in products marketed in several different types of
stages and compared with the set standard in order to achieve packaging material. There is a difference between printing on
better qualitative properties for the final product and improve a mono material with a white background and a multi-layer
its attractiveness on the shelf. Based on the measurement flexible material composed with aluminum or metallized foil.
results, the process stages that provide the highest impact on For the Vegeta Blue spot color, just like for all other world
color deviation were determined. The formula for the Vegeta famous brands, the intended goal is to reach colorimetric
Blue spot color with its allowed tolerance area (ΔE ≤ 1) was values that will be within the desired tolerance area. There
defined as the measurement starting point. are several software solutions which will facilitate the
The results were the measured values for the Vegeta Blue standardization process, but every ordering party should be
spot color. High pressure processing (HPP) causes deviations aware that it is impossible to achieve laboratory conditions in
in the CIE L*a*b* values of spot color and it increases the production plant. It is easy to manage and control
deviation that is a result of the lamination process. The standardization with one supplier, on a defined type of
measurement has shown that the spot colorimetric difference packaging material, with a defined type of printing inks, the
that occurs during the processing phase is lower than the lamination agent that is always the same, printing forms
difference that occurs during the lamination stage. During made with identical parameters and always by the same
high pressure processing, a change in process colors was also manufacturer. All this cannot guarantee a perfect result, but
identified (cyan, magenta, yellow and black), in the colors of any colorimetric difference in such cases has been reduced to
the dish presented as a part of the design. The changes to a minimum. The reality of food companies is unfortunately
process colors were not the subject of the research presented slightly different from the described ideal conditions. All
in this paper, but the Company Anilox measured these values manufacturers, for a variety of reasons, use the services of
against the atlas of the colors applied, specifically to enable several packaging manufacturers and in such situations,
further research. It is exactly the process colors for which standardization, and subsequently the control of colorimetric
deviations resulting from sterilization and high pressure values and other packaging characteristics, become more
processing can be prevented during the graphic preparation complex and difficult in practice.
and print form design stages. The relating research provides a good basis for further
This Paper has defined the Vegeta Blue spot color value research in the area of process color and other packaging
of the colorimetric difference in the CIE L*a*b* color space material deviations, within the area of printing techniques
for printing, lamination and HPP stages. This paper also and processing procedures and other processing methods
proposes guidelines to achieve colorimetric difference within applied to a finished food product.
the previously defined tolerance ranges for individual colors,
and determines those parts of the technological process that 6 REFERENCES
have the greatest effect on color deviations.
Every project relating to the inception and creation of [1] Kotler, P. & Armstrong, G. (2012). Principles of Marketing.
packaging is unique and the same rules do not apply to all 14th Ed. Upper Saddle River, New Jersey: Pearson Prentice
packaging materials. Complexity and a multi-disciplinary Hall.
[2] Govers, P. C. M. & Schoormans, J. P. L. (2005). Product
approach add creativity and challenges to each new project.
personality and its influence on consumer preference. Journal
All stages require excellent correlation and communication of Consumer Marketing, 22(4), 189-197.
among the designer, the ordering party and the packaging https://doi.org/10.1108/07363760510605308
manufacturer. Within the commissioning company itself, the
[3] Young, T. K., Byungjin, M., & Kyung, W. K. (2014) General Authors' contacts:
Characteristics of Packaging Materials for Food System. In
Jung, H.H. (Eds.) Innovations in Food Packaging, 13-20. Renata TOMERLIN, MSc
Podravka d.d. (plc),
Elsevier Ltd. https://doi.org/10.1016/B978-0-12-394601-0.00002-3
Ante Starčevića 32,
[4] Welti-Chanes, J., López-Malo, A., Palou, E., et al. (2005). 48000 Koprivnica, Croatia
Fundamentals and Applications of High Pressure Processing of E-mail: renata.tomerlin@podravka.hr
Foods. Novel Food Processing Technologies. New York: CRC
Press. Mario TOMIŠA, PhD, Assoc. Prof.
[5] HPP Technology: What is High Pressure? (5th February 2018). University North,
Retrived from https://www.hiperbaric.com/en/ high- Trg dr. Žarka Dolinara 1,
48000 Koprivnica, Croatia
pressure E-mail: mario.tomisa@unin.hr
[6] Galić, K., Ciković, N., & Berković, K. (2000). Analiza
ambalažnog materijala. Zagreb, Hrvatska: Hinus d.o.o. Damir VUSIĆ, PhD, Assoc. Prof.
[7] Značenje pojmova iz Pravilnika o ambalaži, Savjetodavna University North,
služba HOK-a, (5th May 2018). Retrived from Jurja Križanića 31b,
http://infos.hok.hr/faq/f_tehnicka_pitanja/f9_zastita_okoline/z 42000 Varaždin, Croatia
nacenje_pojmova_iz_pravilnika_o_ambalazi E-mail: damir.vusic@unin.hr
[8] The Future of Packaging: Long Term Strategic Forecasts to
2028, Smithers Pira. (2018, October 5th). Retrived from
https://www.smitherspira.com/industry-market-reports/
packaging/packaging-long-term-strategic-forecasts-to-2028
[9] Vujković, I. A. (1997). Polimerna i kombinovana ambalaža.
Novi Sad, Srbija: Poli.
[10] Briston, J. H. (1990). Plastics Film. 3rd Ed. Essex, United
Kingdom: Longman Scientific & Technical.
[11] Proposal for a Directive of the European Parliament and of the
Council, (28th May 2018). Retrived from
http://ec.europa.eu/environment/circular-economy/pdf/single-
use_plastics_proposal.pdf
[12] Plastics - the Facts 2016. (18th Dec. 2017). Retrived from
https://www.plasticseurope.org/en/resources/publications/3-
plastics-facts-2016
[13] Bassin, S. B. (1988). Value-Added Packaging Cuts through
Store Clutter, Marketing News, September 26th, pp. 21.
[14] Bolanča, S. (2013). Tisak ambalaže. Zagreb, Hrvatska:
Sveučilišna naklada.
[15] Valdec, D., Miljković, P., & Auguštin, B. (2017). The
influence of printing substrate properties on color
characterization in flexography according to the ISO
specifications, Tehnički glasnik, 11(3), 73-77.
[16] Anyadike, N. (2003). Introduction to Flexible Packaging.
Surrey, UK: Pira International Ltd.
[17] Oldenzijl, R., Gaitens, G. and Dixon, D. (2010). Conduct Radio
Frequencies with Inks, Radio Frequency Identification
Fundamentals and Applications Design Methods and
Solutions, Cristina Turcu (Ed.), ISBN: 978-953-7619-72-5,
InTech, Retrived from: http://www.intechopen.com/books/
radio-frequencyidentification-fundamentals-and-applications-
design-methods-and-solutions/conduct-radio-frequencies-
withinks. https://doi.org/10.5772/7982
[18] Food Safety, High Pressure Processing of Foods, Authority of
Ireland, (September 2015). Retrived from https://www.fsai.ie/
publications_high_pressure_processing/
[19] Kuehni, R. G. (2003). Color space and its divisions: color
order from antiquity to the present. New Jersey: John Wiley &
Sons. https://doi.org/10.1002/0471432261
[20] A Guide to Understanding Color, x-Rite Pantone. (2016).
Retrived from https://www.xrite.com/-/media/xrite/files/
whitepaper_pdfs/l10-001_a_guide_to_understanding_color_
communication/l10-001_understand_color_en.pdf
Gordana JAUKOVIĆ
Abstract: This paper presents advantages and disadvantages of Raman spectroscopy when used for testing banknotes, i.e. in case study I, 1 and 5 dinars banknotes from a
private collection dating from 1876 of numismatic value were analyzed, and in case study II, 1000 dinars banknote from circulation together with its counterfeit. In both cases
banknotes were analyzed by using Raman spectrometer DXR Raman microscope at the Faculty of Physical Chemistry, with a laser of wave-length 780 nm at 12 mW of power.
Specters were recorded under microscope using 50X lens in a range 1800-60 cm-1. Analyses from case study I on both specimens showed spectrum characteristic for fluorescence
appearance, therefore intensive signal of fluorescence covered Raman signal. In case study II, it was established that Raman spectroscopy can reliably detect a counterfeit
banknote. Research, comparisons and analyses regarding original and counterfeit 1000 dinars note were performed on red banknote surfaces.
Keywords: banknotes counterfeiting; forensics; numismatics; printing inks; protective elements; Raman spectroscopy; spectrum
spectrum c) of counterfeit banknote most intensive stripe In region 3 of original banknote (Fig. 6) spectrum that
noticeable is at 1085 cm−1 arising from white pigment – chalk would enable identification of pigments was not received,
[11]. while counterfeit banknote produced very good spectrum of
During recording of spectrum in region 2 of original blue phthalocyanic pigment.
banknote (Fig. 5) intensive fluorescence occurred, so just a In region 4 of original banknote (Fig. 7) TiO2 was
single spectrum was recorded (after correction) with identified in anatase form, spectrum a), naphthol red pigment
noticeable weak stripes at 1586, 1560 and 1364 cm−1 of red (alizarin), spectrum b), while spectrum c) corresponds to blue
pigment (alizarin) [8]. In the same region of counterfeit phthalocyanic pigment. From counterfeit banknote received
banknote green phthalocyanic pigment, spectrum a), was was spectrum of green phthalocyanic pigment, with
identified as well as chalk, spectrum b). pronounced stripes of cellulose at 1120 and 1090 cm−1.
3 CONCLUSION
Acknowledgement
4 REFERENCES
[1] del Hoyo-Melendez, J. M., Gondko, K., Mendys, A., Krol, M.,
Klisinska-Kopacz, A., Sobczyka, J., & Jaworucka-Drath, A.
(2016). A multi-technique approach for detecting and
evaluating material inconsistencies in historical banknotes.
Forensic Science International, 266, 329-337.
https://doi.org/10.1016/j.forsciint.2016.06.018
[2] Badovinac, I. J., Orlic, N., Lofrumento, C., Dobrinic, J., &
Orlic, M. (2010). Spectral analysis of postage stamps and
banknotes from the region of Rijeka in Croatia. Nuc. Instrum.
Meth. A, 619, 487-490. https://doi.org/10.1016/j.nima.2009.10.174
[3] Skenderović-Božičević, M., Gajović, A., & Zjakić, I. (2012).
Identifying a common origin of toner printed counterfeit
banknotes by micro-Raman spectroscopy. Forensic Science
International, 223, 314-320.
https://doi.org/10.1016/j.forsciint.2012.10.007
[4] Itrić, C. & Modrić, D. (2017). Banknote characterization using
the FTIR spectroscopy. Tehnički glasnik, 11(3), 83-88.
[5] Bruna, A., Farinella, G. M., Guarnera, G. C., & Battiato, S.
Figure 7 Raman spectra of region 4 in original and counterfeit banknote (2013). Forgery detection and value identification of euro
banknotes. Sensors, 13(2), 2515-2529.
https://doi.org/10.3390/s130202515
Author’s contact:
Abstract: In the turning process of hard materials with CBN cutting inserts, metalworking fluid is not advisable as it results in poorer quality of finished surface and shortening of
the tool life. Dry machining, otherwise, develops excessive heat, residual stresses in the workpiece material due to thermal stress, and therefore causes possible problems of
achieving tolerances on larger workpieces. During the dry turning of induction hardened raceway on the wind turbine-bearing ring of 1500 mm in diameter, problems with achieving
tolerances occurred. The vortex-cooled air implementation during turning resulted in achieving designing tolerances of the workpiece, with better surface roughness and with less
tool wear. The result was 50% lower cost of cutting inserts and small operating costs associated with air consumption. The vortex-cooled air implementation makes it possible to
extend the tool life of CBN cutting inserts and thus higher cost-effectiveness of machining.
Air vortex is moving toward an adjustable valve at the hot you control the "cold fraction", which is the percentage of the
end that controls the volume of the airflow and the total input of compressed air that exits the cold end of the
temperature existing at the cold end. By adjusting the valve, Vortex Tube.
The adjustable valve leaks smaller part of the swirling air Machining of material 42CrMo4 was experimentally
along the wall of the tube (20% to 40%), and central part of investigated by Sutter et al [13]. They determined that at
air vortex is directed in the opposite direction creating inner cutting speed around 20 m/s, the temperature measured near
vortex through outer air vortex. Inner vortex transfers heat to the tool-chip interface achieved a maximal value of 870 °C
outer vortex near the wall and, with a significant decrease of for 42CrMo4. The increase of the cutting speed from 10 m/s
temperature, air exits on the cold end of the tube. Outer to 65 m/s raises continuously the chips’ temperature and
vortex near tube wall exits on the opposite end with a influence the location of the maximal temperature.
temperature higher than the temperature of inlet air. Vortex
tube has many possible industrial applications and could be Table 1 Workpiece basic data
used as a cooling device at CNC machines, in refrigerators, Description Values
Item name Inner ring
heating processes, etc. High applicability of vortex tube is Number of workpieces 27
based on its simplicity, compactness and the fact that the Material designation 42CrMo4V
system has a small mass and works in quiet mode. These Rolled ring, hardened and tempered at 800-
Material condition
tubes do not have moving parts, so they do not break or wear 900 MPa
which makes them simple for maintenance. Final machining of induction hardened
Machining type
raceway with 56 HRC ±2 HRC for ball ∅50
Application of cold air will decrease the temperature at
the cutting area during machining; it will decrease the
temperature of tool, chip and workpiece due to intense heat
removal by convection which emphasizes the importance of
convection coefficient for tool temperature modelling.
Convection coefficient for cutting fluids based on water is in
the range from 103 to 104 W/(m2K). Convection coefficient
for cold air was investigated for the first time by Liu and
Chou [12] and in simulation, it has values in range 50-5000
W/(m2K), while in the experiment it is about 160 W/(m2K)
with applied air of temperature up to −15 °C and 860
W/(m2K) for cooling with air with temperature up to −25 °C.
• decrease stress in the ring which is caused by thermal Optimal machining parameters according to the
expansion due to higher ring temperatures during recommendations from literature presents a combination of
machining, low feed rate and low depth of cut with higher cutting speed,
• perform dimensions in given tolerances, which is beneficial for reducing cutting temperatures,
• decrease tool costs. machining force and surface roughness (Tab. 4).
All of these requirements should be realized with very Table 4 Basic data for insert and final cutting parameters
limited funding resources. Description Values
RCGX Full Face PCBN grade Insert size 0900700T-
25VM
2.2 Suggested Solution
Insert type
Analysis of literature and taking into account economic
and other limiting factors resulted in the proposal of
application of cold air for cooling of the tool and workpiece.
Installing a cold air nozzle is simple and requires only the CRDCN 3225P 09-A
installation of compressed air hose on the column of the
machine. Accordingly, a SARA vortex tube was purchased. Tool holder
Measuring instrument
4 CONCLUSION
5 REFERENCES
Figure 7 Detail of final machining with spot cooling
[1] Dyck, M. (2007). Beitrag zur Analyse thermisch bedingter
3 DISCUSSION Werkstückdeformationen in Trockenbearbeitungsprozessen.
Dr.-Ing. Dissertation wbk Institute of Production Science,
The following results have been obtained: Karlsruhe.
[2] Klocke, F., Kneer, R., Burghold, M., Deppermann, M., Peng,
• First, it was noticed that there is no characteristic color
B., & Puls, H. (2018). Modelling and Compensation of
of the machined chip. Thermoelastic Workpiece Deformation in Dry Cutting. In:
• There is no characteristic point of burnout chips on Thermal effects in complex machining processes. Lecture
cutting edge due to the fact that cold air decreases Notes in Production Engineering. Springer, Cham,
temperature in point of contact for cutting edge and Switzerland, 63-94. https://doi.org/10.1007/978-3-319-57120-1
raceway. [3] Aydin, M., Karakuzu, C., Uçar, M., Cengiz, A., & Çavuşlu, M.
• During machining without cooling only one side of the A. (2013). Prediction of surface roughness and cutting zone
insert was used for machining of one raceway (1 ring), temperature in dry turning processes of AISI304 stainless steel
and after the application of spot cooling, one side of the using ANFIS with PSO learning. International Journal of
Advanced Manufacturing Technology, 67(1-4), 957-967.
insert performed machining on two raceways (2 rings).
https://doi.org/10.1007/s00170-012-4540-2
That was a surprising result because the insert is [4] Weinert, K., Inasaki, I., Sutherland, J. W., & Wakabayashi, T.
expansive (135 €/piece) and it has only two sides (round (2004). Dry machining and minimum quantity lubrication.
insert which could be set in two positions). The CIRP Annals - Manufacturing Technology, 53(2), 511-537.
durability of the insert was increased for 100% during https://doi.org/10.1016/S0007-8506(07)60027-4
machining of induction hardened raceway with one pass [5] Sharma, V. S., Dogra, M., & Suri, N. M. (2009). Cooling
time of about 1 hour. techniques for improved productivity in turning. International
• The number of workpieces was 27 inner and 27 outer Journal of Machine Tools and Manufacture, 49(6), 435-453.
rings. It means that with spot cooling one insert https://doi.org/10.1016/j.ijmachtools.2008.12.010
[6] Nandy, A. K., Gowrishankar, M. C., & Paul, S. (2009). Some
machined 4 rings instead of only two; therefore, to
studies on high-pressure cooling in turning of Ti-6Al-4V.
machine the rings 14 inserts were spent instead of 28 International Journal of Machine Tools and Manufacture,
inserts without cooling. This allows for significant
Authors’ contacts:
Iryna SOLONENKO
Abstract: The development of road network infrastructure is an important component of the economic development of the European Union. Updating of the road network
contributes to the integration of the economies of countries into a coherent whole. The road network provides the free movement of citizens, the movement of goods and the
effective implementation of various services. The increase in the length of the road network leads to an increase in the financial and material costs necessary to ensure its
maintenance and repair. One of the ways to reduce costs is by strengthening the physic-mechanical and operational characteristics of the pavement due to the widespread use
of cement concrete. The quality of the pavement of cement concrete depends largely on the rational selection of its composition. This allows a significant increase in the durability
of road pavement. The purpose of the research was: the development of recommendations for the rational selection of the composition of the road pavement material of cement
concrete, aimed at upgrading longevity, and taking into account its frost resistance grade. According to the goal, the following tasks were developed: the analyses of the climatic
zones in which the road network of the European Union is located; the development of a research plan, a selection of the response function and influence factors; the study of
physico-mechanical and operational characteristics of the researched material of road pavement; on the basis of the obtained data, the calculation of the complex of experimental-
statistical models, which describe the physico-mechanical and operational characteristics of the road pavement material; on the basis of experimental statistical models, a method
was proposed for selecting the rational compositions of the cement concrete pavement road material depending on the conditions of its application. The results presented in the
article can be used in engineering and scientific practice for the selection of road pavement from cement concrete for highways.
Keywords: cement concrete; experimental statistical models; highway; road pavement; transport corridor; weather conditions
1 INTRODUCTION Table 1 The length of main roads, highways and specific sections of roads in the
EU countries [4]
Ground
The network of European roads has a considerable No. Country
Main-line Hard road
road General
length. They ensure the stable development of the European highway pavement
pavement
Union (EU) countries. European roads are combined in a 1 Austria 2 223 200 000 0 202223
number of transport corridors. The main direction of traffic 2 Belgium 1 763 120 514 33 498 155775
in these corridors: from north to south and from west to east. 3 Bulgaria 801 43 649 440 44890
4 Great Britain 3 557 344 000 54 350 401907
They enable the free movement of citizens and various goods 5 Hungary 1 481 76 075 123 492 201048
and the performance of necessary services [1]. EU roads are 6 Germany 12 917 644 480 0 657397
directly adjacent to the road network of Belarus, Ukraine, 7 Greece 2 311 107 406 9 594 119311
Moldova and Russia. The distinctiveness of the roads of 8 Denmark 1 205 74 558 0 75763
Ukraine and Moldova is the low quality of the road 9 Ireland 1 224 91 145 5 457 97826
10 Spain 16 583 683 175 0 699758
pavement, which is caused by the insufficient allocation of
11 Italy 6 758 487 700 0 485458
funds for their maintenance and repair. 12 Cyprus 254 8 564 4 442 13260
As the analysis of the works [2, 3] shows, the condition 13 Latvia 0 14 707 57 737 72444
of the roads mainly depends on the condition of the road 14 Lithuania 0 13 584 8 242 21 826
pavement. The defects of road pavement (Fig. 1) 15 Luxembourg 152 2 899 0 3051
significantly complicate the effective operation of vehicles, 16 Malta 0 2 704 392 3 096
17 Netherlands 2 808 139 295 0 142103
they reduce their speed and endanger traffic safety. 18 Poland 1 566 292 134 131 863 425563
19 Portugal 2 992 71 294 11 606 85892
20 Romania 806 49 873 34 312 84 991
21 Slovakia 432 38 085 5 676 44193
22 Slovenia 618 38 985 0 39603
23 Finland 863 51 016 27 146 79025
24 France 11 882 1 028 446 0 1040328
Rutting Chuck holes, potholes Cracks 25 Croatia 1 318 26 958 0 28276
26 Czech 1 250 130 671 0 131921
27 Sweden 2 050 579 564 0 581614
28 Estonia 115 10 427 47 985 58 527
Figure 2 Relative construction cost (РΣ), maintenance and repair of roads An analysis of the works [8, 9, 11] showed that the
bitominous with asphalt concrete and cement concrete (a - asphalt concrete, c - durability of the coating for highways made of cement
cement concrete, Pst - construction, Pe - content, Pr - repair, Te – road operation in
time (in years), Pstc = 1,8·Psta)
concrete is significantly affected by the frost resistance of the
material. The resistance of the pavement material to the
As it can be seen from Fig. 3, the total cost of effects of alternating temperatures is important for climatic
construction and maintenance of highways made from zones with unstable negative temperatures (Fig. 2, Tab. 3). In
cement concrete (c), after 10 years of operation, is lower than these zones, throughout the entire winter period, repetitive
for the roads with asphalt concrete (a). Thus, the use of cyclical freeze-thawing of the road pavement is happening.
cement concrete coatings for the European Union countries
reduces the cost of their maintenance and repair.
Table 3 Temperature range in winter time depending on the climatic zones - for each point of the experiment plan, the required
No. of The number of freeze-thaw cycles number of samples was formed in the sizes of
Celsius temperature (°C)
zones for 1 year for 10 years
3 −40 … −34 4 40
0.1×0.1×0.1 meters and 0.07×0.07×0.07 meters [13];
4 −34 … −29 5 50 - the obtained samples (Fig. 4a)) were maintained under
5 −29 … −23 5 50 the standard curing of a normal set for 28 days (t = 20
6 −23 … −18 6 60 °C, W = 80%) [13];
7 −18 … −12 8 80 - on the 28th day, samples of 0.07×0.07×0.07 meters were
8 −12 … −7 12 120
9 −7 … −1 20 200
tested for an abrasion test (LKI-3 device) [14];
10 −1 … +4 55 550 - the samples of 0.1×0.1×0.1 meters were tested for
11 more +4 0 0 compressive strength [13];
- the part of the samples of 0.07×0.07×0.07 and
Table 4 Influencing factors and the range of their changes used in research 0.1×0.1×0.1 meters were tested for frost resistance in the
Impact factors freezer (in the freezing room temperature of −50 °C
х1 – air entraining agent additive Mapeplast РТ-1
[15]);
х2 – polypropylene fiber MAPEFIBRE NS 12/ NS 18
х3 – microsilicasuspension filler - after testing for frost resistance (F50, F 100, F150, F
Range of factors change 200), the samples (Fig. 4b)) were tested for the abrasion
Levels of varying factors x1, l/m3 x2, kg/m3 x3, kg/m3 test and compressive strength [13, 14].
Interval 0,142 0,3 7,5
Upper (+1) 0,285 0,6 15 Table 5 The plan of the experiment and the compositions of the studied concretes
Lower (−1) 0 0 0 The plan of experiment The composition of concrete
РТ-1, Fibre, Microsilica suspension,
No х1 х2 х3
l/m3 kg/m3 kg/m3
The main characteristics of the additives and fillers used: 1 −1 −1 −1 0 0 0
- air entraining agent additive Mapeplast РТ-1 2 −1 1 −1 0 0.60 0
("MAPEI"), its inclusion into the concrete mix provides 3 0 0 −1 0.142 0.30 0
an increase in the amount of the entrapped air (designed 4 1 −1 −1 0.285 0 0
to increase the frost resistance of concrete); 5 1 1 −1 0.285 0.60 0
6 −1 0 0 0 0.30 7.50
- microsilica suspension, particle size 0.1-0.2 microns,
7 0 −1 0 0.142 0 7.50
specific surface area – 20.000 m2/kg (designed to 8 0 0 0 0.142 0.30 7.50
increase the strength of concrete and frost resistance, 9 0 1 0 0.142 0.60 7.50
reduces abrasion and the permeability of concrete, etc.); 10 1 0 0 0.285 0.30 7.50
- polypropylene fibre - Mapefibre NS 12/NS 18 11 −1 −1 1 0 0 15
("MAPEI"), diameter – 0.34 mm, fibre length 12-18 mm, 12 −1 1 1 0 0.60 15
13 0 0 1 0.142 0.30 15
density – 9,1×10-4 kg/m3, tensile strength – 700 MPa, 14 1 −1 1 0.285 0 15
(increases the crack resistance of concrete and makes it 15 1 1 1 0.285 0.60 15
resistant to abrasion).
The experiment plan is presented in Tab. 5. The results of the experiments are presented in Tab. 6.
At each point of the experiment plan, at least three The convenience of the analysis for Tab. 6 is shown in a
experiments were conducted with the subsequent graphical form in the Figs. 5 and 6.
determination of the average value of the measurement As it can be seen from the presented data (Fig. 5), the
result. To eliminate the influence of systematic errors caused compositions used for the manufacture of the samples No. 1,
by external conditions, the order of the experiments was 2, 3, 4 can be used as road pavement on the sections of roads
randomized. that are not affected by the freeze-thaw.
The experiments were conducted in the following Under the influence of the freeze-thaw temperature, the
sequence: compressive strength of these samples decreases:
- the necessary amount of additives and fillers were added - for the sample No. 1 at: F50 – 3.1%, F100 – 8.1%, F150
into the basic composition of the concrete mixture (Tab. – 9.5%, F200 – 14.2%;
4); - for the sample No. 2 at: F50 – 4.6%, F100 – 8.4%, F150
–11.8%, F200 – 15.6%;
- for the sample No. 3 at: F50 – 3.2%, F100 – 6.4%, F150 The compressive strength of the remaining samples
– 9.2%, F200 – 13.2%; subjected to frost resistance tests, even after F200, was more
- for the sample No. 4 it comes down to: F50 – 1.9%, F100 than 50 MPa, which allows its use in all climatic zones.
– 5.4%, F150 – 8%, F200 – 9.3%.
Table 6 The results of the experiments on the compressive strength and abrasion under the cyclic effects of the freeze-thaw test (F)
The results of the experiments
fck.cube (МРа) G (kg/m2)
No fck.cube0 fck.cube50 fck.cube100 fck.cube150 fck.cube200 G0 G50 G100 G150 Gf200
1 2 3 4 5 6 7 8 9 10 11
1 50.25 48.70 46.20 45.50 43.15 0.062 0.065 0.067 0.068 0.071
2 52.60 50.30 48.40 46.70 44.80 0.030 0.033 0.036 0.038 0.039
3 51.80 50.20 48.60 47.20 45.20 0.048 0.050 0.052 0.054 0.057
4 50.10 49.15 47.40 46.10 45.45 0.067 0.069 0.072 0.074 0.076
5 53.40 52.90 52.00 51.10 50.35 0.033 0.036 0.039 0.042 0.045
6 54.60 53.90 52.70 52.00 51.45 0.037 0.040 0.043 0.045 0.048
7 54.20 53.80 52.40 51.20 50.35 0.048 0.053 0.057 0.059 0.061
8 55.20 54.70 53.90 52.00 51.45 0.038 0.040 0.042 0.044 0.047
9 54.70 53.15 52.50 51.70 50.75 0.029 0.032 0.035 0.036 0.039
10 53.75 53.00 52.10 51.40 50.50 0.037 0.039 0.041 0.043 0.046
11 56.50 56.10 55.40 54.90 53.90 0.058 0.060 0.063 0.066 0.069
12 57.80 57.00 56.10 55.20 54.30 0.025 0.028 0.033 0.036 0.039
13 55.45 55.00 54.10 53.25 52.80 0.037 0.039 0.043 0.045 0.048
14 54.90 54.00 53.25 52.30 51.80 0.045 0.048 0.054 0.057 0.062
15 55.00 54.50 54.00 53.85 53.20 0.029 0.031 0.034 0.036 0.038
Figure 5 The results of the testing samples for abrasion before and after freezing
and the thawing test (F0, F50, F 100, F150, F 200)
value of the upper level was equal to +1, and the lower to –1. I
= xi − xi 0
The origin of the coordinates of the influence factors was (3)
transformed and the transition to the normalized value of
each factor was made: The description of the methods for calculating the
coefficients of regression models is beyond the scope of this
article, but if necessary, you can refer to the works [16, 17].
( хi − xi 0 ) The calculated mathematical models describing the
хi =
I (2) change in the compressive strength (fck.cube0, fck.cube200) and
abrasion (G0, G200) of the pavement material before and after
the freezing thawing resistance test of the test samples (F
Where: xi – normalized value; хi – natural value; хi 0 –
200) are presented in Tab. 7.
main level; I – variability interval:
Table 7 Mathematical models describing the change in the compressive strength (fck.cube0, fck.cube200) and abrasion (G0, G200) of the pavement material before and after the freezing
thawing resistance test
No Response function ES-models
1 fck.cube0 (МPа) = 54.49 – 0.46х1 – 0.63х1х3 + 0.75х2 – 0.53х2х3 + 2.15х3 – 0.71x32 (4)
2 G0 (kg/m2) = 0.04 + 0.002x1x2 – 0.002x1x3 – 0.01х2 + 0.002x2x3 – 0.04x3 + 0.004x32 (5)
3 fck.cube200 (МPа) = 51.29 + 0.28х1 – 0.85х1х3 +0.24х2 – 0.31х2х3 +2.69х3 (6)
4 G200 (kg/m2) = 0.05 + 0.001x1x2 – 0.002x1x3 – 0.01х2 + 0.003x22 + 0.001x2x3 – 0.003x3 + 0.006x32 (7)
Figure 7 The example of using the mathematical models No. 1 and 3 for assessing
the resistance of the sample No. 5 to the effect of the freeze-thaw cycles (F) on
compressive strength
Abstract: This article describes the procedure of calculating deflection of rectangular plate using a finite difference method, programmed in Wolfram Mathematica. Homogenous
rectangular plate under uniform pressure is simulated for this paper. In the introduction, basic assumptions are given and the problem is defined. Chapters that follow describe
basic definitions for plate bending, deflection, slope and curvature. The following boundary condition is used in this article: rectangular plate is wedged on one side and simply
supported on three sides. Using finite difference method, linear equation system is given and solved in Wolfram Mathematica. System of equations is built using the mapping
function and solved with solve function. Solutions are given in the graphs. Such obtained solutions are compared to the finite element method solver NastranInCad.
1 ∂ ∂w ∂2 w Ez 1 1
=
− = − 2 (4) =σy +ν
2
(12)
ry ∂y ∂y ∂y 1 −ν ry rx
1 ∂ ∂w ∂2 w
=− = − 2 (5)
rn ∂n ∂n ∂n
1 1 1 1
= cos 2 α − sin 2α + sin 2 α (6)
rn rx rxy ry
where:
∫ σ x z dydz = M x dy (13)
When bending the beam, the assumption is that the beam − h /2
section is unformed and rotates around the neutral axis so that h /2
its position is normal in relation to the elastic line. Axis z is ∫ σ y z dxdz = M y dx (14)
perpendicular to the middle plane and directed downwards − h /2
(Fig. 3).
Elongations of elementary layer abcd are shown in Fig. Substituting (9) and (10):
4:
1 1
z M x D +ν
= (15)
εx = , (9) rx ry
rx
1 1
z M y D +ν
= (16)
εy = . (10) ry rx
ry
Ez 1 1 Eh3
=σx +ν (11) D= (17)
2
1 −ν rx ry 12(1 −ν 2 )
h/2
Qx = ∫ τ xz dz (18)
−h/ 2
h /2 Figure 6 Middle plane of differential element [1]
Qy = ∫ τ yz dz (19)
5 BOUNDARY CONDITIONS
− h /2
Taking into consideration all moments and forces the Assuming the plate is rectangular and the edges are
following equilibrium equation is: parallel to coordinate axes, we set boundary conditions for
the plate.
∂M yx ∂M x
dxdy + dydx − Qx dxdy =
0 (20) 5.1 Wedged Edge
∂y ∂x
If the edge is wedged, deflection on the edge equals zero
Simplified: and tangential plane is corresponding to middle plane of
unloaded plate. For edge x = a
∂M yx ∂M x
+ 0
− Qx = (21) ( w) x = a = 0 , (27)
∂y ∂x
∂w
=0 (28)
Equation of differential element is: ∂x x = a
∂2 w ∂2w ( w) x = a = 0 (29)
− D 2 + ν 2
Mx = (23)
∂ w2 2
∂x ∂y ∂ w
2 + ν 2 0
= (30)
∂x ∂y x = a
∂ w 2 2
∂ w
− D 2 + ν 2
My = (24)
∂y ∂x
dφ φ ( xi + ∆x) − φ ( xi )
= ∆lim (31)
d
x x → 0 ∆x
dφ φi − φi −1 φi − φi −1
=
= (34)
dx i xi − xi −1 ∆x
Central difference:
dφ φi +1 − φi −1 φi +1 − φi −1
=
= (35)
dx i xi +1 − xi −1 2∆x
∂w wi +1, j − wi −1, j
= (36)
∂x 2∆x
Figure 7 Finite difference mesh ∂w wi , j +1 − wi , j −1
= (37)
Fig. 8 shows derivation as a tangent of a function ϕ(x). ∂y 2∆y
We can approximate the tangent of a function with the line ∂2 w wi +1, j − 2 wi , j + wi −1, j
which represents forward difference, backward difference or 2
= (38)
∂x ∆x 2
central difference [5]
Forward difference is given by: ∂2 w wi , j +1 − 2 wi , j + wi , j −1
= (39)
dφ φi +1 − φi φi +1 − φi ∂y 2 ∆y 2
= = (32)
dx i xi +1 − xi ∆x ∂4 w wi + 2, j − 4 wi +1, j + 6 wi , j − 4 wi −1, j + wi − 2, j
4
= (40)
∂x ∆x 4
Forward difference is calculated by:
∂4 w wi , j + 2 − 4 wi , j +1 + 6 wi , j − 4 wi , j −1 + wi , j − 2
= (41)
∂y 4 ∆y 4
7. We solve boundary equations. 10. We make table of coordinates in three dimensional space
(*takes the left boundary conditions, and solves them and plot solution.
to yield the values of w along the left side *) dataPoints=Table[{xmin+i*dx,ymin+j*dy,solutionArr
wleftw=NSolve[Map[#==0&,leftbcw],Table[W[[3, ay[[i+1,j+1]]},{i,0,Length[xgrid]-5},{j,0,Length[ygrid]-
k]],{k,4,Length[ygrid]-2}]]; 5}];
wleft=NSolve[Map[#==0&,leftbc],Table[W[[2,k]], Min[solutionArray]
{k,4,Length[ygrid]-3}]]; ListPlot3D[Flatten[dataPoints,1],AxesLabel->
(*takes the right boundary conditions, and solves {Style[x,Medium,Blue],
them to yield the values of w along the right side *) Style[y,Medium,Blue],Style[u,Medium,Blue]},
wrightw=NSolve[Map[#==0&,rightbcw],Table[W[ PlotRange->All,ColorFunction-
[-3,k]],{k,3,Length[ygrid]-3}]]; >Function[{x,y,z},Hue[z]]]
wrightt=NSolve[Map[#==0&,rightbc],Table[W[[-
2,k]],{k,4,Length[ygrid]-3}]];
7 SOLUTION
wright=Flatten [wrightt/.wrightw,1];
(*takes the bottom boundary conditions, and solves Using Wolfram Mathamatica program, the program code
them to yield the values of w along the bottom side *) from the previous chapters, we obtained maximum deflection
wbottomw=NSolve[Map[#==0&,bottombcw],Tab of 23.874 mm. A three dimensional graph of this solution is
le[W[[k,3]],{k,3,Length[xgrid]-3}]]; given in Fig. 10.
wbottomm=NSolve[Map[#==0&,bottombc],Tabl
e[W[[k,2]],{k,4,Length[xgrid]-3}]];
wbottom=Flatten [wbottomm/.wbottomw,1];
(*takes the top boundary conditions, and solves them
to yield the values of w along the top side *)
wtopw=NSolve[Map[#==0&,topbcw],Table[W[[k,
-3]],{k,4,Length[xgrid]-2}]];
wtopp=NSolve[Map[#==0&,topbc],Table[W[[k,-
2]],{k,4,Length[xgrid]-3}]];
wtop=Flatten [wtopp/.wtopw,1];
(*This is a list of all the boundary values of w*)
boundary1=Flatten[Join[wleft,wright,wtop,wbott
om]];
boundary2=Flatten[Join[wleftw,wrightw,wtopw,wb
ottomw]];
8. We now create a set of equations. We make a table, with Figure 10 Finite difference solution
each entry corresponding to an interior grid point. Each
entry in the table becomes an equation, from the In order to prove that this solution is correct, the
discretized partial differential equation. We use our NastranInCad finite element analysis program was used.
knowledge of the boundary conditions to eliminate the NastranInCad is used as a part of Autodesk Inventor. In the
values of w on the boundary [7, 8]. Inventor, the finite element method has the maximum error
of 23.67 mm. This solution is shown in Fig. 11.
equations=Map[(q0/DD==#)&,Flatten[Table[dwdx4[[i,
j]]+2*dwdx2dy2[[i,j]]+dwdy4[[i,j]],{i,2,xdivisions},{j,2,y
divisions}],1]/.boundary1/.boundary2];
5 REFERENCES
Abstract: Innovations and market changes in warehouse and logistics systems force the adaptation and transformation of the existing business model into a business model
based on modern technology. With the development of the Internet, RFID (radio-frequency identification) technology and sensors, new innovations are being created that allow
the improvement of the existing mode of activity. Implementation of the new technology brings along a number of challenges that organizations must find an adequate response
to. However, warehouse systems are not the only affected by the new technologies. The development of technology and technological innovations enable organizations to develop
sustainability. Sustainable development is imperative due to increasing awareness of the need for environmental protection. The supply chain can also be managed much more
efficiently if sensors that collect information of customer habits and process performance are implemented. Given the accelerating development of Industry 4.0 and the opportunities
offered by newly developed technology, this paper provides an overview of current developments in the implementation of Industry 4.0 technological innovations in logistics.
Keywords: Internet of things; radio-frequency identification; smart technologies; supply chain; warehouses
Chinese economy recognized the importance of this of Industry 4.0 on the logistics system. But there are few that
technology and started investing in its development. [4] provide an overview of all technologies affecting the logistics
Automation of production began in the 1950s, when system. In accordance with this, this chapter provides an
robots were used for the first time in manufacturing for overview of all technologies affecting the logistics system.
simple operations [5]. From 2000 to 2015, the use of
industrial robots increased by 150% [6], whereby the trend of 3.1 Warehouse Systems and Industry 4.0
growth is evident. Automation of the manufacturing process
has often been applied in the automotive industry where in The warehouse system itself contains a large number of
the past 9 of 10 robots were sold to one of the car operations aimed at meeting the requirements of the
manufacturers. Furthermore, the importance of robotics and interested parties. Lee, C.K.M., et al (2017) state that the
automation also comes from the fact that the number of development of industry 4.0 has affected a different way in
robots that are implemented in the todays automotive which past activities within the warehouse are taking place.
industry is 50% of the total number of manufactured robots The fourth industrial revolution at the same time marked the
[7] which can imply a gradual return to human labor, or an use of the IoT technology in warehouses which reflected in
increase in the number of industrial sectors that apply robots the efficiency and performance of the storage systems.
in their business. Traditional WMS (warehouse management system) is
Parallel to the development and implementation of enhanced by the IoT, which exploits the full potential of
robots in industrial production warehouses are implementing RFID technology. With the IoT all activities as well as goods
automated systems that allow faster and easier transportation within the warehouse are controlled, which enables their
of goods [7]. manageability. However, the use of the IoT also entails the
Traditionally, in warehouse, all products have been security issues of such systems as they may often become the
labeled with RFID technology developed during the 70s of target of online attacks. [12] Industry 4.0 does not only
the last century, enabling a revolutionary approach to involve the use of the IoT, but also of other technologies that,
improving warehouse operations and reducing storage costs. through communication as one of the feedback components,
The development of RFID technology is based on IFF have an impact on the improvement of logistic processes.
(Identification friend or foe) technology used during World The IoT does not only affect warehouses and warehouse
War II. [8] Although the purpose of RFID was initially to be operations, but also the entire supply chain. Sensors as an
built into the automotive industry, the banking industry, integral part of this technology collect a large amount of data
security and healthcare, [9] later adaptation and refinement that affects the challenges associated with their storage and
proved to be an ideal storage system. Robotization, management. Machado and Shah emphasize that precisely
automation and internet development transform traditional because of this fact the existing warehouses have to be
warehouses into warehouses that, through the IoT and reorganized. In addition, by deploying intelligent systems to
modern computer support, improve their performance and, the warehouse, the logistic superstructure, which has so far
above all, productivity, efficiency and effectiveness. been managed by employees, becomes autonomous,
Growth of globalization and the strengthening of allowing it to move through the warehouse without or with a
capitalism increases the need for industrial development. slight employee intervention. Furthermore, IoT provides
Industry development results in increasing emissions of immediate insight into inventory stock status. Particular
harmful gases. The first attempts to raise awareness of the importance is the ability to manage the entire supply chain as
issues and challenges related to ecological problems were well as communication between its stakeholders. Real-time
first shown in 1962. In 1987, in the report published by information enables logistical operators to gain insight into
WCED, the concept of sustainable development is mentioned logistics operations, in particular the conditions under which
for the first time in the history. Five years later, in 1992, the freight transport is carried out by means of transport, which
basic environmental goals were defined. In 2015, the UN increases transparency and reduces the possibility of damage
defined the 17 sustainable development goals recommended to goods [13].
by Member States [10], one of which is closely linked to IoT.
[11] 3.2 Warehouses Based on Cyber-Physical System
Past trends have transformed traditional approaches to
business and sustainability management into a modern Innovations that occur within industry 4.0 create new
approach based on technologies developed alongside the technologies that facilitate system management. The Cyber-
development of needs for a different approach and Physical System (CPS) can track and create a virtual copy of
environmental protection relationship. the actual process that can be used to monitor process
performance. In addition, the CPS is very similar to the IoT
3 REVIEW OF RESEARCH concept that allows all components within the system to
communicate with one another i.e. that the physical
The topic of sustainable logistics and the IoT has often components are virtually connected, allowing for cost
been explored and there is a considerable number of papers savings in parallel with increasing efficiency. But, there is
that closely deal with this topic, especially the impact of the one big difference between the CPS and the IoT. The CPS is
IoT on warehouse and logistics systems. However, there are only focused on physical objects while it may be focused on
many of research papers that address the topic of the impact creating network of things that are focused to service
providing. Warehouses based on this technology include sources, energy-efficient building materials which, by their
RFID sensors, Bluetooth technology, Wi-Fi access points, construction, are insulators, improvements in logistics supra-
cameras and robots that are coordinated in the system to structure by improving the forklift, by installing automatic
perform a defined task. The role of a man in such a warehouse closing doors, etc. [18] On the other hand, Amith and
is primarily related to monitoring and reprogramming the Harrison (2013) emphasize that ecological acceptance of the
system if necessary. However, the CPS technology also storage system is ensured through awareness of the impact of
enables inter-machine co-operation between robotic systems, such a system on the environment and suggest a warehouse
thus reducing the need for human work. Such designed design that will maximize the use of natural lighting, water
robotic systems also enable human movement to recognize supply and wastewater management systems, noise reduction
what enables employees to use robotic help in carrying out to the surrounding community and improved temperature
activities. [14] By replacing human labor with automated and control. [19]
robotized systems, or through the implementation of such The transformation of the traditional logistics system
systems as humanitarian aid, the efficiency of the warehouse into the environmentally sustainable results in the
system increases. Activities that are dangerous for a man can improvement of all elements of sustainable development.
be robotized, thus reducing the risk of injury and unlucky Reducing costs increases funding available to improve
cases. ecological acceptability, which also affects the social
component of sustainability. Furthermore, through the
3.3 IoT Conducted Inventory Management optimization of the commissioning route through the WMS,
it affects the performance that the logistics processes
The possibilities of using the IoT technology are growing develop, that is, their efficiency, which ultimately results in
in parallel with innovations. In the context of warehousing greater sustainability. [18]
systems, the IoT allows the connection of previously Many authors point out that logistics is one of the
unreliable physical components to the network, thus fundamental components that needs to be paid particular
managing the storage facility and facilitating it. Particularly attention, especially in large cities where it is timely
problematic components of each warehouse are inventories supplying the necessary resources, crucial to the
and operations associated with their management. development of all processes inherent in cities. [20] Robotics,
Traditional systems in most cases include manual automation, electrification, sensing as well as the application
manipulations of inventory as well as inventory management of the IoT technology, the storage system is managed more
based on human. However, by deploying sensors and linking efficiently, which reduces the negative impact on the
them to the network, sensors can look at real-time inventory, environment.
which greatly facilitates their management. Therefore, the
new warehouse system is based only on IT with human as an 3.5 Future of Warehouse Systems
operator that only controls the process. This approach is
based on RFID technology, which is the basis for deploying By implementing technologies of Industry 4.0, such as
and using this approach. [15] The IoT and RFID technology the IoT, warehouse systems become more efficient and their
not only simplify inventory deployment, but anticipate future performance is improving. By automating warehouse
orders, tracking product durability, temperature, humidity in operations, the warehouse system becomes less dependent
the air, and other parameters that could affect product and responsive to human work variability. By applying
damage. RIFD technology is the foundation for the autonomous systems in transportation, transshipment-
functioning and development of the IoT system. [16] The loading operations, and operations of palletized goods to
application of RFID technology in warehouse management warehouse regal, the flow of goods through warehouses
has a number of advantages that primarily relate to increasing becomes more fluent. However, autonomous systems must
efficiency, accuracy and updating of information, as well as be tailored to the storage environment to minimize the risk of
tracking stocks and inventory losses. [17] By implementing collisions. Existing autonomous systems in a large number of
the IoT and RFID technology in organization’s storage cases, by having a barrier remain standing in place, that is,
systems, they provide competitive advantage based on lower they can only move through previously defined and narrowly
operating costs and greater efficiency of the process. restricted corridors. [21]
A potential solution to such and similar problems of
3.4 Impact of Warehousing on the Environment implementation is the development of artificial intelligence
that is capable of making decisions similar to man-made
Warehousing, just like logistics, is an activity that decisions. Robotic systems with artificial intelligence enable
generates a certain amount of greenhouse gases. Poor energy making simpler decisions and adaptation to new
efficiency of warehouses as well as environmentally circumstances in the environment. [22] As such, artificial
unacceptable heating, air conditioning or lighting systems are intelligence is already present in a number of daily-used
just some of the factors that affect the amount of carbon technologies such as internet browsers as well as systems for
dioxide discharged. In order to increase the ecological predicting earthquakes and weather conditions. [23] The
acceptance of storage processes, Đukić, Česnik and Opetuk Deep Learning concept moves artificial intelligence and its
(2010) propose the implementation of energy-efficient application to a higher level through the ability of such
lighting, the production of electricity from renewable energy systems to learn and improve their knowledge without the
[6] Grau, A., Indri, M., Bello, L. L., & Sauter, T. (2017, October). in the future warehouse. In Proceedings of the 2016 ACM
Industrial robotics in factory automation: From the early stage International Joint Conference on Pervasive and Ubiquitous
to the Internet of Things. In IECON 2017 - 43rd Annual Computing: Adjunct (pp. 1755-1758). ACM.
Conference of the IEEE Industrial Electronics Society (pp. https://doi.org/10.1145/2968219.2968270
6159-6164). IEEE. https://doi.org/10.1109/IECON.2017.8217070 [27] Giannikas, E., Woodall, P., McFarlane, D., & Lu, W. (2017,
[7] Robots in automobile industry. https://www.slideshare.net/ July). The impact of B2C commerce on traditional B2B
NirajRajan/robots-in-automobile-industry-59415113 warehousing. In ISL 2017: 22nd International Symposium on
(10.08.2018) Logistics (pp. 375-383). Nottingham University Business
[8] Ahuja, S. & Potti, P. (2010). An introduction to RFID School.
technology. Communications and Network, 2(3), 183-186. [28] Rozados, I. V. & Tjahjono, B. (2014, December). Big data
https://doi.org/10.4236/cn.2010.23026 analytics in supply chain management: Trends and related
[9] Cardullo, M. (2005). Genesis of the versatile RFID tag. RFID research. In 6th International Conference on Operations and
Journal, 2(1), 13-15. Supply Chain Management, Bali.
[10] Hedenus, F., Persson, M., & Sprei, F. (2016). Sustainable [29] Toshniwal, R., Dastidar, K. G., & Nath, A. (2015). Big data
Development. History, Definition & the Role of the security issues and challenges. International Journal of
Engineer. Chalmers University of Technology. Innovative Research in Advanced Engineering (IJIRAE), 2(2).
[11] The effect of the Internet of Things on sustainability. [30] Terzi, D. S., Terzi, R., & Sagiroglu, S. (2015, December). A
https://www.weforum.org/agenda/2018/01/effect-technology- survey on security and privacy issues in big data. In Internet
sustainability-sdgs-internet-things-iot/ (12.08.2018) Technology and Secured Transactions (ICITST), 2015 10th
[12] Lee, C. K. M., Lv, Y., Ng, K. K. H., Ho, W., & Choy, K. L. International Conference for (pp. 202-207). IEEE.
(2018). Design and application of Internet of things-based https://doi.org/10.1109/ICITST.2015.7412089
warehouse management system for smart [31] Jadeja, Y. & Modi, K. (2012, March). Cloud computing-
logistics. International Journal of Production Research, 56(8), concepts, architecture and challenges. In Computing,
2753-2768. https://doi.org/10.1080/00207543.2017.1394592 Electronics and Electrical Technologies (ICCEET), 2012
[13] Machado, H. & Shah, K. (2016). Internet of Things (IoT) International Conference on (pp. 877-880). IEEE.
impacts on Supply Chain. Retrieved, 19, 2016. https://doi.org/10.1109/ICCEET.2012.6203873
[14] Ding, W. (2013). Study of smart warehouse management [32] Inukollu, V. N., Arsi, S., & Ravuri, S. R. (2014). Security
system based on the IOT. In Intelligence computation and issues associated with big data in cloud computing.
evolutionary computation (pp. 203-207). Springer, Berlin, International Journal of Network Security & Its Applications,
Heidelberg. https://doi.org/10.1007/978-3-642-31656-2_30 6(3), 45. https://doi.org/10.5121/ijnsa.2014.6304
[15] Yerpude, S. & Singhal, T. K. (2018). SMART Warehouse with [33] Roth, M., Klarmann, A., & Franczyk, B. (2013). Future
Internet of Things supported Inventory Management System. logistics-challenges, requirements and solutions for logistics
[16] Jia, X., Feng, Q., Fan, T., & Lei, Q. (2012, April). RFID networks. International Journal of Mechanical, Aerospace,
technology and its applications in Internet of Things (IoT). Industrial, Mechatronic and Manufacturing Engineering,
In Consumer Electronics, Communications and Networks 7(10), 898-903.
(CECNet), 2012 2nd International Conference on (pp. 1282- [34] Stoltz, M. H., Giannikas, V., McFarlane, D., Strachan, J., Um,
1285). IEEE. https://doi.org/10.1109/CECNet.2012.6201508 J., & Srinivasan, R. (2017). Augmented Reality in Warehouse
[17] Krishna, A., Ravinchandra, L., Abdulla, R., & Thang, K. F. Operations: Opportunities and Barriers. IFAC-PapersOnLine,
(2016). Smart Stock Management Control. International 50(1), 12979-12984. https://doi.org/10.1016/j.ifacol.2017.08.1807
Journal of Applied Engineering Research, 11(1), 492-500. [35] Yang, L. & Ye, M. (2014, May). The Design and Development
[18] Đukić, G., Česnik, V., & Opetuk, T. (2010). Order-picking of Intelligent Warehouse Management System based on. NET
methods and technologies for greener warehousing. and Internet of Things. In International Conference on
Strojarstvo, 52(1), 23-31. Logistics Engineering, Management and Computer Science
[19] Amjed, T. W. & Harrison, N. J. (2013). A Model for (LEMCS 2014). Atlantis Press.
sustainable warehousing: from theory to best practices. https://doi.org/10.2991/lemcs-14.2014.98
In Proceedings of the International Decision Sciences Institute [36] Kim, J. Y. & Park, D. J. (2016). Internet-of-Things Based
and Asia Pacific DSI Conference. Approach for Warehouse Management System. International
[20] Rakhmangulov, A., Sladkowski, A., Osintsev, N., & Muravev, journal of Multimedia and ubiquitous Engineering, 11(10),
D. (2017). Green logistics: element of the sustainable 159-166. https://doi.org/10.14257/ijmue.2016.11.10.15
development concept. Part 1. NAŠE MORE, 64(3), 120-126. [37] Liang, C., Chee, K. J., Zou, Y., Zhu, H., Causo, A., Vidas, S.,
https://doi.org/10.17818/NM/2017/3.7 ... & Cheah, C. C. (2015). Automated robot picking system for
[21] Pechanová, Ľ., Stareček, A., Bachár, M., & Caganova, D. e-commerce fulfillment warehouse application. International
(2017). The suggestion for newspapers application federation for the promotion of mechanism and machine
improvement. 1-5. https://doi.org/10.1109/ICETA.2017.8102518 science.
[22] Thamer, Hendrik & Börold, Axel & Yoga Benggolo, Ariandy [38] Mikušová, N., Čujan, Z., & Tomková, E. (2017). Robotization
& Freitag, Michael. (2018). Artificial intelligence in warehouse of Logistics Processes. In MATEC Web of Conferences (Vol.
automation for flexible material handling. 134, p. 00038). EDP Sciences.
[23] See https://ai.google/ https://doi.org/10.1051/matecconf/201713400038
[24] Kim, K. G. (2016). Book Review: Deep Learning. Healthcare
informatics research, 22(4), 351-354.
https://doi.org/10.4258/hir.2016.22.4.351
[25] Brockmann, T. & Godin, P., 1997. Flexibility for the future in
warehouse design. IIE Solutions, 29(7), pp.22-26.
[26] Diete, A., Weiland, L., Sztyler, T., & Stuckenschmidt, H.
(2016, September). Exploring a multi-sensor picking process
Authors’ contacts:
Abstract: The aim of the article is to develop an assessment model for the effect received from the integration of enterprises during the realization of investment projects in the
construction. Having analysed the existing methods of the economic effect evaluation, we came to the conclusion that it is impossible to use them completely to assess the
economic effect of the integration of enterprises during the investment projects in the construction. In the case of an integrated realization of a project in the construction sector, it
is the information and communication links that are combined, and not the enterprises themselves. In addition to that, not all models take into account the synergistic effect. To
solve the above-mentioned problem, we analysed the advantages gained by participants during the realization of the integrated project in construction and the benefits that have
the greatest impact on the synergy effect. The calculation of the effect of the information and communication system management and the innovative nature of the process of its
implementation suggest that this process has a clear investment component. Consequently, the process of enterprise integration should be considered as an investment project
implemented at the level of network association.
Keywords: construction project; information and communication system; integration; synergetic effect
1 INTRODUCTION
The four quadrants of the matrix refer to the varying high. They note that this is the preferred state for
levels of trust and communication. The final state is that of collaboration to occur. The concept of high performance is
integration, where both trust and communication levels are routinely implemented in diverse industries but has received
little attention in the construction industry, where success is characteristics affords ground for affirming or denying the
typically measured using traditional indicators, such as time, existence of a network structure. In our opinion, it is the
cost and quality. Authors suggests that the problem can be interrelationships that can be distinguished from the general
addressed by viewing the construction team as an integrated system of elements and can be used as the basis for
group of participants within a network, rather than as a group identifying the network organizational structures. This is
of participants. The team needs to consist of a cohesive especially important for the network structures of the
network where members focus on building long term investment and construction complex, since they are
relationships that are transferred from activity to activity. characterized by the absence of legal consolidation of
The project delivery method is selected by the owner of integration relationships, such as participation in capital, the
the construction project. That method establishes the transfer of property rights between economic entities.
preconditions for the realization of the project, such as the The effect of network interaction between enterprises has
roles and responsibilities of the participants, including the its own specifics. It is not only the effect of consolidation
legal agreements between the parties by setting the project because enterprises retain some independence (there is no
target and payment method. This tends to be a critical factor saving in management costs, for instant, because the
so if selected correctly, it enables the successful management apparatus of each of the enterprises remains the
implementation of a construction project. Fig. 2 shows three same). At the same time, a systemic dynamic should arise in
levels of integration in construction. the process of interaction and should become the
The aim of the article is to develop the assessment model consequence of it. It means that the coordinated work of
for the effect received from the integration of enterprises independent companies provides an additional synergetic
during project delivery in construction. effect of the system interaction as a whole.
The investigation of the synergy phenomenon began to
2 MATERIALS AND METHODS develop with the activation of integration processes in the
economy, although the very point of the synergetic effect was
Network structures have different configurations described by Aristotle. In the IV century BC, he noted that
including heterogeneous elements that provide the highest "the whole is greater than the sum of its parts" [3]. I. Ansoff
level of collaboration. Such elements in each particular case defined the synergy effect as the possibility of getting the
have their own unique list that provides the energy of better economic effect from the joint work of several
cooperation in each particular situation. The overall business units in comparison with the results of their
economic effect depends not only on the elements the independent activity [4]. In his concept of synergy, he
network consists of, but also on the nature and quality of the focuses on the role the economic base has, on the potential of
links between them. material and non-material assets and their close connection
In modern economic conditions, the decision on the with the potential of the enterprise.
economic agents’ integration considerably allows the change X. Itami considered that the synergy effect received from
of the competitive position of the enterprise. The integration integration is a result of resource utilization efficiency and it
process leads to new subsystem formation and the is growing in an integrated corporation in comparison to their
interconnections between them. At the same time, the use by companies before integration [5]. He identified two
interaction and influence of the subsystems in an integrated types of resources: physical resources (fixed assets,
structure is characterized by the emergence of a new inventories, etc.) and intangible assets (brand, technology
synergetic component, which is considered to be the effect of and know-how, corporate culture, etc.). In his opinion, two
enterprise integration, and is the basis of the goals and types of effects can be formed as a result of integration:
motives of integration. The synergistic effect presence means - complementary effect which refers to the improvement
that the result of the integrated system of enterprises is higher of physical resources usage (the corporation uses its
than the sum of the results of individually operating resource potential in full),
enterprises united in the process of integration. Depending on - synergetic effect which is achieved only by using
synergy effect it can evaluate the effectiveness of the invisible assets.
integration itself.
The efficiency of the integration process depends largely In our opinion, using a network organizational structure
on how effectively the management of various forms of during the delivery process in the construction industry has
interaction between construction enterprises is carried out. the most interesting and at the same time the most difficult
The results of the integrated enterprise's production processes issue to study. It is the issue of information synergy
are not always clear and, in some cases, can lead to an evaluation that is formed due to the optimization of
efficiency decrease, because the enterprises do not fully use information and communication management.
the potential of business units’ interaction and the additional The analysis of theoretical and practical research in the
effects are the results of such potential. If the interaction of field of economic entities integration [6, 7, 8] showed that the
enterprises within the framework of integrated structures efficiency of integration is mainly assessed by minimizing
does not provide the desired effects, it is important to update the costs of the enterprises (they are transactions in most
the methodological approaches that help improve the cases). In this article, only this approach is not sufficient,
mechanism of the integrated structures management. The since it is planned to assess the efficiency of integration not
presence or absence of economic connections with certain only from the position of the interacting entities themselves,
but also from the position of the whole integrated system, The value in Vij – (Vi + Vj) is a synergistic effect. To
within the framework of their interests and its coordination. recoup the integration, this effect should be greater than the
It is expedient to assess the "cumulative effect" of acts of sum (P + E).
initiating integration chains of interaction between the socio- Methods that are used in the framework of the
economic entities during project delivery in construction. comparative approach make it possible to compare the
Such effect can be evaluated only within the framework of a assessed object with similar objects and on this basis, to
synergetic approach to integration evaluation. establish its value. A comparative approach can be used for
Synergy can be classified according to the impact on the both market valuation and the valuation of investment value.
company's cash flows [9, 10] (Fig. 3). However, since there are no two absolutely identical
associations, some average market multipliers are expertly
determined and used to calculate the cost using comparative
approach methods and that reduces the accuracy of the
results.
In practice, the Capital Assets Pricing Method is used
more often and it gives the most objective assessment of the
effect received from the integration. In this approach, the
discounted cash flow method is used to estimate the value of
the company. The very core of the method is to determine the
current value of cash flows that will be received in the future.
In the framework of the Capital Assets Pricing Method, there
are a significant number of models suitable to calculate the
economic effect.
In particular, S. A. Orekhov [15] calculates the effect
received from integration into the holding structure as
follows:
Having analysed the above methods and approaches, we effects that are indicative for the integrated project delivery
can conclude that they cannot be fully used to evaluate the in construction. Taking into account the specifics of this
effect received from enterprise integration during project research, namely, the evaluation of the effect formed after the
delivery in construction, since the vast majority of these information and communication system introduction in the
methods are based on the analysis of enterprise integrated project delivery in construction, we can identify
amalgamation and their cash flows after the takeover. In the the following synergistic effects:
case of integrated project delivery in the construction
industry, it is the information and communication links that I t = ∑ ΔS = S1 + S2 + S3 + S4 + S5 + S6 (6)
are united, but not the enterprises themselves. Additionally,
the synergistic (immaterial) effect is not taken into account
Where: ΔS - the synergy effect, which is formed from the
in all models. Moreover, it is necessary to distinguish more
introduction of the information and communication system
clearly between the terms of economic effect evaluation and
in the integrated project delivery; S1 - the synergetic effect
integration and the synergetic effect which is one of the main
formed from reducing the number of errors and collisions in
factors, but not the only one.
the construction project delivery; S2 - the synergistic effect
from the overall reduction of project duration due to a higher
3 RESULTS AND DISCUSSIONS
coordination of actions; S3 - the synergistic effect of
transaction costs reduction; S4 - the synergetic effect from the
This article proposes, in the framework of the synergetic
cost reduction at the stage of facilities’ management; S5 - the
approach, to take the ratio of the additional income received
synergistic effect of a new value creation; S6 - the synergistic
from the association of enterprises to the costs on this project
effect from the production scale increase.
delivery as a criterion for the economic efficiency of
The next stage is the analysis of the cost structure
enterprise integration programs.
connected with the integrated project delivery. The expenses
The calculation of the economic effect received from
can be divided into two parts, the sustainable (investment)
managing the information and communication system and
and variable (operational) costs.
the innovative character of the process of its implementation
The structure of investment expenses is determined by
shows that this process has a clear investment component.
the task to introduce the integrated model into the network
Therefore, the process of enterprise integration should be
organizational structure of enterprises and is as follows:
considered as an investment project delivered at the level of
a network association. Accordingly, regarding the innovative
process of the information and communication system y
CI=
t ∑ (Ce ∗ Ts bas + Cst ) + CI ce + CI s + CI db +
network model realization, the predictive efficiency of i =1
investments can be expressed by the following economic and (7)
j
mathematical dependence: + ∑ (Cas ∗ Tadapt )
i =1
m
It
∑ (1 + n )t Where: Ce - the average wage rate of an employee who is
i =1 t
Et = m m
(5) having the training; Ts.bas - the time of training, one employee
CI CO coordination and workplace integration; Cst - the average cost
∑ (1 + nt )t + ∑ (1 + nt )t
=i 1 =
t i 1 t of training; CIce - the investment costs for computer
equipment acquisition; CIce - the investment costs for
Where: Et - the efficiency of the information and software acquisition; CIdb - the investment costs for database
communication system operation in the integrated project acquisition; Cas - the average wage rate of the project
delivery over a period of time; It - the total revenue stream manager assistant; Tadapt - the time the assistant needs to
that the integrated enterprise receives from the information adapt in an integrated system; y - the number of employees
and communication system implementation over a period of and workplaces integrated into the system; j - the number of
time; CIt - the investment costs that the enterprise spends to the project manager "assistant".
launch and implement the integrated information and As it is seen from the expression (7), investment costs
communication system; COt - the operational costs to can be divided into the following three groups:
manage the information and communication system of the 1) the costs of training for the employees to learn the new
business entity; t - the analysed period of time; nt - the rate innovative model of network relationships (including
of discounting accepted for calculations; m - the number of both the costs of pay for those who conduct training and
enterprises that integrate. the loss of time during the training period when the
The synergy effect calculation is the most difficult part employees of the enterprise are not engaged in the
in the total revenue stream calculation. It is possible to use performance of their immediate duties);
various factors to assess the synergy, and their sum-total is 2) the expenses for the computer equipment, software and
formed according to the direction and the logic of the database acquisition and adaptation;
research. Therefore, it is a situation similar to the economic 3) the costs to adapt the project manager assistant to the
efficiency integration evaluation, namely the isolation of the integrated system.
The first two groups of expenses are reflected in the COs - the operating costs for software maintenance; COdb is
formula (7) in square brackets and include the training, the operating costs for database maintenance; Tw - the project
computer equipment, software and database acquisition for manager assistant work duration; Ts.add - the time for
all employees participating in the integrated process. The additional training, consultation and formalization of
third term presents the expenses for the project manager knowledge for the employee; x - the number of employees
"assistant". who need additional training and consultations.
Operating costs are given in the following expression: m - the duration of the project, years.
They also form three clear groups:
y j 1) expenses for computer services, software and
CO
= t ∑ (COce + COs + COdb ) + ∑ (Cas ∗ Tw ) + database maintenance;
i 1 =i 1
(8) 2) expenses for the salary of the project manager
x
assistant;
+ ∑ (Ct ∗ Ts add )
3) expenses for the additional training, consultations
i =1
and formalization of knowledge for employees.
Taking into account the expressions (6)-(8), the formula
Where: COce - the operating costs for computer services;
(5) can be expressed by the following economic-
mathematical dependence:
m
S1 + S2 + S3 + S4 + S5 + S6
∑ (1 + nt )t
i =1
Et = (9)
y j y j x
4 CONCLUSIONS 5 REFERENCES
A system of economic and mathematical dependencies [1] Mosey, D. (2009). Early contractor involvement in building
has therefore been formed to assess the effect of enterprise procurement: contracts, partnering and project management,
integration during project delivery in construction. The Chichester: John Wiley & Sons Ltd.
criterion for the efficiency of enterprise integration programs https://doi.org/10.1002/9781444309867
[2] Chinowsky, P. & Songer, A. (2011). Organizational
is the ratio of the additional income received from the management in construction, New York: Spon Press.
association of enterprises to the costs on this project delivery. https://doi.org/10.4324/9780203856109
The additional income is the total synergistic effect that the [3] Aristotel. (2006). Metaphysika (Metaphysics). Moskva: Eksmo
integrated enterprise receives from the information and (in Russian).
communication system implementation. The costs is the [4] Ansoff, I. (2009). Strategicheskiy menedzhment (Strategic
sustainable (investment) and variable (operational) costs on management). Sankt-Peterburg: Piter (in Russian).
this project delivery. The paper proposes how to solve the [5] Asaul, A. N. & Ivanov, S. N. (2008). Snizhenie
issue of separating the share, which is formed as a result of transaktsionnyih zatrat v stroitelstve za schYot optimizatsii
informatsionnogo prostranstva (Decrease in transaction costs
integration processes, from the general effects of the
in construction due to optimization of the information space).
enterprise's economic activity. Finally, the advantages Sankt-Peterburg: ANO IPEV (in Russian).
obtained by the participants of the integrated project in [6] Rowley, J. (2002). Synergy and strategy in e-business.
construction are analyzed and the advantages that have the Marketing Intelligence and Planning, 20(4), 215-222.
greatest impact on the synergistic effect are highlighted. In https://doi.org/10.1108/02634500210431603
our opinion, it is warranted to use projections to make a [7] Franzsen, D., Van der Merwe, R., Buys, S., Bischoff, A.
decision on integration, and retrospective evaluation should Courtney, F., Varachia, A., Cox, E., & Miller, C. (2001).
be used to determine how well the business combination Synergies in the distribution business of Barloworld.
operates. Further studies will be related to the practical Unpublished manuscript, Barloworld Leadership Development
Programme, Gordon Institute of Business Science,
application of the model proposed to evaluate the economic
Johannesburg.
effect received from the integration of enterprises during [8] Naudé, A., Heyns, D., Bester, C., Puig, J. & Tucker, G. (2002).
project delivery in construction. Synergies within Barloworld. Unpublished manuscript,
Authors’ contacts:
while the measuring values and numbers are written results described can be achieved together with the author’s
upright. observations, theoretical statements or measurements.
10pt Preliminary communication contains one or more
Table 1 Table title aligned centre pieces of new scientific information, but without details that
(Style: Arial Narraw, 8pt, Align Centre)
allow recollection as in original scientific papers.
1 2 3 4 5 6
ABC ab ab ab ab ab ab Preliminary communication can give results of an
DEF cd cd cd cd cd cd experimental research, results of a shorter research or
GHI ef ef ef ef ef ef research in progress that is deemed useful for publishing.
10 pt Subject review contains a complete depiction of
Equations in the text are numbered with Arabic conditions and tendencies of a specific domain of theory,
numerals inside the round brackets on the right side of the technology or application. Articles in this category have an
text. Inside the text they are referred to with equation overview character with a critical review and evaluation.
number inside the round brackets i.e. “.... from Eq. (5) Cited literature must be complete enough to allow a good
follows ....” (Create equations with MathType Equation insight and comprehension of the depicted domain.
Editor - some examples are given below). Professional paper can contain a description of an
10pt original solution to a device, assembly or instrument,
1 t0 + t depiction of important practical solutions, and similar.
Favg (t , t0 ) =∫ F (q (τ ), p (τ ) ) dτ , (1) The article need not be related to the original research, but
t t0
10pt it should contains a contribution to an application of known
scientific results and their adaptation to practical needs, so it
α+β α -β
cos α + cos β = 2cos ⋅ cos . (2) presents a contribution to spreading knowledge, etc.
2 2 Outside the mentioned categorization, the Editorial
10pt board of the journal will publish articles of interesting
Variables that are used in equations and also in the text content in a special column. These articles provide
or tables of the article are formatted as italics in the same descriptions of practical implementation and solutions from
font size as the text. the area of production, experiences from device application,
10pt and similar.
10pt
3 WRITING AN ARTICLE
10pt
Article is written in the English language and the
terminology and the measurement system should be
adjusted to legal regulations, standards (ISO 80 000 series)
and the SI international system of units. The article should
Figure 2 The texts under figures be written in third person.
(Style: Arial Narraw, 8pt, Align Centre)
Introduction contains the depiction of the problem and
10pt
an account of important results that come from the articles
Figures and tables that are a part of the article have to
that are listed in the cited literature.
be mentioned inside the text and thus connected to the
Main section of the article can be divided into several
content i.e. „ … as shown in Fig. 1…“ or „data from Tab.
parts or chapters. Mathematical statements that obstruct the
1…“ and similar.
reading of the article should be avoided. Mathematical
10pt
statements that cannot be avoided can be written as one or
2 PRELIMINARY ANNOTATION
more addendums, when needed. It is recommended to use
10pt
an example when an experiment procedure, the use of the
Article that is offered for publication cannot be
work in a concrete situation or an algorithm of the
published beforehand, be it in the same or similar form, and
suggested method must be illustrated. In general, an
it cannot be offered at the same time to a different journal.
analysis should be experimentally confirmed.
Author or authors are solely responsible for the content of
Conclusion is a part of the article where the results are
the article and the authenticity of information and
being given and efficiency of the procedure used is
statements written in the article.
emphasized. Possible procedure and domain constraints
Articles that are accepted for publishing are classified
where the obtained results can be applied should be
into four categories: original scientific papers, preliminary
emphasized.
communications, subject reviews and professional papers.
10pt
Original scientific papers are articles that according to
4 RECAPITULATION ANNOTATION
the reviewer and the editorial board contain original
10pt
theoretical or practical results of research. These articles
In order for the articles to be formatted in the same
need to be written in such a way that based on the
manner as in this template, this document is recommended
information given, the experiment can be repeated and the
for use when writing the article. Finished articles written in
MS Word for Windows and formatted according to this