OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
usermat_solid.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!|| usermat_solid ../engine/source/materials/mat_share/usermat_solid.F
25!||--- called by ------------------------------------------------------
26!|| mmain ../engine/source/materials/mat_share/mmain.F90
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../engine/source/output/message/message.F
29!|| arret ../engine/source/system/arret.F
30!|| epsf2u ../engine/source/materials/mat/mat033/sigeps33.F
31!|| fail_biquad_s ../engine/source/materials/fail/biquad/fail_biquad_s.F
32!|| fail_changchang_s ../engine/source/materials/fail/changchang/fail_changchang_s.F
33!|| fail_cockroft_s ../engine/source/materials/fail/cockroft_latham/fail_cockroft_s.F
34!|| fail_composite_s ../engine/source/materials/fail/composite/fail_composite_s.F90
35!|| fail_emc ../engine/source/materials/fail/emc/fail_emc.F
36!|| fail_energy_s ../engine/source/materials/fail/energy/fail_energy_s.F
37!|| fail_gene1_s ../engine/source/materials/fail/gene1/fail_gene1_s.F
38!|| fail_hashin_s ../engine/source/materials/fail/hashin/fail_hashin_s.F
39!|| fail_hoffman_s ../engine/source/materials/fail/hoffman/fail_hoffman_s.F
40!|| fail_inievo_s ../engine/source/materials/fail/inievo/fail_inievo_s.f
41!|| fail_johnson ../engine/source/materials/fail/johnson_cook/fail_johnson.F
42!|| fail_lemaitre_s ../engine/source/materials/fail/lemaitre/fail_lemaitre_s.F90
43!|| fail_maxstrain_s ../engine/source/materials/fail/max_strain/fail_maxstrain_s.F
44!|| fail_mullins_or_s ../engine/source/materials/fail/mullins_or/fail_mullins_OR_s.F
45!|| fail_orthbiquad_s ../engine/source/materials/fail/orthbiquad/fail_orthbiquad_s.f
46!|| fail_orthenerg_s ../engine/source/materials/fail/orthenerg/fail_orthenerg_s.F
47!|| fail_orthstrain ../engine/source/materials/fail/orthstrain/fail_orthstrain_s.F
48!|| fail_puck_s ../engine/source/materials/fail/puck/fail_puck_s.f
49!|| fail_rtcl_s ../engine/source/materials/fail/rtcl/fail_rtcl_s.F
50!|| fail_sahraei_s ../engine/source/materials/fail/sahraei/fail_sahraei_s.F
51!|| fail_spalling_s ../engine/source/materials/fail/spalling/fail_spalling_s.F90
52!|| fail_syazwan_s ../engine/source/materials/fail/syazwan/fail_syazwan_s.F
53!|| fail_tab2_s ../engine/source/materials/fail/tabulated/fail_tab2_s.F
54!|| fail_tab_old_s ../engine/source/materials/fail/tabulated/fail_tab_old_s.F
55!|| fail_tab_s ../engine/source/materials/fail/tabulated/fail_tab_s.f
56!|| fail_tbutcher_s ../engine/source/materials/fail/tuler_butcher/fail_tbutcher_s.F
57!|| fail_tensstrain_s ../engine/source/materials/fail/tensstrain/fail_tensstrain_s.F
58!|| fail_tsaihill_s ../engine/source/materials/fail/tsaihill/fail_tsaihill_s.F
59!|| fail_tsaiwu_s ../engine/source/materials/fail/tsaiwu/fail_tsaiwu_s.f
60!|| fail_visual_s ../engine/source/materials/fail/visual/fail_visual_s.f
61!|| fail_wierzbicki_s ../engine/source/materials/fail/wierzbicki/fail_wierzbicki_s.F
62!|| fail_wilkins_s ../engine/source/materials/fail/wilkins/fail_wilkins_s.f
63!|| mdtsph ../engine/source/materials/mat_share/mdtsph.F
64!|| mqviscb ../engine/source/materials/mat_share/mqviscb.F
65!|| mreploc ../engine/source/materials/mat_share/mreploc.F
66!|| mrotens ../engine/source/materials/mat_share/mrotens.F
67!|| mrotensns ../engine/source/materials/mat_share/mrotens.F
68!|| nolib_usermat99 ../engine/source/user_interface/nolib_usermat99.F
69!|| nsvisul ../engine/source/materials/mat_share/nsvisul.F
70!|| sreploc3 ../engine/source/elements/solid/solide/sreploc3.F
71!|| startime ../engine/source/system/timer_mod.F90
72!|| stoptime ../engine/source/system/timer_mod.F90
73!|| viscmain ../engine/source/materials/visc/viscmain.F
74!||--- uses -----------------------------------------------------
75!|| dt_mod ../engine/source/modules/dt_mod.F
76!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
77!|| fail_composite_s_mod ../engine/source/materials/fail/composite/fail_composite_s.F90
78!|| fail_lemaitre_s_mod ../engine/source/materials/fail/lemaitre/fail_lemaitre_s.F90
79!|| fail_param_mod ../common_source/modules/mat_elem/fail_param_mod.F90
80!|| fail_spalling_s_mod ../engine/source/materials/fail/spalling/fail_spalling_s.F90
81!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
82!|| law_userso ../engine/source/user_interface/law_userso.F
83!|| mat_elem_mod ../common_source/modules/mat_elem/mat_elem_mod.F90
84!|| matparam_def_mod ../common_source/modules/mat_elem/matparam_def_mod.F90
85!|| message_mod ../engine/share/message_module/message_mod.F
86!|| names_and_titles_mod ../common_source/modules/names_and_titles_mod.F
87!|| nlocal_reg_mod ../common_source/modules/nlocal_reg_mod.F
88!|| table_mod ../engine/share/modules/table_mod.f
89!|| timer_mod ../engine/source/system/timer_mod.F90
90!||====================================================================
91 SUBROUTINE usermat_solid(TIMERS,
92 1 LFT, LLT, NFT, MTN,
93 2 JCVT, PM, OFF, SIG,
94 3 EINT, RHO, QOLD, VOL,
95 4 STRAIN, SIGL, GAMA, UVAR,
96 5 BUFMAT, TF, NPF, STIFN,
97 6 MAT, NGL, NUVAR, DT2T,
98 7 NELTST, ITYPTST, OFFG, GEO,
99 8 PID, EPSD, EL_TEMP, WXX,
100 9 WYY, WZZ, JSPH, MUMAX,
101 A SSP, AIRE, VOLN, VD2,
102 B DELTAX, VIS, D1, D2,
103 C D3, D4, D5, D6,
104 D PNEW, PSH, Q, SSP_EQ,
105 E DVOL, SOLD1, SOLD2, SOLD3,
106 F SOLD4, SOLD5, SOLD6, RX,
107 G RY, RZ, SX, SY,
108 H SZ, TX, TY, TZ,
109 I IPLA, SIGY, DEFP, ISMSTR,
110 J MFXX, MFXY, MFXZ, MFYX,
111 K MFYY, MFYZ, MFZX, MFZY,
112 L MFZZ, IPM, ISORTH, FBUF,
113 M NFAIL, NPG, SIGDD, DXY,
114 N DYX, DYZ, DZY, DZX,
115 O DXZ, FR_WAV, ISROT, V,
116 P VARNL, W, IX, X,
117 Q JTHE, ET, MSSA, DMELS,
118 R IPTR, IPTS, IPTT, TABLE,
119 S FVD2, FDELTAX, FSSP, FQVIS,
120 T IPARG1, IGEO, SIGV,
121 U AL_IMP, SIGNOR, ISTRAIN, NG,
122 V ELBUF_TAB,VBUF, ILAY, VK,
123 W IPARG, BUFVOIS, VDX, VDY,
124 X VDZ, IHET, CONDE, ITASK,
125 Y IEXPAN, VOL_AVG, AMU, EPSTH3,
126 Z EPSTH, SVISC, NEL, ETOTSH,
127 1 ISELECT, TSTAR, MUOLD, AMU2,
128 2 DPDM, RHOREF, RHOSP, NLOC_DMG,
129 3 ITY, JTUR , MAT_ELEM, IDEL7NOK,SVIS,
130 4 DT , glob_therm,damp_buf,idamp_freq_range)
131C-----------------------------------------------
132C M o d u l e s
133C-----------------------------------------------
134 USE names_and_titles_mod, ONLY: ncharline100
135 USE law_userso
136 USE table_mod
137 USE mat_elem_mod
138 USE matparam_def_mod
139 USE fail_param_mod
140 USE message_mod
141 USE nlocal_reg_mod
142 USE elbufdef_mod
143 USE dt_mod
144 use glob_therm_mod
145 USE timer_mod
146 USE fail_lemaitre_s_mod
147 USE fail_spalling_s_mod
148 USE fail_composite_s_mod
149C-----------------------------------------------
150C I m p l i c i t T y p e s
151C-----------------------------------------------
152#include "implicit_f.inc"
153C-----------------------------------------------
154C C o m m o n B l o c k s
155C-----------------------------------------------
156#include "usrplas_c.inc"
157#include "param_c.inc"
158#include "com04_c.inc"
159#include "com08_c.inc"
160#include "com01_c.inc"
161#include "userlib.inc"
162#include "impl1_c.inc"
163#include "timeri_c.inc"
164#include "scr18_c.inc"
165C-----------------------------------------------
166C D u m m y A r g u m e n t s
167C-----------------------------------------------
168 TYPE(timer_), INTENT(INOUT) :: TIMERS
169 INTEGER, INTENT(IN) :: ITY
170 INTEGER, INTENT(IN) :: JTUR
171 INTEGER, INTENT(INOUT) :: IDEL7NOK
172 INTEGER LFT ,LLT ,NFT ,MTN, NUVAR,JCVT,ISORTH,NFAIL,NPG,
173 . j, jsph, ipla,ismstr,neltst, ityptst, nuvarr,ng,ilay,
174 . ix, jthe, isrot,istrain,iptr, ipts, iptt,ihet,iexpan,nel
175 INTEGER NPF(*),MAT(MVSIZ),NGL(MVSIZ),PID(MVSIZ),IPM(NPROPMI,*),
176 . iparg1(*),iparg(nparg,*),igeo(npropgi,*),itask,iselect
177 INTEGER, INTENT(IN) :: IDAMP_FREQ_RANGE
178 TARGET :: ipm
179C
180 my_real
181 . dt2t, bufvois(*)
182 my_real
183 . pm(npropm,*), gama(*),pnew(nel),varnl(nel),
184 . off(*), sig(nel,6), sigl(nel,6), strain(nel,6), eint(*),
185 . rho(*),qold(*), vol(*), bufmat(*), tf(*),vk(*),
186 . offg(*),stifn(*),geo(npropg,*),epsd(*),
187 . wxx(*), wyy(*), wzz(*), mumax(*), ssp(*),
188 . aire(*) ,voln(*) ,vd2(*) ,deltax(*),vis(*) ,
189 . sold1(mvsiz) ,sold2(mvsiz) ,sold3(mvsiz) ,
190 . sold4(mvsiz) ,sold5(mvsiz) ,sold6(mvsiz),
191 . d1(*) ,d2(*),d3(*),d4(*),d5(*),d6(*) ,psh(nel),
192 . q(*),ssp_eq(*),dvol(*),
193 . rx(mvsiz) ,ry(mvsiz) ,rz(mvsiz),
194 . sx(mvsiz) ,sy(mvsiz) ,sz(mvsiz) ,
195 . tx(mvsiz) ,ty(mvsiz) ,tz(mvsiz) , v(*), w(*),
196 . mfxx(*) ,mfxy(*) ,mfxz(*) ,mfyx(*) ,mfyy(*) ,mfyz(*) ,
197 . mfzx(*) ,mfzy(*) ,mfzz(*), fqvis(*),fssp(*),
198 . sigv(nel,6), sigy(*),defp(nel),
199 . dxy(*) ,dyx(*) ,dyz(*) ,dzy(*) ,dzx(*) ,dxz(*), sigdd(nel,6),
200 . fr_wav(*),x(*),et(*), mssa(*), dmels(*),fvd2(*),fdeltax(*),
201 . al_imp(mvsiz),signor(mvsiz,6),vdx(*),vdy(*),vdz(*),
202 . svisc(nel,6),amu(mvsiz),amu2(mvsiz),dpdm(mvsiz),vol_avg(mvsiz),epsth(nel),epsth3(mvsiz),
203 . etotsh(nel,6),tstar(mvsiz),muold(nel),conde(*), rhoref(*), rhosp(*)
204 TARGET :: bufmat,deltax,varnl,defp
205 my_real, DIMENSION(MVSIZ,6), INTENT(INOUT) :: svis
206 my_real, DIMENSION(LLT) :: el_temp
207 my_real, TARGET, DIMENSION(NUVAR*LLT):: uvar
208 my_real, TARGET, DIMENSION(LLT):: scale1,initvarf
209 TYPE (BUF_FAIL_), TARGET :: FBUF
210 TYPE (BUF_VISC_) :: VBUF
211 TYPE (TTABLE), DIMENSION(NTABLE) :: TABLE
212 TYPE (ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
213 TYPE (NLOCAL_STR_) :: NLOC_DMG
214 TYPE (MAT_ELEM_) ,INTENT(INOUT) :: MAT_ELEM
215 TYPE (DT_), INTENT(IN) :: DT
216 TYPE (GLOB_THERM_) ,INTENT(INOUT) :: GLOB_THERM
217 TYPE (BUF_DAMP_RANGE_) ,INTENT(IN):: DAMP_BUF
218C-----------------------------------------------
219C L o c a l V a r i a b l e s
220C-----------------------------------------------
221 INTEGER NV46, NUMEL, INLOC, JSMS
222 INTEGER I,NPAR,IADBUF,NFUNC,NUMTABL,NTABL_FAIL,ISRATE,SIZNUL,MX,IPG,NPTR,NPTS,
223 . IBID,IBIDON1,IBIDON2,IBIDON3,IBIDON4 ,N48,NIX,ILAW_USER,IGTYP,
224 . NVARF,NVARVIS,IR,IRUPT,IMAT,ISVIS,IVISC,NUVARV,ISEQ,IADRUPBUF,IDEV,VISCTYPE
225 my_real
226 .E1,E2,E3,E4,E5,E6,BID1,BID3,Q1,Q2,Q3,SS1,SS2,SS3,SS4,SS5,
227 .SS6,WXXF,WYYF,WZZF,P2,EPSP,DAV,ASRATE,C1(MVSIZ),
228 . EP1(MVSIZ),EP2(MVSIZ),EP3(MVSIZ),E7(MVSIZ),
229 . EP4(2*MVSIZ),EP5(2*MVSIZ),EP6(2*MVSIZ),EINC,
230 . S1(MVSIZ) ,S2(MVSIZ) ,S3(MVSIZ) ,
231 . S4(2*MVSIZ) ,S5(2*MVSIZ) ,S6(2*MVSIZ),
232 . DE1(MVSIZ),DE2(MVSIZ),DE3(MVSIZ),
233 . DE4(2*MVSIZ),DE5(2*MVSIZ),DE6(2*MVSIZ),
234 . SO1(MVSIZ),SO2(MVSIZ),SO3(MVSIZ),
235 . SO4(2*MVSIZ),SO5(2*MVSIZ),SO6(2*MVSIZ),
236 .ES1(MVSIZ),ES2(MVSIZ),ES3(MVSIZ),ES4(MVSIZ),ES5(MVSIZ),ES6(MVSIZ),
237 .SV1(MVSIZ),SV2(MVSIZ),SV3(MVSIZ),SV4(MVSIZ),SV5(MVSIZ),SV6(MVSIZ),
238 .SVO1(MVSIZ),SVO2(MVSIZ),SVO3(MVSIZ),SVO4(MVSIZ),SVO5(MVSIZ),SVO6(MVSIZ),
239 .R11(MVSIZ),R12(MVSIZ),R13(MVSIZ),R21(MVSIZ),R22(MVSIZ),R23(MVSIZ),
240 .R31(MVSIZ),R32(MVSIZ),R33(MVSIZ),EPSP1(MVSIZ),DPLA(MVSIZ),
241 .PAIR(MVSIZ),ES7(MVSIZ),ES8(MVSIZ),ES9(MVSIZ),DEFP0(MVSIZ)
242 my_real FACQ0
243
244 my_real fpsxx(mvsiz),fpsyy(mvsiz),fpszz(mvsiz),fpsxy(mvsiz),
245 * fpsyz(mvsiz),fpszx(mvsiz),fpsyx(mvsiz),fpszy(mvsiz),
246 * fpsxz(mvsiz),bid_arr(mvsiz),
247 * upsxx(mvsiz),upsyy(mvsiz),upszz(mvsiz),upsxy(mvsiz),
248 * upsyz(mvsiz),upsxz(mvsiz)
249 my_real
250 . rho0(mvsiz),bidon1,bidon2,bidon3,bidon4,bidon5,pold
251 my_real tt_local
252 my_real, DIMENSION(NEL), TARGET :: le_max
253C----
254 my_real, DIMENSION(:), POINTER :: uparam,uparamf,uvarf,dfmax,tdel,
255 . yldfac,dam,el_len,el_pla,damini
256 my_real, DIMENSION(MVSIZ) :: el_bfrac,el_tb
257 INTEGER :: SIZBFRAC,SIZTB
258 TYPE(ULAWINTBUF) :: USERBUF
259 TYPE(L_BUFEL_) ,POINTER :: LBUF, LBUF1, LBUF2
260 TYPE(g_bufel_) ,POINTER :: GBUF
261 TYPE (FAIL_PARAM_), POINTER :: FAILPARAM
262 LOGICAL :: LOGICAL_USERL_AVAIL
263 my_real USER_UELR(MVSIZ)
264 INTEGER, DIMENSION(:) ,POINTER :: ITABL_FAIL,IPARAMF
265 INTEGER, DIMENSION(:) ,POINTER :: FLD_IDX,FOFF,IFUNC,ITABLE
266 INTEGER :: DMG_FLAG,LF_DAMMX,NIPAR
267 my_real, DIMENSION(:,:) , ALLOCATABLE :: ISMSTR10_FP
268 my_real ADDITIONAL_FLT_PARAMETERS(MVSIZ)
269 INTEGER ADDITIONAL_INT_PARAMETERS(MVSIZ),LENGTH
270 CHARACTER*256 MDS_LIBNAME
271!
272 CHARACTER OPTION*256
273 INTEGER SIZE
274 CHARACTER(LEN=NCHARLINE100)::IUSER_KEY
275 MY_REAL, DIMENSION(NEL) :: OFF_OLD
276C=======================================================================
277 GBUF => elbuf_tab(ng)%GBUF
278 lbuf => elbuf_tab(ng)%BUFLY(ilay)%LBUF(iptr,ipts,iptt)
279 logical_userl_avail=.false.
280 IF(userl_avail>0) logical_userl_avail=.true.
281 nptr = elbuf_tab(ng)%NPTR
282 npts = elbuf_tab(ng)%NPTS
283 ipg = iptr + ((ipts-1) + (iptt-1)*npts)*nptr
284 igtyp = iparg(38,ng)
285 IF ((igtyp == 20).OR.(igtyp == 21).OR.(igtyp == 22)) THEN
286 ipg = iptr + ((ipts-1) + (ilay-1)*npts)*nptr
287 ENDIF
288 iseq = elbuf_tab(ng)%BUFLY(ilay)%L_SEQ
289 inloc = iparg(78,ng)
290 dmg_flag = elbuf_tab(ng)%BUFLY(ilay)%L_DMGSCL
291 nvarvis = elbuf_tab(ng)%BUFLY(ilay)%NVAR_VISC
292 ! Make sure that the non-local variable increment is positive
293 IF (inloc > 0) THEN
294 DO i = lft,llt
295 varnl(i) = max(varnl(i),zero)
296 lbuf%PLANL(i) = lbuf%PLANL(i) + varnl(i)
297 lbuf%EPSDNL(i) = varnl(i)/max(dt1,em20)
298 ENDDO
299 ENDIF
300 ! needed in mqviscb
301 facq0 = one
302 mx = mat(1)
303 npar = ipm(9,mx)
304 iadbuf = ipm(7,mx)
305 nfunc = ipm(10,mx)
306 numtabl= ipm(226,mx)
307 uparam => bufmat(iadbuf:iadbuf+npar)
308 ifunc => ipm(10+1:10+nfunc,mx)
309 itable => ipm(226+1:226+numtabl,mx)
310c
311 !C
312
313 !c (utilisateurs pour la rupture
314 bidon1 = zero
315 bidon2 = zero
316 bidon3 = zero
317 bidon4 = zero
318 bidon5 = zero
319 bid1 = zero
320 bid3 = zero
321
322 ! (utilisateurs pour la rupture
323 ibid = 0
324 ibidon1 = 0
325 ibidon2 = 0
326 ibidon3 = 0
327 ibidon4 = 0
328 ilaw_user = ipm(217, mx)
329 imat = ipm(1, mx)
330 isvis = igeo(35,pid(1))
331 IF(impl_s >0.AND.idyna == 0) isvis = 0
332c
333 dpla(lft:llt) = zero
334
335 IF (elbuf_tab(ng)%BUFLY(ilay)%L_PLA > 0) THEN
336 defp0(lft:llt) = lbuf%PLA(lft:llt)
337 ELSE
338 defp0(lft:llt) = zero
339 ENDIF
340
341!
342 ! Save old value of OFF flag
343 off_old(1:nel) = off(1:nel)
344
345 IF (gbuf%G_BFRAC > 0) THEN
346 el_bfrac(1:nel) = elbuf_tab(ng)%GBUF%BFRAC(1:nel) ! burn fraction for high explosive
347 sizbfrac = nel
348 ELSE
349 el_bfrac(1:nel) = zero
350 sizbfrac = 0
351 ENDIF
352
353 IF (gbuf%G_TB > 0) THEN
354 el_tb(1:nel) = elbuf_tab(ng)%GBUF%TB(1:nel) ! burning time for high explosive
355 siztb = nel
356 ELSE
357 el_tb(1:nel) = zero
358 siztb = 0
359 ENDIF
360
361 !initial scale factor for yield stress defined per IPG,NPI
362 scale1(lft:llt) = one
363 yldfac => scale1(lft:llt)
364 IF ( elbuf_tab(ng)%BUFLY(ilay)%L_FAC_YLD > 0) THEN
365 yldfac => lbuf%FAC_YLD(1:nel)
366 ENDIF
367c
368 c1(1:mvsiz) = pm(32,mx)
369 rho0(1:mvsiz)= pm( 1,mx)
370 vis(1:mvsiz) = zero
371 ep1(1:mvsiz) = d1(1:mvsiz)*off(1:mvsiz)
372 ep2(1:mvsiz) = d2(1:mvsiz)*off(1:mvsiz)
373 ep3(1:mvsiz) = d3(1:mvsiz)*off(1:mvsiz)
374 ep4(1:mvsiz) = d4(1:mvsiz)*off(1:mvsiz)
375 ep5(1:mvsiz) = d5(1:mvsiz)*off(1:mvsiz)
376 ep6(1:mvsiz) = d6(1:mvsiz)*off(1:mvsiz)
377 sv1(1:mvsiz) = zero
378 sv2(1:mvsiz) = zero
379 sv3(1:mvsiz) = zero
380 sv4(1:mvsiz) = zero
381 sv5(1:mvsiz) = zero
382 sv6(1:mvsiz) = zero
383 e7(1:mvsiz) = zero
384 s1(1:mvsiz) = zero
385 s2(1:mvsiz) = zero
386 s3(1:mvsiz) = zero
387 s4(1:mvsiz) = zero
388 s5(1:mvsiz) = zero
389 s6(1:mvsiz) = zero
390 svo1(1:mvsiz) = zero
391 svo2(1:mvsiz) = zero
392 svo3(1:mvsiz) = zero
393 svo4(1:mvsiz) = zero
394 svo5(1:mvsiz) = zero
395 svo6(1:mvsiz) = zero
396
397 ssp(1:mvsiz)=zero
398
399 if (mat_elem%MAT_PARAM(mat(1))%IVISC == 3) then
400 DO i=1,nel
401 svo1(i) = svisc(i,1)
402 svo2(i) = svisc(i,2)
403 svo3(i) = svisc(i,3)
404 svo4(i) = svisc(i,4)
405 svo5(i) = svisc(i,5)
406 svo6(i) = svisc(i,6)
407 ENDDO
408 ENDIF
409
410 IF (jcvt > 0) THEN
411C---------------------------
412C ISOTROPE ET ORTHOTROPE CONVECTE
413C---------------------------
414 DO i=lft,llt
415 de1(i) = ep1(i)*dt1
416 de2(i) = ep2(i)*dt1
417 de3(i) = ep3(i)*dt1
418 de4(i) = ep4(i)*dt1
419 de5(i) = ep5(i)*dt1
420 de6(i) = ep6(i)*dt1
421 so1(i) = sig(i,1)
422 so2(i) = sig(i,2)
423 so3(i) = sig(i,3)
424 so4(i) = sig(i,4)
425 so5(i) = sig(i,5)
426 so6(i) = sig(i,6)
427 wxx(i)=zero
428 wyy(i)=zero
429 wzz(i)=zero
430 ENDDO
431C
432 ELSEIF (isorth /= 0) THEN
433C---------------------------
434C ORTHOTROPE GLOBAL
435C---------------------------
436 CALL mreploc(
437 1 gama, r11, r12, r13,
438 2 r21, r22, r23, r31,
439 3 r32, r33, rx, ry,
440 4 rz, sx, sy, sz,
441 5 tx, ty, tz, nel,
442 6 jsph)
443 DO i=lft,llt
444 ep4(i) = half*ep4(i)
445 ep5(i) = half*ep5(i)
446 ep6(i) = half*ep6(i)
447 ENDDO
448 CALL mrotens(lft,llt,ep1,ep2,ep3,ep4,ep5,ep6,
449 . r11,r12,r13,
450 . r21,r22,r23,
451 . r31,r32,r33)
452#include "vectorize.inc"
453 DO i=lft,llt
454 ep4(i) = two*ep4(i)
455 ep5(i) = two*ep5(i)
456 ep6(i) = two*ep6(i)
457 de1(i)= ep1(i)*dt1
458 de2(i)= ep2(i)*dt1
459 de3(i)= ep3(i)*dt1
460 de4(i)= ep4(i)*dt1
461 de5(i)= ep5(i)*dt1
462 de6(i)= ep6(i)*dt1
463 so1(i) = sigl(i,1)
464 so2(i) = sigl(i,2)
465 so3(i) = sigl(i,3)
466 so4(i) = sigl(i,4)
467 so5(i) = sigl(i,5)
468 so6(i) = sigl(i,6)
469 wxx(i)=zero
470 wyy(i)=zero
471 wzz(i)=zero
472 ENDDO
473 IF (mat_elem%MAT_PARAM(mat(1))%IVISC == 3) THEN
474 CALL mrotens(lft,llt,svo1,svo2,svo3,svo4,svo5,svo6,
475 . r11,r12,r13,
476 . r21,r22,r23,
477 . r31,r32,r33)
478 ENDIF
479 ELSE
480C---------------------------
481C ISOTROPE GLOBAL
482C---------------------------
483#include "vectorize.inc"
484 DO i=lft,llt
485 de1(i) = ep1(i)*dt1
486 de2(i) = ep2(i)*dt1
487 de3(i) = ep3(i)*dt1
488 de4(i) = ep4(i)*dt1
489 de5(i) = ep5(i)*dt1
490 de6(i) = ep6(i)*dt1
491 so1(i) = sig(i,1)
492 so2(i) = sig(i,2)
493 so3(i) = sig(i,3)
494 so4(i) = sig(i,4)
495 so5(i) = sig(i,5)
496 so6(i) = sig(i,6)
497 ENDDO
498C
499#include "vectorize.inc"
500 DO i=lft,llt
501 wxxf=wxx(i)*off(i)
502 wyyf=wyy(i)*off(i)
503 wzzf=wzz(i)*off(i)
504 q1=strain(i,4)*wzzf
505 q2=strain(i,6)*wyyf
506 q3=strain(i,5)*wxxf
507 ss1=strain(i,1)-q1+q2
508 ss2=strain(i,2)+q1-q3
509 ss3=strain(i,3)-q2+q3
510 ss4=strain(i,4)+2.*wzzf*(strain(i,1)-strain(i,2))+
511 . wyyf*strain(i,5)-wxxf*strain(i,6)
512 ss5=strain(i,5)+2.*wxxf*(strain(i,2)-strain(i,3))+
513 . wzzf*strain(i,6)-wyyf*strain(i,4)
514 ss6=strain(i,6)+2.*wyyf*(strain(i,3)-strain(i,1))+
515 . wxxf*strain(i,4)-wzzf*strain(i,5)
516 strain(i,1)= ss1
517 strain(i,2)= ss2
518 strain(i,3)= ss3
519 strain(i,4)= ss4
520 strain(i,5)= ss5
521 strain(i,6)= ss6
522 ENDDO
523 ENDIF
524c--------------------------
525 IF (istrain > 0) THEN
526 DO i=lft,llt
527 strain(i,1)= strain(i,1) + de1(i)
528 strain(i,2)= strain(i,2) + de2(i)
529 strain(i,3)= strain(i,3) + de3(i)
530 strain(i,4)= strain(i,4) + de4(i)
531 strain(i,5)= strain(i,5) + de5(i)
532 strain(i,6)= strain(i,6) + de6(i)
533 es1(i) = strain(i,1)
534 es2(i) = strain(i,2)
535 es3(i) = strain(i,3)
536 es4(i) = strain(i,4)
537 es5(i) = strain(i,5)
538 es6(i) = strain(i,6)
539 ENDDO
540 ELSE
541 DO i=lft,llt
542 es1(i) = zero
543 es2(i) = zero
544 es3(i) = zero
545 es4(i) = zero
546 es5(i) = zero
547 es6(i) = zero
548 ENDDO
549 ENDIF
550c-------------------------------------------------------------------------
551c TOTAL STRAIN
552c-------------------------------------------------------------------------
553 IF (ismstr == 10.OR.ismstr == 12) THEN
554c [F] = [M_F]+[1]
555c only [B] is used for ISMSTR=10,12 now
556c------ Use uniform [ES] (instead of [MF]) as other laws
557 IF (iselect > 0) THEN
558 DO i=lft,llt
559 es1(i)=etotsh(i,1)*off(i)
560 es2(i)=etotsh(i,2)*off(i)
561 es3(i)=etotsh(i,3)*off(i)
562 es4(i)=etotsh(i,4)*off(i)
563 es6(i)=etotsh(i,6)*off(i)
564 es5(i)=etotsh(i,5)*off(i)
565 END DO
566 ELSE
567 DO i=lft,llt
568 mfxx(i) = off(i)*mfxx(i)
569 mfyy(i) = off(i)*mfyy(i)
570 mfzz(i) = off(i)*mfzz(i)
571 mfxy(i) = off(i)*mfxy(i)
572 mfyx(i) = off(i)*mfyx(i)
573 mfzx(i) = off(i)*mfzx(i)
574 mfxz(i) = off(i)*mfxz(i)
575 mfyz(i) = off(i)*mfyz(i)
576 mfzy(i) = off(i)*mfzy(i)
577 ENDDO
578c Is like 42 ...etc we compute [B]=[F][F]^t - I strain----- e.g. 62,92,88,90
579 DO i=lft,llt
580 es1(i)=mfxx(i)*(two+mfxx(i))+
581 . mfxy(i)*mfxy(i)+mfxz(i)*mfxz(i)
582 es2(i)=mfyy(i)*(two+mfyy(i))+
583 . mfyx(i)*mfyx(i)+mfyz(i)*mfyz(i)
584 es3(i)=mfzz(i)*(two+mfzz(i))+
585 . mfzx(i)*mfzx(i)+mfzy(i)*mfzy(i)
586 es4(i)=mfxy(i)+mfyx(i)+mfxx(i)*mfyx(i)+
587 . mfxy(i)*mfyy(i)+mfxz(i)*mfyz(i)
588 es6(i)=mfxz(i)+mfzx(i)+mfxx(i)*mfzx(i)+
589 . mfxy(i)*mfzy(i)+mfxz(i)*mfzz(i)
590 es5(i)=mfzy(i)+mfyz(i)+mfzx(i)*mfyx(i)+
591 . mfzy(i)*mfyy(i)+mfzz(i)*mfyz(i)
592 ENDDO
593 END IF !(ISELECT>0) THEN
594c
595 IF(idtmin(1)==3.AND.ismstr == 12)THEN
596 DO i=lft,llt
597 IF (offg(i) <=one) cycle
598 es1(i)=mfxx(i)
599 es2(i)=mfyy(i)
600 es3(i)=mfzz(i)
601 es4(i)=half*(mfxy(i)+mfyx(i))
602 es6(i)=half*(mfxz(i)+mfzx(i))
603 es5(i)=half*(mfzy(i)+mfyz(i))
604 ENDDO
605 END IF
606C---------------------------
607C ORTHOTROPE GLOBAL
608C---------------------------
609 IF (jcvt <= 0 .AND. isorth /= 0) THEN
610c
611 CALL mreploc(
612 1 gama, r11, r12, r13,
613 2 r21, r22, r23, r31,
614 3 r32, r33, rx, ry,
615 4 rz, sx, sy, sz,
616 5 tx, ty, tz, nel,
617 6 jsph)
618 CALL mrotens(lft,llt,es1,es2,es3,es4,es5,es6,
619 . r11,r12,r13,
620 . r21,r22,r23,
621 . r31,r32,r33)
622 ENDIF
623 DO i=lft,llt
624 es4(i) = two*es4(i)
625 es5(i) = two*es5(i)
626 es6(i) = two*es6(i)
627 ENDDO
628c-------------------------------
629 ELSEIF (ismstr==11) THEN
630c-------------------------------
631 DO i=lft,llt
632 es1(i)=mfxx(i)
633 es2(i)=mfyy(i)
634 es3(i)=mfzz(i)
635 es4(i)=mfxy(i)+mfyx(i)
636 es6(i)=mfxz(i)+mfzx(i)
637 es5(i)=mfzy(i)+mfyz(i)
638 ENDDO
639 IF (jcvt == 0 .AND. isorth /= 0) THEN
640C---------------------------
641C ORTHOTROPE GLOBAL
642C---------------------------
643 DO i=lft,llt
644 es4(i) = half*es4(i)
645 es5(i) = half*es5(i)
646 es6(i) = half*es6(i)
647 ENDDO
648 CALL mreploc(
649 1 gama, r11, r12, r13,
650 2 r21, r22, r23, r31,
651 3 r32, r33, rx, ry,
652 4 rz, sx, sy, sz,
653 5 tx, ty, tz, nel,
654 6 jsph)
655 CALL mrotens(lft,llt,es1,es2,es3,es4,es5,es6,
656 . r11,r12,r13,
657 . r21,r22,r23,
658 . r31,r32,r33)
659 DO i=lft,llt
660 es4(i) = two*es4(i)
661 es5(i) = two*es5(i)
662 es6(i) = two*es6(i)
663 ENDDO
664 ENDIF
665C---------------------------
666 DO i=lft,llt
667 strain(i,1)= es1(i)
668 strain(i,2)= es2(i)
669 strain(i,3)= es3(i)
670 strain(i,4)= es4(i)
671 strain(i,5)= es5(i)
672 strain(i,6)= es6(i)
673 ENDDO
674c THERM STRESS COMPUTATION-----
675 IF (iexpan > 0) THEN
676 DO i=lft,llt
677 strain(i,1)= strain(i,1)-epsth(i)
678 strain(i,2)= strain(i,2)-epsth(i)
679 strain(i,3)= strain(i,3)-epsth(i)
680 ENDDO
681 ENDIF
682 ENDIF ! ismstr = 10,11,12
683
684c---------------------------------------------------------
685c Strain rate filtering coefficient
686c---------------------------------------------------------
687 israte = ipm(3,mx)
688 IF (israte > 0) THEN
689 asrate = min(one, pm(9,mx)*dt1)
690 ENDIF
691c--------------------------------------------------------
692c END TOTAL STRAIN
693c---------------------------------------------------------
694c
695c--------------------------------------------------------
696c COMPUTE UNDAMAGED EFFECTIVE STRESSES
697c---------------------------------------------------------
698 ! -> Isotropic stress softening
699 IF (dmg_flag == 1) THEN
700 DO i = 1,nel
701 so1(i) = so1(i)/max(lbuf%DMGSCL(i),em20)
702 so2(i) = so2(i)/max(lbuf%DMGSCL(i),em20)
703 so3(i) = so3(i)/max(lbuf%DMGSCL(i),em20)
704 so4(i) = so4(i)/max(lbuf%DMGSCL(i),em20)
705 so5(i) = so5(i)/max(lbuf%DMGSCL(i),em20)
706 so6(i) = so6(i)/max(lbuf%DMGSCL(i),em20)
707 ENDDO
708 ! -> Orthotropic stress softening
709 ELSEIF (dmg_flag == 6) THEN
710 DO i = 1,nel
711 so1(i) = so1(i)/max(lbuf%DMGSCL(i+nel*(1-1)),em20)
712 so2(i) = so2(i)/max(lbuf%DMGSCL(i+nel*(2-1)),em20)
713 so3(i) = so3(i)/max(lbuf%DMGSCL(i+nel*(3-1)),em20)
714 so4(i) = so4(i)/max(lbuf%DMGSCL(i+nel*(4-1)),em20)
715 so5(i) = so5(i)/max(lbuf%DMGSCL(i+nel*(5-1)),em20)
716 so6(i) = so6(i)/max(lbuf%DMGSCL(i+nel*(6-1)),em20)
717 ENDDO
718 ENDIF
719c
720c---------------------------------------------------------
721c SELECT MATERIAL LAW
722c---------------------------------------------------------
723 IF (mtn == 29) THEN
724 u_tagplas(1:nel) = 0
725 IF (logical_userl_avail)THEN
726 tt_local = tt
727 CALL eng_userlib_sigeps( mtn,
728 . llt ,npar,nuvar,nfunc,ifunc,
729 . npf ,tf ,tt_local,dt1,uparam,
730 . rho0,rho ,voln,eint,
731 . ep1 ,ep2 ,ep3 ,ep4 ,ep5 ,ep6 ,
732 . de1 ,de2 ,de3 ,de4 ,de5 ,de6 ,
733 . es1 ,es2 ,es3 ,es4 ,es5 ,es6 ,
734 . so1 ,so2 ,so3 ,so4 ,so5 ,so6 ,
735 . s1 ,s2 ,s3 ,s4 ,s5 ,s6 ,
736 . sv1 ,sv2 ,sv3 ,sv4 ,sv5 ,sv6 ,
737 . ssp ,vis ,uvar,off )
738
739 ELSE
740 ! ----------------
741 ! ERROR to be printed & exit
742 option='/MAT/LAW29 - SOLID '
743 size=len_trim(option)
744 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
745 CALL arret(2)
746 ! ----------------
747!!!
748 ENDIF
749 DO i=lft,llt
750 IF(u_tagplas(i)/=0)THEN
751 sigy(i)=u_yeld(i)
752 defp(i)=u_defp(i)
753 ELSE
754 sigy(i)=ep30
755 defp(i)=zero
756 END IF
757 END DO
758 ELSEIF(mtn == 30)THEN
759 u_tagplas(1:nel) = 0
760 IF (logical_userl_avail)THEN
761 tt_local = tt
762 CALL eng_userlib_sigeps( mtn,
763 . llt ,npar,nuvar,nfunc,ifunc,
764 . npf ,tf ,tt_local,dt1,uparam,
765 . rho0,rho ,voln,eint,
766 . ep1 ,ep2 ,ep3 ,ep4 ,ep5 ,ep6 ,
767 . de1 ,de2 ,de3 ,de4 ,de5 ,de6 ,
768 . es1 ,es2 ,es3 ,es4 ,es5 ,es6 ,
769 . so1 ,so2 ,so3 ,so4 ,so5 ,so6 ,
770 . s1 ,s2 ,s3 ,s4 ,s5 ,s6 ,
771 . sv1 ,sv2 ,sv3 ,sv4 ,sv5 ,sv6 ,
772 . ssp ,vis ,uvar,off )
773
774 ELSE
775 ! ----------------
776 ! ERROR to be printed & exit
777 option='/MAT/LAW30 - SOLID '
778 size=len_trim(option)
779 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
780 CALL arret(2)
781 ! ----------------
782!!!
783 ENDIF
784 DO i=lft,llt
785 IF(u_tagplas(i)/=0)THEN
786 sigy(i)=u_yeld(i)
787 defp(i)=u_defp(i)
788 ELSE
789 sigy(i)=ep30
790 defp(i)=zero
791 END IF
792 END DO
793 ELSEIF(mtn == 31)THEN
794 u_tagplas(1:nel) = 0
795 IF (logical_userl_avail)THEN
796 tt_local = tt
797 CALL eng_userlib_sigeps( mtn,
798 . llt ,npar,nuvar,nfunc,ifunc,
799 . npf ,tf ,tt_local,dt1,uparam,
800 . rho0,rho ,voln,eint,
801 . ep1 ,ep2 ,ep3 ,ep4 ,ep5 ,ep6 ,
802 . de1 ,de2 ,de3 ,de4 ,de5 ,de6 ,
803 . es1 ,es2 ,es3 ,es4 ,es5 ,es6 ,
804 . so1 ,so2 ,so3 ,so4 ,so5 ,so6 ,
805 . s1 ,s2 ,s3 ,s4 ,s5 ,s6 ,
806 . sv1 ,sv2 ,sv3 ,sv4 ,sv5 ,sv6 ,
807 . ssp ,vis ,uvar,off )
808
809 ELSE
810 ! ----------------
811 ! ERROR to be printed & exit
812 option='/MAT/LAW31 - SOLID '
813 size=len_trim(option)
814 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
815 CALL arret(2)
816 ! ----------------
817 ENDIF
818!
819 DO i=lft,llt
820 IF(u_tagplas(i)/=0)THEN
821 sigy(i)=u_yeld(i)
822 defp(i)=u_defp(i)
823 ELSE
824 sigy(i)=ep30
825 defp(i)=zero
826 END IF
827 END DO
828c----------------------------------------
829 ELSEIF(mtn == 99) THEN
830
831 IF(logical_userl_avail)THEN
832
833 IF (ismstr == 10) THEN
834 DO i=lft,llt
835 fpsxx(i) = one+mfxx(i)
836 fpsyy(i) = one+mfyy(i)
837 fpszz(i) = one+mfzz(i)
838 END DO
839 fpsxy = mfxy(1:mvsiz)
840 fpsyz = mfyz(1:mvsiz)
841 fpsxz = mfxz(1:mvsiz)
842 fpsyx = mfyx(1:mvsiz)
843 fpszy = mfzy(1:mvsiz)
844 fpszx = mfzx(1:mvsiz)
845 !------------
846 ! ORTHOTROPE
847 !------------
848 IF (isorth /= 0) THEN
849 CALL mreploc(
850 1 gama, r11, r12, r13,
851 2 r21, r22, r23, r31,
852 3 r32, r33, rx, ry,
853 4 rz, sx, sy, sz,
854 5 tx, ty, tz, nel,
855 6 jsph)
856 CALL mrotensns(lft,llt,
857 . fpsxx,fpsxy,fpsxz,
858 . fpsyx,fpsyy,fpsyz,
859 . fpszx,fpszy,fpszz,
860 . r11,r12,r13,
861 . r21,r22,r23,
862 . r31,r32,r33)
863 END IF ! ISORTH
864
865 CALL epsf2u(
866 1 llt ,fpsxx , fpsxy , fpsxz , fpsyx ,
867 2 fpsyy ,fpsyz , fpszx , fpszy , fpszz ,
868 3 upsxx ,upsyy , upszz , upsxy , upsyz ,
869 4 upsxz )
870
871 ELSE ! ISMSTR /= 10
872 fpsxx(1:mvsiz) = zero
873 fpsyy(1:mvsiz) = zero
874 fpszz(1:mvsiz) = zero
875 fpsxy(1:mvsiz) = zero
876 fpsyz(1:mvsiz) = zero
877 fpsxz(1:mvsiz) = zero
878 fpsyx(1:mvsiz) = zero
879 fpszy(1:mvsiz) = zero
880 fpszx(1:mvsiz) = zero
881 upsxx(1:mvsiz) = zero
882 upsyy(1:mvsiz) = zero
883 upszz(1:mvsiz) = zero
884 upsxy(1:mvsiz) = zero
885 upsyz(1:mvsiz) = zero
886 upsxz(1:mvsiz) = zero
887 END IF ! ISMSTR
888
889 ! Fill Structure in dynamical library
890 CALL eng_userlib_get_law_var(
891 * ncycle, imat,iptr, ipts,iptt,
892 * r11, r12, r13, r21, r22, r23, r31,
893 * r32, r33, so1, so2, so3, so4, so5,
894 * so6, ep1, ep2, ep3, ep4, ep5, ep6,
895 * es1, es2, es3, es4, es5, es6, de1,
896 * de2, de3, de4, de5, de6, rho0, s1,
897 * s2, s3, s4, s5, s6, sv1, sv2,
898 * sv3, sv4, sv5, sv6 )
899
900 IF (dlibtkvers >= 1301501260)THEN
901 siznul=0
902 CALL eng_userlib_get_law_var_2(
903 * fpsxx ,mvsiz ,fpsyy ,mvsiz ,
904 * fpszz ,mvsiz ,fpsxy ,mvsiz ,fpsyz ,mvsiz ,fpszx ,mvsiz,
905 * fpsyx ,mvsiz ,fpszy ,mvsiz ,fpsxz ,mvsiz ,upsxx ,mvsiz,
906 * upsyy ,mvsiz ,upszz ,mvsiz ,upsxy ,mvsiz ,upsyz ,mvsiz,
907 * upsxz ,mvsiz ,el_temp,nel ,dvol ,mvsiz ,el_bfrac,sizbfrac,
908 * el_tb ,siztb ,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
909 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
910 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
911 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
912 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
913 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
914 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
915 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
916 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
917 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
918 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
919 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
920 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
921 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
922 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
923 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
924 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
925 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
926 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul,
927 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr ,siznul)
928 ENDIF
929
930 ! Call SIGEPS99 in dynamical user library
931 tt_local = tt
932 CALL eng_userlib_sigeps99(
933 * llt ,npar ,nuvar ,ilaw_user ,nfunc,
934 * ifunc ,npf ,tf ,tt_local ,dt1,
935 * uparam ,rho ,voln ,eint ,ngl ,
936 * ssp ,vis ,uvar ,off ,sigy ,
937 * defp )
938
939 ! Get back results from structure in user library
940 CALL eng_userlib_set_law_var(
941 * s1, s2, s3, s4, s5, s6,
942 * sv1, sv2, sv3,sv4, sv5, sv6,
943 * dpla )
944
945 IF(dlibtkvers >= 1301501260)THEN
946 siznul=0
947 CALL eng_userlib_set_law_var_2(el_temp,nel,el_bfrac,sizbfrac,
948 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
949 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
950 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
951 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
952 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
953 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
954 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
955 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
956 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
957 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
958 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
959 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
960 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
961 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
962 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
963 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
964 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
965 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
966 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
967 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
968 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
969 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
970 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,
971 * bid_arr,siznul,bid_arr,siznul,bid_arr,siznul,bid_arr,siznul )
972
973 IF (sizbfrac > 0) THEN
974 ! set in element buffer the updated burn fraction for high explosive
975 elbuf_tab(ng)%GBUF%BFRAC(1:nel) = el_bfrac(1:nel)
976 ENDIF
977
978 ENDIF
979!
980 ELSE ! USERL_AVAIL = 0
981!
982 ! ----------------
983 ! ERROR to be printed & exit
984 CALL nolib_usermat99(ilaw_user, iuser_key)
985 option='/MAT/'//iuser_key//' - SOLID'
986 size=len_trim(option)
987 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
988 CALL arret(2)
989 ! ----------------
990 ENDIF ! USERL_AVAIL
991c
992#ifdef DNC
993 ELSE IF (mtn == 200) THEN ! MDS MATERIAL LAW INTERFACE (in construction)
994
995 IF(mds_avail == 1 )THEN
996 ALLOCATE(ismstr10_fp(15,nel))
997
998 CALL sreploc3(
999 1 rx, ry, rz, sx,
1000 2 sy, sz, tx, ty,
1001 3 tz, r11, r21, r31,
1002 4 r12, r22, r32, r13,
1003 5 r23, r33, llt)
1004
1005
1006 IF (ismstr == 10) THEN
1007 DO i=lft,llt
1008 fpsxx(i) = one+mfxx(i)
1009 fpsyy(i) = one+mfyy(i)
1010 fpszz(i) = one+mfzz(i)
1011 END DO
1012 fpsxy = mfxy(1:mvsiz)
1013 fpsyz = mfyz(1:mvsiz)
1014 fpsxz = mfxz(1:mvsiz)
1015 fpsyx = mfyx(1:mvsiz)
1016 fpszy = mfzy(1:mvsiz)
1017 fpszx = mfzx(1:mvsiz)
1018C---------------------------
1019C ORTHOTROPE
1020C---------------------------
1021
1022! ISMSTR10_FP( 1,1) = FPSXX(1:MVSIZ)
1023! ISMSTR10_FP( 2,1) = FPSYY(1:MVSIZ)
1024! ISMSTR10_FP( 3,1) = FPSZZ(1:MVSIZ)
1025! ISMSTR10_FP( 4,1) = FPSXY(1:MVSIZ)
1026! ISMSTR10_FP( 5,1) = FPSYZ(1:MVSIZ)
1027! ISMSTR10_FP( 6,1) = FPSXZ(1:MVSIZ)
1028! ISMSTR10_FP( 7,1) = FPSYX(1:MVSIZ)
1029! ISMSTR10_FP( 8,1) = FPSZY(1:MVSIZ)
1030! ISMSTR10_FP( 9,1) = FPSZX(1:MVSIZ)
1031! ISMSTR10_FP(10,1) = UPSXX(1:MVSIZ)
1032! ISMSTR10_FP(11,1) = UPSYY(1:MVSIZ)
1033! ISMSTR10_FP(12,1) = UPSZZ(1:MVSIZ)
1034! ISMSTR10_FP(13,1) = UPSXY(1:MVSIZ)
1035! ISMSTR10_FP(14,1) = UPSYZ(1:MVSIZ)
1036! ISMSTR10_FP(15,1) = UPSXZ(1:MVSIZ)
1037
1038 IF (isorth /= 0) THEN
1039 CALL mreploc(
1040 1 gama, r11, r12, r13,
1041 2 r21, r22, r23, r31,
1042 3 r32, r33, rx, ry,
1043 4 rz, sx, sy, sz,
1044 5 tx, ty, tz, nel,
1045 6 jsph)
1046 CALL mrotensns(lft,llt,
1047 . ismstr10_fp(1,1),ismstr10_fp(2,1),ismstr10_fp(3,1),
1048 . ismstr10_fp(4,1),ismstr10_fp(5,1),ismstr10_fp(6,1),
1049 . ismstr10_fp(7,1),ismstr10_fp(8,1),ismstr10_fp(9,1),
1050 . r11,r12,r13,
1051 . r21,r22,r23,
1052 . r31,r32,r33)
1053 END IF !(ISORTH /= 0) THEN
1054 CALL epsf2u(
1055 1 llt ,
1056 * ismstr10_fp( 1,1), ismstr10_fp( 4,1), ismstr10_fp( 6,1), ismstr10_fp( 7,1),
1057 * ismstr10_fp( 2,1), ismstr10_fp( 5,1), ismstr10_fp( 9,1), ismstr10_fp( 8,1),ismstr10_fp( 3,1),
1058 * ismstr10_fp(10,1), ismstr10_fp(11,1), ismstr10_fp(12,1), ismstr10_fp(13,1),ismstr10_fp(14,1),
1059 * ismstr10_fp(15,1) )
1060
1061 ELSE ! (ISMSTR == 10)
1062 ismstr10_fp( 1,1) = zero
1063 ismstr10_fp( 2,1) = zero
1064 ismstr10_fp( 3,1) = zero
1065 ismstr10_fp( 4,1) = zero
1066 ismstr10_fp( 5,1) = zero
1067 ismstr10_fp( 6,1) = zero
1068 ismstr10_fp( 7,1) = zero
1069 ismstr10_fp( 8,1) = zero
1070 ismstr10_fp( 9,1) = zero
1071 ismstr10_fp(10,1) = zero
1072 ismstr10_fp(11,1) = zero
1073 ismstr10_fp(12,1) = zero
1074 ismstr10_fp(13,1) = zero
1075 ismstr10_fp(14,1) = zero
1076 ismstr10_fp(15,1) = zero
1077 END if! (ISMSTR == 10)
1078
1079
1080 CALL eng_mds_sigeps_s( ncycle ,imat, ngl,nel ,
1081 * iptr , ipts , iptt , npar ,nuvar ,nfunc ,ifunc ,
1082 * npf , tf , tt , dt1 , uparam , uvar ,
1083 * rho0 , rho , voln , eint , el_temp , amu ,
1084 * ep1 , ep2 , ep3 , ep4 , ep5 , ep6 ,
1085 * de1 , de2 , de3 , de4 , de5 , de6 ,
1086 * es1 , es2 , es3 , es4 , es5 , es6 ,
1087 * so1 , so2 , so3 , so4 , so5 , so6 ,
1088 * s1 , s2 , s3 , s4 , s5 , s6 ,
1089 * sv1 , sv2 , sv3 , sv4 , sv5 , sv6 ,
1090 * ssp , vis , off , r11 , r12 , r13 ,
1091 * r21 , r22 , r23 , r31 , r32 , r33 ,
1092 * ismstr , ismstr10_fp ,defp , dpla , sigy ,
1093 * additional_flt_parameters, additional_int_parameters )
1094C
1095 ELSE
1096 CALL mds_userlib_name_get(mds_libname,length)
1097 CALL ancmsg(msgid=287,anmode=aninfo,
1098 . c1=mds_libname(1:length))
1099 CALL arret(2)
1100 ENDIF
1101#endif
1102 ENDIF ! MTN
1103C=======================================================================
1104C FAILURE MODEL
1105C=======================================================================
1106 IF ((itask==0).AND.(imon_mat==1))CALL startime(timers,121)
1107 IF (nfail > 0) THEN
1108c
1109 ! Failure criterion parameters scaling
1110 IF (inloc > 0) THEN
1111 ! -> Length used for failure criterion parameters scaling is LE_MAX
1112 le_max(1:nel) = nloc_dmg%LE_MAX(mat(1))
1113 el_len => le_max(1:nel)
1114 ELSE
1115 ! -> Length used for failure criterion parameters scaling is LE_MAX
1116 el_len => deltax(1:nel)
1117 ENDIF
1118c
1119 IF (elbuf_tab(ng)%BUFLY(ilay)%L_PLA > 0) THEN
1120C for non-local plastic law
1121 IF (inloc > 0) THEN
1122 DO i=lft,llt
1123 dpla(i) = max(varnl(i),zero)
1124 ENDDO
1125 el_pla => lbuf%PLANL(1:nel)
1126 ELSE
1127C for all plastic law
1128 DO i=lft,llt
1129 dpla(i) = defp(i) - defp0(i)
1130 ENDDO
1131 el_pla => defp(1:nel)
1132 ENDIF
1133 ELSE
1134 DO i=lft,llt
1135 dpla(i) = zero
1136 ENDDO
1137 ENDIF
1138c
1139 DO ir = 1,nfail
1140C------
1141 DO i=lft,llt
1142 dav = (ep1(i)+ep2(i)+ep3(i))*third
1143 e1 = ep1(i) - dav
1144 e2 = ep2(i) - dav
1145 e3 = ep3(i) - dav
1146 e4 = half*ep4(i)
1147 e5 = half*ep5(i)
1148 e6 = half*ep6(i)
1149 epsp1(i) =half*(e1**2+e2**2+e3**2) +e4**2+e5**2+e6**2
1150 epsp1(i) = sqrt(three*epsp1(i))*two_third
1151 ENDDO
1152C------------------------------------------------------------
1153C Recovering non-local plastic strain-rate
1154C------------------------------------------------------------
1155 IF (elbuf_tab(ng)%BUFLY(ilay)%L_PLA > 0) THEN
1156 IF (inloc > 0) THEN
1157 DO i = lft,llt
1158 epsp1(i) = lbuf%EPSDNL(i)
1159 ENDDO
1160 ENDIF
1161 ENDIF
1162c----
1163 uvarf => fbuf%FLOC(ir)%VAR
1164 irupt = fbuf%FLOC(ir)%ILAWF
1165 nvarf = fbuf%FLOC(ir)%NVAR
1166 dfmax => fbuf%FLOC(ir)%DAMMX
1167 damini => fbuf%FLOC(ir)%DAMINI
1168 tdel => fbuf%FLOC(ir)%TDEL
1169 fld_idx => fbuf%FLOC(ir)%INDX
1170 foff => fbuf%FLOC(ir)%OFF
1171 dam => fbuf%FLOC(ir)%DAM
1172 lf_dammx = fbuf%FLOC(ir)%LF_DAMMX
1173c
1174 failparam => mat_elem%MAT_PARAM(mx)%FAIL(ir)
1175 npar = mat_elem%MAT_PARAM(mx)%FAIL(ir)%NUPARAM
1176 nipar = mat_elem%MAT_PARAM(mx)%FAIL(ir)%NIPARAM
1177 uparamf => mat_elem%MAT_PARAM(mx)%FAIL(ir)%UPARAM(1:npar)
1178 iparamf => mat_elem%MAT_PARAM(mx)%FAIL(ir)%IPARAM(1:nipar)
1179 nfunc = mat_elem%MAT_PARAM(mx)%FAIL(ir)%NFUNC
1180 ifunc => mat_elem%MAT_PARAM(mx)%FAIL(ir)%IFUNC(1:nfunc)
1181 ntabl_fail = mat_elem%MAT_PARAM(mx)%FAIL(ir)%NTABLE
1182 itabl_fail => mat_elem%MAT_PARAM(mx)%FAIL(ir)%TABLE(1:ntabl_fail)
1183c----
1184 IF (irupt == 1)THEN
1185C---- Johnson cook
1186 CALL fail_johnson(llt ,npar,nvarf,
1187 2 tt ,dt1 ,uparamf,ngl ,
1188 4 s1 ,s2 ,s3 ,s4 ,s5 ,s6,
1189 5 dpla,epsp1,tstar,uvarf,off,
1190 6 dfmax,tdel )
1191 ELSEIF (irupt == 2) THEN
1192C Tuler Butcher
1193 CALL fail_tbutcher_s(llt ,npar,nvarf,
1194 2 tt ,dt1 ,uparamf,ngl ,
1195 4 s1 ,s2 ,s3 ,s4 ,s5 ,s6,
1196 5 uvarf,off ,dfmax,tdel )
1197 ELSEIF (irupt == 3) THEN
1198C wilkins
1199 CALL fail_wilkins_s(llt ,npar,nuvarr,
1200 2 tt ,dt1 ,uparamf,ngl ,
1201 4 s1 ,s2 ,s3 ,s4 ,s5 ,s6,
1202 5 dpla,uvarf,off,dfmax,tdel )
1203 ELSEIF (irupt == 4)THEN
1204C user1
1205 IF (logical_userl_avail)THEN
1206 user_uelr(1:llt)=gbuf%UELR(1:llt)
1207
1208 tt_local = tt
1209 CALL eng_userlib_flaw(irupt,llt ,npar,nuvarr,nfunc,ifunc,
1210 2 npf ,tf ,tt_local ,dt1 ,uparamf,
1211 3 ngl ,ibidon1,ibidon2,ibidon3 ,ibidon4,
1212 4 ep1 ,ep2 ,ep3 ,ep4 ,ep5 ,ep6 ,
1213 5 es1 ,es2 ,es3 ,es4 ,es5 ,es6 ,
1214 6 s1 ,s2 ,s3 ,s4 ,s5 ,s6 ,
1215 7 defp ,dpla,epsp1,uvarf,off ,
1216 8 deltax,voln,user_uelr,dfmax,bidon5)
1217 gbuf%UELR(1:llt) = user_uelr(1:llt)
1218 ELSE
1219! CALL F04LAW(LLT ,NPAR,NUVARR,NFUNC,IFUNC,
1220! 2 NPF ,TF ,TT ,DT1 ,UPARAMF,
1221! 3 NGL ,IBIDON1,IBIDON2,IBIDON3 ,IBIDON4,
1222! 4 EP1 ,EP2 ,EP3 ,EP4 ,EP5 ,EP6 ,
1223! 5 ES1 ,ES2 ,ES3 ,ES4 ,ES5 ,ES6 ,
1224! 6 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1225! 7 DEFP ,DPLA,EPSP1,UVARF,OFF ,
1226! 8 deltax,voln,gbuf%UELR,dfmax,bidon5)
1227!!!
1228 ! ----------------
1229 ! ERROR to be printed & exit
1230 option='/FAIL/USER1 - SOLID '
1231 size=len_trim(option)
1232 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
1233 CALL arret(2)
1234 ! ----------------
1235!!!
1236 ENDIF
1237 ELSEIF (irupt == 5)THEN
1238C user2
1239 IF (logical_userl_avail)THEN
1240 tt_local = tt
1241 user_uelr(1:llt)=gbuf%UELR(1:llt)
1242 CALL eng_userlib_flaw(irupt,llt ,npar,nuvarr,nfunc,ifunc,
1243 2 npf ,tf ,tt_local ,dt1 ,uparamf,
1244 3 ngl ,ibidon1,ibidon2,ibidon3 ,ibidon4,
1245 4 ep1 ,ep2 ,ep3 ,ep4 ,ep5 ,ep6 ,
1246 5 es1 ,es2 ,es3 ,es4 ,es5 ,es6 ,
1247 6 s1 ,s2 ,s3 ,s4 ,s5 ,s6 ,
1248 7 defp ,dpla,epsp1,uvarf,off ,
1249 8 deltax,voln,user_uelr,dfmax,bidon5)
1250 gbuf%UELR(1:llt) = user_uelr(1:llt)
1251 ELSE
1252! CALL F05LAW(LLT ,NPAR,NUVARR,NFUNC,IFUNC,
1253! 2 NPF ,TF ,TT ,DT1 ,UPARAMF,
1254! 3 NGL ,IBIDON1,IBIDON2,IBIDON3 ,IBIDON4,
1255! 4 EP1 ,EP2 ,EP3 ,EP4 ,EP5 ,EP6 ,
1256! 5 ES1 ,ES2 ,ES3 ,ES4 ,ES5 ,ES6 ,
1257! 6 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1258! 7 DEFP ,DPLA,EPSP1,UVARF,OFF ,
1259! 8 DELTAX ,VOLN,GBUF%UELR,DFMAX,BIDON5)
1260!!!
1261 ! ----------------
1262 ! ERROR to be printed & exit
1263 option='/FAIL/USER2 - SOLID '
1264 size=len_trim(option)
1265 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
1266 CALL arret(2)
1267 ! ----------------
1268!!!
1269 ENDIF
1270 ELSEIF (irupt == 6)THEN
1271C user3
1272 IF (logical_userl_avail)THEN
1273 tt_local = tt
1274 user_uelr(1:llt)=gbuf%UELR(1:llt)
1275 CALL eng_userlib_flaw(irupt,llt ,npar,nuvarr,nfunc,ifunc,
1276 2 npf ,tf ,tt_local ,dt1 ,uparamf,
1277 3 ngl ,ibidon1,ibidon2,ibidon3 ,ibidon4,
1278 4 ep1 ,ep2 ,ep3 ,ep4 ,ep5 ,ep6 ,
1279 5 es1 ,es2 ,es3 ,es4 ,es5 ,es6 ,
1280 6 s1 ,s2 ,s3 ,s4 ,s5 ,s6 ,
1281 7 defp ,dpla,epsp1,uvarf,off ,
1282 8 deltax,voln,user_uelr,dfmax,bidon5)
1283 gbuf%UELR(1:llt) = user_uelr(1:llt)
1284 ELSE
1285! CALL f06law(llt ,npar,nuvarr,nfunc,ifunc,
1286! 2 NPF ,TF ,TT ,DT1 ,UPARAMF,
1287! 3 NGL ,IBIDON1,IBIDON2,IBIDON3 ,IBIDON4,
1288! 4 EP1 ,EP2 ,EP3 ,EP4 ,EP5 ,EP6 ,
1289! 5 ES1 ,ES2 ,ES3 ,ES4 ,ES5 ,ES6 ,
1290! 6 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1291! 7 DEFP ,DPLA,EPSP1,UVARF,OFF ,
1292! 8 DELTAX,VOLN,GBUF%UELR,DFMAX,BIDON5)
1293!!!
1294 ! ----------------
1295 ! ERROR to be printed & exit
1296 option='/fail/user3 - solid '
1297 SIZE=LEN_TRIM(OPTION)
1298 CALL ANCMSG(MSGID=257,C1=OPTION(1:SIZE),ANMODE=ANINFO)
1299 CALL ARRET(2)
1300 ! ----------------
1301!!!
1302 ENDIF
1303 ELSEIF (IRUPT == 7)THEN
1304 ELSEIF (IRUPT == 8)THEN
1305C---- Johnson cook + spalling
1306 CALL FAIL_SPALLING_S(LLT ,NPAR,
1307 2 TT ,UPARAMF,NGL ,
1308 4 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1309 5 DPLA,EPSP1,TSTAR,OFF ,
1310 6 LF_DAMMX ,DFMAX,TDEL ,OFFG ,
1311 7 NIPAR ,IPARAMF,MVSIZ)
1312 ELSEIF (IRUPT == 9)THEN
1313C---- wierzbicki
1314 CALL FAIL_WIERZBICKI_S(LLT ,NPAR,NUVARR,
1315 2 TT ,DT1 ,UPARAMF,NGL ,
1316 4 S1 ,S2 ,S3 ,S4 ,S5 ,S6,
1317 5 DPLA,DEFP,UVARF,OFF ,DFMAX,
1318 6 TDEL )
1319 ELSEIF (IRUPT == 10)THEN
1320C---- strain tension
1321 CALL FAIL_TENSSTRAIN_S(LLT ,NPAR,NVARF,NFUNC,IFUNC ,
1322 2 NPF ,TF ,TT ,DT1 ,UPARAMF,
1323 3 NGL ,DELTAX ,TSTAR ,ISMSTR,
1324 4 ES1 ,ES2 ,ES3 ,ES4 ,ES5 ,ES6 ,
1325 5 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1326 6 EPSP1,UVARF,OFF ,DFMAX,TDEL ,
1327 7 MFXX ,MFXY ,MFXZ ,MFYX ,MFYY ,MFYZ ,
1328 8 MFZX ,MFZY ,MFZZ ,LBUF%DMGSCL)
1329 ELSEIF (IRUPT == 11)THEN
1330C---- energy failure
1331 CALL FAIL_ENERGY_S(
1332 1 NEL ,NPAR ,NVARF ,NFUNC ,IFUNC ,NPF ,
1333 2 TF ,TT ,DT1 ,UPARAMF ,NGL ,EPSP1 ,
1334 3 UVARF ,OFF ,DFMAX ,TDEL ,LBUF%DMGSCL,
1335 4 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1336 5 DE1 ,DE2 ,DE3 ,DE4 ,DE5 ,DE6 )
1337 ELSEIF (IRUPT == 13) THEN
1338C---- chang - chang
1339 CALL FAIL_CHANGCHANG_S(
1340 1 NEL ,NPAR ,NVARF ,UPARAMF ,UVARF ,
1341 2 TT ,IPG ,ILAY ,NPG ,NGL ,
1342 3 LBUF%DMGSCL,DFMAX ,OFF ,LBUF%OFF ,GBUF%NOFF,
1343 4 S1 ,S2 ,S3 ,S4 ,S6 ,
1344 5 TDEL ,LF_DAMMX )
1345 ELSEIF(IRUPT == 14)THEN
1346c --- Hashin failure model
1347 DO I=1,NEL
1348 EPSP1(I) = MAX(ABS(EP1(I)),ABS(EP2(I)),ABS(EP3(I)),EM20)
1349 ENDDO
1350 CALL FAIL_HASHIN_S(
1351 1 NEL ,NVARF ,ILAY ,NPG ,TT ,
1352 2 DT1 ,UPARAMF ,NGL ,OFF ,GBUF%NOFF,
1353 3 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1354 4 UVARF ,NPAR ,DFMAX ,TDEL ,EPSP1 ,LF_DAMMX )
1355 ELSEIF(IRUPT == 16)THEN
1356c --- modified Puck failure model
1357 CALL FAIL_PUCK_S(
1358 1 NEL ,NVARF ,ILAY ,NPG ,TT ,
1359 2 DT1 ,UPARAMF ,NGL ,OFF ,GBUF%NOFF,
1360 3 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1361 4 UVARF ,NPAR ,DFMAX ,LF_DAMMX ,TDEL )
1362 ELSEIF (IRUPT == 23) THEN
1363C---- tabulated failure model
1364 CALL FAIL_TAB_S(
1365 1 LLT ,NVARF ,NPF ,TF ,TT ,
1366 2 UPARAMF ,NGL ,EL_LEN ,
1367 3 S1 ,S2 ,S3 ,S4 ,S5 ,S6,
1368 4 DPLA ,EPSP1 ,TSTAR ,UVARF ,NTABL_FAIL,ITABL_FAIL,
1369 5 OFF ,TABLE ,DFMAX ,TDEL ,NFUNC ,IFUNC )
1370 ELSEIF (IRUPT == 24) THEN
1371c --- Orthotropic strain failure
1372 CALL FAIL_ORTHSTRAIN(
1373 1 LLT ,NPAR ,NVARF ,NFUNC ,IFUNC ,
1374 2 NPF ,TF ,TT ,DT1 ,UPARAMF ,ISMSTR ,
1375 3 EP1 ,EP2 ,EP3 ,EP4 ,EP5 ,EP6 ,
1376 4 ES1 ,ES2 ,ES3 ,ES4 ,ES5 ,ES6 ,
1377 5 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1378 6 UVARF ,OFF ,IPG ,NGL ,DFMAX ,TDEL ,
1379 7 GBUF%UELR,NPG ,DELTAX ,LF_DAMMX )
1380 ELSEIF (IRUPT == 27) THEN
1381c --- extended mohr coulomb failure model
1382 CALL FAIL_EMC(
1383 1 LLT ,NVARF ,TT ,
1384 2 DT1 ,UPARAMF ,NGL ,
1385 3 S1 ,S2 ,S3 ,S4 ,S5 ,S6,
1386 4 DEFP ,DPLA ,EPSP1 ,UVARF ,
1387 5 OFF ,DFMAX ,TDEL )
1388 ELSEIF (IRUPT == 29) THEN
1389c --- MIT Wierzbicki Sahraei electric battery failure
1390 CALL FAIL_SAHRAEI_S(
1391 1 LLT ,NFUNC ,IFUNC ,NPF ,TF ,
1392 2 TT ,NGL ,UPARAMF ,
1393 3 ES1 ,ES2 ,ES3 ,ES4 ,ES5 ,ES6 ,
1394 4 OFF ,DFMAX ,TDEL ,DELTAX ,
1395 5 NVARF ,UVARF )
1396 ELSEIF (IRUPT == 30) THEN
1397c --- Biquadratic failure model
1398 CALL FAIL_BIQUAD_S(
1399 1 LLT ,NPAR ,NVARF ,NFUNC ,IFUNC ,EL_LEN ,
1400 2 NPF ,TF ,TT ,UPARAMF ,TDEL ,
1401 3 NGL ,DPLA ,UVARF ,OFF ,DFMAX ,LBUF%DMGSCL,
1402 4 S1 ,S2 ,S3 ,S4 ,S5 ,S6 )
1403 ELSEIF (IRUPT == 33) THEN
1404c --- Mullins Ogden-Roxburgh damage model
1405 CALL FAIL_MULLINS_OR_S(
1406 1 LLT ,NPAR ,NVARF ,
1407 2 TT ,DT1 ,UPARAMF ,NGL ,
1408 3 S1 ,S2 ,S3 ,S4 ,
1409 4 S5 ,S6 ,UVARF ,OFF ,
1410 5 DFMAX ,LBUF%EPSA)
1411 ELSEIF (IRUPT == 34) THEN
1412C --- Cockroft-Latham failure model
1413 CALL FAIL_COCKROFT_S(LLT ,NPAR ,NVARF,
1414 1 TT ,DT1 ,UPARAMF ,NGL ,
1415 3 EP1 ,EP2 ,EP3 ,EP4 ,EP5 ,EP6 ,
1416 4 ES1 ,ES2 ,ES3 ,ES4 ,ES5 ,ES6 ,
1417 5 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1418 6 DEFP ,DPLA ,EPSP ,UVARF ,OFF ,
1419 7 DFMAX ,TDEL )
1420 ELSEIF (IRUPT == 36) THEN
1421C --- VISUAL failure model
1422 CALL FAIL_VISUAL_S(
1423 . LLT ,NPAR ,NVARF ,TT ,DT1 ,UPARAMF,
1424 . ES1 ,ES2 ,ES3 ,ES4 ,ES5 ,ES6 ,
1425 . S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1426 . UVARF ,OFF ,NGL ,DFMAX ,ISMSTR )
1427c
1428 ELSEIF (IRUPT == 37) THEN
1429c --- tabulated failure model (old, obsolete version)
1430 CALL FAIL_TAB_OLD_S(
1431 1 LLT ,NVARF ,NPF ,TF ,TT ,
1432 2 UPARAMF ,NGL ,EL_LEN ,
1433 3 S1 ,S2 ,S3 ,S4 ,S5 ,S6,
1434 4 DEFP ,DPLA ,EPSP1 ,TSTAR ,UVARF ,
1435 5 OFF ,DFMAX ,TDEL ,
1436 6 NFUNC ,IFUNC )
1437c
1438 ELSEIF (IRUPT == 38) THEN
1439c --- Orthotropic biquadratic failure model
1440 CALL FAIL_ORTHBIQUAD_S(
1441 1 LLT ,NPAR ,NVARF ,NFUNC ,IFUNC ,
1442 2 NPF ,TF ,TT ,DT1 ,UPARAMF ,
1443 3 NGL ,DPLA ,EPSP1 ,UVARF ,OFF ,
1444 4 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1445 5 DFMAX ,TDEL ,EL_LEN )
1446c
1447 ELSEIF (IRUPT == 39) THEN
1448c --- GENE1 failure model
1449 CALL FAIL_GENE1_S(
1450 1 LLT ,NPAR ,NVARF ,NFUNC ,IFUNC ,LBUF%OFF ,
1451 2 NPF ,TF ,TT ,DT1 ,UPARAMF ,IPG ,
1452 3 NGL ,GBUF%DT ,EPSP1 ,UVARF ,OFF ,NPG ,
1453 4 ES1 ,ES2 ,ES3 ,ES4 ,ES5 ,ES6 ,
1454 5 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1455 6 EL_TEMP ,VOLN ,DFMAX ,TDEL ,DELTAX ,TABLE ,
1456 7 IR ,ELBUF_TAB(NG),ILAY ,NTABL_FAIL,ITABL_FAIL,LF_DAMMX,
1457 8 NIPAR ,IPARAMF )
1458c
1459 ELSEIF (IRUPT == 40) THEN
1460c --- RTCL failure model
1461 CALL FAIL_RTCL_S(
1462 1 LLT ,NPAR ,NVARF ,TT ,DT1 ,UPARAMF ,
1463 2 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1464 3 NGL ,DPLA ,UVARF ,OFF ,DFMAX ,TDEL )
1465c
1466 ELSEIF (IRUPT == 41) THEN
1467C---- Tabulated failure model version 2
1468 CALL FAIL_TAB2_S(
1469 1 LLT ,NPAR ,NVARF ,NFUNC ,IFUNC ,
1470 2 NPF ,TABLE ,TF ,TT ,UPARAMF ,
1471 3 NGL ,EL_LEN ,DPLA ,EPSP1 ,UVARF ,
1472 5 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1473 6 EL_TEMP ,OFF ,DFMAX ,TDEL ,LBUF%DMGSCL,
1474 7 GBUF%UELR,IPG ,NPG ,LBUF%OFF ,NTABL_FAIL,ITABL_FAIL,
1475 8 GBUF%NOFF,VOLN )
1476c
1477 ELSEIF (IRUPT == 42) THEN
1478C---- Inievo failure model
1479 CALL FAIL_INIEVO_S(
1480 1 NEL ,NPAR ,NVARF ,
1481 2 TABLE ,NTABL_FAIL,ITABL_FAIL,TT ,UPARAMF ,
1482 3 NGL ,EL_LEN ,DPLA ,EPSP1 ,UVARF ,
1483 5 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1484 6 EL_PLA ,EL_TEMP ,SIGY ,OFF ,DFMAX ,
1485 7 TDEL ,LBUF%DMGSCL,GBUF%UELR,IPG ,NPG ,
1486 8 LBUF%OFF ,DAMINI ,GBUF%VOL ,INLOC )
1487c
1488 ELSEIF (IRUPT == 43) THEN
1489c --- Syazwan failure model
1490 CALL FAIL_SYAZWAN_S(
1491 1 NEL ,UPARAMF ,NPAR ,UVARF ,NVARF ,
1492 2 TT ,NGL ,IPG ,DPLA ,TDEL ,
1493 3 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1494 4 DFMAX ,NFUNC ,IFUNC ,EL_LEN ,OFF ,
1495 5 NPF ,TF ,GBUF%UELR,NPG ,LBUF%OFF)
1496 ELSEIF (IRUPT == 44) THEN
1497C---- Tsai-Wu failure model
1498 CALL FAIL_TSAIWU_S(
1499 1 NEL ,NVARF ,IPG ,ILAY ,NPG ,TT ,
1500 2 DT1 ,UPARAMF ,NGL ,OFF ,LBUF%OFF ,GBUF%NOFF,
1501 3 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1502 4 UVARF ,NPAR ,DFMAX ,LF_DAMMX ,TDEL ,LBUF%DMGSCL)
1503c
1504 ELSEIF (IRUPT == 45) THEN
1505C---- Tsai-Hill failure model
1506 CALL FAIL_TSAIHILL_S(
1507 1 NEL ,NVARF ,IPG ,ILAY ,NPG ,TT ,
1508 2 DT1 ,UPARAMF ,NGL ,OFF ,LBUF%OFF ,GBUF%NOFF,
1509 3 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1510 4 UVARF ,NPAR ,DFMAX ,LF_DAMMX ,TDEL ,LBUF%DMGSCL)
1511c
1512 ELSEIF (IRUPT == 46) THEN
1513C---- Hoffman failure model
1514 CALL FAIL_HOFFMAN_S(
1515 1 NEL ,NVARF ,IPG ,ILAY ,NPG ,TT ,
1516 2 DT1 ,UPARAMF ,NGL ,OFF ,LBUF%OFF ,GBUF%NOFF,
1517 3 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1518 4 UVARF ,NPAR ,DFMAX ,LF_DAMMX ,TDEL ,LBUF%DMGSCL)
1519c
1520 ELSEIF (IRUPT == 47) THEN
1521C---- Maximum strain failure model
1522 CALL FAIL_MAXSTRAIN_S(
1523 1 NEL ,NVARF ,IPG ,ILAY ,NPG ,TT ,
1524 2 DT1 ,UPARAMF ,NGL ,OFF ,LBUF%OFF ,GBUF%NOFF,
1525 3 ES1 ,ES2 ,ES3 ,ES4 ,ES5 ,ES6 ,
1526 4 UVARF ,NPAR ,DFMAX ,LF_DAMMX ,TDEL ,LBUF%DMGSCL)
1527c
1528 ELSEIF (IRUPT == 48) THEN
1529c---- Orthotropic energy failure
1530 CALL FAIL_ORTHENERG_S(
1531 1 NEL ,NPAR ,NVARF ,UPARAMF ,UVARF ,NGL ,
1532 2 NPG ,IPG ,ILAY ,OFF ,LBUF%OFF ,GBUF%NOFF,
1533 3 DE1 ,DE2 ,DE3 ,DE4 ,DE5 ,DE6 ,
1534 5 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1535 6 TT ,TDEL ,DFMAX ,DELTAX ,LBUF%DMGSCL)
1536c
1537 ELSEIF (IRUPT == 50) THEN
1538c---- Lemaitre damage model
1539 CALL FAIL_LEMAITRE_S(
1540 1 NEL ,NPAR ,UPARAMF ,TT,NGL,MAT_ELEM%MAT_PARAM(MX),
1541 2 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1542 3 DPLA ,EL_PLA ,LBUF%OFF ,OFF ,DFMAX ,TDEL ,
1543 4 NIPAR ,IPARAMF ,LBUF%DMGSCL,GBUF%NOFF,NPG )
1544c
1545 ELSEIF (IRUPT == 51) THEN
1546c---- Composite failure model
1547 CALL FAIL_COMPOSITE_S(
1548 1 NEL ,FAILPARAM,NVARF ,UVARF ,TT ,NGL ,
1549 2 IPG ,ILAY ,NPG ,TDEL ,OFF ,LBUF%OFF ,
1550 3 S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1551 4 LBUF%DMGSCL,LF_DAMMX,DFMAX ,GBUF%NOFF)
1552c
1553 ENDIF
1554 ENDDO ! IR = 1,NFAIL
1555 ENDIF ! NFAIL
1556.AND. IF ((ITASK==0)(IMON_MAT==1))CALL STOPTIME(TIMERS,121)
1557C-----------------------------------------------------------------------
1558C viscous stress (/VISC models + Damping frequency range)
1559C-----------------------------------------------------------------------
1560.OR. IF ((MAT_ELEM%MAT_PARAM(MAT(1))%IVISC > 0)(IDAMP_FREQ_RANGE > 0)) THEN
1561 CALL VISCMAIN(MAT_ELEM%MAT_PARAM(MAT(1))%VISC ,NEL ,
1562 . NVARVIS ,VBUF%VAR,RHO0 ,VIS ,SSP ,DT1 ,
1563 . EP1 ,EP2 ,EP3 ,EP4 ,EP5 ,EP6 ,
1564 . DE1 ,DE2 ,DE3 ,DE4 ,DE5 ,DE6 ,
1565 . SVO1 ,SVO2 ,SVO3 ,SVO4 ,SVO5 ,SVO6 ,
1566 . SV1 ,SV2 ,SV3 ,SV4 ,SV5 ,SV6 ,
1567 . MFXX ,MFXY ,MFXZ ,MFYX ,MFYY ,MFYZ ,
1568 . MFZX ,MFZY ,MFZZ ,
1569 . S1 ,S2 ,S3 ,S4 ,S5 ,S6 ,
1570 . DAMP_BUF,IDAMP_FREQ_RANGE,MVSIZ,ET,MAT_ELEM%MAT_PARAM(IMAT)%YOUNG,
1571 . MAT_ELEM%MAT_PARAM(IMAT)%SHEAR,MAT_ELEM%MAT_PARAM(MAT(1))%IVISC)
1572 ENDIF
1573C---------------------------
1574C
1575C calcul of viscosity ( Navier-Stokes)
1576C
1577 IF (ISVIS > 0) THEN
1578 CALL NSVISUL(LLT ,OFF ,RHO ,GEO ,
1579 2 PID ,SSP ,AIRE ,VOLN ,EP1 ,
1580 3 EP2 ,EP3 ,EP4 ,EP5 ,EP6 ,
1581 4 SV1 ,SV2 ,SV3 ,SV4 ,SV5 ,
1582 5 SV6 ,S3 ,ES3 ,RHO0 ,RHOREF )
1583 ENDIF
1584c--------------------------------------------------------
1585c Shooting nodes algorithm activation
1586c--------------------------------------------------------
1587 DO I = 1,NEL
1588.AND. IF ((OFF_OLD(I) > ZERO) (OFF(I) == ZERO)) THEN
1589 IDEL7NOK = 1
1590 ENDIF
1591 ENDDO
1592c--------------------------------------------------------
1593c DAMAGED STRESSES
1594c---------------------------------------------------------
1595 ! -> Isotropic stress softening
1596 IF (DMG_FLAG == 1) THEN
1597 DO I = 1,NEL
1598 S1(I) = S1(I)*LBUF%DMGSCL(I)
1599 S2(I) = S2(I)*LBUF%DMGSCL(I)
1600 S3(I) = S3(I)*LBUF%DMGSCL(I)
1601 S4(I) = S4(I)*LBUF%DMGSCL(I)
1602 S5(I) = S5(I)*LBUF%DMGSCL(I)
1603 S6(I) = S6(I)*LBUF%DMGSCL(I)
1604 ENDDO
1605 ! -> Orthotropic stress softening
1606 ELSEIF (DMG_FLAG == 6) THEN
1607 DO I = 1,NEL
1608 S1(I) = S1(I)*LBUF%DMGSCL(I+NEL*(1-1))
1609 S2(I) = S2(I)*LBUF%DMGSCL(I+NEL*(2-1))
1610 S3(I) = S3(I)*LBUF%DMGSCL(I+NEL*(3-1))
1611 S4(I) = S4(I)*LBUF%DMGSCL(I+NEL*(4-1))
1612 S5(I) = S5(I)*LBUF%DMGSCL(I+NEL*(5-1))
1613 S6(I) = S6(I)*LBUF%DMGSCL(I+NEL*(6-1))
1614 ENDDO
1615 ENDIF
1616c---------------------------------------------------------
1617 IF (ISORTH /= 0) THEN
1618 if (jcvt == 0) then
1619 DO I=1,NEL
1620 SIGL(I,1) = S1(I)*OFF(I)
1621 SIGL(I,2) = S2(I)*OFF(I)
1622 SIGL(I,3) = S3(I)*OFF(I)
1623 SIGL(I,4) = S4(I)*OFF(I)
1624 SIGL(I,5) = S5(I)*OFF(I)
1625 SIGL(I,6) = S6(I)*OFF(I)
1626 SV1(I) = SV1(I)*OFF(I)
1627 SV2(I) = SV2(I)*OFF(I)
1628 SV3(I) = SV3(I)*OFF(I)
1629 SV4(I) = SV4(I)*OFF(I)
1630 SV5(I) = SV5(I)*OFF(I)
1631 SV6(I) = SV6(I)*OFF(I)
1632 ENDDO
1633!
1634 call mrotens(1,nel,
1635 & s1 ,s2 ,s3 , s4 ,s5 ,s6 ,
1636 & r11,r21,r31,
1637 & r12,r22,r32,
1638 & r13,r23,r33)
1639 call mrotens(1,nel,
1640 & sv1 ,sv2 ,sv3 , sv4 ,sv5 ,sv6 ,
1641 & r11,r21,r31,
1642 & r12,r22,r32,
1643 & r13,r23,r33)
1644 else ! jcvt > 0
1645 DO I=1,NEL
1646 SV1(I) = SV1(I)*OFF(I)
1647 SV2(I) = SV2(I)*OFF(I)
1648 SV3(I) = SV3(I)*OFF(I)
1649 SV4(I) = SV4(I)*OFF(I)
1650 SV5(I) = SV5(I)*OFF(I)
1651 SV6(I) = SV6(I)*OFF(I)
1652 ENDDO
1653 end if
1654 END IF
1655!
1656 DO I=1,NEL
1657 SIG(I,1) = S1(I)*OFF(I)
1658 SIG(I,2) = S2(I)*OFF(I)
1659 SIG(I,3) = S3(I)*OFF(I)
1660 SIG(I,4) = S4(I)*OFF(I)
1661 SIG(I,5) = S5(I)*OFF(I)
1662 SIG(I,6) = S6(I)*OFF(I)
1663 SVIS(I,1) = SV1(I)*OFF(I)
1664 SVIS(I,2) = SV2(I)*OFF(I)
1665 SVIS(I,3) = SV3(I)*OFF(I)
1666 SVIS(I,4) = SV4(I)*OFF(I)
1667 SVIS(I,5) = SV5(I)*OFF(I)
1668 SVIS(I,6) = SV6(I)*OFF(I)
1669 ENDDO
1670C
1671C isvis for outp
1672C
1673 VISCTYPE = MAT_ELEM%MAT_PARAM(MAT(1))%IVISC
1674.OR..OR. IF (ISVIS > 0 VISCTYPE == 1 VISCTYPE == 3) THEN
1675 DO I=1,NEL
1676 SIGV(I,1) = SVIS(I,1)
1677 SIGV(I,2) = SVIS(I,2)
1678 SIGV(I,3) = SVIS(I,3)
1679 SIGV(I,4) = SVIS(I,4)
1680 SIGV(I,5) = SVIS(I,5)
1681 SIGV(I,6) = SVIS(I,6)
1682 ENDDO
1683 ENDIF
1684 IF (MAT_ELEM%MAT_PARAM(MAT(1))%IVISC > 0) THEN
1685 DO I=1,NEL
1686 SVISC(I,1) = SVIS(I,1)
1687 SVISC(I,2) = SVIS(I,2)
1688 SVISC(I,3) = SVIS(I,3)
1689 SVISC(I,4) = SVIS(I,4)
1690 SVISC(I,5) = SVIS(I,5)
1691 SVISC(I,6) = SVIS(I,6)
1692 ENDDO
1693 ENDIF
1694C------------------------------------------------------------
1695C Variable to regularize with Non-local
1696C------------------------------------------------------------
1697 IF (ELBUF_TAB(NG)%BUFLY(ILAY)%L_PLA > 0) THEN
1698 IF (INLOC > 0) THEN
1699 DO I=LFT,LLT
1700 IF (OFF(I) == ONE) THEN
1701 VARNL(I) = DEFP(I)
1702 ELSE
1703 VARNL(I) = ZERO
1704 ENDIF
1705 ENDDO
1706 ENDIF
1707 ENDIF
1708C------------------------------------------------------------
1709C DEFINE SOUND SPEED (IN ALL CASE)
1710C DEFINE DYNAMIC VISCOSITY (FOR VISCOUS LAW)
1711C-----------------------
1712 DO I=LFT,LLT
1713 IF (SSP(I) == ZERO) SSP(I)=SQRT(C1(I)/RHO0(I))
1714 ENDDO
1715C-------------------------------------------
1716C BULK VISCOSITY AND TIME STEP COMPUTATION
1717C THIS SUBROUTINE RETURN THE NEW BULK VISCOSITY Q
1718C-----------------------
1719 IF (JSPH == 0)THEN
1720 JSMS = IPARG(52,NG)
1721 CALL MQVISCB(
1722 1 PM, OFF, RHO, BID1,
1723 2 EL_TEMP, SSP, BID3, STIFN,
1724 3 DT2T, NELTST, ITYPTST, AIRE,
1725 4 OFFG, GEO, PID, VOLN,
1726 5 VD2, DELTAX, VIS, D1,
1727 6 D2, D3, PNEW, PSH,
1728 7 MAT, NGL, Q, SSP_EQ,
1729 8 VOL, MSSA, DMELS, IGEO,
1730 9 FACQ0, CONDE, GBUF%DT, GBUF%G_DT,
1731 A IPM, RHOREF, RHOSP, NEL,
1732 B ITY, ISMSTR, JTUR, JTHE,
1733 C JSMS, NPG , glob_therm)
1734!------------------------------------------------
1735 ELSE
1736 CALL MDTSPH(
1737 1 PM, OFF, RHO, BID1,
1738 2 EL_TEMP, BID3, STIFN, DT2T,
1739 3 NELTST, ITYPTST, OFFG, GEO,
1740 4 PID, MUMAX, SSP, VOLN,
1741 5 VD2, DELTAX, VIS, D1,
1742 6 D2, D3, PNEW, PSH,
1743 7 MAT, NGL, Q, SSP_EQ,
1744 8 GBUF%G_DT,GBUF%DT, NEL, ITY,
1745 9 JTUR, JTHE)
1746 ENDIF
1747C-----------
1748 DO I=LFT,LLT
1749 P2 = -(SOLD1(I)+SIG(I,1)+SOLD2(I)+SIG(I,2)+SOLD3(I)+SIG(I,3))* THIRD
1750 E1 = D1(I)*(SOLD1(I)+SIG(I,1) + P2+TWO*SVIS(I,1))
1751 E2 = D2(I)*(SOLD2(I)+SIG(I,2) + P2+TWO*SVIS(I,2))
1752 E3 = D3(I)*(SOLD3(I)+SIG(I,3) + P2+TWO*SVIS(I,3))
1753 E4 = D4(I)*(SOLD4(I)+SIG(I,4) + TWO*SVIS(I,4))
1754 E5 = D5(I)*(SOLD5(I)+SIG(I,5) + TWO*SVIS(I,5))
1755 E6 = D6(I)*(SOLD6(I)+SIG(I,6) + TWO*SVIS(I,6))
1756 EINC = VOL_AVG(I)*DT1 * (E1+E2+E3+E4+E5+E6+E7(I))
1757 . - DVOL(I)*(Q(I) + QOLD(I) + P2)
1758 EINT(I) = EINT(I) + OFF(I)*EINC*HALF
1759 END DO
1760 QOLD(LFT:LLT) = Q(LFT:LLT)
1761
1762C------------------------------------------
1763 RETURN
1764 END SUBROUTINE USERMAT_SOLID
my_real
#define my_real
Definition cppsort.cpp:32
if
if(complex_arithmetic) id
fail_inievo_s
subroutine fail_inievo_s(nel, nuparam, nuvar, table, ntablf, itablf, time, uparam, ngl, aldt, dpla, epsp, uvar, signxx, signyy, signzz, signxy, signyz, signzx, pla, temp, sigy, off, dfmax, tdele, dmg_scale, uelr, ipg, npg, loff, damini, vol, inloc)
Definition fail_inievo_s.F:45
fail_johnson
subroutine fail_johnson(nel, nuparam, nuvar, time, timestep, uparam, ngl, signxx, signyy, signzz, signxy, signyz, signzx, dpla, epsp, tstar, uvar, off, dfmax, tdele)
Definition fail_johnson.F:38
fail_orthbiquad_s
subroutine fail_orthbiquad_s(nel, nuparam, nuvar, mfunc, kfunc, npf, tf, time, timestep, uparam, ngl, dpla, epsp, uvar, off, signxx, signyy, signzz, signxy, signyz, signzx, dfmax, tdele, aldt)
Definition fail_orthbiquad_s.F:40
fail_puck_s
subroutine fail_puck_s(nel, nuvar, ilay, npt0, time, timestep, uparam, ngl, off, noff, signxx, signyy, signzz, signxy, signyz, signzx, uvar, nuparam, dfmax, lf_dammx, tdele)
Definition fail_puck_s.F:37
fail_tab_s
subroutine fail_tab_s(nel, nuvar, npf, tf, time, uparam, ngl, aldt, signxx, signyy, signzz, signxy, signyz, signzx, dpla, epsp, tstar, uvar, ntablf, itablf, off, table, dfmax, tdele, nfunc, ifunc)
Definition fail_tab_s.F:44
fail_tbutcher_s
subroutine fail_tbutcher_s(nel, nuparam, nuvar, time, timestep, uparam, ngl, signxx, signyy, signzz, signxy, signyz, signzx, uvar, off, dfmax, tdele)
Definition fail_tbutcher_s.F:37
fail_tsaiwu_s
subroutine fail_tsaiwu_s(nel, nuvar, ip, ilay, npg, time, timestep, uparam, ngl, off, loff, noff, signxx, signyy, signzz, signxy, signyz, signzx, uvar, nuparam, dfmax, lf_dammx, tdele, dmg_scale)
Definition fail_tsaiwu_s.F:35
fail_visual_s
subroutine fail_visual_s(nel, nuparam, nuvar, time, timestep, uparam, epsxx, epsyy, epszz, epsxy, epsyz, epszx, signxx, signyy, signzz, signxy, signyz, signzx, uvar, off, ngl, dfmax, ismstr)
Definition fail_visual_s.F:36
fail_wilkins_s
subroutine fail_wilkins_s(nel, nuparam, nuvar, time, timestep, uparam, ngl, signxx, signyy, signzz, signxy, signyz, signzx, dpla, uvar, off, dfmax, tdele)
Definition fail_wilkins_s.F:37
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
mreploc
subroutine mreploc(ang, r11, r12, r13, r21, r22, r23, r31, r32, r33, rx, ry, rz, sx, sy, sz, tx, ty, tz, nel, jsph)
Definition mreploc.F:39
mrotensns
subroutine mrotensns(lft, llt, dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, r11, r12, r13, r21, r22, r23, r31, r32, r33)
Definition mrotens.F:96
mrotens
subroutine mrotens(lft, llt, e1, e2, e3, e4, dyz, e6, r11, r12, r13, r21, r22, r23, r31, r32, r33)
Definition mrotens.F:36
dt_mod
Definition dt_mod.F:79
law_userso
Definition law_userso.F:29
message_mod
Definition message_mod.F:1249
names_and_titles_mod
Definition names_and_titles_mod.F:997
names_and_titles_mod::ncharline100
integer, parameter ncharline100
Definition names_and_titles_mod.F:1002
nlocal_reg_mod
Definition nlocal_reg_mod.F:117
table_mod
Definition table_mod.F:112
nolib_usermat99
subroutine nolib_usermat99(ilaw_user, iuser_key)
Definition nolib_usermat99.F:33
epsf2u
subroutine epsf2u(nel, fxx, fxy, fxz, fyx, fyy, fyz, fzx, fzy, fzz, uxx, uyy, uzz, uxy, uyz, uxz)
Definition sigeps33.F:1448
sreploc3
subroutine sreploc3(rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition sreploc3.F:32
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
arret
subroutine arret(nn)
Definition arret.F:87
bidon2
subroutine bidon2(fac_l)
Definition machine.F:58
bidon1
subroutine bidon1(ibid)
Definition machine.F:48
startime
subroutine startime(event, itask)
Definition timer.F:93
tabulated
subroutine tabulated(iflag, nel, pm, off, eint, mu, espe, dvol, df, vnew, mat, psh, pnew, dpdm, dpde, npf, tf)
Definition tabulated.F:32
usermat_solid
subroutine usermat_solid(timers, lft, llt, nft, mtn, jcvt, pm, off, sig, eint, rho, qold, vol, strain, sigl, gama, uvar, bufmat, tf, npf, stifn, mat, ngl, nuvar, dt2t, neltst, ityptst, offg, geo, pid, epsd, el_temp, wxx, wyy, wzz, jsph, mumax, ssp, aire, voln, vd2, deltax, vis, d1, d2, d3, d4, d5, d6, pnew, psh, q, ssp_eq, dvol, sold1, sold2, sold3, sold4, sold5, sold6, rx, ry, rz, sx, sy, sz, tx, ty, tz, ipla, sigy, defp, ismstr, mfxx, mfxy, mfxz, mfyx, mfyy, mfyz, mfzx, mfzy, mfzz, ipm, isorth, fbuf, nfail, npg, sigdd, dxy, dyx, dyz, dzy, dzx, dxz, fr_wav, isrot, v, varnl, w, ix, x, jthe, et, mssa, dmels, iptr, ipts, iptt, table, fvd2, fdeltax, fssp, fqvis, iparg1, igeo, sigv, al_imp, signor, istrain, ng, elbuf_tab, vbuf, ilay, vk, iparg, bufvois, vdx, vdy, vdz, ihet, conde, itask, iexpan, vol_avg, amu, epsth3, epsth, svisc, nel, etotsh, iselect, tstar, muold, amu2, dpdm, rhoref, rhosp, nloc_dmg, ity, jtur, mat_elem, idel7nok, svis, dt, glob_therm, damp_buf, idamp_freq_range)
Definition usermat_solid.F:131