OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_fail_puck.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| hm_read_fail_puck ../starter/source/materials/fail/puck/hm_read_fail_puck.F
25!||--- called by ------------------------------------------------------
26!|| hm_read_fail ../starter/source/materials/fail/hm_read_fail.F
27!||--- calls -----------------------------------------------------
28!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
29!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
30!|| hm_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
31!||--- uses -----------------------------------------------------
32!|| elbuftag_mod ../starter/share/modules1/elbuftag_mod.F
33!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_mod.F
34!|| submodel_mod ../starter/share/modules1/submodel_mod.F
35!||====================================================================
36 SUBROUTINE hm_read_fail_puck(FAIL ,
37 . MAT_ID ,FAIL_ID ,IRUPT ,
38 . LSUBMODEL,UNITAB ,FAIL_TAG )
39C-----------------------------------------------
40c ROUTINE DESCRIPTION :
41c Read Puck failure model parameters
42C-----------------------------------------------
43C M o d u l e s
44C-----------------------------------------------
45 USE fail_param_mod
46 USE unitab_mod
47 USE submodel_mod
48 USE hm_option_read_mod
49 USE elbuftag_mod
50C-----------------------------------------------
51C I m p l i c i t T y p e s
52C-----------------------------------------------
53#include "implicit_f.inc"
54C-----------------------------------------------
55C C o m m o n B l o c k s
56C-----------------------------------------------
57#include "units_c.inc"
58C-----------------------------------------------
59C D u m m y A r g u m e n t s
60C-----------------------------------------------
61 INTEGER ,INTENT(IN) :: FAIL_ID !< failure model ID
62 INTEGER ,INTENT(IN) :: MAT_ID !< material law ID
63 INTEGER ,INTENT(IN) :: IRUPT !< failure model type number
64 TYPE(unit_type_) ,INTENT(IN) :: UNITAB !< table of input units
65 TYPE(submodel_data),INTENT(IN) :: LSUBMODEL(*) !< submodel table
66 TYPE(fail_param_) ,INTENT(INOUT) :: FAIL !< failure model data structure
67 TYPE (FAIL_TAG_) ,INTENT(INOUT) :: FAIL_TAG !< failure model tag for buffer allocation
68C-----------------------------------------------
69C L o c a l V a r i a b l e s
70C-----------------------------------------------
71 INTEGER IFAIL_SH ,ISOLID
72 my_real
73 . sigt1,sigt2,sigt12,sigc1,sigc2,
74 . pn12,pp12,pn22,tmax,fcut,pthkf
75 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
76C--------------------------------------------------
77C S o u r c e L i n e s
78C--------------------------------------------------
79 is_encrypted = .false.
80 is_available = .false.
81
82 CALL hm_option_is_encrypted(is_encrypted)
83
84c---------------------------
85 CALL hm_get_floatv ('Sigma_1t' ,sigt1 ,is_available,lsubmodel,unitab)
86 CALL hm_get_floatv ('Sigma_2t' ,sigt2 ,is_available,lsubmodel,unitab)
87 CALL hm_get_floatv ('Sigma_12' ,sigt12 ,is_available,lsubmodel,unitab)
88 CALL hm_get_floatv ('sigma_1c' ,SIGC1 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
89 CALL HM_GET_FLOATV ('sigma_2c' ,SIGC2 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
90
91 CALL HM_GET_FLOATV ('p12_positive' ,PP12 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
92 CALL HM_GET_FLOATV ('p12_negative' ,PN12 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
93 CALL HM_GET_FLOATV ('p22_negative' ,PN22 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
94 CALL HM_GET_FLOATV ('tau_max' ,TMAX ,IS_AVAILABLE,LSUBMODEL,UNITAB)
95 CALL HM_GET_INTV ('ifail_sh' ,IFAIL_SH ,IS_AVAILABLE,LSUBMODEL)
96 CALL HM_GET_INTV ('ifail_so' ,ISOLID ,IS_AVAILABLE,LSUBMODEL)
97
98 CALL HM_GET_FLOATV ('fcut' ,FCUT ,IS_AVAILABLE,LSUBMODEL,UNITAB)
99c---------------------------
100 IF (SIGT1 == ZERO) SIGT1 = EP20
101 IF (SIGT2 == ZERO) SIGT2 = EP20
102 IF (SIGT12 == ZERO) SIGT12 = EP20
103 IF (SIGC1 == ZERO) SIGC1 = EP20
104 IF (SIGC2 == ZERO) SIGC2 = EP20
105 IF (TMAX == ZERO) TMAX = EP20
106 IF (IFAIL_SH == 0) IFAIL_SH = 1
107 IF (ISOLID == 0) ISOLID = 1
108c
109 IF (IFAIL_SH == 1) THEN
110 PTHKF = EM06
111 ELSEIF (IFAIL_SH == 2) THEN
112 PTHKF = ONE
113 ENDIF
114 FCUT = MAX(ZERO,FCUT)
115c---------------------------
116 FAIL%KEYWORD = 'puck'
117 FAIL%IRUPT = IRUPT
118 FAIL%FAIL_ID = FAIL_ID
119 FAIL%NUPARAM = 12
120 FAIL%NIPARAM = 0
121 FAIL%NUVAR = 14
122 FAIL%NFUNC = 0
123 FAIL%NTABLE = 0
124 FAIL%NMOD = 5
125 FAIL%PTHK = PTHKF
126c
127 ALLOCATE (FAIL%UPARAM(FAIL%NUPARAM))
128 ALLOCATE (FAIL%IPARAM(FAIL%NIPARAM))
129 ALLOCATE (FAIL%IFUNC (FAIL%NFUNC))
130 ALLOCATE (FAIL%TABLE (FAIL%NTABLE))
131c
132 ! Modes of failure
133 FAIL_TAG%LF_DAMMX = FAIL_TAG%LF_DAMMX + FAIL%NMOD
134 ALLOCATE (FAIL%MODE(FAIL%NMOD))
135 FAIL%MODE(1) = "Tensile fiber failure"
136 FAIL%MODE(2) = "Compression fiber failure"
137 FAIL%MODE(3) = "Inter-fiber failure A"
138 FAIL%MODE(4) = "Inter-fiber failure B"
139 FAIL%MODE(5) = "Inter-fiber failure C"
140c
141 FAIL%UPARAM(1) = SIGT1
142 FAIL%UPARAM(2) = SIGT2
143 FAIL%UPARAM(3) = SIGC1
144 FAIL%UPARAM(4) = SIGC2
145 FAIL%UPARAM(5) = SIGT12
146 FAIL%UPARAM(6) = PP12
147 FAIL%UPARAM(7) = PN12
148 FAIL%UPARAM(8) = PN22
149 FAIL%UPARAM(9) = TMAX
150 FAIL%UPARAM(10) = IFAIL_SH
151 FAIL%UPARAM(11) = ISOLID
152 FAIL%UPARAM(12) = FCUT
153c---------------------------
154 IF(IS_ENCRYPTED)THEN
155
156 ELSE
157 WRITE(IOUT, 1000)SIGT1,SIGT2,SIGT12,SIGC1,SIGC2,PP12,PN12,PN22,TMAX,FCUT
158 IF(IFAIL_SH==1)THEN
159 WRITE(IOUT, 1100)
160 ELSEIF(IFAIL_SH==2)THEN
161 WRITE(IOUT, 1200)
162 ENDIF
163 ENDIF
164c-----------
165 1000 FORMAT(
166 & 5X,40H MODIFIED PUCK FAILURE MODEL /,
167 & 5X,40H ------------------------ /,
168 & 5X,'longitudinal tensile strength . . . . . . . =',E12.4/
169 & 5X,'transversal tensile strength . . . . . . . =',E12.4/
170 & 5X,'shear strength . . . . . . . . . . . . . . =',E12.4/
171 & 5X,'longitudinal compressive strength . . . . . =',E12.4/
172 & 5X,'transversal compressive strength . . . . =',E12.4/
173 & 5X,'failure envelop parameter in direction 12 (-) . =',E12.4/
174 & 5X,'failure envelop parameter in direction 12 (+) . =',E12.4/
175 & 5X,'failure envelop parameter in direction 22 (-) . =',E12.4/
176 & 5X,' relaxation time . . . . . . . . . . . . . =',E12.4/
177 & 5X,'stress tensor filtering frequency . . . . . =',E12.4/)
178 1100 FORMAT(
179 & 5X,' shell is deleted IF criteria is reached for one layer ')
180 1200 FORMAT(
181 & 5X,' shell tensor stress is set to zero and shell is deleted',/
182 & 5X,' IF criteria is reached for all layers' )
183c-----------
184 RETURN
185 END
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:418
hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:288
hm_read_fail_puck
subroutine hm_read_fail_puck(fail, mat_id, fail_id, irupt, lsubmodel, unitab, fail_tag)
Definition hm_read_fail_puck.F:39
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
elbuftag_mod
Definition elbuftag_mod.F:230
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