Xodjayorova dilxayot abdiaxatovna deformatsiyalanuvchi muhitdagi silindrik qobiqning buralma tebranishlarini sonli tadqiq etish

§3.2 Doiraviy silindrik qobiqning buralma tebranishlarida
ko’ndalang siljish deformatsiyasi ta’sirini hisobga olish
Doiraviy silindrik qobiq uchun aylanish inersiyasi va ko’ndalang siljish deformatsiyasi ta’sirini hisobga olganda (3.6) va (3.7) tenglamalar sistemasi uchun boshlang’ich va chegaraviy shartlardan foiydalanib ,,Paskal” dasturlash tilida tuzilgan dastur yordamida ko’chishlar va kuchlanishlarning to’rning tugun nuqtalaridagi sonli qiymatlarini aniqlaymiz.
Program Xodjayarova_D;
Uses crt;
const Nz=99; Nt=100; jt=100; kt=2;
type
ff = text;
vec1 = array [0..Nz+1] of real;
var i,j,k,n,it : integer;
u1,u2,u3,v1,v2,v3,x,y,ck,dk,fk,mk,nk,szt : vec1;
l,d,ta,d2,ta2,t,t0,t1,
al1,al0,bet0,bet1,Ab,Bb,r0,ak,gam,eta,gam1,eta1,omega,omega2,
mu,mu1, r1,r2,co,e,nu,aco,bco,cco,kco,dco,qo,
a1,a2,s1,s2,s3,s4,k1,k2,k3,k4 : real;
w1,w2,w3,sigma,tn: array[1..100] of real;
fl,fli : text;
procedure prog3(var y : vec1);
begin
ck[0]:=al1/(al0*d-al1);
dk[0]:=d*Ab/al1;
for i:=1 to Nz do
begin
ck[i]:=1/(mk[i]-nk[i]*ck[i-1]);
dk[i]:=fk[i]-nk[i]*ck[i-1]*dk[i-1];
end;
y[Nz+1]:=(Bb*d+bet1*ck[Nz]*dk[Nz])/(bet0*d+bet1*(ck[Nz]+1));
i:=Nz+1;
while i>=1 do begin i:=i-1; y[i]:=ck[i]*(dk[i]-y[i+1]); end;
end;
function f(t,t1:real): real;
var pi: real;
begin
pi:=3.14159;
if t<=t1 then f:=sin(pi*t/t1) else f:=0;
end;
begin
clrscr;
assign(fl,'BerdZ.dat'); rewrite(fl);
assign(fli,'BerdT.dat'); rewrite(fli);
{Boshlang'ich ma'lumotlar}
l:=10; r0:=1; gam:=0.5; eta:=1-gam; omega:=0.2; omega2:=omega*omega;
mu:=1; mu1:=1.5;
{Boshlang'ich hisoblar}
d:=l/(Nz+1); ta:=0.8*d; d2:=d*d; ta2:=ta*ta; t1:=Nt*ta/4;
{Boshlang'ich shartlar}
for i:=0 to Nz+1 do
begin x[i]:=i*d;
u1[i]:=0; u2[i]:=0; u3[i]:=0;
v1[i]:=0; v2[i]:=0; v3[i]:=0 end;
e:=2.0e11;
nu:=0.3;
mu:=e/(2*(1+nu));
r1:=1;
r2:=1.1;
a1:=1;
a2:=1;
co:=2*(1+nu);
S1:=2*co*(sqr(sqr(r2))-sqr(sqr(r1)))/(sqr(r1*r2)*ln(r2/r1));
S4:=16/(sqr(sqr(r1))+sqr(sqr(r2)));
S3:=(sqr(r2)-sqr(r1))/(sqr(r1*r2)*ln(r2/r1));
S2:=32*(sqr(sqr(r2))-sqr(sqr(r1)))/(sqr(sqr(r1*r2))*ln(r2/r1));
aco:=-e/mu/(a1*co);
bco:=(a2*co+a1*e/mu)/(a1*co);
cco:=S1*S4/(a1*co)/100;
dco:=-S3*e/mu/(a1*co);
kco:=-S2*e/mu/(a1*co);
k1:=bco/d2;
k2:=aco/sqr(d2);
k3:=dco/d2;
k4:=cco;
j:=0;

t:=0; j:=0; k:=0;n:=0;

for it:=1 to Nt do
begin
t:=t+ta;
{Tenglamani yechish}
for i:=0 to Nz+1 do
begin
{Hamma hadlarda ikkiga ajratish}

ak:=-gam/d2;

mk[i]:=(1/ta2+2*gam/d2-cco)/ak;

nk[i]:=1;

qo:=dco/d2*(v2[i+1]-2*v2[i]+v2[i-1])+kco*v2[i];
fk[i]:=(1/ta2*(2*u2[i]-u1[i])+eta/d2*(u1[i+1]-2*u1[i]+u1[i-1])+
cco*u2[i-1]+qo)/ak;

end;
al0:=1; al1:=1; bet0:=1; bet1:=1;

Ab:=f(t,t1); Bb:=0;
prog3(u3);
{Chegaraviy shartlar}
u3[0]:=f(t,t1); u3[Nz+1]:=0;
for i:=0 to Nz+1 do
begin szt[i]:=(u3[i+1]-u3[i-1])/(2*d);
v3[i]:=2*v2[i]-v1[i]+ta2*u3[i];
end;
for i:=0 to Nz+1 do begin
u1[i]:=u2[i]; u2[i]:=u3[i];
v1[i]:=v2[i]; v2[i]:=v3[i]; end;
{Pechat}
j:=j+1; k:=k+1;
if k=kt then begin k:=0; n:=n+1;tn[n]:=t;
w1[n]:=u3[25]; w2[n]:=u3[50]; w3[n]:=u3[75]; end;
if j=jt then
begin
j:=0;
for i:=0 to Nz+1 do
begin
writeln(x[i]:10:5,' ',u3[i]:10:5);
writeln(fl,x[i]:10:5,' ',u3[i]:10:5);
end;
end;
end;
close(fl);
for i:=1 to n do writeln(fli,tn[i]:10:4,w1[i]:10:5,w2[i]:10:5,w3[i]:10:5);
close(fli);
end.

Olingan natijalarni grafik ko’rinishida tasvirlaymiz.

Klassik i=25
Klassik i=50

Ko’ndalang siljish deformatsiyasi i=25
Ko’ndalang siljish deformatsiyasi i=50