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

книги из ГПНТБ / Салтанов, Н. В. Гибкие нити в потоках [монография]

.pdf
Скачиваний:
24
Добавлен:
19.10.2023
Размер:
5.79 Mб
Скачать
(П Д 95)

 

a s

m.’ , 2 >

/ “ *

’ 6 , 2 m ^

' ■ f f c r ' t - l f c F l - r - y f -

% * ( t . 7 ) , \ - ( X , , X 2 ) .

В результате получим систему уравнений

 

Щ

+ х з

 

 

 

о,

(П .196)

 

 

 

 

 

 

 

d

 

 

 

 

 

(П .197)

 

■ f ( f j O + X t f t . r f . f i - 0 .

 

 

 

V

 

*

 

+ ? * ,

 

(П .198)

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d p j

 

 

 

 

Ш .199)

 

 

d X = <“

 

 

 

 

Можно видеть, что соотношения

(П .198) и

(П.199) тождественно удо­

влетворяют уравнению (I I .3 ) .Разрешая

систему уравнений

(П .196)

-

(П .199), можно найти величины

i 7” ( £

) ,

( f )

и J ? ( £ ) .

Далее

с

помощью первого

соотношения (П .195)

можно определить величину Т ( { ) .

^ А с и м п т о т и ч е с к и й

 

м е т о д .

Пусть величины

Х3

и Х^ представим

в виде

 

 

 

 

 

 

 

h

= хзо ( ^ 6 < Р гР )+ Х з 1

 

)

Ш.200)

 

 

 

 

 

 

 

где влияние величин Хз ; is XJ?

на характеристики конфигурации нити

в определенном смысле мало по сравнению

с влиянием величин Я ^ и Х1 д .

В этом случае решение системы уравнений

(П .196)-(П .199)

разыскива­

ем в виде разложений

 

 

 

 

 

 

 

е-б-в+*г+е2 + - >

 

 

 

 

 

(П .201)

 

 

РгЪ оРРл+Ъ г*-',

 

цце индекс у величин справа указывает их порядок по выбранному ма­ лому параметру. Подставим выражения (П .201) в уравнения (П .196) - (П .199). Проводя в получившихся соотношениях разложения и приравни­ вая затем величины одного порядка малости, для решения в нулевом приближении и поправок к нему первого и второго порядков получаем системы уравнений

+хзо(%>&о<Pjo ’Ро) = °>

(П .202)

 

(*>’ 60 ' Рло ’^ о ) ’ 0’

(П .203)

 

70

Ибд

 

Y i + Щ ’

= -+

 

0 £ l2 =

ju n

d $

r

0

( S

. e ^ h o ' K "

(П .204)

(П .205)

*щ)*»(*•*•^ ’&*• ш*206)

d i f t h * 9! ? » ) * *Л ( Z - C o ’ P n , ’ ? ' , ) *

 

+^'w)*10

( * ’ бо Ж ? о ^ ° >

 

 

 

 

 

(П .207)

 

& .

*

Ро ■?,

(П .208)

 

d s

 

 

У Ж

 

 

 

 

 

 

 

dp 1

=

-

(П .209)

 

И

М

 

a

d r

 

 

1 7

 

J

 

а

Ж ) *37 ( * &о Л о Л )+

 

 

^10

j. a

tyo

. j L -

 

d . ir а

z dff„

* ^

'л2 ' Щ

71

 

' П ’ к ^ к ? ' к Л х‘° ( ^ Л ) - . о - ™ »

/rkJл>’

 

 

кг/

 

 

щ *

 

 

 

Ш- 2 П >

a'fc

 

■7Г ,г

(П .212)

Z—A ~ * -------------/ 7T>

* fPo’Pi}

„■*•-«►

 

 

^ ' 2 / u ^ J ! ^ l

2Pd P 2 ~

f + p ;

 

 

 

 

а р i2

(П .213)

 

Отметим, что уравнения для последующих приближений записывают­

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

к гравитационной

составляющей, для величины Х0

имеет место соот­

ношение

__

 

 

%д ~ a ( G ) f n f i * е3 .

(П .214)

71

Напомним, что величина ) есть плавучесть нити. В случае,

когда мал параметр, характеризующий отношение интегрального значе­ ния распределенной силы к характерному значению натяжения^имеет

место

хд = °.

(П .215)

Подставляя выражение (П .214)

