OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
check_active_elem_edge.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!|| check_active_elem_edge ../engine/source/interfaces/interf/check_active_elem_edge.F
25!||--- called by ------------------------------------------------------
26!|| check_edge_state ../engine/source/interfaces/interf/check_edge_state.F
27!|| check_surface_state ../engine/source/interfaces/interf/check_surface_state.F
28!||====================================================================
29 SUBROUTINE check_active_elem_edge( NUMBER_NODE, N1,N2,N3,N4,
30 . DEACTIVATION,GEO,IXS,IXC,
31 . IXT,IXP,IXR,IXTG,IXS10,ADDCNEL,CNEL,
32 . TAG_NODE,TAG_ELEM )
33!$COMMENT
34! CHECK_ACTIVE_ELEM_EDGE description
35! check if a element associated to an edge is active
36! CHECK_ACTIVE_ELEM_EDGE organization
37! loop over the element associated to the N1 node
38! if 1 active element contains N1 & N2, the edge must be kept (for interface 11)
39! if 1 active element contains N1 & N2 & N3 & N4, the surface must be kept (for interface 7)
40!$ENDCOMMENT
41C-----------------------------------------------
42C I m p l i c i t T y p e s
43C-----------------------------------------------
44#include "implicit_f.inc"
45C-----------------------------------------------
46C C o m m o n B l o c k s
47C-----------------------------------------------
48#include "com04_c.inc"
49#include "param_c.inc"
50#include "scr17_c.inc"
51C-----------------------------------------------
52C D u m m y A r g u m e n t s
53C-----------------------------------------------
54 LOGICAL, INTENT(inout) :: DEACTIVATION
55 INTEGER, INTENT(in) :: NUMBER_NODE ! number of node of the edge/surface
56 INTEGER, INTENT(in) :: N1,N2,N3,N4 ! node id : for interface type11 -> N3&N4 are fake node id (ie. =0)
57 INTEGER, DIMENSION(NIXS,NUMELS), INTENT(in) :: IXS ! solid array
58 INTEGER, DIMENSION(NIXC,NUMELC), INTENT(in) :: IXC ! shell array
59 INTEGER, DIMENSION(NIXT,NUMELT), INTENT(in) :: IXT! truss array
60 INTEGER, DIMENSION(NIXP,NUMELP), INTENT(in) :: IXP! beam array
61 INTEGER, DIMENSION(NIXR,NUMELR), INTENT(in) :: IXR! spring array
62 INTEGER, DIMENSION(NIXTG,NUMELTG), INTENT(in) :: IXTG! triangle array
63 INTEGER, DIMENSION(6,NUMELS10), INTENT(in) :: IXS10!< tetra10 data
64 INTEGER, DIMENSION(0:NUMNOD+1), INTENT(in) :: ADDCNEL ! address for the CNEL array
65 my_real, DIMENSION(NPROPG,NUMGEO), INTENT(in) :: geo
66 INTEGER, DIMENSION(0:LCNEL), INTENT(in) :: CNEL ! connectivity node-->element
67 INTEGER, DIMENSION(NUMNOD), INTENT(inout) :: TAG_NODE
68 INTEGER, DIMENSION(NUMELS+NUMELQ+NUMELC+NUMELT+NUMELP+NUMELR+NUMELTG), INTENT(inout) :: TAG_ELEM
69C-----------------------------------------------
70C L o c a l V a r i a b l e s
71C-----------------------------------------------
72 LOGICAL :: STILL_COMPUTE,STILL_ALIVE
73 INTEGER :: I,K
74 INTEGER :: NEXT
75 INTEGER :: ELEM_ID,NODE_ID,OTHER_NODE
76 INTEGER :: OFFSET_SOLID,OFFSET_QUAD,OFFSET_SHELL
77 INTEGER :: OFFSET_TRUSS,OFFSET_BEAM,OFFSET_SPRING
78 INTEGER :: OFFSET_TRIANGLE,OFFSET_UR
79C-----------------------------------------------
80 ! --------------------------
81 offset_solid = 0
82 offset_quad=offset_solid+numels
83 offset_shell=offset_quad+numelq
84 offset_truss=offset_shell+numelc
85 offset_beam=offset_truss+numelt
86 offset_spring=offset_beam+numelp
87 offset_triangle=offset_spring+numelr
88 offset_ur=offset_triangle+numeltg
89 ! --------------------------
90
91 still_compute = .true.
92 next = 0
93 i = addcnel(n1) + next
94 deactivation = .false.
95 still_alive = .false.
96
97 ! ------------------
98 ! check if 1 or more elements is/are associated to
99 ! N1 node
100 ! if no element, the edge can be deactivated (only if
101 ! the remote elements associated to N1/N2 are all deleted)
102 IF((addcnel(n1+1) - addcnel(n1)) ==0) THEN
103 deactivation = .true.
104 still_compute = .false.
105 ENDIF
106 ! ------------------
107
108 DO WHILE( still_compute )
109
110 elem_id = cnel(i)
111 tag_node(n1) = 0
112 tag_node(n2) = 0
113 IF(number_node>2) THEN
114 tag_node(n3) = 0
115 tag_node(n4) = 0
116 ENDIF
117 ! -------------------------
118 ! loop over the active element
119 ! for interface 11 : if still 1 or more active elements with N1 & N2 : need to keep the edge
120 ! else : need to deactivate the edge
121 ! for interface 7 : if still 1 or more active elements with N1 & N2 & N3 & N4 : need to keep the surface
122 ! else : need to deactivate the surface
123 IF(tag_elem(elem_id)>0) THEN
124 ! -----------------
125 ! solid element
126 IF(elem_id<=offset_shell) THEN
127 DO k=2,9
128 node_id = ixs(k,elem_id)
129 tag_node(node_id) = 1
130 ENDDO
131 IF(elem_id>numels8.AND.elem_id<=numels8+numels10) THEN
132 DO k=1,6
133 node_id = ixs10(k,elem_id-numels8)
134 tag_node(node_id) = 1
135 ENDDO
136 ENDIF
137 ELSEIF(elem_id>offset_shell.AND.elem_id<=offset_truss) THEN
138 ! -----------------
139 ! shell element
140 DO k=2,5
141 node_id = ixc(k,elem_id-offset_shell)
142 tag_node(node_id) = 1
143 ENDDO
144 ELSEIF(elem_id>offset_truss.AND.elem_id<=offset_beam) THEN
145 ! -----------------
146 ! truss element
147 DO k=2,3
148 node_id = ixt(k,elem_id-offset_truss)
149 tag_node(node_id) = 1
150 ENDDO
151 ELSEIF(elem_id>offset_beam.AND.elem_id<=offset_spring) THEN
152 ! -----------------
153 ! beam element
154 DO k=2,3
155 node_id = ixp(k,elem_id-offset_beam)
156 tag_node(node_id) = 1
157 ENDDO
158 ELSEIF(elem_id>offset_spring.AND.elem_id<=offset_triangle) THEN
159 ! -----------------
160 ! spring element
161 DO k=2,3
162 node_id = ixr(k,elem_id-offset_spring)
163 tag_node(node_id) = 1
164 ENDDO
165
166 IF(nint(geo(12,ixr(1,elem_id-offset_spring))) == 12) THEN
167 node_id = ixr(4,elem_id-offset_spring)
168 tag_node(node_id) = 1
169 ENDIF
170 ELSEIF(elem_id>offset_triangle.AND.elem_id<=offset_ur) THEN
171 ! -----------------
172 ! triangle element
173 DO k=2,4
174 node_id = ixtg(k,elem_id-offset_triangle)
175 tag_node(node_id) = 1
176 ENDDO
177 ENDIF
178 ! -----------------
179
180 ! -----------------
181 ! for interface 7 : need to check if the element has N3 & N4
182 other_node = 0
183 IF(number_node>2) other_node = tag_node(n3)
184 IF(number_node>3) other_node = other_node + tag_node(n4)
185
186 ! -----------------
187
188 ! -----------------
189 ! check if the element has N1 & N2
190 IF(tag_node(n1)+tag_node(n2)+other_node==number_node) THEN
191 still_compute = .false.
192 still_alive = .true.
193 ENDIF
194 ! -----------------
195
196 ENDIF
197 next = next + 1
198 i = addcnel(n1) + next
199 IF(i>addcnel(n1+1)-1) THEN
200 still_compute = .false.
201 ENDIF
202 ! -------------------------
203 ENDDO
204
205 IF(.NOT.still_alive) THEN
206 deactivation = .true.
207 ENDIF
208
209 RETURN
210 END SUBROUTINE check_active_elem_edge
211
check_active_elem_edge
subroutine check_active_elem_edge(number_node, n1, n2, n3, n4, deactivation, geo, ixs, ixc, ixt, ixp, ixr, ixtg, ixs10, addcnel, cnel, tag_node, tag_elem)
Definition check_active_elem_edge.F:33
my_real
#define my_real
Definition cppsort.cpp:32