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

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

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

Поскольку матрица C ' симметричная, левая часть этого равенства равна

нулю и, следовательно,

'T

'

0

. Для исходных собственных векторов

 

и

1

2

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

условие ортогональности будет иметь вид:

(u 1 )T u 2

1T uT u

2

1T A

2

 

0 .

Собственные векторы ортогональны с матрицей А.

Последовательность решения задачи раскрытия антенны

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

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

2.Решаем упругую систему (число степеней свободы определяется для каждого конкретно примера) и находим собственные частоты и фор-

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

3.Если условие n nдоп выполняется с запасом, часть его можно взять на запас по фиксации, т. е. увеличит угловую скорость штанги в момент фиксации.

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

Далее рассмотрим пример другой возможной эквивалентной схемы ан-

тенного устройства в момент фиксации На практике бывает, что антенное устройство является конструкцией

сравнительно жесткой, тогда в этом случае в качестве эквивалентной схемы механической системы в момент фиксации антенного устройства в рабочее положение можно принять такую, где конструкция устройства – абсолютно твердое тело, а эквивалентная жесткость системы приведена к узлу фиксации в точке С, c {cX , cy , cZ , c x , c y , c z }T , рис. 5.23

111

 

A

 

 

 

 

y

A

 

 

 

 

A

-

ЦМ

 

 

ось

фиксатор

 

С

вращения

 

 

 

 

 

 

z

x

 

 

 

 

КА

 

 

Рис. 5.23 Схема установки антенны Тогда в момент фиксации устройства в рабочее положение, его в пределах

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

боды:

1.Движение центра масс вдоль оси 0Z

2.Вращение относительно оси 0Y

3.Вращение относительно оси 0X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, y }T . Кинетическая энер-

 

