Динамика ротора на комбинированных опорах с центробежными переключателями

Автор: Поляков Роман Николаевич, Савин Леонид Алексеевич, Бондаренко Максим Эдуардович

Журнал: Известия Самарского научного центра Российской академии наук @izvestiya-ssc

Рубрика: Машиностроение и машиноведение

Статья в выпуске: 4-6 т.18, 2016 года.

Бесплатный доступ

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

Ротор, подшипник качения, подшипник скольжения, жесткость, демпфирование, гидродинамическая теория смазки, теория герца

Короткий адрес: https://sciup.org/148204814

IDR: 148204814

Текст научной статьи Динамика ротора на комбинированных опорах с центробежными переключателями

ческих топливных элементов электромобилей, турбодетандеры и различные виды насосов для криогенной техники и нефтехимической промышленности) актуальна задача обеспечения повышенной надежности, что предъявляет более жесткие требования к опорным узлам роторов такого рода машин. И когда ни один из существующих видов подшипников не удовлетворяет техническим требованиям роторной машины, как возможное решение данной проблемы рассматривают вариант совмещения опор различного принципа действия для повышения надежности опорного узла за счет разделения и дублирования функций подшипников. Одним из вариантов является совмещение в едином опорном узле подшипников качения и скольжения по различным пространственным схемам, что позволяет использовать достоинства и исключить недостатки, присущие каждому виду опор на различных режимах работы. Конструкции комбинированных опор и принципы их действия повышают надежность опорного узла за счет разделения и/или дублирования функций подшипников на различных режимах работы, что приводит к выполнению основных критериев работоспособности даже при существенном ужесточении эксплуатационных требований к опорным узлам.

Выделяют два принципиально различных способа построения комбинированных опор: разделение нагрузок или разделение скоростей [1–6] (рис. 1), также существуют варианты, когда одна и та же конструкция может работать как в режиме разделения скоростей, так и в режиме разделения нагрузок [7]. Принцип работы комбинированной опоры, функционирующей по принципу разделения нагрузок (КОРН), заключается в том, что на основном рабочем режиме ПК работает с меньшей нагрузкой по сравнению с его одиночной постановкой. Разгрузка ПК осуществляется за счет возникновения гидродинамических или гидростатических реакций ПС (рис. 1а). Комбинированная опора с разделением скоростей (КОРС) обеспечивает вращение ротора в режиме пуска-останова в ПК, а по мере увеличения скорости вращения и роста реакции ПС, опора качения из работы выключается и передача нагрузки на корпус осуществляется через смазочный слой ПС (рис. 1 б,в).

Для перспективных высокоскоростных роторных машин с линейными скоростями на поверхности цапфы ротора порядка 102 м/с и выше наиболее перспективной является применение комбинированных опор по последовательной схеме расположения относительно оси вала и работающих по принципу разделения скоростей (КОРС). Например, в опоре ротора турбонасоса МК48 для перекачки жидкого кислорода корабля многоразового использования Shuttle [3] применена постоянно включенная КОРС, обеспечивающая стабильную работу роторно-опорного узла при частоте вращения 95000 об/мин (210 м/с) (рисунок 2а). Отечественными учеными была предложена адаптивная КОРС с переключением по давлению [8] с рабочей частотой вращения 10000 об/мин (75 м/с) (рис. 2б). Вместе с тем мож- но говорить о том, что применение комбинированных опор носит экспериментальный характер и отсутствуют единые методики их расчета и проектирования. При этом количество изобретений по классу МПК F16 C21/00 «Комбинации из подшипников скольжения с подшипниками качения» насчитывает тысячи единиц, в которых описываются положительные качества и предпосылки проектирования такого вида опор [9-12].

