1.МАТЕМАТИЧЕСКИЕ СВЕДЕНИЯ

1.1.Матрицы и векторы

Матрицей называется квадратная или прямоугольная таблица, запол-

ненная числами, либо - постоянными или переменными символами.

В основном мы будем рассматривать квадратные матрицы и

матрицы-строки

вида

[a,b,c,d ] =[a b c d ], а также матрицы-столбцы,

T

 

{a b c d} =[a b c d]

 

=[a b c d] , получаемые в результате операции транспо-

нирования матриц-строк и отмечаемые фигурными скобками. Например, ра- диус-вектору r (x, y, z) и любым другим трехкомпонентным векторам, имеющим геометрический образ и правило проектирования на различные

оси, сопоставляется вектор-строка (матрица-строка)

r = [x y z]

или вектор-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

столбец {r} ={x y z} = [x y z] проекций.

 

 

 

 

 

 

 

 

Скалярное произведение двух векторов a и

 

 

вычисляется как произ-

 

b

ведение строки a на столбец {b} =bT :

 

 

 

 

 

 

 

 

 

 

 

 

 

c = a

 

= a{b} = ax bx

+ ay by + az bz .

 

 

 

 

 

 

(1.2)

 

 

 

b

 

 

 

 

 

 

 

Вектор-строки векторного произведения c = a ×

 

= −

 

×a

 

 

 

b

b

 

 

определяется формулой

%

 

 

 

 

 

 

 

 

 

(1.3)

 

 

 

%

 

 

 

 

 

 

 

 

 

 

~

 

 

 

c =ba = −ab .

 

 

 

 

 

 

 

 

 

 

Здесь a — кососимметрическая матрица третьего порядка, вида

 

 

~

0

az

ay

 

 

 

 

 

 

 

 

 

 

(1.4)

 

0

 

ay ;

az 0 ax ; ay ax

0]

a

= − az

ax = [0 az

 

 

 

ax

 

 

 

 

 

 

 

 

 

 

 

 

 

ay

0

 

 

 

 

 

 

 

 

 

 

 

Итак, вектор-строка векторного произведения равна строке b второго сомножителя, умноженной на кососимметрическую матрицу a% первого сомножителя, либо взятому со знаком минус произведению строки a первого со-

множителя и кососимметрической матрицы b% второго сомножителя. В MATLAB вектор-строка c векторного произведения a ×b возвращается функ-

цией cross вида c=cross(a,b).

В результате транспонирования выражений (1.3) и учета свойства ко-

~

~

находим выражения

для векторного

сосимметрической матрицы a

′ = −a ,

произведения в виде вектор-столбца

~

 

 

~

(1.5)

{c} = −a{b} = b{a}.

Квадрат кососимметрической матрицы есть симметрическая матрица

7

%

 

 

a2y az2

axay

axaz

 

 

(1.6)

2

 

ayax

2

 

2

ay az

 

=

a

 

=

ax

az

 

 

 

 

 

a

a

x

a

a

y

a2

a2

 

 

 

 

 

 

z

 

z

 

x

y

 

 

 

=[a2y az2 ,axay ,axaz ;ay ax , ax2 az2 ,ayaz ;az ax ,az ay , ax2 a2y ]

Вэтих выражениях применен формат записи матриц в пакете МatLab, где строки отделяются знаком (;), а элементы в строке отделены пробелами или запятыми.

Библиотека M-функций в механике.

Кососимметрическую матрицу a размера (3 3) можно сформировать в

~ ×

среде МatLab на основе функции cross. А именно, для строки a =[ax ,ay ,az ] создадим M-функцию s33(a), возвращающую (3×3) матрицу вида

function s = s33(a)

(1.7)

s=[ cross(a,[1 0 0]) ; cross(a,[0 1 0]) ;

 

cross(a,[0 0 1])];

 

~2

 

М-функция ss33(a) формирует симметрическую (3×3)-матрицу a

 

function s = ss33(a)

(1.8)

s = s33(a) s33(a);

 

М-функция sp33(a,b) функция двух векторов-строк a , b (угловые скорость и ускорение вращающегося тела: a =ω, b =ε )

function s = sp33(a, b)

(1.9)

s = ss33(a) + s33(b);

 

M-функции для пространственной системы сил.

В кинетостатике и статике при вычислении бивектора W силы

F=[X Y Z] можно применять (3×6)-матрицу радиус-вектора r =[x, y, z], опре-

деляемую М-функцией s36(r) , W=F*s36(r)

 

function s =

s36(r)

(1.10)

s = [s33(r),

eye(3)];

 

где eye(3) = diag[1 1 1] — единичная диагональная матрица третьего порядка.

M-функция координат r = [x y z] точки приложения для силы F = nF e , заданной её модулем и ортом e , где nF=| F | =norm(F), e = [ex ey ez]

function s

= s16(r, e)

(1.11)

s=e*s36(r)

 

 

M-функция координат для силы F = [0 Y Z], вводимой без нуля как, F = [Y Z]

function s = sx26(r)

(1.12)

s=s36(r)

s=s(1, :) = [ ];

8

M-функция координат r = [x y z] силы F = [X 0 Z], вводим как F = [X Z]

function s = sy26(r)

(1.13)

s=s36(r)

s=s(2, :) = [ ];

M-функция координат силы F = [X Y 0], вводим как F = [X Y]]

