Файл: Сергеев В.И. Исследование динамики плоских механизмов с зазорами.pdf

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

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

По указанным выше причинам рассматриваемую задачу будем решать с применением ЭЦВМ. Точность их работы определяется только величиной разрядной сетки машины и используемым ал­ горитмом решения задачи. Так, арифметические операции выпол­ няются машинами типа «Минск-2» и «Минск-32» с точностью порядка ІО-7. Причем подобная точность решения задачи не изменяется со временем и не зависит от времени решения задачи. Формирова­ ние нелинейных функций производится стандартными подпрограм­ мами с такой же высокой точностью. При использовании стан­ дартных подпрограмм решения систем дифференциальных уравне­ ний можно задавать требуемую точность решения. Большим не­ удобством до последнего времени при использовании цифровой вычислительной техники была сложность подготовки задачи к ее решению на машине, т. е. программирование на языке конкретного, типа машины с последующей перфорацией программы и вводом ее в машину. Однако в настоящее время повсеместно находит приме­ нение автоматизация программирования, заключающаяся в напи­ сании программы на алгоритмическом языке типа АЛГОЛ, ФОРТРАН, АКИ-ИҢЖЕНЕРи других. При этом не требуется знать код операций конкретного типа машины. Именно это позволяет широко использовать ЭЦВМ в инженерной практике.

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

Кроме того, важным достоинством ЭЦВМ является простота реализации большого количества сложных логических операций, что необходимо при решении рассматриваемой задачи при переходе от одного вида движения к другому.

Общий моделирующий алгоритм исследования дополнительного движения механизмов с двумя зазорами

Составим общий моделирующий алгоритм [21] исследования ди­ намики плоских механизмов с двумя зазорами в его кинематических парах. Введем следующие операторы (значения индексов / соот­

30

ветствуют порядковым номерам операторов в приводимой ниже блок-схеме программы):

А1 — оператор подготовки исходной информации, осуществляющий перевод констант из десятичной системы счисления в двоич­ ную;

Q i— операторы печати, производящие печать результатов реше­ ния задачи, параметров механизма и начальных условий;

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

Pj — логические операторы, осуществляющие проверку выполне­ ния условий перехода от одного вида движения к другому,

атакже момента окончания решения;

А27— арифметический оператор, вычисляющий некоторую функ­ цию от координат, которая служит для определения наступ­ ления периодического движения механизма с зазором;

Я29— оператор окончания решения.

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

и0, если оно не выполнено.

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

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

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

AlQ* 2- 8’ 13’ 18' 26- 2B/ sQ4P p 7Pel; P8U8' 1319' 26/„Qio

31

Логические операторы P5, Pn , P17, P23 сравнивают текущее значение угла поворота кривошипа а с некоторым его значением а*, при котором возможно наступление периодического движения механизма с зазорами. Операторы Р„ и Р12 проверяют условие R1 0, т. е. регистрируют момент нарушения кинематической связи в паре 1—2, а операторы Р 7, Рв и Р2І проверяют условие R3 > 0. Операторы.Р13, Р25 и Р2а регистрируют момент окончания

Блок-схема программы, составленная в соответствии с написан­ ным алгоритмом, приведена на рис. 2. Оператор Аг переводит заданные параметры исследуемого механизма и начальных условий из десятичной системы счисления в двоичную и передает управле­ ние оператору Q2, который производит печать параметров и началь­

значения старших производных по заданным параметрам механиз­ ма и начальным условиям и интегрирует уравнения движения од­ ним из численных методов, например методом Рунге — Кутта [46]. Как правило, в библиотеке стандартных программ вычислительной машины имеется большое число программ, составленных в соответ­ ствии с методами численного интегрирования обыкновенных диф­ ференциальных уравнений [55]. После интегрирования на первом шаге по значениям полученных координат ух и у3 и их производных Ті и уд определяются величины реакций Рг и Р 3, производится печать этих значений, а также координат ул., Узи их производных или же эти значения передаются в запоминающее устройство для форми­ рования массивов данных, необходимых для построения соответ­ ствующих графиков, затем управление передается оператору Р 5, где проверяется условие а = а*.

а* выбирается равным а„ + 2л/г, т. е. такому значению угла по­

Оба выхода логического оператора Р„ подаются на входы логиче­ ских операторов Р 7 и Р3, проверяющих одно и то же условие R 3 0.

Если окажется, что R1> 0 и R3 > 0, то управление передается снова оператору / 3 и продолжается решение уравнений движения,

32

С4

Рис.

es

о

Подготовка исходной

2 В. И. Сергеев, К. М. Юдин

соответствующих движению механизма с контактом в обоих зазо­ рах на следующем шаге, и т. д. Если же окажется, что Rx ^ О, а Rs !> О, т. е. в паре 1—2 произошел разрыв контакта при сохра­ нении последнего в паре 2—3, то управление передается оператору / 21. Этот оператор решает уравнения движения рассматриваемого механизма при сохранении контакта в паре 2—3 и разрыве в паре

начальные условия и старшие производные и интегрирует получен­

После первого шага интегрирования уравнений движения управ­ ление передается оператору Ргз для проверки условия a = а*.

раняется контакт в паре 2—3 и сохраняются условия бесконтакт­ ного движения в паре 1—2, то управление снова передается опе­ ратору / 21 для дальнейшего интегрирования уравнений движения

IIвида.

Аналогично работают другие части блок-схемы программы.

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

Если же на каком-то этапе решения задачи окажется выполнен­ ным условие а = а*, то управление будет передано арифметиче­ скому оператору Л27 для вычисления некоторой функции | 2 ] от текущих значений координат хг, х3, уг, у 3 или уь у3 и их производ­ ных, значение которой позволяет судить о наступлении периоди­ ческого режима движения.

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

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

34