Передаточные функции группы Асcура второго вида 


Мы поможем в написании ваших работ!



ЗНАЕТЕ ЛИ ВЫ?

Передаточные функции группы Асcура второго вида



Рисунок 5 – Группа Ассура второго вида

Assur2 (Xa, Ya, Xa_1, Ya_1, Xa_2, Ya_2, { входные параметры }

Xn, Yn, Xn_1, Yn_1, Xn_2, Yn_2,

Fing, Fin_1, Fin_2, l, LAB, h, lcc, j, q, Mas: Real;

var FiABg, FiAB_1, FiAB_2: Real), { возвращаемые параметры }

где Xa, Ya, Xa_1, Ya_1, Xa_2, Ya_2 - значения ПФ0, ПФ1, ПФ2 оси вращательной кинематической пары А;

Xn, Yn, Xn_1, Yn_1, Xn_2, Yn_2 - значения ПФ0, ПФ1, ПФ2 точки N, определяющей положение направляющей, по которой перемещается ползун;

Fing, Fin_1, Fin_2 - угловые ПФ0, ПФ1 и ПФ2 направляющей, по которой перемещается ползун;

l - длина участка направляющей, изображаемого на экране, мм;

LAB - длина шатуна, мм;

h - расстояние от оси вращательной КП до оси направляющей, мм;

j - признак сборки, которому следует присвоить значение +1 или -1;

FiABg, FiAB_1, FiAB_2 -угловые ПФ0, ПФ1, ПФ2 шатуна.


Передаточные функции группы Асcура третьего вида

 

Рисунок 6 – Группа Ассура третьего вида

Assur3 (Xa, Ya, Xa_1, Ya_1, Xa_2, Ya_2, { входные параметры }

Xb, Yb, Xb_1, Yb_1, Xb_2, Yb_2,

h1, h2, Mas, lcc, q: Real;

var Fikg, Fik_1, Fik_2: Real), { возвращаемые параметры }

где Xa, Ya, Xa_1, Ya_1, Xa_2, Ya_2,

Xb, Yb, Xb_1, Yb_1, Xb_2, Yb_2 - значения ПФ0, ПФ1, ПФ2 осей внешних вращательных кинематических пар А и В;

h1, h2 - расстояния осей вращательных КП до продольной оси кулисы;

lcc - длина диагонали камня изображаемого на экране;

Fikg, Fik_1, Fik_2 -угловые ПФ0, ПФ1 ПФ2 продольной оси кулисы.


Передаточные функции группы Асcура четвертого вида

 

Рисунок 7 – Группа Ассура четвертого вида

Assur4 (Xn, Yn, Xn_1, Yn_1, Xn_2, Yn_2, { входные параметры }

Fing, Fin_1, Fin_2, hn,

Xm, Ym, Xm_1, Ym_1, Xm_2, Ym_2,

Fimg, Fim_1, Fim_2, hm, Mas, lcc, q: Real;

var Xk, Yk, Xk_1, Yk_1, Xk_2, Yk_2:Real),{ возвращаемые параметры}

 

где Xn, Yn, Xn_1, Yn_1, Xn_2, Yn_2,

Xm, Ym, Xm_1, Ym_1, Xm_2, Ym_2 - значения ПФ точек N и M, определяющих положение направляющих;

Fing, Fin_1, Fin_2 - угловые ПФ продольной направляющей n;

Fimg, Fim_1, Fim_2 - угловые ПФ продольной направляющей m;

hn, hm - расстояния от оси внутренней КП до осей направляющих;

Xk, Yk, Xk_1, Yk_1, Xk_2, Yk_2 - значения ПФ внутренней вращательной кинематической пары.


Передаточные функции группы Асcура пятого вида

 

Рисунок 8 – Группа Ассура пятого вида

Assur5 (Xa, Ya, Xa_1, Ya_1, Xa_2, Ya_2, h, { входные параметры }

Xn, Yn, Xn_1, Yn_1, Xn_2, Yn_2,

Fing, Fin_1, Fin_2, Alfag, Mas, lcc, q: Real;

var Xk, Yk, Xk_1, Yk_1, Xk_2, Yk_2: Real), {возвращаемые параметры}

