ВУЗ:
Составители:
Рубрика:
104
46 pr2= abs((fi(N,1)-ffi(N))/fi(N,1))
if pr2<eps goto 50
ffi(N)=fi(N,1)
s1(i)=s1+1
goto 45
50 dfi(N,1)=df(N,1)-b*(f(N,1)+fi(N,1))
P(N)=-df(N,1)+dfi(N,1)
u(N)=f(N,1)+fi(N,1)
if i>N goto 55
Qs=Qs+m*(df(N,1)+dfi(N,1))/va
goto 57
55 if (-df(N,1)+dfi(N,1))<0 goto 56
Qs=Qs
goto 57
56 Qs=Qs+m*(df(N,1)+dfi(N,1))/va:Qn=Qn+m*(-df(N,1)+dfi(N,1))/va
57 f(N,1)=f(N,0)+df(N,1)*dl:fi(N,1)=fi(N,0)+dfi(N,1)*dl
v(0)=a*(df(0,1)+dfi(0,1))
if i>20*N goto 140
P(0)=-df(0,1)+dfi(0,1)
Qm=V/V0
BQ=Qm+Qs-Qn
if i/z<40 goto 100
print "Относительное время t=";t;" ";
print "Деформация в ударном сечении Po=";P(0);" ";
print "Скор ость ударного сечения V="; V
print " Перемещение ударного сечения uo=";u(0);" ";
print " Перемещение ударной массы xm=";xm;" ";
print " Перемещение нижнего конца сваи un=";u(N)
print "Количество движения ударной массы Qm=";Qm;" ";
print " Количество движения стержня Qs=";Qs;" ";
print "Импульс силы Qn=";Qn;" ";
print " Баланс количества движения BQ=";BQ
print
z=z+1
100 for j=0 to N
f(j,0)=f(j,1)
df(j,0)=df(j,1) fi(j,0)=fi(j,1)
dfi(j,0)=dfi(j,1)
Qs=0
next j; next i
140 print "Относительное время t=";t;" ";
print " Число элементов разбиения N=";N;
print "Жесткость
верхней прокладки наголовника b=";b;
print "Жесткость нижней прокладки наголовника b1=";b1;
print " Скор ость звука в материале a=";a
ee:
END
