CC BY
10
Valery F. Ochkov, Massimiliano Nori
Тheoretica[ questions of computer science, computational mathematics, computer science and cognitive information technologies
УДК 378.147
DOI: 10.25559/SITITO.14.201802.333-343
PLAYING WITH A CHAIN OR PHYSICAL AND MATHEMATICAL INFORMATICS
Valery F. Ochkov1, Massimiliano Nori2
1 National Research University Moscow Power Engineering Institute, Moscow, Russia
2 Saipem S.p.A., San Donato Milanese, Italy
КОМПЬЮТЕРНЫЕ ИГРЫ С ЦЕПОЧКОЙ
ИЛИ ФИЗИКО-МАТЕМАТИЧЕСКАЯ ИНФОРМАТИКА
В.Ф. Очков1, М. Нори2
1 Национальный исследовательский университет «МЭИ», г. Москва, Россия
2 Saipem S.p.A., г. Сан-Донато-Миланезе, Италия
333
© Ochkov V.F., Nori M., 2018
Keywords
Theoretical mechanics; sagging chain; chain function; derivative; integral; potential energy; kinetic energy; Lagrange-Dir-ichlet principle; Newton's law; constrained optimization; system of algebraic equations; computer graphics; animation; Mathcad.
Abstract
The article describes an educational laboratory work within the framework of interdisciplinary connections at the intersection of informatics, mathematics and physics: the study of the sagging of a closed chain with different support points. A technology of computer processing of photos and videos of a physical experiment with the subsequent processing of the media files on the computer is illustrated. We present the n chain number (1.258... - the optimal ratio of the chain length to the distance between its attachment points) and its relation to the previously unexplored problem of the shape of the sag of a closed chain hanging on two nails. The analytic expression for this constant was first found. A new physical and mathematical constant is found - the critical angle of the sagging of a closed chain on the "hangers". The problem of the sagging of a closed chain on a cone is detailed. The applicability of the computer tool "optimization with constraints" for solving problems of theoretical mechanics is investigated. Three main tools of mathematical computer tools are described: numerical mathematics, symbolic mathematics and graphics. The importance of using units of measurement when solving physical problems on a computer is emphasized. The problem of publishing mathematical formulas in articles and books is discussed.
Ключевые слова
Аннотация
Теоретическая механика; провисающая цепь; цепная функция; производная; интеграл; потенциальная энергия; кинетическая энергия; принцип Лагранжа-Дирихле; закон Ньютона; оптимизации с ограничениями; система алгебраических уравнений; компьютерная графика; анимация; Mathcad.
В статье описаны учебные лабораторные работы в русле междисциплинарных связей на стыке информатики, математики и физики: натурное и математическое исследование статического провисания замкнутой цепи с различными точками опоры, а также ее качание как маятника. Описана технология компьютерной обработки фото- и видеосъемки физического эксперимента с цепью с последующим обсчетом медиафайлов на компьютере. Обсуждается цепное число п (1.258... - оптимальное отношение длины цепи к расстоянию между точками ее крепления, находящимися на одном уровне) и его связь с ранее неисследованной задачей о двух формах провисания замкнутой цепи, подвешенной на двух гвоздях (точечных блоках без сцепления). На основе компьютерного вычислительного эксперимента найдена новая физико-математическая константа - критический угол провисания замкнутой цепочки на «плечиках». Детализируется известная задача о провисании замкнутой цепочки на прямом круговом конусе. Проанализированы подходы и сделана попытка численного решения задачи о маятнике-цепи. Исследуется применимость компьютерного инструмента «оптимизация с ограничениями» для решения задач теоретической механики. Дается описание трех основных инструментов математических компьютерных программ: численная математика, символьная математика и графика. Подчеркивается важность использования единиц измерения при решении на компьютере физических задач. Обсуждается проблема публикации математических формул в статьях и книгах.
About the authors:
Valery F. Ochkov, Doctor of Technical Sciences, Professor, National Research University Moscow Power Engineering Institute (14 Krasnokazarmennay St., Moscow 111250, Russia), ORCID: http://orcid.org/0000-0001-6844-5194, ochkov@twt.mpei.ac.ru
Massimiliano Nori, Ph.D. (Chemical Engineering), Process Engineer, Department PRCFA, Saipem S.p.A. (67 Via Martiri di Cefalonia, San Donato Milanese (Milano) 20097, Italy), ORCID: http://orcid.org/0000-0002-3898-9243, massimiliano.nori.home@gmail.com
Vol. 14, no 2. 2018 ISSN 2411-1473 sitito.cs.msu.ru
Modern Information Technologies and IT-Education
334 Теоретические вопросы информатики, прикладной математики, вф 0чков м Нои
компьютерных наук и когнитивно-информационных технологий
Introduction
Before spinners (turntables, twists) appeared, people who did not know what to do with their hands turned a coin, a box of matches or... a chain on their fingers.
Let's also turn the chain, but not on the fingers, but on ... two nails, on "hangers", on a cone ... And we will do this not for idleness, but for scientific and cognitive purposes. In passing, we will discover a couple of new physical and mathematical constants.
Theoretical, implementation and experimental parts
Usually articles of this kind are clearly divided into theoretical, implementation and experimental parts. But in this article it's impossible to do it - they mixed up with each other The author hopes that the reader will understand why this happened, and will take this not quite an ordinary technique.
1. The chain on two nails
In [1] it is described how a chain is suspended at two ends (see Figure 1) and how it is all photographed and sent to the computer, where it is processed appropriately. The fingers of two schoolchildren were used as "nails". In this school class (excellent laboratory work for two classes - physical and computer) at the junction of mathematics, physics and computer science, it was once again shown that the chain is not sagging along the parabola (and this is a very common misconception, and not only for schoolchildren and students), but on the catenary, whose formula [2, 3] will be used in the calculations below.
chain problem.
Take a closed chain1 and hang it not on a finger, but on two nails (on two pins), driven into the wall2 - see Fig. 2. How will it sag? The problem is extremely idealized: the chain itself is "classical" (absolutely flexible and inextensible), and the frictional forces of the chain around the nails are absent. The first condition is easy to be met in an application, while the second is almost impossible, given the fact that the chain with its joints at the attachment points clings fairly tightly to the nails. For all this construction it is necessary, for example, to tap, so that it will come to the expected equilibrium. On nails it is possible to fasten a block with bearing and then to throw a chain on it. This will complicate the geometry of the problem, but it will allow calculations to be compared with an application in the same way as was done in [1].
Fig. 1. Experiment with the sagging chain Рис. 1. Эксперимент с провисающей цепью
Let us consider some modifications that have not yet been investigated for a classical problem of the calculus of variations, the sagging
Fig. 2. Real closed chain suspended on two nails Рис. 2. Реальная замкнутая цепочка, подвешенная на две булавки
1 Circular in the sense that it has neither beginning nor end. Such chains without locks are hung around the neck, if the length of the chain allows it to pass through the head. In the third part of the article, we'll talk about putting the chain on a "neck without a head" - on a cone. By the way, about the locks. They can be regarded as certain point weights on the chain, complicating its calculation and making the task more interesting (the problem of a pendant on the chain, a cabin of the cableway). By "circularity" we can also mean the absolute smoothness of the chain - the absence of friction at the points of its attachment.
2 A sheet of millimeter (graph) paper is attached to the wall, with which you can measure the coordinates of individual points of the sagging chain and compare them with those that were obtained as a result of calculation. Such a work is described in [1]. In our computer age millimeter paper is almost out of use. The author recently raked the rubble in his office and came across rolls of millimeter paper - ordinary, semi-logarithmic and logarithmic in both coordinates. There was an idea to paste it in place of wall paper in the room, where there is the laboratory works with the chain. But this idea of mathematical design of the premises is not new: the floor in the lobby of the headquarters of the Mathematical Association of America (IAA) in Washington is lined with tiles in the form of pentagonal parquet. A search and classification of polygonal parquets is a visual and interesting task of the theory of paving modern combinatorial geometry.
Современные информационные технологии и ИТ-образование
Том 14 № 2 (2018) ISSN 2411-1473 sitito.cs.msu.ru
Тheoretical questions of computer science, computational mathematics, V ry F. Ochkov' M™liano Nori computer science and cognitive information technologies
ЭЭБ
Figure 3 shows the beginning of calculation of the sagging of the closed chain: input of the initial data3 - length of the chain (60 cm) and coordinates of the attachment points of the nails, on which the closed chain is thrown (0-28 and 18-15 cm).
60cm Length of the circular chain
x-| := Ocm y^ := 28cm Left chain fastening point X2 18cm y2 - 15cm Right chain fastening point
such texts are not placed on paper but in electronic form, then they can be "live": the reader is allowed to change the original data and receive a new answer. The example in Fig. 4: author's settlement site for Mathcad Calculation Server technology, with which you can calculate the length of the chain by specifying the coordinates of the points of attachment of its ends h and h2, and also its clearance - the minimum distance from the "floor" h.
¡4>) jg) http://twt.mpei.3c.ru/MCS/Workshtets/chain-S-Eng.xmcd
y(x,h,a,xg) :- h + a-J^cosh^--j - ij Catenary
d x ~ *o"l
y'|x,a,xg) :=—y|x,h,a,xg| -» sinhj-j Derivative ofthe catenary
m„-= 70 — Linear mass of the chain qf -= ^ c cm 1000
A - Jj'x1 - + [y-| - y^ - 22.2cm Distance between the points
Calculation ofthe catenary y(x,x0,h,a):=h + a ^cosh^—— J-1
Ч *\Ш
Linear density of chain
kg/m
PE(hD.hU.aD.aU.x0D.x0u)
V|) <
Уи*
I ,/1 + (y'(x.aD.xœ))2dx Jïi x0 ■= 0-1° Im "" |0-35 I m -|0.10 L "
Л Min. Chain length I := Л2 + (h-, - h2)2 = 0.18m
y(x,hD,aD,x0D)-/l + (y,(x,aD, x0D))2 dx J*1 Given CTOL:- 10" n - 1 X 10" 7
LD Fixing at the left end h1 = y(0,XQ,h,a)
rX2 ,-- Fixing on the right end h2 = y(L-x0-h >a)
y(x,hu,au,xou)4|1 + (y'(x,au,xou)) dx [V ' 0.09 '
J»1 s := Find(xg,S,a) = 0.5187867 m
l-ld I a J 1^0.0324103 J
Chain length
-^y(x,K0,h,a)j
Э-^Уо-Ц) +ïu-(L"lc)]
Fig. 3. Beginning of the Mathcad calculation of the sagging of a closed chain on two
nails
Рис. 3. Начало Mathcad-расчета провисания замкнутой цепочки на двух гвоздях
Important note. The authors do not give individual formulas in this article, but place in it the texts of the solution of problems in the physical and mathematical program Mathcad. At the same time, the authors perfectly understand that this prevents some readers, accustomed to traditional publications with separate numbered formulas, to understand the essence of the tasks being solved. What can you say about this!? First, all the formulas of this article are fairly well read in the texts because Mathcad uses general mathematical notation employing WYSIWIG technology (What You See Is What You Get). Secondly, this style of publication meets the modern world trend, where formulas in books and articles are typed not in the text (mathematical) editors of LaTeX, Word Equation, MathML, etc., but in Maple, Mathematica, Mathcad, etc. By such "living" formulas, you can consider, it is permissible for them to build graphics and create animations, which significantly reduces the likelihood of errors and typos in them. Moreover, more and more often the chapters of books and journal articles are printouts of solving problems in mathematical programs. If
0 ' ^^^^^^^^
Fig. 4. Online-calculation of the sagging chain Рис. 4. Online-расчет провисающей цепи
Figure 4 shows the calculation of the length of the lower part of the sagging closed chain - 51.9 cm. If in this online calculation the value of h is changed from 5.1 to 27.1 cm, the length of the upper short part of the chain is calculated to be equal to 18.1 cm (one millimetre larger than the distance between the attachment points L). In total, this gives 70 cm - this is approximately the length of the chain shown in Fig. 2 after it was removed from the nails and its length was measured straight (using a ruler), and not by an indirect (calculated) method. So the theory is verified by practice!
Figure 3 also shows the user functions - the formula of the catenary4 (y)
3 As a rule, such problems are solved in dimensionless quantities. But the Mathcad package has a unit of measurement tool and it's a sin not to use it. This, in particular, will avoid mistakes of incorrect input of formulas. The user gram-force unit (gf) is introduced into the calculation - one thousandth of the built-in kilogram-force unit (kgf). Our calculations are tuned to a system of physical quantities cm-g-s, so the default values of forces will be output in dynes, which we will change to gram-force.
4 The canonical formula of the catenary, used in paper and electronic mathematical reference books and textbooks, has the form a-cosh(x / a). But in the calculations we will use a non-canonical form of the catenary with two additional parameters h and x0: h + a • cosh(((x-x0) / a) - 1). In this formula, a specific point of this curve is clearly defined - the catenary has a minimum (sagging chain - a > 0) or maximum (arch - a < 0) at the point with coordinates h - x0, which will need to be found in the course of our calculations.
Vol. 14, no 2. 2018 ISSN 2411-1473 sitito.cs.msu.ru
Modern Information Technologies and IT-Education
336
Теоретические вопросы информатики, прикладной математики,
В.Ф. Очков, М. Нори
компьютерных наук и когнитивно-информационных технологий
and its derivative5 (y'). And another function called РЕ with six arguments returns the value of the potential energy sagging on two nails of the closed chain, depending on the parameters of the two segments of the catenary, on which the chain is divided - see Fig. 2. We will conditionally denote these sites by indices as D (down, lower) and U (up, upper). But this division is conventional - the lower part of the chain can appear at the top, and the upper one at the bottom.
The law of mechanics (Lagrange principle [4]) states that the chain sags so that its potential energy becomes minimal. Here looms a typical optimization problem with constraints, that can be solved with the Mathcad built-in function Minimize, along with the Given keyword -see Fig. 5.
Г1(Г
Ни 15
aD 3
аи 1
X0D Ли ^ 10 J0,
11 + У
x.aD.x0D))2<ix+ I + (y'(x,au,*ou))2
J*.
ï(K1.hD.aD.x0Dj=ï1
v(xi'hu'au'x0u} = vi
y(x2,hD,aD,x0D) = y2 y(x2, hu.au, xou) = y2
itov
Да.
mUV
Доо, l Дои.у
:= Minimize ( И hQ ц -p *0D ■ "Oil) =
14.64 1
7.356 10.406 4.403 15.27 10.672 J
4d i y-l + (y'fx, aD,x0D})2dx = 24.13 cm Ly 5
.я.
L - Ln - 35.87 cm
However we should always remember that this answer is approximate, depending on the value of the given accuracy of the numerical solution of the problem.
The value m is the linear mass of the chain. It does not affect the
c
shape of its chain, but it will help us, first, to check the correctness of the calculation shown in Fig. 5 from the point of view of the balance of forces, and secondly, to calculate the forces that stretch the chain in correspondence of its attachment on the two nails. This new problem can be reduced to the numerical solution of an overdetermined system of eight nonlinear algebraic equations with six unknowns describing the "mechanics" of a closed chain sagging on two nails with a minimum potential energy. The calculation uses the Mathcad built-in function Find, which returns the values of unknown equations that convert them into identities with accuracy determined by the method and by the parameters of the numerical solution of the problem.
y1D
~y2D =y1U ~y2U
xD FxU , Giver
Fy1D+ Fy2D= 9 " mc " LD Fy1D=-V,(x1,aD,x0D) ■ FxD
Ру2О=У'(х2-аО-х0о)-рхО
Fy1U + Fy2U=9"mc- LU Fyiu = -y,(x1.au-x(]u) - FxU Fy2U = y'lx2-aU-x0Ul ' FxU
2 2 Fy1D + FxD = F'
y1U
2 c 2 + FxU
y2D
2 2 2 2 " ' fxD = Fy2U + FxU
Fig. 5. Minimizing the potential energy of a hanging closed chain Рис. 5. Минимизация потенциальной энергии висящей замкнутой цепочки
The following calculation steps are executed in Fig. 5:
- input of the initial assumption for the numerical solution of the minimization problem;
- entering constraints after the word Given; these are: the length L of the chain remains constant, and the two sections of the closed chain -the lower (D) and the upper (U) are fixed at the ends at two given points (1 and 2);
- execution of the Minimize function, that returned as answers: the parameters a, h and x0 of the two parts of the catenary (D and U) so that their potential energy PE is minimal and the constraints are satisfied. This answer (the vector of constants) can be transferred to the initial assumptions and assures that the answer does not change.
~y1D "y2D ~y1U ~y2U fxD
rxU
:= F"d(Fy1 D. Fy2D ■ Fy1 U ■ Fy2U • FxD • FxU) =
( 1.54365^ 0.11098 1.74238 0.80299 0.83 945 0.30399
■9f
F1:=2;/Fy1D2+FxD2=3-S14-3f
2 :=
: г/р.лп" + F,
y2D + fxD = 1-694-gf
Fig. 6. Calculation of tension forces of a chain on two nails Рис. 6. Расчет сил натяжения цепочки на двух гвоздях
Physical laws, shown in the calculation in Fig. 6, such as: the sum of the values of the vertical projections of the forces of fastening the two parts of the chain on two nails are equal to the values of the weight of the two parts of the chain; The name variable Fy1D means that it is the vertical (y) projection of the force (F) with which the lower (D) part of the chain pulls the left (1) nail down; the values of the horizontal projections of the forces of fastening the
5 The formula of the derivative of the catenary function is generated by the Mathcad package itself using the symbolic mathematics tool (right arrow operator). Formulas for the length of this curve and for the ordinates of its center of gravity that use a certain integral and are also involved in calculations have been found on the Internet, although they can also be easily derived using the Mathcad symbolic mathematics or "in the mind", based on our own knowledge of the basics of mathematical analysis.
Современные информационные технологии и ИТ-образование
Том 14 № 2 (2018) ISSN 2411-1473 sitito.cs.msu.ru
Valery F. Ochkov, Massimiliano Nori
Тheoretica[ questions of computer science, computational mathematics, computer science and cognitive information technologies
337
parts of the chain on the two nails (FxD and FxU) are equal to each other and are related to the values of the vertical projections of these forces through the values of the derivative of the catenary at these two
points: y'(x^
, y'(x2,
, y'(Xi<
n) and y'(x2, an, x0U);
the tension forces of the lower and upper parts of the chain at the points of their fastening are equal: the chain parts do not overstretch each other and are in equilibrium.
The calculation shown in Fig. 6, is verified. We were convinced that the values of the parameters of the catenary (the variables aD, ay, hD, hy, x0D and x0U), found through minimization of the potential energy, correspond to the values of the forces ensuring the equilibrium conditions of our mechanical system.
But what did the numerical computer experiment show us?! The fact that the closed chain, depending on the location of the two nails can sag in two ways - see Figures 7 and 8. Case 1. The chain is divided into two unequal sections - see Fig. 7. Case 2. The chain does not bifurcate and sags in two equal halves - see Fig. 8.
\ \
\ s
\ L
\ \ LU 1 LU -
\ \
\ \
\ N
V 4 L
\ 4 и
\ 4 4
\
\ J
\ /
\
\ /
\ /
s 1 n
4 /
_2Л = 44-4Q7cm 1 1 1 1 1 1
01234S67B9 101112131415161718
Fig. 7. The closed chain hanging on two nails is divided into two unequal sections Рис. 7. Провисающая на двух гвоздях замкнутая цепочка разделяется на два неравных участка
Fig. 8. The closed chain hanging on two nails is divided into two equal sections Рис. 8. Провисающая на двух гвоздях замкнутая цепочка разделяется на две половинки
At the bottom of Figures 7 and 8, the graphical dependences of the potential energy value of two parts of the closed chain are shown depending on the ratio of their lengths - Ly/L. The following special points are visible on the graphs:
- Fig. 7 - two minima (stable equilibrium) and one maximum (unstable equilibrium - the chain can in principle be of this form with two equal halves (see Figure 8), but at the slightest external disturbance it will slip into the left or right potential wells);
- Fig. 8 - one minimum (stable equilibrium - the chain sags in two equal parts).
Figure 9 shows a Mathcad document that defines a user function that returns the value of the potential energy (PE) of the chain sagging on two nails, depending on the ratio of the length of one part to the total length of the chain (LU/L). In this calculation case the Mathcad package built-in function Find does not return a constant vector as in Fig. 6, but a vector of user-defined functions - dependent on the parameters of the two parts of the catenary, as function of the variable (argument of the function) LU/L. These vector functions named AnsD and
Vol. 14, no 2. 2018 ISSN 2411-1473 sitito.cs.msu.ru
Modern Information Technologies and IT-Education
338
Теоретические вопросы информатики, прикладной математики, Вф Очков М Нои
компьютерных наук и когнитивно-информационных технологий
AnsU are decomposed into separate scalar functions with names hD, aD, etc. This technique allowed us to create a function named PE and then build its graphs, depending on the ratio of the length of the LU to the length L for different coordinates of the points of attachment of the closed chain - see Figures 7 and 8.
yfx1,hD,aD,x0D| = y1 y(x2,hD,aD.x0D| = y2 I J1 + jVjx.aD,x0D|\2 dx =
AnsD(LU/L) := Find(hD,aD,x0D|
hD(LU/L) := AnsD[LU/L)0 aD(LU/L) := AnsD(LU/L)-| x0D(LU/L) := AnsD(LU/L)2
y(x,, hy, аи ,хш) = y, y(x2, hy, ay .xou) = y2 f 2 J^f^yj^f dx = L - L-LU/L
1
AnsU(LU/L) := Find(hy,ay,xou|
hy(LU/L) := AnsU(LU/L)0 ay(LU/L) := AnsU(LU/L)-| xou(LU/L) := AnsU(LU/L)2
|'K2 -
PE[LU/L) := y^.h^LU/L),aD{LU/L),xod{LIMJ)-J 1 + (y(x,aD{LU/L),x0D{LU/L)))2dx ...
I
y(x- hy[LU/L), ay[LU/L) ,XQy[LU/L)|-^ 1 + |y"|x, ay[LU/L) ,XQy[LU/L)jj5 dx
:
Fig. 9. The dependence of the potential energy of the two parts of the sagging closed chain on the ratio of their lengths6 Рис. 9. Создание зависимости потенциальной энергии двух участков провисающей замкнутой цепочки от отношения их длин
To find the boundary separating the two forms of sagging of a closed chain on two nails, we will be helped ... by an interlude.
Interlude. The chain number n
Problem. How much longer should the chain (non-closed, conventional with two ends) be than the distance between the two points of its suspension, which are on the same level, so that the tension at the ends of the chain is minimal?
This is a typical problem of finding a minimum: if you increase the length of the chain, then this force will enhance due to the growth of its vertical projection due to the increase in the weight of the chain. If the chain is shortened, the force of its tension at the attachment points will increase due to the growth of its horizontal projection. When the value of the chain length approaches the value of the distance between the points of its attachment, it ceases to be just a chain and becomes ... a tightly stretched string. The calculation of the string, its vibrations - is a separate very interesting physical and mathematical problem, where the chain, pardon, string cannot be considered inextensible. Here it is necessary to take into account the modulus of elasticity of the string material and other parameters. Figure 10 shows the solution of this problem using Mathcad. To solve the problem, we need to relate the force F and the length of the chain L, i.e. to define a function F(L) - the optimization objective function. In the calculation in Figure 10, two user functions are routinely defined: a catenary whose minimum / maximum falls on the Y axis (the X axis passes through the two points of the chain attachments), and its derivative with respect to x. The problem is reduced to the solution of a system of two non-linear algebraic equations (the
equality of the heights of the attachment points of the chain and the invariance of the value of the chain length] using the Given-Find block, where the Find function returns not two concrete numbers (the numerical solution of the system of two equations], but defines a user function with name Ans (see also this technique in Figure 9]. Then this function-vector (vector because the function returns two numbers] is decomposed into two separate functions a(L) and h(L). Next, the optimization objective function F(L) is formed, along which the graph is plotted and in which the minimum is specified - the point where the derivative value is zero.
y{x .a rh) h+ a■■ COSil - -1 tstenery
yix.a):- ^-yOt.a.hí-^slnhf - | catwuiyi dx
X := OSm
Input Dal«
, , mc-g L ^tiieal fof^je or Uie fuspítnJed chsin « netion of ifc FW(U ---_- length
2
a := 1m
fl := -0.2m Fifit approximation
Given
y{—X ,3 , h) = y(x , Э ,hj = Om Height Equality at the altadiment points
■X _
L= ! 7l+yfe.a)2dx
Ans(Ú Findfa.h)
Chain length
ail) := AnsíOo h(ü Ans(Ll1 Ci(l) atan(f(x,al(l})) Chain angle at Hie right aüitíimerit point
fv(l)
fiúy
Optimization Otyflctiva Function
Y Y V
L := 2-X+ — ,2-Хт- — + —.. 4X 2D 20 300
Un, := root 2-Х+ — .4X = 1.258r
idL 20" J
M
y(MLwiM(L5i*))
0 2
•-9Л
-0.2
2 I
a(Lop,)-56.166 dea
Fig. 10. Determination of the optimal chain length Рис. 10. Определение оптимальной длины цепи
6 In the formula for the potential energy, in Figure 9, the constant mc • g (the linear weight of the chain material - see Figure 5] was removed and other simplifications were made that accelerate the calculation. However it's still quite long. The author, on his rather fast computer, drew the energy curves shown in Fig. 7 and 8, in several tens of minutes. The calculation of the coordinates of each point of this curve requires the numerical solution of a rather complex system of nonlinear algebraic equations with high accuracy. And this task itself is very time-consuming in addition to the other resources of the computer.
Современные информационные технологии и ИТ-образование
Том 14 № 2 (2018) ISSN 2411-1473 sitito.cs.msu.ru
Theoretical questions of computer science, computational mathematics, va[ery R Ochkov, Massimi[iano Nor, computer science and cognitive information technologies
339
The number 1.258 ... (the ratio of the length of the optimal chain to the distance, between the points of its attachment at the same level) can be regarded as some new physical and mathematical constant and give it the name nc (chain number n - catenary) - the ratio of the arc length of the optimal catenary to its "diameter" - the distance between the attachment points. The usual number n is, as everyone knows, the ratio of the arc (the whole circle) to its diameter. It is possible to start a competition - to calculate the maximum number of digits in this constant, similar to how it is done with respect to the "circular" number n: 3.142.
Discussion of this problem with its analytical solution can be found here: https://community.ptc.com/t5/PTC-Mathcad-Questions/Sym-bolic-solution-of-one-optimization-problem/td-p/130281. The author of the analytical solution is Luc Meekes from Netherland - see the Fig. 11.
СЭ- 7 Ф http://www.wolframalpha.com/input/?i=sqr(3)%2B%CE%B6(3)%2F3-7%2F8_- б] Поиск...
sqr(3K(3)/3-7/8
Input:
3 8
Decimal approximation:
1.257736441955408721994025728676355697197800684590546921986...
Fig. 11. Symbolical calculation of the optimal chain length in WolframAlpha site Рис. 11. Символический расчет оптимальной длины цепи на сайте WolframAlpha
The function Z(x)' is the Riemann Zeta function. You might want to check if, with sufficient numerical accuracy, your nc as approximated, reaches this number, or differs from it.
Hence a practical advice. If you want to hang a chain, rope or cable between two posts or walls at the same height, then, in the absence of any special requirements, it will be sufficient to choose a chain length that is greater by about a quarter of the distance between the attachment points, which are naturally on the same level. In this case ... see above. If the points of attachment of the chain are at different levels, then this constant will become different (and be no longer a constant, but a variable instead).
After the publication in the author's book [5] of this research and the emergence of the hope for the discovery of a new physico-mathemat-ical constant, a search was made for the Internet on the key "1.258 catenary", which led to the publication "Optimal form of a sagging cable" [6], in which already appeared the number 1.258. So we are a little late for the priority on this constant, but we consider ourselves to be his co-author. (End of the interlude)
So that's it! Let's sketch a chain on two nails as shown in Fig. 7, and we will gradually reduce its length. A numerical experiment shows that at the moment of the fusion of two sections of the chain (Figure 8), the force applied to the nails acquires a minimum value. Further, this force will only increase. This is a certain boundary between the two forms of sagging of the chain shown in Fig. 7 and 8. A closed chain can be rotated not only on the fingers, but also in a different way: by hanging it on the left and right index fingers and rotating them around the circumference. Figure 12 shows the animation frames of this hand exercise in its mathematical implementation: one point of attachment of the chain is fixed, while the second is rotating in a circle around it. Rather, only one quarter of a circle, but on the other three quarters the result will be identical (symmetrical).
IF Воспроизв
"M1?)"1)
h(j = 12.063067cm aj = 3.287B6 cm xou = 12 413484cm hD= 12.061245cm aD = 3.286013cm x0D = 12.413907cm
IE Воспроизв
hu = 12.731035cm аи = 2.058159 cm xou= 11 181875cm hp = 13.656799 cm aD = 2.982916 cm x0D= 10 419552cm
IE Воспроиз
□
X-X0
H
hu = 13.154966 cm аи= 0 776121cm xou = 7.939224cm hD - 12.936273cm aD = 0.210544 cm x0D = 9.523557cm
Ю Воспроизвести аним... - □ X
■ 1 ГХ_Хэ) Ï
h+ a- cosh - - 1
V 1 a J J
hu=.
зи= .
xou= ■
hD=.
30= .
X0D= ■
► II 9
Fig. 12. Frames of animation of rotation of a closed chain between two nails: the top row is a long chain, the bottom row is a short chain. Рис. 12. Кадры анимации вращения замкнутой цепочки между двумя пальцами, пардон, гвоздями: верхний ряд - длинная цепочка,
нижний ряд укороченная цепочка
Vol. 14, no 2. 2018 ISSN 2411-1473 sitito.cs.msu.ru
Modern Information Technologies and IT-Education
340 Теоретические вопросы информатики, прикладной математики, ВФ Очков М Нои
компьютерных наук и когнитивно-информационных технологий
What conclusions can be drawn by analysing Fig. 12? Firstly, the distance between the points of suspension of the chain is not the only criterion determining the form of its sagging. This criterion is unique for the case when these points are on the same level -see Fig. 10.
Secondly, two frames of animation on the right edge of Fig. 12 once again underline the limits of the numerical solution of the problem: under certain initial conditions, there is no such solution. This is connected with the very essence of numerical methods and with the fact that a vertically sagging chain requires a different formula for calculation, rather than the one shown in the animation frames in Fig. 12.
2. Chain on "hangers"
A hanger, human shoulders designed, is in the form of an obtuse-angled isosceles triangle with a hook at the top. It is hung up with clothes in the wardrobe7.
Some things (women's dresses, for example) often slip off from these "shoulders (hangers)" and fall down. Or vice versa, the straps are assembled at the bottom of the hook - at the top of the triangle. Which, too, is not good in terms of keeping things safe and undisturbed. So our research will have not only purely scientific, but also some applied value! Let's see how this hanger will hold not lingerie underwear, but our closed chain.
Problem. Take the already known closed chain, which is not hung on two nails (Figure 2), but on two straight segments ("on the hangers" -see Figure 13). The forces of friction (a stiction), as in the previous problem, are not taken into account, and the chain is inextensible and completely flexible. How will this chain hang?
Figure 14 shows the protocol for the numerical solution of this problem in Mathcad. The already familiar Minimize function is used, minimizing the potential energy of the closed chain under certain restrictions. The task is similar to the previous one, but the contact points of the chain are not fixed at two points, but "slide" along the segments of the inclined lines - along the "hangers".
It is assumed that the chain on the "hangers" can also hang in an asymmetrical position (obliquely - "a crooked smile"). Therefore, two unknown quantities xL and xR are entered in the calculation: the distances from the vertical axis where the hook is located to the left and right suspension points.
y| x. h. a.xq i := h + a ■ ■ cosh
j x - x0 ^ y'|x. a .xq| := —y{x.h.a.xq| —» sinh.- .
PE{x[_.xpj. hp. hy. sq. ay .xqq .XQyi := _ y(x. hD. aD.x0D)-J 1 + y'(x. aD .x0D;2 dx ... \
" , y(x,hy.ay.xou:,-/fx.ay.xoy^dx \
(*L XR hD hU aD aU xOO X01J := 2 2.5 2 5 7 -1 1 )m Guess values
Given
f"R - f"R -
+ y-(x.aD.x0D;2dx+ J1 - y|x.ay.dx= L
■V -V
y|xL.hD.aD.x0D;,= xLtan:|-.v
( Tt
y|xR.hD.aD.x0D;, = -xRtan:--.i
y|xL.hy.ay.x0y;,= xL tan: | - aj y(xR.hy.ay.xou) =-xR tan: | - |3|
XL *R
aD aU X0D ,K0U
:= Mmimize(PE.xL.xR.hD.hy.aD.ay.x0D.x0u;,=
-1.66502 1.66501 -3.35281 -3.35269 3.03114 3.03188 -0.00002 -0.00002
Fig. 13. Closed chain on hangers Рис. 13. Замкнутая цепочка на "плечиках"
Fig. 14. Calculation of the position of the closed chain on the "hangers" Рис. 14. Расчет положения замкнутой цепочки на «плечиках»
The chain position shown in the photo in Fig. 13, is not an equilibrium one. It is an artificial, determined by the friction force (a stiction) of the chain on the "hangers". Analogy: the bar on the inclined plane (see the centre of Figure 14) is stationary, although it must slide downward, reducing its potential energy. The cause of immobility is the force of friction. And not the frictional force of sliding, which is calculated by a simple formula, but the frictional force of rest, the calculation of which is much more complicated. If these forces would not exist, the chain would either sag with two identical halves symmetrical about the Y axis, or slip off the hangers and sagged vertically,
7 Conditionally all hotels in the world can be divided into two different groups. In the closets of some of them (and these, alas, are a minority) hang ordinary "home" hangers. In others, one might say, not "homemade" but "wild" hangers, which consist of two parts - the "hook", which is tightly tied to the rod in the closet, and actually the hanger itself without a hook, that is necessary to cunningly combine with the hook. This is very inconvenient, but it is done to counter the theft of hangers by hotel guests. The author saw in one hotel also a third type. There the hangers were attached to the closet ... by a chain - the object of our investigation. In the old days, in public place, a chain often attached a mug to a pot of drinking water. Now it is often possible to see in offices a chain fastened to a pen for writing so that visitors do not accidentally or intentionally take it away.
Современные информационные технологии и ИТ-образование
Том 14 № 2 (2018) ISSN 2411-1473 sitito.cs.msu.ru
Theoretical questions of computer science, computational mathematics, Va1ery R Ochkov, Massimi1iano Nori computer science and cognitive information technologies
341
catching the base of the hook of the hanger. And this "either" does not depend on the length of the chain, but is only related to the angle of disclosure of the "hangers". Here there is a qualitative difference between the "hangers" (Figure 13) and the two nails (Figure 2), on which two parts of the closed chain can sag either in two unequal sections (see Figure 7) or symmetrically (Figure 8) to achieve the minimum potential energy. The chain on the "hangers" always hangs in half. It begins, however, to divide into two unequal parts at the moment when the slip begins before slipping off the "hangers".
With the help of a numerical experiment the critical opening angle of the "hangers", in which the chain jumps off the hanger, was determined. To determine this angle, an animation was created, two frames of which are shown in Fig. 15. The animation itself can be seen here https://community.ptc.com/t5/PTC-Mathcad-Questions/Round-chain-on-the-corner/td-p/553624 (site with discussion of the task). With an increase in the angle a (the angle of the "opening" of the hangers) by an amount greater than 50 angular degrees, the chain jumps off the "hangers" and hangs on the hook attachment point - at the apex of the obtuse angle.
■ ' Воспроизвести анимацию — □ X n + л= 100е (\ = 50 " .1=50 7 " ' Воспроизвести анимацию — □ X n + |Э= 103° (\ = 51.5-: .i=51.5-:
► i | ► I |
Fig. 15. Two frames of the animation of the closed chain jump from the hangers Рис. 15. Два кадра анимации соскакивания замкнутой цепочки с плечиков
ей
Catenary y(x,h,a) ?= h + a-^cosh:
L := 7m mr:=70— G := gmr-L = 4.805 N
c m c
a := 1m Guess values
Length, linear mass and weight of chain
y^x, a) := —y(x. h, a) sinh dx
Sver
X = 1m Y — -1m h := -1 rr
tan{a) = — y(X. h. а) = Y
С* ; /
J-x
1 + У(х.а) dx=- y,[*.a) = tai(a)
F[o0 s Fnd(X,Y,h,a) X(o) - F(a)o Y(a) s F[a)i h(o) - F(n)s a[.Y> := FW] rfJ.n;
y(x. h(a). aM) ~|Ï7ytTIiÔ? dx
PEc(ni -
-W
PEV :=
0.5 L
-0.5 L
- - root(PEc(.\} - РЕу.о.ЗМед.ЛМед) - 50.34-d&g Answer
°рг - ~~~
PEcln] .--
«V
Fig. 16. Calculation of the critical angle of a closed chain sagging on the hangers Рис. 16. Расчет критического угла провисания замкнутой цепочки на плечиках
Figure 16 shows the calculation of this critical angle-the numerical determination of the value of a, at which the potential energy of a chain sagging along two identical arcs of a chain line (PEc), is equal to the potential energy of a chain hanging vertically (PEv) and hooked to the top of the "hangers"8. This angle acr was slightly more than 50 degrees, which correlates with the data obtained when viewing the animation (Figure 15).
Most likely, an analytical rather than a numerical approach to solve the chain problem on the "hangers" will give an answer not in angular degrees, but in the form of a fraction n-rc / m, where n and m are natural numbers greater than unity. Mathematicians are always somewhat upset when it is said that a sine of 60 angular degrees9 is equal to 0.866 ... The orthodox mathematician here will say differently - see the second line of calculation in Fig. 17.
You can try to find analytically (symbolically) the value of acr - the values of n and m in the expression n-rc / m, but you can do it in the way shown in Fig. 18 - sorting out the values of n and m in the range from 1 to 30, and remembering those (n = 7, m = 25), which most correspond to the numerical response in Fig. 16.
8 First the catenary parameters X, Y (right attachment point), h (clearance) and a (catenary scaling parameter) are determined as function of the angle a through 4 equations. Then the critical angle is found by solving the equation equating the potential energy of two identical arcs of a chain line (PEc) to the potential energy of a chain hanging vertically (PEv). Conditionally by combining the equation of the arc length and the condition on the derivative we obtain the so-called Whewell equation of the catenary (i.e. L/2=a tan(a)).
9 An air balloon burst from the clouds. From the basket the flyer shouts to a man on the ground: "Excuse me, where are we?" Answer from the man: "You are in a basket of a balloon!" He was a mathematician. Only from a mathematician can we hear an absolutely accurate but completely useless answer. There is a little-known continuation of this old joke. The flyer in response shouts: "You do not understand! Our navigator's batteries are depleted, and we cannot determine our coordinates! ". The mathematician looks at his smartphone and answers: "Zero point ninety-six hundredths radians of northern latitude and zero sixty-four hundredth radians of eastern longitude!" Only mathematicians measure angles in radians, and not in more familiar degrees!
Vol. 14, no 2. 2018 ISSN 2411-1473 sitito.cs.msu.ru
Modern Information Technologies and IT-Education
Теоретические вопросы информатики, прикладной математики, 342 , - В.Ф. Очков, М. Нори
компьютерных наук и когнитивно-информационных технологий
sin (Б0П) = 0.366
sKi)
41 2
Fig. 17. Numerical and symbolic calculation of the sine in Mathcad Рис. 17. Численный и символьный расчет синуса в среде Mathcad
for п = 1.. 50 for ill .г 1 Ы1
■а
if А с Л, А,
50.34deg--тт
m
mm
Л
mm п 1 п ml m
a closed chain can be investigated using the program shown in Fig. 14. The result in Figure 19.
At any angle of inclination of the "hangers", a closed chain wrapped around them does not bifurcate, the angle of connection of the ends of the chain with the "hangers" remains straight, which eliminates the sliding of the chain along the "hangers".
<\- .*•= 75D
t\ - ■ 4l.:>
— тт = 50.4deq 25
Fig. 18. Refinement of the value of the critical angle of the opening of the "hangers" with the chain
Рис. 18. Уточнение значения критического угла раскрытия «плечиков» с цепочкой
The calculation shown in Fig. 18, of course, can be considered a kind of curiosity - an imitation of an analytical approach to the solution of the problem. But the authors emphasize this approach not in vain, since it is very common in our computer age. In some cases, it is very productive, while in others it is counterproductive. We entrust those readers who are strong in analytical transformations, to confirm or refute the assertion that a = 7n / 25.
cr '
The reader can also replace rectilinear hangers by round ones or in the form of a segment of a parabola, a hyperbola or even a segment of a chain line (a catenary) and analyze the behaviour of a closed chain on such "designer hangers." On the chain, you can hang a pendant (point mass) and evaluate the properties of such a mechanical system.
3. Chain on a cone
The behaviour of a closed chain on a circular, straight cone is interesting. The reader can either conduct such experiment by himself, or read the report on this work in [7], the acquaintance with this report moved authors to write this divertissement.
A few words about the chain on the cone. It may, depending on the angle of the opening of the cone, either jump off it, or remain on the cone (the frictional forces are again not taken into account). But there is also a third equilibrium - an unstable one, which is qualitatively displayed on the lower graph of Figure 7. The chain at the slightest external jerk can either jump off the cone, or take a stable equilibrium, enveloping the cone. This was proved by the authors of [4]. It is also interesting to investigate the oblique cone, too, and how a closed chain will behave on it. By the way, indirect - "oblique hangers" with
Fig. 19. Chain on "oblique hangers" Рис. 19. Цепочка на «кривых плечиках»
Conclusions
A chain could not only be worn around neck or twisted on the fingers - it is also possible to conduct interesting physical and mathematical experiments with it - real and on the computer! For educational and scientific purposes.
References
[1] Ochkov V.F., ed. Catenary = physics + mathematics + informatics. Informatics in school. 2018; 3(136): 56-63. (In Russian)
[2] Merkin D.R. Introduction to the mechanics of a flexible thread. M.: Science, 1980. 240 p. (In Russian)
[3] Ochkov V.F., Ochkova N.A. A project of a monument to three mathematicians or MathMetry. Cloud of Science. 2017; 4(4):548-571. Available at: https://elibrary.ru/item.as-p?id=32603280 (accessed 15.04.2018). (In Russian)
[4] Ochkov V.F., Lener F., Chudova Yu.V., Kapitonets V.K., Tarakano-va D.Yu. Physics vs informatics: rope polygon with weights in static, kinematics and dynamics Or Newton vs Lagrange. Cloud of Science. 2017; 4(2):303-336. Available at: https://elibrary. ru/item.asp?id=29461211 (accessed 15.04.2018). (In Russian)
[5] Ochkov V.F., Bogomolova E.P., Ivanov D.A. Physical and mathematical studies with Mathcad and the Internet. Second edition. SPb: Publishing house "Lan", 2018. 560 p. (In Russian)
[6] Wang C.Y. The optimum spanning catenary cable. European Journal of Physics. 2015; 36(2). DOI: 10.1088/01430807/36/2/028001
[7] Zubelevich O.E., Samsonov V.A. Chain on the cone. Collection of scientific and methodical articles "Theoretical mechanics". Issue 30. Publishing house of Moscow State University. 2018, pp. 131-138 (In Russian)
[8] Ochkov V.F., Kolhepp F. Physic and IT: centre of gravity of a black box. Informatics in school. 2017; 7(130):65-70. (In Russian)
[9] Ochkov V.F., Bogomolova E.P. Teaching Mathematics with Mathematical Software. Journal of Humanistic Mathematics. 2015; 5(1):265-285. DOI: 10.5642/jhummath.201501.15.
Современные информационные технологии и ИТ-образование
Том 14 № 2 (2018) ISSN 2411-1473 sitito.cs.msu.ru
Theoretical questions of computer science, computational mathematics, Va[ery F. Ochkov, Massimi[iano Nori computer science and cognitive information technologies
[10] Solodov A., Ochkov V. Differential Models. An Introduction with Mathcad. Springer-Verlag Berlin Heidelberg, 2005. 232 p. DOI: 10.1007/b138129
[11] Korobov V., Ochkov V. Chemical Kinetics with Mathcad and Maple. Springer-Verlag Wien, 2011. 344 p. DOI: 10.1007/9783-7091-0531-3
[12] Behroozi F. A fresh look at the catenary. European Journal of Physics. 2014; 35(5). DOI: 10.1088/0143-0807/35/5/055007
[13] Wang C.Y, Watson L.T. The elastic catenary. International Journal of Mechanical Sciences. 1982; 24(6):349-357. DOI: 10.1016/0020-7403(82)90068-6
[14] Chen J.-S., Li H.-C., Ro W.-C. Slip-through of a heavy elastica on point supports. International Journal of Solids and Structures. 2010; 47(2):261-268. DOI: 10.1016/j.ijsolstr.2009.09.031
[15] Boresi A.P., Schmidt R.J. Engineering Mechanics - Statics. Baker & Taylor Books, North Carolina. 2000. 704 p.
[16] Christensen H.D. Analysis of simply supported elastic beam columns with large deflections. Journal of Aerospace Science. 1962; 29:1112-1121.
[17] Lippmann H., Mahrenholtz O., Johnson W. The heavy elastic strips at large deflections. International Journal of Mechanical Sciences. 1961; 2:294-310.
[18] Parbery R.D. The effect of stiffness on the shape of the elastic catenary. Civil Engineering Transaction. 1976; 18:98-101.
[19] Wang C.Y. A critical review of the heavy elastica. International Journal of Mechanical Sciences. 1986; 28(8):549-559. DOI: 10.1016/0020-7403(86)90052-4
Submitted 15.04.2018; revised 30.05.2018; published online 30.06.2018.
Список использованных источников
[1] ОчковВ.Ф. и др. Цепная линия: физика + математика + информатика // Информатика в школе. 2018. №3(136). С. 56-63.
[2] Меркин Д.Р. Введение в механику гибкой нити. М.: Наука, 1980. 240 с.
[3] Очков В.Ф., Очкова Н.А. Проект памятника трем математикам, или Матметрия // Cloud of Science. 2017. Том 4, № 4. С. 548-571. URL: https://elibrary.ru/item.as-p?id=32603280 (дата обращения: 15.04.2018).
[4] Очков В.Ф., Ленер Ф, Чудова Ю.В., Капитонец В.К., Тараканова Д.Ю. Физика vs информатика: веревочный многоугольник в статике, кинематике и динамике, или Ньютон vs Лагранж // Cloud of Science. 2017. Том 4, № 2. С. 303-336. URL: https://elibrary.ru/item.asp?id=29461211 (дата обращения: 15.04.2018).
[5] Очков В.Ф., Богомолова Е.П., Иванов Д.А. Физико-математические этюды с Mathcad и Интернет. 2 е изд., испр. и доп. СПб.: Издательство «Лань», 2017. 560 с.
[6] Wang C.Y. The optimum spanning catenary cable. European Journal of Physics. 2015; 36(2). DOI: 10.1088/01430807/36/2/028001
[7] Зубелевич О., Самсонов В. Цепь на конусе. Сборник научно-методических статей. Теоретическая механика. Том 30. М.: Изд-во Моск. ун-та, 2018. С. 131-138. .
[8] Очков В.Ф., Кольхепп Ф. Физика и информатика: центр тяжести черного ящика // Информатика в школе. 2017. № 7(130). С. 65-70.
[9] Ochkov V.F., Bogomolova E.P. Teaching Mathematics with Mathematical Software // Journal of Humanistic Mathematics. 2015. Vol. 5, issue 1. Pp. 265-285. DOI: 10.5642/jhum-math.201501.15
[10] Solodov A., Ochkov V. Differential Models. An Introduction with Mathcad. Springer-Verlag Berlin Heidelberg, 2005. 232 p. DOI: 10.1007/b138129
[11] Korobov V., Ochkov V. Chemical Kinetics with Mathcad and Maple. Springer-Verlag Wien, 2011. 344 p. DOI: 10.1007/978-37091-0531-3
[12] Behroozi F. A fresh look at the catenary // European Journal of Physics. 2014. Vol. 35, issue 5. DOI: 10.1088/01430807/35/5/055007
[13] Wang C.Y., Watson L.T. The elastic catenary // International Journal of Mechanical Sciences. 1982.Vol. 24, issue 6. Pp. 349357. DOI: 10.1016/0020-7403(82)90068-6
[14] Chen J.-S., Li H.-C., Ro W.-C. Slip-through of a heavy elastica on point supports // International Journal of Solids and Structures. 2010. Vol. 47, issue 2. Pp. 261-268. DOI: 10.1016/j.ijsol-str.2009.09.031
[15] Boresi A.P., Schmidt R.J. Engineering Mechanics - Statics. Baker & Taylor Books, North Carolina. 2000. 704 p.
[16] Christensen H.D. Analysis of simply supported elastic beam columns with large deflections // Journal of Aerospace Science. 1962. Vol. 29. Pp. 1112-1121.
[17] Lippmann H., Mahrenholtz O., Johnson W. The heavy elastic strips at large deflections // International Journal of Mechanical Sciences. 1961. Vol. 2. Pp. 294-310.
[18] Parbery R.D. The effect of stiffness on the shape of the elastic catenary // Civil Engineering Transaction. 1976. Vol. 18. Pp. 98-101.
[19] Wang C.Y A critical review of the heavy elastica // International Journal of Mechanical Sciences. 1986. Vol. 28, issue 8. Pp. 549-559. DOI: 10.1016/0020-7403(86)90052-4
Поступила 15.04.201; принята в печать 30.05.2018; опубликована онлайн 30.06.2018.
|об авторах:|
Очков Валерий Фёдорович, доктор технических наук, профессор, Национальный исследовательский университет «МЭИ» (111250, Россия, г. Москва, ул. Красноказарменная, д. 14), ORCID: http://orcid.org/0000-0001-6844-5194, ochkov@twt.mpei.ac.ru
Нори Массимилиано, доктор философии по химическому машиностроению, инженер-технолог, Saipem S.p.A. (20097, Италия, г. Сан-Донато-Миланезе, ул. Мартири ди Каталония, д. 67), ORCID: http://orcid.org/0000-0002-3898-9243, massimiliano.nori.home@gmail.com
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted reuse, distribution, and reproduction in any medium provided the original work is properly cited.
Vol. 14, no 2. 2018 ISSN 2411-1473 sitito.cs.msu.ru
Modern Information Technologies and IT-Education
