OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
ebcs_vol2seg.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!|| ebcs_vol2seg ../engine/source/boundary_conditions/ebcs/ebcs_vol2seg.F
25!||--- called by ------------------------------------------------------
26!|| ebcs_main ../engine/source/boundary_conditions/ebcs/ebcs_main.F
27!||--- uses -----------------------------------------------------
28!|| ale_mod ../common_source/modules/ale/ale_mod.F
29!|| segvar_mod ../engine/share/modules/segvar_mod.F
30!||====================================================================
31 SUBROUTINE ebcs_vol2seg(NSEG,SURF_NODES,ISEG,SEGVAR,
32 . A,V,X,VOLMON,FSAV)
33C-----------------------------------------------
34C M o d u l e s
35C-----------------------------------------------
36 USE segvar_mod
37 USE ale_mod
38C-----------------------------------------------
39C I m p l i c i t T y p e s
40C-----------------------------------------------
41#include "implicit_f.inc"
42C-----------------------------------------------
43C D u m m y A r g u m e n t s
44C-----------------------------------------------
45 INTEGER NSEG,ISEG(NSEG),SURF_NODES(NSEG,4)
46 my_real
47 . a(3,*),x(3,*),v(3,*),volmon(*),fsav(*)
48 TYPE(t_segvar) :: SEGVAR
49C-----------------------------------------------
50C L o c a l V a r i a b l e s
51C-----------------------------------------------
52 INTEGER IS,K,KSEG,SEGAD,IRECT(4)
53 my_real orient,
54 . pext,pres,vol,ro,en,
55 . fac,x13,y13,z13,x24,y24,z24,nx,ny,nz,s,
56 . vn1,vn2,vn3,vn4,vel,vn,vplus,vmoins,
57 . gam1,gamrp,vcrt2,fluxm,fluxe,flux,
58 . gv,rv,rsr,pn,enin,roin,enou,roou,pns,fx,fy,fz,
59 . gam,norm,dmin,dmou,hin,hou,pna
60 my_real pna1,pna2,pna3,pna4,rsr1,rsr2,rsr3,rsr4,aire
61C=======================================================================
62C VOLUME MONITORE
63 gam =volmon(1)
64 pext =volmon(3)
65 pres =volmon(12)
66 vol =volmon(16)
67 ro =volmon(20)/vol
68C volumetric Energy
69 en =volmon(13)/vol
70C massic Enthalpy
71 hin = (en+pres)/ro
72
73 gam1 = one/(gam-one)
74 gamrp= half*(gam-one)/gam
75 vcrt2= two*gam*pres/ro/(gam+one)
76
77 flux =zero
78 fluxm=zero
79 fluxe=zero
80 aire =zero
81
82 DO is=1,nseg
83 DO k=1,4
84 irect(k)=surf_nodes(is,k)
85 ENDDO
86 IF(irect(4) == 0) THEN
87 irect(4)=irect(3)
88 fac=third
89 ELSE
90 fac=fourth
91 ENDIF
92C
93C INDICES FOR SEGMENTS, AND OREINTATIONS
94C
95 kseg=abs(iseg(is))
96 orient=float(iseg(is)/kseg)
97
98C NORMAL VECTORS
99 x13=x(1,irect(3))-x(1,irect(1))
100 y13=x(2,irect(3))-x(2,irect(1))
101 z13=x(3,irect(3))-x(3,irect(1))
102 x24=x(1,irect(4))-x(1,irect(2))
103 y24=x(2,irect(4))-x(2,irect(2))
104 z24=x(3,irect(4))-x(3,irect(2))
105 nx=y13*z24-z13*y24
106 ny=z13*x24-x13*z24
107 nz=x13*y24-y13*x24
108 norm=sqrt(nx*nx+ny*ny+nz*nz)
109 s=half*norm
110 IF(norm>em20)THEN
111 norm=orient/norm
112 nx=nx*norm
113 ny=ny*norm
114 nz=nz*norm
115 ENDIF
116
117
118C NORMAL VELOCITIES
119c write(6,*)iseg(is),nx,ny,nz
120 vn1=v(1,irect(1))*nx+v(2,irect(1))*ny+v(3,irect(1))*nz
121 vn2=v(1,irect(2))*nx+v(2,irect(2))*ny+v(3,irect(2))*nz
122 vn3=v(1,irect(3))*nx+v(2,irect(3))*ny+v(3,irect(3))*nz
123 vn4=v(1,irect(4))*nx+v(2,irect(4))*ny+v(3,irect(4))*nz
124 IF(irect(4) == irect(3))THEN
125 vn=fac*(vn1+vn2+vn3)
126 ELSE
127 vn=fac*(vn1+vn2+vn3+vn4)
128 ENDIF
129 vplus =max(vn,zero)
130 vmoins=min(vn,zero)
131C
132C STAGNATION POINT
133C
134 vel=min(vn1,zero)**2
135 vel=min(vel,vcrt2)
136 gv=gamrp*ro*vel/pres
137 IF(gv>em4)THEN
138 rsr = (one-gv)**gam1
139 pna1 = pres*rsr**gam
140 ELSE
141 rv = half*ro*vel
142 rsr = one - rv/(pres*gam)
143 pna1 = pres - rv
144 ENDIF
145 rsr1=rsr
146C
147 vel=min(vn2,zero)**2
148 vel=min(vel,vcrt2)
149 gv=gamrp*ro*vel/pres
150 IF(gv>em4)THEN
151 rsr = (one-gv)**gam1
152 pna2 = pres*rsr**gam
153 ELSE
154 rv = half*ro*vel
155 rsr = one - rv/(pres*gam)
156 pna2 = pres - rv
157 ENDIF
158 rsr2=rsr
159C
160 vel=min(vn3,zero)**2
161 vel=min(vel,vcrt2)
162 gv=gamrp*ro*vel/pres
163 IF(gv>em4)THEN
164 rsr = (one-gv)**gam1
165 pna3 = pres*rsr**gam
166 ELSE
167 rv = half*ro*vel
168 rsr = one - rv/(pres*gam)
169 pna3 = pres - rv
170 ENDIF
171 rsr3=rsr
172C
173 vel=min(vn4,zero)**2
174 vel=min(vel,vcrt2)
175 gv=gamrp*ro*vel/pres
176 IF(gv>em4)THEN
177 rsr = (one-gv)**gam1
178 pna4 = pres*rsr**gam
179 ELSE
180 rv = half*ro*vel
181 rsr = one - rv/(pres*gam)
182 pna4 = pres - rv
183 ENDIF
184 rsr4=rsr
185C
186 IF(irect(4) == irect(3))THEN
187 pna=fac*(pna1+pna2+pna3)
188 rsr=fac*(rsr1+rsr2+rsr3)
189 ELSE
190 pna=fac*(pna1+pna2+pna3+pna4)
191 rsr=fac*(rsr1+rsr2+rsr3+rsr4)
192 ENDIF
193C INCOMING FLOW
194 enin = pna*gam1
195 roin = rsr*ro
196C ROIN =RO
197C ENIN =EN
198C OUTGOING FLOW
199 segad=ale%GLOBAL%NVCONV*(kseg-1)
200 roou = segvar%RHO(kseg)
201 enou = segvar%EINT(kseg)
202
203 segvar%RHO(kseg)=roin
204 segvar%EINT(kseg)=enin
205C
206C MASS FLUX AND ENTHALPY TO AIRBAGS (>0 si rentrant)
207C airbag -> fluid
208C hin = (Ea+Pa)/roa
209C fluide -> airbag
210C hou = (Ef+P)/rof +1/2 v+^2
211C HOU=(ENOU+PNA)/ROOU+0.5*VPLUS**2
212c
213 hou=gam*enou/roou+half*vplus**2
214 dmin=roin*vmoins
215 dmou=roou*vplus
216 aire =aire +s
217 flux =flux +s*vn
218 fluxm=fluxm+s*(dmou+dmin)
219 fluxe=fluxe+s*(hin*dmin+hou*dmou)
220
221C RELATIVE PRESSURE
222 pn = pna1-pext
223 pns=-pn*s*fac
224 fx=pns*nx
225 fy=pns*ny
226 fz=pns*nz
227
228 a(1,irect(1))=a(1,irect(1))+fx
229 a(2,irect(1))=a(2,irect(1))+fy
230 a(3,irect(1))=a(3,irect(1))+fz
231
232 pn = pna2-pext
233 pns=-pn*s*fac
234 fx=pns*nx
235 fy=pns*ny
236 fz=pns*nz
237 a(1,irect(2))=a(1,irect(2))+fx
238 a(2,irect(2))=a(2,irect(2))+fy
239 a(3,irect(2))=a(3,irect(2))+fz
240
241 pn = pna3-pext
242 pns=-pn*s*fac
243 fx=pns*nx
244 fy=pns*ny
245 fz=pns*nz
246 a(1,irect(3))=a(1,irect(3))+fx
247 a(2,irect(3))=a(2,irect(3))+fy
248 a(3,irect(3))=a(3,irect(3))+fz
249
250 IF(irect(4)/=irect(3))THEN
251 pn = pna4-pext
252 pns=-pn*s*fac
253 fx=pns*nx
254 fy=pns*ny
255 fz=pns*nz
256 a(1,irect(4))=a(1,irect(4))+fx
257 a(2,irect(4))=a(2,irect(4))+fy
258 a(3,irect(4))=a(3,irect(4))+fz
259 ENDIF
260
261 ENDDO
262
263 volmon(22)=volmon(22) - fluxe
264 volmon(24)=volmon(24) - fluxm
265 flux=-flux+fsav(9)*fsav(8)
266 fsav(8)=fsav(8) + aire
267 fsav(9)=flux/fsav(8)
268C-----------
269 RETURN
270 END
my_real
#define my_real
Definition cppsort.cpp:32
norm
norm(diag(diag(diag(inv(mat))) -id.SOL), 2) % destroy mumps instance id.JOB
ebcs_vol2seg
subroutine ebcs_vol2seg(nseg, surf_nodes, iseg, segvar, a, v, x, volmon, fsav)
Definition ebcs_vol2seg.F:33
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
ale_mod
Definition ale_mod.F:149
ale_mod::ale
type(ale_) ale
Definition ale_mod.F:249
segvar_mod
Definition segvar_mod.F:61
segvar_mod::t_segvar
Definition segvar_mod.F:69