Моделирование динамических процессов в пневмоцилиндре
Моделирование динамических процессов в пневмоцилиндре
Министерство образования и науки Украины
Донецкий национальный технический университет
Немецкий технический факультет
МехатроникаБАКАЛАВРСКАЯ РАБОТАна тему:«Моделирование динамических процессов впневмоцилиндре»Выполнил Скорик Т.Ю.Руководитель к.т.н. Деркач А.В.Донецк 2006СодержаниеВведение1. Характеристика объекта проектирования и постановка задачи исследования1.1 Схема и циклограмма дифференциального привода 1.2 Процесс наполнения сжатым воздухом рабочей полости привода1.3 Процесс истечения сжатого воздуха из выхлопной полости привода1.4 Динамический расчет дифференциального привода2. Разработка математической модели объекта3. Создание модели внешних нагрузок4. Передаточная функцияЗаключениеСписок литературыВведениеПневматические системы получают все большее распространение в промышленности как наиболее надежные средства автоматизации. Пневмоустройства используют в качестве приводов зажимных и транспортирующих механизмов, для дистанционного управления и регулирования, в контрольно-измерительных приборах, при автоматизации машин и устройств, работающих в агрессивных средах, в условиях пожаро- и взрывоопасности, радиации, а также при значительной вибрации и высоких температурах и т.д. Пневмосистемы распространены в автомобильной промышленности, в самолетостроении, в космонавтике, где они применяются для автоматизации сборочных работ, для управления аварийными системами и т.д. Пневмоустройства используют для управления также в нефтяной, газовой, химической, пищевой промышленности, в горном деле, в строительстве и т.д. Элементы пневмоавтоматики все больше внедряются в медицинские приборы различного назначения (для искусственного дыхания, кровообращения, инъекций и т.д.)Повышение производительности машин применительно к пневматическим системам означает увеличение их быстродействия. Так как пневматические системы являются частью машины, то время их срабатывания входит в общее время рабочего цикла машины и влияет таким образом на производительность машины в целом. Проблема повышения производительности непосредственно связана с разработкой методов динамического исследования машин и выбора их оптимальных параметров и структуры. Эти проблемы составляют основное направление как общей теории машин, так и теории пневматических систем машин. Теория пневматических систем имеет свои характерные черты и особенности. Если в теории механизмов с твердыми звеньями задача динамики состоит в решении уравнения движения последних, то в теории пневматических систем это уравнение всегда решается совместно с уравнениями, характеризующими термодинамические процессы, протекающие в полостях рабочих цилиндров. При составлении этих уравнений используются законы термо- и газодинамики. Так как полости большинства пневмотических устройств постоянно соединены либо с магистралью, либо с атмосферой, то при составлении расчетных уравнений используются методы термодинамики переменной массы газа. Таким образом, в теории пневматичесикх систем кроме методов механики твердого тела применяются методы механики упругой жидкости. В машиностроении получили широкое распространение пневматические системы разных типов. Их можно разделить на следующие основные группы: пневмоприводы, пневматические системы управления и пневматические устройства. Широкое применение пневмоприводов и систем управления объясняется их преимуществами по сравнению с другими средствами автоматизации, в первую очередь надежностью функционирования, которая в современных автоматизированных системах управления играет важную роль. Преимуществом пневмоустройств по сравнению с электрическими исполнительными устройствами является возможность воспроизведения поступательного движения без каких-либо передаточных механизмов. Благодаря этому они получили широкое распространение в тех случаях, когда требуется осуществить возвратно-поступательное движение. Пневмоустройства вращательного движения отличаются от электромоторов меньшими габаритами, нечувствительностью к длительным перегрузкам, простотой регулирования скорости вращения и крутящего момента, полной безопасностью для оператора, но их работа сопровождается большим шумом. По сравнению с гидравлическими пневматические приводы обладают такими преимуществами: их исполнительные устройства имеют большие скорости срабатывания и более низкую стоимость, возвратные линии значительно короче, так как воздух может быть удален в атмосферу из любой точки системы; наличие неограниченного запаса воздуха в качестве рабочего тела также способствует широкому распространению пневмоустройств. Вместе с тем пневматические приводы при равных габаритах с гидравлическими развивают меньшие усилия, что объясняется более высоким давлением жидкости в последних. Пневмоустройства следует применять в тех случаях, когда требуется обеспечить высокие скорости движения рабочего органа при относительно небольших рабочих усилиях. Неизбежные утечки воздуха из системы значительно понижают к.п.д. пневмоустройств. 1. Характеристика объекта проектирования и постановка задачи исследования.1.1 Схема и циклограмма дифференциального приводаСхема одностороннего привода с дифференциальным поршнем 1, площади которого с обеих сторон значительно отличаются друг от друга, приведена на рис. 1. Штоковая полость привода постоянно соединена с магистралью (со сжатым воздухом давлением рм), а поршневая в исходном состоянии соединена с атмосферой (ра), благодаря чему создается начальный перепад давления сжатого воздуха, и поршень 1 занимает крайнее левое положение.После срабатывания управляющего устройства 2 распределитель 3 сообщает с магистралью поршневую полость, давление воздуха в ней повышается, и когда сила давления со стороны этой полости превысит силу сопротивления привода, поршень начинает двигаться вправо (положение, показанное на рисунке). Как только поршень совершит свой рабочий ход s, управляющее устройство переключит распределитель 3, который соединит поршневую полость с атмосферой для того, чтобы поршень 1 мог вернуться в исходное состояние. Дроссели 4 и 5 служат для регулирования скорости движения поршня. Другие обозначения на чертеже будут объяснены позднее.Рис. 1. Схема дифференциального приводаОдносторонние приводы с начальным перепадом давления называют также дифференциальными. Рассмотрим последовательность действия этого привода, т. е. его циклограмму. В циклограмму пневмопривода кроме интервалов времени перемещения и остановки поршня для наглядности включены также интервалы времени изменения давления рl и р2 В полостях рабочего цилиндра в виде диаграммы (рис. 2).Началом циклограммы является момент включения управляющего устройства, время срабатывания которого относится ко времени срабатывания распределителя (см. рис. 1). Процессы открытия отверстий для сжатого воздуха в управляющем устройстве и распределителе могут происходить одновременно с процессом распространения волны сжатого воздуха и заканчиваться в разные моменты времени.В первом приближении предположим, что волна давления воздуха возникает после полного открытия отверстия распределителя, так как у большинства пневмоприводов время открытия распределителя невелико по сравнению с временем всего рабочего цикла. Рис. 2. Циклограмма дифференциального приводаУказанное допущение позволяет отдельно определять интервалы времени t1 (время открытия распределителя) и t2 (время распространения волны давления от распределителя до цилиндра, см. рис. 2). Горизонтальные прямые на верхней диаграмме соответствуют интервалам времени выстоев поршня, а наклонные - интервалам времени его движения. Кривые на нижней диаграмме отражают процессы изменения давлений рl и р2 воздуха соответственно в поршневой и штоковой полостях. Давление в первой из них начинает расти после открытия распределителя и этот процесс продолжается до начала движения поршня (давление рlд, интервал t3). Сумма интервалов t1, t2 и t3 составляет подготовительное время tІ. В период движения поршня tІІ давление в этой полости можетмонотонно увеличиваться (уменьшаться) или колебаться в зависимости от соотношения конструктивных параметров привода и нагрузки. В штоковой полости давление воздуха, равное магистральномy в исходном положении, при движении поршня несколько возрастает, в результате чего часть воздуха перетекает в магистраль, последний процесс может закончиться уже в период остановки поршня (штриховая линия). В поршневой полости после того, как поршень закончит рабочий ход, давление рl обычно уравнивается с магистральными рм или с давлением, требуемым технологическим процессом, которому соответствует время tтехн.После выполнения заданной технологической операции управляющее устройство снова переключается (время tтехн технологической "операции не рассматривается). Затем в той же последовательности начинается обратный ход поршня, при этом вначале давление падает в обеих полостях: в поршневой - за счет соединения с атмосферой, а в штоковой - за счет увеличения ее объема при движении поршня влево. После возвращения поршня в исходное положение давление в поршневой полости падает до атмосферного, а в штоковой - поднимается до магистрального.1.2 Процесс наполнения сжатым воздухом рабочей полости приводаНаполнение сжатым воздухом рабочей полости двигателя объемом V1 (см. рис. 1) происходит из магистрали, баллона или какого-либо другого источника питания, давление в котором принимаем постоянным (рм = const), имея в виду, что на входе в привод находится стабилизатор давления. Термодинамические процессы в пневматических приводах в первом приближении рассматриваем как квазистационарные и протекающие при установившихся режимах истечения. При поступлении сжатого воздуха в полость рабочего цилиндра или при истечении из нее давление в различных точках объема будет неодинаковым. Сначала меняется давление вблизи входного (или выходного отверстия), затем изменение давления постепенно распространяется на весь объем.Учитывая, что выравнивание параметров воздуха происходит достаточно быстро по сравнению с рассматриваемым процессом и большой разницы в их значениях для разных точек объема не наблюдается, в теории пневматических систем процессом выравнивания пренебрегают. Все процессы рассматривают как квазистационарные, т.е. такие, при которых во всех точках объема полости предполагаются одинаковые параметры (давление, температура и плотность). Принимаем, что рабочая полость имеет индекс 1.Считая согласно первому закону термодинамики, что вся подведенная с газом тепловая энергия dQм расходуется на изменение внутренней энергии dU1 и на работу расширения газа dL1, запишем уравнение энергетического балансаdQм = dU1 + dL1 (1.1)Имея в виду, что количество тепловой энергии, поступившей в полость с газом, равно произведению его массы mм на удельную энтальпию (dQM = iMdmM), а внутренняя энергия U1 газа и совершаемая им работа L1 равны соответственно dU1 = d (u lml) и dL1 = р1dV1, представим уравнение (1.1) в следующем виде:iM dmM = u ldml + ml dUl + рl dV1, (1.2) где u l - удельная внутренняя энергия.Выразим в уравнении (1.2) значения энтальпии и внутренней энергии через произведение температуры на теплоемкость соответственно при постоянных давлении ср и объеме сv:сpTMdmM = сvTldml +cvmldTl + рldVl (1.3)Рассматривая воздух как идеальный газ, молекулярными силами сцепления которого можно пренебречь, опишем его состояние с помощью уравнения Клапейрона:рl V1= m1RT1, (1.4)где R - газовая постоянная, R = 287 Дж/(кг, К) для воздуха (при Тм = 290 К).Подставляя в уравнение (1.3) значение m1 dT1, полученное из уравнения (1.4), и полагая в нем и = R, где k - показатель адиабаты, после несложных преобразований получаем следующее выражение:kRT м dmM = V1 dp1 + kp1 dV1.(1.5)Заменим в уравнении (1.6) массу сжатого воздуха dmM, поступающего в полость V1 в течение времени dt, соответствующим значением GM расхода dmM = GM dt и выразим полученное уравнение относительно давления: (1.6)Расход GM воздуха из неограниченного объема (магистрали) определяют чаще всего по формуле Сен-Венана и Ванцеля:, (1.7)где м1 - коэффициент расхода; f1 - площадь входного отверстия; ТМ - температура воздуха в магистрали.Из решения (1.5) совместно с уравнением (1.7) и определяют давление Р1 в полости как функцию времени. Из формулы (1.7) следует, что расход GM является функцией отношения давлений у, в числителе которого всегда находится давление той среды, куда истекает газ, а в знаменателе - давление среды, откуда вытекает этот газ. Представим формулу (1.7) расхода воздуха из магистрали уМ в более удобном виде:, (1.7)где ; ;.При Т м = 293 К расход GM = 0,00912, причем Рм в Па. Чтобы найти максимум расходной функции ер (а), приравняем нулю ее производную , откуда получим, критическое отношение давлений= 0,5282 (при к = 1,4). (1.9)Подставим в уравнение (1.6) значение расхода GM из выражения (1.8), получим уравнение для определения давления в полости наполнения в общем виде: (1.10)Объем V1 рабочей полости равен произведению площади поршня F1 на перемещение поршня х (с учетом его начальной координаты х01), т.е. . Подставим это значение в (1.10), получим, (1.11)где ,V01 - начальный объем рабочей плости,F1 - площадь поршня.Если это уравнение выразить относительно 1, то после интегрирования найдем время наполнения постоянного объема от начального безразмерного давления у11 до заданного давления у12: (1.12)Значения функций Ш1 (у11) и Ш1 (у12) определяем по графику, приведенному на рис.3.Рисунок 3Формулу (1.13) применяют для определения времени наполнения полости постоянного объема до заданного давления, в том числе и для определения времени 3 (см. циклограмму на рис. 3).1.3 Процесс истечения сжатого воздуха из выхлопной полости приводаПри перемещении поршня 1 (см. рис. 1) в дифференциальном приводе давление сжатого воздуха в выхлопной полости 2 может повышаться вследствие уменьшения ее объема. В этом случае происходит истечение воздуха в магистраль. Здесь также может быть применен первый закон термодинамики (1.1), но в этом уравнении следует поставить знак минус в левой части, так как происходит истечение воздуха-dQ2 = dU2 + dL2.(1.14)Соответственно изменяем индекс 1, относящийся к рабочей полости, на индекс 2 выхлопной полости. Далее, выкладки будут аналогичны приведенным в формулах (1.1) - (1.5). Остановимся на последнем выражении-k•R•Т2 dm2 = V2 dР2 + k•Р2 dV2. (1.15) Имея в виду, чтоdm2 = - d () = ,получаем 2 dР2 + k•Р2 d2 = 0 или после интегрирования и потенцирования этого выражения - P2• = const - уравнение адиабаты. Расход воздуха из ограниченного объема V2 в магистраль описывается также формулой Сен-Венана и Ванцеля, однако в ней следует положить ТМ = Т2, Рм =Р2 имея при этом в виду, что все эти величины являются переменными:, (1.16)где при 0,528 < у < 1.Подставив в уравнении (1.15) dm2 = 2 dt и 2 из (1.16), получим уравнение для определения давления в выхлопной полости, соединенной с магистралью:, (1.17) где - площадь поршня со стороны штоковой полости. Температура Т2 в уравнении (1.17) может быть выражена через давление Р2 на основании уравнения адиабаты: .(1.18)Тогда получим следующее уравнение для определения давления воздуха при истечении его из ограниченного объема:(1.19)При обратном ходе подготовительное время t3' будет характеризовать время истечения полости до необходимого давления, определяемого нагрузкой. Однако при обратном ходе эта полость становится выхлопной, соединенной с атмосферой. Давление будет изменяться по уравнению (1.19), в котором следует вместо 1/у2 подставить уа/у2 так как истечение будет происходить в атмосферу пропорционально отношению давлений Ра/Р2 = уа/у2 где уа = Ра/Рм, у2= Р2/РМ.Так как при обратном ходе поршневая полость становится выхлопной, присвоим ей индекс 2:,(1.20)где при 0 < у < 0,528; при 0,528 < у < 1.Расход воздуха при истечении из ограниченного объема полости в атмосферу определяем по формуле (1.16), в которой принимаем РМ = Ра: , (1.21)где .При Тм = 293 К расход 2 = 0,00912. Для определения подготовительного и заключительного времени в уравнение (1.20) следует подставить х = 0, dx = 0. Если полученное уравнение выразить относительно t, то получим после интегрирования время истечения воздуха из постоянного объема в диапазоне изменений давлений у21 (Р21) до у22 (Р22):. (1.22)Значения Ш2 () и Ш2 (), определяем по графику на рис.3, стр. 11. 1.4 Динамический расчет дифференциального приводаУравнение движения поршня дифференциального привода имеет вид: Р, (1.23)где - масса поршня; Р1 - сила вредного сопротивления (трения);Р2 - сила полезного сопротивленияР - результирующая всех сил, приложенных к поршню;; (1.24); (1.25); (1.26). (1.27)Рассмотрим обратный ход поршня. Уравнение его движения при обратном ходе, когда поршневая полость 2 соединяется с атмосферой имеет вид: Р, (1.28)Р = Р1 + Р2 + Рз + Ра (Р1 - Р2),2. Разработка математической модели объектаНа основании дифференциальных уравнений, которые описывают поведение пневмоцилиндра в процессе работы, была составлена динамическая модель пневмоцилиндра. Далее выполняем моделирование (исследование) составленной модели. Нагружаем модель единичным ступенчатым воздействием, который воздействует на поршневую полость.
Таблица 1 Обозначения переменных, используемых в дин.модели
Описание
Обознач.
В схеме
Ед.
Давление в поршневой полости пневмоцилиндра
P1
P1
Па
Давление в поршневой полости пневмоцилиндра
P2
P2
Па
Начальная координата поршня
Х01
Х01
м
Начальный объем рабочей полости
V01
V01
м3
Площадь поршня
F1
F1
м2
Давление в магистрали
PM
Pm
Па
Газовая постоянная
R
R
Температура воздуха в магистрали
TM
Tm
К
Показатель адиабаты
Коэффициент расхода
My1
Площадь входного отверстия
f1
f1
м2
Функция расхода
Fi1
Коэффициент расхода
My2
Площадь входного отверстия
f2
f2
м2
Площадь поршня со стороны штоковой части
F2
F2
м2
Рабочий ход
s
s
м
Конечная координата поршня
Х02
Х02
м
Перемещение поршня
Х
Х
м
Давление в штоковой полости Р2 находится с помощью давления Р1:
Рисунок 2.1 - Подсистема для уравнения нахождения Р2.
Представим каждое уравнение динамической модели в виде схемы:
Рисунок 2.2 - Подсистема для уравнения .
Рисунок 2.3 - Подсистема для уравнения .
Окончательная схема приведена на рисунке 2.4:
Рисунок 2.4 - Схема, описывающая поведение объекта.
Рисунок 2.5 - График изменения суммарных сил, действующих в поршне, от воздействия силы изменяющейся ступенчато
Перемещение поршня
Рисунок 2.6 - График, показывающий перемещение пневмоцилиндра от воздействия всех факторов
3. Создание модели внешних нагрузок
3.1 Скачкообразно изменяющиеся нагрузки
Воздействие нагрузки, меняющейся скачком, соответствует воздействию на пневмоцилиндр однократного усилия, связанного, например, с наполнением поршневой полости цилиндра сжатым воздухом, что повлечет за собой возрастание давления в рабочей полости на определенную величину.
Для имитации будем использовать функцию Step системы MatLab:
Рис 3.1 Схема изменения силы для функции Step
Произведем тестовый расчет поведения системы. Получим:
Рисунок 3.2 - График изменения давления в поршневой полости, от воздействия силы изменяющейся ступенчато
Рисунок 3.3 - График изменения давления в штоковой полости, от воздействия силы изменяющейся ступенчато
4. Передаточная функция
;
Введем преобразователь Лапласа
Принимаем такие значения: масса поршня 3 кг, коэффициент трения 0,14.
Рисунок 4.1 - График нулей и полюсов системы
Как видим полюса системы находятся на одной линии, что означает, что в системе отсутствуют колебания.
Переходный процесс имеет вид:
Рисунок 4.2- График переходного процесса
Рисунок 4.3 - График АЧХ и ФЧХ характеристик
Импульсное воздействие
Рисунок 4.4- График импульсного воздействия на систему
При воздействии импульсной силы на систему амплитуда меняется без скачков, без колебаний, равномерно. Переходный процесс также меняется без колебаний.
Заключение
В результате выполнения данной работы была составлена расчётная схема, а также составлены дифференциальные уравнения, описывающие динамические процессы в пневмоцилиндре. На основании дифференциальных уравнений была разработана модель пневмоцилиндра.
Получены графики, описывающие изменение перемещения поршня в зависимости от изменения наполнения поршневой полости сжатым воздухом.
Решение вышеперечисленных задач было реализовано посредством программного пакета MatLab.
Пневмоприводом называют систему взаимосвязанных пневматических устройств, предназначенных для приведения в движение одного или нескольких твердых тел, входящих в состав машины или механизма. В состав пневмопривода входят: исполнительные устройства (двигатели), предназначенные для преобразования энергии сжатого воздуха в механическую энергию движения рабочих органов машины, которые выполняют заданную технологическую операцию; распределительные устройства, предназначенные для направления потоков воздуха из магистрали в рабочие цилиндры исполнительных устройств и из рабочих цилиндров в атмосферу, и управляющие устройства, предназначенные дл обеспечения последовательности перемещения рабочих органов машин в соответствии с требуемым законом их движения.
В зависимости от структуры исполнительных устройств различают пневмоприводы поршневые, мембранные, сильфонные, пневмомоторы, комбинированные. Наибольшее распространение получили поршневые и мембранные приводы, большинство их образует группу так называемых типовых приводов, у исполнительного устройства которых имеется не более двух полостей. Приводы поступательного и вращательного движения с большим числом полостей у исполнительных устройств образуют сложные приводы.
Различают также приводы дискретного и непрерывного действия.
Типовые приводы в зависимости от вида входящих в них исполнительных устройств разделяются на односторонние и двусторонние. Односторонним пневмоустройством называют устройство, в одну полость которого в период работы подается (или из которой вытекает) сжатый воздух, при этом вторая полость соединена с полостью постоянного давления (магистралью или атмосферой).
Двусторонним пневмоустройством называют устройство, в обе полости которого попеременно поддается сжаты йвоздух, вследствие чего рабочий орган совершает возвратно-поступательное движение. Как односторонние, так и двусторонние устройства разделяются на устройства с начальным перепадом давления воздуха на поршне и без него.
Основными разделами теории пневматических систем следует считать динамический анализ и синтез приводов, структурный анализ и синтез систем управления. Первый из этих разделов построен на базе механики твердого тела и жидкости, второй - на основе кибернетики.
Список использованной литературы
1. Герц Е.В. Динамика пневматических систем и машин. - М.: Машиностроение, 1985. - 256с., ил.
2. Иринг Ю. Проектирование гидравлических и пневматических систем/Пер. со словац.Д.К.Рапопорта.- Л.: Машиностроение. Ленингр. отд-ние, 1983. - 363 с., ил.