Скачиваний:
120
Добавлен:
15.03.2015
Размер:
6.6 Mб
Скачать
k = Pδ ,

находится по формуле типа (8.14)

где P - сила, δ - деформация.

Используя, например, известную формулу [32] для прогиба в центре прямоугольной пластины, свободно опертой по краям при действии в центре сосредоточенной силы P :

δ =

 

 

 

 

4Pa3b3

 

 

 

 

 

,

 

 

π4 D(a2 +b2 )2

 

 

найдем

 

 

 

P

 

 

 

 

π4 D (a2

+b2 )2

 

k

 

 

 

 

 

 

 

 

i

=

 

i

 

=

 

 

 

i

 

 

i

 

i

.

(8.17)

δi

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4ai3bi3

 

Сравнивая формулы (8.16) и (8.17), получим

 

 

 

 

 

m =

mпл.i

 

 

.

 

 

 

 

(8.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

4

 

 

 

 

 

 

 

 

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

D α

2

 

 

 

 

 

 

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

=

 

i

i

.

 

 

 

(8.19)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

4a2ξ

i

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

Пример 8.1. Для уменьшения амплитуд резонансных колебаний двух параллельно расположенных ячеек между ними в центре установлена ДВ. Размеры ячеек 2×100×200 мм и 1×100×200 мм, ДВ-10×10×20 мм. Материал плат - стеклотекстолит с характеристиками: Е=1 106 Па; ρ=2000 кг/м3; ν=0,3. ДВ выполнена из пенополиэтилена с характеристиками: Е=1 106 Па, η=0,5. На каждой ячейке установлены ЭРЭ, равномерно расположенные по площади платы, с общей массой 0,04 кг. Рассчитать коэффициенты передачи в центре каждой ячейки.

Решение. Расчет сосредоточенных параметров модели. Находим цилиндрические жесткости плат:

 

 

 

EH 3

 

2 1010 (2 103)3

 

D =

1

 

=

 

=14,7 Н/м;

 

 

 

 

1

 

12(1 − ν2 )

 

12(1 0,32 )

 

 

 

 

 

D

=

 

EH23

= 2 1010 (2 103)3

=1,84 Н/м.

 

 

2

 

12(1 − ν2 )

 

12(1 0,32 )

 

 

 

 

 

Считая ячейку свободно опертой по краям, найдем

α1 = 9,87(a12 + b12 ) = 9,87(0,122 + 0,112 ) = 49,35.

По формуле (8.19) рассчитаем приведенные значения коэффициентов жесткости

200

k =

D α2

=

14,7 49,352

1 1

 

=111877 Н/м;

1

4a2ξ

 

 

 

4 0,22 2

 

 

 

 

k2 =

D α2

=

1,84 49,352

2

1

=13985Н/м.

 

4a2ξ

 

 

4 0,22 2

Приведенные массы ячеек рассчитаем по формуле (8.18)

m np.i= m4i (i =1, 2) ;

где

m1 = ρabH1 + mexp = 2000 0,2 0,1 2 103 + 0,04 = 0,08 + 0,04 = 0,12 ; m2 = 2000 0,2 0,1 1 103 + 0,04 .

Получим: m = 0,12

= 0,03; m = 0,08

= 0,02.

np.1

4

 

 

 

 

np.2

4

 

 

 

 

 

 

 

 

Коэффициент жесткости ДВ найдем по формуле (8.15)

kä

=

ES

=

1 106 (0,01)

2

= 5000 Н/м.

H

 

0,02

 

 

 

 

 

 

 

Расчет собственных частот колебаний. Используя формулу (8.3), получим

ω

2

111877 +5000

+

13985

+5000

 

±

 

=

2 0,03

2

0,02

 

1,2

 

 

 

 

 

1

111877 +5000

+

13985 +5000

2

 

 

0.5

 

 

4

 

2 0,03

 

 

2 0,02

 

 

 

 

 

 

 

 

 

 

 

 

 

 

=

±

 

111877 13985

 

+111877 5000 +...13985 5000

 

 

 

 

 

 

 

 

 

 

 

 

 

0,03 0,02

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= (2,37 ±1,4) 106.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω01 =1940 рад/с; ω02 = 986рад/с.

 

 

 

Расчет коэффициентов передачи. Проведем вспомогательные расчеты:

 

 

 

kä

= 0,045;

 

kä

 

= 0,36 ; ω012 =

k1

 

= 3,7 106 (рад/с)2;

 

 

 

 

 

 

 

 

k1

 

k2

m1

 

 

 

 

 

 

 

 

 

ω12

 

 

 

 

 

 

 

ω2 =

k2

 

= 0,7 106 (рад/с)2;

 

= 5,4 ;

 

 

 

 

 

m

 

ω2

 

 

 

 

 

02

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

02

 

 

 

 

201

Рис. 8.4. Зависимость коэффициентов передачи ячеек от жесткости ДВ

µ

20

10

0

 

 

 

 

 

 

 

 

 

 

ω2

 

 

 

ω2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

=1,02;

1

 

 

= 5,49;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ω2

 

ω2 ω2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

01

 

 

01

02

 

 

 

 

 

 

 

 

 

 

 

 

ω2

 

 

ω2

 

 

 

ω2

 

 

 

ω4

 

 

 

 

 

 

 

 

2

=

0,26 ;

2

 

=1,39;

02

=

0,19 ;

 

2

 

= 0,36.

 

 

 

 

 

 

 

 

ω2

 

ω2

ω2

 

 

 

ω2

 

 

 

 

ω2

 

 

 

 

 

 

 

 

 

 

01

 

 

02

 

 

01

 

 

 

01

 

02

 

 

 

 

Подставляя найденные значения в формулы (8.9а)-(8.10а), получим

 

 

 

 

 

 

 

 

 

 

 

(1,405 5,4)2 +

 

 

 

 

 

 

0.5

 

 

 

 

 

 

 

 

 

 

 

.

 

µ1

{1,405

1,02[1+

0,36 +5,29(1+ 0,045)]+5,49}2 +

=15,1;

=

 

 

 

 

 

 

+

0,52 0,4052

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

..

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ 0,25[0,405

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

(0,36 1,02 + 0,045 5,4)]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1,405 0,26)2 +

 

 

 

 

 

 

