qdefo2.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"
#include "com08_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	qdefo2 (v, w, y1, y2, y3, y4, z1, z2, z3, z4, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, py1, py2, pz1, pz2, wyz, dyz, dzy, eyy, ezz, ett, eyz, eyt, ezt, rx, ry, rz, sx, sy, sz, tx, ty, tz, voln, aire, airem, nc1, nc2, nc3, nc4, nel)

Function/Subroutine Documentation

qdefo2()

subroutine qdefo2	(		v,
			w,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4,
			vy1,
			vy2,
			vy3,
			vy4,
			vz1,
			vz2,
			vz3,
			vz4,
			py1,
			py2,
			pz1,
			pz2,
			wyz,
			dyz,
			dzy,
			eyy,
			ezz,
			ett,
			eyz,
			eyt,
			ezt,
			rx,
			ry,
			rz,
			sx,
			sy,
			sz,
			tx,
			ty,
			tz,
			voln,
			aire,
			airem,
		integer, dimension(*)	nc1,
		integer, dimension(*)	nc2,
		integer, dimension(*)	nc3,
		integer, dimension(*)	nc4,
		integer, intent(in)	nel )

Definition at line 29 of file qdefo2.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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      my_real
     .   v(3,*), w(3,*),
     .   y1(*),y2(*),y3(*),y4(*),z1(*),z2(*),z3(*),z4(*),
     .   vy1(*),vy2(*),vy3(*),vy4(*),vz1(*),vz2(*),vz3(*),vz4(*),
     .   py1(*), py2(*), pz1(*), pz2(*),
     .   wyz(*), dyz(*), dzy(*),
     .   eyy(*),ezz(*), ett(*), eyz(*), eyt(*), ezt(*),
     .   rx(*),ry(*), rz(*), sx(*), sy(*), sz(*), tx(*), ty(*), tz(*),
     .   voln(*),aire(*),airem(*)
C
      INTEGER NC1(*), NC2(*), NC3(*), NC4(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com08_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   vy13(mvsiz), vy24(mvsiz), vz13(mvsiz),vz24(mvsiz),
     .   ym1(mvsiz), ym2(mvsiz), ym3(mvsiz), ym4(mvsiz),
     .   zm1(mvsiz), zm2(mvsiz), zm3(mvsiz), zm4(mvsiz), dth,
     .   yavg(mvsiz),a1(mvsiz) , a2(mvsiz)
C-----------------------------------------------
      dth=half*dt1
      
      DO i=1,nel
        !from qforc2 subroutine : &W  = &V  
        !from bforc2 subroutine : &W /= &V  
        vy1(i)=w(2,nc1(i))
        vz1(i)=w(3,nc1(i))
        vy2(i)=w(2,nc2(i))
        vz2(i)=w(3,nc2(i))
        vy3(i)=w(2,nc3(i))
        vz3(i)=w(3,nc3(i))
        vy4(i)=w(2,nc4(i))
        vz4(i)=w(3,nc4(i))
      ENDDO
 
      DO i=1,nel
        ym1(i)=y1(i)-dth*vy1(i)
        ym2(i)=y2(i)-dth*vy2(i)
        ym3(i)=y3(i)-dth*vy3(i)
        ym4(i)=y4(i)-dth*vy4(i)
        zm1(i)=z1(i)-dth*vz1(i)
        zm2(i)=z2(i)-dth*vz2(i)
        zm3(i)=z3(i)-dth*vz3(i)
        zm4(i)=z4(i)-dth*vz4(i)
      ENDDO
 
      DO i=1,nel
        vy1(i)=v(2,nc1(i))
        vz1(i)=v(3,nc1(i))
        vy2(i)=v(2,nc2(i))
        vz2(i)=v(3,nc2(i))
        vy3(i)=v(2,nc3(i))
        vz3(i)=v(3,nc3(i))
        vy4(i)=v(2,nc4(i))
        vz4(i)=v(3,nc4(i))
      ENDDO
 
      DO i=1,nel
        py1(i)=half*(zm2(i)-zm4(i))
        py2(i)=half*(zm3(i)-zm1(i))
        pz1(i)=half*(ym4(i)-ym2(i))
        pz2(i)=half*(ym1(i)-ym3(i))      
      ENDDO
 
      DO i=1,nel
        a1(i) =ym2(i)*(zm3(i)-zm4(i))+ym3(i)*(zm4(i)-zm2(i))+ym4(i)*(zm2(i)-zm3(i))
        a2(i) =ym2(i)*(zm4(i)-zm1(i))+ym4(i)*(zm1(i)-zm2(i))+ym1(i)*(zm2(i)-zm4(i))
        airem(i)=(a1(i)+a2(i))*half
        !AUTRE FORMULE AIREM(I)=2.*(PZ2(I)*PY1(I)-PZ1(I)*PY2(I))
        yavg(i) =ym1(i)+ym2(i)+ym3(i)+ym4(i)
      ENDDO
 
      DO i=1,nel
        vy13(i)=vy1(i)-vy3(i)
        vy24(i)=vy2(i)-vy4(i)
        vz13(i)=vz1(i)-vz3(i)
        vz24(i)=vz2(i)-vz4(i)
       ENDDO
 
      DO i=1,nel
        IF(airem(i)>zero) THEN
          eyy(i)=(py1(i)*vy13(i)+py2(i)*vy24(i))/airem(i)
          ezz(i)=(pz1(i)*vz13(i)+pz2(i)*vz24(i))/airem(i)
          ett(i)=zero
          dzy(i)=(py1(i)*vz13(i)+py2(i)*vz24(i))/airem(i)
          dyz(i)=(pz1(i)*vy13(i)+pz2(i)*vy24(i))/airem(i)
          eyt(i)=zero
          ezt(i)=zero   
        ELSE
          eyy(i)=zero
          ezz(i)=zero
          ett(i)=zero
          dzy(i)=zero
          dyz(i)=zero
          eyt(i)=zero
          ezt(i)=zero
        END IF
      ENDDO
 
      IF(n2d==1) THEN
        DO i=1,nel
          IF (yavg(i)<=zero) cycle
          ett(i)=(vy1(i)+vy2(i)+vy3(i)+vy4(i))/yavg(i)
        ENDDO
      ENDIF
 
      DO i=1,nel
        eyz(i)= dzy(i)+dyz(i)
        wyz(i)=half*dt1*(dzy(i)-dyz(i))
      ENDDO
C---------------------------------------
C     DERIVEES POUR LE CALCUL DES FORCES
C---------------------------------------
      DO i=1,nel
        py1(i)=half*(z2(i)-z4(i))
        py2(i)=half*(z3(i)-z1(i))
        pz1(i)=half*(y4(i)-y2(i))
        pz2(i)=half*(y1(i)-y3(i))      
      ENDDO
C---------------------------------------
C     REPERE LOCAL POUR MATERIAUX NON ISOTROPES
C---------------------------------------
      DO i=1,nel
        rx(i)=one
        ry(i)=zero
        rz(i)=zero
        sx(i)=zero    
        sy(i)=-pz1(i)-pz2(i)
        sz(i)=+py1(i)+py2(i)
        tx(i)=zero
        ty(i)=+pz1(i)-pz2(i)
        tz(i)=-py1(i)+py2(i)
      ENDDO
 
      RETURN