networkforc.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	networkforc (vl, tnl, imat, konl, mi, ntmat_, shcon, nshcon, rhcon, nrhcon)

Function/Subroutine Documentation

networkforc()

subroutine networkforc	(	real*8, dimension(0:mi(2),20)	vl,
		real*8, dimension(9)	tnl,
		integer	imat,
		integer, dimension(20)	konl,
		integer, dimension(*)	mi,
		integer	ntmat_,
		real*8, dimension(0:3,ntmat_,*)	shcon,
		integer, dimension(*)	nshcon,
		real*8, dimension(0:1,ntmat_,*)	rhcon,
		integer, dimension(*)	nrhcon
	)

!
!     calculates the concentrated flux of a generic networkelement
!     element label: D + blank
!
      implicit none
!
      integer konl(20),mi(*),imat,nshcon(*),nrhcon(*),ntmat_
!
      real*8 vl(0:mi(2),20),tnl(9),gastemp,shcon(0:3,ntmat_,*),
     &  cp,r,dvi,rhcon(0:1,ntmat_,*),rho
!
      gastemp=(vl(0,1)+vl(0,3))/2.d0
!
      call materialdata_tg(imat,ntmat_,gastemp,shcon,nshcon,cp,r,
     &  dvi,rhcon,nrhcon,rho)
!
!     internal force = - external force
!
      if(vl(1,2).gt.0.d0) then
         tnl(3)=cp*(vl(0,3)-vl(0,1))*vl(1,2)
      else
         tnl(1)=-cp*(vl(0,1)-vl(0,3))*vl(1,2)
      endif
!
      return