Программа 35. Задача 49. Кручение

[Maple OLE 2.0 Object]

>    restart;

>    n:=3:#Число участков

>    d:=0.09,0.14,0.08,0.07: # Диаметры, м

>    G:=8*10^10:             # Модуль сдвига

>    M0:=30,-10,20:          # Моменты, кНм

>    l:=1.1,1.2,0.8:         # Длины, м

>    for i to 2 do           # Моменты инерции, м4

>     J[i]:=evalf(Pi*d[i]^4/32);

>    od;

J[1] := .6441246688e-5

J[2] := .3771481981e-4

>    J[3]:=evalf(Pi*(d[3]^4-d[4]^4)/32);

J[3] := .1664062359e-5

>    for i to n do # Моменты в сечениях

>     M[i]:=add(M0[k],k=i..n);

>    od;

M[1] := 40

M[2] := 10

M[3] := 20

>    for i to n do # Углы относительные

>     ft[i]:=M[i]*1000*l[i]/G/J[i];

>    od;

ft[1] := .8538719702e-1

ft[2] := .3977216403e-2

ft[3] := .1201878036

>    for i to n do # Углы абсолютные

>     f[i]:=add(ft[k],k=1..i);

>    od;

f[1] := .8538719702e-1

f[2] := .8936441342e-1

f[3] := .2095522170

>    for i to n do

>     W[i]:=J[i]/(d[i]/2); #Моменты сопротивления м3

>     tau[i]:=M[i]*1e3*d[i]/J[i]/2/1e6;#Кас.напр.,MPa

>    od;

W[1] := .1431388153e-3

tau[1] := 279.4490084

W[2] := .5387831402e-3

tau[2] := 18.56034322

W[3] := .4160155898e-4

tau[3] := 480.7512144

>