szsigpara.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/.
!||====================================================================
!||    szsigpara       ../engine/source/elements/solid/solidez/szsigpara.F
!||--- called by ------------------------------------------------------
!||    strs_tenscor3   ../engine/source/output/h3d/h3d_results/strs_tenscor3.F
!||    tensgps3        ../engine/source/output/anim/generate/tensor6.F
!||    tensgps_skin    ../engine/source/output/anim/generate/tensor6.F
!||====================================================================
      SUBROUTINE szsigpara(JR0  ,JS0  ,JT0  ,FHOUR  ,SIG0 ,
     .                     SIG  ,KSI  ,ETA  ,ZETA   ,NU   ,
     .                     NEL  ,I  )
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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NEL,I
      my_real
     .   fhour(nel,3,4),jr0(*),js0(*),jt0(*) ,
     .   sig(mvsiz,6) ,sig0(nel,6),nu(*),ksi  ,eta  ,zeta
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real
     .   coef,nu1,nu2,nu3,sxx,syy,szz,sxy,sxz,syz
C---JR0=-X1+X2+X3-X4-X5+X6+X7-X8 (local coordinates)     
C---JS0=-Y1-Y2+Y3+Y4-Y5-Y6+Y7+Y8     
C---JT0=-Z1-Z2-Z3-Z4+Z5+Z6+Z7+Z8
C---NU : Poisson's coefficient     
C-----------------------------------------------
        coef=twenty4
c        DO I=LFT,LLT
         nu1=one/(one-nu(i))
         nu2=nu(i)*nu1
         nu3=one+nu(i)
         sxx=nu1*(fhour(i,1,2)*zeta+fhour(i,1,3)*eta)+
     .       nu2*(fhour(i,2,1)*zeta+fhour(i,3,1)*eta)+         
     .       nu3*fhour(i,1,4)*zeta*eta         
         syy=nu1*(fhour(i,2,1)*zeta+fhour(i,2,3)*ksi)+
     .       nu2*(fhour(i,1,2)*zeta+fhour(i,3,2)*ksi)+         
     .       nu3*fhour(i,2,4)*zeta*ksi         
         szz=nu1*(fhour(i,3,1)*eta+fhour(i,3,2)*ksi)+
     .       nu2*(fhour(i,1,3)*eta+fhour(i,2,3)*ksi)+         
     .       nu3*fhour(i,3,4)*eta*ksi         
         sig(i,1)=sig0(i,1)+two*sxx*twenty4
         sig(i,2)=sig0(i,2)+two*syy*twenty4
         sig(i,3)=sig0(i,3)+two*szz*twenty4
c        ENDDO
C        
c        DO I=LFT,LLT
         sxy=zeta*(jr0(i)*fhour(i,1,1)/js0(i)+
     .             js0(i)*fhour(i,2,2)/jr0(i))         
         sxz=eta*(jr0(i)*fhour(i,1,1)/jt0(i)+
     .            jt0(i)*fhour(i,3,3)/jr0(i))
         syz=ksi*(js0(i)*fhour(i,2,2)/jt0(i)+
     .            jt0(i)*fhour(i,3,3)/js0(i))
         sig(i,4)=sig0(i,4)+sxy*twenty4         
         sig(i,5)=sig0(i,5)+syz*twenty4
         sig(i,6)=sig0(i,6)+sxz*twenty4
c        ENDDO
C
      RETURN
      END