Физика пневматического удара в клапанах и трубопроводных системах

Разберите физику пневматического удара, распространение волн давления, ограничения формулы Жуковского для воздуха и безопасные методы измерения переходных процессов клапана.

Поделиться
David Li, Главный технический консультант, Бепто Пневматик

Об авторе

David Li

Главный технический консультант

Я David, Главный технический консультант в Бепто Пневматик. Помогаю проверять требования к пневматическим продуктам, применение компонентов и технические детали до расчета предложения.

Статьи автораDavid@bepto.com

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

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

Ключевые выводы

  • При температуре около 20°C небольшое возмущение в идеальном сухом воздухе распространяется со скоростью примерно 343 m/s (NASA Glenn).
  • Сопоставляйте длительность срабатывания клапана с 2L/a{2L/a} для предварительной оценки, а не для предсказания пикового давления.
  • Формула Жуковского для жидкости не является универсальной формулой для сжатого воздуха.
  • Перед выбором решения используйте абсолютные параметры состояния и синхронные динамические измерения.

Что такое пневматический удар в системе?

Пневматический удар — это переходный процесс течения газа, а не задача о стационарном падении давления. Быстрое закрытие клапана — лишь один из возможных инициирующих факторов. Быстрое открытие, перемещение обратного клапана, соединение объемов с разным давлением, реакция регулятора и резкое изменение потребления также могут запускать возмущения давления и скорости (Rao и Eswaran, 1993).

Газ сжимаем, поэтому во время распространения возмущения могут одновременно меняться давление, плотность, температура и скорость. Поэтому задача отличается от расчета потери давления на клапане при постоянных рабочих условиях. ISO 6358-1 прямо относится к испытаниям стационарного испытаний расхода пневматических компонентов, а ISO 6358-3 оценивает стационарные расходные характеристики соединенных компонентов и трубопроводов (ISO 6358-1, 2013; ISO 6358-3, 2014). Звуковая проводимость и критическое отношение давлений полезны для подбора пропускной способности, но сами по себе не предсказывают историю переходного давления.

Во время диагностики используйте более точные названия неисправностей:

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

Первое событие информативнее самого громкого места. Рама, панели и плохо закрепленные трубы могут передавать звук события, начавшегося в другом месте. По нашему опыту, временное совмещение сигнала клапана и кривой давления — самый быстрый способ проверить ошибочно определенное место возникновения. Если конструкция скрывает начало события, добавьте запись звука или ускорения.

В отчете о пневматическом ударе следует указать изменившееся граничное условие. «Труба ударила» — это наблюдение. «Клапан V3 за 18 ms уменьшил соединение между работающим ответвлением и запертым объемом после клапана» — это проверяемая гипотеза о переходном процессе.

С какой скоростью возмущение давления распространяется в воздухе?

Для небольшого возмущения в идеальном газе первой оценкой для предварительного анализа служит местная скорость звука. NASA Glenn выводит эту зависимость из уравнений сохранения и изоэнтропического допущения о малом возмущении (NASA Glenn, вывод скорости звука, 2021):

a=γRsTa = \sqrt{\gamma R_s T}

где:

  • aa — скорость распространения малой волны в газе, в m/s
  • γ\gamma — отношение теплоемкостей
  • RsR_s — удельная газовая постоянная, в J/(kg·K)
  • TT — абсолютная температура, в K

Для сухого воздуха при температуре около 20°C, если принять γ1.4\gamma \approx 1.4, Rs287.05 J/(kg⋅K)R_s \approx 287.05\ \text{J/(kg·K)}, и T293.15 KT \approx 293.15\ \text{K}, получаем:

a343 m/sa \approx 343\ \text{m/s}

Это число оценивает момент прихода небольшого возмущения давления. Это не означает, что рабочий поток к приводу проходит по линии со скоростью 343 m/s, и не подтверждает величину пика давления. Реальное поведение также зависит от начального состояния потока, теплообмена, податливости труб, трения, ограничений, ответвлений, движения клапана и величины изменения давления.

Полезны две временные величины:

t1=Lat_{1} = \frac{L}{a}
tr=2Lat_{r} = \frac{2L}{a}

