13 element projections

[halleluja]

Hi everyone,

i want to know how i can project 13 value to nodes in the case of 6-nodes element mesh. Therefore i have progammed a 6-nodes element with 7 gauss points subroutine elmt49, see elmt49.f in attachment.
For projection first do local least squares on each element to get values to element nodes, then average.
For simplicity i set from the 4th projected value to the 13th projected  value zero.
At last i plot cont 1, but the output is not smooth comparing to the case if only 7 values is projected, see pictures in attachment. But i thought the output of cases that 13 projected values and 7 projected values should be the same.

P.S. i use feap8.5 and in attachment you find the elmt49.f, inputfile Imesh49 and two pictures of the plot.


Best

[Prof. R.L. Taylor]

The default projection for 6-node triangles is a linear projection to the vertex nodes and then an average to the midside nodes.  See (I think) ./elements/solid2d/slcn2d.f

[halleluja]

The default projection for 6-node triangles is a linear projection to the vertex nodes and then an average to the midside nodes.  See (I think) ./elements/solid2d/slcn2d.f

yes, that i have done in my code elmt49.f as following

Code: [Select]

c       from slcn2d.f line 108 to 122
c     do l = 1,lint
          do j = 1,3
            xg = el(j,l)*el(4,l)
            do i = 1,3
              mat(i,j) = mat(i,j) + el(i,l)*xg
            end do ! i
            do i = 1,nes
              sigg(j,i) = sigg(j,i) + valp(i)*xg
            end do ! i
          end do ! j
c       end do ! l
c---------------plot------------c
400   continue

      if (isw.eq.8) then     

        iste=nes

c       from slcn2d.f line 123 to 144
      call invert(mat,3,3)
        do i = 1,nes
          do j = 1,3
            gsig(j,i) = mat(j,1)*sigg(1,i)
     &                + mat(j,2)*sigg(2,i)
     &                + mat(j,3)*sigg(3,i)
c            geps(j,i) = mat(j,1)*epsg(1,i)
c     &                + mat(j,2)*epsg(2,i)
c     &                + mat(j,3)*epsg(3,i)
          end do ! j

c         Project stresses

          valplot(1,i) =  gsig(1,i)                         
          valplot(2,i) =  gsig(2,i)
          valplot(3,i) =  gsig(3,i)
          valplot(4,i) = (gsig(1,i) + gsig(2,i))*0.5d0
          valplot(5,i) = (gsig(2,i) + gsig(3,i))*0.5d0
          valplot(6,i) = (gsig(3,i) + gsig(1,i))*0.5d0
        end do ! i