hm__read__fail__chang_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_chang       ../starter/source/materials/fail/changchang/hm_read_fail_chang.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       -----------------------------------------------------

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

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

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

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


      SUBROUTINE hm_read_fail_chang(FAIL ,

     .           MAT_ID   ,FAIL_ID  ,IRUPT    ,

     .           LSUBMODEL,UNITAB   ,FAIL_TAG )

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

C   ROUTINE DESCRIPTION :

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

C   READ FAILURE CHANG

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

      USE elbuftag_mod

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      "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 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

      TYPE (FAIL_TAG_)    ,INTENT(INOUT) :: FAIL_TAG   !< failure model tag for buffer allocation

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

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

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

      INTEGER :: IFAIL_SH,FAILIP

      my_real :: SIGT1,SIGT2,SIGT12,SIGC1,SIGC2,BETA,TMAX,PTHKF

      LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED

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

C   S o u r c e   L i n e s

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

      is_encrypted = .false.

      is_available = .false.


      CALL hm_option_is_encrypted(is_encrypted)


      CALL hm_get_floatv    ('Sigma_1t'   ,sigt1      ,is_available,lsubmodel,unitab)

      CALL hm_get_floatv    ('Sigma_2t'   ,sigt2      ,is_available,lsubmodel,unitab)

      CALL hm_get_floatv    ('Sigma_12'   ,sigt12     ,is_available,lsubmodel,unitab)

      CALL hm_get_floatv    ('Sigma_1c'   ,sigc1      ,is_available,lsubmodel,unitab)

      CALL hm_get_floatv    ('Sigma_2c'   ,sigc2      ,is_available,lsubmodel,unitab)


      CALL hm_get_floatv    ('Beta'       ,beta       ,is_available,lsubmodel,unitab)

      CALL hm_get_floatv    ('Tau_max'    ,tmax       ,is_available,lsubmodel,unitab)

      CALL hm_get_intv      ('Ifail_sh'   ,ifail_sh     ,is_available,lsubmodel)

      CALL hm_get_intv      ('Failip'     ,failip     ,is_available,lsubmodel)

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


      IF (sigt1  == zero) sigt1  = infinity

      IF (sigt2  == zero) sigt2  = infinity

      IF (sigt12 == zero) sigt12 = infinity

      IF (sigc1  == zero) sigc1  = infinity

      IF (sigc2  == zero) sigc2  = infinity

      IF (tmax   == zero) tmax   = infinity

      IF (ifail_sh == 0)    ifail_sh = 1

      IF (failip == 0)    failip = 1

c

      IF (ifail_sh == 1) THEN        ! matrix or fiber failure

        pthkf = em06

      ELSEIF (ifail_sh == 2) THEN    ! matrix or fiber failure

        pthkf = one

      ELSEIF (ifail_sh == 3) THEN    ! fiber only failure

        pthkf = em06

      ELSEIF (ifail_sh == 4) THEN    ! fiber only failure

        pthkf = one

      END IF

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

      fail%KEYWORD = 'CHANG'

      fail%IRUPT   = irupt

      fail%FAIL_ID = fail_id

      fail%NUPARAM = 9

      fail%NIPARAM = 0

      fail%NUVAR   = 1

      fail%NFUNC   = 0

      fail%NTABLE  = 0

      fail%NMOD    = 4

      fail%PTHK    = pthkf

c

      ALLOCATE (fail%UPARAM(fail%NUPARAM))

      ALLOCATE (fail%IPARAM(fail%NIPARAM))

      ALLOCATE (fail%IFUNC (fail%NFUNC))

      ALLOCATE (fail%TABLE (fail%NTABLE))

c

      ! Modes of failure

      fail_tag%LF_DAMMX = fail_tag%LF_DAMMX + fail%NMOD

      ALLOCATE (fail%MODE(fail%NMOD))

      fail%MODE(1) = "Failure index tension dir. 1"

      fail%MODE(2) = "Failure index compression dir. 1"

      fail%MODE(3) = "Failure index tension dir. 2"

      fail%MODE(4) = "Failure index compression dir. 2"

c

      fail%UPARAM(1) = sigt1

      fail%UPARAM(2) = sigt2

      fail%UPARAM(3) = sigt12

      fail%UPARAM(4) = sigc1

      fail%UPARAM(5) = sigc2

      fail%UPARAM(6) = beta

      fail%UPARAM(7) = tmax

      fail%UPARAM(8) = ifail_sh

      fail%UPARAM(9) = failip

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

      IF(is_encrypted)THEN


      ELSE

       WRITE(iout, 1000)sigt1,sigt2,sigt12,sigc1,sigc2,beta,tmax

       IF(ifail_sh==1)THEN

        WRITE(iout, 1100)

       ELSEIF(ifail_sh==2)THEN

        WRITE(iout, 1200)

       ENDIF

       WRITE(iout,1300) failip

      ENDIF

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

 1000 FORMAT(

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

     & 5x,'              CHANG CHANG FAILURE MODEL                  ',/,

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

     & 5x,'LONGITUDINAL TENSILE STRENGTH SIGMA1_T . . . . . . . . .=',1pg20.13/

     & 5x,'TRANSVERSAL  TENSILE STRENGTH SIGMA2_T . . . . . . . . .=',1pg20.13/

     & 5x,'SHEAR STRENGTH SIGMA_12. . . . . . . . . . . . . . . . .=',1pg20.13/

     & 5x,'LONGITUDINAL COMPRESSIVE STRENGTH SIGMA1_C . . . . . . .=',1pg20.13/

     & 5x,'TRANSVERSAL  COMPRESSIVE STRENGTH SIGMA2_C . . . . . . .=',1pg20.13/

     & 5x,'SHEAR SCALING FACTOR BETA. . . . . . . . . . . . . . . .=',1pg20.13/

     & 5x,'RELAXATION TIME TAU_MAX. . . . . . . . . . . . . . . . .=',1pg20.13/)

 1100 FORMAT(

     & 5x,'SHELL IS DELETED IF CRITERIA IS REACHED FOR ONE LAYER    ')

 1200 FORMAT(

     & 5x,'SHELL TENSOR STRESS IS SET TO ZERO AND SHELL IS DELETED  ',/

     & 5x,'IF CRITERIA IS REACHED FOR ALL LAYERS' )

 1300 FORMAT(

     & 5x,'NUMBER OF FAILED INTG. POINT PRIOR TO SOLID DELETION . .=',i5/)

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

      RETURN


      END

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

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

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

hm_read_fail_chang
subroutine hm_read_fail_chang(fail, mat_id, fail_id, irupt, lsubmodel, unitab, fail_tag)
Definition hm_read_fail_chang.F:39

elbuftag_mod
Definition elbuftag_mod.F:227

hm_option_read_mod
Definition hm_option_read_mod.F:408

submodel_mod
Definition submodel_mod.F:712

unitab_mod
Definition unitab_mod.F:513