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

                                                                 169                                                                                                                            170

                                                                                                                                                 Нахождение узловых перемещений:

j1 := 1 .. neq                      qtd         := 0              qtd           := 0.5⋅ P               qtd        := P
                                           j1                             14                                  16                         i := 1 .. 6                     ue          := ud
                                                                                                                                                                              i, k            ( MI k , i)
qtd        := P             qtd            := P                       qtd        := P                 qtd        := 0.5⋅ P
      18                            20                                      22                              24
                                                                                                                                         rz :=      i1 ← 0
qtd        := 0.5⋅ P             qtd        := P                 qtd           := P                   qtd        := P                               for k1 ∈ 1 .. nuz
      26                               28                                 30                                32
                                                                                                                                                         for j1 ∈ 1 .. nsu
qtd        := P                  qtd        := 0.5⋅ P                 qtd        := 0.5⋅ P             qtd        := P                                      i1 ← i1 + 1
      34                               36                                   38                               40
                                                                                                                                                            rz          ← ud
                                                                                                                                                              j1 , k1          i1
qtd        := P                    qtd          := 0.75⋅ P            qtd        := 0.5⋅ P             qtd        := 0.25⋅ P                        rz
      42                                   44                               46                               48

                                                                                                                                                 Построение линии прогибов верхней кромки панели:
qtd        := 0.5⋅ P                   qtd         := P                   qtd         := P             qtd        := 0.5⋅ P
      50                                      52                                54                           56
                                                                                                                                         rz         = 32.238                     rz            = 29.789     rz        = 24.556
                                                                                                                                           2 , 29                                    2 , 25                  2 , 19
qtd        := 0.25⋅ P            qtd        := 0.5⋅ P                 qtd        := 0.5⋅ P             qtd        := 0.25⋅ P
      58                               60                                   62                               64                          rz         = 17.509                     rz         = 8.945         rz       =0
                                                                                                                                           2 , 13                                    2, 7                    2, 1

                                                                                                                                         w := rz                                 w := rz                    w := rz
          Учет граничных условий:                                                                                                          1         2 , 29                          2          2 , 25       3        2 , 19
j := 1 .. neq                              i := 1 .. rows ( nsd )
                                                                                                                                         w := rz                                 w := rz                     w := rz
                                                                                                                                           4         2 , 13                          5          2, 7             6        2, 1
qtd := for i ∈ 1 .. rows ( nsd )
      j
                                                                                 qtd
                                                                                        ( nsdi) := sgl ( nsdi , nsdi) ⋅ dis i
                 qtd ← qtd − sgl                              ⋅ dis
                      j                j           j , nsdi           i



sgl              :=       sgl                   if nsd            j             sgl              :=   sgl               if nsd       j
      nsdi , j                  ( nsdi , j)               i                           j , nsdi              j , nsdi             i
                          0 otherwise                                                                 0 otherwise