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

UptoLike

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

63
4
I:
=
N
-
1,1,
-
1
Q
=
C[I]
J:
=
I
+
1,N
Q:
=
Q
-
X[J]
*
A[I,J]
X[I]
=
Q/A[I,J]
=
Печать X[I]
Вычисление
вектора M[I],
I:=1
,
k
Конец
Печать вектора
перемещений
Z
64
2.1.3. Программа для расчета рам на языке Turbo Pascal
(*РАСЧЕТ СТАТИЧЕСКИ НЕОПРЕДЕЛИМОЙ РАМЫ М. ПЕ-
РЕМ. В МАТРИЧНОЙ ФОРМЕ*)
PROGRAM RAMA_MP;
CONST
(*ВВОД ИСХОДНЫХ ДАННЫХ*)
(*КОЛИЧЕСТВО НЕИЗВЕСТНЫХ И СЕЧЕНИЙ*)
N=2; M=9; (**)
TYPE
MASS = ARRAY[1..M, 1..N] OF REAL;
MASS1= ARRAY[1..M, 1..M] OF REAL;
MASS2= ARRAY[1..N, 1..M] OF REAL;
MASS3= ARRAY[1..N, 1..N] OF REAL;
MASS4= ARRAY[1..M] OF REAL;
MASS5= ARRAY[1..N] OF REAL;
VAR
B1:MASS; F:MASS1; BT1,C:MASS2; D:MASS3;
S0P,S0P1,SP,CS0P:MASS4; X:MASS5;
I,J,L,KK:INTEGER; Q:REAL;
(*МАТРИЦА ЕДИНИЧНЫХ МОМЕНТОВ*)
CONST
KB1:MASS=
(( -0.667, 0.667), (**)
( 0.333, 0 ), (**)
( 0.333, 0 ), (**)
( 1.333, -0.667), (**)
( -1, 0 ), (**)
( -0.5, 0 ), (**)
( 0, 0 ), (**)
( 0, 0 ), (**)
( 0, -0.333)); (**)
      63                                                        64

                                    2.1.3. Программа для расчета рам на языке Turbo Pascal
           4
                                    (*РАСЧЕТ СТАТИЧЕСКИ НЕОПРЕДЕЛИМОЙ РАМЫ М. ПЕ-
                                    РЕМ. В МАТРИЧНОЙ ФОРМЕ*)
   I:=N-1,1,-1
                                    PROGRAM RAMA_MP;
      Q=C[I]                        CONST
                                          (*ВВОД ИСХОДНЫХ ДАННЫХ*)
     J:=I+1,N
                                    (*КОЛИЧЕСТВО НЕИЗВЕСТНЫХ И СЕЧЕНИЙ*)
Q:=Q-X[J]*A[I,J]                                    N=2; M=9;                (**)
                                    TYPE
                                           MASS = ARRAY[1..M, 1..N] OF REAL;
                                           MASS1= ARRAY[1..M, 1..M] OF REAL;
   X[I]=Q/A[I,J]                           MASS2= ARRAY[1..N, 1..M] OF REAL;
                                           MASS3= ARRAY[1..N, 1..N] OF REAL;
                                           MASS4= ARRAY[1..M] OF REAL;
    I:=1,N                                 MASS5= ARRAY[1..N] OF REAL;
                                    VAR
                   Печать вектора          B1:MASS; F:MASS1; BT1,C:MASS2; D:MASS3;
   Печать X[I]     перемещений        S0P,S0P1,SP,CS0P:MASS4; X:MASS5;
                   Z
                                           I,J,L,KK:INTEGER; Q:REAL;
                                    (*МАТРИЦА ЕДИНИЧНЫХ МОМЕНТОВ*)
                                    CONST
  Вычисление                        KB1:MASS=
  вектора M[I],                                   (( -0.667, 0.667), (**)
      I:=1,k                                       ( 0.333, 0      ), (**)
                                                   ( 0.333, 0      ), (**)
      Конец                                        ( 1.333, -0.667), (**)
                                                   ( -1,      0    ), (**)
                                                   ( -0.5,    0    ), (**)
                                                   ( 0,       0     ), (**)
                                                   ( 0,       0     ), (**)
                                                   ( 0,      -0.333)); (**)