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

Go to the source code of this file.

Functions/Subroutines

subroutine qinit2 (elbuf_str, ms, ixq, pm, x, detonators, geo, veul, ale_connectivity, iparg, dtelem, sigi, igeo, nel, skew, msq, ipart, ipartq, ipm, nsigs, wma, ptquad, bufmat, npf, tf, ipargg, iloadp, facload, partsav, v)

Function/Subroutine Documentation

◆ qinit2()

subroutine qinit2 ( type(elbuf_struct_), target elbuf_str,
ms,
integer, dimension(nixq,*) ixq,
pm,
x,
type(detonators_struct_) detonators,
geo,
veul,
type(t_ale_connectivity), intent(inout) ale_connectivity,
integer, dimension(*) iparg,
dtelem,
sigi,
integer, dimension(npropgi,*) igeo,
integer nel,
skew,
msq,
integer, dimension(lipart1,*) ipart,
integer, dimension(*) ipartq,
integer, dimension(npropmi,*) ipm,
integer nsigs,
wma,
integer, dimension(*) ptquad,
bufmat,
integer, dimension(*) npf,
tf,
integer, dimension(*) ipargg,
integer, dimension(sizloadp,*), intent(in) iloadp,
dimension(lfacload,*), intent(in) facload,
partsav,
v )

Definition at line 45 of file qinit2.F.

