Ankilov
.pdfб) Задана интенсивность теплового потока через торцевое сечение x a:
Ka y (a) qa , |
Ka |
K(a). |
В частности, |
если стержень теплоизолирован при |
||
|
0 . |
|
|
|
|
|
x a , то y (a) |
x a |
|
|
|
|
|
в) На конце |
имеет место |
теплообмен с окружающей |
средой |
|||
известной температуры Ta : |
|
|
|
|
||
|
|
|
|
|
Ka K(a). |
|
|
|
Ka y (a) αa y(a) Ta , |
|
|||
Здесь a – коэффициент |
теплообмена |
на |
конце x a. Последнее |
условие |
(условие Ньютона) означает, что тепловой поток, передаваемый в единицу времени с единицы площади поверхности в окружающую среду, пропорционален разности температур поверхности тела и окружающей среды.
Аналогичные краевые условия могут быть заданы и на правом конце стержня при x b . Например, условие теплообмена при x b имеет вид
Kb y (b) b ( y(b) Tb ).
Втаблице 1.1 приведены возможные варианты краевых условий для
определения стационарного распределения температуры в стержне согласно уравнению (1.7).
Таблица 1.1
Варианты краевых условий для уравнения (1.7)
x |
x a |
|
|
x b |
||
№ |
|
|
||||
|
|
|
|
|
|
|
1 |
y Ta |
|
|
y Tb |
||
2 |
y T |
|
|
Kb y qb |
||
|
a |
|
|
|
|
|
3 |
y Ta |
Kb y b ( y Tb ) |
||||
4 |
Ka y qa |
|
|
y T |
||
|
|
|
|
|
|
b |
5 |
Ka y qa |
|
|
Kb y qb |
||
6 |
Ka y qa |
K |
b |
y |
b |
( y T ) |
|
|
|
|
b |
||
7 |
Ka y a ( y Ta ) |
|
|
y Tb |
||
8 |
Ka y a ( y Ta ) |
|
Kb y qb |
|||
9 |
Ka y a ( y Ta ) |
K |
b |
y |
b |
( y T ) |
|
|
|
|
b |
Напомним еще раз используемые в таблице 1.1 обозначения: Ka K (a), Kb K (b) – коэффициенты теплопроводности;
αа, αb – коэффициенты теплообмена на левом и правом концах стержня
соответственно;
Tа, Tb – температуры, которые поддерживаются на концах стержня при x a и
при x b ;
qа,qb – интенсивности тепловых потоков при x a и при x b .
Очевидно, что все приведенные в таблице 1.1 варианты краевых условий можно записать в виде
a |
|
|
|
(a) a |
|
|
|
y(a) a y |
|
(1.8) |
|||
|
0 |
1 |
|
|
2 |
|
b0 y(b) b1 y (b) b2 , |
|
11
при соответствующем выборе значений коэффициентов ai , bi . Например, для первого варианта условий из таблицы 1.1 имеем
a0 1, |
a1 0, |
|
a2 Ta , |
|
b0 1, |
b2 0, |
|
b2 Tb ; |
|
а для девятого – |
|
|
|
|
a0 a , |
a1 Ka , |
a2 aTa , |
||
b0 a , |
b1 Kb , |
b2 bTb ; |
Таким образом, математическая задача одномерной стационарной теплопроводности формулируется следующим образом: требуется найти функцию y(x) , удовлетворяющую на отрезке a,b обыкновенному линейному
дифференциальному уравнению (1.7), а на концах отрезка – граничным условиям (1.8).
1.3. Постановка начально-краевой задачи нестационарной одномерной теплопроводности
В разделе 1.2 рассмотрена краевая задача для одномерного стационарного уравнения теплопроводности (1.7), которая представляет собой краевую задачу для обыкновенного дифференциального уравнения второго порядка. В случае нестационарной теплопроводности к краевым (граничным) условиям (1.8) добавляется начальное условие в некоторый начальный момент времени t t0
(обычно t 0)
u(x,t0 ) (x) , |
(1.9) |
и говорят, что задана начально-краевая задача для уравнения параболического типа (1.4).
1.4. Постановка краевых задач двухмерной стационарной теплопроводности
Согласно (1.6) стационарное (установившееся во времени) распределение теплового поля в пластине описывается уравнением
|
u |
|
|
u |
|
|
||
|
K |
|
|
|
K |
|
F(x, y) . |
(1.10) |
|
|
|||||||
x |
x |
|
|
|
|
|
|
|
|
y |
y |
|
|
При решении краевых задач для уравнения эллиптического типа (1.10) наиболее часто используются три типа краевых условий.
а) Краевая задача с граничными условиями первого рода (первая краевая задача).
Требуется найти решение уравнения (1.10) в некоторой области D, принимающее на границе этой области заданные значения. Т. е. нужно найти
12
(применительно к рассматриваемой задаче) стационарное распределение температуры внутри области, если задана температура на границе этой области
u |
|
|
g(x, y) . |
(1.11) |
|
||||
|
|
D |
|
|
Здесь D – граница области D, g(x, y) – известная функция.
б) Краевая задача с граничными условиями второго рода (вторая краевая задача).
Требуется найти решение уравнения (1.10) в некоторой области, на
границе которой задана внешняя нормальная производная |
u |
(т. е. на границе |
||||||||
задана интенсивность теплового потока) |
|
|
|
n |
|
|||||
|
|
|
|
|
||||||
|
|
K |
u |
|
q(x, y) или |
u |
|
q(x, y) 0, |
|
(1.12) |
|
|
|
|
|
||||||
|
|
|
n |
|
|
n |
|
|
|
|
|
q |
|
|
|
D |
|
|
D |
|
|
где q |
. Здесь |
– граница области D, |
|
q(x, y) – интенсивность теплового |
||||||
|
|
|||||||||
|
K |
D |
|
|
|
|
|
|
|
|
|
|
|
|
q 0, то тепловой поток направлен наружу, а если |
||||||
потока. |
При этом, если |
|
q 0 , то тепловой поток направлен внутрь области. При q 0 |
имеем условие |
|||||||||||||||||||
теплоизоляции u |
0 . |
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
в) Краевая задача с граничными условиями третьего рода (третья краевая |
||||||||||||||||||||
задача). |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
Требуется найти решение уравнения (1.9) в некоторой области, которое |
||||||||||||||||||||
удовлетворяет на границе условию |
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
K u |
|
|
u |
|
T или |
u |
|
|
|
|
u |
|
|
T , |
(1.13) |
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
n |
|
|
|
n |
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
D |
|
|
|
|
|
|
D |
|
|
|||
|
|
|
|
|
|
|
D |
|
|
|
|
|
D |
|
|
|
|
|
|
|
где |
|
|
. Здесь |
|
|
|
– граница области D, |
на которой задан теплообмен с |
||||||||||||
|
|
|
||||||||||||||||||
|
|
|
||||||||||||||||||
|
|
|
K |
|
|
D |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
T ; – |
|
||
окружающей средой, |
температура которой |
равна |
коэффициент |
|||||||||||||||||
теплообмена. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Если на различных частях границы D заданы условия различного рода, то такие условия и соответствующие им задачи называют смешанными.
1.5. Вывод уравнений поперечных колебаний струны
Рассмотрим тонкую гибкую упругую нить (струну), которая в положении равновесия занимает отрезок a,b оси Ox и концы которой закреплены.
Полагая струну тонкой, пренебрегаем весом струны по сравнению с внутренними силами натяжения и внешней нагрузкой. Полагая струну гибкой, считаем, что внутренние усилия, возникающие в струне, направлены по касательной к мгновенному профилю в каждой точке, т. е. струна не
13
сопротивляется изгибу. Предполагаем также, что внешние силы лежат в вертикальной плоскости, в которой совершают колебания точки струны.
Рассмотрим элемент струны между точками x и x dx (рис. 1.1) и обозначим смещение точек струны через u(x,t) , а длину элемента струны через
ds .
Тогда
|
x dx |
|
u 2 |
|
1 |
|
|
|
2 |
||||
ds |
|
1 |
|
|
|
dx , |
|
x |
|
|
x |
|
|
откуда, предполагая смещение струны u(x,t) малыми настолько, что
|
u |
2 |
|
|
x |
1, |
(1.14) |
получаем ds dx , т. е. в пределах принятой точности удлинения участков струны в процессе колебаний не происходит. Следовательно, согласно закону Гука величина натяжения в каждой точке струны не меняется со временем и является функцией только x , т. е. T T (x) .
Рис. 1.1. Иллюстрация к выводу уравнения колебаний струны
Запишем условия динамического равновесия элемента струны, на который действуют в плоскости Oxu силы натяжения T1 T (x), T2 T (x dx), внешняя
распределенная по длине дуги с линейной плотностью F(x,t) поперечная сила и сила инерции, направленная вдоль оси Ou .
Проектируя силы на ось Ox , получаем |
|
T (x dx) cos( 2 ) T (x) cos( 1 ) 0. |
(1.15) |
Так как, согласно тождествам тригонометрии и геометрического смысла производной,
cos( 2 ) |
|
1 |
|
|
|
1 |
|
|
, |
||
|
tg2 ( 2 ) |
|
u(x |
dx,t) 2 |
|||||||
1 |
|
|
|
|
|||||||
|
|
|
|
|
|
|
1 |
x |
|
|
|
|
|
|
|
|
|
|
|
|
(1.16) |
||
|
|
|
1 |
|
|
|
|
1 |
|
|
|
cos( 1) |
|
|
|
|
|
|
|
, |
|
||
|
tg2 ( 1) |
|
|
|
|
||||||
1 |
|
|
u(x,t) 2 |
|
|
||||||
|
|
|
|
|
|
|
1 |
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
14
то, учитывая условие (1.15), из (1.16) получим T (x dx) T (x) . Откуда, в силу
произвольности выбора точек x |
и x dx , следует, что величина натяжения не |
||||||||||
зависит и от x , т. е. является постоянной, |
T (x) T0 |
const . |
|
|
|||||||
Проектируя теперь все силы на ось Ou , получаем |
F(z,t)dz, |
(1.17) |
|||||||||
|
(z) |
2 |
u(z,t) dz T0 sin( 2 ) T0 sin( 1 ) |
||||||||
x dx |
|
|
|
|
|
|
x dx |
|
|
||
x |
|
|
t 2 |
|
|
|
x |
|
|
|
|
где (x) – линейная плотность струны. |
|
|
|
|
|
|
|||||
Аналогично формулам (1.16) устанавливаем |
|
|
|
|
|||||||
|
|
|
|
|
tg( 2 ) |
|
u(x dx,t) |
|
|
||
|
sin( 2 ) |
|
|
|
x |
|
, |
|
|||
|
|
tg2 ( 2 ) |
u(x dx,t) 2 |
|
|||||||
|
|
|
1 |
|
|
|
|||||
|
|
|
|
|
|
|
1 |
x |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
tg( 1 ) |
|
|
|
u(x,t) |
|
||
sin( 1 ) |
|
|
|
|
|
|
x |
|
, |
|
|
tg2 ( 1 ) |
|
|
u(x,t) 2 |
||||||
|
1 |
1 |
|
|
||||||
|
|
|
|
|
|
x |
|
|
||
|
|
|
|
|
|
|
|
|
|
|
откуда, согласно условию (1.14), имеем |
|
|
|
|
|
|
||||
sin( 2 ) |
u(x dx,t) |
, |
|
sin( 1 ) |
u(x,t) . |
|||||
|
|
|
x |
|
|
|
|
|
x |
|
Теперь, применяя для входящих в формулу (1.17) интегралов теорему о среднем, а для ux (x dx,t) – формулу Тейлора первого порядка с остатком в
форме Пеано, получаем
|
|
|
2 |
u( 1 |
,t) |
|
|
2 |
u(x,t) |
|
|
|
||
|
( 1 ) |
|
|
|
|
|
F( 2 ,t)dx, |
|||||||
|
|
|
t |
2 |
|
dx T0 |
|
|
x |
2 |
dx o(dx) |
|||
|
|
|
|
|
|
|
|
|
x, x dx . |
|
|
|||
где 1 |
и 2 принадлежат отрезку |
Почленно деля последнее |
равенство на dx и осуществляя предельный переход при dx 0 , получаем уравнение колебания струны следующего вида:
|
|
(x) |
2u(x,t) |
T |
2u(x,t) |
F(x,t) . |
|
(1.18) |
|||
|
|
t 2 |
|
|
x2 |
|
|||||
|
|
|
|
|
0 |
|
|
|
|
||
Если струна дополнительно по всей длине связана с вязкоупругим |
|||||||||||
основанием, то для описания ее колебания можно получить уравнение |
|||||||||||
(x) |
2u(x,t) |
T0 |
2u(x,t) |
|
|
|
|
t |
|
|
|
t 2 |
x2 |
(x,t) u(x,t) Q(x,t, )u(x, )d |
|||||||||
|
|
|
|
|
|
0 |
|
(1.19) |
|||
|
|
|
(x,t) u(x,t) |
|
|
|
|||||
|
|
|
F(x,t), |
|
|
||||||
где (x,t), (x,t) – |
|
|
|
|
t |
|
|
|
|
|
|
коэффициенты |
жесткости и демпфирования |
основания; |
Q(x,t, ) – ядро релаксации, учитывающее изменение с течением времени физико-механических свойств материала основания (т. е. его старение).
15
Заметим, что при выводе уравнения (1.19) предполагалось, что реакция основания пропорциональна его деформации (модель Винклера).
В статических задачах профиль струны u u(x) определяется, согласно
(1.12), решением уравнения |
(x) |
F(x) |
|
|
|
|
|
u |
. |
|
|
|
(1.20) |
||
T |
u T |
|
|
|
|||
|
0 |
0 |
|
|
|
|
|
1.6. Вывод уравнений продольных и крутильных колебаний стержня |
|||||||
Для вязкоупругого тела при одномерном растяжении (сжатии) связь между |
|||||||
деформацией (относительным удлинением) (x,t) |
и напряжением |
(x,t) |
|||||
представляется формулой |
t |
|
|
|
|
|
|
|
|
|
|
(x,t) |
|
|
|
|
|
|
|
(x,t) |
, |
(1.21) |
|
(x,t) E(x,t) (x,t) |
R(x,t, ) (x, )d |
t |
|||||
|
0 |
|
|
|
|
|
|
где E – модуль упругости; R – ядро релаксации, учитывающее старение |
|||||||
материала тела; – коэффициент внутреннего трения. Заметим, если |
R 0 и |
0, то получаем закон Гука для упругого тела.
Рассмотрим элемент стержня (рис. 1.2), заключенный между поперечным сечением с координатами x и x dx .
Рис. 1.2 Иллюстрация к выводу уравнения продольных колебаний стержня
В сечении « x » на элемент действует сила N (x,t) (x,t)S(x), где S(x) – площадь сечения, в сечении « x dx » – сила N (x dx,t) (x dx,t)S(x dx)) .
Предполагая, что на стержень действует внешняя нагрузка, распределенная по длине стержня с объемной плоскостью F(x,t) , аналогично выводу уравнения
(1.18) получаем уравнение продольных колебаний стержня следующего вида:
(x)S(x) |
2u(x,t) |
|
|
S(x) (x,t) S(x)F(x,t), |
(1.22) |
|
t 2 |
x |
|||||
|
|
|
|
где (x) – объемная плотность материала стержня; u(x,t) – продольное
смещение сечения стержня с координатой x в момент времени t от положения, которое занимало это сечение, когда стержень находился в ненапряженном состоянии.
Учитывая, что |
u(x dx,t) u(x,t) |
|
u(x,t) |
|
(x,t) lim |
|
, |
||
dx 0 |
dx |
|
x |
|
и подставляя (1.21) в (1.22), имеем |
|
|
|
|
16
|
|
2 |
u(x,t) |
|
|
|
|
|
u(x,t) |
t |
u |
|||
|
|
|
|
|
|
|||||||||
(x)S(x) |
|
|
|
2 |
|
|
|
|
|
|
|
|
R(x,t, ) |
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
t |
|
|
S(x)E(x,t) |
|
x |
|
||||||
|
|
|
|
|
x |
|
|
|
|
0 |
|
|||
|
|
|
|
|
|
|
|
|
2 |
u(x,t) |
|
|
|
|
|
|
|
|
|
(x,t)S(x) |
|
|
S(x)F(x,t). |
|
|||||
|
|
|
|
|
x t |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
(x, )dx
Если боковая поверхность стержня скреплена с вязкоупругим основанием (модель Винклера), то приходим к следующему уравнению:
|
|
2 |
u(x,t) |
|
|
|
|
|
|
|
u(x,t) |
t |
u(x, ) |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|||||||||||
(x)S(x) |
|
|
|
2 |
|
|
|
|
|
|
|
|
|
R(x,t, ) |
|
|
d |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
t |
|
|
|
S(x)E(x,t) |
x |
x |
|
|||||||||||
|
|
|
|
|
|
|
|
x |
|
|
|
0 |
|
|
|
|
||||
|
|
|
|
|
|
2 |
u(x,t) |
|
|
|
|
t |
|
|
|
|
||||
(x,t)S(x) |
|
|
|
|
|
|
|
(1.23) |
||||||||||||
|
x t |
|
|
(x,t) u(x,t) Q(x,t, )u(x, )d |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|||||
|
|
|
|
|
|
|
|
(x,t) |
u(x,t) S(x)F (x,t), |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
t |
|
|
|
|
|
|
|
где (x,t), (x,t) – коэффициенты жесткости и демпфирования основания; Q(x,t, ) – ядро релаксации основания. Заметим, что форма записи уравнения (1.23) не изменится, если считать S и зависящими от времени t .
Статические продольные смещения u(x) сечений стержня определяются, согласно (1.23), решением уравнения
|
|
|
|
|
|
(x)u S(x)F(x). |
(1.24) |
||
S(x)E(x)u |
Для вязкоупругого стержня, находящегося в состоянии кручения (рис. 1.3), связь между напряжением , вызванным сдвигом образующей на угол , и
этим углом может быть представлена формулой
|
|
t |
R(x,t, ) (x, )d |
|
(x,t) |
(x,t) |
, |
(1.25) |
||
(x,t) G(x,t) (x,t) |
|
|
|
|||||||
|
|
|
|
|
|
t |
|
|
||
|
|
0 |
|
|
|
|
|
|
|
|
где G – модуль сдвига; R – ядро релаксации стержня; |
|
– коэффициент |
||||||||
внутреннего трения. |
|
|
|
|
|
|
|
|
|
|
Заметим, если R 0, 0 , |
то получаем |
известный |
закон сдвига для |
|||||||
упругого тела. |
|
|
|
|
|
|
|
|
|
|
Если обозначить через u(x,t) |
|
угол поворота сечения с координатой |
x , то |
|||||||
(см. рис. 1.3) из равенства rdu dx , имеем |
|
|
|
|
|
|
|
|||
|
|
r u . |
|
|
|
|
|
|
(1.26) |
|
|
|
|
x |
|
|
|
|
|
|
|
Рис.1.3. Иллюстрация к выводу уравнения крутильных колебаний стержня
17
Крутящий момент M (x,t) , действующий в сечении S стержня, соответствующем координате x , определяется формулой
M (x,t) r ds .
S
Отсюда, используя выражения (1.25), (1.26), получаем
|
u(x,t) |
t |
u(x, ) |
|
|
|
2 |
u(x,t) |
|
|
M (x,t) J0 |
R(x,t, ) |
(x,t)J0 (x) |
|
, (1.27) |
||||||
(x)G(x,t) |
x |
x |
d |
|
x t |
|||||
|
|
0 |
|
|
|
|
где J0 r 2dS – полярный момент инерции сечения.
S
Рассмотрим элемент стержня, заключенный между поперечными сечениями с координатами x и x dx (рис. 1.3). В сечении « x » действует крутящий момент M (x,t) , в сечении « x dx » – M (x dx,t) . Предполагая, что
на стержень действует крутящий момент внешних сил, распределенный по длине стержня с линейной плотностью F(x,t) , из уравнения динамического
равновесия получаем
( 1 )J0 ( 1 ) 2u(t 21,t) M (xx,t) dx o(dx) F( 2 ,t)dx,
где – плотность стержня; 1 и 2 – принадлежат x, x dx . Откуда
аналогично уравнению (1.18) получаем уравнение крутильных колебаний стержня
|
|
|
|
(x)J0 |
(x) |
2u(x,t) |
|
M (x,t) |
F(x,t), |
|
|
|||||
|
|
|
|
|
t 2 |
|
x |
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
которое, с учетом (1.27), принимает вид |
|
|
|
|
|
|
||||||||||
|
|
2 |
u(x,t) |
|
|
|
|
|
|
|
u(x,t) |
t |
u(x, ) |
|
||
|
|
|
|
|
|
|
|
|||||||||
(x)J0 (x) |
|
|
|
2 |
|
|
|
|
|
(x)G(x,t) |
|
R(x,t, ) |
|
d |
||
|
|
|
|
|
|
|
|
|||||||||
|
|
t |
|
|
J0 |
x |
x |
|||||||||
|
|
|
|
|
x |
|
|
|
|
0 |
|
(x,t)J0 (x) 2u(x,t) F(x,t).
Если боковая поверхность стержня скреплена с вязкоупругим основанием (модель Винклера), то для описания крутильных колебаний приходим к уравнению
|
|
2 |
u(x,t) |
|
|
|
|
u(x,t) |
t |
u(x, ) |
|
|
|||
|
|
|
|
|
|
|
|||||||||
(x)J0 (x) |
|
|
2 |
|
|
|
|
|
(x)G(x,t) |
x |
R(x,t, ) |
x |
d |
|
|
|
|
|
|
|
|
||||||||||
|
|
|
x J0 |
|
|||||||||||
|
|
|
t |
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
2 |
u(x,t) |
|
|
|
t |
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||||||
(x,t)J0 (x) |
|
|
|
|
|
|
(x,t) u(x,t) |
Q(x,t, )u(x, )d |
|
(1.28) |
|||||
|
x t |
||||||||||||||
|
|
|
|
|
|
|
|
0 |
|
|
|
(x,t) u(x,t) F(x,t),
t
где , ,Q – коэффициенты жесткости, демпфирования и ядро релаксации основания.
18
Заметим, что форма записи уравнения (1.27) не изменится, если считать и J0 функциями двух переменных x и t .
Статические углы поворота u(x) сечений стержня при кручении определяются, согласно (1.28), решением уравнения
|
|
|
|
|
|
(x)u(x) F(x) . |
(1.29) |
||
J0 (x)G(x)u |
1.7.Постановка статических краевых задач для струны и стержня
Встатическом варианте профиль струны, продольные и угловые перемещения сечений стержня, согласно (1.20), (1.24) и (1.29), определяется решением уравнения
L( y) (K (x) y ) |
|
(x) y g(x), |
(1.30) |
|
|
|
где y(x) u(x); a x b ;
K (x) T0 , g(x) F(x), если рассматривается задача (1.20);
K (x) S(x)E(x), g(x) S(x)F(x), если – задача (1.24);
K (x) J0 (x)G(x), g(x) F(x), если – задача (1.29).
Перечислим основные типы граничных условий при x a для уравнений
(1.20), (1.24), (1.29) в обозначениях уравнения (1.30).
а) y(a) 0; это условие соответствует жесткому закреплению левого конца
струны и стержня.
б) K (a) y (a) qa ; это условие соответствует заданию на левом конце стержня продольной силы N(a) qa для задачи (1.24) и заданию крутящего момента M (a) qa в случае задачи (1.29). В частности, если левый конец свободен, то qa 0.
в) K(a) y (a) a y(a) ; это условие соответствует упругому закреплению левого конца стержня, когда qa a y(a) ( qa или равно N(a), или – M (a) ), гдеa – соответствующий задаче (1.24) или (1.29) коэффициент закрепления.
Аналогичные краевые условия могут быть заданы и на правом конце струны или стержня при x b . Очевидно, что все возможные варианты краевых условий для уравнения (1.30) можно получить из условий (1.8) при соответствующем выборе значений коэффициентов ai ,bi .
Таким образом, рассматриваемые статические краевые задачи для струны и стержня математически формулируется так же, как и задача стационарной теплопроводности из раздела 1.2.
1.8. Краевые задачи в теории колебаний струн и стержней
Предположим, что геометрические и прочностные характеристики упругих тел (струны, стержня) и оснований, на которые они опираются, зависят только от x , и запишем уравнения движения без учета демпфирования и старения.
19
Уравнение продольных колебаний струны |
|
|
|||||||||||||||
|
(x) |
2u(x,t) |
|
T |
|
2u(x,t) |
(x)u(x,t) F(x,t). |
(1.31) |
|||||||||
|
|
t 2 |
|
|
|
x2 |
|
||||||||||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|||||
Уравнение продольных колебаний стержня |
|
|
|||||||||||||||
|
2u(x,t) |
|
|
|
|
|
|
u(x,t) |
|
|
|||||||
(x)S(x) |
|
t |
2 |
|
|
|
|
S(x)E(x) |
|
x |
|
(x)u(x,t) S(x)F(x,t). |
(1.32) |
||||
|
|
|
|
||||||||||||||
|
|
|
|
|
x |
|
|
|
|
|
|
||||||
Уравнение крутильных колебаний стержня |
|
|
|||||||||||||||
|
|
2u(x,t) |
|
|
|
|
|
|
|
u(x,t) |
(1.33) |
||||||
(x)J 0 (x) |
|
t |
2 |
|
|
|
|
J 0 |
(x)G(x) |
|
x |
(x)u(x,t) F (x,t). |
|||||
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
Уравнения (1.31)–(1.33) являются уравнениями гиперболического типа. Рассмотрим гармонические колебания упругих тел. В этом случае решение
уравнений (1.31)–(1.33) и приложенную внешнюю нагрузку F(x,t) представим в виде:
u(x,t) u* (x)sin( t ), |
F(x,t) F* (x)sin( t ), |
|
(1.34) |
где (частота колебаний) и – постоянные. Тогда для u* (x) y(x) |
получим |
||
уравнение (1.30), в котором F(x) следует заменить на F* (x) , а (x) |
– на * (x) , |
||
где |
|
|
|
* (x) (x) (x) 2 соответствует уравнению (1.31),
* (x) (x) (x)S(x) 2 – уравнению (1.32),
* (x) (x) (x)J0 (x) 2 – уравнению (1.33).
Приведем основные типы граничных условий при x a .
а) u(x,t) a (t); это условие соответствует движению левого конца струны или стержня по закону a (t) .
б) K (a) u(a,t) qa (t); это условие соответствует заданию на левом конце
x
стержня продольной силы N (a,t) qa (t) для задачи (1.32) и заданию крутящего момента M (a,t) qa (t) в случае задачи (1.33). В частности, если левый конец свободен, то qa 0 .
в) K (a) |
u(a,t) |
a u(a,t) a (t) ; это условие соответствует упругому |
|
x |
|
закреплению левого сечения стержня, движущегося (вращающегося) по закону
a (t) .
Предполагая функции a (t), a (t),qa (t) периодическими во времени,
аналогично (1.34) положим |
|
|
a (t) a0 sin( t ), |
a (t) a0 sin( t ), |
qa (t) qa0 sin( t ), |
где a0 , a0 ,qa0 – постоянные. Тогда для u* (x) y(x) будем иметь граничные условия следующего вида:
а) y(a) a0 ;
20