Здесь t1t_1 — время прихода в одну сторону на эффективной длине LL, а trt_r — простое время прохождения туда и обратно. В линии длиной 20 m при принятой скорости 343 m/s первый приход происходит примерно через 58 ms, а идеализированное прохождение туда и обратно занимает около 117 ms.

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

Сравнение времени Инженерная интерпретация Чего это не доказывает
срабатывание клапана намного медленнее, чем 2L/a{2L/a} за время события информация о давлении может многократно пройти по линии что выброс давления невозможен
срабатывание клапана сопоставимо с 2L/a{2L/a} движение границы и возврат волны могут сильно взаимодействовать конкретный коэффициент усиления
срабатывание клапана намного быстрее, чем 2L/a{2L/a} рассматривать событие как возможный быстрый переходный процесс что применимы уравнения гидравлического удара для жидкости

Это сравнение служит ориентиром для предварительного анализа, а не критерием «годен/не годен». Короткое ответвление, соединенное с большой камерой, может вести себя иначе, чем прямая труба такой же длины. Обратный клапан, отскакивающий после срабатывания, может создать несколько изменений границы вместо одного четкого закрытия.

Почему нужны абсолютные давление и температура?

Газовые уравнения требуют термодинамических параметров состояния, поэтому используйте абсолютное давление и абсолютную температуру. NASA Glenn приводит уравнение идеального газа в форме через плотность и отмечает ограниченность модели идеального газа (NASA Glenn, уравнение состояния):

pabs=ρRsTp_{abs} = \rho R_s T

Манометрическое давление измеряется относительно местного атмосферного давления. Перед оценкой плотности газа переведите его в абсолютное давление:

pabs=pg+patmp_{abs} = p_g + p_{atm}

Например, 6 bar(g) соответствуют примерно 7.0 bar(a), если местное атмосферное давление составляет около 1.0 bar. При 20°C соответствующая плотность идеального газа составляет примерно 8.3 kg/m³, а не значение, полученное при подстановке 6 bar непосредственно как абсолютного давления. В реальном расчете эти округленные значения следует заменить местным атмосферным давлением и фактической температурой газа.

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

Записывайте каждое давление в модели как bar(a), bar(g), kPa(a) или kPa(g). Неуказанный тип «6 bar» может привести к существенной ошибке плотности еще до запуска расчета переходного процесса.

Почему уравнение Жуковского не предсказывает пневматический удар?

Известное соотношение Жуковского получено для переходных процессов в жидкости и не должно представляться как универсальное уравнение максимального давления для сжатого воздуха. При анализе переходного процесса в сжимаемой среде рассчитывают изменение давления, скорости и плотности по заданной сети (Rao и Eswaran, 1993). Обычно соотношение для жидкости записывают так:

Δp=ρaΔv\Delta p = \rho a\,\Delta v

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

Rao и Eswaran применили метод характеристик к сложным сетям со сжимаемой средой и явно смоделировали резкие операции с клапаном. В более поздних исследованиях газопроводов сравнивались метод характеристик и метод конечных объемов для медленных и быстрых переходных процессов (Koo, 2022; Koo, 2022).

Уравнение для жидкости все же можно привести в инженерном обзоре как напоминание о важности быстрого изменения импульса. Но без вывода, соответствующего модели газа, граничным условиям и размеру возмущения, его нельзя называть «максимальным давлением пневматического удара».

Метод расчета Подходящее применение Главное ограничение
проводимость и критическое отношение давлений по ISO 6358 подбор компонентов по стационарному расходу не модель временной истории волны давления
скорость звука в идеальном газе первая оценка прихода волны допущения о малом возмущении и свойствах среды
L/aL/a и 2L/a{2L/a} временная оценка сопоставление геометрии со временем срабатывания клапана не уравнение амплитуды
соотношение Жуковского для жидкости оценка гидравлического удара при указанных допущениях не универсальная формула пика для сжатого воздуха
одномерная сжимаемая модель методом характеристик или конечных объемов история давления, скорости и плотности в заданной сети требует корректных моделей компонентов, границ и трения
переходный расчет сжимаемого потока CFD или взаимодействия жидкости и конструкции сложная локальная геометрия или связь с конструкцией требует валидации и значительно большего объема исходных данных

Ни один обоснованный расчет не начинается с правила «обычное давление умножить на пять». Пик может оказаться выбросом, провалом, колебанием или последовательностью локальных событий. Ответ зависит от конкретной системы.

