OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
updmat.F File Reference
#include "implicit_f.inc"
#include "units_c.inc"
#include "param_c.inc"
#include "scr17_c.inc"
#include "scr19_c.inc"
#include "com04_c.inc"
#include "tabsiz_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine updmat (bufmat, pm, ipm, table, func_id, npc, pld, sensors, nloc_dmg, mlaw_tag, mat_param)

Function/Subroutine Documentation

◆ updmat()

subroutine updmat ( target bufmat,
pm,
integer, dimension(npropmi,nummat), intent(inout), target ipm,
type(ttable), dimension(ntable), intent(inout) table,
integer, dimension(nfunct) func_id,
integer, dimension(snpc) npc,
pld,
type (sensors_), intent(in) sensors,
type (nlocal_str_) nloc_dmg,
type(mlaw_tag_), dimension(nummat), target mlaw_tag,
type(matparam_struct_), dimension(nummat), intent(inout), target mat_param )

Definition at line 74 of file updmat.F.

77C-----------------------------------------------
78C M o d u l e s
79C-----------------------------------------------
80 USE table_mod
81 USE nlocal_reg_mod
82 USE elbuftag_mod
83 USE message_mod
84 USE matparam_def_mod
85 USE sensor_mod
86 USE law19_upd_mod , ONLY : law19_upd
87 USE law81_upd_mod , ONLY : law81_upd
88 USE law87_upd_mod , ONLY : law87_upd
89 USE law129_upd_mod , ONLY : law129_upd
90 USE law133_upd_mod , ONLY : law133_upd
91 USE law163_upd_mod , ONLY : law163_upd
92 USE law190_upd_mod , ONLY : law190_upd
93 USE visc_param_mod
94 USE names_and_titles_mod , ONLY : nchartitle
95 USE eosfun_usr2sys_mod , ONLY : eosfun_usr2sys
96C-----------------------------------------------
97C I m p l i c i t T y p e s
98C-----------------------------------------------
99#include "implicit_f.inc"
100C-----------------------------------------------
101C C o m m o n B l o c k s
102C-----------------------------------------------
103#include "units_c.inc"
104#include "param_c.inc"
105#include "scr17_c.inc"
106#include "scr19_c.inc"
107#include "com04_c.inc"
108#include "tabsiz_c.inc"
109C-----------------------------------------------
110C D u m m y A r g u m e n t s
111C-----------------------------------------------
112 INTEGER NPC(SNPC), FUNC_ID(NFUNCT)
113 INTEGER, DIMENSION(NPROPMI,NUMMAT) ,INTENT(INOUT) :: IPM
114 my_real pm(npropm,nummat), pld(stf),bufmat(sbufmat)
115 TARGET ipm,bufmat
116 TYPE(MLAW_TAG_), TARGET, DIMENSION(NUMMAT) :: MLAW_TAG
117 TYPE(TTABLE), DIMENSION(NTABLE) ,INTENT(INOUT) :: TABLE
118 TYPE (NLOCAL_STR_) :: NLOC_DMG
119 TYPE (SENSORS_) ,INTENT(IN) :: SENSORS
120 TYPE(MATPARAM_STRUCT_) ,DIMENSION(NUMMAT), INTENT(INOUT) ,TARGET :: MAT_PARAM
121C-----------------------------------------------
122C L o c a l V a r i a b l e s
123C-----------------------------------------------
124 INTEGER IMAT,ILAW,IRUP,MAT_ID,IADD,IFC,FLAG_HE,NFUNC,NFUNCL,NUMTABL,
125 . NPARAM,NUPARF,NFAIL,IR,FLAG_FIT,IFAILG,
126 . NPRONY,J,ITYPE_VISC,FORM
127 my_real :: gama_inf
128 INTEGER, DIMENSION(NUMMAT) :: MULLINS,GURSON
129 INTEGER, DIMENSION(MAXTAB) :: TABLE_ID
130 INTEGER, DIMENSION(:), POINTER :: IFUNC,ITABLE
131 my_real, DIMENSION(:), POINTER :: uparam,uparf
132 TYPE(MLAW_TAG_), POINTER :: MTAG
133 CHARACTER(LEN=NCHARTITLE)::TITR, TITLE
134 TYPE(MATPARAM_STRUCT_) , POINTER :: MATPARAM
135 TYPE (VISC_PARAM_) , POINTER :: VISC
136 INTEGER :: NFUNC_EOS
137C=======================================================================
138c TRAITEMENT OF MATERIAL PARAMETERS AND FUNCTIONS
139c-----------------------------------------------------------------------
140 mullins(1:nummat) = 0
141 gurson(1:nummat) = 0
142 DO imat=1,nummat
143 ilaw = mat_param(imat)%ILAW
144 nfail = mat_param(imat)%NFAIL
145 IF (nfail > 0) THEN
146 mtag => mlaw_tag(imat)
147 DO ir = 1,nfail
148 irup = mat_param(imat)%FAIL(ir)%IRUPT
149 IF (irup == 33 .AND. ilaw /=100 .AND. ilaw /=95 .AND. ilaw/=92) THEN
150 CALL ancmsg(msgid=1769, msgtype=msgerror, anmode=aninfo_blind_1,
151 . i1=imat,
152 . i2=ilaw)
153 ENDIF
154 IF (irup == 33 .AND. (ilaw ==100 .OR. ilaw ==95 .OR. ilaw==92)) THEN
155 mullins(imat) = ir
156 ELSEIF (irup == 35 .and. ilaw==104) THEN
157 IF (gurson(imat) == 0) THEN
158 gurson(imat) = ir
159 ELSE
160 ! Error message
161 CALL ancmsg(msgid=1650,msgtype=msgerror,
162 . anmode=aninfo_blind_1,i1=imat)
163 ENDIF
164 ELSE IF (irup == 14) THEN ! Hashin
165 mtag%G_DMG = 1
166 ELSE IF (irup == 28) THEN ! Alter
167 mtag%G_DMG = 2
168 ! -> Isotropic stress softening
169 ELSE IF (irup == 10 .OR. irup == 11 .OR. irup == 13 .OR. irup == 30 .OR.
170 . irup == 41 .OR. irup == 42 .OR. irup == 44 .OR. irup == 45 .OR.
171 . irup == 46 .OR. irup == 47 .OR. irup == 50 .OR. irup == 51) THEN
172 mtag%L_DMGSCL = 1
173 mtag%G_DMGSCL = 1
174 ! -> Orthotropic stress softening
175 ELSE IF (irup == 48) THEN
176 mtag%G_DMGSCL = 6
177 mtag%L_DMGSCL = 6
178 ENDIF
179 ENDDO
180 ENDIF
181 ENDDO
182c-----------------------------------------------------------------------
183
184
185 DO imat=1,nummat
186C
187 mat_id = mat_param(imat)%MAT_ID
188 ilaw = mat_param(imat)%ILAW
189
190 !EoS
191 nfunc_eos = mat_param(imat)%EOS%NFUNC
192 IF(nfunc_eos > 0)THEN
193 title = mat_param(imat)%EOS%title
194 CALL eosfun_usr2sys(title,mat_id,nfunc_eos,mat_param(imat)%EOS%FUNC,func_id, nfunct )
195 ENDIF
196
197c
198 IF (ilaw < 29.and.ilaw/=19.or.ilaw==32.or.ilaw==43.or.ilaw==49.or.ilaw==52.or.
199 . ilaw==57.or.ilaw==59.or.ilaw==60.or.ilaw==65.or.ilaw==75.or.
200 . ilaw==78 ) cycle
201c
202 CALL fretitl2(titr,ipm(npropmi-ltitr+1,imat),ltitr)
203c
204 iadd = max(1,ipm(7,imat))
205 nparam = min(ipm(9,imat),sbufmat)
206 uparam => bufmat(iadd:iadd+nparam-1)
207 !/visc/prony
208 itype_visc = mat_param(imat)%IVISC
209 gama_inf = one
210 IF (itype_visc == 2) THEN
211 visc => mat_param(imat)%VISC
212 nprony = visc%IPARAM(1)
213 form = visc%IPARAM(2)
214 IF(form == 1) THEN
215 gama_inf = one
216 DO j=1,nprony
217 gama_inf = gama_inf - visc%UPARAM(j)
218 ENDDO
219 gama_inf = max(em20, gama_inf)
220 ENDIF
221 ENDIF
222 gama_inf = min(gama_inf, one)
223c
224c-----------------------------------------------------------------------
225c Replace Func ID -> function system number in IFUNC table
226c-----------------------------------------------------------------------
227c
228 nfunc = ipm(10 ,imat)
229 numtabl = ipm(226,imat)
230 ifunc => ipm(10+1:10+nfunc,imat)
231 itable => ipm(226+1:226+numtabl,imat)
232c
233 IF (nfunc > 0) THEN
234 CALL matfun_usr2sys(titr,mat_id,nfunc,ifunc,func_id )
235 ENDIF
236
237 IF (ilaw /= 73 .and. ilaw /= 74 .and. ilaw /= 80) THEN
238 IF (numtabl > 0 ) THEN
239 table_id(1:numtabl) = ipm(226+1:226+numtabl,imat)
240 CALL mattab_usr2sys(titr,mat_id,ntable,table,numtabl,itable)
241 ENDIF
242 ENDIF
243
244c-----------------------------------------------------------------------
245 SELECT CASE (ilaw)
246c-------------------------
247 CASE (19)
248 CALL law19_upd(mat_param(imat),sensors)
249c-------------------------
250 CASE (36)
251 mtag => mlaw_tag(imat)
252 CALL law36_upd(iout ,titr ,mat_id ,nparam,uparam ,
253 . nfunc ,ifunc ,func_id,npc ,pld ,
254 . mtag ,nfunct )
255c-------------------------
256 CASE (42)
257 CALL law42_upd(mat_param(imat),iout,titr,mat_id,pm(1,imat),gama_inf)
258c-------------------------
259 CASE (58)
260 nfuncl = ipm(6,imat)
261 ifunc => ipm(10+1:10+nfunc+nfuncl,imat)
262 CALL law58_upd(mat_param(imat),titr,npc,pld,
263 . nfunc,nfuncl,ifunc,mat_id,func_id,
264 . pm(1,imat),sensors)
265c-------------------------
266 CASE (62)
267 CALL law62_upd(iout,titr,mat_id,nparam,uparam, pm(1,imat),gama_inf)
268c------------------
269 CASE (69)
270 CALL law69_upd(iout, titr , mat_id, uparam, nfunc,nfunct,
271 . ifunc, func_id , npc , pld ,
272 . pm(1,imat),ipm(1,imat),gama_inf)
273c
274c------------------
275 CASE (70)
276c------------------
277 CALL law70_upd(mat_param(imat) ,titr ,mat_id ,nparam ,uparam ,
278 . nfunc ,ifunc ,npc ,pld ,iout ,
279 . nfunct ,func_id ,npropm ,pm(1,imat))
280c------------------
281 CASE (76)
282 CALL law76_upd(iout ,titr ,mat_id ,nparam ,mat_param(imat) ,
283 . uparam ,numtabl ,itable ,table ,nfunc ,
284 . ifunc ,npc ,pld )
285c
286c------------------
287 CASE (77)
288 CALL law77_upd(titr ,mat_id ,nparam ,mat_param(imat) ,
289 . uparam ,nfunc ,ifunc ,npc ,pld )
290c
291c------------------
292 CASE (81)
293 CALL law81_upd(mat_param(imat) ,nfunc ,ifunc ,npc ,snpc ,
294 . pld ,stf ,pm(1,imat),npropm,iout ,
295 . mat_id ,titr )
296c
297c------------------
298 CASE (87) !BARLAT 2D
299 flag_fit = mat_param(imat)%IPARAM(3)
300 IF (flag_fit == 1) THEN
301 CALL law87_upd(iout ,titr ,mat_id ,mat_param(imat) )
302 ENDIF
303c------------------
304 CASE (90)
305c
306 CALL law90_upd(iout,titr,mat_id ,uparam,
307 . ipm(1,imat),func_id,npc,pld,pm(1,imat),
308 . nfunct)
309c------------------
310 CASE (92) ! Arruda-Boyce hyperelastic law
311c------------------
312c
313 uparam(13) = mullins(imat)
314 ifc = ipm(11,imat)
315 IF (ifc > 0) THEN
316 CALL law92_upd(iout ,titr ,mat_id ,uparam ,nfunc ,nfunct,
317 . ifunc ,func_id ,npc ,pld ,pm(1,imat),ipm(1,imat))
318 ENDIF
319c------------------
320 CASE (95) ! Bergstrom-Boyce visco-elastic law
321 uparam(21) = mullins(imat)
322c------------------
323 CASE (100) ! parallel network model
324c
325 flag_he = uparam(2)
326 uparam(3) = mullins(imat)
327 IF (flag_he == 2 )THEN
328 ifc = ipm(11,imat)
329 IF( ifc /=0 ) THEN
330 CALL law100_upd_ab(iout, titr , mat_id, uparam, nfunc,
331 . ifunc, func_id , npc , pld , pm(1,imat))
332 ENDIF
333 ELSEIF(flag_he == 13 )THEN
334 CALL law100_upd_nht(iout, titr , mat_id, uparam, nfunc,
335 . ifunc, func_id , npc , pld , pm(1,imat))
336 ENDIF
337c------------------
338 CASE (104) ! Drucker-Gurson
339c
340 ir = gurson(imat)
341 IF (ir > 0) THEN
342 ifailg = 1
343 nuparf = mat_param(imat)%FAIL(ir)%NUPARAM
344 uparf => mat_param(imat)%FAIL(ir)%UPARAM(1:nuparf)
345c
346 CALL law104_upd(ifailg ,nparam,nuparf,uparam,uparf ,
347 . nloc_dmg,imat ,mlaw_tag(imat),ipm ,
348 . mat_param(imat))
349 ENDIF
350c------------------
351 CASE (108) ! Spring material------------------
352 CALL law108_upd(iout,titr,uparam,npc,pld,
353 . nfunc,ifunc,mat_id,func_id,
354 . pm(1,imat))
355c------------------
356 CASE (109) ! ------------------
357 IF (mat_param(imat)%NLOC > 0) THEN
358 mlaw_tag(imat)%NUVAR = 1
359 ENDIF
360c------------------
361 CASE (111) ! Marlow hyperelastic law
362c
363 !! UPARAM(13) = MULLINS(IMAT)
364 ifc = ipm(11,imat)
365 IF (ifc > 0) THEN
366 CALL law111_upd(iout ,titr ,mat_id ,uparam ,nfunc ,
367 . ifunc ,func_id ,npc ,pld ,pm(1,imat),ipm(1,imat))
368 ENDIF
369c------------------
370 CASE (112) ! Paperboard material------------------
371 CALL law112_upd(titr ,uparam,npc ,pld ,
372 . numtabl,itable,table,mat_id)
373c------------------
374 CASE (113) ! Spring material------------------
375 CALL law113_upd(iout,titr,uparam,npc,pld,
376 . nfunc,ifunc,mat_id,func_id,
377 . pm(1,imat))
378
379c------------------
380 CASE (114) ! Spring material------------------
381 CALL law114_upd(iout,titr,uparam,npc,pld,
382 . nfunc,ifunc,mat_id,func_id,
383 . pm(1,imat))
384c------------------
385 CASE (119) ! fabric material for seatbelts
386 CALL law119_upd(nparam ,numtabl ,itable ,table ,table_id,
387 . uparam,pm(1,imat),titr ,mat_id )
388c------------------
389 CASE (120) ! TAPO material for adhesives
390 CALL law120_upd(nparam ,numtabl ,itable ,table ,uparam)
391c------------------
392 CASE (129)
393 CALL law129_upd(mat_param(imat),sensors)
394c------------------
395 CASE (133)
396 CALL law133_upd(mat_param(imat),pm(1,imat),npropm)
397c-------------------------
398 CASE (158)
399 nfuncl = ipm(6,imat)
400 ifunc => ipm(10+1:10+nfunc+nfuncl,imat)
401 CALL law158_upd(mat_param(imat),titr,npc,pld,
402 . nfunc,nfuncl,ifunc,mat_id,func_id,
403 . pm(1,imat),sensors)
404c------------------
405 CASE (163)
406 CALL law163_upd(mat_param(imat),pm(1,imat),npropm)
407c------------------
408 CASE (190) !Path DEPENDE FAAM (DUBOIS)
409 CALL law190_upd(mat_param(imat) ,numtabl ,itable ,table ,pm(1,imat),
410 . npropm, ntable)
411c------------------
412c
413 END SELECT ! (ILAW)
414 pm(107,imat) = two*max(pm(32,imat),pm(100,imat)) ! used for contact stif of Icontrol>0
415 IF (ilaw==62) pm(107,imat) = hundred*pm(107,imat)
416c------------------
417 ENDDO ! imat=1,nummat
418c-----------
419 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
law100_upd_ab
subroutine law100_upd_ab(iout, titr, mat_id, uparam, nfunc, ifunc, func_id, npc, pld, pm)
Definition law100_upd.F:106
law100_upd_nht
subroutine law100_upd_nht(iout, titr, mat_id, uparam, nfunc, ifunc, func_id, npc, pld, pm)
Definition law100_upd.F:35
law104_upd
subroutine law104_upd(ifailg, nuparam, nuparf, uparam, uparf, nloc_dmg, imat, mlaw_tag, ipm, matparam)
Definition law104_upd.F:34
law108_upd
subroutine law108_upd(iout, titr, uparam, npc, pld, nfunc, ifunc, mat_id, func_id, pm)
Definition law108_upd.F:36
law111_upd
subroutine law111_upd(iout, titr, mat_id, uparam, nfunc, ifunc, func_id, npc, pld, pm, ipm)
Definition law111_upd.F:38
law112_upd
subroutine law112_upd(titr, uparam, npc, pld, numtabl, itable, table, mat_id)
Definition law112_upd.F:35
law113_upd
subroutine law113_upd(iout, titr, uparam, npc, pld, nfunc, ifunc, mat_id, func_id, pm)
Definition law113_upd.F:36
law114_upd
subroutine law114_upd(iout, titr, uparam, npc, pld, nfunc, ifunc, mat_id, func_id, pm)
Definition law114_upd.F:36
law119_upd
subroutine law119_upd(nuparam, numtabl, itable, table, table_id, uparam, pm, titr, mat_id)
Definition law119_upd.F:36
law120_upd
subroutine law120_upd(nuparam, numtabl, itable, table, uparam)
Definition law120_upd.F:31
law158_upd
subroutine law158_upd(mat_param, titr, npc, pld, nfunc, nfunl, ifunc, mat_id, func_id, pm, sensors)
Definition law158_upd.F:38
law36_upd
subroutine law36_upd(iout, titr, mat_id, nuparam, uparam, nfunc, ifunc, func_id, npc, pld, mtag, nfunct)
Definition law36_upd.F:38
law42_upd
subroutine law42_upd(mat_param, iout, titr, mat_id, pm, gama_inf)
Definition law42_upd.F:32
law58_upd
subroutine law58_upd(mat_param, titr, npc, pld, nfunc, nfunl, ifunc, mat_id, func_id, pm, sensors)
Definition law58_upd.F:40
law62_upd
subroutine law62_upd(iout, titr, mat_id, nuparam, uparam, pm, gama_inf)
Definition law62_upd.F:33
law69_upd
subroutine law69_upd(iout, titr, mat_id, uparam, nfunc, nfunct, ifunc, func_id, npc, pld, pm, ipm, gama_inf)
Definition law69_upd.F:38
law70_upd
subroutine law70_upd(mat_param, titr, mat_id, nuparam, uparam, nfunc, ifunc, npc, pld, iout, nfunct, func_id, npropm, pm)
Definition law70_upd.F:38
law76_upd
subroutine law76_upd(iout, titr, mat_id, nuparam, matparam, uparam, numtabl, itable, table, nfunc, ifunc, npc, pld)
Definition law76_upd.F:37
law77_upd
subroutine law77_upd(titr, mat_id, nuparam, mat_param, uparam, nfunc, ifunc, npc, pld)
Definition law77_upd.F:36
law90_upd
subroutine law90_upd(iout, titr, mat_id, uparam, ipm, func_id, npc, pld, pm, nfunct)
Definition law90_upd.F:36
law92_upd
subroutine law92_upd(iout, titr, mat_id, uparam, nfunc, nfunct, ifunc, func_id, npc, pld, pm, ipm)
Definition law92_upd.F:37
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
matfun_usr2sys
subroutine matfun_usr2sys(titr, mat_id, nfunc, ifunc, func_id)
Definition matfun_usr2sys.F:36
mattab_usr2sys
subroutine mattab_usr2sys(titr, mat_id, ntable, table, ntabl, itable)
Definition mattab_usr2sys.F:39
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
nlocal_reg_mod
Definition nlocal_reg_mod.F:119
table_mod
Definition table_mod.F:113
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
fretitl2
subroutine fretitl2(titr, iasc, l)
Definition freform.F:799