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

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

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

Еще

Гидравлическая сила, статический анализ, электр о гидравлическая рулевая машина

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

IDR: 143178834

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

В качестве исполнительных органов рулевых трактов систем управления вектором тяги (СУВТ) жидкостных ракетных двигателей космических разгонных блоков (РБ) широкое распространение получили автономные однокаскадные электрогидравлические рулевые машины (РМ) [1].

Поскольку РМ РБ чаще работают в режиме позиционирования, а не слежения, к таким РМ, как правило, не предъявляются динамические требования, а предъявляются только требования обеспечения статических характеристик, к коим относятся скоростные и силовые (или моментные) характеристики.

Автономные однокаскадные РМ могут выполняться по схемам с двух-и четырёхдроссельным электрогидрав-лическими усилителями (ЭГУ), при этом четырёхдроссельные ЭГУ могут иметь как положительное (нулевое), так и отрицательное перекрытие дроссельных окон золотниковыми плунжерами, тогда как двухдроссельные ЭГУ выполняются, как правило, с отрицательным перекрытием дроссельных окон золотниковыми плунжерами. Схемы таких рулевых машин представлены на рис. 1.

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

а)                                 б)                                 в)

Рис. 1. Электрогидравлические схемы РМ: а — с двухдроссельным ЭГУ с отрицательным перекрытием; б — с четырёхдроссельным ЭГУ с отрицательным перекрытием; в — с четырёхдроссельным ЭГУ с нулевым (положительным) перекрытием. 1 — поляризованное реле (электромеханический преобразователь); 2 — коромысло; 3 — плоская нагрузочная пружина; 4 — золотниковый плунжер; 5 — предохранительный клапан; 6 — электродвигатель; 7 — трёхшестерённый насос; 8 — трубопровод; 9 — поршень; 10 — силовой гидроцилиндр

В основу исследования составляющих гидравлических сил РМ были положены:

  • •    результаты исследований гидродинамических сил, действующих на золотниковые плунжеры [2–6];

  • •    итерационные методы статического анализа РМ [7, 8].

  • •    результаты исследования влияния конструкционных и регулировочных параметров ЭГУ на статические характеристики РМ [9, 10];

  • •    результаты оптимизации РМ по критерию электропотребления [11];

В процессе проведения статического анализа [7, 8] исследования влияния конструкционных и регулировочных параметров ЭГУ на статические характеристики [9, 10] и оптимизации РМ [11] было выявлено, что силовая характеристика РМ с четырёхдроссельным ЭГУ с отрицательным перекрытием дроссельных окон золотниковыми плунжерами не является однозначной. Вычислительные эксперименты показали, что при заторможенном поршне РМ при малейшем смещении коромысла ЭГУ (имитация силового воздействия) в ту или другую сторону золотниковые плунжеры сваливаются под действием гидравлической силы (суммы гидростатической силы и стационарной составляющей гидродинамической силы) на значительную величину в сторону смещения (сторону силового воздействия).

Это явление иллюстрируется рис. 2, 3. И в том, и в другом случае зависимости гидравлической силы от перемещения золотниковых плунжеров при заторможенном поршне РМ F г = f ( X з) не начинаются из начала координат, что указывает на статическую неустойчивость такого типа РМ.

Рис. 2. Первый вариант зависимости F г = f ( X з) РМ с четырёхдроссельным ЭГУ при заторможенном поршне

Соответствующие этим смещениям золотниковых плунжеров ( X з) силовые характеристики представлены на рис. 4.

Таким образом, в результате проведённых исследований было установлено, что нейтральное положение золотниковых плунжеров ЭГУ РМ (Xз = 0) является точкой бифуркации [12], а точки, в которые смещаются золотниковые плунжеры при незначительных внешних воздействиях, являются аттракторами [12]. При этом составляющие четырёхдроссельного ЭГУ простейшие трёхходовые золотниковые гидрораспределители с отрицательным перекрытием дроссельных окон золотником таким свойством, согласно работе [13], не обладают, и их зависимости Fг = f(Xз) начинаются из начала координат.

Рис. 3. Второй вариант зависимости F г = f ( X з) РМ с четырёхдроссельным ЭГУ при заторможенном поршне

Рис. 4. Зависимости усилия на поршне F п = f ( 1 к ) РМ с четырёхдроссельным ЭГУ: 1 — первый вариант; 2 — второй вариант

