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

UptoLike

ВУЗ: 

ВСГУТУ | Улан-Удэ

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

Рубрика: 

65
(*МАТРИЦА ПОДАТЛИВОСТЕЙ*)
KF:MASS1=
((0.5, 0.25, 0, 0, 0, 0, 0, 0, 0 ), (**)
(0.25, 0.5, 0, 0, 0, 0, 0, 0, 0 ), (**)
(0, 0, 0.5, 0.25, 0, 0, 0, 0, 0 ), (**)
(0, 0, 0.25, 0.5, 0, 0, 0, 0, 0 ), (**)
(0, 0, 0, 0, 0.5, 0, 0, 0, 0 ), (**)
(0, 0, 0, 0, 0, 2, 0, 0, 0 ), (**)
(0, 0, 0, 0, 0, 0, 0.5, 0, 0 ), (**)
(0, 0, 0, 0, 0, 0, 0, 1, 0.5), (**)
(0, 0, 0, 0, 0, 0, 0, 0.5, 1 )); (**)
(*МАТРИЦА ГРУЗОВЫХ МОМЕНТОВ*)
KS0P:MASS4= (1.125,-1.125,-1.125,1.125, 2.25, -1.125, 0, 0, 0);
(**)
(*МАТРИЦА ГРУЗОВЫХ МОМЕНТОВ В
СТАТ.ОПРЕД.СИСТЕМЕ*)
KS0P1:MASS4=(4.5, 0, 0, 0, 0, -2.25, 0, 0, 0); (**)
BEGIN
(*ВВОД МАТРИЦЫ ЕДИНИЧНЫХ МОМЕНТОВ*)
B1:=KB1;
WRITELN('МАТРИЦА ЕДИНИЧНЫХ МОМЕНТОВ');
FOR I:=1 TO M DO BEGIN
FOR J:=1 TO N DO WRITE(' ', B1[I,J]:11:6); WRITELN
END;
WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ <ENTER>');
READLN;READLN;
(*ВВОД МАТРИЦЫ ПОДАТЛИВОСТЕЙ*)
F:=KF;
WRITELN('МАТРИЦА ПОДАТЛИВОСТЕЙ');
FOR I:=1 TO M DO BEGIN
FOR J:=1 TO M DO WRITE(' ', F[I,J]:6:3); WRITELN
END;
WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ <ENTER>');
READLN;READLN;
(*ВВОД МАТРИЦЫ ГРУЗОВЫХ МОМЕНТОВ*)
S0P:=KS0P;
66
WRITELN('МАТРИЦА ГРУЗОВЫХ МОМЕНТОВ');
FOR I:=1 TO M DO WRITE(' ', S0P[I]:6:3); WRITELN;
WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ <ENTER>');
READLN;READLN;
(*ВВОД МАТРИЦЫ ГРУЗОВЫХ МОМЕНТОВ*)
S0P1:=KS0P1;
WRITELN('МАТРИЦА ГРУЗОВЫХ МОМЕНТОВ В
СТАТ.ОПРЕД.СИСТЕМЕ');
FOR I:=1 TO M DO WRITE(' ', S0P1[I]:6:3); WRITELN;
WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ <ENTER>');
READLN;READLN;
FOR I:=1 TO M DO FOR J:=1 TO N DO BT1[J,I]:=B1[I,J];
(*ПЕРЕМНОЖЕНИЕ МАТРИЦ*)
FOR I:=1 TO N DO
FOR J:=1 TO M DO
FOR L:=1 TO M DO
C[I,J]:=C[I,J]+BT1[I,L]*F[L,J];
FOR I:=1 TO N DO
BEGIN
FOR J:=1 TO N DO
FOR L:=1 TO M DO
D[I,J]:=D[I,J]+C[I,L]*B1[L,J];
FOR L:=1 TO M DO
CS0P[I]:=CS0P[I]+C[I,L]*S0P1[L]
END;
(*РЕШЕНИЕ СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ
УРАВНЕНИЙ*)
FOR I:=1 TO N-1 DO
FOR J:=I+1 TO N DO
BEGIN
D[J,I]:=-D[J,I]/D[I,I];
FOR KK:=I+1 TO N DO
D[J,KK]:=D[J,KK]+D[J,I]*D[I,KK];
CS0P[J]:=CS0P[J]+D[J,I]*CS0P[I]
END; 
                               65                                                           66

(*МАТРИЦА ПОДАТЛИВОСТЕЙ*)
KF:MASS1=                                                          WRITELN('МАТРИЦА ГРУЗОВЫХ МОМЕНТОВ');
        ((0.5, 0.25, 0, 0, 0, 0, 0, 0, 0 ),             (**)       FOR I:=1 TO M DO WRITE(' ', S0P[I]:6:3); WRITELN;
         (0.25, 0.5, 0, 0, 0, 0, 0, 0, 0 ),             (**)       WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ ');
         (0,    0,   0.5, 0.25, 0, 0, 0, 0, 0 ),         (**)     READLN;READLN;
         (0, 0,      0.25, 0.5, 0, 0, 0, 0, 0 ),         (**)     (*ВВОД МАТРИЦЫ ГРУЗОВЫХ МОМЕНТОВ*)
         (0, 0,      0,    0, 0.5, 0, 0, 0, 0 ),         (**)      S0P1:=KS0P1;
         (0, 0,      0,    0, 0, 2, 0, 0, 0 ),           (**)      WRITELN('МАТРИЦА ГРУЗОВЫХ МОМЕНТОВ В
         (0, 0,      0,    0, 0, 0, 0.5, 0, 0 ),          (**)    СТАТ.ОПРЕД.СИСТЕМЕ');
         (0, 0,      0,    0, 0, 0, 0, 1, 0.5),           (**)     FOR I:=1 TO M DO WRITE(' ', S0P1[I]:6:3); WRITELN;
         (0, 0,      0,    0, 0, 0, 0, 0.5, 1 ));         (**)     WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ ');
(*МАТРИЦА ГРУЗОВЫХ МОМЕНТОВ*)                                     READLN;READLN;
KS0P:MASS4= (1.125,-1.125,-1.125,1.125, 2.25, -1.125, 0, 0, 0);   FOR I:=1 TO M DO FOR J:=1 TO N DO BT1[J,I]:=B1[I,J];
(**)                                                              (*ПЕРЕМНОЖЕНИЕ МАТРИЦ*)
(*МАТРИЦА ГРУЗОВЫХ МОМЕНТОВ В                                      FOR I:=1 TO N DO
СТАТ.ОПРЕД.СИСТЕМЕ*)                                                FOR J:=1 TO M DO
KS0P1:MASS4=(4.5, 0, 0, 0, 0, -2.25, 0, 0, 0);    (**)                    FOR L:=1 TO M DO
BEGIN                                                                           C[I,J]:=C[I,J]+BT1[I,L]*F[L,J];
(*ВВОД МАТРИЦЫ ЕДИНИЧНЫХ МОМЕНТОВ*)                                FOR I:=1 TO N DO
 B1:=KB1;                                                           BEGIN
 WRITELN('МАТРИЦА ЕДИНИЧНЫХ МОМЕНТОВ');                                   FOR J:=1 TO N DO
 FOR I:=1 TO M DO BEGIN                                                         FOR L:=1 TO M DO
        FOR J:=1 TO N DO WRITE(' ', B1[I,J]:11:6); WRITELN                       D[I,J]:=D[I,J]+C[I,L]*B1[L,J];
END;                                                                      FOR L:=1 TO M DO
 WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ ');                                    CS0P[I]:=CS0P[I]+C[I,L]*S0P1[L]
READLN;READLN;                                                      END;
(*ВВОД МАТРИЦЫ ПОДАТЛИВОСТЕЙ*)                                    (*РЕШЕНИЕ СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ
 F:=KF;                                                           УРАВНЕНИЙ*)
 WRITELN('МАТРИЦА ПОДАТЛИВОСТЕЙ');                                 FOR I:=1 TO N-1 DO
 FOR I:=1 TO M DO BEGIN                                             FOR J:=I+1 TO N DO
        FOR J:=1 TO M DO WRITE(' ', F[I,J]:6:3); WRITELN                  BEGIN
END;                                                                            D[J,I]:=-D[J,I]/D[I,I];
 WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ ');                                    FOR KK:=I+1 TO N DO
READLN;READLN;                                                                    D[J,KK]:=D[J,KK]+D[J,I]*D[I,KK];
(*ВВОД МАТРИЦЫ ГРУЗОВЫХ МОМЕНТОВ*)                                              CS0P[J]:=CS0P[J]+D[J,I]*CS0P[I]
 S0P:=KS0P;                                                               END;

Страницы