hgforce.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	hgforce (fn, elas, a, gs, vl, mi, konl)

Function/Subroutine Documentation

hgforce()

subroutine hgforce	(	real*8, dimension(0:mi(2),*)	fn,
		real*8, dimension(1)	elas,
		real*8	a,
		real*8, dimension(8,4)	gs,
		real*8, dimension(0:mi(2),20)	vl,
		integer, dimension(*)	mi,
		integer, dimension(20)	konl
	)

!
!     hourglass control forces for 8-node solid mean strain element
!
!     Reference: Flanagan, D.P., Belytschko, T.; "Uniform  strain hexahedron
!     and quadrilateral with orthogonal Hourglass control". Int. J. Num.
!     Meth. Engg., Vol. 17, 679-706, 1981. 
!
!     author: Otto-Ernst Bernhardi
!
      implicit none
      integer i,j,k,mi(*),konl(20)
      real*8 gs(8,4),a,elas(1),ahr
      real*8 vl(0:mi(2),20), fn(0:mi(2),*)
      real*8 hglf(3,4)
!
      ahr=elas(1)*a
c     write(6,*) "force:", ahr
!
      do i=1,3
         do k=1,4    
            hglf(i,k)=0.0d0
            do j=1,8
               hglf(i,k)=hglf(i,k)+gs(j,k)*vl(i,j)
            enddo
            hglf(i,k)=hglf(i,k)*ahr
         enddo
      enddo
      do i=1,3
         do j=1,8
            do k=1,4
               fn(i,konl(j))=fn(i,konl(j))+hglf(i,k)*gs(j,k)
            enddo
         enddo
      enddo
      return