
книги из ГПНТБ / Цейтлин Г.М. Аэродинамика и динамика полета самолета с ТРД учебник
.pdfВозвратившись к выводу уравнения (2.2-1), нетрудно устано вить, что правая часть записанной выше пропорции есть отношение
. |
A- |
D |
V |
произведении |
ат—^-, |
которые в соответствии с известным из ме |
ханики принципом Даламбера можно назвать «силами инерции», препятствующими изменению скорости элементарных воздушных частиц в двух динамически подобных явлениях. Переходя к про изводной пропорции
|
.„ |
|
2 |
dVt |
dV2 |
||
—j—dbi |
u . a - 3 |
— d b 3 |
мы видим, что в двух явлениях, подобных по силам вязкого трения, отношения «сил инерции» к силам трения должны быть одинако выми. Производя сокращения с учетом соотношений, вытекающих
|
|
|
|
|
|
|
IS |
|
S |
из наличия геометрического и кинематического подобия (j^1 |
— |
j ^ - ; |
|||||||
V] |
V2 |
^ |
и заменяя |
элементарные |
1 |
к |
ли |
||
•jy-— |
"^т/")> |
отрезки нормалей |
|||||||
ниям |
тока |
dti\ |
и dn2 |
пропорциональными |
им характерными |
линей |
|||
ными размерами 1\ и h, |
получаем |
|
|
|
|||||
|
|
|
|
|
PiVi^i |
__ |
|
|
(26) |
|
|
|
|
|
(A, |
(J.2 |
|
|
> |
Безразмерное отношение «сил инерции» к силам вязкого |
трения |
||||||||
|
|
|
|
|
R e ~ i £ U - £ |
|
(2.7) |
||
называют |
ч и с л о м |
Р е й н о л ь д е а. |
|
|
|
||||
Согласно уравнению |
(2.6) в |
явлениях, |
подобных по силам |
тре |
ния, числа Рейнольдса одинаковы. Повторяя рассуждения в обрат ной последовательности, можно убедиться, что равенство чисел Рейнольдса (Rei = Re2) является не только необходимым, но и до статочным условием подобия аэродинамических явлений по силам трения. Другими словами, число Рейнольдса является критерием подобия аэродинамических явлений по силам трения.
Чем меньше число Рейнольдса, тем больше силы трения, по нуждающие частицы к изменению скорости, по сравнению с силами инерции, препятствующими изменению скорости. При малых чи слах Re силы трения преобладают над силами инерции и оказы вают существенное влияние на характер всего потока. При боль
ших числах |
Re |
(соответствующих |
реальным скоростям |
полета) |
||
вязкость воздуха |
проявляется только |
в непосредственной |
близости |
|||
от |
поверхности |
тела. Коэффициент |
сопротивления трения тела |
|||
сх?р |
— ~~^ |
с увеличением числа |
Re |
уменьшается. |
|
Теперь предположим, что в двух динамически подобных явле ниях действуют только силы давления. Так как элементарная сила
50
давления dP = pdS, то в соответствии с уравнением (2.2-1) можно записать
Разделяя параметры, относящиеся к натуре и модели явления,
и учитывая, что в силу геометрического подобия |
= -ф-, |
полу- |
|
чаем |
|
|
|
Pi |
Pi |
|
|
Отношение — можно определить из формулы |
скорости |
звука |
a2 = kRT= k ~ : -£p = JL^ Тогда полученное выше равенство, вы ражающее условие подобия явлений по силам давления, прини
мает вид |
4 ( т. е. Mj = Ма . |
«1 |
|
1 |
2 |
Таким образом, в явлениях, подобных по силам давления, чи сла М должны быть одинаковыми. Иначе говоря, число М являет ся критерием подобия аэродинамических явлений по силам дав ления.
Из сказанного очевидно, что динамическое подобие аэродинами ческих явлений достигается соблюдением геометрического и кине матического подобия, а также равенством чисел Re и М. При на личии этих условий все одноименные аэродинамические коэффици енты одинаковы. Если изменять критерии подобия, то, естественно, будут изменяться и аэродинамические коэффициенты. Другими сло вами, безразмерные коэффициенты аэродинамических сил и момен тов являются функциями безразмерных критериев подобия.
Зависимости аэродинамических коэффициентов от критериев геометрического, кинематического и динамического подобия назы вают а э р о д и н а м и ч е с к и м и х а р а к т е р и с т и к а м и т е л . Исследование аэродинамических характеристик крыла и других частей самолета, а также самолета в целом — основной предмет прикладной авиационной аэродинамики.
§ 2.4. Пограничный слой при малых числах М
Поверхность тела, как бы хорошо она ни была обработана, имеет неровности — выступы и впадины, размеры которых намного больше размеров молекул воздуха. Это и понятно. При самой тон кой обработке поверхности с нее снимаются частицы материала, состоящие, по крайней мере, из нескольких молекул.
51
Пусть около такой поверхности движется воздушный поток с некоторой скоростью V (рис. 2.4). Так как кроме движения всей воздушной массы относительно тела внутри нее имеется еще хаоти ческое (тепловое) движение молекул, то скорость любой молекулы можно представить векторной суммой
VH = V+V.
где Ут — скорость теплового движения данной молекулы в данное мгновение.
За счет скорости VT молекулы отклоняются от общего направле ния движения воздушной массы и сталкиваются с неровностями стенки. Так как элементарные участки поверхности, о которые уда
Рис. 2.4. К объяснению образования погра ничного слоя
ряются молекулы, расположены совершенно случайно относительно векторов VM, то молекулы после столкновения со стенкой отлетают в случайных направлениях. При этом упорядоченное движение мо лекул разрушается, превращаясь в хаотическое, кинетическая энер гия воздушной массы переходит в тепло, а скорость направленного движения воздуха непосредственно на поверхности тела падает до нуля.
Поскольку на поверхности тела скорость воздуха равна нулю, а на некотором расстоянии от этой поверхности поток движется со скоростью V, то в результате проявления вязкости воздуха здесь образуется переходный слой (см. § 1.5).
Слой воздуха, непосредственно прилегающий к поверхности те ла, в котором действуют силы вязкого трения и скорость постепен
но нарастает |
от нуля до |
скорости внешнего потока, называют по |
г р а н и ч н ы м |
с л о е м . |
|
Общие закономерности, действующие в пограничном слое, удобнее всего рассмотреть на примере тонкой плоской пластинки, параллельной линиям тока невозмущенного потока (рис. 2.5). Око ло такой пластинки внешний поток во всех точках имеет одинако вые скорости и давления, чем исключается воздействие на погра ничный слой различных дополнительных факторов.
Как уже говорилось, проявление вязкости воздуха определяется
52
числом Re. Для исследования состояния пограничного слоя в раз личных его точках обычно пользуются местным числом Рейнольдса:
при подсчете которого в качестве характерного линейного размера используется расстояние х данной точки от передней кромки пла стинки.
Если число М невелико, можно считать, что плотность и темпе ратура воздуха в пограничном слое постоянны и, следовательно, не меняется коэффициент кинематической вязкости v. Так как ско-
X
X
Рис. 2.5. Пластинка, |
параллельная ли |
Рис. 2.6. Пограничный слой пластинки |
ниям |
тока |
|
рость на внешней границе пограничного слоя пластинки тоже по
стоянна и равна скорости невозмущенного |
потока, |
то |
при У» — |
||||
= const число Re* пропорционально |
координате |
х, |
т. е. увеличи |
||||
вается по линейному |
закону |
с удалением от |
передней кромки. |
||||
Вблизи передней |
кромки |
число Rex |
мало. Это значит, |
что силы |
|||
трения велики по сравнению |
с инерционными |
силами. |
Поэтому |
скорость быстро нарастает вдоль нормали к поверхности и уже на
небольшом расстоянии от нее достигает |
скорости внешнего пото |
ка — толщина пограничного слоя 8 здесь |
мала. |
Структура пограничного слоя на передней части пластинки ла минарная (слоистая) — перемешивания воздуха между смежными элементарными струйками практически нет.
С удалением от носка пластинки число Re* постепенно увеличи вается и силы трения, вызывающие изменения скорости по толщи не пограничного слоя, играют все меньшую роль по сравнению с инерционными силами, препятствующими этим изменениям. Следо вательно, нарастание скорости с удалением от поверхности замед ляется, а толщина пограничного слоя увеличивается.
Выделив в пограничном слое произвольную воздушную частицу (например, в виде кубика, рис. 2.6) и представив себе профиль скоростей, легко обнаружить, что сторона частицы, расположенная дальше от стенки, движется быстрее, чем сторона, расположенная
53
ближе к стенке. Это приводит к деформации |
и вращению |
частиц, |
|
т. е. к образованию вихрей. |
|
|
|
С увеличением местного |
числа Рейнольдса |
(с удалением |
от но |
ска пластинки) количество, интенсивность и размеры вихрей воз
растают и они начинают взаимодействовать между собой. В ко |
|
нечном счете это приводит к |
разрушению ламинарной структуры |
и переходу ее в турбулентную |
(завихренную), характеризующуюся |
интенсивным перемешиванием воздуха по всей толщине погранич ного слоя.
Протяженность |
Ахт |
переходного |
(от ламинарного |
к турбулент |
|||||
ному |
состоянию) |
участка пограничного |
слоя |
сравнительно |
неве |
||||
лика, |
и при решении |
практических |
задач |
этот участок |
обычно |
||||
|
|
|
сводят |
к так называемой |
точке |
пере |
|||
|
|
|
хода Т. |
|
|
|
|
|
|
|
|
|
Характерно, что разрушение |
лами |
|||||
|
|
|
нарной |
структуры |
пограничного слоя |
||||
|
|
|
происходит |
при определенном |
значе |
||||
|
|
|
нии местного числа Рейнольдса, кото |
||||||
|
|
|
рое называют |
критическим: |
|
||||
|
|
|
|
Re,кр • |
|
|
|
(2.9) |
Величина критического числа Рей нольдса тем больше, чем лучше обра ботана поверхность пластинки и чем меньше турбулентность (завихрен
ность) внешнего потока. В среднем можно считать ReK P =105 . Если известны параметры невозмущенного потока, то, зная ве
личину ReK P , легко определить протяженность ламинарного участка пограничного слоя, равную координате точки перехода:
Минимально возможная скорость полета современного самоле
та составляет примерно |
200 км/ч = 55 |
м/с. Такая |
скорость дости |
||||||||
гается на |
посадке, |
т. |
е. |
непосредственно |
у |
земли, |
где |
v = |
|||
|
|
|
|
|
|
105.] 4 5 . Ю - 6 |
|
|
|||
= 1,45-10~5 |
м 2 / с |
Для |
этого |
случая хт — |
'-^ |
|
= 0,027. Из |
||||
рассмотренного |
примера |
видно, что даже при минимальной |
ско |
||||||||
рости значение координаты х т намного меньше |
|
протяженности |
|||||||||
любой из основных частей самолета |
в направлении |
оси |
Ох. Это |
||||||||
значит, что в полете на большей части |
внешних поверхностей |
само |
|||||||||
лета пограничный слой имеет турбулентную структуру. |
|
|
|||||||||
В отличие от ламинарного пограничного слоя, в котором |
пере |
||||||||||
нос количества движения по нормали к линиям |
тока |
осуществляет |
ся только за счет хаотического движения молекул, в турбулентном слое благодаря интенсивным вихревым движениям вдоль нормали перемещаются не только отдельные молекулы, но и сравнительно большие воздушные массы. Это приводит к тому, что с удалением
54
непосредственно от стенки (в нижней части пограничного слоя) скорость нарастает значительно быстрее. По той же причине затор моженные воздушные частицы выносятся значительно дальше от стенки в воздушный поток, вследствие чего толщина турбулентного пограничного слоя больше. Так как интенсивность вихревых дви жений мало зависит от местного числа Рейнольдса, то увеличение темпа нарастания скорости в нижней части турбулентного погра ничного слоя с увеличением числа Rex происходит значительно медленнее, чем в ламинарном пограничном слое. Сравнение про филей скоростей в турбулентном и ламинарном пограничных слоях схематично показано на рис. 2.7.
§ 2.5. Нагрев воздуха в пограничном слое. Кинетический
нагрев конструкции летательного аппарата
Как уже говорилось, в пограничном слое воздух тормозится до полной остановки непосредственно около поверхности тела. По скольку кинетическая энергия воздуха при этом переходит в тепло, температура в нижней части пограничного слоя повышается.
Если считать процесс торможения воздуха в пограничном слое адиабатным, то его температуру на поверхности тела можно опре делить по формуле (1.23-2):
Т 0 п = Тт + 6,2Тт№т.
Прирост температуры пропорционален числу М^:
&Tu = Tn-Tm=0,2rm№m.
В реальных условиях главным образом за счет конвективного отвода тепла от нижней области пограничного слоя в менее на
гретую верхнюю область и внешний поток часть тепла |
рассеивает |
||||||||||
ся и действительное |
повышение |
температуры АГ |
оказывается |
||||||||
меньшим, |
чем при отсутствии |
теплообмена. |
Отношение |
г— |
|||||||
|
|
|
|
|
|
|
|
|
|
|
ал |
называют |
к о э ф ф и ц и е н т о м |
|
в о с с т а н о в л е н и я |
т е м п е |
|||||||
р а т у р ы . |
Как показывают |
теоретические |
и экспериментальные |
||||||||
исследования, для ламинарного пограничного слоя г=0,85 |
и для |
||||||||||
турбулентного г=0,9. |
Поэтому |
повышение |
температуры |
в |
нижней |
||||||
части пограничного слоя можно |
определять |
|
по формуле |
|
|
||||||
|
ЬТ = 0,2гТ |
М2 в |
0,1757 |
М 2 . |
|
|
(2.10) |
||||
|
|
СС |
00 |
' |
ОС |
|
ОС |
|
\ |
/ |
|
График |
этой зависимости |
показан |
на рис. 2.8. |
|
|
|
От воздуха тепло передается конструкции летательного аппа рата (кинетический нагрев конструкции) и через некоторое время она прогревается до температуры нижней части пограничного слоя. В связи с этим уже в современном самолетостроении возникает ряд новых сложных технических проблем. При температурах 60—80° С, что соответствует числам М.= 1,21,3, теряет свою проч ность и становится непригодным для фонарей кабин органическое
55
стекло. В длительном полете с числом М = 2-г-2,5 обшивка самолета прогревается до 200—300° С и алюминиевые сплавы становятся не пригодными для ее изготовления. При еще больших числах М вы ходят на предел прочности титановые сплавы, а затем и жаропроч ные стали, возникает необходимость в отыскании новых, более жа ропрочных конструкционных материалов.
Не менее сложной оказывается задача охлаждения кабин эки пажа, отсеков с топливными баками и различным оборудованием, чувствительным к высоким температурам. Дело в том, что при
высоких |
скоростях |
полета системы охлаждения, основанные на |
||||||||||||||
|
|
|
|
непосредственном |
использовании |
холодного |
||||||||||
|
|
|
|
забортного воздуха, непригодны, ибо (в ре |
||||||||||||
|
|
|
|
зультате |
торможения) |
в |
воздухозаборни |
|||||||||
|
|
|
|
ках и радиаторах |
воздух |
нагревается |
|
так |
||||||||
|
|
|
|
же, как и в пограничном |
слое. Приходится |
|||||||||||
|
|
|
|
применять |
холодильные |
установки, |
осно |
|||||||||
|
|
|
|
ванные |
на |
механическом |
срабатывании |
|||||||||
|
|
|
|
теплоперепада |
|
(турбохолодильники) |
и |
на |
||||||||
|
|
|
|
испарении |
специальных |
жидкостей. |
|
Та |
||||||||
|
|
|
|
кие |
установки |
|
имеют |
|
определенные |
га |
||||||
|
|
|
|
бариты |
и |
вес, |
потребляют |
значительные |
||||||||
|
|
|
|
мощности. |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
Одним |
из |
наиболее |
кардинальных |
ре- |
||||||||
и |
j |
g j |
Мм шений |
проблем, |
связанных с |
кинетическим |
||||||||||
|
|
|
|
нагревом |
конструкции |
летательного |
аппа- |
|||||||||
Рис. 2.8. Повышение тем- |
рата, |
является |
увеличение высоты полета, |
|||||||||||||
пературы |
|
пограничного |
В в е р |
х н и |
х |
С Л О Я х атмосферы, |
на |
высотах |
||||||||
|
|
|
|
60—80 км и более, воздух сильно |
разрежен, |
|||||||||||
|
|
|
|
в связи с чем коэффициент |
теплопередачи |
|||||||||||
от пограничного слоя к стенке |
столь |
мал, |
что |
даже |
длитель |
|||||||||||
ный полет с числами Л4^8 - М0 не вызовет опасного нагрева |
|
кон |
||||||||||||||
струкции. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Предполагается, что, применяя специальные теплоизолирующие покрытия и системы охлаждения кратковременного действия, в ко торых пары охлаждающей жидкости выводились бы непосредствен но в поток без конденсации для повторного использования, можно будет обеспечить полет с числами М ~ 4 ч - 6 и на значительно мень ших высотах (10—20 км). Разумеется, безопасная продолжитель ность такого полета строго ограничена.
Нагрев воздуха в пограничном слое приводит к увеличению коэффициента динамической вязкости р.. Казалось бы, что за счет этого с увеличением числа М т скорость должна быстрее нара стать вдоль нормалей к поверхности тела и, следовательно, тол щина пограничного слоя должна уменьшиться. Однако это не так.
Поскольку давление по толщине пограничного |
слоя не меняется |
|
и остается равным давлению во внешнем потоке |
(для тонкой пла |
|
стинки— давление рк), то с повышением температуры |
происхо |
|
дит изобарное расширение воздуха. Средние значения |
расстояния |
56
между молекулами и длины их свободного пробега увеличивают ся, то же изменение скорости происходит на большем отрезке нор
мали, и толщина |
пограничного слоя |
увеличивается. |
|
§ 2.6. Сопротивление трения |
|
Разделим пограничный слой пластинки на множество предель |
||
но тонких слоев |
(рис. 2.9). Так как |
скорости в смежных слоях |
различны, между ними действуют силы вязкого трения, равные
количеству |
движения, переносимому |
за единицу времени |
из слоя |
|
в слой. Слой / увлекает слой 2, слой |
2 — слой 3 и т. д. В |
конеч- |
||
|
V |
|
|
|
——Ч |
' |
|
|
|
1 |
* - / |
2 |
|
|
|
—) |
3 |
|
|
|
•У' |
••• |
|
|
|
п-1 |
|
|
|
|
|
п |
|
|
Рис. 2.9. К |
пояснению |
природы сопро |
Рис. 2.10. К определению коэффи |
|
тивления трения |
циента трения |
|
ном счете силы трения замыкаются на поверхности тела в форме
непрерывной |
бомбардировки |
молекулами |
неровностей зтой |
по |
|||
верхности. |
|
|
|
|
|
|
|
Используем |
введенное в |
§ 2.4 представление о |
скорости |
мо |
|||
лекулы как о результирующей скоростей |
хаотического |
движе |
|||||
ния VT и направленного движения V (VM = |
Vr•+ V). |
|
|
|
|||
Если за |
счет |
скорости VT |
вероятность |
и средняя |
сила |
ударов |
молекул на элементарных участках неровностей поверхности тела, обращенных вперед и назад, совершенно одинаковы, то наличие скорости V увеличивает вероятность и среднюю силу ударов моле кул на участках поверхности, обращенных вперед. В результате
суммарная сила воздействия молекул воздуха |
на тело |
приобре |
тает составляющую, направленную назад (в |
сторону |
движения |
потока). Эта составляющая, численно равная секундной потере количества движения воздуха в пограничном слое, и является со
противлением |
трения. |
|
|
||
|
Выделим на поверхности пластинки, имеющей постоянную ши |
||||
рину / и длину b |
(рис. 2.10), |
элементарную площадку |
dSTp=ldx. |
||
В |
соответствии |
с |
формулой |
(1.11) сила трения, действующая на |
|
эту |
площадку, |
будет |
|
|
57
где ( - ^ f ) 0 H р.о — соответственно градиент скорости по нормали и коэффициент динамической вязкости воздуха непосредственно на
поверхности пластинки. |
|
|
|
|
|
|
||
По формуле |
(2.3), |
полученной |
на |
основании |
теории |
подобия, |
||
ту же силу можно представить в виде |
|
|
|
|
||||
|
|
dQrp = |
c'fqjdx. |
|
|
|||
Вошедший в последнюю формулу |
коэффициент |
|
||||||
|
C r ^ = ? : h |
|
w . |
|
( 2 Л 1 ) |
|||
называют м е с т н ы м |
(для данного |
элемента трущейся |
поверхно |
|||||
сти) к о э ф ф и ц и е н т о м |
т р е н и я . |
|
|
|
||||
Общее значение коэффициента трения для данного тела в дан |
||||||||
ном потоке определяется |
осреднением |
местного |
коэффициента c'f |
|||||
по всей трущейся |
о воздух |
поверхности |
ST p : |
|
|
|||
|
|
|
5 |
|
|
ь |
|
|
При малых числах М м коэффициент динамической вязкости воздуха цо= Iх„ —const и изменение коэффициента c'f вдоль пла-
f |
dV\ |
стинки определяется изменением градиента |
y-j^j- |
Как было показано, при ламинарном пограничном слое увели чение местного числа Re* приводит к более плавному нарастанию скорости по толщине пограничного слоя, т. е. к уменьшению гра-
диента |
соответственно |
коэффициент |
cf с удалением от |
передней кромки пластинки уменьшается. |
|
||
При определении общего числа Рейнольдса, |
характеризующего |
||
обтекание |
данного тела данным |
потоком, в |
качестве линейной |
базы используется какой-либо характерный размер тела. Для пря моугольной пластинки таким размером является ее длина Ь в на правлении оси Ох. Таким образом,
R e = - = - .
Отношение местного числа Re* к общему числу Re есть отно сительная координата данной точки:
Re* х
Не Ь
Очевидно, что при Re<ReK p ламинарная структура погранич-* ного слоя сохраняется на всей пластинке. В этом случае выполне-
58
ние интегрирования в уравнении (2.12) (приводим без доказатель ства) дает результат
1.328
Если Re>ReK P , то на пластинке имеются и ламинарный, и тур булентный участки пограничного слоя. Турбулизация погранич ного слоя сопровождается интенсивным увеличением градиента
("Ж")о' |
3 |
с л е д о в |
а т е л ь |
н о > и |
местного коэффициента трения с'}, |
|||||||
С |
дальнейшим |
увеличением |
числа Re* (свыше |
ReK p ) |
коэффи |
|||||||
циент |
c'j постепенно |
падает, |
но медленнее, |
чем |
в |
ламинарном |
||||||
слое. |
|
|
|
|
|
|
|
|
|
|
|
|
|
При |
больших |
и средних |
скоростях |
полета |
Re^Re l t p . В |
этом |
|||||
случае |
относительная |
протяженность |
ламинарного |
участка |
хт = |
|||||||
|
х |
Re Kp |
|
|
|
|
|
|
|
|
|
|
= |
— = |
—j7^- невелика и можно считать, что весь |
пограничный |
|||||||||
слой имеет |
турбулентную структуру. Тогда |
(приводим |
тоже без |
доказательства) для осредненного коэффициента трения пластинки получается выражение
0.455 |
|
(2.14) |
с = Ц « » |
|
|
(lg Re)2 |
|
|
При малых скоростях полета, особенно |
на больших |
высотах, |
где коэффициент кинематической вязкости |
v = - t - с и л ь н о |
возра |
стает за счет падения плотности воздуха, числа Re не намного превышают ReK p . В этих условиях ламинарный участок погранич ного слоя по своей протяженности соизмерим с турбулентным и формула (2.14) дает завышенные значения сх
При оценке влияния сжимаемости воздуха на коэффициент трения необходимо иметь в виду, что в реальных условиях полета применительно к основным частям самолета, имеющим значитель ную длину в направлении линий тока, большим числам М„ со ответствуют и большие числа Re^>ReK p . Кроме того, при боль ших числах Мм сопротивление трения обычно невелико по сравне нию с сопротивлением давления. Поэтому при оценке влияния сжимаемости воздуха в практических авиационных задачах допу стимо считать пограничный слой турбулентным и пользоваться
полученной для |
этих условий |
приближенной |
формулой |
||||
|
|
Cf |
~ С/нс ^ + о д е ш ^ ) 0 , 6 7 " * |
|
|||
Здесь |
с / н с — коэффициент |
трения |
пластинки, определяемый по |
||||
формуле |
(2.14) |
при действительном |
числе |
Re |
без учета сжимае |
||
мости воздуха (в несжимаемом потоке). |
|
|
|||||
Поверхности |
частей |
самолета обычно |
имеют небольшую кри |
||||
визну в |
направлении |
линий |
тока. |
Поэтому |
по коэффициентам |
59