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

Толстых-Уравнматем-физики

.pdf
Скачиваний:
153
Добавлен:
27.03.2016
Размер:
5.07 Mб
Скачать

u x,t

sin x at sin x at

 

1

 

x at

 

 

 

 

 

dy

 

 

 

 

 

 

 

 

 

 

2a

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

x at

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

sinxcosat

y

 

x at .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2a

 

 

 

x at

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Таким образом, форма струны определяется функцией

 

 

 

 

u x,t sin xcosat t.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u x,

 

sinxcos

 

 

 

 

 

 

 

 

u x,

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

2a

 

2a

2 2a

 

 

 

 

 

2a

 

 

то есть в момент

времени t

 

 

 

 

струна

будет

параллельна оси Ox

 

2a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(рис. 9).

Рис. 9

Задачи для самостоятельного решения

1.Неограниченная струна возбуждена начальным локальным отклонением вида, изображенного на рисунке 10. Начальные скорости отсутствуют.

Изобразить положения струны в моменты времени tk

kc

,

k

 

(см.

1,5

 

 

2a

 

 

 

задачу 6.1). Найти решение методом Даламбера.

 

 

 

2. Начертить профиль неограниченной струны для различных моментов времени, если начальное отклонение ее от положения равновесия равно нулю: f x 0, а начальная скорость постоянна F x v0 на участкес; 2c и равна нулю вне этого участка (ср. задачу 6.1).

3.Неограниченная струна возбуждена локальным начальным отклонением, имеющим форму параболы (рис. 11). Найти закон u x,t движения точек

51

струны в любой момент времени, если начальные скорости точек струны равны нулю.

4. Методом Даламбера найти решение волнового уравнения ut"2 a2u"x2

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

а) a 2, u x,0 0, u't

t 0

x;

 

 

 

 

 

б) a 1, u x,0 еx, u't

 

 

4x;

 

 

 

 

 

 

 

t 0

 

 

 

 

 

в) a 3, u x,0 х х 2 ,

ut'

t 0

еx.

 

 

 

 

 

 

Рис. 10

Рис. 11

5.Решить задачи 1, 2 для полуограниченной струны 0 х при условии, что конец x 0 закреплен.

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

0,

0 x l,

2l x ,

 

 

 

 

x

 

 

 

 

 

 

 

u x,0

 

 

 

 

 

 

 

 

sin

 

,

l x 2l.

 

 

l

 

 

 

 

 

 

 

l

 

l

 

Определить графически форму струны в моменты времени t

,

,

 

 

 

 

 

 

4a

a

5l , 3l , 7l , 9l , предполагая, что начальные скорости отсутствуют. 4a 2a 4a 4a

7.Линейные вибрационные системы

сраспределенными параметрами

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

52

виброзащитных систем. Основным элементом виброзащитной системы является амортизатор, а его наиболее существенной частью – упругий элемент. В результате внутреннего трения в упругом элементе происходит демпфирование колебаний. В ряде конструкций амортизаторов применяют специальные демпфирующие устройства для рассеивания энергии колебаний. Выбор модели упругого элемента определяется спектром частот вибрационного воздействия. Чем они выше, тем более сложной оказывается модель упругого амортизатора. Если частоты возбуждения больше собственных частот амортизатора, его необходимо рассматривать в виде системы с распределенными параметрами. В этом случае продольные упругие колебания амортизатора описываются волновым уравнением относительно смещения u x,t плоскости его поперечного сечения:

 

2u

2

u

c

2

 

2u

,

(7.1)

t2

t

 

x2

 

 

 

 

 

где коэффициент распределенных потерь; с – скорость распространения смещения вдоль амортизатора

 

 

 

 

 

 

 

 

 

 

 

с

 

 

Е

,

 

 

 

 

 

 

 

 

 

 

 

(7.2)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Е – модуль упругости, – плотность материала амортизатора.

Пользуясь методом Фурье, решение (7.1) ищем в виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

u x,t X x T t .

 

 

 

 

 

 

 

 

 

 

(7.3)

Подставляя (7.3) в (7.1), получим XT

 

