sfor_visn8.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	sfor_visn8 (vc, fld, tol_v, ifc1, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, for_t1, for_t2, for_t3, for_t4, for_t5, for_t6, for_t7, for_t8, ifctl, stif, mu, nel, e_distor, dt1)

Function/Subroutine Documentation

sfor_visn8()

subroutine sfor_visn8	(	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)	vx5,
		intent(in)	vx6,
		intent(in)	vx7,
		intent(in)	vx8,
		intent(in)	vy1,
		intent(in)	vy2,
		intent(in)	vy3,
		intent(in)	vy4,
		intent(in)	vy5,
		intent(in)	vy6,
		intent(in)	vy7,
		intent(in)	vy8,
		intent(in)	vz1,
		intent(in)	vz2,
		intent(in)	vz3,
		intent(in)	vz4,
		intent(in)	vz5,
		intent(in)	vz6,
		intent(in)	vz7,
		intent(in)	vz8,
		intent(inout)	for_t1,
		intent(inout)	for_t2,
		intent(inout)	for_t3,
		intent(inout)	for_t4,
		intent(inout)	for_t5,
		intent(inout)	for_t6,
		intent(inout)	for_t7,
		intent(inout)	for_t8,
		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_visn8.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,
     5                vx5,     vx6,     vx7,     vx8,
     6                vy1,     vy2,     vy3,     vy4,
     7                vy5,     vy6,     vy7,     vy8,
     8                vz1,     vz2,     vz3,     vz4,
     9                vz5,     vz6,     vz7,     vz8
      my_real, DIMENSION(MVSIZ,3), INTENT (IN)  :: vc
      my_real, DIMENSION(MVSIZ,3), INTENT (INOUT) :: 
     .                     for_t1, for_t2, for_t3, for_t4,
     .                     for_t5, for_t6, for_t7, for_t8
      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,J,IFCONT
C                                                                     12
      my_real
     .   vre(6),vxc,vyc,vzc,fx,fy,fz,fac,vnj(8),vl,tol_v2,
     .   vxmax,vymax,vzmax,v2max,vc2,lam_min,lam_min2,lam_max
C----------------------------
        tol_v2 = tol_v*tol_v
        ifctl = 0
        ifcont = 1 ! no more for contact condition due to P/on issue
         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)
           vnj(5) = vx5(i)*vx5(i) + vy5(i)*vy5(i) + vz5(i)*vz5(i)
           vnj(6) = vx6(i)*vx6(i) + vy6(i)*vy6(i) + vz6(i)*vz6(i)
           vnj(7) = vx7(i)*vx7(i) + vy7(i)*vy7(i) + vz7(i)*vz7(i)
           vnj(8) = vx8(i)*vx8(i) + vy8(i)*vy8(i) + vz8(i)*vz8(i)
           v2max = max(vnj(1),vnj(2),vnj(3),vnj(4))
           v2max = max(v2max,vnj(5),vnj(6),vnj(7),vnj(8))
           IF (v2max > vl) ifc1(i) = 1
           IF (ifc1(i) > 0)   ifctl=1
