OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
multifluid_init2t.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_init2t ../starter/source/multifluid/multifluid_init2t.F
25!||--- called by ------------------------------------------------------
26!|| initia ../starter/source/elements/initia/initia.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| c3veok3 ../starter/source/elements/sh3n/coque3n/c3veok3.F
30!|| m5in2t ../starter/source/initial_conditions/detonation/m5in2t.F
31!|| matini ../starter/source/materials/mat_share/matini.F
32!||--- uses -----------------------------------------------------
33!|| detonators_mod ../starter/share/modules1/detonators_mod.F
34!|| message_mod ../starter/share/message_module/message_mod.F
35!||====================================================================
36 SUBROUTINE multifluid_init2t(ELBUF_STR, NEL, NSIGS, NVC, IPARG, IXTG, ALE_CONNECTIVITY,
37 . IGEO, IPART, IPARTTG, IPM, PTSH3N, NPF, ILOADP,
38 . XGRID, PM, GEO, SIGI, SKEW, TF, BUFMAT, FACLOAD, MULTI_FVM,
39 . ERROR_THROWN, DETONATORS,MAT_PARAM)
40C-----------------------------------------------
41C M o d u l e s
42C-----------------------------------------------
43 USE elbufdef_mod
44 USE multi_fvm_mod
45 USE message_mod
46 USE detonators_mod
47 USE ale_connectivity_mod
48 USE matparam_def_mod, ONLY : matparam_struct_
49 use element_mod , only : nixtg
50C-----------------------------------------------
51C I m p l i c i t T y p e s
52C-----------------------------------------------
53#include "implicit_f.inc"
54C-----------------------------------------------
55C G l o b a l P a r a m e t e r s
56C-----------------------------------------------
57! MVSIZ
58#include "mvsiz_p.inc"
59C-----------------------------------------------
60C C o m m o n B l o c k s
61C-----------------------------------------------
62! NIXTG, NPROPMI, NPROPM, LFACLOAD
63#include "param_c.inc"
64! JEUL, NFT
65#include "vect01_c.inc"
66! NUMSH3N
67#include "scry_c.inc"
68! LIPART1
69#include "scr17_c.inc"
70! SIZLOADP
71#include "com04_c.inc"
72! N2D
73#include "com01_c.inc"
74C-----------------------------------------------
75C D u m m y A r g u m e n t s
76C-----------------------------------------------
77 TYPE(elbuf_struct_), INTENT(IN), TARGET :: ELBUF_STR
78 INTEGER, INTENT(IN) :: NEL, NSIGS, IXTG(NIXTG, *),
79 . IGEO(*), PTSH3N(*), NPF(*), ILOADP(SIZLOADP, *),
80 . ipart(lipart1, *), iparttg(*), ipm(npropmi, *)
81 INTEGER, INTENT(INOUT) :: IPARG(*)
82 INTEGER, INTENT(OUT) :: NVC
83 my_real, INTENT(IN) :: XGRID(3, *), FACLOAD(LFACLOAD, *)
84 my_real,INTENT(INOUT) :: pm(npropm, *)
85 my_real, INTENT(INOUT) :: geo(*), sigi(nsigs, *), skew(lskew, *),
86 . tf(*), bufmat(*)
87 TYPE(multi_fvm_struct) :: MULTI_FVM
88 LOGICAL :: ERROR_THROWN
89 TYPE(DETONATORS_STRUCT_) DETONATORS
90 TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
91 TYPE(matparam_struct_) ,DIMENSION(NUMMAT) ,INTENT(IN) :: MAT_PARAM
92C-----------------------------------------------
93C L o c a l V a r i a b l e s
94C-----------------------------------------------
95 TYPE(l_bufel_) ,POINTER :: LBUF
96 TYPE(G_BUFEL_) ,POINTER :: GBUF
97 TYPE(BUF_MAT_) ,POINTER :: MBUF
98 INTEGER :: NLAY, NGL(MVSIZ), MAT(MVSIZ), IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ)
99 my_real ::
100 . y1(mvsiz), z1(mvsiz),
101 . y2(mvsiz), z2(mvsiz),
102 . y3(mvsiz), z3(mvsiz),
103 . nx(mvsiz), ny(mvsiz), nz(mvsiz),
104 . lgth1(mvsiz), lgth2(mvsiz), lgth3(mvsiz),pres,vfrac
105 my_real :: tempel(nel)
106 INTEGER :: II, I, IP, IBID, ILAY, MATLAW
107C-----------------------------------------------
108C P r e - C o n d i t i o n s
109C-----------------------------------------------
110 IF(n2d>0 .AND. mtn /=151) THEN
111 CALL ancmsg(msgid=1646,
112 . msgtype=msgerror,
113 . anmode=aninfo_blind_1)
114 RETURN
115 ENDIF
116C-----------------------------------------------
117C B e g i n n i n g o f s u b r o u t i n e
118C-----------------------------------------------
119 tempel(:) = zero
120C Global buffer
121 gbuf => elbuf_str%GBUF
122C Number of layers ( = number of material in law 151)
123 nlay = elbuf_str%NLAY
124C Gather nodale coordinates, compute surface (and volume)
125 DO ii = 1, nel
126 i = ii + nft
127 ix1(ii) = ixtg(1 + 1, i)
128 ix2(ii) = ixtg(1 + 2, i)
129 ix3(ii) = ixtg(1 + 3, i)
130 y1(ii) = xgrid(2, ixtg(1 + 1, i))
131 z1(ii) = xgrid(3, ixtg(1 + 1, i))
132 y2(ii) = xgrid(2, ixtg(1 + 2, i))
133 z2(ii) = xgrid(3, ixtg(1 + 2, i))
134 y3(ii) = xgrid(2, ixtg(1 + 3, i))
135 z3(ii) = xgrid(3, ixtg(1 + 3, i))
136 ngl(ii) = ixtg(6, i)
137 nx(ii) = half * ((y2(ii) - y1(ii)) * (z3(ii) - z1(ii)) -
138 . (z2(ii) - z1(ii)) * (y3(ii) - y1(ii)))
139 ny(ii) = zero
140 nz(ii) = zero
141 gbuf%AREA(ii) = sqrt(nx(ii) * nx(ii) + ny(ii) * ny(ii) + nz(ii) * nz(ii))
142 IF (n2d /= 1) THEN
143 gbuf%VOL(ii) = gbuf%AREA(ii)
144 ELSE
145 gbuf%VOL(ii) = (y1(ii) + y2(ii) + y3(ii)) * (
146 . y1(ii) * (z2(ii) - z3(ii)) +
147 . y2(ii) * (z3(ii) - z1(ii)) +
148 . y3(ii) * (z1(ii) - z2(ii))) * one_over_6
149 ENDIF
150 ! TODO(DC) Calcul volume axi
151 lgth1(ii) = sqrt((y2(ii) - y1(ii)) * (y2(ii) - y1(ii)) +
152 . (z2(ii) - z1(ii)) * (z2(ii) - z1(ii)))
153 lgth2(ii) = sqrt((y3(ii) - y2(ii)) * (y3(ii) - y2(ii)) +
154 . (z3(ii) - z2(ii)) * (z3(ii) - z2(ii)))
155 lgth3(ii) = sqrt((y1(ii) - y3(ii)) * (y1(ii) - y3(ii)) +
156 . (z1(ii) - z3(ii)) * (z1(ii) - z3(ii)))
157 gbuf%DELTAX(ii) = gbuf%AREA(ii) / max(lgth1(ii), lgth2(ii), lgth3(ii))
158 ENDDO
159
160 CALL c3veok3(nvc ,ix1 ,ix2 ,ix3 )
161
162 pm(104,ixtg(1, 1 + nft)) = zero !global pressure
163
164C Loop over the materials
165 DO ilay = 1, nlay
166C Layer buffer
167 lbuf => elbuf_str%BUFLY(ilay)%LBUF(1,1,1)
168 mbuf => elbuf_str%BUFLY(ilay)%MAT(1,1,1)
169 DO ii = 1, nel
170C Material
171 mat(ii) = mat_param( ixtg(1,ii+nft) )%MULTIMAT%MID(ilay)
172C Partial volumes
173 lbuf%VOL(ii) = mat_param( ixtg(1,ii+nft) )%MULTIMAT%VFRAC(ilay) * gbuf%VOL(ii)
174 ENDDO
175C Material initialization
176 ip = 1
177 ibid = 0
178 CALL matini(pm, ixtg, nixtg, xgrid,
179 . geo, ale_connectivity, detonators, iparg,
180 . sigi, nel, skew, igeo,
181 . ipart,iparttg,
182 . mat, ipm, nsigs, numsh3n, ptsh3n,
183 . ip, ngl, npf, tf, bufmat,
184 . gbuf, lbuf, mbuf, elbuf_str, iloadp,
185 . facload, gbuf%DELTAX,tempel,mat_param )
186
187 vfrac = mat_param( ixtg(1,1+nft) )%MULTIMAT%VFRAC(ilay)
188 pres = pm(104, mat_param( ixtg(1,1+nft) )%MULTIMAT%MID(ilay) )
189 pm(104,ixtg(1, 1 + nft)) = pm(104,ixtg(1, 1 + nft)) + vfrac * pres !global pressure
190
191 matlaw = ipm(2, mat(1))
192 IF (matlaw == 5) THEN
193! JWL MAT
194 IF (.NOT. error_thrown) THEN
195 IF (pm(44, mat(1)) == zero) THEN
196 CALL ancmsg(msgid = 1623, msgtype = msgerror, anmode = aninfo,
197 . i1 = ipm(1, ixtg(1, 1 + nft)), i2 = ipm(1, mat(1)))
198 ENDIF
199 error_thrown = .true.
200 ENDIF
201 CALL m5in2t(pm, mat, ipm(1, ixtg(1,1+nft)), detonators, lbuf%TB, xgrid, ixtg, nixtg)
202 ENDIF
203 ENDDO
204
205 IF (nlay > 1) THEN
206
207C Mass globalization
208 gbuf%RHO(1:nel)=zero
209 DO ilay = 1, nlay
210 lbuf => elbuf_str%BUFLY(ilay)%LBUF(1,1,1)
211 DO ii = 1, nel
212 gbuf%RHO(ii) = gbuf%RHO(ii) + lbuf%RHO(ii) * mat_param( ixtg(1,ii+nft) )%MULTIMAT%VFRAC(ilay)
213 ENDDO
214 ENDDO
215
216C Temperature globalization. We must solve later T such as e+p/rho=integral(Cp_global(T),dT)
217 gbuf%TEMP(1:nel)=zero
218 DO ilay = 1, nlay
219 lbuf => elbuf_str%BUFLY(ilay)%LBUF(1,1,1)
220 DO ii = 1, nel
221 gbuf%TEMP(ii) = gbuf%TEMP(ii) + lbuf%TEMP(ii) *
222 . mat_param( ixtg(1,ii+nft) )%MULTIMAT%VFRAC(ilay)*lbuf%RHO(ii)/gbuf%RHO(ii) !volfrac*densfrac=massfrac
223 ENDDO
224 ENDDO
225
226 ENDIF
227 END SUBROUTINE multifluid_init2t
c3veok3
subroutine c3veok3(nvc, ix1, ix2, ix3)
Definition c3veok3.F:36
m5in2t
subroutine m5in2t(pm, mat, m151_id, detonators, tb, x, ix, nix)
Definition m5in2t.F:36
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, mat_param)
Definition matini.F:83
multifluid_init2t
subroutine multifluid_init2t(elbuf_str, nel, nsigs, nvc, iparg, ixtg, ale_connectivity, igeo, ipart, iparttg, ipm, ptsh3n, npf, iloadp, xgrid, pm, geo, sigi, skew, tf, bufmat, facload, multi_fvm, error_thrown, detonators, mat_param)
Definition multifluid_init2t.F:40
ale_connectivity_mod
Definition ale_connectivity_mod.F:225
detonators_mod
Definition detonators_mod.F:80
message_mod
Definition message_mod.F:1257
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