s8e__sig__a_8F_source.html

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/.

!||====================================================================

!||    s8e_sig_a   ../engine/source/elements/solid/solide8e/s8e_sig_a.F

!||--- called by ------------------------------------------------------

!||    s8eforc3    ../engine/source/elements/solid/solide8e/s8eforc3.F

!||====================================================================


      SUBROUTINE s8e_sig_a(

     1   C1,      G,       OFF,     BXY1,

     2   BXY2,    BXY3,    BXY4,    BXY5,

     3   BXY6,    BXY7,    BXY8,    BYX1,

     4   BYX2,    BYX3,    BYX4,    BYX5,

     5   BYX6,    BYX7,    BYX8,    BXZ1,

     6   BXZ2,    BXZ3,    BXZ4,    BXZ5,

     7   BXZ6,    BXZ7,    BXZ8,    BZX1,

     8   BZX2,    BZX3,    BZX4,    BZX5,

     9   BZX6,    BZX7,    BZX8,    BYZ1,

     A   BYZ2,    BYZ3,    BYZ4,    BYZ5,

     B   BYZ6,    BYZ7,    BYZ8,    BZY1,

     C   BZY2,    BZY3,    BZY4,    BZY5,

     D   BZY6,    BZY7,    BZY8,    VX1,

     E   VX2,     VX3,     VX4,     VX5,

     F   VX6,     VX7,     VX8,     VY1,

     G   VY2,     VY3,     VY4,     VY5,

     H   VY6,     VY7,     VY8,     VZ1,

     I   VZ2,     VZ3,     VZ4,     VZ5,

     J   VZ6,     VZ7,     VZ8,     PIJ,

     K   SIG_A,   EINT,    VOL,     D1,

     L   D2,      D3,      SIG,     S1,

     M   S2,      S3,      NEL,     OFFG)

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      "com08_c.inc"

C-----------------------------------------------

C   D u m m y   A r g u m e n t s

C-----------------------------------------------

      INTEGER NEL

      my_real

     .   VX1(*), VX2(*), VX3(*), VX4(*), VX5(*), VX6(*), VX7(*), VX8(*),

     .   VY1(*), VY2(*), VY3(*), VY4(*), VY5(*), VY6(*), VY7(*), VY8(*),

     .   VZ1(*), VZ2(*), VZ3(*), VZ4(*), VZ5(*), VZ6(*), VZ7(*), VZ8(*),

     .   BXY1(*),BXY2(*),BXY3(*),BXY4(*),

     .   BXY5(*),BXY6(*),BXY7(*),BXY8(*),

     .   BYX1(*),BYX2(*),BYX3(*),BYX4(*),

     .   BYX5(*),BYX6(*),BYX7(*),BYX8(*),

     .   BXZ1(*),BXZ2(*),BXZ3(*),BXZ4(*),

     .   BXZ5(*),BXZ6(*),BXZ7(*),BXZ8(*),

     .   BZX1(*),BZX2(*),BZX3(*),BZX4(*),

     .   BZX5(*),BZX6(*),BZX7(*),BZX8(*),

     .   BYZ1(*),BYZ2(*),BYZ3(*),BYZ4(*),

     .   BYZ5(*),BYZ6(*),BYZ7(*),BYZ8(*),

     .   BZY1(*),BZY2(*),BZY3(*),BZY4(*),

     .   BZY5(*),BZY6(*),BZY7(*),BZY8(*),

     .   PIJ(NEL,*),C1,G,SIG_A(MVSIZ,3),OFF(*),EINT(*),

     .   VOL(*),D1(*),D2(*),D3(*),SIG(NEL,6),S1(*),S2(*),S3(*),OFFG(*)

C-----------------------------------------------

C   L o c a l   V a r i a b l e s

C-----------------------------------------------

      INTEGER I,J,K

      my_real

     .   dc1,dg2,dxx(mvsiz),dyy(mvsiz),dzz(mvsiz),ekk,dav,

     .   dt05,e1,e2,e3,de,p2

C PIJ(NEL,25-27) are used call stress with v-terms

