bulkigeo3.F File Reference

#include "implicit_f.inc"
#include "vect01_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	bulkige3 (volu, nctrl, mat, pm, vnige, bnige, px, py, pz, nctrlmax)

Function/Subroutine Documentation

bulkige3()

subroutine bulkige3	(		volu,
		integer	nctrl,
		integer, dimension(*)	mat,
			pm,
			vnige,
			bnige,
		integer	px,
		integer	py,
		integer	pz,
		integer	nctrlmax )

Definition at line 28 of file bulkigeo3.F.

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      "vect01_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NCTRL, NCTRLMAX, MAT(*), PX,PY,PZ
C     REAL
      my_real
     .   vnige(nctrlmax,*),bnige(nctrlmax,*),volu(*),
     .   pm(npropm,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, N
C     REAL
      my_real
     .   av1,cv1
C-----------------------------------------------------------------------
C
      IF(nctrl==0) RETURN
C 
C     isogeometric elements :
      DO i=lft,llt
        av1 = volu(i)/nctrl
        cv1 = pm(32,mat(i))*av1
        DO n=1,nctrl
          vnige(n,i)=av1
          bnige(n,i)=cv1
        ENDDO
      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
      RETURN