ig3dfint.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	ig3dfint (pm, mxt, kxig3d, sig, nctrl, matb, fx, fy, fz, vol, btdba, dba, ssp_eq, stig, nel, nft)

Function/Subroutine Documentation

ig3dfint()

subroutine ig3dfint	(		pm,
		integer, dimension(*)	mxt,
		integer, dimension(nixig3d,*)	kxig3d,
			sig,
		integer	nctrl,
		dimension(3*nctrl,mvsiz)	matb,
			fx,
			fy,
			fz,
			vol,
			btdba,
			dba,
			ssp_eq,
			stig,
		integer	nel,
		integer, intent(in)	nft )

Definition at line 28 of file ig3dfint.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NFT
      INTEGER NCTRL, KXIG3D(NIXIG3D,*), MXT(*),NEL
      my_real
     .   sig(nel,6), vol(*), pm(npropm,*)
      my_real,
     .  DIMENSION(3*NCTRL,MVSIZ) :: matb
      my_real
     .  fx(mvsiz,*),fy(mvsiz,*),fz(mvsiz,*),btdba(3*nctrl,mvsiz),
     .  dba(6,mvsiz),stig(mvsiz,*),ssp_eq(*),sumx,sumy,sumz
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, INCTRL
      my_real
     .   aa,wip(mvsiz),stin(mvsiz,nctrl)
C-----------------------------------------------
C            _
C           /   t
C     F = _/ [B]  sig dvol
C
C                   NEL
C   [B] = [dNi/dx  dNk/dx  ..  ] \  3
C         [dNi/dy  dNk/dy  ..  ] |  *
C         [dNi/dz  dNk/dz  ..  ] |  N
C         [dNj/dx  dNl/dx  ..  ] |  C
C         [dNj/dy  dNl/dy  ..  ] |  T
C         [dNj/dz  dNl/dz  ..  ] |  R
C         [  ..      ..    ..  ] /  L
C
      DO i=1,nel
       sumx=zero 
       sumy=zero
       sumz=zero
       DO j=1,nctrl
        IF( j <= kxig3d(3,i+nft) ) THEN
         inctrl = (j-1)*3
         sumx=sumx+abs(matb(inctrl+1,i))
         sumy=sumy+abs(matb(inctrl+2,i))
         sumz=sumz+abs(matb(inctrl+3,i))
        ENDIF
       ENDDO
       DO j=1,nctrl
        IF( j <= kxig3d(3,i+nft) ) THEN
         inctrl = (j-1)*3
         fx(i,j)=fx(i,j)-vol(i)*(matb(inctrl+1,i)*sig(i,1)
     +                         + matb(inctrl+2,i)*sig(i,4)  
     +                         + matb(inctrl+3,i)*sig(i,6))
 
         fy(i,j)=fy(i,j)-vol(i)*(matb(inctrl+1,i)*sig(i,4)
     +                         + matb(inctrl+2,i)*sig(i,2)
     +                         + matb(inctrl+3,i)*sig(i,5))
 
         fz(i,j)=fz(i,j)-vol(i)*(matb(inctrl+1,i)*sig(i,6)
     +                         + matb(inctrl+2,i)*sig(i,5) 
     +                         + matb(inctrl+3,i)*sig(i,3))
         stin(i,j)= half*vol(i)*
     .      (abs(matb(inctrl+1,i))*sumx 
     .     + abs(matb(inctrl+2,i))*sumy 
     .     + abs(matb(inctrl+3,i))*sumz)
         aa=pm(89,mxt(i))*ssp_eq(i)*ssp_eq(i)
         stig(i,j)=stig(i,j)+stin(i,j)*aa
        ENDIF
       ENDDO
      ENDDO
 
 
      RETURN