или (П .215) в уравнения для последую­

щих приближений (П .202) - (П .205), (П .206) -

(П .209) и т . д . , можно

убедиться в том, что решение получающихся при этом уравнений сво­

дится к последовательному вычислению квадратур.

Рассмотрим для про­

стоты

случай

(П .215).

Решение

системы уравнений

(П .202) -

(П .205)

здесь

тлеет

вид

 

 

 

 

_

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

?d =&o

'

 

 

 

 

 

 

 

 

 

 

 

 

- *°0 *

 

 

 

(

? - ? * )

, K +

f i a

-

к в) ,

 

(п ,2 1 6 )

где

ff'fi,

fig ,

ff" ,

„О „

->0

постоянные

интегрирования.

Из уравнения

г;и

и

р

“-

(П .206)

с

учетом выражения

(П .215) получаем

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(II. 217)

где

 

sro

-

 

 

 

 

 

 

t

 

величины

_

,

/>

и

опреде­

J j

постоянная интегрирования,

о,

ляются соотношениями

(П .216).

 

 

 

 

 

 

 

 

 

 

Из уравнения

(П .207) с

учетом выражения (П .215)

получаем

 

 

 

 

 

 

А

 

 

 

О

 

J0 •>„

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

( u f -

постоянная интегрирования,

величины

&0 , f i Q, &

, р

ж

.определяются согласно соотношениям Ш.2Т6) и (П.217).

 

 

 

 

Зная величину р 1

 

, с

помощью соотношений

(П .208)

и

(П .209)

находим

 

 

 

 

 

f

I h + fi

 

 

Г

 

 

 

 

 

 

 

 

 

 

 

б?. = в ° i

f

о

 

И fi

 

 

 

 

 

 

 

 

 

 

 

Л

 

 

 

 

 

 

 

 

(П .219)

 

 

 

 

 

1

1

Л

p z

 

 

 

 

 

 

 

 

где

и р

-

постоянные

интегрирования,

величины fi Q и

р 1

опре­

деляются

 

согласно соотношениям (П. 2Т6) и (П. 218).

 

 

 

 

 

 

Таким образом, показано, что в случае

(П .215) отыскание

попра­

вок первого

порядка к

решению

(П .216)

сводится к квадратурам. Да­

лее

исходя из системы уравнений (П .210) -

(П .213) можно убедиться

в том, что отыскание

поправок второго

порядка к

решению (П .216)

также

сводится к квадратурам

и т .д .

 

 

 

 

 

 

 

 

 

К р а е в ы е

 

у с л о в и я

д л я

п о с л е д у ю щ и х

п р и б л и ж е н и й .

 

Получение краевых условий для последующих

приближений в рассматриваемом

случае

проиллюстрируем на примере

72

следующих условий:

 

 

 

 

 

 

 

 

 

 

 

, Р / Х, 0 = ° ’ &/ К. 0 = ° -

 

 

(П .220)

н.

 

 

 

 

 

 

 

 

 

 

* ( Р ) / К. } ’

 

 

 

где Ф (р )~ приложенная к

концу нити сосредоточенная сила, являющая­

ся заданной функцией своих аргументов.

 

 

 

 

 

 

 

При этом в качестве

характерной длины

/

выбрана у 3

-

коорди­

ната концевой точки нити. Используя обозначения

(П .195),

а

также

обозначение

-*

 

 

 

 

 

 

 

 

 

 

ФЛ

3 (

ф, -

V

,

 

 

 

(П .221)

условия (П .220) преобразуем к

виду

 

 

 

 

 

 

 

 

 

P i/x = o

° ' в/ х - о ~ 0.

 

 

 

 

Подставим выражения

(П .201) для

,

&

и

f t

в условия

(П .222).

Проводя разложения и приравнивая величины одного порядка малости,

краевые условия для решения в нулевом приближении и поправок к не­ му первого и второго порядков малости запишем в виде

 

Ф±о1х=о 0 >

^ о /х = 7 ~

(П .223)

 

 

 

 

=- / S(Ъ о * )]Ь , - ( W b H ?(?ю •

 

 

pn h ; * o ~ 0 ' е’! $ = 0 = О.

 

 

Г

W /f.r

(fijo •*)]L , fJ

ш ' г 2 4 >

 

?/r./ '* {[$ 12 • ajsio * 2

 

7

'

& о ? г + % Ъ * % & / г ~ Г

/ (П .225)

Ч /? п ж * 1 (ъ -ж ,Т ]ъ

■ )

 

