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

       Условно диск разделим на 3 участка (рис.14.1,в). Число шагов интегрирования на
участках: n1 = 16 , n2 = 96 , n3 = 6 . Шаг интегрирования на каждом участке:
                                       (r − r )    (r − r )     (r − r )
                                 h1 = 3 1 ; h2 = 4 3 ; h3 = 2 4 .
                                          n1          n2           n3
       Внутри шага печати должно умещаться целое число шагов интегрирования.
Программа расчета диска приведена в п.14.3.1. В ней предусмотрен учет изменения
температуры, толщины диска и модуля упругости материала. В программе приняты
следующие обозначения:
       Y(2) – вектор состояния;
       XS(4) – массив радиусов, соответствующих границам участков;
       NN(3) – массив чисел шагов интегрирования на участках;
       KR(3) – массив чисел шагов вывода результатов на участках;
       Y1(2),P0(2),P1(2),P2(2),P3(2) рабочие массивы для подпрограммы RK;
       NP – число шагов интегрирования в шаге вывода результатов;
       KP – число точек вывода результатов на участке;
       х – переменная интегрирования;
       СР – коэффициент Пуассона;
       Р(2) – вектор правых частей;
       RN – число оборотов диска в минуту;
       RO – плотность материала диска;
       OM – угловая скорость;
       N - порядок системы уравнений;
       Y(N) – вектор неизвестных;
       H – шаг интегрирования;
       NI – номер интегрирования;
       NS – номер участка;
       TH(X) - толщина диска в зависимости от радиуса.
       Эпюры u, σ r , σ t постоянные по результатам машинного счета представлены на
рис.14.1,а.

                      14.3.1. Программа расчета вращающихся дисков

      Program DISK;
      Uses Crt;
      Label 1;
      Const
      N=2; RN=3000; RO=8E3; AL=0; Cp=0.3;
      Type
      Wect=array[1..4] of real;
      Var
      Ns,Ni,Ip,I,Np,Kp,:integer;
      H,x,Tx,Thx,Ex,Ur,T2,Sig1,Sig2,C,C2:real;
      Y,XS:wect;
      NN,KR:array[1..3] of integer;
      {Вычисление температуры в зависимости от радиуса}
      Function T(x:real):real;
             Begin
             If Ns=1 then T:150+892.9*(x-0.019);
             If Ns=2 then T:165+3311*(x-0.0358);
             IfNs=3 then T:=510+6667*(x-0.140);T:=0;
             End;



                                         254