OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
multifluid_init2.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "vect01_c.inc"
#include "scr17_c.inc"
#include "scry_c.inc"
#include "com04_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine multifluid_init2 (nel, nsigs, iparg, ixq, ipm, ale_connectivity, igeo, ipart, ipartq, npf, ptquad, iloadp, x, pm, geo, sigi, skew, tf, bufmat, facload, elbuf_str, error_thrown, detonators, mat_param)

Function/Subroutine Documentation

◆ multifluid_init2()

subroutine multifluid_init2 ( integer, intent(in) nel,
integer, intent(in) nsigs,
integer, dimension(*), intent(inout) iparg,
integer, dimension(nixq, *), intent(in) ixq,
integer, dimension(npropmi, nummat), intent(in) ipm,
type(t_ale_connectivity), intent(inout) ale_connectivity,
integer, dimension(*), intent(in) igeo,
integer, dimension(lipart1, *), intent(in) ipart,
integer, dimension(*), intent(in) ipartq,
integer, dimension(*), intent(in) npf,
integer, dimension(*), intent(in) ptquad,
integer, dimension(sizloadp, *), intent(in) iloadp,
dimension(*), intent(in) x,
dimension(npropm, nummat), intent(inout) pm,
dimension(*), intent(inout) geo,
dimension(nsigs, *), intent(inout) sigi,
dimension(lskew, *), intent(inout) skew,
dimension(*), intent(inout) tf,
dimension(*), intent(inout) bufmat,
dimension(lfacload, *), intent(in) facload,
type(elbuf_struct_), intent(in), target elbuf_str,
logical error_thrown,
type(detonators_struct_) detonators,
type(matparam_struct_), dimension(nummat), intent(in) mat_param )

Definition at line 38 of file multifluid_init2.F.