где Xa, Ya, Xa_1, Ya_1, Xa_2, Ya_2 - значения ПФ оси вращательной КП;

h - расстояние от оси вращательной КП до оси внутренней направляющей;

Xn, Yn, Xn_1, Yn_1, Xn_2, Yn_2 - значения ПФ точки N, определяющей положение внешней направляющей;

Fing, Fin_1, Fin_2 - угловые ПФ оси внешней направляющей;

Alfag - угол между осями направляющих;

Xk, Yk, Xk_1, Yk_1, Xk_2, Yk_2 - значения ПФ внешней поступательной кинематической пары.


Пример проектирования гиперболической лемнискаты Баутса

 

Рассмотрим фрагмент программы, в котором определяются передаточные функции кулисно-рычажного механизма Артоболевского для воспроизведения гиперболической лемнискаты Баутса, изображенного на рис. 9.

Рисунок 9 – Кулисно-рычажный механизм Артоболевского

 

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

Формула строения механизма

I(0, 1) → II(2, 3)1 → II(4, 5)3 → II(6, 7)2 .

Рассмотрим план расчета механизма, основываясьна его формуле строения. Вначале, используя процедуру Kriv, следует определить параметры точки А кривошипа. После этого, используя процедуру Assur1, следует определить передаточные функции звеньев 2 и 3, а затем с помощью процедуры Point - передаточные функции точки В. На следующем этапе с помощью процедуры Assur3 следует определить передаточные функции кулисы 4, после чего с помощью процедуры Assur2 определить передаточные функции шатуна 6. В заключение с помощью процедуры Point следует определить параметры точки D и с помощью процедуры Track получить ее траекторию - лемнискату Баутса. В процессе разработки этого этапа программы следует использовать описания процедур, приведенные в конце исходной программы, в соответствии с рекомендациями. С целью упрощения замены формальных параметров фактическими в описании используемых процедур целесообразно вместо переменных, величина которых заведомо равна нулю (например, ПФ1 и ПФ2 неподвижных точек и неподвижных направляющих), подставлять переменную z=0 или просто 0.

В качестве одного из результатов работы этого этапа программы на рис. 10 представлена копия экрана с изображением механизма и траектории точки D.

 

 

Рисунок 10 – Копия экрана с изображением механизма и участком

траектории камня

 


Ниже приводится фрагмент кода программы для создания лемнискаты Баутса.

 

FiOAg0:=0; { Начальный угол поворота кривошипа }

z:=0;

For i:=1 To m+1 Do Begin { Цикл определения ПФ }

ClearDevice; { Очистка графического экрана }

jk:=1;

FiOAg[i]:=FiOAg0+360/(m+1)*(i-1)*jk; { Угол поворота кривошипа }

 

{ КРИВОШИП }

Xo:=0; { Координаты оси вращения кривошипа }

Yo:=0;

LOA:=500; { Длина кривошипа }

q:=1;

Mas:=8; { Масштабный коэффициент }

Kriv (Xo, Yo, LOA, FiOAg[i], q, jk, Mas,

Xa[i], Ya[i], Xa_1[i], Ya_1[i], Xa_2[i], Ya_2[i]);

 

{ ГРУППА 1-ГО ВИДА }

Xc:=1000; { Координаты опоры – точки С }

Yс:=0;

LAB:=1000; { Длины звеньев группы }

LCB:=LOA;

If Ya[i]<=0 Then j1:=1 Else j1:=-1;

{ Это условие приходится вводить, чтобы четырехзвенник ОАВС

всегда был антипараллелограммом }

Assur1 (Xa[i], Ya[i], Xa_l[i], Ya_l[i], Xa_2[i], Ya[i]2,

Xc, Yс, z, z, z, LAB, LCB, j1, q, Mas,

FiABg[i], FiCBg[i], FiAB_1[i], FiCB_l[i], FiAB_2[i], FiCB_2[i]);

 

FiABr[i]:=FiABg[i]/180*pi; { Значения углов в радианах }

FiCBr[i]:=FiCBg[i]/180*pi;

 

{ ОПРЕДЕЛЯЕМ ТОЧКУ В }

r:=5;