function s = sz26(r)

(1.14)

s=s36(r)

s=s(3, :) = [ ];

M-функция координат для силы F = [X 0 0] ,вводим как F = X , без нулей

function s = syz16(r)

(1.15)

s=s36(r) s=s(1, :) ;

M-функция координат для силы F = [0 Y 0], вводим как F = Y

function s = sxz16(r)

(1.16)

s=s36(r) s=s(2, :) ;

M-функция координат для силы F = [0 0 Z], вводим как F = Z

function s = sxy16(r)

(1.17)

s=s36(r) s=s(3, :);

M-функция для момента пары сил с моментом M = M e = nMe , e =[ex ey ez ]– вектор-строка орта e . M = nM*e, бивектор пары W = nM*sp16(e)

function s

= sp16(e)

(1.18)

s

= [ex

ey

ez

0 0 0]

 

%

или s

= [e,

zeros(1, 3)];

 

M-функции для плоской системы сил, расположенной в плоскости Oxy

Вычисление (2×3)-матрицы координат для случая плоской системы сил, когда r=[x y], F = [X Y], W = F*s23(r) - бивектор силы F

function s

=

s23(r)

(1.19)

s = [0 -1;

1

0]*[r;[0 –1]; [1 0]]’;

 

Бивектор скалярного момента пары сил M Z = ± | M |, W = nM [±1,0,0]

function s

=

sp13(M)

(1.20)

s = M * [1

0

0];

 

Бивектор единичной силы e , заданной строкой e = [ex ey], при r = [x y] function s = s13(e)

s = e * s23(r);

9

1.2.Конкатенация матриц и линейных функций

Вертикальная конкатенация (сцепление) двух матричных выражений, содержащих один и тот же столбцевой множитель {X}, вида A1{X} и A2{X}, где A1, A2 — прямоугольные матрицы с одинаковым числом столбцов, приводит к сцепленному по вертикали выражению A{X}, где матрица A = [A1; A2]

— результат вертикальной конкатенации матриц A1 и A2.

Горизонтальная конкатенация двух выражений XB1 и XB2, где X есть n-строка, B1 и B2 — прямоугольные матрицы с n-строками, n1 и n2-столбцами, приводит к объединенному выражению XB, где B = [B1, B2] — nx(n1+n2)- матрица, результат горизонтальной конкатенации матриц.

Любую сумму произведений вида ax +by +cz + ku можно представить как произведение вектора-строки (матрицы-строки) [a,b,c,k] образованной из первых сомножителей и вектора-столбца [x, y, z,u] ' вторых сомножителей

или наоборот, — строки вторых сомножителей, умноженной на столбец первых сомножителей. Например, возьмем три суммы произведений и каждую из них представим в виде произведения строки на столбец:

a1x +b1 y +c1z + k1u =[x, y, z,u][a1,b1,c1,k1 ]'

(1.21)

a2 x +b2 y +c2 z + k2u =[x, y, z,u][a2 ,b2 ,c2 ,k2 ]'

a3 x +b3 y + c3 z + k3u =[x, y, z,u][a3 ,b3 ,c3 ,k3 ]'

Три линейных алгебраических выражения объединяем по горизонтали в одно матричное строчное выражение VA, где V=[ x y z u], A=[a1 a2 a3 ; b1 b2 b3 ; c1 c2 c3 ; k1 k2 k3]. Либо сцепляем вертикально в столбцевое выражение A΄{V}