0.5

 

 

 

 

 

 

 

 

..

 

µ2

{1,405 1,39[1+

0,045 + 0,19(1+ 0,36)]+ 0,36}2 +

 

=4,55.

=

 

 

 

 

 

 

+

0,25 0,4052

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

+ 0,25[0,405 (0,36 0,26 + 0,405 1,39)]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= 0,08 + 0,04 = 0,12

 

 

 

 

 

 

 

 

 

 

На рис. 8.4 приведены зависимости

 

 

 

 

 

 

 

 

 

коэффициентов передачи µ1 и µ2 от жест-

µ1

 

кости демпфера kд. Видно, что, подбирая

µ2

величину kд, можно получить минимальное

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

значение µ1 или µ2 . Наиболее часто тре-

 

 

 

 

 

 

 

 

 

буется уменьшить амплитуды резонансных

 

 

 

 

 

 

 

 

 

колебаний обеих

ячеек, т.е.

µ1 = µ2 . В

1 kдопт 3 4 kд. 4,

данном конкретном примере это достигается

 

 

 

 

 

 

 

 

 

при k

Ä

=1,8 10

4 и µ = µ

2

= 4,8 .

 

 

 

 

 

 

 

 

 

 

 

1

 

Динамический метод.1) Приве-

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

1) Способ определения приведенных параметров механической системы: Патент №2006717 РФ, F16F15/22/. Долгов Г.Ф., Талицкий Е.Н. (РФ).

202

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

При кинематическом возбуждении уравнение движения системы с одной степенью свободы (рис. 8.5)

) && &&

)

)

(8.20)

m(Ζ −ΖΑ ) + kZ (1

+ jη) = 0,

где знак « » означает принадлежность к системе с сосредоточенными параметрами (ССП).

Считая ускорение основания Z&&Α изменяющимся по гармоническому

закону Z&&Α (t)= Z&&Αe jωt , из (8.20) найдем:

