Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Синтез_мех_систем

.pdf
Скачиваний:
43
Добавлен:
16.03.2015
Размер:
2.5 Mб
Скачать

труда. Можно также и обеспечить заданное ограничение на время поворота,

которое, обычно, не является жестким.

Однако, учитывая, что антенна после приведения ее в рабочее положение должна выполнять свою задачу – быть работоспособной то она не должна быть подвержена большим ударным нагрузкам, действующим на неѐ при фиксации штанги в рабочем положении. Чтобы оценить эти нагрузки не дос-

таточно использовать только уравнения механики твердого тела, необходимо использовать аппарат механики сплошных сред.

Задача проектирования системы раскрытия и фиксации антенны в рабо-

чем положении будет заключаться в подборе энергопривода и, возможно,

жесткостных характеристик стержня и узла фиксации, обеспечивающих две главные задачи – гарантированную фиксацию стержня в рабочее положение и не превышение при этом допустимых ударных нагрузок на конструкцию антенны.

Во многих случаях область допустимых угловых скоростей, при которых происходит фиксация, бывает достаточно узкой, когда обеспечиваются эти требования.

Используют обычно электромеханические, гидромеханические или про-

сто пружинные энергоприводы. Далее предполагаем, что для раскрытия ан-

тенны используется пружинный привод, как наиболее надежный, дешевый и простой.

Очевидно, что на этапе формирования структурной схемы системы рас-

крытия необходимо иметь сравнительно простой аппарат анализа работоспо-

Асобности предлагаемых схем. Причем, как обычно бывает,

 

анализ надо проводить за сравнительно короткий срок.

В

Поэтому на этой стадии проектирования, задачу поведе-

ц.м.

ния упругой системы с бесконечным числом степеней сво-

 

 

боды, заменяют на задачу поведения так называемых приве-

 

денных, или эквивалентно-приведенных систем с конечным числом степеней свободы.

101

Например, груз, подвешенный к неподвижной точке А на пружине АВ,

если учитывать распределенную массу пружины представляет собой систему с бесконечным числом степеней свободы. Но когда масса груза значительно превышает массу пружины, при нахождении наименьшей (основной) частоты колебаний без большой погрешности можно пренебречь массой пружины,

сохраняя ее свойства упругости. Если, кроме того, предположить, что груз совершает прямолинейные колебания, то рассматриваемая система обраща-

ется в приведенную систему с одной степенью свободы. Для определения движения такой системы достаточно найти только одну величину как функ-

цию времени – отклонение ц.м. груза от положения равновесия.

Рассмотрим более сложный пример. При определении основных частот поперечных колебаний горизонтальной балки, шарнирно закрепленной на 2х

опорах, массу, распределенную по тому или иному закону вдоль длины бал-

ки, заменяют несколькими сосредоточенными массами, величина и положе-

ние которых определяется особым способом.

q1

q2

q3

q4

m1

m2

 

m4

 

 

 

 

m3

Рассматривая малые колебания балки, пренебрегают незначительными горизонтальными смещениями масс. При таких условиях для определения

поперечных перемещений приведенной системы достаточно найти в функ-

102

ции времени только четыре величины q1, q2, q3, q4. Ими будут определятся прогибы отдельных масс при колебаниях (скорости, ускорения и, как следст-

вие, перегрузки), а вместе с тем и форма упругой линии балки в наибольшем ее отклонении. Здесь мы имеем случай замены упругой системы эквивалент-

ной системой с 4мя степенями свободы.

Математическим аппаратом решения таких эквивалентных систем являет-

ся теория малых колебаний систем с конечным числом степеней свободы и уравнения Лагранжа второго рода в прямой и обратной форме. Их еще назы-

вают уравнениями в усилиях и в перемещениях, соответственно. Для боль-

шинства колеблющихся систем более простым является подход с использо-

ванием уравнений в в усилиях, что требует вычисление матрицы жесткости,

но имеется много случаев, когда удобнеее противоположный подход, когда построенине матрицы жесткости представляет некоторые затруднения, в то время как матрица податливости вычисляется сравнительно простио. Полу-

чим эти уравнения.

5.4.2 Уравнения Лагранжа IIго рода

Известно, что в основе теории малых колебаний систем с конечным чис-

лом степеней свободы лежит понятие кинетической и потенциальной энер-

гии, которые являются определенно – положительными квадратичными фор-

мами от обобщенных скоростей и координат, соответственно. Эти энергии в общем виде имеют вид.

 

 

 

 

T

 

 

q

A q 2T (q), q

C q 2T (q) ,

T

 

 

 

 

 

где: А, С – квадратичные, симметрические матрицы с постоянными коэффи-

