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

Go to the source code of this file.

Functions/Subroutines

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

Function/Subroutine Documentation

◆ hm_read_mat83()

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

Definition at line 39 of file hm_read_mat83.F.

43C-----------------------------------------------
44C ROUTINE DESCRIPTION :
45C ===================
46C READ MAT LAW83 WITH HM READER
47C-----------------------------------------------
48C DUMMY ARGUMENTS DESCRIPTION:
49C ===================
50C
51C NAME DESCRIPTION
52C
53C IPM MATERIAL ARRAY(INTEGER)
54C PM MATERIAL ARRAY(REAL)
55C UNITAB UNITS ARRAY
56C ID MATERIAL ID(INTEGER)
57C TITR MATERIAL TITLE
58C LSUBMODEL SUBMODEL STRUCTURE
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 "param_c.inc"
76#include "units_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 INTEGER, INTENT(INOUT) :: IFUNC(MAXFUNC),NFUNC,MAXFUNC,MAXUPARAM,NUPARAM,NUVAR,ISRATE
82 INTEGER, INTENT(IN) :: MAT_ID
83 my_real, INTENT(INOUT) :: pm(npropm),parmat(100),uparam(maxuparam)
84 CHARACTER(LEN=NCHARTITLE) ,INTENT(IN) :: TITR
85 TYPE(SUBMODEL_DATA),INTENT(IN) :: LSUBMODEL(NSUBMOD)
86 TYPE(MLAW_TAG_), INTENT(INOUT) :: MTAG
87 TYPE(MATPARAM_STRUCT_) ,INTENT(INOUT) :: MATPARAM
88C-----------------------------------------------
89C L o c a l V a r i a b l e s
90C-----------------------------------------------
91 INTEGER :: I,J,NRATE,IFILTR,IDYIELD,IFUNN,IFUNT,
92 . RHOFLAG,ICOMP,IPLAS,VP,ILAW
93 my_real :: youngt,youngc,g,young,fcut,xfac,yfac,xscale,rn,rs,beta,alpha,
94 . a1,a2,aa,e0,emax,epsmax,rho0,rhor,yfac_unit,
95 . xfac_unit,xscale_unit,rn_unit,rs_unit
96 LOGICAL :: IS_AVAILABLE,IS_ENCRYPTED
97!=======================================================================
98 ilaw = 83
99!
100 is_encrypted = .false.
101 is_available = .false.
102!
103 CALL hm_option_is_encrypted(is_encrypted)
104card1
105 CALL hm_get_floatv('MAT_RHO' ,rho0 ,is_available, lsubmodel, unitab)
106 CALL hm_get_floatv('Refer_Rho' ,rhor ,is_available, lsubmodel, unitab)
107card2
108 CALL hm_get_floatv('MAT_E' ,youngt ,is_available, lsubmodel, unitab)
109 CALL hm_get_floatv('MAT_G' ,g ,is_available, lsubmodel, unitab)
110 CALL hm_get_intv ('MAT_IMASS' ,rhoflag ,is_available, lsubmodel)
111 CALL hm_get_intv ('COMP_OPT' ,icomp ,is_available, lsubmodel)
112 CALL hm_get_floatv('MAT_ECOMP' ,youngc ,is_available, lsubmodel, unitab)
113card3
114 CALL hm_get_intv ('FUN_A1' ,idyield ,is_available ,lsubmodel)
115 CALL hm_get_floatv('FScale11' ,yfac ,is_available, lsubmodel, unitab)
116 CALL hm_get_floatv('FScale22' ,xfac ,is_available, lsubmodel, unitab)
117 CALL hm_get_floatv('MAT_ALPHA' ,alpha ,is_available, lsubmodel, unitab)
118 CALL hm_get_floatv('MAT_Beta' ,beta ,is_available, lsubmodel, unitab)
119card4
120 CALL hm_get_floatv('MAT_R00' ,rn ,is_available, lsubmodel, unitab)
121 CALL hm_get_floatv('MAT_R45' ,rs ,is_available, lsubmodel, unitab)
122 CALL hm_get_intv ('Fsmooth' ,israte ,is_available ,lsubmodel)
123 CALL hm_get_floatv('Fcut' ,fcut ,is_available, lsubmodel, unitab)
124card5
125 CALL hm_get_intv ('FUN_A2' ,ifunn ,is_available ,lsubmodel)
126 CALL hm_get_intv ('FUN_A3' ,ifunt ,is_available ,lsubmodel)
127 CALL hm_get_floatv('FScale33' ,xscale ,is_available, lsubmodel, unitab)
128c-------------------------------------------------------------------------------------
129 CALL hm_get_floatv_dim('FScale11' ,yfac_unit ,is_available, lsubmodel, unitab)
130 CALL hm_get_floatv_dim('FScale22' ,xfac_unit ,is_available, lsubmodel, unitab)
131 CALL hm_get_floatv_dim('FScale33' ,xscale_unit ,is_available, lsubmodel, unitab)
132 CALL hm_get_floatv_dim('MAT_R00' ,rn_unit ,is_available, lsubmodel, unitab)
133 CALL hm_get_floatv_dim('MAT_R45' ,rs_unit ,is_available, lsubmodel, unitab)
134c-------------------------------------------------------------------------------------
135 IF (rhor == zero) rhor=rho0
136 pm(1) = rhor
137 pm(89)= rho0
138!
139 IF (xfac == zero) xfac = one*xfac_unit
140 IF (beta == zero) beta = two
141 IF (rn == zero) rn = one*rn_unit
142 IF (rs == zero) rs = one*rs_unit
143 iplas = 0 ! hidden
144 IF (iplas == 0) iplas = 2
145 IF (israte == 0) THEN
146 vp = 0
147 ELSE
148 vp = 1
149 israte = 1
150 END IF
151 IF (fcut == zero) fcut = 10000.0d0*unitab%FAC_T_WORK
152C---
153 nfunc = 3
154 ifunc(1) = ifunn
155 ifunc(2) = ifunt
156 ifunc(3) = idyield
157C---
158 IF (ifunc(3) == 0) THEN
159 CALL ancmsg(msgid=126, msgtype=msgerror, anmode=aninfo_blind_1,
160 . i1=mat_id,
161 . c1=titr,
162 . i2=ifunc(3))
163 ENDIF
164
165 IF (yfac == zero) yfac = one*yfac_unit
166 IF (xscale == zero) xscale = one*xscale_unit
167 IF( g == zero) g = youngt
168 IF (youngc == zero) youngc = youngt
169 IF (icomp /= 1) icomp = 0
170 young = max(youngt,youngc)
171 alpha = min(alpha, one)
172C----------------
173 uparam(1) = youngt
174 uparam(2) = alpha
175 uparam(3) = beta
176 uparam(4) = yfac
177 uparam(5) = xscale
178 uparam(6) = rn
179 uparam(7) = rs
180 uparam(8) = xfac
181 uparam(9) = rhoflag
182 uparam(10)= iplas
183 uparam(11) = g
184 uparam(12) = icomp
185 uparam(13) = youngc
186 uparam(14) = vp ! total or plastic strain rate flag
187C----------------
188 nuparam = 14
189 nuvar = 1
190C----------------
191 parmat(1) = young/three
192 parmat(2) = young
193 parmat(5) = fcut
194C Formulation for solid elements time step computation.
195 parmat(16) = 2
196 parmat(17) = one ! (ONE - TWO*NU)/(ONE - NU), NU=0
197!
198 mtag%G_GAMA = 9
199 mtag%G_EPSD = 1
200 mtag%G_PLA = 1 !global
201 mtag%L_EPE = 3
202 mtag%L_PLA = 1
203 mtag%L_EPSD = 1
204 mtag%L_DMG = 1
205 mtag%G_DMG = 1
206c-----------------
207 ! MATPARAM keywords
208 CALL init_mat_keyword(matparam,"HOOK")
209 ! Properties compatibility
210 CALL init_mat_keyword(matparam,"SOLID_COHESIVE")
211C----------------
212 WRITE(iout,1100) trim(titr),mat_id,83
213 WRITE(iout,1000)
214 IF (is_encrypted) THEN
215 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA'
216 ELSE
217 WRITE(iout,1200) rho0
218 WRITE(iout,1300) youngt,youngc,g,icomp,rhoflag,idyield,yfac,xfac,alpha,beta
219 WRITE(iout,1400)ifunn,ifunt,xscale,rn,rs,israte,fcut
220 IF (vp ==0) THEN
221 WRITE(iout,1500)
222 ELSE
223 WRITE(iout,1600)
224 END IF
225 ENDIF
226C-----------
227 RETURN
228C-----------
229 1000 FORMAT(
230 & 5x,' CONNECTION MATERIAL LAW 83 ',/,
231 & 5x,' -------------------------- ',/)
232 1100 FORMAT(/
233 & 5x,a,/,
234 & 5x,'MATERIAL NUMBER . . . . . . . . . . . . . . . . . .=',i10/,
235 & 5x,'MATERIAL LAW. . . . . . . . . . . . . . . . . . . .=',i10/)
236 1200 FORMAT(
237 & 5x,'INITIAL DENSITY . . . . . . . . . . . . . . . . . .=',1pg20.13)
238 1300 FORMAT(
239 & 5x,'YOUNG MODULUS PER THICKNESS UNIT IN TENSION . . . .=',1pg20.13/,
240 & 5x,'YOUNG MODULUS PER THICKNESS UNIT IN COMPRESSION . .=',1pg20.13/,
241 & 5x,'SHEAR MODULUS PER THICKNESS UNIT . . . . . . . . .=',1pg20.13/,
242 & 5x,'COMPRESSION BEHAVIOR (FLAG) . . . . . . . . . . . .=',i10/,
243 & 5x,' = 0 => ELASTO-PLASTIC '/,
244 & 5x,' = 1 => LINEAR ELASTIC '/,
245 & 5x,'DENSITY UNIT FLAG . . . . . . . . . . . . . . . . .=',i10/,
246 & 5x,' = 0 => VOLUME DENSITY '/,
247 & 5x,' = 1 => SURFACE DENSITY '/,
248 & 5x,'YIELD STRESS FUNCTION . . . . . . . . . . . . . . .=',i10/,
249 & 5x,'SCALE FACTOR FOR YIELD FUNCTION . . . . . . . . . .=',1pg20.13/,
250 & 5x,'SCALE FACTOR FOR YIELD FUNCTION ABSCISSA . . . . .=',1pg20.13/,
251 & 5x,'PARAMETER ALPHA IN YIELD FUNCTION . . . . . . . . .=',1pg20.13/,
252 & 5x,'PARAMETER BETA IN YIELD FUNCTION . . . . . . . . .=',1pg20.13)
253 1400 FORMAT(
254 & 5x,'STRAIN RATE FUNCTION IN NORMAL DIRECTION. . . . . .=',i10/,
255 & 5x,'STRAIN RATE FUNCTION IN TANGENT DIRECTION . . . . .=',i10/,
256 & 5x,'SCALE FACTOR FOR STRAIN RATE IN FUNCTION . . . . .=',1pg20.13/,
257 & 5x,'RN VARIABLE . . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
258 & 5x,'RS VARIABLE . . . . . . . . . . . . . . . . . . . .=',1pg20.13/,
259 & 5x,'STRAIN RATE FILTERING FLAG . . . . . . . . . . . .=',i10/,
260 & 5x,'CUT FREQ FOR STRAIN RATE FILTERING . . . . . . . .=',1pg20.13)
261 1500 FORMAT(5x,'USING TOTAL STRAIN RATE',/)
262 1600 FORMAT(5x,'USING PLASTIC STRAIN RATE',/)
263C--------
my_real
#define my_real
Definition cppsort.cpp:32
alpha
#define alpha
Definition eval.h:35
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
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