2 XT

 

c

2

 

 

 

 

 

 

 

 

X T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 T

 

2

2

T 0,

T

2

T

X

 

 

 

 

 

T

 

 

 

c

 

 

 

 

 

 

2

const 0

 

 

 

 

 

 

 

 

 

 

 

 

 

c2T

 

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

X

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

 

 

X 0,

 

 

(рассмотреть самостоятельно возможность const > 0).

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

 

t

Acos

 

2

 

2

t Bsin

 

2

 

2

t ,

 

T e

 

 

 

 

 

(7.4)

 

 

 

 

 

 

 

 

 

 

 

 

 

X Ccos x Dsin x,

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

53

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

Рассмотрим сначала свободные колебания амортизатора.

Если концы амортизатора свободны, то граничные условия имеют

вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

u

 

 

 

 

0,

u

 

 

 

 

 

0,

 

(7.5)

 

 

 

 

 

x

 

 

 

x

 

 

 

 

 

 

 

x 0

 

 

 

 

 

 

 

 

x l

 

 

 

 

 

 

 

 

 

 

 

 

 

где l – длина амортизатора.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как u x,t 0, то, учитывая (7.3), имеем X 0 X

l 0.

Из (7.4) X Csin x Dcos x.

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

D = 0, а для

нахождения имеем характеристическое уравнение sin l 0, откуда

 

n

n

,

 

n 1,2,...

 

(7.6)

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

и значения собственных частот амортизатора

 

 

 

nc

 

n

 

 

 

 

 

 

 

 

 

n nc

 

 

 

 

E

,

n 1,2,...

 

(7.7)

 

l

 

 

 

 

 

 

 

l

 

 

 

 

 

 

 

 

 

Таким образом, решение (7.1) с учетом (7.4), (7.3), соответствующее свободным колебаниям амортизатора, определяется рядом:

 

 

nx

 

 

 

 

 

 

u x,

t e t cos

(an cos

n2 2

t bn sin

n2 2 t).

(7.8)

 

 

n 1

l

 

 

 

 

 

Коэффициенты an, bn определяются начальными условиями.

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

(рис. 12).

Закон колебаний описывается уравнением (7.1) при 0.

Пусть М – амортизируемая масса, s – площадь поперечного сечения амортизатора, m sl – масса амортизатора. Тогда граничные условия имеют вид:

u

 

x 0

0,

M

2u

Es

u

,

x l.

(7.9)

 

 

 

2

x

Рис. 12

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

54

Первое условие, примененное к (7.4), приводит к равенству C 0, а из второго при 0, получим

u TX Acos t Bsin t sin x, u T cos x,

 

 

2u

 

 

 

 

x

 

 

 

2TX,

M 2 sin l Es cos l,

 

 

t2

2

 

 

E

 

M sin l scos l,

 

c2 2

2

l

 

 

 

 

 

 

 

 

 

tg sl m .

M M

Итак,

Рис. 13

tg

1

,

 

M

,

(7.10)

 

 

m

l

l

 

m

,

 

Es

.

(7.11)

 

 

 

 

c

 

 

l

 

Решение трансцендентного уравнения (7.10) изображено на рисунке 13 и имеет вид:

n 1 n 2n 1

 

,

n 1,2,...

 

2

 

 

(7.12)

или

n

n 1 n,

0 n

 

 

,

lim n

0.

(7.13)

 

 

 

 

2

 

n

 

 

Если велико (легкий амортизатор нагружен тяжелой массой), то из (7.10), (7.13) в первом приближении

1

 

1

 

.

(7.14)

 

 

 

1

 

 

 

 

 

 

 

 

Действительно, при 0 справедливы соотношения эквивалентности:

tg , 1 tg 2.

Отсюда частота собственных колебаний первой формы

1

 

c

1

 

 

 

 

.

 

1

 

 

1

 

 

 

 

(7.15)

l

 

M

 

 

 

 

m

 

 

55

Более точное приближение получим, если в (7.10) положим

 

 

 

 

 

 

 

 

 

1

