OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
fail_fld_xfem.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!|| fail_fld_xfem ../engine/source/materials/fail/fld/fail_fld_xfem.F
25!||--- called by ------------------------------------------------------
26!|| mulawc ../engine/source/materials/mat_share/mulawc.F90
27!|| usermat_shell ../engine/source/materials/mat_share/usermat_shell.F
28!||--- calls -----------------------------------------------------
29!|| finter ../engine/source/tools/curve/finter.F
30!|| finterfld ../engine/source/materials/fail/fld/fail_fld_c.F
31!||====================================================================
32 SUBROUTINE fail_fld_xfem(
33 1 NEL ,NUPARAM ,NUVAR ,NFUNC ,IFUNC ,
34 2 NPF ,TF ,TIME ,UPARAM ,
35 3 NGL ,IPT ,NPTOT ,SSP ,TENS ,
36 4 SIGNXX ,SIGNYY ,SIGNXY ,SIGNYZ ,SIGNZX ,
37 5 EPSXX ,EPSYY ,EPSXY ,EPSYZ ,EPSZX ,
38 6 UVAR ,NOFF ,OFF ,LF_DAMMX ,
39 7 ELCRKINI ,IXFEM ,IXEL ,ILAY ,IPTT ,
40 8 OFFL ,NPTT ,NPTTF ,DFMAX ,TDEL ,
41 9 DAM ,FLD_IDX ,NIPARAM ,IPARAM ,PLA ,
42 A DEPSXX ,DEPSYY ,DEPSXY ,DT1 )
43C-----------------------------------------------
44c FLD failure model for XFEM
45C-----------------------------------------------
46C I m p l i c i t T y p e s
47C-----------------------------------------------
48#include "implicit_f.inc"
49C-----------------------------------------------
50C G l o b a l P a r a m e t e r s
51C-----------------------------------------------
52#include "units_c.inc"
53#include "comlock.inc"
54#include "com_xfem1.inc"
55C---------+---------+---+---+--------------------------------------------
56C VAR | SIZE |TYP| RW| DEFINITION
57C---------+---------+---+---+--------------------------------------------
58C NEL | 1 | I | R | SIZE OF THE ELEMENT GROUP NEL
59C NUPARAM | 1 | I | R | SIZE OF THE USER PARAMETER ARRAY
60C UPARAM | NUPARAM | F | R | USER MATERIAL PARAMETER ARRAY
61C NUVAR | 1 | I | R | NUMBER OF USER ELEMENT VARIABLES
62C UVAR |NEL*NUVAR| F |R/W| USER ELEMENT VARIABLE ARRAY
63C---------+---------+---+---+--------------------------------------------
64C NFUNC | 1 | I | R | NUMBER FUNCTION USED FOR THIS USER LAW not used
65C IFUNC | NFUNC | I | R | FUNCTION INDEX
66C NPF | * | I | R | FUNCTION ARRAY
67C TF | * | F | R | FUNCTION ARRAY
68C---------+---------+---+---+--------------------------------------------
69C TIME | 1 | F | R | CURRENT TIME
70C---------+---------+---+---+--------------------------------------------
71C EPSXX | NEL | F | R | STRAIN XX
72C EPSYY | NEL | F | R | STRAIN YY
73C ... | | | |
74C---------+---------+---+---+--------------------------------------------
75C OFF | NEL | F |R/W| DELETED ELEMENT FLAG (=1. ON, =0. OFF)
76C---------+---------+---+---+--------------------------------------------
77C ILAY CURRENT LAYER
78C IPTT CURRENT INTEGRATION POINT IN THE LAYER (FOR OUTPUT ONLY)
79C---------+---------+---+---+--------------------------------------------
80C I N P U T A r g u m e n t s
81C-----------------------------------------------
82 INTEGER NEL,NUPARAM,NUVAR,NFUNC,IPT,IPTT,NPTOT,IXFEM,IXEL,ILAY,NPTT
83 INTEGER ELCRKINI(NXLAYMAX,NEL),NOFF(NEL),IFUNC(NFUNC),NGL(NEL),FLD_IDX(NEL)
84 my_real
85 . TIME,UPARAM(NUPARAM),SIGNXX(NEL),SIGNYY(NEL),
86 . SIGNZZ(NEL),SIGNXY(NEL),SIGNYZ(NEL),SIGNZX(NEL),SSP(NEL),
87 . EPSXX(NEL),EPSYY(NEL),EPSXY(NEL),EPSYZ(NEL),EPSZX(NEL),
88 . TENS(NEL,5),DFMAX(NEL,LF_DAMMX),TDEL(NEL),NPTTF(NEL),DAM(NEL)
89 INTEGER, INTENT(IN) :: NIPARAM, LF_DAMMX
90 INTEGER, DIMENSION(NIPARAM), INTENT(IN) :: IPARAM
91 my_real, DIMENSION(NEL), INTENT(IN) :: PLA,DEPSXX,DEPSYY,DEPSXY
92 my_real, INTENT(IN) :: dt1
93C-----------------------------------------------
94C I N P U T O U T P U T A r g u m e n t s
95C-----------------------------------------------
96 my_real uvar(nel,nuvar),off(nel),offl(nel)
97C-----------------------------------------------
98C VARIABLES FOR FUNCTION INTERPOLATION
99C-----------------------------------------------
100 INTEGER NPF(*)
101 my_real finter , finterfld ,tf(*)
102 EXTERNAL finter
103C Y = FINTER(IFUNC(J),X,NPF,TF,DYDX)
104C Y : y = f(x)
105C X : x
106C DYDX : f'(x) = dy/dx
107C IFUNC(J): FUNCTION INDEX
108C J : FIRST(J=1), SECOND(J=2) .. FUNCTION USED FOR THIS LAW
109C NPF,TF : FUNCTION PARAMETER
110C-----------------------------------------------
111C L o c a l V a r i a b l e s
112C-----------------------------------------------
113 INTEGER I,II,J,K,L,IENG,IAD,NINDX,NINDXP,ISHELL,LAYXFEM,LENF
114 INTEGER INDX(NEL),INDXP(NEL),RFLAG(NEL),RFLAGP(NEL)
115 my_real
116 . RANI,R1,R2,DADV,SIGMAX,S1,S2,SS,Q,DYDX,E12
117 my_real, ALLOCATABLE, DIMENSION(:) :: XF
118 my_real
119 . TBA(NEL),TBK(NEL),SIGR(NEL),EMAJ(NEL),EMIN(NEL),
120 . EM(NEL),EM_DCRIT(NEL),EM_DADV(NEL),DEMAJ(NEL),DEMIN(NEL),
121 . BETA(NEL),ALPHA,FCUT
122 CHARACTER (LEN=3) :: XCHAR
123C=======================================================================
124 iad = (ipt-1)*nel
125 xchar = ' '
126 DO i=1,nel
127 rflag(i) = 0
128 rflagp(i)= 0
129 indxp(i) = 0
130 indx(i) = 0
131 ENDDO
132 nindx = 0
133 nindxp = 0
134c---
135 !=================================================================
136 ! - INITIALISATION OF COMPUTATION ON TIME STEP
137 !=================================================================
138 ! Recovering failure criterion parameters
139 ! -> Integer parameter
140 ishell = iparam(1)
141 ieng = iparam(3)
142 ! -> Real parameters
143 dadv = uparam(2)
144 rani = uparam(3)
145 fcut = uparam(5)
146 IF (uparam(6) > zero) THEN
147 alpha = uparam(6)
148 ELSE
149 alpha = two*pi*fcut*dt1/(one + two*pi*fcut*dt1)
150 ENDIF
151c
152 layxfem = ixfem
153 IF (ixfem == 1 .and. ishell == 1) ishell=2
154c
155 DO i=1,nel
156 em_dcrit(i)= zero
157 em_dadv(i) = zero
158 END DO
159c
160 !=================================================================
161 ! - MINOR AND MAJOR (TRUE) STRAIN DEFORMATION
162 !=================================================================
163 DO i = 1,nel
164 ! For classic linear formulation
165 e12= half*epsxy(i)
166 s1 = half*(epsxx(i) + epsyy(i))
167 s2 = half*(epsxx(i) - epsyy(i))
168 q = sqrt(s2**2 + e12**2)
169 emaj(i) = s1 + q
170 emin(i) = s1 - q
171 IF (emin(i) >= emaj(i)) THEN
172 ss = emin(i)
173 emin(i) = emaj(i)
174 emaj(i) = ss
175 ENDIF
176 ! For non-linear path formulation
177 e12 = half*depsxy(i)
178 s1 = half*(depsxx(i) + depsyy(i))
179 s2 = half*(depsxx(i) - depsyy(i))
180 q = sqrt(s2**2 + e12**2)
181 demaj(i) = s1 + q
182 demin(i) = s1 - q
183 ! Strain ratio filtering value for non-linear path formulation
184 demaj(i) = alpha*demaj(i) + (one - alpha)*uvar(i,2)
185 demin(i) = alpha*demin(i) + (one - alpha)*uvar(i,3)
186 beta(i) = demin(i)/sign(max(abs(demaj(i)),em20),demaj(i))
187 uvar(i,2) = demaj(i)
188 uvar(i,3) = demin(i)
189 IF (ieng == 2) THEN
190 dfmax(i,4) = beta(i)
191 ENDIF
192 ENDDO
193c
194 !=================================================================
195 ! FAILURE MAJOR STRAIN FROM INPUT CURVE AND DAMAGE RATIO
196 !=================================================================
197 ! -> Engineering strains input
198 IF (ieng == 1) THEN ! transform input fld curve to true strain
199 ii = npf(ifunc(1))
200 lenf = npf(ifunc(1)+ 1) - npf(ifunc(1))
201 ALLOCATE(xf(lenf))
202 DO i = 1,lenf
203 xf(i) = log(tf(ii + i-1) + one)
204 ENDDO
205c
206 DO i = 1,nel
207 em(i) = finterfld(emin(i),lenf,xf)
208 dam(i) = emaj(i) / em(i)
209 dfmax(i,2) = dam(i)
210 dfmax(i,1) = min(one, dam(i))
211 ENDDO
212 ! -> True strain input
213 ELSE
214 ! -> Classic formulation
215 IF (ieng == 0) THEN
216 DO i = 1,nel
217 em(i) = finter(ifunc(1),emin(i),npf,tf,dydx)
218 dam(i) = emaj(i) / em(i)
219 dfmax(i,2) = dam(i)
220 dfmax(i,1) = min(one, dam(i))
221 ENDDO
222 ! -> Non-linear path formulation
223 ELSEIF (ieng == 2) THEN
224 DO i = 1,nel
225 em(i) = finter(ifunc(1),emin(i),npf,tf,dydx)
226 dam(i) = max(pla(i) / em(i),dam(i))
227 dfmax(i,2) = dam(i)
228 dfmax(i,1) = min(one, dam(i))
229 ENDDO
230 ENDIF
231 ENDIF
232c
233 !=================================================================
234 ! FLD ZONE INDEX CALCULATION FOR ANIM OUTPUT
235 !=================================================================
236 r1 = two*em02
237 r2 = rani/(rani+one)
238 IF (ieng < 2) THEN
239 DO i = 1,nel
240 IF (emaj(i) >= em(i)) THEN
241 fld_idx(i) = 6 ! zone 6 = failure
242 ELSEIF (emaj(i) >= em(i) - em01) THEN
243 fld_idx(i) = 5 ! zone 5 = margin to fail
244 ELSEIF (emaj(i)**2 + emin(i)**2 < r1**2) THEN
245 fld_idx(i) = 1 ! zone 1 = radius 0.02
246 ELSEIF (emaj(i) >= abs(emin(i))) THEN
247 fld_idx(i) = 4 ! zone 4 = safe (45 deg line)
248 ELSEIF (emaj(i) >= r2*abs(emin(i))) THEN
249 fld_idx(i) = 3 ! zone 3 = angle atan(r/(1+r)) - compression
250 ELSE
251 fld_idx(i) = 2 ! zone 2 - high wrinkle tendency
252 ENDIF
253 dfmax(i,3) = fld_idx(i)
254 ENDDO
255 ELSE
256 DO i = 1,nel
257 IF (pla(i) >= em(i)) THEN
258 fld_idx(i) = max(6,fld_idx(i)) ! zone 6 = failure
259 ELSEIF (pla(i) >= em(i) - em01) THEN
260 fld_idx(i) = max(5,fld_idx(i)) ! zone 5 = margin to fail
261 ELSEIF (pla(i)**2 + beta(i)**2 < r1**2) THEN
262 fld_idx(i) = max(1,fld_idx(i)) ! zone 1 = radius 0.02
263 ELSEIF (pla(i) >= abs(beta(i))) THEN
264 fld_idx(i) = max(4,fld_idx(i)) ! zone 4 = safe (45 deg line)
265 ELSEIF (pla(i) >= r2*abs(beta(i))) THEN
266 fld_idx(i) = max(3,fld_idx(i)) ! zone 3 = angle atan(r/(1+r)) - compression
267 ELSE
268 fld_idx(i) = max(2,fld_idx(i)) ! zone 2 - high wrinkle tendency
269 ENDIF
270 dfmax(i,3) = fld_idx(i)
271 ENDDO
272 ENDIF
273c--------------------------------------------------------------------
274 DO i = 1,nel
275 tens(i,1) = signxx(i)
276 tens(i,2) = signyy(i)
277 tens(i,3) = signxy(i)
278 tens(i,4) = signyz(i)
279 tens(i,5) = signzx(i)
280 em_dcrit(i)= em(i)
281 IF (nfunc > 1 .and. ifunc(2) > 0) THEN
282 em_dadv(i) = dadv*finter(ifunc(2),emin(i),npf,tf,dydx)
283 ELSE
284 em_dadv(i) = dadv*em_dcrit(i)
285 ENDIF
286 END DO
287C
288 IF (ixel > 0) THEN ! testing phantom elements
289 IF (ixel == 1) THEN
290 xchar = '1st'
291 ELSEIF (ixel == 2) THEN
292 xchar = '2nd'
293 ELSEIF (ixel == 3) THEN
294 xchar = '3rd'
295 ENDIF
296 ELSE
297 xchar = 'standard'
298 ENDIF
299c--------------------------
300 SELECT CASE (layxfem)
301c---------------
302 CASE (1) ! multilayer XFEM
303c---------------
304 DO i = 1,nel
305 IF (ieng == 2) emaj(i) = pla(i)
306 IF (off(i) == one) THEN
307 IF (uvar(i,1) == zero) THEN
308 IF (ixel == 0) THEN
309 IF (elcrkini(ilay,i)==0 .and. emaj(i)>=em_dcrit(i)) THEN
310 nindxp = nindxp+1
311 indxp(nindxp) = i
312 uvar(i,1) = one
313 offl(i) = zero
314 noff(i) = noff(i) + 1
315 npttf(i) = npttf(i) + one
316 rflagp(i) = 1
317 IF (int(npttf(i)) == nptt) THEN
318 nindx = nindx+1
319 indx(nindx) = i
320 elcrkini(ilay,i) = -1 ! one layer failed (by initiation)
321 rflag(i) = 1
322 ENDIF
323 IF (noff(i) == nptot) THEN
324 off(i) = four_over_5
325 tdel(i)= time
326 ENDIF
327 ELSEIF (elcrkini(ilay,i)==2 .and. emaj(i)>=em_dadv(i)) THEN
328 nindxp = nindxp + 1
329 indxp(nindxp) = i
330 uvar(i,1) = one
331 offl(i) = zero
332 noff(i) = noff(i) + 1
333 npttf(i) = npttf(i)+ one
334 rflagp(i) = 1
335 IF (int(npttf(i)) == nptt) THEN
336 nindx = nindx+1
337 indx(nindx) = i
338 elcrkini(ilay,i) = 1 ! one layer failed (by advancing)
339 rflag(i) = -1
340 ENDIF
341 IF (noff(i) == nptot) THEN
342 off(i) = four_over_5
343 tdel(i)= time
344 ENDIF
345 ENDIF
346 ELSEIF (emaj(i) >= em_dcrit(i)) THEN ! IXEL > 0
347 nindxp = nindxp + 1
348 indxp(nindxp) = i
349 uvar(i,1) = one
350 offl(i) = zero
351 npttf(i) = npttf(i) + one
352 rflagp(i) = 1
353 IF (int(npttf(i)) == nptt) THEN
354 nindx = nindx+1
355 indx(nindx) = i
356 off(i) = four_over_5
357 tdel(i) = time
358 rflag(i) = 3
359 ENDIF
360 ENDIF ! IXEL == 0
361 ELSEIF (uvar(i,1) == one) THEN
362 signxx(i) = zero
363 signyy(i) = zero
364 signxy(i) = zero
365 signyz(i) = zero
366 signzx(i) = zero
367 ENDIF ! UVAR
368 ENDIF ! IF (OFF(I) == ONE)
369 ENDDO ! DO I=1,NEL
370 IF (nindxp > 0) THEN
371 DO j=1,nindxp
372 i = indxp(j)
373 signxx(i) = zero
374 signyy(i) = zero
375 signxy(i) = zero
376 signyz(i) = zero
377 signzx(i) = zero
378#include "lockon.inc"
379 IF (rflagp(i) == 1) WRITE(iout, 3800)ngl(i),ilay,iptt
380 IF (rflagp(i) == 1) WRITE(istdo,3900)ngl(i),ilay,iptt,time
381#include "lockoff.inc"
382 ENDDO
383 ENDIF ! NINDXP > 0
384 IF (nindx > 0) THEN
385 DO j=1,nindx
386 i = indx(j)
387#include "lockon.inc"
388c initialization std element
389 IF (rflag(i) == 1) WRITE(iout ,3000) ngl(i),ilay
390 IF (rflag(i) == 1) WRITE(istdo,3100) ngl(i),ilay,time
391c advancement std element
392 IF (rflag(i) == -1) WRITE(iout ,3200) ngl(i),ilay
393 IF (rflag(i) == -1) WRITE(istdo,3300) ngl(i),ilay,time
394c delete phantom element
395 IF (rflag(i) == 3) WRITE(iout, 3400) xchar,ngl(i),ilay
396 IF (rflag(i) == 3) WRITE(istdo,3500) xchar,ngl(i),ilay,time
397#include "lockoff.inc"
398 ENDDO
399 ENDIF ! NINDX > 0
400c---------------
401 CASE (2) ! monolayer XFEM
402c---------------
403 IF (ishell == 1) THEN
404c
405 IF (layxfem > 0) THEN
406 DO i = 1,nel
407 IF (ieng == 2) emaj(i) = pla(i)
408 IF (off(i) == one) THEN
409 IF (ixel == 0) THEN ! not cut yet
410 IF (elcrkini(ilay,i) == 0 .and.
411 . emaj(i) >= em_dcrit(i)) THEN
412 elcrkini(ilay,i) = -1 ! ready to start
413 off(i) = four_over_5
414 nindx=nindx+1
415 indx(nindx)=i
416 rflag(i) = 1
417 tdel(i)= time
418 ELSEIF (elcrkini(ilay,i) == 2 .and.
419 . emaj(i) >= em_dadv(i)) THEN
420 elcrkini(ilay,i) = 1 ! ready to advance
421 off(i) = four_over_5
422 nindx=nindx+1
423 indx(nindx)=i
424 rflag(i) = -1
425 tdel(i)= time
426 ENDIF
427 ELSE IF (emaj(i) >= em_dcrit(i)) then! IXEL > 0
428 off(i) = four_over_5
429 nindx=nindx+1
430 indx(nindx)=i
431 rflag(i) = 3
432 ENDIF ! IXEL
433 ENDIF
434 ENDDO
435 ENDIF ! LAYXFEM
436C
437 IF (nindx > 0) THEN
438 DO j=1,nindx
439 i=indx(j)
440#include "lockon.inc"
441 WRITE(iout, 4800)ngl(i),iptt
442 WRITE(istdo,4900)ngl(i),iptt,time
443c initialization std element
444 IF (rflag(i) == 1) WRITE(iout ,4000) ngl(i)
445 IF (rflag(i) == 1) WRITE(istdo,4100) ngl(i),time
446c advancement std element
447 IF (rflag(i) == -1) WRITE(iout, 4200) ngl(i)
448 IF (rflag(i) == -1) WRITE(istdo,4300) ngl(i),time
449c delete phantom
450 IF (rflag(i) == 3) WRITE(iout, 4400) xchar,ngl(i)
451 IF (rflag(i) == 3) WRITE(istdo,4500) xchar,ngl(i),time
452#include "lockoff.inc"
453 ENDDO
454 ENDIF
455c
456 ELSEIF (ishell == 2)THEN
457c
458 DO i = 1,nel
459 IF (ieng == 2) emaj(i) = pla(i)
460 IF (off(i) == one) THEN
461 IF (uvar(i,1) == zero) THEN
462 IF (ixel == 0) THEN ! testing standard element
463 IF (elcrkini(ilay,i) == 0 .and. emaj(i) >= em_dcrit(i)) THEN
464 uvar(i,1) = one
465 nindx=nindx+1
466 indx(nindx)=i
467 nindxp = nindxp + 1
468 indxp(nindxp) = i
469 rflagp(i) = 1
470 noff(i) = noff(i) + 1
471 IF (noff(i) == nptot) THEN
472 elcrkini(ilay,i) = -1 ! one layer failed (by initiation)
473 off(i) = four_over_5
474 rflag(i) = 1
475 tdel(i)= time
476 ENDIF
477 ELSEIF (elcrkini(ilay,i) == 2 .and. emaj(i) >= em_dadv(i)) THEN
478 uvar(i,1) = one
479 nindx=nindx+1
480 indx(nindx)=i
481 nindxp = nindxp + 1
482 indxp(nindxp) = i
483 rflagp(i) = 1
484 noff(i) = noff(i) + 1
485 IF (noff(i) == nptot) THEN
486 off(i) = four_over_5
487 elcrkini(ilay,i) = 1 ! one layer failed (by advancing)
488 rflag(i) = -1
489 tdel(i)= time
490 ENDIF
491 ENDIF
492 ELSEIF (emaj(i) >= em_dcrit(i)) THEN ! IXEL > 0
493 uvar(i,1) = one
494 nindx = nindx+1
495 indx(nindx) = i
496 nindxp = nindxp + 1
497 indxp(nindxp) = i
498 rflagp(i) = 1
499 noff(i) = noff(i) + 1
500 IF (noff(i) == nptot) THEN
501 off(i) = four_over_5
502 rflag(i) = 3
503 ENDIF
504 ENDIF ! IXEL
505 ELSEIF (uvar(i,1) == one) THEN
506 signxx(i) = zero
507 signyy(i) = zero
508 signxy(i) = zero
509 signyz(i) = zero
510 signzx(i) = zero
511 ENDIF ! UVAR
512 ENDIF
513 ENDDO
514c
515 IF (nindxp > 0) THEN
516 DO j=1,nindxp
517 i = indxp(j)
518 signxx(i) = zero
519 signyy(i) = zero
520 signxy(i) = zero
521 signyz(i) = zero
522 signzx(i) = zero
523#include "lockon.inc"
524 IF (rflagp(i) == 1) WRITE(iout, 4800)ngl(i),iptt
525 IF (rflagp(i) == 1) WRITE(istdo,4900)ngl(i),iptt,time
526#include "lockoff.inc"
527 ENDDO
528 ENDIF ! NINDXP > 0
529 IF (nindx > 0) THEN
530 DO j=1,nindx
531 i = indx(j)
532#include "lockon.inc"
533c initialization std element
534 IF (rflag(i) == 1) WRITE(iout ,4000) ngl(i)
535 IF (rflag(i) == 1) WRITE(istdo,4100) ngl(i),time
536c advancement std element
537 IF (rflag(i) == -1) WRITE(iout, 4200) ngl(i)
538 IF (rflag(i) == -1) WRITE(istdo,4300) ngl(i),time
539c delete phantom
540 IF (rflag(i) == 3) WRITE(iout, 4400) xchar,ngl(i)
541 IF (rflag(i) == 3) WRITE(istdo,4500) xchar,ngl(i),time
542#include "lockoff.inc"
543 ENDDO
544 ENDIF ! nindx > 0
545 ENDIF ! ISHELL
546c-----------------
547 END SELECT ! LAYXFEM
548c
549 IF (ieng == 1) DEALLOCATE(xf)
550c-----------------------------------------------
551 2000 FORMAT(1x,'FOR SHELL ELEMENT (FLD)',i10,1x,'LAYER',i3,':',/
552 . 1x, 'STRESS TENSOR SET TO ZERO')
553 2100 FORMAT(1x,'FOR SHELL ELEMENT (FLD)',i10,1x,'LAYER',i3,':',/,
554 . 1x, 'STRESS TENSOR SET TO ZERO',1x,'AT TIME :',1pe20.13)
555 2400 FORMAT(1x,1pg20.13,' % OF THICKNESS OF SHELL BROKEN ')
556 2410 FORMAT(1x,1pg20.13,' % OF THICKNESS OF SHELL ',i10,' BROKEN ')
557 2500 FORMAT(1x,' LOWER SKIN -> UPPER SKIN ')
558 2600 FORMAT(1x,' UPPER SKIN -> LOWER SKIN ')
559c--- multilayer xfem
560 3000 FORMAT(1x,'CRACK INITIALIZATION IN SHELL ELEMENT',i10,1x,'LAYER',i3)
561 3100 FORMAT(1x,'CRACK INITIALIZATION IN SHELL ELEMENT',i10,1x,'LAYER',i3,/
562 . 1x,'AT TIME :',1pe12.4)
563 3200 FORMAT(1x,'CRACK ADVANCEMENT IN SHELL ELEMENT',i10,' LAYER',i3)
564 3300 FORMAT(1x,'CRACK ADVANCEMENT IN SHELL ELEMENT',i10,' LAYER',i3/
565 . 1x,'AT TIME :',1pe12.4)
566 3400 FORMAT(1x,'DELETE ',a4,' PHANTOM ELEMENT, SHELL ID=',i10,' LAYER',i3)
567 3500 FORMAT(1x,'DELETE ',a4,' PHANTOM ELEMENT, SHELL ID=',i10,' LAYER',i3,/
568 . 1x,'AT TIME :',1pe12.4)
569 3800 FORMAT(1x,'FLD FAILURE IN SHELL',i10,1x,'LAYER',i3,1x,'INT POINT',i2)
570 3900 FORMAT(1x,'FLD FAILURE IN SHELL',i10,1x,'LAYER',i3,1x,'INT POINT',i2,/
571 . 1x,'AT TIME :',1pe12.4)
572c--- monolayer xfem
573 4000 FORMAT(1x,'CRACK INITIALIZATION IN SHELL ELEMENT',i10)
574 4100 FORMAT(1x,'CRACK INITIALIZATION IN SHELL ELEMENT',i10,/
575 . 1x,'AT TIME :',1pe12.4)
576 4200 FORMAT(1x,'CRACK ADVANCEMENT IN SHELL ELEMENT',i10)
577 4300 FORMAT(1x,'CRACK ADVANCEMENT IN SHELL ELEMENT',i10,/
578 . 1x,'AT TIME :',1pe12.4)
579 4400 FORMAT(1x,'DELETE ',a4,' PHANTOM ELEMENT, SHELL ID=',i10)
580 4500 FORMAT(1x,'DELETE ',a4,' PHANTOM ELEMENT, SHELL ID=',i10,/
581 . 1x,'AT TIME :',1pe12.4)
582 4800 FORMAT(1x,'FLD FAILURE IN SHELL',i10,1x,'INT POINT',i2)
583 4900 FORMAT(1x,'FLD FAILURE IN SHELL',i10,1x,'INT POINT',i2,1x,'AT TIME :',1pe12.4)
584c-----------------------------------------------
585 RETURN
586 END
fail_fld_xfem
subroutine fail_fld_xfem(nel, nuparam, nuvar, nfunc, ifunc, npf, tf, time, uparam, ngl, ipt, nptot, ssp, tens, signxx, signyy, signxy, signyz, signzx, epsxx, epsyy, epsxy, epsyz, epszx, uvar, noff, off, lf_dammx, elcrkini, ixfem, ixel, ilay, iptt, offl, nptt, npttf, dfmax, tdel, dam, fld_idx, niparam, iparam, pla, depsxx, depsyy, depsxy, dt1)
Definition fail_fld_xfem.F:43
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21