Механизм работы КОРС предполагает следующие направления улучшения рабочих характеристик опорного узла:

  • 1)    уменьшение повреждающих моментов, действующих на втулку подшипника скольжения во время переходных режимов, что ведет к повышению его ресурса;

  • 2)    уменьшение центробежных нагрузок в подшипнике качения за счет разделения скоростей;

  • 3)    отстранение от режимов самовозбуждаю-щихся колебаний, характерных движению роторов в подшипниках скольжения: во время неустойчивого движения ротора опасность существенных повреждений (например, заклинивания цапфы во втулке ПС) практически исключается за счет возможности дополнительного поворота в ПК;

  • 4)    повышение возможностей обеспечения необходимых динамических характеристик опорного узла для обеспечения прохождения резонансных режимов с допустимыми амплитудами колебаний;


    Рис. 1. Схемы комбинированных опор



    a)


    б)

    Рис. 2. Конструкции комбинированных опор


  • 5)    повышение надежности опорного узла за счет повышенной работоспособности опорного узла в аварийных ситуациях и при неустойчивой работе за счет дублирования функций ПК и ПС;

  • 6)    отстройка от режимов турбулентности за счет адаптивного перераспределения скоростей между подшипником качения и скольжения.

Цель работы: обоснование положительных эффектов для КОРС с центробежным переключателем в виде упругого лепестка теоретическими и экспериментальными методами. В работе рассматривается комбинированная опора [13], схема которой представлена на рис. 3. Комбинированная опора работает следующим образом:

  • 1)    в момент пуска машины передача нагрузки на корпус 1 происходит через подшипник качения 2 и упругие лепестки 4, которые за счет предварительной деформации центрируют вал 3;

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

  • 3)    при достижении определенной частоты

вращения (назовем ее частотой переключения ю пер ) когда центробежные силы станут настолько велики, что превысят силы упругости от предварительной деформации произойдет полное расцепление лепестков и наружного кольца ПК, а в образовавшемся зазоре возникнет гидродинамическая реакция; нагрузка с вала 3 на корпус 1 передается последовательно через ПС и ПК, наружное кольцо ПК вращается с частотой меньшей, чем частота вращения ротора.

МЕТОДИКА ИССЛЕДОВАНИЙ

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

Математическая модель радиально-упорного подшипника качения (рис. 4) основывается на

Рис. 3. Схема комбинированной опоры с центробежным переключателем

Рис. 4. Схема шарикового радиально-упорного подшипника качения (а) и эпюра распределения нагрузок между телами качения (б)

б)

решении контактной задачи теории упругости [14, 15] и строилась при следующих допущениях:

  • . рассматривается подшипник с идеальными геометрическими параметрами: отсутствуют микронеровности, отклонения от формы всех деталей ПК;

. материалы контактирующих тел однородны по своим механическим свойствам и абсолютно упругие;

. нагрузка приложена перпендикулярно поверхности контакта;

. деформации в контакте малы по сравнению с размерами площадки;

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

. Реакция со стороны тел качения в зависи- мости от смещения центра внутреннего кольца в направлении координат X, Y, Z

R™ = Kx^^ ^(tgpsmpcosy0 + cospcosyj + tgy0cos/?sinyj) ^cospcosyj ;

R™ = Ky^ ^(tg/?sin/? s iny0 + ctgy0cos/?cosyy + cos^sinyy) ^cos^sinyy;

R™ = Kz’k (sin/? + ctg/?cos/?cosy0cosy, + ctg/?cos/?siny0sinyy) 2sin/j'.

где                – коэффициент пропор- k З^Р '2

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

k

Н/м3/2, здесь     - величина являющаяся ото бражением геометрии контакта и определяемая в зависимости от вспомогательной величины

;

- сумма главных кривизн поверхностей соприкосновения тел в месте их первоначального контакта;

- коэффициент, характеризующий упругие свойства материала, здесь                - соответствен но, коэффициенты Пуассона и модули упругости материалов контактирующих тел;

- деформация от предварительного натяга.

Математическое описание сил трения основывается на полуэмпирической методике [16]

p 0.55

