Файл: Козобков, А. А. Электрическое моделирование вибраций трубопроводов.pdf
ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 24.10.2024
Просмотров: 98
Скачиваний: 0
20 |
T |
1 |
|
|
(110) |
E l d*X
Ho dx‘l
полагая, что E J ( x ) = E / =const, |а0(д:)= ч0—const, т. e., что тру бопровод однороден.
Уравнение (ПО) распадается на два самостоятельных урав нения
(111)
( 112)
Для определения собственных частот трубопровода, описы ваемого уравнением (23), из уравнений (111) и (112) нас будет интересовать лишь уравнение (112) для одной s-той собственной формы.
Очевидно, что если уравнению (112) будет удовлетворять ре шение Xs, то ему будет удовлетворять и решение уравнения
uP uq |
будем иметь |
(108). Тогда, обозначая — |
EJ
:(п з)
Учитывая уравнения (65), (66) и (113), обозначим началь ные значения производных или значения производных в сечении
х=0:
Учитывая начальные значения производных, найдем реше ние уравнения (113) в виде системы линейных гиперболо-три гонометрических уравнений
—yfiskVx + 0оs$x+ M 0s~ + Qos - -xk- ; |
(ИЗ) |
||
M s = y 0tE J ^ l / x + |
B0sE J k V x + |
M 0sS x + Q 0t Ц - ; |
(116) |
< 1 ,= У о ' Е Л * Т х + |
% , Е Л ? и я + |
М 0£ У л + 4 ъ р „ |
(117) |
48
где Тх, Vx, Ux — известные функции Крылова:
с |
cos кх + |
ch кх |
х ~ |
2 |
( 1 1 8 ) |
’ |
Л‘
„ |
, |
Г n |
|
, |
sin кх + sh кх |
|
|
|
T x= k |
j |
Sxdx = |
-----------------; |
|
|
|||
|
|
О |
|
|
|
|
|
|
и х= & j |
[ |
' ^ ^ |
= - соз^ 2+с11*- |
; |
(120) |
|||
|
|
6 |
6 |
|
|
|
|
|
V д. = k3 | |
j1^ Sxdx dx d x = ~ 81n ** + sh l!X ; |
(121) |
||||||
|
|
обо |
|
|
|
|
||
Sx = |
Id | |
j‘ |
j‘ |
j‘ s xdx dx dx dx = cos |
kx . |
(122) |
||
|
|
0 |
0 |
0 |
0 |
|
|
|
Полученная система уравнений (114)—(117) и соотношения (118)—(122) позволяют определить значение s-той собственной частоты однородной трубы, которую необходимо знать, чтобы определить частотную погрешность моделирования.
Для нахождения спектра собственных частот трубы со5 необ ходимо знать граничные условия. Пусть закрепления будут шар нирными, т. е.
МД0)=Л45(/)=0; z/s(0 )= ys(Z)=0,
где Msl и ys(l)—момент и перемещение для s-той собственной формы на конце трубы длиной I. Подставляя выражения для граничных условий в соотношение (114) и (116), после некото рых преобразований получим выражение для спектра собствен ных частот
Раскрывая это выражение для Тх и Vx по формулам (119) и (121), получим выражение для спектра собственных частот участка трубопровода длиной I с шарнирным закреплением
СО |
(123) |
5 |
|
Для нахождения частотной погрешности |
необходимо |
еще найти значения s собственных частот a Ms системы, описы
ваемой уравнением (52). Для этого потребуется определение s первых собственных частот системы, описываемой п—1 диффе ренциальными уравнениями вида
49
Ho
EJ
HO
EJ
HO'
EJ
1s*toho 1 dt2
dfl
dt2
54£2
Дд-1 54£3
Дл-4 <И£4
Дл-1
|
|
|
( 1 2 4 ) |
НО |
<1-Ут-2 _ |
В;1£ш 2 |
|
EJ |
dfl |
Дл-1 |
|
НО |
Ц2г/„_2 _ |
В-1£"—2 |
|
EJ |
dfl |
Д-v'i |
|
— с шарнирным закреплением |
|
(125) |
|
Уо = Уп= 0\ M0 = Mn = Q. |
|||
Получение аналитического решения системы уравнении (124) |
|||
с граничными условиями |
(125) в |
конечном |
виде не представ |
ляется возможным, поэтому воспользуемся несколько иным ме тодом нахождения собственных частот системы, описываемой уравнением вида (52).
Предположим, что s-тая собственная форма системы, описы ваемой уравнением (52), совпадает с s-той известной формой системы, описываемой уравнением (23), что возможно при до статочно большом числе разбиений п. Тогда по известной собст венной форме найдем частоту <mms для этой формы. Найдем вы ражение для s-той собственной формы. Для этого подставим выражения для шарнирных граничных условий в систему урав нений (II4) —(117). После преобразований получим
|
|
|
у ( х ) = — 5o£_sinyfex. |
|
|
(126) |
||||
|
|
|
|
1 |
EJk? |
|
|
|
|
|
|
Следовательно, |
учитывая |
сделанное выше |
предположение, |
||||||
можно записать с учетом выражения (126) |
|
|
|
|||||||
|
|
Ут+Р=У{хт + Lxp)= - |
sin (xm+ |
Дхр) k, |
|
(127) |
||||
где р — любое |
целое |
число, удовлетворяющее |
условию |
|||||||
т + |р |^ / г . Перепишем |
выражение |
(52-) с учетом уравнения |
||||||||
(127). Тогда получим |
|
|
|
|
|
|
||||
H-o |
diy, |
kl |
—y« |
E J |
sin k - |
(in-j-2) —4 sin — (m-j- 1)Д- |
||||
dft |
Sin--- in — |
Дх-l |
||||||||
|
n |
|
n |
|
n |
EJ |
|
|||
-(-6 sin — m — 4sin — (m — 1)-)- sin — (m — 2) = — Уm |
P. |
|||||||||
Да"1 |
||||||||||
|
|
n |
n |
|
n |
|
|
(128) |
||
|
|
|
|
|
|
|
|
|
50
Как известно из теории дифференциальных уравнений, зна чения собственных частот системы, описываемой уравнениями вида (128), удовлетворяют следующему выражению:
------- —------- 3. (129) ki
(л0Д x'l sin — т
п
Подставляя уравнения (129) и (123) в выражение (107) для частотной погрешности, после неко торых преобразований получим
и- . 2 |
|
|
0,6 |
Sin- |
|
. l i |
|
= 1 |
18 |
0,5 |
|
±_ |
324’ |
|
2
(130) 0,4
где q — величина, характеризую щая число звеньев модели (или чис ло разбиений) на одну полуволну
ч=- П-Л
Аналогичным |
образом |
могут |
быть получены выражения |
частот |
|
ных погрешностей |
для второго и |
|
третьего приближений |
|
|
sin ns |
|
|
0“," = 1 - |
2 п |
X |
ns |
||
|
2л |
|
|
X- |
I |
|
2 |
|
|
, |
|
|
1 —— sin2 |
|
|
V |
3 |
0,3
$ 0 SI
0,1
ЕГ
0,1
0 1 2 3 4 5 6 2L-H q~S
Рис. 23. Графики расчетной за висимости частотной погрешно сти от числа разбиений для моделей трех приближений
ns
-----
2 п
/
sin “X— |
|
|
|
|
|
|
|
2 |
п |
|
|
|
|
|
|
J I I = I |
/ |
|
|
|
|
|
|
ns |
2 |
( |
ns |
3 |
ns \ |
||
V ~2п |
1 — — |
[ sin2 -— + — |
sin4 — |
) |
|||
3 \ |
2n |
10 |
2n ) |
||||
На рис. 23 изображены зависимости |
Sm= 8 Ш(— ) для |
всех |
|||||
трех приближений. Из графика видно, |
|
s |
f \ Ч |
I |
числе |
||
|
что при заданном |
51
(4—6 ) разбиений на одну полуволну ошибка моделирования по спектру собственных частот для первого приближения не превы сит 8 —2%. Увеличение точности моделирования при использова нии звеньев (см. рис. 1 2 ) может быть достигнуто при помощи уве личения их числа. Теперь мы можем получить аналитическое вы ражение, которое ' связывает число разбиений (или число звеньев) на половину длины волны в функции заданной погреш ности. Для этого решим биквадратное уравнение (130)
q = b V
Получим, что число разбиений- |
|
|
|
п = — |
---- |
1 |
(131) |
3 |
/ 1 - / 1 - 4 8 Ил1 |
|
|
Очевидно, что если высшая гармоническая составляющая |
|||
процесса пом, где г — номер высшей гармонической |
составляю |
||
щей, и, если первая собственная |
частота моделируемой трубы |
||
соь то на частоте гсом по |
длине |
трубы уложится s —----- полу- |
|
|
|
|
ш, |
волн, следовательно, число разбиений п, исходя из выражения (131), составит
где 8 Ш ^ 1 — частотная погрешность на частоте гым.
Если в трубопроводной системе имеется несколько участков разной длины, то число разбиений по формуле (132) выбирается по самому короткому участку, так как спектр собственных ча стот более длинных участков ниже, чем у короткого, следова тельно, частотная погрешность моделирования других участков при той же длине участка разбиения будет заведомо ниже.
Получение выражений, аналогичных выражению (132), для электрических моделей поперечных колебаний трубы второго и третьего приближений представляется довольно сложным. Однако анализ кривых (см. рис. 23) и экспериментальные иссле дования показывают, что при заданной погрешности моделиро вания по частоте число звеньев второго приближения оказы вается на 25% меньше, а число звеньев третьего приближения в два раза меньше числа звеньев первого приближения.
§ 5. СРАВНЕНИЯ ЗВЕНЬЕВ ТРЕХ ПОРЯДКОВ ПРИБЛИЖЕНИЯ
Сравним модели всех трех порядков приближения при одном и том же числе моделирующих звеньев. Критерием для оценки эффективности звеньев послужит частотная погрешность модели
52