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

862

.pdf
Скачиваний:
2
Добавлен:
09.01.2024
Размер:
6.93 Mб
Скачать

91

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

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

 

 

 

 

 

 

 

 

 

d

d 0

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

Vd

 

 

nd

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

V 2

 

 

d qnd

 

 

 

 

 

 

 

 

 

 

 

dt

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

fd

 

 

.

(4.47)

 

 

 

 

 

 

nVd Ad

 

 

 

 

 

 

 

Так как в качестве объемных сил в данном случае являются силы воздействия электромагнитного поля, то согласно (4.41)

f E

1

( j H ) .

(4.48)

e c

Кроме того, под А подразумевается приток энергии только за счет электромагнитного поля

A E e V j E .

(4.49)

С учетом этого систему уравнений сплошной среды в дифференциальной форме можно получить в виде:

 

 

 

dV

 

 

dt

 

 

 

 

d

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt divV 0

 

 

 

 

 

 

 

 

 

 

ˆ

1

 

 

 

divP eE

c

( j H)

.

(4.50)

 

 

 

 

 

 

 

 

V 2

 

 

ˆ

 

 

 

 

 

 

 

 

 

 

 

div q div(P V )

E j

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

К сожалению, данная система не является еще замкнутой.

Эти уравнения позволяют определить все механические характеристики, если известно электромагнитное поле. Для этого необходимо либо задавать плотность тока как функцию координат времени, либо связывать плотность тока с остальными величинами, входящими в уравнения.

В этом случае уравнение механики, уравнение Максвелла и выражение для плотности тока должны рассматриваться совместно.

Соотношения, связывающие плотность тока с остальными величинами, называется обобщенным законом Ома [4.16], например, в такой форме

jnc E 1V H eV ,

c

где nc – величина, называемая проводимостью среды и зависиящая от свойств среды;

92

eV – конвективные токи, возникающие в частично ионизированном газе. После ряда некоторых упрощений в этом случае можно получить си-

стему уравнений магнитной гидродинамики

 

 

 

 

 

 

 

d

divV 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dV

 

 

 

 

 

 

1

 

 

 

 

 

 

 

qrad

Р

rot (H H) div ˆ

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

4

 

 

 

 

 

 

 

 

dS

 

 

 

 

 

V

 

 

 

 

 

T

divq

(rotH )2

 

 

 

 

 

 

m

 

 

,

(4.51)

 

 

dt

4

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dH

rot (V H ) rot (Vm rotH )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

divH 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где div ˆ

ik e

;

 

 

 

 

( , P, H,

Ve

)

 

iR

iR

 

 

 

xk

i

 

 

 

 

 

 

 

 

xm

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q q( , P, H,

T

) ;

 

 

 

Ф

1

 

 

(

Vi

Vk ) ;

T T ( , P)

x

 

 

 

 

 

ik

 

 

m

 

 

 

 

 

 

2

 

 

x

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

k

i

 

О схожести уравнений гидродинамики и электрогидродинамики говорят даже сравнительно общие характеристики среды.

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

магнитную вязкость [4.18].

Аналогично числу Рейнольдса движущего потока введено понятие

магнитного числа Рейнольдcа.

Если, например, Rm > 1, то в уравнениях движения можно пренебречь электрической силой по сравнению с магнитной.

4.3 Модель процесса пневмотранспортирования семян

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

Негативное влияние семяпровода на равномерность распределения уже давно было выявлено В.Е. Комаристовым и затем многими другими исследователями.

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

93

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

Преимуществами таких аппаратов является универсальность, отсутствие на сеялке (посевной секции) вращающихся частей (цепей, звездочек, дисков), использование пневмоавтоматики для целей регулирования процесса и контроля за качеством высева.

Немаловажным является опыт решения аналогичных задач в теории контейнерного трубопроводного транспорта, в частности, организации пневмопочты с очень сложными системами трубопроводов, соответствующих произвольному профилю трассы, путевыми сбросами, эжекторами, байпасными устройствами при высокой плотности элементов [4.13], [4.14].

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

Устройство, обеспечивающее выброс семян из пневмопровода с заданной скоростью и снижение отрицательного воздействия воздуха в посевной борозде изучено в Воронежском ГАУ [4.19] сотрудниками кафедры сельскохозяйственных машин А.П. Тарасенко, В.В. Василенко, В.А. Гребцовым и др.

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

Рис.4.10 Схема высева с пневмотранспортированием семян:

