Компьютерная обработка информации (в пакете MathCAD). Мартьянова А.Е. - 146 стр.

UptoLike

ВУЗ: 

АГТУ | Астрахань

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

Рубрика: 

Мартьянова А.Е. Компьютерная обработка информации
146
0.5 0 0.5 1 1.5
0
1
узлы КЭ
центры масс КЭ
Пластина из двух конечных элементов
Расчет удвоенной площади конечных элементов -
форм. (7.5в)
∆2n
ie
1
1
1
Xn
ie 1,
Xn
ie 2,
Xn
ie 3,
Yn
ie 1,
Yn
ie 2,
Yn
ie 3,
⎛
⎜
⎜
⎜
⎝
⎞
⎟
⎠
:=
Матрицы связи деформаций и узловых смещений -
форм. (7.12)
Bn
ie j,
1
∆2n
ie
bn
ie j,
0
cn
ie j,
0
cn
ie j,
bn
ie j,
⎛
⎜
⎜
⎜
⎝
⎞
⎟
⎠
⋅:=
BBn ie( ) stack Bn
ie 1,
T
Bn
ie 2,
T
, Bn
ie 3,
T
,
()()
T
:=
Матрица жесткости конечного элемента - форм. (7.13)
K_
ie
BBn ie()
T
D_N⋅ BBn ie()⋅ t⋅
∆
2n
ie
2
⋅:=
Составление глобальной матрицы жесткости (7.14)
ii 1 2 nn_point
⋅..:= jj 1 2 nn_point⋅..:= K
ii jj,
0:=
K
M_
ie ii,
M_
ie jj,
,
()
K
M_
ie ii,
M_
ie jj,
,
()
K_
ie
()
ii jj,
+:=
146                              Мартьянова А.Е. Компьютерная обработка информации

                           Пластина из двух конечных элементов

                   1



                   0

                       0.5             0          0.5            1             1.5
                                     узлы КЭ
                                     центры масс КЭ

      Расчет удвоенной площади конечных элементов -
      форм. (7.5в)

               ⎛ 1 Xnie , 1 Ynie , 1 ⎞
               ⎜
      ∆2nie := ⎜ 1 Xnie , 2 Ynie , 2 ⎟
               ⎜ 1 Xn
               ⎝     ie , 3 Ynie , 3 ⎠

      Матрицы связи деформаций и узловых смещений -
      форм. (7.12)

                        ⎛ bn ie , j   0 ⎞
                    1   ⎜
      Bnie , j :=      ⋅⎜   0       cnie , j ⎟
                  ∆2nie ⎜
                        ⎝ cnie , j bn ie , j ⎠

                       (         (                                   ))
                                                                          T
                                           T        T            T
      BBn( ie) := stack Bnie , 1 , Bnie , 2 , Bnie , 3

      Матрица жесткости конечного элемента - форм. (7.13)

                             T                      ∆2nie
      K_ie := BBn( ie) ⋅ D_N⋅ BBn( ie) ⋅ t⋅
                                                        2

      Составление глобальной матрицы жесткости (7.14)
       ii := 1 .. 2⋅ nn_point           jj := 1 .. 2⋅ nn_point                Kii , jj := 0

       K(M_ie, ii, M_ie, jj) := K(M_ie, ii, M_ie, jj) + ( K_ie)
                                                                ii , jj

Страницы