Программа 36. Задача 50
> | restart; n:=4: # Число участков |
> | m:=0,10,-18,15:# Моменты (кН) |
> | L:=2,4,3,5: # Длины участков(м) |
> | dm:=0.09,0.18,0.25,0.32:#Диаметры (м) |
> | G:=80e9: # Модуль сдвига(Па) |
> | M[0]:=-Ma: f[0]:=0: x[0]:=0: |
> | for i to n do |
> | Ji:=Pi*dm[i]^4/32: |
> | Wi:=2*Ji/dm[i]: |
> | M[i]:=M[i-1]+m[i]*1e3: |
> | ffi:=M[i]*L[i]/G/Ji; |
> | f[i]:=f[i-1]+ffi: |
> | x[i]:=x[i-1]+L[i]: |
> | tau[i]:=M[i]/Wi/1e6;# МПа |
> | od: |
> | Ma:=solve(f[n],Ma);#Момент на опоре (Нм) |
> | plot([seq([x[i],f[i]],i=0..n)],title=`fi`); |
> | plot([seq(op([[x[i],tau[i+1]],[x[i+1],tau[i+1]]]),i=0..n-1)],title=`tau`); |
> | plot([seq(op([[x[i],M[i+1]],[x[i+1],M[i+1]]]),i=0..n-1)],title=`M`); |
> | evalf(seq(tau[i],i=1..n)); # касательные напряжения (МПа) |
> | evalf(seq(f[i]*1000,i=1..n-1));# углы*1000 |
> | evalf(seq(M[i],i=1..n)); # моменты (Нм) |