OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_mat28.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_mat28 ../starter/source/materials/mat/mat028/hm_read_mat28.f
25!||--- called by ------------------------------------------------------
26!|| hm_read_mat ../starter/source/materials/mat/hm_read_mat.F90
27!||--- calls -----------------------------------------------------
28!|| hm_get_floatv ../starter/source/devtools/hm_reader/hm_get_floatv.F
29!|| hm_get_floatv_dim ../starter/source/devtools/hm_reader/hm_get_floatv_dim.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_mat28(UPARAM ,MAXUPARAM,NUPARAM ,ISRATE , IMATVIS ,
39 . NUVAR ,IFUNC ,MAXFUNC ,NFUNC , PARMAT ,
40 . UNITAB ,MAT_ID ,TITR ,MTAG , LSUBMODEL,
41 . PM ,MATPARAM )
42C-----------------------------------------------
43C D e s c r i p t i o n
44C-----------------------------------------------
45C READ MAT LAW28 WITH HM READER ( TO BE COMPLETED )
46C
47C DUMMY ARGUMENTS DESCRIPTION:
48C ===================
49C
50C NAME DESCRIPTION
51C
52C IPM MATERIAL ARRAY(INTEGER)
53C PM MATERIAL ARRAY(REAL)
54C UNITAB UNITS ARRAY
55C MAT_ID MATERIAL ID(INTEGER)
56C TITR MATERIAL TITLE
57C LSUBMODEL SUBMODEL STRUCTURE
58C
59C-----------------------------------------------
60C M o d u l e s
61C-----------------------------------------------
62 USE unitab_mod
63 USE elbuftag_mod
64 USE message_mod
65 USE submodel_mod
66 USE matparam_def_mod
67 USE names_and_titles_mod , ONLY : nchartitle
68C-----------------------------------------------
69C I m p l i c i t T y p e s
70C-----------------------------------------------
71#include "implicit_f.inc"
72C-----------------------------------------------
73C C o m m o n B l o c k s
74C-----------------------------------------------
75#include "units_c.inc"
76#include "param_c.inc"
77C-----------------------------------------------
78C D u m m y A r g u m e n t s
79C-----------------------------------------------
80 TYPE (UNIT_TYPE_),INTENT(IN) ::UNITAB
81 my_real, DIMENSION(NPROPM) ,INTENT(INOUT) :: PM
82 my_real, DIMENSION(100) ,INTENT(INOUT) :: PARMAT
83 my_real, DIMENSION(MAXUPARAM) ,INTENT(INOUT) :: uparam
84 INTEGER, DIMENSION(MAXFUNC) ,INTENT(INOUT) :: IFUNC
85 INTEGER, INTENT(INOUT) :: ISRATE,IMATVIS,NFUNC,MAXFUNC,MAXUPARAM,NUPARAM,NUVAR
86 TYPE(mlaw_tag_),INTENT(INOUT) :: MTAG
87 INTEGER,INTENT(IN) :: MAT_ID
88 CHARACTER(LEN=NCHARTITLE) ,INTENT(IN) :: TITR
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 my_real
95 . e11,e22,e33,g12,g23,g31,emx11,emx22,emx33,emx12,emx23,emx31,
96 . fac1,fac2,fac3,fac4,fac5,fac6,dmin,dmax
97 INTEGER I11,I22,I33,I12,I23,I31,IF1,IF2,IMAT
98 my_real :: RHOR,RHO0,FAC_UNIT
99 LOGICAL :: IS_ENCRYPTED,IS_AVAILABLE
100C-----------------------------------------------
101C S o u r c e L i n e s
102C-----------------------------------------------
103 is_encrypted = .false.
104 is_available = .false.
105 israte=0
106 imatvis= 0
107 imat = 28
108
109
110 CALL hm_option_is_encrypted(is_encrypted)
111
112 CALL hm_get_floatv('MAT_RHO' ,rho0 ,is_available, lsubmodel, unitab)
113 CALL hm_get_floatv('Refer_Rho' ,rhor ,is_available, lsubmodel, unitab)
114
115 CALL hm_get_floatv('MAT_EA' ,e11 ,is_available, lsubmodel, unitab)
116 CALL hm_get_floatv('MAT_EB' ,e22 ,is_available, lsubmodel, unitab)
117 CALL hm_get_floatv('MAT_EC' ,e33 ,is_available, lsubmodel, unitab)
118
119 CALL hm_get_floatv('MAT_GAB' ,g12 ,is_available, lsubmodel, unitab)
120 CALL hm_get_floatv('MAT_GBC' ,g23 ,is_available, lsubmodel, unitab)
121 CALL hm_get_floatv('MAT_GCA' ,g31 ,is_available, lsubmodel, unitab)
122
123 CALL hm_get_intv ('FUN_A1' ,i11 ,is_available, lsubmodel)
124 CALL hm_get_intv ('FUN_B1' ,i22 ,is_available, lsubmodel)
125 CALL hm_get_intv ('FUN_A2' ,i33 ,is_available, lsubmodel)
126 CALL hm_get_intv ('Gflag' ,if1 ,is_available, lsubmodel)
127 CALL hm_get_floatv('FScale11' ,fac1 ,is_available, lsubmodel, unitab)
128 CALL hm_get_floatv('FScale22' ,fac2 ,is_available, lsubmodel, unitab)
129 CALL hm_get_floatv('FScale33' ,fac3 ,is_available, lsubmodel, unitab)
130
131 CALL hm_get_floatv('MAT_EPSR1' ,emx11 ,is_available, lsubmodel, unitab)
132 CALL hm_get_floatv('MAT_EPSR2' ,emx22 ,is_available, lsubmodel, unitab)
133 CALL hm_get_floatv('MAT_EPSR3' ,emx33 ,is_available, lsubmodel, unitab)
134
135 CALL hm_get_intv ('FUN_A3' ,i12 ,is_available, lsubmodel)
136 CALL hm_get_intv ('FUN_B3' ,i23 ,is_available, lsubmodel)
137 CALL hm_get_intv ('FUN_A4' ,i31 ,is_available, lsubmodel)
138 CALL hm_get_intv ('Vflag' ,if2 ,is_available, lsubmodel)
139 CALL hm_get_floatv('FScale12' ,fac4 ,is_available, lsubmodel, unitab)
140 CALL hm_get_floatv('FScale23' ,fac5 ,is_available, lsubmodel, unitab)
141 CALL hm_get_floatv('FScale13' ,fac6 ,is_available, lsubmodel, unitab)
142
143 CALL hm_get_floatv('MAT_EPSR4' ,emx12 ,is_available, lsubmodel, unitab)
144 CALL hm_get_floatv('MAT_EPSR5' ,emx23 ,is_available, lsubmodel, unitab)
145 CALL hm_get_floatv('MAT_EPSR6' ,emx31 ,is_available, lsubmodel, unitab)
146
147 CALL hm_get_floatv_dim('FScale11' ,fac_unit ,is_available, lsubmodel, unitab)
148
149
150C -------------------------
151 IF (fac1 == zero) fac1 = one * fac_unit
152 IF (fac2 == zero) fac2 = one * fac_unit
153 IF (fac3 == zero) fac3 = one * fac_unit
154 IF (fac4 == zero) fac4 = one * fac_unit
155 IF (fac5 == zero) fac5 = one * fac_unit
156 IF (fac6 == zero) fac6 = one * fac_unit
157 uparam(1)=e11
158 uparam(2)=e22
159 uparam(3)=e33
160 uparam(4)=g12
161 uparam(5)=g23
162 uparam(6)=g31
163 uparam(7)=if1
164 uparam(8)=if2
165 uparam(9) =emx11
166 uparam(10)=emx22
167 uparam(11)=emx33
168 uparam(12)=emx12
169 uparam(13)=emx23
170 uparam(14)=emx31
171 IF(uparam(9) ==zero)uparam(9) =infinity
172 IF(uparam(10)==zero)uparam(10)=infinity
173 IF(uparam(11)==zero)uparam(11)=infinity
174 IF(uparam(12)==zero)uparam(12)=infinity
175 IF(uparam(13)==zero)uparam(13)=infinity
176 IF(uparam(14)==zero)uparam(14)=infinity
177 uparam(15)=fac1
178 uparam(16)=fac2
179 uparam(17)=fac3
180 uparam(18)=fac4
181 uparam(19)=fac5
182 uparam(20)=fac6
183 nuparam = 20
184C
185 parmat(1) = max(e11,e22,e33,g12,g23,g31)
186C
187 ifunc(1)=i11
188 ifunc(2)=i22
189 ifunc(3)=i33
190 ifunc(4)=i12
191 ifunc(5)=i23
192 ifunc(6)=i31
193 nfunc = 6
194C
195 nuvar = 6
196CCC Formulation for solid elements time step computation.
197 dmin = min(e11*e22, e22*e33,e11*e33)
198 dmax = max(e11,e22,e33)
199 parmat(16) = 1
200 parmat(17) = dmin/dmax/dmax
201
202 mtag%G_PLA = 1
203 mtag%L_PLA = 1
204c
205 ! MATPARAM keywords
206 CALL init_mat_keyword(matparam,"HOOK")
207 CALL init_mat_keyword(matparam,"ORTHOTROPIC")
208c
209 ! Properties compatibility
210 CALL init_mat_keyword(matparam,"SOLID_ORTHOTROPIC")
211 CALL init_mat_keyword(matparam,"SPH")
212c
213 IF(rhor==zero)rhor=rho0
214 pm(01)=rhor
215 pm(89)=rho0
216c--------------------------------------------------
217 CALL init_mat_keyword(matparam,"COMPRESSIBLE")
218 CALL init_mat_keyword(matparam,"SMALL_STRAIN")
219c--------------------------------------------------
220 WRITE(iout,1001) trim(titr),mat_id,imat
221 WRITE(iout,1000)
222
223 IF(is_encrypted)THEN
224 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA'
225 ELSE
226 WRITE(iout,1002) rho0
227 WRITE(iout,1100)e11,e22,e33,g12,g23,g31,i11,i22,i33,i12,i23,i31,fac1,fac2,fac3,fac4,fac5,fac6
228 IF(if1+if2/=0)WRITE(iout,1200)if1,if2
229 IF(emx11+emx22+emx33+emx12+emx23+emx31/=0)WRITE(iout,1300)emx11,emx22,emx33,emx12,emx23,emx31
230 ENDIF
231C
232 RETURN
233 1000 FORMAT(
234 & 5x,40h honeycomb law ,/,
235 & 5x,40h ------------- ,//)
236 1001 FORMAT(/
237 & 5x,a,/,
238 & 5x,'MATERIAL NUMBER . . . . . . . . . . . .=',i10/,
239 & 5x,'MATERIAL LAW. . . . . . . . . . . . . .=',i10/)
240 1002 FORMAT(
241 & 5x,'INITIAL DENSITY . . . . . . . . . . . .=',1pg20.13/)
242 1100 FORMAT(
243 & 5x,'E11 . . . . . . . . . . . . . . . . . .=',1pg20.13/
244 & 5x,'E22 . . . . . . . . . . . . . . . . . .=',1pg20.13/
245 & 5x,'E33 . . . . . . . . . . . . . . . . . .=',1pg20.13/
246 & 5x,'G12 . . . . . . . . . . . . . . . . . .=',1pg20.13/
247 & 5x,'G23 . . . . . . . . . . . . . . . . . .=',1pg20.13/
248 & 5x,'G31 . . . . . . . . . . . . . . . . . .=',1pg20.13/
249 & 5x,'YIELD STRESS 11 FUNCTION NUMBER . . . .=',i10/
250 & 5x,'YIELD STRESS 22 FUNCTION NUMBER . . . .=',i10/
251 & 5x,'YIELD STRESS 33 FUNCTION NUMBER . . . .=',i10/
252 & 5x,'YIELD STRESS 12 FUNCTION NUMBER . . . .=',i10/
253 & 5x,'YIELD STRESS 23 FUNCTION NUMBER . . . .=',i10/
254 & 5x,'YIELD STRESS 31 FUNCTION NUMBER . . . .=',i10/
255 & 5x,'STRESS 11 FUNCTION SCALE FACTOR . . . .=',1pg20.13/
256 & 5x,'STRESS 22 FUNCTION SCALE FACTOR . . . .=',1pg20.13/
257 & 5x,'STRESS 33 FUNCTION SCALE FACTOR . . . .=',1pg20.13/
258 & 5x,'STRESS 12 FUNCTION SCALE FACTOR . . . .=',1pg20.13/
259 & 5x,'STRESS 23 FUNCTION SCALE FACTOR . . . .=',1pg20.13/
260 & 5x,'STRESS 31 FUNCTION SCALE FACTOR . . . .=',1pg20.13/)
261 1200 FORMAT(
262 & 5x,'YIELD FUNCTION 11,22,33 FLAG . . . . .=',i10/
263 & 5x,'YIELD FUNCTION 12,23,31 FLAG . . . . .=',i10/)
264 1300 FORMAT(
265 & 5x,'TENSION FAILURE STRAIN 11 . . . . . . .=',1pg20.13/
266 & 5x,'TENSION FAILURE STRAIN 22 . . . . . . .=',1pg20.13/
267 & 5x,'TENSION FAILURE STRAIN 33 . . . . . . .=',1pg20.13/
268 & 5x,'SHEAR FAILURE STRAIN 12 . . . . . . . .=',1pg20.13/
269 & 5x,'SHEAR FAILURE STRAIN 23 . . . . . . . .=',1pg20.13/
270 & 5x,'SHEAR FAILURE STRAIN 31 . . . . . . . .=',1pg20.13//)
271
272
273 END
274
275
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_mat28
subroutine hm_read_mat28(uparam, maxuparam, nuparam, israte, imatvis, nuvar, ifunc, maxfunc, nfunc, parmat, unitab, mat_id, titr, mtag, lsubmodel, pm, matparam)
Definition hm_read_mat28.F:42
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
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
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