(t)

 

 

 

)

 

ω2 Ζ

 

 

 

 

Z (t) = −

)

− ω2

Α

)

.

 

(8.21)

 

ω2

+ jηω2

 

 

 

0

 

 

 

0

 

 

 

Перемещение ячейки в области резонанса в точке с координатами

X,Y

 

ω2 ΖΑ (t)Kif (x, y)

 

 

Ζ(x, y,t)= −

,

(8.22)

ω02if

− ω2 + jηω02if

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

Сравнивая выражения (8.21) и (8.22), нетрудно заметить, что при равенстве частот ω0if и ω0 есть два пути приведения КМП:

1. Считать, что

Ζ(x, y,t)= Ζ(t) .

Тогда

(ω02if −ω2 ) [1Kif (x, y)]+ jηωoif2

)

η =

 

.

 

 

jKif (x, y)ω02if

а)

б)

Рис. 8.5. Модель ячейки

203

2. Принимать, что η = η. В этом случае:

 

Ζ(x, y,t)= Ζ(t)Kif (x, y),

(8.23)

где Ζ(t) находится по формуле (8.21), а коэффициент формы колебанийKif (x, y) - методами, изложенными в гл. 3.

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

Приведение массы и жесткости ячейки. С о б с т в е н н а я ф о р м а к о л е б а н и й я ч е й к и и з в е с т н а. При приведении распределенных параметров к сосредоточенным необходимо выполнение условия: энергия ССП должна равняться энергии системы с распределенными параметрами (СРП) [3].

Кинетическая (Tс) и потенциальная (Пc) энергии ССП равны:

 

1 ) )

2

 

 

1

) )

2

 

 

Tc =

2 mZ

 

(t) ;

Пс =

2

ΚZ

 

(t) .

(8.24)

После подстановки выражения (8.21), получим

 

 

 

 

 

 

) 4

2

Α(t)

 

 

 

 

T = 1

 

 

 

mω

Z

 

;

 

(8.25)

)

 

 

 

))

2

 

c

 

2

 

 

 

 

 

 

 

 

 

(ω2

−ω2 + jηω2 )

 

 

 

 

 

 

0

 

 

 

0

 

 

 

 

П

c

=

1

 

 

Κω4Z 2Α(t)

 

 

.

(8.26)

)

 

 

))

 

2

 

 

 

 

 

 

 

 

 

 

2 (ω2

− ω2 + jηω2 )

 

 

 

 

 

 

 

 

0

 

 

 

0

 

 

 

Кинетическая энергия СРП описывается выражением

T= 1 ∫∫m0 Z 2 (x, y,t)dx .

2 S

После подстановки выражения для перемещения ячейки Z (x, y,t)

при кинематическом возбуждении (3.55), (3.59) и ряда преобразований получим

 

 

 

 

 

 

 

∫∫

if

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 m

 

 

 

 

 

 

W

(x, y)ds

 

 

 

T =

ω4Ζ2

(t)

 

 

S

 

 

 

 

 

 

 

. (8.27)

 

 

 

 

 

 

 

 

 

2

 

2 0

 

Α

 

(ω2

 

 

 

2

 

W 2

 

 

 

 

 

 

 

−ω2 + jηω2 )

 

∫∫

(x, y)ds

 

 

 

 

 

 

oif

 

 

 

oif

 

if

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

Потенциальную

энергию

пластины

(3.37)

после

подстановки

Z (x, y, t) можно привести к виду: 204

 

D

ω4 ΖΑ2 (t) ∫∫Wif

(x, y)ds 2

 

 

 

 

(S )

 

 

2

П =

 

 

 

 

 

 

 

∫∫[W′′(x)+W′′(y)] ds. (8.28)

2 (ω02if

−ω2 + jηω02if )2 ∫∫Wif2 (x, y)ds

2

 

 

S

 

 

 

 

 

S

 

 

 

Приравнивая выражения (8.25) и (8.27), а затем (8.26) и (8.28), найдем, что сосредоточенная масса

 

 

 

∫∫Wif

2

 

 

)

 

 

(x, y)ds

 

 

 

S

 

 

 

