CC BY
33
2015
НАУЧНЫЙ ВЕСТНИК МГТУ ГА
№ 217
УДК 66-5
ВОПРОСЫ ОПТИМИЗАЦИИ РАЗМЕРОВ ЭЛАСТИЧНЫХ РЕЗЕРВУАРОВ
ДЛЯ ПОЛЕВЫХ СКЛАДОВ ГОРЮЧЕГО
Ю.Н. РЫБАКОВ, С.В. ЛАРИОНОВ, С.И. ЧИРИКОВ, Р.И. КЮННАП
В статье изложены теоретические подходы к вопросам оптимизации геометрических параметров эластичных резервуаров (ЭР) для нефтепродуктов с учётом общей длины швов на верхней и нижней поверхностях. Предложены алгоритм, математические модели и критерий оптимизации, позволяющие существенно (до 20%) повысить эффективность раскроя конструкционного материала. На основе полученных зависимостей приведены формулы для определения оптимальных высоты, длины и ширины эластичных резервуаров. Приведены расчеты оптимальных размеров на примере эластичного резервуара вместимостью 250 м3.
Ключевые слова: эластичные резервуары, оптимизация геометрических размеров, метод Виета-Кардано, критерий оптимизации, конструкционный материал, эффективность.
Вопросам определения оптимальных параметров резервуаров посвящено значительное количество работ различных авторов. Классическая постановка и метод решения задачи нахождения оптимальных размеров цилиндрических резервуаров принадлежит В.Г. Шухову [1]. В дальнейшем различными авторами были рассмотрены и решены аналогичные задачи как для случая постоянной, так и переменной толщины стенки резервуара. Большинство полученных решений относилось к вопросу определения оптимальных параметров стальных цилиндрических резервуаров. Основной проблемой в этих работах являлось нахождение отношения высоты резервуара к его диаметру.
В случае нахождения оптимальных параметров эластичных резервуаров для полевых складов горючего вместимостью 300 - 15 000 м3 (рис. 1) возникает аналогичная задача. Также необходимо определить оптимальные соотношения между геометрическими размерами резервуара. Однако, если в случае цилиндрической формы резервуара задача нахождения его оптимальных размеров сводилась лишь к определению отношения между его высотой и диаметром, то для эластичного резервуара эта задача является более сложной.
б
Рис. 1. Полевой склад горючего, вместимостью 6000 м : а - участок хранения полевого склада горючего; б - раскрой и сварка оболочки эластичного резервуара ЭР-250
Эластичный резервуар [2-5] представляет собой геометрическое тело сложной формы. Для эффективного решения задачи нахождения его оптимальных размеров необходимо аппроксимировать его форму, представив ее в виде более простых геометрических тел. Для этого представим эластичный резервуар в виде параллелепипеда с присоединенными к каждой его боковой грани полуцилиндрическими поверхностями. Таким образом, объем рассматриваемого тела, аппроксимирующего эластичный резервуар, будет складываться из объема параллелепипеда, объемов цилиндров с высотами, равными длине и ширине
а
параллелепипеда, и радиусами, равными половине высоты параллелепипеда, а также объемов четырех цилиндров с тем же радиусом и высотами, равными высоте параллелепипеда. Таким образом, объем эластичного резервуара равен
У = Уи + УЦ1+Уц2+4Уц3,
где УП — объем параллелепипеда УП = 2Я-Ь-Н; Уц1,Уц2,Уц3 — объемы цилиндров: Уц1 = жЯ2 -Ь; Уц2 = жЯ2-Н; Уц3 = жЯ2-2Я = 2жЯ3, Ь - длина эластичного резервуара; Н - его ширина; Я -радиус цилиндрической поверхности, равный половине высоты резервуара.
Отношение между шириной и длиной параллелепипеда будем считать заданной,
определяемой особенностью эксплуатации, обозначим ее ^ = кО
Тогда
Обозначим отношение длины и ширины эластичного резервуара как ^ = к, значит Н = кЬ.
Уп = 2Я-Ь -кЬ = кЬ2-2Я; Уц2 = жЯ2-кЬ.
В итоге объем тела, аппроксимирующего эластичный резервуар, равен:
У = к\}-2Я + жЯ2-1 + жЯ2-кЬ + 4-2жЯ3;
У = к12-2Я + жЯ2-Ь-(1 + к) + 8жЯ3. (1)
Для дальнейшего решения задачи оптимизации размеров ЭР выразим длину эластичного резервуара через его объем.
Приведя выражение для объема тела (1), аппроксимирующего эластичный резервуар, к виду уравнения относительно длины Ь, получим
]}-2кЯ + Ь-жЯ2(1 + к) + (8жЯ3 — У) = 0. (2)
Решая (2) относительно Ь, получим
_ —жЯ2(1 + к)+ V8кЯУ + л2Я4(1 + к)2 — 64жкЯ4
1= 4Ш "
Обозначим 8кЯУ + ж2Я4(1 + к)2 — 64жкЯ4 = Б*, Э* > 0, тогда
т -жЯ2(1+к)+4о^
Ь =--—--. (3)
Величиной, минимизирующей стоимость эластичного резервуара, является минимум площади поверхности фигуры, аппроксимирующей эластичный резервуар. Ее величина равна
5 = 5П + 5ц1 + Бц2 + 46цз,
2
где 5П = 2ЬН = 2Ь-кЬ = 2кЪ 5ц1 = 2жЯ-Ь; Бц2 = 2жЯ-Н = 2жЯ-кЬ; Бц3 = 2жЯ-2Я = 4жЯ2.
Подставив БП, $ц1,5ц2,5ц3 в общую формулу площади поверхности, получим:
5 = 2кЬ2 + 2жЯ-Ь + 2жЯ-кЬ + 4-4жЯ2;
5 = 2кЬ2 + 2жЯЬ(1 + к) + 16жЯ2. (4)
В соответствии с поставленной целью основным критерием оптимизации будет
Б^тт. (5)
Для выполнения условия (5) сначала подставим в уравнение (4) найденное значение Ь из (3), затем найдем производную площади поверхности по радиусу цилиндров (—) и приравняем ее к
йЯ
нулю.
Итак, 5 = Ь2-2к + Ь-2—К(1 + к) + 16пИ2;
£ _ -—Я2(1+к)+
= 4кЯ '
где й* = Я-8кУ + И4-[л2(1 + к)2 _ 64—к];
= 16к2Я2 = 8к2Я ' ( )
Подставляя найденные Ь и Ь2, получим
п2Я2(1+к)2 + я(1+к)^Р* + 2
4Я 4 к 4 к
V „2 (л Л т?(1+к)2\ 1 1 V , —1+к)
или 5 = К2-(14—---—) + К
5 = К2.(14П _ «М2) + R-1.1L + ф*)1 (7)
V 4 к / 4 4 к
Найдем производную по Я
или
^ = 2К.(14п _ *(1+к)2) _ К-2 .V +1
йЯ V 4к ) 4 2 йЯ 4к
й^_2/? (^1471 — —2(1+к)2\__й£* —1+к)
йЯ = V 4к ) 4Я2 йЯ
2 йЯ 8кЛг' Производная Б* по Я равна
НП*
— = 1-8кУ + 4К3-[—2(1 + к)2 _ 64лк] = 8кУ + 4—2К3(1 + к)2 _ 256—К2
йЯ
Тогда
йБ „ л2(1+к)2\ V , л(1+к)[2к1+л2Я3(1+к)2-64лкЯ3]
— = К( 28л---—— )--- + ——, гА (8)
йЯ V 2к ) 4Я2 2кл18кЯ1+л2Я4(1+к)2-64лкЯ4 У
2к^8кЯ1+л2Я4(1+к)2-64лкЯ4 ' л2(1+к)2 \ V л(1+к)[2кУ+тЛ2Я3(1+к)2-64лкЯ3]
Приравниваем производную к нулю
„ (л2(1+к)2 __ \ , V л(1+к)[2к1+л2Я3(1+к)2-64лкЯ3]
К(——-_28л)+ — = ——, А (9)
V 2к ) 4Я2 2к^81кЯУ+^Я4(1+к:)2-64—кЯ4 У
Преобразовав последнее выражение, получим пусть л2 + 2л2к + л2к2 _ 56— = а
2—1+к№Я^-8—Я5'] = V8кКУ + _ 8лк)
Упростим выражение и получим
4К9[1792—3к3 + 32п2к2(35а + 16) + 16лка(а + 20) _ а3] + +4К6[896л2к2 + 8—(а _ 1) _ а _ 6а2]кУ + К3[232— _ 29а]к2У2 _ 8к3У3 = 0. Если представить К3 = Р, то получим кубическое уравнение вида
4?3[1792л3к3 + 32л2к2(35а + 16) + 16лка(а + 20) _ а3] + +4Р2[896л2к2 + 8—(а _ 1) _ а _ 6а2]кУ + Р[232лк _ 29а]к2У2 _ 8к3У3 = 0
или
Р3-А + Р2-В + РС + в = 0, (10)
где А = 4(1792л3к3 + 32л2к2(35а + 16) + 16лка(а + 20) _ а3); В = 4(896л2к2 + 8—(а _ 1) _ а _ 6а2)кУ; С = (232— _ 29а)к2У2;
в = -8к3У3.
Из решения уравнения (10) относительно Я3 = Р будет определяться оптимальное значение И и далее по формуле (6) можно определить L - длину эластичного резервуара и Н= кЬ - его ширину.
Расчет стоимости эластичного резервуара с учетом швов
Введем следующие переменные:
2 2 а - стоимость 1 м эластичного материала, руб./м ;
в - стоимость 1 м сварного шва эластичного материала, руб./м;
I - общая длина швов во всем резервуаре;
£ - площадь поверхности резервуара;
С - общая стоимость резервуара с учетом швов.
Критерий оптимизации
С = — тт. (11)
Количество швов на верхней (нижней) поверхности резервуара —,
где Н - ширина резервуара;
к - расстояние между швами.
ь
А
V
к
пЯ ~2
н
пЯ ~2
Н
к = т
Н = кЬ
Общая длина швов на верхней и нижней поверхностях резервуара (с учетом боковой поверхности длиной L)
„т Н+пЯ
2Ь--.
К
Общая длина швов на боковой поверхности длиной Н (по ширине)
2пЯ •
К'
где пЯ - длина одного шва на боковой поверхности. Тогда общая длина всех швов резервуара
1 г Н+ЯИ , « х-. Н ,
1 = 21--+ 2пЯ • - . (12)
К К
Общий вид целевой функции с учетом стоимости швов
С = (*• (2кЬ2 + 2кЯЬ(1 + к) + 16кИ2) + ¡3 • (2Ь • + 2пЯ • (13)
Для выполнения условия (11) найдем производную функции (13) стоимости резервуара по
ЛС
радиусу цилиндров (—) и приравняем ее к нулю.
йС йБ , „ й1
— = а---+ В • — . (14)
ая ая ^ ая к '
Так как — была рассчитана в пункте 4, то посчитаем отдельно —.
Преобразуем I = — + - • (1 + к) и вспомним, что по формуле (4)
5 = 1? -2к + 1-2пЯ(1 + к) + 16кИ2. Следовательно, I = -- 16пК =-(5- 16пИ2).
ННН
И = -(?1-32пК). (15)
йИ НКйК ) у '
Подставим результаты (10) и (15) в формулу (14) и приравняем производную к нулю:
ас ¿я / в\ в
йР. йЯ V К К
аь /аь \ и
— •а = 1—-32пР)^.
ая ш ) к
Упростив последнее выражение, получим
4Я9[х2(3х + 2,5Рпк - 64пк) + (х - 0,5Рпк + 64пк)(х + 0,5Рпк)(3х - 0,5Рпк + 64пк)] + +4И6[4х(3х + 0,5Рпк) + (х - 0,5Рпк + 64пк)(1 -2х- Рпк)]кУ + +И3[35х + 2,5Рпк - 64пк]к2У2 - 8к3У3 = 0. Если представить И3 = Р, то получим кубическое уравнение вида
Р3-А' + Р2В' + РС' + С = 0, (16)
где А = 4(х2(3х + 2,5Рпк - 64пк) + (х - 0,5Рпк + 64пк)(х + 0,5Рпк)(3х - 0,5Рпк + 64пк)); В' = 4(4х(3х + 0,5Рпк) + (х - 0,5Рпк + 64пк)(1 - 2х - Рпк))кУ; С = (35х + 2,5Рпк - 64пк)к2У2; С = -8к3У3.
Пример расчета заданных объемов резервуаров и их цены с учетом швов
«-» 3
Найдем оптимальные размеры эластичного резервуара объемом 250 м , где отношение длины и
2 2 ширины к=0,75. Стоимость 1 м эластичного материала а = 10 000 руб./м . Стоимость 1 м шва
эластичного материала в = 3000 руб./м.
Расстояние между швами И = 1,5 м
2ва+% 28-10000 + 4-3™
■л __К ___ го о/:
Р = 3000 = 52,36
10 000 --пг"
К 1,5
Найдем значения коэффициентов:
А' =1310691,00; В' = -5205783,52; С' = -5740443,37; С = -52734375,00. Подставим данные коэффициенты в уравнение (16) и найдем решение по методу Виета-Кардано [6]
1310691 х 3 - 5205783 х 2 - 5740443 х - 52734375 = 0.
Коэффициенты:
а = -3.9717851118227; Ь = -4.3797073452095; с = -40.234025411024; _ (а2 - 3Ъ) _ ((-3.9717851118227)2 - 3 • (-4,3797073452095)) _
У = 9 = 9 = 3,2127>
(2а3 - 9аЪ + 27с)
И =--—-- = -25,3368.
54
Следовательно, по методу Виета-Кардано уравнение имеет один действительный корень (общий случай) или два (вырожденный). Кроме действительного корня имеются два комплексно-сопряженных:
x i = 5.88; x 2 = -0.954 + i х (2.436); x з = -0.954 - i х (2.436).
Учитываем только действительный корень Р = 5,88.
R = 0,87 м; 2R = 1,74 м; L = 11,09 м; Я = 9,92 м; 5 = 257,08 м2.
Длина всех швов резервуара I = 62,81 м.
Минимальная стоимость резервуара заданного объема
С = а • 5 + $ • I = 10 000 • 257,08 + 3 000 • 62,81 = 2 759 230 руб.
Выводы
1. Разработана математическая модель для определения оптимальных размеров ЭР, базирующаяся на минимизации площади его поверхности. На основе полученных зависимостей приведены формулы для определения оптимальных соотношений между высотой, длиной и шириной ЭР.
2. Приведена модель для определения оптимальных размеров ЭР с учетом стоимости его швов. Показано, что учет стоимости швов ЭР оказывает влияние на величины их оптимальных размеров.
3. Получены формулы для определения оптимальных соотношений между геометрическими параметрами ЭР, учитывающие влияние стоимости его швов. Рассчитан пример определения оптимальных размеров ЭР объемом 250 м2.
ЛИТЕРАТУРА
1. Шухов В.Г. Искусство конструкции / Райнер Грефе, Оттмар Перчи, Ф.В. Шухов, М.М. Гаппоев и др. М.: Мир, 1994. 192 с.
2. Патент РФ 2304553. Мобильный эластичный резервуар для нефтепродуктов / Рыбаков Ю.Н., Харламова О.Д., Паталах И.И., Абрамов М.Д. // Заявл. 07.04.2006. Опубл. 20.08.2007. Бюл. № 23.
3. О принятии на снабжение полевого склада горючего ПСГ-600 вместимостью 600 м3: приказ Министра обороны Российской Федерации № 132 от 02.03.2015.
4. ТУ 7981-023-08151164-2013. Полевой склад горючего ПСГ-1000. Технические условия. М.: ФАУ «25 ГосНИИ химмотологии Минобороны России», 2013.
5. Исследование по созданию опытно-промышленного производства передвижных эластичных резервуаров: отчет о НИР. М.: ФАУ «25 ГосНИИ Минобороны России», 2010.
6. Бронштейн И.Н., Семендяев К.А. Справочник по математике. 7-е изд., стереотипное. М.: Государственное издательство технико-теоретической литературы, 1967. С. 138-139.
PROBLEMS OF SIZE OPTIMIZATION OF FLEXIBLE TANKS FOR FUEL FARMS
Ribakov Yu.N., Larionov S.V., Chirikov S.I., Kyunnap R.I.
The article outlines the theoretical approaches to the optimization of the geometrical parameters of flexible tanks for petroleum products taking into account the total length of the stitches on the top and bottom surfaces. The algorithm, mathematical models and optimization criterion to significantly (up to 20%) increase the efficiency of cutting the material of construction are proposed. On the basis of these curves are the formulas for determining the optimal height, length and width of the flexible container. The calculations are for the optimal sizes for example, flexible container with a capacity of 250 m3.
Keywords: flexible tanks, optimization of geometrical dimensions, Vieta-Cardano method, optimization criterion, construction material, efficiency.
REFERENCES
1. Shukhov V.G. Iskusstvo konstruktsii. Rayner Grefe, Ottmar Perchi, F.V. Shukhov, M.M. Gappoyev i dr. M.: Mir. 1994. 192 p. (In Russian).
2. Patent RF 2304553. Mobil'nyy elastichnyy rezervuar dlya nefteproduktov. Rybakov Y.N., Kharlamova O.D., Patalakh I.I., Abramov M.D. Declared 07.04.2006. Published 20.08.2007. Bulletin № 23. (In Russian).
3. O prinyatii na snabzheniye polevogo sklada goryuchego PSG-600 vmestimost'yu 600 m3: prikaz Ministra oborony Rossiyskoy Federatsii № 132 ot 02.03.2015. (In Russian).
4. TU 7981-023-08151164-2013. Polevoy sklad goryuchego PSG-1000. Tekhnicheskiye usloviya. M.: FAU «25 GosNII khimmotologii Minoborony Rossii». 2013. (In Russian).
5. Issledovaniye po sozdaniyu opytno-promyshlennogo proizvodstva peredvizhnykh elastichnykh rezervuarov: otchet o NIR. M.: FAU «25 GosNII Minoborony Rossii». 2010. (In Russian).
6. Bronshteyn I.N., Semendyayev K.A. Spravochnik po matematike. 7-e izd, stereotipnoye. M.: Gosudarstvennoye izdatel'stvo tekhniko-teoreticheskoy literatury. 1967. Pp. 138-139. (In Russian).
Сведения об авторах
Рыбаков Юрий Николаевич, 1961 г.р., окончил МИНХ и ГП им. И.М. Губкина (1983), кандидат технических наук, заслуженный изобретатель Российской Федерации, лауреат стипендии Президента РФ «Работникам оборонно-промышленного комплекса за выдающиеся заслуги в создании военной и специальной техники», начальник 23 отдела «25 ГосНИИ химмотологии Минобороны России», автор более 200 научных работ, область научных интересов - нефтепродуктообеспечение, полимерные материалы.
Ларионов Сергей Владимирович, 1956 г.р., окончил МИНХ и ГП им. И.М. Губкина (1978), кандидат технических наук, доцент кафедры нефтепродуктообеспечения и газоснабжения РГУ нефти и газа им. И.М. Губкина, автор более 150 научных работ, область научных интересов - технические средства нефтепродуктообеспечения и газообеспечения.
Чириков Сергей Игоревич, 1990 г.р., окончил МАТИ (2013), младший научный сотрудник 23 отдела «25 ГосНИИ химмотологии Минобороны России», автор 5 научных работ, область научных интересов -полимерные материалы.
Кюннап Роман Игоревич, 1989 г.р., окончил УВВТУ (ВИ) (2011), младший научный сотрудник 23 отдела «25 ГосНИИ химмотологии Минобороны России», автор 7 научных работ, область научных интересов - технические средства нефтепродуктообеспечения.
