sfor_n2stria.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/.
!||====================================================================
!||    sfor_n2stria   ../engine/source/elements/solid/solide4/sfor_n2stria.F
!||--- called by ------------------------------------------------------
!||    s4for_distor   ../engine/source/elements/solid/solide4/s4for_distor.F
!||    s6for_distor   ../engine/source/elements/thickshell/solide6c/s6for_distor.F90
!||--- calls      -----------------------------------------------------
!||    sfor_n2s3      ../engine/source/elements/solid/solide4/sfor_n2s3.F
!||====================================================================
      SUBROUTINE sfor_n2stria(XC,      YC,     ZC, 
     .                        X1,      X2,     X3,
     .                        Y1,      Y2,     Y3,
     .                        Z1,      Z2,     Z3,
     .                       VCX,     VCY,    VCZ,
     .                       VX1,     VX2,    VX3,
     .                       VY1,     VY2,    VY3,
     .                       VZ1,     VZ2,    VZ3,
     .                     FOR_T1, FOR_T2, FOR_T3,
     .                     FORC_N,  STIF0,   STIF, 
     .                     FQMAX , PENMIN, PENREF,
     .                         LL,  IFCTL,   NEL ,
     .                   E_DISTOR,   DT1)
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, INTENT (IN)  :: NEL
      INTEGER, INTENT (OUT) :: IFCTL
      my_real,                   INTENT (IN) :: FQMAX,DT1
      my_real, DIMENSION(MVSIZ), INTENT (IN) :: STIF0,LL,
     .                        X1,      X2,     X3,
     .                        Y1,      Y2,     Y3,
     .                        Z1,      Z2,     Z3,
     .                        XC,      YC,     ZC,
     .                       VCX,     VCY,    VCZ,
     .                       VX1,     VX2,    VX3,
     .                       VY1,     VY2,    VY3,
     .                       VZ1,     VZ2,    VZ3,
     .                     penmin, penref
      my_real, DIMENSION(MVSIZ), INTENT (INOUT) :: stif
      my_real, DIMENSION(MVSIZ,3), INTENT (INOUT) :: forc_n,
     .                           for_t1, for_t2, for_t3
      my_real, DIMENSION(NEL), INTENT (INOUT) :: e_distor
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,IFCTL1,IFC2(MVSIZ)
C                                                                     12
      my_real
     .   fx,fy,fz,fkt(mvsiz)
C----------------------------
        ifc2(1:nel) = 1
        fkt(1:nel) = one
        DO i=1,nel
          IF (stif0(i)>zero) cycle
          ifc2(i) = 0
        ENDDO
         CALL sfor_n2s3(xc ,     yc ,    zc , forc_n,     
     .                   x1,      y1,     z1, for_t1, 
     .                   x2,      y2,     z2, for_t2,
     .                   x3,      y3,     z3, for_t3,
     .                stif0,    ifc2,  ifctl, penmin,
     .               penref,     fkt,    ll ,  fqmax,
     .                  vx1,     vx2,    vx3,  vcx  ,
     .                  vy1,     vy2,    vy3,  vcy  ,
     .                  vz1,     vz2,    vz3,  vcz  ,
     .                  nel,e_distor,    dt1)
        DO i=1,nel
          IF (stif0(i)==zero) cycle
          stif(i) =max(stif(i),fkt(i)*stif0(i))
        ENDDO
C
      RETURN
      END