Анализ волнового уравнения и расчет собственных чисел и собственных функций для колебаний давления в трубе при наличии осевого градиента температуры
Рефераты >> Физика >> Анализ волнового уравнения и расчет собственных чисел и собственных функций для колебаний давления в трубе при наличии осевого градиента температуры

Наш интеграл имеет вид:

(55)

Определим а0 из условия симметрии:

p(0)=-p(1) (56)

Подставляя условие (56) в уравнение (48), получим:

или

(57)

Отсюда

(58)

Перепишем (55) в следующем виде:

(59)

Подведем итоги и запишем вид функционала (45) после всех проведенных вычислений:

(60)

Получили квадратичную форму коэффициентов аi. Теперь получим собственные функции и собственные значения при ε=0,1. В этом случае (60) примет вид:

(61)

Исследуем функционал на минимум. Для этого найдем производные по а1, а3, а5.

(62)

На минимуме производные (62) должны обращаться в нуль, поэтому получаем систему однородных уравнений относительно а1, а3, а5.

(63)

Чтобы система (63) имела решение, ее определитель должен быть равен нулю.

(64)

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

(65)

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

Раскрывая определитель, получим:

, 1, 2, 3,…6. (66)

Вычислим Сi:

Сложив все Сi, получим уравнение:

(67)

Разделив уравнение (67) на коэффициент при последнем члене, получим:

(68)

Алгебраическое уравнение вида (68) можно записать как

, (69)

где

Делая в уравнении (69) замену неизвестного

, , (70)

получаем, так называемое, приведенное уравнение

, (71)

где

Вычислим дискриминант:

Здесь , поэтому уравнение (71) имеет три действительных решения. Положим

, ,

, .

Тогда решениями приведенного уравнения будут:

(72)

Из (70) найдем :

Отсюда

(73)

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

Для определения чисел а1, а3 , а5 подставим значения в систему уравнений (63), получим:

(74)

Разделим каждое уравнение на коэффициент при последнем члене:

Теперь рассмотрим уравнения попарно, найдем коэффициенты а1, а3, а5 для каждой пары уравнений, а затем проверим точность, подставив решения каждой пары в неиспользованное уравнение.

1)

Проверим точность:

2)

Проверим точность:

3)

Проверим точность:

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

Коэффициент а1 определим из условия нормировки (42):

Вычисления остальных коэффициентов дают:

(из (57))


Страница: