Руководство к решению задач по механике материалов и конструкций. Егодуров Г.С - 222 стр.

       WP – полярный момент сопротивления сечения Wp;
       M – внутренний крутящий момент T;
       M1 – внешний крутящий момент M;
       M2 – интенсивность внешнего распределенного крутящего момента te;
       l1, l 2 - длина соответственно первого и второго участков стержня;
       l - вся длина стержня;
       D1 – диаметр стержня d1;
       D2 – диаметр стержня d2;
       H – шаг интегрирования ∆z;
       N – число шагов интегрирования n;
       d – диаметр текущего поперечного сечения стержня d;
            Результаты расчета представлены графически на рис. 11.1,б.

            11.1.2 Программа расчета стержня переменной жесткости на кручение

Program VAL; {Вал переменной жесткости}
Uses CRT;
Const kl=50;
Var i,n,n1:integer;
        l,ll,l2,d,d1,d2,M1,M2,G,z,h,Ip,zl:real;
        T,F,K,M,Wp,ZZ:array[1..k1] of real;
        Otlet:char;
Begin
Clrscr;
Writeln(‘Размер вала’);
Writeln(‘l(m),l1(m),l2(m),d1(m),d2(m)’);
Readln(l,l1,l2,d1,d2);
Writeln(‘Нагрузка M1(H*m),M2(H*m)’);
Readln(M1,M2);
Writeln(‘Модуль сдвига G(Па’);
Readln(G);
Writeln(‘Число шагов интегрирования n’);
Readln(n);
Writeln(‘Результаты расчета’);
Writeln;
Writeln(‘z mm MH*m T MПа                 F рад   Wp см3);
Writeln;
h:=l/n;
n1:=n+1; {Количество узлов}
For i:=1 to n1 do
Beqin
z:=h*(i-1);
ZZ[i]:=z;
If (z>=0) and (z<=l1) then M[i]:=-M1;
If (z>l1) and (z<=l2) then M[i]:=-M1+M2*(z-l1);
If (z>l2) and (z<=l) then M[i]:=-M1+M2*(l2-l1);
If (z>=0) and (z<=l1) then d:=d1;
If (z>l1) and (z<=l2) then d:=((d2-d1)*z+(d1*l2-d2*l1))/(l2-l1);
If (z>l2) and (z<=l) then M[i]:=-d:=d2;
Ip:=3.14*exp(4*ln(d))/32;
Wp[i]:=2*Ip/d;
T[i]:=M[i]/Wp[i];



                                                222