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

UptoLike

ВУЗ: 

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

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

Рубрика: 

151
Составление матрицы индексов
:
IT
T
1
2
3
4
4
5
2
4
4
6
=
IT
1
2
3
4
4
5
2
4
4
6
T
:=
nue
2
:=
nse
6
:=
neq 1
8
:=
nel
5
:=
I 1 neq..:=
J 1 neq
.
.:=
SGL
IJ,
0
:=
QGL
I
0
:=
i1 1 ne
l
..:=
j1 1
3
..:=
k1 1
3
..:=
IG
i1 k1,
IT
i1 1,
1
(
)
3
⋅
k
1
:=
JG
i1 k1,
IT
i1 2,
1
(
)
3
⋅
k
1
:=
e1ne
l
..:=
i 1 nse
.
.:=
j 1 nse..:=
seg e i, j,( ) seg e()
i
j
,
:=
Q1 e i,() Qe()
i
:=
SGL
MI
ei,
MI
ej,
,
()
SGL
MI
ei,
MI
ej,
,
()
seg e i
,
j
,
()+:=
QGL
MI
ei,
()
QGL
MI
ei,
()
Q1 e i,()
:=
QGL
T
12345
1
0 -2.5 -2.083 0 -2.5
=
152
Учет граничных условий
Обновление матрицы системы в соответствии с заданными
граничными условиями:
i1 1 ndis
.
.:=
nsd
1
1:=
nsd
2
2
:=
nsd
3
4
:=
nsd
4
5
:=
nsd
5
8
:=
nsd
6
1
1
:=
nsd
7
1
4
:=
nsd
8
1
6
:=
nsd
9
1
7
:=
nsd
10
1
8
:=
QGL
j1 nsd
i1
←
k1 nsd
i1
←
QGL
k1
QGL
k1
SGL
k1 j1,
dis
i1
⋅−←
k1 k1 1+←
k1 neq≤if
l1 1 neq..∈for
QGL
j1
dis
i1
←
i1 1 ndis
.
.
∈
for
QGL
:=
QGL
T
12345
1
0 0 -2.083 0 0
=
neq 1
8
:=
k1 1 neq
.
.:=
SGL
nsd
i1
k1,
if k1 nsd
i1
1
,
0
,
(
)
:=
dis
i1
0
:=
                                                       151                                                                                                           152

                                                                                                                     Учет граничных условий


                                                                                                             nsd := 1 nsd := 2                     nsd := 4               nsd := 5          nsd := 8
                                                                                                                 1               2                       3                   4                        5
                                                                                                             nsd := 11 nsd := 14 nsd := 16 nsd := 17                                        nsd            := 18
                                                                                                                 6                   7                       8               9                        10
                                                                                                             i1 := 1 .. ndis             dis := 0
                                                                                                                                              i1

       Составление матрицы индексов:                                                                          QGL:=          for i1 ∈ 1 .. ndis
                  T                                                                                                              j1 ← nsd
                                                                                                                                                   i1
      1 3 4 2 4                                                          T 1 3 4 2 4
                                                                        IT = 
IT := 
      2 4 5 4 6                                                            2 4 5 4 6                                         k1 ← nsd
                                                                                                                                                   i1
                                                                                                                                 for l1 ∈ 1 .. neq
nue := 2             nse := 6                neq := 18                  nel := 5
                                                                                                                                     if k1 ≤ neq
I := 1 .. neq                  J := 1 .. neq                SGL         := 0             QGL := 0                                          QGL           ← QGL            − SGL          ⋅ dis
                                                                 I, J                           I                                                  k1                k1          k1 , j1         i1
                                                                                                                                           k1 ← k1 + 1
      i1 := 1 .. nel                   j1 := 1 .. 3         k1 := 1 .. 3
                                                                                                                                 QGL ← dis
                                                                                                                                          j1             i1
IG
   i1 , k1      (
             := IT
                    i1 , 1         )
                             − 1 ⋅ 3 + k1                      JG
                                                                 i1 , k1
                                                                           := IT(   i1 , 2      )
                                                                                             − 1 ⋅ 3 + k1                    QGL


                                                                                                                  T                  1              2                3           4               5
                                                                                                             QGL =
e := 1 .. nel                i := 1 .. nse                  j := 1 .. nse                                                1                0                  0   -2.083              0                0

seg ( e , i, j) := seg ( e)                           Q1( e , i) := Q( e)                                        Обновление матрицы системы в соответствии с заданными
                               i, j                                         i
                                                                                                            граничными условиями:
SGL
     ( MI e, i , MI e, j) := SGL( MI e, i , MI e, j) + seg ( e , i, j)                                       neq := 18
                                                                                                             k1 := 1 .. neq
QGL
      ( MI e, i) := QGL( MI e , i) + Q1( e , i)
                        1                2              3           4               5
                                                                                                             SGL
                                                                                                                 nsd i1 , k1         (
                                                                                                                               := if k1        nsd , 1 , 0
                                                                                                                                                        i1       )
     T
QGL =
              1                0         -2.5         -2.083               0         -2.5

Страницы