cdkfint3.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/.
!||====================================================================
!||    cdkfint3   ../engine/source/elements/sh3n/coquedk/cdkfint3.F
!||--- called by ------------------------------------------------------
!||    cdkforc3   ../engine/source/elements/sh3n/coquedk/cdkforc3.F
!||====================================================================
          SUBROUTINE cdkfint3(JFT,JLT,VOL,THK0,FOR,MOM,
     1                      PX2,PY2,PX3,PY3,
     2                      BZ1,BZ2,BZ3,
     3                      BRX1,BRX2,BRX3,BRY1,BRY2,BRY3,
     4                      F11,F12,F13,F21,F22,F23,F32,F33,
     5                      M11,M12,M13,M21,M22,M23,
     6                      NEL)
C--------------------------------------------------------------------------------------------------
C-----------------------------------------------
C   I M P L I C I T   T Y P E S
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D U M M Y   A R G U M E N T S
C-----------------------------------------------
      INTEGER JFT,JLT,NG,NEL
      MY_REAL 
     .   FOR(NEL,5),MOM(NEL,3), 
     .   PX2(*), PY2(*), PX3(*), PY3(*), 
     .   BRX1(MVSIZ,3),BRX2(MVSIZ,3),BRX3(MVSIZ,3),BRY1(MVSIZ,3),BRY2(MVSIZ,3),BRY3(MVSIZ,3),
     .   BZ1(MVSIZ,2),BZ2(MVSIZ,2),BZ3(MVSIZ,2),
     .   f11(*), f12(*), f13(*),
     .   f21(*), f22(*), f23(*), f32(*), f33(*),
     .   m11(*), m12(*), m13(*),
     .   m21(*), m22(*), m23(*),
     .   thk0(*),vol(*)
C-----------------------------------------------
C   L O C A L   V A R I A B L E S
C-----------------------------------------------
      INTEGER I
      MY_REAL 
     .    C2,FX2,FX3,FY2,FY3
C--------------------------------------------------------------------------------------------------
      DO I=jft,jlt
          c2=vol(i)
        fx2 =c2*(px2(i)*for(i,1)+py2(i)*for(i,3))
        fy2 =c2*(py2(i)*for(i,2)+px2(i)*for(i,3))
        fx3 =c2*(px3(i)*for(i,1)+py3(i)*for(i,3))
        fy3 =c2*(py3(i)*for(i,2)+px3(i)*for(i,3))
        f12(i)=f12(i)+fx2
        f22(i)=f22(i)+fy2
        f13(i)=f13(i)+fx3
        f23(i)=f23(i)+fy3
C        F11(I)=F11(I)-FX2-FX3
C        F21(I)=F21(I)-FY2-FY3
      ENDDO     
      DO i=jft,jlt
          c2=vol(i)*thk0(i)
        f32(i)=f32(i)+c2*(bz1(i,1)*mom(i,1)+bz2(i,1)*mom(i,2)+
     .                    bz3(i,1)*mom(i,3))
        f33(i)=f33(i)+c2*(bz1(i,2)*mom(i,1)+bz2(i,2)*mom(i,2)+
     .                    bz3(i,2)*mom(i,3))
        m11(i)=m11(i)+c2*(brx1(i,1)*mom(i,1)+brx2(i,1)*mom(i,2)+
     .                    brx3(i,1)*mom(i,3))
        m21(i)=m21(i)+c2*(bry1(i,1)*mom(i,1)+bry2(i,1)*mom(i,2)+
     .                    bry3(i,1)*mom(i,3))
        m12(i)=m12(i)+c2*(brx1(i,2)*mom(i,1)+brx2(i,2)*mom(i,2)+
     .                    brx3(i,2)*mom(i,3))
        m22(i)=m22(i)+c2*(bry1(i,2)*mom(i,1)+bry2(i,2)*mom(i,2)+
     .                    bry3(i,2)*mom(i,3))
        m13(i)=m13(i)+c2*(brx1(i,3)*mom(i,1)+brx2(i,3)*mom(i,2)+
     .                    brx3(i,3)*mom(i,3))
        m23(i)=m23(i)+c2*(bry1(i,3)*mom(i,1)+bry2(i,3)*mom(i,2)+
     .                    bry3(i,3)*mom(i,3))
      ENDDO     
C    
      RETURN
      END