ВУЗ:
Составители:
Рубрика:
103
Q(j,2)=2/(L(j)+1)
next j
d2f(0,0)=va*b
P(0)=-df(0,0)+dfi(0,0)
u(0)=f(0,0)+fi(0,0)
xm=u(0)-(-df(0,0)+dfi(0,0))/b
u(N)=f(N,0)+fi(N,0)
v=a*(df(0,0)+dfi(0,0))-a*(-d2f(0,0)+d2fi(0,0))/b
print "Относительное время t=";t;" ";
print "Деформация в ударном сечении Po=";P(0);" ";
print "Скор ость ударного сечения V=";V
print " Перемещение ударного сечения uo=";u(0);" ";
print " Перемещение ударной массы xm=";xm;" ";
print " Перемещение нижнего конца сваи un=";u(N)
40 for i=1 to 50*N
t=i*dt
s=1
s1=1
if b/(N*m)>4 goto 42
Pt=0
tm=0
goto 43
42 q=(1-4*N*m/b)^.5:q1=q*b*i/N:q2=(q-1)*b*i/(2*N)
Pt=exp(q2)*(1-exp(-q1))/q
tm=(-N*log((1-q)/(1+q)))/(q*b)
if (xm-u(0))>0 goto 43
v=v:xm=xm+v*dl/a
dfi(0,1)=dfi(1,0)
df(0,1)=dfi(0,1) vm=vm
v1=v1
P(0)=0
u(0)=u(0)+(df(0,1)+dfi(0,1))*dl
goto 44
43 xm=xm+v*dl/a:dfi(0,1)=dfi(1,0):df(0,1)=b*(xm-u(0))+dfi(0,1) u(0)=u(0)+(df(0,1)+dfi(0,1))*dl
Vm=Vm-b*a*m*(xm-u(0))
v1=v1+a*m*(-df(0,1)+dfi(0,1))
v=(v+v1)/2
44 for j=1 to N-1
df(j,1)=Q(j,1)*df(j-1,0)+R(j,2)*dfi(j+1,0)
f(j,1)=f(j,0)+df(j,1)*dl
dfi(j,1)=R(j,1)*df(j-1,0)+Q(j,2)*dfi(j+1,0)
fi(j,1)=fi(j,0)+dfi(j,1)*dl
P(j)=-df(j,1)+dfi(j,1)
v(j)=a*(df(j-1,0)+dfi(j,1))
Qs=Qs+m*(df(j-1,1)+dfi(j,1))/va
next j
df(N,1)=df(N-1,0)
ffi(N)=fi(N,0)
45 fi(N,1)=fi(n,0)+df(N,1)*dl-b1*(f(N,1)+ffi(N))*dl
if abs(fi(N,1))>0 goto 46
goto 50