циентами инерции (аik) и жесткости (cik) соответственно.

Поскольку считается, что система малых колебаний консервативна, то

Е=Т+П=const и отсюда получают уравнения Лагранжа второго рода.

d

(T

 

 

П )

 

0, или в скалярной фоме :

dt

 

 

 

 

 

 

 

 

 

 

 

d

(

 

T

)

 

П

0, i 1,2,3...n

 

 

 

 

 

 

dt

 

 

q

i

q

i

 

 

 

 

 

 

 

103

где: Q

П

`- обобщенная сила.

 

 

i

qi

 

 

 

 

 

Обобщенную силу, трактуя как n – мерный вектор можно записать и так -

 

 

 

Q

Cq.

 

 

 

 

 

 

 

Обозначим

P - вектор обобщенной силы, который нужно приложить к

системе, чтобы удержать ее в равновесом положении, когда обобщенные ко-

 

 

 

 

 

0 .

 

 

 

 

ординаты не раны нулю, т.е. q

0 , следовательно Pq 0

Для этого должны выполняться условия равновесия:

 

 

 

 

 

 

 

Q

P 0, откуда следуе что P

Q, или P

Cq (5.4.1)

Поскольку определитель матрицы С не равен нулю, то систему (5.4.1)

можно разрешить относительно q

 

 

C 1

 

q

P,

(5.4.2)

Это уравнение представляет собой обобщенный закон Гука для системы с n степенями свободы. Его также можно рассматривать как условие равнове-

 

 

 

сия между силой P и восстанавливающей силой Q в так называемой обрат-

ной форме.

- это матрица коэффициентов aik, которые называются стати-

ческими коэффициентами влияния. Эта матрица, так же как и матрица С симметрическая, aik= aki, что отражает принцип взаимности перемещений.

Если в уравнении (5.4.2) положить все силы кроме Рк равными нулю, а Рк=1,

то получим что:

qi

ik , i

1...n

 

 

 

 

 

 

q1

 

11

12 .

1,n 1

1n

 

P1

 

 

 

 

q2

 

21

21 .

21

2 n

 

P2

 

.

 

.

.

.

.

.

 

.

 

qn 1

 

an 1,1

n 1, 2 .

n 1,n 1

n 1,n

 

Pn 1

 

qn

 

an1

n , 2 .

n ,n 1

nn

 

Pn

104

Это значит, что коэффициенты aik определяют конфигурацию равновесно-

го положения системы, когда к ней будет приложена только одна единичная сила Pk. Иначе говоря, aik – это изменение координат qi от единичной силы,

соответствующей координате qk. Поскольку уравнение (5.4.2) выражает ус-

ловие равновесия между внешними силами P и восстанавливающими сила-

ми Q , возникающими при отклонении системы из положения устойчивого равновесия, в качестве этих сил можно взять обобщенные силы инерции, т. е.

 

d

 

T

 

 

 

 

 

 

 

 

 

 

Pi

 

 

 

,i

1...n, или F

Aq. (5.4.3)

dt

 

qi

 

 

 

 

 

 

 

 

 

 

Если учесть условие равновесия Q

P 0 и подставив сюда

P через

 

 

 

 

 

силы инерции P

 

 

 

Aq , получим принцип Даламбера:

 

 

 

 

 

 

 

Q 0 , отсюда следует уравнение Лагранжа второго рода.

Aq

Подставив (5.4.3) в систему (5.4.2) получим уравнение движения в обрат-

ной форме по отношению к уравнению Лагранжа второго рода.

 

 

 

 

q

 

A, q

 

Aq, D

Dq. (5.4.4)

Необходимо отметить, что даже если кинетическая энергия выражается в канонической форме, матрица D уже не будет симметрической.

Чтобы решить систему (5.4.4) необходимо найти собственные формы,

частоты и коэффициенты влияния, как видно, потенциальная энергия систе-

мы нам здесь не понадобится и еѐ не надо определять.

Конкретный вид решения получим применительно к нашей задаче – кон-

сольной балки с сосредоточенной массой на свободном конце.

Пусть штанга представляет собой стержень с равномерно распределенной массой и изгибной жесткостью по длине (EJ=const, (l)= o=const) и сосре-

доточенной массой m3 на конце.

105

o

 

x

l

m3

 

Для простоты массу штанги представим двумя сосредоточенными масса-

ми m

o l

, m

 

o l

 

2

 

1

 

1

 

1

 

 

m1

m2

 

x

l/4

 

l/4

m3

 

 

Тогда эта механическая система будет иметь три степени свободы

 

{q , q

, q

}T .

 

 

 

q

Кинетическую энергию запишем в канонической форме

 

 

1

 

2

 

 

 

3

 

 

 

 

T

3

 

1

m q

2

, а матрица В, B=D будет иметь вид:

 

 

 

 

 

 

i 1

2

i

 

 

 

i

 

 

 

 

 

 

 

 

 

 

 

 

m1

0

0

 

 

 

 

 

 

 

 

 

 

 

B

 

 

 

ij

 

 

 

0

m2

0

, i, j 1,2,3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

m3

 

Поскольку систем консервативна, решение уравнения (5.4.4) ищется в

гармоничской форме.

 

 

)

 

q

sin( t

(5.4.5)

106

 

 

 

 

получим 2 B

 

 

Подставив (5.4.5)

в (5.4.4)

 

, или, вводя единичную

матрицу E и обозначив p2

1

 

, получим:

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

(B p2 E)

0,

или B

 

p2 .

(5.4.6)

 

 

 

 

 

Это задача о собственных векторах - и числах - p2 матрицы B.

 

 

 

 

 

Чтобы собственный вектор был ненулевым,

0 , должен быть равен

нулю определитель:

 

 

 

B p2 E

 

(5.4.7)

 

 

 

 

 

 

 

 

 

Раскрывая его, получим кубическое уравнение относительно p2. На прак-

тике, для подобного рода задач корни различные и положительные. Их мож-

но отыскать любым приближенным способом. В результате получим 3 собст-

венные частоты 1 , 2 , 3 , 1

2

3 .

 

 

 

 

Подставив в (5.4.6) найдем 3 собственных вектора

(1) ,

( 2 ) ,

( 3) .

Поскольку определитель системы (5.4.7) равен нулю, а ранг матрицы двум (корни векового уравнения различны) собственные векторы будут оп-

ределены с точностью до произвольной постоянной.

В скалярной форме уравнение (5.4.6) имеет вид

(b

 

p2 )

1

b

2

b

3

11

 

 

12

13

b

1

(b

 

p2 )

2

b

3

21

22

 

 

23

b

1

b

2

(b

 

p2 )

3

31

32

33

 

 

0

0 (5.4.8)

0

Отбросим

 

последнее

уравнение, из оставшихся первых 2х, определим

1 , 2 через

 

3 , значение которого произвольно.

Обозначив миноры определителя J31, J32, J33, получим:

 

J31

 

 

,

 

J32

 

3 , 3 1

1

J33

3

2

J33

 

 

 

 

 

107

Итак, для каждого значения 1 , 2 , 3 , и получим три вектора

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1) ,

( 2 ) ,

 

( 3) . Ими будут определяться собственные формы колебаний, каж-

дая из которых соответствует собственной частоте.

 

 

 

 

 

Обычно, формы приводят к безразмерному виду, т.е. нормируют.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Нормирующий множитель (в нашем случае) N

2

2

2

 

1

1

2 .

Безразмерные формы соответственно равны:

 

 

 

 

u

 

 

1

,u

 

2

,u

 

3

, u 2

u 2

u 2

1

 

 

 

 

1

 

 

2

 

3

 

 

 

 

 

 

 

N

 

N

N

1

2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Теперь можно записать общий интеграл уравнения (5.4.4)

 

3

 

 

 

 

 

 

 

 

 

 

t

 

), A ,

 

const.

(5.4.9)

q

 

A u ( i ) sin(

i

i

i

 

i 1

i

 

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку q(0)

 

0 все

i

0. Постоянные Ai найдутся из с учетом на-

чальных условий по скоростям:

q

(0)

 

 

l

,

q

(0)

 

 

3l

,

q

(0)

 

l

(5.4.10)

k 4

 

 

1

 

 

 

2

 

 

k

4

 

3

 

k

 

 

k - угловая скорость штанги в момент фиксации. Аi можно определить из линейного уравнения:

 

A

1

 

u

1

u

2

u

3

 

 

1

 

 

1

 

1

 

1

 

 

 

 

 

1

 

2

 

3

 

 

A2

 

, u

 

 

 

 

q(0) u

2

u2

u2

u2

(5.4.11)

 

A

3

 

u

1

u

2

u

3

 

 

3

 

 

3

 

3

 

3

 

После определения Ai, i=1,2,3 можно будет найти перегрузку в ц.м. ан-

тенны и сравнить еѐ с допустимой:

n

1

max

 

3 A u(i )

2

sin(

 

t)

 

n

(5.4.12)

 

 

 

 

i

i

 

g

t

 

i 3

 

 

 

доп

 

 

 

i 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

5.4.3 Определение коэффициентов влияния условия ортогональности

Определение основывается на понятии упругой линии балки, находя-

ik

щейся под действием заданных сил и моментов в статическом равновесии. Из курса сопротивления материалов известно, что если балка под действием со-

108

средоточенных сил и моментов находится в равновесии, то уравнение еѐ уп-

ругой линии выглядит следующим образом, [2]:

EJy(x) EJ o x EJfo

M

(x a)2

P

(x b)3

2!

3!

 

 

 

 

 

 

 

 

 

y

 

M

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

b

 

 

P

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При защемленном левом конце прогиб

fo и угол

o равны нулю.

Для нашего примера для определения

ik приложим поочередно единич-

ную силу в точках сосредоточенных масс mi , i 1, 2, 3

и в соответствии с урав-

нением упругой линии вычислим перемещения в этих точках. Заметим, так-

же, что матрица перемещений коэффициентов влияния симметричная, т.е.

aij=aji.

y

 

 

 

М 0

m1

m2

x

 

 

m3

 

l1

 

 

 

 

R

 

l2

P=1

 

 

 

 

 

l3

 

Опорная реакция в заделке и момент при приложении единичной силы в точке mi, i=1,2,3 ,будут равны Ri=mP, Mi=Pl.

Перемещения в точек mi при приложении едеиничной силы в этих точках (т.е. a11, a22, a33), будут равны:

109

 

Pl l2

 

Pl3

 

2Pl3

 

 

l3

 

i i

 

i

 

i

;

a

 

i

 

ii

2EJ 6EJ 6EJ

ii

 

3EJ

 

 

 

Поскольку матрица ||

ij

|| симметричная, ее коэффициенты ik , i k удоб-

нее вычислять следующим образом. Приложив единичную в т. m3 вычислим

13 и 23, и приложив единичную силу в т. m2 вычислим 12. Получим:

 

M (l 0)2

P(l 0)3

Pl l 2

 

 

Pl3

 

 

l 2

(3l l )

 

 

 

 

1

 

 

 

1

 

 

 

3

1

 

 

 

 

1

 

1

3

1

 

 

 

 

13

2 EJ

 

 

 

 

6 EJ

 

 

 

 

2 EJ

 

6 EJ

 

 

 

6 EJ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M (l

2

0)2

 

 

P(l

2

0)3

 

 

l2 (3l l

2

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

23

2 EJ

 

 

 

 

6 EJ

 

 

 

6 EJ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M (l 0)2

P(l 0)3

P l

2

(l 2

0)2

 

 

P(l 0)3

 

l 2

(3l

2

l )

 

1

 

 

 

1

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

1

 

1

12

2 EJ

 

 

 

 

2 EJ

 

 

 

 

 

 

2 EJ

 

 

 

 

 

 

 

 

6 EJ

 

 

6 EJ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Проверку полученных результатов (возможных ошибок при вычислениях)

произведем по условиям ортогональности собственных форм. Из теории ли-

нейных преобразований известно, что если матрица линейного преобразова-

 

 

 

 

ния симметричная, а векторы x и

y собственные и соответствуют собствен-

ным значениям p1

и p2 , p1 p2 , то имеет место условие ортогональности этих

 

 

0 . Для нашего случая условие ортогональности получает-

векторов, т.е. x

y

ся следующим образом. Матрицу А в уравнении (5.4.4) представим в виде

1 1

произведения двух матриц A uT u , uT u A2 , где A2 - диагональная матрица

с диагональными элементами, равными квадратному корню из соответст-

вующих элементов матрицы

А. Тогда уравнение (5.4.4) запишется в виде

 

u

 

 

.

Умножив

 

его

на

матрицу

u и

 

 

 

q

uq

 

обозначения q '

uq , получим

 

T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

где матрица C '

u

u - симметричная, т.к. она получена преобразо-

q '

C 'q

 

 

 

 

 

 

 

 

 

T

 

 

 

 

 

 

 

ванием подобия. Для двух собственных векторов ' ,

' и соответствующих

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

2

 

им

собственным

значениям

p2

, p2

, p2

p2

преобразования,

аналогичного

 

 

 

 

 

 

 

 

 

 

 

1

2

1

2

 

 

 

(5.4.6) будем иметь

 

 

 

 

 

 

 

 

 

 

 

 

c ' '

 

p2 ' , c '

'

p2 ' .

 

 

 

 

 

 

 

 

 

1

 

 

1

1

2

2

2

 

 

 

 

 

 

 

 

 

 

Умножив первое равенство на

' , а второе на ' и вычитая первое из вто-

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

1

 

 

рого получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

'T C '

'

 

'T C ' '

( p2

p2 ) 'T ' .

 

 

 

 

 

 

 

1

 

2

 

2

1

 

2

1

1

2

 

 

 

 

 

 

110