scdefo3.F File Reference

#include "implicit_f.inc"
#include "com08_c.inc"
#include "impl1_c.inc"
#include "scr05_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	scdefo3 (dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, d4, d5, d6, dcxx, dcxy, dcxz, dcyx, dcyy, dcyz, dczx, dczy, dczz, dc4, dc5, dc6, dhxx, dhxy, dhxz, dhyx, dhyy, dhyz, dhzx, dhzy, dhzz, dh4, dh5, dh6, zi, wi, vzl, vol, volg, volo, ddhv, sig, sigzm, volm, usb, eint, off, offg, dti, offs, dvc, vol0dp, voldp, ipres, nel)

Function/Subroutine Documentation

scdefo3()

subroutine scdefo3	(		dxx,
			dxy,
			dxz,
			dyx,
			dyy,
			dyz,
			dzx,
			dzy,
			dzz,
			d4,
			d5,
			d6,
			dcxx,
			dcxy,
			dcxz,
			dcyx,
			dcyy,
			dcyz,
			dczx,
			dczy,
			dczz,
			dc4,
			dc5,
			dc6,
			dhxx,
			dhxy,
			dhxz,
			dhyx,
			dhyy,
			dhyz,
			dhzx,
			dhzy,
			dhzz,
			dh4,
			dh5,
			dh6,
			zi,
			wi,
			vzl,
			vol,
			volg,
			volo,
			ddhv,
			sig,
			sigzm,
			volm,
			usb,
			eint,
			off,
			offg,
			dti,
			offs,
			dvc,
		double precision, dimension(*)	vol0dp,
		double precision, dimension(*)	voldp,
		integer, intent(in)	ipres,
		integer, intent(in)	nel )

Definition at line 29 of file scdefo3.F.

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      "com08_c.inc"
#include      "impl1_c.inc"
#include      "scr05_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL,IPRES
C     REAL
      my_real
     .   dxx(*), dxy(*), dxz(*), dyx(*), dyy(*),dyz(*), 
     .   dzx(*), dzy(*), dzz(*), d4(*), d5(*), d6(*),
     .   dcxx(*), dcxy(*),dcxz(*), dcyx(*), dcyy(*),dcyz(*), 
     .   dczx(*), dczy(*), dczz(*), dc4(*), dc5(*), dc6(*),
     .   dhxx(*), dhxy(*),dhxz(*), dhyx(*), dhyy(*),dhyz(*), 
     .   dhzx(*), dhzy(*), dhzz(*), 
     .   dh4(*), dh5(*), dh6(*), ddhv(*),
     .   zi,wi,vzl(*),vol(*),volg(*),volo(*),
     .   sig(nel,6),sigzm(*),volm(*),usb(*),eint(*),off(*),offg(*),dti,
     .   offs(*),dvc(*)
      DOUBLE PRECISION 
     .   VOL0DP(*),VOLDP(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  I
      my_real
     .  dv,dvz,ddv,tol,dt1d2,dt1d1
C=======================================================================
      tol=one-em20
      DO i=1,nel
        voldp(i)=half*wi*(volg(i)+vzl(i)*zi)
        vol(i)=voldp(i)
        off(i)=offg(i)
        IF (vol(i)<=zero) THEN
           vol(i)= em20
           off(i)= zero
        ELSEIF(off(i)==zero.OR.offs(i)==two.OR.ismdisp>0)THEN
           voldp(i)=max(em20,voldp(i))
           vol(i)=max(em20,vol(i))
        ENDIF
      ENDDO
      IF (ipres==1) THEN
        DO i=1,nel
C  ------due to the assumed strain terms----
          dv = ddhv(i)*dt1*zi
C  ------due to the sigzz traitement----
          dvz = -(sig(i,3)-sigzm(i)/max(em20,volm(i)))*usb(i)
          ddv = (dv+dvz)*off(i)
          dvc(i) = ddv
          IF (ddv>tol) THEN
           vol(i)= em20
           off(i)= zero
           ddv = zero
          ELSEIF(off(i)==zero.OR.offs(i)==two.OR.ismdisp>0)THEN
           ddv = zero
          END IF
C         
          IF (iresp==1) vol0dp(i) = vol0dp(i)*(one- ddv)
          IF(ismdisp>0.OR.offs(i)==two) dvc(i) = dvz
C      
          volo(i) = volo(i)*(one- ddv)
          eint(i) = eint(i)/(one- ddv)
        ENDDO
      ELSE
        DO i=1,nel
          dvc(i) = zero
        ENDDO
      END IF
C   +++ -----partie non constante------
      DO i=1,nel
        dxx(i)=dcxx(i)+zi*dhxx(i)
        dyy(i)=dcyy(i)+zi*dhyy(i)
        dzz(i)=dczz(i)+zi*dhzz(i)
        dxy(i)=dcxy(i)+zi*dhxy(i)
        dyx(i)=dcyx(i)+zi*dhyx(i)
        dzx(i)=dczx(i)+zi*dhzx(i)
        dzy(i)=dczy(i)+zi*dhzy(i)
        dxz(i)=dcxz(i)+zi*dhxz(i)
        dyz(i)=dcyz(i)+zi*dhyz(i)
c        D4(I)=DC4(I)+ZI*DH4(I)
c        D5(I)=DC5(I)+ZI*DH5(I)
c        D6(I)=DC6(I)+ZI*DH6(I)
      ENDDO
C-----------
      dt1d1=dt1
      IF (ismdisp>0.AND.iscau==0) dt1d1=zero
      dt1d2=half*dt1d1
       DO i=1,nel
        d4(i)   = dxy(i)+dyx(i)
     .           -dt1d1*(dxx(i)*dxy(i)+dyx(i)*dyy(i)+dzx(i)*dzy(i))
        d5(i)   = dyz(i)+dzy(i)
     .           -dt1d1*(dyy(i)*dyz(i)+dzy(i)*dzz(i)+dxy(i)*dxz(i))
        d6(i)   = dxz(i)+dzx(i)
     .           -dt1d1*(dzz(i)*dzx(i)+dxz(i)*dxx(i)+dyz(i)*dyx(i))
        dxx(i)  = dxx(i)
     .           -dt1d2*(dxx(i)*dxx(i)+dyx(i)*dyx(i)+dzx(i)*dzx(i))
        dyy(i)  = dyy(i)
     .           -dt1d2*(dyy(i)*dyy(i)+dzy(i)*dzy(i)+dxy(i)*dxy(i))
        dzz(i)  = dzz(i)
     .           -dt1d2*(dzz(i)*dzz(i)+dxz(i)*dxz(i)+dyz(i)*dyz(i))
       ENDDO
C       
      RETURN