s8ejacit3.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/.
!||====================================================================
!||    s8ejacit3   ../engine/source/elements/solid/solide8e/s8ejacit3.F
!||--- called by ------------------------------------------------------
!||    s8eforc3    ../engine/source/elements/solid/solide8e/s8eforc3.F
!||====================================================================
      SUBROUTINE s8ejacit3(
     1   VOLDP,   KSI,     ETA,     ZETA,
     2   WI,      AJP1,    AJP2,    AJP3,
     3   AJP4,    AJP5,    AJP6,    AJP7,
     4   AJP8,    AJP9,    AJI1,    AJI2,
     5   AJI3,    AJI4,    AJI5,    AJI6,
     6   AJI7,    AJI8,    AJI9,    JAC_I,
     7   NEL,     ISMSTR)
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.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   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: ISMSTR
C     REAL
      my_real
     .   KSI,ETA,ZETA,WI,
     .   AJP1(*),AJP2(*),AJP3(*),
     .   AJP4(*),AJP5(*),AJP6(*),
     .   ajp7(*),ajp8(*),ajp9(*),
     .   aji1(mvsiz), aji2(mvsiz), aji3(mvsiz),
     .   aji4(mvsiz), aji5(mvsiz), aji6(mvsiz),
     .   aji7(mvsiz), aji8(mvsiz), aji9(mvsiz),
     .   jac_i(10,mvsiz)
      DOUBLE PRECISION 
     .   VOLDP(MVSIZ)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J ,ICOR
C     REAL
C                                                                     12
      my_real
     .   det(mvsiz) ,dett(mvsiz) ,  
     .   aj1(mvsiz),aj2(mvsiz),aj3(mvsiz),
     .   aj4(mvsiz),aj5(mvsiz),aj6(mvsiz),
     .   aj7(mvsiz),aj8(mvsiz),aj9(mvsiz),
     .   jac_59_68(mvsiz), jac_67_49(mvsiz), jac_48_57(mvsiz),
     .   jac_38_29(mvsiz), jac_19_37(mvsiz), jac_27_18(mvsiz),
     .   jac_26_35(mvsiz), jac_34_16(mvsiz), jac_15_24(mvsiz),
     .   aj12(mvsiz), aj45(mvsiz), aj78(mvsiz)
      DOUBLE PRECISION 
     .   DETDP
C-----------------------------------------------
      IF (ISMSTR==11) then
       aj1(1:nel)=ajp1(1:nel)
       aj2(1:nel)=ajp2(1:nel)
       aj3(1:nel)=ajp3(1:nel)
       aj4(1:nel)=ajp4(1:nel)
       aj5(1:nel)=ajp5(1:nel)
       aj6(1:nel)=ajp6(1:nel)
       aj7(1:nel)=ajp7(1:nel)
       aj8(1:nel)=ajp8(1:nel)
       aj9(1:nel)=ajp9(1:nel)
 
C     JACOBIAN
C
      DO i=1,nel
       jac_59_68(i)=aj5(i)*aj9(i)-aj6(i)*aj8(i)
       jac_67_49(i)=aj6(i)*aj7(i)-aj4(i)*aj9(i)
       jac_38_29(i)=(-aj2(i)*aj9(i)+aj3(i)*aj8(i))
       jac_19_37(i)=( aj1(i)*aj9(i)-aj3(i)*aj7(i))
       jac_27_18(i)=(-aj1(i)*aj8(i)+aj2(i)*aj7(i))
       jac_26_35(i)=( aj2(i)*aj6(i)-aj3(i)*aj5(i))
       jac_34_16(i)=(-aj1(i)*aj6(i)+aj3(i)*aj4(i))
       jac_15_24(i)=( aj1(i)*aj5(i)-aj2(i)*aj4(i))
       jac_48_57(i)=aj4(i)*aj8(i)-aj5(i)*aj7(i)
      ENDDO
C
      DO i=1,nel
      detdp=one_over_512*(aj1(i)*jac_59_68(i)+aj2(i)*jac_67_49(i)+aj3(i)*jac_48_57(i))
      det(i) = detdp
      voldp(i)= wi*detdp
      dett(i)=one_over_512/det(i)
      ENDDO
C
      DO i=1,nel
       aji1(i)=dett(i)*jac_59_68(i)
       aji4(i)=dett(i)*jac_67_49(i)
       aji7(i)=dett(i)*jac_48_57(i)
       aji2(i)=dett(i)*jac_38_29(i)
       aji5(i)=dett(i)*jac_19_37(i)
       aji8(i)=dett(i)*jac_27_18(i)
       aji3(i)=dett(i)*jac_26_35(i)
       aji6(i)=dett(i)*jac_34_16(i)
       aji9(i)=dett(i)*jac_15_24(i)
      ENDDO
C---------don't forget this for Ismstr11 and R,S,T
c      DO I=1,NEL
c       IF (IDEG(I)==0) DELTAX(I)=MIN(DELTAX(I),HUNDRED28*DET(I)*SMAX(I))
c      ENDDO
C------ ISMSTR =10
      ELSE
      DO i=1,nel
       aji1(i)=jac_i(1,i)
       aji4(i)=jac_i(4,i)
       aji7(i)=jac_i(7,i)
       aji2(i)=jac_i(2,i)
       aji5(i)=jac_i(5,i)
       aji8(i)=jac_i(8,i)
       aji3(i)=jac_i(3,i)
       aji6(i)=jac_i(6,i)
       aji9(i)=jac_i(9,i)
       voldp(i) =jac_i(10,i)
      END DO
      END IF !(ISMSTR==11)
C
      RETURN
      END