szstraingps.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    szstraingps      ../engine/source/elements/solid/solidez/szstraingps.F
!||--- called by ------------------------------------------------------
!||    gpsstrain_skin   ../engine/source/output/anim/generate/tensgpstrain.F
!||    strn_tenscor3    ../engine/source/output/h3d/h3d_results/h3d_strn_tenscor3.F
!||    tensgpstrain     ../engine/source/output/anim/generate/tensgpstrain.F
!||====================================================================
      SUBROUTINE szstraingps(
     1   STRAIN,  STR_PI,  STRHG,   NEL)
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL
      my_real
     .   strain(nel,6) ,str_pi(mvsiz,6,8),strhg(nel,3,6)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,  J,IR,IS,IT,IP
      my_real
     .   ksi,eta,zeta,coef,f1,f2,f3,f1x,f2y,f3z
      my_real
     .  a_gauss(9,9)
      DATA a_gauss / 
     1 0.               ,0.               ,0.               ,
     1 0.               ,0.               ,0.               ,
     1 0.               ,0.               ,0.               ,
     2 -.577350269189626,0.577350269189626,0.               ,
     2 0.               ,0.               ,0.               ,
     2 0.               ,0.               ,0.               ,
     3 -.774596669241483,0.               ,0.774596669241483,
     3 0.               ,0.               ,0.               ,
     3 0.               ,0.               ,0.               ,
     4 -.861136311594053,-.339981043584856,0.339981043584856,
     4 0.861136311594053,0.               ,0.               ,
     4 0.               ,0.               ,0.               ,
     5 -.906179845938664,-.538469310105683,0.               ,
     5 0.538469310105683,0.906179845938664,0.               ,
     5 0.               ,0.               ,0.               ,
     6 -.932469514203152,-.661209386466265,-.238619186083197,
     6 0.238619186083197,0.661209386466265,0.932469514203152,
     6 0.               ,0.               ,0.               ,
     7 -.949107912342759,-.741531185599394,-.405845151377397,
     7 0.               ,0.405845151377397,0.741531185599394,
     7 0.949107912342759,0.               ,0.               ,
     8 -.960289856497536,-.796666477413627,-.525532409916329,
     8 -.183434642495650,0.183434642495650,0.525532409916329,
     8 0.796666477413627,0.960289856497536,0.               ,
     9 -.968160239507626,-.836031107326636,-.613371432700590,
     9 -.324253423403809,0.               ,0.324253423403809,
     9 0.613371432700590,0.836031107326636,0.968160239507626/
      
C-----------------------------------------------
      DO ir=1,2
       DO is=1,2
        DO it=1,2
         ip = ir + ( (is-1) + (it-1)*2 )*2
         ksi = a_gauss(ir,2)
         eta = a_gauss(is,2)
         zeta = a_gauss(it,2)
         f1 = eta*zeta
         f2 = ksi*zeta
         f3 = eta*ksi
C
         DO i=1,nel
          str_pi(i,1,ip)=strain(i,1)+
     +                   strhg(i,1,1)*ksi+strhg(i,1,2)*eta+strhg(i,1,3)*zeta+
     +                   strhg(i,1,5)*f1+strhg(i,2,6)*f2+strhg(i,3,6)*f3
          str_pi(i,2,ip)=strain(i,2)+
     +                   strhg(i,2,1)*ksi+strhg(i,2,2)*eta+strhg(i,2,3)*zeta+
     +                   strhg(i,2,5)*f2+strhg(i,1,6)*f1+strhg(i,3,6)*f3
          str_pi(i,3,ip)=strain(i,3)+
     +                   strhg(i,3,1)*ksi+strhg(i,3,2)*eta+strhg(i,3,3)*zeta+
     +                   strhg(i,3,5)*f3+strhg(i,1,6)*f1+strhg(i,2,6)*f2
          str_pi(i,4,ip)=strain(i,4)+strhg(i,1,4)*zeta
          str_pi(i,6,ip)=strain(i,5)+strhg(i,2,4)*eta
          str_pi(i,5,ip)=strain(i,6)+strhg(i,3,4)*ksi
         END DO
        ENDDO 
       ENDDO 
      ENDDO
C      
      RETURN
      END