m2lawt.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "com08_c.inc"
#include "units_c.inc"
#include "scr17_c.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

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

Function/Subroutine Documentation

m2lawt()

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

Definition at line 28 of file m2lawt.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.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"
#include      "units_c.inc"
#include      "scr17_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),NGL(MVSIZ)
      my_real
     .   pm(npropm,*),geo(npropg,*),off(*),for(*),eint(*),
     .   area(*),al0(*),pla(*),sti(*),eps(mvsiz),al(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER ICC(MVSIZ)
      INTEGER I,J
      my_real
     .   ym(mvsiz),gap(mvsiz),
     .   ca(mvsiz),cb(mvsiz),cn(mvsiz),ymax(mvsiz),
     .   epmx(mvsiz),yld(mvsiz),ff(mvsiz),hh(mvsiz),
     .   aa(mvsiz),pr(mvsiz),yma
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      my_real
     .   cvmgt
C-----------------------------------------------      
      DO i=1,nel
        ym(i)  =pm(20,mat(i))
        pr(i)  =pm(21,mat(i))
        ca(i)  =pm(38,mat(i))
        cb(i)  =pm(39,mat(i))
        cn(i)  =pm(40,mat(i))
        epmx(i)=pm(41,mat(i))
        ymax(i)=pm(42,mat(i))
        icc(i) =nint(pm(49,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)=yma / al(i)
      ENDDO
c
      DO i=1,nel
        yld(i)=ca(i)+cb(i)*pla(i)**cn(i)
      ENDDO
C-----------------------
C     MODULE ECROUISSAGE
C-----------------------
      DO i=1,nel
        IF (cn(i) == one) THEN
          hh(i)= cb(i)
        ELSE
          IF (pla(i) /= zero) THEN
            hh(i)= cb(i)*cn(i)/pla(i)**(one -cn(i))
          ELSE
            hh(i)=zero
          ENDIF
        ENDIF
      ENDDO
c
      DO i=1,nel
        aa(i)  = (ym(i)+hh(i))*area(i)
        yld(i) = min(yld(i),ymax(i))
        ff(i)  = abs(for(i))-yld(i)*area(i)
        ff(i)  = max(zero,ff(i))
      ENDDO
c
      DO i=1,nel
        pla(i)=pla(i)+ff(i)/aa(i)
      ENDDO
c
      DO i=1,nel
        for(i)=cvmgt(sign(yld(i)*area(i),for(i)),for(i),ff(i) > zero)
      ENDDO
C--------------------------------
C     TEST DE RUPTURE DUCTILE
C-------------------------------
      DO i=1,nel
        IF (off(i) < em01) off(i)=zero
        IF (off(i) < one) off(i)=off(i)*four_over_5
      ENDDO 
c
      DO i=1,nel
        IF (off(i) < one)      cycle
        IF (pla(i) < epmx(i)) cycle
        off(i)=off(i)*four_over_5
        idel7nok = 1
      ENDDO
c
      DO i=1,nel
        IF (off(i) /= four_over_5) cycle
#include "lockon.inc"
        WRITE(iout,1000)  ngl(i)
        WRITE(istdo,1100) ngl(i),tt
#include "lockoff.inc"
      ENDDO
c
      DO i=1,nel
        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
 1000 FORMAT(1x,'-- RUPTURE OF TRUSS ELEMENT NUMBER ',i10)
 1100 FORMAT(1x,'-- RUPTURE OF TRUSS ELEMENT :',i10,' AT TIME :',g11.4)
C-----------------------------------------------
      RETURN