linvec10.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	linvec10 (vec, konl, vecl, istart, iend, shp)

Function/Subroutine Documentation

linvec10()

subroutine linvec10	(	real*8, dimension(istart:iend,*)	vec,
		integer, dimension(*)	konl,
		real*8, dimension(3)	vecl,
		integer	istart,
		integer	iend,
		real*8, dimension(4,*)	shp
	)

!
!     calculates a linear approximation to the quadratic interpolation
!     of the temperatures in a C3D10 element. A
!     quadratic interpolation of the temperatures leads to quadratic
!     thermal stresses, which cannot be handled by the elements 
!     displacement functions (which lead to linear stresses). Thus,
!     the temperatures are approximated by a linear function.
!
      implicit none
!
      integer konl(*),j,istart,iend
!
      real*8 vec(istart:iend,*),vecl(3),shp(4,*)
!
      do j=1,3
         vecl(j)=
     &    (shp(4,1)+(shp(4,5)+shp(4,7)+shp(4,8))/2.d0)*vec(j,konl(1))
     &    +(shp(4,2)+(shp(4,5)+shp(4,6)+shp(4,9))/2.d0)*vec(j,konl(2))
     &    +(shp(4,3)+(shp(4,6)+shp(4,7)+shp(4,10))/2.d0)*vec(j,konl(3))
     &    +(shp(4,4)+(shp(4,8)+shp(4,9)+shp(4,10))/2.d0)*vec(j,konl(4))
      enddo
!     
      return