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

       После определения постоянных интегрирования С1, С2 вектор состояния в
произвольном сечении r вычисляют по формуле (15.10). Решения                      y 1 ( r ), y0 ( r )
интегрированием уравнения (15.8), причем, при определении y1 ( r ) следует положить
g = 0 , либо находить эти два решения параллельно интегрированием следующего
дифференциального уравнения:
                                                 dY                                   (15.11)
                                                     = FY + G ,
                                                  dr
где Y ( r ) = ( y1 ( r ), y0 ( r )) матрица частных решений; G = ( 0, g ) - матрица, первый столбец
которого нулевой, а вторым столбцом является вектор g .
      Заметим, что вектор частного решения y0 ( r ) - это второй столбец матрицы Y.
Поэтому в программных обозначениях компоненты этого вектора имеют вторым индексом
цифру 2 (т.е. y12, y22, y32).
      В приведенной программе интегрируется матричное уравнение (15.11).

                                      15.2. Описание программы

       Для численного расчета пластины методом начальных параметров разработана
программа (см. п. 15.2.1), позволяющая рассчитать пластину с одной кольцевой опорой при
осесимметричной поперечной нагрузке.
       При расчете пластина разбивается на участки. Участки отличаются друг от друга либо
различным законом изменения толщины, либо различным законом изменения поперечной
силы, либо и тем и другим.

       15.2.1. Программа расчета осесимметрично нагруженной одной кольцевой опорой

       Program KRUG_PLAST;
       Uses Crt;
       Const
       CP=0.3; E=2E11; NU=2; N=3; M=2;
       TYPE
       Matr:array[1..3, 1..3] of real;
       Amatr:array[1..50, 1..3, 1..2] of real;
       Var
       NU1, NS : integer;
       CCP, H, RR: real;
       NN, KR : array[1..N] of integer;
       RG: array[1..NU+1] of real;
       {Функция для вычисления толщины}
       Function DD (R;real):real;
       Begin
                DD:=E*EXP(3/LN(TOL(NS,R)))/(12/ccP);
       End;
       {Процедура общего вида начальных значений векторов}
       Procedure NUS (Var Y:Matr);
       Var
       i,j: integer;
       Begin
                For i:=1 to 3 do
                For j:=1 to 2 Y[I,j]:=0;
                Y[2,1]:=1;


                                                261