OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_fail_hashin.F File Reference
#include "implicit_f.inc"
#include "units_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine hm_read_fail_hashin (fail, mat_id, fail_id, irupt, frwave, lsubmodel, unitab, fail_tag)

Function/Subroutine Documentation

◆ hm_read_fail_hashin()

subroutine hm_read_fail_hashin ( type(fail_param_), intent(inout) fail,
integer, intent(in) mat_id,
integer, intent(in) fail_id,
integer, intent(in) irupt,
integer, intent(inout) frwave,
type(submodel_data), dimension(*), intent(in) lsubmodel,
type(unit_type_), intent(in) unitab,
type (fail_tag_), intent(inout) fail_tag )
Parameters
[in]fail_idfailure model ID
[in]mat_idmaterial law ID
[in]iruptfailure model type number
[in]unitabtable of input units
[in]lsubmodelsubmodel table
[in,out]frwavefailure model type number
[in,out]failfailure model data structure
[in,out]fail_tagfailure model tag for buffer allocation

Definition at line 36 of file hm_read_fail_hashin.F.

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-----------
my_real
#define my_real
Definition cppsort.cpp:32
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
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
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