ini_outmax.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/.
!||====================================================================
!||    ini_tmax       ../engine/source/output/ini_outmax.F
!||--- called by ------------------------------------------------------
!||    resol_init     ../engine/source/engine/resol_init.F
!||--- calls      -----------------------------------------------------
!||    ini_tmnorm2    ../engine/source/output/outmaxsubr.F
!||    my_barrier     ../engine/source/system/machine.F
!||    sig2d2pin      ../engine/source/output/outmaxsubr.F
!||    sig3d2pin      ../engine/source/output/outmaxsubr.F
!||    upd_tmax       ../engine/source/output/upd_outmax.F
!||--- uses       -----------------------------------------------------
!||    elbufdef_mod   ../common_source/modules/mat_elem/elbufdef_mod.F90
!||    element_mod    ../common_source/modules/elements/element_mod.F90
!||    outmax_mod     ../common_source/modules/outmax_mod.F
!||    stack_mod      ../engine/share/modules/stack_mod.F
!||====================================================================
      SUBROUTINE ini_tmax(ELBUF_TAB,IPARG   ,GEO     ,PM   ,
     .           IXS  ,IXS10   ,IXS16   ,IXS20   ,IXQ     ,
     .           IXC  ,IXTG   ,IXT    ,IXP     ,IXR     ,
     .           X  ,D       ,V       ,IAD_ELEM,FR_ELEM ,
     .             WEIGHT ,IPM   ,IGEO    ,STACK   ,ITASK   )
C-----------------------------------------------
C   M o d u l e s
C----------------------------------------------- 
      USE elbufdef_mod            
      USE outmax_mod
      USE stack_mod
      use element_mod , only : nixs,nixq,nixc,nixp,nixr,nixt,nixtg
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      "com01_c.inc"
#include      "com04_c.inc"
#include      "com08_c.inc"
#include      "param_c.inc"
#include      "tabsiz_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN):: ITASK
      INTEGER, DIMENSION(NPARG,NGROUP) ,INTENT(IN):: IPARG
      INTEGER, DIMENSION(NIXC,NUMELC) ,INTENT(IN):: IXC
      INTEGER, DIMENSION(NIXTG,NUMELTG) ,INTENT(IN):: IXTG
      INTEGER, DIMENSION(NIXS,NUMELS) ,INTENT(IN):: IXS
      INTEGER, DIMENSION(NIXQ,NUMELQ) ,INTENT(IN):: IXQ
      INTEGER, DIMENSION(NIXT,NUMELT) ,INTENT(IN):: IXT
      INTEGER, DIMENSION(NIXP,NUMELP) ,INTENT(IN):: IXP
      INTEGER, DIMENSION(NIXR,NUMELR) ,INTENT(IN):: IXR
      INTEGER, DIMENSION(6,NUMELS10) ,INTENT(IN):: IXS10
      INTEGER, DIMENSION(8,NUMELS16) ,INTENT(IN):: IXS16
      INTEGER, DIMENSION(12,NUMELS20) ,INTENT(IN):: IXS20
      INTEGER, DIMENSION(2,NSPMD+1)   ,INTENT(IN):: IAD_ELEM
      INTEGER, DIMENSION(SFR_ELEM)     ,INTENT(IN):: FR_ELEM
      INTEGER, DIMENSION(NUMNOD)     ,INTENT(IN):: WEIGHT
      INTEGER, DIMENSION(NPROPMI,NUMMAT) ,INTENT(IN):: IPM
      INTEGER, DIMENSION(NPROPGI,NUMGEO) ,INTENT(IN):: IGEO
      my_real, DIMENSION(NPROPG,NUMGEO) ,INTENT(IN):: geo
      my_real, DIMENSION(NPROPM,NUMMAT) ,INTENT(IN):: pm
      my_real, DIMENSION(3,NUMNOD) ,INTENT(IN):: x,d,v
      TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
      TYPE (STACK_PLY) :: STACK
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NG, ITY, NEL
      TYPE(G_BUFEL_)  ,POINTER :: GBUF     
 
 
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------     
c-----doing TT=0 initialization instead of Starter; elementary first      
       IF(itask == 0) nmax_e = 0
       IF (tt==zero) THEN
         IF (itask==0) THEN
           DO ng=1,ngroup
             ity=iparg(5,ng)
             nel=iparg(2,ng)
             IF (ity /= 1 .AND. ity /= 3 .AND. ity /= 7) cycle
             gbuf => elbuf_tab(ng)%GBUF
             nmax_e = nmax_e + gbuf%G_TM_YIELD
             nmax_e = nmax_e + gbuf%G_TM_SEQ  
             nmax_e = nmax_e + gbuf%G_TM_EINT 
             nmax_e = nmax_e + gbuf%G_TM_DMG  
             nmax_e = nmax_e + gbuf%G_TM_SIG  
             nmax_e = nmax_e + gbuf%G_TM_STRA 
           END DO
           nmax_n = lmax_dis+lmax_vel+lmax_nsig+lmax_nstra
           IF (nmax_n>0) THEN
             IF (lmax_dis>0) tm_pdis = zero
             IF (lmax_vel>0) tm_pvel = zero
             IF (lmax_nsig>0) THEN
               tm_pnsig(1:numnod) = zero
               tm_pnsig(numnod+1:2*numnod) = ep30
             END IF
             IF (lmax_nstra>0) THEN
               tm_pnstra(1:numnod) = zero
               tm_pnstra(numnod+1:2*numnod) = ep30
             END IF
           END IF
         END IF !(ITASK==0) THEN
