Программа 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); |
> | R0:=Norm(R,Euclidean); |
Главный момент
> | M:=add(CrossProduct(A[i],F[i]),i=1..N); |
> | 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); |
> | Цвет:=[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),Вид); |
> | 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); |