Краевые условия для поправок к решению более высокого порядка малос­ ти имеют аналогичный вид. Пусть

 

 

 

Ф = ф(к) , х= 7 ( с , е,]и)

(П .226)

Удовлетворяя

с

помощью выражений (П .216) и (П .217) -

(П .219) усло­

виям (П .223)

и

(П .224),

определяем в явном ввде входящие в указан­

ные выражения постоянные

интегрирования. В результате

решение в ну-

73

левом приближении и поправки к этому решению первого порядка малое

ти запишем в

виде

 

~

 

 

 

 

 

(П .227)

 

 

У> = / Хи « . ео ’ ? о )

 

 

 

 

 

 

(П .228)

Можно у к азать

ряд

представляющих определенный ин терес

сл у ч аев ,

к о гд а интегралы ,

входящие в соотношения

( П .2 2 8 ) , выражаются через

элементарные функции. Очевидно, одним из

таких случаев

я в л я ется

случай % = %(ju).

 

 

 

Полученные выше решения удобны, когда краевые условия заданы

при фиксированных

значениях координаты

X - Если краевые условия

заданы при фиксированных значениях какой -нибудь другой переменной,

то ц елесообразн о и сп ользовать последнюю в к ач еств е независимой.

В связи с этим отметим, что вариант рассм отренного выше асимптоти­

ческого м етода

с

текущей длиной нити

в к ач еств е

независимой

перег

менной изложен

в

р аботах / 1 1 8 ,

1 2 0 7 .

При этом

эффективность

и

точ­

н о сть м етода проиллюстрированы

на

конкретных

примерах.

 

 

Выше рассмотрены случаи наличия

в зад ач е

равновесия гибкой ни­

ти в потоке малого п арам етра,

к о гд а

отыскание

эффективного

прибли­

женного решения

 

своди тся к квадратурам . Е ст ест в ен н о ,

когда

не

уд а­

е т с я ввести удобный малый парам етр,

могут применяться др уги е

мето­

ды получения решений исходной системы уравнений

(численные,

графи­

ч еск и е , методы

аналогий и д р у г и е ).

При этом ,

в

ч астн о сти , могут

о к а за т ь ся удобными формы уравнений

равн овесия

гибких

нитей

в

пото­

к а х , полученные

 

выше.

 

 

 

 

 

 

 

 

Г л а в а Ш

НЕЛИНЕЙНАЯ МОДЕЛЬ ДИНАМИКИ НИТЕЙ В ПОТОКАХ

I . Основные уравнения

Пусть эффекты, связанные с растяжимостью нити, не существенны.

Тогда имеет

место условие нерастяжимости нити следующего вида:

 

(Ш.1)

где г и s

г- соответственно текущие радиус-вектор и длина нити.

Силовой

баланс для единицы длины нити представим так:

 

(Ш.2)

 

 

 

 

(Ш .З )

где t - время,

Т - натяжение,

m

и w - соответственно

масса и

плавучесть единицы длины нити,

Рл

- вес к -го сосредоточенного

груза, sk - s

- координата к

-го

сосредоточенного груза,

N - чис­

ло сосредоточенных грузов,

L* - характерная длина,

у - ускорение

силы тяжести, ]~е

- распределенная сила, обусловленная упругими

свойствами нити,

- распределенная сила гидродинамического воз­

действия потока на нить, /

- функция Дирака /907,

нормированная

следующим образом:

st *o

 

 

где <р - вектор упругого усилия в нити.

Поскольку продольная"составляющая вектора усилия в нити"(натя­

жения Г ) выделена,

(Ш.6)

75

(ШЛ2)

В пренебрежении инерцией вращения нити и сосредоточенных гру­

зов и распределенным моментом внешних сил уравнение моментов для

нити запишем в виде /52, 77,

1177

 

 

 

 

 

 

 

 

т

+

ds

е

=0

 

(ш .7)

 

М

- момент

ds

 

'

 

 

 

где

напряжений в

нити.

 

 

 

 

 

Пусть момент

напряжений

М

пропорционален кривизне

нити

в

рассматриваемой точке и направлен

по

бинормали к ней

 

 

 

 

 

 

 

ds

