wjmscr.for

采用由FORTRAN编写,边界元前处理程序!

   	Program main	
	dimension PR(2),CON(2),DSMIN2(2),DSMAX2(2),E(2),G(2)
	dimension NODBS(2),NUMBS(2),NODPS(2),NUMOS(2),
     1	NODCS(2),NEDPS(2),NEDCS(2)
 	COMMON/S8/MS(2),MD(6),NS(2),ND(2)
      dimension XN(2,3000),YN(2,3000),ZN(2,3000),kod(2,3000) 
	DIMENSION TITLE(20),corx(8),cory(8),corz(8),mtyp(2)
	DIMENSION NTTYP(2,3000),NUTYP(2,3000)
	dimension Nord(2,3000,8),Kodt(2,3000),Kode(2,3000,8)
!	integer  NUMTS,NUMUS,MODEX,ITERL,MULTI,node1,ment1,node2,ment2
	
	open(7,file='WWW.dat',status='old')
      open(100,file='drawz2502b.scr',status='old')
	open(101,file='drawg2502b.scr',status='old')

	READ(7,1) (TITLE(I),I=1,20)
 	READ(7,*) MODEX,ITERL,MULTI,node1,ment1,node2,ment2

	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)
	DO 60 J=1,NODBS(I)
	READ(7,*) N,KOD(I,N),XN(I,N),YN(I,N),ZN(I,N)
60	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)
  	write(*,*) N,(NORD(I,N,L),L=1,4),KODT(I,N),(KODE(I,N,L),L=1,4),
     #            NTTYP(I,N),NUTYP(I,N)

176   CONTINUE
	DO 150 K=1,numbs(2)
 	DO 210 M=1,8
	
! Calculate values of node coordinates after distorded.
	corx(m)=xn(2,nord(2,k,m))
      cory(m)=yn(2,nord(2,k,m))
	corz(m)=zn(2,nord(2,k,m))
-----------------------------------------------------
! 写绘制锥套变形后图形的命令文件(*.scr)						 1125
!	write(11 ,960) k,((ul(k,m,i),i=1,3),m=1,4)
	write(100,970) (corx(m),cory(m),corz(m),m=1,4)	

!     corx(m) = xn(1,nord(1,k,m))
!	cory(m) = yn(1,nord(1,k,m))
!   	corz(m) = zn(1,nord(1,k,m))

!-----------------------------------------------------
! 写绘制轴承座变形后图形的命令文件(*.scr)
! 	write(11 ,960) k,((ul(k,m,i),i=1,3),m=1,4)
!	write(101,970) (corx(m),cory(m),corz(m),m=1,4)
210   CONTINUE

!------------------------------------------------------
150   CONTINUE



280   CONTINUE
1	FORMAT(20A4)
		


!960	Format(1x,i4,6(3x,e12.6)/5x,6(3x,e12.6))	
970   format(7H3DFACE ,SP,e12.6,1h,,SP,e12.6,1h,,SP,e12.6,1x,sp,e12.6,  
     #  ',',sp,e12.6,',',sp,e12.6,1x,sp,e12.6,',',sp,e12.6,',',sp,e12.6,     
     #  1x,sp,e12.6,',',sp,e12.6,',',sp,e12.6,' ')

      
 	CLOSE(7)
	Close(100)	
	Close(101)
	END