OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
alefvm_tfext.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!|| alefvm_tfext ../engine/source/ale/alefvm/alefvm_tfext.F
25!||--- called by ------------------------------------------------------
26!|| sforc3 ../engine/source/elements/solid/solide/sforc3.F
27!||--- uses -----------------------------------------------------
28!|| ale_connectivity_mod ../common_source/modules/ale/ale_connectivity_mod.f
29!|| alefvm_mod ../common_source/modules/ale/alefvm_mod.F
30!|| element_mod ../common_source/modules/elements/element_mod.F90
31!||====================================================================
32 SUBROUTINE alefvm_tfext (
33 1 IXS, NV46, SIG, ALE_CONNECT,
34 2 N1X, N2X, N3X, N4X, N5X, N6X ,
35 3 N1Y, N2Y, N3Y, N4Y, N5Y, N6Y ,
36 4 N1Z, N2Z, N3Z, N4Z, N5Z, N6Z ,
37 5 X , IPM, NALE,
38 6 W , P , NEL, WFEXT)
39C-----------------------------------------------
40C D e s c r i p t i o n
41C-----------------------------------------------
42C 'alefvm' is related to a collocated scheme (built from FVM and based on Godunov scheme)
43C which was temporarily introduced for experimental option /INTER/TYPE22 (FSI coupling with cut cell method)
44C This cut cell method is not completed, abandoned, and is not an official option.
45C There is no other use for this scheme which is automatically enabled when /INTER/TYPE22 is defined (INT22>0 => IALEFVM=1).
46C
47C This subroutine is treating an uncut cell.
48C-----------------------------------------------
49C M o d u l e s
50C-----------------------------------------------
51 USE alefvm_mod
52 USE ale_connectivity_mod
53 use element_mod , only : nixs
54C-----------------------------------------------
55C I m p l i c i t T y p e s
56C-----------------------------------------------
57#include "implicit_f.inc"
58C-----------------------------------------------
59C G l o b a l P a r a m e t e r s
60C-----------------------------------------------
61#include "mvsiz_p.inc"
62C-----------------------------------------------
63C C o m m o n B l o c k s
64C-----------------------------------------------
65#include "vect01_c.inc"
66#include "com01_c.inc"
67#include "com06_c.inc"
68#include "com08_c.inc"
69#include "param_c.inc"
70#include "comlock.inc"
71C-----------------------------------------------
72C D e s c r i p t i o n
73C-----------------------------------------------
74C This subroutines computes External force
75C work for internal energy computation.
76C
77CC WARNING THIS SUBROUTINE CAN BE OPTIMIZED (TEMPORARY ONE)
78C
79C-----------------------------------------------
80C D u m m y A r g u m e n t s
81C-----------------------------------------------
82 INTEGER :: IXS(NIXS,*),NV46,IPM(NPROPMI,*),NALE(*),NEL
83 my_real :: SIG(NEL,6),X(3,*),W(3,*),P(MVSIZ)
84 my_real :: N1X(*), N2X(*), N3X(*), N4X(*), N5X(*), N6X(*),
85 . N1Y(*), N2Y(*), N3Y(*), N4Y(*), N5Y(*), N6Y(*),
86 . N1Z(*), N2Z(*), N3Z(*), N4Z(*), N5Z(*), N6Z(*)
87 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
88 DOUBLE PRECISION,INTENT(INOUT) :: WFEXT
89C-----------------------------------------------
90C L o c a l V a r i a b l e s
91C-----------------------------------------------
92 INTEGER :: I, II, IV, J, MT, IALEFVM_FLG, ICF(4,6),IX(4)
93 INTEGER :: NC1(MVSIZ),NC2(MVSIZ),NC3(MVSIZ),NC4(MVSIZ),NC5(MVSIZ),NC6(MVSIZ),NC7(MVSIZ),NC8(MVSIZ)
94 my_real :: x1(mvsiz), x2(mvsiz), x3(mvsiz) , x4(mvsiz), x5(mvsiz), x6(mvsiz), x7(mvsiz), x8(mvsiz),
95 . y1(mvsiz), y2(mvsiz), y3(mvsiz) , y4(mvsiz), y5(mvsiz), y6(mvsiz), y7(mvsiz), y8(mvsiz),
96 . z1(mvsiz), z2(mvsiz), z3(mvsiz) , z4(mvsiz), z5(mvsiz), z6(mvsiz), z7(mvsiz), z8(mvsiz),
97 . swn(6) , wface(3,6,mvsiz), wfext_add, wfextt
98
99 DATA icf/1,4,3,2,3,4,8,7,5,6,7,8,1,2,6,5,2,3,7,6,1,5,8,4/
100 INTEGER :: IAD2
101C-----------------------------------------------
102C P r e - C o n d i t i o n s
103C-----------------------------------------------
104 IF(ALEFVM_Param%IEnabled==0)return
105 mt = ixs(1,nft+lft)
106 ialefvm_flg = ipm(251,mt)
107 IF(ialefvm_flg <= 1)RETURN
108
109 IF(jale == 0)RETURN
110C-----------------------------------------------
111C S o u r c e L i n e s
112C-----------------------------------------------
113
114 !-------------------------------------------------------------!
115 ! NORMAL VECTOR FOR ALE !
116 !-------------------------------------------------------------!
117 IF(jale==1)THEN
118 DO i=lft,llt
119 ii = i + nft
120 !---8 local node numbers NC1 TO NC8 for solid element I ---!
121 nc1(i)=ixs(2,ii)
122 nc2(i)=ixs(3,ii)
123 nc3(i)=ixs(4,ii)
124 nc4(i)=ixs(5,ii)
125 nc5(i)=ixs(6,ii)
126 nc6(i)=ixs(7,ii)
127 nc7(i)=ixs(8,ii)
128 nc8(i)=ixs(9,ii)
129 !
130 !---Coordinates of the 8 nodes
131 x1(i)=x(1,nc1(i))
132 y1(i)=x(2,nc1(i))
133 z1(i)=x(3,nc1(i))
134 !
135 x2(i)=x(1,nc2(i))
136 y2(i)=x(2,nc2(i))
137 z2(i)=x(3,nc2(i))
138 !
139 x3(i)=x(1,nc3(i))
140 y3(i)=x(2,nc3(i))
141 z3(i)=x(3,nc3(i))
142 !
143 x4(i)=x(1,nc4(i))
144 y4(i)=x(2,nc4(i))
145 z4(i)=x(3,nc4(i))
146 !
147 x5(i)=x(1,nc5(i))
148 y5(i)=x(2,nc5(i))
149 z5(i)=x(3,nc5(i))
150 !
151 x6(i)=x(1,nc6(i))
152 y6(i)=x(2,nc6(i))
153 z6(i)=x(3,nc6(i))
154 !
155 x7(i)=x(1,nc7(i))
156 y7(i)=x(2,nc7(i))
157 z7(i)=x(3,nc7(i))
158 !
159 x8(i)=x(1,nc8(i))
160 y8(i)=x(2,nc8(i))
161 z8(i)=x(3,nc8(i))
162 ENDDO
163 DO i=lft,llt
164 ! Face-1
165 n1x(i)=(y3(i)-y1(i))*(z2(i)-z4(i)) - (z3(i)-z1(i))*(y2(i)-y4(i))
166 n1y(i)=(z3(i)-z1(i))*(x2(i)-x4(i)) - (x3(i)-x1(i))*(z2(i)-z4(i))
167 n1z(i)=(x3(i)-x1(i))*(y2(i)-y4(i)) - (y3(i)-y1(i))*(x2(i)-x4(i))
168 ! Face-2
169 n2x(i)=(y7(i)-y4(i))*(z3(i)-z8(i)) - (z7(i)-z4(i))*(y3(i)-y8(i))
170 n2y(i)=(z7(i)-z4(i))*(x3(i)-x8(i)) - (x7(i)-x4(i))*(z3(i)-z8(i))
171 n2z(i)=(x7(i)-x4(i))*(y3(i)-y8(i)) - (y7(i)-y4(i))*(x3(i)-x8(i))
172 ! Face-3
173 n3x(i)=(y6(i)-y8(i))*(z7(i)-z5(i)) - (z6(i)-z8(i))*(y7(i)-y5(i))
174 n3y(i)=(z6(i)-z8(i))*(x7(i)-x5(i)) - (x6(i)-x8(i))*(z7(i)-z5(i))
175 n3z(i)=(x6(i)-x8(i))*(y7(i)-y5(i)) - (y6(i)-y8(i))*(x7(i)-x5(i))
176 ! Face-4
177 n4x(i)=(y2(i)-y5(i))*(z6(i)-z1(i)) - (z2(i)-z5(i))*(y6(i)-y1(i))
178 n4y(i)=(z2(i)-z5(i))*(x6(i)-x1(i)) - (x2(i)-x5(i))*(z6(i)-z1(i))
179 n4z(i)=(x2(i)-x5(i))*(y6(i)-y1(i)) - (y2(i)-y5(i))*(x6(i)-x1(i))
180 ! Face-5
181 n5x(i)=(y7(i)-y2(i))*(z6(i)-z3(i)) - (z7(i)-z2(i))*(y6(i)-y3(i))
182 n5y(i)=(z7(i)-z2(i))*(x6(i)-x3(i)) - (x7(i)-x2(i))*(z6(i)-z3(i))
183 n5z(i)=(x7(i)-x2(i))*(y6(i)-y3(i)) - (y7(i)-y2(i))*(x6(i)-x3(i))
184 ! Face-6
185 n6x(i)=(y8(i)-y1(i))*(z4(i)-z5(i)) - (z8(i)-z1(i))*(y4(i)-y5(i))
186 n6y(i)=(z8(i)-z1(i))*(x4(i)-x5(i)) - (x8(i)-x1(i))*(z4(i)-z5(i))
187 n6z(i)=(x8(i)-x1(i))*(y4(i)-y5(i)) - (y8(i)-y1(i))*(x4(i)-x5(i))
188 ENDDO
189
190 DO i=lft,llt
191 wface(1,1,i) = fourth*(w(1,nc1(i))+w(1,nc2(i))+w(1,nc3(i))+w(1,nc4(i)))
192 wface(2,1,i) = fourth*(w(2,nc1(i))+w(2,nc2(i))+w(2,nc3(i))+w(2,nc4(i)))
193 wface(3,1,i) = fourth*(w(3,nc1(i))+w(3,nc2(i))+w(3,nc3(i))+w(3,nc4(i)))
194
195 wface(1,2,i) = fourth*(w(1,nc3(i))+w(1,nc4(i))+w(1,nc7(i))+w(1,nc8(i)))
196 wface(2,2,i) = fourth*(w(2,nc3(i))+w(2,nc4(i))+w(2,nc7(i))+w(2,nc8(i)))
197 wface(3,2,i) = fourth*(w(3,nc3(i))+w(3,nc4(i))+w(3,nc7(i))+w(3,nc8(i)))
198
199 wface(1,3,i) = fourth*(w(1,nc5(i))+w(1,nc6(i))+w(1,nc7(i))+w(1,nc8(i)))
200 wface(2,3,i) = fourth*(w(2,nc5(i))+w(2,nc6(i))+w(2,nc7(i))+w(2,nc8(i)))
201 wface(3,3,i) = fourth*(w(3,nc5(i))+w(3,nc6(i))+w(3,nc7(i))+w(3,nc8(i)))
202
203 wface(1,4,i) = fourth*(w(1,nc1(i))+w(1,nc2(i))+w(1,nc5(i))+w(1,nc6(i)))
204 wface(2,4,i) = fourth*(w(2,nc1(i))+w(2,nc2(i))+w(2,nc5(i))+w(2,nc6(i)))
205 wface(3,4,i) = fourth*(w(3,nc1(i))+w(3,nc2(i))+w(3,nc5(i))+w(3,nc6(i)))
206
207 wface(1,5,i) = fourth*(w(1,nc2(i))+w(1,nc3(i))+w(1,nc6(i))+w(1,nc7(i)))
208 wface(2,5,i) = fourth*(w(2,nc2(i))+w(2,nc3(i))+w(2,nc6(i))+w(2,nc7(i)))
209 wface(3,5,i) = fourth*(w(3,nc2(i))+w(3,nc3(i))+w(3,nc6(i))+w(3,nc7(i)))
210
211 wface(1,6,i) = fourth*(w(1,nc1(i))+w(1,nc4(i))+w(1,nc5(i))+w(1,nc8(i)))
212 wface(2,6,i) = fourth*(w(2,nc1(i))+w(2,nc4(i))+w(2,nc5(i))+w(2,nc8(i)))
213 wface(3,6,i) = fourth*(w(3,nc1(i))+w(3,nc4(i))+w(3,nc5(i))+w(3,nc8(i)))
214 ENDDO
215
216 !-------------------------------------------------------------!
217 ! NORMAL VECTOR FOR EULER !
218 !-------------------------------------------------------------!
219 !ELSE
220 ! DO i=lft,llt
221 ! II = I+NFT
222 ! N1X(I)=VEUL(14,II)
223 ! N2X(I)=VEUL(15,II)
224 ! N3X(I)=VEUL(16,II)
225 ! N4X(I)=VEUL(17,II)
226 ! N5X(I)=VEUL(18,II)
227 ! N6X(I)=VEUL(19,II)
228 ! N1Y(I)=VEUL(20,II)
229 ! N2Y(I)=VEUL(21,II)
230 ! N3Y(I)=VEUL(22,II)
231 ! N4Y(I)=VEUL(23,II)
232 ! N5Y(I)=VEUL(24,II)
233 ! N6Y(I)=VEUL(25,II)
234 ! N1Z(I)=VEUL(26,II)
235 ! N2Z(I)=VEUL(27,II)
236 ! N3Z(I)=VEUL(28,II)
237 ! N4Z(I)=VEUL(29,II)
238 ! N5Z(I)=VEUL(30,II)
239 ! N6Z(I)=VEUL(31,II)
240 ! ENDDO
241 ! DO I=LFT,LLT
242 ! Wface(1:3,1:6,I) = ZERO
243 ! ENDDO
244
245 ENDIF
246
247 !-------------------------------------------------------------!
248 ! EXTERNAL FORCE WORK !
249 !-------------------------------------------------------------!
250 alefvm_buffer%WFEXT_CELL(lft+nft:llt+nft) = zero
251 wfextt = zero
252 DO i=lft,llt
253 ii = i + nft
254 iad2 = ale_connect%ee_connect%iad_connect(ii)
255 ! < W, 2S.n> = 2S.wn
256 swn(1) = wface(1,1,i)*n1x(i) + wface(2,1,i)*n1y(i) + wface(3,1,i)*n1z(i)
257 swn(2) = wface(1,2,i)*n2x(i) + wface(2,2,i)*n2y(i) + wface(3,2,i)*n2z(i)
258 swn(3) = wface(1,3,i)*n3x(i) + wface(2,3,i)*n3y(i) + wface(3,3,i)*n3z(i)
259 swn(4) = wface(1,4,i)*n4x(i) + wface(2,4,i)*n4y(i) + wface(3,4,i)*n4z(i)
260 swn(5) = wface(1,5,i)*n5x(i) + wface(2,5,i)*n5y(i) + wface(3,5,i)*n5z(i)
261 swn(6) = wface(1,6,i)*n6x(i) + wface(2,6,i)*n6y(i) + wface(3,6,i)*n6z(i)
262 wfext_add = zero
263 Do j=1,nv46
264 iv = ale_connect%ee_connect%connected(iad2 + j - 1)
265 !IF(IV /= 0)CYCLE !domain border only
266 ix(1) = ixs(icf(1,j)+1,ii)
267 ix(2) = ixs(icf(2,j)+1,ii)
268 ix(3) = ixs(icf(3,j)+1,ii)
269 ix(4) = ixs(icf(4,j)+1,ii)
270 IF(sum(iabs(nale(ix(1:4))))==4)cycle !4 lagrangian nodes only
271 !dt * P*S*wn
272 wfext_add = wfext_add - dt1*third*(sig(i,1)+sig(i,2)+sig(i,3))*swn(j)*half !Swn=2S.w.n
273 ENDDO
274 alefvm_buffer%WFEXT_CELL(ii) = wfext_add
275 wfextt = wfextt + wfext_add
276 ENDDO
277
278#include "lockon.inc"
279 wfext = wfext + wfextt
280#include "lockoff.inc"
281
282
283 !DEBUG-OUTPUT---------------!
284 if(alefvm_param%IOUTP_WFEXT == 1)then
285 IF(wfextt /= zero)THEN
286!#!include "lockon.inc"
287 print *, " |----alefvm_tfext.F------|"
288 print *, " | THREAD INFORMATION |"
289 print *, " |------------------------|"
290 print *, " NCYCLE =", ncycle
291 do i=lft,llt
292 ii = nft + i
293 if(alefvm_buffer%WFEXT_CELL(ii)==zero)cycle
294 write(*,fmt='(A,I10,A,F30.16,A,F30.16,A,6F30.16)') " brique=", ixs(11,ii)," Wfext=",
295 . alefvm_buffer%WFEXT_CELL(ii), "P(I)=", p(i), "SWn(1:6)=",swn(1:6)
296 enddo
297 ENDIF
298!#!include "lockoff.inc"
299 endif
300 !-----------------------------------------!
301
302 RETURN
303 END
alefvm_tfext
subroutine alefvm_tfext(ixs, nv46, sig, ale_connect, n1x, n2x, n3x, n4x, n5x, n6x, n1y, n2y, n3y, n4y, n5y, n6y, n1z, n2z, n3z, n4z, n5z, n6z, x, ipm, nale, w, p, nel, wfext)
Definition alefvm_tfext.F:39
ale_connectivity_mod
Definition ale_connectivity_mod.F:225
alefvm_mod
Definition alefvm_mod.F:88
alefvm_mod::alefvm_buffer
type(alefvm_buffer_), target alefvm_buffer
Definition alefvm_mod.F:120