C-----------------         
         CALL my_barrier()
C-----------------         
         CALL upd_tmax(elbuf_tab,iparg   ,geo     ,pm   ,
     .           ixs  ,ixs10   ,ixs16   ,ixs20   ,ixq     ,
     .           ixc  ,ixtg   ,ixt    ,ixp     ,ixr     ,
     .           x  ,d       ,v       ,iad_elem,fr_elem ,
     .             weight ,ipm   ,igeo    ,stack   ,itask   )
C
       ELSEIF (itask==0) THEN
         DO ng=1,ngroup
           ity=iparg(5,ng)
           nel=iparg(2,ng)
           IF (ity /= 1 .AND. ity /= 3 .AND. ity /= 7) cycle
           gbuf => elbuf_tab(ng)%GBUF
           nmax_e = nmax_e + gbuf%G_TM_YIELD
           nmax_e = nmax_e + gbuf%G_TM_SEQ  
           nmax_e = nmax_e + gbuf%G_TM_EINT 
           nmax_e = nmax_e + gbuf%G_TM_DMG  
           nmax_e = nmax_e + gbuf%G_TM_SIG  
           nmax_e = nmax_e + gbuf%G_TM_STRA 
          SELECT CASE (ity)
            CASE(1)
               IF (gbuf%G_TM_SIG>0) THEN
                 CALL sig3d2pin(gbuf%TM_SIG1,gbuf%TM_SIG3,gbuf%TM_PSIG,nel)
               END IF
               IF (gbuf%G_TM_STRA>0) THEN
                 CALL sig3d2pin(gbuf%TM_STRA1,gbuf%TM_STRA3,gbuf%TM_PSTRA,nel)
               END IF
            CASE(3,7)
              IF (gbuf%G_TM_SIG>0) THEN
                CALL sig2d2pin(gbuf%TM_SIG1,gbuf%TM_SIG3,gbuf%TM_PSIG,nel)
              END IF
              IF (gbuf%G_TM_STRA>0) THEN
                CALL sig2d2pin(gbuf%TM_STRA1,gbuf%TM_STRA3,gbuf%TM_PSTRA,nel)
              END IF
          END SELECT
         END DO
         nmax_n = lmax_dis+lmax_vel+lmax_nsig+lmax_nstra
         IF (nmax_n>0) THEN
           IF (lmax_dis>0) CALL ini_tmnorm2(tm_dis,tm_pdis,numnod)
           IF (lmax_vel>0) CALL ini_tmnorm2(tm_vel,tm_pvel,numnod)
           IF (lmax_nsig>0)CALL sig3d2pin(tm_nsig1,tm_nsig3,tm_pnsig,numnod)
           IF (lmax_nstra>0)CALL sig3d2pin(tm_nstra1,tm_nstra3,tm_pnstra,numnod)
         END IF
       END IF !(TT==ZERO) THEN
 
      RETURN
      END