c3fcum3.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/.
!||====================================================================
!||    c3fcum3       ../engine/source/elements/sh3n/coque3n/c3fcum3.F
!||--- called by ------------------------------------------------------
!||    c3forc3       ../engine/source/elements/sh3n/coque3n/c3forc3.F
!||    c3forc3_crk   ../engine/source/elements/xfem/c3forc3_crk.F
!||====================================================================
      SUBROUTINE c3fcum3(JFT,JLT,F,
     .             F11,F12,F13,F21,F22,
     .             F23,F31,F32,F33,FZERO,
     .             E1X      ,E1Y      ,E1Z     ,E2X      ,                 
     .             E2Y      ,E2Z      ,E3X     ,E3Y      ,E3Z      )       
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 JFT, JLT
C     REAL
      my_real F(3,*),FZERO(3,4,*),
     .     E1X(MVSIZ), E1Y(MVSIZ), E1Z(MVSIZ),
     .     e2x(mvsiz), e2y(mvsiz), e2z(mvsiz),
     .     e3x(mvsiz), e3y(mvsiz), e3z(mvsiz),
     .     f11(mvsiz), f12(mvsiz), f13(mvsiz),
     .     f21(mvsiz), f22(mvsiz), f23(mvsiz),
     .     f31(mvsiz), f32(mvsiz), f33(mvsiz)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com08_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J
      my_real G11, G12, G13, G21, G22, G23, G31, G32, G33,
     .     TEMP1,TEMP2,TEMP3
C-----------------------------------------------
C               POUR ISIGI = 5 OU -5
C     CHANGEMENT DE REPERE DE FZERO POUR TT=DT1
C          REPERE GLOBAL ---> REPERE LOCAL
C-----------------------------------------------
      IF(iabs(isigi)==5.AND.ncycle==1.AND.irun==1) THEN
        DO i=jft,jlt
          DO j=1,3
            temp1=fzero(1,j,i)
            temp2=fzero(2,j,i)
            temp3=fzero(3,j,i)
            fzero(1,j,i)=temp1*e1x(i)+temp2*e1y(i)+temp3*e1z(i)
            fzero(2,j,i)=temp1*e2x(i)+temp2*e2y(i)+temp3*e2z(i)
            fzero(3,j,i)=temp1*e3x(i)+temp2*e3y(i)+temp3*e3z(i)
          ENDDO
        ENDDO
      ENDIF
C-----------------------------------------------
C               POUR ISIGI = 5 OU -5
C     PRISE EN COMPTE DES FORCES D EQUILIBRAGE
C      DANS L EVALUATION DES EFFORTS INTERNES
C-----------------------------------------------
      IF((iabs(isigi)==5.AND.ncycle>=1.AND.irun==1).OR.
     .   (iabs(isigi)==5.AND.irun>1)) THEN
        DO i=jft,jlt
          g11 = f11(i)+fzero(1,1,i)
          g21 = f21(i)+fzero(2,1,i)
          g31 = f31(i)+fzero(3,1,i)
          f11(i) = e1x(i)*g11 + e2x(i)*g21 + e3x(i)*g31
          f21(i) = e1y(i)*g11 + e2y(i)*g21 + e3y(i)*g31
          f31(i) = e1z(i)*g11 + e2z(i)*g21 + e3z(i)*g31
          g12 = f12(i)+fzero(1,2,i)
          g22 = f22(i)+fzero(2,2,i)
          g32 = f32(i)+fzero(3,2,i)
          f12(i) = e1x(i)*g12 + e2x(i)*g22 + e3x(i)*g32
          f22(i) = e1y(i)*g12 + e2y(i)*g22 + e3y(i)*g32
          f32(i) = e1z(i)*g12 + e2z(i)*g22 + e3z(i)*g32
          g13 = f13(i)+fzero(1,3,i)
          g23 = f23(i)+fzero(2,3,i)
          g33 = f33(i)+fzero(3,3,i)
          f13(i) = e1x(i)*g13 + e2x(i)*g23 + e3x(i)*g33
          f23(i) = e1y(i)*g13 + e2y(i)*g23 + e3y(i)*g33
          f33(i) = e1z(i)*g13 + e2z(i)*g23 + e3z(i)*g33
        ENDDO
      ENDIF
      IF(iabs(isigi)/=5.OR.tt==0) THEN
        DO i=jft,jlt
            g11 = f11(i)
            g21 = f21(i)
            g31 = f31(i)
            f11(i) = e1x(i)*g11 + e2x(i)*g21 + e3x(i)*g31
            f21(i) = e1y(i)*g11 + e2y(i)*g21 + e3y(i)*g31
            f31(i) = e1z(i)*g11 + e2z(i)*g21 + e3z(i)*g31
            g12 = f12(i)
            g22 = f22(i)
            g32 = f32(i)
            f12(i) = e1x(i)*g12 + e2x(i)*g22 + e3x(i)*g32
            f22(i) = e1y(i)*g12 + e2y(i)*g22 + e3y(i)*g32
            f32(i) = e1z(i)*g12 + e2z(i)*g22 + e3z(i)*g32
            g13 = f13(i)
            g23 = f23(i)
            g33 = f33(i)
            f13(i) = e1x(i)*g13 + e2x(i)*g23 + e3x(i)*g33
            f23(i) = e1y(i)*g13 + e2y(i)*g23 + e3y(i)*g33
            f33(i) = e1z(i)*g13 + e2z(i)*g23 + e3z(i)*g33
        ENDDO
      ENDIF
C
      RETURN
      END