52C-----------------------------------------------
53C M o d u l e s
54C-----------------------------------------------
55 USE elbufdef_mod
56 USE message_mod
57 USE detonators_mod
58 USE ale_connectivity_mod
59 USE names_and_titles_mod , ONLY : nchartitle
60C-----------------------------------------------
61C I m p l i c i t T y p e s
62C-----------------------------------------------
63#include "implicit_f.inc"
64C-----------------------------------------------
65C G l o b a l P a r a m e t e r s
66C-----------------------------------------------
67#include "mvsiz_p.inc"
68C-----------------------------------------------
69C C o m m o n B l o c k s
70C-----------------------------------------------
71#include "com04_c.inc"
72#include "param_c.inc"
73#include "scry_c.inc"
74#include "vect01_c.inc"
75#include "scr17_c.inc"
76C-----------------------------------------------
77C D u m m y A r g u m e n t s
78C-----------------------------------------------
79 INTEGER IXQ(NIXQ,*), IPARG(*),IGEO(NPROPGI,*),
80 . NEL,IPART(LIPART1,*),IPARTQ(*),IPM(NPROPMI,*), PTQUAD(*),
81 . NSIGS, NPF(*),IPARGG(*)
82 my_real
83 . ms(*), pm(npropm,*), x(*), geo(npropg,*),
84 . veul(10,*), dtelem(*),sigi(nsigs,*),skew(lskew,*),
85 . msq(*), bufmat(*), tf(*),wma(*),partsav(20,*),v(*)
86 TYPE(ELBUF_STRUCT_), TARGET :: ELBUF_STR
87 INTEGER,INTENT(IN) :: ILOADP(SIZLOADP,*)
88 my_real,INTENT(IN) :: facload(lfacload,*)
89 TYPE(DETONATORS_STRUCT_)::DETONATORS
90 TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
91C-----------------------------------------------
92C L o c a l V a r i a b l e s
93C-----------------------------------------------
94 INTEGER MAT(MVSIZ), PID(MVSIZ), NGL(MVSIZ)
95 INTEGER NF1, I, IGTYP, IP, IBID, IPID1
96 my_real e1y(mvsiz),e1z(mvsiz),e2y(mvsiz),e2z(mvsiz),
97 . bid(1), dtx(mvsiz),
98 . sy(mvsiz) ,sz(mvsiz) ,ty(mvsiz) ,tz(mvsiz)
99 my_real fv
100 CHARACTER(LEN=NCHARTITLE)::TITR
101 TYPE(G_BUFEL_) ,POINTER :: GBUF
102 TYPE(L_BUFEL_), POINTER :: LBUF
103 TYPE(BUF_MAT_) ,POINTER :: MBUF
104 my_real y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
105 . z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
106 . aire(mvsiz), deltax(mvsiz)
107 my_real :: tempel(nel)
108 INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ)
109C-----------------------------------------------
110C S o u r c e L i n e s
111C=======================================================================
112 gbuf => elbuf_str%GBUF
113 lbuf => elbuf_str%BUFLY(1)%LBUF(1,1,1)
114 mbuf => elbuf_str%BUFLY(1)%MAT(1,1,1)
115c
116 igtyp = iparg(38)
117 jcvt = iparg(37)
118 ibid = 0
119 bid = zero
120 tempel(:) = zero
121 nf1 = nft+1
122c
123 IF(jcvt == 0)THEN
124 CALL qcoor2(x,ixq(1,nf1),ngl,mat,
125 . pid, ix1, ix2, ix3, ix4,
126 . y1, y2, y3, y4,
127 . z1, z2, z3, z4,
128 . sy, sz, ty, tz)
129 ELSE
130 CALL qrcoor2(x ,ixq(1,nf1),ngl ,mat ,
131 . pid, ix1, ix2, ix3, ix4,
132 . y1, y2, y3, y4,
133 . z1, z2, z3, z4,
134 . sy, sz, ty, tz,
135 . e1y, e1z, e2y, e2z)
136 END IF
137 IF (igtyp == 6) CALL qmorth2(pid ,geo ,igeo ,gbuf%GAMA, nel,
138 . sy ,sz ,ty ,tz ,
139 . e1y ,e1z , e2y, e2z)
140 CALL qvoli2(gbuf%VOL,ixq(1,nf1),
141 . ngl, aire,
142 . y1, y2, y3, y4,
143 . z1, z2, z3 ,z4)
144 CALL qdlen2(iparg(63),
145 . aire, deltax,
146 . y1, y2, y3, y4,
147 . z1, z2, z3, z4)
148 IF(jeul/=0) CALL edlen2(veul(1,nf1), aire, deltax)
149C
150 ip=0
151 CALL matini(pm ,ixq ,nixq ,x ,
152 . geo ,ale_connectivity ,detonators,iparg ,
153 . sigi ,nel ,skew ,igeo ,
154 . ipart ,ipartq ,
155 . mat ,ipm ,nsigs ,numquad ,ptquad ,
156 . ip ,ngl ,npf ,tf ,bufmat ,
157 . gbuf ,lbuf ,mbuf ,elbuf_str ,iloadp ,
158 . facload, deltax ,tempel )
159C----------------------------------------
160C INITIALISATION DE LA THERMIQUE ET TURBULENCE
161C----------------------------------------
162 IF (jthe/=0) CALL atheri(mat,pm,gbuf%TEMP)
163 IF (jtur/=0) CALL aturi2(ipargg ,gbuf%RHO,pm,ixq,x,
164 . gbuf%RK,gbuf%RE, aire)
165C------------------------------------------
166C INITIALISATION DE LA MATRICE DE MASSE
167C------------------------------------------
168 IF (jlag+jale+jeul/=0)
169 . CALL qmasi2(pm,mat,ms,gbuf%VOL,msq(nf1),wma,ipartq(nft+1),partsav,
170 . ix1, ix2, ix3, ix4 ,x ,v)
171C-------------------------------------------
172C CALCUL DES DT ELEMENTAIRES
173C-------------------------------------------
174 CALL dtmain(geo ,pm ,ipm ,pid ,mat ,fv ,
175 . gbuf%EINT ,gbuf%TEMP ,gbuf%DELTAX ,gbuf%RK ,gbuf%RE ,bufmat, deltax, aire,
176 . gbuf%VOL, dtx, igeo,igtyp)
177c
178 DO i=lft,llt
179 IF(ixq(6,i+nft)/=0) THEN
180 IF (igtyp/=0 .AND. igtyp/=6 .AND.
181 . igtyp/=14.AND.igtyp/=15)THEN
182 ipid1=ixq(nixq-1,i+nft)
183 CALL fretitl2(titr,igeo(npropgi-ltitr+1,ipid1),ltitr)
184 CALL ancmsg(msgid=226,
185 . msgtype=msgerror,
186 . anmode=aninfo_blind_1,
187 . i1=igeo(1,ipid1),
188 . c1=titr,
189 . i2=igtyp)
190 ENDIF
191 ENDIF
192 dtelem(nft+i)=dtx(i)
193 ENDDO
194C-----------
195 RETURN
atheri
subroutine atheri(mat, pm, temp)
Definition atheri.F:41
aturi2
subroutine aturi2(iparg, rho, pm, ix, x, rk, re, aire)
Definition aturi2.F:32
my_real
#define my_real
Definition cppsort.cpp: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
edlen2
subroutine edlen2(veul, aire, deltax)
Definition edlen2.F:31
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
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
qmasi2
subroutine qmasi2(pm, mat, ms, vol, msq, wma, ipart, partsav, ix1, ix2, ix3, ix4, x, v)
Definition qmasi2.F:33
qmorth2
subroutine qmorth2(pid, geo, igeo, gama, nel, ry, rz, sy, sz, e1y, e1z, e2y, e2z)
Definition qmorth2.F:37
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
qrcoor2
subroutine qrcoor2(x, ixq, ngl, mxt, pid, ix1, ix2, ix3, ix4, y1, y2, y3, y4, z1, z2, z3, z4, sy, sz, ty, tz, e1y, e1z, e2y, e2z)
Definition qrcoor2.F:34
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