Раздел 6. Инженерная экология и смежные вопросы. Кинетостатический расчет шестизвенного рычажного механизма
аналитическим методом
к.т.н. доц. Иванов В. А.
Университет машиностроения 8(499)267-12-00, tmirl94Kajmail.ru
Аннотация. В статье рассмотрена возможность аналитического определения реакций связей многозвенных рычажных механизмов 2-го класса 2-го порядка на основе метода кинетостатики. Выведены аналитические зависимости для определения реакций кинематических пар кулисно-рычажного 6-тизвенного механизма. Приведены примеры графиков, полученных с использованием расчетной математической программы МаЛСас! 14.
Ключевые слова: аналитический расчет, реакции кинематических пар, кинетостатика, многозвенный рычажный механизм.
Современный уровень развития машиностроения требует высокой точности и большого объема выполняемых расчетов. В настоящее время в кинетостатическом расчете применяют, как правило, графоаналитический метод, который, наряду с наглядностью, доступностью, простотой, имеют существенный недостаток: он становится затруднительным, если требуется провести большой объем однообразных построений и вычислений. В настоящей работе приведены аналитические зависимости для определения реакций, возникающих в кинематических парах шестизвенного кулисно-рычажного механизма 2-го класса с использованием метода кинетостатики, показан порядок и последовательность выполнения расчетов, дается анализ знаком членов, входящих в расчетные зависимости. Вычисления удобно выполнять с помощью современных прикладных расчетных программ для ЭВМ, например, МаЛСас!. Ранее в работах [1,2] были выведены аналогичные аналитические зависимости для расчета силовых характеристик рычажных механизмов насоса и пресса, использованные студентами при выполнении курсовых работ с применением расчетных математических программ МаШСас! 14.
Рассмотрим схему шестизвенного рычажного механизма, имеющего два выходных зве-
Рисунок 1 - Схема механизма
При выполнении силового расчета важно правильно выбрать направление вращения кривошипа.
Преодоление силы полезного сопротивления происходит только на прямом (рабочем) ходу выходного звена, в промежутке между двумя крайними положениями механизма, ограничивающими траекторию точки звена, совершающего возвратные движения (в нашем случае движение точки Б звена 6, а также точки С звена 4).
Примем движение штока 6 на рабочем ходу в направлении из крайнего нижнего в
крайнее верхнее его положение, а движение ползуна 4 на рабочем ходу - из крайнего левого в крайнее правое его положение.
Для обеспечения наибольшей производительности необходимо, чтобы время рабочего хода составляло большую долю в полном времени кинематического цикла механизма, чем время обратного хода, т.е. угол фрх поворота кривошипа на рабочем участке движения выходного звена должен быть больше угла фхх поворота кривошипа при обратном движении выходного звена.
Для ползуна 4 , как видно из рисунка 1, углы поворота кривошипа на рабочем и на холостом ходу одинаковы, так как крайними в кривошипно-ползунном механизме ABC будут положения при ф=0° и при ф=180°. Для штока 6 крайние положения могут быть установлены только после анализа его скоростных характеристик. В крайних положениях механизма ABCDE скорость выходного звена 6 равна нулю. Это дает нам возможность, например, по предварительно построенному графику Уб(ф) определить крайние положения кривошипа, а по ним и углы поворота его на прямом и холостом ходу штока 6. Эти углы не одинаковы. Направление вращения надо выбрать так, чтобы при повороте кривошипа из одного крайнего положения в другое на больший угол происходило движение штока 6 из крайнего нижнего в крайнее верхнее положение. При этом одновременно должно быть обеспечено частичное перекрытие рабочего хода ползуна 4.
Очевидно, что в рассматриваемом механизме вращение кривошипа будет происходить по часовой стрелке.
Для определения реакций, возникающих в кинематических парах, воспользуемся методом кинетостатики. Для этого рассмотрим равновесие структурных групп, отсоединив их предварительно от механизма, т.е. нарушив связи, которые звенья группы образуют с другими звеньями механизма, и для каждой группы составим уравнения равновесия - кинетостатики.
Силовой расчет структурной группы 5-6
Структурная группа 5-6 (рисунок 2) включает в себя два шарнирно-связанных звена -камень 5 и шток 6. Эти звенья в механизме образуют поступательные пары: соответственно D с шатуном и Е с неподвижной направляющей. Элементы поступательных пар на схеме группы показаны пунктирными линиями.
Реакция в поступательной паре имеет два неизвестных параметра: величину и точку приложения. Известна линия ее действия - по нормали к поверхности соприкосновения звеньев. Приложим условно реакцию r35 (силу со стороны шатуна 3 на камень 5) на расстоянии h5 0T центра шарнира D перпендикулярно оси шатуна 3, а реакцию R06 (силу со стороны направляющей на шток 6) - на расстоянии Ьб от центра шарнира D перпендикулярно оси штока 6. Эти расстояния являются неизвестными параметрами, характеризующими положение точек приложения сил, поэтому и силы r35 и r06 предварительно показаны пунктиром. Следовательно, необходимо составить четыре уравнения равновесия.
Расчет всегда начинают с составления уравнений моментов относительно точки D внутреннего шарнира для каждого звена отдельно.
Уравнение моментов для звена 5:
ZMd(5) = 0, R35-h5 = 0. (1)
Так как r35 в общем случае не равна нулю, то hs = 0, т е- R35 проходит через точку D. Покажем её сплошной линией.
Уравнение моментов для звена 6:
ZMD(6) = 0, R06-h6 = 0. (2)
Так как r06 в общем случае не равна нулю, то Ьб = 0,т.е. r06 проходит через точку D. Покажем её сплошной линией.
YFC6 Ro:6 и
Рисунок 2 - Силовой расчет структурной группы 5-6
Точки приложения реакций поступательных пар найдены. Теперь определяем их величины. Для этого составляем уравнения проекций сил, приложенных к звеньям группы, на оси прямоугольной системы координат, начало которой разместим в точке D. Положительные абсциссы будем откладывать вправо, а положительные ординаты - вверх.
Z Пр. Fi (X) = 0, -R06 + R35-sinB = 0. (3)
x Пр. Fi (Y) = 0, -Fc6-G6 + Fh6 + R35-cosb = 0. (4)
Здесь Fc6 ~~ сила технологического сопротивления штоку 6; G6 - сила тяжести звена 6, Fh6 ~ сила инерции звена 6, в - угол наклона шатуна 3.
Масса камня 5 не учитывается, поэтому Gs = О и Fms = О-Решая совместно уравнения проекций сил, получим:
r35 = (fee + g6" FH6)/cosb , (5)
R06 = (Fc6 + G6 ~ Риб) ' sine/cOSB . (6)
Fee и G6 всегда направлены вниз по оси Y. При выводе формул было принято это направление, поэтому знак перед ними меняться не будет. Fh6 может изменить свое направление. При выводе формул было принято её положительное направление (ускорение звена 6 отрицательно). Значит, знак перед Fh6 изменится на отрицательный, если ускорение звена 6 имеет положительные значения. Кроме того, надо помнить, что Fc6 ~~ сила технологических сопротивлений механизма действует лишь на рабочем ходу звена 6, т.е. при yD < 0 Fc6 =0-
Для определения реакции во внутреннем шарнире D достаточно рассмотреть равновесие камня 5. Тогда R65 = -R35.
Силовой расчет структурной группы 3-4 Структурная группа 3-4 (рисунок 3) включает в себя два звена, связанных шарниром С: шатун 3 и ползун 4. Внешними парами являются вращательная пара В и поступательная пара С. Расчет реакций связей этой группы начинаем с уравнений моментов, действующих на каждое звено отдельно, относительно точки С внутреннего шарнира. Предварительно разложим полную реакцию шарнира В r23 (усилие со стороны кривошипа 2 на шатун 3) на две
составляющие: нормальную Rj3, направленную по оси шатуна 3, и тангенциальную Ri,3, направленную перпендикулярно оси шатуна 3.
Уравнение моментов, действующих на звено 3: ZMc(3) = 0, или
~ R23 • 1вс + R53 • icd + Миз + fh3 • Ьз + Fm • Ьз + G3 • Ьз = 0 • Из уравнения моментов получаем:
Rf23 = (R53 • Icd + Миз + F-из • h3r + Fh3 • hf + G3 ■ h3B)/lBc • (7)
Здесь 1BC - длина шатуна 3; r53 - усилие, передаваемое со стороны камня 5 на шатун
3, в соответствии с третьим законом Ньютона r53 = - r35; Icd = V(Xc~Xd)2 + Yd - плечо силы r53; mh3 = is3'e3 - момент сил инерции; is3 - момент инерции звена 3; Риз = тз'а;зз ~~ проекция силы инерции звена 3 на ось Y; F^3 = тз' a"s3 - проекция силы инерции звена 3 на ось X; hf = lcs3' sinB - плечо силы ; h3 = lcs3' C0SB - плечо силы и силы тяжести G3; lcs3 - часть длины шатуна, определяющая положение его центра тяжести.
Рисунок 3 - Силовой расчет структурной группы 3-4
Формула (7) составлена с учетом направлений сил, показанных на рисунке 3. Большинство сил, входящих в равенство (7), меняют свое направление за цикл работы механизма.
Усилие R53 со стороны камня 5 на шатун 3 неизменно по направлению, если только на холостом ходу (при у0 < 0) не окажется, что Биб > Об (при рИ6 > - Об)
Направление Миз зависит от £3. Из анализа графика £3 (ф) кривошипно-ползунного механизма [3] следует, что при 0<ср<180° £3 направлено против часовой стрелки, тогда Миз - п0 часовой стрелке. При 180° < ф <360° £3 направлено по часовой стрелке, а, значит, Миз - против часовой стрелки. Таким образом, при расчетах по формуле (7) надо учитывать условие: при 0<ф<180° знак перед Миз отрицательный, при 180°<ф<360° знак перед Миз положительный.
Знаки перед Биз и Риз устанавливают следующим образом.
Знак перед Биз положительный, если а^з >(-) и Р>(-) или «Э^зи
Знак перед Биз отрицательный, если а^з' а или а^з <(-)' а Р>(-)-
Знак перед Биз положителен, если а^з >(-) на всем диапазоне изменения р.
Знак перед Биз отрицателен при а^з <(-)-
Направление Сз неизменно за весь цикл работы механизма.
Уравнение моментов, действующих на ползун 4, можно не составлять, так как анализ сил, приложенных к ползуну 4, с учетом ранее выполненных пояснений для группы 5-6, позволяет сделать вывод о том, что реакция R04 со стороны направляющей на ползун 4 будет проходить через точку С.
Переходим к определению величин нормальной составляющей R23 реакции шарнира В и реакции в поступательной паре R04 . Для этого составляем уравнения проекций сил, действующих на звенья группы 3-4, на оси прямоугольной системы координат. Положительное 256 Известия МГТУ «МАМИ» № 2(14), 2012, т. 4
Раздел 6. Инженерная экология и смежные вопросы, направление оси X - вправо, положительное направление оси Y - вверх.
Уравнение проекций сил на ось X: r23 • cosb + r23 • sine - r53 • sine - Fm - Fh4 ~ Fc4 = 0-
Отсюда:
R23 =(R53-sinB-Rt23-sinB +Fh3 + Fh4 + Fc4)/cosb . (8)
Уравнение проекций сил на ось Y:
- r23 • sine + r123 • cosb - r53 • cosb - Fh3 " G3" G4 + Ro4 = 0.
Отсюда:
R04 = R53' cosb + r23 • sine - r23 • cosb + Fro + G3 + G4 • (9)
В формулах (8) и (9) Fh4 ~ сила инерции ползуна 4, a Fc4 ~ сила технологических сопротивлений ползуну 4.
Fh4 = Ш4 • as4 > где as4 = ас •
В зависимостях (8) и (9) сумма проекций сил алгебраическая. Анализ знака перед r53
был дан ранее. Сила r23 имеет знак, полученный в результате расчета по формуле (7). На рабочем ходу (при yD > 0) знак перед r23 всегда положителен.
Знак перед Fh3 положительный, если а^ >0. Если а^з т0 знак перед Fh3 изменится на отрицательный.
Знак перед Fm положительный, если а<зз >(-)- Отрицательный знак будет при а<зз <(-)-
Знак перед Fh4 положительный, если as4>0-
Знак перед G3 и G4 не изменяется.
Сила Fc4 неизменна по направлению. Направление ее действия устанавливается по заданному графику, но всегда противоположно направлению движения ползуна 4.
Реакция r23 входит в равенство (9) с тем знаком, который она получит после подсчета по уравнению (8).
Полную реакцию в шарнире В рассчитываем по теореме Пифагора:
R23 = A/(R23)2 + (R23)2-
Для определения реакции r34 в шарнире С необходимо разорвать связь ползуна 4 и шатуна 3 и рассмотреть равновесие ползуна.
Векторное уравнение сил имеет вид:
R04 + G4 + Fc4 + FH4 + R34 = т.е. решение его, как это видно из рисунка 3, будет представлять собой прямоугольный треугольник, в котором r34 является гипотенузой.
Тогда по теореме Пифагора найдем
R34 = V(Ro4+G4)2 + (FH4-FC4)2 • (10)
Силовой расчет входного звена
Расчетная схема кривошипа представлена на рисунке 4.
Неизвестными факторами являются реакция в шарнире A R02 (сила, действующая со стороны стойки на кривошип) и уравновешивающий момент Мур, который нужно приложить к входному звену для уравновешивания всех сил, действующих на звенья механизма.
Составляем уравнение моментов относительно точки А.
R32 • г • sin(<3 - в) + R^ • г • cos(<? - в) + G2 • 1AS2 • cos (Р -Мур =0.
Отсюда:
Myp = R?2-r-sin(^-B)+Ri2-r-cos(^-B) + G2TAS2-cos^ (11)
Л 02
____« X
Рисунок 4 - Силовой расчет входного звена
Здесь Яз2- составляющая реакции со стороны шатуна на кривошип, принятая для удобства расчетов равной - К23 ; К32 - составляющая реакции со стороны шатуна на кривошип, принятая равной - я23; г - длина кривошипа; 1А82 - расстояние до центра масс кривошипа; Ог ~~ сила тяжести кривошипа.
Остальные составляющие равенства показаны на рисунке 4.
Переходим к определению реакции в шарнире А Я02 . Реакция имеет два неизвестных параметра: величину и направление.
Представим ее в виде двух составляющих и Яд2 - проекций Я02 на оси X и У прямоугольной системы координат с началом координат в точке А, положительным направлением оси X вправо, а оси У - вверх.
Уравнение проекций сил на ось X: Яо2 - • совв - я32 • вшв - рИ2 • со$<р = 0, откуда
Ы02 =Ы32 ' С08В + Р-32 ' 81ПВ + РИ2 ' С0$<р. (12)
Уравнение проекций сил на ось У:
Я02 + Я32 ' вШВ " Я32 ' С08В + Ри2 ' ЫПф " С2 = 0,
откуда:
1^-02 =Яз2'С08в-к-З2 ' ^ " ри2 • втр + С2
(13)
^зз(ф) 2.75« 10
4001
240 — "7 80-/-
-80«—< -240 —
МурСФ)
-400 -560 -720-
-т-
- 1.04x10-~ 1.2к!0:
г8—
к
м-
Рисунок 5 - График изменения реакции в паре шатун (звено 3) - ползун (звено 5)
Рисунок 6 - График изменения уравновешивающего момента на кривошипе
Раздел 6. Инженерная экология и смежные вопросы. Здесь рИ2 = тг • и! • 1а82 ~ сила инерции звена 2. Полная реакция шарнира А
Я02 = д/^Мког)2
Результаты расчетов могут быть представлены в виде таблиц, графиков или годографов усилий в кинематических парах. Примеры графиков представлены на рисунках 5 и 6.
Проверка с помощью графо-аналитического метода показала хорошее совпадение результатов.
Литература
1. Иванов В.А. Кинематический и силовой расчет рычажного механизма насоса аналитическим методом: Методические указания. - М.: МГУИЭ, 2003. - 28 с.
2. Иванов В.А. Расчет шестизвенного механизма пресса аналитическим методом: Методические указания. - М.: МГУИЭ, 2008. - 24 с.
3. Артоболевский И.И. Теория механизмов и машин: Учебник для втузов. - М.: Наука, 1988. - 640 с.
Экспертная система по выбору технологического процесса переработки
природного газа
д.т.н. проф. Кольцова Э.М., к.т.н. Филиппова Е.Б., к.т.н. доц. Зубов Д.В.
РХТУ им. Д. И. Менделеева, Университет машиностроения
Аннотация. Разработана экспертная система, позволяющая получить информацию о полезности интересующего месторождения, сделать выбор оптимального способа переработки природного газа и утилизации попутного нефтяного газа непосредственно в местах добычи, что является весьма актуальным.
Ключевые слова: экспертная система, база данных, переработка природного газа
Как известно, Россия занимает первое место в мире по запасам углеводородного сырья и объемам его добычи. При этом следует отметить тот факт, что более 20 % от всего объема извлекаемых углеводородных газов сжигается сейчас на факеле, либо, что еще хуже, просто выбрасывается в атмосферу.
Начисная с 2005 года правительство РФ приняло ряд законопроектов, обязывающих недропользователей довести степень использования добываемого природного газа до 95 % , т.е. сократить потери до 5%.
Однако удаленность большинства месторождений от газоприемной системы «Газпрома» и несоответствие добываемого природного газа требованиям ряда ГОСТов на его качество не дают возможности непосредственно подавать этот газ в транспортную систему. Поэтому в газовой отрасли возникла задача переработки углеводородных газов на самом месторождении в полезные, а главное, легко транспортируемые продукты, что является наиболее экономически выгодным решением.
В такой ситуации создание экспертной системы (ЭС), позволяющей быстро сделать выбор оптимального способа переработки природного газа (Ш) и утилизации попутного нефтяного газа (ПНГ) на месторождениях, непосредственно в местах добычи, является весьма актуальным.
Такая экспертная система создавалась на основе обширной информационной системы, включающей ряд баз данных. Разработка ЭС состояла из следующих этапов:
• сбор, анализ и систематизация информации о месторождениях России и СНГ;
• сбор, анализ и систематизация информации о современных технологиях переработки ПГ и утилизации ПНГ;