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

Go to the source code of this file.

Functions/Subroutines

subroutine hm_read_mat102 (uparam, maxuparam, nuparam, israte, nuvar, nfunc, parmat, unitab, mat_id, titr, mtag, lsubmodel, pm, ipm, matparam)

Function/Subroutine Documentation

◆ hm_read_mat102()

subroutine hm_read_mat102 ( real(kind=wp), dimension(maxuparam), intent(inout) uparam,
integer, intent(inout) maxuparam,
integer, intent(inout) nuparam,
integer, intent(inout) israte,
integer, intent(inout) nuvar,
integer, intent(inout) nfunc,
real(kind=wp), dimension(100), 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(nsubmod), intent(in) lsubmodel,
real(kind=wp), dimension(npropm), intent(inout) pm,
integer, dimension(npropmi), intent(inout) ipm,
type(matparam_struct_), intent(inout) matparam )

Definition at line 38 of file hm_read_mat102.F.

42C-----------------------------------------------
43C D e s c r i p t i o n
44C-----------------------------------------------
45C READ MAT LAW102 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 , ONLY : submodel_data, nsubmod
66 USE matparam_def_mod
67 USE names_and_titles_mod , ONLY : nchartitle, ncharkey, ncharfield
68 USE constant_mod , ONLY : pi, hundred80, four, zero, infinity, nine, one, six, three, two
69 USE precision_mod , ONLY : wp
70
71C-----------------------------------------------
72C I m p l i c i t T y p e s
73C-----------------------------------------------
74 IMPLICIT NONE
75C-----------------------------------------------
76C C o m m o n B l o c k s
77C-----------------------------------------------
78#include "units_c.inc"
79#include "param_c.inc"
80C-----------------------------------------------
81C D u m m y A r g u m e n t s
82C-----------------------------------------------
83 TYPE (UNIT_TYPE_),INTENT(IN) :: UNITAB
84 REAL(KIND=wp), DIMENSION(NPROPM),INTENT(INOUT) :: pm
85 REAL(KIND=wp), DIMENSION(100),INTENT(INOUT) :: parmat
86 REAL(KIND=wp), DIMENSION(MAXUPARAM) ,INTENT(INOUT) :: uparam
87 INTEGER, DIMENSION(NPROPMI),INTENT(INOUT) :: IPM
88 INTEGER, INTENT(INOUT) :: ISRATE,NFUNC,MAXUPARAM,NUPARAM,NUVAR
89 TYPE(MLAW_TAG_),INTENT(INOUT) :: MTAG
90 INTEGER,INTENT(IN) :: MAT_ID
91 CHARACTER(LEN=NCHARTITLE) ,INTENT(IN) :: TITR
92 TYPE(SUBMODEL_DATA),INTENT(IN) :: LSUBMODEL(NSUBMOD)
93 TYPE(MATPARAM_STRUCT_) ,INTENT(INOUT) :: MATPARAM
94C-----------------------------------------------
95C L o c a l V a r i a b l e s
96C-----------------------------------------------
97 REAL(KIND=wp) :: e,nu,c,pstar,amax,g, delta,stifint,pmin,phi_deg
98 REAL(KIND=wp) :: a0,a1,a2,rho0,rhor
99 INTEGER :: IFORM
100 REAL(KIND=wp) :: phi,k,alpha
101 LOGICAL :: IS_ENCRYPTED,IS_AVAILABLE
102 CHARACTER*64 :: CHAIN
103C-----------------------------------------------
104C S o u r c e L i n e s
105C-----------------------------------------------
106 is_encrypted = .false.
107 is_available = .false.
108 k = 0.0d0
109
110 CALL hm_option_is_encrypted(is_encrypted)
111
112 CALL hm_get_intv ('IFORM',iform ,is_available, lsubmodel)
113
114 CALL hm_get_floatv('MAT_RHO' ,rho0 ,is_available, lsubmodel, unitab)
115
116 CALL hm_get_floatv('MAT_E' ,e ,is_available, lsubmodel, unitab)
117 CALL hm_get_floatv('MAT_NU' ,nu ,is_available, lsubmodel, unitab)
118
119 CALL hm_get_floatv('MAT102_C' ,c ,is_available, lsubmodel, unitab)
120 CALL hm_get_floatv('MAT102_PHI' ,phi ,is_available, lsubmodel, unitab)
121 CALL hm_get_floatv('MAT102_AMAX' ,amax ,is_available, lsubmodel, unitab)
122
123 CALL hm_get_floatv('MAT102_PMIN' ,pmin ,is_available, lsubmodel, unitab)
124
125 ! C : Mohr-Coulomb cohesion (intercept) (in Pa)
126 ! PHI : Mohr-Coulomb angle of internal friction (angle, degree)
127 ! Iform : Formulation Flag (circumscribed,inscribed,middle)
128
129C-----------------------------------------------
130C UNITS
131C-----------------------------------------------
132 phi_deg = phi
133 phi = phi*pi/hundred80
134C-----------------------------------------------
135C DEFAULTS
136C-----------------------------------------------
137 IF(iform<=0 .OR. iform>=4)iform=2
138C-----------------------------------------------
139C YIELD PARAMETERS
140C-----------------------------------------------
141 g=e/two/(one+nu)
142 SELECT CASE(iform)
143 CASE(1)
144 k = six*c*cos(phi)/sqrt(three)/(three-sin(phi))
145 alpha = two*sin(phi)/sqrt(three)/(three-sin(phi))
146 CASE(2)
147 k = six*c*cos(phi)/sqrt(three)/(three+sin(phi))
148 alpha = two*sin(phi)/sqrt(three)/(three+sin(phi))
149 CASE(3)
150 k = three*c*cos(phi)/sqrt(nine+three*sin(phi)*sin(phi))
151 alpha = sin(phi)/sqrt(nine+three*sin(phi)*sin(phi))
152 END SELECT
153 a0 = k*k
154 a1 = six*k*alpha
155 a2 = nine*alpha*alpha
156
157 IF(e<=zero)THEN
158 chain='YOUNG MODULUS MUST BE DEFINED '
159 CALL ancmsg(msgid=829, msgtype=msgerror, anmode=aninfo, i1=10, i2=mat_id, c1='ERROR', c2=titr, c3=chain)
160 ENDIF
161
162 IF(nu<=zero)THEN
163 chain='POISSON RATIO MUST BE DEFINED '
164 CALL ancmsg(msgid=829, msgtype=msgerror, anmode=aninfo, i1=10, i2=mat_id, c1='ERROR', c2=titr, c3=chain)
165 ENDIF
166
167 pstar = -infinity
168 IF(a2==zero .AND. a1/=zero)THEN
169 pstar=-a0/a1
170 ELSEIF(a2/=zero)THEN
171 delta = a1*a1-four*a0*a2
172 !If intersection with the axis
173 IF(delta >= zero)THEN
174 delta=sqrt(delta)
175 pstar = (-a1+delta)/two/a2
176
177 ELSE
178 pstar = -a1/two/a2 ! extremum
179 ! chain='FITTED YIELD FUNCTION HAS NO ROOT. CHECK INPUT PARAMETER '
180 ! CALL ANCMSG( MSGID=829, MSGTYPE=MSGWARNING, ANMODE=ANINFO, I1=10, I2=MAT_ID, C1='WARNING', C2=TITR, C3=chain)
181 ENDIF
182 ELSE
183 !do nothing let user do what he wants
184 pstar = -infinity
185 ENDIF
186C
187 IF(amax==zero) amax = infinity
188 IF(pmin==zero) pmin =-infinity
189
190 rhor=rho0
191 pm(1) = rhor
192 pm(89)= rho0
193 pm(37)= pmin
194
195 israte=0
196
197C-----------------------------------------------
198C USER MATERIAL PARAMETERS DEFINITION
199C-----------------------------------------------
200 nuparam = 11
201 uparam(1) = c
202 uparam(2) = phi
203 uparam(3) = pstar
204 uparam(4) = a0
205 uparam(5) = a1
206 uparam(6) = a2
207 uparam(7) = amax
208 uparam(8) = g
209 uparam(9) = iform
210 uparam(10)= e
211 uparam(11)= nu
212 nuvar = 0
213 nfunc = 0
214 stifint = e
215 parmat(1) = stifint/three !
216 parmat(2) = stifint
217 parmat(3) = nu
218c---------------------------
219 mtag%G_PLA = 1
220 mtag%L_PLA = 1
221!---------------------------------------------------------
222 !MATPARAM keywords
223 matparam%IEOS = 18 ! Linear EOS is used by default
224 ipm(4) = 18
225 pm(32) = e / three/(one - two*nu) ! Bulk modulus for default EOS
226!---------------------------------------------------------
227
228 CALL init_mat_keyword(matparam,"ELASTO_PLASTIC")
229
230 ! Material compatibility with /EOS option
231 CALL init_mat_keyword(matparam,"EOS")
232
233 ! EOS/Thermo keyword for pressure treatment in elements
234 CALL init_mat_keyword(matparam,"HYDRO_EOS")
235
236 ! Properties compatibility
237 CALL init_mat_keyword(matparam,"SOLID_ISOTROPIC")
238 CALL init_mat_keyword(matparam,"SPH")
239C-------------------------------------------------
240C LISTING OUTPUT
241C-------------------------------------------------
242 WRITE(iout,1001) trim(titr),mat_id,102
243 WRITE(iout,1000)
244 IF(is_encrypted)THEN
245 WRITE(iout,'(5X,A,//)')'CONFIDENTIAL DATA'
246 ELSE
247 WRITE(iout,1002)rho0,rhor
248 WRITE(iout,1100)e,nu,c,phi_deg,pmin
249 WRITE(iout,1200)iform
250 SELECT CASE(iform)
251 CASE(1)
252 WRITE(iout,1201)
253 CASE(2)
254 WRITE(iout,1202)
255 CASE(3)
256 WRITE(iout,1203)
257 CASE(4)
258 WRITE(iout,1204)
259 END SELECT
260 WRITE(iout,1300)a0,a1,a2,amax,pstar
261 ENDIF
262C
263 1000 FORMAT(
264 & 5x,' EXTENDED DRUCKER-PRAGER MATERIAL (DPRAG2) ',/,
265 & 5x,' ----------------------------------------- ')
266 1001 FORMAT(/
267 & 5x,a,/,
268 & 5x,'MATERIAL NUMBER . . . . . . . . . . . .=',i10/,
269 & 5x,'MATERIAL LAW. . . . . . . . . . . . . .=',i10/)
270 1002 FORMAT(
271 & 5x,'INITIAL DENSITY . . . . . . . . . . . .=',1pg20.13/,
272 & 5x,'REFERENCE DENSITY . . . . . . . . . . .=',1pg20.13/)
273 1100 FORMAT(
274 & 5x,'YOUNG MODULUS . . . . . . . . . . . . .=',1pg20.13/
275 & 5x,'POISSON RATIO . . . . . . . . . . . . .=',1pg20.13/
276 & 5x,'COHESION. . . . . . . . . . . . . . . .=',1pg20.13/
277 & 5x,'ANGLE OF INTERNAL FRICTION. . . . . . .=',1pg20.13/
278 & 5x,'MINIMUM PRESSURE. . . . . . . . . . . .=',1pg20.13)
279 1200 FORMAT(
280 & 5x,'DRUCKER-PRAGER MATERIAL CRITERION DEFINED FROM MOHR-COULOMB PARAMETERS',/,
281 & 5x,'FORMULATION FLAG. . . . . . . . . . . .=',i10)
282 1201 FORMAT(
283 & 5x,'-> CIRCUMBSCRIBED CRITERIA')
284 1202 FORMAT(
285 & 5x,'-> MIDDLE CRITERIA')
286 1203 FORMAT(
287 & 5x,'-> INSCRIBED CRITERIA')
288 1204 FORMAT(
289 & 5x,'-> ORIGINAL MOHR-COULOMB CRITERIA')
290 1300 FORMAT(
291 & 5x,'PARAMETERS USED TO DEFINE CRITERIA',/,
292 & 5x,'A0. . . . . . . . . . . . . . . . . . .=',1pg20.13/
293 & 5x,'A1. . . . . . . . . . . . . . . . . . .=',1pg20.13/
294 & 5x,'A2. . . . . . . . . . . . . . . . . . .=',1pg20.13/
295 & 5x,'AMAX. . . . . . . . . . . . . . . . . .=',1pg20.13/
296 & 5x,'YIELD FUNCTION PRESSURE ROOT. . . . . .=',1pg20.13//)
297
298 RETURN
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:418
hm_get_intv
subroutine hm_get_intv(name, ival, is_available, lsubmodel)
Definition hm_get_intv.F:528
hm_option_is_encrypted
subroutine hm_option_is_encrypted(is_encrypted)
Definition hm_option_is_encrypted.F:288
init_mat_keyword
subroutine init_mat_keyword(matparam, keyword)
Definition init_mat_keyword.F:153
constant_mod
Definition constant_mod.F:355
elbuftag_mod
Definition elbuftag_mod.F:230
message_mod
Definition message_mod.F:1257
names_and_titles_mod
Definition names_and_titles_mod.F:1002
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1008
names_and_titles_mod::ncharkey
integer, parameter ncharkey
Definition names_and_titles_mod.F:1005
names_and_titles_mod::ncharfield
integer, parameter ncharfield
Definition names_and_titles_mod.F:1010
submodel_mod
Definition submodel_mod.F:717
submodel_mod::nsubmod
integer nsubmod
Definition submodel_mod.F:725
unitab_mod
Definition unitab_mod.F:517
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:895
submodel_mod::submodel_data
Definition submodel_mod.F:728