bulkfakeigeo3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	bulkfakeigeo3 (elbuf_tab, iparg, pm, kxig3d, igrsurf, stifint)

Function/Subroutine Documentation

bulkfakeigeo3()

subroutine bulkfakeigeo3	(	type(elbuf_struct_), dimension(ngroup), target	elbuf_tab,
		integer, dimension(nparg,ngroup)	iparg,
			pm,
		integer, dimension(nixig3d,*)	kxig3d,
		type (surf_), dimension(nsurf)	igrsurf,
			stifint )

Definition at line 29 of file bulkfakeigeo3.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_mod
      USE groupdef_mod
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"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IPARG(NPARG,NGROUP),
     .        KXIG3D(NIXIG3D,*)
C     REAL
      my_real
     .   pm(npropm,*),stifint(*)
      TYPE(ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
      TYPE (SURF_)   , DIMENSION(NSURF)   :: IGRSURF
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: J, K, N, NG, IEL, ITEL, NUMFAKECP, ITSURF
      my_real :: av1,cv1, vol(numelig3d), vnige(numfakenodigeo), bnige(numfakenodigeo)
      TYPE(G_BUFEL_) ,POINTER :: GBUF 
C-----------------------------------------------------------------------
C 
C     isogeometric elements :
c
      DO ng=1,ngroup
       IF (iparg(5,ng) == 101) THEN
        DO itel=1,iparg(2,ng)
          gbuf => elbuf_tab(ng)%GBUF
          vol(itel+iparg(3,ng)) = gbuf%VOL(itel)
        ENDDO
       ENDIF
      ENDDO
 
      vnige = 0
      bnige = 0
c
      DO itsurf=1,nsurf
        IF(igrsurf(itsurf)%NSEG_IGE>0) THEN !we are on an isogeometric surface
          DO j=1,igrsurf(itsurf)%NSEG_IGE
            iel=igrsurf(itsurf)%ELEM_IGE(j)
            av1 = vol(iel)/56
            cv1 = pm(32,kxig3d(1,iel))*av1
            DO k=1,4
              numfakecp=igrsurf(itsurf)%NODES_IGE(j,k) 
              vnige(numfakecp-numnod)=vnige(numfakecp-numnod)+av1
              bnige(numfakecp-numnod)=bnige(numfakecp-numnod)+cv1
            ENDDO
          ENDDO
        ENDIF
      ENDDO
c
      DO n=1,numfakenodigeo
        bnige(n)=bnige(n)/max(em30,vnige(n))
        stifint(numnod+n)=stifint(numnod+n)+bnige(n)*(2.*vnige(n))**third
      ENDDO
 
Cnote: the volume can be distributed differently depending on whether the points are at the corners, in the middle
Cor on an edge, with px, py, pz you can find the place of each one
CThis is for 20-node elements, there are 3/192 * volume at corners and 14/192 * volume at the mid-nodes
C
Cnote: it seems that numfakenodigeo is too large here and so the arrays are too large and end with 0
 
      RETURN