m = m

 

,

(8.29)

0

 

∫∫Wif2 (x, y)ds

 

 

 

 

S

а сосредоточенная жесткость

 

 

2

 

)

∫∫Wif (x, y)ds

 

S

 

 

Κ = D

 

.

(8.30)

∫∫[Wif′′(x) +Wif′′( y)]2ds

S

Рассмотренный метод приведения легко применим для однородных прямоугольных пластин с классическими способами типа жесткого защемления, свободного опирания или незакрепленного края. Собственная форма колебаний W(x,y) таких пластин описывается выражением (3.39).

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

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

а)

б)

Рис. 8.6. Модель ячейки с дополнительным инерционным элементом

но моделирует колебания ячейки в точке с координатами X,Y, то максимальные кинетические (Тмакс) и потенциальные (Пмакс) энергии колебаний ячейки должны быть равны максимальной кинетической и потенциальной

205

энергии системы с одной степенью свободы [3], т.е.

 

ΚΖ2

)

ω2

 

 

 

2

 

 

Пмакс =

макс

; Тмакс =

макс

0

,

(8.31)

2

2

 

 

 

 

 

 

где Ζмакс = Ζ(x, y,t) / Кфг (x, y) ;

Если на ячейку в точку с координатами (x, y) установить дополнительный инерционный элемент массой m2 (рис. 8.6, а), то кинетическая энергия такой системы опишется формулой

 

 

 

 

 

 

)

ω2

m Ζ2

(x, y)ω2

T

= Т

 

+ Т

 

=

2

 

 

макс

+

г

,

смакс

 

пмакс

 

гмакс

 

2

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

где Ζмакс = Ζ(x, y,ω) / kфг (x, y) ;

Тпмакс – максимальная кинетическая энергия ячейки; Тгмакс – максимальная кинетическая энергия дополнительного инер-

ционного элемента; ω– круговая резонансная частота ячейки с дополнительным инер-

ционным элементом.

Если масса дополнительного инерционного элемента mГ много меньше сосредоточенной массы ячейки, то Кф(x,y)фг(x,y). Потенциальная энергия ячейки при установке точечного дополнительного инерционного элемента остается без изменения. Учитывая, что при резонансе Пмаксмакс для ячейки без дополнительного инерционного элемента и Псмакссмакс для ячейки с дополнительным инерционным элементом, правомерно записать систему уравнений

 

)

2

 

 

 

 

 

 

 

 

 

 

 

 

К mω = 0,

 

 

 

 

)

)

2

2

K

2

2

.

К mω

= m

ф

(x, y)ω

 

 

г

 

 

Решение системы относительно неизвестных K и m дает выражения

Κ) =

mгKф2

(x, y) f02 f02г 4π2

;

(8.32)

 

f02 f02г

 

 

 

 

 

 

 

 

 

)

m K 2

(x, y) f

2

 

 

 

г

ф

 

 

 

m =

 

 

 

 

,

 

(8.33)

 

f

02 f02г

 

 

 

 

 

 

 

 

 

где f0 и f– собственные частоты колебаний ячейки и ячейки с дополнительным инерционным элементом соответственно.

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

Определение сосредоточенной жесткости демпфирующей вставки.

206

Физическая модель конструкции, содержащая упругий элемент, установленный между основаниями ячейки и корпуса, может быть представлена в виде системы с одной степенью свободы (рис. 8.7). Потенциальная энергия конструкций с упругим элементом Пс складывается из потенциальной энергии ячейки П, рассчитываемой по выражению (8.31), и потенциальной энергии упругого элемента Пд, рассчитываемой по формуле [3]

= ΚДΖ2 (x, y,t)

Д2 .

 

kД

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а)

 

 

 

 

б)

Рис. 8.7. Модель ячейки с дополнительным упругим элементом

После подстановки перемещения ячейки в точке крепления упругого

элемента, учитывая (8.24), получим

(x, y)]2

[Κ) + Κ)ДkфД (x, y)]Ζ)2

ПС = П + ПД =

Κ)Ζ)2

 

ΚД [Ζ)kфД

 

 

+

 

 

=

 

