repla3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	repla3 (xk, dpx, tf, npf, iecrou, ifunc, ifv, epla, nel)

Function/Subroutine Documentation

repla3()

subroutine repla3	(		xk,
			dpx,
			tf,
		integer, dimension(*)	npf,
		integer, dimension(mvsiz)	iecrou,
		integer, dimension(mvsiz)	ifunc,
		integer, dimension(mvsiz)	ifv,
			epla,
		integer, intent(in)	nel )

Definition at line 35 of file repla3.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-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER NPF(*),IECROU(MVSIZ),IFV(MVSIZ)
C     REAL
      my_real
     .   xk(*), dpx(*), tf(*),epla(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IFUNC(MVSIZ), I, NPOINT
C     REAL
      my_real
     .    xki, fp
C-----------------------------------------------
      DO i=1,nel
        IF(ifunc(i)/=0.AND.iecrou(i)==3)THEN
          npoint=(npf(ifunc(i)+1)-npf(ifunc(i)))/2
C-------------------------------------
C       ELASTO PLASTIQUE (TENSION COMPRESSION COUPLEE 6 DOF COUPLES)
C-------------------------------------
          CALL interp(tf(npf(ifunc(i))),dpx(i),npoint,fp,xki)
          fp = fp * xk(i) / max( em20, (xk(i) - xki) )
          fp = max(abs(fp), em20)
          dpx(i) = dpx(i) + epla(i) / fp
        ENDIF
      ENDDO
C
      RETURN