Так проявляется потенцированное действие четырёхдроссельного ЭГУ или синергизм его рабочего процесса [14], т. е. результирующий эффект общего действия (совместное воздействие факторов обтекания золотниковых плунжеров потоками жидкости через дроссельные окна и вращения гильз (полых осей трёхшес-терённого насоса)) превосходит сумму эффектов отдельного действия указанных факторов [15, 16].

При этом силовая характеристика РМ с двухдроссельным ЭГУ с отрицательным перекрытием дроссельных окон золотниковыми плунжерами, представленная на рис. 5, как и соответствующая ей зависимость гидравлической силы от перемещения золотниковых плунжеров F г = f ( X з) (рис. 6) в аналогичных условиях демонстрирует однозначный характер [9], т. е. существование только одного вида характеристики. Кроме этого, из графика на рис. 6 видно, что зависимость F г = f ( X з) при заторможенном поршне РМ начинается из начала координат, что указывает на статическую устойчивость такого типа РМ.

Рис. 5. Зависимость F п = f ( I к ) РМ с двухдроссельным ЭГУ

Рис. 6. Зависимость F г = f ( X з) РМ с двухдроссельным

ЭГУ при заторможенном поршне

Это даёт право утверждать, что составляющие двухдроссельного ЭГУ простейшие двухходовые золотниковые гидрораспределители оказывают аддитивное действие [17, 18], т. е. результирующий эффект общего действия (воздействия факторов обтекания золотниковых плунжеров потоками жидкости через дроссельные окна и вращения гильз (полых осей трёхшестерённого насоса)) равен сумме отдельных эффектов указанных факторов.

При статических режимах работы РМ гидравлические силы F г, действующие на золотниковые плунжеры ЭГУ обоих типов, слагаются из двух составляющих [7, 8]:

  • •    гидростатической силы F гc;

  • •    стационарной составляющей гидродинамической силы F гдс, т. е.

F г = F гc + F гдс .

Гидростатические силы описываются уравнениями:

  • •    (1) — для двухдроссельного ЭГУ [7]

F гc = S т1( p г1 р г2) + S т2( р з1 р з2),        (1)

где S т1, S т2 — площади внешней и внутренней торцевых поверхностей золотникового плунжера; p г1, р г2 — давления в полостях начала каналов гильз; р з1, р з2 — давления в камерах золотниковых плунжеров;

  • •    (2) — для четырёхдроссельного ЭГУ [8]

F гc = S 12( pz 1 рz 2 pz 3 + рz 4) +

+ S 34( p г1 р г2),                             (2)

где S 12, S 34 — площади наливных и сливных торцевых поверхностей пояска золотникового плунжера; pz 1, рz 2, pz 3, рz 4 — давления в началах и концах полостей между штоками золотниковых плунжеров и гильзами.

Стационарные составляющие гидродинамических сил описываются уравнениями:

  • •    (3) — для двухдроссельного ЭГУ [7]

    n


    F гдс = [2 S 2 кг ( p г1 р г3 )]/ ζ кг. i + λ кг1 ( l кг / d кг )


    г - 1


    [2 S 2 кз ( p г1 р з1 )] /


    n

    ζ кз. i + λ кз1( l кз / d кз)

    г -1


    – [2 S 2кг( p г2 р г4)] /


    n

    ζ кг. i + λ кг2( l кг/ d кг)

    г - 1


    +


    [2 S 2кз( p г2 р з2)] /


    n

    ζ кз. i + λ кз2( l кз / d кз)

    г - 1


    +


+ [2 m озµ2оз1 S оз р з1cos(θ/2)] / εоз1 – [2 m озµ2оз2 S оз р з2cos(θ/2)] / εоз2

