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

UptoLike

ВУЗ: 

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

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

Рубрика: 

167
Формирование матрицы индексов степеней свободы:
Формирование матрицы жесткости элемента k
nug
40 1,
3
1
:=
nug
40 2,
2
7
:=
nug
40 3,
3
0
:=
nug
40 1,
3
1
:=
nug
40 2,
2
7
:=
nug
40 3,
3
0
:=
nug
42 1,
3
2
:=
nug
42 2,
2
8
:=
nug
42 3,
3
1
:=
i435
0
..:=
j1
3
..:=
nug
ij,
0
:=
k1ne
l
..:=
i1
3
.
.:=
k1ne
l
..:=
i1
3
.
.:=
MI
k2i⋅ 1−,
2 nug
ki,
⋅ 1
−
:=
MI
k2i⋅,
2 nug
k
i
,
⋅
:=
x cuz
1
〈〉
:=
y cuz
2
〉
:=
ak()
x
nug
k3,
()
x
nug
k2,
()
−
x
nug
k1,
()
x
nug
k3,
()
−
x
nug
k2,
()
x
nug
k1,
()
−
:=
bk()
y
nug
k2,
()
y
nug
k3,
()
−
y
nug
k3,
()
y
nug
k1,
()
−
y
nug
k1,
()
y
nug
k2,
()
−
:=
A
k
0.5
1
1
1
x
nug
k1,
()
x
nug
k2,
()
x
nug
k3,
()
y
nug
k1,
()
y
nug
k2,
()
y
nug
k3,
()
:=
A
1
510
3−
×=
168
Матрица деформаций элемента под номером k:
Матрица упругости элемента:
Матрица жесткости элемента k:
Формирование глобальной матрицы жесткости и правой
части системы уравнений
i 1 neq
.
.:=
j 1 neq
.
.:=
sgl
ij,
0
:=
i1
6
.
.:=
j1
6
.
.:=
k1ne
l
.
.:=
                                                        167                                                                                             168


nug            := 31                         nug            := 27                   nug                 := 30
      40 , 1                                       40 , 2                                      40 , 3                   Матрица деформаций элемента под номером k:

nug            := 31                         nug            := 27                   nug                 := 30
      40 , 1                                       40 , 2                                      40 , 3


nug            := 32                         nug            := 28                   nug                 := 31
      42 , 1                                       42 , 2                                      42 , 3


i := 43.. 50                                 j := 1 .. 3                            nug               := 0
                                                                                               i, j                     Матрица упругости элемента:
        Формирование матрицы индексов степеней свободы:
k := 1 .. nel                  i := 1 .. 3


k := 1 .. nel                  i := 1 .. 3


MI                 := 2⋅ nug          −1                    MI              := 2⋅ nug
     k , 2⋅ i− 1               k, i                              k , 2⋅ i               k, i
        Формирование матрицы жесткости элемента k
        〈1〉                        〈2〉
x := cuz                   y := cuz


         x( nugk , 3) − x( nugk , 2)                                        y ( nugk , 2) − y( nugk , 3)            Матрица жесткости элемента k:
                                                                                                          
         x
a( k) := ( nugk , 1)   − x                                                   y             − y             
                         ( nugk , 3)                              b ( k) := ( nugk , 3)
                                                                                                ( nugk , 1) 
         x nug − x nug                                                      y nug − y nug 
         ( k , 2) ( k , 1)                                                  ( k , 1) ( k , 2)                     Формирование глобальной матрицы жесткости и правой
                                                                                                                 части системы уравнений
         1            x
                        ( nugk , 1) y( nugk , 1) 
                                                                                                                i := 1 .. neq       j := 1 .. neq             sgl          := 0
                                                                                                                                                                     i, j
A := 0.5 1                                       
                       x               y                                                  −3
 k                     ( nugk , 2) ( nugk , 2)                      A = 5 × 10
                                                                            1
                                                                                                                 i := 1 .. 6     j := 1 .. 6         k := 1 .. nel
         1
         
                       x
                        ( nugk , 3) y( nugk , 3) 

Страницы