Мельник В.М., Карачун В.В.

Національний технічний університет України «КПІ»

ЛІНІЙНО-ПРУЖНА МОДЕЛЬ ДІЇ ЗВУКОВОЇ ХВИЛІ НА ПЛАСТИНУ

 

При нестаціонарній пружній взаємодії, диференціальне рівняння збуреного згинного руху пластини можна записати у вигляді:

,                          (1)

де ;  – час.

Праву частину рівняння задамо наступним чином

,                                      (2)

де .

Розв’язок рівняння (1) шукаємо у вигляді:

.                                             (3)

Підставимо (3) у вихідне рівняння (1):

,                                        (4)

де .

Відшукуємо наближені розв’язки у формі:

.                                          (5)

Стовпець

                                               (6)

підлягає визначенню.

Підстановка (5) в рівняння (4) приводить до наближеної рівності

,                                    (7)

для  якої стовпець (6) вважаємо за найбільш слушний у тому розумінні, що проекції лівої і правої частин виразу (7) на лінійну оболонку  образів координатних функцій були би рівними.

Помножимо обидві частини рівності (7) на

.

Одержуємо:

.                       (8)

Матриця Грама образів координатних функцій має вигляд:

.

Матрицю Грама координатних функцій  по енергії оператора . Запишемо наступними  чином–

,                                     (9)

.                 (10)

Знаючи матрицю Грама образів координатних функцій, а також матрицю Грама  координатних функцій по енергії оператора  (10), систему (8) можна навести інакше –

.                                                  (11)

Якщо , то ця система однозначно розв’язувана, тобто

.                                              (12)

У зв’язку з тим, що матриця  неособлива (невироджена), тому

,

де  – одинична матриця. Отже,

;

.

Але  є многочлен ступеня  відносно . І якщо  виявиться додатним коренем цього рівняння, тоді система (11) може стати нерозв’язуваною.

 Проаналізуємо цю тезу докладніше. Припустимо для простоти, що рівняння однорідне –

,                                             (13)

де , при однорідних граничних умовах.

Ненульові розв’язки шукаємо у вигляді

,                                             (14)

де множник  задовольняє граничним умовам (13). Підставляючи (9.41) в рівняння (13), одержуємо:

;   ,    де .                 (15)

.                                                   (16)

Наближений розв’язок задачі (15) шукаємо у вигляді

                                   (17)

де  – координатні функції.

Таким чином, після підстановки виразу (17) у співвідношення (13), отримуємо наближену рівність –

                                         (18)

з невідомим стовпцем .

Як і в попередньому, найкращим вважається стовпець, при якому проекції лівої і правої частини виразу (18) на лінійну оболонку  образів координатних функцій дорівнюють одна одній, тобто

.

Стовпець  не повинен бути нульовим, тому

.

З додатної визначенності оператора  походить, що  і, отже,  постають додатними коренями многочлена . В нашій задачі вони являються наближеними власними числами оператора .  

Приймаючи до уваги, що , з виразу (16) з’ясовуємо:

.

Це означає, що у прийнятому наближенні власні частоти обчислюються за формулами –

.

Отже, для того, щоб задача визначення вимушених коливань, задача (4), була розв’язуваною, необхідно вилучити співпадання частот щільності акустичного випромінювання і власних. Іншими словами, запобігти частотному резонансу.

Якщо обмежитися тільки першими трьома координатними функціями

тоді матриця  виявиться діагональною

.                         (19)

Нулі її визначника будуть спостерігатися за умови, що

; 

з одночасним виконанням співвідношення

.