– [2nоµ2ос1Sос1(рр1 – рг1)cos(βос1)] / εос1 + [2nоµ2ос1Sос2(рр2 – рг2)cos(βос2)] / εос2, где βос1, βос2 — углы истечения потоков рабочей жидкости в сечениях сегментных дроссельных окон золотникового гидрораспределителя; nо — количество дроссельных окон; θ — угол истечения потоков рабочей жидкости в отверстиях золотниковых плунжеров; mоз — количество отверстий в золотниковом плунжере; Sоз — площадь проходного сечения отверстия в золотниковом плунжере; µоз1, µоз2 — коэффициенты расхода отверстий nn в золотниковом плунжере; ∑ζкг.i, ∑ζкз.i — суммы коэффициентов местных гидрав- лических сопротивлений, обусловленных изменениями параметров русла каналов гильз и каналов золотниковых плунжеров, соответственно; λкг1, λкг2, λкз1, λкз2 — коэффициенты гидравлических потерь на трение по длине каналов гильз и каналов золотниковых плунжеров, соответственно; lкг, lкз — длины каналов гильз и каналов золотниковых плунжеров, соответственно; dкг, dкз — диаметры каналов гильз и каналов золотниковых плунжеров, соответственно; Sкг, Sкз — площади проходных сечений каналов гильз и каналов золотников, соответственно;

  • •    (4) — для четырёхдроссельного ЭГУ [8]

    F гдс = µ2 о1 S 2 о1 ( p z 3 p p1 ){[2 n о cos(β о1 )] /( S о1 ε о1 ) – [(2 n 2 о )/ S 12 ](ρ z 1. z 3 z 3.p1 )} –


    – µ2 о3 S 2 о3 ( p p1 p г1 ){[2 n о cos(β о3 )] /( S о3 ε о3 ) – [(2 n 2 о )/ S 34 ](ρ г1.г3 р1.г1 )} –


    – µ2 о2 S 2 о2 ( p z 4 p р2 ){[2 n о cos(β о2 )] /( S о2 ε о2 ) – [(2 n 2 о )/ S 12 ](ρ z 2. z 4 z 4.р2 )} +



    + µ2 о4 S 2 о4 ( p p2 p г2 ){[2 n о cos(β о4 )]/( S о4 ε о4 ) – [(2 n 2 о )/ S 34 ](ρ г2.г4 р2.г1 )},


где βо1, βо2, βо3, βо4 — углы истечения потоков рабочей жидкости в сечениях дроссельных окон золотникового гидрораспределителя; S о1, S о2, S о3, S о4 — площади проходных сечений дроссельных окон; µо1, µо2, µо3, µо4 — коэффициенты расхода дроссельных окон; εо1, εо2, εо3, εо4 — коэффициенты сжатия потоков в дроссельных окнах; n о — количество дроссельных окон; ρ z 1. z 3, ρ z 2. z 4 — средние значения плотности жидкости в полостях между штоками золотниковых плунжеров и гильзами; ρ z 3.p1, ρ z 4.p2 — средние значения плотности жидкости в наливных дроссельных окнах; ρр1.г1, ρр2.г1 — средние значения плотности жидкости в сливных дроссельных окнах; ρг1.г3, ρг2.г4 — средние значения плотности жидкости в каналах гильз; p p1, p p2 — давления в рабочих полостях ЭГУ РМ.

Постановка задачи

Для изучения возможности устранения статической неустойчивости РМ с четырёхдроссельным ЭГУ необходимо выяснить её причины.

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

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

результаты вычислительных экспериментов

Исследования составляющих гидравлических сил, действующих на золотниковые плунжеры РМ с двух- и четырёхдроссельным ЭГУ с отрицательным перекрытием дроссельных окон золотниковыми плунжерами, проводились методом статического анализа [7, 8], в результате которого были определены зависимости F гс = f ( X з) и F гдс = f ( X з) при заторможенном штоке РМ при расчётах силовых характеристик.

Исследования проводились с использованием математических моделей оптимизированных РМ [11].

Результаты расчётов гидростатических сил, действующих на золотниковые плунжеры РМ с двух- и четырёхдроссельным ЭГУ, в зависимости от

Рис. 7. Зависимость F гс = f ( X з) РМ с двухдроссельным ЭГУ при заторможенном поршне

Рис. 8. Первый вариант зависимости F гс = f ( X з) РМ с четырёхдроссельным ЭГУ при заторможенном поршне

Результаты расчётов стационарных составляющих гидродинамических сил, действующих на золотниковые плунжеры РМ с двухдроссельным ЭГУ, представлены на рис. 10.

Результаты расчётов стационарных составляющих гидродинамических сил, действующих на золотниковые плунжеры РМ с четырёхдроссельным ЭГУ, представлены на рис. 11, 12.

Рис. 9. Второй вариант зависимости F= f ( Х з ) РМ с четырёхдроссельным ЭГУ при заторможенном поршне

Рис. 12. Второй вариант зависимости F _дс = f ( X з) РМ с четырёхдроссельным ЭГУ при заторможенном поршне

Рис. 10. Зависимость F гд = f ( X з) РМ с двухдроссельным

ЭГУ при заторможенном поршне

Рис. 11. Первый вариант зависимости F гдс = f ( X з) РМ с четырёхдроссельным ЭГУ при заторможенном поршне

