OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
multifluid_init3.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!|| multifluid_init3 ../starter/source/multifluid/multifluid_init3.F
25!||--- called by ------------------------------------------------------
26!|| initia ../starter/source/elements/initia/initia.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| atheri ../starter/source/ale/atheri.F
30!|| aturi3 ../starter/source/ale/ale3d/aturi3.F
31!|| dtmain ../starter/source/materials/time_step/dtmain.F
32!|| fretitl2 ../starter/source/starter/freform.F
33!|| m5in3 ../starter/source/initial_conditions/detonation/m5in3.F
34!|| matini ../starter/source/materials/mat_share/matini.F
35!|| scoor3 ../starter/source/elements/solid/solide/scoor3.F
36!|| sderi3 ../starter/source/elements/solid/solide/sderi3.F
37!|| sdlen3 ../starter/source/elements/solid/solide/sdlen3.F
38!|| sjac_i ../starter/source/elements/solid/solide/sderi3.F
39!|| smass3 ../starter/source/elements/solid/solide/smass3.F
40!|| smorth3 ../starter/source/elements/solid/solide/smorth3.F
41!|| srcoor3 ../starter/source/elements/solid/solide/srcoor3.F
42!|| sveok3 ../starter/source/elements/solid/solide/sveok3.F
43!||--- uses -----------------------------------------------------
44!|| detonators_mod ../starter/share/modules1/detonators_mod.F
45!|| message_mod ../starter/share/message_module/message_mod.F
46!||====================================================================
47 SUBROUTINE multifluid_init3(ELBUF_STR,MAS ,IXS ,PM ,X ,
48 . GEO ,ALE_CONNECTIVITY ,IPARG_GR,
49 . DTELEM ,SIGI ,NEL ,SKEW ,IGEO ,
50 . STIFN ,PARTSAV ,V ,IPARTS ,MSS ,
51 . IPART ,SIGSP ,NG ,IPARG ,GLOB_THERM,
52 . NSIGI ,MSNF ,NVC ,MSSF ,IPM ,
53 . IUSER ,NSIGS ,VOLNOD ,BVOLNOD ,VNS ,
54 . BNS ,IN ,VR ,INS ,WMA ,
55 . PTSOL ,BUFMAT ,MCP ,MCPS ,TEMP ,
56 . XREFS ,NPF ,TF ,MSSA ,STRSGLOB,
57 . STRAGLOB ,FAIL_INI ,SPBUF ,KXSP ,IPARTSP ,
58 . NOD2SP ,SOL2SPH ,IRST ,ILOADP ,FACLOAD,
59 . MULTI_FVM, ERROR_THROWN,DETONATORS,MAT_PARAM)
60C-----------------------------------------------
61C M o d u l e s
62C-----------------------------------------------
63 USE elbufdef_mod
64 USE message_mod
65 USE multi_fvm_mod
66 USE detonators_mod
67 USE ale_connectivity_mod
68 USE names_and_titles_mod , ONLY : nchartitle
69 USE matparam_def_mod, ONLY : matparam_struct_
70 use glob_therm_mod
71C-----------------------------------------------
72C I m p l i c i t T y p e s
73C-----------------------------------------------
74#include "implicit_f.inc"
75C-----------------------------------------------
76C G l o b a l P a r a m e t e r s
77C-----------------------------------------------
78#include "mvsiz_p.inc"
79C-----------------------------------------------
80C C o m m o n B l o c k s
81C-----------------------------------------------
82#include "com01_c.inc"
83#include "com04_c.inc"
84#include "param_c.inc"
85#include "scr03_c.inc"
86#include "scr17_c.inc"
87#include "scry_c.inc"
88#include "sphcom.inc"
89#include "vect01_c.inc"
90C-----------------------------------------------
91C D u m m y A r g u m e n t s
92C-----------------------------------------------
93 INTEGER IXS(NIXS,*),IPARG(NPARG,NGROUP),
94 . IPARG_GR(NPARG),IPARTS(*),IGEO(NPROPGI,*),
95 . IPM(NPROPMI,NUMMAT),IPART(LIPART1,*),PTSOL(*),
96 . NG, NSIGI ,NVC,NEL,IUSER, NSIGS, NPF(*),
97 . STRSGLOB(*),STRAGLOB(*),FAIL_INI(*),
98 . KXSP(NISP,*), IPARTSP(*), NOD2SP(*), SOL2SPH(2,*), IRST(3,*)
99 my_real
100 . MAS(*), PM(NPROPM,NUMMAT), X(3,*), GEO(NPROPG,*),
101 . DTELEM(*),SIGI(NSIGS,*),SKEW(LSKEW,*),STIFN(*),
102 . PARTSAV(20,*), V(3, *), MSS(8,*),
103 . SIGSP(NSIGI,*),MSNF(*), MSSF(8,*), WMA(*),
104 . VOLNOD(*), BVOLNOD(*), VNS(8,*), BNS(8,*),
105 . in(*),vr(*), ins(8,*),bufmat(*),
106 . mcp(*), mcps(8,*), temp(*),
107 . xrefs(8,3,*), tf(*), mssa(*),
108 . spbuf(nspbuf,*)
109 TYPE(elbuf_struct_), TARGET :: ELBUF_STR
110 INTEGER,INTENT(IN) :: ILOADP(SIZLOADP,*)
111 my_real,INTENT(IN) :: facload(lfacload,*)
112 TYPE(multi_fvm_struct) :: MULTI_FVM
113 LOGICAL :: ERROR_THROWN
114 TYPE(DETONATORS_STRUCT_) :: DETONATORS
115 TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
116 TYPE(MATPARAM_STRUCT_) ,DIMENSION(NUMMAT) ,INTENT(IN) :: MAT_PARAM
117 type (glob_therm_) ,intent(in) :: glob_therm
118C-----------------------------------------------
119C L o c a l V a r i a b l e s
120C-----------------------------------------------
121 INTEGER I,II,J, NF1, IBID, JHBE, IREP,IGTYP, NUVAR, IDEF,
122 . ipt,lvloc,ipid1,nptr,npts,nptt,nlay,
123 . nsphdir, l_pla
124 INTEGER MAT(MVSIZ), PID(MVSIZ), NGL(MVSIZ),
125 . ix1(mvsiz),ix2(mvsiz),ix3(mvsiz),ix4(mvsiz),
126 . ix5(mvsiz),ix6(mvsiz),ix7(mvsiz),ix8(mvsiz)
127 my_real
128 . volu(mvsiz),dtx(mvsiz),
129 . x1(mvsiz),x2(mvsiz),x3(mvsiz),x4(mvsiz),x5(mvsiz),x6(mvsiz),
130 . x7(mvsiz),x8(mvsiz),y1(mvsiz),y2(mvsiz),y3(mvsiz),y4(mvsiz),
131 . y5(mvsiz),y6(mvsiz),y7(mvsiz),y8(mvsiz),z1(mvsiz),z2(mvsiz),
132 . z3(mvsiz),z4(mvsiz),z5(mvsiz),z6(mvsiz),z7(mvsiz),z8(mvsiz),
133 . rx(mvsiz) ,ry(mvsiz) ,rz(mvsiz) ,sx(mvsiz) ,
134 . sy(mvsiz) ,sz(mvsiz) ,tx(mvsiz) ,ty(mvsiz) ,tz(mvsiz) ,
135 . e1x(mvsiz),e1y(mvsiz),e1z(mvsiz),
136 . e2x(mvsiz),e2y(mvsiz),e2z(mvsiz),
137 . e3x(mvsiz),e3y(mvsiz),e3z(mvsiz),
138 . f1x(mvsiz) ,f1y(mvsiz) ,f1z(mvsiz) ,
139 . f2x(mvsiz) ,f2y(mvsiz) ,f2z(mvsiz),rhocp(mvsiz),temp0(mvsiz),
140 . px1(mvsiz),px2(mvsiz),px3(mvsiz),px4(mvsiz),
141 . py1(mvsiz),py2(mvsiz),py3(mvsiz),py4(mvsiz),
142 . pz1(mvsiz),pz2(mvsiz),pz3(mvsiz),pz4(mvsiz),
143 . rhof(mvsiz),alpha(mvsiz), deltax(mvsiz), aire(mvsiz), dummy, pres, vfrac
144 my_real :: tempel(nel)
145 my_real
146 . bid, fv, sti
147 double precision
148 . xd1(mvsiz), xd2(mvsiz), xd3(mvsiz), xd4(mvsiz),
149 . xd5(mvsiz), xd6(mvsiz), xd7(mvsiz), xd8(mvsiz),
150 . yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz),
151 . yd5(mvsiz), yd6(mvsiz), yd7(mvsiz), yd8(mvsiz),
152 . zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz),
153 . zd5(mvsiz), zd6(mvsiz), zd7(mvsiz), zd8(mvsiz),voldp(mvsiz)
154 INTEGER :: ILAY, MATLAW
155C-----------------------------------------------
156 CHARACTER(LEN=NCHARTITLE) :: TITR1
157 parameter(lvloc = 51)
158C-----------------------------------------------
159 TYPE(l_bufel_) ,POINTER :: LBUF
160 TYPE(G_BUFEL_) ,POINTER :: GBUF
161 TYPE(BUF_MAT_) ,POINTER :: MBUF
162 TYPE(BUF_LAY_) ,POINTER :: BUFLY
163C-----------------------------------------------
164C S o u r c e L i n e s
165C=======================================================================
166C Global buffer
167 gbuf => elbuf_str%GBUF
168C Number of layer
169 nlay = elbuf_str%NLAY
170C Current group starts at NF1
171 nf1=nft+1
172C ???
173 jhbe = iparg_gr(23)
174 irep = iparg_gr(35)
175 jcvt = iparg_gr(37)
176 igtyp = iparg_gr(38)
177 IF (jcvt==1.AND.isorth/=0) jcvt=2
178 idef = 0 ! initialization flag for the total strain
179 bid = zero
180 ibid = 0
181 nptr = elbuf_str%NPTR
182 npts = elbuf_str%NPTS
183 nptt = elbuf_str%NPTT
184
185 tempel(:) = zero
186 DO i=1,nel
187 rhocp(i) = pm(69,ixs(1,nft+i))
188 temp0(i) = pm(79,ixs(1,nft+i))
189C for air + foam
190 rhof(i) = pm(192,ixs(1,nft+i))
191 alpha(i) = pm(193,ixs(1,nft+i))
192 ENDDO
193C-----JAC_I [J]^-1 is calculated in global system
194 IF (ismstr==10.OR.ismstr==12) THEN
195 CALL scoor3(x,xrefs(1,1,nf1),ixs(1,nf1), geo ,mat ,pid ,ngl ,
196 . ix1 ,ix2 ,ix3 ,ix4 ,ix5 ,ix6 ,ix7 ,ix8 ,
197 . x1 ,x2 ,x3 ,x4 ,x5 ,x6 ,x7 ,x8 ,
198 . y1 ,y2 ,y3 ,y4 ,y5 ,y6 ,y7 ,y8 ,
199 . z1 ,z2 ,z3 ,z4 ,z5 ,z6 ,z7 ,z8 ,
200 . rx ,ry ,rz ,sx ,sy ,sz ,tx ,ty ,tz ,
201 . e1x ,e1y ,e1z ,e2x ,e2y ,e2z ,e3x ,e3y ,e3z ,
202 . f1x ,f1y ,f1z ,f2x ,f2y ,f2z ,temp0,temp ,glob_therm%NINTEMP,
203 . xd1 ,xd2 ,xd3 ,xd4 ,xd5 ,xd6 ,xd7 ,xd8 ,
204 . yd1 ,yd2 ,yd3 ,yd4 ,yd5 ,yd6 ,yd7 ,yd8 ,
205 . zd1 ,zd2 ,zd3 ,zd4 ,zd5 ,zd6 ,zd7 ,zd8 )
206 CALL sjac_i(
207 . x1 ,x2 ,x3 ,x4 ,x5 ,x6 ,x7 ,x8 ,
208 . y1 ,y2 ,y3 ,y4 ,y5 ,y6 ,y7 ,y8 ,
209 . z1 ,z2 ,z3 ,z4 ,z5 ,z6 ,z7 ,z8 ,
210 . gbuf%JAC_I ,nel)
211 END IF
212C Orthotropy wrt reference geometry
213 IF (jcvt == 0) THEN
214 CALL scoor3(x,xrefs(1,1,nf1),ixs(1,nf1),geo ,mat ,pid ,ngl ,
215 . ix1 ,ix2 ,ix3 ,ix4 ,ix5 ,ix6 ,ix7 ,ix8 ,
216 . x1 ,x2 ,x3 ,x4 ,x5 ,x6 ,x7 ,x8 ,
217 . y1 ,y2 ,y3 ,y4 ,y5 ,y6 ,y7 ,y8 ,
218 . z1 ,z2 ,z3 ,z4 ,z5 ,z6 ,z7 ,z8 ,
219 . rx ,ry ,rz ,sx ,sy ,sz ,tx ,ty ,tz ,
220 . e1x ,e1y ,e1z ,e2x ,e2y ,e2z ,e3x ,e3y ,e3z ,
221 . f1x ,f1y ,f1z ,f2x ,f2y ,f2z ,temp0,temp ,glob_therm%NINTEMP,
222 . xd1 ,xd2 ,xd3 ,xd4 ,xd5 ,xd6 ,xd7 ,xd8 ,
223 . yd1 ,yd2 ,yd3 ,yd4 ,yd5 ,yd6 ,yd7 ,yd8 ,
224 . zd1 ,zd2 ,zd3 ,zd4 ,zd5 ,zd6 ,zd7 ,zd8 )
225 ELSE
226 CALL srcoor3(x,xrefs(1,1,nf1),ixs(1,nf1),geo ,mat ,pid ,ngl ,jhbe ,
227 . ix1 ,ix2 ,ix3 ,ix4 ,ix5 ,ix6 ,ix7 ,ix8 ,
228 . x1 ,x2 ,x3 ,x4 ,x5 ,x6 ,x7 ,x8 ,
229 . y1 ,y2 ,y3 ,y4 ,y5 ,y6 ,y7 ,y8 ,
230 . z1 ,z2 ,z3 ,z4 ,z5 ,z6 ,z7 ,z8 ,
231 . rx ,ry ,rz ,sx ,sy ,sz ,tx ,ty ,tz ,
232 . e1x ,e1y ,e1z ,e2x ,e2y ,e2z ,e3x ,e3y ,e3z ,
233 . f1x ,f1y ,f1z ,f2x ,f2y ,f2z ,temp0,temp ,glob_therm%NINTEMP,
234 . xd1 ,xd2 ,xd3 ,xd4 ,xd5 ,xd6 ,xd7 ,xd8 ,
235 . yd1 ,yd2 ,yd3 ,yd4 ,yd5 ,yd6 ,yd7 ,yd8 ,
236 . zd1 ,zd2 ,zd3 ,zd4 ,zd5 ,zd6 ,zd7 ,zd8 )
237
238 ENDIF
239C
240C Orthotropy
241 IF (isorth == 1)
242 . CALL smorth3(pid ,geo ,igeo ,skew ,irep ,gbuf%GAMA ,
243 . rx ,ry ,rz ,sx ,sy ,sz ,tx ,ty ,tz ,
244 . e1x ,e1y ,e1z ,e2x ,e2y ,e2z ,e3x ,e3y ,e3z ,
245 . f1x ,f1y ,f1z ,f2x ,f2y ,f2z ,nsigi,sigsp,nsigs,
246 . sigi ,ixs ,x ,jhbe ,ptsol,nel ,iparg_gr(28))
247C
248C-----------
249 CALL sveok3(nvc,8, ix1, ix2, ix3, ix4, ix5, ix6, ix7, ix8)
250
251 CALL sderi3(
252 . gbuf%VOL ,dummy ,geo ,igeo ,
253 . xd1 ,xd2 ,xd3 ,xd4 ,xd5 ,xd6 ,xd7 ,xd8 ,
254 . yd1 ,yd2 ,yd3 ,yd4 ,yd5 ,yd6 ,yd7 ,yd8 ,
255 . zd1 ,zd2 ,zd3 ,zd4 ,zd5 ,zd6 ,zd7 ,zd8 ,
256 . rx ,ry ,rz ,sx ,sy ,sz ,ngl ,pid ,
257 . px1 ,px2 ,px3 ,px4 ,py1 ,py2 ,py3 ,py4 ,
258 . pz1 ,pz2 ,pz3 ,pz4 ,volu ,voldp,nel ,jeul ,
259 . nxref,imulti_fvm )
260 CALL sdlen3(
261 . x1 ,x2 ,x3 ,x4 ,x5 ,x6 ,x7 ,x8 ,
262 . y1 ,y2 ,y3 ,y4 ,y5 ,y6 ,y7 ,y8 ,
263 . z1 ,z2 ,z3 ,z4 ,z5 ,z6 ,z7 ,z8,
264 . deltax, volu)
265
266 gbuf%RHO(:) = zero
267 pm(104,ixs(1, 1 + nft)) = zero !global pressure
268
269 DO ilay = 1, nlay
270 lbuf => elbuf_str%BUFLY(ilay)%LBUF(1,1,1)
271 mbuf => elbuf_str%BUFLY(ilay)%MAT(1,1,1)
272 bufly => elbuf_str%BUFLY(ilay)
273 nuvar = elbuf_str%BUFLY(ilay)%NVAR_MAT
274 l_pla = elbuf_str%BUFLY(ilay)%L_PLA
275 DO i = 1, nel
276 mat(i) = mat_param( ixs(1,i+nft) )%MULTIMAT%MID(ilay)
277C Fill partial volumes
278 lbuf%VOL(i) = mat_param( ixs(1,i+nft) )%MULTIMAT%VFRAC(ilay) * gbuf%VOL(i)
279 ENDDO
280
281! PT=1 => material initializations in matini are not fully compatible with Isolid=12 (lm)
282 ipt=1
283 CALL matini(pm ,ixs ,nixs ,x ,
284 2 geo ,ale_connectivity ,detonators,iparg_gr ,
285 3 sigi ,nel ,skew ,igeo ,
286 4 ipart ,iparts ,
287 5 mat ,ipm ,nsigs ,numsol ,ptsol ,
288 6 ipt ,ngl ,npf ,tf ,bufmat ,
289 7 gbuf ,lbuf ,mbuf ,elbuf_str ,iloadp ,
290 8 facload, deltax ,tempel)
291
292 vfrac = mat_param( ixs(1,1+nft) )%MULTIMAT%VFRAC(ilay)
293 pres = pm(104, mat_param( ixs(1,1+nft) )%MULTIMAT%MID(ilay))
294 pm(104,ixs(1, 1 + nft)) = pm(104,ixs(1, 1 + nft)) + vfrac * pres !global pressure
295
296 matlaw = ipm(2, mat(1))
297 IF (matlaw == 5) THEN
298! JWL material
299 IF (.NOT. error_thrown) THEN
300 IF (pm(44, mat(1)) == zero) THEN
301 CALL ancmsg(msgid = 1623, msgtype = msgerror, anmode = aninfo,
302 . i1 = ipm(1, ixs(1, 1 + nft)), i2 = ipm(1, mat(1)))
303 ENDIF
304 error_thrown = .true.
305 ENDIF
306 CALL m5in3(pm, mat, ipm(1, ixs(1,1+nft)), detonators, lbuf%TB, iparg, x, ixs, nixs)
307 ENDIF
308 IF (matlaw == 6) THEN
309 IF (pm(24, mat(1)) > zero) THEN
310 multi_fvm%NS_DIFF = .true.
311 ENDIF
312 ENDIF
313
314 ENDDO ! ILAY = 1, NLAY
315
316 IF (nlay > 1) THEN
317
318C Mass globalization
319 DO ilay = 1, nlay
320 lbuf => elbuf_str%BUFLY(ilay)%LBUF(1,1,1)
321 DO i = 1, nel
322 gbuf%RHO(i) = gbuf%RHO(i) + lbuf%RHO(i) * mat_param( ixs(1,i+nft) )%MULTIMAT%VFRAC(ilay)
323 ENDDO
324 ENDDO
325
326C Temperature globalization. We must solve later T such as e+p/rho=integral(Cp_global(T),dT)
327 gbuf%TEMP(1:nel)=zero
328 DO ilay = 1, nlay
329 lbuf => elbuf_str%BUFLY(ilay)%LBUF(1,1,1)
330 DO i = 1, nel
331 gbuf%TEMP(i) = gbuf%TEMP(i) + lbuf%TEMP(i) * mat_param( ixs(1,i+nft) )%MULTIMAT%VFRAC(ilay)*lbuf%RHO(i)/gbuf%RHO(i) !volfrac*densfrac=massfrac
332 ENDDO
333 ENDDO
334
335 ENDIF
336C----------------------------------------
337C Thermal and turbulence initialization
338C----------------------------------------
339 IF(jthe /=0) CALL atheri(mat,pm ,gbuf%TEMP)
340 IF(jtur /=0) CALL aturi3(iparg ,gbuf%RHO,pm,ixs,x,
341 . gbuf%RK ,gbuf%RE,volu)
342C----------------------------------------
343C Masses initialization
344C----------------------------------------
345 IF(jlag+jale+jeul/=0) THEN
346 CALL smass3(
347 . gbuf%RHO ,mas ,partsav ,x ,v ,
348 . iparts(nf1),mss(1,nf1) ,volu ,
349 . msnf ,mssf(1,nf1),in ,
350 . vr ,ins(1,nf1) ,wma ,rhocp ,mcp ,
351 . mcps(1,nf1),mssa ,rhof ,alpha ,gbuf%FILL,
352 . ix1, ix2, ix3, ix4, ix5, ix6, ix7, ix8)
353 ENDIF
354C------------------------------------------
355C Element time steps
356C------------------------------------------
357 aire(:) = zero
358 CALL dtmain(geo , pm , ipm , pid , mat , fv ,
359 . gbuf%EINT, gbuf%TEMP, gbuf%DELTAX, gbuf%RK, gbuf%RE, bufmat, deltax, aire, volu, dtx, igeo,igtyp)
360C
361 DO i=1,nel
362 IF(ixs(10,i+nft)/=0) THEN
363 IF( igtyp/=0 .AND.igtyp/=6 .AND. igtyp/=14
364 . .AND.igtyp/=15.AND. igtyp/=29) THEN
365 ipid1=ixs(nixs-1,i+nft)
366 CALL fretitl2(titr1,igeo(npropgi-ltitr+1,ipid1),ltitr)
367 CALL ancmsg(msgid=226,
368 . msgtype=msgerror,
369 . anmode=aninfo_blind_1,
370 . i1=igeo(1,ipid1),
371 . c1=titr1,
372 . i2=igtyp)
373 ENDIF
374 ENDIF
375 dtelem(nft+i)=dtx(i)
376C
377C STI = 0.25 * RHO * VOL / (DT*DT)
378 sti = fourth * gbuf%FILL(i) * gbuf%RHO(i) * volu(i) /
379 . max(em20,dtx(i)*dtx(i))
380 stifn(ixs(2,i+nft))=stifn(ixs(2,i+nft))+sti
381 stifn(ixs(3,i+nft))=stifn(ixs(3,i+nft))+sti
382 stifn(ixs(4,i+nft))=stifn(ixs(4,i+nft))+sti
383 stifn(ixs(5,i+nft))=stifn(ixs(5,i+nft))+sti
384 stifn(ixs(6,i+nft))=stifn(ixs(6,i+nft))+sti
385 stifn(ixs(7,i+nft))=stifn(ixs(7,i+nft))+sti
386 stifn(ixs(8,i+nft))=stifn(ixs(8,i+nft))+sti
387 stifn(ixs(9,i+nft))=stifn(ixs(9,i+nft))+sti
388 ENDDO
389 RETURN
390 END SUBROUTINE multifluid_init3
atheri
subroutine atheri(mat, pm, temp)
Definition atheri.F:41
aturi3
subroutine aturi3(iparg, rho, pm, ix, x, rk, re, volu)
Definition aturi3.F:32
dtmain
subroutine dtmain(geo, pm, ipm, pid, mat, fv, eint, temp, deltax, rk, re, bufmat, ddeltax, aire, vol, dtx, igeo, igtyp)
Definition dtmain.F:67
alpha
#define alpha
Definition eval.h:35
m5in3
subroutine m5in3(pm, mat, m151_id, detonators, tb, iparg, x, ix, nix)
Definition m5in3.F:39
max
#define max(a, b)
Definition macros.h:21
matini
subroutine matini(pm, ix, nix, x, geo, ale_connectivity, detonators, iparg, sigi, nel, skew, igeo, ipart, ipartel, mat, ipm, nsig, nums, pt, ipt, ngl, npf, tf, bufmat, gbuf, lbuf, mbuf, elbuf_str, iloadp, facload, ddeltax, tempel)
Definition matini.F:81
multifluid_init3
subroutine multifluid_init3(elbuf_str, mas, ixs, pm, x, geo, ale_connectivity, iparg_gr, dtelem, sigi, nel, skew, igeo, stifn, partsav, v, iparts, mss, ipart, sigsp, ng, iparg, glob_therm, nsigi, msnf, nvc, mssf, ipm, iuser, nsigs, volnod, bvolnod, vns, bns, in, vr, ins, wma, ptsol, bufmat, mcp, mcps, temp, xrefs, npf, tf, mssa, strsglob, straglob, fail_ini, spbuf, kxsp, ipartsp, nod2sp, sol2sph, irst, iloadp, facload, multi_fvm, error_thrown, detonators, mat_param)
Definition multifluid_init3.F:60
ale_connectivity_mod
Definition ale_connectivity_mod.F:223
detonators_mod
Definition detonators_mod.F:79
message_mod
Definition message_mod.F:1249
names_and_titles_mod
Definition names_and_titles_mod.F:997
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1003
smass3
subroutine smass3(rho, ms, partsav, x, v, ipart, mss, volu, msnf, mssf, in, vr, ins, wma, rhocp, mcp, mcps, mssa, rhof, frac, fill, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8)
Definition smass3.F:44
smorth3
subroutine smorth3(pid, geo, igeo, skew, irep, gama, rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, f1x, f1y, f1z, f2x, f2y, f2z, nsigi, sigsp, nsigs, sigi, ixs, x, jhbe, pt, nel, isolnod)
Definition smorth3.F:43
scoor3
subroutine scoor3(x, xrefs, ixs, geo, mxt, ngeo, ngl, ix1, ix2, ix3, ix4, ix5, ix6, ix7, ix8, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, f1x, f1y, f1z, f2x, f2y, f2z, temp0, temp, nintemp, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8)
Definition scoor3.F:52
sderi3
subroutine sderi3(vol, veul, geo, igeo, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8, jac1, jac2, jac3, jac4, jac5, jac6, ngl, ngeo, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, det, voldp, nel, jeul, nxref, imulti_fvm)
Definition sderi3.F:44
sjac_i
subroutine sjac_i(x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, jac_i, nel)
Definition sderi3.F:264
sdlen3
subroutine sdlen3(x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, deltax, voln)
Definition sdlen3.F:41
srcoor3
subroutine srcoor3(x, xrefs, ixs, geo, mxt, ngeo, ngl, jhbe, ix1, ix2, ix3, ix4, ix5, ix6, ix7, ix8, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, f1x, f1y, f1z, f2x, f2y, f2z, temp0, temp, nintemp, xd1, xd2, xd3, xd4, xd5, xd6, xd7, xd8, yd1, yd2, yd3, yd4, yd5, yd6, yd7, yd8, zd1, zd2, zd3, zd4, zd5, zd6, zd7, zd8)
Definition srcoor3.F:52
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889
fretitl2
subroutine fretitl2(titr, iasc, l)
Definition freform.F:804
sveok3
subroutine sveok3(nvc, nod, ix1, ix2, ix3, ix4, ix5, ix6, ix7, ix8)
Definition sveok3.F:33