1 – высевающий аппарат; 2 – семяпровод; 3 – выбросное устройство; 4– сошник

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

При высеве семян кукурузы оптимальными оказались семяпроводы с внутренним диаметром 20 мм, выполненными по эллиптической или параболической кривой.

Диаметр витка выбросного устройства в этом случае должен равняться 200 мм. Коэффициент живого сечения радиальных отверстий необ-

94

ходимо изменять в пределах 0,06…0,25. Оптимальной скоростью воздушного потока оказались величины 28…32 м/с.

Скорость выброса семян кукурузы в этом устройстве может регулироваться в пределах 1,6…3,5 м/с.

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

Большее нарушение равномерности исходного потока характерно для малой нормы высева. Коэффициент вариации интервалов времени вылета отдельных семян увеличивается на 3…11%. Дробление семян не превышает

0,3%.

Дисперсионный анализ технологического процесса показал, что при равенстве скорости выброса семян с рабочей скоростью посевного агрегата относительное влияние пневмосемяпровода на точность распределения семян составляет около 60%, относительное влияние перекатывания семян по дну борозды не превышает 10%.

Увеличение длины семяпровода снижает равномерность распределения семян и при длине трубок, равных 0,74 м, высев находится уже на уровне серийных пунктирных сеялок (σс = 10…13 см).

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

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

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

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

Уравнение, предложенное автором в работе [4.19], не учитывают ряда существенных факторов. В них не учитывается величина напора воздуха (давления) ни в начальной точке, ни в любом из последующих положений. В уравнения не входят даже такой существенный параметр, как внутренний диаметр трубопровода.

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

95

4.3.1 Основные уравнения одномерного течения воздуха в семяпроводе (трубе)

Как и в случае изучения течения воздуха около семян, в основе описания движения воздуха в семяпроводе могут быть использованы уравнения механики сплошной среды [4.13], выражающие законы сохранения массы газа (4.3) приводящей к уравнению неразрывности потока (4.4), количества движения (4.5), изменения количества движения (4.6) и энергии

(4.8).

В интегральной форме применительно к движению по трубе эти

уравнения могут быть представлены как

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 0

 

 

 

 

(4.52)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

Vd

 

Pnd

 

Fd

 

 

(4.53)

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

 

(U

V 2

)d

 

 

P Vd

 

F Vd

 

q d ,

(4.54)

 

dt

 

 

 

 

 

 

 

2

 

 

 

 

 

n

 

n

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

A

 

 

 

 

A

 

 

где η(dη) –произвольный объем сплошной среды, состоящий из одних и тех же частиц: А()- поверхность, ограничивающая этот объем, ρ - плотность газа, V - скорость газа, Рn – вектор напряжений на площадке с нормалью n, F – плотность массовых сил , U – внутренняя энергия единицы объема газа, qn –нормальная составляющая вектора потока тепла.

Преобразуя поверхностные интегралы в формулах (4.52), (4.54) в интегралы по объему η с помощью формулы Гаусса–Остроградского, можно получить:

 

 

 

 

 

 

d

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

d 0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

d divPd Fd

(4.55)

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

( U

V 2

)d div(PV )d divqd F Vd ,

(4.56)

dt

2

 

 

 

 

 

 

 

 

 

 

 

 

где Р - тензор напряжений в газовой среде;

Р = -РI + T (здесь Р – давление газа; I – единичный тензор; Т – тензор касательных напряжений).

Используя известную формулу дифференцирования по времени интеграла с переменными пределами интегрирования

d

A d (

A

divAV)d ,

dt

t

 

 

 

 

 

 

можно уравнения (4.52), (4.55), (4.56) привести к виду:

( div V )d 0;

t

96

 

V

div( V V PI ) d divTd Fd ;

 

 

 

t

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V

2

 

 

 

P

 

V

2

 

 

 

 

 

( U

 

) div

 

U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

2

 

 

 

 

 

 

2

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

div(VT q) F V d

где ρV V –диадное произведение векторов.

Так как область интегрирования η выбрана произвольно, приведенные выше соотношения равносильны дифференциальным уравнениям неразрывности:

 

 

 

 

 

div V 0 ;

(4.57)

 

 

 

 

t

движения

 

 

 

 

 

 

 

 

 

 

 

div( V V PI ) divT F;

(4.58)

 

 

t

 

энергии

 

 

 

 

 

 

 

 

 

 

V 2

 

 

V 2

 

t

( U

 

) div[( i

 

 

)V ] div(TV q) F;

(4.59)

