sfor_visn4.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sfor_visn4 (vc, fld, tol_v, ifc1, vx1, vx2, vx3, vx4, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, for_t1, for_t2, for_t3, for_t4, ifctl, stif, mu, nel, e_distor, dt1)

Function/Subroutine Documentation

sfor_visn4()

subroutine sfor_visn4	(	intent(in)	vc,
		intent(in)	fld,
		intent(in)	tol_v,
		integer, dimension(mvsiz), intent(inout)	ifc1,
		intent(in)	vx1,
		intent(in)	vx2,
		intent(in)	vx3,
		intent(in)	vx4,
		intent(in)	vy1,
		intent(in)	vy2,
		intent(in)	vy3,
		intent(in)	vy4,
		intent(in)	vz1,
		intent(in)	vz2,
		intent(in)	vz3,
		intent(in)	vz4,
		intent(inout)	for_t1,
		intent(inout)	for_t2,
		intent(inout)	for_t3,
		intent(inout)	for_t4,
		integer, intent(out)	ifctl,
		intent(inout)	stif,
		intent(in)	mu,
		integer, intent(in)	nel,
		intent(inout)	e_distor,
		intent(in)	dt1 )

Definition at line 28 of file sfor_visn4.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
      INTEGER, INTENT (OUT) :: IFCTL
      INTEGER, DIMENSION(MVSIZ),INTENT (INOUT) :: IFC1
      my_real, DIMENSION(MVSIZ), INTENT (INOUT) :: stif
      my_real, DIMENSION(MVSIZ), INTENT (IN) :: fld ,
     4                vx1,     vx2,     vx3,     vx4,
     6                vy1,     vy2,     vy3,     vy4,
     8                vz1,     vz2,     vz3,     vz4
      my_real, DIMENSION(MVSIZ,3), INTENT (IN)  :: vc
      my_real, DIMENSION(MVSIZ,3), INTENT (INOUT) :: 
     .                     for_t1, for_t2, for_t3, for_t4
      my_real, INTENT (IN) ::  tol_v,mu,dt1
      my_real, DIMENSION(NEL), INTENT (INOUT)  :: e_distor
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C                                                                     12
      my_real
     .   fac,vnj(4),vl,tol_v2,
     .   v2max,vc2,vref(mvsiz,3),vr2
C----------------------------
        tol_v2 = tol_v*tol_v
        ifctl = 0
         DO i=1,nel
           vc2 = vc(i,1)*vc(i,1)+vc(i,2)*vc(i,2)+vc(i,3)*vc(i,3)
           IF (vc2 <em20.OR.stif(i)==zero) cycle
           vl = tol_v2*vc2
           vnj(1) = vx1(i)*vx1(i) + vy1(i)*vy1(i) + vz1(i)*vz1(i)
           vnj(2) = vx2(i)*vx2(i) + vy2(i)*vy2(i) + vz2(i)*vz2(i)
           vnj(3) = vx3(i)*vx3(i) + vy3(i)*vy3(i) + vz3(i)*vz3(i)
           vnj(4) = vx4(i)*vx4(i) + vy4(i)*vy4(i) + vz4(i)*vz4(i)
           v2max = max(vnj(1),vnj(2),vnj(3),vnj(4))
           vr2 = min(vnj(1),vnj(2),vnj(3),vnj(4))
           vl = tol_v2*vr2
           IF (vr2 == vnj(1)) THEN 
              vref(i,1) = vx1(i)
              vref(i,2) = vy1(i)
              vref(i,3) = vz1(i) 
            ELSE IF (vr2 == vnj(2)) THEN
              vref(i,1) = vx2(i)
              vref(i,2) = vy2(i)
              vref(i,3) = vz2(i) 
            ELSE IF (vr2 == vnj(3)) THEN
              vref(i,1) = vx3(i)
              vref(i,2) = vy3(i)
              vref(i,3) = vz3(i) 
            ELSE 
              vref(i,1) = vx4(i)
              vref(i,2) = vy4(i)
              vref(i,3) = vz4(i)
            END IF
           IF (v2max > vl) ifc1(i) = 1
           IF (ifc1(i) > 0)   ifctl=1
         END DO
C
       IF (ifctl==1) THEN
         fac = one + two*mu
         DO i=1,nel
           IF (ifc1(i)==0) cycle
           for_t1(i,1) = for_t1(i,1) - fld(i)*(vx1(i)-vref(i,1))
           for_t1(i,2) = for_t1(i,2) - fld(i)*(vy1(i)-vref(i,2))
           for_t1(i,3) = for_t1(i,3) - fld(i)*(vz1(i)-vref(i,3))
           for_t2(i,1) = for_t2(i,1) - fld(i)*(vx2(i)-vref(i,1))
           for_t2(i,2) = for_t2(i,2) - fld(i)*(vy2(i)-vref(i,2))
           for_t2(i,3) = for_t2(i,3) - fld(i)*(vz2(i)-vref(i,3))
           for_t3(i,1) = for_t3(i,1) - fld(i)*(vx3(i)-vref(i,1))
           for_t3(i,2) = for_t3(i,2) - fld(i)*(vy3(i)-vref(i,2))
           for_t3(i,3) = for_t3(i,3) - fld(i)*(vz3(i)-vref(i,3))
           for_t4(i,1) = for_t4(i,1) - fld(i)*(vx4(i)-vref(i,1))
           for_t4(i,2) = for_t4(i,2) - fld(i)*(vy4(i)-vref(i,2))
           for_t4(i,3) = for_t4(i,3) - fld(i)*(vz4(i)-vref(i,3))
           stif(i)      = fac*stif(i)
           e_distor(i)=e_distor(i)- dt1*(for_t1(i,1)*(vx1(i)-vref(i,1))+
     .                                   for_t1(i,2)*(vy1(i)-vref(i,2))+
     .                                   for_t1(i,3)*(vz1(i)-vref(i,3))+
     .                                   for_t2(i,1)*(vx2(i)-vref(i,1))+
     .                                   for_t2(i,2)*(vy2(i)-vref(i,2))+
     .                                   for_t2(i,3)*(vz2(i)-vref(i,3))+
     .                                   for_t3(i,1)*(vx3(i)-vref(i,1))+
     .                                   for_t3(i,2)*(vy3(i)-vref(i,2))+
     .                                   for_t3(i,3)*(vz3(i)-vref(i,3))+
     .                                   for_t4(i,1)*(vx4(i)-vref(i,1))+
     .                                   for_t4(i,2)*(vy4(i)-vref(i,2))+
     .                                   for_t4(i,3)*(vz4(i)-vref(i,3)))
         ENDDO
       END IF 
C
      RETURN