OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
multifluid_init3t.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_init3t ../starter/source/multifluid/multifluid_init3t.F
25!||--- called by ------------------------------------------------------
26!|| initia ../starter/source/elements/initia/initia.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| m5in3 ../starter/source/initial_conditions/detonation/m5in3.F
30!|| matini ../starter/source/materials/mat_share/matini.F
31!|| s4coor3 ../starter/source/elements/solid/solide4/s4coor3.F
32!|| s4deri3 ../starter/source/elements/solid/solide4/s4deri3.F
33!|| s4mass3 ../starter/source/elements/solid/solide4/s4mass3.F
34!||--- uses -----------------------------------------------------
35!|| defaults_mod ../starter/source/modules/defaults_mod.F90
36!|| detonators_mod ../starter/share/modules1/detonators_mod.F
37!|| message_mod ../starter/share/message_module/message_mod.F
38!||====================================================================
39 SUBROUTINE multifluid_init3t(ELBUF_STR,
40 . NEL, NSIGS, NSIGI, IXS, IGEO, IPM, IPARG, ALE_CONNECTIVITY, IPARTS, PTSOL, NPF, IPART, ILOADP,
41 . XREFS, GEO, PM, FACLOAD, TF, SKEW, SIGI, BUFMAT, X,
42 . WMA, PARTSAV, MAS, V, MSS, MSSF, MSSA, MSNF, MCPS, ERROR_THROWN, DETONATORS,DEFAULTS,
43 . MAT_PARAM,NINTEMP)
44C-----------------------------------------------
45C M o d u l e s
46C-----------------------------------------------
47 USE elbufdef_mod
48 USE message_mod
49 USE detonators_mod
50 USE ale_connectivity_mod
51 USE defaults_mod
52 USE matparam_def_mod, ONLY : matparam_struct_
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! NIXS, LVEUL
66#include "param_c.inc"
67! NFT
68#include "vect01_c.inc"
69! NUMSOL
70#include "scry_c.inc"
71! LIPART1
72#include "scr17_c.inc"
73! SIZLOADP
74#include "com04_c.inc"
75C-----------------------------------------------
76C D u m m y A r g u m e n t s
77C-----------------------------------------------
78 TYPE(elbuf_struct_), INTENT(IN), TARGET :: ELBUF_STR
79 INTEGER, INTENT(IN) :: NEL, NSIGS, NSIGI, IGEO(NPROPGI, *), IPM(NPROPMI, *),
80 . IPARTS(*), PTSOL(*), NPF(*), IPART(LIPART1, *), ILOADP(SIZLOADP, *)
81 INTEGER, INTENT(INOUT) :: IPARG(*), IXS(NIXS,*)
82 INTEGER, INTENT(IN) :: NINTEMP
83 my_real, INTENT(IN) :: x(3, *), geo(npropg, *),
84 . facload(lfacload, *), tf(*), skew(lskew, *), sigi(nsigi, *), bufmat(*)
85 my_real, INTENT(INOUT) :: xrefs(8, 3, *)
86 my_real,INTENT(INOUT) :: pm(npropm, *)
87 my_real, INTENT(INOUT) :: wma(*), partsav(20, *), mas(*), v(*),
88 . msnf(*), mcps(8, *), mssf(8, *), mss(8, *), mssa(*)
89 LOGICAL :: ERROR_THROWN
90 TYPE(detonators_struct_) DETONATORS
91 TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
92 TYPE(DEFAULTS_), INTENT(IN) :: DEFAULTS
93 TYPE(matparam_struct_) ,DIMENSION(NUMMAT) ,INTENT(IN) :: MAT_PARAM
94C-----------------------------------------------
95C L o c a l V a r i a b l e s
96C-----------------------------------------------
97 TYPE(l_bufel_) ,POINTER :: LBUF
98 TYPE(G_BUFEL_) ,POINTER :: GBUF
99 TYPE(BUF_MAT_) ,POINTER :: MBUF
100 INTEGER :: ILAY, NLAY, PID(MVSIZ), NGL(MVSIZ), MAT(MVSIZ),
101 . ix1(mvsiz), ix2(mvsiz), ix3(mvsiz), ix4(mvsiz)
102 my_real :: x1(mvsiz), y1(mvsiz), z1(mvsiz),
103 . x2(mvsiz), y2(mvsiz), z2(mvsiz),
104 . x3(mvsiz), y3(mvsiz), z3(mvsiz),
105 . x4(mvsiz), y4(mvsiz), z4(mvsiz),
106 . rx(mvsiz),ry(mvsiz),rz(mvsiz),
107 . sx(mvsiz),sy(mvsiz),sz(mvsiz),
108 . tx(mvsiz),ty(mvsiz),tz(mvsiz),
109 . px1(mvsiz),px2(mvsiz),px3(mvsiz),px4(mvsiz),
110 . py1(mvsiz),py2(mvsiz),py3(mvsiz),py4(mvsiz),
111 . pz1(mvsiz),pz2(mvsiz),pz3(mvsiz),pz4(mvsiz),
112 . volu(mvsiz), bid(mvsiz), dummy, pres,vfrac
113 my_real :: tempel(nel)
114 INTEGER :: II, IP, IBID, MATLAW,IMAS_DS
115 DOUBLE PRECISION
116 . VOLDP(MVSIZ)
117C-----------------------------------------------
118C S o u r c e L i n e s
119C-----------------------------------------------
120 gbuf => elbuf_str%GBUF
121C Number of layers ( = number of material in law 151)
122 nlay = elbuf_str%NLAY
123 imas_ds = defaults%SOLID%IMAS
124C Gather coordinates
125 CALL s4coor3(x, xrefs(1, 1, nft + 1), ixs(1, nft + 1), ngl,
126 . mat, pid, ix1, ix2, ix3, ix4,
127 . x1, x2, x3, x4,
128 . y1, y2, y3, y4,
129 . z1, z2, z3, z4)
130C Volume
131 CALL s4deri3(gbuf%VOL, dummy, geo, igeo,
132 . rx, ry, rz, sx, sy, sz, tx, ty, tz,
133 . x1, x2, x3, x4,
134 . y1, y2, y3, y4,
135 . z1, z2, z3, z4,
136 . px1, px2, px3, px4,
137 . py1, py2, py3, py4,
138 . pz1, pz2, pz3, pz4, gbuf%JAC_I,
139 . gbuf%DELTAX, volu, ngl, pid, mat,
140 . pm ,voldp)
141
142 tempel(:) = zero
143 pm(104,ixs(1, 1 + nft)) = zero !global pressure
144
145C Loop over the materials
146 DO ilay = 1, nlay
147C Layer buffer
148 lbuf => elbuf_str%BUFLY(ilay)%LBUF(1,1,1)
149 mbuf => elbuf_str%BUFLY(ilay)%MAT(1,1,1)
150 DO ii = 1, nel
151C Material
152 mat(ii) = mat_param( ixs(1,ii+nft) )%MULTIMAT%MID(ilay)
153C Partial volumes
154 lbuf%VOL(ii) = mat_param( ixs(1,ii+nft) )%MULTIMAT%VFRAC(ilay) * gbuf%VOL(ii)
155 lbuf%VOL0DP(ii) = mat_param( ixs(1,ii+nft) )%MULTIMAT%VFRAC(ilay) * voldp(ii)
156 ENDDO
157C Material initialization
158 ip = 1
159 ibid = 0
160 CALL matini(pm, ixs, nixs, x,
161 . geo, ale_connectivity, detonators, iparg,
162 . sigi, nel, skew, igeo,
163 . ipart,iparts,
164 . mat, ipm, nsigs, numsol, ptsol,
165 . ip, ngl,npf, tf, bufmat,
166 . gbuf, lbuf, mbuf, elbuf_str, iloadp,
167 . facload, gbuf%DELTAX,tempel, mat_param )
168
169 vfrac = mat_param( ixs(1,1+nft) )%MULTIMAT%VFRAC(ilay)
170 pres = pm(104, mat_param( ixs(1,1+nft) )%MULTIMAT%MID(ilay))
171 pm(104,ixs(1, 1 + nft)) = pm(104,ixs(1, 1 + nft)) + vfrac * pres !global pressure
172
173 matlaw = ipm(2, mat(1))
174 IF (matlaw == 5) THEN
175! JWL MAT
176 IF (.NOT. error_thrown) THEN
177 IF (pm(44, mat(1)) == zero) THEN
178 CALL ancmsg(msgid = 1623, msgtype = msgerror, anmode = aninfo,
179 . i1 = ipm(1, ixs(1, 1 + nft)), i2 = ipm(1, mat(1)))
180 ENDIF
181 error_thrown = .true.
182 ENDIF
183 CALL m5in3(pm, mat, ipm(1, ixs(1,1+nft)), detonators, lbuf%TB, iparg, x, ixs, nixs)
184 ENDIF
185 ENDDO
186
187 IF (nlay > 1) THEN
188
189C Mass globalization
190
191 DO ii = 1, nel
192 gbuf%RHO(ii) = zero
193 ENDDO
194 DO ilay = 1, nlay
195 lbuf => elbuf_str%BUFLY(ilay)%LBUF(1,1,1)
196 DO ii = 1, nel
197 gbuf%RHO(ii) = gbuf%RHO(ii) + lbuf%RHO(ii) * mat_param( ixs(1,ii+nft) )%MULTIMAT%VFRAC(ilay)
198 ENDDO
199 ENDDO
200
201C Temperature globalization. We must solve later T such as e+p/rho=integral(Cp_global(T),dT)
202 gbuf%TEMP(1:nel)=zero
203 DO ilay = 1, nlay
204 lbuf => elbuf_str%BUFLY(ilay)%LBUF(1,1,1)
205 DO ii = 1, nel
206 gbuf%TEMP(ii) = gbuf%TEMP(ii) + lbuf%TEMP(ii) *
207 . mat_param( ixs(1,ii+nft) )%MULTIMAT%VFRAC(ilay) *lbuf%RHO(ii)/gbuf%RHO(ii) !volfrac*densfrac=massfrac
208 ENDDO
209 ENDDO
210
211 ENDIF
212 CALL s4mass3(
213 1 gbuf%RHO ,mas ,partsav,x ,v,
214 2 iparts(nft + 1),mss(1,nft + 1),msnf ,mssf(1,nft + 1),wma,
215 3 bid ,bid ,mcps(1,nft + 1),bid,bid ,
216 4 mssa ,ix1 ,ix2 ,ix3 ,ix4 ,
217 5 gbuf%FILL, gbuf%VOL ,imas_ds, nintemp)
218 END SUBROUTINE multifluid_init3t
my_real
#define my_real
Definition cppsort.cpp:32
m5in3
subroutine m5in3(pm, mat, m151_id, detonators, tb, iparg, x, ix, nix)
Definition m5in3.F:39
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, mat_param)
Definition matini.F:83
multifluid_init3t
subroutine multifluid_init3t(elbuf_str, nel, nsigs, nsigi, ixs, igeo, ipm, iparg, ale_connectivity, iparts, ptsol, npf, ipart, iloadp, xrefs, geo, pm, facload, tf, skew, sigi, bufmat, x, wma, partsav, mas, v, mss, mssf, mssa, msnf, mcps, error_thrown, detonators, defaults, mat_param, nintemp)
Definition multifluid_init3t.F:44
ale_connectivity_mod
Definition ale_connectivity_mod.F:225
detonators_mod
Definition detonators_mod.F:80
message_mod
Definition message_mod.F:1257
s4mass3
subroutine s4mass3(rho, ms, partsav, x, v, ipart, mss, msnf, mssf, wma, rhocp, mcp, mcps, temp0, temp, mssa, ix1, ix2, ix3, ix4, fill, volu, imas_ds, nintemp)
Definition s4mass3.F:41
s4coor3
subroutine s4coor3(x, xrefs, ixs, ngl, mxt, ngeo, ix1, ix2, ix3, ix4, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4)
Definition s4coor3.F:40
s4deri3
subroutine s4deri3(vol, veul, geo, igeo, rx, ry, rz, sx, sy, sz, tx, ty, tz, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, jac_i, deltax, det, ngl, ngeo, mxt, pm, voldp)
Definition s4deri3.F:47
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:895
ale_connectivity_mod::t_ale_connectivity
Definition ale_connectivity_mod.F:256
detonators_mod::detonators_struct_
Definition detonators_mod.F:155