extrapolate_u1.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	extrapolate_u1 (yi, yn, ipkon, inum, kon, lakon, nfield, nk, ne, mi, ndim, orab, ielorien, co, iorienloc, cflag, vold, force, ielmat, thicke, ielprop, prop, i)

Function/Subroutine Documentation

extrapolate_u1()

subroutine extrapolate_u1	(	real*8, dimension(ndim,mi(1),*), intent(in)	yi,
		real*8, dimension(nfield,*), intent(inout)	yn,
		integer, dimension(*), intent(in)	ipkon,
		integer, dimension(*), intent(in)	inum,
		integer, dimension(*), intent(in)	kon,
		character*8, dimension(*), intent(in)	lakon,
		integer, intent(in)	nfield,
		integer, intent(in)	nk,
		integer, intent(in)	ne,
		integer, dimension(*), intent(in)	mi,
		integer, intent(in)	ndim,
		real*8, dimension(7,*), intent(in)	orab,
		integer, dimension(mi(3),*), intent(in)	ielorien,
		real*8, dimension(3,*), intent(in)	co,
		integer, intent(in)	iorienloc,
		character*1, intent(in)	cflag,
		real*8, dimension(0:mi(2),*), intent(in)	vold,
		logical, intent(in)	force,
		integer, dimension(mi(3),*), intent(in)	ielmat,
		real*8, dimension(mi(3),*), intent(in)	thicke,
		integer, dimension(*), intent(in)	ielprop,
		real*8, dimension(*), intent(in)	prop,
		integer, intent(in)	i
	)

!
!     extrapolates field values at the integration points to the 
!     nodes for user element i of type u1
!
!
!     INPUT:
!
!     yi(ndim,mi(1),*)   value of the variables at the integration
!                        points
!     ipkon(i)           points to the location in field kon preceding
!                        the topology of element i
!     inum(i)            < 0: node i is a network node
!                        > 0: node i is a structural node
!                        =0: node i is not used
!     kon(*)             contains the topology of all elements. The
!                        topology of element i starts at kon(ipkon(i)+1)
!                        and continues until all nodes are covered. The
!                        number of nodes depends on the element label
!     lakon(i)           contains the label of element i
!     nfield             number of variables to be extrapolated
!     nk                 maximum node number in the mesh
!     ne                 maximum element number in the mesh
!     ndim               number of variables in the integration point
!                        field to be extrapolated
!     orab(7,*)          description of all local coordinate systems.
!                        (cf. List of variables and their meaning in the
!                        User's manual)
!     ielorien(i)        orientation in element i
!     co(1..3,i)         global coordinates of node i
!     iorienloc          0: extrapolated variables requested in global 
!                           coordinates
!                        1: extrapolated variables requested in local
!                           coordinates
!     cflag (char*1)     I: interpolate 3D results onto 1D/2D
!                        E: store extrapolated 1D/2D results
!                        M: store 1D section forces
!                        blank: any other case
!     vold(j,i)          value of variable j in node i at the end
!                        of the previous iteration   
!     force              logical variable; if true the values to
!                        be extrapolated are force values; important for
!                        interpolation from 3D expanded structures on the
!                        original 1D/2D structure: forces across the
!                        expansion have to be summed, not interpolated
!     ielmat(i)          material of element i
!     thicke(j,i)        thickness of layer j in node i
!     ielprop(i)         properties for element i are stored in
!                        prop(ielprop(i)+1),prop(ielprop(i)+2),....
!                        (number of properties depends on the type of
!                        element)
!     prop               property field
!     i                  number of the element for which the extrapolation
!                        is to be performed
!
!     OUTPUT:
!
!     yn(nfield,*)       value of the variables at the nodes
!
      implicit none
!
      logical force
!
      character*1 cflag
      character*8 lakon(*)
!
      integer ipkon(*),inum(*),kon(*),mi(*),ne,nfield,nk,i,ndim,
     &  iorienloc,ielorien(mi(3),*),ielmat(mi(3),*),ielprop(*)
!
      real*8 yi(ndim,mi(1),*),yn(nfield,*),orab(7,*),co(3,*),prop(*),
     &  vold(0:mi(2),*),thicke(mi(3),*)
!
      intent(in) yi,ipkon,inum,kon,lakon,nfield,nk,
     &  ne,mi,ndim,orab,ielorien,co,iorienloc,cflag,
     &  vold,force,ielmat,thicke,ielprop,prop,i
!
      intent(inout) yn
!
!     START OF THIS SUBROUTINE
!
      integer indexe,j,k,node
!
      if(iorienloc.ne.0) then
         write(*,*) '*ERROR in extrapolate_u1'
         write(*,*) '       no local orientation for variables'
         write(*,*) '       belonging to this type of element'
         write(*,*) '       allowed'
         call exit(201)
      endif
!
      if(nfield.eq.6) then
         indexe=ipkon(i)
         do j=1,2
            node=kon(indexe+j)
            do k=1,nfield
               yn(k,node)=yn(k,node)+yi(k,1,i)
            enddo
         enddo
      else
         write(*,*) '*ERROR in extrapolate_u1'
         write(*,*) '       extropolation for element of type u1'
         write(*,*) '       is only coded for fields with 6'
         write(*,*) '       entries'
         call exit(201)
      endif
!
!     END OF THIS SUBROUTINE
!
      return