thermexpp.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	thermexppi (elbuf_str, nel, npt, mat, pid, pm, geo, al, eth, off, for, eint)
subroutine	thermexppg (nel, mat, pid, pm, geo, off, eth, for, eint)

Function/Subroutine Documentation

thermexppg()

subroutine thermexppg	(	integer	nel,
		integer, dimension(mvsiz)	mat,
		integer, dimension(mvsiz)	pid,
			pm,
			geo,
			off,
			eth,
			for,
			eint )

Definition at line 123 of file thermexpp.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"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL,MAT(MVSIZ),PID(MVSIZ)
C     REAL
      my_real
     .   pm(npropm,*), for(nel,3),  eint(nel,2), geo(npropg,*),
     .   off(*),al(mvsiz),eth(nel)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J
C     REAL
      my_real
     .   ym,a1
C-----------------------------------------------
      DO i=1,nel
        ym  =pm(20,mat(i))
        a1  =geo(1,pid(i))
        for(i,1) = (for(i,1)  - eth(i)*a1*ym )*off(i)   
!!        EINT(I,1) = EINT(I,1) - FOR(I,1)*ETH(I)*AL(I)*HALF*OFF(I)
      ENDDO             
C-----------------------------------------------
      RETURN

thermexppi()

subroutine thermexppi	(	type (elbuf_struct_), target	elbuf_str,
		integer	nel,
		integer	npt,
		integer, dimension(mvsiz)	mat,
		integer, dimension(mvsiz)	pid,
			pm,
			geo,
			al,
			eth,
			off,
			for,
			eint )

Definition at line 30 of file thermexpp.F.

C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_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"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NPT,NEL
      INTEGER MAT(MVSIZ),PID(MVSIZ)
C     REAL
      my_real
     .   pm(npropm,*), geo(npropg,*),for(nel,3), eint(nel,2),
     .   off(*), al(*),eth(nel)
C
      TYPE (ELBUF_STRUCT_), TARGET :: ELBUF_STR
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER ICC(MVSIZ),IRTY(MVSIZ),INDX(MVSIZ)
      INTEGER I, J, J1, J2, NPIF, MX, IPT, IPY,IPZ,IPA,NINDX,
     .        IR,IS,ILAYER,II(3)
C     REAL
      my_real
     .  e, signxx,dfxx,apt
C
      TYPE(L_BUFEL_),POINTER :: LBUF  
      TYPE(BUF_LAY_) ,POINTER :: BUFLY
C=======================================================================
!   
   
C-------------------------------------
C     DEBUT DE BOUCLE SUR POINTS INTEGRATION
C--------------------------------------
C             
      ipa  = 400  
      DO i=1,3
        ii(i) = nel*(i-1)
      ENDDO 
      DO ipt= 1,npt
        ilayer=1
        ir = 1
        is = 1
        lbuf => elbuf_str%BUFLY(ilayer)%LBUF(ir,is,ipt)
        bufly => elbuf_str%BUFLY(ilayer)
c---    Total strain
C 
        IF (bufly%L_STRA > 0) THEN
          DO i= 1,nel
            lbuf%STRA(ii(1)+i) = lbuf%STRA(ii(1)+i) + eth(i)
          ENDDO
        ENDIF        
C
C---    Contraintes elastiques
C
        DO i =1 ,nel             
          mx = mat(i)                            
          e    = pm(20,mx)                         
          signxx =  -e*eth(i)  
C       FORCES ET MOMENTS
          apt = geo(ipa+ipt,pid(i)) 
          dfxx = apt*signxx
          for(i,1) = for(i,1) + dfxx     
C                                                    
          lbuf%SIG(ii(1)+i) = lbuf%SIG(ii(1)+i) + signxx
        ENDDO                                                       
C-------------------------------------
C       FIN DE BOUCLE SUR POINT INTEGRATION
C-------------------------------------
       ENDDO   
C-----------
      RETURN