OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_mat10.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_mat10 ../starter/source/materials/mat/mat010/hm_read_mat10.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_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
30!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
31!|| hm_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
32!|| init_mat_keyword ../starter/source/materials/mat/init_mat_keyword.F
33!||--- uses -----------------------------------------------------
34!|| elbuftag_mod ../starter/share/modules1/elbuftag_mod.F
35!|| message_mod ../starter/share/message_module/message_mod.F
36!|| reader_old_mod ../starter/share/modules1/reader_old_mod.F90
37!|| submodel_mod ../starter/share/modules1/submodel_mod.F
38!||====================================================================
39 SUBROUTINE hm_read_mat10(
40 . LSUBMODEL ,MTAG ,UNITAB ,IPM ,PM ,
41 . MAT_ID ,TITR ,EOS_TAG ,ISRATE ,MATPARAM)
42C-----------------------------------------------
43C M o d u l e s
44C-----------------------------------------------
45 USE unitab_mod
46 USE elbuftag_mod !MAXEOS
47 USE message_mod
48 USE submodel_mod , ONLY : nsubmod, submodel_data
49 USE matparam_def_mod
50 USE names_and_titles_mod , ONLY : nchartitle
51 USE reader_old_mod , ONLY : irec
52C-----------------------------------------------
53C ROUTINE DESCRIPTION :
54C ===================
55C READ MAT LAW21 WITH HM READER
56C-----------------------------------------------
57C DUMMY ARGUMENTS DESCRIPTION:
58C ===================
59C UNITAB UNITS ARRAY
60C MAT_ID MATERIAL ID(INTEGER)
61C TITR MATERIAL TITLE
62C LSUBMODEL SUBMODEL STRUCTURE
63C-----------------------------------------------
64C I m p l i c i t T y p e s
65C-----------------------------------------------
66#include "implicit_f.inc"
67C-----------------------------------------------
68C C o m m o n B l o c k s
69C-----------------------------------------------
70#include "scr17_c.inc"
71#include "units_c.inc"
72#include "param_c.inc"
73C-----------------------------------------------
74C D u m m y A r g u m e n t s
75C-----------------------------------------------
76 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
77 INTEGER, INTENT(IN) :: MAT_ID
78 INTEGER, INTENT(INOUT) :: ISRATE
79 INTEGER, DIMENSION(NPROPMI) ,INTENT(INOUT) :: IPM
80 CHARACTER(LEN=NCHARTITLE) ,INTENT(IN) :: TITR
81 my_real, DIMENSION(NPROPM) ,INTENT(INOUT) :: pm
82 TYPE(submodel_data), DIMENSION(NSUBMOD),INTENT(IN) :: LSUBMODEL
83 TYPE(mlaw_tag_), INTENT(INOUT) :: MTAG
84 TYPE(eos_tag_) , TARGET, DIMENSION(0:MAXEOS) ,INTENT(INOUT) :: EOS_TAG
85 TYPE(matparam_struct_),INTENT(INOUT) :: MATPARAM
86C-----------------------------------------------
87C L o c a l V a r i a b l e s
88C-----------------------------------------------
89 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
90 INTEGER :: I,IFUNC,COUNT,ID,IFORM,IEOS
91 my_real :: rho0,rhor,e,nu,a0,a1,a2,amx,c0,c1,c2,c3,fac_y
92 my_real :: pmin,pext,xmumx,bunl,e0,g,pstar,psh,delta,det,young,anu
93 CHARACTER*64 :: chain,chain2
94C-----------------------------------------------
95C S o u r c e L i n e s
96C-----------------------------------------------
97 pstar = -infinity
98 c0 = zero
99 c1 = zero
100 c2 = zero
101 c3 = zero
102 pext = zero
103 psh = zero
104 pmin = zero
105 bunl = zero
106 xmumx = zero
107 israte = 0
108 is_encrypted = .false.
109 is_available = .false.
110!---
111 CALL hm_option_is_encrypted(is_encrypted)
112!
113 CALL hm_get_intv ('Line_count' ,count ,is_available,lsubmodel)
114 IF(count==0)count=3 ! format 4.* and older
115 count=count+2
116
117 CALL hm_get_floatv('MAT_RHO' ,rho0 ,is_available, lsubmodel, unitab)
118 CALL hm_get_floatv('Refer_Rho' ,rhor ,is_available, lsubmodel, unitab)
119 CALL hm_get_floatv('MAT_E' ,young ,is_available, lsubmodel, unitab)
120 CALL hm_get_floatv('MAT_NU' ,anu ,is_available, lsubmodel, unitab)
121 CALL hm_get_floatv('MAT_A0' ,a0 ,is_available, lsubmodel, unitab)
122 CALL hm_get_floatv('MAT_A1' ,a1 ,is_available, lsubmodel, unitab)
123 CALL hm_get_floatv('MAT_A2' ,a2 ,is_available, lsubmodel, unitab)
124 CALL hm_get_floatv('MAT_AMAX' ,amx ,is_available, lsubmodel, unitab)
125 CALL hm_get_floatv('EOS_COM_C0' ,c0 ,is_available, lsubmodel, unitab)
126 CALL hm_get_floatv('EOS_COM_C1' ,c1 ,is_available, lsubmodel, unitab)
127 CALL hm_get_floatv('EOS_COM_C2' ,c2 ,is_available, lsubmodel, unitab)
128 CALL hm_get_floatv('EOS_COM_C3' ,c3 ,is_available, lsubmodel, unitab)
129 CALL hm_get_floatv('EOS_COM_B' ,bunl ,is_available, lsubmodel, unitab)
130 CALL hm_get_floatv('eos_com_mue_max' ,XMUMX ,IS_AVAILABLE, LSUBMODEL, UNITAB)
131
132 CALL HM_GET_FLOATV('mat_pc' ,PMIN ,IS_AVAILABLE, LSUBMODEL, UNITAB)
133 CALL HM_GET_FLOATV('pext' ,PEXT ,IS_AVAILABLE, LSUBMODEL, UNITAB)
134
135!---
136 IF (RHOR == ZERO) RHOR=RHO0
137 PM(1) = RHOR
138 PM(89)= RHO0
139!
140 IF(COUNT == 5) THEN
141 IFORM= 1 !old format (embedded EoS : C0,C1,C2, ... parameters)
142 IEOS = 13
143 EOS_TAG(IEOS)%L_MU = 1
144 MATPARAM%IEOS = IEOS
145 ELSE
146 IFORM=2 !new format (Eos defined from /EOS option)
147 IEOS = 18 ! linear EOS is used by default
148 IREC = IREC+1
149 ENDIF
150!
151 MATPARAM%IEOS = IEOS
152 IPM(4) = IEOS ! keep this temporarily for output and fluid sections
153
154 E0 = ZERO
155
156 IF(C1 == ZERO)C1=THIRD*YOUNG/(ONE-TWO*ANU) !bulk modulus
157
158 IF(YOUNG <= ZERO)THEN
159 chain='young modulus must be defined '
160 chain2=''
161 CALL ANCMSG(MSGID=829, MSGTYPE=MSGERROR, ANMODE=ANINFO, I1=10, I2=MAT_ID, C1='error', C2=TITR, C3=chain, C4=chain2)
162 ENDIF
163
164 IF(ANU <= ZERO)THEN
165 chain='poisson ratio must be defined '
166 chain2=''
167 CALL ANCMSG(MSGID=829,MSGTYPE=MSGERROR,ANMODE=ANINFO,I1=10,I2=MAT_ID,C1='error',C2=TITR,C3=chain,C4=chain2)
168 ENDIF
169
170.AND. IF(A1 < ZERO A2 == ZERO)THEN
171 chain ='inverted yield surface because a1 is negative. '
172 chain2='check a0,a1,a2 yield parameters '
173 CALL ANCMSG(MSGID=829,MSGTYPE=MSGWARNING,ANMODE=ANINFO,I1=10,I2=MAT_ID,C1='warning',C2=TITR,C3=chain,C4=chain2)
174 ENDIF
175
176 IF(A2 < ZERO)THEN
177 chain ='unexpected yield surface : a2 is negative '
178 chain2='check a0,a1,a2 yield parameters '
179 CALL ANCMSG(MSGID=829,MSGTYPE=MSGWARNING,ANMODE=ANINFO,I1=10,I2=MAT_ID,C1='warning',C2=TITR,C3=chain,C4=chain2)
180 ENDIF
181.AND. IF(IFORM == 1 C1 <= ZERO) THEN
182 chain='tensile bulk modulus c1 is lower or equal to 0. '
183 chain2=''
184 CALL ANCMSG(MSGID=829,MSGTYPE=MSGERROR,ANMODE=ANINFO,I1=10,I2=MAT_ID,C1='error',C2=TITR,C3=chain,C4=chain2)
185 END IF
186
187.AND. IF(A2==ZEROA1/=ZERO)THEN !(A2=A1=ZERO => error message)
188 PSTAR=-A0/A1
189 ELSEIF(A2 /= ZERO)THEN
190 DELTA = A1*A1-FOUR*A0*A2
191 !Si intersection avec l'axe
192 IF(delta >= zero)THEN
193 delta=sqrt(delta)
194 pstar = (-a1+delta)/two/a2
195
196 ELSE
197 pstar = -a1/two/a2 !extremum
198 chain ='YIELD SURFACE J2(P)=A0+A1.P+A2.P^2 HAS NO INTERSECTION WITH '
199 chain2='PRESSURE AXIS. ASSUMING SURFACE CLOSURE AT P= '
200 WRITE(chain2(46:59),fmt=('(E13.6)')) pstar
201 CALL ancmsg( msgid=829,msgtype=msgwarning,anmode=aninfo,i1=10,i2=mat_id,c1='WARNING',c2=titr,c3=chain,c4=chain2)
202 ENDIF
203 ELSE
204 !do nothing let user do what he wants
205 ENDIF
206C
207 IF(amx == zero) amx = ep20
208 IF(pmin == zero) pmin =-infinity
209 !-----------------------------------------------------------------------------!
210 ! THESE FOLLOWING CHECKS ARE RELEVANT ONLY IN CASE OF EMBEDDED COMPACTION EOS !
211 !-----------------------------------------------------------------------------!
212
213 IF(iform == 1)THEN
214 IF(xmumx == zero.AND.bunl /= zero)THEN
215 chain= 'MISSING MUMAX VALUE IS AUTOMATICALLY ESTIMATED FROM BUNL VALUE. '
216 chain2=''
217 CALL ancmsg(msgid=829,msgtype=msgwarning,anmode=aninfo,i1=10,i2=mat_id,c1='WARNING',c2=titr,c3=chain,c4=chain2)
218 IF(c3 == zero)THEN
219 IF(c2 == zero)THEN
220 xmumx=ep20
221 ELSE
222 xmumx=(bunl-c1)/(two*c2)
223 ENDIF
224 ELSE
225 det=sqrt(c2**2 + three*c3*(bunl-c1))
226 xmumx=(det-c2)/(three*c3)
227 ENDIF
228 ENDIF
229
230 IF(xmumx /= zero.AND.bunl == zero)THEN
231 chain= 'MISSING BUNL VALUE IS AUTOMATICALLY ESTIMATED FROM MUMAX '
232 chain2=''
233 CALL ancmsg(msgid=829,msgtype=msgwarning,anmode=aninfo,i1=10,i2=mat_id,c1='WARNING',c2=titr,c3=chain,c4=chain2)
234 IF(c3 == zero)THEN
235 IF(c2 == zero)THEN
236 bunl = c1
237 ELSE
238 bunl = c1 + two*c2*xmumx
239 ENDIF
240 ELSE
241 bunl = c1 + two*c2*xmumx + three*c3*c3*xmumx**two
242 ENDIF
243 ENDIF
244 IF(bunl == zero) bunl = c1
245 ENDIF !IFORM == 1
246
247 g=young/(two*(one + anu))
248
249 !--- STORAGE OF MATERIAL PARAMETERS IN MATERIAL BUFFER ---!
250 pm(20)=young
251 pm(21)=anu
252 pm(22)=g
253 pm(23)=e0
254 pm(31)=c0
255 pm(32)=c1
256 pm(33)=c2
257 pm(34)=c3
258 pm(49)=c0
259 pm(35)=zero
260 pm(36)=zero
261 pm(37)=pmin
262 pm(38)=a0
263 pm(39)=a1
264 pm(40)=a2
265 pm(41)=amx
266 pm(43)=pext
267 pm(44)=pstar
268 pm(45)=bunl
269 pm(46)=xmumx
270 pm(47)=zero ! MUMIN
271 pm(48)=iform
272 pm(79)=three100
273 pm(88)=psh
274 pm(104)=c0
275 pm(105) = two*g/(c1+four_over_3*g) ! =(1-2*Nu)/(1-Nu)
276
277 !Formulation for solid elements time step computation.
278 ipm(252)= 2
279
280C---- Buffer Size for specific element buffer allocations
281
282 mtag%G_PLA = 1
283 mtag%G_EPSQ = 1
284 mtag%G_MU = 1 !unloading history (volumetric plastic strain)
285
286 mtag%L_PLA = 1 ! /TH(VPLA)
287 mtag%L_EPSQ = 1 ! /TH(EPSP)
288 mtag%L_MU = 1 ! unloading history (volumetric plastic strain)
289
290 IF(count == 5) THEN
291 !embedded compaction eos (old format)
292 matparam%eos%nuparam = 3
293 matparam%eos%niparam = 1
294 matparam%eos%nfunc = 0
295 matparam%eos%ntable = 0
296 call matparam%eos%construct() !allocations
297 matparam%eos%uparam(1) = xmumx
298 matparam%eos%uparam(2) = zero
299 matparam%eos%uparam(3) = bunl
300 matparam%eos%iparam(1) = iform
301 ENDIF
302
303 ! Material compatibility with /EOS option
304 CALL init_mat_keyword(matparam,"EOS")
305
306 ! EOS/Thermo keyword for pressure treatment in elements
307 CALL init_mat_keyword(matparam,"HYDRO_EOS")
308
309 CALL init_mat_keyword(matparam,"ELASTO_PLASTIC")
310
311 ! Properties compatibility
312 CALL init_mat_keyword(matparam,"SOLID_ISOTROPIC")
313 CALL init_mat_keyword(matparam,"SPH")
314
315 !--- OUTPUT IN STARTER LISTING FILE ---!
316 WRITE(iout,1100) trim(titr),mat_id,10
317 WRITE(iout,1000)
318 IF(is_encrypted)THEN
319 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA'
320 ELSE
321 WRITE(iout,1200)rho0,rhor
322 WRITE(iout,1300)young,anu,g
323 WRITE(iout,1400)a0,a1,a2,amx
324 IF(iform == 1)THEN
325 WRITE(iout,1500)c0,c1,c2,c3,bunl,xmumx,pmin,pstar
326 ELSE
327 WRITE(iout,1501)bunl,xmumx,pmin,pstar
328 ENDIF
329 ENDIF
330C-----------
331 RETURN
332 1000 FORMAT(
333 & 5x,40h drucker-prager(law10) ,/,
334 & 5x,40h ---------------------- ,//)
335 1100 FORMAT(/
336 & 5x,a,/,
337 & 5x,'MATERIAL NUMBER . . . . . . . . . . . .=',i10/,
338 & 5x,'MATERIAL LAW. . . . . . . . . . . . . .=',i10/)
339 1200 FORMAT(
340 & 5x,'INITIAL DENSITY . . . . . . . . . . . .=',1pg20.13/,
341 & 5x,'REFERENCE DENSITY . . . . . . . . . . .=',1pg20.13/)
342 1300 FORMAT(
343 & 5x,40hyoung'S MODULUS . . . . . . . . . . . .=,E12.4/,
344 & 5X,40HPOISSON's ratio . . . . . . . . . . . .=,e12.4/,
345 & 5x,40hshear modulus . . . . . . . . . . . . .=,e12.4//)
346 1400 FORMAT(
347 & 5x,40hyield coefficient a0. . . . . . . . . .=,e12.4/,
348 & 5x,40hyield coefficient a1. . . . . . . . . .=,e12.4/,
349 & 5x,40hyield coefficient a2. . . . . . . . . .=,e12.4/,
350 & 5x,40ha-max . . . . . . . . . . . . . . . . .=,e12.4//)
351 1500 FORMAT(
352 & 5x,40hc0. . . . . . . . . . . . . . . . . . .=,e12.4/,
353 & 5x,40hc1. . . . . . . . . . . . . . . . . . .=,e12.4/,
354 & 5x,40hc2. . . . . . . . . . . . . . . . . . .=,e12.4/,
355 & 5x,40hc3. . . . . . . . . . . . . . . . . . .=,e12.4/,
356 & 5x,40hunloading bulk. . . . . . . . . . . . .=,e12.4/,
357 & 5x,40hmax volumic compression . . . . . . . .=,e12.4/,
358 & 5x,40hfracture pressure . . . . . . . . . . .=,e12.4/,
359 & 5x,40hyield surface pressure root . . . . . .=,e12.4//)
360 1501 FORMAT(
361 & 5x,40hunloading bulk. . . . . . . . . . . . .=,e12.4/,
362 & 5x,40hmax volumic compression . . . . . . . .=,e12.4/,
363 & 5x,40hfracture pressure . . . . . . . . . . .=,e12.4/,
364 & 5x,40hyield surface pressure root . . . . . .=,e12.4//)
365C-----------
366 RETURN
367 END
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
hm_read_mat10
subroutine hm_read_mat10(lsubmodel, mtag, unitab, ipm, pm, mat_id, titr, eos_tag, israte, matparam)
Definition hm_read_mat10.F:42
init_mat_keyword
subroutine init_mat_keyword(matparam, keyword)
Definition init_mat_keyword.F:150
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
submodel_mod::nsubmod
integer nsubmod
Definition submodel_mod.F:720
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::eos_tag_
Definition elbuftag_mod.F:353
elbuftag_mod::mlaw_tag_
Definition elbuftag_mod.F:232
submodel_mod::submodel_data
Definition submodel_mod.F:723