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

UptoLike

ВУЗ: 

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

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

Рубрика: 

83
(*МАТРИЦА ПОДАТЛИВОСТЕЙ DФ*)
KD:MASS1=
(( 5.64, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),
( 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),
( 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),
( 0, 0, 0,5.64, 0, 0, 0, 0, 0, 0, 0, 0, 0),
( 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0),
( 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0),
( 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0),
( 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0),
( 0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0),
( 0, 0, 0, 0, 0, 0, 0, 0, 0,5.64, 0, 0, 0),
( 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0),
( 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,5.64, 0),
( 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4));
(*ВЕКТОР ПРОДОЛЬНЫХ СИЛ ОТ ВНЕШНЕЙ НАГРУЗКИ В
ОСНОВНОЙ ФЕРМЕ*)
KNP:MASS4=( -17.7,-15,-15,-17.7,12.5,12.5,12.5,12.5,5,3.53,-
5,3.53,5);
BEGIN
(*ВВОД ВЕКТОРА ПРОДОЛЬНЫХ СИЛ ОТ ЕДИНИЧНОЙ
НАГРУЗКИ*)
N0:=KN0;
WRITELN('ВЕКТОР ПРОДОЛЬНЫХ СИЛ ОТ ЕД.НАГРУЗКИ
В ОСН.ФЕРМЕ');
FOR I:=1 TO K DO BEGIN
FOR J:=1 TO N DO WRITE(' ', N0[I,J]:6:2); WRITELN
END;
WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ <ENTER>');
READLN;READLN;
(*ВВОД МАТРИЦЫ ПОДАТЛИВОСТЕЙ*)
D:=KD;
WRITELN('МАТРИЦА ПОДАТЛИВОСТЕЙ');
FOR I:=1 TO K DO BEGIN
FOR J:=1 TO K DO WRITE(' ', D[I,J]:6:2); WRITELN
END;
84
WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ <ENTER>');
READLN;READLN;
(*ВВОД ВЕКТОРА ГРУЗОВЫХ ПРОДОЛЬНЫХ СИЛ*)
NP:=KNP;
WRITELN('ВЕКТОР ПРОДОЛЬНЫХ СИЛ ОТ ВНЕШНЕЙ НА-
ГРУЗКИ В ОСН.ФЕРМЕ');
FOR I:=1 TO K DO WRITE(' ', NP[I]:6:2); WRITELN;
WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ <ENTER>');
READLN;READLN;
FOR I:=1 TO K DO FOR J:=1 TO N DO NT[J,I]:=N0[I,J];
(*ПЕРЕМНОЖЕНИЕ МАТРИЦ*)
FOR I:=1 TO N DO
FOR J:=1 TO K DO
FOR L:=1 TO K DO
H[I,J]:=H[I,J]+NT[I,L]*D[L,J];
FOR I:=1 TO N DO
BEGIN
FOR J:=1 TO N DO
FOR L:=1 TO K DO
A[I,J]:=A[I,J]+H[I,L]*N0[L,J];
FOR L:=1 TO K DO
C[I]:=C[I]-H[I,L]*NP[L]
END;
(*РЕШЕНИЕ СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ
УРАВНЕНИЙ*)
FOR I:=1 TO N-1 DO
FOR J:=I+1 TO N DO
BEGIN
A[J,I]:=-A[J,I]/A[I,I];
FOR KK:=I+1 TO N DO
A[J,KK]:=A[J,KK]+A[J,I]*A[I,KK];
C[J]:=C[J]+A[J,I]*C[I]
END;
X[N]:=C[N]/A[N,N];
FOR I:=N-1 DOWNTO 1 DO
BEGIN 
                             83                                                         84

(*МАТРИЦА ПОДАТЛИВОСТЕЙ DФ*)                                    WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ ');
KD:MASS1=                                                      READLN;READLN;
        (( 5.64, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),          (*ВВОД ВЕКТОРА ГРУЗОВЫХ ПРОДОЛЬНЫХ СИЛ*)
         ( 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),              NP:=KNP;
         ( 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),              WRITELN('ВЕКТОР ПРОДОЛЬНЫХ СИЛ ОТ ВНЕШНЕЙ НА-
         ( 0, 0, 0,5.64, 0, 0, 0, 0, 0, 0, 0, 0, 0),           ГРУЗКИ В ОСН.ФЕРМЕ');
         ( 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0, 0),              FOR I:=1 TO K DO WRITE(' ', NP[I]:6:2); WRITELN;
         ( 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0, 0),              WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ ');
         ( 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0, 0),             READLN;READLN;
         ( 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0, 0),              FOR I:=1 TO K DO FOR J:=1 TO N DO NT[J,I]:=N0[I,J];
         ( 0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0, 0, 0),             (*ПЕРЕМНОЖЕНИЕ МАТРИЦ*)
         ( 0, 0, 0, 0, 0, 0, 0, 0, 0,5.64, 0, 0, 0),            FOR I:=1 TO N DO
         ( 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4, 0, 0),               FOR J:=1 TO K DO
         ( 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,5.64, 0),                   FOR L:=1 TO K DO
         ( 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 4));                          H[I,J]:=H[I,J]+NT[I,L]*D[L,J];
(*ВЕКТОР ПРОДОЛЬНЫХ СИЛ ОТ ВНЕШНЕЙ НАГРУЗКИ В                   FOR I:=1 TO N DO
ОСНОВНОЙ ФЕРМЕ*)                                                 BEGIN
KNP:MASS4=( -17.7,-15,-15,-17.7,12.5,12.5,12.5,12.5,5,3.53,-           FOR J:=1 TO N DO
5,3.53,5);                                                                   FOR L:=1 TO K DO
BEGIN                                                                         A[I,J]:=A[I,J]+H[I,L]*N0[L,J];
(*ВВОД ВЕКТОРА ПРОДОЛЬНЫХ СИЛ ОТ ЕДИНИЧНОЙ                             FOR L:=1 TO K DO
НАГРУЗКИ*)                                                                   C[I]:=C[I]-H[I,L]*NP[L]
 N0:=KN0;                                                        END;
 WRITELN('ВЕКТОР ПРОДОЛЬНЫХ СИЛ ОТ ЕД.НАГРУЗКИ                 (*РЕШЕНИЕ СИСТЕМЫ ЛИНЕЙНЫХ АЛГЕБРАИЧЕСКИХ
В ОСН.ФЕРМЕ');                                                 УРАВНЕНИЙ*)
 FOR I:=1 TO K DO BEGIN                                         FOR I:=1 TO N-1 DO
        FOR J:=1 TO N DO WRITE(' ', N0[I,J]:6:2); WRITELN        FOR J:=I+1 TO N DO
END;                                                                   BEGIN
 WRITELN('ДЛЯ ПРОДОЛЖЕНИЯ НАЖМИТЕ ');                                 A[J,I]:=-A[J,I]/A[I,I];
READLN;READLN;                                                               FOR KK:=I+1 TO N DO
(*ВВОД МАТРИЦЫ ПОДАТЛИВОСТЕЙ*)                                                 A[J,KK]:=A[J,KK]+A[J,I]*A[I,KK];
 D:=KD;                                                                      C[J]:=C[J]+A[J,I]*C[I]
 WRITELN('МАТРИЦА ПОДАТЛИВОСТЕЙ');                                     END;
 FOR I:=1 TO K DO BEGIN                                         X[N]:=C[N]/A[N,N];
        FOR J:=1 TO K DO WRITE(' ', D[I,J]:6:2); WRITELN        FOR I:=N-1 DOWNTO 1 DO
END;                                                                  BEGIN

Страницы