updatecont2d.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	updatecont2d (koncont, ncont, co, cg, straight)

Function/Subroutine Documentation

updatecont2d()

subroutine updatecont2d	(	integer, dimension(3,*)	koncont,
		integer	ncont,
		real*8, dimension(2,*)	co,
		real*8, dimension(2,*)	cg,
		real*8, dimension(9,*)	straight
	)

!
!     update geometric date of the contact master surface triangulation
!
      implicit none
!
      integer koncont(3,*),ncont,i,j,k,node
!
      real*8 co(2,*),cg(2,*),straight(9,*),col(2,3)
!
      do i=1,ncont
         do j=1,3
            node=koncont(j,i)
            do k=1,2
               col(k,j)=co(k,node)
            enddo
         enddo
!
!        center of gravity of the triangles
!
         do k=1,2
            cg(k,i)=col(k,1)
         enddo
         do j=2,3
            do k=1,2
               cg(k,i)=cg(k,i)+col(k,j)
            enddo
         enddo
         do k=1,2
            cg(k,i)=cg(k,i)/3.d0
         enddo
!
!        calculating the equation of the triangle plane and the planes
!        perpendicular on it and through the triangle edges
!
         call straighteq2d(col,straight(1,i))
!
      enddo
!
      return