при {V}=V΄.

Сумма произведений строчных матриц и согласованных с ними по умножению прямоугольных матриц объединяется в произведение удлиненной вектор строки и блочной матрицы, полученной в результате вертикальной конкатенации прямоугольных матриц. Например 1x3 – вектор-строка вида:

[x, y, z]A33 +[u, v]B23 + wC13 =[x, y, z, u, v, w][A33; B23; C13 ],

(1.22)

где вторые матричные сомножители имеют размеры (3x3), (2x3), (1x3).

1.3.Перепроецирование векторов и тензоров второго ранга

Пусть имеется две декартовы системы координат: Oxnynzn и Oxsyszs, первую назовем — n-системой или исходной системой, вторую — повернутой системой или — s-системой, с отсчетом углов от исходной системы. Расмотрим в них один и тот же вектор r , он определяется строками проекций, которые обозначим соответственно rn=[xn yn zn] и rs=[xs ys zs].

Верна матричная формула перепроецирования вектора r из n-системы в s-систему, связывающая его проекции на оси s-системы, n-системы:

10

rs = rnCsn

(1.23)

Здесь Csn — матрица поворота s-системы от n-системы, элементами которой являются направляющие косинусы осей повернутой системы, т.е. скалярные произведения соответствующих ортов осей,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

is in

 

 

js in

 

ks in

in

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.24)

 

 

 

 

 

 

 

 

js jn

 

 

 

 

 

 

 

 

[is

js ks ]

Csn = is jn

 

ks jn

jn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

js kn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

is kn

ks kn

kn

 

 

 

 

 

 

 

 

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

Мнемоническое правило: в формуле (1.23) индексы s, n повторяются в том же порядке, а в произведении матриц (1.24) – в обратном порядке.

Общая формула перепроектирования любого вектора a

из n-системы на

s-систему и обратно

(1.25)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

as = anCsn , an = asCns при Сns =[is js ks ]' [in jn kn ] = Csn

 

Примем, например, за исходную систему Ox1y1z1 и за повернутую — “нулевую” систему Oxyz. Тогда по формуле (1.23) находим строку r=[x y z]

проекций вектора r

на “нулевую” систему:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

r = rC

при

С

=[i

 

j

 

 

 

 

 

 

 

 

 

j k ]

 

 

k ] [i

 

1

01

 

01

1

 

1

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

И наоборот, перепроецирование r

из системы Oxyz в Ox1y1z1 имеет вид:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(1.26)

r1 = rC10

при

С10 =[i

 

j

 

 

 

 

[i1

 

j1 k1]

 

 

k ]

 

 

Пусть например, имеется четыре системы координат, из них вторая Ox2y2z2 повернута от первой Ox1y1z1 на угол ψ вокруг Oz1 , третья от второй повернута на угол θ вокруг Ox2 , четвертая от третьей на угол ϕ вокруг Oz3.

Векторы-строки вектора r на эти оси обозначим rS = [xS yS zS ] , s =1,2,3,4 . По формулам (1.24), (1.23) находим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cψ sψ

0

 

 

 

 

 

 

i1

 

 

 

 

 

 

 

 

 

 

 

 

 

C

 

(ψ, z) =

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= sψ

cψ

0

,

r

= rC

21

 

j

 

 

i

 

j

 

k

2

 

 

 

1

 

2

2

 

 

 

 

 

 

 

2

1 21

 

 

 

 

k

 

 

 

 

 

 

 

 

 

 

 

 

0

0

1

 

 

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

i

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

C32 (θ, x) = j2

 

i3

 

 

 

 

 

 

 

 

 

 

 

k2

 

 

 

 

 

 

cϕ sϕ

C

43

(ϕ, z) = sϕ cϕ

 

 

 

0

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

0

0

 

 

 

 

(1.27)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j3

 

 

=

 

cθ sθ

, r3 = r2C32

k3

0

 

 

 

 

 

 

 

 

0

 

sθ

cθ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

, r

= r C

43

= rC C

32

C

43

,cψ cosψ, sψ sinψ

 

 

4

 

3

1

21

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Отсюда находим матрицу поворота C41 четвертой системы от первой системы и матрицу C14 обратного поворота.

11