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

                          MT:=DZ[i]*CCP*S[2]/RZ[i]+CP*MR;
                          TETA:=S[2];
                          Writeln(RZ[i]:10,``,NZ[i]:10,``.QRZ[i]:10,``,
                                 W:10,``,MR:10,``,MT:10,``,TETA:10);
                      End;
                   End;
                   Begin {Основная программа}
                   KR[1]:4; KR[2]:=8;
                   NN[1]:=16; NN[2]:=80;
                   RG[1]:=0.02; RG[2]:=0.04; RG[3]:=0.125;
                   NU1:=NU+1;
                   CCP:=1-CP*CP;
                   PLAST;
                   Readln;
                   END.
       Программа состоит из девяти программных единиц.
       В основную программу вводятся следующие данные о пластине:
       NU – количество участков;
       NU1 – количество границ участков (NU1=NU+1);
       RG(NU1) – массив радиусов (в метрах) границ участков, включая внутренний и
внешний радиусы пластины;
       NN(NU) – массив числа шагов интегрирования на каждом участке;
       KR(NU) – массив числа шагов печати по участкам; при этом каждый шаг печати
должен включать целое число шагов интегрирования, поэтому каждый элемент массива NN
должен нацело делиться на соответствующий элемент массива KR (при вводе исходных
данных это условие нужно выполнить).
       Основная программа обращается к подпрограмме PLAST, реализующей алгоритм
расчета пластины методом начальных параметров. В описательной части подпрограммы
PLAST задаются следующие массивы:
       Y –матрица решений Y = ( y1 , y0 ) ;
       Y1,P0,P1,P2,P3 – рабочие массивы подпрограммы интегрирования размерностью 3х2;
       RZ,HZ,DZ,QRZ – одномерные массивы, в которых запоминаются значения радиуса,
толщины, жесткости на изгиб и поперечной силы в точках вывода информации;
       Z – трехмерный массив, предназначенный для запоминания матрицы решений в тех
же точках;
       S – одномерный рабочий массив, используемый при обработке результатов.
       При этом максимально возможное число точек вывода не должно превышать 50, в
противном случае нужно изменить размеры следующих массивов: RZ,HZ,DZ,QRZ,Z.
       Вначале вводят значения следующих величин:
       M, N – число столбцов и строк матрицы Y (M=2; N=3);
       K, K1 – начальные значения счетчиков точек вывода результатов и общего числа
точек вывода.
       После обращения в подпрограмме NUS(Y), формирующей начальное значение
матрицы Y ( r1 ) , где r1 - левая граница интервала интегрирования, вычисляется общее число
точек вывода результатов счета.
       Далее имеем три вложенных цикла:
                   1) внешний цикл – по номеру участка NS,
                   2) промежуточный цикл – по номеру точки вывода IP,
                   3) внутренний цикл – по номеру шага интегрирования IK.
       Во внешнем цикле вычисляются значение радиуса R внутреннего контура участка,
шаг интегрирования H, число шагов интегрирования NP в шаге печати I и число точек
вывода результатов KP на участке; затем осуществляется переход к промежуточному циклу.


                                            264