s8sav12.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/.
!||====================================================================
!||    s8sav12    ../engine/source/elements/solid/solide/s8sav12.F
!||--- called by ------------------------------------------------------
!||    s8eforc3   ../engine/source/elements/solid/solide8e/s8eforc3.F
!||    s8zforc3   ../engine/source/elements/solid/solide8z/s8zforc3.F
!||    sforc3     ../engine/source/elements/solid/solide/sforc3.F
!||    szforc3    ../engine/source/elements/solid/solidez/szforc3.f
!||====================================================================
      SUBROUTINE s8sav12(
     1   OFFG,    OFFG0,   SAV,     X,
     2   XDP,     NC1,     NC2,     NC3,
     3   NC4,     NC5,     NC6,     NC7,
     4   NC8,     NEL,     JCVT)
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   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr05_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: JCVT
      INTEGER NEL
C     REAL
      INTEGER 
     .        nc1(*),nc2(*),nc3(*),nc4(*),nc5(*),nc6(*),nc7(*),nc8(*)
      my_real
     .   offg(*),offg0(*),x(3,*)
      DOUBLE PRECISION 
     .   XDP(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I ,IH    
      DOUBLE PRECISION 
     .   XD1(MVSIZ), XD2(MVSIZ), XD3(MVSIZ), XD4(MVSIZ),
     .   XD5(MVSIZ), XD6(MVSIZ), XD7(MVSIZ), XD8(MVSIZ),
     .   yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz),
     .   yd5(mvsiz), yd6(mvsiz), yd7(mvsiz), yd8(mvsiz),
     .   zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz),
     .   zd5(mvsiz), zd6(mvsiz), zd7(mvsiz), zd8(mvsiz), sav(nel,21)
C=======================================================================
C-----Always in global system and fixed then (not like Ismstr=11
C-----for JCVT>0 Fixed SAV is rotated w/ local system not like Ismstr=1
        ih=0
        DO i=1,nel
         IF(offg(i)/=offg0(i).AND.abs(offg(i)) > one ) ih=1
          END DO
          IF (ih==0) RETURN
C
      IF(iresp==1)THEN 
        DO i=1,nel
          xd1(i)=xdp(1,nc1(i))  
          yd1(i)=xdp(2,nc1(i))  
          zd1(i)=xdp(3,nc1(i))  
          xd2(i)=xdp(1,nc2(i))  
          yd2(i)=xdp(2,nc2(i))  
          zd2(i)=xdp(3,nc2(i))  
          xd3(i)=xdp(1,nc3(i))  
          yd3(i)=xdp(2,nc3(i))  
          zd3(i)=xdp(3,nc3(i))  
          xd4(i)=xdp(1,nc4(i))  
          yd4(i)=xdp(2,nc4(i))  
          zd4(i)=xdp(3,nc4(i))  
          xd5(i)=xdp(1,nc5(i))  
          yd5(i)=xdp(2,nc5(i))  
          zd5(i)=xdp(3,nc5(i))  
          xd6(i)=xdp(1,nc6(i))  
          yd6(i)=xdp(2,nc6(i))  
          zd6(i)=xdp(3,nc6(i))  
          xd7(i)=xdp(1,nc7(i))  
          yd7(i)=xdp(2,nc7(i))  
          zd7(i)=xdp(3,nc7(i))  
          xd8(i)=xdp(1,nc8(i))  
          yd8(i)=xdp(2,nc8(i))  
          zd8(i)=xdp(3,nc8(i))              
        ENDDO
      ELSE
        DO i=1,nel
          xd1(i)=x(1,nc1(i))  
          yd1(i)=x(2,nc1(i))  
          zd1(i)=x(3,nc1(i))  
          xd2(i)=x(1,nc2(i))  
          yd2(i)=x(2,nc2(i))  
          zd2(i)=x(3,nc2(i))  
          xd3(i)=x(1,nc3(i))  
          yd3(i)=x(2,nc3(i))  
          zd3(i)=x(3,nc3(i))  
          xd4(i)=x(1,nc4(i))  
          yd4(i)=x(2,nc4(i))  
          zd4(i)=x(3,nc4(i))  
          xd5(i)=x(1,nc5(i))  
          yd5(i)=x(2,nc5(i))  
          zd5(i)=x(3,nc5(i))  
          xd6(i)=x(1,nc6(i))  
          yd6(i)=x(2,nc6(i))  
          zd6(i)=x(3,nc6(i))  
          xd7(i)=x(1,nc7(i))  
          yd7(i)=x(2,nc7(i))  
          zd7(i)=x(3,nc7(i))  
          xd8(i)=x(1,nc8(i))  
          yd8(i)=x(2,nc8(i))  
          zd8(i)=x(3,nc8(i))  
        ENDDO      
      ENDIF
C   -----note: format is different to Ismstr=10,11
        DO i=1,nel
          IF(offg(i)/=offg0(i).AND.abs(offg(i)) > one )THEN
           sav(i,1)=xd1(i)-xd8(i)
           sav(i,2)=yd1(i)-yd8(i)
           sav(i,3)=zd1(i)-zd8(i)
           sav(i,4)=xd2(i)-xd8(i)
           sav(i,5)=yd2(i)-yd8(i)
           sav(i,6)=zd2(i)-zd8(i)
           sav(i,7)=xd3(i)-xd8(i)
           sav(i,8)=yd3(i)-yd8(i)
           sav(i,9)=zd3(i)-zd8(i)
           sav(i,10)=xd4(i)-xd8(i)
           sav(i,11)=yd4(i)-yd8(i)
           sav(i,12)=zd4(i)-zd8(i)
           sav(i,13)=xd5(i)-xd8(i)
           sav(i,14)=yd5(i)-yd8(i)
           sav(i,15)=zd5(i)-zd8(i)
           sav(i,16)=xd6(i)-xd8(i)
           sav(i,17)=yd6(i)-yd8(i)
           sav(i,18)=zd6(i)-zd8(i)
           sav(i,19)=xd7(i)-xd8(i)
           sav(i,20)=yd7(i)-yd8(i)
           sav(i,21)=zd7(i)-zd8(i)
          ENDIF
        ENDDO
C-----------
      RETURN
      END