Аналитика элементарных модулей
Автор: Семенов Б.П.
Журнал: Известия Самарского научного центра Российской академии наук @izvestiya-ssc
Рубрика: Механика и машиностроение
Статья в выпуске: 2 т.4, 2002 года.
Бесплатный доступ
Аналитические решения задачи о функциях и производных первого и второго порядков функций векторных модулей могут быть использованы в системах управления, в частности, робототехнических систем.
Короткий адрес: https://sciup.org/148197713
IDR: 148197713 | УДК: 621.01
Analytics of elementary modules
Analytical decisions of a task on functions and derivatives of the first and second orders of functions of vector modules represent practical interest, in particular, at the analysis and synthesis various on structure and purpose of mechanisms.
Текст научной статьи Аналитика элементарных модулей
Самарский государственный аэрокосмический университет
Аналитические решения задачи о функциях и производных первого и второго порядков функций векторных модулей могут быть использованы в системах управления, в частности, робототехнических систем.
Модуль векторной модели - векторный замкнутый контур или минимальная по составу система взаимо связанных векторных замкнутых контуров с определенностью между функциями и аргументами.
Под определённостью между функцией и аргументом будем понимать возможность вычисления значений функций контура по заданным значениям его аргументов и констант.
Все многообразие векторных модулей можно разделить на одноконтурные (простые) и многоконтурные (с обратными связями).
Элементарный векторный модуль - одноконтурный модуль с числом функций, равным числу измерений контура.
Элементарные плоские и пространственные векторные модули, приведены в табл.1. Практический интерес для определения частных производных представляют приведённые в этой таблице основные аргументы одно- и двухаргументных модулей.
Функции элементарных плоских модулей
Аналитические решения для функций элементарных m-векторных плоских модулей табл.1 получим из условия замкнутости:
∑m pi ⋅cos ϕi = 0, (1) i=1
m
∑ pi⋅sin ϕi = 0. (2)
i=1
Используем вариант индексации векторов: i = 1, j = 2 и k = 3.
Модуль Пл1 (функции p1 и ϕ 1 ). Преобразуем (1) и (2)
p1 ⋅ cos ϕ1 = -X,
p1 ⋅ sin ϕ1 = -Y, где
X= ∑ m pi ⋅ cos ϕ i,Y= ∑ m pi ⋅ sin ϕ i, i=2 i=2
тогда p1 = VX2 +Y2 (3)
и
ϕ 1 = sign(- ⋅ Ya)rccos(-pX1/) . (4)
Модуль Пл2 (функции p1 и p2 ). Преобразуем (1) и (2)
p1 ⋅ cosϕ(1 - ϕ2 )+ p2 + X = 0 , (5)
p1 ⋅ sin( ϕ 1 - ϕ 2)+Y=0, (6)
где
m
X=∑ pi ⋅cos( ϕi -ϕ2), i=3
Рис.1. Геометрия модуля Пл2
Таблица 1. Элементарные векторные модули
|
s m |
Св И И |
hQ О s |
Обобщённый список функций |
a Ch О П У S T u к s |
Основные аргументы |
||||||||||||||||||
|
двухаргументных модулей |
|||||||||||||||||||||||
|
одноаргументных модулей |
|||||||||||||||||||||||
|
u 1 |
u 2 |
u 3 |
ri |
α i |
β i |
rj |
α j |
β j |
rk |
α k |
β k |
r i1 |
α i1 |
β i1 |
r i2 |
α i2 |
β i2 |
||||||
|
Пространственные |
|||||||||||||||||||||||
|
I |
I |
Пр 1 |
ri |
α i |
β i |
1 |
x |
x |
x |
x |
x |
x |
|||||||||||
|
II |
2 |
Пр 2 |
ri |
α i |
rj |
1 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
||||||||
|
Пр 3 |
ri |
α i |
α j |
2 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
||||||||||
|
Пр 4 |
ri |
α i |
β j |
2 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
||||||||||
|
3 |
Пр 5 |
ri |
β i |
rj |
2 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
|||||||||
|
Пр 6 |
ri |
β i |
α j |
4 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
||||||||||
|
Пр 7 |
ri |
β i |
β j |
2 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
||||||||||
|
4 |
Пр 8 |
α i |
β i |
rj |
2 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
|||||||||
|
Пр 9 |
α i |
β i |
α j |
2 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
||||||||||
|
Пр 10 |
α i |
β i |
β j |
2 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
||||||||||
|
III |
5 |
Пр 11 |
ri |
rj |
rk |
1 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
|||||
|
Пр 12 |
ri |
rj |
α k |
2 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
|||||||
|
Пр 13 |
ri |
rj |
β k |
2 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
|||||||
|
6 |
Пр 14 |
ri |
α j |
α k |
4 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
||||||
|
Пр 15 |
ri |
α j |
β k |
4 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
|||||||
|
7 |
Пр 16 |
α i |
α j |
α k |
8 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
||||||
|
Пр 17 |
α i |
α j |
β k |
8 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
|||||||
|
8 |
Пр 18 |
β i |
β j |
rk |
2 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
||||||
|
Пр 19 |
β i |
β j |
α k |
4 |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
||||||
|
Пр 20 |
β i |
β j |
β k |
∞ |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
x |
|||||||
|
Плоские |
|||||||||||||||||||||||
|
u 1 |
u 2 |
pi |
ϕ i |
pj |
ϕ j |
pil |
ϕ il |
pi2 |
ϕ i2 |
||||||||||||||
|
I |
1 |
Пл 1 |
pi |
ϕ i |
1 |
x |
x |
x |
x |
||||||||||||||
|
II |
2 |
Пл 2 |
1 |
x |
x |
x |
x |
x |
x |
||||||||||||||
|
Пл 3 |
2 |
x |
x |
x |
x |
x |
x |
||||||||||||||||
|
3 |
Пл 4 |
2 |
x |
x |
x |
x |
x |
x |
|||||||||||||||
Примечание: i 1 и i2 - индексы векторов, не имеющих параметров-функций.
Y = j pi ■ sin(фi-ф2), p1=-Y/sin>1-^2) (7)
i=3 затем из уравнения (5)
т.е. перейдем к системе координат X2O2Y2
(рис.1), тогда из уравнения (6) получим p 2 = -X-p 2 ⋅ cos( ϕ 1 - ϕ 2 ) . (8)
Модуль Пл3 (функции p1 и ϕ 2 ). В системе координат X1O1Y1 (рис.2),
X= ∑ m pi ⋅ cos ϕ (i - ϕ 1) , i=3
-
Y= ∑ m pi ⋅ sin ϕ (i- ϕ 1) .
i=3
Первый вариант решения (V=1) принимаем при | ϕ 2 1 | > π /2.
Пусть k = 1 при V = 1 и k = -1 при V = 2.
В треугольнике АВО1 ∆p= p22 -Y 2 , тогда p1 =-X+k ⋅ ∆p, (9)
ϕ 2 = sign(- ⋅ Ya)rccos( ⋅ ∆ kp/p2) + ϕ 1. (10)
Модуль Пл4 (функции ϕ 1 и ϕ 2 ). Предварительно определим параметры суммарного вектора p1-2 = p1 + p2 .
p 1-2 = VX2+Y2 ,
ϕ 1-2 =sign(- ⋅ Ya)rccos(-pX1/-2) , где
X= ∑ m pi ⋅ cos ϕ i, Y= ∑ m pi ⋅ si ϕ ni . i=3 i=3
Перейдем к системе координат X1О1Y1 (рис.3).
В треугольнике со сторонами a, b, c , и противолежащими соответственно углами А, В, С:
А = 2 ⋅ arccos л/ p ⋅ (p - a)/bc
Рис.2. Геометрия модуля Пл3
Рис.3. Геометрия модуля Пл4
и B=2 ⋅ arccosp ⋅ (p-b)/a , где p=0,5⋅ (a + b + c)
или при a = p1, b = p2 и с = p1-2 (рис.3)
А =2 ⋅ arccoss ⋅ (p-p1)/p2 ,
B=2 ⋅ arccoss⋅ (p-p2)/p1 , где p=0,5⋅ (p1 + p2 + p1-2) , s= p/p1-2 .
Принимаем первый вариант решения
-
V=1 при ϕ11-2 ≥ 0.
Пусть k = 1 при V = 1 и k = -1 при V = 2, тогда
ϕi= ϕ1-2 + ∆ϕi, (11)
где
i = 1 или 2, ∆ϕ1 =k⋅ B , ∆ϕ2 =-k⋅ A .
Используя отображения пространственного векторного контура на вспомогательной плоскости, представим решение задачи о функциях большинства элементарных пространственных модулей табл.1 в виде следующих параметрических формул.
Модуль Пр1 (функции r1, α 2, β 1 ) .
p4 =∑m ricosαicosβi , j4 = 0, i=2
p5 =∑m ricosαisinβi , j5 = 90°, i=2 5
m
p2 =∑risinαi , j = 90°.
i = 2 2
Решение единственное.
(V=1) - при β21-2 ≥
π
.
Максимальное число решений ‒ 2.
Модуль Пр2 (функции r 1 , α 1 , r 2 ). m
Модуль Пр7 (функции r1, β 1, β 2 ) .
r 2 =
-
∑ricosαisin(βi -β1)
i = 3
cosα2 ⋅ sin(β2 - β1)
r1
m
∑ri sinαi i=2 ,
sinα1
r 2
m
(1, 2, …, m) Пp1( r 1 , α 1 , β 1 )
p1
p2
= r1 cosα1
= r2 cosα2
m
p3 =∑ricosαicosβi , j3 = 0, i=3
Решение единственное.
Модуль Пр3 (функции r1 , α 1 , α 2 ).
∑m ri cosαi sin(βi -β1)
α 2 = ± arccos i = 3
r 2 ⋅ sin( β 2 - β 1)
α 2
m
(1, 2, …, m) Пp1( r 1 , α 1 , β 1 )
Принимаем: первый вариант решения ( V =1) - при α 2 ≥ 0 .
Максимальное число решений ‒ 2.
Модуль Пр4 (функции r1, α 1, β 2 ).
m
p4 = ∑ ricosαisinβi , j4 = 90°.
i = 3
|
m =4 (1, 2, 3, 4) |
|
|
Пл4( ϕ 1 , ϕ 2 ) |
|
|
V |
|
β 1 = ϕ 1
β 2 = ϕ 2
Принимаем: первый вариант решения ( V =1) - при β 1 2 ≥ 0
Максимальное число решений ‒ 2.
Модуль Пр8 (функции α1, β1, r2). ϕ2 = α2
p3 =∑m ricosαicos(βi-β2), ϕ3 = 0, i=3
m
-∑ri cosαi sin(βi - β1)
β = arcsin i=3 +
r2 ⋅ cosα2
+ β 1
β2 = β или
β 2 = 180 ° - β
m
p4 =∑risinαi, ϕ4 = 0, i=3
p5 = r1 - p4
Принимаем: первый вариант решения
π
( V =1) - при ϕ 5 ≥ 2 .
Модуль Пр18 (функции β 1, β 2, r3 ).
-
r 3
m
r1 sinα1 + r2 sinα2 + ∑ri sinαi i=4
sinα3
Максимальное число решений ‒ 2.
Модуль Пр9 (функции α 1, β 1, α 2 ).
m2 = m
(1, 2, …, m) Пp7( r 1 , β 1 , β 2 )
p2 = r2,
p3 =∑m ricosαicos(βi-β2), j3 = 0, i=3
-
p4 = ∑m ri sinαi , j4 = 0, i=3
p5 = r1 - p4
Принимаем: первый вариант решения ( V =1) ‒ при β 11 - 2 ≥ 0.
Максимальное число решений ‒ 2.
Модуль Пр19 (функции β 1, β 2, α 3 ).
m
-
- ( r 1 sin α 1 + r 2 sin α 2 + ∑ ri sin α i )
α = arcsin i = 4
r 3
α3 = α или
= 180 ° - α
α 3
m 2 = m
(1, 2, …, m) Пp7( r 1 , β 1 , β 2 )
Принимаем: первый вариант решения
Принимаем: первый вариант решения ( V =1) - при α 2 ≥ 0 .
Максимальное число решений ‒ 2.
Модуля Пр10 (функции α 1, β 1, β 2 ).
m
-
- r sin α
ii
α1 = arcsin i=2
r1
α 1
m2 = m
(1, 2, …, m) Пp7( r 1 , β 1 , β 2 )
π
(V=1) ‒ при β11-2 ≥0 и α3 ≥ 2 .
Максимальное число решений ‒ 4.
Прямые аналитические решения для функций модулей Пр12 и Пр13 получены ранее методом Гаусса [1-3], для модулей Пр6, Пр11, Пр14-Пр17, применим метод обращения, при котором одну из функций исходного модуля используем в качестве аргумента обращённого .
Производные первого порядка функций элементарных модулей
Для одноаргументной функции
Принимаем: первый вариант решения ( V =1) - при β 11 - 2 ≥ 0.
Максимальное число решений ‒ 2.
dU X) = dU dx dt ∂x dt ,
для многоаргументной функции
dU(x1,
..
dt
., xk ) = ∑ k ∂ U dx j , j = 1 ∂ xj dt
т.е. для определения полной производной функции необходимо определить все её частные производные первого порядка.
Дифференцируем уравнения (1) и (2)
QU,A,.
= , dxjA где
Таблица 2. Коэффициенты a. и свободные s члены b1s матриц формулы (16) для плоских модулей
d 2 U ( x 1,..., xk ) dt 2
A =
a
a
a
a,
,
A i =
b 1 1
a 2 1
“
a
,
A 2 =
a
a
b 1 1
b 1 2
.
k
=1
j = 1
Для пространственного модуля (i = 1, 2, 3 и s = 1, 2, 3) используем формулу (16) и данные табл.3, при этом
d 2 U ( dXj ) 2 d U d 2 Xj
+ dx 2j dt dx dt2
j V 7 j
+ 21 1 d 2 U dx- dx j~1 -2dj+i dx- dx -2 dt dt
|
a 1 1 |
a 2 1 |
a 3 1 |
“ b i i |
a 2 i |
a 3 1 |
||
|
A= |
a 1 2 |
a 2 2 |
a 3 2 |
, A 1 = |
" Ь 1 2 |
a 2 2 |
a 3 2 |
|
a 1 3 |
a 2 3 |
a 3 3 |
“ Ь 1 з |
a 2 3 |
a 3 3 |
|
a 1 1 |
“ b 1 i |
a 3 1 |
a i i |
a 2 i |
" b 1 i |
||
|
A 2 = |
a 1 2 |
“ Ь 1 2 |
a 3 2 |
, A = |
a 1 2 |
a 2 2 |
- bx 1 2 |
|
a 1 3 |
“ Ь 1 з |
a 3 3 |
a 1 3 |
a 2 3 |
— b 1 3 |
Производные второго порядка функций элементарных модулей
Дифференцируем выражения (12) и (13)
т.е. для определения полной производной второго порядка функции необходимо определить все частные производные первого и второго порядков этой функции, как чистые так и смешанные.
Дифференцируем систему уравнений (14) и (15)
a
d U 1
dx . • dx . j 1 j 2
d U 2 + b 2 = 0, (19)
dx . •dx .
j1
a
d 2U
-----1— + a 2 dx • dx 2 dx • dx j1 j2 j1
= 0, (20)
Таблица 3. Коэффициенты a i и свободные члены b1 s матриц формулы (16) для пространственных модулей
|
U i a>\ |
r k |
a k |
P k |
x k b 1 |
|
a i 1 |
cos a k sin ^ k |
-r k sin a k cos P k |
-r k cos a k sin P k |
b 1 1 |
|
a i 2 |
sin a k |
r k sin a k |
0 |
b 1 2 |
|
a i 3 |
cos a k sin ^ k |
-r k sin a k sin P k |
-r k cos a k cos P k |
b 1 3 |
∂α i ∂ b 1
Таблица 4. Частные производные S и S матриц формулы (17) для плоских модулей
∂x ∂x j2 j2
|
x j 1 |
p k |
ϕ k |
x j 1 |
||||||
|
x j 2 ∂α iS ∂ x j 2 |
p k |
ϕ k |
p l |
ϕ l |
p k |
ϕ k |
p l |
ϕ l |
xj 2 ∂ b 1 S ∂ x j 2 |
|
∂ a i 1 ∂ x j 2 |
0 |
-sin ϕ k |
0 |
0 |
-sin ϕ k |
-p k cos ϕ k |
0 |
0 |
∂ b 11 ∂ x j 2 |
|
∂ a i 2 ∂ x j 2 |
0 |
cos ϕ k |
0 |
0 |
cos ϕ k |
-p k sin ϕ k |
0 |
0 |
∂ b 12 ∂ x j 2 |
тогда при i = 1, 2; s = 1, 2; j1 = 1, …, k; j2 = 1, …, k, используя данные табл.4,
где
∂ 2 U i = ∆ i
∂ xj 1 ∂ xj 2 ∆
b 2 1
= ∑3 ∂αi1 ⋅ ∂Ui + ∂b11 i=1 ∂xj2 ∂xj1 ∂xj2
где
|
∆= |
α 1 1 |
α 2 1 |
, |
- b 21 |
α 21 |
, |
α 1 1 |
- b 21 |
|
α 1 2 |
α 2 2 |
1 |
- b 22 |
α 22 |
2 |
α 1 2 |
- b 22 |
b 2 = ∑ 3 ∂α i 2 ⋅ ∂ U i + ∂ b 12
2 i=1 ∂xj2 ∂xj1 ∂xj2 ,
b2
при
b =∑2 ∂αi1 ⋅∂Ui + ∂b11
21 i = 1 ∂ xj 2 ∂ xj 1 ∂ xj 2
∑ 3 ∂α i 3 i = 1 ∂ xj 2
∂Ui + ∂b12 ∂xj1 ∂xj2 .
Для каждого одно- и двухаргументно-
b =∑2 ∂αi2 ⋅ ∂Ui + ∂b12
22 i=1 ∂xj2 ∂xj1 ∂xj2 .
го модуля табл.1 на основе приведённых выше общих матричных решений могут быть получены аналитические выражения,
Для пространственного модуля (i = 1, 2, 3 и s = 1, 2, 3) используем данные табл.5 и формулу (17), при этом
минимизированнные по числу вычислительных операций.
Аналитическое и программное обеспечение абстрактных векторных модулей, ис-
∆=
∆ 2
α
α
α
α
α
α
α 2 1
α
-
b 2 1
α 2 1
α
α
α
∆
1 =
-
,
b
α
α
α
α
-
b
α
α
,
пользованное при решении различных научных и инженерных задач теории механизмов, может быть рекомендовано к широкому практическому применению, в частности, применительно к системам управления, например робототехнических систем.
- b21
α
αα
1 1 2 1
-
b
α
∆=
-
b
α
,
α
α
α
α
- b21
- b22
- b23