fsigtrini.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	fsigtrini (fxbelm, iparg, x, pm, ixt, geo, fxbmod, fxbsig, nelt)

Function/Subroutine Documentation

fsigtrini()

subroutine fsigtrini	(	integer, dimension(*)	fxbelm,
		integer, dimension(nparg,*)	iparg,
			x,
			pm,
		integer, dimension(nixt,*)	ixt,
			geo,
			fxbmod,
			fxbsig,
		integer	nelt )

Definition at line 29 of file fsigtrini.F.

      use element_mod , only : nixt
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 FXBELM(*), IPARG(NPARG,*), IXT(NIXT,*), NELT
      my_real
     .        fxbsig(*), x(3,*), pm(npropm,*), fxbmod(*),
     .        geo(npropg,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IG, OFFSET, LAST, NFT, NFS, I, NG, IEL,
     .        MAT, PROP, N1, N2
      my_real
     .        area, e, x1, y1, z1, x2, y2, z2, l0, d11, d12, d13,
     .        d21, d22, d23, xc1, yc1, zc1, xc2, yc2, zc2, lc
C-----------------------------------------------
C
      DO ig=1,nelt,mvsiz
         offset=ig-1
         last=min(mvsiz,nelt-offset)
         nft=offset*7
         nfs=offset*2
         DO i=1,last
            ng=fxbelm(nft+7*(i-1)+1)
            iel=iparg(3,ng)+fxbelm(nft+7*(i-1)+2)
            mat=ixt(1,iel)
            prop=ixt(4,iel)
            area=geo(1,prop)
            e=pm(20,mat)
            x1=x(1,ixt(2,iel))
            y1=x(2,ixt(2,iel))
            z1=x(3,ixt(2,iel))
            x2=x(1,ixt(3,iel))
            y2=x(2,ixt(3,iel))
            z2=x(3,ixt(3,iel))
            l0=sqrt((x2-x1)**2+(y2-y1)**2+(z2-z1)**2)
            n1=fxbelm(nft+7*(i-1)+3)
            n2=fxbelm(nft+7*(i-1)+4)
            d11=fxbmod(6*(n1-1)+1)
            d12=fxbmod(6*(n1-1)+2)
            d13=fxbmod(6*(n1-1)+3)
            d21=fxbmod(6*(n2-1)+1)
            d22=fxbmod(6*(n2-1)+2)
            d23=fxbmod(6*(n2-1)+3)
            xc1=x1+d11
            yc1=y1+d12
            zc1=z1+d13
            xc2=x2+d21
            yc2=y2+d22
            zc2=z2+d23
            lc=sqrt((xc2-xc1)**2+(yc2-yc1)**2+(zc2-zc1)**2)
C
            fxbsig(nfs+2*(i-1)+1)=e*area*(lc-l0)
            fxbsig(nfs+2*(i-1)+2)=half*e*area*(lc-l0)**2
         ENDDO
      ENDDO
C
      RETURN