hm__read__fail__visual_8F_source.html

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/.

!||====================================================================

!||    hm_read_fail_visual      ../starter/source/materials/fail/visual/hm_read_fail_visual.F

!||--- called by ------------------------------------------------------

!||    hm_read_fail             ../starter/source/materials/fail/hm_read_fail.F

!||--- calls      -----------------------------------------------------

!||    hm_get_floatv            ../starter/source/devtools/hm_reader/hm_get_floatv.F

!||    hm_get_intv              ../starter/source/devtools/hm_reader/hm_get_intv.F

!||    hm_option_is_encrypted   ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F

!||--- uses       -----------------------------------------------------

!||    hm_option_read_mod       ../starter/share/modules1/hm_option_read_mod.F

!||    submodel_mod             ../starter/share/modules1/submodel_mod.F

!||====================================================================


       SUBROUTINE hm_read_fail_visual(

     .            FAIL    ,MAT_ID   ,FAIL_ID  ,IRUPT    ,

     .            LSUBMODEL,UNITAB   )

C-----------------------------------------------

c    ROUTINE DESCRIPTION :

c    Read Visual failure estimation model parameters

C      Only visualization by damage of stress or strain is between to limited values

C            (no element deletion)

c     ! User variable # 1, to store the previous damage value

c     ! User variable # 2, to store the previous stress or strain value (for EMA filtering)

      ! What it should be is like this:

      ! Value = USER_INPUT * 2 * Pi * DT1

      ! Alpha = Value / (Value + 1)

      ! Actual_filtered_stress = Alpha * actual_Stress + (1-Alpha) * previous_filtered_stress

      ! UPARAM(4) = EMA  * TWO * 3.14159265

C-----------------------------------------------

C   M o d u l e s

C-----------------------------------------------

      USE fail_param_mod

      USE unitab_mod

      USE submodel_mod

      USE hm_option_read_mod

C---------------------------------------------------------------------

C I m p l i c i t T y p e s

C---------------------------------------------------------------------

#include  "implicit_f.inc"

#include  "units_c.inc"

C-----------------------------------------------

C D u m m y A r g u m e n t s

C-----------------------------------------------

      INTEGER            ,INTENT(IN) :: FAIL_ID         ! failure model ID

      INTEGER            ,INTENT(IN) :: MAT_ID          ! material law ID

      INTEGER            ,INTENT(IN) :: IRUPT           ! failure model type number

      TYPE(unit_type_)   ,INTENT(IN) :: UNITAB          ! table of input units

      TYPE(submodel_data),INTENT(IN) :: LSUBMODEL(*)    ! submodel table

      TYPE(fail_param_)  ,INTENT(INOUT) :: FAIL         ! failure model data structure

C-----------------------------------------------

C L o c a l V a r i a b l e s

C-----------------------------------------------

       INTEGER TYPE_MAX,F_FLAG,STRDEF,NUVAR

       my_real C_MIN , C_MAX,EMA,BUTTER

       LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED

C=======================================================================

      is_encrypted = .false.

      is_available = .false.


      CALL hm_option_is_encrypted(is_encrypted)

C--------------------------------------------------

       type_max  = 0

       c_min     = zero

       c_max     = zero

       butter    = zero

c!CARD1

      CALL hm_get_intv   ('Type'                ,type_max  ,is_available,lsubmodel)

      CALL hm_get_floatv ('C_min'               ,c_min     ,is_available,lsubmodel,unitab)

      CALL hm_get_floatv ('C_max'               ,c_max     ,is_available,lsubmodel,unitab)

      CALL hm_get_floatv ('Alpha_Exp_Avg_Filter',ema       ,is_available,lsubmodel,unitab)

      CALL hm_get_floatv ('f_cutoff'            ,ema       ,is_available,lsubmodel,unitab)

      CALL hm_get_intv   ('f_flag'              ,f_flag    ,is_available,lsubmodel)

      CALL hm_get_intv   ('Strdef'              ,strdef    ,is_available,lsubmodel)

c-----------------------------------------------------------------------

      IF (f_flag == 0 .OR. f_flag == 1 ) THEN

         f_flag = 1

         IF (ema == zero) ema = one

         nuvar  = 5

         ! User variable # 1: damage value

         ! User variable # 2: stress or strain value in the previous step

         ! User variable # 3-5: Exx, Exy, Exz for the failure of beam elements type3

      ELSE

         butter = ema

         ema    = one

         nuvar  = 14

         ! User variable # 1: damage value

         ! User variable # 2: stress or strain value in the previous step

         ! User variable # 3-11: Storage values for the Butterworth filter

         ! User variable # 12-14: Exx, Exy, Exz for the failure of beam elements type3

      ENDIF


      IF (type_max == 0) type_max = 1

      IF (strdef == 0)   strdef   = 1

