m1lawt.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "com08_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	m1lawt (pm, geo, off, for, eint, area, al0, sti, mat, mgm, eps, al, nel)

Function/Subroutine Documentation

m1lawt()

subroutine m1lawt	(		pm,
			geo,
			off,
			for,
			eint,
			area,
			al0,
			sti,
		integer, dimension(mvsiz)	mat,
		integer, dimension(mvsiz)	mgm,
			eps,
			al,
		integer, intent(in)	nel )

Definition at line 28 of file m1lawt.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
#include      "com08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER MAT(MVSIZ),MGM(MVSIZ)
      my_real
     .   pm(npropm,*),geo(npropg,*),off(*),for(*),eint(*),
     .   area(*),al0(*),sti(*),eps(mvsiz),al(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J
      my_real
     .   pr(mvsiz),ym(mvsiz),gap(mvsiz),yma
C-----------------------------------------------      
      DO i=1,nel
        ym(i)  = pm(20,mat(i))
        pr(i)  = pm(21,mat(i))
        gap(i) = geo(2,mgm(i))
      ENDDO
c
      DO i=1,nel
        IF(gap(i) > zero .AND. al(i) <= (al0(i)-gap(i))) off(i)=one
      ENDDO
c
      DO i=1,nel
        eint(i)=eint(i)+for(i)*eps(i)*al(i)*dt1*half
      ENDDO
c
      DO i=1,nel
        area(i)=area(i)*(one - two*pr(i)*eps(i)*dt1*off(i))
      ENDDO
c
      DO i=1,nel
        yma    = ym(i)*area(i)
        for(i) = for(i)+yma*eps(i)*dt1
        sti(i) = zero
        IF (al(i) /= zero)
     .  sti(i) = yma / al(i)
        sti(i) = sti(i)*off(i)
        for(i) = for(i)*off(i)
      ENDDO
c
      DO i=1,nel
        eint(i)=eint(i)+for(i)*eps(i)*al(i)*dt1*half
      ENDDO
C-----------------------------------------------
      RETURN