OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_mat40.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_mat40 ../starter/source/materials/mat/mat040/hm_read_mat40.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_mat40(UPARAM ,MAXUPARAM, NUPARAM ,NUVAR ,
39 . MAXFUNC ,NFUNC , STIFINT ,UNITAB ,MAT_ID ,
40 . MTAG ,TITR , LSUBMODEL ,PM ,IMATVIS,
41 . MATPARAM )
42C-----------------------------------------------
43C M o d u l e s
44C-----------------------------------------------
45 USE unitab_mod
46 USE elbuftag_mod
47 USE message_mod
48 USE submodel_mod
49 USE matparam_def_mod
50 USE hm_option_read_mod
51 USE names_and_titles_mod , ONLY : nchartitle
52C-----------------------------------------------
53C I m p l i c i t T y p e s
54C-----------------------------------------------
55#include "implicit_f.inc"
56C-----------------------------------------------
57C C o m m o n B l o c k s
58C-----------------------------------------------
59#include "units_c.inc"
60#include "param_c.inc"
61C-----------------------------------------------
62C D u m m y A r g u m e n t s
63C-----------------------------------------------
64 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
65 my_real, DIMENSION(NPROPM), INTENT(INOUT) :: PM
66 my_real, DIMENSION(100), INTENT(INOUT) :: STIFINT
67 my_real, DIMENSION(MAXUPARAM), INTENT(INOUT) :: uparam
68 INTEGER, INTENT(INOUT) :: NFUNC, NUPARAM, NUVAR, IMATVIS
69 TYPE(mlaw_tag_), INTENT(INOUT) :: MTAG
70 INTEGER, INTENT(IN) :: MAT_ID, MAXFUNC, MAXUPARAM
71 CHARACTER(LEN=NCHARTITLE),INTENT(IN) :: TITR
72 TYPE(submodel_data), DIMENSION(NSUBMOD), INTENT(IN) :: LSUBMODEL
73 TYPE(matparam_struct_) ,INTENT(INOUT) :: MATPARAM
74C-----------------------------------------------
75C L o c a l V a r i a b l e s
76C-----------------------------------------------
77 my_real :: ak, g0, g1, g2, g3, g4, g5, gt, beta1,
78 . beta2, beta3, beta4, beta5, nu1, nu2,
79 . astas, bstas, vmisk, fac_l, fac_t, fac_m, fac_c,
80 . rho0, rhor
81 LOGICAL :: IS_ENCRYPTED, IS_AVAILABLE
82C-----------------------------------------------
83C S o u r c e L i n e s
84C-----------------------------------------------
85! Parameter initialization
86 nuvar = 40
87 nfunc = 0
88 imatvis = 1
89
90 is_encrypted = .false.
91 is_available = .false.
92
93! Check input encryption
94 CALL hm_option_is_encrypted(is_encrypted)
95! Initial and reference density
96 CALL hm_get_floatv('MAT_RHO' ,rho0 ,is_available, lsubmodel, unitab)
97 CALL hm_get_floatv('Refer_Rho' ,rhor ,is_available, lsubmodel, unitab)
98 IF (rhor == zero) THEN
99 rhor = rho0
100 ENDIF
101 pm(1) = rhor
102 pm(89) = rho0
103! Line 1
104 CALL hm_get_floatv('MAT_BULK', ak, is_available, lsubmodel, unitab)
105 CALL hm_get_floatv('MAT_GI', g0, is_available, lsubmodel, unitab)
106 CALL hm_get_floatv('Astass', astas, is_available, lsubmodel, unitab)
107 CALL hm_get_floatv('Bstass', bstas, is_available, lsubmodel, unitab)
108 CALL hm_get_floatv('Kvm', vmisk, is_available, lsubmodel, unitab)
109! Line 2
110 CALL hm_get_floatv('MAT_G0', g1, is_available, lsubmodel, unitab)
111 CALL hm_get_floatv('MAT_G2', g2, is_available, lsubmodel, unitab)
112 CALL hm_get_floatv('MAT_G3', g3, is_available, lsubmodel, unitab)
113 CALL hm_get_floatv('MAT_G4', g4, is_available, lsubmodel, unitab)
114 CALL hm_get_floatv('MAT_G5', g5, is_available, lsubmodel, unitab)
115! Line 3
116 CALL hm_get_floatv('MAT_DECAY', beta1, is_available, lsubmodel, unitab)
117 CALL hm_get_floatv('MAT_DECAY2', beta2, is_available, lsubmodel, unitab)
118 CALL hm_get_floatv('MAT_DECAY3', beta3, is_available, lsubmodel, unitab)
119 CALL hm_get_floatv('MAT_DECAY4', beta4, is_available, lsubmodel, unitab)
120 CALL hm_get_floatv('MAT_DECAY5', beta5, is_available, lsubmodel, unitab)
121
122 IF (astas <= em20) astas = infinity
123 IF (bstas <= em20) bstas = infinity
124 IF (vmisk <= em20) vmisk = infinity
125 nu1 = (three * ak - two * g0) / (two * g0 + six * ak)
126 gt = g0 + g1 + g2 + g3 + g4 + g5
127 nu2 = (three * ak - two * gt) / (two * gt + six * ak)
128 IF (nu1 < zero .OR. nu1 >= half) THEN
129 CALL ancmsg(msgid = 49,
130 . msgtype = msgerror,
131 . anmode = aninfo,
132 . r1 = nu1,
133 . i1 = mat_id,
134 . c1 = titr)
135 ENDIF
136 IF (nu2 < zero .OR. nu2 >= half) THEN
137 CALL ancmsg(msgid = 49,
138 . msgtype = msgerror,
139 . anmode = aninfo,
140 . r1 = nu2,
141 . i1 = mat_id,
142 . c1 = titr)
143 ENDIF
144 nuparam = 15
145 IF(nuparam > maxuparam)THEN
146 CALL ancmsg(msgid = 309,
147 . msgtype = msgerror,
148 . anmode = aninfo,
149 . i1 = mat_id,
150 . c1 = titr,
151 . i2 = nuparam,
152 . i3 = maxuparam)
153 ELSE
154 uparam(1) = ak
155 uparam(2) = g0
156 uparam(3) = g1
157 uparam(4) = g2
158 uparam(5) = g3
159 uparam(6) = g4
160 uparam(7) = g5
161 uparam(8) = max(beta1, em20)
162 uparam(9) = max(beta2, em20)
163 uparam(10) = max(beta3, em20)
164 uparam(11) = max(beta4, em20)
165 uparam(12) = max(beta5, em20)
166 uparam(13) = astas
167 uparam(14) = bstas
168 uparam(15) = vmisk
169 ENDIF
170
171 stifint(1) = ak
172
173! Formulation for solid elements time step computation.
174 stifint(16) = 2
175 stifint(17) = two * g0 / (ak + four_over_3 * g0)
176c-----------------
177 IF (nu1 >= 0.49 .or. nu2 >= 0.49) THEN
178 CALL init_mat_keyword(matparam,"INCOMPRESSIBLE")
179 ELSE
180 CALL init_mat_keyword(matparam,"COMPRESSIBLE")
181 END IF
182 CALL init_mat_keyword(matparam,"HOOK")
183 ! Properties compatibility
184 CALL init_mat_keyword(matparam,"SOLID_ISOTROPIC")
185 CALL init_mat_keyword(matparam,"SPH")
186c-----------------
187 WRITE(iout, 800) trim(titr), mat_id, 40
188 WRITE(iout,1000)
189 IF(is_encrypted)THEN
190 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA'
191 ELSE
192 WRITE(iout, 850) rho0
193 WRITE(iout,1100)ak,g0,g1,g2,g3,g4,g5,
194 . beta1,beta2,beta3,beta4,beta5,
195 . astas,bstas,vmisk
196 ENDIF
197C
198 800 FORMAT(/
199 & 5x,a,/,
200 & 5x,'MATERIAL NUMBER. . . . . . . . . . . . . . .=',i10/,
201 & 5x,'MATERIAL LAW . . . . . . . . . . . . . . . .=',i10/)
202 850 FORMAT(
203 & 5x,'INITIAL DENSITY . . . . . . . . . . . . . .=',1pg20.13/)
204 1000 FORMAT(
205 & 5x,' MAXWELL VISCO-ELASTIC LAW ',/,
206 & 5x,' ------------------------- ',//)
207 1100 FORMAT(
208 & 5x,'BULK MODULUS . . . . . . . . . . . . .=',1pg20.13/
209 & 5x,'LONG TIME SHEAR MODULUS . . . . . . . .=',1pg20.13/
210 & 5x,'SHEAR MODULUS 1 . . . . . . . . . . . .=',1pg20.13/
211 & 5x,'SHEAR MODULUS 2 . . . . . . . . . . . .=',1pg20.13/
212 & 5x,'SHEAR MODULUS 3 . . . . . . . . . . . .=',1pg20.13/
213 & 5x,'SHEAR MODULUS 4 . . . . . . . . . . . .=',1pg20.13/
214 & 5x,'SHEAR MODULUS 5 . . . . . . . . . . . .=',1pg20.13/
215 & 5x,'DECAY CONSTANT 1 . . . . . . . . . . .=',1pg20.13/
216 & 5x,'DECAY CONSTANT 2 . . . . . . . . . . .=',1pg20.13/
217 & 5x,'DECAY CONSTANT 3 . . . . . . . . . . .=',1pg20.13/
218 & 5x,'DECAY CONSTANT 4 . . . . . . . . . . .=',1pg20.13/
219 & 5x,'DECAY CONSTANT 5 . . . . . . . . . . .=',1pg20.13/
220 & 5x,'STASSI A COEFFICIENT . . . . . . . . .=',1pg20.13/
221 & 5x,'STASSI B COEFFICIENT . . . . . . . . .=',1pg20.13/
222 & 5x,'K VON MISES COEFFICIENT . . . . . . . =',1pg20.13//)
223C
224 RETURN
225 END
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:415
hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:284
hm_read_mat40
subroutine hm_read_mat40(uparam, maxuparam, nuparam, nuvar, maxfunc, nfunc, stifint, unitab, mat_id, mtag, titr, lsubmodel, pm, imatvis, matparam)
Definition hm_read_mat40.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
hm_option_read_mod
Definition hm_option_read_mod.F:408
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