2025
.pdf81
COMMON /M11/KM2,km3 COMMON /MK2/IPRIN,dm,ik,ik1
COMMON /M1/G,Q,E,L,XAQ,XT,R,XAD,QBEG COMMON /M3/FD0,FF0,FQ0,M0,RYD,RYQ COMMON /M4/ID0,IQ0
COMMON /M10/N,N1,N2,IgR,HQ,Q2
AM1=M0
M0=AM1-AMOB(2.)+am3(1.)
IF (Q.EQ.QBEG) GOTO 1112
IF (M0.LE.0) GO TO 1112
IF (IPRIN.EQ.1) GOTO 1112
IF (DM.LE.M0) GOTO 1112
PRINT 1111
1111 FORMAT(' MAX='$) IPRIN=1
PRINT 43,L,B1,Q,AM1,M0,KU1,U2,KM2,km3,g,r 1112 CONTINUE
IF (IGR.EQ.0) then
IK=0
IF (RAS(n1,1).LE.0.) IK=1
IK1=0
DO 505 I7=1,N
505IF (RAS(i7,1).LE.0.) IK1=1 IGR=1
else
DO 35 I=1,N
IF (RAS(i,1).LE.0.) then
IF (IK1.NE.1) then
* устойчивость
82
PRINT 40
PRINT 43,L,B2,Q2,DM1,DM,U1,U2,KM2,km3,g,r * самораскачивание
PRINT 31
PRINT 43,L,B1,Q,AM1,M0,KU1,U2,KM2,km3,g,r IK1=1
endif GOTO 36 endif
35CONTINUE
IF (IK1.EQ.1) then
* самораскачивание
PRINT 31
PRINT 43,L,B2,Q2,DM1,DM,U1,U2,KM2,km3,g,r * устойчивость
PRINT 40
PRINT 43,L,B1,Q,AM1,M0,KU1,U2,KM2,km3,g,r IK1=0
endif
36CONTINUE
IF (RAS(n1,1).LE.0.) then IF (IK.NE.1) then
IF (IK1.EQ.1) then
* самораскачивание
PRINT 31 else
* устойчивость
PRINT 40 endif
83
PRINT 43,L,B2,Q2,DM1,DM,U1,U2,KM2,km3,g,r
* сползание
PRINT 38
PRINT 43,L,B1,Q,AM1,M0,KU1,U2,KM2,km3,g,r IF (IK.EQ.0) return 1
IK=0 endif
elseif (IK.NE.0) then
* сползание
PRINT 38
PRINT 43,L,B2,Q2,DM1,DM,U1,U2,KM2,km3,g,r
* устойчивость
PRINT 40
PRINT 43,L,B1,Q,AM1,M0,KU1,U2,KM2,km3,g,r ik = 0
endif
endif RC2=RAS(n,1) RC1=RAS(n1,1) DM1=AM1 DM=M0 U1=KU1 B2=B1
Q2=Q
40 FORMAT(' УСТОЙЧИВОСТЬ')
38 FORMAT(' СПОЛЗАНИЕ')
31 FORMAT(' САМОРАСКАЧИВАНИЕ')
43 FORMAT(1X,'L=',F7.4,' B0=',F7.2,' Q=',F6.2, *' M01=',F8.4,' M0=',F8.4,' U1=',F6.2,' U2=',F6.2,
84
*' KM2=',F8.2,' KM3=',f8.2,' G=',f8.2,' R=',f8.2) RETURN
END
subroutine lever(n,a,fu)
double precision m(6),a(6,6),b(6,6),c(6,6),fu(7),s,r do1i=1,n
do2j=1, n 2 b(i,j) = 0.0
1b(i,i) = 1.0 do3 i= 1,n
s = 0.0 do4 kk=1,n
do6 ii=1,n c(kk,ii) = 0.0 do6 j=1,n
6 |
c(kk,ii) = c(kk,ii) + b(kk,j)*a(j,ii) |
4 |
s = s + c(kk,kk) |
|
do7 ii=1,n |
|
do7 j=1,n |
7 |
b(ii,j) = c(ii,j) |
3m(i) = s fu(1) = 1. fu(2) = -m(1) do8 i=2,n
kk = i
s = 0.0 do9 j=1, i-1
kk = kk - 1
85
9 |
s = s + m(kk)*fu(j+1) |
|
r = i |
8fu(i+1) = -(m(i) + s)/r do20i=1,3
j = n + 2 - i r = fu(i)
fu(i) = fu(j)
20fu(j) = r return end