c-----------------------------------------------------------------------

      fail%KEYWORD = 'VISUAL'

      fail%IRUPT   = irupt

      fail%FAIL_ID = fail_id

      fail%NUPARAM = 7

      fail%NIPARAM = 0

      fail%NUVAR   = nuvar

      fail%NFUNC   = 0

      fail%NTABLE  = 0

      fail%NMOD    = 0

c

      ALLOCATE (fail%UPARAM(fail%NUPARAM))

      ALLOCATE (fail%IPARAM(fail%NIPARAM))

      ALLOCATE (fail%IFUNC (fail%NFUNC))

      ALLOCATE (fail%TABLE (fail%NTABLE))

c

      fail%UPARAM(1) = real(type_max) + 1e-4

      fail%UPARAM(2) = c_min

      fail%UPARAM(3) = c_max

      fail%UPARAM(4) = ema

      fail%UPARAM(5) = butter

      fail%UPARAM(6) = real(f_flag) +0.00001

      fail%UPARAM(7) = real(strdef)

c-----------------------------------------------------------------------

      IF (is_encrypted) THEN

        WRITE (iout,'(5X,A,//)') 'CONFIDENTIAL DATA'

      ELSE

       IF (type_max == 1) THEN

         WRITE(iout, 1000) c_min, c_max

         IF (f_flag ==1 .AND. ema /= one) THEN

           WRITE(iout, 2010) ema

          ELSEIF(f_flag ==2)THEN

           WRITE(iout, 2020) butter

          ELSE

           WRITE(iout, 2030)

         ENDIF

       ENDIF

       IF (type_max == 2) THEN

         WRITE(iout, 2000) c_min, c_max, strdef

         IF (f_flag ==1 .AND. ema /= one) THEN

           WRITE(iout, 2010) ema

         ELSEIF(f_flag ==2)THEN

           WRITE(iout, 2020) butter

         ELSE

           WRITE(iout, 2030)

         ENDIF

       ENDIF

      ENDIF

c-----------------------------------------------------------------------

 1000 FORMAT(

     & 5x,'         VISUAL FAILURE MODEL           ',/

     & 5x,' -------------------------------------- ',/

     & 5x,' TYPE IS stress ',/

     & 5X,' minimum 1st principal stress. . . . . =',E12.4,/

     & 5X,' maximum 1st principal stress. . . . . =',E12.4,/)


 2000 FORMAT(

     & 5X,'         visual failure model           ',/

     & 5X,' -------------------------------------- ',/

     & 5X,' TYPE IS strain ',/

     & 5X,' minimum 1st principal strain. . . . . =',E12.4/

     & 5X,' maximum 1st principal strain. . . . . =',E12.4/

     & 5X,' failure strain input TYPE . . . . . . =',I10/

     & 5X,' =1 (strain measure defined in property)',/

     & 5X,' =2 (engineering strain) . . . . . . . .',/

     & 5X,' =3 (true strain). . . . . . . . . . . .')


 2010 FORMAT(

     & 5X,' exponential moving average filter used:',/

     & 5X,' f-coefficient . . . . . . . . . . . . =',E12.4,/)

 2020 FORMAT(

     & 5X,' 4-pole butterworth filter used:',/

     & 5X,' f-coefficient . . . . . . . . . . . . =',E12.4,/)


 2030 FORMAT(

     & 5X,' no filtering used:',/)


      RETURN


      END

butterworth
subroutine butterworth(dt, freq, x2, x1, x, fx2, fx1, fx)
Definition butterworth.F:31

hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:418

hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:528

hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:288

hm_read_fail_visual
subroutine hm_read_fail_visual(fail, mat_id, fail_id, irupt, lsubmodel, unitab)
Definition hm_read_fail_visual.F:38

hm_option_read_mod
Definition hm_option_read_mod.F:412

submodel_mod
Definition submodel_mod.F:717

unitab_mod
Definition unitab_mod.F:517

submodel_mod::submodel_data
Definition submodel_mod.F:728

unitab_mod::unit_type_
Definition unitab_mod.F:524