cdk6defo3.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/.
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
!||====================================================================
!||    cdk6defo3   ../engine/source/elements/sh3n/coquedk6/cdk6defo3.F
!||--- called by ------------------------------------------------------
!||    cdk6forc3   ../engine/source/elements/sh3n/coquedk6/cdk6forc3.F
!||====================================================================
      SUBROUTINE cdk6defo3(JFT,JLT,VLX,VLY,VLZ,VZ4,VZ5,VZ6,
     1                     PX2,PY2,PX3,PY3,PB1,PB2,PB3,VDEF)
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
      my_real
     .   px2(*),py2(*),px3(*),py3(*),pb1(mvsiz,3,3),pb2(mvsiz,3,3),pb3(mvsiz,3,6),
     .   vlx(mvsiz,2),vly(mvsiz,2),vlz(mvsiz,6),vz4(mvsiz,3),vz5(mvsiz,3),vz6(mvsiz,3),vdef(mvsiz,8)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real a5,a6,a4
C-----------------------------------------------
C
      DO i=jft,jlt
       vdef(i,1)=px2(i)*vlx(i,1)+px3(i)*vlx(i,2)
       vdef(i,2)=py2(i)*vly(i,1)+py3(i)*vly(i,2)
       vdef(i,3)=px2(i)*vly(i,1)+px3(i)*vly(i,2)
     .          +py2(i)*vlx(i,1)+py3(i)*vlx(i,2)
       vdef(i,4)=zero
       vdef(i,5)=zero
      ENDDO
C
      DO i=jft,jlt
       a5 = pb3(i,1,4)*vz5(i,1)+pb3(i,1,5)*vz5(i,2)
     .     +pb3(i,1,6)*vz5(i,3)
       a6 = pb3(i,2,4)*vz6(i,1)+pb3(i,2,5)*vz6(i,2)
     .     +pb3(i,2,6)*vz6(i,3)
       a4 = pb3(i,3,4)*vz4(i,1)+pb3(i,3,5)*vz4(i,2)
     .     +pb3(i,3,6)*vz4(i,3)
       vdef(i,6)=pb1(i,1,1)*vlz(i,1)+pb1(i,1,2)*vlz(i,2)
     .          +pb1(i,1,3)*vlz(i,3)+pb2(i,1,1)*vlz(i,4)
     .          +pb2(i,1,2)*vlz(i,5)+pb2(i,1,3)*vlz(i,6)
     .          +pb3(i,1,1)*a5+pb3(i,1,2)*a6+pb3(i,1,3)*a4
       vdef(i,7)=pb1(i,2,1)*vlz(i,1)+pb1(i,2,2)*vlz(i,2)
     .          +pb1(i,2,3)*vlz(i,3)+pb2(i,2,1)*vlz(i,4)
     .          +pb2(i,2,2)*vlz(i,5)+pb2(i,2,3)*vlz(i,6)
     .          +pb3(i,2,1)*a5+pb3(i,2,2)*a6+pb3(i,2,3)*a4
       vdef(i,8)=pb1(i,3,1)*vlz(i,1)+pb1(i,3,2)*vlz(i,2)
     .          +pb1(i,3,3)*vlz(i,3)+pb2(i,3,1)*vlz(i,4)
     .          +pb2(i,3,2)*vlz(i,5)+pb2(i,3,3)*vlz(i,6)
     .          +pb3(i,3,1)*a5+pb3(i,3,2)*a6+pb3(i,3,3)*a4
      ENDDO
C
      RETURN
      END