OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
upxfem1.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!|| upxfem1 ../engine/source/elements/xfem/upxfem1.F
25!||--- called by ------------------------------------------------------
26!|| resol ../engine/source/engine/resol.F
27!||--- calls -----------------------------------------------------
28!|| startimeg ../engine/source/system/timer.F
29!|| stoptimeg ../engine/source/system/timer.F
30!|| upenric1_n3 ../engine/source/elements/xfem/upenric1_n3.F
31!|| upenric1_n4 ../engine/source/elements/xfem/upenric1_n4.F
32!|| upenric2_n3 ../engine/source/elements/xfem/upenric2_n3.F
33!|| upenric2_n4 ../engine/source/elements/xfem/upenric2_n4.F
34!||--- uses -----------------------------------------------------
35!|| crackxfem_mod ../engine/share/modules/crackxfem_mod.F
36!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
37!|| element_mod ../common_source/modules/elements/element_mod.F90
38!||====================================================================
39 SUBROUTINE upxfem1(XFEM_TAB,
40 . IPARG ,IXC ,NGROUC ,IGROUC ,IXTG ,
41 . IADC_CRK ,IEL_CRK ,INOD_CRK,ELCUTC ,NODEDGE ,
42 . ENRTAG ,CRKEDGE ,XEDGE4N ,XEDGE3N )
43C-----------------------------------------------
44C M o d u l e s
45C-----------------------------------------------
46 USE crackxfem_mod
47 USE elbufdef_mod
48 use element_mod , only : nixc,nixtg
49C-----------------------------------------------
50C I m p l i c i t T y p e s
51C-----------------------------------------------
52#include "implicit_f.inc"
53C-----------------------------------------------
54C C o m m o n B l o c k s
55C-----------------------------------------------
56#include "param_c.inc"
57#include "com01_c.inc"
58#include "com04_c.inc"
59#include "com_xfem1.inc"
60#include "task_c.inc"
61#include "vect01_c.inc"
62C-----------------------------------------------
63C D u m m y A r g u m e n t s
64C-----------------------------------------------
65 INTEGER IPARG(NPARG,*),IXC(NIXC,*),NGROUC,IGROUC(*),IEL_CRK(*),
66 . ELCUTC(2,*),IADC_CRK(*),INOD_CRK(*),IXTG(NIXTG,*),
67 . NODEDGE(2,*),ENRTAG(NUMNOD,*),XEDGE4N(4,*),XEDGE3N(3,*)
68C
69 TYPE (ELBUF_STRUCT_), TARGET ,DIMENSION(NGROUP,NXEL) :: XFEM_TAB
70 TYPE (XFEM_EDGE_) , DIMENSION(*) :: CRKEDGE
71C-----------------------------------------------
72C L o c a l V a r i a b l e s
73C-----------------------------------------------
74 INTEGER IG,NG,JFT,JLT,NEL,ITG1,ITG2,NXLAY,OFF,IXFEM
75C=======================================================================
76C check cracks connections
77c---------------------------
78C SMP dynamic parallel loop
79!$OMP DO SCHEDULE(DYNAMIC,1)
80c---------------------------
81 DO ig = 1, ngrouc
82 ng = igrouc(ig)
83 off = iparg(8,ng)
84 ixfem = iparg(54,ng) * iparg(70,ng) ! IXFEM * ACTIVEXFEM
85 IF (off == 0 .and. ixfem > 0) THEN
86 IF (iddw > 0) CALL startimeg(ng)
87C---
88 ity = iparg(5,ng)
89 nel = iparg(2,ng)
90 nft = iparg(3,ng)
91 nxlay = iparg(59,ng)
92 lft = 1
93 llt = min(nvsiz,nel)
94 jft=lft
95 jlt=llt
96c
97C--- crack connection on an edge (enrs) : iecut=3
98 IF (ity == 3) THEN
99 CALL upenric1_n4(iparg ,ixc ,nft ,jft ,jlt ,
100 . elcutc ,iadc_crk,iel_crk,inod_crk ,nxlay ,
101 . nodedge,enrtag ,crkedge ,xedge4n )
102 ELSEIF (ity == 7) THEN
103 itg1 = 1+numelc
104 itg2 = 1+4*ecrkxfec
105 CALL upenric1_n3(
106 . iparg ,ixtg ,nft ,jft ,jlt ,
107 . elcutc(1,itg1),iadc_crk(itg2),iel_crk(itg1),inod_crk,nxlay,
108 . nodedge ,enrtag,crkedge ,xedge3n )
109 ENDIF
110C---
111 IF (iddw > 0) CALL stoptimeg(ng)
112 ENDIF
113 ENDDO
114!$OMP END DO
115C-----------------------------------------------
116C
117C SMP dynamic parallel loop
118C
119!$OMP DO SCHEDULE(DYNAMIC,1)
120 DO ig = 1, ngrouc
121 ng = igrouc(ig)
122 off = iparg(8,ng)
123 ixfem = iparg(54,ng) * iparg(70,ng) ! IXFEM * ACTIVEXFEM
124 IF (off == 0 .and. ixfem > 0) THEN
125 IF (iddw > 0) CALL startimeg(ng)
126C---
127 ity = iparg(5,ng)
128 nel = iparg(2,ng)
129 nft = iparg(3,ng)
130 nxlay = iparg(59,ng)
131 lft = 1
132 llt = min(nvsiz,nel)
133 jft=lft
134 jlt=llt
135C---
136 IF (ity == 3) THEN
137 CALL upenric2_n4(iparg ,ixc ,nft ,jft ,jlt ,
138 . elcutc ,iadc_crk,iel_crk,inod_crk,enrtag,
139 . nxlay ,crkedge ,xedge4n )
140 ELSEIF (ity == 7) THEN
141 itg1 = 1+numelc
142 itg2 = 1+4*ecrkxfec
143 CALL upenric2_n3(
144 . iparg ,ixtg ,nft ,jft ,jlt ,
145 . elcutc(1,itg1),iadc_crk(itg2),iel_crk(itg1),inod_crk,enrtag,
146 . nxlay ,crkedge ,xedge3n )
147 ENDIF
148C---
149 IF (iddw > 0) CALL stoptimeg(ng)
150 ENDIF
151 ENDDO
152!$OMP END DO
153c-------------
154 RETURN
155 END
startimeg
subroutine startimeg(ng)
Definition timer.F:1371
stoptimeg
subroutine stoptimeg(ng)
Definition timer.F:1419
min
#define min(a, b)
Definition macros.h:20
crackxfem_mod
Definition crackxfem_mod.F:121
upenric1_n3
subroutine upenric1_n3(iparg, ixtg, nft, jft, jlt, elcutc, iad_crktg, iel_crktg, inod_crk, nxlay, nodedge, enrtag, crkedge, xedge3n)
Definition upenric1_n3.F:35
upenric1_n4
subroutine upenric1_n4(iparg, ixc, nft, jft, jlt, elcutc, iadc_crk, iel_crk, inod_crk, nxlay, nodedge, enrtag, crkedge, xedge4n)
Definition upenric1_n4.F:34
upenric2_n3
subroutine upenric2_n3(iparg, ixtg, nft, jft, jlt, elcutc, iad_crktg, iel_crktg, inod_crk, enrtag, nxlay, crkedge, xedge3n)
Definition upenric2_n3.F:34
upenric2_n4
subroutine upenric2_n4(iparg, ixc, nft, jft, jlt, elcutc, iadc_crk, iel_crk, inod_crk, enrtag, nxlay, crkedge, xedge4n)
Definition upenric2_n4.F:34
upxfem1
subroutine upxfem1(xfem_tab, iparg, ixc, ngrouc, igrouc, ixtg, iadc_crk, iel_crk, inod_crk, elcutc, nodedge, enrtag, crkedge, xedge4n, xedge3n)
Definition upxfem1.F:43