OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_mat74.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_mat74 ../starter/source/materials/mat/mat074/hm_read_mat74.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_floatv_dim ../starter/source/devtools/hm_reader/hm_get_floatv_dim.F
31!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
32!|| hm_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
33!|| init_mat_keyword ../starter/source/materials/mat/init_mat_keyword.F
34!||--- uses -----------------------------------------------------
35!|| elbuftag_mod ../starter/share/modules1/elbuftag_mod.F
36!|| message_mod ../starter/share/message_module/message_mod.F
37!|| submodel_mod ../starter/share/modules1/submodel_mod.F
38!||====================================================================
39 SUBROUTINE hm_read_mat74(UPARAM ,MAXUPARAM,NUPARAM ,ISRATE , IMATVIS ,
40 . NUVAR ,IFUNC ,MAXFUNC ,NFUNC , PARMAT ,
41 . UNITAB ,MAT_ID ,TITR ,MTAG , LSUBMODEL,
42 . ITABLE ,MAXTABL ,NTABLE ,PM ,IPM ,
43 . MATPARAM)
44C-----------------------------------------------
45C D e s c r i p t i o n
46C-----------------------------------------------
47C READ MAT LAW74 WITH HM READER
48C
49C DUMMY ARGUMENTS DESCRIPTION:
50C ===================
51C
52C NAME DESCRIPTION
53C
54C PM MATERIAL ARRAY(REAL)
55C UNITAB UNITS ARRAY
56C MAT_ID MATERIAL ID(INTEGER)
57C TITR MATERIAL TITLE
58C LSUBMODEL SUBMODEL STRUCTURE
59C
60C-----------------------------------------------
61C M o d u l e s
62C-----------------------------------------------
63 USE unitab_mod
64 USE elbuftag_mod
65 USE message_mod
66 USE submodel_mod
67 USE matparam_def_mod
68 USE names_and_titles_mod , ONLY : nchartitle
69C-----------------------------------------------
70C I m p l i c i t T y p e s
71C-----------------------------------------------
72#include "implicit_f.inc"
73C-----------------------------------------------
74C C o m m o n B l o c k s
75C-----------------------------------------------
76#include "units_c.inc"
77#include "param_c.inc"
78C-----------------------------------------------
79C D u m m y A r g u m e n t s
80C-----------------------------------------------
81 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
82 my_real, INTENT(INOUT) :: PM(NPROPM),PARMAT(100),UPARAM(MAXUPARAM)
83 INTEGER, INTENT(INOUT) :: IPM(NPROPMI),ISRATE,IFUNC(MAXFUNC),NFUNC,MAXFUNC,MAXUPARAM,NUPARAM, NUVAR,IMATVIS
84 . ,NTABLE,ITABLE(MAXTABL)
85 TYPE(mlaw_tag_),INTENT(INOUT) :: MTAG
86 INTEGER,INTENT(IN) :: MAT_ID,MAXTABL
87 CHARACTER(LEN=NCHARTITLE) ,INTENT(IN) :: TITR
88
89 TYPE(submodel_data),INTENT(IN) :: LSUBMODEL(*)
90 TYPE(matparam_struct_) ,INTENT(INOUT) :: MATPARAM
91C-----------------------------------------------
92C L o c a l V a r i a b l e s
93C-----------------------------------------------
94 INTEGER ILAW
95 my_real
96 . E,NU,G,C1,EPSMAX,EPSR1,EPSR2,X2FAC,YFAC,
97 . S11Y,S22Y,S33Y,S12Y,S23Y,S31Y,
98 . fisokin,
99 . t0, rhocp, ff, gg, hh, ll, mm, nn, fcut,
100 . einf,ce,fac_pres,fac_strain_rate
101 INTEGER IFUNCE, OPTE
102 my_real :: RHO0, RHOR
103 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
104C-----------------------------------------------
105C S o u r c e L i n e s
106C-----------------------------------------------
107 is_encrypted = .false.
108 is_available = .false.
109 mtag%G_EPSD = 1
110 mtag%G_PLA = 1
111 mtag%G_SEQ = 1
112 mtag%L_EPSD = 1
113 mtag%L_PLA = 1
114 mtag%L_TEMP = 1
115!
116 israte=0
117 imatvis=0
118 ilaw=74
119C
120 CALL hm_option_is_encrypted(is_encrypted)
121 !line+1
122 CALL hm_get_floatv('MAT_RHO' ,rho0 ,is_available, lsubmodel, unitab)
123 !line-2
124 CALL hm_get_floatv('MAT_E' ,e ,is_available, lsubmodel, unitab)
125 CALL hm_get_floatv('MAT_NU' ,nu ,is_available, lsubmodel, unitab)
126 CALL hm_get_floatv('MAT_EPS' ,epsmax ,is_available, lsubmodel, unitab)
127 CALL hm_get_floatv('MAT_EPST1' ,epsr1 ,is_available, lsubmodel, unitab)
128 CALL hm_get_floatv('MAT_EPST2' ,epsr2 ,is_available, lsubmodel, unitab)
129 !line-3
130 CALL hm_get_intv ('Yr_fun' ,ifunce ,is_available, lsubmodel)
131 CALL hm_get_floatv('MAT_EFIB' ,einf ,is_available, lsubmodel, unitab)
132 CALL hm_get_floatv('MAT_C' ,ce ,is_available, lsubmodel, unitab)
133 !line-4
134 CALL hm_get_intv ('Fsmooth' ,israte ,is_available, lsubmodel)
135 CALL hm_get_floatv('MAT_HARD' ,fisokin ,is_available, lsubmodel, unitab)
136 CALL hm_get_floatv('Fcut' ,fcut ,is_available, lsubmodel, unitab)
137 !line-5
138 CALL hm_get_floatv('MAT_SIGT1' ,s11y ,is_available, lsubmodel, unitab)
139 CALL hm_get_floatv('MAT_SIGT2' ,s22y ,is_available, lsubmodel, unitab)
140 CALL hm_get_floatv('MAT_SIGT3' ,s33y ,is_available, lsubmodel, unitab)
141 !line-6
142 CALL hm_get_floatv('MAT_SIGYT1' ,s12y ,is_available, lsubmodel, unitab)
143 CALL hm_get_floatv('MAT_SIGYT2' ,s23y ,is_available, lsubmodel, unitab)
144 CALL hm_get_floatv('MAT_SIGYT3' ,s31y ,is_available, lsubmodel, unitab)
145 !line-7
146 CALL hm_get_intv ('FUN_A1' ,itable(1) ,is_available, lsubmodel)
147 CALL hm_get_floatv('MAT_FScale' ,yfac ,is_available, lsubmodel, unitab)
148 CALL hm_get_floatv('MAT_PScale' ,x2fac ,is_available, lsubmodel, unitab)
149 !line-8
150 CALL hm_get_floatv('T_Initial' ,t0 ,is_available, lsubmodel, unitab)
151 CALL hm_get_floatv('MAT_SPHEAT' ,rhocp ,is_available, lsubmodel, unitab)
152 !unit
153 CALL hm_get_floatv_dim('MAT_FScale' ,fac_pres ,is_available, lsubmodel, unitab)
154 CALL hm_get_floatv_dim('MAT_PScale' ,fac_strain_rate ,is_available, lsubmodel, unitab)
155 !========== DEFAULT VALUES=============!
156 rhor = rho0
157 pm(1) =rhor
158 pm(89)=rho0
159!!
160 IF(yfac == zero)yfac=fac_pres
161 IF(x2fac == zero)x2fac=fac_strain_rate
162C
163C tab(1st entry = plastic strain, 2nd entry = strain rate, 3rd entry = temperature)
164C
165 IF(epsmax == zero) epsmax=infinity
166 IF(epsr1 == zero) epsr1=infinity
167 IF(epsr2 == zero) epsr2=two*infinity
168 IF(t0 == zero) t0=twohundred93
169C
170 ntable=1
171 uparam(1)=fisokin
172 uparam(2)=e
173 uparam(3)=e/(one-nu*nu)
174 uparam(4)=nu*uparam(3)
175 g=half*e/(1.+nu)
176 uparam(5)=g
177 uparam(6)=nu
178
179C----------- FF...NN --------
180 IF (s11y <= zero) THEN
181 CALL ancmsg(msgid=822,
182 . msgtype=msgerror,
183 . anmode=aninfo,
184 . i1=mat_id,
185 . c1=titr,
186 . c2='S11Y')
187 ENDIF
188 IF (s22y <= zero) THEN
189 CALL ancmsg(msgid=822,
190 . msgtype=msgerror,
191 . anmode=aninfo,
192 . i1=mat_id,
193 . c1=titr,
194 . c2='S22Y')
195 ENDIF
196 IF (s33y <= zero) THEN
197 CALL ancmsg(msgid=822,
198 . msgtype=msgerror,
199 . anmode=aninfo,
200 . i1=mat_id,
201 . c1=titr,
202 . c2='S33Y')
203 ENDIF
204 IF (s12y <= zero) THEN
205 CALL ancmsg(msgid=822,
206 . msgtype=msgerror,
207 . anmode=aninfo,
208 . i1=mat_id,
209 . c1=titr,
210 . c2='S12Y')
211 ENDIF
212 IF (s23y <= zero) THEN
213 CALL ancmsg(msgid=822,
214 . msgtype=msgerror,
215 . anmode=aninfo,
216 . i1=mat_id,
217 . c1=titr,
218 . c2='S23Y')
219 ENDIF
220 IF (s31y <= zero) THEN
221 CALL ancmsg(msgid=822,
222 . msgtype=msgerror,
223 . anmode=aninfo,
224 . i1=mat_id,
225 . c1=titr,
226 . c2='s31y')
227 ENDIF
228C
229 FF=HALF*(ONE/S22Y**2+ONE/S33Y**2-ONE/S11Y**2)
230 GG=HALF*(ONE/S11Y**2+ONE/S33Y**2-ONE/S22Y**2)
231 HH=HALF*(ONE/S11Y**2+ONE/S22Y**2-ONE/S33Y**2)
232 LL=HALF/S23Y**2
233 MM=HALF/S31Y**2
234 NN=HALF/S12Y**2
235 UPARAM(7)=FF
236 UPARAM(8)=GG
237 UPARAM(9)=HH
238 UPARAM(10)=LL
239 UPARAM(11)=MM
240 UPARAM(12)=NN
241 UPARAM(13)=ONE/X2FAC
242 UPARAM(14)=YFAC
243 UPARAM(15)=EPSMAX
244 UPARAM(16)=EPSR1
245 UPARAM(17)=EPSR2
246 UPARAM(18)=TWO*G
247 UPARAM(19)=THREE*G
248 C1=E/THREE/(ONE - TWO*NU)
249 UPARAM(20)=C1
250 UPARAM(21)=C1+ FOUR_OVER_3*G
251 UPARAM(22)=T0
252 IF(RHOCP==ZERO)THEN
253 UPARAM(23)=ZERO
254 ELSE
255 UPARAM(23)=ONE/RHOCP
256 END IF
257c--------
258 OPTE = 0
259 IF (IFUNCE > 0 )OPTE = 1
260 NFUNC = 1
261 IFUNC(NFUNC) = IFUNCE
262 UPARAM(24) = NFUNC
263 UPARAM(25) = OPTE
264 UPARAM(26) = EINF
265 UPARAM(27) = CE
266c-------
267 NUPARAM =27
268C
269 ISRATE = 1
270 PARMAT(1)=C1
271 PARMAT(2)=E
272 PARMAT(3)=NU
273 PARMAT(4) = ISRATE
274 PARMAT(5) = FCUT
275C
276 NUVAR = 10
277CC Formulation for solid elements time step computation.
278 PARMAT(16) = 2
279 PARMAT(17) = (ONE - TWO*NU)/(ONE - NU) ! == TWO*G/(C1+FOUR_OVER_3*G)
280C
281 ! activate heat source calculation in material
282 MATPARAM%HEAT_FLAG = 1
283!
284 ! MATPARAM keywords
285 CALL INIT_MAT_KEYWORD(MATPARAM,"ORTHOTROPIC")
286C
287 ! Properties compatibility
288 CALL INIT_MAT_KEYWORD(MATPARAM,"SOLID_ORTHOTROPIC")
289C
290 WRITE(IOUT,1001) TRIM(TITR),MAT_ID,ILAW
291 WRITE(IOUT,1000)
292 IF(IS_ENCRYPTED)THEN
293 WRITE(IOUT,'(5x,a,//)')'confidential data'
294 ELSE
295 WRITE(IOUT,1002)RHO0
296 WRITE(IOUT,1100)E,NU,G,S11Y,S22Y,S33Y,S12Y,S23Y,S31Y,FISOKIN,
297 . ISRATE,FCUT
298 WRITE(IOUT,1300)EPSMAX,EPSR1,EPSR2
299 WRITE(IOUT,1200)ITABLE(1),X2FAC,YFAC
300 WRITE(IOUT,1400)T0,RHOCP,IFUNCE,EINF,CE
301 WRITE(IOUT,*)' '
302 ENDIF
303C
304 RETURN
305 1000 FORMAT(
306 & 5X,50H TABULATED HILL ORTHOTROPIC PLASTICITY FOR SOLIDS,/,
307 & 5X,50H ------------------------------------------------,//)
308 1001 FORMAT(
309 & 5X,A,/,
310 & 5X,'material number . . . . . . . . . . . .=',I10/,
311 & 5X,'material law. . . . . . . . . . . . . .=',I10/)
312 1002 FORMAT(
313 & 5X,'initial density . . . . . . . . . . . .=',1PG20.13/)
314 1100 FORMAT(
315 & 5X,'young modulus . . . . . . . . . . . . .=',1PG20.13/
316 & 5X,'poisson ratio . . . . . . . . . . . . .=',1PG20.13/
317 & 5X,'shear modulus . . . . . . . . . . . . .=',1PG20.13/
318 & 5X,'yield in direction 1. . . . . . . . . .=',1PG20.13/
319 & 5X,'yield in direction 2. . . . . . . . . .=',1PG20.13/
320 & 5X,'yield in direction 3. . . . . . . . . .=',1PG20.13/
321 & 5X,'yield in shear direction 12 . . . . . .=',1PG20.13/
322 & 5X,'yield in shear direction 23 . . . . . .=',1PG20.13/
323 & 5X,'yield in shear direction 31 . . . . . .=',1PG20.13/
324 & 5X,'iso-kinematic hardenning factor . . . .=',1PG20.13/
325 & 5X,'smooth strain rate option . . . . . . .=',I10/
326 & 5X,'strain rate cutting frequency . . . . .=',1PG20.13/)
327 1200 FORMAT(
328 & 5X,'yield stress table number . . . . . . .=',I10/
329 & 5X,'2nd entry(strain rate) scale factor. .=',1PG20.13/
330 & 5X,'yield scale factor. . . . . . . . . . .=',1PG20.13)
331 1300 FORMAT(
332 & 5X,'maximum plastic strain. . . . . . . . .=',1PG20.13/
333 & 5X,'tensile failure strain 1. . . . . . . .=',1PG20.13/
334 & 5X,'tensile failure strain 2. . . . . . . .=',1PG20.13)
335 1400 FORMAT(
336 & 5X,'initial temperature . . . . . . . . . .=',1PG20.13/
337 & 5X,'heat capacity per volume unit . . . . .=',1PG20.13/
338 & 5X,'young modulus scale factor FUNCTION . .=',I10/
339 & 5X,'young modulus einf. . . . . . . . . . .=',1PG20.13/
340 & 5X,'parameter ce. . . . . . . . . . . . . .=',1PG20.13)
341 END
342
343
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:418
hm_get_floatv_dim
subroutine hm_get_floatv_dim(name, dim_fac, is_available, lsubmodel, unitab)
Definition hm_get_floatv_dim.F:161
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_mat74
subroutine hm_read_mat74(uparam, maxuparam, nuparam, israte, imatvis, nuvar, ifunc, maxfunc, nfunc, parmat, unitab, mat_id, titr, mtag, lsubmodel, itable, maxtabl, ntable, pm, ipm, matparam)
Definition hm_read_mat74.F:44
elbuftag_mod
Definition elbuftag_mod.F:230
message_mod
Definition message_mod.F:1257
names_and_titles_mod
Definition names_and_titles_mod.F:1002
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1008
submodel_mod
Definition submodel_mod.F:717
unitab_mod
Definition unitab_mod.F:517
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:895
elbuftag_mod::mlaw_tag_
Definition elbuftag_mod.F:235
submodel_mod::submodel_data
Definition submodel_mod.F:728