M™ = M(F) + M((n) = 0S-l(T6\-^\ X \L0/

X (3Fa - 0.1Fr)D0 + 2 • 10-10O)2/3D03 , (2)

где – эквивалентная нагрузка на ПК; Fr, Fa - радиальная и осевая нагрузка, приведенные к подшипнику; X, Y - коэффициенты, учитывающие разное повреждающее действие радиальной и осевой нагрузки, опре- деляются в зависимости от соотношения Fr /C0, здесь C0 - статическая грузоподъемность ПК; V - коэффициент вращения; KБ, KТ - коэффициенты безопасности и температуры; v - кинематическая вязкость смазки; n – частота вращения.

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

Определение сил гладкого гидродинамиче- ского и гидростатодинамического с точечными камерами подшипников (рис. 3) основывается на интегрировании поля давлений p(x, z) по опорной поверхности на основании решения уравнения Рейнольдса

, (3)

Эх LpKx SxJ dz LpKz dzj Эх             r \ x где - функция радиального зазора; p, ц - плотность и динамическая вязкость смазки; U, V – скорость смазочного материала на поверхности цапфы в окружном и радиальном направлении.

Реакции и момент трения подшипника сколь- жения

L nltD

;

о J

L

И лО

;

о J

- ,пс       D rL rnD ,, ,

;

ip         2

h dp pU

2 ox h

Реакция ПС находится из решения уравнения Рейнольдса [17, 18], расчетная схем гидростатодинамического подшипника с точечными питающими камерами представлена на рис. 5а, на рис. 5б представлен пример расчета поля давлений.

Температура в смазочном слое определяется на основании численного решения уравнения баланса энергий

ph-Cp

QT dT (U dt + dx \2

/i2 dp\

12pKx dx/

dT / h? dp\ dz \12pKz dz /

£p + ^ ^E + y2^

dt 2 dx h

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

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

Рис. 5. Схема радиального гидростатодинамического подшипника (а) и расчетное поле давлений (б)

б)

строения 3-х мерных поверхностей зависимости вязкости и плотности от температуры и давления по табличным данным [19].

Для водорода (H 2 ), температура T =10 – 40 К:

  • -    динамическая вязкость Па х с;

  • -    плотность , кг/м 3 .

Для кислорода (О 2 ), температура T= 90 – 130 К:

  • -    динамическая вязкость Па х с;

  • -    плотность кг/м 3 .

Математическая модель лепестка базировалась на конечно-элементном анализе упругой криволинейной балки в системе автоматизированного проектирования APM Win Machine [20]. Расчетная схема лепестка представлена на рисунке 6а. Схема нагружения лепестка (рисунок 6б) менялась согласно имитационной модели работы КОРС: на этапе разгона действуют силы упругости от начальной деформации лепестка и центробежные силы; на основном режиме работы действуют сила упругости лепестка, центробежные и гидродинамические силы. Суммарная деформация лепестка определялась на основании принципа суперпозиции от каждой действующей силы.

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

Грузоподъемность КОРС:

  • -    на этапе разгона РУкорс = И^пк;

    Рис. 6. Конечно-элементная схема упругого лепестка (а) и схема нагружения (б)



    - на переходном режиме РИкорс = РИУЛ, здесь 1Уул - грузоподъемность упругого ле-

  • пестка;
  • -    на основном режиме И/корс = Wnc.

Момент трения (коэффициент трения, потери мощности на трение):

  • - на этапе разгона М^РС = М™;

    - на переходном режиме М^РС = Мул, здесь Мул - момент трения скольжения при проскальзывании упругого лепестка относительно

наружного кольца ПК;

- на основном режиме М^®РС = М^рс.

Расход смазочного материала через комби- нированную опору:

, где QnK, QnC - расход смазочного материала через ПК и ПС, соответственно;

– коэффициенты гидравлического сопротивления ПС и ПК, соответственно.

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

Температурный режим комбинированной опоры:

  • -    на этапе разгона tK0PC = Спк, здесь Спк -температура ПК;

  • -    на переходном режиме tK0PC = tnK + (ул , здесь tYn - температура нагрева лепестка от проскальзывания поверхностей упругого лепестка относительно наружного кольца ПК;

  • -    на основном режиме tK0PC = tnc.

Жесткость комбинированной опоры:

Г^КОРС _ КпКк™

  • -    на этапе разгона К = Пкул , здесь /<пк- жесткость ПК, КУЛ - жесткость упругого лепестка;

корс кпску;1

  • -    на переходном режиме К = кпсул , КОРС _                         +

здесь К      кпсул , - жесткость ПС;

rzKOPC КПСКУЛ

  • -    на основном режиме К = Псул

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

г> КОРС   о УЛ        о УЛ

  • -    на этапе разгона В = В , здесь - демпфирование упругого лепестка, которое определяется внутреннем трением в материале и конструкционным трением;

  • -    на переходном и основном режимах Вкорс _ руд + Впс, здесь Впс - демпфирование ПС.

Коэффициенты жесткости и демпфирования для ПК и ПС определяются в квазистационарной линейной постановке [21, 22].

РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ

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

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

  • -    диаметр вала d= 40 мм;

  • -    частота вращения 0...10000 об/мин;

  • -    смазочная среда - жидкий водород;

  • -    толщина лепестка 5 =100 мкм;

  • -    ширина лепестка b =72 мм;

  • -    длина лепестка l =18 мм;

  • -    радиус лепестка r =41 мм;

  • -    давление подачи смазочного материала: р0 =2,5 атм.

Особенность работы многолепесткового подшипника в КОРС состоит в том, что на основном режиме на него действует поле гидродинамического давления, с одной стороны, и поле центробежных сил с другой стороны (рис. 7а). Это создает предпосылки для управления зазором лепесткового подшипника на основном режиме, например, путем применения утяжелителей, которые могут быть размещены на нерабочей стороне лепестка (рис. 7б). Также можно снизить чувствительность лепестка к перепадам гидродинамического давления в результате колебаний ротора.

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

На рис. 9а представлена имитационная схема деформации лепестка в зависимости от скорости вращения, которая наглядно показывает как изменяется зазор лепесткового подшипника на переходном режиме. На рис. 9б представлена зависимость жесткости КОРС с центробежными

Рис . 7. Сложнонагруженное состояние лепестка в КОРС

‘ 1 б)

Рис. 8. Динамическая модель КОРС

б)

