lect_05 (966005), страница 2
Текст из файла (страница 2)
При заданных значениях функций для n+1 равноотстоящих значений аргумента
квадратурные формулы Ньютона-Котеса имеют вид:
правило трапеций для n шагов
правило трапеций для n=1
правило Симпсона для n=2
правило Уэддля для n=6
При вычислениях на ЭВМ используют программы, имеющиеся в каталоге конкретной машины (например, QTFG или QSF).
При графическом определении интеграла подынтегральная функция задается графиком. Для примера рассмотрим определение угла поворота выходного звена по заданной кривой
, полученной экспериментально.
График угловой скорости изображается в декартовых координатах с учетом числовых значений масштабов: угловой скорости
и времени
. Промежуток времени от
до
, делится на такое количество интервалов
, которое позволяет считать, что на каждом малом промежутке времени
движение можно принять равномерным.
Эти промежутки времени, отмеченные на рис. 5.7 , а точками 0, 1, 2, 3, 4, не обязательно должны быть равными.
В каждом интервале времени, например от до
можно приближенно считать, что
т. е. можно принять, что площадь криволинейной трапеции равновелика площади прямоугольника высотой и основанием
.
Концы средних ординат для каждого интервала проецируют на ось ординат и соединяют найденные точки 1', 2', 3' , ... , i' с точкой D, которая ограничивает слева выбранный отрезок интегрирования OD длиной К, мм (рис. 5.7, a).
Лучи D1', D2', D3', ... , проведенные через точку D, образуют углы D1’, с положительным направлением оси х, причем
.
На искомом графике (рис. 5.7, б) проводят линии 01", 1"2" , 2"3" , ... , параллельные в пределах соответствующих интервалов лучам Dl', D2' , D3', ... . Первый отрезок 01" проводят через начало координат 0, следующие отрезки соответственно через точку 1", затем через точку 2" и т. д. Эти линии наклонены относительно положительного направления оси х под углами
соответственно, т. е.
Отрезки на графиках связаны с соответствующими физическими параметрами с помощью масштабов соотношениями:
Рис. 5.7
Приравнивая правые части написанных выше соотношений для тангенса угла
Откуда масштаб искомого графика:
Графическое и численное дифференцирование.
Графическое дифференцирование начинают с построения графика функции по заданным значениям. При экспериментальном исследовании такой график получают с помощью самопишущих приборов. Далее проводят касательные к кривой в фиксированных положениях и вычисляют значения производной по тангенсу угла, образованного касательной с осью абсцисс.
На рис. 5.8, а изображена кривая
полученная экспериментально на установке (рис. 5.6). Определение углового ускорения (искомой функции) проводят графическим дифференцированием по соотношению:
Тангенс угла наклона касательной к кривой
в некоторой точке i представляют в виде отношения отрезков
, где К – выбранный отрезок интегрирования (рис. 5.8, б)
После подстановки этого соотношения в соотношение (5.19) получают
где - ордината искового графика углового ускорения;
— масштаб искомого графика ; единицы СИ:[
] = мм; [
] = мм/(рад
с-2).
График функции строят по найденным значениям ординат для ряда позиций. Точки на кривой соединяют от руки плавной линией, а затем обводят с помощью лекала.
Графическое дифференцирование рассмотренным методом касательных имеет относительно низкую точность. Более высокую точность получают при графическом дифференцировании методом хорд (рис. 5.8, в и г).
Н
а заданной кривой отмечают ряд точек 1", 2", 3", которые соединяют хордами, т.е. заменяют заданную кривую ломаной линией. Принимают следующее, допущение: угол наклона касательных в точках, расположенных посередине каждого участка кривой, равен углу наклона соответствующей хорды. Это допущение вносит некоторую погрешность, но она относится только к данной точке. Эти погрешности не суммируются, что обеспечивает приемлемую точность метода.
Рис (5.8)
Остальные построения аналогичны ранее описанным при графическом дифференцировании методом касательных. Выбирают отрезок (мм); проводят лучи, наклоненные под углами
до пересечения с осью ординат в точках 1', 2', 3' ... , которые переносят на ординаты, проведенные в середине каждого из интервалов. Полученные точки 1*, 2*, 3* являются точками искомой функции
.
Масштабы по осям координат при этом методе построения связаны таким же соотношением (5.21), которое было выведено для случая графического дифференцирования методом касательных.
Дифференцирование функции f(x), заданной (либо вычисленной) в виде массива чисел, выполняют методом численного дифференцирования с применением ЭВМ.
Чем меньше шаг в массиве чисел, тем точнее можно вычислить значение производной функции в этом интервале
Можно пользоваться также выражением
При численном дифференцировании используют интерполяционные формулы, которые сопоставляют заданные значения какой-либо величины с функцией известного .класса, зависящей от нескольких параметров, выбранную так, чтобы при заданных значениях аргумента (в узлах интерполяции) значения функции совпадали с заданными значениями величины, т. е. чтобы график функции проходил через заданные точки. Численное дифференцирование чувствительно к ошибкам, вызванным неточностью исходных данных. Для функции у(х), заданной таблицей разностей для равноотстоящих значений аргумента с шагом , используют следующие соотношения для вычисления аргумента и производных:
При разработке прикладных программ для численного дифференцирования на ЭВМ используют интерполяционные формулы Стирлинга, Бесселя, Ньютона и др.
Метод преобразования координат.
Применение ЭВМ для кинематического анализа механизмов связано с разработкой соответствующих алгоритмов и программ расчёта. Наиболее просто такие алгоритмы реализуются с использованием уравнений преобразования координат в матричной форме записи необходимых операций вычисления.
При этом методе выбирают некоторое число систем координат, достаточное для математического описания геометрической формы звеньев и относительного движения звеньев в каждой кинематической паре. Число систем координат определяется числом элементов звеньев, образующих кинематические пары. Неподвижная система координат связана со стойкой. В каждой кинематической паре выбирают две системы координат (способ 1) или одну систему координат (способ 2). При 1-м способе две системы координат относятся к элементам пары звеньев, образующих эту пару. При втором способе каждой кинематической паре соответствует прямоугольная система координат, одна из осей которой связана с характерными признаками звена, например осевой линией. Для примера на рис. 5.9, а показаны координатные оси
(или
) четырехзвенной открытой кинематической цепи из звеньев 1, 2, 3, 4, моделирующей структуру руки человека (рис. 5.9, б). Ось
направляют вдоль оси пары, а ось
дополняет правую систему координат
.
Начало координат каждой i-й локальной координатной системы совмещают с той кинематической парой, которой данное звено соединено с предыдущим звеном. Для плоских механизмов оси параллельны между собой, так как они перпендикулярны базовой плоскости, в которой рассматривается движение звеньев плоского механизма.
П
Рис 5.9
ереход от i-й локальной координатной системы к другой (i+1) системе определяется уравнениями преобразования декартовых прямоугольных координат, в общем случае – переноса и поворота координатных осей.Применение метода будет продемонстрированною на конкретном примере в лекции 25 «Манипуляторные роботы».
Контрольные вопросы к лекции N5:
-
Как построить графически функцию положения механизма и её производные?
-
Как рассчитать масштабы кинематических диаграмм?
-
Как определить величину и направление угловых скоростей и угловых ускорений звеньев?
-
Как аналитически определить функцию положения, передаточные функции скорости и ускорения ползуна кривошипного механизма?
-
В чём заключается преимущества и недостатки аналитического и графического методов кинематического анализа?