OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_mat63.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_mat63 ../starter/source/materials/mat/mat063/hm_read_mat63.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_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
31!|| init_mat_keyword ../starter/source/materials/mat/init_mat_keyword.F
32!||--- uses -----------------------------------------------------
33!|| elbuftag_mod ../starter/share/modules1/elbuftag_mod.F
34!|| hm_option_read_mod ../starter/share/modules1/hm_option_read_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_mat63(UPARAM ,MAXUPARAM,NUPARAM ,NUVAR ,
39 . PARMAT ,MAT_ID ,PM ,TITR ,
40 . UNITAB ,LSUBMODEL,MTAG ,MATPARAM )
41C-----------------------------------------------
42C ROUTINE DESCRIPTION :
43C ===================
44C READ MAT LAW63 WITH HM READER
45C-----------------------------------------------
46C DUMMY ARGUMENTS DESCRIPTION:
47C ===================
48C UNITAB UNITS ARRAY
49C MAT_ID MATERIAL ID(INTEGER)
50C TITR MATERIAL TITLE
51C LSUBMODEL SUBMODEL STRUCTURE
52C-----------------------------------------------
53C M o d u l e s
54C-----------------------------------------------
55 USE elbuftag_mod
56 USE message_mod
57 USE submodel_mod
58 USE matparam_def_mod
59 USE hm_option_read_mod
60 USE unitab_mod
61 USE names_and_titles_mod , ONLY : nchartitle
62C-----------------------------------------------
63C I m p l i c i t T y p e s
64C-----------------------------------------------
65#include "implicit_f.inc"
66C-----------------------------------------------
67C C o m m o n B l o c k s
68C-----------------------------------------------
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(IN) :: MAT_ID,MAXUPARAM
75 TYPE(UNIT_TYPE_), INTENT(IN) :: UNITAB
76 my_real, DIMENSION(NPROPM) ,INTENT(INOUT) :: pm
77 CHARACTER(LEN=NCHARTITLE) ,INTENT(IN) :: TITR
78 INTEGER, INTENT(INOUT) :: NUPARAM,NUVAR
79 my_real, DIMENSION(MAXUPARAM) ,INTENT(INOUT) :: uparam
80 my_real, DIMENSION(100),INTENT(INOUT) :: parmat
81 TYPE(submodel_data), DIMENSION(*),INTENT(IN) :: LSUBMODEL
82 TYPE(mlaw_tag_), INTENT(INOUT) :: MTAG
83 TYPE(matparam_struct_) ,INTENT(INOUT) :: MATPARAM
84C-----------------------------------------------
85C L o c a l V a r i a b l e s
86C-----------------------------------------------
87 my_real
88 . e,nu,cp,a, b, q,c,d,p, ahs,bhs,m,n,k1,k2,dh,vm0,eps0,
89 . c1,g,temp,hl,coef,rho0,rhor
90 INTEGER ILAW
91 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
92C=======================================================================
93 is_encrypted = .false.
94 is_available = .false.
95c
96 ilaw = 63
97c------------------------------------------
98 CALL hm_option_is_encrypted(is_encrypted)
99c------------------------------------------
100C
101 CALL hm_get_floatv('MAT_RHO' ,rho0 ,is_available, lsubmodel, unitab)
102 CALL hm_get_floatv('Refer_Rho',rhor ,is_available, lsubmodel, unitab)
103Card1
104 CALL hm_get_floatv('MAT_E' ,e ,is_available, lsubmodel, unitab)
105 CALL hm_get_floatv('MAT_NU' ,nu ,is_available, lsubmodel, unitab)
106 CALL hm_get_floatv('CP' ,cp ,is_available, lsubmodel, unitab)
107Card2
108 CALL hm_get_floatv('MAT_A' ,a ,is_available, lsubmodel, unitab)
109 CALL hm_get_floatv('MAT_B' ,b ,is_available, lsubmodel, unitab)
110 CALL hm_get_floatv('Q' ,q ,is_available, lsubmodel, unitab)
111 CALL hm_get_floatv('MAT_C' ,c ,is_available, lsubmodel, unitab)
112 CALL hm_get_floatv('MAT_D' ,d ,is_available, lsubmodel, unitab)
113Card3
114 CALL hm_get_floatv('MAT_P0' ,p ,is_available, lsubmodel, unitab)
115 CALL hm_get_floatv('A_HS' ,ahs ,is_available, lsubmodel, unitab)
116 CALL hm_get_floatv('B_HS' ,bhs ,is_available, lsubmodel, unitab)
117 CALL hm_get_floatv('MAT_M' ,m ,is_available, lsubmodel, unitab)
118 CALL hm_get_floatv('MAT_N' ,n ,is_available, lsubmodel, unitab)
119Card4
120 CALL hm_get_floatv('K1' ,k1 ,is_available, lsubmodel, unitab)
121 CALL hm_get_floatv('K2' ,k2 ,is_available, lsubmodel, unitab)
122 CALL hm_get_floatv('delta_H' ,dh ,is_available, lsubmodel, unitab)
123 CALL hm_get_floatv('Kvm' ,vm0 ,is_available, lsubmodel, unitab)
124 CALL hm_get_floatv('Epsilon_0',eps0 ,is_available, lsubmodel, unitab)
125Card5
126 CALL hm_get_floatv('MAT_T0' ,temp ,is_available, lsubmodel, unitab)
127 CALL hm_get_floatv('MAT_HL' ,hl ,is_available, lsubmodel, unitab)
128 CALL hm_get_floatv('MAT_ETA' ,coef ,is_available, lsubmodel, unitab)
129C
130c-----------------------------------------------------------------------
131c Default values
132c-----------------------------------------------------------------------
133 IF(coef <= zero .OR. coef > one) coef = one
134 IF(vm0 == zero .OR. vm0 >= one) vm0=em20
135 IF (cp == zero) THEN
136 cp = zero
137 ELSE
138 cp = one/cp
139 ENDIF
140C
141 IF (b == 0) b = -one
142 IF ((b+1)/b>=p) THEN
143 CALL ancmsg(msgid=301,
144 . msgtype=msgerror,
145 . anmode=aninfo,
146 . i1=63,
147 . i2=mat_id,
148 . c1=titr)
149 END IF
150c-----------------------------------------------------------------------
151 g = half*e/(one+nu)
152 c1 = e/three/(one - two*nu)
153c-----------------------------------------------------------------------
154 uparam(1) = e
155 uparam(2) = e/(one - nu*nu)
156 uparam(3) = nu*uparam(2)
157 uparam(4) = g
158 uparam(5) = nu
159 uparam(6) = a
160 uparam(7) = b
161 uparam(8) = q
162 uparam(9) = c
163 uparam(10) = d
164 uparam(11) = p
165 uparam(12) = ahs
166 uparam(13) = bhs
167 uparam(14) = m
168 uparam(15) = n
169 uparam(16) = k1
170 uparam(17) = k2
171 uparam(18) = dh
172 uparam(19) = vm0
173 uparam(20) = eps0
174 uparam(21) = cp
175 uparam(22) = temp
176 uparam(23) = hl
177 uparam(24) = coef
178c-----------------------------------------------------------------------
179 nuparam = 24
180 nuvar = 4
181c-----------------------------------------------------------------------
182 parmat(1) = c1
183 parmat(2) = e
184 parmat(3) = nu
185C Formulation for solid elements time step computation.
186 parmat(16) = 2
187 parmat(17) = (one - two*nu)/(one - nu) ! == TWO*G/(C1+FOUR_OVER_3*G)
188c-----------------------------------------------------------------------
189 pm(1) = rhor
190 pm(89) = rho0
191 mtag%G_PLA = 1
192 mtag%L_PLA = 1
193c-----------------------------------------------------------------------
194 ! Properties compatibility
195 CALL init_mat_keyword(matparam,"shell_isotropic")
196c-----------------------------------------------------------------------
197c--------------------------------------------------
198c Starter output
199c--------------------------------------------------
200 WRITE(IOUT,1000) TRIM(TITR),MAT_ID,ILAW
201 WRITE(IOUT,1100)
202 IF (IS_ENCRYPTED) THEN
203 WRITE(IOUT,'(5X,A,//)')'CONFIDENTIAL DATA'
204 ELSE
205 WRITE(IOUT,1200) RHO0
206 WRITE(IOUT,1300) E,NU,ONE/CP
207 WRITE(IOUT,1400) A,B,Q,C,D
208 WRITE(IOUT,1500) P,AHS,BHS,M,N
209 WRITE(IOUT,1600) K1,K2,DH,VM0,EPS0
210 WRITE(IOUT,1700) TEMP,HL,COEF
211 ENDIF
212c-----------------------------------------------------------------------
213 1000 FORMAT(/
214 & 5X,A,/,
215 & 5X,'MATERIAL NUMBER. . . . . . . . . . . . . . .=',I10/,
216 & 5X,'MATERIAL LAW . . . . . . . . . . . . . . . .=',I10/)
217 1100 FORMAT
218 &(5X,'MATERIAL : TRIP-STEELS PLASTIC LAW',/,
219 & 5X,'-----------------------------------',/)
220 1200 FORMAT(
221 & 5X,'INITIAL DENSITY . . . . . . . . . . . . . .=',1PG20.13/)
222 1300 FORMAT(
223 & 5X,'YOUNG''S MODULUS . . . . . . . . . . . . . .=',1PG20.13/,
224 & 5X,'POISSON''S RATIO . . . . . . . . . . . . . .=',1PG20.13/,
225 & 5X,'THERMAL MASSIC CAPACITY . . . . . . . . . .=',1PG20.13/)
226 1400 FORMAT(
227 & 5X,'PARAMETER 1 (A) . . . . . . . . . . . . . .=',1PG20.13/,
228 & 5X,'PARAMETER 2 (B) . . . . . . . . . . . . . .=',1PG20.13/,
229 & 5X,'PARAMETER 3 (Q) . . . . . . . . . . . . . .=',1PG20.13/,
230 & 5X,'PARAMETER 4 (C) . . . . . . . . . . . . . .=',1PG20.13/,
231 & 5X,'PARAMETER 5 (D) . . . . . . . . . . . . . .=',1PG20.13/)
232 1500 FORMAT(
233 & 5X,'PARAMETER 6 (P) . . . . . . . . . . . . . .=',1PG20.13/,
234 & 5X,'PARAMETER 7 (AHS) . . . . . . . . . . . . .=',1PG20.13/,
235 & 5X,'PARAMETER 8 (BHS) . . . . . . . . . . . . .=',1PG20.13/,
236 & 5X,'PARAMETER 9 (M) . . . . . . . . . . . . . .=',1PG20.13/,
237 & 5X,'PARAMETER 10 (N) . . . . . . . . . . . . . .=',1PG20.13/)
238 1600 FORMAT(
239 & 5X,'PARAMETER 11 (K1) . . . . . . . . . . . . .=',1PG20.13/,
240 & 5X,'PARAMETER 12 (K2). . . . . . . . . . . . . .=',1PG20.13/,
241 & 5X,'PARAMETER 13 (DH). . . . . . . . . . . . . .=',1PG20.13/,
242 & 5X,'INITIAL MARTENSITE VOLUME FRACTION . . . . .=',1PG20.13/,
243 & 5X,'INITIAL PLASTIC STRAIN . . . . . . . . . . .=',1PG20.13/)
244 1700 FORMAT(
245 & 5X,'INITIAL TEMPERATURE . . . . . . . . . . . .=',1PG20.13/
246 & 5X,'LATENT HEAT OF MARTENSITE . . . . . . . . .=',1PG20.13/
247 & 5X,'TAYLOR-QUINNEY COEFFICIENT . . . . . . . . .=',1PG20.13/)
248 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:418
hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:288
hm_read_mat63
subroutine hm_read_mat63(uparam, maxuparam, nuparam, nuvar, parmat, mat_id, pm, titr, unitab, lsubmodel, mtag, matparam)
Definition hm_read_mat63.F:41
init_mat_keyword
subroutine init_mat_keyword(matparam, keyword)
Definition init_mat_keyword.F:153
elbuftag_mod
Definition elbuftag_mod.F:230
hm_option_read_mod
Definition hm_option_read_mod.F:412
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
starter
program starter
Definition starter.F:39
elbuftag_mod::mlaw_tag_
Definition elbuftag_mod.F:235
submodel_mod::submodel_data
Definition submodel_mod.F:728