3

 

 

 

 

 

 

 

1

2

 

 

 

 

 

 

tg

 

 

1

 

 

.

 

 

 

 

 

3

3

 

 

 

 

 

1

1

 

 

1

 

 

 

 

1

 

 

1

 

Тогда 2

 

1

2

 

tg

1

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

1

 

3 1

 

1 1

 

 

1

 

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

откуда второе приближение

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

1

 

 

 

 

 

2

1

(7.16)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

3

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и второе приближение частоты собственных колебаний первой формы:

(2)

 

c

 

(2)

 

 

 

 

 

 

 

 

 

 

.

(7.17)

l

 

 

1

 

1

 

M

m

 

 

 

 

 

 

 

 

 

 

 

3

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

Если мало (амортизатор нагружен легкой массой), можно считать

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

 

 

 

 

 

 

 

 

tg

 

 

 

 

 

 

 

 

 

ctg

2

 

 

 

 

2

 

2

 

2

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

,

0,

 

 

 

 

 

 

 

 

 

 

 

tg

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

откуда

 

 

 

 

 

 

– малая величина,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2 1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

1 ,

 

 

 

 

(7.18)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

 

 

1

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

M

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

.

(7.19)

 

 

 

 

 

 

 

 

 

 

 

m

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

1

 

 

 

 

 

 

 

 

 

m

 

 

m

 

 

56

При 0 из (7.10) получим характеристическое уравнение

cos cos l 0,

(7.20)

соответствующее случаю, когда один конец амортизатора закреплен, а второй свободен. Из (7.20)

n

 

2n 1

,

 

n 0,1,2,...

(7.21)

 

 

 

 

 

 

2l

 

 

n

nc

2n 1 c

,

n 0,1,2,...

(7.22)

 

 

 

 

 

2l

 

 

При произвольном значении для в первом приближении можно брать средние значения между крайними пределами, определяемыми формулами (7.14) и (7.18), соответствующими M , M 0.

Параметр 2 определяется как волновое число, где длина c

продольной волны, бегущей по амортизатору. Значения волнового числа, определяемые формулами (7.6) и (7.21), соответствуют волновому резонансу, когда на длине амортизатора укладывается целое число полуволн. На частотах, соответствующих указанным значениям (7.7) и (7.22), виброизоляция отсутствует. Вблизи этих частот она может быть даже отрицательной из-за возникновения значительных продольных колебаний амортизатора. Следует выбирать такие соотношения параметров системы (рис. 12), чтобы избежать подобной ситуации. Аналогично можно рассмотреть вынужденные колебания амортизатора

(см. [19]).

Заметим, что математические модели (см. [2]) подвижной рамы грузового вагона, кузова вагона с подвижной хребтовой балкой и т.п. подобны рассмотренным в этом параграфе примерам.

8.Задача о распространении тепла в стержне

иее решение методом Фурье. Интеграл Пуассона

8.1.Решение уравнения теплопроводности методом Фурье

Рассмотрим тонкий однородный стержень длины l, расположенный между точками x = 0, x = l оси .

Боковая поверхность стержня предполагается теплонепроницаемой, сечение настолько мало, что температуру во всех точках сечения можно

57

считать постоянной. В начальный момент времени t 0 дано распределение температуры вдоль стержня, характеризуемое функцией f(x). Указан также тепловой режим на концах стержня. Задача состоит в определении температуры точек стержня как функции абсциссы x и времени t : u u x,t . Можно показать (см., например, [11], гл. XVIII, § 4; [4], гл. 5),

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

удовлетворяет уравнению теплопроводности

u

a2

2u

,

(8.1)

t

x2

 

 

 

если F x,t плотность тепловых источников, то

u

a2

2u

f x,t ,

f x,t

F x,t

,

t

x2

 

 

 

 

c

a2 k , k – коэффициент теплопроводности стержня; с, ρ – c

соответственно, теплоемкость и плотность вещества стержня. В пространстве уравнение теплопроводности имеет вид:

u

2

2

 

2u

2u

2u

 

a u a

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

2

t

 

 

x

y