Из рассмотрения графиков F гс = f ( X з) рис. 7–9 следует, что гидростатические силы, действующие на торцы поясков золотниковых плунжеров РМ с четырёхдроссельным ЭГУ, ожидаемо на порядок больше гидростатических сил, действующих на торцевые поверхности золотниковых плунжеров РМ с двухдроссельным ЭГУ.

Из рассмотрения графиков F гдс = f ( X з) рис. 11, 12 видно, что зависимости стационарной составляющей гидродинамической силы от перемещения золотниковых плунжеров при заторможенном поршне РМ с четырёхдроссельным ЭГУ не начинаются из начала координат и практически повторяют зависимости F г = f ( X з), представленные на рис. 2, 3, в то время как аналогичная зависимость при заторможенном поршне РМ с двухдроссельным ЭГУ (рис. 10) из начала координат начинается и практически повторяет зависимость, представленную на рис. 6.

Сравнительный анализ значений стационарных составляющих гидродинамических сил (рис. 10–12) с гидростатическими силами, представленный на рис. 7–9, показывает, что гидростатические силы пренебрежимо малы по сравнению со стационарными составляющими гидродинамических сил и поэтому не могут быть причиной статической неустойчивости РМ.

Причиной статической неустойчивости РМ с четырёхдроссельным ЭГУ с отрицательным перекрытием дроссельных окон золотниковыми плунжерами является стационарная составляющая гидродинамической силы.

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

Анализ формул (3) и (4), а также то, что углы истечения потоков в дроссельных окнах зависят от фактора вращения гильзы [7, 8], показывают, что найти аналитическое решение возможности уменьшения стационарной составляющей гидродинамической силы практически невозможно. Проведённые исследования влияния конструкционных и регулировочных параметров ЭГУ на статические характеристики РМ [10] показали, что повысить статическую устойчивость РМ с четырёхдроссельным ЭГУ возможно увеличением относительного коэффициента упругости плоской нагрузочной пружины узла управления ЭГУ РМ и одновременным увеличением относительной ширины дроссельных отверстий в гильзе при сохранении начальных площадей дроссельных окон. Однако лабораторная отработка РМ показала, что в этом случае для обеспечения требуемых статических характеристик РМ значительно уменьшается величина перемещения золотниковых плунжеров. И эта величина становится соизмерима с величиной колебаний золотниковых плунжеров в условиях воздействия вибраций, ударов и линейных ускорений в процессе эксплуатации, что делает РМ практически неуправляемой.

Полученный результат показывает, что стационарная составляющая гидродинамической силы, действующая на золотниковые плунжеры РМ с четырёхдроссельным ЭГУ, является дестабилизирующим фактором, который может приводить к возникновению автоколебаний в рулевых трактах системы управления вектором тяги жидкостных ракетных двигателей РБ.

Так, в рулевых трактах РБ типа Д и ДМ долгое время применялась РМ с четырёхдроссельным ЭГУ с отрицательным перекрытием дроссельных окон золотниковыми плунжерами, всё это время в них периодически возникали автоколебания, которые с переменным успехом пытались устранить либо перенастройкой, либо заменой РМ. Кроме этого, автоколебания выходного штока РМ иногда наблюдались и на испытательном стенде РМ даже при отсутствии замыкания РМ электрической обратной связью.

Окончательно избавить рулевые тракты РБ от автоколебаний удалось только после внедрения на РБ ДМ– SL оптимизированной РМ с двухдроссельным ЭГУ [11] с отрицательным перекрытием дроссельных окон золотниковыми плунжерами. А затем такие РМ стали применять на РБ ДМ– SL Б и на новом варианте РБ типа ДМ.

Применение разработанной высокостабильной и высокоустойчивой РМ с двухдроссельным ЭГУ с отрицательным перекрытием дроссельных окон золотниковыми плунжерами в новых перспективных разработках РБ должно обеспечить устойчивость рулевых трактов систем управления вектором тяги жидкостных ракетных двигателей.

заключение

В результате проведённых исследований установлено, что причиной статической неустойчивости РМ с четырёхдроссельным ЭГУ с отрицательным перекрытием дроссельных окон золотниковыми плунжерами является стационарная составляющая гидродинамической силы, действующая на золотниковые плунжеры ЭГУ РМ.

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

Для исключения возникновения автоколебаний в рулевых трактах СУВТ жидкостных ракетных двигателей космических РБ рекомендуется не применять в качестве исполнительных органов СУВТ РМ с четырёхдроссельным ЭГУ с отрицательным перекрытием дроссельных окон золотниковыми плунжерами.

