OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
brezo2.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!|| brezo2 ../engine/source/ale/ale2d/brezo2.F
25!||--- called by ------------------------------------------------------
26!|| arezon ../engine/source/ale/arezon.F
27!||--- uses -----------------------------------------------------
28!|| ale_connectivity_mod ../common_source/modules/ale/ale_connectivity_mod.F
29!||====================================================================
30 SUBROUTINE brezo2(ALE_CONNECT,VAR,PHI,FLUX,VOL,BHOLE,NM)
31C-----------------------------------------------
32C M o d u l e s
33C-----------------------------------------------
34 USE ale_connectivity_mod
35C-----------------------------------------------
36C I m p l i c i t T y p e s
37C-----------------------------------------------
38#include "implicit_f.inc"
39C-----------------------------------------------
40C C o m m o n B l o c k s
41C-----------------------------------------------
42#include "com01_c.inc"
43#include "com08_c.inc"
44#include "vect01_c.inc"
45C-----------------------------------------------
46C D u m m y A r g u m e n t s
47C-----------------------------------------------
48 INTEGER BHOLE(NMULT,*),NM
49 my_real var(*), phi(*), flux(4,*), vol(*)
50 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
51C-----------------------------------------------
52C L o c a l V a r i a b l e s
53C-----------------------------------------------
54 INTEGER I, IE, IV1, IV2, IV3, IV4, IAD2
55 my_real voln
56C-----------------------------------------------
57C S o u r c e L i n e s
58C-----------------------------------------------
59 IF(jmult == 0)THEN
60 DO i=lft,llt
61 ie =nft+i
62 iad2 = ale_connect%ee_connect%iad_connect(ie)
63 iv1=ale_connect%ee_connect%connected(iad2 + 1 - 1)
64 iv2=ale_connect%ee_connect%connected(iad2 + 2 - 1)
65 iv3=ale_connect%ee_connect%connected(iad2 + 3 - 1)
66 iv4=ale_connect%ee_connect%connected(iad2 + 4 - 1)
67 IF(iv1 <= zero)THEN
68 iv1=ie
69 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv1)) /= abs(bhole(nm,ie)))THEN
70 iv1=ie
71 ENDIF
72 IF(iv2 <= 0)THEN
73 iv2=ie
74 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv2)) /= abs(bhole(nm,ie)))THEN
75 iv2=ie
76 ENDIF
77 IF(iv3 <= 0)THEN
78 iv3=ie
79 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv3)) /= abs(bhole(nm,ie)))THEN
80 iv3=ie
81 ENDIF
82 IF(iv4 <= 0)THEN
83 iv4=ie
84 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv4)) /= abs(bhole(nm,ie)))THEN
85 iv4=ie
86 ENDIF
87 var(i) = var(i) + half * dt1 * (phi(ie)*
88 2 (flux(1,i)+flux(2,i)+flux(3,i)+flux(4,i))
89 3 -phi(iv1)*flux(1,i)-phi(iv2)*flux(2,i)-phi(iv3)*flux(3,i)
90 4 -phi(iv4)*flux(4,i)) / vol(i)
91 ENDDO !nexi i
92 ELSE
93 DO i=lft,llt
94 ie =nft+i
95 iad2 = ale_connect%ee_connect%iad_connect(ie)
96 iv1=ale_connect%ee_connect%connected(iad2 + 1 - 1)
97 iv2=ale_connect%ee_connect%connected(iad2 + 2 - 1)
98 iv3=ale_connect%ee_connect%connected(iad2 + 3 - 1)
99 iv4=ale_connect%ee_connect%connected(iad2 + 4 - 1)
100 IF(iv1 <= zero)THEN
101 iv1=ie
102 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv1)) /= abs(bhole(nm,ie)))THEN
103 iv1=ie
104 ENDIF
105 IF(iv2 <= zero)THEN
106 iv2=ie
107 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv2)) /= abs(bhole(nm,ie)))THEN
108 iv2=ie
109 ENDIF
110 IF(iv3 <= zero)THEN
111 iv3=ie
112 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv3)) /= abs(bhole(nm,ie)))THEN
113 iv3=ie
114 ENDIF
115 IF(iv4 <= zero)THEN
116 iv4=ie
117 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv4)) /= abs(bhole(nm,ie)))THEN
118 iv4=ie
119 ENDIF
120 voln=vol(i)-dt1 * (flux(1,i)+flux(2,i)+flux(3,i)+flux(4,i))
121 IF (voln > em15) THEN
122 var(i) = var(i) + half * dt1 * (phi(ie)*
123 2 (flux(1,i)+flux(2,i)+flux(3,i)+flux(4,i))
124 3 -phi(iv1)*flux(1,i)-phi(iv2)*flux(2,i)-phi(iv3)*flux(3,i)
125 4 -phi(iv4)*flux(4,i)) / max(em15,voln)
126 ENDIF
127 ENDDO
128 ENDIF
129 RETURN
130 END
131!||====================================================================
132!|| brezo2_sig ../engine/source/ale/ale2d/brezo2.F
133!||--- uses -----------------------------------------------------
134!|| ale_connectivity_mod ../common_source/modules/ale/ale_connectivity_mod.F
135!||====================================================================
136 SUBROUTINE brezo2_sig(ALE_CONNECT,SIG,NV,PHI,NDIM,FLUX,VOL,BHOLE,NM,NEL)
137 USE ale_connectivity_mod
138C-----------------------------------------------
139C I m p l i c i t T y p e s
140C-----------------------------------------------
141#include "implicit_f.inc"
142C-----------------------------------------------
143C C o m m o n B l o c k s
144C-----------------------------------------------
145#include "com01_c.inc"
146#include "com08_c.inc"
147#include "vect01_c.inc"
148C-----------------------------------------------
149C D u m m y A r g u m e n t s
150C-----------------------------------------------
151 INTEGER NDIM,NV,NEL
152 INTEGER BHOLE(NMULT,*),NM
153 my_real sig(nel,ndim), phi(*), flux(4,*), vol(*)
154 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
155C-----------------------------------------------
156C L o c a l V a r i a b l e s
157C-----------------------------------------------
158 INTEGER I, IE, IV1, IV2, IV3, IV4, IAD2
159 my_real voln
160C-----------------------------------------------
161C S o u r c e L i n e s
162C-----------------------------------------------
163 IF(jmult == 0)THEN
164 DO i=lft,llt
165 ie =nft+i
166 iad2 = ale_connect%ee_connect%iad_connect(ie)
167 iv1=ale_connect%ee_connect%connected(iad2 + 1 - 1)
168 iv2=ale_connect%ee_connect%connected(iad2 + 2 - 1)
169 iv3=ale_connect%ee_connect%connected(iad2 + 3 - 1)
170 iv4=ale_connect%ee_connect%connected(iad2 + 4 - 1)
171
172 IF(iv1 <= zero)THEN
173 iv1=ie
174 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv1)) /= abs(bhole(nm,ie)))THEN
175 iv1=ie
176 ENDIF
177 IF(iv2 <= 0)THEN
178 iv2=ie
179 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv2)) /= abs(bhole(nm,ie)))THEN
180 iv2=ie
181 ENDIF
182 IF(iv3 <= 0)THEN
183 iv3=ie
184 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv3)) /= abs(bhole(nm,ie)))THEN
185 iv3=ie
186 ENDIF
187 IF(iv4 <= 0)THEN
188 iv4=ie
189 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv4)) /= abs(bhole(nm,ie)))THEN
190 iv4=ie
191 ENDIF
192 sig(i,nv) = sig(i,nv) + half * dt1 * (
193 . phi(ie)*(flux(1,i)+flux(2,i)+flux(3,i)+flux(4,i))
194 . -phi(iv1)*flux(1,i)-phi(iv2)*flux(2,i)-phi(iv3)*flux(3,i)-phi(iv4)*flux(4,i)
195 . ) / vol(i)
196 ENDDO
197 ELSE
198 DO i=lft,llt
199 ie =nft+i
200 iad2 = ale_connect%ee_connect%iad_connect(ie)
201 iv1=ale_connect%ee_connect%connected(iad2 + 1 - 1)
202 iv2=ale_connect%ee_connect%connected(iad2 + 2 - 1)
203 iv3=ale_connect%ee_connect%connected(iad2 + 3 - 1)
204 iv4=ale_connect%ee_connect%connected(iad2 + 4 - 1)
205 IF(iv1 <= zero)THEN
206 iv1=ie
207 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv1)) /= abs(bhole(nm,ie)))THEN
208 iv1=ie
209 ENDIF
210 IF(iv2 <= zero)THEN
211 iv2=ie
212 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv2)) /= abs(bhole(nm,ie)))THEN
213 iv2=ie
214 ENDIF
215 IF(iv3 <= zero)THEN
216 iv3=ie
217 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv3)) /= abs(bhole(nm,ie)))THEN
218 iv3=ie
219 ENDIF
220 IF(iv4 <= zero)THEN
221 iv4=ie
222 ELSEIF(bhole(nm,ie) < zero .AND.abs(bhole(nm,iv4)) /= abs(bhole(nm,ie)))THEN
223 iv4=ie
224 ENDIF
225 voln=vol(i)-dt1 * (flux(1,i)+flux(2,i)+flux(3,i)+flux(4,i))
226 IF (voln > em15) THEN
227 sig(i,nv) = sig(i,nv) + half * dt1 * (phi(ie)* (flux(1,i)+flux(2,i)+flux(3,i)+flux(4,i))
228 . -phi(iv1)*flux(1,i)-phi(iv2)*flux(2,i)-phi(iv3)*flux(3,i)-phi(iv4)*flux(4,i)) / max(em15,voln)
229 ENDIF
230 ENDDO !next I=LFT,LLT
231 ENDIF
232C-----------------------------------------------
233 RETURN
234 END
brezo2
subroutine brezo2(ale_connect, var, phi, flux, vol, bhole, nm)
Definition brezo2.F:31
brezo2_sig
subroutine brezo2_sig(ale_connect, sig, nv, phi, ndim, flux, vol, bhole, nm, nel)
Definition brezo2.F:137
my_real
#define my_real
Definition cppsort.cpp:32
max
#define max(a, b)
Definition macros.h:21
ale_connectivity_mod
Definition ale_connectivity_mod.F:223
ale_connectivity_mod::t_ale_connectivity
Definition ale_connectivity_mod.F:254