книги из ГПНТБ / Зак М.А. Неклассические проблемы механики сплошных сред
.pdfимеет разрыв вторых производных от радиус-вектора по коорди натам, т. е. разрыв формы, то через некоторый интервал време ни может возникнуть неоднозначность характеристической ско рости распространения этой поверхности вследствие появления разрыва нормали к ней. Если же исходная поверхность разрыва в начальный момент дважды непрерывно дифференцируема, то неоднозначности характеристической скорости ее перемещения возникнуть не может, так как дифференциальные уравнения дви жения этой поверхности линейны относительно производных dnr/dqin при п > 1, и, следовательно, «самозарождение» разры вов dnr/dqin в результате накопления разрывов dn+1 r/dqin+l не возможно.
§ 4. Границы применимости классических уравнений движения нити
1. Характеристические скорости распространения волн в нити.
Будем исходить из принципа Гамильтона в форме (II—1—8) и после введения неопределенных множителей Лагранжа полу чим
А9а |
dr* |
dr* |
|
|
р V • О V — тъ |
F-8 r* 1dqdt = 0, |
|||
Л9. |
dq |
dq |
|
|
F = |
F' + F". |
(Ill—4 -1 ) |
||
|
Здесь F', F" — соответственно мертвые и следящие силы. Представим движение нити в виде суммы возмущенного и не
возмущенного движений: г= г0+Дг, причем возмущения будем считать малыми. Варьируя только возмущенное движение и учи тывая преобразования гл. II, § 2, п. 1, вместо (1) будем иметь
A 9i
+ Д F" 8 Д г* dq| d t— 0.
Подынтегральное выражение первого внутреннего интеграла и подстановка равны нулю в силу уравнений невозмущенного дви жения. Поэтому
А9а |
|
|
д Дг |
|
д Дг* |
, |
d Дг |
— Д |
Т |
8 |
|||
dt |
~д<Г |
dq |
‘ |
|||
А9. |
|
|
|
|
|
|
+ AF"8Ar*J dqdt = 0. |
|
|
(Ill—4—2) |
|||
70
Положим, что возмущенное движение начинается со скачкооб разного изменения производных <9r/dt и dr/dq. Учитывая кинема тические соотношения на фронте разрыва (I—2—47), получим
{dAr/df]=—7,[dAr/dq], б [<?Дг*/д/]=^- Аб [d^r*/dq], где А,— скорость распространения разрыва по нити. Полагая, что вариации б [■••] обращаются в нуль на концах участка q\, qo и в моменты време ни t\ п t% получим из (2)
о |
Аг 1. 2 |
Д/Т дг |
|
|
|
|
|
dq |
|
|
|
|
|
Положим, что |
|
|
|
<7i |
|
|
|
m = o , |
|
|
(III—4—3) |
||
Тогда |
|
|
|
|||
|
|
|
q\ |
|
|
|
|
|
|
d Дг |
|
(III—4 -4 ) |
|
|
(ol2- |
Т0) |
j‘AF"rf<7 . |
|||
|
~д<Г |
|||||
Из (3) следует, что |
|
|
<и |
|
|
|
|
|
|
|
|
||
(р)2 _ Т0) <?Дг ' |
|
|
д Д г ' |
■ь°=о, |
(III—4—5) |
|
|
dq |
’ = |
0, (РХ2 — |
70) dq |
||
так как |
|
|
|
|
|
|
|
|
[ J aF " ^ ] = 0. |
|
(Ill—4 -6 ) |
||
|
|
|
<7i |
|
|
|
Для существования ненулевых разрывов [дДг/дд]-п° и [5Дг/^]-Ь°, где п°, Ь° — нормаль и бинормаль невозмущенной нити, необхо димо, чтобы рА,2=7'°, т. е.
Ч * = ± V t • (III 4 7)
Таким образом, Ai, 2 в (7) представляет собой скорость распрост
ранения поперечной волны в нити. |
|
|
(3) имеет место ра |
|
Положим, что при сохранении условия |
||||
венство [дДг/д^] = [т]= 0. Тогда нормаль п |
может претерпевать |
|||
разрыв, направленный перпендикулярно касательной |
||||
’0'[п] = 4 г |
дг Дг |
•т° = 0 |
||
dq2 |
||||
при условии, что кривизна формы нити |
Йз |
непрерывна, т. е. |
||
[й3]= 0. Таким образом, рассматриваемый разрыв состоит в том, что соприкасающаяся плоскость поворачивается на некоторый конечный угол ф относительно касательной т.
Используя уравнение движения нити в виде
<Э3г |
д / ~ дг |
и учитывая, что
F" = - P '. j |
F"2 |
J »3г_ _F% _ |
дЧ |
U3 |
dq3 ' й3 |
dq ^ dg3 ’ |
получаем
(p),2 - T°) |
г аз 4г |
. ь * — |
F"> |
‘ д3Дг ' |
•b* |
(III—4—8) |
|
dq3 |
Ь ~ |
Q, |
dq- |
|
|
Здесь F2" — проекция следящей силы F" на нормаль n, b* — пер пендикуляр к биссектрисе угла между нормалями до и после разрыва.
Из (8) получается выражение для скорости
Ч * = ± ] / r-L (r> + - ^ - ) |
(III—4—9) |
распространения волны «закручивания» нити, переносящей раз рывы угла поворота соприкасающейся плоскости.
.Положим, наконец, что [Г]=^0. Тогда из (4) с учетом (5) по лучим pA,2[ft] = [7т], где dr/dq = x-f. Следовательно,
(о)2/ ' — Т') т' = (Pl2f — Т") 1", (III—4—10)
причем [/]=/"— |Т] = Г"—Т', |
[т]= т"—г'. Пусть т'= т"= т, т. е. |
|
[т]= 0. Тогда |
|
|
= |
• |
(Ill—4—11) |
Очевидно, что Яб,б является скоростью распространения продоль ной волны растяжения. Ограничиваясь случаем малых возмуще ний растяжений [f] и предполагая материал нити идеально упру гим, получаем
V . = ± / - f - 4 г ■ |
(Ш -4 -12) |
Если в (10) положить [т]=^=0, то рЯ2/'—T' = pXzf"—Т"—0, и снова приходим к формуле (12).
Итак, в упруго растяжимой нити возможны лишь три системы волн и соответственно три характеристические скорости их рас пространения; эти скорости определяются формулами (7), (9) и (12). Заметим, что при выводе скоростей (7) и (9) не делалось предположений о малости переносимых разрывов, а при выводе скорости (12) малость переносимых разрывов была существенна.
Если нить движется по поверхности, то формулы (7) и (12) для скоростей сохраняются. Скорость (9) сохраняет свой вид лишь тогда, когда нить располагается по геодезической кривой поверхности. В противном случае
= ± У т ( Т° + -%-) • |
(Ш -4 -13) |
Формула (13) отличается от (9) тем, что вместо полной кривиз ны нити й 3 здесь фигурирует геодезическая кривизна Qr, a Р2*
72
является проекцией F" не на главную нормаль, как в (9), а на нормаль к нити, лежащую в касательной к поверхности плоско сти и образующую острый угол с главной нормалью нити.
Отметим особенности распространения волн разрывов в гиб ком шланге. Сохраняя обозначения гл. II, § 3, п. 5 и учитывая, что
и' = и' |
дг |
|
д Дг ' |
а дг* |
+ и'о |
д Дг* 1\ |
|
|
+ и'о dq |
dt |
dq J/ |
||
|
д Д г ! |
' д Дг* ' |
|
д Дг* |
|
|
+ Р dt |
dt |
|
dq |
|
|
|
в предположении несжимаемости жидкости [ц'] = 0, получим
|рХ2 + р' (X— ц'о)3|
откуда следуют три системы волн разрывов и соответствующие им три характеристические скорости
|
^1.2 |
'U'o , |
Г |
Т*—Р* |
РР До |
(III—4—14) |
|
4 - d' — |
р + р ' |
(Р + Р')2 ’ |
|||
|
|
Р + Р |
|
|
|
|
*3,( ■ |
|
т*—р *+ F'у а 3 |
(й г к ” |
' (ш - 4- 15) |
||
Р“> ± |
] / |
р + |
р' |
|||
|
р + |
|
|
|
||
|
|
к5,6- |
±Vd iP)-dT/df |
(III—4—16) |
||
для волн поперечного разрыва, волн «закручивания» и волн рас
тяжения стенок шланга. |
|
волн сжа |
Если жидкость сжимаема, то появляется система |
||
тия, распространяющаяся в случае малых Ди' со скоростью |
||
) Л .8 - « о '+ / - f - |
> |
(HI—4—17) |
характерной для соответствующего газа |
(см. § 6). |
|
2. Критерии структурной неустойчивости движения нити. Сопо ставим результаты предыдущего пункта с результатами § 2. Ес ли искать малые возмущения Дг в виде, аналогичном (2—3):
Дг = Дг0 exp Х0г
то придем к характеристическим соотношениям (5) относительновеличин [dtkTjdq\-n и [дДг/д<?]-Ь, т. е. относительно величин «из ломов» нити. Напомним, что при выводе этих соотношений не делалось предположений относительно малости изломов нити. Поэтому результаты § 2, п. 1 полностью могут быть отнесенцг к данному случаю, и, следовательно, при условии К О , вытекаю щему из (7), имеет место оценка
д Дг
dq Пу J
д 04
Дг
-О
, со
оо при t ->■О,
73
причем
д Дг |
|
д Дг |
d q |
110 |
У |
d q |
||
|
|
1 = 0 |
Таким образом, сжатая нить становится нигде не дифферен цируемой, хотя остается непрерывной. Ее форма принимает пи лообразный профиль с «зубцами» бесконечно большой частоты и бесконечно малой высоты. При этом углы при вершине зубцов могут быть достаточно малыми вплоть до нуля.
Интересно отметить, что в классе таких функций нить может «свернуться» в точечный комок без нарушения условия нерастяжимости. Однако для практических целей удобнее перейти к ре шениям, являющимся в некотором смысле обобщенными и со держащимся в классе достаточно гладких функций. С этой целью заключим мысленно нить в воображаемую бесконечно тон кую трубочку, которая будет полностью «вмещать» нить, оста ваясь непрерывно дифференцируемой кривой, ввиду бесконечно малой высоты «зубцов» нити; другими словами, нить и трубочка будут совпадать « в среднем» на любом сколь угодно малом, но конечном участке. При этом следует, конечно, отказаться от условия |dr/ds| =1, если радиус-вектором г фиксировать поло жения точек трубочки, так как расстояния между точками нити, отсчитываемые вдоль трубочки, будут уменьшаться, и, значит, | dr/ds | < 1. Поскольку это неравенство означает отсутствие гео метрической связи, реакцией которой является натяжение, то на тяжение Т < 0, имеющее место в модели нити, не войдет в урав
нения движения модели трубочки. |
Следовательно, на участке |
|
ГСО вместо обычных уравнений движения |
нити следует рас |
|
сматривать уравнение |
|
|
РОО ^ |
1== F . |
(III—4—18) |
Если F = F (s, t), то из этого уравнения сразу получается закон движения точек трубочки
|
t |
t |
r(s, f) = r0(s) + v0(s)-t + |
j |
| [ f (s, t)dt]^dt, s, < s < s2. |
|
о |
b |
Здесь r0 (s), v (s) — координаты и скорость точек нити в тот мо мент, когда на участке s i< s < s 2 возникает отрицательное натя жение.
Уравнение (18) справедливо до тех пор, пока |d r/d s|< l. На
тех участках, |
где | dr/ds | > 1 |
(неравенство имеет место в случае |
растяжимой |
нити), следует |
переходить к обычным уравнениям |
движения нити. Однако для отыскания истинного (а не обоб щенного) решения задачи о движении нити в области сжатия следует воспользоваться непосредственной минимизацией функ ционала, например, в форме (II—1—4) при условии непрерывно-
74
сти (но не обязательно дифференцируемости) функции г (q,t). Подчеркнем, что исследование движения нити в области сжатия представляет определенный практический интерес, так как сжа тое состояние нити не противоречит ни одной из аксиом механи ки и часто реализуется при работе тех или иных конструций. Так, например, если одну из точек закрепления весомой нити, подве шенной в двух точках и имеющей форму цепной линии, начать резко перемещать вниз (удар по опоре), то возникнет область, для которой Г<0. Точно те же явления, что и со сжатой нитью, происходят с гибким шлангом при условии, вытекающем из (14):
Т* < Р |
РР |
|
Р + р' |
|
|
|
|
|
или |
|
(III—4—19) |
9Т < 9'р + 99'щ '\ |
||
Как видно, это условие более слабое, чем для нити, и выполня ется даже при растяжении шланга, если достаточно велики дав ление жидкости в нем и скорость протекания этой жидкости. Для исследования гибкого шланга в области (19) необходимо перей ти к непосредственной минимизации функционала (II—1—4).
Итак, границами применения дифференциальных уравнений движения нити (и соответствующих им интегральных уравнений в форме (1)) являются неравенства Т< 0 или неравенства (19), которые, кстати говоря, не зависят от свойств материала нити.
3. Критерии неустойчивости формы нити. В рассуждени предыдущего пункта была использована лишь характеристиче ская скорость распространения поперечных волн. Обратимся те перь к скорости волны «закручивания» (9) или (15). При выво де этих скоростей также нигде не оговаривалась малость разры вов угла закручивания ср, так что относительно этого угла прохо дят рассуждения предыдущего пункта, и, следовательно, при вы полнении условий
Т < |
т* < р* — |
РР'«о'3 |
(III—4—20) |
|
|
р + р' |
|
соответственно для нити и гибкого |
шланга возникают всюду |
||
плотные разрывы угла закручивания ср. Таким образом, при вы полнении условий (20) становится недетерминированной форма нити или шланга (в частности, кручение формы, а следователь но, и положение осей естественного трехгранника); поэтому ста новится недетерминированной и сама сила F", если ее задать про екциями на оси естественного трехгранника. Для отыскания истинного решения следует обратиться к прямой минимизации функционала (II—1—7).
Итак, дополнительными границами применимости дифферен циальных уравнений движения нити (и соответствующих им ин
75
тегральных уравнений) являются условия (20), которые не зави сят от свойств материала нити.
Характеристические скорости распространения продольных волн (12) и .(16) не являются специфическими для нити, и их изучение будет проведено в последующих параграфах при изу чении упругого твердого тела и идеального газа.
4. Случай нулевой характеристической ркорости. При иссл довании характеристических скоростей нити (7) и (9) рассмат ривались случаи, когда выполняется точное неравенство Г<0. То же относилось и к гибкому шлангу (неравенства (19)). Можно показать [11], что при выполнении неравенств противоположного знака соответствующие дифференциальные уравнения малых возмущений корректны в классе дифференцируемых функций, п возникающие в силу начальных условий локальные разрывы со ответствующих параметров остаются локалвными и в дальней шем, имея при этом вещественную (характеристическую) ско рость распространения. Однако остается неясным предельный случай, когда характеристическая скорость оказывается равной нулю. С целью исследования такого случая рассмотрим один интересный пример.
Рассмотрим нерастяжимую гибкую нить, имеющую в невоз мущенном состоянии малую кривизну и малое кручение формы. Исследуем вопрос о единственности н устойчивости решения за дачи о малых возмущениях для случая, когда один конец нити свободен [14].
Будем исходить из векторного уравнения динамики нити pd2rfdt2 = d/ds (Т dr/ds)-{-F, |dr/<3s| = l. Здесь р — плотность, г — радиус-вектор точек нити, Т — натяжение, F — внешние силы, s — дуговая координата, t — время.
Для малых возмущений Дг и АТ следует (dr0/ds) • (dAr/ds) =0,
рд2Дг dt2 — д ds (Т0д Ar ds) д ds (Д Т drn ds). (а)
Здесь T0(s), г0(s )— невозмущенные значения Т и г. Разложим возмущение Дг по натуральным осям т0=т(1), п0=т<2>, &0= tt3>, со
ответствующим форме нити в невозмущенном состоянии: Дг= = Дгг--т«.
Положим вначале, что форма нити в невозмущенном состоя нии близка к плоской, т. е. величины dbQ/ds, d2b2/ds2, d20r/ds2, dn°/ds имеют тот же порядок малости, что и возмущения Дг, АТ. Тогда, умножая (а) скалярно на Ьо, получим
(б)
Таким образом, бинормальные малые возмущения могут изу чаться независимо от остальных.
Если форма нити в начальном состоянии имеет еще и малую кривизну, т. е. величины dxQ/ds имеют тотже порядок малости,
76
что и возмущения А г, |
А Т, |
то, умножая |
(а) |
скалярно на п 0 и то, |
||||
получим |
|
|
|
|
|
|
|
|
<Э3 Дг3 |
д { |
у, |
д Дг2 \ |
d 2 Д/'i |
|
d Д Т |
|
, . |
= |
|
|
Р - д ё - = -д Г -- |
|
М |
|||
Следовательно, в последнем случае все |
компоненты |
малых |
||||||
возмущений могут изучаться независимо один от другого. |
и |
|||||||
Рассмотрим вначале уравнение (б), полагая, что |
р (s) |
|||||||
To(s) — известные функции, |
характеризующие |
невозмущенное |
||||||
состояние нити. Если нить |
|
закреплена |
в |
двух |
точках |
(s=0, |
||
s = l) или движение ее в этих точках задано кинематически, |
то |
|||||||
существует единственное решение, удовлетворяющее начальным и граничным условиям [12]:
Ar3(s, 0) ='? (s), д Д/'з dt |/=0= J» (s), Ar3(0, f) = |
\i.l (t), |
Дг3(/, = |
(г) |
Пусть в точке s — l нить имеет свободный конец |
|
Го|8=г=0. |
(д) |
Тогда второе граничное условие в (г) отпадает. Покажем, что и в этом случае единственность решения обеспечивается, правда, при некоторых ограничениях.
Следуя [12], допустим, что существуют два решения рассмат
риваемой задачи Ar3' (s, t) и Ar3" |
(s, t) и рассмотрим их разность |
||||||||
w (s, t). |
|
t) |
удовлетворяет однородному уравнению с од |
||||||
Функция w (s, |
|||||||||
нородным дополнительным условием |
|
|
|
|
|||||
= |
|
|
^ (s, 0) = |
0, |
дщ'д£\/=0= |
0 , |
w(0, t) = 0 . |
||
Рассмотрим энергию |
i |
|
|
|
|
|
|
||
|
|
— |
(dw dsf + р(dw dt)2} ds. |
|
|
||||
|
|
0 |
|
|
|||||
Очевидно, |
что |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
i |
|
|
|
|
|
|
|
|
dE (t)/dt= j |
{T0dwds (д'-w ds dt) + pdwldt(d2w dt2)j ds. |
||||||||
При этом |
о |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
d2w |
ds |
T |
dw |
|
d_ |
T |
dw |
ds. |
|
|
||||||||
о |
dsdt |
0 |
ds |
|
ds |
0 |
ds |
||
|
|
|
|
|
|
|
|
|
|
Подстановка обращается в |
нуль |
в |
силу граничного |
условия |
|||||
w (0, t) =0 |
и условия |
(д). Следовательно, |
|
|
|
||||
77
I
dE d t — ^ dw,c^[ p d'lw:df- — djds (TQdwtds)l ds = 0.
0 |
|
|
Учитывая начальные условия, получаем |
|
|
|
i |
|
£(£) = const = .£(0) = |
Jjro^TO'^s)2 + |
ds = 0. (e) |
В отличие от [12] из (е) здесь еще не следует, что w (s, t) = 0, так как имеют место более слабые неравенства, чем в [12], а именно
р (s )> 0 при |
Т0 ( s ) > 0 при 0 ^ s < /, To\s=i=0. |
И только если искать решение в классе функций, в котором про изводная dw/ds непрерывна в замкнутом промежутке
из (е) вытекает единственность решения.
Такое своеобразие поставленной задачи связано с тем, что исходное уравнение является гиперболическим в открытом про межутке 0 ^ 5 < /, но вырождается в параболическое в закрытом промежутке
Исследуем устойчивость решения в |
смысле |
корректности |
||||
в постановке задачи с начальными условиями. Пусть |
||||||
<Р(s) > 0 при 0 < |
s°2 < s < s0! < I, |
<р (s) = |
0 |
при |
s \ |
< s < I, (ж) |
i(s) = |
0 при 0 < s < / , |
ji.1(z^) = |
0 |
при |
ty-0. |
|
Тогда |
|
|
|
|
|
|
^ > = Т 1 М ^ Г + Ч ^ ) > = Т | И ^ Г +
так как |
|
-I- Т, |
|
)"I ds = |
Е0 = |
const > 0, |
(а) |
|
|
|
|
|
|
|
|
||
d E ( t ) |
гр |
д Дг3 |
d Агз |
С ( d* |
Дг3 |
_d_ |
ds = 0. |
|
d t |
0 |
d s |
dt |
J l ^ |
ds |
|||
|
||||||||
|
|
|
|
о |
|
|
|
|
Здесь si и S2 — дуговые |
координаты переднего |
и заднего |
||||||
фронтов волны'разрыва производных d2Ar3/dt2 и d2Ar3/ds2, при чем Si = Si°, S2=S2° при ^=0.
•Очевидно, что Дг3= 0 при s<si и s> S 2-
Из дифференциальных уравнений характеристик
dSi/dt= + Y Т0/р, ds2jdt= ± Y Т01р
находим уравнения характеристик, проходящих через si0 и S20:
78
<i-n |
di| |
при 0 < s,, s2 < |
l- |
|
Ут0Ш(г>) ’ |
лI Y 7 oCnj/pCO |
|||
|
|
|||
|
|
|
При s\,i = l в силу (д) имеет место особое решение, совпада ющее для обеих характеристик. Могут представиться два слу чая.
1°. Несобственный интеграл
|
|
|
ti |
|
rfT) |
|
|
|
|
|
IV |
|
|
(и) |
|||
|
|
T0(-q)/p (.т)) |
|
|||||
|
|
|
|
|||||
сходится при т] —>■I, т. |
е. совпадение характеристик происходит |
|||||||
при конечном t='t*. Тогда |
|
($2—Si)—>-0 при |
Но в этом слу |
|||||
чае, как следует из (з), |
|
|
|
|
|
|
|
|
дАг3 |
|
|
д Аг3 |
—> оо |
при S -*■ I |
|||
Р dt |
+ Т0 |
Os |
|
|||||
и ввиду ограниченности р и Т0 с учетом |
(д) |
приходим к оценке |
||||||
|
д&г3 |
-> |
со |
|
при S |
А |
(к)- |
|
|
dt |
|
|
|||||
которая имеет место при сколь угодно малых начальных усло виях.
Эта оценка свидетельствует о некорректности постановки за
дачи в замкнутом интервале |
т. е. о неустойчивости ре |
шения вблизи свободного конца. |
|
Сформулированный результат получен при частных началь ных условиях (ж). Однако, пользуясь принципом суперпозиции для исходного линейного уравнения и добавляя к произвольным начальным условиям условия (ж), которые могут сколь угодно мало отличаться от нулевых, придем к той же оценке.
Рассмотренная задача получилась в результате линеариза ции исходной физической задачи, поэтому математическую фор мулировку неустойчивости следует ослабить и записывать в ви де
■д^[3 ■> 8 > 0 при ®(s) -» 0.
Физическим проявлением отмеченной выше неустойчивости будет резкое повышение скоростей точек нити вблизи свободно го конца (щелчок кнута). Сходимость интеграла (и) имеет ме сто в большинстве практически важных случаев, в частности, в случае, когда однородная весомая нить имеет прямолинейную
невозмущенную форму. |
расходится при |
2°. Пусть несобственный интеграл (и) |
|
т. е. совпадение характеристик происходит |
при £*->оо. Тогда |
любое возмущение, возникшее в промежутке О<As<A в конечном интервале времени не достигает свободного конца, и любое воз
79
