cbavisc.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	cbavisc (jft, jlt, vdef, amu, off, shf, nu, rho, ssp, area, thk, for, mom, npt, mtn, ipartc, evis, dt1, nel)
subroutine	cbavisnp1 (jft, jlt, vxyz, rxyz, vcore, amu, off, rho, ssp, area, thk, g, dt1, vf, ipartc, evis, kfts)

Function/Subroutine Documentation

cbavisc()

subroutine cbavisc	(	integer	jft,
		integer	jlt,
			vdef,
			amu,
			off,
			shf,
			nu,
			rho,
			ssp,
			area,
			thk,
			for,
			mom,
		integer	npt,
		integer	mtn,
		integer, dimension(*)	ipartc,
			evis,
			dt1,
		integer	nel )

Definition at line 30 of file cbavisc.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "param_c.inc"
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT, JLT ,NPT,MTN,MX,IPARTC(*),NEL
      my_real
     .   vdef(mvsiz,8),off(*),shf(*),
     .   for(nel,5),mom(nel,3),amu(*),rho(*),ssp(*),
     .   area(*),thk(*),nu(*),evis(npsav,*),dt1
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J
      my_real
     .   dtinv,visc(mvsiz),gg,gs,fac,fx,fy,fxy,dv(mvsiz)
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
C
        IF (npt == 1.AND.mtn /= 1) THEN
          fac =zero
        ELSE
          fac =one
        ENDIF
        DO i=jft,jlt
          visc(i) =onep414*off(i)*amu(i)*rho(i)*ssp(i)*sqrt(area(i))
          gg =half/(one+nu(i))
          fx =visc(i)*(vdef(i,1)+nu(i)*vdef(i,2))
          fy =visc(i)*(vdef(i,2)+nu(i)*vdef(i,1))
          fxy =visc(i)*vdef(i,3) *gg
          for(i,1)=for(i,1)+fx
          for(i,2)=for(i,2)+fy
          for(i,3)=for(i,3)+fxy
          dv(i) =area(i)*thk(i)*dt1
          mx = ipartc(i)
          evis(8,mx)=evis(8,mx)+(fx*vdef(i,1)+fy*vdef(i,2))*dv(i)
        ENDDO
        IF (npt /= 1) THEN
         DO i=jft,jlt
          visc(i)    =zep3*thk(i)*visc(i)
          gg =half/(one+nu(i))
          fx =visc(i)*(vdef(i,6)+nu(i)*vdef(i,7))
          fy =visc(i)*(vdef(i,7)+nu(i)*vdef(i,6))
          fxy =visc(i)* vdef(i,8) *gg
          mom(i,1)=mom(i,1)+fx
          mom(i,2)=mom(i,2)+fy
          mom(i,3)=mom(i,3)+fxy
          dv(i) =dv(i)*thk(i)
          mx = ipartc(i)
          evis(8,mx)=evis(8,mx) 
     .             +(fx*vdef(i,6)+fy*vdef(i,7)+fxy*vdef(i,8))*dv(i)
         ENDDO
        ENDIF
C
      RETURN

cbavisnp1()

subroutine cbavisnp1	(	integer	jft,
		integer	jlt,
			vxyz,
			rxyz,
			vcore,
			amu,
			off,
			rho,
			ssp,
			area,
			thk,
			g,
			dt1,
			vf,
		integer, dimension(*)	ipartc,
			evis,
		integer	kfts )

Definition at line 100 of file cbavisc.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "param_c.inc"
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT, JLT ,IPARTC(*)
      my_real
     .   off(*),vf(mvsiz,3,4),g(*),
     .   vcore(mvsiz,3,4),vxyz(mvsiz,3,4),rxyz(mvsiz,2,4),
     .   amu(*),rho(*),ssp(*),area(*),thk(*),evis(npsav,*),dt1
      INTEGER KFTS
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,MX, IC, II, J, JST(MVSIZ+1)
        my_real 
     .  c2,hsura,mx23,my23,mx34,my34,mx13r,my13r,
     .  sc6_v,sc5_v,cxz_v,cyz_v,ss3_v,hvl,rx13,ry13,rx24,ry24,
     .  rxsom,rysom,tesy(mvsiz),vhxz(mvsiz),vhyz(mvsiz)
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
       DO i=jft,jlt
        rx13 =rxyz(i,1,1)
        rx24 =rxyz(i,1,2)
        rxsom=rxyz(i,1,4)
        ry13 =rxyz(i,2,1)
        ry24 =rxyz(i,2,2)
        rysom=rxyz(i,2,4)
        my13r=vcore(i,3,4)*rxsom
        mx13r=vcore(i,2,4)*rysom
        vhxz(i)=vxyz(i,3,3)-(my13r-vcore(i,2,3)*(rx13+rx24)
     1                     -mx13r+vcore(i,1,3)*(ry13+ry24))*fourth
        vhyz(i)=vxyz(i,3,3)-(my13r-vcore(i,1,4)*(rx13-rx24)
     1                     -mx13r+vcore(i,3,3)*(ry13-ry24))*fourth
       ENDDO
C
       DO i=jft,jlt
         mx23=vcore(i,1,3)
         mx34=vcore(i,3,3)
         my23=vcore(i,2,3)
         my34=vcore(i,1,4)
         c2 =one_over_12*g(i)*rho(i)*area(i)
         hvl = two1000*amu(i)*sqrt(c2)*off(i)
         cxz_v=(my34*my34+mx34*mx34)*hvl
         cyz_v=(my23*my23+mx23*mx23)*hvl
         hsura=thk(i)/area(i)
         sc5_v=cxz_v*vhxz(i)*hsura
         sc6_v=cyz_v*vhyz(i)*hsura
         ss3_v = sc5_v+sc6_v
C--------NOEUDS 1,3-------
         vf(i,3,3)=vf(i,3,3)+ss3_v
C--------NOEUDS 2,4-------
         vf(i,3,4)=vf(i,3,4)-ss3_v
         tesy(i)= (sc5_v*vhxz(i)+sc6_v*vhyz(i))*dt1
       ENDDO 
C
       mx = ipartc(jft)
       DO i=jft,jlt
         evis(8,mx)=evis(8,mx) + tesy(i)
       ENDDO
C
      RETURN