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

UptoLike

ВУЗ: 

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

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

Рубрика: 

165
Матрица координат узлов и глобальные номера узлов ансамбля
элементов:
nx
6
:=
ny
6
:=
nx1
4
:=
ny1
4
:=
cuz k1 0←
sx hx
1
−←
sy hy
1
−←
sx sx hx
i1
+←
k1 k1 1+←
sy sy hy
j1
+←
cuz
k1 1,
sx←
cuz
k1 2,
sy←
j1 1 ny..∈for
i1 nx1≤if
sy hy
1
−←
sx sx hx
i1
+←
k1 k1 1+←
sy sy hy
j1
+←
cuz
k1 1,
sx←
cuz
k1 2,
sy←
j1 1 ny1..∈for
i1 nx1>if
i1 1 nx..∈for
cuz
:=
cuz
12
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
00
0 0.1
0 0.2
0 0.3
0 0.4
0 0.5
0.1 0
0.1 0.1
0.1 0.2
0.1 0.3
0.1 0.4
0.1 0.5
0.2 0
0.2 0.1
0.2 0.2
0.2 0.3
0.2 0.4
0.2 0.5
0.3 0
0.3 0.1
0.3 0.2
0.3 0.3
0.3 0.4
=
166
Генерация глобальных номеров узлов элементов ансамбля:
nug
37 1,
2
5
:=
nug
37 2,
2
9
:=
nug
37 3,
2
6
:=
nug
38 1,
3
0
:=
nug
38 2,
2
6
:=
nug
38 3,
2
9
:=
nug
39 1,
2
6
:=
nug
39 2,
3
0
:=
nug
39 3,
2
7
:=
                                              165                                                                  166

Матрица координат узлов и глобальные номера узлов ансамбля                      Генерация глобальных номеров узлов элементов ансамбля:
элементов:
 nx := 6           ny := 6              nx1 := 4           ny1 := 4


 cuz :=    k1 ← 0
                                                                1       2
           sx ← −hx
                      1                                    1        0       0
           for i1 ∈ 1 .. nx                                2        0   0.1
                 if i1 ≤ nx1                               3        0   0.2
                    sy ← −hy                               4        0   0.3
                                  1
                                                           5        0   0.4
                    sx ← sx + hx
                                       i1                  6        0   0.5
                     for j1 ∈ 1 .. ny                      7    0.1         0
                          k1 ← k1 + 1                      8    0.1     0.1

                          sy ← sy + hy                     9    0.1     0.2
                                              j1
                                                           10   0.1     0.3
                          cuz          ← sx
                              k1 , 1                       11   0.1     0.4
                                                   cuz =
                          cuz          ← sy                12   0.1     0.5
                              k1 , 2
                                                           13   0.2         0
                 if i1 > nx1
                                                           14   0.2     0.1
                    sy ← −hy                               15   0.2     0.2
                                  1
                    sx ← sx + hx                           16   0.2     0.3
                                       i1
                                                           17   0.2     0.4
                     for j1 ∈ 1 .. ny1
                                                           18   0.2     0.5
                          k1 ← k1 + 1                      19   0.3         0    nug            := 25   nug            := 29   nug            := 26
                                                                                       37 , 1                 37 , 2                 37 , 3
                          sy ← sy + hy                     20   0.3     0.1
                                              j1
                          cuz          ← sx                21   0.3     0.2
                              k1 , 1                                             nug            := 30   nug            := 26   nug            := 29
                                                           22   0.3     0.3            38 , 1                 38 , 2                 38 , 3
                          cuz          ← sy
                              k1 , 2                       23   0.3     0.4
                                                                                 nug            := 26   nug            := 30   nug            := 27
           cuz                                                                         39 , 1                 39 , 2                 39 , 3

Страницы