Программа 8. Задача 18

>    restart:#Инварианты приведения системы сил

>    with(LinearAlgebra):with(plots):N:=4:

   Проекции  векторов сил

>    F[1]:=Vector([0,0,1]):

>    F[2]:=Vector([0,-3,0]):

>    F[3]:=Vector([0,5,0]):

>    F[4]:=Vector([2,0,0]):

   Координаты точек приложения

>    A[1]:=Vector([0,3,0]):

>    A[2]:=Vector([1,3,0]):

>    A[3],A[4]:=Vector([0,0,4])$2:

>    A[4]:=A[3]:

>    #ro:=Vector([1,-1,0]):for i to N do A[i]:=A[i]-ro;od:#Для проверки

Главный вектор

>    R:=add(F[i],i=1..N);

R := Vector(%id = 152056508)

>    R0:=Norm(R,Euclidean);

Главный момент

>    M:=add(CrossProduct(A[i],F[i]),i=1..N);

R0 := 3

M := Vector(%id = 152057876)

>    M0:=Norm(M,Euclidean);

Скалярный инвариант

>    Inv:=DotProduct(R,M);

Минимальный момент

>    Mmin:=Inv/R0;

Шаг винта

>    p:=Mmin/R0;

>    V:=Vector([x,y,z]);

Уравнение центральной винтовой оси

>    AX:=evalm(M-CrossProduct(V,R)-R*p);

M0 := 362^(1/2)

Inv := -21

Mmin := -7

p := -7/3

V := Vector(%id = 152058036)

AX := vector([-37/3-y+2*z, 38/3-2*z+x, -2/3-2*x+2*y])

>    Цвет:=[red,blue,black,green]:

>    Стрелка:=width=0.05,head_length=0.2,color=Цвет[i],length=2:

>    Вид:=scaling=constrained,axes=normal,orientation=[65,76],

>    axesfont=[HELVETICA,6],tickmarks=[3,3,4]:

>    display(seq(arrow(A[i],F[i],Стрелка),i=1..N),Вид);

[Maple Plot]

>    for j to 90 do

>    Вид:=axes=normal,orientation=[65+j,76]:

>    P[j]:=display(seq(arrow(A[i],F[i]),i=1..N),Вид):
od:

>    display(seq(P[i],i=1..90),insequence=true);

[Maple Plot]