Что происходит, когда волна достигает клапана, ответвления или замкнутого объема?

Граница отражает, пропускает или рассеивает часть возмущения в зависимости от ее динамического импеданса и состояния. Модели сжимаемых сетей представляют трубы, узлы, ответвления и меняющиеся во времени клапаны как связанные задачи с граничными условиями (Rao и Eswaran, 1993). В закрытом конце локальный массовый расход стремится к нулю, а узел распределяет возмущение между соединенными путями.

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

В недавнем исследовании взаимодействия потока и конструкции для линии природного газа газовой турбины пневматический удар при открытии и закрытии клапана моделировали совместным расчетом потока и конструкции. Результаты относятся к той системе высокого давления, но показывают важный общий принцип: закон движения клапана, геометрию, историю волны давления и отклик конструкции необходимо рассматривать совместно (Frontiers in Energy Research, 2026).

Временная схема переходного процесса и измерений в пневмосистеме Вертикальная инженерная схема показывает изменение границы клапана, распространение волны, отражение на границах сети и синхронное измерение давления. Боковое сравнение сопоставляет время срабатывания клапана со временем акустического прохождения туда и обратно. Следите за событием, а не только за шумом 1. Изменяется границаплощадь прохода клапана, положение обратного клапана или соединенный объем меняются со временем 2. Распространяются возмущения давления и скоростипервая оценка: a = √(γRₛT)время прихода не равно времени полного заполнения камеры 3. Границы отражают и пропускаютответвления, замкнутые объемы, ресиверы, регуляторы, клапаны и фитингикоэффициенты зависят от фактического динамического состояния 4. Датчики фиксируют историю во временисинхронизируйте давление до клапана, давление после клапана и состояние клапанапроверьте динамику преобразователя, монтаж, диапазон и временную базу Измеренное время срабатывания клапанаtᵥ по данным о положении или площади прохода Акустический масштаб прохождения туда и обратноtᵣ = 2L/a Сопоставьте два времени при планировании испытания; не используйте это сравнение как множитель пикового давления.
Обоснованный анализ переходного процесса связывает меняющуюся границу клапана с распространением волны, отражениями в сети и синхронными измерениями. Сравнение временных масштабов помогает спланировать испытание, но не рассчитывает пиковое давление. Источники: вывод скорости звука NASA Glenn и модель сжимаемого трубопровода Rao и Eswaran.

Какие исходные данные нужны для обоснованной модели переходного процесса?

Модель настолько надежна, насколько надежны ее начальное состояние, граничные функции, данные компонентов и проверочная осциллограмма. Dorao и Fernandino описывают динамику газопровода как одномерное переходное течение сжимаемого газа, управляемое изменением потребления и регулирующими устройствами (журнале Natural Gas Science and Engineering, 2011). Заводская пневматическая линия меньше, но к ней применима та же дисциплина моделирования.

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

  1. Начальное состояние газа: состав газа, абсолютное давление, температура и начальный массовый расход или скорость.
  2. Геометрия сети: фактические внутренние диаметры и длины, перепад высоты при необходимости, фитинги, ответвления, тупиковые участки, камеры и ресиверы.
  3. Свойства труб и трубок: материал, толщина стенки, закрепление, податливость, температурный предел и номинальное рабочее давление.
  4. История граничного поведения клапана: эффективная площадь прохода или коэффициент во времени, работа пилотной части, ход золотника или тарелки, динамика обратного клапана и утечки.
  5. Модели компонентов: регуляторы, фильтры, глушители, быстроразъемные соединения, камеры цилиндров, предохранительные устройства и выпускные границы.
  6. Тепловые допущения и допущения о трении: изотермический, адиабатический режим или учет теплообмена; стационарное или нестационарное трение; взаимодействие со стенкой.
  7. Данные валидации: синхронные кривые давления в точках, позволяющих различить источник, распространение и отражение.

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

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

Как измерять переходный процесс давления в пневмосистеме?

Используйте подходящие для динамики датчики давления по обе стороны предполагаемого источника и синхронизируйте измерения с состоянием клапана. NIST отмечает, что процессы, зависящие от времени, и волны давления требуют динамических измерений давления, а в некоторых применениях важен отклик на уровне микросекунд (NIST Pressure/Vacuum Calibrations, updated 2026). Это не означает, что для каждого заводского испытания нужна микросекундная система. Это означает, что отклик прибора должен быть обоснован относительно события, которое требуется разрешить.

