CC BY
40
Математика
УДК 517.95:532.5 + 532.511 DOI: 10.14529/mmph180401
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ДВИЖЕНИЯ ПОРШНЯ
ПОД ВОЗДЕЙСТВИЕМ ГОРЯЩЕГО ГАЗА С УЧЕТОМ ЗАЗОРА
МЕЖДУ ПОРШНЕМ И ТРУБОЙ
А.В. Геренштейн, Н.С. Мидоночева
Южно-Уральский государственный университет, г. Челябинск, Российская Федерация E-mail: [email protected]
Рассматриваются математические модели пневматической системы, состоящей из трубки, закрытой с одной стороны и открытой с другой. В трубке находится поршень, ограничивающий некоторый объем сжатого газа. Для нахождения параметров движения поршня под действием давления расширяющегося газа строится математическая модель системы несколькими способами: с помощью обыкновенных дифференциальных уравнений и с помощью уравнений в частных производных. В последнюю включаются такие уравнения, как уравнение движения, уравнение неразрывности и уравнение сохранения энергии, т. е. уравнения газовой динамики. Кроме того, определяются соответствующие краевые условия. При этом учитывается возможный нагрев газа и возможные потери некоторого объема газа сквозь имеющийся зазор между цилиндром и поршнем. Все уравнения, входящие в состав математической модели, приводятся к безразмерной форме. Для выполнения расчетов используются методы конечных разностей и характеристик, при которых все частные производные в уравнениях заменяются конечными разностями в узлах некоторой сетки. По имеющемуся шаблону находится приближенное значение каждого уравнения в каждом узле сетки по пространству, затем происходит переход на следующий временной слой. Расчеты выполняются до тех пор, пока поршень не достиг открытого конца трубы или до тех пор, пока поршень не начал замедляться. Затем проводится сравнение результатов, полученных с помощью рассматриваемых методов, по критериям быстродействия и точности, а также даются рекомендации относительно целесообразности использования каждого метода построения математической модели.
Ключевые слова: математическая модель; сжатый газ; пневматическая система.
Введение
Существующие методы расчета параметров пневматических систем не являются полностью приемлемыми. Связано это с тем, что при расчетах математических моделей точными аналитическими методами можно решить лишь некоторые простейшие задачи. Для решения более сложных задач применяются различные численные методы. Численный метод позволяет получить лишь приближенное решение задачи. Методы сравниваются между собой по критериям точности и быстродействия. Следовательно, некоторые методы решений являются более
предпочтительными, а некоторые -менее предпочтительными.
Необходимо исследовать каждый метод и принять решение о его предпочтительности.
В данной работе рассматривается класс систем «поршень-труба» (рис. 1), предназначенных для создания
Газ Поршень Труба
Рис. 1. Система «Труба-поршень» с зазором
ускоренного движения поршня. В общем виде подобного рода системы состоят из двух основных частей: герметичного с одной стороны цилиндра, который содержит сжатый и подогреваемый газ, и подвижного поршня, который под действием давления сжатого газа приводится в движение внутри цилиндра. В этом случае скорость поршня до достижения им конца цилиндра (до достижения максимальной скорости) имеет принципиальное значение (рис. 1) [1-10].
Кроме того, необходимо учитывать возможный нагрев газа в процессе движения поршня, а также возможный зазор между цилиндром и поршнем.
Цель работы - сравнить результаты, полученные с помощью рассматриваемых методов по критериям быстродействия и точности, и дать рекомендации относительно целесообразности использования каждого метода.
1. Обозначения
и = и (х,() - смещение частицы смеси в момент времени ^. р = 1 + Эи/Эх - деформация (удлинение) частицы смеси. V = Эи/Э^ - скорость частицы смеси. Т - температура газа (Кельвин). ср - удельная теплоемкость газа при постоянном давлении. су - удельная теплоемкость газа при постоянном объеме. р - плотность частицы смеси в момент времени ^. р0 - плотность частицы смеси в момент времени ^ = 0, р/3 = р0 . р - давление в частице газа в момент времени ^. Я = ср-су - универсальная газовая постоянная. у = cp|cv - отношение удельных теплоемкостей. q - удельная мощность тепла, выделяемого сгорающим порохом. qs -удельная мощность тепла, выделяемого внешним источником. Р - площадь сечения трубы (и поршня). Р - давление на поршень справа. 10 - начальная длина участка трубы, занятой смесью.
М - масса поршня. V - доля массы твердой фазы (пороха, V = V(7) - заданная функция времени).
2. Математическая модель задачи, построенная с помощью уравнений в частных производных
1.1. С учетом нагрева газа
В основе данной математической модели лежат уравнения газовой динамики [6-10]:
1.1.1. Уравнение неразрывности (совместности)
Имеет место:
(1)
(2)
1.1.3. Уравнение энергии
(3)
Учитывая уравнения совместности, движения и Клапейрона, получим
(4)
или, дифференцируя произведение, деля на рв и приводя подобные, получим
gM+1 = r (qs )(r-1) b dt p dt pb
Геренштейн А.В., Математическая модель движения поршня под воздействием
Мидоночева Н.С. горящего газа с учетом зазора между поршнем и трубой
1.1.4. Начальные и краевые условия
В начальный момент времени (при t = 0) в любом сечении плотность, давление и температура одна и та же, а скорость и деформация равны нулю, поэтому (рис. 2)
р( x,0 ) = р0, p (x,0) = p0, v (x,0) = 0,ß( x,0) = 1. На левом конце (x = 0 ) газ прилегает вплотную к закрытому концу трубы и неподвижен:
v (0, t ) = 0. (5)
Рис. 2. Пневматическая система
На правом конце (х = ¡0 ) поршень движется под действием давления газа и противодавления. Поэтому при х = 10 имеем
Р = Р(1о,X), V = V(1о,X),р = р(¡о,X) и Щ =(р -Р)Р. (6)
1.2. С учетом зазора между цилиндром и поршнем
Рассмотрим пневматическую систему, изображенную на рис. 1. В данной системе поршень неплотно прилегает к стенкам трубки, и возникает промежуток (зазор) между поршнем и трубкой (например, при нарезном канале ствола, поршне нестандартной толщины и др.).
Заметим, что в начальный момент времени газ между поршнем и трубкой не просачивается, эта проблема возникает лишь после начала движения поршня. Имеем, что при движении поршня сквозь зазор в единицу времени просачивается некоторый объем газа. При этом уравнения математической модели не меняются, корректируются лишь методы расчета.
Также необходимо отметить, что расчет параметров движения поршня происходит в одномерном пространстве, и поэтому не учитывается характер течения газа сквозь зазор, а учитываются лишь потери некоторых объемов газа. Так, некоторый объем газа, вытекая через зазор, перестает участвовать в продвижении поршня вперед, и, следовательно, параметры движения поршня меняются.
3. Математическая модель задачи, построенная с помощью обыкновенных дифференциальных уравнений
X - время, и - координата поршня (и(0) = 1), q = q(X) - мощность тепла, выделяемая
единицей массы горящего пороха, V - скорость поршня (V(0) = 0), р - давление, р(0) = 1,су -удельная теплоемкость при постоянном давлении, у - отношение удельных теплоемкостей, т -сумма масс газа и пороха, М - масса поршня, ¡л = т/(уЫ), р1 - противодавление на поршне. а2 = ур/р = урР1р0 - квадрат скорости звука, V - начальная доля массы еще несгоревшего пороха, X - показатель скорости сгорания пороха, Р - площадь сечения трубы.
(7)
M^V = F(Р -Pi).
При равномерном расширении газ рассматривается как единое целое, поэтому применим условие (7) и выразим р через квадрат скорости звука. При этом система уравнений (1)-(3), (6) сводится к уравнению:
v
2 , ( 2 a v 1 a
= r(r-1)q-(g-1)—-— i1 - me1)-Pi r(r-1)
с краевыми условиями:
V (0, t ) = 0,
Mvt = ( p - P) F.
(8) (9)
4. Методы вычислений
1.1. Приведение к безразмерному виду
Обозначим а0 = ^ур0/р0 (исходная скорость звука).
Положим
Р
—. х =—, V =-, Р =
I
. = ta0 x = x v = p = p m = PsFk p= P q = qlor(r-1) = qsior(r-1)
i
an
3
0 1 а0 Р0 МГ Р0 а0
Теперь под г,х,V,р,д,дв будем подразумевать г',х',у', р',д',дя . В безразмерных переменных получим уравнение совместности:
= Эу Э' Эх
^ qs =-
a
Уравнение движения Уравнение энергии
dv 1 dp _ — + —— = 0. dt r dx
rdb+1 dp=qv+qs
/ dt p dt p/
(10) (11) (12) (13)
Начальные условия (при ' = 0)
p (х,0 ) = 1, v (x,0) = 0,Ь( x,0 ) = 1. Условие на дне трубы (x = 0 )
v (0, ' ) = 0.
Условие на поршне примет вид (при х = 1):
ЭУ="( Р -п). ....
После решения задачи (8)-(13) надо параметры ',х,у,р заменить на г',х',у',р', а затем эти «штрихованные» параметры с помощью соотношений (12) снова заменить на исходные параметры ', х, у, а . Эти манипуляции выполняются исключительно для упрощения записей. Замечание. Обозначим
'д (*>(*)+дs (*)
^ (t )=р
0 p(t)/(t) В уравнении (3) перейдем к квадратурам. Получим
-dt.
p = p0
f / Y / J
g (t)
Произведя некоторую перегруппировку сомножителей, получим такое представление [1-5]
'■д (*>(*)+дs (*)
г
p = J" 0
/(t)
-dt.
1.2. Метод конечных разностей
2
Геренштейн А.В., Мидоночева Н.С.
Математическая модель движения поршня под воздействием горящего газа с учетом зазора между поршнем и трубой
В соответствии с шаблоном (рис. 3) получим следующее соотношение для уравнения совместности:
4+1
ь =ь+
>/-1
2И
ст.
(ч+И
Для уравнения движения:
-1 - р,±1 т г г 2Н у Уравнение энергии:
'я(т)у(т)±Ч* (т)
(ил Ь (ч) Ь
Рис. 3. Шаблон
(H-l.il
Р=V
Ь(Т)
ёт.
1 2 3
И
1.3. Учет зазора между поршнем и трубкой
Для выполнения расчетов с учетом зазора между поршнем и трубкой так же строится сетка по методу конечных разностей, то есть полость трубки разбивается на участки, как показано на рис. 4
Затем на первом шаге работы метода считаем, что прилегающий к поршню слой под номером п полностью вытек сквозь зазор между поршнем и трубкой. Тогда необходимо скорректировать положение поршня, и, поскольку реальное продвижение поршня в направлении
незакрытого конца трубки меньше рассчитанного, то мы «сдвигаем» поршень в противоположном направлении. То есть новое, скорректированное положение поршня более соответствует действительности.
Затем наступает очередь слоя п-1, затем п—2 и т. д.
Слои газа Поршень
Рис. 4. Схема расчета параметров движения поршня
о
5. Сравнение математических моделей
Сравнение математических моделей, в основе которых лежат обыкновенные ДУ и уравнения газовой динамики с учетом нагрева газа (подогрев 2-10+7 Дж/кг): решение, полученное с помощью обыкновенных дифференциальных уравнений: время 0,00240744 с, деформация 11,5002, путь 0,500009 м, скорость 235,899 м/с, давление 1,65846 атм. Среднее время работы метода 0,06 с.
С помощью уравнений газовой динамики: время 0.00240771 с, деформация 11,5789, путь 0,500003 м, скорость 235,87 м/сек, давление 1,64565 атм. Среднее время работы метода 7 с.
Разность скоростей: 0,02879 м/с, что в процентах составляет: 0,012206 %. Разность путей: 0,006263, что в процентах: 0,001253 % (рис. 5).
Полученные данные показывают, что метод обыкновенных дифференциальных уравнений позволяет получить решение задачи в несколько раз быстрее (в 117 раз для данного примера). При этом метод конечных разностей более предпочтителен, чем метод характеристик.
6. Практическая значимость
Применяемые способы
математического моделирования
применяются к расчету параметров пневматических систем с поршнем.
Выводы
- В ходе исследования были рассмотрены методы построения моделей «труба-поршень» и выявлены их недостатки;
—
/А
/
// 'Л
/< /1 Шаг по в ремеми 1
-Точное аналитическое решение
~ Обыкновенные дифференциальные уравнения Уравнения газовой динамики
Рис. 5. График зависимости скорости поршня от времени
- построены математические модели движения поршня в трубе с учетом давления газа при равномерном и неравномерном расширении газа, с учетом нагрева газа и его проникновением в зазор между поршнем и трубой;
- сравнительный анализ математических моделей позволил сделать следующие выводы: решение, полученное в случае уравнений газовой динамики, является более точным, нежели в случае обыкновенных дифференциальных уравнений: отклонения от точных аналитических решений в простых задачах составляют не более 0,4-1,1 % (тогда как в случае обыкновенных дифференциальных уравнений до 1,3 %);
- но при рассмотрении нагрева газа метод обыкновенных дифференциальных уравнений позволяет получить решение почти в 117 раз быстрее, нежели метод уравнений газовой динамики при незначительной разнице в результатах. При этом метод конечных разностей более предпочтителен, чем метод характеристик.
Литература
1. Клейман, Я.З. Скорость звука в смесях, содержащих взвешенные частицы / Я.З. Клейман // Акустический журнал. - 1961. - Т. 7, № 2. - С. 262-264.
2. Клейман, Я.З. О распространении сильных разрывов в многокомпонентной среде / Я.З. Клейман // Прикл. математика и механика. - 1958. - Т. 22, № 2. - С. 197-205.
3. Клейман Я.З. Некоторые особенности движения смесей / Я.З. Клейман // Акустический журнал. - 1959. - Т. 5, № 2. - С. 157-165.
4. Рахматулин, Х.А. Основы газодинамики взаимопроникающих движений сжимаемых сред / Х.А. Рахматулин // Прикладная математика и механика. - 1956. - Т. 20, № 2. - С. 184-195.
5. Голубятников, А.Н. Интегральные неравенства в задачах газовой динамики / А.Н. Голубятников // Аэромеханика и газовая динамика. - 2001. - № 1. - С. 74-81.
6. Голубятников, А.Н. К оптимизации ускорения тела в классе движений толкающего газа с однородной деформацией / А.Н. Голубятников, Н.Е. Леонтьев // Аэромеханика и газовая динамика. - 2001.- № 2.- С. 27-34.
7. Кушнер, Е.Н. Нормальные формы некоторых уравнений газовой динамики / Е.Н. Куш-нер // Научный вестник Московского государственного технического университета гражданской авиации. - 2013. - № 194. - С. 20-23.
8. Баутин, С.П. Одно точное стационарное решение системы уравнений газовой динамики / С.П. Баутин, А.Г. Обухов // Известия высших учебных заведений. Нефть и газ. - 2013. - № 4. -С.81-86.
9. Рылов, А.И. Функциональная зависимость между законами сохранения газовой динамики, отвечающими разделению переменных / А.И. Рылов // Доклады Академии наук. - 2014. - Т. 454, № 6. - С. 647-650.
10. Газодинамические основы внутренней баллистики / С.А. Бетехтин, А.М. Виницкий, М.С. Горохов и др. - М.: Оборонгиз, 1957. - 384 с.
11. Жаровцев, В.В. Математическое моделирование и оптимальное проектирование баллистических установок / В.В. Жаровцев, Л.В. Комаровский, Е.И. Погорелов. - Томск: Издательство Томского университета, 1989. - 253 с.
12. Корнер, Дж. Внутренняя баллистика орудий / Дж. Корнер. - М.: Издательство иностранной литературы, 1953. - 462 с.
13. Геренштейн, А.В. Математическая модель движения поршня в трубе при действии давления газа / А.В. Геренштейн, Н.С. Кастрюлина // XVII Международная научно-практическая конференция: «Научное обозрение физико-математических и технических наук в XXI веке». - 2015. - №5(17). - С. 134-138.
Поступила в редакцию 11 октября 2017 г.
Геренштейн А.В., Мидоночева Н.С.
Математическая модель движения поршня под воздействием горящего газа с учетом зазора между поршнем и трубой
Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2018, vol. 10, no. 4, pp. 5-12
DOI: 10.14529/mmph180401
MATHEMATICAL SIMULATION OF PNEUMATIC SYSTEM WITH A CLEARANCE BETWEEN THE PISTON AND THE PIPE
A.V. Herreinstein, N.S. Midonocheva
South Ural State University, Chelyabinsk, Russian Federation E-mail: [email protected]
The article regards mathematical models of pneumatic system, consisting of a pipe which is closed from one end and open from the other. In the pipe, there is a piston that limits some volume of compressed gas. In order to determine parameters of the piston's motion under the pressure of expanding gas, mathematical model of the system gets constructed using several methods: with the use of ordinary differential equations, and with the use of partial differential equations. The last method includes such equations as motion equation, continuity equation, and energy conservation equation, i.e. the equations of gas dynamics. Besides, corresponding boundary conditions get determined. At that, possible heating of the gas and probable loss of some volume of the gas through an existing clearance between the cylinder and the piston are taken into account. All equations included into the mathematical model get reduced to the dimensionless form. Methods of finite differences and characteristics are used for calculations, at which all partial derivatives in equations get replaced with finite differences in nodes of a grid. By the existing template, approximate value of each equation gets determined in each node of the grid by the space; then a transition to the next temporal layer takes place. Calculations are being performed either until the piston reaches the open end of the tube or until the piston started to slow down. After that, results obtained with the use of methods under consideration get compared by the criteria of fast operation and accuracy, and recommendations regarding advisability of using each method for construction of a mathematical model get provided.
Keywords: mathematical model; compressed gas; pneumatic system.
References
1. Kleyman Ya.Z. Akusticheskiy zhurnal, 1961, Vol. 7, no. 2, pp. 262-264. (in Russ.).
2. Kleyman Ya.Z. Prikladnaya matematika i mekhanika, 1958, Vol. 22, no. 2, pp. 197-205. (in Russ).
3. Kleyman Ya.Z. Akusticheskiy zhurnal, 1959, Vol. 5, № 2, pp. 157-165. (in Russ.).
4. Rakhmatulin Kh.A. Prikladnaya matematika i mekhanika, 1956, Vol. 20, no. 2, pp. 184-195. (in Russ).
5. Golubyatnikov A.N. Aeromekhanika i gazovaya dinamika, 2001, no.1, pp. 74-81. (in Russ.).
6. Golubyatnikov A.N., Leont'ev N.E. Aeromekhanika i gazovaya dinamika, 2001, no. 2, pp. 27-34. (in Russ.).
7. Kushner E.N. Normal FORMS of Some Gas Dynamics Equations. Civil Aviation High Technologies, 2013, no. 194, pp. 20-23. (in Russ.).
8. Bautin S.P., Obukhov A.G. A single exact stationary solution of the gas dynamics equations system. Proceedings of Higher Educational Establishments. Oil and Gas, 2013, no. 4, pp. 81-86. (in Russ).
9. Rylov A.I. Doklady Akademii nauk, 2014, Vol. 454, no. 6, pp. 647-650. DOI: 10.7868/S0869565214060073
10. Betekhtin S.A., Vinitskiy A.M., Gorokhov M.S. et al. Gazodinamicheskie osnovy vnutrenney ballistiki (Gas dynamic fundamentals of internal ballistics). Moscow, Oborongiz Publ., 1957, 384 p. (in Russ).
11. Zharovtsev V.V., Komarovskiy L.V., Pogorelov E.I. Matematicheskoe modelirovanie i opti-mal'noe proektirovanie ballisticheskikh ustanovok (Mathematical modeling and optimal design of ballistic installations). Tomsk, Izdatel'stvo Tomskogo universiteta Publ., 1989, 253 p. (in Russ.).
12. Korner, Dzh. Vnutrennyaya ballistika orudiy (Theory of Interior Ballistics of Guns). Moscow, Izdatel'stvo inostrannoy literatury Publ., 1953, 462 p. (in Russ.). [Corner, J., Theory of Interior Ballistics of Guns, John Wiley and Sons, Inc., New York, 1950, 443 p.]
13. Gerenshteyn A.V., Kastryulina N.S. Matematicheskaya model' dvizheniya porshnya v trube pri deystvii davleniya gaza (Mathematical model of piston motion in a pipe under the action of gas pressure). XVII Mezhdunarodnaya nauchno-prakticheskaya konferentsiya: "Nauchnoe obozrenie fiziko-matematicheskikh i tekhnicheskikh nauk v XXI veke " (Proc. XVII International Scientific and Practical Conference: "Scientific Review of Physical, Mathematical and Technical Sciences in the 21st Century"), 2015, no. 5(17), pp. 134-138. (in Russ.).
Received October 11, 2017