.

2

2

 

2

 

 

 

 

 

 

Аналогично получим выражение для потенциальной энергии системы с одной степенью свободы

ПС = (Κ + ΚД )Ζ2 ,

2

где Ζ находится из выражения (8.21). Если КД << Κ, то kфд (x,y)kф(x,y). Условия равенства энергии будет выполняться только в случае, когда

ΚД = ΚДkф2 .

(8.34)

Жесткость упругой вставки находится из геометрических размеров и модуля упругости 1-го рода Е материалов вставки по формуле

КД =

ЕS

,

(8.35)

H

 

 

 

где S и H – площадь поперечного сечения вставки и ее высота соответственно.

207

Рис. 8.8. Физическая модель n - ячеечной конструкции

8.3. КОНСТРУКЦИИ ЭЛЕКТРОННЫХ СРЕДСТВ С n ЯЧЕЙКАМИ1)

Если в блоке РЭС установлено n ячеек с n-1 ДВ, физическую модель можно представить в виде, показанном на рис. 8.8.

Колебания ячеек в этом случае описываются системой n уравнений

) &&

)

 

 

(Z

 

 

&&

 

 

)

 

 

 

P

+ P

= 0;

mZ

1

K

 

A

Z

1

)(1 + jη )

 

 

1

 

 

 

 

1

 

 

Д0

 

 

Д1

 

 

 

) &)&

)

 

 

(Z

 

&)&

)

 

 

 

) P

 

+ P

 

= 0;

mZ

2

K

2

A

Z

 

2

)(1 + jη

2

 

 

 

 

 

 

 

 

 

 

 

 

Д1

 

Д2

 

 

) &)&

)

 

 

(Z

 

&)&

 

)

 

 

) P

 

+ P

 

= 0;

mZ

3

K

3

A

Z

3

)(1 + jη

3

 

 

 

 

 

 

 

 

 

 

 

Д2

 

 

Д3

 

 

 

− − − − − − − − − − − − − − − − − − − − − − −

 

) &)&

)

 

 

(Z

 

 

&)&

 

 

)

)

P

 

+ P

 

= 0;

mZ

i

K

i

A

Z

i

 

)(1 + jη

 

 

 

 

 

 

 

 

 

i

 

 

 

Дi-1

 

 

Дi

 

 

 

− − − − − − − − − − − − − − − − − − − − − − −

 

) &)&

)

 

 

(Z

 

&)&

)

 

 

) P

 

 

+ P

 

 

= 0;

mZ

n

K

n

A

Z

n

)(1 + jη

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Дn-1

Дn

 

где n- количество ячеек в блоке;

}

i=1,n -текущий номер ячейки;

Zi -перемещение инерционных элементов;

ZA- перемещение основания;

PДi -силы, действующие со стороны i-й ДВ. Указанные силы находятся по выраже-

нию

PÄi = KÄi (Zi Zi+1)(1 + jηÄi )

при 1 i n 1, а PД0 = РДn = 0.

Таким образом, существует система уравнений, позволяющая рассчитать колебания n ячеек с ДВ, у которых известны сосредоточенные параметры.

Если блок РЭА содержит двух- и трехъячеечные конструкции с ДВ, то система n уравнений разбивается на несколько подсистем из двух и трех уравнений. Так, подсистема уравнений, описавших колебания трехъячеечной конструкции с ДВ при кинематическом возбуждении, будет выглядеть

1) Написан совместно с Долговым Г.Ф.

208

 

 

 

 

 

 

 

) &&

 

 

 

 

 

 

)

 

 

(Z

 

 

 

 

&&

 

 

 

 

 

 

 

)

 

 

 

P

 

+ P

 

 

= 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mZ

1

 

K

 

A

Z

)(1 + jη )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

Д0

 

 

 

 

Д1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

) &)&

 

 

 

 

 

 

)

 

 

(Z

 

 

 

 

 

&)&

 

 

 

 

 

 

 

)

 

 

 

) P

 

 

 

+ P

 

 

= 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mZ

2

K

2

A

Z

2

)(1 + jη

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Д1

 

 

Д2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

) &)&

 

 

 

 

 

 

)

 

 