Рис. 9. Динамический изгиб упругого лепестка (а) и коэффициенты жесткости (КОРС) (б)

переключателями от угловой скорости. Анализируя представленные результаты, можно сказать, что на начальном этапе эквивалентная жесткость опорного узла растет, что связано с тем, что под действием центробежных сил упругий переключатель становится жестче. В тот момент, когда происходит отрыв лепестка от поверхности подшипника качения жесткость КОРС определяется жесткостью смазочного слоя, а так как в этот момент зазор весьма мал, то и жесткость смазочного слоя велика. С дальнейшем увеличением скорости зазор увеличивается – жесткость смазочного слоя и, соответственно, всей КОРС падает, а с дальнейшим ростом скорости растет.

Особенность жесткостных характеристик КОРС также заключается в том, что эквивалентная жесткость всего опорного узла не может быть больше жесткости самого податливого из ее элементов. То есть с увеличением угловой скорости жесткостная характеристика будет асимптотически приближаться к этому значению.

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

Используя полученные значения коэффициентов жесткости для проведения модального анализа [23] можно построить диаграмму Кэмпбелла (рис. 10). Скачок коэффициентов жесткости КОРС на этапе разгона приводит к скачку собственных частот, что приводит к появлению точек пересечения с частотой возбуждения и возникновению резонансных явлений на переходном режиме, на этапе останова происходят те же явления. Данная особенность требует более тщательного анализа поведения ротора на КОРС на переходных режимах и учета этого явления при выборе характеристик совмещаемых подшипников.

ВЫВОДЫ

Проведенный анализ результатов позволяет сделать вывод о существенных преимуществах применения комбинированных опорных узлов в агрегатах с тяжелыми условиями для работы подшипников: позволяет повысить ресурс опорного узла, повысить грузоподъемность, быстроход- ность, улучшить динамические характеристики с позиций уменьшения уровня вибраций и отстранения от режимов критических частот. Но достижение перечисленных положительных эффектов КО по сравнению с одиночной постановкой ПК или ПС возможно лишь в случае правильного выбора геометрических параметров опорных узлов при заданных рабочих (частота вращения, внешняя нагрузка, дисбаланс ротора, рабочая среда, температура, давление). Например, ПК для КОРС, выбранный только с позиций удовлетворения заданной геометрии конструкции может не обеспечить необходимые динамические характеристики всего опорного узла, выбор гидродинамического подшипника вместо гидростатического при низкой вязкости смазочного материала может не обеспечить режим разделения скоростей и т.п.

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

