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,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: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
subroutine hm_read_fail(mat_param, nummat, maxfail, fail_tag, ntable, table, failwave, nloc_dmg, unitab, lsubmodel)
Definition hm_read_fail.F:92
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
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
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
starter
program starter
Definition starter.F:39
submodel_mod::submodel_data
Definition submodel_mod.F:728