Os2

.

(Ш.8)

 

 

 

 

 

 

 

 

где

£ -

эффективный модуль Юнга нити, I

- момент инерции попереч­

ного

сечения нити относительно диаметра.

Подставляя выражение

(Ш.8)

в соотношение (Ш.7) и разрешая получившееся уравнение относительно

вектора упругого усилия

в нити

ф ,

находим

 

 

 

 

Ф = - £ Г ^ + Л ^ ,

(W-9)

 

А -

 

е

d s 3

ds

 

где

произвольный множитель.

 

 

 

 

Подставляя выражение

(Ш.9) в условие (Ш.6) и учитывая уравне­

ние

(Ш .1),

определяем величину

А

. В результате

получаем

 

 

 

9е = £1

d s [ ds

d s 3J

(ШЛО)

 

 

 

 

Отметим, что если конфигурация нити представляет собой дугу окруж­

ности,

то из выражения (ШЛО) получается естественное

следствие

= 0

/827.

Таким образом, при сделанных допущениях распределен­

ную силу

/ е

, обусловленную упругими свойствами нити,

определяют

выражения (Ш.5) и (ШЛО).

 

 

 

 

Проведем детализацию гидродинамической составляющей распреде­

ленной

силы,

действующей

на нить. Представим ее в следующем виде:

 

 

 

 

 

J/,~ f a +/ п *i f * f v '

(шли

 

j

-

 

 

/ п -

 

где

сила инерции^присоединенной массы жидкости,

сила

сопротивления формы, /

- сила сопротивления трения,

-

боко­

вая сила, обусловленная срывом вихрей.

 

 

 

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

щейся распределенной инерционной нагрузки /977. Скорость движения распределенной нагрузки отождествим с касательной к нити состав­ ляющей w относительной скорости потока и нити,

dr ds

76

где

V

- скорость потока ^жидкости. (Всюду в настоящей главе по­

ток

считается однородным,

V = c o n s t .}

 

 

 

 

В результате для силы

/

запишем

 

 

 

 

 

J 0

 

йо ) I Й-, ■•//У

} '

1.13)

 

м -

 

 

tfy //(*??

г im

 

где

присоединенная масса

единицы длины нити.

 

 

 

Выражение (ШЛЗ) для силы, обусловленной инерцией присоединен­

ной массы жидкости, представляет собой определенное обобщение фор­

мулы, рассмотренной в работах /8 0 , 8 1 ,

180, 184,

18§7. Силу /~п при­

мем пропорциональной квадрату нормальной к нити составляющей отно­

сительной

скорости

 

 

H r

 

 

/ = _ //& £ ? * ? .

(ill.14)

 

Jn

2

i s

d t / i s J J i s /\ i t J i s j

 

 

где Kn -

коэффициент сопротивления формы нити, i - ее диаметр,

р- плотность жидкости.

Силу / представляется естественным принять в следующем виде:

лК / p d j/r _ A E j^ y . И ?-

(Ш.15)

4

где К - коэффициент сопротивления трения нити.

Как представляется, аппроксимационные соотношения для боковой силы / , обусловленной срывом вихрей, в настоящее время разра­

ботаны еще

недостаточно /43] . Поэтому

конкретизацию выражения для

величины

целесообразно проводить в

каждом специальном случае

особ’о . Отметим, что предложенное разбиение силы р^ на указанные

выше составляющие носит в определенной степени условный характер.

Тем не менее, исходя из существующих в настоящее время представле­

ний, его можно считать целесообразным. В подтверждение этого отме­

тим, например, работы / 8 , 180,

182, 184, 1867. В работе /1 8 2 7 , в

частности, рассматривается

вопрос

о выборе аппроксимирующего экс­

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

выражения для силы, действующей на колеблю­

щийся в жидкости цилиндр (величин

Kn ,

М и т . п . ) .

 

Подставляя выражения (111.3),

(Ш .5),

(Ш Л О ),(Ш .II)

.1 3 )

11.15) в соотношение

(Ш .2),

приходим к уравнению

 

1 г р

 

i t 1

 

(т дГ\-

 

a - t J J _

 

i s !

i s ) ~

 

m,rip/ ^

 

 

 

 

dr d3r l )

 

> + {■ 1

r /j 3

 

is l i s t~ d s * i s 3] !

 

LX A*7

 

 

 

 

 

 

