книги / Теория вертолета. Кн. 1
.pdf262 Глава 5
диске винта в продольном направлении. Так как изменение ин дуктивной скорости в поперечном направлении влияет на коэф фициент Pic не столь существенно, теория, основанная на рав номерном распределении индуктивных скоростей, позволяет вполне надежно рассчитать продольный наклон ПКЛ. Теорети
ческие зависимости махового движения от общего шага |
и угла |
||||||||||
|
|
|
|
наклона |
|
вала |
винта |
||||
|
|
|
|
также |
вполне |
надеж |
|||||
|
|
|
|
ны. |
Таким |
образом, |
|||||
|
|
|
|
расхождение |
теорети |
||||||
|
|
|
|
ческих |
|
и |
эксперимен |
||||
|
|
|
|
тальных |
|
результатов |
|||||
|
|
|
|
обусловлено |
главным |
||||||
|
|
|
|
образом |
|
неравномер |
|||||
|
|
|
|
ностью |
|
распределения |
|||||
|
|
|
|
индуктивных |
скоростей |
||||||
|
|
|
|
при |
полете |
вперед. |
|||||
|
|
|
|
Харрис |
|
сопоставил |
ве |
||||
|
|
|
|
личины pls, полученные |
|||||||
Рис. 5.39. Сопоставление |
экспериментальных |
в эксперименте при ма |
|||||||||
лых |
р, |
|
Ст/а = |
0,08 и |
|||||||
величин поперечного наклона конуса лопастей |
|
||||||||||
с величинами, рассчитанными |
при Ст1о = 0,08, |
апкл |
= |
|
Г, |
|
с |
величи |
|||
а п к л = 1° и различных распределениях индук |
нами, |
|
рассчитанными |
||||||||
тивных скоростей [Н.48]. |
|
|
по |
различным |
тео |
||||||
J — равномерное распределение; |
2 —линейное распре |
риям: |
1) |
с |
равномер |
||||||
деление при 6 x ss2 ii; 3 —линейное распределение при |
ным |
|
распределением |
||||||||
kx=tg (х/2); 4—по |
вихревой |
теории для равномерно |
|
||||||||
нагруженного диска; 5—по |
вихревой теории для не |
индуктивных |
скоро |
||||||||
равномерно нагруженного диска; |
6 — распределение, |
стей, |
2) |
с |
линейным |
||||||
"индуцируемое следом заданной |
формы; О —экспери |
||||||||||
ментальные точки. Рисунок воспроизведен с разреше |
распределением |
|
при |
||||||||
ния Ф. Харриса и |
Американского вертолетного обще |
|
|||||||||
ства. |
|
|
|
kx = |
2р, 3) |
с линейным |
|||||
A* = tg(x/2) |
[С.78], 4) |
|
распределением |
|
при |
||||||
по вихревой теории для равномерно на |
груженного диска [С.35], 5) по вихревой теории для неравно мерно нагруженного диска [Н.84] и, наконец, 6) по теории, в которой распределение индуктивных скоростей было найдено численно по схеме следа, состоящего из дискретных вихрей заранее заданной формы (результаты вихревой теории винта рассмотрены в гл. 4). На рис. 5.39 видно, что учет неравномер ности распределения индуктивных скоростей действительно при водит к более правильной оценке Pis как по величине, так и по характеру изменения в зависимости от р, но все-таки попереч ный наклон ПКЛ получается существенно заниженным. Это на-
•водит на мысль, что теория несущего винта не позволяет пока надежно рассчитать даже первые гармоники неравномерного
распределения индуктивной скорости, которые определяют на клон ПКЛ. Харрис высказал предположение, что основной при чиной остающихся расхождений является использование схемы
Полет вперед И |
263 |
следа заранее заданной формы. Индуцируемые следом ско рости сильно зависят от положения концевых вихрей, а на режиме малых р след претерпевает значительные дефор мации вблизи диска из-за взаимной индукции вихрей. Таким образом, для повышения надежности расчетов потребуется уточ нить и сделать более эффективными методы определения фор мы следа и индуктивных скоростей (дальнейшее обсуждение способов расчета формы следа и неравномерного распределения индуктивных скоростей см. в гл. 13).
Петерс и Ормистон [Р.55] распространили методы расчета установившегося махового движения на бесшарнирные винты и исследовали влияние различных элементов расчетной схемы на получаемое решение. В результате исследования они сделали следующие выводы относительно выбора расчетной схемы при анализе махового движения и нагрузок лопастей. Для надеж ного расчета n-й гармоники махового движения анализ должен охватывать все гармоники до от-й, где от = п при 0 < ц < 0,4 и от = п + 1 при 0,4 ^ р, <С 1,0. Зону обратного обтекания следует
учитывать |
только при |
р > 0,6, неоперенную |
часть |
лопасти — |
только при р > 1,0, а концевые потери всегда |
важны. Сжимае |
|||
мость воздуха имеет существенное значение, но при A4I>9O<;0,9 |
||||
достаточна простая поправка, получаемая для |
гэфф = |
0,75. При |
||
Mi, до > 0 ,9 |
необходимо |
учитывать изменение |
числа |
Маха по |
радиусу и азимуту. Схема эквивалентной пружины и относа не вполне удовлетворительна при расчете формы изгиба бесшарнирного винта; гораздо предпочтительнее использовать реаль ные формы упругой консольно закрепленной лопасти. Для на дежного расчета нагрузок и движения лопастей нужно учиты
вать |
лишь |
одну форму |
при 0 < р < 0,6, две формы при 0,6 < |
< р |
< 1,2 |
и три формы |
при 1,2 < р < 1,6. Эти выводы приме |
нимы также к шарнирным винтам, так как шарнирно подвешен ную лопасть можно рассматривать как предельный случай кон сольно закрепленной гибкой лопасти.
В качестве дополнительной литературы о работе несущего винта при полете вперед можно рекомендовать: [К-45, М.162,
G.86, |
G.88, |
S.64, |
|
G.90, |
|
W.46 — W.48, |
W.50, W.52 — W.55, |
В.52, |
||||||||||||
S.34, S.169, |
W.56 — W.60, В.35, P.70, |
В.136, |
В.2, В.З, |
Н.112, |
К.49, |
|||||||||||||||
P.85, Н.178, S.119, S.120, |
В.53, P.56, |
S.62, G.123,M.l 14,МЛ 15, А.12, |
||||||||||||||||||
G.133, |
L.79, |
F.5, |
М.168, |
W.104, |
D.51, |
С.21, |
S.170, |
S.181, |
S.192, |
|||||||||||
S.193, |
D.48, |
Р.27 — P.30, |
P.32, P.35, P.37, А.17, Т.28, В.122, С.37, |
|||||||||||||||||
G. 63, |
J.64, |
J.67, |
Н.61, |
М.6, |
М.48, |
G.64, |
R.3, |
Y.16, |
Y.19, |
Н.181, |
J.9, |
|||||||||
P.90, |
А.42, |
М.11, |
|
S.218, |
Е.5, |
J.2, |
N.24, |
Н.46, |
Н.47, |
М.59, |
N.23, |
|||||||||
С.73, |
М.13, |
В. 127, |
P.4, |
|
Y.4, Y.5, L.2, R.69, С.7, С.8, J.18—J.20, |
|||||||||||||||
L.50, |
D.47, |
Н.138, |
L.51, |
1.12, |
В.152, |
V.6, |
V.7, |
D.65, |
Н.177, |
В.85, |
||||||||||
H.82, |
Н.83, |
К.П, S.83, L.23, В.18, С.9, D.66, 1.11, S.82, S.142, |
||||||||||||||||||
S.151, S.161, S.201, Т.57]. Интерференции следа несущего винта |
||||||||||||||||||||
до стенками |
аэродинамической |
трубы |
й поправкам |
на |
влияние |
264 |
Глава 5 |
трубы посвящена следующая литература: [G.16, Н.75, Н.76, Н.79, Н.80, R.10, R.l 1, L.102/L.44, L.101, R. 12]. Экспериментальные ис следования аэродинамических нагрузок несущего винта описы ваются в работах: [М.108, R.l, R.4, М.101, В.166, D.89, Н.17, Н.18, Н.181, S.19,—S.21, S.15, Н.28, D.94, R.5, R.6, R.8, P.93, Р.96, В.29, В.94, F.18, В.116, Т.68, W.13, В.56, D.4, R.59, V.4, G.106, S.96, S.95, К.52, М. 151, S.87].
6
Аэродинамический расчет вертолета
Аэродинамический расчет вертолета сводится в основном к определению потребной и располагаемой мощностей в рассмат риваемом диапазоне режимов полета. Данные о мощности мо гут быть затем преобразованы в такие величины, как скоро подъемность, потолок, дальность и максимальная скорость, ко торые определяют летно-технические характеристики вертолета. Потребную мощность можно представить суммой четырех час тей: 1) индуктивной мощности, затрачиваемой на создание силы тяги винта, 2) профильной мощности, необходимой для враще ния винта в воздухе, 3) затрат мощности на преодоление вред ного сопротивления, т. е. на продвижение вертолета в воздухе, и 4) затрат мощности на набор высоты, т. е. на изменение по тенциальной энергии вертолета. На режиме висения для преодо ления вредного сопротивления мощность не затрачивается, а ин дуктивная мощность составляет 60 4- 70% общих затрат. С уве личением скорости полета индуктивная мощность уменьшается, профильная слегка возрастает, а мощность, затрачиваемая на вредное сопротивление, увеличивается вплоть до того, что ста новится доминирующей при больших скоростях. Таким образом, потребная мощность велика на висении вследствие больших ин дуктивных затрат при приемлемой нагрузке на диск (хотя винт и малонагруженный), далее она сначала уменьшается с ростом скорости полета в результате уменьшения индуктивной мощ ности, а затем снова увеличивается, так как при больших ско ростях велика мощность, затрачиваемая на преодоление вред ного сопротивления. Потребная мощность минимальна прибли зительно в середине диапазона скоростей вертолета.
Главной задачей в анализе характеристик вертолета являет ся расчет нагрузок и мощности несущего винта. Методы такого расчета изложены в предшествующих главах. Существует два основных подхода к расчету аэродинамических характеристик несущего винта: метод тяг и метод мощностей. При использова нии первого метода интегрируют элементарные силы, действую щие в сечениях лопастей, и получают результирующие силы и аэродинамический крутящий момент несущего винта. Для этого нужно знать индуктивные скорости и движение лопастей, по ко торым находят распределение углов атаки. Затем из условий равновесия сил и моментов определяют балансировочные углы,
266 |
Глава 6 |
а по ним — силы |
и положение вертолета в пространстве, кото |
рые требуются для сохранения заданного режима полета. Дру гой способ состоит в том, что аэродинамические характеристики винта вычисляют для некоторого диапазона сил тяги и поло жений вертолета в пространстве, а затем строят серии графи ков характеристик. В гл. 5 были выведены выражения силы тяги и крутящего момента винта через элементарные силы лопасти. Метод тяг, даже в его простейшем варианте, сложен, и потому он лучше подходит для численного решения. Он также весьма удобен при использовании наиболее усовершенствованных рас четных схем обтекания винта. С помощью метода тяг можно также построить только графики аэродинамических характери стик несущего винта, а затем использовать их для определения характеристик всего вертолета.
Второй метод расчета характеристик вертолета состоит в том, что потребную мощность выражают через определяемые по от дельности затраты энергии на вертолете. В гл. 5 условие ба ланса энергии было получено из условия равновесия сил и тем самым показано, что оба метода эквивалентны: получаемые ими результаты совпадают, если совпадают исходные предположе ния. По нескольким причинам метод мощностей удобнее для выполнения стандартных расчетов характеристик. Во-первых, равновесие сил, действующих на вертолет в продольной пло скости, уже было рассмотрено, так что мощность можно нахо дить сразу, без необходимости определять - балансировочные углы. Во-вторых, мощности, затрачиваемые на преодоление вред ного сопротивления и на набор высоты, вычисляются по про стым и в то же время точным формулам. Индуктивную же и профильную мощности можно определять отдельно, и примене ние соответствующих приближенных выражений не вызывает затруднений. Если использовать простейшие приближенные вы ражения, то метод мощностей позволяет рассчитать характера; стики быстро и с. приемлемой погрешностью, вследствие чего он очень удобен для расчетов на предварительной стадии проекти рования. Для более обстоятельного анализа характеристик нужны уточненные формулы индуктивной и профильной мощ ностей, применение которых снова потребует расчета распре деления углов атаки. Таким образом, численные методы тяг и мощностей даже с вычислительной точки зрения эквивалентны, хотя разделение всей требуемой мощности на индуктивную, про фильную, мощность на преодоление вредного сопротивления и мощность на набор высоты полезно и при численном решении для интерпретации результатов.
Итак, аэродинамический расчет вертолета выполняют, во обще говоря, тремя способами: 1) методом мощностей с исполь зованием весьма простых выражений для индуктивной и про фильной мощностей, 2) с помощью графиков характеристик,
Аэродинамический расчет вертолета |
2 6 7 |
построенных на базе аналитических или численных решений, получаемых обычно методом тяг, и 3) численно с использова нием расчетной схемы, выбираемой в зависимости от целей ана лиза и располагаемых возможностей. Численный анализ с ис пользованием ЦВМ повсеместно вошел в практику расчетов на заключительной стадии проектирования. В настоящей главе при ведена сводка методов аэродинамического расчета вертолета на режимах висения и полета вперед. Подробные выводы формул и описания методов расчета даны в различных местах, в част ности в гл. 2 и 5, а в гл. 14 изложены численные методы. Затем рассмотрены величины, описывающие летно-технические харак теристики вертолета, в том числе потребную мощность, макси мальную скорость, скорость набора высоты и спуска, потолок, дальность и продолжительность полета. Требование, чтобы вер толет имел определенные характеристики, в широком смысле означает способность наиболее эффективным образом выпол нить конкретную задачу, а возможность рассчитать эти харак теристики нужна не только для определения эксплуатационных возможностей вертолета, но и для исследования круга задач, которые он может выполнить. Глава завершается обсуждением описанных в литературе методов расчета характеристик и по лученных результатов.
6.1. АЭРОДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ В ВЕРТИКАЛЬНОМ ПОЛЕТЕ
6.1.1. ПОТРЕБНАЯ МОЩНОСТЬ НА РЕЖИМАХ ВИСЕНИЯ И ВЕРТИКАЛЬНОГО ПОЛЕТА
Коэффициент мощности, потребной для вертикального по
лета, был определен в виде СР — СР. + Ср0 + |
Срс, где |
|
R |
1 |
|
CPl = ^ Я(- dCf, |
СРо = ^ (огс^0/2) г3 dr, |
Ср0 = ксСр |
Га |
О |
|
соответственно коэффициенты индуктивной и профильной мощ ностей, а также мощности, затрачиваемой на набор высоты. В размерной форме это соотношение будет следующим:
BR
P = Pt + Po + Pc= \ (V + v)dT + pA(QRfCPo. Го
Индуктивная мощность Р ,— это рассеиваемая в единицу вре мени в следе несущего винта энергия, которую винт сообщает воздуху, отбрасывая его вниз. Реакция воздуха является подъем ной силой винта. В гл. 2 и 3 с использованием импульсной тео рии были получены простейшие формулы для индуктивной мощ*
268 Глава 6
ности. Путем введения эмпирических поправок в формулы для индуктивной скорости на тех режимах обтекания, к которым импульсная теория неприменима, для суммы индуктивной мощ ности и мощности, затрачиваемой на набор высоты в вер тикальном полете, был получен универсальный график зависи мости P /P &= ( V -\-v)/vв от V/vB, где vl = T/(2pA) (рис. 3.8). На
висении |
индуктивная |
скорость определяется соотношением |
v = kvB= |
k д/77(2рЛ), |
где k — эмпирический коэффициент, учи |
тывающий дополнительные потери мощности (в основном кон цевые потери и потери, обусловленные неравномерностью про текания); в типичных случаях k = 1,1 ч- 1,2 (см. разд. 3.1.3.1). Для вертикального полета несколько лучшую оценку индуктив ной мощности можно получить по элементно-импульсной тео рии (см. разд. 2.5).
Профильная мощность Р0— это энергия, рассеиваемая в еди ницу времени вследствие сопротивления лопастей при их дви жении в вязком воздухе. Грубая оценка коэффициента Ср0 была получена при использовании среднего по лопасти коэффициента сопротивления: CPo — acdj8. Для более точной оценки этого
коэффициента нужно интегрировать элементарные сопротивле ния по радиусу лопасти, используя реальные распределения уг лов атаки и чисел Маха по лопасти. Такие расчеты могут быть выполнены на базе элементно-импульсной теории.
Характеристики несущего винта на режиме висения можно представить в виде поляры винта — графика зависимости СР от Ст■Для построения поляры обычно используют приближенную
формулу |
_ |
|
СР= |
kCr2/ V 2 |
+ ocdjs |
или |
|
|
Р = kT V Т/2рА + 9A (QRf (<rcdJ8)-
Для построения поляры винта в вертикальном полете исполь зуется формула Ср = (А,г + Я,с) Ст+ <rcdo/8, причем при малых
скоростях набора высоты или снижения Я,/ + Хс » Яв + Яс/2. Теория элемента лопасти позволяет найти общий шаг. Для
лопастей с линейной круткой и постоянной хордой при равно мерном протекании справедливо соотношение
00,75 == бСу/сго 4 “ (3/2) К.
При численном определении характеристик винта общий шаг фактически является параметром, по которому вычисляют Ср и Ст, получая затем поляру.
Эффективность работы несущего винта на висении измеряют коэффициентом совершенства
м = г Л/ Т Ш / Р = сГ Ы Ъ С р,
Аэродинамический расчет вертолета |
269 |
который выражает отношение индуктивной и профильной мощ ностей. В типичных случаях профильная мощность составляет приблизительно 30% полной, а реальная индуктивная мощ ность— около 10%, так что М ~ 0,6.
6.1.2.НАБОР ВЫСОТЫ И СНИЖЕНИЕ
Вразд. 3.3 была получена следующая формула для скорости набора высоты по вертикали при заданном избытке мощности:
V = (АР К ) (2ив + AP/T)/{Vv + АР/Т).
Здесь по-прежнему и^ = 7'/2рЛ, а А Р — избыток мощности по
сравнению с той, которая нужна для висения. Так как эта фор мула получена по импульсной теории, она применима и при ма лых скоростях снижения. Если скорость набора высоты или снижения мала, то формула упрощается: V ~ 2&Р/Т. Таким об разом, уменьшение индуктивной скорости вследствие увеличе ния потока сквозь винт при наборе высоты удваивает эффек тивность использования заданного избытка мощности.
В гл. 3 были рассмотрены способы расчета скорости верти кального спуска при авторотации винта на режиме турбулент ного следа с помощью универсальной кривой индуктивной мощ
ности. Авторотация |
реального винта определяется условием |
р = 7(1/ + ц )+ Р о = |
0 или (V + v)/и в = —Ро/Рв, откуда и на |
ходят скорость V/vB спуска. В разд. 3.2 была получена прибли женная формула, основанная на том, что на режиме турбулент ного следа кривая индуктивной мощности аппроксимируется прямой линией, а профильная мощность считается такой же, как на висении. Эта формула имеет вид
VlvB~ - [1,71 + 0,29(l/A f-/e)].
6.1.3. РАСПОЛАГАЕМАЯ МОЩНОСТЬ
Располагаемая мощность определяется характеристиками силовой установки. Обычно мощность силовой установки сни жается с ростом высоты и температуры, а также в какой-то мере зависит от скорости полета. Поэтому при расчете характеристик вертолета важны изменения располагаемой мощности. Следует также учитывать потери мощности в силовой установке и в трансмиссии, включая потери в редукторе и в системе охлаж дения двигателя, а также мощность, затрачиваемую на привод вспомогательных агрегатов, таких, как генераторы и насосы гид росистемы. Часто все эти потери выражают посредством общего к. п. д. г), т. е. считают, что общая потребная мощность в (1/ri) раз больше той, которая требуется для вращения винта:
1
270 |
Глава 6 |
При типичных потерях |
в двигателе и системе привода т] = |
= 0,91-г 0,96.
Кроме затрат мощности на отдельный несущий винт имеются еще дополнительные потери. Потери на аэродинамическую ин терференцию несущих винтов и винта с фюзеляжем составляют значительную часть располагаемой мощности, особенно у вер толетов продольной схемы. У вертолетов одновинтовой схемы нужно учитывать также потери на рулевой винт. Расчет харак теристик рулевого винта осложнен тем, что этот винт работает в следе несущего винта и фюзеляжа. Интерференция уменьшает эффективность рулевого винта; особенно увеличиваются его на грузки и вибрации. При маневрировании по рысканию рулевой винт может даже попасть в режим вихревого кольца, вследствие чего ухудшается управление и значительно усиливаются вибра ции. Характеристики рулевого винта можно рассчитать, учиты вая, что его сила тяги задана аэродинамическим моментом не сущего винта, т. е. Тр. „ = Q/lP. в, где 1Р. „ — плечо рулевого винта относительно вала несущего винта. Так как потребная мощность рулевого винта составляет малую часть общей мощности, а по тери на интерференцию нужно как-то оценить, часто прибе гают к весьма приближенным формулам. Потери на интерферен цию между частями вертолета и потери на рулевой винт можно также учесть в общем к.п.д. rj. При этом нужно рассчитать только затраты мощности на несущий винт, а полная потребная мощность определяется умножением этих з'атрат на коэффи циент 1/т]. Если принять в расчет потери в силовой установке и в трансмиссии, а также потери на интерференцию и рулевой винт, то на режиме висения в типичном случае ri составляет 0,80 -н 0,87. При полете вперед ц, как правило, больше, по скольку потери на интерференцию и на рулевой винт умень шаются.
6.2. АЭРОДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ПРИ ПОЛЕТЕ ВПЕРЕД
6.2.1.ПОТРЕБНАЯ МОЩНОСТЬ ПРИ ПОЛЕТЕ ВПЕРЕД
Вгл. 5 коэффициент аэродинамического крутящего момента
несущего |
винта был |
получен в виде интеграла от лежащих |
в плоскости диска составляющих элементарных сил лопасти: |
||
|
|
1 |
|
СР= CQ = |
^ (1/2) arU2(ci sin cp + cd cos <p) dr, |
|
|
6 |
где C/2 = |
Mj + Hp и cp = arctg (мр/«г). Слагаемое, содержащее Ci, |
представляет собой ускоряющий момент, а слагаемое, содержа щее Са, — замедляющий. Отсюда было выведено уравнение 6а«
Аэродинамический расчет вертолета |
271 |
ланса мощностей при полете вертолета вперед
Ср = СР. + Сро + сРвр + г .р с,
где коэффициенты индуктивной и профильной мощностей, а также мощностей, затрачиваемых на преодоление вредного со противления и набор высоты, вычисляются соответственно по формулам
|
R |
|
|
CPi |
$ Яг dCT, |
Срвр = |
DV/[pA (Я/?)3], |
|
Го |
|
|
Cp0 = |
4 + ^c v |
Cpc = |
v cw / [ p A w n |
Напомним, что при выводе этого уравнения предположение о малости углов не использовано (см. разд. 5.4 и 5.18). При по лете вперед в баланс мощностей входит мощность PBP^=DV, затрачиваемая на преодоление сопротивления D вертолета, дви жущегося в воздухе.
В гл. 4 было выведено уравнение баланса мощностей Р = = Т( V sin а + v) = Pi -f РВр + Рс для случая полета вперед вер толета с идеальным винтом (без профильных потерь).На рис.4.4 представлено решение этого уравнения, полученное по элемент но-импульсной теории и экспериментальным данным в виде за висимости Р/Рв = (V'sin.a + v) /цв от V cosa/yBи V sin a /v B. По импульсной теории индуктивный коэффициент протекания при полете вперед равен
|
|
Я, = Сг/2 V IA2 + Я2, |
|
||
где Я = |
Я |
, p. tg а. Для |
любых |
скоростей полета, |
кроме самьце |
малых, |
выражение Я, |
хорошо |
аппроксимируется формулой |
||
%{~Ст/2ц |
(см. разд. 4.1.1). Эта формула очень |
полезна, тан |
как в ней не фигурирует скорость набора высоты или снижения. Если еще ввести эмпирическую поправку k, то для коэффициента индуктивной мощности при полете вперед получается выраже ние Cpt = XICT — £C|/2p, т. е. Pi = kV'I'lpA V.
Коэффициент профильной мощности несущего винта был оп ределен в разд. 5.12 в виде
Ср«= SKfa/2)(4 +
о
Здесь учтено влияние зоны обратного обтекания, радиального течения и радиальной силы сопротивления. Если использовать средний коэффициент сопротивления сечений, то приходим к приближенной формуле
Cpa’±=acdAl + 4 ,6и2)/81