Файл: Козобков, А. А. Электрическое моделирование вибраций трубопроводов.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

 

 

 

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