objective_stress_dx.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	objective_stress_dx (nodeset, istartset, iendset, ialset, nk, idesvarc, iobject, dgdx, ndesi, nobject, stn, dstn, objectset, g0)

Function/Subroutine Documentation

objective_stress_dx()

subroutine objective_stress_dx	(	integer	nodeset,
		integer, dimension(*)	istartset,
		integer, dimension(*)	iendset,
		integer, dimension(*)	ialset,
		integer	nk,
		integer	idesvarc,
		integer	iobject,
		real*8, dimension(ndesi,nobject)	dgdx,
		integer	ndesi,
		integer	nobject,
		real*8, dimension(6,*)	stn,
		real*8, dimension(6,*)	dstn,
		character*81, dimension(4,*)	objectset,
		real*8, dimension(nobject)	g0
	)

!
!     calculates derivative of the sum of the square of the 
!     von Mises stress of a node set w.r.t. the coordinates of the mesh
!
      implicit none
!
      character*81 objectset(4,*)
!
      integer nk,istartset(*),iendset(*),ialset(*),nodeset,
     &  idesvarc,iobject,j,k,ndesi,nobject,idesvar
!
      real*8 dgdx(ndesi,nobject),stn(6,*),dstn(6,*),g0(nobject),
     &  rho,xstress,svm,dsvmdx,p
!
      idesvar=idesvarc+1
!
!     reading rho and the mean stress for the Kreisselmeier-Steinhauser
!     function
!
      read(objectset(2,iobject)(41:60),'(f20.0)') rho
      read(objectset(2,iobject)(61:80),'(f20.0)') xstress
!
!     check for the existence of a set, else take the complete mesh
!
      if(nodeset.eq.0) then
         do j=1,nk
            p=-(stn(1,j)+stn(2,j)+stn(3,j))/3.d0
            svm=dsqrt(1.5d0*
     &        ((stn(1,j)+p)**2+(stn(2,j)+p)**2+(stn(3,j)+p)**2+
     &         2.d0*(stn(4,j)**2+stn(5,j)**2+stn(6,j)**2)))
            dsvmdx=
     &           (2.d0*stn(1,j)-stn(2,j)-stn(3,j))*dstn(1,j)+
     &           (2.d0*stn(2,j)-stn(1,j)-stn(3,j))*dstn(2,j)+
     &           (2.d0*stn(3,j)-stn(1,j)-stn(2,j))*dstn(3,j)+
     &           6.d0*(stn(4,j)*dstn(4,j)+stn(5,j)*dstn(5,j)+
     &                 stn(6,j))
            dgdx(idesvar,iobject)=dgdx(idesvar,iobject)+
     &           dexp(rho*svm/xstress)*dsvmdx
         enddo
         dgdx(idesvar,iobject)=dgdx(idesvar,iobject)*
     &         dexp(-rho*g0(iobject))/xstress
      else
         do j=istartset(nodeset),iendset(nodeset)
            if(ialset(j).gt.0) then
               k=ialset(j)
               p=-(stn(1,k)+stn(2,k)+stn(3,k))/3.d0
               svm=dsqrt(1.5d0*
     &              ((stn(1,k)+p)**2+(stn(2,k)+p)**2+(stn(3,k)+p)**2+
     &              2.d0*(stn(4,k)**2+stn(5,k)**2+stn(6,k)**2)))
               dsvmdx=
     &              (2.d0*stn(1,k)-stn(2,k)-stn(3,k))*dstn(1,k)+
     &              (2.d0*stn(2,k)-stn(1,k)-stn(3,k))*dstn(2,k)+
     &              (2.d0*stn(3,k)-stn(1,k)-stn(2,k))*dstn(3,k)+
     &              6.d0*(stn(4,k)*dstn(4,k)+stn(5,k)*dstn(5,k)+
     &              stn(6,k))
               dgdx(idesvar,iobject)=dgdx(idesvar,iobject)+
     &              dexp(rho*svm/xstress)*dsvmdx
            else
               k=ialset(j-2)
               do
                  k=k-ialset(j)
                  if(k.ge.ialset(j-1)) exit
                  p=-(stn(1,k)+stn(2,k)+stn(3,k))/3.d0
                  svm=dsqrt(1.5d0*
     &                 ((stn(1,k)+p)**2+(stn(2,k)+p)**2+(stn(3,k)+p)**2+
     &                 2.d0*(stn(4,k)**2+stn(5,k)**2+stn(6,k)**2)))
                  dsvmdx=
     &                 (2.d0*stn(1,k)-stn(2,k)-stn(3,k))*dstn(1,k)+
     &                 (2.d0*stn(2,k)-stn(1,k)-stn(3,k))*dstn(2,k)+
     &                 (2.d0*stn(3,k)-stn(1,k)-stn(2,k))*dstn(3,k)+
     &                 6.d0*(stn(4,k)*dstn(4,k)+stn(5,k)*dstn(5,k)+
     &                 stn(6,k))
                  dgdx(idesvar,iobject)=dgdx(idesvar,iobject)+
     &                 dexp(rho*svm/xstress)*dsvmdx
               enddo
            endif
         enddo
         dgdx(idesvar,iobject)=dgdx(idesvar,iobject)*
     &         dexp(-rho*g0(iobject))/xstress
      endif
!     
      return