C-----------------------------------------------

       IF (tt==zero) pij(1:nel,25:27)=zero

C----1,compute Dii with v-terms, compute incremental stress

       DO i=1,nel

               IF (offg(i) >one) cycle

           dyy(i) =bxy1(i)*vx1(i)+bxy2(i)*vx2(i)

     +              +bxy3(i)*vx3(i)+bxy4(i)*vx4(i)

     +              +bxy5(i)*vx5(i)+bxy6(i)*vx6(i)

     +              +bxy7(i)*vx7(i)+bxy8(i)*vx8(i)

     +              +bzy1(i)*vz1(i)+bzy2(i)*vz2(i)

     +              +bzy3(i)*vz3(i)+bzy4(i)*vz4(i)

     +              +bzy5(i)*vz5(i)+bzy6(i)*vz6(i)

     +              +bzy7(i)*vz7(i)+bzy8(i)*vz8(i)

           dzz(i) =bxz1(i)*vx1(i)+bxz2(i)*vx2(i)

     +              +bxz3(i)*vx3(i)+bxz4(i)*vx4(i)

     +              +bxz5(i)*vx5(i)+bxz6(i)*vx6(i)

     +              +bxz7(i)*vx7(i)+bxz8(i)*vx8(i)

     +              +byz1(i)*vy1(i)+byz2(i)*vy2(i)

     +              +byz3(i)*vy3(i)+byz4(i)*vy4(i)

     +              +byz5(i)*vy5(i)+byz6(i)*vy6(i)

     +              +byz7(i)*vy7(i)+byz8(i)*vy8(i)

           dxx(i) =byx1(i)*vy1(i)+byx2(i)*vy2(i)

     +              +byx3(i)*vy3(i)+byx4(i)*vy4(i)

     +              +byx5(i)*vy5(i)+byx6(i)*vy6(i)

     +              +byx7(i)*vy7(i)+byx8(i)*vy8(i)

     +              +bzx1(i)*vz1(i)+bzx2(i)*vz2(i)

     +              +bzx3(i)*vz3(i)+bzx4(i)*vz4(i)

     +              +bzx5(i)*vz5(i)+bzx6(i)*vz6(i)

     +              +bzx7(i)*vz7(i)+bzx8(i)*vz8(i)

       ENDDO

       dc1=c1*dt1

       dg2=two*g*dt1

       sig_a(1:nel,1:3)=pij(1:nel,25:27)

       DO i=1,nel

               IF (offg(i) >one) cycle

         ekk=(dxx(i)+dyy(i)+dzz(i))*dc1

c       IF (OFFG(I) >ONE) EKK=ZERO

         dav=-third*(dxx(i)+dyy(i)+dzz(i))

         pij(i,25)=pij(i,25)*off(i)+ekk+dg2*(dxx(i)+dav)

         pij(i,26)=pij(i,26)*off(i)+ekk+dg2*(dyy(i)+dav)

         pij(i,27)=pij(i,27)*off(i)+ekk+dg2*(dzz(i)+dav)

       END DO

       dt05 =half*dt1

      DO i=1,nel

              IF (offg(i) >one) cycle

       e1=(dxx(i)+d1(i))*(sig_a(i,1)+pij(i,25))+dxx(i)*(s1(i)+sig(i,1))

       e2=(dyy(i)+d2(i))*(sig_a(i,2)+pij(i,26))+dyy(i)*(s2(i)+sig(i,2))

       e3=(dzz(i)+d3(i))*(sig_a(i,3)+pij(i,27))+dzz(i)*(s3(i)+sig(i,3))

       de= (e1+e2+e3)*dt05*off(i)

       eint(i)=eint(i)+de

       sig_a(i,1:3)=pij(i,25:27)

      ENDDO

C-----------

      RETURN


      END

s8e_sig_a
subroutine s8e_sig_a(c1, g, off, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, pij, sig_a, eint, vol, d1, d2, d3, sig, s1, s2, s3, nel, offg)
Definition s8e_sig_a.F:51