Исследование напряженно-деформированного состояния однородной цилиндрической оболочки при действии бокового давления, страница 3

Из графика прогиба срединной поверхности (рис.1) можно сказать, что максимальный прогиб будет отрицателен. А максимальное угловое напряжение (рис.4) отрицательно и меняет свой знак.

Модель задачи в Ansys.

Сравнение результатов.

Также задача была смоделирована в конечно разностном пакете Ansys. Результаты мы сравнили с результатами полученными в программе Maple.

Перемещение верхней лицевой поверхности в Ansys (рис.7) и в программе Maple (рас.2), практически полностью совпадают численно и знаками. В нижней части цилиндра перемещение будет максимальное, а в верхней минимальное.

Аналогично можно сказать, что перемещение нижней лицевой поверхности в Ansys (рис.6) совпадает с перемещением в Maple (рис.3). Только теперь в нижней части цилиндра перемещение будет минимальное, а в верхней максимальное.

Модели прогиба срединной поверхности (рис.7) также аналогична результатам полученным в программе Maple (рис.1). Здесь хорошо видно закрепление на концах и отрицательный прогиб. И напряжение в программе Ansys (рис.8), совпадает с результатами полученными в Maple (рис.4).

Вывод.

В работе был попробован метод коечных элементов, но он не дал результатов, так как матрица имеет плохую обусловленность. Поэтому был применён метод Бубнова- Галёркина. Были определены перемещения которые имеют максимальные значения: max=0,00015, min=-0,00015 на верхней лицевой поверхности, max=0,0004, min=-0,0004 на нижней лицевой поверхности, прогиб имеет максимальное значение -0,003, а угловое напряжение- . Также полученные результаты сравнили с результатами полученными в программе Ansys. Полученные результаты практически совпадают, что свидетельствует о возможности применения выбранного метода для решения задач механики деформируемого твёрдого тела.

Приложение.

Maple.

> restart: with(orthopoly): R:=0.45;h:=0.1;L:=2;nu:=0.3;E:=1.60*10^11;sigma:=-200*10^6;

>

>

> de[1]:=diff(N11(s),s);

> de[2]:=diff(Q1(s),s)+(R+h/2)*sigma-N22(s);

> de[3]:=diff(M11(s),s)-Q1(s);

> de[4]:=E/(12(1-nu^2)*R)*(h^2*diff(psi(s),s)+12*diff(u(s),s)+12*nu*w(s))-N11(s);

> de[5]:=E*h^3/(12*(1-nu^2)*R)*(diff(u(s),s)+R*diff(psi(s),s))-M11(s);

> de[6]:=E/(1-nu^2)*(w(s)*(h/R)+h*diff(u(s),s))-N22(s);

> W:=rhs(isolate(de[6]=0,w(s)));

> de[7]:=E*h/(2*(1+nu))*(psi(s)+diff(w(s),s))-Q1(s);

> N11(s):=sum(P(i,s/L)*d[i],i=1..3)+c[5];

> u(s):=sum(P(i,s/L)*A[i],i=1..3)+c[1];

> psi(s):=sum(P(i,s/L)*B[i],i=1..3)+c[2];

> w(s):=sum(P(i,s/L)*C[i],i=1..3)+c[3];

> M11(s):=sum(P(i,s/L)*e[i],i=1..3)+c[4];

> N22(s):=sum(P(i,s/L)*f[i],i=1..3);

> Q1(s):=sum(P(i,s/L)*g[i],i=1..3)+c[7];

> de[1]:=eval(de[1]);

> eq[1]:=int(de[1]*P(1,s/L),s=0..L)=0;

> eq[2]:=int(de[1]*P(2,s/L),s=0..L)=0;

> eq[3]:=int(de[1]*P(3,s/L),s=0..L)=0;

>

> de[2]:=eval(de[2]);

>

> eq[4]:=int(de[2]*P(1,s/L),s=0..L)=0;

> eq[5]:=int(de[2]*P(2,s/L),s=0..L)=0;

> eq[6]:=int(de[2]*P(3,s/L),s=0..L)=0;

>

> de[3]:=eval(de[3]);

> eq[7]:=int(de[3]*P(1,s/L),s=0..L)=0;

> eq[8]:=int(de[3]*P(2,s/L),s=0..L)=0;

> eq[9]:=int(de[3]*P(3,s/L),s=0..L)=0;

>

> de[4]:=eval(de[4]);

> eq[10]:=int(de[4]*P(1,s/L),s=0..L)=0;

> eq[11]:=int(de[4]*P(2,s/L),s=0..L)=0;