OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s16init3.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!|| s16init3 ../starter/source/elements/thickshell/solide16/s16init3.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!|| dtmain ../starter/source/materials/time_step/dtmain.F
31!|| failini ../starter/source/elements/solid/solide/failini.F
32!|| fretitl2 ../starter/source/starter/freform.f
33!|| matini ../starter/source/materials/mat_share/matini.F
34!|| s16coor3 ../starter/source/elements/thickshell/solide16/s16coor3.F
35!|| s16deri3 ../starter/source/elements/thickshell/solide16/s16deri3.F
36!|| s16mass3 ../starter/source/elements/thickshell/solide16/s16mass3.F
37!|| s16rst ../starter/source/elements/thickshell/solide16/s16deri3.F
38!|| s20msi ../starter/source/elements/solid/solide20/s20mass3.F
39!|| s20temp ../starter/source/elements/solid/solide20/s20temp.F90
40!|| sbulk3 ../starter/source/elements/solid/solide/sbulk3.F
41!|| sigin20b ../starter/source/elements/solid/solide20/s20mass3.F
42!||--- uses -----------------------------------------------------
43!|| detonators_mod ../starter/share/modules1/detonators_mod.F
44!|| message_mod ../starter/share/message_module/message_mod.F
45!|| s20temp_mod ../starter/source/elements/solid/solide20/s20temp.F90
46!||====================================================================
47 SUBROUTINE s16init3(
48 1 ELBUF_STR ,MAS ,IXS ,PM ,X ,
49 2 DETONATORS,GEO ,VEUL ,ALE_CONNECTIVITY ,IPARG ,
50 3 DTELEM ,SIGI ,NEL ,SKEW ,IGEO ,
51 4 STIFN ,PARTSAV ,V ,IPARTS ,MSS ,
52 5 IXS16 ,IPART ,MSSX ,SIGSP ,NSIGI ,
53 6 IPM ,IUSER ,NSIGS ,VOLNOD ,BVOLNOD ,
54 7 VNS ,BNS ,VNSX ,BNSX ,PTSOL ,
55 8 BUFMAT ,MCP ,MCPS ,MCPSX ,TEMP ,
56 9 NPF ,TF ,STRSGLOB,STRAGLOB,FAIL_INI ,
57 A ILOADP ,FACLOAD ,RNOISE ,PERTURB ,MAT_PARAM,
58 B GLOB_THERM)
59C-----------------------------------------------
60C M o d u l e s
61C-----------------------------------------------
62 USE elbufdef_mod
63 USE message_mod
64 USE detonators_mod
65 USE ale_connectivity_mod
66 USE matparam_def_mod
67 USE names_and_titles_mod , ONLY : nchartitle
68 use glob_therm_mod
69 use s20temp_mod
70C-----------------------------------------------
71C I m p l i c i t T y p e s
72C-----------------------------------------------
73#include "implicit_f.inc"
74C-----------------------------------------------
75C G l o b a l P a r a m e t e r s
76C-----------------------------------------------
77#include "mvsiz_p.inc"
78C-----------------------------------------------
79C C o m m o n B l o c k s
80C-----------------------------------------------
81#include "com04_c.inc"
82#include "param_c.inc"
83#include "scr12_c.inc"
84#include "scr17_c.inc"
85#include "scry_c.inc"
86#include "vect01_c.inc"
87C-----------------------------------------------
88C D u m m y A r g u m e n t s
89C-----------------------------------------------
90 INTEGER IXS(NIXS,*), IPARG(*),IPARTS(*),IGEO(NPROPGI,*),
91 . IXS16(8,*), IPART(LIPART1,*),IPM(NPROPMI,*), PTSOL(*),
92 . NPF(*),STRSGLOB(*),STRAGLOB(*),FAIL_INI(*),PERTURB(NPERTURB)
93 INTEGER NEL,NSIGI,IUSER,NSIGS
94 my_real
95 . MAS(*), PM(NPROPM,*), X(*), GEO(NPROPG,*),
96 . VEUL(LVEUL,*), DTELEM(*),SIGI(NSIGS,*),SKEW(LSKEW,*),STIFN(*),
97 . PARTSAV(20,*), V(*), MSS(8,*), MSSX(12,*), SIGSP(NSIGI, *),
98 . VOLNOD(*), BVOLNOD(*), VNS(8,*), BNS(8,*),
99 . VNSX(12,*), BNSX(12,*),BUFMAT(*),RNOISE(NPERTURB,*),
100 . MCP(*), MCPS(8,*),MCPSX(12,*), TEMP(*), TF(*)
101 TYPE(elbuf_struct_), TARGET :: ELBUF_STR
102 INTEGER,INTENT(IN) :: ILOADP(SIZLOADP,*)
103 my_real,INTENT(IN) :: FACLOAD(LFACLOAD,*)
104 TYPE(detonators_struct_)::DETONATORS
105 TYPE(t_ale_connectivity), INTENT(INOUT) :: ALE_CONNECTIVITY
106 TYPE (MATPARAM_STRUCT_) ,DIMENSION(NUMMAT) ,INTENT(INOUT) :: MAT_PARAM
107 type (glob_therm_) ,intent(in) :: glob_therm
108C-----------------------------------------------
109C L o c a l V a r i a b l e s
110C-----------------------------------------------
111 INTEGER NF1,IBID,I,IGTYP,IP,NF2,NPTR,NPTS,NPTT,NLAY,IL,IR,IS,IT,
112 . n, nuvar,iint, ncc,idef,jhbe,ipid1,l_pla,l_sigb
113 INTEGER NC(MVSIZ,16),MAT(MVSIZ), PID(MVSIZ), NGL(MVSIZ),RBID(1)
114 INTEGER ,PARAMETER :: NPE=16
115 CHARACTER(LEN=NCHARTITLE)::TITR1
116 my_real
117 . bid, fv,wi,aa,bb
118 my_real
119 . mass(mvsiz),
120 . volp(mvsiz,5), sti(mvsiz),deltax(mvsiz),deltax2(mvsiz),
121 . xx(mvsiz,16), yy(mvsiz,16), zz(mvsiz,16),
122 . vx(mvsiz,16), vy(mvsiz,16), vz(mvsiz,16),
123 . px(mvsiz,16), py(mvsiz,16), pz(mvsiz,16),
124 . rx(mvsiz),ry(mvsiz),rz(mvsiz),
125 . sx(mvsiz),sy(mvsiz),sz(mvsiz),volg(mvsiz),
126 . tx(mvsiz),ty(mvsiz),tz(mvsiz),ul(mvsiz,16),
127 . ni(mvsiz,16),dnidr(mvsiz,16),dnids(mvsiz,16),dnidt(mvsiz,16),
128 . dtx(mvsiz),stin(mvsiz,16), rhocp(mvsiz),temp0(mvsiz), aire(mvsiz)
129 my_real :: tempel(nel)
130 TYPE(l_bufel_) ,POINTER :: LBUF
131 TYPE(G_BUFEL_) ,POINTER :: GBUF
132 TYPE(BUF_MAT_) ,POINTER :: MBUF
133C-----------------------------------------------
134 my_real
135 . W_GAUSS(9,9),A_GAUSS(9,9),W_LOBATTO(9,9),A_LOBATTO(9,9),
136 . W_NEWTON(9,9),A_NEWTON(9,9)
137C-----------------------------------------------
138 DATA W_GAUSS /
139 1 2. ,0. ,0. ,
140 1 0. ,0. ,0. ,
141 1 0. ,0. ,0. ,
142 2 1. ,1. ,0. ,
143 2 0. ,0. ,0. ,
144 2 0. ,0. ,0. ,
145 3 0.555555555555556,0.888888888888889,0.555555555555556,
146 3 0. ,0. ,0. ,
147 3 0. ,0. ,0. ,
148 4 0.347854845137454,0.652145154862546,0.652145154862546,
149 4 0.347854845137454,0. ,0. ,
150 4 0. ,0. ,0. ,
151 5 0.236926885056189,0.478628670499366,0.568888888888889,
152 5 0.478628670499366,0.236926885056189,0. ,
153 5 0. ,0. ,0. ,
154 6 0.171324492379170,0.360761573048139,0.467913934572691,
155 6 0.467913934572691,0.360761573048139,0.171324492379170,
156 6 0. ,0. ,0. ,
157 7 0.129484966168870,0.279705391489277,0.381830050505119,
158 7 0.417959183673469,0.381830050505119,0.279705391489277,
159 7 0.129484966168870,0. ,0. ,
160 8 0.101228536290376,0.222381034453374,0.313706645877887,
161 8 0.362683783378362,0.362683783378362,0.313706645877887,
162 8 0.222381034453374,0.101228536290376,0. ,
163 9 0.081274388361574,0.180648160694857,0.260610696402935,
164 9 0.312347077040003,0.330239355001260,0.312347077040003,
165 9 0.260610696402935,0.180648160694857,0.081274388361574/
166 DATA a_gauss /
167 1 0. ,0. ,0. ,
168 1 0. ,0. ,0. ,
169 1 0. ,0. ,0. ,
170 2 -.577350269189626,0.577350269189626,0. ,
171 2 0. ,0. ,0. ,
172 2 0. ,0. ,0. ,
173 3 -.774596669241483,0. ,0.774596669241483,
174 3 0. ,0. ,0. ,
175 3 0. ,0. ,0. ,
176 4 -.861136311594053,-.339981043584856,0.339981043584856,
177 4 0.861136311594053,0. ,0. ,
178 4 0. ,0. ,0. ,
179 5 -.906179845938664,-.538469310105683,0. ,
180 5 0.538469310105683,0.906179845938664,0. ,
181 5 0. ,0. ,0. ,
182 6 -.932469514203152,-.661209386466265,-.238619186083197,
183 6 0.238619186083197,0.661209386466265,0.932469514203152,
184 6 0. ,0. ,0. ,
185 7 -.949107912342759,-.741531185599394,-.405845151377397,
186 7 0. ,0.405845151377397,0.741531185599394,
187 7 0.949107912342759,0. ,0. ,
188 8 -.960289856497536,-.796666477413627,-.525532409916329,
189 8 -.183434642495650,0.183434642495650,0.525532409916329,
190 8 0.796666477413627,0.960289856497536,0. ,
191 9 -.968160239507626,-.836031107326636,-.613371432700590,
192 9 -.324253423403809,0. ,0.324253423403809,
193 9 0.613371432700590,0.836031107326636,0.968160239507626/
194C-----
195 DATA w_lobatto /
196 1 2. ,0. ,0. ,
197 1 0. ,0. ,0. ,
198 1 0. ,0. ,0. ,
199 2 1. ,1. ,0. ,
200 2 0. ,0. ,0. ,
201 2 0. ,0. ,0. ,
202 3 0.333333333333333,1.333333333333333,0.333333333333333,
203 3 0. ,0. ,0. ,
204 3 0. ,0. ,0. ,
205 4 0.166666666666667,0.833333333333333,0.833333333333333,
206 4 0.166666666666667,0. ,0. ,
207 4 0. ,0. ,0. ,
208 5 0.1 ,0.544444444444444,0.711111111111111,
209 5 0.544444444444444,0.1 ,0. ,
210 5 0. ,0. ,0. ,
211 6 0.066666666666667,0.37847496 ,0.55485838 ,
212 6 0.55485838 ,0.37847496 ,0.066666666666667,
213 6 0. ,0. ,0. ,
214 7 0.04761904 ,0.27682604 ,0.43174538 ,
215 7 0.48761904 ,0.43174538 ,0.27682604 ,
216 7 0.04761904 ,0. ,0. ,
217 8 0.03571428 ,0.21070422 ,0.34112270 ,
218 8 0.41245880 ,0.41245880 ,0.34112270 ,
219 8 0.21070422 ,0.03571428 ,0. ,
220 9 0.027777777777778,0.1654953616 ,0.2745387126 ,
221 9 0.3464285110 ,0.3715192744 ,0.3464285110 ,
222 9 0.2745387126 ,0.1654953616 ,0.027777777777778/
223 DATA a_lobatto /
224 1 0. ,0. ,0. ,
225 1 0. ,0. ,0. ,
226 1 0. ,0. ,0. ,
227 2 -1. ,1. ,0. ,
228 2 0. ,0. ,0. ,
229 2 0. ,0. ,0. ,
230 3 -1. ,0. ,1. ,
231 3 0. ,0. ,0. ,
232 3 0. ,0. ,0. ,
233 4 -1. ,-.44721360 ,0.44721360 ,
234 4 1. ,0. ,0. ,
235 4 0. ,0. ,0. ,
236 5 -1. ,-.65465367 ,0. ,
237 5 0.65465367 , 1. ,0. ,
238 5 0. ,0. ,0. ,
239 6 -1. ,-.76505532 ,-.28523152 ,
240 6 0.28523152 ,0.76505532 , 1. ,
241 6 0. ,0. ,0. ,
242 7 -1. ,-.83022390 ,-.46884879 ,
243 7 0. ,0.46884879 ,0.83022390 ,
244 7 1. ,0. ,0. ,
245 8 -1. ,-.87174015 ,-.59170018 ,
246 8 -.20929922 ,0.20929922 ,0.59170018 ,
247 8 0.87174015 , 1. ,0. ,
248 9 -1. ,-.8997579954 ,-.6771862795 ,
249 9 -.3631174638 ,0. ,0.3631174638 ,
250 9 0.6771862795 ,0.8997579954 , 1. /
251C-----
252C-----
253 DATA w_newton /
254 1 2. ,0. ,0. ,
255 1 0. ,0. ,0. ,
256 1 0. ,0. ,0. ,
257 2 1. ,1. ,0. ,
258 2 0. ,0. ,0. ,
259 2 0. ,0. ,0. ,
260 3 0.5 ,1. ,0.5 ,
261 3 0. ,0. ,0. ,
262 3 0. ,0. ,0. ,
263 4 0.166666666666667,0.833333333333333,0.833333333333333,
264 4 0.166666666666667,0. ,0. ,
265 4 0. ,0. ,0. ,
266 5 0.25 ,0.5 ,0.5 ,
267 5 0.5 ,0.25 ,0. ,
268 5 0. ,0. ,0. ,
269 6 0.066666666666667,0.37847496 ,0.55485838 ,
270 6 0.55485838 ,0.37847496 ,0.066666666666667,
271 6 0. ,0. ,0. ,
272 7 0.04761904 ,0.27682604 ,0.43174538 ,
273 7 0.48761904 ,0.43174538 ,0.27682604 ,
274 7 0.04761904 ,0. ,0. ,
275 8 0.03571428 ,0.21070422 ,0.34112270 ,
276 8 0.41245880 ,0.41245880 ,0.34112270 ,
277 8 0.21070422 ,0.03571428 ,0. ,
278 9 0.027777777777778,0.1654953616 ,0.2745387126 ,
279 9 0.3464285110 ,0.3715192744 ,0.3464285110 ,
280 9 0.2745387126 ,0.1654953616 ,0.027777777777778/
281 DATA a_newton /
282 1 0. ,0. ,0. ,
283 1 0. ,0. ,0. ,
284 1 0. ,0. ,0. ,
285 2 -1. ,1. ,0. ,
286 2 0. ,0. ,0. ,
287 2 0. ,0. ,0. ,
288 3 -1. ,0. ,1. ,
289 3 0. ,0. ,0. ,
290 3 0. ,0. ,0. ,
291 4 -1. ,-.44721360 ,0.44721360 ,
292 4 1. ,0. ,0. ,
293 4 0. ,0. ,0. ,
294 5 -1. ,-.5 ,0. ,
295 5 0.5 , 1. ,0. ,
296 5 0. ,0. ,0. ,
297 6 -1. ,-.76505532 ,-.28523152 ,
298 6 0.28523152 ,0.76505532 , 1. ,
299 6 0. ,0. ,0. ,
300 7 -1. ,-.83022390 ,-.46884879 ,
301 7 0. ,0.46884879 ,0.83022390 ,
302 7 1. ,0. ,0. ,
303 8 -1. ,-.87174015 ,-.59170018 ,
304 8 -.20929922 ,0.20929922 ,0.59170018 ,
305 8 0.87174015 , 1. ,0. ,
306 9 -1. ,-.8997579954 ,-.6771862795 ,
307 9 -.3631174638 ,0. ,0.3631174638 ,
308 9 0.6771862795 ,0.8997579954 , 1. /
309C-----------------------------------------------
310C S o u r c e L i n e s
311C=======================================================================
312 gbuf => elbuf_str%GBUF
313 nptr = elbuf_str%NPTR
314 npts = elbuf_str%NPTS
315 nptt = elbuf_str%NPTT
316 nlay = elbuf_str%NLAY
317c
318 jhbe = iparg(23)
319 iint = iparg(36)
320 igtyp = iparg(38)
321 idef = 0
322 nf1=nft+1
323 nf2=nf1-(numels8+numels10+numels20)
324 ibid = 0
325 rbid = zero
326C
327 DO i=lft,llt
328 rhocp(i) = pm(69,ixs(1,nft+i))
329 temp0(i) = pm(79,ixs(1,nft+i))
330 ENDDO
331C
332 CALL s16coor3(
333 1 x ,v ,ixs(1,nf1) ,ixs16(1,nf2),xx ,
334 2 yy ,zz ,vx ,vy ,vz ,
335 3 nc ,ngl ,mat ,pid ,mass ,
336 4 dtelem(nft+1),sti ,gbuf%SIG ,gbuf%EINT ,gbuf%RHO,
337 5 gbuf%QVIS ,temp0 ,temp ,nel ,glob_therm%NINTEMP)
338C
339 DO n=1,16
340 DO i=lft,llt
341 ul(i,n) = zero
342 ENDDO
343 ENDDO
344 DO i=lft,llt
345 volg(i) = zero
346 ENDDO
347C
348C----------------------------------------
349C INITIALISATION DE LA THERMIQUE
350C----------------------------------------
351 IF(jthe /=0) CALL atheri(mat,pm,gbuf%TEMP)
352C-----------------------------
353C POINTS D' INTEGRATION
354C-----------------------------
355 is = 1
356 DO it=1,nptt
357 DO ir=1,nptr
358 DO il=1,nlay
359c
360 lbuf => elbuf_str%BUFLY(il)%LBUF(ir,is,it)
361 mbuf => elbuf_str%BUFLY(il)%MAT(ir,is,it)
362 l_pla = elbuf_str%BUFLY(il)%L_PLA
363 l_sigb = elbuf_str%BUFLY(il)%L_SIGB
364 ip = ir + ( (il-1) + (it-1)*nlay )*nptr
365C
366 IF (iint == 1) THEN
367C--- GAUSS
368 wi = w_gauss(ir,nptr)*w_gauss(il,nlay)*w_gauss(it,nptt)
369c
370 CALL s16rst(
371 1 a_gauss(ir,nptr),a_gauss(il,nlay),a_gauss(it,nptt),ni ,
372 2 dnidr ,dnids ,dnidt )
373
374 CALL s16deri3(ngl,lbuf%OFF,
375 1 a_gauss(ir,nptr),a_gauss(il,nlay),a_gauss(it,nptt),wi,
376 2 dnidr ,dnids ,dnidt ,rx ,ry ,rz ,
377 3 sx ,sy ,sz ,tx ,ty ,tz ,
378 4 xx ,yy ,zz ,px ,py ,pz ,
379 5 lbuf%VOL,deltax ,stin ,ni ,volg ,ul ,lbuf%VOL0DP)
380 ELSEIF (iint == 2) THEN
381C--- LOBATO
382 wi = w_gauss(ir,nptr)*w_lobatto(il,nlay)*w_gauss(it,nptt)
383c
384 CALL s16rst(
385 1 a_gauss(ir,nptr),a_lobatto(il,nlay),a_gauss(it,nptt),ni ,
386 2 dnidr ,dnids ,dnidt )
387
388 CALL s16deri3(ngl,lbuf%OFF,
389 1 a_gauss(ir,nptr),a_lobatto(il,nlay),a_gauss(it,nptt),wi,
390 2 dnidr ,dnids ,dnidt ,rx ,ry ,rz ,
391 3 sx ,sy ,sz ,tx ,ty ,tz ,
392 4 xx ,yy ,zz ,px ,py ,pz ,
393 5 lbuf%VOL,deltax ,stin ,ni ,volg ,ul ,lbuf%VOL0DP )
394 ENDIF
395c
396 IF (jthe == 0 .and. glob_therm%NINTEMP > 0) THEN
397 CALL s20temp(nel,numnod,mvsiz,npe, nc,ni(1,ip),temp,tempel)
398 ELSE
399 tempel(1:nel) = temp0(1:nel)
400 ENDIF
401!
402 CALL matini(pm ,ixs ,nixs ,x ,
403 . geo ,ale_connectivity ,detonators ,iparg ,
404 . sigi ,nel ,skew ,igeo(1,1) ,
405 . ipart ,iparts ,
406 . mat ,ipm ,nsigs ,numsol ,ptsol ,
407 . ip ,ngl ,npf ,tf ,bufmat ,
408 . gbuf ,lbuf ,mbuf ,elbuf_str ,iloadp ,
409 . facload, deltax ,tempel )
410C----------------------------------------
411C INITIALISATION DES MASSES
412C----------------------------------------
413 CALL s20msi(lbuf%RHO ,mass ,lbuf%VOL ,dtelem(nft+1),sti ,
414 . lbuf%OFF ,lbuf%SIG ,lbuf%EINT ,dtx ,nel ,
415 . gbuf%OFF ,gbuf%SIG ,gbuf%EINT ,gbuf%RHO ,wi/eight)
416C
417C-------------------------------------------
418C INITIALISATION DES CONTRAINTES + ENERGIE+ DENSITE + PLA
419C---------------------------------------------
420 IF (mtn >= 28) THEN
421 nuvar = ipm(8,ixs(1,nft+1))
422 idef =1
423 ELSE
424 nuvar = 0
425 IF(mtn == 14 .OR. mtn == 12)THEN
426 idef =1
427 ELSEIF(mtn == 24)THEN
428 idef =1
429 ELSEIF(istrain == 1)THEN
430 IF(mtn == 1)THEN
431 idef =1
432 ELSEIF(mtn == 2)THEN
433 idef =1
434 ELSEIF(mtn == 4)THEN
435 idef =1
436 ELSEIF(mtn == 3.OR.mtn == 6.OR.mtn ==10.OR.
437 . mtn == 21.OR.mtn == 22.OR.
438 . mtn == 23.OR.mtn == 49)THEN
439 idef =1
440 ENDIF
441 ENDIF
442 ENDIF
443 CALL sigin20b(lbuf%SIG,pm ,lbuf%VOL,sigsp ,
444 . sigi ,lbuf%EINT,lbuf%RHO,mbuf%VAR,lbuf%STRA,
445 . ixs ,nixs ,nsigi ,ip ,nuvar ,
446 . nel ,iuser ,idef ,nsigs ,strsglob ,
447 . straglob,jhbe ,igtyp ,x ,lbuf%GAMA,
448 . mat ,lbuf%PLA ,l_pla ,ptsol ,lbuf%SIGB,
449 . l_sigb ,ipm ,bufmat ,lbuf%VOL0DP)
450 ENDDO
451 ENDDO
452 ENDDO
453C----------
454 DO i=lft,llt
455 aa = max(ul(i,1),ul(i,2),ul(i,3),ul(i,4),
456 . ul(i,5),ul(i,6),ul(i,7),ul(i,8))
457 bb = max(ul(i,9) ,ul(i,10),ul(i,11),ul(i,12),ul(i,13),ul(i,14),
458 . ul(i,15),ul(i,16))
459 deltax2(i) = aa/max(aa,bb)
460 aa = aa*thirty2
461 bb = bb*thirty2*third
462 deltax(i) = sqrt(two*volg(i)/max(aa,bb))
463 gbuf%VOL(i) = volg(i)
464 ENDDO
465C
466 aire(:) = zero
467 CALL dtmain(geo ,pm ,ipm ,pid ,mat ,fv ,
468 . gbuf%EINT ,gbuf%TEMP ,gbuf%DELTAX ,gbuf%RK ,gbuf%RE ,bufmat, deltax, aire,
469 . gbuf%VOL, dtx, igeo,igtyp)
470c
471 CALL s16mass3(
472 1 mass ,mas,partsav,iparts(nf1),mss(1,nf1),volg,
473 2 xx ,yy ,zz ,vx ,vy ,vz ,
474 3 nc ,sti,stifn ,deltax2 ,gbuf%RHO ,dtx ,
475 4 dtelem(nft+1),mssx(1,nf1),rhocp, mcp, mcps(1,nf1) ,
476 5 mcpsx(1,nf1) ,gbuf%FILL )
477C----------------------------------------
478c Failure model initialisation
479C----------------------------------------
480 CALL failini(elbuf_str,nptr,npts,nptt,nlay,
481 . ipm,sigsp,nsigi,fail_ini ,
482 . sigi,nsigs,ixs,nixs,ptsol,
483 . rnoise,perturb,mat_param)
484C------------------------------------------
485C assemblage des Volumes nodaux et Modules nodaux
486C (pour rigidites d'interface)
487C------------------------------------------
488 IF (i7stifs /= 0) THEN
489 ncc=16
490 CALL sbulk3(volg ,nc ,ncc,mat,pm ,
491 2 volnod,bvolnod,vns(1,nf1),bns(1,nf1),
492 3 vnsx(1,nf1),bnsx(1,nf1) ,gbuf%FILL)
493 ENDIF
494C------------------------------------------
495 DO i=lft,llt
496 IF(ixs(10,i+nft) /= 0) THEN
497 IF (igtyp/=0 .AND. igtyp /= 14 .AND. igtyp/=20 .AND.
498 . igtyp/=21) THEN
499 ipid1=ixs(nixs-1,i+nft)
500 CALL fretitl2(titr1,igeo(npropgi-ltitr+1,ipid1),ltitr)
501 CALL ancmsg(msgid=226,
502 . msgtype=msgerror,
503 . anmode=aninfo_blind_1,
504 . i1=igeo(1,ipid1),
505 . c1=titr1,
506 . i2=igtyp)
507 ENDIF
508 ENDIF
509 ENDDO
510C-----------
511 RETURN
512 END
atheri
subroutine atheri(mat, pm, temp)
Definition atheri.F:41
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
freform
subroutine freform(irunn, irfl, irfe, h3d_data, flag_cst_ams, dynain_data, sensors, dt, output, glob_therm)
Definition freform.F:88
failini
subroutine failini(elbuf_str, nptr, npts, nptt, nlay, ipm, sigsp, nsigi, fail_ini, sigi, nsigs, ix, nix, pt, rnoise, perturb, mat_param)
Definition failini.F:43
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
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
s16coor3
subroutine s16coor3(x, v, ixs, ixs16, xx, yy, zz, vx, vy, vz, nc, ngl, mxt, ngeo, mass, dtelem, sti, sigg, eintg, rhog, qg, temp0, temp, nel, nintemp)
Definition s16coor3.F:38
s16init3
subroutine s16init3(elbuf_str, mas, ixs, pm, x, detonators, geo, veul, ale_connectivity, iparg, dtelem, sigi, nel, skew, igeo, stifn, partsav, v, iparts, mss, ixs16, ipart, mssx, sigsp, nsigi, ipm, iuser, nsigs, volnod, bvolnod, vns, bns, vnsx, bnsx, ptsol, bufmat, mcp, mcps, mcpsx, temp, npf, tf, strsglob, straglob, fail_ini, iloadp, facload, rnoise, perturb, mat_param, glob_therm)
Definition s16init3.F:59
s16mass3
subroutine s16mass3(mass, ms, partsav, ipart, mss, volg, xx, yy, zz, vx, vy, vz, nc, sti, stifn, deltax2, rho, dtx, dtelem, mssx, rhocp, mcp, mcps, mcpsx, fill)
Definition s16mass3.F:34
s20msi
subroutine s20msi(rho, mass, volu, dtelem, sti, off, sig, eint, dtx, nel, offg, sigg, eintg, rhog, wip)
Definition s20mass3.F:281
sigin20b
subroutine sigin20b(sig, pm, vol, sigsp, sigi, eint, rho, uvar, eps, ix, nix, nsigi, ipt, nuvar, nel, iuser, idef, nsigs, strsglob, straglob, jhbe, igtyp, x, bufgama, mat, epsp, l_pla, pt, sigb, l_sigb, ipm, bufmat, voldp)
Definition s20mass3.F:350
sbulk3
subroutine sbulk3(volu, nc, nnc, mat, pm, volnod, bvolnod, vns, bns, vnsx, bnsx, fill)
Definition sbulk3.F:42
s16rst
subroutine s16rst(r, s, t, ni, dnidr, dnids, dnidt)
Definition s16deri3.F:495
s16deri3
subroutine s16deri3(ngl, off, r, s, t, w, dnidr, dnids, dnidt, dxdr, dydr, dzdr, dxds, dyds, dzds, dxdt, dydt, dzdt, xx, yy, zz, px, py, pz, vol, deltax, kxx, ni, volg, ul, voldp)
Definition s16deri3.F:37
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
starter
program starter
Definition starter.F:39
ale_connectivity_mod::t_ale_connectivity
Definition ale_connectivity_mod.F:254
detonators_mod::detonators_struct_
Definition detonators_mod.F:154