Point (Xa[i], Ya[i], Xa_1[i], Ya_l[i], Xа_2[i], Ya_2[i],

FiABg[i], FiAB_1[i], FiAB_2[i],

LAB, z, q, Mas, r,

Xb[i], Yb[i], Xb_1[i], Yb_1[i], Xb_2[i], Yb_2[i]);

 

{ ГРУППА 3-ГО ВИДА }

lcc: = 100; { Длина диагонали камня }

Assur3 (Xb[i], Yb[i], Xb_1[i], Yb_1[i], Xb_2[i],Yb_2[i],

Xо, Yо, z, z, z, z, z, z, Mas, lcc, q,

Fikg[i], Fik_1[i], Fik_2[i]);

Fikr[i]: = Fikg[i]/180*pi; { Значения углов в радианах }

 

{ ГРУППА 2-ГО ВИДА }

LAD: =LOA; { Длины звеньев группы }

j2: =1; { Признак сборки группы }

{ Здесь опущено условие изменения признака сборки группы 2-го вида, требуемое для получения всей лемнискаты }

Assur2 (Xa[i], Ya[i], Xa_1[i], Ya_l[i], Xа_2[i], Ya_2[i],

Xo, Yо, z, z. z, z,

Fikg[i], Fik_1[i], Fik_2[i],

Z, LAD, z, lcc, j2, q, Mas,

FiADg[i], FiAD_1[i], FiAD_2[i]);

FiADr[i]: = FiADg[i]/180*pi; { Значения углов в радианах }

 

{ ОПРЕДЕЛЯЕМ ТОЧКУ D }

Point (Xa[i], Ya[i], Xa_1[i], Ya_l[i], Xа_2[i], Ya_2[i],

FiADg[i], FiAD_1[i], FiAD_2[i], LAD, z, q, Mas, r,

Xd[i], Yd[i], Xd_1[i], Yd_1[i],Xd_2[i], Yd_2[i]);

k:=i;

Track(k, q, Mas, Xd, Yd); { Рисуем траекторию точки D }

End; { Конец цикла }


7 КОНТРОЛЬ ОПРЕДЕЛЕНИЯ ПЕРЕДАТОЧНЫХ ФУНКЦИЙ

 

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

Количественная оценка ПФ0 может быть выполнена в результате построения кинематической схемы в масштабе методом засечек для одного или нескольких положений начального звена. Рассмотрим порядок построения кинематической схемы на примере механизма, изображенного на рис. 9.

Известны координаты точки С и длины звеньев lOA, lAB, lCB, lAD. Сначала следует выбрать стандартную величину масштаба (например, 2:1, 1:1, 1:2 и т.д.) и вычислить масштабные значения координат точки С и длин звеньев. Выбрав на чертеже начало координат (точку О), следует по координатам xC и yC построить точку С. После этого для некоторого угла поворота кривошипа определить положение точки А. Угол, определяющий положение кривошипа, следует отсчитывать от положительного направления оси х против часовой стрелки. При этом угол кривошипа должен быть определен по формуле

, (1)

где в качестве i принимается одно из значений от 1 до 48.

В качестве начального угла кривошипа целесообразно принять угол, соответствующий началу рабочего хода, о чем более подробно смотри в начале раздела 8. После этого раствором циркуля, равным масштабным значениям длин звеньев lAB и lCB, сделать засечки из центров А и С и на их пересечении получить точку В. Углы, определяющие положение звеньев АВ и СВ, нужно измерить транспортиром, отсчитывая от прямой, проведенной параллельно оси х через начальную точку в обозначении звена, против часовой стрелки. После этого следует провести звено 4 и определить его угол наклона jOB. Далее необходимо из точки А сделать засечку радиусом lAD на оси звена 4, получить точку D и измерить угол jAD. После этого необходимо определить координаты всех характерных точек механизма (осей вращательных кинематических пар и центров масс звеньев). Результаты графических построении следует сравнить с результатами, полученными на базе аналитических методов в процессе работы программы. Для этого необходимо использовать процедуру Grafik.

Описание процедуры:

Grafik (m: Integer; OK, Y1, Y2, Y3: Massiv; S1, S2, S3: String),

где OK – формальное имя обобщенной координаты, например, угла поворота кривошипа, являющегося абсциссой представляемых графиков;

m - число значений элементов массивов, представляемых на экране;

