cnloc_matini.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	cnloc_matini (elbuf_str, nel, ipm, bufmat, time, varnl, failure)

Function/Subroutine Documentation

cnloc_matini()

subroutine cnloc_matini	(	type(elbuf_struct_), target	elbuf_str,
		integer	nel,
		integer, dimension(npropmi,*)	ipm,
		target	bufmat,
			time,
			varnl,
		logical	failure )

Definition at line 35 of file cnloc_matini.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_mod      
      USE message_mod
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-----------------------------------------------
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL
      TYPE(ELBUF_STRUCT_), TARGET :: ELBUF_STR
      my_real
     .        varnl(nel,*),bufmat(*),time
      TARGET :: bufmat
      INTEGER IPM(NPROPMI,*)
      LOGICAL :: FAILURE
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,K,IR,IS,IT,IADBUF,NPTT,NPTR,NPTS,
     .        IJ(5),IJ1,IJ2,IJ3,IJ4,IJ5,IPG,
     .        ILAW,IMAT,NUVAR,NUPARAM
      TYPE(BUF_LAY_), POINTER :: BUFLY
      TYPE(L_BUFEL_), POINTER :: LBUF     
      TYPE(BUF_MAT_), POINTER :: MBUF
      TYPE(G_BUFEL_), POINTER :: GBUF
      my_real, 
     .  DIMENSION(:), POINTER :: uparam,uvar,thkn
C=======================================================================
      ! Buffer of the first layer
      bufly => elbuf_str%BUFLY(1)
      ! General buffer of the elements of the group
      gbuf  => elbuf_str%GBUF
      ! Global thickness of the shell
      thkn  => gbuf%THK(1:nel)
C
      ! Number of integration points in the thickness
      nptr  = elbuf_str%NPTR
      npts  = elbuf_str%NPTS
      nptt  = bufly%NPTT
      ilaw  = bufly%ILAW
      imat  = bufly%IMAT
      nuvar = bufly%NVAR_MAT
C
      ! Material properties and parameters
      iadbuf  = max(1,ipm(7,imat))
      nuparam = ipm(9,imat) 
      uparam  => bufmat(iadbuf:iadbuf+nuparam)       
C
      DO k=1,5
        ij(k) = nel*(k-1)
      ENDDO
      ij1 = ij(1) + 1
      ij2 = ij(2) + 1
      ij3 = ij(3) + 1
      ij4 = ij(4) + 1
      ij5 = ij(5) + 1
C
      ! Loop over integration points in the plane
      DO ir = 1,nptr
        DO is = 1,npts
          ! Current Gauss point
          ipg  = (is-1)*elbuf_str%NPTR + ir
          ! Loop over thickness integration points
          DO it = 1,nptt       
            ! Integration point buffer
            lbuf => bufly%LBUF(ir,is,it)
            ! Interval variable buffer
            uvar => bufly%MAT(ir,is,it)%VAR
            ! Non-local plastic strain update
            IF (gbuf%G_PLANL > 0) THEN 
              DO i = 1,nel 
                varnl(i,it)   = max(varnl(i,it),zero)
                lbuf%PLANL(i) = lbuf%PLANL(i) + varnl(i,it)
              ENDDO 
            ENDIF
            ! Select the material initialization routine
            SELECT CASE (ilaw)
              CASE(104)
                CALL cnloc_mat104_ini(nel      ,ipg      ,it       ,nuparam  ,nuvar    ,uparam   ,
     .                                uvar     ,lbuf%PLA ,lbuf%OFF ,lbuf%THK ,gbuf%OFF ,lbuf%SIG(ij1),
     .                                lbuf%SIG(ij2),lbuf%SIG(ij3)  ,lbuf%SIG(ij4)      ,lbuf%SIG(ij5),
     .                                thkn     ,lbuf%DMG ,nptr     ,npts     ,nptt     ,bufly        ,
     .                                time     ,varnl(1,it),failure)                  
            END SELECT
            ! Plastic strain to regularize
            IF (gbuf%G_PLANL > 0) THEN 
              DO i = 1,nel 
                varnl(i,it) = lbuf%PLA(i)
              ENDDO
            ENDIF
          ENDDO
        ENDDO
      ENDDO
C-----------