-'” *

4

4

1

(

4

r

p L

<fsi i \ i t

Л

" dsJ

J isJ

 

 

+//nJp/_ I7l J ( 7

 

 

 

Щ

г t i s / (

 

 

И / i s Л i s 11

i t j

p $

/ ( ? . * ) .

 

i r

l

(111.16)

isj*

 

77

Многоточие в выражении для силы Д . означает, что эта сила может зависеть также от производных радиуса-вектора Т по s , t и т .п .

Уравнение движения нити (Ш.16) удобно при получении условий "сшив­ ки ", которым должны удовлетворять решения в местах расположения сосредоточенных грузов, при выделении критериев подобия и в некото­ рых других случаях. В промежутках между сосредоточенными грузами оно несколько упрощается и принимает следующий вид:

dr Гдг З37

m ~—3t^ -Ts ( Ti ) =w7- + E I^ds I ds I вs * ds3jw -

 

Wi

T 3s,

 

 

 

 

 

 

 

 

dt

Щ

 

 

. «•*>

 

Ч

 

 

 

Пусть величина

д

как функция времени,

текущей длины и производ­

ных радиуса-вектора по

t

и s

задана.

Тогда получаем, что систе­

ма четырех нелинейных уравнений в частных производных

(Ш.1)

и

(Ш.17) служит для

определения четырех

величин у 7, у 2 ,

у а и

Т .

(Напомним, что

здесь

величины

у 1, у 2

и

у ч е с т ь

проекции радиуса-

вектора элемента длины нити г

на оси декартовой системы коорди­

н ат). Умножая уравнение

(Ш.17)

скалярно

на единичный вектор каса­

тельной к нити

(

Зг /

ds

)

и делая естественное

допущение о

том,

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

Ё 1 ds

(Ш.18)

Первое слагаемое справа в уравнении (Ш.18) описывает изменение на­

тяжения вдоль

нити

за счет сил инерции нити, второе - за

счет пла­

вучести

нити

и третье - за счет

гидродинамической силы сопротивле­

ния трения нити.

 

 

 

 

 

Введем в рассмотрение скорости текущего элемента нити

v

в

( v, ,

v , v3 )

и единичный вектор касательной к нити

m

г

в

(n t j,

тг , т} )

с помощью соотношений

 

 

 

 

 

 

3 7

V ,

(И1.19)

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

I E

'•

(Ш.20)

 

 

 

 

ds

 

 

 

Пусть для вектора т справедливо представление (П .П ). Тогда урав­ нение (Ш.1) удовлетворяется тождественно. Подставляя выражения

78

(Ш.19) и (L1I.20) в уравнения (Ш.17) и (Ш Л 8), соответственно запи-

~

 

 

 

 

 

 

 

я - 22)

Система уравнений,

состоящая из соотношения (П .И )

и уравнении

(Ш.19) - (1D.2I),

служит для

определения девяти

величин у 1 , y z, y 3, vJt

,1 $ , Т , ы. и

f .

Она может оказаться более

удобной по сравнению

с системой (Ш.1) и

(Ш .17),

например,

при проведении численных рас­

четов. Уравнение (Ш .22)(или

(Ш .18))

удобно при получении различных

оценок для натяжения. В некоторых случаях

гол удобнее

воспользовать­

ся вместо одной из

компонент уравнения (Ш .21).

 

 

 

Полагая в приведенных выше уравнениях

скорость

Г

= 0 , полу­

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

находящейся в покоящейся

жидкости.

 

 

 

 

 

 

 

 

 

2 .

Некоторые предельные

случаи

 

 

С л у ч а й

м а л ы х

у г л о в

м е ж д у

 

н и т ь ю

ис к о р о с т ь ю п о т о к а . Далее за исключением специ-

ально оговоренных случаев будем считать, что ось

у 7 направлена

ВДОЛЬ потока,

*

*

 

 

 

о Т

(Ш.23)

Пусть углы между нитью и скоростью потока малы,

т .е .

г я se~7 + у2 Т2 v- у3Т3 ,

(Ш.24)

где

 

 

 

 

 

« и

(Ш.25)

Пусть также выполнены

следующие неравенства:

 

 

 

 

(Ш.26)

 

м

ду3

(Ш.27)

 

a .

I s V « 7 .

 

79

Соседние файлы в папке книги из ГПНТБ