dmasanif.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "param_c.inc"
#include "task_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	dmasanif (elbuf_tab, x, d, geo, iparg, ixt, ixp, ixr, mas, pm, el2fa, nbf)

Function/Subroutine Documentation

dmasanif()

subroutine dmasanif	(	type (elbuf_struct_), dimension(ngroup), target	elbuf_tab,
			x,
			d,
			geo,
		integer, dimension(nparg,*)	iparg,
		integer, dimension(nixt,*)	ixt,
		integer, dimension(nixp,*)	ixp,
		integer, dimension(nixr,*)	ixr,
			mas,
			pm,
		integer, dimension(*)	el2fa,
		integer	nbf )

Definition at line 31 of file dmasanif.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_mod
      use element_mod , only : nixt,nixr,nixp
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "param_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real :: mas(*) ,pm(npropm,*),geo(npropg,*),x(3,*), d(3,*)
      INTEGER IPARG(NPARG,*),IXT(NIXT,*),IXP(NIXP,*),IXR(NIXR,*),EL2FA(*),NBF
C
      TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real
     .   off, p, vonm2, vonm, s1, s2, s12, s3, VALUE,thk0,a0,al0,
     .   rho0,xx1,xx2,xx3,yy1,yy2,yy3,zz1,zz2,zz3
      INTEGER I, NG, NEL, NFT, IAD, ITY, LFT, 
     .        IADD, N, J, LLT, MLW, NB1, NB2, NB3, NB4, NB5,NB6, 
     .        NB7, NB8, NB9, NB10, NB11, NB12, NB13, NB14, NB15, NB16,
     .        ISTRAIN,NN, K1, K2,JTURB,MT,JALE, IMID, IALEL,IPID,
     .        N1,N2,N3,N4,
     .        NN1,NN2,NN3,NN4,NN5,NN6,NN7,NN8,NN9,NN10,
     .        OFFSET,NEL_OLD,ITY_OLD,NFT_FA,N_FA,
     .        NUVAR
 
C
      TYPE(G_BUFEL_) ,POINTER :: GBUF
C-----------------------------------------------
      nn1 = 1
      nn2 = 1
      nn3 = 1
      nn4 = nn3 
      nn5 = nn4 
      nn6 = nn5 
      nn7 = nn6 + numelt
      nn8 = nn7 + numelp
      nn9 = nn8 + numelr
      nn10= nn9
C-----------------------------------------------
      nel_old = 0
      ity_old = 0
      DO ng=1,ngroup
        mlw   =iparg(1,ng)
        nel   =iparg(2,ng)
        ity   =iparg(5,ng)
        gbuf => elbuf_tab(ng)%GBUF
        IF (ispmd == 0) THEN
          IF (ity /= ity_old) THEN
            nel_old = 0
            ity_old= ity
          ENDIF
          nft_fa = nel_old
          nel_old = nel_old + nel
        ENDIF
        nft   =iparg(3,ng)
        iad   =iparg(4,ng)
        lft=1
        llt=nel
        IF (ispmd == 0) THEN
          nft_fa = nel_old - nel 
        ELSE
          nft_fa = nft
        ENDIF
C-----------------------------------------------
C       TRUSS
C-----------------------------------------------
        IF (ity == 4) THEN
          DO i=lft,llt
            n = i + nft
            n_fa = i + nft_fa
            rho0 = pm(1,ixt(1,n))
            a0 = geo(1,ixt(4,n))
            n1 = ixt(2,n)
            n2 = ixt(3,n)
            xx1 = x(1,n2)-d(1,n2)-x(1,n1)+d(1,n1)
            yy1 = x(2,n2)-d(2,n2)-x(2,n1)+d(2,n1)
            zz1 = x(3,n2)-d(3,n2)-x(3,n1)+d(3,n1)
            al0  = sqrt(xx1*xx1 + yy1*yy1 + zz1*zz1)
            mas(el2fa(nn6+n_fa)) = rho0*al0*a0
          ENDDO
C-----------------------------------------------
C       POUTRES
C-----------------------------------------------
        ELSEIF (ity == 5) THEN
          DO i=lft,llt
            n = i + nft
            n_fa = i + nft_fa
            rho0 = pm(1,ixp(1,n))
            a0 = geo(1,ixp(5,n))
            n1 = ixp(2,n)
            n2 = ixp(3,n)
            xx1 = x(1,n2)-d(1,n2)-x(1,n1)+d(1,n1)
            yy1 = x(2,n2)-d(2,n2)-x(2,n1)+d(2,n1)
            zz1 = x(3,n2)-d(3,n2)-x(3,n1)+d(3,n1)
            al0  = sqrt(xx1*xx1 + yy1*yy1 + zz1*zz1)
            mas(el2fa(nn7+n_fa)) = rho0*al0*a0
          ENDDO
C-----------------------------------------------
C       RESSORTS
C-----------------------------------------------
        ELSEIF (ity == 6) THEN
          IF(mlw==3)THEN
            DO i=lft,llt
              n = i + nft
              n_fa = i + nft_fa
              mas(el2fa(nn8+n_fa))   = half*geo(1,ixr(1,n))
              mas(el2fa(nn8+n_fa)+1) = half*geo(1,ixr(1,n))
            ENDDO
          ELSEIF (mlw == 5) THEN
c           NB1 =IAD - 1
c           NB2 =NB1  + NEL
c           NUVAR =  NINT(GEO(25,IXR(1,1+NFT)))
c           NB3 =NB2  + 3*NEL
c           NB4 =NB3  + NEL
c           NB5 =NB4  + 3*NEL
c           NB6 =NB5  
c           NB7 =NB6  
c           NB8 =NB7  
c           NB9 =NB8  + 3*NEL
c           NB10=NB9  + 3*NEL
c           NB11=NB10 
c           NB12=NB11 
c           NB13=NB12 
c           NB14=NB13 
c           NB15 = NB14 + 3*NEL
c           NB16 = NB15 + NUVAR*NEL
            DO i=lft,llt
              n = i + nft
              n_fa = i + nft_fa
              mas(el2fa(nn8+n_fa)) = gbuf%MASS(i)
cc            MAS(EL2FA(NN8+N_FA)) = BUFEL(NB16+I)
            ENDDO
          ELSE
            DO i=lft,llt
              n = i + nft
              n_fa = i + nft_fa
              mas(el2fa(nn8+n_fa)) = geo(1,ixr(1,n))
            ENDDO
          ENDIF ! IF(MLW)
        ENDIF ! IF (ITY)
C-----------------------------------------------
C       end of loop over offsets
C-----------------------------------------------
      ENDDO
C-----------------------------------------------
      RETURN