OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
forintc.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!|| forintc ../engine/source/elements/forintc.F
25!||--- called by ------------------------------------------------------
26!|| resol ../engine/source/engine/resol.F
27!||--- calls -----------------------------------------------------
28!|| c3forc3 ../engine/source/elements/sh3n/coque3n/c3forc3.F
29!|| c3forc3_crk ../engine/source/elements/xfem/c3forc3_crk.F
30!|| cbaforc3 ../engine/source/elements/shell/coqueba/cbaforc3.F
31!|| cdk6forc3 ../engine/source/elements/sh3n/coquedk6/cdk6forc3.F
32!|| cdkforc3 ../engine/source/elements/sh3n/coquedk/cdkforc3.F
33!|| cforc3 ../engine/source/elements/shell/coque/cforc3.F
34!|| cforc3_crk ../engine/source/elements/xfem/cforc3_crk.F
35!|| czforc3 ../engine/source/elements/shell/coquez/czforc3.F
36!|| czforc3_crk ../engine/source/elements/xfem/czforc3_crk.F
37!|| seatbelt_reduction_factor ../engine/source/tools/seatbelts/seatbelt_reduction_factor.F
38!|| section_3n ../engine/source/tools/sect/section_3n.F
39!|| section_c ../engine/source/tools/sect/section_c.F
40!|| sensor_energy_part ../engine/source/tools/sensor/sensor_energy_part.F
41!|| startimeg ../engine/source/system/timer.F
42!|| stoptimeg ../engine/source/system/timer.F
43!||--- uses -----------------------------------------------------
44!|| crackxfem_mod ../engine/share/modules/crackxfem_mod.F
45!|| drape_mod ../engine/share/modules/drape_mod.F
46!|| dt_mod ../engine/source/modules/dt_mod.F
47!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.f90
48!|| element_mod ../common_source/modules/elements/element_mod.F90
49!|| failwave_mod ../common_source/modules/failwave_mod.F
50!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
51!|| groupdef_mod ../common_source/modules/groupdef_mod.F
52!|| mat_elem_mod ../common_source/modules/mat_elem/mat_elem_mod.F90
53!|| nlocal_reg_mod ../common_source/modules/nlocal_reg_mod.F
54!|| output_mod ../common_source/modules/output/output_mod.f90
55!|| pinchtype_mod ../common_source/modules/pinchtype_mod.F
56!|| rupture_mod ../engine/share/modules/rupture_mod.F
57!|| sensor_mod ../common_source/modules/sensor_mod.F90
58!|| stack_mod ../engine/share/modules/stack_mod.F
59!|| table_mod ../engine/share/modules/table_mod.F
60!|| timer_mod ../engine/source/system/timer_mod.f90
61!||====================================================================
62 SUBROUTINE forintc(TIMERS,
63 1 PM ,GEO ,X ,A ,AR ,
64 2 V ,VR ,MS ,IN ,NLOC_DMG ,
65 3 WA ,STIFN ,STIFR ,FSKY ,CRKSKY ,
66 4 TF ,BUFMAT ,PARTSAV ,D ,MAT_ELEM ,
67 5 DR ,EANI ,TANI ,FANI ,
68 6 FSAV ,SENSORS ,SKEW ,FAILWAVE ,
69 7 DT2T ,THKE ,BUFGEO ,IADC ,IADTG ,
70 8 IPARG ,NPC ,IXC ,IXTG ,NELTST ,
71 A IPARI ,ITYPTST ,NSTRF ,
72 B IPART ,IPARTC ,IPARTTG ,SECFCUM ,
73 D FSAVD ,GROUP_PARAM_TAB,
74 F FZERO ,IXTG1 ,IADTG1 ,IGEO ,IPM ,
75 G MADFAIL ,XSEC ,ITASK ,MCP ,
76 H TEMP ,FTHE ,FTHESKY ,
77 I MS_PLY ,ZI_PLY ,INOD_PXFEM,XEDGE4N ,XEDGE3N ,
78 J IEL_PXFEM ,IADC_PXFEM,IGROUC ,NGROUC ,GRESAV ,
79 K GRTH ,IGRTH ,MSTG ,DMELTG ,MSC ,
80 L DMELC ,TABLE ,KNOD2ELC ,PTG ,MSZ2 ,
81 M INOD_CRK ,IEL_CRK ,IADC_CRK ,ELCUTC ,NODENR ,
82 N IBORDNODE ,NODEDGE ,CRKNODIAD ,ELBUF_TAB ,
83 O XFEM_TAB ,CONDN ,CONDNSKY ,CRKEDGE ,
84 P STACK ,ITAB ,GLOB_THERM,
85 S DRAPE_SH4N ,DRAPE_SH3N ,SUBSET ,XDP ,VPINCH,
86 T APINCH ,STIFPINCH ,DRAPEG,OUTPUT ,DT ,
87 Y SNPC , STF, USERL_AVAIL,MAXFUNC,
88 U SBUFMAT)
89C-----------------------------------------------
90C M o d u l e s
91C-----------------------------------------------
92 USE timer_mod
93 USE table_mod
94 USE mat_elem_mod
95 USE crackxfem_mod
96 USE stack_mod
97 USE failwave_mod
98 USE nlocal_reg_mod
99 USE groupdef_mod
100 USE pinchtype_mod
101 USE sensor_mod
102 USE rupture_mod
103 USE drape_mod
104 USE output_mod
105 USE elbufdef_mod
106 USE dt_mod
107 use glob_therm_mod
108 use element_mod , only : nixc,nixtg
109C-----------------------------------------------
110C I m p l i c i t T y p e s
111C-----------------------------------------------
112#include "implicit_f.inc"
113#include "comlock.inc"
114C-----------------------------------------------
115C G l o b a l P a r a m e t e r s
116C-----------------------------------------------
117#include "mvsiz_p.inc"
118#include "param_c.inc"
119C-----------------------------------------------
120C C o m m o n B l o c k s
121C-----------------------------------------------
122#include "com01_c.inc"
123#include "com04_c.inc"
124#include "com06_c.inc"
125#include "com08_c.inc"
126#include "com_xfem1.inc"
127#include "vect01_c.inc"
128#include "scr06_c.inc"
129#include "scr07_c.inc"
130#include "scr17_c.inc"
131#include "task_c.inc"
132#include "couple_c.inc"
133#include "impl1_c.inc"
134#include "stati_c.inc"
135C-----------------------------------------------
136C D u m m y A r g u m e n t s
137C-----------------------------------------------
138 TYPE(timer_), INTENT(INOUT) :: TIMERS
139 INTEGER,INTENT(IN) :: USERL_AVAIL ! Flag for User libraries availability
140 INTEGER,INTENT(IN) :: MAXFUNC ! Maximum number of functions
141 INTEGER,INTENT(IN) :: SBUFMAT ! Size of Bufmat
142 INTEGER,INTENT(IN) :: STF ! Size of TF
143 INTEGER,INTENT(IN) :: SNPC ! Size of NPC
144 INTEGER IXC(NIXC,*), IXTG(NIXTG,*), IGEO(NPROPGI,*), IPM(NPROPMI,*),
145 . NPC(*), IPARG(NPARG,*), IPARI(NPARI,*),
146 . NSTRF(*), IPART(LIPART1,*), IPARTC(*), IPARTTG(*),
147 . IADC(4,*), IADTG(3,*),NELTST,
148 . ITYPTST,IXTG1(4,*),XEDGE4N(4,*),XEDGE3N(3,*),
149 . IADTG1(3,*),MADFAIL(*),ITASK,
150 . INOD_PXFEM(*),IEL_PXFEM(*) ,IADC_PXFEM(4,*), IGROUC(*),
151 . NGROUC,GRTH(*),IGRTH(*),KNOD2ELC(*),
152 . INOD_CRK(*),IEL_CRK(*),IADC_CRK(*),ELCUTC(2,*),
153 . NODENR(*),IBORDNODE(*),NODEDGE(2,*),CRKNODIAD(*),
154 . ITAB(*)
155C REAL
156 my_real
157 . X(3,*) ,D(3,*) ,V(3,*) ,VR(3,*),
158 . MS(*) ,IN(*) ,PM(NPROPM,*),SKEW(LSKEW,*),
159 . GEO(NPROPG,*),BUFMAT(*) ,TF(STF) ,FSAV(NTHVKI,*) ,
160 . WA(*), THKE(*),
161 . A(3,*) ,AR(3,*) ,FANI(3,*) ,PARTSAV(NPSAV,*) ,
162 . STIFN(*) ,STIFR(*),FSKY(*) ,
163 . DR(3,*) ,TANI(*),EANI(*),
164 . bufgeo(*) ,dt2t, secfcum(7,numnod,nsect),
165 . fsavd(nthvki,*),
166 . fzero(3,4,(numelc+numeltg)),xsec(4,3,nsect),
167 . mcp(*),temp(*),fthe(*),fthesky(*),
168 . ms_ply(*), zi_ply(*),gresav(*),
169 . mstg(*), dmeltg(*), msc(*), dmelc(*),condn(*),condnsky(*),
170 . ptg(3,*),msz2(*),apinch(3,*),stifpinch(*),vpinch(3,*)
171 DOUBLE PRECISION :: XDP(3,*)
172 TYPE (TTABLE) TABLE(*)
173 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP) :: ELBUF_TAB
174 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP,NXEL) :: XFEM_TAB
175 TYPE (XFEM_EDGE_) , DIMENSION(*) :: CRKEDGE
176 TYPE (XFEM_SKY_) , DIMENSION(*) :: CRKSKY
177 TYPE (STACK_PLY) :: STACK
178 TYPE (FAILWAVE_STR_) :: FAILWAVE
179 TYPE (NLOCAL_STR_) :: NLOC_DMG
180 TYPE (GROUP_PARAM_) , DIMENSION(NGROUP) :: GROUP_PARAM_TAB
181 TYPE (SUBSET_) , DIMENSION(NSUBS) :: SUBSET
182 TYPE (DRAPE_) :: DRAPE_SH4N(NUMELC_DRAPE),DRAPE_SH3N(NUMELTG_DRAPE)
183 TYPE (DRAPEG_) :: DRAPEG
184 TYPE (SENSORS_) ,INTENT(INOUT) ,TARGET :: SENSORS
185 TYPE(OUTPUT_), INTENT(INOUT) :: OUTPUT !< output structure
186 TYPE (MAT_ELEM_) ,INTENT(INOUT) :: MAT_ELEM
187 TYPE (DT_) ,INTENT(IN) :: DT
188 type (glob_therm_) ,intent(inout) :: glob_therm
189C-----------------------------------------------
190C L o c a l V a r i a b l e s
191C-----------------------------------------------
192 DOUBLE PRECISION, POINTER :: pFBSAV6
193 INTEGER INDXOF(MVSIZ)
194 INTEGER I,II,J,N, NG, NVC, MLW, JFT, JLT,ISOLNOD,ITHK,IPLA,
195 . K2, NF1,IPRI, OFFSET,TYP,
196 . K0, K5, K9, NSG, NEL, KFTS,IOFC, ISTRA,
197 . icnod,nft1,isens_energy,
198 .
199 .
200 . iexpan, ig,itg1,itg2,itg3,nlevxf,iseatbelt
201 INTEGER IPARSENS,ISECT
202 INTEGER ISH3N,ISHPLYXFEM,IXFEM
203 INTEGER IXEL,ACTIFXFEM,ISUBSTACK
204 my_real
205 . fx(mvsiz,20),fy(mvsiz,20),fz(mvsiz,20),
206 . mx(mvsiz,4),my(mvsiz,4),mz(mvsiz,4)
207C
208 my_real
209 . uxint_mean(nxel*nxlaymax,mvsiz),uyint_mean(nxel*nxlaymax,mvsiz),
210 . uzint_mean(nxel*nxlaymax,mvsiz)
211C======================================================================|
212 itg1 = 1+4*ecrkxfec
213 itg2 = 2*numelc
214 itg3 = 1+numelc
215 ipri = 0
216 IF(mod(ncycle,iabs(ncpri))==0.OR.tt >= output%TH%THIS.OR.
217 + mdess /= 0.OR.tt >= output%TH%THIS1(1).OR.tt >= output%TH%THIS1(2)
218 + .OR.tt >= output%TH%THIS1(3) .OR.tt >= output%TH%THIS1(4).OR.tt >= output%TH%THIS1(5)
219 + .OR.tt >= output%TH%THIS1(6) .OR.tt >= output%TH%THIS1(7).OR.tt >= output%TH%THIS1(8)
220 + .OR.tt >= output%TH%THIS1(9) .OR.nth /= 0.OR.nanim/=0
221 + .OR.tt >= tabfis(1).OR.tt >= tabfis(2)
222 + .OR. tt >= tabfis(3).OR.tt >= tabfis(4).OR.tt >= tabfis(5)
223 + .OR. tt >= tabfis(6).OR.tt >= tabfis(7).OR.tt >= tabfis(8)
224 + .OR. tt >= tabfis(9).OR. tt >= tabfis(10).OR. istat==3) ipri=1
225!
226 isens_energy = 0
227 DO i=1,sensors%NSENSOR
228 typ = sensors%SENSOR_TAB(i)%TYPE
229 IF (typ == 14) isens_energy = 1 ! save internal/kinetic energy (PARTSAV)
230 ENDDO ! DO I=1,NSENSOR
231C
232C Dynamic parallel loop SMP
233C
234!$OMP DO SCHEDULE(DYNAMIC,1)
235 DO ig = 1, ngrouc
236C
237 ng = igrouc(ig)
238 sensors%NGR_SENSOR(itask+1) = ng
239 ngr_fail_wind = ng
240 IF (iparg(1, ng) == 151) THEN ! Bypass law 151
241 cycle
242 ENDIF
243C---------temporarily used to avoid pass KTBUF_STR everywhere
244 ng_imp = ng
245 IF(iparg(8,ng)==1)GOTO 250
246 ity =iparg(5,ng)
247C IF(ITY/=3.AND.ITY/=7)GOTO 250
248 IF (iddw>0) CALL startimeg(ng)
249 offset = 0
250 mlw = iparg(1,ng)
251C MLW= 0 ----> void
252C MLW = 13 ----> rigid material
253 IF (mlw == 0 .OR. mlw == 13) GOTO 250
254C---
255 nel = iparg(2,ng)
256 nft = iparg(3,ng)
257 npt = iparg(6,ng)
258 jale = iparg(7,ng)
259 ismstr = iparg(9,ng)
260 nsg = iparg(10,ng)
261 jeul = iparg(11,ng)
262 jtur = iparg(12,ng)
263 jthe = iparg(13,ng)
264 jlag = iparg(14,ng)
265 istra = iparg(44,ng)
266 nvc = iparg(19,ng)
267 jmult = iparg(20,ng)
268 jhbe = iparg(23,ng)
269 ish3n = iparg(23,ng)
270 jivf = iparg(24,ng)
271 jpor = iparg(27,ng)
272 ithk = iparg(28,ng)
273 isolnod = iparg(28,ng)
274 ipla = iparg(29,ng)
275 icnod = iparg(11,ng)
276 irep = iparg(35,ng)
277 iint = iparg(36,ng)
278 jcvt = iparg(37,ng)
279 igtyp = iparg(38,ng)
280 isorth = iparg(42,ng)
281 isorthg = isorth
282 israt = iparg(40,ng)
283 isrot = iparg(41,ng)
284 ifailure= iparg(43,ng)
285 kfts = iparg(30,ng)
286 jclose = iparg(33,ng)
287 icsen = iparg(39,ng)
288 iexpan = iparg(49,ng)
289 ishplyxfem = iparg(50,ng)
290 igre = iparg(51,ng)
291 jsms = iparg(52,ng)
292 ixfem = iparg(54,ng)
293 nlevxf = iparg(65,ng)
294 actifxfem=iparg(70,ng)
295 isubstack=iparg(71,ng)
296 iseatbelt=iparg(91,ng)
297 lft = 1
298 llt = min(nvsiz,nel)
299 mtn = mlw
300 jft=lft
301 jlt=llt
302 nf1 = nft+1
303 jsph=0
304C----6---------------------------------------------------------------7---------8
305 IF(ity==3)THEN
306 iofc = 0
307C
308C-- Reduction factor computation for seatbelts in sliprings
309 IF ((nslipring > 0).AND.(iseatbelt==1)) THEN
310 CALL seatbelt_reduction_factor(elbuf_tab(ng),jft,jlt,ixc,nel,x,nft)
311 ENDIF
312C
313 IF (isens_energy == 1 .AND. ipri == 0)
314 . CALL sensor_energy_part(ipartc(nf1) ,subset ,ipri,sensors,itask)
315 IF (jhbe >= 11.AND.jhbe <= 19) THEN
316
317 CALL cbaforc3(timers,
318 1 elbuf_tab(ng), jft, jlt, nft,
319 2 npt, ipari, mlw, ipri,
320 3 ithk, neltst, ityptst, itab,
321 4 mat_elem, istra, ipla, tt,
322 5 dt1, dt2t, pm, geo,
323 6 partsav, ixc(1,nf1), failwave, bufmat,
324 7 tf, npc, iadc(1,nf1), x,
325 8 d, dr, v, vr,
326 9 a, ar, stifn, stifr,
327 a fsky, tani, offset, eani,
328 b indxof, ipartc(nf1), thke(nf1), nvc,
329 c iofc, jhbe, fx(1,1), fx(1,2),
330 d fx(1,3), fx(1,4), fy(1,1), fy(1,2),
331 e fy(1,3), fy(1,4), fz(1,1), fz(1,2),
332 f fz(1,3), fz(1,4), mx(1,1), mx(1,2),
333 g mx(1,3), mx(1,4), my(1,1), my(1,2),
334 h my(1,3), my(1,4), mz(1,1), mz(1,2),
335 i mz(1,3), mz(1,4), kfts, ismstr,
336 j igeo, group_param_tab(ng),ipm, ifailure,
337 k itask, jthe, temp, fthe,
338 l fthesky, iexpan, ishplyxfem, ms,
339 m in, ms_ply, zi_ply, inod_pxfem,
340 n iel_pxfem, iadc_pxfem, gresav, grth,
341 o igrth(nf1), msc(nf1), dmelc(nf1), jsms,
342 p table, iparg(1,ng), sensors, msz2,
343 q condn, condnsky, isubstack, stack,
344 r drape_sh4n, nel, nloc_dmg, vpinch,
345 s apinch, stifpinch, drapeg%INDX_SH4N, igre,
346 t jtur, dt, ncycle, snpc,
347 y stf, glob_therm, nxlaymax, idel7nok,
348 u userl_avail, maxfunc, sbufmat ,
349 x ipart , lipart1)
350
351 ELSEIF (jhbe >= 21 .AND. jhbe <= 29) THEN
352
353 CALL czforc3(timers,
354 1 elbuf_tab(ng), jft, jlt, nft,
355 2 npt, itab, mlw,
356 3 ipri, ithk, neltst, istra,
357 4 ipla, dt1, dt2t,
358 5 pm, geo, partsav, ixc(1,nf1),
359 6 ityptst, bufmat, tf, npc,
360 7 iadc(1,nf1), failwave, x,
361 8 dr, v, vr, a,
362 9 ar, stifn, stifr, fsky,
363 a tani, indxof, ismstr,
364 b group_param_tab(ng),ipartc(nf1), thke(nf1), nvc,
365 c iofc, jhbe, fx(1,1), fx(1,2),
366 d fx(1,3), fx(1,4), fy(1,1), fy(1,2),
367 e fy(1,3), fy(1,4), fz(1,1), fz(1,2),
368 f fz(1,3), fz(1,4), mx(1,1), mx(1,2),
369 g mx(1,3), mx(1,4), my(1,1), my(1,2),
370 h my(1,3), my(1,4), mz(1,1), mz(1,2),
371 i mz(1,3), mz(1,4), kfts, fzero(1,1,nf1),
372 j igeo, ipm, ifailure, itask,
373 k jthe, temp, fthe, fthesky,
374 l iexpan, gresav, grth, igrth(nf1),
375 m xedge4n, msc(nf1), dmelc(nf1), jsms,
376 n table, iparg(1,ng), mat_elem , ixfem,
377 o knod2elc, sensors, elcutc(1,nf1), inod_crk,
378 p iel_crk, nodenr, iadc_crk,
379 q nodedge, crknodiad, condn, condnsky,
380 r stack, isubstack, xfem_tab(ng,1:nxel),crkedge,
381 s drape_sh4n, nel, nloc_dmg, drapeg%INDX_SH4N,
382 t igre, jtur, dt , ncycle,
383 y snpc, stf, glob_therm, idel7nok,
384 u userl_avail, maxfunc, sbufmat ,
385 x ipart , lipart1)
386
387 IF (icrack3d > 0 .AND. ixfem > 0 .AND. actifxfem > 0) THEN
388 DO ixel=1,nxel
389 CALL czforc3_crk(timers, xfem_tab(ng,ixel),
390 1 jft ,jlt ,nft ,ityptst ,
391 2 ipari ,mlw ,ipri ,ithk ,neltst ,
392 3 istra ,ipla ,tt ,dt1 ,dt2t ,
393 4 pm ,geo ,partsav ,ixc(1,nf1),group_param_tab(ng),
394 5 bufmat ,tf ,npc ,iadc(1,nf1),failwave ,
395 6 x ,d ,dr ,v ,vr ,
396 7 a ,ar ,stifn ,stifr ,fsky ,
397 8 tani ,offset ,eani ,indxof ,
398 9 ipartc(nf1),thke(nf1) ,nvc ,iofc ,jhbe ,
399 a fx(1,1) ,fx(1,2) ,fx(1,3) ,fx(1,4) ,fy(1,1) ,
400 b fy(1,2) ,fy(1,3) ,fy(1,4) ,fz(1,1) ,fz(1,2) ,
401 c fz(1,3) ,fz(1,4) ,mx(1,1) ,mx(1,2) ,mx(1,3) ,
402 d mx(1,4) ,my(1,1) ,my(1,2) ,my(1,3) ,my(1,4) ,
403 e mz(1,1) ,mz(1,2) ,mz(1,3) ,mz(1,4),
404 f kfts ,fzero(1,1,nf1),ismstr,mat_elem ,
405 i igeo ,ipm ,ifailure ,itask ,jthe ,
406 j temp ,fthe ,fthesky ,iexpan ,gresav ,
407 k grth ,igrth(nf1) ,msc(nf1) ,dmelc(nf1),jsms ,
408 l table ,iparg(1,ng),ixfem ,inod_crk ,iel_crk ,
409 m iadc_crk ,elcutc(1,nf1),crksky ,
410 n sensors ,ixel ,
411 o isubstack ,uxint_mean ,uyint_mean,uzint_mean,nlevxf ,
412 p nodedge ,crkedge ,stack ,drape_sh4n ,nloc_dmg,drapeg%INDX_SH4N,igre,
413 * dt ,ncycle ,snpc , stf ,glob_therm ,
414 * idel7nok ,userl_avail,maxfunc ,sbufmat ,
415 x ipart , lipart1)
416 ENDDO
417 ENDIF
418c
419 ELSE
420 CALL cforc3(timers,
421 1 elbuf_tab(ng), jft, jlt, pm,
422 2 ixc(1,nf1), x, a, ar,
423 3 v, vr, failwave, nvc,
424 4 mlw, geo, tf, npc,
425 5 bufmat, partsav, dt2t, neltst,
426 6 ityptst, stifn, stifr, fsky,
427 7 iadc(1,nf1), itab, d, dr,
428 8 tani, offset, eani, fx(1,1),
429 9 fx(1,2), fx(1,3), fx(1,4), fy(1,1),
430 a fy(1,2), fy(1,3), fy(1,4), fz(1,1),
431 b fz(1,2), fz(1,3), fz(1,4), mx(1,1),
432 c mx(1,2), mx(1,3), mx(1,4), my(1,1),
433 d my(1,2), my(1,3), my(1,4), mz(1,1),
434 e mz(1,2), mz(1,3), mz(1,4), indxof,
435 f ipartc(nf1), thke(nf1), group_param_tab(ng),mat_elem,
436 g nel, istra, jhbe, ithk,
437 h iofc, ipla, nft, ismstr,
438 i npt, kfts, fzero(1,1,nf1), igeo,
439 j ipm, ifailure, itask, jthe,
440 k temp, fthe, fthesky, iexpan,
441 l gresav, grth, xedge4n, igrth(nf1),
442 m msc(nf1), dmelc(nf1), jsms, table,
443 n iparg(1,ng), ixfem, knod2elc, sensors ,
444 o elcutc(1,nf1), inod_crk, iel_crk, ibordnode,
445 p nodenr, iadc_crk, nodedge, crknodiad,
446 q condn, condnsky, stack, isubstack,
447 r xfem_tab, ng ,crkedge, drape_sh4n, ipri,
448 s nloc_dmg, drapeg%INDX_SH4N, igre, jtur,
449 t output, dt, snpc, stf ,
450 u glob_therm, userl_avail, maxfunc, sbufmat ,
451 x ipart)
452c
453 IF (icrack3d > 0 .AND. ixfem > 0 .AND. actifxfem > 0) THEN
454 DO ixel=1,nxel
455 CALL cforc3_crk(timers, xfem_tab(ng,ixel),
456 1 jft ,jlt ,pm ,ixc(1,nf1),x ,
457 2 a ,ar ,v ,vr ,failwave ,
458 3 nvc ,mlw ,geo ,tf ,npc ,
459 4 bufmat ,partsav ,dt2t ,neltst ,ityptst ,
460 5 stifn ,stifr ,fsky ,crksky ,iadc(1,nf1),
461 6 d ,dr ,tani ,offset ,eani ,
462 7 fx(1,1) ,fx(1,2) ,fx(1,3) ,fx(1,4) ,fy(1,1) ,
463 8 fy(1,2) ,fy(1,3) ,fy(1,4) ,fz(1,1) ,fz(1,2) ,
464 9 fz(1,3) ,fz(1,4) ,mx(1,1) ,mx(1,2) ,mx(1,3) ,
465 a mx(1,4) ,my(1,1) ,my(1,2) ,my(1,3) ,my(1,4) ,
466 b mz(1,1) ,mz(1,2) ,mz(1,3) ,mz(1,4) ,indxof ,
467 c ipartc(nf1),thke(nf1),group_param_tab(ng),mat_elem ,
468 f nel ,istra ,jhbe ,kfts ,
469 g ithk ,iofc ,ipla ,nft ,ismstr ,
470 h fzero(1,1,nf1),igeo ,ipm ,ifailure ,itask ,
471 i jthe ,temp , fthe ,fthesky ,iexpan ,
472 j gresav ,grth ,
473 k igrth(nf1) ,msc(nf1) ,dmelc(nf1) ,jsms ,table ,
474 l iparg(1,ng) ,ixfem ,inod_crk ,iel_crk ,iadc_crk ,
475 m elcutc(1,nf1),
476 n sensors,ixel ,stack ,
477 o isubstack ,uxint_mean ,uyint_mean ,uzint_mean,nlevxf ,
478 p nodedge ,crkedge ,drape_sh4n ,ipri ,nloc_dmg ,
479 q drapeg%INDX_SH4N,igre ,dt ,snpc , stf ,
480 r glob_therm ,userl_avail, maxfunc,sbufmat ,
481 x ipart )
482 ENDDO
483 END IF
484
485 ENDIF
486!----------------------------------------------------------------------------
487 IF(nsect>0)THEN
488 k0=nstrf(25)
489 n=ninter+nrwall+nrbody
490 DO i=1,nsect
491 n=n+1
492 k2=k0+30+nstrf(k0+14)
493 k5=k0+30+nstrf(k0+14)+nstrf(k0+6)
494 . +2*nstrf(k0+7)+2*nstrf(k0+8)
495 iparsens=0
496 isect=0
497 IF (sensors%STABSEN > 0) isect=sensors%TABSENSOR(i+1)-sensors%TABSENSOR(i)
498 NULLIFY(pfbsav6)
499 IF(isect/=0) THEN
500 iparsens=1
501 pfbsav6 => sensors%FSAV(1,1,isect)
502 ENDIF
503C
504 CALL section_c(jft,jlt,nft,nstrf(k0+9),nstrf(k0+3),
505 2 nstrf(k0+4),nstrf(k0+5),nstrf(k5),x,v,vr,fsav(1,n),
506 3 ixc ,fani(1,1+2*(i-1)), secfcum(1,1,i) ,
507 4 fx ,fy ,fz ,mx ,my ,mz ,
508 5 nstrf(k0),nstrf(k0+14),nstrf(k0+26),nstrf(k0+6),
509 6 nstrf(k2),ms,
510 8 xsec(1,1,i) ,pfbsav6,iparsens)
511 k0 = nstrf(k0+24)
512 ENDDO
513 ENDIF
514 IF(nexmad/=0.AND.iofc/=0)THEN
515 imadfsh4=1
516 DO j=1,iofc
517 ii=indxof(j)+jft-1+nft
518 madfail(ii)=1
519 ENDDO
520 ENDIF
521C----6---------------------------------------------------------------7--
522 ELSEIF(ity==7)THEN
523 iofc = 0
524 IF (isens_energy == 1 .AND. ipri == 0)
525 . CALL sensor_energy_part(iparttg(nf1) ,subset ,ipri,sensors,itask)
526 IF (icnod == 6) THEN
527 nft1 = nf1-numeltg+numeltg6
528 CALL cdk6forc3(timers,
529 1 elbuf_tab(ng), jft, jlt, pm,
530 2 ixtg(1,nf1), x, a, ar,
531 3 v, vr, failwave, nvc,
532 4 mlw, geo, tf, npc,
533 5 bufmat, partsav, dt2t, neltst,
534 6 ityptst, stifn, stifr, fsky,
535 7 iadtg(1,nf1), group_param_tab(ng), tani(1+6*numelc),offset,
536 8 iparttg(nf1), thke(numelc+nf1), fx(1,1), fx(1,2),
537 9 fx(1,3), fy(1,1), fy(1,2), fy(1,3),
538 a fz(1,1), fz(1,2), fz(1,3), fx(1,4),
539 b fx(1,5), fx(1,6), fy(1,4), fy(1,5),
540 c fy(1,6), fz(1,4), fz(1,5), fz(1,6),
541 d mat_elem, nel, istra, ish3n,
542 e ithk, iofc, ipla, nft,
543 f ismstr, npt, kfts, ixtg1(1,nft1),
544 g iadtg1(1,nft1), igeo, ipm, ifailure,
545 h iexpan, gresav, grth, igrth(numelc+numelt+numelp+numelr+nf1),
546 i mstg(nf1), dmeltg(nf1), jsms, table,
547 j iparg(1,ng), sensors, ptg(1,nf1), jthe,
548 k condn, condnsky, isubstack, stack,
549 l itask, drape_sh3n, ipri, nloc_dmg,
550 m drapeg%INDX_SH3N,igre, jtur, dt ,
551 n ncycle, snpc, stf , glob_therm ,
552 o nxlaymax, idel7nok, userl_avail, maxfunc,
553 p sbufmat ,
554 x ipart , lipart1)
555 ELSE
556
557 IF (ish3n == 30) THEN
558 CALL cdkforc3(timers,
559 1 elbuf_tab(ng), jft, jlt, pm,
560 2 ixtg(1,nf1), x, a, ar,
561 3 v, vr, failwave, nvc,
562 4 mlw, geo, tf, npc,
563 5 bufmat, partsav, dt2t, neltst,
564 6 ityptst, stifn, stifr, fsky,
565 7 iadtg(1,nf1), itab, tani(1+6*numelc),iparttg(nf1),
566 8 thke(numelc+nf1), group_param_tab(ng), fx(1,1), fx(1,2),
567 9 fx(1,3), fy(1,1), fy(1,2), fy(1,3),
568 a fz(1,1), fz(1,2), fz(1,3), mx(1,1),
569 b mx(1,2), mx(1,3), my(1,1), my(1,2),
570 c my(1,3), mz(1,1), mz(1,2), mz(1,3),
571 d mat_elem, nel, istra, ish3n,
572 e ithk, iofc, ipla, nft,
573 f ismstr, npt, kfts, igeo,
574 g ipm, ifailure, gresav, grth,
575 h igrth(numelc+numelt+numelp+numelr+nf1),mstg(nf1), dmeltg(nf1), jsms,
576 i table, iparg(1,ng), sensors, ptg(1,nf1),
577 j jthe, condn, condnsky, isubstack,
578 k stack, itask, drape_sh3n, ipri,
579 l nloc_dmg, drapeg%INDX_SH3N, igre, jtur,
580 m dt, ncycle, snpc, stf,
581 m glob_therm , nxlaymax, idel7nok, userl_avail,
582 m maxfunc, sbufmat ,
583 x ipart , lipart1)
584
585 ELSE
586 CALL c3forc3(timers,
587 1 elbuf_tab(ng), jft, jlt, pm,
588 2 ixtg(1,nf1), x, a, ar,
589 3 v, vr, failwave, nvc,
590 4 mlw, geo, tf, npc,
591 5 bufmat, partsav, dt2t, neltst,
592 6 ityptst, stifn, stifr, fsky,
593 7 iadtg(1,nf1), itab, tani(1+6*numelc), offset,
594 8 iparttg(nf1), thke(numelc+nf1), fx(1,1), fx(1,2),
595 9 fx(1,3), fy(1,1), fy(1,2), fy(1,3),
596 a fz(1,1), fz(1,2), fz(1,3), mx(1,1),
597 b mx(1,2), mx(1,3), my(1,1), my(1,2),
598 c my(1,3), mz(1,1), mz(1,2), mz(1,3),
599 d group_param_tab(ng), mat_elem, nel, istra,
600 e ish3n, xedge3n, ithk, iofc,
601 f ipla, nft, ismstr, npt,
602 g kfts, fzero(1,1,nf1+numelc), igeo, ipm,
603 h ifailure, itask, jthe, temp,
604 i fthe, fthesky, iexpan, gresav,
605 j grth, igrth(numelc+numelt+numelp+numelr+nf1),mstg(nf1), dmeltg(nf1),
606 k jsms, table, iparg(1,ng), ixfem,
607 l sensors, ptg(1,nf1), ibordnode, elcutc(1,nf1+itg3-1),
608 m inod_crk, iel_crk(itg3), nodenr, iadc_crk(itg1),
609 n nodedge, crknodiad, knod2elc, condn,
610 o condnsky, stack, isubstack, xfem_tab(ng,1:nxel),
611 p crkedge, drape_sh3n, ipri, nloc_dmg,
612 q xdp, drapeg%INDX_SH3N, igre, jtur,
613 r dt, snpc, stf, glob_therm,
614 s idel7nok, userl_avail, maxfunc, sbufmat ,
615 x ipart , lipart1 )
616
617 IF (icrack3d > 0 .AND. ixfem > 0 .AND. actifxfem > 0) THEN
618 DO ixel=1,nxel
619 CALL c3forc3_crk(timers,
620 1 xfem_tab(ng,ixel), jft, jlt, pm,
621 2 ixtg(1,nf1), x, a, ar,
622 3 v, vr, failwave, nvc,
623 4 mlw, geo, tf, npc,
624 5 bufmat, partsav, dt2t, neltst,
625 6 ityptst, stifn, stifr, fsky,
626 7 crksky, iadtg(1,nf1), tani(1+6*numelc), offset,
627 8 iparttg(nf1), thke(numelc+nf1), fx(1,1), fx(1,2),
628 9 fx(1,3), fy(1,1), fy(1,2), fy(1,3),
629 a fz(1,1), fz(1,2), fz(1,3), mx(1,1),
630 b mx(1,2), mx(1,3), my(1,1), my(1,2),
631 c my(1,3), mz(1,1), mz(1,2), mz(1,3),
632 d kfts, group_param_tab(ng), mat_elem, nel,
633 e istra, ish3n, ithk, iofc,
634 f ipla, nft, ismstr, fzero(1,1,nf1+numelc),
635 g igeo, ipm, ifailure, itask,
636 h jthe, temp, fthe, fthesky,
637 i iexpan, gresav, grth, igrth(numelc+numelt+numelp+numelr+nf1),
638 j mstg(nf1), dmeltg(nf1), jsms, table,
639 k iparg(1,ng), sensors, ptg(1,nf1), ixfem,
640 l inod_crk, iel_crk(itg3), iadc_crk(itg1), elcutc(1,nf1+itg3-1),
641 m ixel, stack, isubstack, uxint_mean,
642 n uyint_mean, uzint_mean, nlevxf, nodedge,
643 o crkedge, drape_sh3n, ipri, nloc_dmg,
644 p drapeg%INDX_SH3N, igre, dt , ncycle ,
645 q snpc, stf, glob_therm , idel7nok,
646 a userl_avail, maxfunc, sbufmat ,
647 x ipart , lipart1 )
648 ENDDO
649 ENDIF
650 ENDIF
651 ENDIF
652!----------------------------------------------------------------------------------------
653 IF(nsect>0)THEN
654 k0=nstrf(25)
655 n=ninter+nrwall+nrbody
656 DO i=1,nsect
657 n=n+1
658 k2=k0+30+nstrf(k0+14)
659 k9=k0+30+nstrf(k0+14)+nstrf(k0+6)
660 1 +2*nstrf(k0+7)+2*nstrf(k0+8)+2*nstrf(k0+9)
661 2 +2*nstrf(k0+10)+2*nstrf(k0+11)+2*nstrf(k0+12)
662 iparsens=0
663 isect=0
664 IF (sensors%STABSEN > 0) isect=sensors%TABSENSOR(i+1)-sensors%TABSENSOR(i)
665 NULLIFY(pfbsav6)
666 IF(isect/=0) THEN
667 iparsens=1
668 pfbsav6 => sensors%FSAV(1,1,isect)
669 ENDIF
670C
671 CALL section_3n(jft,jlt,nft,nstrf(k0+13),nstrf(k0+3),
672 2 nstrf(k0+4),nstrf(k0+5),nstrf(k9),x,v,vr,fsav(1,n),
673 3 ixtg ,fani(1,1+2*(i-1)), secfcum(1,1,i) ,
674 4 fx ,fy ,fz ,mx ,my, mz,
675 5 nstrf(k0),nstrf(k0+14),nstrf(k0+26),nstrf(k0+6),
676 6 nstrf(k2),ms,
677 8 xsec(1,1,i) , pfbsav6,iparsens)
678 k0=nstrf(k0+24)
679 ENDDO
680 ENDIF
681 ENDIF
682 IF (iddw>0) CALL stoptimeg(ng)
683
684 250 CONTINUE
685 END DO
686!$OMP END DO
687C----6---------------------------------------------------------------7---------8
688 RETURN
689 END
690
691
c3forc3
subroutine c3forc3(timers, elbuf_str, jft, jlt, pm, ixtg, x, f, m, v, r, failwave, nvc, mtn, geo, tf, npf, bufmat, pmsav, dt2t, neltst, ityptst, stifn, stifr, fsky, iadtg, itab, epsdot, offset, iparttg, thke, f11, f12, f13, f21, f22, f23, f31, f32, f33, m11, m12, m13, m21, m22, m23, m31, m32, m33, group_param, mat_elem, nel, istrain, ish3n, xedge3n, ithk, iofc, ipla, nft, ismstr, npt, kfts, fzero, igeo, ipm, ifailure, itask, jthe, temp, fthe, fthesky, iexpan, gresav, grth, igrth, mstg, dmeltg, jsms, table, iparg, ixfem, sensors, ptg, ibordnode, elcutc, inod_crk, iel_crk, nodenr, iadtg_crk, nodedge, crknodiad, knod2elc, condn, condnsky, stack, isubstack, xfem_str, crkedge, drape_sh3n, ipri, nloc_dmg, xdp, indx_drape, igre, jtur, dt, snpc, stf, glob_therm, idel7nok, userl_avail, maxfunc, sbufmat, ipart, lipart1)
Definition c3forc3.F:114
c3forc3_crk
subroutine c3forc3_crk(timers, xfem_str, jft, jlt, pm, ixtg, x, f, m, v, r, failwave, nvc, mtn, geo, tf, npf, bufmat, pmsav, dt2t, neltst, ityptst, stifn, stifr, fsky, crksky, iadtg, epsdot, offset, iparttg, thke, f11, f12, f13, f21, f22, f23, f31, f32, f33, m11, m12, m13, m21, m22, m23, m31, m32, m33, kfts, group_param, mat_elem, nel, istrain, ish3n, ithk, iofc, ipla, nft, ismstr, fzero, igeo, ipm, ifailure, itask, jthe, temp, fthe, fthesky, iexpan, gresav, grth, igrth, mstg, dmeltg, jsms, table, iparg, sensors, ptg, ixfem, inod_crk, iel_crk, iadtg_crk, elcutc, ixel, stack, isubstack, uxint_mean, uyint_mean, uzint_mean, nlevxf, nodedge, crkedge, drape_sh3n, ipri, nloc_dmg, indx_drape, igre, dt, ncycle, snpc, stf, glob_therm, idel7nok, userl_avail, maxfunc, sbufmat, ipart, lipart1)
Definition c3forc3_crk.F:96
cbaforc3
subroutine cbaforc3(timers, elbuf_str, jft, jlt, nft, npt, ipari, mtn, ipri, ithk, neltst, ityptst, itab, mat_elem, istrain, ipla, tt, dt1, dt2t, pm, geo, partsav, ixc, failwave, bufmat, tf, npf, iadc, x, d, dr, v, vr, f, m, stifn, stifr, fsky, tani, offset, eani, indxof, ipartc, thke, nvc, iofc, ihbe, f11, f12, f13, f14, f21, f22, f23, f24, f31, f32, f33, f34, m11, m12, m13, m14, m21, m22, m23, m24, m31, m32, m33, m34, kfts, ismstr, igeo, group_param, ipm, ifailure, itask, jthe, temp, fthe, fthesky, iexpan, ishplyxfem, ms, in, ms_ply, zi_ply, inod_pxfem, iel_pxfem, iadc_pxfem, gresav, grth, igrth, msc, dmelc, jsms, table, iparg, sensors, msz2, condn, condnsky, isubstack, stack, drape_sh4n, nel, nloc_dmg, vpinch, fpinch, stifpinch, indx_drape, igre, jtur, dt, ncycle, snpc, stf, glob_therm, nxlaymax, idel7nok, userl_avail, maxfunc, sbufmat, ipart, lipart1)
Definition cbaforc3.F:131
cdk6forc3
subroutine cdk6forc3(timers, elbuf_str, jft, jlt, pm, ixtg, x, f, m, v, r, failwave, nvc, mtn, geo, tf, npf, bufmat, pmsav, dt2t, neltst, ityptst, stifn, stifr, fsky, iadtg, group_param, epsdot, offset, iparttg, thke, f11, f12, f13, f21, f22, f23, f31, f32, f33, f14, f15, f16, f24, f25, f26, f34, f35, f36, mat_elem, nel, istrain, ihbe, ithk, iofc, ipla, nft, ismstr, npt, kfts, ixtg1, iadtg1, igeo, ipm, ifailure, iexpan, gresav, grth, igrth, mstg, dmeltg, jsms, table, iparg, sensors, ptg, jthe, condn, condnsky, isubstack, stack, itask, drape_sh3n, ipri, nloc_dmg, indx_drape, igre, jtur, dt, ncycle, snpc, stf, glob_therm, nxlaymax, idel7nok, userl_avail, maxfunc, sbufmat, ipart, lipart1)
Definition cdk6forc3.F:85
cdkforc3
subroutine cdkforc3(timers, elbuf_str, jft, jlt, pm, ixtg, x, f, m, v, r, failwave, nvc, mtn, geo, tf, npf, bufmat, pmsav, dt2t, neltst, ityptst, stifn, stifr, fsky, iadtg, itab, epsdot, iparttg, thke, group_param, f11, f12, f13, f21, f22, f23, f31, f32, f33, m11, m12, m13, m21, m22, m23, m31, m32, m33, mat_elem, nel, istrain, ihbe, ithk, iofc, ipla, nft, ismstr, npt, kfts, igeo, ipm, ifailure, gresav, grth, igrth, mstg, dmeltg, jsms, table, iparg, sensors, ptg, jthe, condn, condnsky, isubstack, stack, itask, drape_sh3n, ipri, nloc_dmg, indx_drape, igre, jtur, dt, ncycle, snpc, stf, glob_therm, nxlaymax, idel7nok, userl_avail, maxfunc, sbufmat, ipart, lipart1)
Definition cdkforc3.F:88
cforc3
subroutine cforc3(timers, elbuf_str, jft, jlt, pm, ixc, x, f, m, v, vr, failwave, nvc, mtn, geo, tf, npf, bufmat, partsav, dt2t, neltst, ityptst, stifn, stifr, fsky, iadc, itab, d, dr, tani, offset, eani, f11, f12, f13, f14, f21, f22, f23, f24, f31, f32, f33, f34, m11, m12, m13, m14, m21, m22, m23, m24, m31, m32, m33, m34, indxof, ipartc, thke, group_param, mat_elem, nel, istrain, ihbe, ithk, iofc, ipla, nft, ismstr, npt, kfts, fzero, igeo, ipm, ifailure, itask, jthe, temp, fthe, fthesky, iexpan, gresav, grth, xedge4n, igrth, msc, dmelc, jsms, table, iparg, ixfem, knod2elc, sensors, elcutc, inod_crk, iel_crk, ibordnode, nodenr, iadc_crk, nodedge, crknodiad, condn, condnsky, stack, isubstack, xfem_str, ig, crkedge, drape_sh4n, ipri, nloc_dmg, indx_drape, igre, jtur, output, dt, snpc, stf, glob_therm, userl_avail, maxfunc, sbufmat, ipart)
Definition cforc3.F:113
cforc3_crk
subroutine cforc3_crk(timers, xfem_str, jft, jlt, pm, ixc, x, f, m, v, vr, failwave, nvc, mtn, geo, tf, npf, bufmat, partsav, dt2t, neltst, ityptst, stifn, stifr, fsky, crksky, iadc, d, dr, tani, offset, eani, f11, f12, f13, f14, f21, f22, f23, f24, f31, f32, f33, f34, m11, m12, m13, m14, m21, m22, m23, m24, m31, m32, m33, m34, indxof, ipartc, thke, group_param, mat_elem, nel, istrain, ihbe, kfts, ithk, iofc, ipla, nft, ismstr, fzero, igeo, ipm, ifailure, itask, jthe, temp, fthe, fthesky, iexpan, gresav, grth, igrth, msc, dmelc, jsms, table, iparg, ixfem, inod_crk, iel_crk, iadc_crk, elcutc, sensors, ixel, stack, isubstack, uxint_mean, uyint_mean, uzint_mean, nlevxf, nodedge, crkedge, drape_sh4n, ipri, nloc_dmg, indx_drape, igre, dt, snpc, stf, glob_therm, userl_avail, maxfunc, sbufmat, ipart)
Definition cforc3_crk.F:91
czforc3
subroutine czforc3(timers, elbuf_str, jft, jlt, nft, npt, itab, mtn, ipri, ithk, neltst, istrain, ipla, dt1, dt2t, pm, geo, partsav, ixc, ityptst, bufmat, tf, npf, iadc, failwave, x, dr, v, vr, f, m, stifn, stifr, fsky, tani, indxof, ismstr, group_param, ipartc, thke, nvc, iofc, ihbe, f11, f12, f13, f14, f21, f22, f23, f24, f31, f32, f33, f34, m11, m12, m13, m14, m21, m22, m23, m24, m31, m32, m33, m34, kfts, fzero, igeo, ipm, ifailure, itask, jthe, temp, fthe, fthesky, iexpan, gresav, grth, igrth, xedge4n, msc, dmelc, jsms, table, iparg, mat_elem, ixfem, knod2elc, sensors, elcutc, inod_crk, iel_crk, nodenr, iadc_crk, nodedge, crknodiad, condn, condnsky, stack, isubstack, xfem_str, crkedge, drape_sh4n, nel, nloc_dmg, indx_drape, igre, jtur, dt, ncycle, snpc, stf, glob_therm, idel7nok, userl_avail, maxfunc, sbufmat, ipart, lipart1)
Definition czforc3.F:116
czforc3_crk
subroutine czforc3_crk(timers, xfem_str, jft, jlt, nft, ityptst, ipari, mtn, ipri, ithk, neltst, istrain, ipla, tt, dt1, dt2t, pm, geo, partsav, ixc, group_param, bufmat, tf, npf, iadc, failwave, x, d, dr, v, vr, f, m, stifn, stifr, fsky, tani, offset, eani, indxof, ipartc, thke, nvc, iofc, ihbe, f11, f12, f13, f14, f21, f22, f23, f24, f31, f32, f33, f34, m11, m12, m13, m14, m21, m22, m23, m24, m31, m32, m33, m34, kfts, fzero, ismstr, mat_elem, igeo, ipm, ifailure, itask, jthe, temp, fthe, fthesky, iexpan, gresav, grth, igrth, msc, dmelc, jsms, table, iparg, ixfem, inod_crk, iel_crk, iadc_crk, elcutc, crksky, sensors, ixel, isubstack, uxint_mean, uyint_mean, uzint_mean, nlevxf, nodedge, crkedge, stack, drape_sh4n, nloc_dmg, indx_drape, igre, dt, ncycle, snpc, stf, glob_therm, idel7nok, userl_avail, maxfunc, sbufmat, ipart, lipart1)
Definition czforc3_crk.F:91
startimeg
subroutine startimeg(ng)
Definition timer.F:1371
stoptimeg
subroutine stoptimeg(ng)
Definition timer.F:1419
forintc
subroutine forintc(timers, pm, geo, x, a, ar, v, vr, ms, in, nloc_dmg, wa, stifn, stifr, fsky, crksky, tf, bufmat, partsav, d, mat_elem, dr, eani, tani, fani, fsav, sensors, skew, failwave, dt2t, thke, bufgeo, iadc, iadtg, iparg, npc, ixc, ixtg, neltst, ipari, ityptst, nstrf, ipart, ipartc, iparttg, secfcum, fsavd, group_param_tab, fzero, ixtg1, iadtg1, igeo, ipm, madfail, xsec, itask, mcp, temp, fthe, fthesky, ms_ply, zi_ply, inod_pxfem, xedge4n, xedge3n, iel_pxfem, iadc_pxfem, igrouc, ngrouc, gresav, grth, igrth, mstg, dmeltg, msc, dmelc, table, knod2elc, ptg, msz2, inod_crk, iel_crk, iadc_crk, elcutc, nodenr, ibordnode, nodedge, crknodiad, elbuf_tab, xfem_tab, condn, condnsky, crkedge, stack, itab, glob_therm, drape_sh4n, drape_sh3n, subset, xdp, vpinch, apinch, stifpinch, drapeg, output, dt, snpc, stf, userl_avail, maxfunc, sbufmat)
Definition forintc.F:89
min
#define min(a, b)
Definition macros.h:20
crackxfem_mod
Definition crackxfem_mod.F:121
drape_mod
Definition drape_mod.F:54
dt_mod
Definition dt_mod.F:80
failwave_mod
Definition failwave_mod.F:61
groupdef_mod
Definition groupdef_mod.F:668
nlocal_reg_mod
Definition nlocal_reg_mod.F:119
pinchtype_mod
Definition pinchtype_mod.F:69
rupture_mod
Definition rupture_mod.F:30
rupture_mod::ngr_fail_wind
integer ngr_fail_wind
Definition rupture_mod.F:41
stack_mod
Definition stack_mod.F:75
table_mod
Definition table_mod.F:113
seatbelt_reduction_factor
subroutine seatbelt_reduction_factor(elbuf_str, jft, jlt, ixc, nel, x, nft)
Definition seatbelt_reduction_factor.F:32
section_3n
subroutine section_3n(lft, llt, nft, nseg, n1, n2, n3, nstrf, x, v, vr, fsav, ixtg, fopta, secfcum, fx, fy, fz, mx, my, mz, type, nsint, ifram, nnod, nod, ms, xsec, fbsav6, iparsens)
Definition section_3n.F:40
section_c
subroutine section_c(lft, llt, nft, nseg, n1, n2, n3, nstrf, x, v, vr, fsav, ixc, fopta, secfcum, fx, fy, fz, mx, my, mz, type, nsint, ifram, nnod, nod, ms, xsec, fbsav6, iparsens)
Definition section_c.F:40
sensor_energy_part
subroutine sensor_energy_part(ipart, subset, ipri, sensors, itask)
Definition sensor_energy_part.F:33