s4coork.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/.
!||====================================================================
!||    s4coork       ../engine/source/elements/solid/solide4/s4coork.F
!||--- called by ------------------------------------------------------
!||    s4ke3         ../engine/source/elements/solid/solide4/s4ke3.F
!||--- uses       -----------------------------------------------------
!||    element_mod   ../common_source/modules/elements/element_mod.F90
!||====================================================================
      SUBROUTINE s4coork(
     1   X,       IXS,     X1,      X2,
     2   X3,      X4,      Y1,      Y2,
     3   Y3,      Y4,      Z1,      Z2,
     4   Z3,      Z4,      OFFG,    OFF,
     5   SAV,     NC1,     NC2,     NC3,
     6   NC4,     NGL,     MXT,     NGEO,
     7   K11,     K12,     K13,     K14,
     8   K22,     K23,     K24,     K33,
     9   K34,     K44,     NEL,     ISMSTR)
      use element_mod , only : nixs
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "scr18_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: ISMSTR
      INTEGER NEL
      my_real
     .   X(3,*),
     .   X1(*), X2(*), X3(*), X4(*),  
     .   Y1(*), Y2(*), Y3(*), Y4(*),  
     .   Z1(*), Z2(*), Z3(*), Z4(*), 
     .   OFFG(*), OFF(*)
      DOUBLE PRECISION
     .  sav(nel,9)
      my_real
     .   k11(9,*)   ,k12(9,*)   ,k13(9,*)  ,k14(9,*)    ,k22(9,*)  ,
     .   k23(9,*)   ,k24(9,*)   ,k33(9,*)   ,k34(9,*)   ,k44(9,*)     
      INTEGER NC1(*), NC2(*), NC3(*), NC4(*),MXT(*), NGL(*),NGEO(*)
      INTEGER IXS(NIXS,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,MXT_1
C     REAL
C-----------------------------------------------
C
      mxt_1 = ixs(1,1)
      DO i=1,nel
      ngeo(i)=ixs(10,i)
      ngl(i)=ixs(11,i)
      mxt(i)=mxt_1
      nc1(i)=ixs(2,i)
      nc2(i)=ixs(4,i)
      nc3(i)=ixs(7,i)
      nc4(i)=ixs(6,i)
      ENDDO
C----------------------------
C     NODAL COORDINATES     |
C----------------------------
      DO i=1,nel
        x1(i)=x(1,nc1(i))
        y1(i)=x(2,nc1(i))
        z1(i)=x(3,nc1(i))
        x2(i)=x(1,nc2(i))
        y2(i)=x(2,nc2(i))
        z2(i)=x(3,nc2(i))
        x3(i)=x(1,nc3(i))
        y3(i)=x(2,nc3(i))
        z3(i)=x(3,nc3(i))
        x4(i)=x(1,nc4(i))
        y4(i)=x(2,nc4(i))
        z4(i)=x(3,nc4(i))
        off(i) = min(one,abs(offg(i)))
      ENDDO
C-----------
      IF(ismstr==1.OR.(ismstr==2.AND.idtmin(1)==3))THEN
       DO i=1,nel
       IF(abs(offg(i))>one)THEN
        x1(i)=sav(i,1)
        y1(i)=sav(i,2)
        z1(i)=sav(i,3)
        x2(i)=sav(i,4)
        y2(i)=sav(i,5)
        z2(i)=sav(i,6)
        x3(i)=sav(i,7)
        y3(i)=sav(i,8)
        z3(i)=sav(i,9)
        x4(i)=zero
        y4(i)=zero
        z4(i)=zero
        off(i) = abs(offg(i))-one
       ELSE
        off(i) = offg(i)
       ENDIF
       ENDDO
C
      ENDIF
C-----------
      DO j=1,9      
      DO i=1,nel
        k11(j,i)=zero
        k12(j,i)=zero
        k13(j,i)=zero
        k14(j,i)=zero
        k22(j,i)=zero
        k23(j,i)=zero
        k24(j,i)=zero
        k33(j,i)=zero
        k34(j,i)=zero
        k44(j,i)=zero
      ENDDO
      ENDDO
C-----------
      RETURN
      END