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

UptoLike

ВУЗ: 

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

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

Рубрика: 

121
defor[iel]:=0.0;
sigma[iel]:=0.0;
for j1:=1 to nue do
begin{2}
if j1=1 then bb[iel,j1]:=-1/dl[iel]
else bb[iel,j1]:=1/dl[iel];
defor[iel]:=defor[iel]+bb[iel,j1]*rz[iel,j1]
end;{2}
sigma[iel]:=defor[iel]*ee[iel]
end;{1}
for j1:=1 to nel do write(j1:3,')',sigma[j1]:5:4);
writeln;
end;
Begin
clrscr;
dl:=kdl;
nug:=knug;
nsd:=knsd;
dis:=kdis; ar:=kar; ee:=kee;
for i:=1 to nuz do
begin
rez[i]:=0.0;
for j:=1 to nuz do sgl[i,j]:=0.0;
end;
for iel:=1 to nel do
begin
MEL;
MGL;
end;
PRAV;
GRAN;
SISTEM;
STRESS;
end.
122
3.5. Расчет рам методом конечных элементов
Матрица жесткости балочного элемента конструкции
Рассмотрим расчетную схему линейно-упругой рамы в глобаль-
ных (общей для всей системы) осях координат xyz (рис.3.9,а).
Рис. 3.9
Разобъем ось этой рамы на m частей (конечных элемен-
тов), соединенных между собой в n узлах (рис.3.9,а). Каждому
элементу под номером е поставим в соответствие систему ло-
кальных осей координат x
/
y
/
z
/
(рис.3.9,б). Рассмотрим в плоско-
сти x
/
y
/
деформацию поперечного изгиба элемента. На концах
этого элемента укажем векторы узловых перемещений
[
]
T
e
uuuuu
4321
)(
=
ρ
и узловых усилий
[
]
T
e
RRRRR
4321
)(
=
ρ
. Нумерация и положительные направ-
ления компонентов этих векторов показаны на рис.3.9,б. Связь 
                                  121                                              122

   defor[iel]:=0.0;                                        3.5. Расчет рам методом конечных элементов
   sigma[iel]:=0.0;
   for j1:=1 to nue do                                Матрица жесткости балочного элемента конструкции
     begin{2}                                         Рассмотрим расчетную схему линейно-упругой рамы в глобаль-
      if j1=1 then bb[iel,j1]:=-1/dl[iel]             ных (общей для всей системы) осях координат xyz (рис.3.9,а).
            else bb[iel,j1]:=1/dl[iel];
      defor[iel]:=defor[iel]+bb[iel,j1]*rz[iel,j1]
     end;{2}
      sigma[iel]:=defor[iel]*ee[iel]
 end;{1}
 for j1:=1 to nel do write(j1:3,')',sigma[j1]:5:4);
 writeln;
end;
      Begin
      clrscr;
      dl:=kdl;
      nug:=knug;
      nsd:=knsd;
      dis:=kdis; ar:=kar; ee:=kee;
      for i:=1 to nuz do
      begin
        rez[i]:=0.0;
        for j:=1 to nuz do sgl[i,j]:=0.0;
      end;                                                                       Рис. 3.9
      for iel:=1 to nel do
        begin                                                  Разобъем ось этой рамы на m частей (конечных элемен-
          MEL;                                        тов), соединенных между собой в n узлах (рис.3.9,а). Каждому
          MGL;                                        элементу под номером е поставим в соответствие систему ло-
        end;                                          кальных осей координат x/y/z/ (рис.3.9,б). Рассмотрим в плоско-
          PRAV;                                       сти x/y/ деформацию поперечного изгиба элемента. На концах
          GRAN;                                       этого элемента укажем векторы узловых перемещений
          SISTEM;                                      ρ
                                                      u ( e ) = [u1 u 2 u 3 u 4 ]
                                                                                 T
                                                                                      и         узловых       усилий
          STRESS;                                      ρ( e )
                                                      R = [R1 R2 R3 R4 ] . Нумерация и положительные направ-
                                                                                   T
      end.
                                                      ления компонентов этих векторов показаны на рис.3.9,б. Связь

Страницы