OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_mat42.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_mat42 ../starter/source/materials/mat/mat042/hm_read_mat42.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_float_array_index ../starter/source/devtools/hm_reader/hm_get_float_array_index.F
30!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
31!|| hm_get_floatv_dim ../starter/source/devtools/hm_reader/hm_get_floatv_dim.F
32!|| hm_get_intv ../starter/source/devtools/hm_reader/hm_get_intv.F
33!|| hm_option_is_encrypted ../starter/source/devtools/hm_reader/hm_option_is_encrypted.F
34!|| init_mat_keyword ../starter/source/materials/mat/init_mat_keyword.F
35!||--- uses -----------------------------------------------------
36!|| message_mod ../starter/share/message_module/message_mod.F
37!|| submodel_mod ../starter/share/modules1/submodel_mod.F
38!||====================================================================
39 SUBROUTINE hm_read_mat42(
40 . MATPARAM ,NUVAR ,MAXFUNC ,NFUNC ,IFUNC ,
41 . PARMAT ,IMATVIS ,UNITAB ,LSUBMODEL,ID ,
42 . TITR ,PM )
43C-----------------------------------------------
44C M o d u l e s
45C-----------------------------------------------
46 USE unitab_mod
47 USE message_mod
48 USE submodel_mod
49 USE matparam_def_mod
50 USE names_and_titles_mod , ONLY : nchartitle
51C-----------------------------------------------
52C ROUTINE DESCRIPTION :
53C ===================
54C READ MAT LAW42 WITH HM READER
55C-----------------------------------------------
56C I m p l i c i t T y p e s
57C-----------------------------------------------
58#include "implicit_f.inc"
59C-----------------------------------------------
60C C o m m o n B l o c k s
61C-----------------------------------------------
62#include "units_c.inc"
63#include "param_c.inc"
64C-----------------------------------------------
65C D u m m y A r g u m e n t s
66C-----------------------------------------------
67 INTEGER, INTENT(IN) :: ID
68 INTEGER, INTENT(IN) :: MAXFUNC
69 INTEGER, INTENT(INOUT) :: NUVAR
70 INTEGER, INTENT(INOUT) :: NFUNC
71 INTEGER, INTENT(INOUT) :: IMATVIS
72 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
73 my_real, DIMENSION(NPROPM) ,INTENT(INOUT) :: pm
74 CHARACTER(LEN=NCHARTITLE) ,INTENT(IN) :: TITR
75 INTEGER, DIMENSION(MAXFUNC) ,INTENT(INOUT) :: IFUNC
76 my_real, DIMENSION(100),INTENT(INOUT) :: parmat
77 TYPE(submodel_data), DIMENSION(*),INTENT(IN) :: LSUBMODEL
78 TYPE(matparam_struct_) ,INTENT(INOUT) :: MATPARAM
79C-----------------------------------------------
80C L o c a l V a r i a b l e s
81C-----------------------------------------------
82 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
83 INTEGER :: I,IFORM,IFUNC_BLK,NPRONY,NUPARAM,NIPARAM,NORDER,ILAW
84 my_real :: RHO0,RHOR,BULK,SMAX,GS,MU0,NU,GVMAX,C1,ZEP495,FSCAL,FSCAL_UNIT
85 my_real, DIMENSION(10) :: mu,al
86 my_real, DIMENSION(100) :: gi,taux
87C=======================================================================
88 is_encrypted = .false.
89 is_available = .false.
90 ilaw = 42
91c------------------------------------------
92 CALL hm_option_is_encrypted(is_encrypted)
93c------------------------------------------
94c
95!card1 - 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!card2 - Poisson's ratior - Bulk definition - Flags
99 CALL hm_get_floatv('MAT_NU' ,nu ,is_available, lsubmodel, unitab)
100 CALL hm_get_floatv('MAT_SIGCUT' ,smax ,is_available, lsubmodel, unitab)
101 CALL hm_get_intv ('FUN_BULK' ,ifunc_blk,is_available,lsubmodel)
102 CALL hm_get_floatv('MAT_FScale' ,fscal ,is_available, lsubmodel, unitab)
103 CALL hm_get_intv ('ORDER' ,nprony ,is_available,lsubmodel)
104 CALL hm_get_intv ('IFORM' ,iform ,is_available,lsubmodel)
105!card3 - Shear hyperelastic modulus parameters
106 CALL hm_get_floatv('MAT_MUE1' ,mu(1) ,is_available, lsubmodel, unitab)
107 CALL hm_get_floatv('MAT_MUE2' ,mu(2) ,is_available, lsubmodel, unitab)
108 CALL hm_get_floatv('MAT_MUE3' ,mu(3) ,is_available, lsubmodel, unitab)
109 CALL hm_get_floatv('MAT_MUE4' ,mu(4) ,is_available, lsubmodel, unitab)
110 CALL hm_get_floatv('MAT_MUE5' ,mu(5) ,is_available, lsubmodel, unitab)
111!card4 - Shear hyperelastic modulus parameters
112 CALL hm_get_floatv('MAT_MUE6' ,mu(6) ,is_available, lsubmodel, unitab)
113 CALL hm_get_floatv('MAT_MUE7' ,mu(7) ,is_available, lsubmodel, unitab)
114 CALL hm_get_floatv('MAT_MUE8' ,mu(8) ,is_available, lsubmodel, unitab)
115 CALL hm_get_floatv('MAT_MUE9' ,mu(9) ,is_available, lsubmodel, unitab)
116 CALL hm_get_floatv('MAT_MUE10' ,mu(10) ,is_available, lsubmodel, unitab)
117!card5 - Material exponents
118 CALL hm_get_floatv('MAT_ALPHA11',al(1) ,is_available, lsubmodel, unitab)
119 CALL hm_get_floatv('MAT_ALPHA22',al(2) ,is_available, lsubmodel, unitab)
120 CALL hm_get_floatv('MAT_ALPHA33',al(3) ,is_available, lsubmodel, unitab)
121 CALL hm_get_floatv('MAT_ALPHA44',al(4) ,is_available, lsubmodel, unitab)
122 CALL hm_get_floatv('MAT_ALPHA55',al(5) ,is_available, lsubmodel, unitab)
123!card6 - Material exponents
124 CALL hm_get_floatv('MAT_ALPHA6 ',al(6) ,is_available, lsubmodel, unitab)
125 CALL hm_get_floatv('MAT_ALPHA7 ',al(7) ,is_available, lsubmodel, unitab)
126 CALL hm_get_floatv('MAT_ALPHA8 ',al(8) ,is_available, lsubmodel, unitab)
127 CALL hm_get_floatv('MAT_ALPHA9 ',al(9) ,is_available, lsubmodel, unitab)
128 CALL hm_get_floatv('MAT_ALPHA10',al(10) ,is_available, lsubmodel, unitab)
129!card7 - Prony series
130 IF (nprony > 0) THEN
131 DO i=1,nprony
132 CALL hm_get_float_array_index('gamma_arr',GI(I) ,I,IS_AVAILABLE,LSUBMODEL,UNITAB)
133 ENDDO
134 DO I=1,NPRONY
135 CALL HM_GET_FLOAT_ARRAY_INDEX('tau_arr' ,TAUX(I),I,IS_AVAILABLE,LSUBMODEL,UNITAB)
136 ENDDO
137 ENDIF
138c
139c-----------------------------
140c Check and default values
141c-----------------------------
142c
143 ! Reference density
144 IF (RHOR == ZERO) RHOR = RHO0
145c
146 ! Poisson's ratio
147 zep495 = zep4 + nine*em02 + five*em03
148 IF (nu == zero ) nu = zep495
149c
150 ! Tensile stress limit
151 IF (smax <= zero) smax = ep20
152c
153 ! Scale factor for bulk tabulation
154 CALL hm_get_floatv_dim('MAT_FScale',fscal_unit ,is_available, lsubmodel, unitab)
155 IF (fscal == zero) fscal = one*fscal_unit
156c
157 ! Check strain energy density formulation flag
158 IF (iform == 0) iform = 1
159 iform = min(2,iform)
160 iform = max(1,iform)
161c
162 ! Viscosity parameters
163 imatvis = 1
164 IF (nprony > 0 ) THEN
165 imatvis = 3
166 ENDIF
167 ! Ground state viscous shear modulus
168 gvmax = zero
169 DO i=1,nprony
170 gvmax = gvmax + gi(i)
171 ENDDO
172!------------------------------------------------------------
173! Check Ogden parameter order (skip zero values of mu)
174!------------------------------------------------------------
175 norder = 0
176 DO i = 1,10
177 IF (mu(i) /= zero) THEN
178 norder = norder + 1
179 mu(norder) = mu(i)
180 al(norder) = al(i)
181 END IF
182 END DO
183!
184 gs = zero
185 DO i = 1,norder
186 gs = gs + mu(i)*al(i)
187 END DO
188 IF (gs <= zero) THEN
189 CALL ancmsg(msgid=828, msgtype=msgerror, anmode=aninfo,
190 . i1=id,
191 . c1=titr)
192 END IF
193c
194 ! Initial shear and bulk modulus
195 mu0 = gs/two
196 bulk = gs*(one+nu)/max(em20,three*(one-two*nu))
197c--------------------------
198c Filling buffer tables
199c--------------------------
200 nfunc = 1 ! Number of functions
201 nuvar = 12 + 6*nprony ! Number of user variables
202 niparam = 3
203 nuparam = 24 + 2*nprony ! Number of material parameters
204!
205 matparam%NIPARAM = niparam
206 matparam%NUPARAM = nuparam
207 ALLOCATE (matparam%IPARAM(niparam))
208 ALLOCATE (matparam%UPARAM(nuparam))
209c
210 ! Material parameters
211 matparam%IPARAM(1) = norder
212 matparam%IPARAM(2) = nprony
213 matparam%IPARAM(3) = iform
214
215 matparam%UPARAM(1) = mu(1)
216 matparam%UPARAM(2) = mu(2)
217 matparam%UPARAM(3) = mu(3)
218 matparam%UPARAM(4) = mu(4)
219 matparam%UPARAM(5) = mu(5)
220 matparam%UPARAM(6) = mu(6)
221 matparam%UPARAM(7) = mu(7)
222 matparam%UPARAM(8) = mu(8)
223 matparam%UPARAM(9) = mu(9)
224 matparam%UPARAM(10) = mu(10)
225 matparam%UPARAM(11) = al(1)
226 matparam%UPARAM(12) = al(2)
227 matparam%UPARAM(13) = al(3)
228 matparam%UPARAM(14) = al(4)
229 matparam%UPARAM(15) = al(5)
230 matparam%UPARAM(16) = al(6)
231 matparam%UPARAM(17) = al(7)
232 matparam%UPARAM(18) = al(8)
233 matparam%UPARAM(19) = al(9)
234 matparam%UPARAM(20) = al(10)
235 matparam%UPARAM(21) = bulk
236 matparam%UPARAM(22) = nu
237 matparam%UPARAM(23) = smax
238 matparam%UPARAM(24) = fscal
239 IF (nprony > 0) THEN
240 DO i=1,nprony
241 matparam%UPARAM(24 + i) = gi(i)
242 matparam%UPARAM(24 + nprony + i) = taux(i)
243 ENDDO
244 ENDIF
245c
246 ! Function IDs
247 ifunc(1) = ifunc_blk
248c
249 ! PARMAT table
250 parmat(1) = gs
251 parmat(2) = gs*(one+nu)
252 parmat(3) = nu
253 parmat(6) = bulk
254 c1 = third*gs*(one + nu)/(one - two*nu)
255 parmat(16) = 2
256 parmat(17) = gs/(c1 + two_third*gs)
257c
258 ! PM /matparam common parameters
259
260 matparam%RHO = rhor
261 matparam%RHO0 = rho0
262 pm(100)= bulk
263c
264 ! MATPARAM keywords
265 CALL init_mat_keyword(matparam,"INCOMPRESSIBLE")
266 CALL init_mat_keyword(matparam,"TOTAL")
267 CALL init_mat_keyword(matparam,"HOOK")
268 ! Properties compatibility
269 CALL init_mat_keyword(matparam,"SHELL_ISOTROPIC")
270 CALL init_mat_keyword(matparam,"SOLID_ISOTROPIC")
271 CALL init_mat_keyword(matparam,"SPH")
272c
273c--------------------------
274c Parameters printout
275c--------------------------
276 WRITE(iout,1100) trim(titr),id,42
277 WRITE(iout,1000)
278 IF (is_encrypted) THEN
279 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA'
280 ELSE
281 WRITE(iout,1200) rho0
282 WRITE(iout,1300) nu,smax,ifunc_blk,fscal,iform
283 WRITE(iout,1400) norder
284 DO i = 1,norder
285 WRITE(iout,1500) i,mu(i),al(i)
286 ENDDO
287 WRITE(iout,1600) mu0,bulk
288 IF (nprony > 0) THEN
289 WRITE(iout,1700) nprony
290 DO i = 1, nprony
291 WRITE(iout,1800) i,gi(i),taux(i)
292 ENDDO
293 ENDIF
294 ENDIF
295C-----------------
296 RETURN
297C-----------------
298 1000 FORMAT(
299 & 5x,'------------------------------------------',/
300 & 5x,' MATERIAL MODEL : GREEN-ELASTIC (OGDEN) ',/,
301 & 5x,'------------------------------------------',/)
302 1100 FORMAT(/
303 & 5x,a,/,
304 & 5x,'MATERIAL NUMBER . . . . . . . . . . . . . . . . .=',i10/,
305 & 5x,'MATERIAL LAW. . . . . . . . . . . . . . . . . . .=',i10/)
306 1200 FORMAT(
307 & 5x,'INITIAL DENSITY . . . . . . . . . . . . . . . . .=',1pg20.13/)
308 1300 FORMAT(
309 & 5x,'POISSON RATIO NU. . . . . . . . . . . . . . . . .=',1pg20.13/,
310 & 5x,'CUT-OFF STRESS IN TENSION SIGMA CUT . . . . . . .=',1pg20.13/,
311 & 5x,'BULK FUNCTION NUMBER ID . . . . . . . . . . . . .=',i10/,
312 & 5x,'SCALE FACTOR FOR BULK FUNCTION. . . . . . . . . .=',1pg20.13/,
313 & 5x,'INCOMP. FORMULATION FLAG IFORM. . . . . . . . . .=',i2/,
314 & 5x,' IFORM = 1: STANDARD STRAIN ENERGY DENSITY (DEFAULT)',/,
315 & 5x,' IFORM = 2: MODIFIED STRAIN ENERGY DENSITY ',/)
316 1400 FORMAT(
317 & 5x,'OGDEN MODEL PARAMETERS: ',/,
318 & 5x,'------------------------ ',/,
319 & 5x,' ',/,
320 & 5x,'NUMBER OF TERMS IN OGDEN MODEL. . . . . . . . . .=',i8/)
321 1500 FORMAT(
322 & 5x,'OGDEN TERM NUMBER # ' ,i8/,
323 & 5x,'ground shear modulus. . . . . . . . . . . . . . .=',1PG20.13/,
324 & 5X,'exponent alpha. . . . . . . . . . . . . . . . . .=',1PG20.13/)
325 1600 FORMAT(
326 & 5X,'initial shear modulus . . . . . . . . . . . . . .=',1PG20.13/,
327 & 5X,'bulk modulus . . . . . . . . . . . . . . . . . .=',1PG20.13/)
328 1700 FORMAT(
329 & 5X,'prony series parameters: ',/,
330 & 5X,'------------------------ ',/,
331 & 5X,' ',/,
332 & 5X,'number of terms in prony series m . . . . . . . .=',I8/)
333 1800 FORMAT(
334 & 5X,'prony term number # ' ,I8/,
335 & 5X,'shear stiffness . . . . . . . . . . . . . . . . .=',1pg20.13/,
336 & 5x,'RELAXATION TIME . . . . . . . . . . . . . . . . .=',1pg20.13/)
337C-----------------
338 END
my_real
#define my_real
Definition cppsort.cpp:32
alpha
#define alpha
Definition eval.h:35
hm_get_float_array_index
subroutine hm_get_float_array_index(name, rval, index, is_available, lsubmodel, unitab)
Definition hm_get_float_array_index.F:105
hm_get_floatv
subroutine hm_get_floatv(name, rval, is_available, lsubmodel, unitab)
Definition hm_get_floatv.F:415
hm_get_floatv_dim
subroutine hm_get_floatv_dim(name, dim_fac, is_available, lsubmodel, unitab)
Definition hm_get_floatv_dim.F:159
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_mat42
subroutine hm_read_mat42(matparam, nuvar, maxfunc, nfunc, ifunc, parmat, imatvis, unitab, lsubmodel, id, titr, pm)
Definition hm_read_mat42.F:43
init_mat_keyword
subroutine init_mat_keyword(matparam, keyword)
Definition init_mat_keyword.F:150
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
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
submodel_mod::submodel_data
Definition submodel_mod.F:723