mmat_h1.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	mmat_h1 (hh, hh1, fac, icp, ipredu, nel, mtn, ismstr, jhbe)

Function/Subroutine Documentation

mmat_h1()

subroutine mmat_h1	(		hh,
			hh1,
			fac,
		integer	icp,
		integer	ipredu,
		integer, intent(in)	nel,
		integer, intent(in)	mtn,
		integer, intent(in)	ismstr,
		integer, intent(in)	jhbe )

Definition at line 30 of file mmat_h1.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-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: MTN
      INTEGER, INTENT(IN) :: ISMSTR
      INTEGER, INTENT(IN) :: JHBE
      INTEGER ICP,IPREDU
C     REAL
      my_real
     .   hh(2,*),hh1(2,*),fac(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,IFP
C     REAL
C-----------------------------------------------
       ipredu=0
       IF (icp==1.OR.(icp==2.AND.jhbe==17)) ipredu=1
C-----exclude nonlinear pressure laws (don't use p=c1*rho/(rho0-1)--hyperelastic first       
       IF (mtn==42.OR.mtn==62.OR.mtn==82) ipredu=0
       ifp = 0
       IF (icp==1.AND.jhbe==17.AND.ismstr==10) ifp=1
       IF (ifp==1) ipredu=1
C       
        IF (ipredu == 0) THEN
          DO i=1,nel
           hh1(1,i)= hh(1,i)
           hh1(2,i) =hh(2,i)
          ENDDO
        ELSEIF (icp==2.AND.jhbe==17.OR.ifp==1) THEN
C-----------constant pressure-----  
         DO i=1,nel
          hh1(1,i)=hh(2,i)*(one-fac(i))-two_third*hh(2,i)*fac(i)
          hh1(2,i) = hh(2,i)
         ENDDO
        ELSEIF (icp == 1) THEN
C-----------constant pressure-----  
         DO i=1,nel
          hh1(1,i)=-two_third*hh(2,i)
          hh1(2,i) = hh(2,i)
         ENDDO
        END IF !(IPREDU == 0) THEN
C
      RETURN