soltosph_hour.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/.
!||====================================================================
!||    sig_heph1     ../engine/source/elements/sph/soltosph_hour.F
!||--- called by ------------------------------------------------------
!||    soltosphp     ../engine/source/elements/sph/soltosph.F
!||--- uses       -----------------------------------------------------
!||    element_mod   ../common_source/modules/elements/element_mod.F90
!||====================================================================
      SUBROUTINE sig_heph1(
     1   JR0,     JS0,     JT0,     GSIG,
     2   FHOUR,   SIG_HEPH,PM,      IXS,
     3   II,      NEL,     LFT,     LLT)
      use element_mod , only : nixs
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"
#include      "param_c.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, INTENT(INOUT) :: LFT
      INTEGER, INTENT(INOUT) :: LLT
      INTEGER IXS(NIXS,*),II(6),NEL
      my_real
     .   gsig(*),fhour(nel,3,4),jr0(*),js0(*),jt0(*) ,
     .   sig_heph(mvsiz,6,7),pm(npropm,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   NU1,NU2,NU3
C-----------------------------------------------
 
C-----------------------------------------------
C
C       SIG_HEPH(,1) -> mean stress
C       SIG_HEPH(,2) -> hourglass stress component proportional to ZETA
C       SIG_HEPH(,3) -> hourglass stress component proportional to ETA 
C       SIG_HEPH(,4) -> hourglass stress component proportional to KSI 
C       SIG_HEPH(,5) -> hourglass stress component proportional to ZETA*ETA
C       SIG_HEPH(,6) -> hourglass stress component proportional to ZETA*KSI 
C       SIG_HEPH(,7) -> hourglass stress component proportional to ETA*KSI    
C
        DO i=lft,llt
           nu1=one/(one-pm(21,ixs(1,i)))
           nu2=pm(21,ixs(1,i))*nu1
           nu3=one+pm(21,ixs(1,i))
C
           sig_heph(i,1,1)=gsig(ii(1)+i)
           sig_heph(i,2,1)=gsig(ii(2)+i)
           sig_heph(i,3,1)=gsig(ii(3)+i)
           sig_heph(i,4,1)=gsig(ii(4)+i)
           sig_heph(i,5,1)=gsig(ii(5)+i)
           sig_heph(i,6,1)=gsig(ii(6)+i)
C
           sig_heph(i,1,2)=2*twenty4*(nu1*fhour(i,1,2)+nu2*fhour(i,2,1))
           sig_heph(i,2,2)=2*twenty4*(nu1*fhour(i,2,1)+nu2*fhour(i,1,2))
           sig_heph(i,4,2)=twenty4*(jr0(i)*fhour(i,1,1)/js0(i)+js0(i)*fhour(i,2,2)/jr0(i))
           sig_heph(i,1,3)=2*twenty4*(nu1*fhour(i,1,3)+nu2*fhour(i,3,1))
           sig_heph(i,3,3)=2*twenty4*(nu1*fhour(i,3,1)+nu2*fhour(i,1,3))
           sig_heph(i,6,3)=twenty4*(jr0(i)*fhour(i,1,1)/jt0(i)+jt0(i)*fhour(i,3,3)/jr0(i))
           sig_heph(i,2,4)=2*twenty4*(nu1*fhour(i,2,3)+nu2*fhour(i,3,2))
           sig_heph(i,3,4)=2*twenty4*(nu1*fhour(i,3,2)+nu2*fhour(i,2,3))
           sig_heph(i,5,4)=twenty4*(js0(i)*fhour(i,2,2)/jt0(i)+jt0(i)*fhour(i,3,3)/js0(i))
           sig_heph(i,1,5)=2*twenty4*nu3*fhour(i,1,4)        
           sig_heph(i,2,6)=2*twenty4*nu3*fhour(i,2,4)
           sig_heph(i,3,7)=2*twenty4*nu3*fhour(i,3,4)
        END DO
C
      RETURN
      END
C