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

UptoLike

ВУЗ: 

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

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

Рубрика: 

43
Рис.1.8
да
B1
j
:=1
,
nel
i:=1
,
nuz
SC[i,j]=1
нет
SZ[2*i-1,j]:=
α
[j,1];
SZ[2*i-1,j]:=α[j,2];
SC[i,j]=-1
SZ[2*i-1,j]:=-
α
[j,1];
SZ[2*i-1,j]:=-α[j,2];
да
нет
Вво
д
Р
Получение матрицы
P
S и
векто
р
а
Q
ρ
Решение СЛАУ методом Гаус-
са с выделением главного
элемента
Печать вектора
N
ρ
конец
В
i:=1,nuz2
j
:=1,nel
SZ[i,j]:=0.0
В1
Составление матрицы
*
44
1.2.3. Программа для расчета ферм на алгоритмиче-
ском языке Турбо Паскаль
Program ferma;
uses Crt;
label 1;
const nel=17; {число стержней фермы}
nuz=10; {число узлов фермы}
nuz2=20; {удвоенное число узлов}
nopr=3; {число уравнений, которые нужно удалить}
type mas1=array[1..nuz,1..nel] of integer;
mas2=array[1..nel,1..nuz] of integer;
mas3=array[1..nuz2,1..nel] of real;
mas4=array[1..nel,1..nel] of real;
mas5=array[1..nel,1..2] of real;
mas6=array[1..nuz,1..2] of real;
mas7=array[1..nel] of real;
mas8=array[1..nuz2] of real;
mas9=array[1..nopr] of integer;
var i,j,k,i1:integer;
(* sc-структурная матрица; st-трансп.структ. матрица;
Dl-вектор длин стержней; pr-вектор проекций;
ALFA-вектор направляющих косинусов;
c-вектор координат узлов; p-вектор внешней нагрузки;
sz-прямоуг.матрица s*; n-вектор номеров строк,
которые нужно удалить из sz и p; sp-матрица СЛАУ;
Q-вектор правой части СЛАУ; *)
sc:mas1; st:mas2; dl,Q,b,x:mas7; pr,ALFA:mas5;
sz:mas3; c:mas6; sp,a,a1:mas4; p:mas8;
n:mas9;
const
(* задание структурной матрицы *)
ksc:mas1=(( 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),
(-1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),
( 0,-1,-1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),
( 0, 0, 0,-1,-1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0), 
                             43                                                                          44


           В                                  B1                              1.2.3. Программа для расчета ферм на алгоритмиче-
                                                                      ском языке Турбо Паскаль
                                           j:=1,nel                   Program ferma;
       i:=1,nuz2
                                                                       uses Crt;
                                           i:=1,nuz                    label 1;
       j:=1,nel
                                                                         const nel=17; {число стержней фермы}
                                                                нет              nuz=10; {число узлов фермы}
      SZ[i,j]:=0.0                        SC[i,j]=1
                                                                                 nuz2=20; {удвоенное число узлов}
                                         да                                      nopr=3; {число уравнений, которые нужно удалить}
                                    SZ[2*i-1,j]:=α[j,1];               type mas1=array[1..nuz,1..nel] of integer;
                                    SZ[2*i-1,j]:=α[j,2];                         mas2=array[1..nel,1..nuz] of integer;
          В1                                                                     mas3=array[1..nuz2,1..nel] of real;
                                                                                 mas4=array[1..nel,1..nel] of real;
Составление матрицы                                                              mas5=array[1..nel,1..2] of real;
                              нет
         S*                               SC[i,j]=-1                             mas6=array[1..nuz,1..2] of real;
                                                                                 mas7=array[1..nel] of real;
                                         да
                                                                                 mas8=array[1..nuz2] of real;
                                    SZ[2*i-1,j]:=-α[j,1];                        mas9=array[1..nopr] of integer;
                                    SZ[2*i-1,j]:=-α[j,2];
                                                                       var i,j,k,i1:integer;
                                                                                 (* sc-структурная матрица; st-трансп.структ. матрица;
                                                                                  Dl-вектор длин стержней; pr-вектор проекций;
                                                                                  ALFA-вектор направляющих косинусов;
                                           Ввод Р                                 c-вектор координат узлов; p-вектор внешней нагрузки;
                                                                                  sz-прямоуг.матрица s*; n-вектор номеров строк,
                              Получение матрицы        SP   и                     которые нужно удалить из sz и p; sp-матрица СЛАУ;
                                      ρ                                           Q-вектор правой части СЛАУ; *)
                              вектора Q                                           sc:mas1; st:mas2; dl,Q,b,x:mas7; pr,ALFA:mas5;
                                                                                  sz:mas3; c:mas6; sp,a,a1:mas4; p:mas8;
                              Решение СЛАУ методом Гаус-
                              са с выделением главного
                                                                                  n:mas9;
                              элемента                                const
                                                                      (* задание структурной матрицы *)
                                                       ρ                     ksc:mas1=(( 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),
                     конец          Печать вектора     N                                  (-1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),
                                                                                          ( 0,-1,-1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),
                                                                                          ( 0, 0, 0,-1,-1, 0, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0),
                             Рис.1.8

Страницы