Можно сделать вывод о том, что комбинирование подшипников качения и скольжения дают

Рис. 10. Характерная диаграмма Кэмпбелла для ротора на КОРС

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

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

Список литературы Динамика ротора на комбинированных опорах с центробежными переключателями

  • Ханович М.Г. Опоры жидкостного трения и комбинированные. Л.: Машгиз, 1960. 272 с.
  • Понькин В.Н., Горюнов Л.В., Такмовцев В.В. Совмещенные опоры быстроходных турбомашин, принципы конструирования и экспериментальное исследование: препринт. Казань: Изд-во КГТУ, 2003. 62 с.
  • Hannum N.P. Nielson C.E. The performance and Application of High Speed Long Life LH2 Hybrid Bearing for Reusable Rocket Engine Turbomachinery//(NASA TM-83417) AIAA №83-1389, 1983. 26 p.
  • Butner M.F., Murphy B.T. Space shuttle main engine long-life bearings. Final report//NASA-CR-179455, Rockwell International Corp., 1986. 163 p.
  • Stek J.B. Design of Turbines for Rocket Engines, with Turbiness for a Spase Shuttke as Examples//Ingenieur. 1971. Vol. 83. № 43. P. 119 -130.
  • Найпен, Скиббе, Хемрок. Оптимальное распределение скоростей в сериесном комбинированном подшипнике//Проблемы трения и смазки. 1973. № 1. С. 83 -89.
  • Поляков Р.Н., Соломин О.В. Комбинированная опора//Патент России №2243425. 2004. Бюл. №36.
  • Демидович В.М., Горюнов Л.В., Клюшкин А.П., Якимов Н.А. Комбинированная опора//Патент СССР №1035311. Опубл. 15.08.1983.
  • Лентов В.В. Комбинированная опора//Патент России №2079015. Опубл. 10.05.1997.
  • Савин Л.А. Комбинированная опора//Патент России №2073801. 1997. Бюл. №5.
  • Демидович В.М., Горюнов Л.В., Клюшкин А.П., Якимов Н.А. Комбинированная опора//Патент СССР №1935311.1983. Бюл. №30.
  • Поляков Р.Н., Соломин О.В. Комбинированная опора//Патент России 2243425. 2004. Бюл.№36.
  • Поляков Р.Н., Савин Л.А., Соломин О.В., Пугачев А.О. Комбинированная опора//Патент России №2228470. 2004. Бюл.13.
  • Демидов С.П. Теория упругости: учебник для вузов. М.: Высш. школа, 1979. 432 с.,
  • Wensing J.A. On the dynamics of ball bearing: PhD thesis//University of Twente, Enschede, The Netherlands, 1998. 183 p.
  • Перель Л.Я., Филатов А.А. Подшипники качения: Расчет, проектирование и обслуживание опор: Справочник. 2-е изд., перераб. и доп. М.: Машиностроение, 1992. 608 с.
  • Коровчинский М.В. Теоретические основы работы подшипников скольжения. М.: Машгиз, 1959. 404 с.
  • Савин Л.А., Соломин О.В. Моделирование роторных систем с подшипниками жидкостного трения: монография. М.: Машиностроение-1, 2006. 334 с.
  • Варгафтик Н.Б. Справочник по теплофизическим свойствам газов и жидкостей. М.: Наука, 1972. 720 с.
  • НТЦ АПМ -CAD/CAE -оптимальные решения в машиностроении . URL: http://www.apm.ru (дата обращения 12.10.2016)
  • Лунд. Разработка понятия динамических коэффициентов радиальных подшипников жидкостного трения//Проблемы трения и смазки. 1987. №1. С. 40-44.
  • Кельзон А.С., Циманский Ю.П., Яковлев В.И. Динамика роторов в упругих опорах. М.: Наука, 1982. 280 с.
  • Krämer E. Dynamics of Rotors and Foundations. Springer-Verlag, 1993. 384 p.
Еще
Статья научная