force_imp.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/.
! wrapper for implicit calls , cannot be the included in the module
! because the calls does not respect Fortran standard and will cause
! compiler errors
!||====================================================================
!||    force_imp       ../engine/source/loads/general/force_imp.F
!||--- called by ------------------------------------------------------
!||    dyna_ina        ../engine/source/implicit/imp_dyna.F
!||    dyna_wex        ../engine/source/implicit/imp_dyna.F
!||    imp_chkm        ../engine/source/implicit/imp_solv.F
!||    imp_solv        ../engine/source/implicit/imp_solv.F
!||--- calls      -----------------------------------------------------
!||    force           ../engine/source/loads/general/force.F90
!||--- uses       -----------------------------------------------------
!||    force_mod       ../engine/source/loads/general/force.F90
!||    h3d_mod         ../engine/share/modules/h3d_mod.F
!||    pinchtype_mod   ../common_source/modules/pinchtype_mod.F
!||    sensor_mod      ../common_source/modules/sensor_mod.F90
!||    skew_mod        ../common_source/modules/skew_mod.F90
!||    th_surf_mod     ../common_source/modules/interfaces/th_surf_mod.F
!||====================================================================
      SUBROUTINE force_imp(IB      ,FAC     ,SNPC    ,NPC        ,TF         ,
     2                     A       ,V       ,X       ,SKEWS      ,
     3                     AR      ,VR      ,NSENSOR ,SENSOR_TAB ,WFEXC      ,
     4                     IADC    ,FSKY    ,FEXT    ,H3D_DATA   ,CPTREAC    ,
     5                     FTHREAC ,NODREAC ,TH_SURF ,
     6                     DPL0CLD ,VEL0CLD ,D       ,DR         ,NCONLD     ,
     7                     NUMNOD  ,NFUNCT  ,STF     ,WFEXT)
C-----------------------------------------------
C   M o d u l e s
C----------------------------------------------- 
      USE h3d_mod
      USE pinchtype_mod 
      USE sensor_mod
      USE th_surf_mod , ONLY : th_surf_
      USE force_mod   , ONLY : force
      USE skew_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
!#include      "com04_c.inc"
#include      "com06_c.inc"
#include      "com08_c.inc"
#include      "impl1_c.inc"
#include      "scr05_c.inc"
#include      "scr14_c.inc"
#include      "scr16_c.inc"
#include      "parit_c.inc"
C-----------------------------------------------,
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: SNPC
      INTEGER ,INTENT(IN) :: NSENSOR,NCONLD,NUMNOD,NFUNCT
      INTEGER ,INTENT(IN) :: STF
      INTEGER NPC(SNPC),CPTREAC,NODREAC(*)
      INTEGER IB(NIBCLD,*)
      INTEGER IADC(4,*)
      my_real fac(lfaccld,*),tf(*),a(3,*),v(3,*),ar(3,*),vr(3,*),x(3,*),wfexc,fsky(8,lsky),fext(3,*),fthreac(6,*)
      TYPE(h3d_database) :: H3D_DATA
      TYPE (SENSOR_STR_) ,DIMENSION(NSENSOR) :: SENSOR_TAB
      TYPE(skew_),INTENT(INOUT) :: SKEWS
      TYPE (TH_SURF_) , INTENT(INOUT) :: TH_SURF
      my_real, INTENT(IN) ::  dpl0cld(6,nconld),vel0cld(6,nconld)
      my_real, INTENT(IN) ::  d(3,numnod), dr(3,numnod)
      DOUBLE PRECISION,INTENT(INOUT) :: WFEXT
C-----------------------------------------------
      call force( nibcld  ,ib      ,lfaccld    ,fac      ,snpc      ,
     2            npc     ,stf     ,tf         ,a        ,v         ,
     3            x       ,skews   ,ar         ,
     4            vr      ,nsensor ,sensor_tab ,wfexc    ,iadc      ,
     5            lsky    ,fsky    ,fext       ,h3d_data ,cptreac   ,
     6            fthreac ,nodreac ,th_surf    ,
     7            dpl0cld ,vel0cld ,d          ,dr       ,nconld    ,
     8            numnod  ,nfunct  , anim_v    ,outp_v   ,
     9            iparit  ,tt      ,dt1        ,n2d      ,wfext     ,
     a            impl_s) 
 
      END SUBROUTINE