книги / Статика в задачах биомеханики
..pdfзов руки является тяговый протез. Тяговый (механический) протез управляется с помощью тяг и полностью контролируется усилиями самого человека. Наиболее распространенные варианты тягового управления – от движения плеча, локтя, лучезапястного сустава, проксимальной фаланги пальца. Сильная сторона такого механизма – возможность контролировать усилие. При выполнении хвата пользователь сам определяет силу сжатия, скорость и может почувствовать сопротивление, когда кисть / крюк упирается в предмет и не может продолжать сдавливание. Механизм прост в управлении и обслуживании. Слабая сторона − это ограниченность силы хвата возможностями самого человека. При длительной практике ношения пользователи демонстрируют возможность манипуляции различными предметами – могут писать, пользоваться ножом и вилкой, играть втеннис, зажигатьспички, переноситьгрузит.п.
Несмотря на значительный прогресс в области протезирования, совершенный учеными и медиками за последние пару веков, остается достаточно много проблем, связанных с высокой технологичностью и сложностью подобного оборудования. Для создания протеза необходимы обширная материально-технологическая база и целая команда специалистов из различных областей: ученых, врачей, инженеров, конструкторов, программистов.
Для развития навыков и умений ставить и решать задачи проектирования протезов как руки, так и нижней конечности человека рассмотрим близкую и имеющую самостоятельное значение задачу о расчете устойчивости мобильного антропоморфного робота-манипулятора.
Расчет на устойчивость имеет важное значение при проектировании мобильных грузоподъемных механизмов. Под устойчивым понимается такое положение манипулятора, при котором любое допустимое нагружение не может привести к потере его равновесия или к опрокидыванию.
Большое влияние на устойчивость оказывает распределение массы между различными частями конструкции. Для повышения общей устойчивости робота-манипулятора часто увеличивают
41
вес его центральной части, что в большинстве случаев приводит к необоснованному перерасходу металла. Таким образом, возникает противоречивая ситуация. С одной стороны, для повышения устойчивости необходимо увеличивать вес некоторых частей машины, а с другой стороны, необходимо добиваться уменьшения обшей металлоемкости конструкции. Это возможно только путем постановки и решения соответствующей оптимальной задачи, для чего необходимо в достаточной мере владеть аппаратом теоретической механики и знать методы оптимизации.
Рассмотрим методику решения подобных задач на примере минимизации веса промышленного робота-манипулятора, широко применяемого на различных промышленных предприятиях при выполнении таких операций, как сварка, пайка, сверление, манипуляции с грузами. Это, на наш взгляд, позволит глубже разобраться в различных подходах к проблеме оптимизации конструкции антропоморфных роботов-манипуляторов, приобретающей все более важное значение для современного производства.
Условие задачи
Промышленный робот-манипулятор (рис. 1.1) размещен на горизонтально-ориентированной мобильной платформе. Расстояние между опорами мобильной платформы равно DE.
Центр тяжести мобильной платформы и основания манипулятора (общим весом РB) находится в точке B и лежит в плоскости ее симметрии (т.е. DB = BE).
Робот-манипулятор удерживает груз весом PG. Центр тяжести груза расположен в точке G. Наибольший груз, который может поднять манипулятор, имеет вес РGmax. Угол подъема груза (α) может изменяться впределах0°–80°.
Центр тяжести звеньев манипулятора (весом РC) находится в точке C. Конструкция манипулятора такова, что при подъеме груза угол наклона звеньев манипулятора β изменяется пропорционально углу α:
β = α + 20°. |
(1.136) |
42
Рис. 1.1. Промышленный робот-манипулятор
Центр тяжести противовеса (весом РA) находится в точке A. Расстояние между опорой D и проекцией центра тяжести противовеса на горизонтальную плоскость обозначим LA.
Все необходимые параметры робота-манипулятора для всех вариантов заданий приведены в табл. 1.1.
Таблица 1.1
Параметры робота-манипулятора
№ |
BC, |
BG, |
DE, |
РC, |
РGmax, |
Способ выполнения задания |
δ, |
|||
п/п |
м |
м |
м |
кН |
кН |
I |
|
II |
кН |
|
LA, м |
Lmin, м |
Lmax, м |
||||||||
1 |
1,3 |
3,4 |
0,8 |
21,7 |
9,0 |
0,8 |
0,5 |
|
1,1 |
14,3 |
2 |
1,4 |
3,6 |
0,8 |
21,7 |
9,0 |
0,8 |
0,5 |
|
1,1 |
14,3 |
3 |
1,5 |
3,8 |
0,9 |
21,7 |
9,0 |
0,8 |
0,5 |
|
1,1 |
14,3 |
4 |
1,6 |
4,0 |
0,9 |
21,7 |
9,0 |
0,8 |
0,5 |
|
1,1 |
14,3 |
5 |
1,7 |
4,2 |
1,0 |
21,7 |
9,0 |
0,8 |
0,5 |
|
1,1 |
14,3 |
6 |
1,8 |
4,4 |
1,0 |
43,4 |
18,0 |
1,0 |
0,3 |
|
1,7 |
14,3 |
7 |
1,9 |
4,6 |
1,1 |
43,4 |
18,0 |
1,0 |
0,3 |
|
1,7 |
14,3 |
8 |
2,0 |
4,8 |
1,1 |
43,4 |
18,0 |
1,0 |
0,3 |
|
1,7 |
14,3 |
9 |
2,1 |
5,0 |
1,2 |
43,4 |
18,0 |
1,0 |
0,3 |
|
1,7 |
14,3 |
|
|
|
|
|
|
|
|
|
|
43 |
Окончание табл. 1.1
№ |
BC, |
BG, |
DE, |
РC, |
РGmax, |
Способ выполнения задания |
δ, |
|||
п/п |
м |
м |
м |
кН |
кН |
I |
|
II |
кН |
|
LA, м |
Lmin, м |
Lmax, м |
||||||||
10 |
2,2 |
5,2 |
1,2 |
43,4 |
18,0 |
1,0 |
0,3 |
|
1,7 |
14,3 |
11 |
0,9 |
2,2 |
0,6 |
5,9 |
6,0 |
0,4 |
0,1 |
|
0,7 |
4,6 |
12 |
0,9 |
2,3 |
0,6 |
5,9 |
6,0 |
0,4 |
0,1 |
|
0,7 |
4,6 |
13 |
1,0 |
2,4 |
0,7 |
5,9 |
6,0 |
0,4 |
0,1 |
|
0,7 |
4,6 |
14 |
1,0 |
2,5 |
0,7 |
5,9 |
6,0 |
0,4 |
0,1 |
|
0,7 |
4,6 |
15 |
1,1 |
2,6 |
0,8 |
5,9 |
6,0 |
0,4 |
0,1 |
|
0,7 |
4,6 |
16 |
1,1 |
2,7 |
0,8 |
3,7 |
12,0 |
0,5 |
0,2 |
|
0,8 |
4,6 |
17 |
1,2 |
2,8 |
0,9 |
3,7 |
12,0 |
0,5 |
0,2 |
|
0,8 |
4,6 |
18 |
1,2 |
2,9 |
0,9 |
3,7 |
12,0 |
0,5 |
0,2 |
|
0,8 |
4,6 |
19 |
1,3 |
3,0 |
1,0 |
3,7 |
12,0 |
0,5 |
0,2 |
|
0,8 |
4,6 |
20 |
1,3 |
3,1 |
1,0 |
3,7 |
12,0 |
0,5 |
0,2 |
|
0,8 |
4,6 |
21 |
0,7 |
1,8 |
0,3 |
0,8 |
0,1 |
0,2 |
0,1 |
|
0,3 |
0,4 |
22 |
0,6 |
1,8 |
0,3 |
0,8 |
0,1 |
0,2 |
0,1 |
|
0,3 |
0,4 |
23 |
0,7 |
1,9 |
0,3 |
0,8 |
0,1 |
0,2 |
0,1 |
|
0,3 |
0,4 |
24 |
0,6 |
1,9 |
0,4 |
0,8 |
0,1 |
0,2 |
0,1 |
|
0,3 |
0,4 |
25 |
0,7 |
2,0 |
0,4 |
0,8 |
0,1 |
0,2 |
0,1 |
|
0,3 |
0,4 |
26 |
0,6 |
2,0 |
0,4 |
0,4 |
0,2 |
0,3 |
0,2 |
|
0,4 |
0,4 |
27 |
0,7 |
2,1 |
0,5 |
0,4 |
0,2 |
0,3 |
0,2 |
|
0,4 |
0,4 |
28 |
0,6 |
2,1 |
0,5 |
0,4 |
0,2 |
0,3 |
0,2 |
|
0,4 |
0,4 |
29 |
0,7 |
2,2 |
0,5 |
0,4 |
0,2 |
0,3 |
0,2 |
|
0,4 |
0,4 |
30 |
0,6 |
2,2 |
0,6 |
0,4 |
0,2 |
0,3 |
0,2 |
|
0,4 |
0,4 |
Требуется определить вес основания манипулятора (РB), положение (расстояние LA) и вес (РА) противовеса, при которых общий вес робота-манипулятора будет минимальным. Считается, что размеры и вес звеньев манипулятора неизменны, а действием динамических сил можно пренебречь. При этом должны выполняться следующие условия устойчивости манипулятора:
ND ≥ δ,
NE ≥ δ, (1.137)РG [0, РGmax],
α [0°, +80°],
где ND, NE – реакции опоры в точках D и E, направленные по нормали к опорной поверхности; δ − некоторая положительная величина,
44
характеризующая заданный запас на устойчивость мобильного ро- бота-манипулятора; РG – вес груза; α – уголподъема груза.
Условия (1.137) должны выполнятся одновременно. Очевидно, что при максимальном весе поднимаемого груза
РGmax и угле α, равном 0°, вес противовеса РA также должен быть наибольшим. Однако вес противовеса не должен быть слишком велик, иначе робот-манипулятор может опрокинуться в тот момент, когда он находится без груза, а звенья манипулятора подняты в верхнее положение. Понятно, что в данном случае необходимы постановка и решение задачи оптимизации.
Постановка задачи оптимизации
Исходя из условий задачи, в качестве критерия оптимизации выберем общий вес манипулятора (Робщ):
Робщ = РA + РB + РC. |
(1.138) |
Теперь данную задачу можно сформулировать как задачу минимизации с ограничениями: необходимо найти вес основания (РB), вес (РA) противовеса и положение (LA) противовеса, доставляющие минимум следующему критерию оптимизации:
min Pобщ{LA, РA, РB}. |
(1.139) |
LA ,PB ,PA
При ограничениях (1.137), обеспечивающих устойчивость конструкции, и ограничениях, накладываемых на искомые величины:
РB ≥ 0 кH, |
|
РA ≥ 0 кH, |
(1.140) |
Lmin ≤ LA ≤ Lmax.
Искомые величины LA, РB, РA называют параметрами оптимизации. Как следует из теории оптимизации [25, 41], задача относится к задаче линейного программирования в том случае, если
45
минимизируемая функция и ограничения являются линейными относительно параметров оптимизации. В противном случае задача является задачей нелинейного программирования.
Для определения условий выполнения ограничений (1.137) необходимо составить уравнения моментов всех сил, действующих на робот-манипулятор, относительно точек D и E. Нетрудно заметить, что, поскольку в эти уравнения войдут произведения параметров оптимизации, это делает ограничения (1.137) нелинейными. Поставленная задача оптимизации (1.137)–(1.140) в этом случае является задачей нелинейного программирования. Однако если расстояние LA принимает заданное значение, то данная оптимальная задача превращается в задачу линейного программирования. Это может быть использовано при ее решении.
Метод решения задачи оптимизации
В данном учебном пособии предусмотрено несколько способов решения поставленной оптимальной задачи.
Первый способ выполнения расчетной работы предполагает, что положение противовеса известно и расстояние LA берется из таблицы. В этом случае оптимальная задача (1.137)–(1.140) является задачей линейного программирования и ее решение получается аналитически, без привлечения ЭВМ. Как известно [25], основным методом решения подобных задач является симплекс-метод. Идея этого метода состоит в построении допустимого многогранника и упорядоченном исследовании его вершин. Под допустимым многогранником в этом случае понимается область допустимых решений оптимальной задачи, образуемая пересечением m множеств. Каждое из них определяется соответствующимнеравенством:
ai1x1 + ... + ainxn ≤ bi, |
(1.141) |
и представляет собой полупространство, лежащее по одну сторону от гиперплоскости:
ai1x1 + ... + ainxn = bi. |
(1.142) |
46
В выражениях (1.141), (1.142) под xj, j = 1,n, понимается со-
вокупность параметров оптимизации; aij, bi (i = 1,m , j = 1,n ) есть
заданные константы.
Из теории линейного программирования известно [25], что решение оптимальной задачи находится в одной из вершин такого многогранника.
Применение симплекс-метода рассмотрим на примере данной задачи (1.137)–(1.140) при фиксированном значении величины LA. В этом случае останется только два параметра оптимизации РA и РB. В соответствии с уравнениями (1.137), ограничения
задачи можно записать в виде |
|
a11PA + a12PB ≥ b1, |
|
a21PA + a22PB ≥ b2, |
(1.143) |
PA ≥ 0, |
|
PB ≥ 0. |
|
Рис. 1.2. Допустимый многогранник в пространстве переменных РA и РB
В соответствии с системой неравенств (1.143), нетрудно построить допустимый многогранник в пространстве переменных РA и РB. Пример построения такого многогранника показан на рис. 1.2. Здесь множество допустимых решений задачи неограни-
47
ченно. Из рисунка видно, что решением задачи может быть только точка Q1 или точка Q2. Окончательный выбор решения зависит от угла наклона прямой a21PA + a22PB = b2 к оси абсцисс. В данном случае оптимальное решение задачи находится в точке Q1.
При выполнении расчетной работы вторым способом расстояние LA считается искомой величиной и решение задачи усложняется. В этом случае аналитическое решение затруднено и предполагается использование ЭВМ.
Решение задачи вторым способом заключается в том, что решение задачи нелинейного программирования заменяется решением последовательности задач линейного программирования при различных фиксированных значениях параметра LA. В этом случае отрезок [Lmin, Lmax] разбивается на N участков. Далее задача линейного программирования решается N раз с помощью сим- плекс-метода. На каждом шаге определяются оптимальные РAi и РBi и анализируется величина критерия оптимальности. То решение (РAK, РBK, LAK), при котором общий вес манипулятора будет наименьшим, считается оптимальным. Точность решения задачи в данном случае зависит от числа разбиений отрезка [Lmin, Lmax], поэтому здесь выгодно использование ЭВМ.
Для этого необходимо разработать алгоритм решения, основанный на многократном использовании симплекс-метода. Как было показано выше, решение задачи линейного программирования при фиксированном значении параметра LA находится в точке пересечения прямых:
a11PA + a12PB = b1, |
(1.144) |
a21PA + a22PB = b2, |
(1.145) |
где коэффициенты aij (i, j = 1,2) зависят от расстояния LA. Для решения задачи необходимо аналитически разрешить систему двух алгебраических уравнений относительно весов РA и РB как функций от параметра LA. Иными словами, решение системы уравнений (1.144)–(1.145) необходимо записать в виде
48
РA = (b1 − a12(a11b2 − a21b1)/(a11a22 − a21a12)) / a11, |
(1.146) |
РB = (a11b2 − a21b1)/(a11a22 − a21a12). |
(1.147) |
Отметим, что функции РB(LA) и РA(LA) лучше представить в явном виде, который зависит от исходных данных задачи.
Далее, организовав перебор параметра LA [Lmin, Lmax], можно легко найти оптимальное решение задачи. Конечно, последнюю операцию, требующую больших вычислений, удобно поручить ЭВМ. С этой целью предлагается составить небольшую программу, реализующую приведенный выше алгоритм, и получить решение задачи с помощью ЭВМ.
Порядок решения задачи
1.Ознакомиться с методами решения задач статики твердого тела и методами оптимизации по приведенному списку литературы.
2.Начертить расчетную схему промышленного робота-ма- нипулятора с указанием всех размеров, сил и реакции связей, действующих на элементы конструкции. Исходные данные задачи взять из приведенной таблицы (см. табл. 1.1) согласно выданному варианту и способу выполнения расчетной работы.
3.Из условий равновесия [13, 29, 33, 43] найти зависимо-
сти значений давления робота-манипулятора на опоры ND и NE от параметров PA, PB, PC, РGmax, LA. Получить условия устойчи-
вости робота-манипулятора (1.137), не забывая, что при его работе может изменяться как угол подъема α, так и вес поднимаемого груза РG.
4. Исследовать устойчивость робота-манипулятора в зависимости от веса поднимаемого груза РG [0, РGmax]. При выполнении расчетной работы первым способом величина LA берется из табл. 1.1. Затем необходимо составить систему ограничений (1.139), построить допустимый многогранник и найти оптимальное решение.
49
5.При выполнении расчетной работы вторым способом для исследования устойчивости робота-манипулятора необходимо:
а) составить систему уравнений (1.144)–(1.145), выразив ко-
эффициенты aij и bi, i, j = 1, 2, через переменную величину LA и другие параметры задачи: DE, РGmax, PC, δ;
б) cоставить для ЭВМ программу поиска оптимальных зна-
чений PA, PB и LA. При этом числе разбиений заданного отрезка [Lmin, Lmax] выбирается произвольно, но не менее 20.
в) провести расчеты на ЭВМ и получить оптимальные значения параметров PA, PB и LA.
6.Для всех способов выполнения расчетной работы проверить выполнениеусловийустойчивостиробота-манипулятора (1.137) при найденных значениях параметровоптимизации.
7.Составить отчет о работе, обязательно включающий постановку задачи оптимизации, расчетную схему, описание метода решения оптимальной задачи, анализ полученных результатов и необходимые выводы.
50