Разработанная высокостабильная и высокоустойчивая РМ с двухдроссельным ЭГУ с отрицательным перекрытием дроссельных окон золотниковыми плунжерами обеспечит устойчивость рулевых трактов систем управления вектором тяги жидкостных ракетных двигателей в новых перспективных разработках разгонных блоков.

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

  • Kudryavtsev V.V., Stepan G.A., Shutenko V.I., Chertok B.E. The rocket steering actuators // IAC'94 International Aerospace congress. Theory, Applications, Technologies. Abstracts. August 15-19, 1994, Moscow, Russia.
  • Попов Д.Н., Сосновский Н.Г., Княжанский А.А. Идентификация гидродинамических характеристик золотникового распределителя // Наука и образование: Науч. изд. МГТУ им. Н.Э. Баумана. 2015. № 05. С. 26-39.
  • Солман М.И., Попов Д.Н. Компьютерное исследование и расчёт гидродинамических нагрузок на золотник // Наука и образование: Науч. изд. МГТУ им. Н.Э. Баумана. 2012. № 10. С. 79-92.
  • Попов Д.Н., Сосновский Н.Г., Сиу-хин М.В. Компьютерное моделирование гидросистем с типовыми нелинейными характеристиками // Наука и образование: Науч. изд. МГТУ им. Н.Э. Баумана. 2017. № 07. С. 63-74.
  • Кишкевич П.Н. Статический и динамический расчёт гидро- и пневмо-распределителей: Уч.-метод. пос. для студ. спец. 1-36 01 07 «Гидропневмо-системы мобильных и технологических машин» / П.Н. Кишкевич, М.И. Жилевич, П.Р. Бартош. Минск: БНТУ, 2012. 80 с.
  • Даршт Я.А. Гидропривод и средства автоматики: Уч. пос. Ковров: ФГБОУ ВПО «КГТА им. В.А. Дегтярева», 2012. 108 с.
  • Белоногов О.Б. Итерационные методы статического анализа двухдрос-сельной электрогидравлической рулевой машины ракетных блоков // Космическая техника и технологии. 2018. № 2(21). С. 93-105.
  • Белоногов О.Б. Итерационные методы статического анализа четырёх-дроссельной электрогидравлической рулевой машины ракетных блоков // Космическая техника и технологии. 2019. № 2(25). C. 115-126.
  • Белоногов О.Б. Исследование влияния конструкционных и регулировочных параметров двухдроссельного электрогидравлического усилителя на статические характеристики рулевой машины ракетного блока / / Известия РАН. Энергетика. 2020. № 1. С. 40-49.
  • Belonogov O.B. Investigation of the influence of structural and adjusting parameters of a quadruple -orifices electro hydraulic amplifier on the static characteristics of rocket pack steering actuator // Herald of the Bauman Moscow state technical university. Series mechanical engineering. 2019. № 6. P. 32-40.
  • Belonogov O.B. Scheme-parametric optimization by the criterion of electric consumption of autonomous one-rock electrohydraulic steering actuators of space rocket packs // Journal of Computer and Systems Sciences International. 2020. V. 59. № 3. P. 406-414.
  • Малинецкий Г.Г., Потапов А.Б. Современные проблемы нелинейной динамики. М.: Эдиториал УРСС, 2000. 336 с.
  • Taft C.K., Twill J.P. An Analysis of the Three-Way Underlepped Hydraulic Spool Servovalve // Transaction of ASME. Journal of Dynamic sistems, measurment, and control. June 1978. V. 100. P. 117-123.
  • Данилов Ю.А. Лекции по нелинейной динамике. Элементарное введение: Уч. пос. Изд. 2-е, испр. М.: КомКнига, 2006. 208 с.
  • Черногор А.Ф. О нелинейности в природе и науке. Харьков: ХНУ имени В.Н. Каразина, 2009. 386 с.
  • Черногор А.Ф. Нелшшна радюф1зика, шдручник / 3-те вид доп. I перероб. Харшв: ХНУ 1меш В.Н. Каразша, 2016. 208 с.
  • Аддитивность // Большая Российская Энциклопедия. М.: Науч. изд-во БРЭ, 2005. Т. 1. С. 221.
  • Мелких А.В., Повзнер А.А. Физика нелинейных явлений: Уч. пос. Екатеринбург: ЧГТУ-УПИ, 2009. 144 с.
Еще
Статья научная