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

UptoLike

ВУЗ: 

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

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

Рубрика: 

61
A = M
ed
T
·B·M
ed
– матрица реакций R
r
;
C = M
ed
T
·B·M
/
P
– вектор реактивных усилий в доп. связях
R
P
.
I = 1
,
N
J = 1
,k
L = 1
,k
H
[
I
,
J
]
=H
[
I
,
J
]
+MT
[
I
,
L
]
*B
[
L
,
J
]
I = 1
,
N
J = 1
,k
L = 1
,k
A
[
I
,
J
]
=A
[
I
,
J
]
+H
[
I
,
L
]
*M0
[
L
,
J
]
2 3
1
Выч. H=M
ed
T
·B
62
L
=
1,k
=
+
*
3
I
=
1,N
-
1
J
=
I
+
1,N
A[J,I]=-A[J,I]/A[I,I]
Kk
=
I
+
1,N
A[J,kk]
=
A[J,kk]
+
A[J,I]
*
A[I,kk]
C[J]=C[J]+A[J,I]*C[I]
4
X[N]:
=
C[N]/A[N,N]
Обратный ход
метода Гаусса
Реш. СЛАУ
                             61                                                          62


                                     1         Выч. H=MedT·B
                                                                                                                     3
                                  I = 1,N                                              L=1,k

                                                                               C[I]=C[I]+H[I,L]*MP[L]
                                   J = 1,k

                                  L = 1,k
                                                                   Реш. СЛАУ
                                                                                       I=1,N-1
                        H[I,J]=H[I,J]+MT[I,L]*B[L,J]
                                                                                       J=I+1,N

                                                                                  A[J,I]=-A[J,I]/A[I,I]

                                  I = 1,N
                                                                                      Kk=I+1,N
                                   J = 1,k                                 A[J,kk]=A[J,kk]+A[J,I]*A[I,kk]

                                  L = 1,k                                       C[J]=C[J]+A[J,I]*C[I]

                        A[I,J]=A[I,J]+H[I,L]*M0[L,J]



                                                                                   X[N]:=C[N]/A[N,N]
                                     2                         3
                                                                                                          Обратный ход
      A = MedT·B·Med – матрица реакций Rr;                                                    4           метода Гаусса
      C = MedT·B·M/P – вектор реактивных усилий в доп. связях
RP.

Страницы