-□ □-
Методом планування експерименту отрима-на математична модель, яка описуе вплив про-цесу наповнення-пластифжацп на технологiчнi властивостi шкри. Виконана максимшна згорт-ка основних параметрiв якостi для переходу до однокритершльнои оптимiзацii. Розроблено обчис-лювальний модуль в середовищ^ Visual basic for Application для автоматизации розрахунку оптимального складу дисперсноИ системи методом Хука-Дживса. Отримат оптимальт значення концентрацш iнгредiентiв наповнювально-пласти-фжуючоН композицй i технологiчних властивостей шкряного матершлу
Ключовi слова: метод Хука-Дживса, оптимiза-ця процесу наповнення-пластифжацп, шкряний
матерiал, технологiчнi властивостi
□-□
Методом планирования эксперимента получена математическая модель, описывающая влияние процесса наполнения-пластификации на технологические свойства кожи. Выполнена мак-симинная свертка основных параметров качества для перехода к однокритериальной оптимизации. Разработан вычислительный модуль в среде Visual basic for Application для автоматизации расчета оптимального состава дисперсной системы методом Хука-Дживса. Получены оптимальные значения концентраций ингредиентов наполни-тельно-пластифицирующей композиции и технологических свойств кожевенного материала
Ключевые слова: метод Хука-Дживса, оптимизация процесса наполнения-пластификации, кожевенный материал, технологические свойства -□ □-
1. Introduction
Development of innovative technologies of the leather and fur skin materials formation requires application of a wide range of new synthetic and natural reactants [1, 2] in the filling-plasticizing process. The efficient implementation of the process requires creating conditions for a deep diffusion of chemical particles in a microfibrillar collagen structure of the dermis. This can be performed due to the high dispersion of plasticizing reactants of hydrophilic and hydrophobic nature to improve features of elastic polyfunctional leather materials [3, 4]. Taking into account the wide range of interactions both between components of a filling-plasticizing system and the system component interaction with collagen of semi-finished structured material during the tanning process, it is necessary to create conditions for the diffusion of reactants. In this connection, during first interaction minimal interaction of the system components is provided.
Optimization of the chemical composition of the above-mentioned systems is an important aspect of the problem. Usage of natural materials such as natural fats of animal or vegetable origin, natural aluminosilicates that are subject to modification depending on the stage of the technological process, is a long-term research area [4]. This is an important point both in terms of obtaining aggregation-stable colloidal
©
UDC 675.026.2
|DOI: 10.15587/1729-4061.2016.79459|
OPTIMIZATION OF LEATHER MATERIAL FILLING-PLASTICIZING PROCESS USING DISPERSE SYSTEM
A. Danylkovych
Doctor of Technical Sciences, Professor Department of Biotechnology, leather and fur Kyiv National University of Technologies and Design Nemirovich-Danchenko str., 2, Kyiv, Ukraine, 01011
Е-mail: [email protected] O. Sanginova PhD, Associate Professor* E-mail: [email protected] V. Chervinskyi Postgraduate Student* E-mail: [email protected] *Department of Cybernetic Chemical Technological Process National Technical University of Ukraine "Igor Sikorsky Kyiv Polytechnic Institute" Peremohy ave., 37, Kyiv, Ukraine, 03056
systems used in leather production technology and their diffusing capacity in the semi-finished leather.
The complex component composition of the filler-plasticiz-ing dispersed system and multiple stages of the structure and properties of the leather material formation process requires the use of multi-criteria optimization approach. Such an approach promotes the rational use of leather raw materials and chemical reagents. In this case, minimax optimization criteria convolution of filling-plasticizing process has been used. This method allows improving the worst indicator among the technological properties of the material. After such a convolution of criteria, single-criteria optimization problem has been solved by the method of Hooke-Jeeves. This method consists of two main stages - the exploring search and pattern search (accelerating search). Acceleration search allows finding the global optimum of the function briefly.
2. Literature review and problem statement
The analysis of publications devoted to the development of innovative technologies of leather and fur skin materials formation showed that the authors mainly use the of experiment planning methods to obtain a mathematical description of the mentioned processes. In [5] it is shown that
the simulation model and optimal area of the skin chrome tanning process were obtained by using Mathcad software where the noncomposite symmetric quasi-optimal plan was implemented. Mentioned software is not typically used at Ukrainian leather factories. Thus, such solution cannot be used to correct the optimum area of the process in operation mode under changing production conditions.
Plackett-Burman plan in the paper [6] is used to decrease the number of factors from 7 to 4. The survey of the full factorial experiment to study the process of chrome tanning with the use of modern polymer material - a derivative of maleic acid is shown in [7]. The adequate simulation model that describes the dependencies of semi-finished leather indicators on the consumption of basic materials and processing time is obtained. The rational conditions of tanning are identified. It should be noted that the Harrington's desirability function used for the formation of the generalized criterion of optimization is performed in the papers [6, 7]. With this approach, the formation of optimization criterion depends on the experience of the researcher [8]. The application of SWOT-analysis methods [9] in the marketing strategies area, connecting the factors of threats and opportunities, allows us to solve the problematic issues. However, it contains subjective evaluation of the company's management. It may adversely affect the activity of the enterprise in terms of dynamism and uncertainty.
The approach proposed in [5-8] to solve the optimization problem is associated with the complexity of the for-malization and requires a significant amount of calculations. In papers [10, 11] it is shown that in some cases the wrong selection of driving directions to the optimum does not lead to the desired result. The choice of the optimization method and the optimality criterion formation method have a decisive impact on the quality of the obtained solutions. Using the modified method of Hooke-Jeeves [12] allows to solve constrained optimization problems as well as it allows to optimize the investigated processes with a higher degree and achieve the optimal area more rapidly.
3. Purpose and objectives of the study
The purpose of this paper is an improvement of basic technological properties of leather material formed by fill-ing-plasticizing process optimization using the disperse system ingredients.
In accordance with the set goal, the following research tasks are identified:
- obtaining simulation models that adequately describe the effect of the filling-plasticizing process using dispersion system components on the technological properties of the leather material;
- development of software module allowing automation of the calculation of the required components concentration of the dispersed system;
- optimization of components concentration of the fill-ing-plasticizing process using the developed software module.
4. Materials and research methods of dispersed system
In accordance with the set purpose, in the optimization of the leather material-filling-plasticizing process, semifinished leather is used. It is obtained from cattle raw material using chrome tanning method and slicing the material
to 1.5-1.6 mm thickness. The finely aluminosilicate (FA) modified by sodium format is used as a filler of semi-finished leather. The dispersing mixture of natural fats and oils: beef tallow, sunflower oil, and fish oil in a ratio of 1:1.75:2.25 by weight % [13] is used for plasticizing the semi-finished product. The nonionic surface active agent is used to emulsify the fats.
The process of filling-plasticizing of semi-finished leather is carried out in the experimental section of PJSC Chinbar (Kiev, Ukraine) at the ratio of water:semi-finished product 1:1 and temperature of 55-60 °C. The concentrations of the plasticizer, surfactant and aluminosilicate filler were varied, respectively, at a rate: 49.6-140.4; 4.4-31.6; 15.0-30.0 (g/l).
The effectiveness of the filling-plasticizing process is determined by the concentration of the following substances which influence the technological properties of the materials, g/l:
- plasticizer (Xj);
- modified filler BAC (X2);
- nonionic surfactant (X3).
The technological effect of the ingredients of filling-plas-ticizing dispersion system application on the leather material obtained by the control technology was evaluated by the following quality parameters:
- output area of the semi-finished product, % (yi);
- the total content of fatty substances in the leather material, % (y2);
- hardness of material, cH (y3).
The output area of the leather material was calculated as the area ratio of the received prototypes of the finished material to the finished material obtained by the current technology. The total content of fatty substances in the finished leather was determined by the content of fat in the sample. The extraction in an organic solvent mixture after the hydrolytic destruction by alkali solution, followed by the extraction fatty acids with ether was used [14]. The hardness was measured by the PZU-12M equipment after drying and humidification processes. After that, the samples were cured in a desiccator at 20-22 °C and 16 % humidity for 24 hours.
The central rotatable second-order composite plan was designed to conduct the experiment. The plan is based on the full factorial experiment, complemented by experiments in the center of the plan and asterisk points [15, 16] which allow to evaluate the dispersion of reproducibility. The area of the experiment is described by the following model (1).
n n-1 n n
y = bo + Xb,x, +IIbljx1xj + ^b11x2, (1)
i=1 1=1 j=i+1 1=1 1*j
where xi, xj - factors, n - the total amount of factors bi, bij, bii - coefficients.
This model includes the first-order coefficients as well as pair and quadratic interaction effects. The results of the experiment are shown in Table 1.
Factors X1, x2, x3 take the zero value at the concentration of (g/l): plasticizer - 90, modified filler - 18, nonionic surfactant - 22.5 respectively, and their variation interval - 30; 8.1 and 4.5.
The transition from the technological values Xi to the encoded xi is carried out using the formula (2).
where Xi0 - factor zero value, AX¡ - variation interval.
Table 1
Central rotatable second-order composite plan with coded values of factors
The coefficients of the simulation model (1) were calculated using the least squares method in a matrix form (3)
№
10
11
12
13
14
15
16
17
18
19
20
xo
X1
1.68
1.68
0
X2
1.68
1.68
X3
1.68
1.68
0
X1-X2 X1-X3
1
0
1
0
X2-X3
1
0
x1
2.8224
2.8224
0
x2
0
2.8224
2.8224
0
x3
0
2.8224
2.8224
0
B = (xtx)-1 (xty).
(3)
Reproducibility dispersion using the experimental points in the center of the plan is defined as follows (4)
1
(4)
The value of the simulation model coefficients was tested by the Student's criteria and obtained by the formula (5)
tp_ >tT(Î0, q).
(5)
The criterion impact was considered significant if its estimated value exceeded a tabulated point. After the elimination of non-significant coefficients, an adequacy of the model was determined by the Fisher criterion. For this purpose two dispersions were calculated: the residual dispersion and dispersion of adequacy by the formulas (6) and (7).
S (yu - yj2
Si = ^
N -1
(6)
Note: the encoded of the plan; -1.68 levels of the plan
values -1, 0, 1 correspond to the lowest, zero and top-levels 1.68 - correspond to the asterisk points on the lowest and top
S2 -
ad _
S„ - S„
(7)
5. Simulation model of the filling-plasticizing process of semi-finished leather
According to the experimental plan, to study the efficiency of the effect of input variables on the technological properties of the leather material, 20 lots were collected. For this reason, the method of alternating halves [14] was used for semi-finished chrome-tanned leather, 3 pieces in each lot. The values obtained for the technological properties of semi-finished chrome-tanned leather after the implementation of the filling-plasticizing process are given in Table 2.
The Fisher criterion is calculated by the formula (8).
S2
F = -ad
c q2 • -
(8)
Table 2
Experiment Results
№ 1 2 3 4 5 6 7 8 9 10
Y1 105.6 98.7 103.8 96.5 101 95.1 98.4 93.2 104.7 93.7
Y2 7.3 5.9 6.7 3.9 5.6 4.8 5.1 3.7 6.7 4.7
Y3 24.5 28.5 23 29 22.5 24.5 32 37.5 25.5 31
Y1 104.3 100.9 99.3 93.6 105.9 105.2 105.1 104.3 105.6 105.4
Y2 7.8 5.7 4.3 2.8 7.5 6.4 7.2 6.1 7.4 6.9
Y3 17 19 33 43 17.5 20 18 22 19 21
The tabulated point of the Fisher criterion was taken according to the degree of freedom of dispersions of adequacy, dispersion and value level q, FT (fad, fo, q).
As the elimination of insignificant coefficients from the obtained model led to inadequate model FP>FT, the final model includes all the factors (Table 3).
In this case, the calculated value of the Fisher's criterion was less than the tabulated point FP<FT, for the specified degrees of freedom (Table 4). Thus, the resulting model is adequate.
Thus, the mathematical model adequately describes the process of filling-plasticizing of semi-finished leather. The following mathematical model (9) was used to optimize the composition of the dispersed system.
y1 = 105,2504 +3,208549 ^+1,041569 x2+ +1,940374x3 - 2,146144x2 -0,9414786 x2 -3,120475x2; y2 = 6,888165+0,7152907^+0 ,5663697x2+ +0,5218105x3 - 0,398291x2 - 1,160054x3; y3 = 19,61442 - 1,959721x1 - 1,821939x2 - 2 ,07405x3+ +2,9375x^3+2,87359x2 - 0,7580667x2 + 6,328101x2.
1
1
1
2
1
1
3
1
1
4
1
5
1
1
1
6
1
1
7
1
1
8
1
1
1
9
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
0
ad
Coefficients of the simulation model
bii bo b1 b2 b3 b23 bn b22 b33
Yi 105.2504 3.208548 1.041569 1.940373 -0.06251 -2.14613 -0.94148 -3.12048
Y2 6.911947 0.715291 0.56637 0.521811 0.1251 -0.40115 -0.02911 -1.16292
Y3 19.61442 -1.95972 -1.82194 -2.07405 2.9375 2.873593 -0.75807 6.328099
Note: the impact of factors on the technological properties of the material is insignificant
Table 4
Determination of the models adequacy
Model Y1 Y2 Y3
Dispersion
- residual; 0.12902 0.03702 2.63670
- adequacy 0.22691 0.19585 4.15753
Fisher criterion
- calculated; 0.4315 0.1165 0.8669
- tabulated 4.9503 4.7351 4.7351
Table 3 until the optimum is found with a given accuracy [17, 18].
As the obtained model includes three optimization parameters and each of them contains all three factors, the software module in the Visual Basic for Application was developed to automate the calculation. It allows to specify the starting point, step size, and the factors boundaries. The main window of the module is shown in Fig. 1.
A user places initial data on the MS Excel sheet (Fig. 2), and then can run the calculation of the simulation model (9).
As a consequence of solving the multicriteria optimization task using the software module, the factors optimal values of the leather material filling-plasticizing process were obtained. The 115.267, 22.492, 23.883 (g/l) concentration of plasticizer, modified filler and nonionic surfactant should be set respectively.
6. Optimization of the filler-plasticizer dispersion system composition
Using the simulation model (9) as an objective function specifies the solution of the multiobjective constrained optimization problem. Firstly, we need to obtain a compromise solution to solve this problem. Namely, to find values of the factors where the technological properties parameters of the material take the same values in the given range of values. The methods of the main criterion, additive, multiplicative convolution of criteria, generalized Harrington's desirability function [12, 17] can be used. In addition, the papers that describe the alternative compromise optimizing modern methods [18] are known. However, the mentioned methods are not applicable to the solution of chemical engineering problems.
In this paper, the maximin convolution of optimization criteria was used as it provides the best value of the worst criterion (10).
Fig. 1. Main window of software module
3 * □
1
Coefficients bO bl b2 b3 bl2 b23 bl3 bll b22 b33
yl 105,2504 3,208548 1,041569 1,940373 0 0 0 -2,14613 -0,94148 -3,12048
y- 6,911947 0,715291 0,56637 0,521811 0 0 0 -0,40115 0 -1,16292
y3 19,61442 -1,95972 -1,82194 -2,07405 0 0 2,9375 2,873593 -0,75807 6,328099
Y(x1,x2,x3) = max min yi(x1,x2,x3).
(10)
Fig. 2. Coefficients of simulation model in MSExcel
Thus, the maximum leather material output area y1=107.0393 %, the material hardness y3=18.986 sN and the total content of fatty substances in the material y2=7.565 % are provided.
According to the formula (10), the minimum value of each parameter yi is calculated, and the corresponding factors are determined. Optimum values of factors are defined as the maximum of the minimum value of the output variables. The method of Hooke-Jeeves, which is referred to the zero-order direct methods [12, 19], is used to solve this task. The Hooke-Jeeves method algorithm consists of two stages - the exploring search and the pattern search. The Hooke-Jeeves method specifies the starting point near the center of the test range. The starting point can be set as a random value. After that, the exploring search was conducted in the vicinity of the starting point, which allows to determine the direction of the function steepest ascent. Then, the accelerating search was carried out in the identified direction. If the function value ceased to increase, the search is conducted around the new point. If the exploring search failed in all directions and if the predetermined calculation accuracy is not reached, the step is narrowed down. The search is carried out
7. Analysis of the semi-finished leather filling-plasticizing process simulation model
Presented graphical diagrams (Fig. 3, a, b) demonstrate the technological properties of leather materials and the composition of the dispersed system relationships. Namely, the increase in the plasticizer concentration in the system leads to increase in the material output area up to the saturation point 115.27 g/l. Further increase of the plasticizer content decreased the output area. The reason for this is that the mobility of the collagen dermis structure fibers increased up to the critical plasticizer concentration in the material during its deformation in different directions. After water removing, the deformation degree of semi-finished leather structure is stabilised.
Influence of the modified filler and a nonionic surfactant concentration on the leather material output area is symbatical as compared to the plasticizer influence. However, the nonionic
surfactant provides the plasticizer and the filler particles diffusion to the product structure and their uniform distribution in the material by layers. Effective implementation of the process of filling- plasticizing of semi-finished leather can be carried out at the optimal concentration of surfactants in the disperse system. Insufficient amount of surfactants in the system does not provide effective plasticizer particles dispersion. However, the plasticizer overdosage increases the disperse system stability. In the former case, the plasticizer is deposited in the surface layer of the material. In the latter case, the plasticizer does not provide an active diffusion that leads to increased material shrinkage. Effect of filler on the material output area is due to its reduction resulted from thickness increase at a concentration above the critical value.
The major impact on the leather hardness is caused by nonionic surfactant concentration (Fig. 3, e, f). The increase of concentration from 15 to 24.95 g/l leads to an increase of material elasticity that is caused by aggregate stability increase of disperse systems and reduction of diffusion of the plasticizer and filler particles in the leather structure. In this connection, the increase of surfactant concentration above the critical value leads to a decrease of the plasticizing effect of fattening substances in the material structure. At the same time, the hardness of the leather material is increased.
To evaluate the influence of each factor on the output area of leather material W(xi), the single-factor dependencies for x1, x2, x3 respectively were obtained by the formula (11).
W = (b + b1JxJ)x1 + ba • x?.
(11)
For this purpose, each plot contains only analyzed factor index effects and the other factors are fixed at zero level. The diagrams of obtained dependencies (Fig. 4) have a similar nature.
Fig. 4. Dependency for ingredients concentration output area, g/l: 1 — plasticizer, 2 — filler, 3 — nonionic surfactants
e f
Fig. 3. Single-factor dependency for technological properties of leather material y! (a, b), y2 (c, d), y3 (e, f) in minimum and maximum concentration area, g/l: 1 — plasticizer, 2 — modified filler, 3 — nonionic surfactant: a, c, e — valley region; b, d, f — peak region
The concentration increase of the plasticizer from 49.6 to 119.05 g/l leads to an increase of the total content of fatty substances in the skin (Fig. 3, c, d). Further increase of fatty substances concentration in the system leads to a reduction of the total content of the plasticizer in the material. As a result, the degree of emulsification is reduced and the content in the disperse system is increased. An increase of the modified filler concentration due to the increased diffusion capacity of the hydrophilic filler leads to an increase of fatty substances content in the leather material.
Thus, the leather material output area is mainly affected by the concentration of plasticizer and a nonionic surfactant. The stabilization of the filler concentration at zero level provides the optimal area of the process (Fig. 5).
Fig. 5. Optimal area of semi-finished leather filling-plasticizing process
Fig. 5 shows that the optimal area of the filling-plasticizing process of semi-finished chrome-tanned leather produced from cattle raw material. The optimal area is located within the plas-ticizer and the nonionic surfactant concentration 84.2-123.1 and 21.8-26.6 g/l, respectively. The maximum output area of the leather material is achieved.
b
a
c
8. Discussion of the study results of the leather filling-plasticizing process
The obtained optimal area of the process of filling-plasti-cizing of the semi-finished leather allows to increase technological properties of the product, to decrease the stiffness of the leather material and to maximize the output area of the semi-finished product.
The output area of leather material with a minimum hardness of 19.0 sN is increased as compared to the control technology by 7.0 % using 115.3, 22.5, 23.9 concentrations in the studied solution (g/l) of plasticizer, modified filler and nonionic surfactant, respectively. Wherein fatty substances content in the finished material is 7.6 %, that complies with DSTU 311595 "Leather for Garments. General Specifications".
The developed software module can be used to successfully optimize technological processes of forming various composite materials.
9. Conclusions
1. The second-order mathematical model is obtained on the basis of experimental data at the PJSC Chinbar (Ukraine). The
model describes the semi-finished leather filling-plasticizing process adequately. The diagrams of single-factor dependencies for simulation models are obtained, and the influence of disperse system ingredients on basic technological properties of the leather material are analysed.
2. In order to convert multi-criteria task into single-criterion task using maximin convolution and Hooke-Jeeves method to solve the optimization task, a software module in the VBA environment (Visual Basic for Application) was developed. The developed module implements the considered computational schemes of a compromise solution. As a result, the optimum values of the factors that influence the semi-finished leather filling-plasticizing process conditions are found.
3. The optimum values of the concentration parameters of disperse system ingredients using multi-criteria constrained optimization of the filling- plasticizing process are found. The output area of leather material with a minimum hardness of 19.0 sN is increased as compared to the control technology by 7.0 % with 115.3, 22.5, 23.9 concentrations in solution (g/l) of plasticizer, modified filler and a nonionic surfactant, respectively. In addition, the content of fatty substances in the finished material is 7.6 % that meets the standard. The developed software module can successfully be applied for the optimization of technological processes of composite materials formation.
References
1. Danylkovych, A. G. Tekhnolohiya i materialy vyrobnytstva shkiry [Text] / A. G. Danylkovych, O. R. Mokrousova, O. A. Okhmat; A. G. Danylkovych (Ed.). - Kyiv: Feniks, 2009. - 580 p.
2. Grigor'ev, B. S. Himicheskie materialy i tehnologichsekie obrabotki pushno-mehovogo i ovchinno-shubnogo syr'ja [Text] / B. S. Grigor'ev. -Moscow: OOO «Al'kon Prestizh», 2004. - 48 p.
3. Innovacijni tehnologii virobnictva shkirjanih i hutrovih materialiv ta virobiv [Text] / A. G. Danylkovycha (Ed.). - Kyiv: Feniks, 2012. - 344 p.
4. Grishhenko, I. M. Polifunkcional'ni shkirjani materiali [Text] / I. M. Grishhenko, A. G. Danylkovych, O. R. Mokrousova; A. G. Danylkovych (Ed.). - Kyiv: Feniks, 2013. - 268 p.
5. Nikonova, A. Application of advanced polymeric compounds for development of leather production [Text] / A. Nikonova, O. Andreyeva, L. Maistrenko // IOP Conference Series: Materials Science and Engineering. - 2016. - Vol. 111, Issue 1. - P. 012024. doi: 10.1088/1757-899x/111/1/012024
6. Gomes, C. S. Optimizing adsorption parameters in tannery-dye-containing effluent treatment with leather shaving waste [Text] / C. S. Gomes, J. S. Piccin, M. Gutterres // Process Safety and Environmental Protection. - 2016. - Vol. 99. - P. 98-106. doi: 10.1016/j.psep.2015.10.013
7. Majstrenko, L. A. Doslidzhennya procesu xromovogo dublennya iz zastosuvannyam polimernogo materialu - poxidnoyi maleyinovoyi kysloty [Text] / L. A. Majstrenko, O. A. Andreyeva, M. K. Kolyada // Tehnologiyi ta dyzajn. - 2013. - Vol. 4.
8. Onem, E. Optimization of tannin isolation from acorn and application in leather processing [Text] / E. Onem, G. Gulumser, S. Akay, O. Yesil-Celiktas // Industrial Crops and Products. - 2014. - Vol. 53. - P. 16-22. doi: 10.1016/j.indcrop.2013.12.014
9. Solomenko, O. E. SWOT analysis as the tool of marketing strategies determination [Text] / O. E. Solomenko, Yu. V. Chernyuk // Visnyk KNUDT. - 2011. - Vol. 2, Issue 58. - P. 7-11.
10. Kuznetsov, A. S. Flexible Technology Options Assessment in the Case of Complex Criteria [Text] / A. S. Kuznetsov // Vestnik SIBGUTI. - 2012. - Vol. 1. - P. 13-22.
11. Lutsenko, I. Systems engineering of optimal control I. Synthesis of the structure of the technological product conversion system (part 1) [Text] / I. Lutsenko // Eastern-European Journal of Enterprise Technologies. - 2014. - Vol. 6, Issue 2 (72). - P. 28-37. doi: 10.15587/1729-4061.2014.28724
12. Branovitskaja, S. V. Application of Modified Hooke-Jeeves Method for Multi-Objective Optimization [Text] / S. V. Branovitskaja, S. G. Bondarenko, A. G. Danilkovich, O. V. Sanginova, V. A. Chervinsky // Naukivi visti NTUU "KPI". - 2014. - Vol. 1, Issue 26. - P. 24-30.
13. Patent na KM 23194, Ukraina, № u200613912 [Text] / Mokrousova O. R., Olijnik M. M. - Sklad dlja obrobki shkirjanogo napivfabrikatu. declareted: 27.12.2007; published: 10.05.2007, Bul. 6.
14. Golovteeva, A. A. Laboratornyj praktykum po khimii i tekhnologii ko i mecha [Text] / A. A. Golovteeva, D. A. Kucidi, L. B. Sankin. -Moscow: Legka i pishchevaya promyshlenost, 1982. - 312 p.
15. Adler, Ju. P. Planirovanie jeksperimenta pri poiske optimal'nyh uslovij [Text] / Ju. P. Adler, E. V. Markova, Ju. V. Granovskij. - Moscow: Nauka, 1986. - 278 p.
16. Montgomeri, D. K. Planirovanie jeksperimenta i analiz dannyh [Text] / D. K. Montgomeri. - Leningrad: Sudostroenie, 1980. - 384 p.
17. Ehrgott, M. Multicriteria optimization [Text] / M. Ehrgott. - Springer, 2005. - 328 p. doi: 10.1007/3-540-27659-9
18. Novakova, A. V. Metod kompromissnoj suboptimizacii dlya resheniya mnogokriterial'nyh zadach [Text] / A. V. Novakova, A. M. Davidenko, M. D. Kats // Eastern European Journal of Enterprise Technologies. - 2009. - Vol. 5, Issue 4 (41). - P. 39-44. - Available at: http:// journals.uran.ua/eejet/article/view/23963/21500
19. Hamdy, A. T. Operation research an introduction [Text] / A. T. Hamdy. - Prentice Hall, 2005. - 898 p.