OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_mat93.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_mat93 ../starter/source/materials/mat/mat093/hm_read_mat93.F
25!||--- called by ------------------------------------------------------
26!|| hm_read_mat ../starter/source/materials/mat/hm_read_mat.F90
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| hm_get_float_array_index ../starter/source/devtools/hm_reader/hm_get_float_array_index.F
30!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
31!|| hm_get_floatv_dim ../starter/source/devtools/hm_reader/hm_get_floatv_dim.F
32!|| hm_get_int_array_index ../starter/source/devtools/hm_reader/hm_get_int_array_index.F
33!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
34!|| hm_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
35!|| init_mat_keyword ../starter/source/materials/mat/init_mat_keyword.F
36!||--- uses -----------------------------------------------------
37!|| elbuftag_mod ../starter/share/modules1/elbuftag_mod.F
38!|| message_mod ../starter/share/message_module/message_mod.F
39!|| submodel_mod ../starter/share/modules1/submodel_mod.F
40!||====================================================================
41 SUBROUTINE hm_read_mat93(UPARAM ,MAXUPARAM,NUPARAM ,ISRATE , IMATVIS ,
42 . NUVAR ,IFUNC ,MAXFUNC ,NFUNC , PARMAT ,
43 . UNITAB ,MAT_ID ,TITR ,MTAG , LSUBMODEL,
44 . PM ,IPM ,MATPARAM ,NVARTMP )
45C-----------------------------------------------
46C D e s c r i p t i o n
47C-----------------------------------------------
48C READ MAT LAW93 WITH HM READER ( TO BE COMPLETED )
49C
50C DUMMY ARGUMENTS DESCRIPTION:
51C ===================
52C
53C NAME DESCRIPTION
54C
55C PM MATERIAL ARRAY(REAL)
56C UNITAB UNITS ARRAY
57C MAT_ID MATERIAL ID(INTEGER)
58C TITR MATERIAL TITLE
59C LSUBMODEL SUBMODEL STRUCTURE
60C
61C-----------------------------------------------
62C M o d u l e s
63C-----------------------------------------------
64 USE unitab_mod
65 USE elbuftag_mod
66 USE message_mod
67 USE submodel_mod
68 USE matparam_def_mod
69 USE names_and_titles_mod , ONLY : nchartitle
70C-----------------------------------------------
71C I m p l i c i t T y p e s
72C-----------------------------------------------
73#include "implicit_f.inc"
74C-----------------------------------------------
75C C o m m o n B l o c k s
76C-----------------------------------------------
77#include "units_c.inc"
78#include "param_c.inc"
79C-----------------------------------------------
80C D u m m y A r g u m e n t s
81C-----------------------------------------------
82 TYPE (UNIT_TYPE_),INTENT(IN) :: UNITAB
83 my_real, INTENT(INOUT) :: PM(NPROPM),PARMAT(100),UPARAM(MAXUPARAM)
84 INTEGER, INTENT(INOUT) :: IPM(NPROPMI),ISRATE,IFUNC(MAXFUNC),NFUNC,
85 . maxfunc,maxuparam,nuparam,nuvar,imatvis ,
86 . nvartmp
87 TYPE(mlaw_tag_),INTENT(INOUT) :: MTAG
88 INTEGER,INTENT(IN) :: MAT_ID
89 CHARACTER(LEN=NCHARTITLE) ,INTENT(IN) :: TITR
90 TYPE(submodel_data),INTENT(IN) :: LSUBMODEL(*)
91 TYPE(matparam_struct_) ,INTENT(INOUT) :: MATPARAM
92C-----------------------------------------------
93C L o c a l V a r i a b l e s
94C-----------------------------------------------
95 INTEGER J,NRATE,I,ILAW,VP
96 my_real
97 . E11,E22,E33,NU12,NU23,NU13,G12,G13,G23,QR1,QR2,CR1,CR2,
98 . sigy,r11,r22,r33,r12,r13,r23,a1,a2,a3,hh,ff,gg,ll,mm,nn,
99 . d11,d22,d33,d12,d13,d23,a11,a22,a12,c11,c22,c33,c12,c13,
100 . c23,nu21,nu31,nu32,detc,fac,yfac(100),rate(100),dmin,dmax,
101 . yscale_unit,rho0,rhor,asrate
102 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
103C-----------------------------------------------
104C S o u r c e L i n e s
105C-----------------------------------------------
106 is_encrypted = .false.
107 is_available = .false.
108 ilaw = 93
109c
110c-----------------------------------------------
111 CALL hm_option_is_encrypted(is_encrypted)
112c-----------------------------------------------
113c
114card1 - Density
115 CALL hm_get_floatv('MAT_RHO',rho0 ,is_available, lsubmodel, unitab)
116card2 - Orthotropic elastic parameters
117 CALL hm_get_floatv('LAW93_E11' ,e11 ,is_available, lsubmodel, unitab)
118 CALL hm_get_floatv('LAW93_E22' ,e22 ,is_available, lsubmodel, unitab)
119 CALL hm_get_floatv('LAW93_E33' ,e33 ,is_available, lsubmodel, unitab)
120 CALL hm_get_floatv('LAW93_G12' ,g12 ,is_available, lsubmodel, unitab)
121 CALL hm_get_floatv('LAW93_Nu12',nu12 ,is_available, lsubmodel, unitab)
122card3 - Orthotropic elastic parameters
123 CALL hm_get_floatv('LAW93_G13' ,g13 ,is_available, lsubmodel, unitab)
124 CALL hm_get_floatv('LAW93_G23' ,g23 ,is_available, lsubmodel, unitab)
125 CALL hm_get_floatv('LAW93_Nu13',nu13 ,is_available, lsubmodel, unitab)
126 CALL hm_get_floatv('LAW93_Nu23',nu23 ,is_available, lsubmodel, unitab)
127card4 - Number of curves for each rate
128 CALL hm_get_intv ('LAW93_NL' ,nrate ,is_available, lsubmodel)
129 CALL hm_get_floatv('FCUT' ,asrate ,is_available, lsubmodel, unitab)
130 CALL hm_get_intv ('VP' ,vp ,is_available, lsubmodel)
131card5 - Curves parameters
132 IF (nrate > 0) THEN
133 DO i=1,nrate
134 CALL hm_get_int_array_index ('LAW93_arr1',ifunc(i) ,i ,is_available, lsubmodel)
135 CALL hm_get_float_array_index('LAW93_arr2',yfac(i) ,i ,is_available, lsubmodel, unitab)
136 CALL hm_get_float_array_index('LAW93_arr3',rate(i) ,i ,is_available, lsubmodel, unitab)
137 IF (yfac(i) == zero) THEN
138 CALL hm_get_floatv_dim('LAW93_arr2' ,yscale_unit ,is_available, lsubmodel, unitab)
139 yfac(i) = one * yscale_unit
140 ENDIF
141 ENDDO
142 ENDIF
143card6 - Continuous hardening yield stress of Voce combination
144 CALL hm_get_floatv('LAW93_Sigma_y',sigy ,is_available, lsubmodel, unitab)
145 CALL hm_get_floatv('LAW93_QR1' ,qr1 ,is_available, lsubmodel, unitab)
146 CALL hm_get_floatv('LAW93_CR1' ,cr1 ,is_available, lsubmodel, unitab)
147 CALL hm_get_floatv('LAW93_QR2' ,qr2 ,is_available, lsubmodel, unitab)
148 CALL hm_get_floatv('LAW93_CR2' ,cr2 ,is_available, lsubmodel, unitab)
149card7 - Hill parameters (normalized yield stresses)
150 CALL hm_get_floatv('LAW93_R11' ,r11 ,is_available, lsubmodel, unitab)
151 CALL hm_get_floatv('LAW93_R22' ,r22 ,is_available, lsubmodel, unitab)
152 CALL hm_get_floatv('LAW93_R12' ,r12 ,is_available, lsubmodel, unitab)
153card8 - Hill parameters (normalized yield stresses)
154 CALL hm_get_floatv('LAW93_R33' ,r33 ,is_available, lsubmodel, unitab)
155 CALL hm_get_floatv('LAW93_R13' ,r13 ,is_available, lsubmodel, unitab)
156 CALL hm_get_floatv('LAW93_R23' ,r23 ,is_available, lsubmodel, unitab)
157c
158 !========== DEFAULT VALUES=============!
159C
160 ! Default value for functions
161 nfunc = nrate
162 IF (nrate > 1) THEN
163 IF (rate(1) == zero) THEN
164 nfunc = nrate
165 ELSE
166 nfunc = nrate + 1
167 DO j = nrate,1,-1
168 ifunc(j+1) = ifunc(j)
169 rate(j+1) = rate(j)
170 yfac(j+1) = yfac(j)
171 ENDDO
172 rate(1) = zero
173 ENDIF
174 ENDIF
175C
176 ! Yield stresses
177 IF(sigy == zero) sigy = infinity
178 IF(r11 == zero) r11 = one
179 IF(r22 == zero) r22 = one
180 IF(r33 == zero) r33 = one
181 IF(r12 == zero) r12 = one
182 IF(r23 == zero) r23 = one
183 IF(r13 == zero) r13 = one
184C
185 ! Young modulus
186 IF (e22 == zero) e22 = e11
187 IF (e33 == zero) e33 = e22
188 ! Shear modulus
189 IF (g13 == zero) g13 = g12
190 IF (g23 == zero) g23 = g12
191 ! Remaining Poisson's ratio
192 nu21 = nu12*e22/e11
193 nu31 = nu13*e33/e11
194 nu32 = nu23*e33/e22
195 !check stability
196 ! checking poisson's ratio
197 if(nu12*nu21 >= one ) then
198 call ancmsg(msgid=3068,
199 . msgtype=msgerror,
200 . anmode=aninfo_blind_2,
201 . i1=mat_id ,
202 . c1=titr)
203 else if(nu13*nu31 >= one ) then
204 call ancmsg(msgid=3069,
205 . msgtype=msgerror,
206 . anmode=aninfo_blind_2,
207 . i1=mat_id ,
208 . c1=titr)
209 else if(nu23*nu32 >= one ) then
210 call ancmsg(msgid=3070,
211 . msgtype=msgerror,
212 . anmode=aninfo_blind_2,
213 . i1=mat_id ,
214 . c1=titr)
215 endif
216C
217 ! Hill coefficient
218 a1 = one/r11/r11
219 a2 = one/r22/r22
220 a3 = one/r33/r33
221 ff = half*(a2 + a3 - a1)
222 gg = half*(a3 + a1 - a2)
223 hh = half*(a1 + a2 - a3)
224 ll = three_half/r23/r23
225 mm = three_half/r13/r13
226 nn = three_half/r12/r12
227C
228 ! Elasticity matrix for 2D plane stress
229 fac = one/(one - nu12*nu21)
230 a11 = e11*fac
231 a12 = nu21*a11
232 a22 = e22*fac
233 ! Compliance matrix for 3D
234 c11 = one/e11
235 c22 = one/e22
236 c33 = one/e33
237 c12 =-nu12/e11
238 c13 =-nu31/e33
239 c23 =-nu23/e22
240 ! Checking input
241 detc= c11*c22*c33-c11*c23*c23-c12*c12*c33+c12*c13*c23
242 + +c13*c12*c23-c13*c22*c13
243 IF(detc<=zero) THEN
244 CALL ancmsg(msgid=307,
245 . msgtype=msgerror,
246 . anmode=aninfo,
247 . i1=mat_id,
248 . c1=titr)
249 ENDIF
250 ! 3D elastic matrix
251 d11 = (c22*c33-c23*c23)/detc
252 d12 =-(c12*c33-c13*c23)/detc
253 d13 = (c12*c23-c13*c22)/detc
254 d22 = (c11*c33-c13*c13)/detc
255 d23 =-(c11*c23-c13*c12)/detc
256 d33 = (c11*c22-c12*c12)/detc
257 dmin = min(d11*d22 -d12**2, d11*d33 - d13**2, d22*d33 - d23**2 )
258 dmax = max(d11,d22,d33)
259C
260 ! Strain-rate filtering
261 IF (nfunc > 1) THEN
262 israte = 1
263 IF (vp == 0) vp = 2
264 IF (vp == 1) THEN
265 asrate = 1.0d4*unitab%FAC_T_WORK
266 ELSE
267 IF (asrate == zero) asrate = 1.0d4*unitab%FAC_T_WORK
268 ENDIF
269 ELSE
270 israte = 0
271 asrate = zero
272 ENDIF
273C
274 ! PM table
275 rhor = zero
276 pm(1) = rhor
277 pm(89) = rho0
278C
279 ! PARMAT table
280 parmat(1) = max(a11,a22,d11,d22,d33)
281 parmat(2) = max(e11,e22,e33)
282 parmat(3) = max(nu12,nu13,nu23)
283 parmat(4) = israte
284 parmat(5) = asrate
285 parmat(16) = 1
286 parmat(17) = dmin/dmax/dmax
287C
288 ! MTAG variable activation
289 mtag%G_PLA = 1
290 mtag%L_PLA = 1
291 mtag%G_SEQ = 1
292 mtag%L_SEQ = 1
293 mtag%G_EPSD = 1
294 mtag%L_EPSD = 1
295C
296 ! MATPARAM parameters
297 CALL init_mat_keyword(matparam ,"ELASTO_PLASTIC")
298 CALL init_mat_keyword(matparam ,"INCREMENTAL" )
299 CALL init_mat_keyword(matparam ,"LARGE_STRAIN" )
300 CALL init_mat_keyword(matparam ,"HOOK")
301 CALL init_mat_keyword(matparam,"ORTHOTROPIC")
302C
303 ! Properties compatibility
304 CALL init_mat_keyword(matparam,"SHELL_ORTHOTROPIC")
305 CALL init_mat_keyword(matparam,"SOLID_ORTHOTROPIC")
306 CALL init_mat_keyword(matparam,"SPH")
307C
308 ! No viscosity
309 imatvis = 0
310C
311 ! Number of user variable
312 IF ((nrate > 1).AND.(vp /= 2)) THEN
313 nuvar = 1
314 ELSE
315 nuvar = 0
316 ENDIF
317
318 ! Number of material parameter
319 nuparam = 30 + 2*nfunc
320 nvartmp = nfunc
321c
322 ! Filling the parameter table
323 ! -> Elastic parameters
324 uparam(1) = a11
325 uparam(2) = a22
326 uparam(3) = a12
327 uparam(4) = d11
328 uparam(5) = d12
329 uparam(6) = d13
330 uparam(7) = d22
331 uparam(8) = d23
332 uparam(9) = d33
333 uparam(10) = g12
334 uparam(11) = g13
335 uparam(12) = g23
336 uparam(13) = e11
337 uparam(14) = e22
338 uparam(15) = e33
339 uparam(16) = nu12
340 uparam(17) = nu13
341 uparam(18) = nu23
342 ! -> Yield criterion parameters
343 uparam(19) = ff
344 uparam(20) = gg
345 uparam(21) = hh
346 uparam(22) = ll
347 uparam(23) = mm
348 uparam(24) = nn
349 ! -> Continuous hardening parameters
350 uparam(25) = sigy
351 uparam(26) = qr1
352 uparam(27) = cr1
353 uparam(28) = qr2
354 uparam(29) = cr2
355 ! -> Strain-rate computation flag
356 uparam(30) = vp
357 ! -> Tabulated hardening parameters
358 IF (nfunc > 0) THEN
359 DO j=1,nfunc
360 uparam(30 + j) = rate(j)
361 uparam(30 + nfunc + j) = yfac(j)
362 ENDDO
363 ENDIF
364c
365c--------------------------
366c Parameters printout
367c--------------------------
368 WRITE(iout,1001) trim(titr),mat_id,ilaw
369 WRITE(iout,1000)
370 IF(is_encrypted)THEN
371 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA'
372 ELSE
373 WRITE(iout,1002) rho0
374 WRITE(iout,1300) e11,e22,e33,g12,g13,g23,nu12,nu13,nu23
375 IF (nrate == 0) THEN
376 WRITE(iout,1450)
377 WRITE(iout,1400) sigy,qr1,cr1,qr2,cr2
378 ELSE
379 WRITE(iout,1550)
380 DO j=1,nfunc
381 WRITE(iout,1500) ifunc(j),yfac(j),rate(j)
382 ENDDO
383 IF (nrate > 1) THEN
384 WRITE(iout,1575) asrate,vp
385 ENDIF
386 ENDIF
387 WRITE(iout,1600) r11,r22,r33,r12,r13,r23
388 ENDIF
389C-----------------------------------------------------------------
390 1000 FORMAT(
391 & 5x,' ORTHOTROPIC ELASTIC + HILL CRITERION '/,
392 & 5x,' ------------------------------------ '//)
393 1001 FORMAT(
394 & 5x,a,/,
395 & 5x,'MATERIAL NUMBER . . . . . . . . . . .=',i10/,
396 & 5x,'MATERIAL LAW. . . . . . . . . . . . .=',i10/)
397 1002 FORMAT(
398 & 5x,'INITIAL DENSITY . . . . . . . . . . .=',1pg20.13/)
399 1300 FORMAT(
400 & 5x,'YOUNG MODULUS IN 11 DIRECTION . . . .=',1pg20.13/,
401 & 5x,'YOUNG MODULUS IN 22 DIRECTION . . . .=',1pg20.13/,
402 & 5x,'YOUNG MODULUS IN 33 DIRECTION . . . .=',1pg20.13/,
403 & 5x,'SHEAR MODULUS IN 12 DIRECTION . . . .=',1pg20.13/,
404 & 5x,'SHEAR MODULUS IN 13 DIRECTION . . . .=',1pg20.13/,
405 & 5x,'SHEAR MODULUS IN 23 DIRECTION . . . .=',1pg20.13/,
406 & 5x,'POISSON RATIO 12. . . . . . . . . . .=',1pg20.13/,
407 & 5x,'POISSON RATIO 13. . . . . . . . . . .=',1pg20.13/,
408 & 5x,'POISSON RATIO 23. . . . . . . . . . .=',1pg20.13//)
409 1550 FORMAT(
410 & 5x,'--------------------------------------'/,
411 & 5x,'TABULATED YIELD STRESS '/,
412 & 5x,'--------------------------------------'//)
413 1500 FORMAT(
414 & 5x,'YIELD STRESS FUNCTION NUMBER. . . . .=',i10/,
415 & 5x,'YIELD SCALE FACTOR. . . . . . . . . .=',1pg20.13/,
416 & 5x,'STRAIN RATE . . . . . . . . . . . . .=',1pg20.13/)
417 1575 FORMAT(
418 & 5x,'STRAIN RATE CUTTING FREQUENCY . . . .=',1pg20.13/
419 & 5x,'STRAIN RATE CHOICE FLAG . . . . . . .=',i10/
420 & 5x,' VP=1 EQUIVALENT PLASTIC STRAIN RATE'/
421 & 5x,' VP=2 TOTAL STRAIN RATE (DEFAULT)'/
422 & 5x,' VP=3 DEVIATORIC STRAIN RATE'/)
423 1450 FORMAT(
424 & 5x,'--------------------------------------'/,
425 & 5x,'CONTINUOUS YIELD STRESS '/,
426 & 5x,'--------------------------------------'//)
427 1400 FORMAT(
428 & 5x,'INITIAL YIELD STRESS. . . . . . . . .=',1pg20.13/,
429 & 5x,'PARAMETER QR1 OF HARDENING . . . . .=',1pg20.13/,
430 & 5x,'PARAMETER CR1 OF HARDENING . . . . .=',1pg20.13/,
431 & 5x,'PARAMETER QR2 OF HARDENING . . . . .=',1pg20.13/,
432 & 5x,'PARAMETER CR2 OF HARDENING . . . . .=',1pg20.13/,
433 & 5x,'REFERENCE STRAIN. . . . . . . . . . .=',1pg20.13//)
434 1600 FORMAT(
435 & 5x,'RATIO YIELD PARAMETER R11 . . . . . .=',1pg20.13/,
436 & 5x,'RATIO YIELD PARAMETER R22 . . . . . .=',1pg20.13/,
437 & 5x,'RATIO YIELD PARAMETER R33 . . . . . .=',1pg20.13/,
438 & 5x,'RATIO YIELD PARAMETER R12 . . . . . .=',1pg20.13/,
439 & 5x,'RATIO YIELD PARAMETER R13 . . . . . .=',1pg20.13/,
440 & 5x,'RATIO YIELD PARAMETER R23 . . . . . .=',1pg20.13/)
441 RETURN
442 END
hm_get_float_array_index
subroutine hm_get_float_array_index(name, rval, index, is_available, lsubmodel, unitab)
Definition hm_get_float_array_index.F:105
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:415
hm_get_floatv_dim
subroutine hm_get_floatv_dim(name, dim_fac, is_available, lsubmodel, unitab)
Definition hm_get_floatv_dim.F:159
hm_get_int_array_index
subroutine hm_get_int_array_index(name, ival, index, is_available, lsubmodel)
Definition hm_get_int_array_index.F:164
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_mat93
subroutine hm_read_mat93(uparam, maxuparam, nuparam, israte, imatvis, nuvar, ifunc, maxfunc, nfunc, parmat, unitab, mat_id, titr, mtag, lsubmodel, pm, ipm, matparam, nvartmp)
Definition hm_read_mat93.F:45
init_mat_keyword
subroutine init_mat_keyword(matparam, keyword)
Definition init_mat_keyword.F:150
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
message_mod
Definition message_mod.F:1249
names_and_titles_mod
Definition names_and_titles_mod.F:997
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1003
submodel_mod
Definition submodel_mod.F:712
unitab_mod
Definition unitab_mod.F:513
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889
elbuftag_mod::mlaw_tag_
Definition elbuftag_mod.F:232
submodel_mod::submodel_data
Definition submodel_mod.F:723