Машиностроение
УДК 539.3; 532.5; 517.9
О НЕКОТОРЫХ МАТЕМАТИЧЕСКИХ МОДЕЛЯХ МЕХАНИЧЕСКОЙ СИСТЕМЫ «ТРУБОПРОВОД - ДАТЧИК ДАВЛЕНИЯ»1
П.А. Вельмисов, Ю.В. Покладова
Ульяновский государственный технический университет 432027, Россия, г. Ульяновск, ул. Северный Венец, 32 E-mail: pokladovau@inbox.ru
Рассматриваются математические модели механической системы, включающей в себя трубопровод с рабочей средой и датчик, составной частью которого является упругий элемент. Датчик предназначен для определения давления рабочей среды (например, на выходе из камеры сгорания двигателя), закон изменения которого считается заданным. Получены дифференциальные уравнения, описывающие колебания упругого элемента, и на их основе проведен численный эксперимент по исследованию динамики элемента.
Ключевые слова: трубопровод, датчик давления, деформация, упругий элемент, интег-ро-дифференциальные уравнения, динамика.
Введение
Независимо от принципа преобразования все датчики давления в той или иной степени критичны к воздействию широкого диапазона температур и повышенных уровней виброускорений. Размещение датчика давления непосредственно на корпусе двигателя принципиально обеспечивает более высокую достоверность измерения, но, как правило, сопровождается воздействием на датчики давления высоких температур и виброускорений. Причем во многих реальных случаях (например, при взлете и посадке аппарата) температура воздействия как на авиационный, так и на ракетный двигатель носит нестационарный характер. В результате датчики подвергаются воздействию нестационарной температуры, повышенных виброускорений, что приводит к погрешности измерений, а в ряде случаев к разрушению упругого чувствительного элемента датчика.
В связи с вышеизложенным, возникает задача проектирования механической системы «трубопровод - датчик давления». В системе датчик расположен на некотором расстоянии от двигателя и соединен с ним с помощью трубопровода, что позволяет ослабить воздействие высоких температур и виброускорений. Задача состоит в получении уравнений, связывающих закон изменения рабочей среды на входе в трубопровод (на выходе из камеры сгорания двигателя) и деформацию упругого элемента датчика, и предназначенных по величине деформации элемента рассчитать давление в двигателе. Математические модели системы «трубопровод - датчик дав-
1 Работа выполнена в рамках реализации ФЦП «Научные и научно-педагогические кадры инновационной России» на 2009-2013 годы (НК-177П, ГК №П1183).
Петр Александрович Вельмисов - д.ф.-м.н., профессор.
Юлия Валерьевна Покладова - к.ф.-м.н.
ления» рассматривались в [1-5]. Здесь предложены новые модели системы «трубопровод - датчик давления», отличающиеся или уравнениями, описывающими колебания упругого элемента датчика (нелинейные модели), или геометрией трубопровода.
1. Математические модели механической системы «трубопровод - датчик давления»
Рассматриваются плоские модели механической системы «трубопровод - датчик давления» для трубопровода конечной длины (см. рис.1,а), бесконечно длинного трубопровода (см. рис.1,б) с датчиком, закрепленным на боковой стенке, и бесконечно длинного трубопровода с датчиком, расположенным на стенке полости трубопровода (см. рис.1,в).
Рис. 1. Схема механической системы:
1 - двигатель, 2 - трубопровод, 3 - датчик, 4 - рабочая среда, 5 - пластина (упругий элемент датчика).
Рассмотрим задачу о динамике упругого элемента датчика давления рабочей среды, расположенного на боковой стенке трубопровода конечной длины (рис.1,а).
В линейной постановке уравнения и граничные условия имеют следующий вид:
Фхх + Ф уу = 0 (Х У) е в = {(X, у): 0 < х < Xo, О < у < Уо}; (1)
Фу (х,0, t) = 0, х е (0, х0); (2)
Фх(0, y, t) = 0 y е (о, yo); (3)
Фу (х, jy0, t) = 0, х е (0,a) u (b,х0); (4)
Фу(х,y0,t) = со(х,t), х е (a,b); (5)
P-pФt (xo,У, t) = р*(У t) У е (0,J0); (6)
L(o) = Meó + Dra"" + Nb' + acó "" + pcó + ую =
= P - р0(хt)-PФt Уo, 0. х е (ab)- (7)
Здесь (1) - уравнение Лапласа, описывающее движение рабочей среды (жидкости) в трубопроводе; (2)-(5) - условия непротекания среды через соответствующие границы; условие (6) задает закон изменения давления на входе в трубопровод; (7) -уравнение динамики пластины; ф( х, у, t)- потенциал скорости среды; ю( х, t) - прогиб упругого элемента; Р*(у, t) - закон распределения давления среды в сечении х = х0 (на выходе из двигателя); Р0(х, t) - распределенная внешняя нагрузка, действующая на упругий элемент; хо, Уо- продольный и поперечный размеры трубопровода; а, Ь - координаты концов упругого элемента; Р - давление рабочей среды в трубопроводе в состоянии покоя; р - плотность среды; М - погонная масса; Б -изгибная жесткость; N - сжимающее (растягивающее) усилие; а, р - коэффициенты внутреннего и внешнего демпфирования; у - коэффициент жесткости основания; точка и штрих, также как и индексы t и у снизу, обозначают частную производную по t и у соответственно.
Если продольный размер трубопровода хо значительно превосходит его поперечный размер уо, то можно считать хо =+<х> (рис.1,б). Для модели, изображенной на рис. 1,б, линеаризованные уравнения и граничные условия имеют следующий вид:
= 0, (х,у) е G = {(х,y): 0 < х <+<х>,0 < y < У0}; (8)
ФУ (х,0, t) = 0, х е (0,+<х>); (9)
ФУ (х, y0, t) = 0, х е (0, a) u (b,+<x>); (10)
Фх(0, y t) = ^ y е (0, y0); (11)
Ф y (х, y0, t) = ó (х, t), х е (a, b); (12)
lim (ф2 +ф2у)= 0 ; х (13)
lim (р -pФt )= p*(t); (14)
х ^+^
L(o) = Meó + Da'"' + Nb' + acó '"' + pd) + ую =
= P - P0(х, t) -рфt (х, У0, t), х е (a, b). (15)
Здесь P.(t) - закон изменения давления на входе в трубопровод; Р0(х,t) - распределенная внешняя нагрузка.
В линейной постановке, соответствующей малым прогибам упругого элемента и малым возмущениям потенциала скорости рабочей среды, математическая модель, геометрическая схема которой изображена на рис.1,в, определяется следующими уравнениями и граничными условиями:
Фхх + Фуу = ^ (х у) е G0; (16)
Фу (х ,0, t) = 0, х е (0, да); (17)
Фу (х, У0, t) = 0, х е (1, да); (18)
Фх(0, уt) = ^ у е (0, H); (19)
Фх(1, у t) = ^ у е (Уo, H); (20)
ФУ (х, H, t) = 0 , х е (0, a) u (b, 1); (21)
ФУ (х, H, t) = a (х, t), х е (a, b); (22)
lim (p* -pфt(x, y t)) = p(t), y е (0, У0 ); (23)
х ^-+<»
L(a) = Meó + Da'"' + Na' + acó "" + pcó + ya =
= P* -рФt (х, H, t) - P0(х, t), х е (a, b). (24)
Постановка задач (1) - (7), (8) - (15), (16) - (24) соответствует линейной теории аэрогидроупругости, когда динамика жидкости (газа), а также динамика чувствительного элемента датчика описываются линейными уравнениями.
Линейные математические модели системы «трубопровод - датчик давления» рассматривались в [1 - 5]. В данной статье авторами предлагаются нелинейные модели, которые определяются приведенными выше уравнениями (1) - (6), (8) - (14), (16) - (23), при этом в уравнениях (7), (15), (24), описывающих динамику пластины, дифференциальный оператор Х(ю) заменяется следующим:
Da"
L(a) = Meó
' 1 + (о')2 Полагая прогибы малыми, заменим
+ Na + ace — 5со + pee + seo + ya + цш .
1 (1 - -3 (a')2 ). В результате
Л2 ~
1 + (a')
получим
3
L(a) = Mió + Da” -2Do''''(o')2 -9DaW' -3D(o'')3 + Na'' + acá'''' -
- 5co ' + peo + ya + цш3 + sa 3.
(25)
2. Сведение задачи к исследованию уравнения для функции деформации упругого элемента
В случае трубопровода конечной длины (рис.1,а) решение аэрогидродинамиче-ской задачи осуществлялось на основе метода Фурье, при этом потенциал скорости рабочей среды был представлен в виде ряда.
В результате получено уравнение, связывающее закон изменения давления Р.(у, t) рабочей среды на входе в трубопровод (х = х0) и функцию прогиба ю(х,t) упругого элемента датчика, расположенного в сечении х = 0 :
Цга) = -Р0(х,¿) + — ^Р*(у,t)dy-—I^(^пх)С^ п^о)^со(х,I)cos(Xnx)dx+
У0 0 х° п=1 '“п
2 Х"1 I г \ ' +—^ск(у пх)- )
У0 П=1 Ск(хпх0)
,С0^ пУ0)
У0
| Р*( у Г )cos(V пУ^У,
(26)
пп
(2п - 1)п
у0 2х0
Это уравнение выражает связь между прогибом ю(х, t) упругого элемента датчика и законом изменения давления Р*( у, t) на входе в трубопровод.
В случае бесконечно длинного трубопровода (рис.1,б, рис.1,в) на основе методов теории функций комплексного переменного (с помощью интеграла Кристофеля-Шварца, формул Шварца и Сохоцкого) получено уравнение, связывающее закон изменения давления Р^)на входе в трубопровод и деформацию ю(х,t) упругого элемента датчика.
Для модели, изображенной на рис.1,б, имеем
ь
Дю) = Р.^) -Р0(х,0 + — Геб(х,t)ln
П ■>
. и . пх
СП------------------СП —
У0
У0
ds,
(27)
а для модели, изображенной на рис.1,в -
Дю) = Р^) - Р0 (х, t) — — Геб (т, t) 1п п ■>
^(т) - п
^(т) Ч( х)
где
^( х) - функция, обратная к функции х(^) = С0
dт , х е (а,Ь), (28)
, £, е [0,1]. Концы
х(1 - х) п - х
интервала (а,Р) определяются из условий: х(а) = а , х(Р) = Ь. Параметры С0, т , п выражаются через параметры I, У0, Н .
а
0
а
Ь
а
0
3. Решение уравнения для деформации
Согласно методу Галеркина, решение ю( х, t) уравнений (26), (27), (28) ищется в
т
виде ю(х,t) = 1®к^)gk(х), где ^к(х)}^ - полная на [а,Ь] система базисных
к=1
функций, удовлетворяющих граничным условиям, соответствующим условиям закрепления пластины.
Например, в случае жесткого защемления концов пластины (ю = 0, ю' = 0, при х = а, х = Ь) решение уравнений (26), (27), (28) будем искать в
т
виде ю( х, t) = юк (0Ук(x), где
к=1
Vк(х) = сМ^к(х - а)) - с°^к (х - а)) -
- ^к (Ь- а))- ^кГ В ^к (X - а)) - ™(цк (X - а))),
к(Ь - а)) - яп(Цк(Ь - а))
при этом ^к находятся из уравнения Л(|а,к (Ь - а))соэ(цк (Ь - а)) = 1 (к = 1,..,т ).
ь
Функции {ук (х)}к=1 ортогональны на [а,Ь], т.е. |у г (х)у^ (х)йх = 0 при г Ф ] . Мож-
а
Ь
но показать
, что | у 2 (х)йх = Ь - а .
Для шарнирного закрепления концов пластины (ю = 0,ю" = 0 при х = а, х = Ь)
%к
можно положить ю(х,г) = ^ак (г)sinрк(х - а), где Рк =
к=1
Ь - а
Из условия ортогональности невязки уравнений (26), (27), (28) к системе базисных функций ^ (х)}т=1 получим систему из т обыкновенных дифференциальных уравнений для Юк(I)(к = 1,..,т ).
Например, для модели, изображенной на рис.1,в, получим систему из т обыкновенных дифференциальных уравнений для Юк (I) (шарнирное закрепление концов пластины, а = 0, Ь = I, Xк =пк/
где
^ Аік а к(г) + Віа і (£) + С{юг- (г) + ^(£) = 0,
к=1
Ак =-
1п
і (і М (|
0 V 0
Аіі = М- + ^ И 2
^(г)- п
£(г) Ч(х)
і (і
Л
8Ш Xк гйг
8т Х{хйх, і Ф к.
£п
к - і
1п
^(г) - п
^(г) Ч(х)
Л
sin Х; г<іг
sin Хіхйх, і = к ,
о V 0 - • у
Ві = ^(аХ4 +р), Сі = -М,2 +у),
2
(г) = | Р0 (х, г) 8ш Х^хйх + Р(г)
1 -(-1)
о
і, к = 1,.., т.
Хі
а
4. Численное моделирование
Для всех рассмотренных моделей задача Коши для систем из т обыкновенных дифференциальных уравнений решается с помощью системы МаШешайса 7.0. Проведено численное моделирование на ЭВМ динамики упругого элемента датчика в зависимости от закона изменения давления в двигателе. Исследовалась деформация элемента как функция времени (в фиксированных точках элемента) и как функция координаты (в фиксированные моменты времени) для различных параметров механической системы.
Пример. Рассмотрим модель, изображенную на рис.1,а (жесткое закрепление концов упругого элемента, нелинейная модель). Рабочая среда - вода (р = 103 кг/ 3),
/ м
пластина изготовлена из алюминия (Е = 7 • 1010 Па, у = 0,36, ро = 2,7 • 103 кг/мг , h = 1,5 40-3 м).
Для значений параметров х0 = 0,3м , у0 = 0,015м, а = 0 м, Ь = 0,015м, М = 4,05 к/ 2, Б = 22,619 Дж, N = 103 Н , а = 0,5 мк/ , Р = 0,3 к/ 2 , у = 0,2 к/ 22,
/м /с * / м с / м с
N1 = 15.
5 = 0'
8 = 0 кг^с
о(х,0) = 0 м, а> (х,0) = 0-7
Р0 (х, г) = 0 Па, Р*(у, г) = 105 (10 + cos(10г))Па получено решение для функции
о(х, г) = ^Юк (г)у к (х) в точке х* =
а + Ь
(рис.2):
к=1
¿у (*-,/) 0.0001
о.ооооз о.оооо в
0.00004
0.00002
кЛМЛЛМЛММЛЛЛ/
1
10
2 4 6 3
Рис. 2
Прогибы упругого элемента в фиксированные моменты времени г = г0 для нелинейной модели представлены на рис.3.
о.оос
0.0002
0.0001
0.0000
*□ = 0,005
<и(Х,Ід) ^ = 0,01
¿¿{Х.Ід}
Рис. 3
На рис.2, 3 представлены графики, полученные при численном моделировании на ЭВМ в случае нелинейной модели для трубопровода конечной длины с упругим элементом на боковой стенке (рис.1,а). Линейная модель рассматривалась в [1, 5]. Графики для линейной и нелинейной моделей для одних и тех же значений параметров механической системы существенно различаются (в частности, в нелинейной модели происходит увеличение деформации). Можно сделать вывод, что учет нели-
4
м4с2
2
нейных членов в уравнении, описывающем динамику пластины, имеет важное значение при исследовании динамики чувствительного элемента.
Заключение
Развитие авиационной техники требует постоянного совершенствования и разработки новых типов первичных преобразователей, в частности, датчиков давления. Поэтому естественным образом возникает актуальная задача разработки специальных методов исследования динамики и устойчивости упругих элементов датчиков давления, взаимодействующих с жидкостью. Предложенные авторами новые модели, методика решения задач аэрогидроупругости и их численная реализация существенно позволяют дополнить базу современного проектирования датчиков давления и усовершенствовать ее.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Vel’misov P. A., Pokladova Yu. V. An investigation of mathematical models of a mechanical system “Pipeline-Pressure Sensor” // Romai Journal. Pite§ti, Romania. 2005. v.2, №1. - Р.214-219.
2. Vel’misov P. A., Pokladova Yu. V. Investigation of dynamics of an elastic element of a pressure sensor // Applications of Mathematics in Engineering and Economics. - Soft trade, Sofia, Bulgaria. 2006. - P.51-57.
3. Анкилов А. В., Вельмисов П.А., Горбоконенко В.Д., Покладова Ю.В. Математическое моделирование механической системы «трубопровод-датчик давления» / Ульяновск: УлГТУ, 2008. 188с.
4. Анкилов А. В., Вельмисов П.А., Покладова Ю.В. Математические модели механической системы «трубопровод-датчик давления» // Вестник Саратовского государственного технического университета. 2007. №3. С.7-14.
5. Вельмисов П.А., Горбоконенко В.Д., Решетников Ю.А. Математическое моделирование механической системы «трубопровод - датчик давления» // Датчики и системы. 2003. №6(49). С.12-15.
6. Вельмисов П.А., Покладова Ю.В. О некоторых математических моделях механической системы «трубопровод - датчик давления» // Proceeding of the international conference “Education, sience and economics at universities. Integration to internatibional education area” (Poland, Plock, 20.09.10 -25.09.10). - Plock: NOVUM, 2010. - P.492-499.
Статья поступила в редакцию 24 февраля 2011 г.
UDC 539.3; 532.5; 517.9
ON SOME MATHEMATICAL MODELS OF MECHANICAL SYSTEM “PIPELINE - PRESSURE SENSOR” P.A. Velmisov, Yu. V. Pokladova
Ulyanovsk State Technical University
32, Severny Venets st., Ulyanovsk, 432027
In this work the mathematical models of the mechanical system including the pipeline with a working medium and the sensor with the elastic element are considered. The sensor is intended for determination ofpressure of a working medium (for example, on an output from the camera of combustion of the engine) law of change of which is considered set. The differential equations describing oscillations of an elastic element are received, and on their basis numerical experiment on research of dynamics of an element is made.
Keywords: pipeline, pressure sensor, deformation, elastic element, integro-differential equations, dynamics.
Petr A. Velmisov - Doctor of Physics and Mathematics, Professor. Yulia V. Pokladova - Candidate of Physics and Mathematics.