Руководство к решению задач по механике твердого деформируемого тела матричными методами. Бундаев В.В. - 41 стр.

UptoLike

Составители: 

81
Рис.2.5
2.2.3. Программа для расчета статически
неопределимых ферм
B
Решение СЛАУ
CXA
ρ
ρ
=
печать вектора
реакций
X
ρ
i=1,k
C[i]:=0
j=1,n
C[i]:=C[i]+N0[i,j]*X[j]
N[i]:=C[i]+NP[i]
выво
д
N
ρ
конец
XNC
ед
ρ
ρ
=
Pед
NXNN
ρ
ρ
ρ
+=
N
ρ
вектор результи-
рующих продольных сил
82
В соответствии с блок-схемой, изображенной на рис.2.5
составляем программу на языке Turbo Pascal.
PROGRAM FERMA_MS;
CONST
(*ВВОД ИСХОДНЫХ ДАННЫХ*)
(*КОЛИЧЕСТВО СТОЛБЦОВ И СТЕРЖНЕЙ*)
N=1; K=13;
TYPE
MASS = ARRAY[1..K, 1..N] OF REAL;
MASS1= ARRAY[1..K, 1..K] OF REAL;
MASS2= ARRAY[1..N, 1..K] OF REAL;
MASS3= ARRAY[1..N, 1..N] OF REAL;
MASS4= ARRAY[1..K] OF REAL;
MASS5= ARRAY[1..N] OF REAL;
VAR
N0:MASS; D:MASS1; NT,H:MASS2; A:MASS3;
NP,N,C:MASS4; X:MASS5;
I,J,L,KK:INTEGER; Q:REAL;
(*ВЕКТОР ПРОДОЛЬНЫХ СИЛ ОТ ЕДИНИЧНОЙ НАГРУЗКИ
В ОСН. ФЕРМЕ*)
CONST
KN0:MASS=((-0.707),
(-1 ),
(-1 ),
(-0.707),
( 0.5),
( 0.5),
( 0.5),
( 0.5),
( 0 ),
( 0.707),
( 0 ),
(0.707),
(0));
                                81                                                             82


                 B
                                                                       В соответствии с блок-схемой, изображенной на рис.2.5
                                                                  составляем программу на языке Turbo Pascal.
           Решение СЛАУ
              ρ ρ
            A⋅X =C                                                PROGRAM FERMA_MS;
                                                                  CONST
                                                                         (*ВВОД ИСХОДНЫХ ДАННЫХ*)
           печать вектора
                      ρ                                           (*КОЛИЧЕСТВО СТОЛБЦОВ И СТЕРЖНЕЙ*)
           реакций    X                                                         N=1; K=13;
                                                                  TYPE
                                                                        MASS = ARRAY[1..K, 1..N] OF REAL;
                 i=1,k                                                  MASS1= ARRAY[1..K, 1..K] OF REAL;
                                                                        MASS2= ARRAY[1..N, 1..K] OF REAL;
                C[i]:=0                                                 MASS3= ARRAY[1..N, 1..N] OF REAL;
                                                                        MASS4= ARRAY[1..K] OF REAL;
                 j=1,n                                                  MASS5= ARRAY[1..N] OF REAL;
                                          ρ        ρ              VAR
      C[i]:=C[i]+N0[i,j]*X[j]             C = N ед X                    N0:MASS; D:MASS1; NT,H:MASS2; A:MASS3;
                                                                  NP,N,C:MASS4; X:MASS5;
                                                                        I,J,L,KK:INTEGER; Q:REAL;
                                                                  (*ВЕКТОР ПРОДОЛЬНЫХ СИЛ ОТ ЕДИНИЧНОЙ НАГРУЗКИ
                                          ρ        ρ ρ
                                          N = N ед X + N P        В ОСН. ФЕРМЕ*)
        N[i]:=C[i]+NP[i]
                                                                  CONST
                      ρ                   ρ                         KN0:MASS=((-0.707),
              вывод   N                   N  – вектор результи-                       (-1 ),
                                          рующих продольных сил                       (-1 ),
                                                                                      (-0.707),
                                                                                      ( 0.5),
                 конец                                                                ( 0.5),
                                                                                      ( 0.5),
                                                                                      ( 0.5),
                                Рис.2.5                                               ( 0 ),
                                                                                      ( 0.707),
                                                                                      ( 0 ),
2.2.3. Программа для расчета статически                                               (0.707),
неопределимых ферм                                                                    (0));