szstrainhg.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/.
 
!||====================================================================
!||    szstrainhg   ../engine/source/elements/solid/solidez/szstrainhg.F
!||--- called by ------------------------------------------------------
!||    szhour3      ../engine/source/elements/solid/solidez/szhour3.F
!||    szhour3_or   ../engine/source/elements/solid/solidez/szhour3_or.F
!||    szhour_ctl   ../engine/source/elements/solid/solidez/szhour_ctl.F
!||====================================================================
      SUBROUTINE szstrainhg(
     1   JR_1,    JS_1,    JT_1,    STRHG,
     2   NEL,     HGX1,    HGX2,    HGX3,
     3   HGX4,    HGY1,    HGY2,    HGY3,
     4   HGY4,    HGZ1,    HGZ2,    HGZ3,
     5   HGZ4,    NU,      NU1)
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-----------------------------------------------
#include      "com08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL
      my_real
     .   HGX1(MVSIZ), HGX2(MVSIZ), HGX3(MVSIZ), HGX4(MVSIZ),
     .   HGY1(MVSIZ), HGY2(MVSIZ), HGY3(MVSIZ), HGY4(MVSIZ),
     .   HGZ1(MVSIZ), HGZ2(MVSIZ), HGZ3(MVSIZ), HGZ4(MVSIZ),
     .   jr_1(*),js_1(*),jt_1(*) ,strhg(nel,3,6),nu(*),nu1(*)
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
     .   HGX(MVSIZ,4), HGY(MVSIZ,4), HGZ(MVSIZ,4),
     .   KSI,ETA,ZETA,COEF,NJR_1(MVSIZ),NJS_1(MVSIZ),NJT_1(MVSIZ),
     .   JSR_1(MVSIZ),JRS_1(MVSIZ),JTR_1(MVSIZ),JRT_1(MVSIZ),JTS_1(MVSIZ),
     .   jst_1(mvsiz),f1,f2,f3,f1x,f2y,f3z
C-----------------------------------------------
      coef = eight
      DO i=1,nel
       njr_1(i) = coef*jr_1(i)*dt1
       njs_1(i) = coef*js_1(i)*dt1
       njt_1(i) = coef*jt_1(i)*dt1
      ENDDO
      DO i=1,nel
       hgx(i,1) = njr_1(i)*hgx1(i)
       hgx(i,2) = njr_1(i)*hgx2(i)
       hgx(i,3) = njr_1(i)*hgx3(i)
       hgx(i,4) = njr_1(i)*hgx4(i)
       hgy(i,1) = njs_1(i)*hgy1(i)
       hgy(i,2) = njs_1(i)*hgy2(i)
       hgy(i,3) = njs_1(i)*hgy3(i)
       hgy(i,4) = njs_1(i)*hgy4(i)
       hgz(i,1) = njt_1(i)*hgz1(i)
       hgz(i,2) = njt_1(i)*hgz2(i)
       hgz(i,3) = njt_1(i)*hgz3(i)
       hgz(i,4) = njt_1(i)*hgz4(i)
      ENDDO
C------ STRHG(NEL,3,j) : J=1,3: term ksi,eta,zeta of Eii; J=4: shear;   
C------                  J=5: mode 4 HGi4; J=6: nu*HGi4;   
         DO i=1,nel
          f1x = hgx(i,4)
          f2y = hgy(i,4)
          f3z = hgz(i,4)
          strhg(i,1,1)= strhg(i,1,1)-nu1(i)*(hgy(i,3)+hgz(i,2))
          strhg(i,1,2)= strhg(i,1,2)+hgx(i,3)
          strhg(i,1,3)= strhg(i,1,3)+hgx(i,2)
          strhg(i,1,4)= strhg(i,1,4)+hgx1(i)*njs_1(i)+hgy2(i)*njr_1(i)
          strhg(i,1,5)= strhg(i,1,5)+f1x
          strhg(i,1,6)= strhg(i,1,6)-nu(i)*f1x
          strhg(i,2,1)= strhg(i,2,1)+hgy(i,3)
          strhg(i,2,2)= strhg(i,2,2)-nu1(i)*(hgx(i,3)+hgz(i,1))
          strhg(i,2,3)= strhg(i,2,3)+hgy(i,1)
          strhg(i,2,4)= strhg(i,2,4)+hgx1(i)*njt_1(i)+hgz3(i)*njr_1(i)
          strhg(i,2,5)= strhg(i,2,5)+f2y
          strhg(i,2,6)= strhg(i,2,6)-nu(i)*f2y
          strhg(i,3,1)= strhg(i,3,1)+hgz(i,2)
          strhg(i,3,2)= strhg(i,3,2)+hgz(i,1)
          strhg(i,3,3)= strhg(i,3,3)-nu1(i)*(hgx(i,2)+hgy(i,1))
          strhg(i,3,4)= strhg(i,3,4)+hgy2(i)*njt_1(i)+hgz3(i)*njs_1(i)
          strhg(i,3,5)= strhg(i,3,5)+f3z
          strhg(i,3,6)= strhg(i,3,6)-nu(i)*f3z
         END DO
C      
      RETURN
      END