(Z

 

 

 

 

 

&)&

 

 

 

 

 

 

 

)

 

 

) P

 

 

+ P

 

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

mZ

3

K

3

A

Z

3

)(1 + jη

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Д2

 

 

 

Д3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

PД1 = KД1 (Z1 Z2 )(1 + jηД1 ),

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

PД2

= KД2 (Z2 Z3 )(1 + jηД2 ).

 

 

 

 

 

 

 

 

 

 

 

 

Так как свободное колебания масс описывается гармоническими

 

функциями Zi

= Zi (t) = Zi e j(ωt−ϕ) , получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

ω

2 )

 

 

 

) )

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

) )

 

 

 

 

 

 

 

)

 

 

 

 

) )

 

 

 

 

)

 

)

 

 

m

Z

 

 

+ K Z (1

+ jη )

K Z

 

 

(1 + jη )

+ K Z (1 + jη

Ä1

 

 

1

 

 

1

 

 

1 1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1 A

 

 

 

 

 

 

1

 

 

 

 

Ä1 1

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

)

 

 

 

 

 

 

)

 

 

 

 

 

 

 

= 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

KÄ1Z2 (1 + jηÄ1)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

)

 

 

) )

 

 

 

 

 

 

 

 

)

 

 

 

)

 

 

) )

 

 

 

 

)

 

 

 

 

) )

 

 

 

 

)

 

 

)

m

 

ω2Z

2

+ K

2

Z

2

(1 + jη

2

 

K

2

Z

A

 

(1 + jη

2

) + K Z

 

(1 + jη

Ä1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ä1 2

 

 

 

 

 

 

 

 

)

 

 

 

)

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

)

 

 

 

)

 

 

 

 

 

 

 

)

 

 

 

 

)

 

 

 

)

 

 

 

 

 

)

 

 

 

 

 

 

 

 

KÄ1Z2 (1 + jηÄ1)

+ KÄ2Z2 (1 + jηÄ2 ) KÄ2Z3

(1 + jηÄ3) = 0,

 

 

)

 

 

 

2 )

 

 

) )

 

 

 

 

 

 

 

 

)

 

 

)

 

 

) )

 

 

 

 

 

)

 

 

 

 

 

)

 

)

 

 

 

 

)

 

 

 

)

m

ω Z

3

+ K

3

Z

3

(1 + jK

3

K

3

Z

A

(1 + jη

 

) + K

Ä2

Z

3

(1 + jη

Ä2

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

)

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

KÄ2Z2

 

(1 + jηÄ2 ) = 0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

После подстановки в эту систему уравнений выражения Kд = Kдkф

и деления всех членов на ZА(t) получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

+ K

Φ

 

(1 + jη ) K

(1 + jη ) + K

 

k

2

Φ

 

(1 + jη

 

)

 

 

 

 

 

m ω2Φ

1

 

Д1

1

Д1

 

 

 

 

 

1

 

 

 

 

 

1

1

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

Φ1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

KД1kΦ21Φ1 (1 + jηД1 ) = 0,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

m2ω2Φ2 + K2Φ2 (1 + jη2 ) K2 (1 + jη2 ) KД1kΦ21Φ1 (1 + jηД2 )

 

 

KД1kΦ2 2Φ2 (1 + jηД1 ) + KД2 kΦ2

2Φ2 (1 + jηД1 ) KД1kΦ2 3Φ3 (1 + jηД2 ) = 0,

)

 

 

 

 

 

)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

)

 

 

(1 + jη3 ) KД1kΦ2 3Φ3 (1 + jηД2 )

 

 

m3ω2Φ3

+ K3Φ3

(1 + jη3 ) K3

 

 

KД2 kΦ2 2Φ2 (1 + jηД2 ) = 0.

где Фi = Zi(t)/ZA(t) при i=1, 2, 3.

Получилась система уравнений с тремя неизвестными Фi, которая легко решается на ЭВМ при любом значении частоты ω.

КП ячейки в точке с координатами (х,у) рассчитывается по формуле

µi = Фik(x,y)k(x,y)+1 .

209