43C-----------------------------------------------
44C M o d u l e s
45C-----------------------------------------------
46 USE elbufdef_mod
47 USE message_mod
48 USE detonators_mod
49 USE ale_connectivity_mod
50 USE matparam_def_mod, ONLY : matparam_struct_
51C-----------------------------------------------
52C I m p l i c i t T y p e s
53C-----------------------------------------------
54#include "implicit_f.inc"
55C-----------------------------------------------
56C G l o b a l P a r a m e t e r s
57C-----------------------------------------------
58! MVSIZ
59#include "mvsiz_p.inc"
60C-----------------------------------------------
61C C o m m o n B l o c k s
62C-----------------------------------------------
63! NIXQ, NPROPMI, NPROPM, LFACLOAD
64#include "param_c.inc"
65! JEUL, NFT
66#include "vect01_c.inc"
67! LIPART1
68#include "scr17_c.inc"
69! NUMQUAD
70#include "scry_c.inc"
71! SIZLOADP
72#include "com04_c.inc"
73C-----------------------------------------------
74C D u m m y A r g u m e n t s
75C-----------------------------------------------
76 INTEGER, INTENT(IN) :: NEL, NSIGS,IXQ(NIXQ, *), IPM(NPROPMI, NUMMAT),
77 . IGEO(*), IPART(LIPART1, *), IPARTQ(*), PTQUAD(*),
78 . NPF(*), ILOADP(SIZLOADP, *)
79 INTEGER, INTENT(INOUT) :: IPARG(*)
80 TYPE(ELBUF_STRUCT_), INTENT(IN), TARGET :: ELBUF_STR
81 my_real, INTENT(IN) :: x(*), facload(lfacload, *)
82 my_real, INTENT(INOUT) :: pm(npropm, nummat)
83 my_real, INTENT(INOUT) :: geo(*), sigi(nsigs, *),
84 . skew(lskew, *), tf(*), bufmat(*)
85 LOGICAL :: ERROR_THROWN
86 TYPE(DETONATORS_STRUCT_) :: DETONATORS
87 TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
88 TYPE(MATPARAM_STRUCT_) ,DIMENSION(NUMMAT) ,INTENT(IN) :: MAT_PARAM
89C-----------------------------------------------
90C L o c a l V a r i a b l e s
91C-----------------------------------------------
92 INTEGER :: ILAY, NLAY, II, IP, IBID, MATLAW
93 INTEGER :: NGL(MVSIZ), MAT(MVSIZ), PID(MVSIZ)
94 INTEGER :: IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ)
95 my_real :: y1(mvsiz), y2(mvsiz), y3(mvsiz),y4(mvsiz),
96 . z1(mvsiz), z2(mvsiz), z3(mvsiz),z4(mvsiz),
97 . sy(mvsiz), sz(mvsiz), ty(mvsiz),tz(mvsiz),pres,vfrac
98 my_real :: tempel(nel)
99 TYPE(L_BUFEL_) ,POINTER :: LBUF
100 TYPE(G_BUFEL_) ,POINTER :: GBUF
101 TYPE(BUF_MAT_) ,POINTER :: MBUF
102C-----------------------------------------------
103C B e g i n n i n g o f s u b r o u t i n e
104C-----------------------------------------------
105C Global buffer
106 gbuf => elbuf_str%GBUF
107C Number of layers ( = number of material in law 151)
108 nlay = elbuf_str%NLAY
109C Gather coordinates, material Id and so on
110 CALL qcoor2(x, ixq(1, nft + 1), ngl, mat, pid,
111 . ix1, ix2, ix3, ix4,
112 . y1, y2, y3, y4,
113 . z1, z2, z3, z4,
114 . sy, sz, ty, tz)
115C Volume, area computation
116 CALL qvoli2(gbuf%VOL, ixq(1, nft + 1), ngl, gbuf%AREA,
117 . y1, y2, y3, y4,
118 . z1, z2, z3, z4)
119C Compute deltax
120 IF (jeul /= 0) THEN
121 CALL qdlen2(iparg(63),
122 . gbuf%AREA, gbuf%DELTAX,
123 . y1, y2, y3, y4,
124 . z1, z2, z3, z4)
125 ENDIF
126 tempel(:) = zero
127 pm(104,ixq(1, 1 + nft)) = zero !global pressure
128
129C Loop over the materials
130 DO ilay = 1, nlay
131C Layer buffer
132 lbuf => elbuf_str%BUFLY(ilay)%LBUF(1,1,1)
133 mbuf => elbuf_str%BUFLY(ilay)%MAT(1,1,1)
134 DO ii = 1, nel
135C Material
136 mat(ii) = mat_param( ixq(1,ii+nft) )%MULTIMAT%MID(ilay)
137C Partial volumes
138 lbuf%VOL(ii) =mat_param( ixq(1,ii+nft) )%MULTIMAT%VFRAC(ilay) * gbuf%VOL(ii)
139 ENDDO
140C Material initialization
141 ip = 1
142 ibid = 0
143 CALL matini(pm, ixq, nixq, x,
144 . geo, ale_connectivity, detonators, iparg,
145 . sigi, nel, skew, igeo,
146 . ipart,ipartq,
147 . mat, ipm, nsigs, numquad, ptquad,
148 . ip, ngl,npf, tf, bufmat,
149 . gbuf, lbuf, mbuf, elbuf_str, iloadp,
150 . facload, gbuf%DELTAX,tempel)
151
152 vfrac = mat_param( ixq(1,1+nft) )%MULTIMAT%VFRAC(ilay)
153 pres = pm(104, mat_param( ixq(1,1+nft) )%MULTIMAT%MID(ilay) )
154 pm(104,ixq(1, 1 + nft)) = pm(104,ixq(1, 1 + nft)) + vfrac * pres !global pressure
155
156 matlaw = ipm(2, mat(1))
157 IF (matlaw == 5) THEN
158! JWL MAT - layer detonation times
159 IF (.NOT. error_thrown) THEN
160 IF (pm(44, mat(1)) == zero) THEN
161 CALL ancmsg(msgid = 1623, msgtype = msgerror, anmode = aninfo,
162 . i1 = ipm(1, ixq(1, 1 + nft)), i2 = ipm(1, mat(1)))
163 ENDIF
164 error_thrown = .true.
165 ENDIF
166 CALL m5in2(pm, mat, ipm(1, ixq(1,1+nft)), detonators, lbuf%TB, x, ixq, nixq)
167 ENDIF
168 ENDDO
169
170 IF (nlay > 1) THEN
171
172C Mass globalization
173 gbuf%RHO(1:nel)=zero
174 DO ilay = 1, nlay
175 lbuf => elbuf_str%BUFLY(ilay)%LBUF(1,1,1)
176 DO ii = lft, llt
177 gbuf%RHO(ii) = gbuf%RHO(ii) + lbuf%RHO(ii) * mat_param( ixq(1,ii+nft) )%MULTIMAT%VFRAC(ilay)
178 ENDDO
179 ENDDO
180
181C Temperature globalization. We must solve later T such as e+p/rho=integral(Cp_global(T),dT)
182 gbuf%TEMP(1:nel)=zero
183 DO ilay = 1, nlay
184 lbuf => elbuf_str%BUFLY(ilay)%LBUF(1,1,1)
185 DO ii = 1, nel
186 gbuf%TEMP(ii) = gbuf%TEMP(ii) + lbuf%TEMP(ii) * mat_param( ixq(1,ii+nft) )%MULTIMAT%VFRAC(ilay) *
187 . lbuf%RHO(ii)/gbuf%RHO(ii) !volfrac*densfrac=massfrac
188 ENDDO
189 ENDDO
190
191 ENDIF
192
my_real
#define my_real
Definition cppsort.cpp:32
m5in2
subroutine m5in2(pm, mat, m151_id, detonators, tb, x, ix, nix)
Definition m5in2.F:40
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
ale_connectivity_mod
Definition ale_connectivity_mod.F:223
detonators_mod
Definition detonators_mod.F:79
message_mod
Definition message_mod.F:1249
qvoli2
subroutine qvoli2(volu, ixq, ngl, aire, y1, y2, y3, y4, z1, z2, z3, z4)
Definition qvoli2.F:40
qcoor2
subroutine qcoor2(x, ixq, ngl, mxt, pid, ix1, ix2, ix3, ix4, y1, y2, y3, y4, z1, z2, z3, z4, sy, sz, ty, tz)
Definition qcoor2.F:37
qdlen2
subroutine qdlen2(iparg, aire, deltax, y1, y2, y3, y4, z1, z2, z3, z4)
Definition qdlen2.F:39
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