2

2

 

 

 

 

 

где i = U + P/ρ – энтальпия среды.

Уравнение (4.59) можно преобразовать с помощью уравнений (4.58) и неразрывности (4.57), исключив из него производную от кинетической энергии. При этом получается уравнение

 

dU

 

P d

T V divq ,

(4.60)

 

 

 

 

dt

dt

 

 

 

 

которое часто называют [4.14] уравнением притока тепла.

Система дифференциальных уравнений одномерного течения газа в трубе

Пусть ось ОХ будет направлена вдоль трубы, а остальные оси расположены в плоскости, перпендикулярной к оси трубы.

Одномерные уравнения трубопроводной газовой динамики, т.е. такие уравнения, в которых все функции представляют собой среднее по сечению значения соответствующих величин и зависят только от координаты Х и времени t, получают из системы уравнений (4.57)…(4.59) осреднением по площади поперечного сечения трубки.

Для этого используют операцию

A(x,t) A(x,y,z,t

4

 

A(x, y, z,t)d

,

d

2

 

 

s

 

 

 

 

 

где S – площадь сечения трубы; d –диаметр трубы.

При осреднении обычно применяют допущения о том, что сжимаемость газа много больше упругости трубы, плотность газа его давление,

97

внутренняя энергия и температура мало изменяются по сечению трубки, скорость газа на внутренней поверхности трубы равна нулю [4.14].

Осредненные по сечению трубы уравнения имеют вид:

 

 

 

 

 

 

 

Q

0

 

 

 

 

 

 

 

 

 

 

 

 

t

x

 

 

 

 

 

 

 

 

Q

 

 

P V 2

 

4

 

 

F ;

t

 

 

 

 

 

 

 

x

 

 

 

 

d

 

 

в

 

x

 

du

 

P d

 

 

 

4

 

 

 

qx

 

 

 

 

 

 

 

T V

 

qn

 

dt

dt

 

d

x

Q V

, (4.61)

;

где Q(x,t) –расход газа; V - средняя скорость газа; ηв - среднее касательное напряжение на внутренней поверхности трубки; β – коэффициент, учитывающий неравномерность распределения скоростей газа по сечению трубы (для развитых турбулентных течений β ≈ 1); qn – осредненный по периметру сечения тепловой поток; qх – тепловой поток вдоль оси трубы.

Гидравлическое приближение. Для замыкания системы уравнений (4.61) необходимо ввести дополнительные допущения, чтобы компенсировать потерянную при осреднении информацию. Эти допущения составляют гидравлическое приближение и выражаются в следующем.

1.Касательное напряжение на внутренней поверхности трубопровода связано со средними параметрами течения формулой Дарси

U

 

1

 

V

 

V

 

 

,

(4.62)

 

 

 

 

 

 

 

 

 

 

d

d

2

 

 

 

 

 

 

 

 

 

в которой безразличный коэффициент гидравлического сопротивления λ имеет тот же вид, как и для стационарного движения (гипотеза квазистационарности).

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

В соответствии с гипотезой квазистационарности скорость диссипации механической энергии движения газа за счет сил внутреннего трения имеет вид:

T V

1 V 3

.

(4.63)

d

 

2

 

 

 

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

M VC 1.

Это допущение позволяет не учитывать в уравнении движения конвективный перенос количества движения, т.е. величину βρV2, по сравнению с колебаниями давления.

98

Согласно формуле Н.Е. Жуковского, колебания давления δρ можно

оценить равенством

с Q ,

а колебания величины βρV2 равенством

( V 2 ) V Q ,

где δQ – колебания расхода газа.

Отношение колебаний рассмотренных величин, составляет

( V 2 ) V ,C

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

3.Тепловой поток вдоль оси трубы мал по сравнению с потоком тепла в перпендикулярном направлении и его можно не учитывать.

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

P RTг ;

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

U CV Tг const ;

i CpTг const ,

где R - газовая постоянная; CV, Cp – удельные теплоемкости изохорного и изобарного процессов; Тг - температура газа.

С учетом сделанных допущений систему уравнений (4.61) и добавочные условия для U и P можно представить в виде

 

 

 

 

 

 

 

Q 0;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

 

 

P

 

 

 

 

 

 

V

 

V

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

;

 

 

 

 

 

 

t

 

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(4.64)

 

du

 

 

P

 

d

 

1

 

 

V

3

 

4qn

 

 

 

 

 

 

 

 

;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dt

 

 

 

 

dt

 

d

 

 

2

 

 

 

 

 

 

d

 

 

P RT ;

 

u C T const;

 

 

Q V

 

 

 

г

 

 

 

 

 

V г

 

 

 

 

 

 

 

 

 

 

 

 

 

 

При выводе этих уравнений принято, что проекция массовых сил, действующих на газ, равна нулю. Система уравнений (4.64) является уже замкнутой, т.е. позволяет определить шесть неизвестных величин

, Р,V,Q,Тг ,U .

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

Допущение о малости числа Маха приводит к уравнению

 

 

di

 

 

d

 

 

4qn

;

 

dt

dt

 

d

 

 

 

 

 

 

 

 

или

 

dTв

 

P

 

 

4qn

 

 

Cn dt

t

 

 

d

 

 

 

 

 

 

99

Раскрывая в последней формуле полную производную температуры по времени, находим

C

 

Tв QC

Tв

Р

 

4qn

 

t

 

 

р

t

p x

 

d

Полная система уравнений для определения функций (x,t), P(x,t), Q(x,t) и Tв(x,t) будет иметь вид:

 

 

 

Q 0;

 

 

 

 

 

 

 

 

t

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Q

 

Р

 

 

Q

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

;

 

 

 

 

t

x

 

 

 

 

 

 

 

 

 

2d

 

 

 

 

 

 

 

T

 

 

Tв

р

 

 

 

4qn

.

C

 

QC

 

;

 

 

 

 

 

 

 

P

t

 

 

p x

t

 

 

 

d

 

 

 

 

P RTb

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

четырех неизвестных

(4.65)

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

Возможна еще одна упрощенная форма уравнений (4.65), описывающих движение газа в условиях, близких к изотермическим [4.13].

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

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

При этом колебания температуры определяются по формуле:

dTв

 

1

P

.

(4.66)

dt

Ср

t

 

 

Исключив из уравнения неразрывности плотность, с помощью уравнения состояния получают:

P

R

Tв PT Q

0 .

t

 

t

в x

 

После подстановки в это уравнение частной производной температуры по времени, определяемую по формуле (4.66), можно найти

P

(1

R

) RT

Q

0 ,

t

 

 

 

C

p

в x

 

 

 

 

 

 

 

или, учитывая формулу Майера R = CP – CV,

 

 

P

RT

Q

0

,

 

t

 

 

в x

 

 

где СP / СV - показатель адиабаты.

100

Выражение RTв представляет собой адиабатическую скорость

звука в газе:

C 2 ( dPd )s RTb ,

где производную определяют при постоянной энтропии S.

Система уравнений (4.65) для рассматриваемого случая примет вид:

 

P

С

2

Q

0

 

 

 

 

 

t

 

x

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.

(4.67)

Q P

C

2

 

Q

 

Q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

x

 

2 d

 

P

 

 

 

 

 

 

 

 

 

 

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

Система уравнений (4.67) относительно двух функций – давления P(x,t) и расхода Q(x,t) – может быть использована в качестве исходной для расчета движения газа по трубам в изотермических условиях.

4.3.2 Математическая модель движения семян по семяпроводу

Пусть семяпровод представляет собой трубку малого диаметра с произвольным профилем, движение газа в котором (при отсутствии семян) описывается системой уравнений (4.67).

Рассмотрим участок семяпровода, расположенного между сечениями х=0 и х=L, где L – длина участка, отсчитываемого вдоль оси трубки х.

Пусть в некоторый момент времени на этом участке находится одно или несколько семян. Движение каждого семени описывается неизвестной и подлежащей определению функцией x j (t) (здесь j=1,2…N), которую

будем называть траекторией j-го семени.

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

 

dV

 

 

 

 

m j

j

P ( j ) P ( j ) S m j g sin

( j ) Fcj ,

(4.68)

dt

 

 

 

d j

V

 

,

 

(4.69)

 

 

 

 

j

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

где mj - масса j-го семени; Vj - скорость j-го семени;

P ( j ) и P ( j ) - дав-

ление газа за и перед семенем; g - ускорение свободного падения; α(х) - угол наклона оси семяпровода к горизонтали в той точке трассы x j (t) ,

в которой находится j-е семя в данный момент времени; Fcj - сила сопротивления, вызванная трением.

Уравнения (4.68) и (4.69) незамкнуты, и их нельзя решить отдельно, так как в каждое из них входят неизвестные P ( j ) и P ( j ) , значения которых определяются в процессе решения газодинамических уравнений.

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