📄 libearingfilm.for
字号:
30 CONTINUE
K=0
40 DO 60 I=1,M
DO 50 J=1,N
PT(I,J)=P(I,J)
P(I,J)=(1.0-WP)*P(I,J)+WP*(H3(I,J)*P(I,J+1)+H4(I,J)*P(I+1,J)
1 +H5(I,J)*P(I,J-1)+H6(I,J)*P(I-1,J)-H7(I,J))
IF(P(I,J).LT.0)THEN
P(I,J)=0
ENDIF
diff(i,j)=abs(pt(i,j)-p(i,j))
S0=S0+DIFF(I,J)
S1=S1+P(I,J)
50 continue
60 CONTINUE
DP1=S0/S1
! write(*,*)dp1
IF(DP1.GT.1E-3)THEN
K=K+1
GOTO 40
ENDIF
DO 80 I=1,M
DO 70 J=1,N
PS(I,J)=P(I,J)
70 CONTINUE
80 CONTINUE
return
END
subroutine cznl(m,n,x,y,pwj,p,w,kw,o)
dimension p(0:m+1,0:n+1),kw1(m,n),o1(m,n),fi(m,0:n)
double precision p,kw1,o1,fi
common pii,xdjx,oilnd
real w,kw
DO 2 I=1,M
DO 1 J=0,N
FI(I,J)=2*PIi/3.0+X*(J-1)-PWJ
1 CONTINUE
2 CONTINUE
KW=0
O=0
KW1(I,J)=0
O1(I,J)=0
DO 220 I=1,M
DO 200 J=1,N
KW1(I,J)=abs(P(I,J)*COS(FI(I,J)+PWJ-pii)*X*Y)
O1(I,J)=abs(P(I,J)*SIN(FI(I,J)+PWJ)*X*Y)
KW=KW+KW1(I,J)
O=O+O1(I,J)
200 CONTINUE
220 CONTINUE
w=sqrt(kw**2.0+o**2.0)
return
end
subroutine input(numbe,nodal,nttyp,nutyp,U0, mtyp,xn,yn,zn,msi,
# cosbex,cosbey,cosbez,kod,nord,kodt,kode,mors,morp,
# morc,merp,merc,mord,mort,node,MNB)
COMMON/S1/PI,PR(2),PR1(2),PR2(2),PR3(2),PR4(2),CON(2),DSMIN2(2),
1 DSMAX2(2),E(2),G(2)
COMMON/S2/AXX(8),AXY(8),AXZ(8),AYX(8),AYY(8),AYZ(8),
1 AZX(8),AZY(8),AZZ(8),BXX(8),BXY(8),BXZ(8),
2 BYX(8),BYY(8),BYZ(8),BZX(8),BZY(8),BZZ(8)
COMMON/S4/HJACOB,COSAX,COSAY,COSAZ
COMMON/S6/FJACOB,COSBX,COSBY,COSBZ
COMMON/S7/NODBS(2),NUMBS(2),NODPS(2),NODCS(2),NEDPS(2),NEDCS(2)
COMMON/S8/MS(2),MD(6),NS(2),ND(2)
dimension COSBEX(2,numbe,3),COSBEY(2,numbe,3),COSBEZ(2,numbe,3)
DIMENSION NTTYP(2,msi/3),NUTYP(2,msi/3),NODAL(8),FN(8),XK(8),
1 YK(8),ZK(8)
common/s30/BTX(2000,8),BTY(2000,8),BTZ(2000,8),RTX(2,2000,8),
1 RTY(2,2000,8),RTZ(2,2000,8),RUX(2,2000,8),RUY(2,2000,8),
1 RUZ(2,2000,8)
dimension Nord(2,msi/3,8),KODT(2,msi/3),Kode(2,msi/3,8)
dimension kod(2,msi/3)
dimension Mors(2,numbe),Morp(2,numbe),Morc(2,numbe),Merp(2,numbe),
# Merc(2,numbe)
dimension Mord(2,msi/3),Mort(2,msi/3),Node(2,msi/3,8)
dimension XN(2,msi/3),YN(2,msi/3),ZN(2,msi/3)
DIMENSION NUMOS(2),MTYP(2),U0(numbe)
DIMENSION MNB(2,3000)
PI=4.*ATAN(1.E+0)
C DEFINE LOCATION,SIZES,DRENTATIONS,BOUNDARY CONDITIONS
C CONTACT COUPLES OF BOUNDARY ELEMENTS
DO 280 I=1,2
READ(7,*) NODBS(I),NUMBS(I),NUMOS(I),MTYP(I),NUMTS,NUMUS,
1 DSMIN2(I),DSMAX2(I),PR(I),E(I)
CON(I)=1./(8.*PI*(1.-PR(I)))
G(I)=0.5*E(I)/(1.+PR(I))
PR1(I)=1.-2.*PR(I)
PR2(I)=1.-4.*PR(I)
PR3(I)=3.-3.*PR(I)
PR4(I)=3.-4.*PR(I)
DO 60 J=1,NODBS(I)
READ(7,*) N,KOD(I,N),XN(I,N),YN(I,N),ZN(I,N)
60 CONTINUE
write(*,*)"wjm"
c DO 601 J=1,NODBS(I)
c READ(7,*) N
c601 CONTINUE
IF(NUMTS.LE.0) GOTO 80
DO 70 N=1,NUMTS
READ (7,*) M,(BTX(M,J),BTY(M,J),BTZ(M,J),J=1,2)
70 READ(7,*) (BTX(M,J),BTY(M,J),BTZ(M,J),J=3,4)
80 IF(NUMUS.LE.0) GOTO 100
DO 90 N=1,NUMUS
READ(7,*) M,(RUX(I,M,J),RUY(I,M,J),RUZ(I,M,J),J=1,2)
90 READ(7,*) (RUX(I,M,J),RUY(I,M,J),RUZ(I,M,J),J=3,4)
100 continue
DO 176 M=1,NUMBS(I)
READ(7,*) N,(NORD(I,N,L),L=1,4),KODT(I,N),(KODE(I,N,L),L=1,4),
# NTTYP(I,N),NUTYP(I,N)
IF(NTTYP(I,N).LE.0.OR.NTTYP(I,N).GT.NUMTS) GOTO 120
DO 110 J=1,8
RTX(I,NTTYP(I,N),J)=0.5*BTX(NTTYP(I,N),J)/G(I)
RTY(I,NTTYP(I,N),J)=0.5*BTY(NTTYP(I,N),J)/G(I)
110 RTZ(I,NTTYP(I,N),J)=0.5*BTZ(NTTYP(I,N),J)/G(I)
! WRITE(99,*) (RTX(I,NTTYP(I,N),J),J=1,8)
GOTO 140
120 NTTYP(I,N)=NUMTS+1
DO 130 J=1,8
RTX(I,NTTYP(I,N),J)=0.
RTY(I,NTTYP(I,N),J)=0.
130 RTZ(I,NTTYP(I,N),J)=0.
140 IF(NUTYP(I,N).GT.0.AND.NUTYP(I,N).LE.NUMUS) GOTO 176
NUTYP(I,N)=NUMUS+1
DO 150 J=1,8
RUX(I,NUTYP(I,N),J)=0.
RUY(I,NUTYP(I,N),J)=0.
150 RUZ(I,NUTYP(I,N),J)=0.
176 CONTINUE
WRITE(*,*)"HJL"
c DO k=1,NUMBS(I)
c DO J=1,8
c IF (NORD(I,k,J).EQ.0) GOTO 211
c KODE(I,k,J)=KOD(I,NORD(I,k,J))
c211 CONTINUE
c ENDDO
c ENDDO
CALL ASSUM(I,NUMBE,MULTI,kod,nord,kodt,mors,morp,
# morc,merp,merc,mord,mort,node,msi,MNB)
IF(DSMIN2(I).NE.0..AND.DSMAX2(I).NE.0.) GOTO 200
XK(1)=(XN(I,2)-XN(I,1))*(XN(I,2)-XN(I,1))+(YN(I,2)-YN(I,1))*
1 (YN(I,2)-YN(I,1))+(ZN(I,2)-ZN(I,1))*(ZN(I,2)-ZN(I,1))
DO 190 M=3,NODBS(I)
XK(2)=(XN(I,M)-XN(I,1))*(XN(I,M)-XN(I,1))+(YN(I,M)-YN(I,1))*
1 (YN(I,M)-YN(I,1))+(ZN(I,M)-ZN(I,1))*(ZN(I,M)-ZN(I,1))
190 IF(XK(2).GT.XK(1)) XK(1)=XK(2)
DSMIN2(I)=0.146*XK(1)
DSMAX2(I)=0.382*XK(1)
200 continue !WRITE(99,18) DSMIN2(I),DSMAX2(I)
! WRITE(99,19)
DO 210 N=1,NODBS(I)
! WRITE(99,20) N,KOD(I,N),XN(I,N),YN(I,N),ZN(I,N)
210 WRITE(*,20) N,KOD(I,N),XN(I,N),YN(I,N),ZN(I,N)
IF(NUMTS.LE.0) GOTO 220
! WRITE(99,21)
DO 215 M=1,NUMTS
! WRITE(99,22) M,(BTX(M,J),BTY(M,J),BTZ(M,J),J=1,4)
WRITE(*,22)M,(BTX(M,J),BTY(M,J),BTZ(M,J),J=1,4)
215 continue !WRITE(99,23) (BTX(M,J),BTY(M,J),BTZ(M,J),J=5,8)
220 IF(NUMUS.LE.0) GOTO 230
! wRITE(99,24)
DO 225 M=1,NUMUS
! WRITE(99,22) M,(RUX(I,M,J),RUY(I,M,J),RUZ(I,M,J),J=1,4)
225 continue !WRITE(99,23) (RUX(I,M,J),RUY(I,M,J),RUZ(I,M,J),J=5,8)
230 continue !WRITE(99,25)
DO 240 M=1,NUMBS(I)
! WRITE(*,26)M,KODT(I,M),(NORD(I,M,J),J=1,8),(KODE(I,M,J),J=1,8)
240 continue !WRITE(99,26) M,KODT(I,M),(NORD(I,M,J),J=1,8),(KODE(I,M,J),J=1,8)
! WRITE(99,27)
DO 280 M=1,NUMBS(I)
DO 260 J=1,8
NODAL(J)=NORD(I,M,J)
IF(NODAL(J).EQ.0) GOTO 250
FN(J)=0.5
XK(J)=XN(I,NODAL(J))
YK(J)=YN(I,NODAL(J))
ZK(J)=ZN(I,NODAL(J))
GOTO 260
250 FN(J)=0.
XK(J)=0.
YK(J)=0.
ZK(J)=0.
260 CONTINUE
FN(1)=0.25-0.5*(FN(5)+FN(8))
FN(2)=0.25-0.5*(FN(5)+FN(6))
FN(3)=0.25-0.5*(FN(6)+FN(7))
FN(4)=0.25-0.5*(FN(7)+FN(8))
XM=FN(1)*XK(1)+FN(2)*XK(2)+FN(3)*XK(3)+FN(4)*XK(4)+
1 FN(5)*XK(5)+FN(6)*XK(6)+FN(7)*XK(7)+FN(8)*XK(8)
YM=FN(1)*YK(1)+FN(2)*YK(2)+FN(3)*YK(3)+FN(4)*YK(4)+
1 FN(5)*YK(5)+FN(6)*YK(6)+FN(7)*YK(7)+FN(8)*YK(8)
ZM=FN(1)*ZK(1)+FN(2)*ZK(2)+FN(3)*ZK(3)+FN(4)*ZK(4)+
1 FN(5)*ZK(5)+FN(6)*ZK(6)+FN(7)*ZK(7)+FN(8)*ZK(8)
CALL JACOB(NODAL,0.0,0.0,XK,YK,ZK,0)
ANGLX=180./PI*ACOS(COSBX)
ANGLY=180./PI*ACOS(COSBY)
ANGLZ=180./PI*ACOS(COSBZ)
NBT=NTTYP(I,M)
NBU=NUTYP(I,M)
IF(NBT.GT.NUMTS) NBT=0
IF(NBU.GT.NUMUS) NBU=0
280 continue !WRITE(99,28) M,XM,YM,ZM,ANGLX,ANGLY,ANGLZ,NBT,NBU
DO 300 LM=1,NUMBE
c if(LM.LE.140.OR.LM.GT.NUMBE-140)THEN
c U0(LM)=-0.6000
c ELSEIF((LM.GT.140.AND.LM.LE.210).OR.(LM.GT.NUMBE-210.AND.LM.LE
c #.NUMBE-140))THEN
c U0(LM)=-1.8000
c ELSEc U0(LM)=-1.80000
c ENDIF
U0(LM)=0.01
300 CONTINUE
CALL LOCALE(NUMBE,xn,yn,zn,msi,cosbex,cosbey,cosbez,
# nord,mors,kodt)
!----COMPUTE INFLUENCE COEFFICIENTS AND UP SYSTEM OF ALGEBRIC EQUATIONS
38 Print *, 'End checking data !'
2 FORMAT(4I5)
11 FORMAT(6I5,4F10.0)
13 FORMAT(2I5,3F10.0,I5)
14 FORMAT(I5,6F10.0,/,F15.0,5F10.0)
15 FORMAT(F15.0,5F10.0)
16 FORMAT(9I5,1X,8I2,4I3)
18 FORMAT(//,35H INTEGRAL CALCULATION PARAMETER.,//,9H SQUARE ,
1 22HOF MINMUN DISTANCE = ,F12.4,//,26H SQUARE OF MAXMUM DISTAN,
2 5HCE = ,F12.4)
19 FORMAT(///,/,21H BOUNDARY NODES DATA,//18H NODE KODE ,
1 43H X_COORDINATE Y_COORDINATE Z_COORDINATE,/)
20 FORMAT(2I7,3F15.4)
21 FORMAT(///,/,29H INITIAL BOUNDARY STRESS DATA,/,12H TYPE ,
1 4HTX_1,6X,4HTY_1,6X,4HTZ_1,6X,4HTX_2,6X,4HTY_2,6X,4HTZ_2,6X,2HTX,
2 2H_3,6X,4HTY_3,6X,4HTZ_3,6X,4HTX_4,6X,4HTY_4,6X,4HTZ_4,/,12X,1HT,
3 3HX_5,6X,4HTY_5,6X,4HTZ_5,6X,4HTX_6,6X,4HTY_6,6X,4HTZ_6,6X,3HTX_,
4 1H7,6X,4HTY_7,6X,4HTZ_7,6X,4HTX_8,6X,4HTY_8,6X,4HTZ_8,6X,/)
22 FORMAT(I5,2X,12F10.2)
23 FORMAT(7X,12F10.2)
24 FORMAT(//,36H INITIAL BOUNDARY DISPLACEMENT DATA.,/, 9H TYPE ,
1 3X,
1 4HUX_1,6X,4HUY_1,6X,4HUZ_1,6X,4HUX_2,6X,4HUY_2,6X,4HUZ_2,6X,
2 4HUX_3,6X,4HUY_3,6X,4HUZ_3,6X,4HUX_4,6X,4HUY_4,6X,4HUZ_4,/,12X,
3 4HUX_5,6X,4HUY_5,6X,4HUZ_5,6X,4HUX_6,6X,4HUY_6,6X,4HUZ_6,6X,
4 4HUX_7,6X,4HUY_7,6X,4HUZ_7,6X,4HUX_8,6X,4HUY_8,6X,4HUZ_8,6X,/)
25 FORMAT(///,/,23H BOUNDARY ELEMENT DATA.,//,17H ELEMENT CODE EL ,
1 60HNODE_1 NODE_2 NODE_3 NODE_4 NODE_5 NODE_6 NODE_7 NODE_8 KODE,
2 51H_1 KODE_2 KODE_3 KODE_4 KODE_5 KODE_6 KODE_7 KODE_8,/)
26 FORMAT(I6,I8,16I7)
27 FORMAT(///,22H BOUNDARY ELEMENT DATA,//,20H ELEM X_(CENTER) Y_(,
1 60HCENTER) Z_(CENTER) ANGLE(N,X) ANGLE(N,Y) ANGLE(N,Z) CASE(T) ,
2 7HCASE(U),/)
28 FORMAT(I4,3F10.4,4X,3F10.4,2X,I5,I8)
CLOSE(7)
return
end subroutine
subroutine matrix(iterl,multi,numbe,nodal,nttyp,nutyp,U0,mtyp,
# msi,nsi,ksi,xn,yn,zn,cosbex,cosbey,cosbez,kod,nord,kodt,
# kode,mors,morp,morc,merp,merc,mord,mort,node,MNB,film,TFORCE)
COMMON/S1/PI,PR(2),PR1(2),PR2(2),PR3(2),PR4(2),CON(2),DSMIN2(2),
1 DSMAX2(2),E(2),G(2)
COMMON/S2/AXX(8),AXY(8),AXZ(8),AYX(8),AYY(8),AYZ(8),
1 AZX(8),AZY(8),AZZ(8),BXX(8),BXY(8),BXZ(8),
2 BYX(8),BYY(8),BYZ(8),BZX(8),BZY(8),BZZ(8)
COMMON/S4/HJACOB,COSAX,COSAY,COSAZ
COMMON/S6/FJACOB,COSBX,COSBY,COSBZ
COMMON/S7/NODBS(2),NUMBS(2),NODPS(2),NODCS(2),NEDPS(2),NEDCS(2)
COMMON/S8/MS(2),MD(6),NS(2),ND(2)
dimension COSBEX(2,numbe,3),COSBEY(2,numbe,3),COSBEZ(2,numbe,3)
DIMENSION NTTYP(2,msi/3),NUTYP(2,msi/3),NODAL(8),XK(8),YK(8),ZK(8)
common/s30/BTX(2000,8),BTY(2000,8),BTZ(2000,8),RTX(2,2000,8),
1 RTY(2,2000,8),RTZ(2,2000,8),RUX(2,2000,8),RUY(2,2000,8),
1 RUZ(2,2000,8)
dimension A(msi,nsi),P(msi),AC(2,3*numbe,ksi),PC(2,3*numbe)
DIMENSION A1(msi,nsi),p1(msi)
dimension XN(2,msi/3),YN(2,msi/3),ZN(2,msi/3),
1film(9,19),TFORCE(9,19)
Dimension SUMX(msi),SUMY(msi),SUMZ(msi),U0(numbe),mtyp(2)
dimension Nord(2,msi/3,8),KODT(2,msi/3),Kode(2,msi/3,8)
dimension kod(2,msi/3)
dimension Mors(2,numbe),Morp(2,numbe),Morc(2,numbe),Merp(2,numbe),
# Merc(2,numbe)
dimension Mord(2,msi/3),Mort(2,msi/3),Node(2,msi/3,8)
DIMENSION MNB(2,3000)
DOUBLE PRECISION film,TFORCE,SUMX
INTEGER S
Do 988 i=1,2
! WRITE(99,1001) I
WRITE(*,1001) I
WRITE(*,12) I,MTYP(I),PR(I),E(I)
c write(*,1002) nodbs(i),numbs(i),i,ms(i),i,ns(i),i,nd(i),numbe
1001 FORMAT(1X,5HBODY=,I2)
12 FORMAT(5X,19HInput initial data./5x,18H Body of contact : ,
# I4,/5x,19H Type of the body: ,I3,/5x,18H Poisson's ratio: ,F5.2,
# /5x,18H Young's modulus: ,E11.4,/)
1002 format(6x,10Hnodbs: ,i5/6x,10Helements: ,i5/6x,3hMs(,i2,3h) =,
# i7/6x,3hNs(,i2,3h) =,i7/6x,3hNd(,i2,3h) =,i7/6x,8hNumbe =,i6)
311 DO 310 K=1,MS(I)
P(K)=0.
SUMX(K)=0.
SUMY(K)=0.
SUMZ(K)=0.
DO 310 L=1,NS(I)
310 A(K,L)=0.
DO 600 J=1,NODBS(I)
WRITE(*,*) J
! WRITE(99,*) J
JZ=3*J
JY=JZ-1
JX=JZ-2
DO 460 K=1,NUMBS(I)
CALL INITL(2)
L=0
N=NTTYP(I,K)
S=NUTYP(I,K)
DO 330 M=1,8
NODAL(M)=NORD(I,K,M)
IF(NODAL(M).EQ.0) GOTO 320
IF(NODAL(M).EQ.J) L=M
XK(M)=XN(I,NODAL(M))
YK(M)=YN(I,NODAL(M))
ZK(M)=ZN(I,NODAL(M))
GOTO 330
320 XK(M)=0.
YK(M)=0.
ZK(M)=0.
330 CONTINUE
c PRINT *,'i=',i,'k=',k,'j=',j
IF(L.NE.0) GOTO 340
CALL COEFF(I,NODAL,XN(I,J),YN(I,J),ZN(I,J),XK,YK,ZK)
GOTO 350
340 CALL IMPLE(I,NODAL,L,XN(I,J),YN(I,J),ZN(I,J),XK,YK,ZK)
350 DO 460 M=1,8
IF(NORD(I,K,M).EQ.0) GOTO 450
L=MORD(I,NORD(I,K,M))
KZ=3*L
KY=KZ-1
KX=KZ-2
GOTO (365,368,370,372,375,378,380,382,385,388,390,392,395,398,
1 400,402,405,408,410,412,415,418,420,422,425,430,435,440,
1 440),KODE(I,K,M)
365 A(JX,KX)=A(JX,KX)+AXX(M)
A(JX,KY)=A(JX,KY)+AXY(M)
A(JX,KZ)=A(JX,KZ)+AXZ(M)
A(JY,KX)=A(JY,KX)+AYX(M)
A(JY,KY)=A(JY,KY)+AYY(M)
A(JY,KZ)=A(JY,KZ)+AYZ(M)
A(JZ,KX)=A(JZ,KX)+AZX(M)
A(JZ,KY)=A(JZ,KY)+AZY(M)
A(JZ,KZ)=A(JZ,KZ)+AZZ(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXY(M)*RTY(I,N,M)+BXZ(M)*RTZ(I,N,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYY(M)*RTY(I,N,M)+BYZ(M)*RTZ(I,N,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZY(M)*RTY(I,N,M)+BZZ(M)*RTZ(I,N,M)
GOTO 450
368 A(JX,KX)=A(JX,KX)-BXX(M)
A(JX,KY)=A(JX,KY)-BXY(M)
A(JX,KZ)=A(JX,KZ)-BXZ(M)
A(JY,KX)=A(JY,KX)-BYX(M)
A(JY,KY)=A(JY,KY)-BYY(M)
A(JY,KZ)=A(JY,KZ)-BYZ(M)
A(JZ,KX)=A(JZ,KX)-BZX(M)
A(JZ,KY)=A(JZ,KY)-BZY(M)
A(JZ,KZ)=A(JZ,KZ)-BZZ(M)
P(JX)=P(JX)-AXX(M)*RUX(I,S,M)-AXY(M)*RUY(I,S,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)-AYX(M)*RUX(I,S,M)-AYY(M)*RUY(I,S,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)-AZX(M)*RUX(I,S,M)-AZY(M)*RUY(I,S,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
370 A(JX,KX)=A(JX,KX)+AXX(M)
A(JX,KY)=A(JX,KY)+AXY(M)
A(JX,KZ)=A(JX,KZ)-BXZ(M)
A(JY,KX)=A(JY,KX)+AYX(M)
A(JY,KY)=A(JY,KY)+AYY(M)
A(JY,KZ)=A(JY,KZ)-BYZ(M)
A(JZ,KX)=A(JZ,KX)+AZX(M)
A(JZ,KY)=A(JZ,KY)+AZY(M)
A(JZ,KZ)=A(JZ,KZ)-BZZ(M)
P(JX)=P(JX)+BXX(M)*RTX(I,N,M)+BXY(M)*RTY(I,N,M)-AXZ(M)*RUZ(I,S,M)
P(JY)=P(JY)+BYX(M)*RTX(I,N,M)+BYY(M)*RTY(I,N,M)-AYZ(M)*RUZ(I,S,M)
P(JZ)=P(JZ)+BZX(M)*RTX(I,N,M)+BZY(M)*RTY(I,N,M)-AZZ(M)*RUZ(I,S,M)
GOTO 450
372 A(JX,KX)=A(JX,KX)-BXX(M)
A(JX,KY)=A(JX,KY)-BXY(M)
A(JX,KZ)=A(JX,KZ)+AXZ(M)
A(JY,KX)=A(JY,KX)-BYX(M)
A(JY,KY)=A(JY,KY)-BYY(M)
A(JY,KZ)=A(JY,KZ)+AYZ(M)
A(JZ,KX)=A(JZ,KX)-BZX(M)
A(JZ,KY)=A(JZ,KY)-BZY(M)
A(JZ,KZ)=A(JZ,KZ)+AZZ(M)
P(JX)=P(JX)+BXZ(M)*RTZ(I,N,M)-AXX(M)*RUX(I,S,M)-AXY(M)*RUY(I,S,M)
P(JY)=P(JY)+BYZ(M)*RTZ(I,N,M)-AYX(M)*RUX(I,S,M)-AYY(M)*RUY(I,S,M)
P(JZ)=P(JZ)+BZZ(M)*RTZ(I,N,M)-AZX(M)*RUX(I,S,M)-AZY(M)*RUY(I,S,M)
GOTO 450
375 A(JX,KX)=A(JX,KX)-BXX(M)
A(JX,KY)=A(JX,KY)+AXY(M)
A(JX,KZ)=A(JX,KZ)+AXZ(M)
A(JY,KX)=A(JY,KX)-BYX(M)
A(JY,KY)=A(JY,KY)+AYY(M)
A(JY,KZ)=A(JY,KZ)+AYZ(M)
A(JZ,KX)=A(JZ,KX)-BZX(M)
A(JZ,KY)=A(JZ,KY)+AZY(M)
A(JZ,KZ)=A(JZ,KZ)+AZZ(M)
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -