OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_fail_hashin.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_hashin ../starter/source/materials/fail/hashin/hm_read_fail_hashin.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_hashin(FAIL ,
37 . MAT_ID ,FAIL_ID ,IRUPT ,FRWAVE ,
38 . LSUBMODEL,UNITAB ,FAIL_TAG)
39C-----------------------------------------------
40c ROUTINE DESCRIPTION :
41c Read Hashin 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 INTEGER ,INTENT(INOUT) :: FRWAVE !< failure model type number
67 TYPE(fail_param_) ,INTENT(INOUT) :: FAIL !< failure model data structure
68 TYPE (FAIL_TAG_) ,INTENT(INOUT) :: FAIL_TAG !< failure model tag for buffer allocation
69C-----------------------------------------------
70C L o c a l V a r i a b l e s
71C-----------------------------------------------
72 INTEGER :: IFAIL_SH,IMODEL,ISOLID,DMG_FLAG,TMOD
73 my_real :: SIGT1,SIGT2,SIGT3,SIGT12,SIGC1,SIGC2,SOFT1,SOFT2,SOFT3,
74 . fsig12,msig12,msig23,msig13,angle,sdelam,csig,pthkf,
75 . tmax,ratio,tcut,epsp_ref,epsp_min
76 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
77C=======================================================================
78 is_encrypted = .false.
79 is_available = .false.
80c
81 epsp_min = zero
82 soft1 = one
83 soft2 = one
84 soft3 = one
85C--------------------------------------------------
86C check crypting
87C--------------------------------------------------
88c
89 CALL hm_option_is_encrypted(is_encrypted)
90c
91C--------------------------------------------------
92C Read DATA
93C--------------------------------------------------
94 CALL hm_get_intv ('Iform' ,imodel ,is_available,lsubmodel)
95 CALL hm_get_intv ('Ifail_sh' ,ifail_sh ,is_available,lsubmodel)
96 CALL hm_get_intv ('Ifail_so' ,isolid ,is_available,lsubmodel)
97 CALL hm_get_floatv ('ratio' ,ratio ,is_available,lsubmodel,unitab)
98 CALL hm_get_intv ('I_Dam' ,dmg_flag ,is_available,lsubmodel)
99 CALL hm_get_intv ('Imod' ,tmod ,is_available,lsubmodel)
100 CALL hm_get_intv ('I_frwave' ,frwave ,is_available,lsubmodel)
101 CALL hm_get_floatv ('Epsp_min' ,epsp_min ,is_available,lsubmodel,unitab)
102c
103 CALL hm_get_floatv ('Sigma_1t' ,sigt1 ,is_available,lsubmodel,unitab)
104 CALL hm_get_floatv ('Sigma_2t' ,sigt2 ,is_available,lsubmodel,unitab)
105 CALL hm_get_floatv ('Sigma_3t' ,sigt3 ,is_available,lsubmodel,unitab)
106 CALL hm_get_floatv ('sigma_1c' ,SIGC1 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
107 CALL HM_GET_FLOATV ('sigma_2c' ,SIGC2 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
108c
109 CALL HM_GET_FLOATV ('sigma_c' ,CSIG ,IS_AVAILABLE,LSUBMODEL,UNITAB)
110 CALL HM_GET_FLOATV ('sigma_12f' ,FSIG12 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
111 CALL HM_GET_FLOATV ('sigma_12m' ,MSIG12 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
112 CALL HM_GET_FLOATV ('sigma_23m' ,MSIG23 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
113 CALL HM_GET_FLOATV ('sigma_13m' ,MSIG13 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
114c
115 CALL HM_GET_FLOATV ('phi' ,ANGLE ,IS_AVAILABLE,LSUBMODEL,UNITAB)
116 CALL HM_GET_FLOATV ('sdel' ,SDELAM ,IS_AVAILABLE,LSUBMODEL,UNITAB)
117 CALL HM_GET_FLOATV ('tau_max' ,TMAX ,IS_AVAILABLE,LSUBMODEL,UNITAB)
118 CALL HM_GET_FLOATV ('epsp_ref' ,EPSP_REF ,IS_AVAILABLE,LSUBMODEL,UNITAB)
119 CALL HM_GET_FLOATV ('tcut' ,TCUT ,IS_AVAILABLE,LSUBMODEL,UNITAB)
120c
121 CALL HM_GET_FLOATV ('soft' ,SOFT1 ,IS_AVAILABLE,LSUBMODEL,UNITAB)
122c--------------------------------------------------
123c Default values
124c--------------------------------------------------
125 IF (FRWAVE == 0) FRWAVE = 1
126 IF (FRWAVE > 1) FRWAVE = 1
127 IF (SIGT1 == ZERO) SIGT1 = EP20
128 IF (SIGT2 == ZERO) SIGT2 = EP20
129 IF (SIGT3 == ZERO) SIGT3 = EP20
130 IF (SIGC1 == ZERO) SIGC1 = EP20
131 IF (SIGC2 == ZERO) SIGC2 = EP20
132 IF (CSIG == ZERO) CSIG = EP20
133 IF (FSIG12 == ZERO) FSIG12 = EP20
134 IF (MSIG12 == ZERO) MSIG12 = EP20
135 IF (MSIG13 == ZERO) MSIG13 = EP20
136 IF (MSIG23 == ZERO) MSIG23 = EP20
137 IF (TMAX == ZERO) TMAX = EP20
138 IF (TCUT == ZERO) TCUT = TMAX
139 IF (EPSP_REF == ZERO) EPSP_REF = EM20
140 IF (SDELAM == ZERO) SDELAM = ONE
141 IF (IMODEL == 0) IMODEL = 1
142 IF (IFAIL_SH == 0) IFAIL_SH = 1
143 IF (ISOLID == 0) ISOLID = 1
144 IF (RATIO == ZERO) RATIO = ONE - EM06
145 RATIO = MIN(RATIO, ONE)
146 RATIO = MAX(RATIO,-ONE)
147 IF (DMG_FLAG == 0) DMG_FLAG = 1 ! DMG_FLAG = 1 => damage reduces forces after criteria is met
148 IF (DMG_FLAG == 2) DMG_FLAG = 0 ! DMG_FLAG = 2->0 => damage reduces saved stress (old method)
149 SOFT1 = MIN(ONE,MAX(ZERO,SOFT1))
150c
151 IF (IFAIL_SH == 1) THEN
152 PTHKF = EM06
153 ELSEIF (IFAIL_SH == 2) THEN
154 PTHKF = RATIO
155 ELSE ! IFAIL_SH = 3
156 PTHKF = ONE - EM06
157 ENDIF
158c--------------------------------------------------
159 FAIL%KEYWORD = 'hashin'
160 FAIL%IRUPT = IRUPT
161 FAIL%FAIL_ID = FAIL_ID
162 FAIL%NUPARAM = 26
163 FAIL%NIPARAM = 0
164 FAIL%NUVAR = 12
165 FAIL%NFUNC = 0
166 FAIL%NTABLE = 0
167 FAIL%PTHK = PTHKF
168c
169 ALLOCATE (FAIL%UPARAM(FAIL%NUPARAM))
170 ALLOCATE (FAIL%IPARAM(FAIL%NIPARAM))
171 ALLOCATE (FAIL%IFUNC (FAIL%NFUNC))
172 ALLOCATE (FAIL%TABLE (FAIL%NTABLE))
173c
174 ! Modes of failure
175 ! -> Uni-directional lamina model
176 IF (IMODEL == 1) THEN
177 FAIL%NMOD = 5
178 ALLOCATE (FAIL%MODE(FAIL%NMOD))
179 FAIL%MODE(1) = "Tensile/Shear fiber"
180 FAIL%MODE(2) = "Compression fiber"
181 FAIL%MODE(3) = "Crush"
182 FAIL%MODE(4) = "Matrix failure"
183 FAIL%MODE(5) = "Delamination"
184 ! -> Fabric lamina model
185 ELSE
186 FAIL%NMOD = 7
187 ALLOCATE (FAIL%MODE(FAIL%NMOD))
188 FAIL%MODE(1) = "Tensile/Shear fiber 1"
189 FAIL%MODE(2) = "Tensile/Shear fiber 2"
190 FAIL%MODE(3) = "Compression fiber 1"
191 FAIL%MODE(4) = "Compression fiber 2"
192 FAIL%MODE(5) = "Crush"
193 FAIL%MODE(6) = "Shear failure matrix"
194 FAIL%MODE(7) = "Matrix failure"
195 ENDIF
196 FAIL_TAG%LF_DAMMX = FAIL_TAG%LF_DAMMX + FAIL%NMOD
197c
198 FAIL%UPARAM(1) = IMODEL
199 FAIL%UPARAM(2) = SIGT1
200 FAIL%UPARAM(3) = SIGT2
201 FAIL%UPARAM(4) = SIGT3
202 FAIL%UPARAM(5) = SIGC1
203 FAIL%UPARAM(6) = SIGC2
204 FAIL%UPARAM(7) = CSIG
205 FAIL%UPARAM(8) = FSIG12
206 FAIL%UPARAM(9) = MSIG12
207 FAIL%UPARAM(10) = MSIG13
208 FAIL%UPARAM(11) = MSIG23
209 FAIL%UPARAM(12) = ANGLE*PI/HUNDRED80 ! degre ---> radian
210 FAIL%UPARAM(13) = SDELAM
211 FAIL%UPARAM(14) = TMAX
212 FAIL%UPARAM(15) = IFAIL_SH
213 FAIL%UPARAM(16) = ISOLID
214 FAIL%UPARAM(17) = RATIO
215 FAIL%UPARAM(18) = DMG_FLAG
216 FAIL%UPARAM(19) = TMOD
217 FAIL%UPARAM(20) = EPSP_REF
218 FAIL%UPARAM(21) = TCUT
219 FAIL%UPARAM(22) = EPSP_MIN
220 FAIL%UPARAM(23) = 0 !not used
221 FAIL%UPARAM(24) = SOFT1
222 FAIL%UPARAM(25) = SOFT2
223 FAIL%UPARAM(26) = SOFT3
224c-----------------------------------------------------
225 IF(IS_ENCRYPTED)THEN
226 WRITE(IOUT, 1000)
227 ELSE
228 WRITE(IOUT, 1100) MAT_ID,FAIL_ID,SIGT1,SIGT2,SIGT3,SIGC1,SIGC2,
229 . CSIG,FSIG12,MSIG12,MSIG23,MSIG13,
230 . ANGLE,SDELAM, IMODEL,RATIO,TMAX,
231 . TMOD,EPSP_REF,TCUT,EPSP_MIN,FRWAVE
232 IF (FRWAVE > 1) WRITE(IOUT, 1200) SOFT1
233 ENDIF
234c-----------
235 RETURN
236c--------------------------------------------------
237 1000 FORMAT(
238 & 5X,' crypted DATA in failure model '/,
239 & 5X,' ----------------------------- '/)
240 1100 FORMAT(
241 & 5X,' hashin failure model ',/
242 & 5X,' ------------------------ ',/
243 & 5X,' mat_id =',I10/
244 & 5X,' fail_id =',i10/
245 & 5x,'Longitudinal tensile strength. . . . . . . . =',1pg20.13/
246 & 5x,'Transversal tensile strength . . . . . . . . =',1pg20.13/
247 & 5x,'Through thickness strength. . . . . . . . . =',1pg20.13/
248 & 5x,'Longitudinal compressive strength. . . . . . =',1pg20.13/
249 & 5x,'Transversal compressive strength. . . . . . =',1pg20.13/
250 & 5x,'Crush strength . . . . . . . . . . . . . =',1pg20.13/
251 & 5x,'Fiber Mode shear strength . . . . . . . . . =',1pg20.13/
252 & 5x,'Matrix mode shear strength (plane 12 ) . . . =',1pg20.13/
253 & 5x,'Matrix mode shear strength (plane 23 ) . . . =',1pg20.13/
254 & 5x,'Matrix mode shear strength (plane 13 ) . . . =',1pg20.13/
255 & 5x,'Coulomb friction angle and delamination' /,
256 & 5x, ' failure < 90° =',1pg20.13/
257 & 5x,'Scale factor for delamination criteria . . . =',1pg20.13/
258 & 5x,'Material failure model: . . . . . . . . . . =',i10 /
259 & 15x, '1: Unidirectional layer model' ,/
260 & 15x, '2: Fabric layer model ' ,/
261 & 5x,'Element deletion parameter ratio. . . . . . .=',1pg20.13,/,
262 & 5x,' > 0.0 : Fraction of failed thickness ',/,
263 & 5x,' < 0.0 : Fraction of failed intg. points or layers',/,
264 & 5x,'Relaxation time . . . . . . : . . . . . . .=',1pg20.13 /
265 & 5x,'Relaxation time calculation : . . . . . . .=',i10 /
266 & 15x, '1 : Based on timestep' ,/
267 & 15x, 'Other value : Constant ' ,/
268 & 5x,'Ref strain rate. . . . . . . . . . . . . . . =',1pg20.13/
269 & 5x,'Strain rate cutoff period. . . . . . . . . . =',1pg20.13/
270 & 5x,'Low strain rate limit . . .. . . . . . . . . =',1pg20.13/
271 & 5x,'Front Wave : . . . . . . . . . . . . . . . =',i10 /
272 & 15x, '1: No frontwave propagation ',/
273 & 15x, '2: Isotropic propagation ',/)
274 1200 FORMAT(
275 & 5x,'Softening factor . . . . . . . . . . . . . .=',1pg20.13//)
276c-----------
277 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_hashin
subroutine hm_read_fail_hashin(fail, mat_id, fail_id, irupt, frwave, lsubmodel, unitab, fail_tag)
Definition hm_read_fail_hashin.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
submodel_mod::submodel_data
Definition submodel_mod.F:723
unitab_mod::unit_type_
Definition unitab_mod.F:520