边界元程序输入模板.doc

C COMPUTER POTENTLAL OR ITS DERIVIAL OF THE BOUNDARY AND INITIAL POINTS C IMPLICIT REAL*8(A-H,O-Z) COMMON/S1/PI COMMON/S5/GI(6,3),OME(6,3) COMMON/S6/FJACOB,COSBX,COSBY,COSBZ COMMON/S3/YY(242) COMMON/S71/MTYPE,NORD(240,4) COMMON/S72/X(242),y(242),z(242) common/s10/ICOD(240,4),KCD(242),FI(242),TK(242),TQ(240,4) TKT(242),TQT(242),ML(242) COMMON/S11/AT(4),BT(4),SUM(1) COMMON/SO/XX(29282) COMMON/INTE/L,CX(20),CY(20),CZ(20),SOL(20) DIMENSION XK(4),YK(4),ZK(4),ND(4),EQU(242),C(1,242),INDEX(242) CALL INPUT DO 1001 J=1,NODBS XP=X(J) YP=Y(J) ZP=Z(J) SUM(1)=0 YY(J)=0 DO 1002 K=1,NODBS EQU(K)=0. 1002 C(1,K)=0 DO 1003 K=1,NUMBS LL=0 CALL INITL DO 1004=1,4 IF(NORD(K,I).EQ.J)LL=1 ND(I)=I XK(I)=X(NORD(K,I)) YK(I)=Y(NORD(K,I)) ZK(I)=Z(NORD(K,I)) 1004 CONTINUE IF(LL.NE.0)GOTO 1006 CALL COFT(ND,XP,YP,ZP,XK,YK,ZK) GOTO 1015 1006 CALL IMPT(LL,ND,XP,YP,ZP,XK,YK,ZK) 1015 CONTINUE DO 1003 I=1,4 KI=NORD(K,I) GOTO(1375,1380,1445),ICOD(K,I) 1375 C(1,KI)=C(1,KI)-BT(I) YY(J)=YY(J)-AT(I)*TK(K,I) GOTO 1016 1480 C(1,KI)=C(1,KI)+AT(I) YY(J)=YY(J)+BT(I)*TQ(K,I) GOTO 1016 1445 YY(J)=YY(J)-AT(I)*TK(K,I)+BT(I)*TQ(K,I) 1016 SUM(1)=SUM(1)-AT(I) 1003 CONTINUE SUM(1)=SUM(1)+DBLE(MTYPE) IF(KCD(J).EQ.1)GOTO 1020 YY(J)=YY(J)-SUM(1)*FI(J) GOTO 1018 1020 (1,J)=SUM(1) 1018 DO 1019 I1=1,NODBS 1019 EQU(I1)=C(1,I1) CALL SOLVER(EQU,YY,J,NODBS,INDEX,1,1,0D-6) WRITE*(,)`BEM3 CALL SOLVER NUMBS YY`,J,YY(J) 1001 CONTINUE DO 1021 KI=1,NUMBS DO 1021 KJ=1,4 IJ=NORD(KI,KJ) GOTO (1025,1022,1021),ICOD(KI,KJ) 1025 TQ(KI,KJ)=YY(IJ) 1021 CONTINUE DO 1301 NJ=1,NODBS TKT(NJ)=0 TQT(NJ)=0 ML(NJ)=0 DO 1300 KI=1,SUMDS DO 1300 KJ=1,4 IF(NJ.EQ.NORD(KI,KJ))THEN ML(NJ)=ML(NJ)+1 TKT(NJ)=TKT(NJ)+TK(KI,KJ) TQT(NJ)=TQT(NJ)=TQ(KI,KJ