Y1, Y2, Y3 - формальные имена массивов, графики которых будут представлены на экране;

S1, S2, S3 - формальные строковые переменные, служащие для идентификации графиков на экране.

Поскольку все три графика строятся в одинаковом масштабе, целесообразно на экран выводить массивы, максимальные величины которых имеют одинаковый порядок, например передаточные функции точек механизма. В связи с этим при выводе на экран передаточных функций, характеризующих вращательное движение звеньев, углы следует выводить в радианах. Первый график представляется на экране зеленой линией, второй - красной и третий ‑ желтой. Кроме того, на экран выводится вертикальная линия, которая с помощью клавиш “” и “®” может перемещаться влево и вправо и занимать положение, соответствующее любому из m положений кривошипа. При этом на экран выводятся численные значения обобщенной координаты (например, угла, определяющего положение кривошипа), а также значения выводимых функций. На рис. 11 представлена копия экрана с изображением графиков для абсциссы точки D механизма. Здесь сплошная линия соответствует зеленой линии, штриховая – красной, а пунктирная – желтой. Установив вертикальную линию в положение, соответствующее построенной кинематической схеме, можно сравнить расчетное значение координаты со значением, полученным графическим методом.

Рисунок 11 – Графики передаточных функций камня

 

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

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

xA, yA, jAB, jCB, , , , , xB, yB,

jOB, , , jAD, , , xD, yD.

При этом нужно добиваться правильности анализируемого параметра и только после этого переходить к анализу следующего. В случае правильности ПФ0 рассматриваемого параметра с помощью выведенных на экран графиков можно дать качественную оценку правильности ПФ1 и ПФ2. Рассмотрим это на примере графиков, представленных на рис. 11.

Эти графики связаны между собой операцией дифференцирования. При этом на интервале, где функция xD (сплошная линия) убывает, ее производная (штриховая линия) отрицательна, а на интервале возрастания x D производная положительна. Точке экстремума графика хD соответствует значение = 0. Аналогичную связь можно обнаружить и между графиками и (пунктирная линия). Следует также обратить внимание на тот факт, что абсцисса, где = 0, соответствует точке перегиба на графике хD.

Для количественной оценки правильности передаточных функций следует использовать процедуру TestPF. В ее основу положен метод численного дифференцирования. Значения первой и второй производных функции в какой-либо точке можно определить, зная ее значение в этой и трех соседних точках справа и слева, взятых с одинаковым шагом Dj, по формулам:

(2)

. (3)

Метод численного дифференцирования является менее точным, чем аналитический, однако расчеты показывают, что для достаточно плавных функций производные, найденные этими методами, отличаются весьма незначительно. На этом и основано применение рассматриваемой процедуры которая предусматривает ввод ПФ0, ПФ1 и ПФ2 какого-либо элемента механизма, расчет ПФ1 и ПФ2 по формулам (2) и (3) и сравнение передаточных функций, определенных двумя способами. Если разница этих величин не превышает 0,1%, аналитические значения передаточных функций можно считать правильными. В противном случае на экран выводятся значения передаточных функций длятех точек, где их разница превышает 0,1%.

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

В некоторых случаях может наблюдаться разница значений передаточных функций, найденных двумя методами, в семи последовательных точках. Это возможно, когда функция претерпевает разрыв, например, для графика угла звена механизма, где величина разрыва первого рода может равняться периоду тригонометрической функции, используемой для определения угла наклона звена. Этот разрыв можно наблюдать на графике ПФ0, после чего следует принять меры по его устранению. Значения же ПФ1 и ПФ2 следует считать правильными.

Описание процедуры:

TestPF (m: Integer; Yf, Yf1, Yf2: Massiv; StrYf: String),

где m - число значений элементов массива;

Yf, Yf1, Yf2 - формальные имена передаточных функций, определенных аналитически;

StrYf - формальная строковая переменная, идентифицирующая тестируемый параметр. Шаг численного дифференцирования .

При тестировании углов, определяющих положение звеньев, значение ПФ0 следует вводить в радианах.

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


8 ОПРЕДЕЛЕНИЕ ПАРАМЕТРОВ ДИНАМИЧЕСКОЙ

МОДЕЛИ МЕХАНИЗМА

 

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

При определении угла кривошипа, соответствующего началу рабочего хода, следует иметь в виду, что выходное звено занимает при этом крайнее положение и его ПФ1 равна нулю. Методика определения такого угла кривошипа следующая. В первом приближении его следует принять по графику ПФ1 равным значению левой границы интервала, на котором ПФ1 меняет знак. Присвоив это значение переменной FiOA0g и уменьшив шаг dFi, следует найти новое значение начального угла кривошипа. Такие уточнения следует продолжить до тех пор, пока величина ПФ1 на левой границе интервала изменения знака не окажется по модулю меньшей 0,0005. Это значение обобщенной координаты и будет соответствовать началу рабочего хода. Присвоив после этого шагу dFi первоначальное значение, следует продолжить разработку программы.

При последующей работе над программой с целью сокращения времениее отладки параметру q можно присвоить значение 0 (на экран не будет выводиться изображение механизма) и с помощью скобок { } отключить процедуры Grafik и TestPF.

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

, (4)

где , - проекции внешней силы, приложенной в точке Ki -го звена;

, - ПФ1 точки Ki;

Мi - внешний момент, приложенный к i -му звену;

- угловая ПФ1 i -гo звена;

- проекция силы тяжести i -гo звена;

- ПФ1 центра масс i -гo звена.

Суммирование проводится по всем подвижным звеньям механизма, кроме кривошипа.

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

В качестве примера рассмотрим задание силы технологического сопротивления для механизма, у которого выходное звено (ползун), обозначаемое буквой D, совершает возвратно-поступательное движение параллельно оси абсцисс. Рабочий ход совершается при движении ползуна справа налево. График зависимости силы сопротивления от положения ползуна на рабочем ходу показан на рис. 12. При холостом ходе сила сопротивления равна нулю.

Рисунок 12 – График силы сопротивления

На этом графике:

- рабочий ход ползуна, ;

- текущее значение перемещения ползуна, ;

- максимальное значение силы технологического сопротивления, задаваемое в исходных данных на проектирование.

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

при £ 0,2: =0;

при 0,2 < £ 0,4: = ( -0,2);

при 0,4< £ 0,7: = ;

пpи 0,7< £ 1,0: = (1- )/0,3;

при холостом ходе: = 0.

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

For i:=1 to 49 Do Begin

if XD[i]>XDmax Then XDmax:=XD[i];

if XD[i]<XDmin Then XDmin:=XD[i];

End; {For}

SDmax:=XDmax-XDmin;

Pcmax:=10000;

For i:=1 to 49 Do Begin

S_[i]]:=(XDmax-XD[i])/SDmax;

if S_[i]<=0.2 Then Pcx[i]:=0;

if (S_[i]>0.2) And (S_[i]<=0.4) Then Pсх[i]:=Pcmax/0.2*(s_[i]-0.2);

if (S_[i]>0.4) And (S_[i]<=0.7) Then Pcx[i]:=Pcmax;

if (S_[i]>0.7) Then Pcх[i]:=Pcmax/0.3*(1-s_[i]);

if Xd_1[i]>0 Then Pcx[i]:=0;

End; {For}

Для проверки правильности задания силы сопротивления на экран следует вывести график = f(OK). При этом наклонные прямолинейные участки на графике (см. рис. 12) вероятнее всего станут криволинейными.

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

; (5)

, (6)

где mi, - масса и момент инерции i -гo звена;

, , , - ПФ1 и ПФ2 центра масс i -гo звена;

, - угловые ПФ1 и ПФ2 i -гo звена.

В исходной программе эти формулы приведены, и в них лишь следует внести коррективы, учитывающие особенности рассчитываемого механизма. Полученные функции и следует с помощью процедуры Grafik вывести на экран и проанализировать, а затем протестировать с помощью процедуры TestPF. При использовании процедуры Grafik в качестве 3-й функции удобно использовать пустой массив zz [i]=0, который следует объявить в программе. При использовании процедуры TestPF вместо массива ПФ2 также можно ввести пустой массив и не тестировать вторую производную.



Поделиться:


Последнее изменение этой страницы: 2017-02-19; просмотров: 243; Нарушение авторского права страницы; Мы поможем в написании вашей работы!

infopedia.su Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Обратная связь - 52.14.121.242 (0.082 с.)