Используйте следующую последовательность полевых измерений:

  1. Перед установкой или перемещением датчиков выполните процедуру предприятия по контролю опасной энергии.
  2. Выбирайте преобразователи с подходящим абсолютным или манометрическим отсчетом, диапазоном, перегрузочной способностью, полосой пропускания или временем нарастания, совместимостью со средой и подтверждением калибровки.
  3. Делайте соединительные каналы короткими и фиксируйте каждый адаптер или испытательный шланг, добавляющий мертвый объем.
  4. По возможности установите один датчик до, а другой после предполагаемого клапана, ограничения или ответвления.
  5. Записывайте команду клапану и, если возможно, положение золотника, тарелки или привода по тем же часам.
  6. Выбирайте частоту сбора по полосе датчика и длительности события, а не по универсальному правилу количества отсчетов в секунду.
  7. До изменения одного контролируемого условия запишите несколько штатных циклов.
  8. Сравнивайте между каналами время прихода, форму давления, затухание и повторяемость.

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

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

Как снизить риск пневматического удара, не создав другую опасность?

Устраните подтвержденную причину, сохранив требуемую реакцию машины на безопасность. ISO 4414 рассматривает опасности пневматического оборудования на уровне системы, а не сертифицирует универсальное исправление отдельного компонента (ISO 4414, подтверждено в 2021 году). Решение может включать выбор и размещение клапана, объем трубок, управляемое повышение давления, пропускную способность выпуска, крепление, локальный ресивер или другую последовательность управления.

Начинайте с механизма, который показывает осциллограмма:

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

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

Не замедляйте аварийный или защитный выпуск только ради улучшения кривой давления. ISO 4414 касается существенных опасностей в пневматических системах машин и требует учитывать безопасность на уровне системы (ISO 4414, подтверждено в 2021 году). В США стандарт OSHA 29 CFR 1910.147 включает пневматическую энергию и требует во время соответствующего обслуживания снять, отключить, удержать или иным образом обезопасить опасную запасенную или остаточную энергию (OSHA 1910.147). Требуемое поведение при остановке, изоляции, выпуске, сбросе и проверке имеет приоритет над общим правилом защиты от удара.

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

Контрольный список инженерной проверки

Перед утверждением изменения конструкции убедитесь, что проверка отвечает на следующие вопросы:

  • Является ли событие переходным процессом в сухом газе, ударом конденсата, ударом привода, ограничением выпуска или стационарным падением давления?
  • Какая граница изменилась первой и какие данные подтверждают момент изменения?
  • Обозначены ли во всех расчетах состояния газа абсолютное давление и абсолютная температура?
  • Известна или измерена ли история эффективной площади прохода клапана?
  • Использовались ли L/aL/a и 2L/a{2L/a} использовались ли только как временные масштабы?
  • Исключена ли формула Жуковского для жидкости как универсальный расчет пика в газе?
  • Достаточны ли отклик датчика, диапазон, монтаж, калибровка и синхронизация часов?
  • Сохраняет ли предлагаемое исправление требуемое безопасное поведение при остановке, выпуске, изоляции и повторном запуске?
  • Проверены ли номиналы труб, клапанов, фитингов, сосудов и приводов по подтвержденному результату переходного расчета и применимым правилам?
  • Проведена ли наладка измененной системы во всем реалистичном диапазоне давлений, температур и циклов?

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

Заключение

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

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

Часто задаваемые вопросы о пневматическом ударе в системах

Является ли пневматический удар тем же, что и гидравлический удар в жидкости?

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

Можно ли по уравнению Жуковского рассчитать максимальное давление пневматического удара?

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

С какой скоростью распространяется волна давления при пневматическом ударе?

Для небольшого возмущения в идеальном сухом воздухе при температуре около 20°C первая оценка составляет около 343 m/s. Фактические распространение и затухание зависят от состояния газа, размера возмущения, податливости линии, ограничений и геометрии. Это значение предсказывает первый приход, а не полное заполнение или пиковое давление.

Предотвратит ли клапан плавного запуска пневматический удар?

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

Что следует измерять при испытании на пневматический удар?

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

Источники и технические справочные материалы