OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
hm_read_mat28.F File Reference
#include "implicit_f.inc"
#include "units_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine hm_read_mat28 (uparam, maxuparam, nuparam, israte, imatvis, nuvar, ifunc, maxfunc, nfunc, parmat, unitab, mat_id, titr, mtag, lsubmodel, pm, matparam)

Function/Subroutine Documentation

◆ hm_read_mat28()

subroutine hm_read_mat28 ( intent(inout) uparam,
integer, intent(inout) maxuparam,
integer, intent(inout) nuparam,
integer, intent(inout) israte,
integer, intent(inout) imatvis,
integer, intent(inout) nuvar,
integer, dimension(maxfunc), intent(inout) ifunc,
integer, intent(inout) maxfunc,
integer, intent(inout) nfunc,
intent(inout) parmat,
type (unit_type_), intent(in) unitab,
integer, intent(in) mat_id,
character(len=nchartitle), intent(in) titr,
type(mlaw_tag_), intent(inout) mtag,
type(submodel_data), dimension(*), intent(in) lsubmodel,
intent(inout) pm,
type(matparam_struct_), intent(inout) matparam )

Definition at line 38 of file hm_read_mat28.F.

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
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: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
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