!           IF (V2MAX > FOUR*VC2) IFCONT=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)-vc(i,1))
           for_t1(i,2) = for_t1(i,2) - fld(i)*(vy1(i)-vc(i,2))
           for_t1(i,3) = for_t1(i,3) - fld(i)*(vz1(i)-vc(i,3))
           for_t2(i,1) = for_t2(i,1) - fld(i)*(vx2(i)-vc(i,1))
           for_t2(i,2) = for_t2(i,2) - fld(i)*(vy2(i)-vc(i,2))
           for_t2(i,3) = for_t2(i,3) - fld(i)*(vz2(i)-vc(i,3))
           for_t3(i,1) = for_t3(i,1) - fld(i)*(vx3(i)-vc(i,1))
           for_t3(i,2) = for_t3(i,2) - fld(i)*(vy3(i)-vc(i,2))
           for_t3(i,3) = for_t3(i,3) - fld(i)*(vz3(i)-vc(i,3))
           for_t4(i,1) = for_t4(i,1) - fld(i)*(vx4(i)-vc(i,1))
           for_t4(i,2) = for_t4(i,2) - fld(i)*(vy4(i)-vc(i,2))
           for_t4(i,3) = for_t4(i,3) - fld(i)*(vz4(i)-vc(i,3))
           for_t5(i,1) = for_t5(i,1) - fld(i)*(vx5(i)-vc(i,1))
           for_t5(i,2) = for_t5(i,2) - fld(i)*(vy5(i)-vc(i,2))
           for_t5(i,3) = for_t5(i,3) - fld(i)*(vz5(i)-vc(i,3))
           for_t6(i,1) = for_t6(i,1) - fld(i)*(vx6(i)-vc(i,1))
           for_t6(i,2) = for_t6(i,2) - fld(i)*(vy6(i)-vc(i,2))
           for_t6(i,3) = for_t6(i,3) - fld(i)*(vz6(i)-vc(i,3))
           for_t7(i,1) = for_t7(i,1) - fld(i)*(vx7(i)-vc(i,1))
           for_t7(i,2) = for_t7(i,2) - fld(i)*(vy7(i)-vc(i,2))
           for_t7(i,3) = for_t7(i,3) - fld(i)*(vz7(i)-vc(i,3))
           for_t8(i,1) = for_t8(i,1) - fld(i)*(vx8(i)-vc(i,1))
           for_t8(i,2) = for_t8(i,2) - fld(i)*(vy8(i)-vc(i,2))
           for_t8(i,3) = for_t8(i,3) - fld(i)*(vz8(i)-vc(i,3))
           stif(i)      = fac*stif(i)
           e_distor(i)=e_distor(i) -dt1*(for_t1(i,1)*(vx1(i)-vc(i,1))+
     .                                   for_t1(i,2)*(vy1(i)-vc(i,2))+
     .                                   for_t1(i,3)*(vz1(i)-vc(i,3))+
     .                                   for_t2(i,1)*(vx2(i)-vc(i,1))+
     .                                   for_t2(i,2)*(vy2(i)-vc(i,2))+
     .                                   for_t2(i,3)*(vz2(i)-vc(i,3))+
     .                                   for_t3(i,1)*(vx3(i)-vc(i,1))+
     .                                   for_t3(i,2)*(vy3(i)-vc(i,2))+
     .                                   for_t3(i,3)*(vz3(i)-vc(i,3))+
     .                                   for_t4(i,1)*(vx4(i)-vc(i,1))+
     .                                   for_t4(i,2)*(vy4(i)-vc(i,2))+
     .                                   for_t4(i,3)*(vz4(i)-vc(i,3))+
     .                                   for_t5(i,1)*(vx5(i)-vc(i,1))+
     .                                   for_t5(i,2)*(vy5(i)-vc(i,2))+
     .                                   for_t5(i,3)*(vz5(i)-vc(i,3))+
     .                                   for_t6(i,1)*(vx6(i)-vc(i,1))+
     .                                   for_t6(i,2)*(vy6(i)-vc(i,2))+
     .                                   for_t6(i,3)*(vz6(i)-vc(i,3))+
     .                                   for_t7(i,1)*(vx7(i)-vc(i,1))+
     .                                   for_t7(i,2)*(vy7(i)-vc(i,2))+
     .                                   for_t7(i,3)*(vz7(i)-vc(i,3))+
     .                                   for_t8(i,1)*(vx8(i)-vc(i,1))+
     .                                   for_t8(i,2)*(vy8(i)-vc(i,2))+
     .                                   for_t8(i,3)*(vz8(i)-vc(i,3)))
         ENDDO
       END IF 
!
        ifctl = ifcont  ! used for self-contact compute less strict : 2 times
C
      RETURN