z

.

 

 

 

 

 

 

 

 

 

Сначала рассмотрим простейший случай, когда концы стержня поддерживаются при температуре, равной нулю.

Таким образом, нужно найти решение уравнения (8.1),

удовлетворяющее начальному условию

u x,0 f x

(8.2)

и граничным условиям

u l,t 0.

 

u 0,t 0,

(8.3)

В более общем случае первая краевая задача для уравнения теплопроводности имеет граничные условия

u 0,t 1 t , u l,t 2 t .

(8.3')

Согласно методу Фурье, решение ищем аналогично п. 4, в виде:

u x,t X x T t .

58

Подстановка такого вида решения в (8.2) дает

 

2

X

 

 

T

X

2

const 0

 

 

 

X

X x T t a

 

x T t a2T

 

 

T a2 2T 0 T Ce a2 2t ,

 

X 2 X 0

 

 

X Acos x Bsin x,

A2 B2 0;

граничные условия имеют вид: u 0,t X 0 T t 0

X 0 0 A 0; u l,t X l T t 0 X l 0

 

sin l 0 n

n

,

n 1,2,...

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

l

 

 

 

 

 

 

Полагая n

и BC bn , получим ряд частных решений:

 

 

 

 

 

 

a n

2

 

nx

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

u x,t bne

l

 

 

sin

,

(8.4)

 

 

 

 

 

 

n 1

 

 

 

 

 

 

 

 

 

 

 

l

 

причем (см. (8.2))

 

nx

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f x bn sin

,

 

 

 

0 x l ,

 

 

 

 

 

 

 

 

 

 

 

n 1

 

l

 

 

 

 

 

 

 

откуда, если функция f(x) разложима в ряд Фурье,

 

 

 

2 l

nx

 

 

 

 

 

 

 

 

bn

 

f x sin

 

 

dx,

n 1,2,...

(8.5)

 

l

 

l

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Аналогично решается задача при других граничных условиях. Например, если на концах стержня поддерживается постоянная температура, отличная от нуля:

u 0,t A,

u l,t B,

(8.6)

то решение

 

 

u x,t v x,t u* x,t ,

(*)

где u* x,t удовлетворяет данным граничным условиям, например,

59

u* x,t A

B A

x.

(8.7)

 

 

l

 

Тогда v x,t удовлетворяет задаче типа (8.1), (8.3), (8.2):

 

v x,0 f x u* x,0 .

(8.2')

то есть ищется по формулам (8.4), (8.5).

Таким образом, если в концах стержня поддерживается постоянная температура (отличная от нуля), то температурный режим стержня определяется функцией:

 

 

 

 

B A

 

 

a n 2

 

 

 

 

nx

 

u x,t A

 

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

x bne l

 

 

sin

,

 

 

 

 

 

 

 

 

 

l

 

 

n 1

 

 

 

 

 

 

 

 

l

 

 

2 l

 

 

 

 

B A

 

 

nx

 

 

 

bn

 

 

f (x) A

 

 

 

x sin

 

 

 

dx

l

l

 

l

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

l

 

nx

 

 

 

2

A 1 n B .

 

 

 

 

 

 

f x sin

dx

 

 

 

 

l

 

 

 

 

 

 

 

0

 

 

 

l

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.8)

(8.9)

При t распределение температуры внутри стержня стремится к линейному u* x,t u* x .

Задача 8.1. Дан тонкий однородный стержень длины l с теплоизолированной поверхностью, начальная температура которого равна

cx l x

f x l2 . Концы стержня поддерживаются при температуре равной нулю. Определить температуру стержня в любой момент времени.

Решение. Требуется решить задачу типа (8.1), (8.2), (8.3) при

f x

cx l x

 

. Решение ищется согласно (8.4), (8.5):

l2

 

 

 

 

 

 

 

 

 

 

 

8.4

a n 2

nx

 

 

 

 

 

 

 

t

 

 

 

 

 

 

 

 

 

u x,t bne

l

 

sin

.

 

 

 

 

 

 

 

n 1

 

 

 

 

l

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

60

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]