КАТЕГОРИИ: Архитектура-(3434)Астрономия-(809)Биология-(7483)Биотехнологии-(1457)Военное дело-(14632)Высокие технологии-(1363)География-(913)Геология-(1438)Государство-(451)Демография-(1065)Дом-(47672)Журналистика и СМИ-(912)Изобретательство-(14524)Иностранные языки-(4268)Информатика-(17799)Искусство-(1338)История-(13644)Компьютеры-(11121)Косметика-(55)Кулинария-(373)Культура-(8427)Лингвистика-(374)Литература-(1642)Маркетинг-(23702)Математика-(16968)Машиностроение-(1700)Медицина-(12668)Менеджмент-(24684)Механика-(15423)Науковедение-(506)Образование-(11852)Охрана труда-(3308)Педагогика-(5571)Полиграфия-(1312)Политика-(7869)Право-(5454)Приборостроение-(1369)Программирование-(2801)Производство-(97182)Промышленность-(8706)Психология-(18388)Религия-(3217)Связь-(10668)Сельское хозяйство-(299)Социология-(6455)Спорт-(42831)Строительство-(4793)Торговля-(5050)Транспорт-(2929)Туризм-(1568)Физика-(3942)Философия-(17015)Финансы-(26596)Химия-(22929)Экология-(12095)Экономика-(9961)Электроника-(8441)Электротехника-(4623)Энергетика-(12629)Юриспруденция-(1492)Ядерная техника-(1748) |
Ol — V2X 56 страница
§ 92. Косая ударная волна Рассмотрим стационарную ударную волну, отказавшись при этом от подразумевавшегося везде выше выбора системы координат, в которой скорость газа направлена перпендикулярно к данному элементу поверхности волны. Линии тока могут пересекать поверхность такой ударной волны наклонно, причем пересечение сопровождается преломлением линий тока. Касательная составляющая скорости газа не меняется при прохождении через ударную волну, а нормальная составляющая согласно (87,4) падает: Vu = V2t, f l/i > V2n. Поэтому ясно, что при прохождении через ударную волну линии тока приближаются к ней (как это показано на рис. 63). Таким образом, преломление линий тока на ударной волне происходит всегда в определенном направлении. Выберем направление скорости vi газа перед ударной волной в качестве оси х, и пусть ср — угол между поверхностью разрыва и осью х (рис. 63). Возможные значения угла ср ограничены лишь условием, чтобы нормальная составляющая скорости vi ') Выражение (91,7) для произвольной (не политропной) среды и его связь с условиями гофрировочной неустойчивости ударных волн указаны С. Г. Су гаком и В. Е. Фортовым, А. Л, Ни (1981). превышала скорость звука с\. Поскольку Vi„ = v sin ф, то отсюда следует, что ф может иметь произвольные значения в интервале между я/2 и углом Маха оц: а, < ф < л/2, sin cci = c\/v\ = 1/Mi.
Движение позади ударной волны может быть как до-, так и сверхзвуковым (меньше скорости звука с2 должна быть лишь нормальная компонента скорости); движение же перед ударной волной — непременно сверхзвуковое. Если движение газа по обе стороны от ударной волны является сверхзвуковым, то все возмущения могут распространяться вдоль ее поверхности лишь в ту сторону, ку-х да направлена касательная к ней составляющая скорости газа. В этом смысле можно говорить о «направлении» ударной волны и различать по отношению к какому-либо месту «исходящие» из него и «приходящие» волны (подобно тому как мы это уже делали для характеристик, вокруг которых движение всегда является сверхзвуковым; см. § 82). Если же движение позади ударной волны является дозвуковым, то понятие о ее направлении теряет, строго говоря, смысл, так как возмущения могут распространяться вдоль ее поверхности во все стороны. Выведем соотношение, связывающее друг с другом две компоненты скорости газа после его прохождения через косую ударную волну; при этом будем предполагать газ политропным. Непрерывность касательной к волне составляющей скорости означает, что v\ cos ф = v2x cos ф + viu sin ф, или tgop-. (92,1)
Далее, воспользуемся формулой (89,6); в этой формуле Vi и i>2 обозначают нормальные к плоскости ударной волны составляющие скорости и должны быть теперь заменены на vi sin ф и v%x sin ф — V2y cos ф, так что имеем: уц gjnФ — v2y cos.cp _ у — 1 "-2 2с1
(92,2) t), sin ф Из двух написанных соотношений можно исключить угол ф. После простых преобразований получим следующую формулу, определяющую связь между v2x и v2y (при заданных v\ и С\)
2—G- -^)-^-v Vly = (Vl - V2xf Y+ 1 0| Этой формуле можно придать более изящный вид, если ввести в нее критическую скорость. Согласно уравнению Бернулли и определению критической скорости имеем: 2 v (ср. задачу 1 § 89), откуда Т + с:
■v2 + - Y+ 1 1 Y+ 1 Введя эту величину в (92,3), получим: (92,4)
Чу = (»1 - 4xf
Y-+ 1 viv2X -«; ■»1 - vl°ix + с1
(92,5) Уравнение (92,5) называют уравнением ударной поляры {A. Busemann, 1931). На рис. 64 изображен график этой зависимости; это есть кривая третьего порядка (так называемая \ строфоида или декартов лист). Она пересекает ось абсцисс в точках Р и Q (рис. 64), соответствующих значениям v2x = c;/rj и v2x = vii). Проведя из начала о координат луч (ОВ на рис. 64) под углом х к оси абсцисс по длине его отрезка до точки пересечения с кривой ударной поляры, мы определяем скорость газа за скачком, поворачивающим поток на угол %. Таких точек пересечения имеется две (А и В), т. е. заданному значению х отвечают две различные ударные волны. Направление ударной волны тоже может быть ') От точки Q, являющейся двойной точкой кривой, строфоида в действительности продолжается еще. в виде двух уходящих к бесконечным \v2y\ ветвей (не изображенных на рис. 64) с общей вертикальной асимптотой Ds* —«?M+2«»,/(Y + 1) Однако точки этих ветвей не имеют физическою смысла: они дали бы для v2x, vzy значения, при которых v2nlv\n > 1, что невозможно. сразу определено графически по этой же диаграмме — оно определяется перпендикуляром, опущенным из начала координат на прямую, проведенную из точки Q соответственно через точку В или А (на рис. 64 изображен угол ср для волны, соответствующей точке В). При уменьшении % точка А приближается к точке Р, отвечающей прямому (ср = л/2) скачку с v2 = cl/vi- Точка же В приближается при этом к точке Q, причем интенсивность ударной волны (скачок скорости в ней) стремится к нулю; в пределе, в самой точке Q, угол ср равен, как и следовало, углу Маха ат (угол наклона касательной к поляре к оси абсцисс в этой точке равен я/2 + ai). Из диаграммы ударной поляры сразу можно вывести важное заключение, что угол отклонения % потока в ударной волне не может превышать некоторого максимального значения %тах, соответствующего луч>, проведенному из точки О касательно к кривой, хтах является, конечно, функцией числа Mi = v\/c\\ мы не приводим ее здесь ввиду ее громоздкости. При Mi = 1 имеем Хтах = 0, а при возрастании Mi угол хтах монотонно растет и при Mi->-сю стремится к конечному пределу. Легко рассмотреть оба предельных случая. Если скорость V\ близка к с», то вместе с ней близка к с» и скорость v2, а угол % мал; уравнение ударной поляры (92,5) можно тогда приближенно переписать в виде') X2 = JT3i (у 1 - у2)2 (»1 + Щ - 2с.) (92,6> (ввиду малости угла % здесь положено v2x та v2, v2y ж с*%). Отсюда элементарным путем найдем2): = lVY±I^ _ Л3/2= 2?/2 (М, _ П3/2 m 7у В обратном предельном случае, при Mi -*- со, ударная поляра вырождается в окружность
Легко видеть, что при этом Хтах = arcsi-n(l/y). (92,8)
На рис. 65 изображен график зависимости Хтах от Mi для воздуха (у = 1,4); горизонтальный пунктирный отрезок показывает предельное значение %тах(°°) = 45,6° (верхняя кривая на рисунке— аналогичный график для обтекания конуса; см. § 113).
Окружность V2 — с* пересекает ось абсцисс между точками Р и Q (рис. 64) и поэтому делит ударную поляру на две,части, соответствующие до- и сверхзвуковым скоростям газа позади разрыва. Точка пересечения окружности V2 = с* с полярой лежит пра- во* вее точки С, но очень близко к ней; поэтому весь участок РС соответствует переходам к дозвуковым скоростям, а участок CQ (за исключением лишь очень небольшого участка вблизи точки С) — переходам к сверхзвуковым скоростям. /, Изменения давления и плотно-
и Ф получаются из формул (89,6—7) просто путем замены в них Mi на Mi sin ф:
2(М2 sin2 ф — l)] (Y — 1)М2 sin* ф-г-2 Эти отношения монотонно возрастают при увеличении угла ф от значения ф = а\ (когда р2/р\ == Рг/pi = 1) до я/2, т. е. по мере перемещения по ударной поляре от точки Q к точке Р. Приведем еще, для справок, формулу, выражающую угол поворота % скорости через число Mi и угол ф: ctg % = tg ф £ (Y+l)Mf 2(М28ш2ф- l) (92,11) и формулу, определяющую число М2 = v2/c2 по Mi и ф|
2 + (Y-I)Mf (Y-1) 2 + (Y-l)M2sin2? (92,12) (при ф = я/2 последнее выражение переходит в (89,9)). Две ударные волны, определяемые ударной полярой для заданного угла поворота скорости, называют волнами слабого и сильного семейства. Ударная волна сильного семейства (участок РС поляры) обладает большей интенсивностью (большим отношением p2/pi), образует больший угол ф с направлением скорости vi и превращает течение из сверх- в дозвуковое. Волна же слабого семейства (участок QC поляры) обладает § 93. Ширина ударных волн Мы говорили до сих пор об ударных волнах как о геометрических поверхностях, не обладающих толщиной. Рассмотрим теперь вопрос о структуре реальных физических поверхностей разрыва. Мы увидим, что ударные волны с небольшими скачками величин представляют собой в действительности переходные слои конечной толщины, уменьшающейся при увеличении величины скачков. Если же скачки величин в ударной волне не малы, то, действительно, разрыв происходит настолько резко, что в макроскопической теории не имеет смысла говорить о его толщине. Для определения структуры и толщины переходного слоя надо учесть вязкость и теплопроводность газа, влиянием которых мы до сих пор пренебрегали. Соотношения (85,1—3) на ударной волне были получены из условий постоянства потоков вещества, импульса и энергии. Если рассматривать поверхность разрыва как слой конечной толщины, то эти условия надо писать не в виде равенства соответствующих величин по обе стороны разрыва, а в виде их постоянства вдоль всей толщины разрывного слоя. Первое из этих условий (85,1) не меняется: pv =е / = const. (93,1) В двух же других условиях надо учесть дополнительные потоки импульса и энергии, обусловленные внутренним трением и теплопроводностью. Плотность потока импульса (вдоль оси х), обусловленного внутренним трением, определяется компонентой — а'хх вязкого тензора напряжений; согласно общему выражению (15,3) для этого тензора имеем:
Условие (85,2) приобретает теперь вид1) p + pf2-(4n + S)-57 = const. Как и в § 85, введем вместо скорости v удельный объем V согласно v = jV. Постоянную же в правой стороне равенства выразим через предельные значения величин на большом расстоянии впереди ударной волны (сторона /). Тогда написанное ') Положительное направление оси к совпадает с направлением движения газа через неподвижную ударную волну. Если перейти к системе отсчета, в которой неподвижен газ перед ударной волной, то сама ударная волна будет двигаться в отрицательном направлении оси х. условие примет вид Р - Pi + j2 (V - V,) - (-i г, + С) j = о. (93,2) Далее, плотность потока энергии, обусловленного теплопроводностью, есть —кдТ/дх. Поток же энергии, связанный с внутренним трением, есть / (4. „\ dv - °Xivt = - °xxv--{T^ + 4v-dx-- Таким образом, условие (85,3) напишется в виде (I v2 \ (4 I Л dv dT i
или, снова введя v = jV и выразив const через величины с индексом 1: w-Wi + ^(v2-vb-i(U+^^-i^^- 03.3) Мы будем рассматривать здесь ударные волны, в которых все величины испытывают лишь малый скачок. Тогда и все разности V—V\, р — Pi ит. п. между значениями величин внутри переходного слоя и вне его тоже малы. Из получающихся ниже соотношений видно, что 1/6 (где б — ширина разрыва) есть величина первого порядка малости по р2 — р\. Поэтому дифференцирование по х увеличивает порядок малости на единицу (так, производная dp/dx — величина второго порядка). Умножим уравнение (93,2) на (V+Vi)/2 и вычтем его из уравнения (93,3). Тогда получим: (^-^,)-l(p-p,)(K + Ki) = yg- (93,4) (здесь опущен член, содержащий (V — Vi)dV/dx, являющийся малой величиной третьего порядка). Разложим выражение в левой стороне (93,4) по степеням р — рх и s — s\, выбрав давление и энтропию в качестве основных независимых переменных. Члены первого и второго порядка по р — рх в этом разложении выпадают (ср. вычисления при выводе формулы (86,1)) и, опустив члены более высокого порядка, получим просто T(s — s\). Производную же dT'/dx пишем в виде dx \др)s dx ' V. ds Jp dx Член с производной ds/dx можно опустить как малую величину третьего порядка (см. ниже), и в результате находим формулу, выражающую функцию s(x) через функцию р(х): 'с—>-t(v).*- (9М) Обратим внимание на то, что разность s — Si внутри переходного слоя оказывается величиной второго порядка малости, между тем как полный скачок s2 — si является (как было показано в § 86) величиной третьего порядка по сравнению со скачком давления р2— рь Это связано с тем, что (как будет показано ниже) давление р(х) меняется в переходном слое монотонно от одного предельного значения pi до другого р2; энтропия же s(x), определяясь производной dp/dx, проходит через максимум, достигая наибольшего значения внутри переходного слоя. Уравнение, определяющее функцию р(х), можно было бы получить путем аналогичного разложения уравнений (93,2—3) и их комбинирования друг с другом. Мы, однако, изберем другой, более поучительный способ, позволяющий более ясно понять происхождение различных членов в уравнении. В § 79 было показано, что монохроматическое слабое возмущение состояния газа (звуковая волна) затухает по мере своего распространения с декрементом, пропорциональным квадрату частоты: у = асо2; положительный коэффициент а выражается через коэффициенты вязкости и теплопроводности согласно формуле (79,6). Там же было показано, что это затухание может быть описано (для произвольной плоской звуковой волны) введением дополнительного члена в линеаризованное уравнение движения — см. (79,9). Заменив в этом уравнении вторую производную по времени второй производной по координате и изменив знак перед производной др'/дх (что отвечает распространению волны в отрицательном направлении оси х1)), запишем его в виде с дх —ас дх2 • ^'Ь> где р' — переменная часть давления. Для учета слабой нелинейности надо ввести в это уравнение член вида р'др'/дх:
Коэффициент <хр в нелинейном члене определяется путем соответствующего разложения гидродинамических уравнений идеальной (без диссипации) жидкости и оказывается равным
(см. задачу) 2). Уравнение (93,7) описывает распространение возмущений в слабо диссипирующей, слабо нелинейной среде. В применении к слабой ударной волне оно описывает ее распространение в системе отсчета, в которой невозмущенный газ (перед волной) неподвижен. Требуется найти решение со стационарным (не зависящим от времени) профилем, в котором вдали от волны, при #->-±оо, давление принимает заданные значения р2 и р\\ разность р2 — р\ есть скачок давления в разрыве1). Волна со стационарным профилем описывается решением вида p'(x,t)^p'(x + Vlt), ■ (93,9) где v\ — скорость распространения такой волны. Подстановка в (93,7) приводит к уравнению ± [(vx - с) р' - % р'2 _ ас3 = 0, % = х + vd, первый интеграл которого: acS-^ = -^p'2 + (Vi — с)//+ const. (93,10) Квадратный трехчлен в правой стороне равенства должен обращаться в ноль при значениях р', отвечающих предельным условиям на бесконечностях, где производная dp'/d% обращается в ноль. Эти значения равны р2— рх и 0 если условиться отсчитывать р' от невозмущенного давления р\ перед волной. Это значит, что указанный трехчлен может быть представлен в виде ~^\р' ~(р2-рхЛр', причем константа v\ выражается через р\ и р2 согласно oi = с + -^ (ft-/>.)■ (93,11) Для самого же давления р уравнение (93,10) принимает вид ас -ж = — -т(р — Р\)(р -PJ-
(93,7) к виду ди, ди дги ■ЬТ + и-Щ = *-аё' (93,7а) в котором его называют уравнением Бюргерса (J. М. Burgers, 1940). ■) Мы увидим в дальнейшем (§ 102), что в отсутствии диссипации эффекты нелинейности приводят к искажению профиля волны по мере ее распространения — постепенному возрастанию крутизны фронта волны. В свою очередь, это возрастание приводит к усилению диссипативных эффектов, стремящихся уменьшить крутизну профиля (т. е. уменьшить градиенты меняющихся величин). Именно взаимная компенсация этих противоположных тенденций приводит к возможности распространения волн со стационарным профилем в нелинейной диссипативиой среде. Решение этого уравнения, удовлетворяющее требуемым условиям есть _ Pl + Р2, Р2 — Pl.. (Р2 — Pl)(*+PlQ V 2 ~Г 2 1 4ас3/ар Этим решается поставленная задача. Возвратившись снова к системе отсчета, в которой ударная волна покоится, напишем формулу, определяющую ход изменения давления в ней в виде р_ £l±EL= P-ZZlth-J. {93,12) где 8aV2 б== (Р2-Р1) (d2V/dp%- (93,13) Практически все изменение давления от р\ до р2 происходит на расстоянии ~б— ширине ударной волны. Мы видим, что ширина волны уменьшается с увеличением ее интенсивности — скачка давления р2— Р\ '). Для хода изменения энтропии внутри разрыва имеем из (93,5) и (93,12): *~s'=l&aVT(4R(-0"),(Р*~Р>>2 Отсюда видно, что энтропия меняется не монотонно, а имеет максимум внутри ударной волны (при х — 0). При л: = ±оо эта формула дает одинаковые значения s = sr. это связано с тем, что полное изменение энтропии s2 — S\ является величиной третьего порядка по р2 — pi (ср. (86,1)), в то время как s — si — второго. Формула (93,12) применима количественно только при достаточно малых разностях р2— Р\. Однако качественно мы можем применить формулу (93,13) для определения порядка величины ширины ударной волны и в тех случаях, когда разность р2 — Р\ — порядка величины самих давлений ри р2. Скорость звука в гггг — порядка величины тепловой скорости v молекул. Кинематическая же вязкость, как известно из кинетической теории газов, v ~ lv ~ 1с, где / — длина свободного пробега молекул. Поэтому а ~ 1/с2 (оценка члена с теплопроводностью дает то же самое). Наконец, (d2V/dp2)s ~ V/p2 и pV ~ с2. Внося эти выражения в (93,13), получаем: б ~ I. (93,15) ') Для ударной волны, распространяющейся в смеси, определенный вклад в ее ширину возникает также и от процессов диффузии в переходном слое. Вычисление этого вклада см. Дьяков С. П. — ЖЭТФ, 1954, т. 27, с. 283 Упомянем также, что ударные волны слабой интенсивности остаются устойчивыми по отношению к поперечной модуляции (ср. примечание на стр. 477) и ппи учете их диссипативной структуры; см. Спектор М. Д. — Письма ЖЭТФ, 1983, т. 35, с. 181. Таким образом, ширина ударных волн большой интенсивности оказывается порядка величины длины свободного пробега молекул газа1). Но в макроскопической газодинамике, трактующей газ как сплошную среду, длина свободного пробега должна рассматриваться как равная нулю. Поэтому, строго говоря, чисто газодинамические методы непригодны для исследования внутренней структуры ударных волн большой интенсивности. Задачи 1. Определить коэффициент нелинейности а в уравнении (93,7) для распространения звуковых волн в газе. Решение. Точные гидродинамические уравнения одномерного движения идеального (без диссипации) газа: dv dv 1 dp др д... dt дх p dx dt dx ' Произведем их разложение с учетом членов второго порядка малости. Для этого полагаем р = Ро+Р', p==po+ ^l+^l(0)s. <2) Члены второго порядка в уравнениях можно упростить, приведя их всех к одинаковому виду — содержащему произведение р'др'/dx. Для этого замечаем, что для волны, распространяющейся в отрицательном направлении оси х (со скоростью с) дифференцирование по t эквивалентно дифференцированию по х/с; при этом v — —р'/сро- После всех этих замен получим из (1) и (2) следующие уравнения: dt р дх w dv, 1 dp' (dW\, dp' ~dx^~pc^~dr = Cp\dp^)sP dx <4> (индекс 0 у постоянных равновесных значений величин опускаем); здесь использовано также равенство
(V = 1/р —удельный объем). Дифференцируя уравнения (3) и (5) соответственно по х и по t и вычтя одно из другого, получим 421/ \ Я, Яр' Ч (\ d д \(\ д. д \, 2 г(d2V\ d (, дх С той же точностью заменяем в левой стороне этого уравнения д/дх + + d/cdt -*■ 2д/дх. Наконец, вычеркнув с обоих сторон дифференцирования по х и сравнив получившееся уравнение с (93,7), найдем для ар значение (93,8). Уравнение для скорости v можно получить непосредственно из (93,7), не повторяя заново вычислений, подобных произведенным выше. Действительно, сумма членов первого порядка в левой стороне (93,7) содержит оператор [) Сильная ударная волна сопровождается значительным увеличением температуры; под / надо понимать длину пробега, соответствующую некоторой средней температуре газа в волне. d/dt— сд/дх, который надо рассматривать как малый первого порядка: он обращает в ноль функцию р'(х, t) в ее линейном приближении. Поэтому мы получим уравнение для функции v (х, t) в требуемом приближении, просто заменив в (93,7) р' согласно линейному соотношению р' = —реи:
dv dv, dv, d2v где \ dp2)s v 2V3 \ dp Величина a„ безразмерна; для политропного газа av = (у + 1)/2. 2. Путем нелинейной подстановки привести уравнение Бюргерса (93,7а) к виду линейного уравнения теплопроводности (Е. Hopj, 1950). Решение. Подстановкой в(С,0 = -2ц-|^-1пф(Е,/) (1) уравнение (93,7а) приводится к виду
откуда Зф сЭ2ф df (t) dt *" dt,2 * dt • где посредством df/dt обозначена произвольная функция t. Переобозначением ф->-фе' (не меняющим искомой функции и(£, г)) это уравнение преобра- -^- = 11-^- (3) dt Ц dl2 • (d) Решение этого уравнения с начальным условием ф(£, 0) = фо(£) дается формулой (51,3): оо Ф (С, г) = 2(яц0-1/2 J Фо Ю ехр {- — оо Начальная же функция фо(£) связана с начальным значением искомой функции и\%, I) равенством I In ф„(С)= - Л-J щ(j)ц (5) о (выбор нижнего предела в интеграле произволен).
§ 94. Ударные волны в релаксирующей среде К значительному расширению ударной волны может привести наличие в газе сравнительно медленно протекающих релаксационных процессов — медленно протекающие химические реакции, замедленная передача энергии между различными
п. (Я. Б. Зельдович, степенями 1946) >). Пусть т — порядок величины времени релаксации. Как на- если интенсивность волны превышает определенный предел, то структура волны усложняется, в чем можно убедиться следующим образом. На рис. 67 сплошной линией изображена ударная адиабата, проведенная через заданную начальную точку /, в предположении полной равновесности конечных состояний газа; наклон касательной к этой кривой в точке / определяется «равновесной» скоростью звука, которую мы обозначали в § 81 посредством с0. Пунктиром же изображена ударная адиабата, проведенная через ту же точку /, в предположении, что релаксационные процессы «заморожены» и не происходят вовсе; наклон касательной к этой кривой в точке / определяется значением скорости звука, которое было обозначено в § 81 как сж. Если скорость ударной волны такова, что с0 < v\ < с», то хорда 12 расположена так, как указано на рис. 67 нижним отрезком. В этом случае мы получим простое расширение ударной волны, причем все промежуточные состояния между начальным состоянием / и конечным 2 изображаются в плоскости р, V точками на отрезке 12. Это следует из того, что (при пренебрежении обычными вязкостью и теплопроводностью) все последовательно проходимые газом состояния удовлетворяют уравнениям сохранения вещества pv = / = const и сохранения импульса р -f- j2V = const (ср. подробнее аналогичные соображения в § 129). Если же vi > с», то хорда занимает положение 11'2'. Все точки, лежащие на ее отрезке между точками / и вообще не соответствуют каким-либо реальным состояниям газа; первой
Дата добавления: 2014-11-07; Просмотров: 396; Нарушение авторских прав?; Мы поможем в написании вашей работы! Нам важно ваше мнение! Был ли полезен опубликованный материал? Да | Нет |