Обобщенными координатами будут q {z, X

гия системы будет равна:

 

 

 

 

 

T

1

( mz2

J

 

2

J

2

) ;

 

 

 

2

x

x

y y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а потенциальная:

 

 

 

 

 

 

 

П

 

1

( с

 

z 2

 

с

 

2

с

2

);

 

 

 

 

z

 

 

x

y

 

 

 

2

 

c

 

 

x

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

zc z

yc

x

xc

y .

Найдем матрицу

жесткостей. Поскольку

 

П

 

, получим:

 

 

 

 

 

 

Cq

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

112

 

 

 

cZ (z yC X

 

xC y )

 

 

 

 

 

 

 

 

 

 

 

 

 

Cq

cZ yC (z yC X

xC y ) c x X

 

 

 

 

cZ xC (z yC X

xC y ) c y y

 

 

cZ

 

 

yC cZ

 

xC cZ

 

 

z

 

 

 

 

 

 

 

c y

(c

x

c y2 )

 

c y x

 

 

x

 

Z

C

 

Z C

 

Z

C C

 

 

 

c x

( c x y )

c

y

c

x2

 

 

y

 

Z

C

 

Z

C C

 

 

Z C

 

 

Теперь имеем уравнение Лагранжа в стандартной форме:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cq

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Aq

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя подстановку (5.4.5), получим характеристическое уравнение

|

2 A

C |

 

0, или

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(m

2

 

c11 )

 

 

c12

 

 

c13

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c21

 

 

(J X

2

 

 

c22 )

c23

 

0 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

c31

 

 

 

 

c32

2

c33 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(J y

 

 

 

 

 

решая которое найдем собственные частоты

1< 2< 3,

а затем соответст-

вующие

собственные

векторы i ,i

 

1,2,3 и

общее решение

исходного

уравнения в виде (5.4.9).

 

 

 

 

 

 

 

 

 

 

Поскольку в начальный момент времени все обобщенные координаты

равны нулю, qi(0)=0, все ai

должны быть также нулевыми. Постоянные ci

найдутся из условий

 

 

 

 

 

 

 

 

 

 

 

 

 

 

z

0

 

 

u

1

u

2

u 3

 

c

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

3

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

2

3

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

c2

 

 

 

 

 

 

 

 

 

 

 

 

u2

u2

u2

 

2

, т.е. q(0) ux .

 

 

 

 

 

 

y 0

 

 

u

1

u

2

u 3

 

c

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

3

3

 

 

 

 

 

 

 

 

 

 

Здесь столбцы матрицы и являются столбцами нормированных собствен-

ных векторов.

Если

умножить это

 

равенство слева

на u-1,

получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u 1q(0)

x . Определив таким образом координатный столбец x , определим

и все ci. Будем иметь ci=xi/

 

i, i=1,2,3.

 

 

 

 

 

 

 

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

убедиться в ортогональности собственных векторов i ,i

1,2,3 (или нор-

мированных векторов u ( i ) ).

113

Теперь можно приступить к определению перегрузок в центре масс ан-

тенны. Если обозначить rA - радиус вектор центра масс системы относитель-

но ц.м. антенного устройства, то будем иметь следующее векторное равенст-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

во aA

a

 

rA

 

 

 

(

 

 

rA ), где a - ускорение ц.м. антенного устройст-

ва, а

- вектор угловой скорости вращения устройства. Учитывая, что

 

 

 

 

 

(

 

 

 

 

 

(

 

 

2

 

 

 

 

 

:

 

 

 

 

 

 

r )

 

 

r )

r

 

, запишем равенство для a

A

 

 

 

 

 

 

 

A

 

 

 

 

 

A

A

 

 

 

 

 

 

 

 

xA

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

xA

 

 

 

 

 

 

 

0

 

i

j

k

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

yA

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

x

y

0

 

 

y

(

x xA

y yA )

yA

(

x2

y2 )

 

 

 

zA

 

z

 

xA

yA

0

 

0

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

nx

ny

nz

1

g max(xA (t))

t

1

g max( yA (t))

t

1

g max(zA (t))

t

nXдоп

nдоп

(5.4.13)

y

 

nzдоп ,

Если эти неравенства не удовлетворяются, необходимо рассмотреть возмож-

ность уменьшения жесткосей системы и энергетики раскрытия антенного устройства.

5.4.4 Возможные особенности задачи

Решая уравнения Лагранжа (пусть для определенности в обратной форме

 

 

 

 

 

 

 

cq

0 ) с использованием способа Эйлера в виде q

 

sin( t

) ,

Aq

 

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

вания:

 

A

 

 

,

1/

2

(5.4.14)

 

 

 

 

Далее, решая характеристическое уравнение,

 

 

A E

 

0

 

 

(5.4.15)

 

 

 

 

и ссылаясь

на

опыт

практики, сказали, что все собственные числа

i , i 1...n , полученные из решения уравнения (5.4.14) будут вещественны,

114

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

нения будет иметь вид:

q

n B (i ) sin(

i

t

i

) ,

(5.4.16)

i

 

 

 

i 1

где i - столбцы собственных векторов уравнения (1), а i , i const.

В теории линейных преобразований доказывается, что если вещественная

 

 

матрица линейного преобразования Ax

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

числа вещественны (не обязательно различны), а собственные векторы, соот-

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

бой. Если же еще эта матрица удовлетворяет критерию Сильвестра (все глав-

ные миноры положительные), то собственные значения – положительны.

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

 

 

ние подобия Ex

x . Здесь все собственные числа равны между собой и лю-

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

После всего этого , возвращаясь к нашей матрице A мы видим, что она несимметричная, хотя каждая из матиц , A симметрична. Это, должно вы-

звать сомнение в том, что решение (5.4.16) справедливо всегда. Исследуем этот вопрос.

Заменим преобразование (5.4.14) ему эквивалентным. Для этого матрицу

А (в нашем случае она диагональная), запишем в виде A uT u , где u – диа-

гональная матрица

 

 

m1

 

0

 

0

 

 

u

0

 

m2

0

 

.

 

0

 

0

 

m3

 

Тогда будем иметь:

 

 

 

 

uT u

.

(5.4.17)

115

Теперь умножим равенство (5.4.17) на u слева: u

uT u

u , u

', и

u uT '

'

(5.4.18)

 

Преобразование (5.4.18) симметричное, и, следовательно, собственные

значения

i

будут вещественны, а если учесть, что

- положительна и пре-

 

 

 

образование u uT - преобразование подобия, следовательно матрица u uT

также положительна, и собственные значения вещественны и положительны.

Таким образом, почти доказано, что решение (5.4.16) справедливо. Да-

лее рассмотрим случай кратных корней.

 

 

Из (5.4.18) следует, что если i

j , то

'i

' j 0 . Следовательно пер-

воначальные собственные векторы

будут

ортогональны с «весом», т.е.

i A j 0 .

 

 

 

В общем случае, решая характеристическое уравнение (5.4.15), среди его

корней могут оказаться и кратные. Пусть например k

k 1 . Тогда, в этом

случае будут равны и частоты

k

k 1

соответствующих колебаний.

 

 

 

Из теории обыкновенных дифференциальных уравнений известно, что

решение, соответствующие корню второй кратности, будет иметь вид:

( k

k t)sin(

k

t

k

) .

 

 

 

 

На этом основании Лагранж, когда исследовал малые движения и писал свою знаменитую «Аналитическую механику» полагал, что в этом случае движение будет неустойчивым т.к. при t , член k t также неограни-

ченно возрастает.

Эта ошибка держалась долгое время (более 50 лет) и именно благодаря авторитету этого выдающегося ученого. И только уже после его смерти Вей-

ерштрасс обнаружил, что это не так. В середине 19-го века он доказал, что в этом случае имеет место тождественность 2х уравнений системы (5.4.14) и,

следовательно, корень характеристического уравнения (5.4.15) фактически остается первой кратности и тогда отпадает решение, соответствующее крат-

ному корню а вместе с тем и заключение Лагранжа. Он пришел к этому вы-

116

воду на основании теории квадратичных форм, которыми являются кинети-

ческая и потенциальная энергии системы. Если повторить рассуждения Вей-

ерштрасса, то оно выглядит следующим образом. Две квадратичные формы

 

1 T

 

1 T

 

 

T

 

q Aq, П

 

q Cq

из которых хотя бы одна определенно-

2

2

положительна (в данном случае это кинетическая энергия) одним линейным преобразованием могут быть приведены к каноническому виду. Если гово-

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

скании главных осей поверхности второго порядка или, например, главных направлений тензора инерции.

Известно, что одну форму – определенно-положительную можно привес-

ти к каноническому виду многими способами. Один из них связан с заменой

исходного базиса e 0 . Т. е., если найти собственные векторы линейного пре-

 

 

 

образования Ax

x

и из них построить новый базис e ', то в этом новом

базисе¥ матрица А станет диагональной, диагональные элементы которой равны собственным числам i , i 1...n , а форма примет вид

 

1

 

1

0

0

 

 

 

 

 

 

T

T

... ... ...

 

.

2

 

 

 

 

0

0

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

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

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

 

 

 

 

 

 

 

 

 

 

 

K

 

, где К – верхнедиагональная матрица вида:

щую подстановку: q

 

 

 

k12 ...

k1n

 

 

 

1

 

 

K

0

1 ...

k2 n

,

 

 

0

0

1 ...

 

 

 

0

0

0

1

 

 

 

 

 

 

 

 

 

117

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

p1

0

0

 

 

 

 

в этом случае

T

T

K

T

AK

 

, A' K

T

AK ,

A'

 

0

...

0

.

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

0

pn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Поскольку

A' также как и

A симметрична и

det K

 

 

det K T

1, все

главные

 

миноры

 

 

матриц

A'

 

 

и

 

A

 

равны.

 

 

Следовательно

p

 

,

 

a , p p

 

a11

 

a12

 

 

 

, p

 

2

, p p p

 

, p

 

 

3

и т.д.

 

 

 

 

 

 

 

 

 

 

 

 

 

1

1

 

 

 

 

 

 

 

2

 

 

3

 

 

 

 

1

 

11

1

2

 

a21

 

a22

 

 

2

 

1

2

3

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

2

 

 

 

Чтобы матрица

A' была диагональной,

остается последовательно подоб-

рать коэффициенты Kij.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Наконец, способ, заключающийся в замене матрицы А на произведении 2х

треугольных матриц (если А вещественна и положительна – это возможно),

т.е. A uT u , u верхнедиагональная матрица.

Пусть Вейерштрасс, также как и мы сейчас воспользовался эти способом.

 

Тогда

 

 

 

 

 

 

 

uq,

 

 

1 T

 

1 T

 

T

 

T

 

q

Aq

 

q

u

 

u q . Если ввести новые переменные

2

2

 

 

 

 

 

uq, q

u 1

, qT

T (u 1 )T , получим :

 

1

 

 

 

 

1

T

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

 

 

T

 

1

T

 

1

 

*

 

1

T

 

1

 

T

 

 

E

 

, П

 

q

C q

 

 

(u

 

)

C u

 

, C

 

(u

 

)

C u

 

..

2

 

 

2

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Матрица С* является также как и С симметричной матрицей, т.к. она по-

лучена преобразованием подобия. А раз так, то для нее можно найти собст-

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

ском виде. Матрица Е так и остается прежней, а матрица C* преобразуется в диагональную с коэффициентами, равными собственным значениям преобра-

 

 

 

зования C* x

x . Действительно, пусть исходный базис

e 0 . Запишем век-

торы базиса в виде матрицы – строки:

 

0

 

0

 

0

 

, где еѐ элементы – векто-

e

,e

,...e

 

 

1

2

n

 

 

 

 

 

 

 

 

 

 

 

 

 

ры. Тогда в этом «старом» базисе «новый» базис e '

из нормированных соб-

ственных векторов можно представить так:

 

 

 

 

 

 

 

 

 

Г ,

e'

,e'

2

,...e'

n

 

e 0

,e 0

,...e 0

 

1

 

 

 

1

2

n

 

118

или проще: e' e 0 Г , где Г – матрица, столбцами которой являются коорди-

натные столбцы «новых» векторов в «старом» базисе – направляющие коси-

нусы.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Любой вектор x в этих базисах будет иметь свой столбец. Пусть в базисе

 

0

 

 

 

 

 

 

{

1 , 2 ,...

T

 

 

{ '1 ,

'2 ,...

T

 

 

 

 

e

 

это

 

 

n } ,

а в базисе e '

'

'n } .

Поскольку

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

это

 

e

0

,e

0

,...e

0

 

 

...

 

 

 

 

 

 

 

 

0 Г

 

 

 

 

 

 

 

 

 

 

 

2

n

 

 

 

 

e

и вместе с этим x

e'

' , то учитывая, что e '

e

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

будем

 

иметь

 

связь

между столбцами

e 0

e'

',

или

e 0

e 0 Г

' ,

т.е.

Г' .

Теперь имеем

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

T

T

E

 

 

'T

Г

T

 

'

и

 

2

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

T

 

 

T

*

 

 

 

1

 

T

 

 

**

 

П

 

'

 

 

Г

 

C

Г '

 

 

'

 

C

 

'.

2

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Преобразование Г это преобразование поворота, т.е. преобразование ор-

тогональное. Итак, в новом базисе (базисе из собственных векторов) имеем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

..

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

T

 

'

 

E

', т.к.

Г

 

EГ E, П

 

 

C

 

 

.. ..

..

-

1

 

T

 

Т

1

'

T

**

', где C

**

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

0 ..

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

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

Значит все i - положительные числа.

Теперь уравнения колебаний в этих новых координатах (их называют

нормальными) будут иметь простейший вид.

 

 

 

E ' C**

' 0

(5.4.19)

119

Эти колебания называются главными. Как видно из (5.4.19) эти уравнения определяют независимые друг от друга движения в скалярной форме

'

k

k

'

k

0, k 1...n

 

 

 

 

 

 

 

 

 

 

 

Их решение имеет вид: 'k

k sin( k t

k ), k

 

 

k . Таким обра-

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

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

Отсюда с особенной ясностью следует заключение Вейерштрасса – если несколько корней характеристического уравнения между собой равны (на-

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

ковые частоты. При этом имеет место только явление унисона и не происхо-

дит возрастание координат. Ясно также, что в системе уравнений (5.4.19) два

будут тождественны и, следовательно в системе уравнений (5.4.14).

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

зор инерции к каноническому виду, т.е. когда моменты инерции твердого те-

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

нейного преобразования (два момента инерции) равны между собой.

В заключение убедимся, что частоты главных колебаний, найденные из

(5.4.19) это именно те частоты, которые мы бы определили из исходного

уравнения (5.4.14).

 

 

 

 

 

Имеем C**

2 , или

Г(u 1 )T Cu 1 Г

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

u 1 Г слева получим u 1 ГГ Т (u 1 )T Cu 1 Г

2u 1 Г , или т.к. '

u 1 Г

1

 

'

2

'. Еще

умножим последнее

равенство на А

получим

A C

 

C '

 

 

2 A ' или

'

A '

2 . Это равенство с точностью до обозначения

120