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 (mat_param, geo, off, for, eint, area, al0, pla, sti, mgm, ngl, eps, al, nel)

Function/Subroutine Documentation

m2lawt()

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

Definition at line 30 of file m2lawt.F.

!-----------------------------------------------
!   m o d u l e s
!-----------------------------------------------
      use matparam_def_mod
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 MGM(MVSIZ),NGL(MVSIZ)
      my_real
     .   geo(npropg,*),off(*),for(*),eint(*),
     .   area(*),al0(*),pla(*),sti(*),eps(mvsiz),al(mvsiz)
      type (matparam_struct_)  ,intent(in)    :: mat_param
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      INTEGER ICC(MVSIZ)
      my_real yma
      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)
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      my_real :: cvmgt
C-----------------------------------------------      
      !J-C formulation only
      
      DO i=1,nel
        ym(i)  = mat_param%young
        pr(i)  = mat_param%nu
        icc(i) = mat_param%iparam(2)
        ca(i)  = mat_param%uparam(1)
        cb(i)  = mat_param%uparam(2)
        cn(i)  = mat_param%uparam(3)
        epmx(i)= mat_param%uparam(4)
        ymax(i)= mat_param%uparam(5)
        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     ductile failure test
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