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 I, J, K, N, NG, IEL, ITEL, NUMFAKECP,
     .        OFFSET_IGBUFSSG, OFFSET_NUMFAKENOD, ITSURF
C     REAL
      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 !! on est sur une surface isogeometrique
          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
 
CC NB : on peut repartir le volume differemment selon si les points sont dans les coins, au milieu
cC      ou sur une arrete, avec PX,PY,PZ on peut retrouver la place de chacun
Cc pour les elements 20 neouds, il y a 3/192 * vol aux coins et 14/192 * vol aux points milieux
C
CC  NB: il me semble que numfakenodigeo est trop grand ici et donc les tableaux sont trop grands et se terminent par 0
 
      RETURN