OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_mat03.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_mat03 ../starter/source/materials/mat/mat003/hm_read_mat03.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!|| submodel_mod ../starter/share/modules1/submodel_mod.F
37!||====================================================================
38 SUBROUTINE hm_read_mat03(
39 . LSUBMODEL ,MTAG ,UNITAB ,IPM ,PM ,
40 . MAT_ID ,TITR ,ISRATE ,MATPARAM)
41C-----------------------------------------------
42C M o d u l e s
43C-----------------------------------------------
44 USE elbuftag_mod !MAXEOS
45 USE message_mod
46 USE submodel_mod
47 USE matparam_def_mod
48 USE unitab_mod
49 USE names_and_titles_mod , ONLY : nchartitle
50C-----------------------------------------------
51C ROUTINE DESCRIPTION :
52C ===================
53C READ MAT LAW03 WITH HM READER
54C-----------------------------------------------
55C DUMMY ARGUMENTS DESCRIPTION:
56C ===================
57C UNITAB UNITS ARRAY
58C MAT_ID MATERIAL ID(INTEGER)
59C TITR MATERIAL TITLE
60C LSUBMODEL SUBMODEL STRUCTURE
61C-----------------------------------------------
62C I m p l i c i t T y p e s
63C-----------------------------------------------
64#include "implicit_f.inc"
65C-----------------------------------------------
66C C o m m o n B l o c k s
67C-----------------------------------------------
68#include "scr03_c.inc"
69#include "units_c.inc"
70#include "param_c.inc"
71C-----------------------------------------------
72C D u m m y A r g u m e n t s
73C-----------------------------------------------
74 INTEGER, INTENT(INOUT) :: MAT_ID
75 INTEGER, INTENT(INOUT) :: ISRATE
76 INTEGER, DIMENSION(NPROPMI) ,INTENT(INOUT) :: IPM
77 CHARACTER(LEN=NCHARTITLE) ,INTENT(IN) :: TITR
78 my_real, DIMENSION(NPROPM) ,INTENT(INOUT) :: pm
79 TYPE(UNIT_TYPE_) ,INTENT(IN) :: UNITAB
80 TYPE(submodel_data), DIMENSION(NSUBMOD),INTENT(IN) :: LSUBMODEL
81 TYPE(mlaw_tag_), INTENT(INOUT) :: MTAG
82 TYPE(matparam_struct_),INTENT(INOUT) :: MATPARAM
83C-----------------------------------------------
84C L o c a l V a r i a b l e s
85C-----------------------------------------------
86 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED,IS_EOS, EOS_EMBEDDED
87 INTEGER JTUR, I, EMPTY_LINE,COUNT,IEOS,LAW6_OPT
88 my_real
89 . young, anu, ca, cb, cn, epsm, sigm, bulk,
90 . pmin, e0, g, fac_l, fac_t, fac_m, fac_c,
91 . bid2 ,bid3 ,bid4 ,bid5 ,
92 . bid2_,bid3_,bid4_,bid5_,
93 . c0,c1,c2,c3,c4,c5,psh,rho0,rhor
94C-----------------------------------------------
95C S o u r c e L i n e s
96C-----------------------------------------------
97 count = 0
98 is_eos=.false.
99 ieos=0
100 eos_embedded=.false.
101 bulk=zero
102 is_encrypted = .false.
103 is_available = .false.
104!---
105 CALL hm_option_is_encrypted(is_encrypted)
106
107 CALL hm_get_floatv('MAT_RHO' ,rho0 ,is_available, lsubmodel, unitab)
108 CALL hm_get_floatv('Refer_Rho' ,rhor ,is_available, lsubmodel, unitab)
109
110 CALL hm_get_floatv('MAT_E' ,young ,is_available, lsubmodel, unitab)
111 CALL hm_get_floatv('MAT_NU' ,anu ,is_available, lsubmodel, unitab)
112
113 CALL hm_get_floatv('MAT_SIGY' ,ca ,is_available, lsubmodel, unitab)
114 CALL hm_get_floatv('MAT_BETA' ,cb ,is_available, lsubmodel, unitab)
115 CALL hm_get_floatv('MAT_HARD' ,cn ,is_available, lsubmodel, unitab)
116 CALL hm_get_floatv('MAT_EPS' ,epsm ,is_available, lsubmodel, unitab)
117 CALL hm_get_floatv('MAT_SIG' ,sigm ,is_available, lsubmodel, unitab)
118
119 CALL hm_get_floatv('MAT_PC' ,pmin ,is_available, lsubmodel, unitab)
120
121 CALL hm_get_intv ('Line_count' ,count ,is_available, lsubmodel)
122
123 IF(invers_init>=2018)THEN
124 eos_embedded=.false.
125 ELSEIF(invers_init>=110)THEN
126 IF(count==3)eos_embedded=.true.
127 ELSE
128 eos_embedded=.true.
129 ENDIF
130
131 !----------------------------------------------------
132 !CHECK IF OLD CARD FORMAT 100 IS USED (OBSOLETE)
133 ! MAKES IT COMPATIBLE + WARNING
134 !----------------------------------------------------
135 IF(eos_embedded)THEN ! POLYNOMIAL EOS
136 IF(invers_init>=140) THEN
137 CALL ancmsg(msgid=1072,
138 . msgtype=msgerror,
139 . anmode=aninfo,
140 . i1 = mat_id,
141 . c1 = titr)
142 RETURN
143 ENDIF
144 CALL hm_get_floatv('MAT_C0' ,c0 ,is_available, lsubmodel, unitab)
145 CALL hm_get_floatv('MAT_C1' ,c1 ,is_available, lsubmodel, unitab)
146 CALL hm_get_floatv('MAT_C2' ,c2 ,is_available, lsubmodel, unitab)
147 CALL hm_get_floatv('MAT_C3' ,c3 ,is_available, lsubmodel, unitab)
148 CALL hm_get_floatv('MAT_C4' ,c4 ,is_available, lsubmodel, unitab)
149 CALL hm_get_floatv('MAT_C5' ,c5 ,is_available, lsubmodel, unitab)
150 CALL hm_get_floatv('MAT_E0' ,e0 ,is_available, lsubmodel, unitab)
151 CALL hm_get_floatv('MAT_PSH' ,psh ,is_available, lsubmodel, unitab)
152 bulk = c1
153 pm(23)= e0
154 pm(31)= c0-psh
155 pm(32)= c1
156 pm(33)= c2
157 pm(34)= c3
158 pm(35)= c4
159 pm(36)= c5
160 pm(88)= psh
161 ieos = 1 ! POLYNOMIAL EOS
162 ELSE
163 ieos = 18 ! linear EOS is used by default
164 ENDIF
165!
166 matparam%IEOS = ieos
167 ipm(4) = ieos ! keep this temporarily for output and fluid sections
168
169 IF (pmin==zero) pmin=-ep20
170 IF (rhor==zero) rhor=rho0
171 pm(1) = rhor
172 pm(89)= rho0
173 IF(cn==zero.OR.cn==one) cn = onep0001
174 IF(epsm==zero) epsm = ep20
175 IF(sigm==zero) sigm = ep20
176 g = young/(two*(one +anu))
177 IF(.NOT.eos_embedded)bulk = young/(three*(one - two*anu))
178 IF(.NOT.eos_embedded)e0 = zero
179
180 pm(20)=young
181 pm(21)=anu
182 pm(22)=g
183 pm(23)=e0
184 pm(32)=bulk
185 pm(37)=pmin
186 pm(38)=ca
187 pm(39)=cb
188 pm(40)=cn
189 pm(41)=epsm
190 pm(42)=sigm
191 pm(79)=three100
192 !Formulation for solid elements time step computation.
193 ipm(252)= 2
194 pm(105) = (one -two*anu)/(one - anu) ! TWO*G/(BULK + FOUR_OVER_3*G)
195
196 mtag%G_PLA = 1
197 mtag%L_PLA = 1
198 mtag%G_TEMP = 1
199 mtag%L_TEMP = 1
200
201 CALL init_mat_keyword(matparam,"ELASTO_PLASTIC")
202 CALL init_mat_keyword(matparam,"INCREMENTAL")
203 CALL init_mat_keyword(matparam,"LARGE_STRAIN")
204 ! Material compatibility with /EOS
205 CALL init_mat_keyword(matparam,"EOS")
206
207 ! EOS/Thermo keyword
208 CALL init_mat_keyword(matparam,"HYDRO_EOS")
209
210 ! Properties compatibility
211 CALL init_mat_keyword(matparam,"SOLID_ISOTROPIC")
212 CALL init_mat_keyword(matparam,"SPH")
213
214 WRITE(iout,2001) titr,mat_id,3
215 WRITE(iout,1000)
216
217 IF(is_encrypted)THEN
218 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA'
219 ELSE
220 WRITE(iout,2002)rho0,rhor
221 WRITE(iout,1300)young,anu,g,bulk
222 WRITE(iout,1400)ca,cb,cn,epsm,sigm
223 WRITE(iout,1500)pmin
224 ENDIF
225 !--------------------------------------------------!
226 ! OLD FORMAT 100 (OBSOLETE) !
227 !--------------------------------------------------!
228 IF(eos_embedded)THEN !
229 WRITE(iout,2000) !
230 IF(is_encrypted)THEN !
231 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA' !
232 ELSE !
233 WRITE(iout,2500)c0,c1,c2,c3,c4,c5,e0,pm(1),psh !
234 ENDIF !
235 ENDIF !
236 !--------------------------------------------------!
237C-----------
238 RETURN
239 1000 FORMAT(
240 & 5x,' ELASTIC PLASTIC HYDRODYNAMIC ',/,
241 & 5x,' ---------------------------- ',//)
242 1300 FORMAT(
243 & 5x,'YOUNG',1h','s modulus . . . . . . . . . . . .=',1PG20.13/,
244 & 5X,'poisson',1H','S RATIO . . . . . . . . . . . .=',1pg20.13/,
245 & 5x,'SHEAR MODULUS . . . . . . . . . . . . . .=',1pg20.13/,
246 & 5x,'BULK MODULUS . . . . . . . . . . . . . .=',1pg20.13//)
247 1400 FORMAT(
248 & 5x,'YIELD COEFFICIENT CA. . . . . . . . . .=',1pg20.13/,
249 & 5x,'YIELD COEFFICIENT CB. . . . . . . . . .=',1pg20.13/,
250 & 5x,'YIELD COEFFICIENT CN. . . . . . . . . .=',1pg20.13/,
251 & 5x,'EPS-MAX . . . . . . . . . . . . . . . .=',1pg20.13/,
252 & 5x,'SIG-MAX . . . . . . . . . . . . . . . .=',1pg20.13//)
253 1500 FORMAT(
254 & 5x,'PRESSURE CUTOFF . . . . . . . . . . . .=',1pg20.13//)
255 2000 FORMAT(
256 & 5x,' POLYNOMIAL EOS ',/,
257 & 5x,' -------------- ',/,
258 & 5x,' INPUT 100 IS OBSOLETE ',/,
259 & 5x,' IT CAN BE ERASED BY /EOS CARD ',/)
260 2500 FORMAT(
261 & 5x,'C0. . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
262 & 5x,'C1. . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
263 & 5x,'C2. . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
264 & 5x,'C3. . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
265 & 5x,'C4. . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
266 & 5x,'C5. . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
267 & 5x,'INITIAL INTERNAL ENERGY PER UNIT VOLUME .=',1pg20.13/,
268 & 5x,'REFERENCE DENSITY . . . . . . . . . . . .=',1pg20.13/,
269 & 5x,'PRESSURE SHIFT. . . . . . . . . . . . . .=',1pg20.13//)
270 2001 FORMAT(/
271 & 5x,a,/,
272 & 5x,' MATERIAL NUMBER . . . . . . . . . . . .=',i10/,
273 & 5x,' MATERIAL LAW. . . . . . . . . . . . . .=',i10/)
274 2002 FORMAT(
275 & 5x,' INITIAL DENSITY . . . . . . . . . . . .=',1pg20.13/,
276 & 5x,' REFERENCE DENSITY . . . . . . . . . . .=',1pg20.13/)
277
278 RETURN
279 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_mat03
subroutine hm_read_mat03(lsubmodel, mtag, unitab, ipm, pm, mat_id, titr, israte, matparam)
Definition hm_read_mat03.F:41
init_mat_keyword
subroutine init_mat_keyword(matparam, keyword)
Definition init_mat_keyword.F:150
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
starter
program starter
Definition starter.F:39
elbuftag_mod::mlaw_tag_
Definition elbuftag_mod.F:232
submodel_mod::submodel_data
Definition submodel_mod.F:723