OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s4forc3.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!|| s4forc3 ../engine/source/elements/solid/solide4/s4forc3.F
25!||--- called by ------------------------------------------------------
26!|| alemain ../engine/source/ale/alemain.F
27!|| forint ../engine/source/elements/forint.F
28!||--- calls -----------------------------------------------------
29!|| a4mass3 ../engine/source/elements/solid/solide4/a4mass3.F
30!|| a4mass3p ../engine/source/elements/solid/solide4/a4mass3p.F
31!|| a4momt3 ../engine/source/elements/solid/solide4/a4momt3.F
32!|| boltst ../engine/source/elements/solid/solide/boltst.F
33!|| check_off_ale ../engine/source/elements/solid/solide/check_off_ale.f
34!|| e4pxle3 ../engine/source/elements/solid/solide4/e4pxle3.F
35!|| mmain ../engine/source/materials/mat_share/mmain.F90
36!|| s11fx3 ../engine/source/elements/solid/solide/s11fx3.F
37!|| s4bilan ../engine/source/elements/solid/solide4/s4bilan.F
38!|| s4coor3 ../engine/source/elements/solid/solide4/s4coor3.F
39!|| s4cumu3 ../engine/source/elements/solid/solide4/s4cumu3.F
40!|| s4cumu3p ../engine/source/elements/solid/solide4/s4cumu3p.F
41!|| s4defo3 ../engine/source/elements/solid/solide4/s4defo3.F
42!|| s4defot3 ../engine/source/elements/solid/solide4/s4defot3.F
43!|| s4deri3 ../engine/source/elements/solid/solide4/s4deri3.F
44!|| s4derit3 ../engine/source/elements/solid/solide4/s4derit3.F
45!|| s4derito3 ../engine/source/elements/solid/solide4/s4derito3.F
46!|| s4fillopt ../engine/source/elements/solid/solide4/s4fillopt.F
47!|| s4fint3 ../engine/source/elements/solid/solide4/s4fint3.F
48!|| s4fint_reg ../engine/source/elements/solid/solide4/s4fint_reg.F
49!|| s4for_distor ../engine/source/elements/solid/solide4/s4for_distor.F
50!|| s4malla3 ../engine/source/elements/solid/solide4/s4mall3.F
51!|| s4rcoor12 ../engine/source/elements/solid/solide4/s4rcoor12.F
52!|| s4sav12 ../engine/source/elements/solid/solide4/s4sav12.F
53!|| s4sav3 ../engine/source/elements/solid/solide4/s4sav3.f
54!|| s4therm ../engine/source/elements/solid/solide4/s4therm.f
55!|| s4upd11t12 ../engine/source/elements/solid/solide4/s4upd11t12.F
56!|| s4voln_m ../engine/source/elements/solid/solide4_sfem/s4voln_m.F
57!|| sdistor_ini ../engine/source/elements/solid/solide/sdistror_ini.F90
58!|| sgcoor3 ../engine/source/elements/solid/solide/sgcoor3.f
59!|| sgeodel3 ../engine/source/elements/solid/solide/sgeodel3.F
60!|| smallb3 ../engine/source/elements/solid/solide/smallb3.F
61!|| smallgeo3 ../engine/source/elements/solid/solide/smallgeo3.F
62!|| sordef12 ../engine/source/elements/solid/solidez/sordef12.F
63!|| sordeft12 ../engine/source/elements/solid/solidez/sordeft12.F
64!|| sorthdir3 ../engine/source/elements/solid/solide/sorthdir3.F
65!|| sreploc3 ../engine/source/elements/solid/solide/sreploc3.F
66!|| srho3 ../engine/source/elements/solid/solide/srho3.f
67!|| srota3 ../engine/source/elements/solid/solide/srota3.F
68!|| sroto12_sig ../engine/source/elements/solid/solidez/sroto12_sig.F
69!|| sstra3 ../engine/source/elements/solid/solide/sstra3.F
70!|| startime ../engine/source/system/timer_mod.F90
71!|| stoptime ../engine/source/system/timer_mod.F90
72!||--- uses -----------------------------------------------------
73!|| ale_connectivity_mod ../common_source/modules/ale/ale_connectivity_mod.F
74!|| ale_mod ../common_source/modules/ale/ale_mod.F
75!|| dt_mod ../engine/source/modules/dt_mod.F
76!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
77!|| element_mod ../common_source/modules/elements/element_mod.F90
78!|| glob_therm_mod ../common_source/modules/mat_elem/glob_therm_mod.F90
79!|| mat_elem_mod ../common_source/modules/mat_elem/mat_elem_mod.F90
80!|| mmain_mod ../engine/source/materials/mat_share/mmain.F90
81!|| nlocal_reg_mod ../common_source/modules/nlocal_reg_mod.f
82!|| output_mod ../common_source/modules/output/output_mod.f90
83!|| sdistor_ini_mod ../engine/source/elements/solid/solide/sdistror_ini.F90
84!|| sensor_mod ../common_source/modules/sensor_mod.F90
85!|| table_mod ../engine/share/modules/table_mod.F
86!|| timer_mod ../engine/source/system/timer_mod.F90
87!||====================================================================
88 SUBROUTINE s4forc3(TIMERS, OUTPUT,
89 1 ELBUF_TAB, NG, PM, GEO,
90 2 IXS, X, A, V,
91 3 MS, W, FLUX, FLU1,
92 4 VEUL, FV, ALE_CONNECT,IPARG,
93 5 TF, NPF, BUFMAT, PARTSAV,
94 6 NLOC_DMG, DT2T, NELTST, ITYPTST,
95 7 STIFN, FSKY, IADS, OFFSET,
96 8 EANI, IPARTS, F11, F21,
97 9 F31, F12, F22, F32,
98 A F13, F23, F33, F14,
99 B F24, F34, NEL, FSKYM,
100 C MSNF, IPM, IGEO, BUFVOIS,
101 D ISTRAIN, ITASK, TEMP, FTHE,
102 E FTHESKY, IEXPAN, GRESAV, GRTH,
103 F IGRTH, MSSA, DMELS, TABLE,
104 G XDP, SFEM_NODVAR, VOLN, CONDN,
105 H CONDNSKY, D, SENSORS, IOUTPRT,
106 I MAT_ELEM, H3D_STRAIN, DT, IDEL7NOK,NSVOIS,SZ_BUFVOIS,
107 * SNPC, STF, SBUFMAT, SVIS,IDTMINS,IRESP,
108 * IDEL7NG, MAXFUNC, USERL_AVAIL, GLOB_THERM,
109 V impl_s, idyna, S_SFEM_NODVAR)
110C-----------------------------------------------
111C M o d u l e s
112C-----------------------------------------------
113 USE timer_mod
114 USE output_mod, only : output_
115 USE mmain_mod
116 USE table_mod
117 USE mat_elem_mod
118 USE nlocal_reg_mod
119 USE ale_connectivity_mod
120 USE sensor_mod
121 USE ale_mod
122 USE dt_mod
123 USE elbufdef_mod
124 USE sdistor_ini_mod, ONLY : sdistor_ini
125 use glob_therm_mod
126 USE sensor_mod
127 use element_mod , only : nixs
128C-----------------------------------------------
129C I m p l i c i t T y p e s
130C-----------------------------------------------
131#include "implicit_f.inc"
132C-----------------------------------------------
133C G l o b a l P a r a m e t e r s
134C-----------------------------------------------
135#include "mvsiz_p.inc"
136C-----------------------------------------------
137C C o m m o n B l o c k s
138C-----------------------------------------------
139#include "com01_c.inc"
140#include "com04_c.inc"
141#include "com08_c.inc"
142#include "vect01_c.inc"
143#include "scr06_c.inc"
144#include "parit_c.inc"
145#include "param_c.inc"
146#include "timeri_c.inc"
147#include "scr18_c.inc"
148C-----------------------------------------------
149C D u m m y A r g u m e n t s
150C-----------------------------------------------
151 TYPE(timer_) ,INTENT(INOUT) :: TIMERS
152 TYPE(OUTPUT_), INTENT(INOUT) :: OUTPUT
153 INTEGER, INTENT(IN) :: S_SFEM_NODVAR
154 INTEGER, INTENT(IN) :: SNPC
155 INTEGER, INTENT(IN) :: STF
156 INTEGER, INTENT(IN) :: NSVOIS
157 INTEGER, INTENT(IN) :: SZ_BUFVOIS
158 INTEGER, INTENT(IN) :: SBUFMAT
159 INTEGER, INTENT(INOUT) :: IDEL7NOK
160 INTEGER, INTENT(IN) :: IDTMINS
161 INTEGER ,INTENT(IN) :: IRESP
162 INTEGER ,INTENT(IN) :: MAXFUNC
163 INTEGER ,INTENT(IN) :: IDEL7NG
164 INTEGER, INTENT(IN) :: IMPL_S
165 INTEGER, INTENT(IN) :: IDYNA
166 INTEGER, INTENT(IN) :: USERL_AVAIL
167 INTEGER IXS(NIXS,*), IPARG(NPARG,NGROUP), NPF(*),IADS(8,*),
168 . IPARTS(*),IPM(NPROPMI,*),IGEO(NPROPGI,*),ITASK,
169 . GRTH(*),IGRTH(*),IOUTPRT
170C
171 INTEGER NELTST,ITYPTST,OFFSET,NEL,NG, ISTRAIN,
172 . iexpan,h3d_strain
173
174 double precision
175 . xdp(3,*)
176
177 my_real
178 . dt2t
179 my_real
180 . pm(npropm,*), geo(npropg,*), x(3,*), a(*), v(*), ms(*), w(*), flux(6,*),
181 . flu1(*), veul(*), fv(*), tf(*), bufmat(*),
182 . partsav(*),stifn(*), fsky(*),eani(*), fskym(*),
183 . f11(mvsiz),f21(mvsiz),f31(mvsiz),
184 . f12(mvsiz),f22(mvsiz),f32(mvsiz),
185 . f13(mvsiz),f23(mvsiz),f33(mvsiz),
186 . f14(mvsiz),f24(mvsiz),f34(mvsiz),d(*),
187 . temp(*), fthe(*), fthesky(*),gresav(*), mssa(*), dmels(*), voln(mvsiz)
188 my_real msnf(*),sfem_nodvar(s_sfem_nodvar),condn(*),condnsky(*),bufvois(6,*)
189 my_real, DIMENSION(MVSIZ,6), INTENT(INOUT) :: svis
190 TYPE(ttable) TABLE(*)
191 TYPE (ELBUF_STRUCT_), TARGET, DIMENSION(NGROUP) :: ELBUF_TAB
192 TYPE (NLOCAL_STR_) , TARGET :: NLOC_DMG
193 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
194 TYPE (MAT_ELEM_) ,INTENT(INOUT) :: MAT_ELEM
195 TYPE (SENSORS_) , INTENT(INOUT) :: SENSORS
196 TYPE(DT_), INTENT(INOUT) :: DT
197 type (glob_therm_) ,intent(inout) :: glob_therm
198C-----------------------------------------------
199C L o c a l V a r i a b l e s
200C-----------------------------------------------
201 INTEGER I,IBID,NF1,IFLAG,IPTR,IPTS,IPTT,ILAY,IMAT
202 INTEGER IBIDON(1),ITET,IP
203 my_real SUM
204 INTEGER MXT(MVSIZ),NGL(MVSIZ),NGEO(MVSIZ)
205 my_real
206 . VD2(MVSIZ) , DVOL(MVSIZ),DELTAX(MVSIZ),
207 . VIS(MVSIZ) , QVIS(MVSIZ), CXX(MVSIZ) ,
208 . S1(MVSIZ) , S2(MVSIZ) , S3(MVSIZ) ,
209 . s4(mvsiz) , s5(mvsiz) , s6(mvsiz) ,
210 . b1(mvsiz) , b2(mvsiz) , b3(mvsiz) ,
211 . b4(mvsiz) , b5(mvsiz) , b6(mvsiz) ,
212 . dxx(mvsiz) , dyy(mvsiz) , dzz(mvsiz) ,
213 . d4(mvsiz) , d5(mvsiz) , d6(mvsiz) ,
214 . rx(mvsiz) , ry(mvsiz) , rz(mvsiz) ,
215 . sx(mvsiz) , sy(mvsiz) , sz(mvsiz) ,
216 . tx(mvsiz) , ty(mvsiz) , tz(mvsiz) ,
217 . vdx(mvsiz), vdy(mvsiz), vdz(mvsiz),ssp_eq(mvsiz),
218 . conde(mvsiz),divde(mvsiz)
219
220 DOUBLE PRECISION
221 . X0(MVSIZ,4),Y0(MVSIZ,4),Z0(MVSIZ,4),
222 . xd1(mvsiz), xd2(mvsiz), xd3(mvsiz), xd4(mvsiz),
223 . yd1(mvsiz), yd2(mvsiz), yd3(mvsiz), yd4(mvsiz),
224 . zd1(mvsiz), zd2(mvsiz), zd3(mvsiz), zd4(mvsiz),voldp(mvsiz)
225
226 ! variables used as arguments by materials
227 my_real
228 . sti(mvsiz), gama(mvsiz,6),
229 . wxx(mvsiz) , wyy(mvsiz) , wzz(mvsiz),aire(mvsiz)
230 ! variables used as arguments by materials (sph only)
231 my_real
232 . muvoid(mvsiz)
233
234 ! Variables void MMAIN
235 my_real
236 . sigy(mvsiz),et(mvsiz),r3_free(mvsiz)
237 ! variables used only in solid routines (as arguments)
238 INTEGER NC1(MVSIZ), NC2(MVSIZ), NC3(MVSIZ), NC4(MVSIZ)
239 my_real
240 . OFF(MVSIZ) , RHOO(MVSIZ),OFFG0(MVSIZ) ,
241 . x1(mvsiz), x2(mvsiz), x3(mvsiz), x4(mvsiz),
242 . y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
243 . z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
244 . vx1(mvsiz),vx2(mvsiz),vx3(mvsiz),vx4(mvsiz),
245 . vy1(mvsiz),vy2(mvsiz),vy3(mvsiz),vy4(mvsiz),
246 . vz1(mvsiz),vz2(mvsiz),vz3(mvsiz),vz4(mvsiz),
247 . px1(mvsiz),px2(mvsiz),px3(mvsiz),px4(mvsiz),
248 . py1(mvsiz),py2(mvsiz),py3(mvsiz),py4(mvsiz),
249 . pz1(mvsiz),pz2(mvsiz),pz3(mvsiz),pz4(mvsiz),
250 . vdx1(mvsiz),vdx2(mvsiz),vdx3(mvsiz),vdx4(mvsiz),
251 . vdy1(mvsiz),vdy2(mvsiz),vdy3(mvsiz),vdy4(mvsiz),
252 . vdz1(mvsiz),vdz2(mvsiz),vdz3(mvsiz),vdz4(mvsiz),
253 . dxy(mvsiz),dyx(mvsiz),
254 . dyz(mvsiz),dzy(mvsiz),
255 . dzx(mvsiz),dxz(mvsiz),
256 . e1x(mvsiz) , e1y(mvsiz) , e1z(mvsiz) ,
257 . e2x(mvsiz) , e2y(mvsiz) , e2z(mvsiz) ,
258 . e3x(mvsiz) , e3y(mvsiz) , e3z(mvsiz) ,rho_0,
259 . tempel(mvsiz), them(mvsiz,4) , die(mvsiz)
260
261 my_real
262 . vx0(mvsiz,4),vy0(mvsiz,4),vz0(mvsiz,4),
263 . mfxx(mvsiz),mfxy(mvsiz),mfyx(mvsiz),
264 . mfyy(mvsiz),mfyz(mvsiz),mfzy(mvsiz),
265 . mfzz(mvsiz),mfzx(mvsiz),mfxz(mvsiz),bid(mvsiz),amu(mvsiz),
266 . fld(mvsiz),sti_c(mvsiz),ll(mvsiz),mu,fqmax
267 my_real,
268 . DIMENSION(:), POINTER :: eint
269 my_real, dimension(mvsiz) :: fheat
270
271 TYPE(g_bufel_) ,POINTER :: GBUF
272 TYPE(L_BUFEL_) ,POINTER :: LBUF
273C----- Variables for non-local computation
274 INTEGER INOD(4),IPOS(4), L_NLOC, INLOC,ISM12_11
275 my_real, DIMENSION(:), ALLOCATABLE :: VAR_REG
276 my_real, DIMENSION(:), POINTER :: DNL
277C
278
279C-----
280c Flag Bolt Preloading
281 INTEGER IBOLTP,NBPRELD,ISCTL,ISTAB(MVSIZ),PID,SZ_IX
282 my_real,
283 . DIMENSION(:), POINTER :: bpreld
284C-----------------------------------------------
285C S o u r c e L i n e s
286C=======================================================================
287 gbuf => elbuf_tab(ng)%GBUF
288 lbuf => elbuf_tab(ng)%BUFLY(1)%LBUF(1,1,1)
289C-----special traitment when switching Ismstr from 12 to 11
290 ism12_11 = elbuf_tab(ng)%BUFLY(1)%L_SIGL
291 ibid = 0
292 ibidon(1) = 0
293
294 tempel(:) = zero
295 fheat(:) = zero
296C
297 sz_ix=numelq+numels+nsvois
298 iboltp = iparg(72,ng)
299 inloc = iparg(78,ng)
300 ALLOCATE(var_reg(nel))
301 nbpreld = gbuf%G_BPRELD
302 bpreld =>gbuf%BPRELD(1:nbpreld*nel)
303C-----------
304 nf1 = nft+1
305 pid = ixs(10,nf1)
306 igtyp = igeo(11,pid)
307 isctl = igeo(97,pid)
308c
309 CALL s4coor3(
310 1 x, ixs(1,nf1),v, w,
311 2 x1, x2, x3, x4,
312 3 y1, y2, y3, y4,
313 4 z1, z2, z3, z4,
314 5 vx1, vx2, vx3, vx4,
315 6 vy1, vy2, vy3, vy4,
316 7 vz1, vz2, vz3, vz4,
317 8 vdx1, vdx2, vdx3, vdx4,
318 9 vdy1, vdy2, vdy3, vdy4,
319 a vdz1, vdz2, vdz3, vdz4,
320 b vdx, vdy, vdz, vd2,
321 c vis, gbuf%OFF, off, gbuf%SMSTR,
322 d gbuf%RHO, rhoo, nc1, nc2,
323 e nc3, nc4, ngl, mxt,
324 f ngeo, f11, f21, f31,
325 g f12, f22, f32, f13,
326 h f23, f33, f14, f24,
327 i f34, xd1, xd2, xd3,
328 j xd4, yd1, yd2, yd3,
329 k yd4, zd1, zd2, zd3,
330 l zd4, xdp, nel, jale,
331 m ismstr, jeul, jlag)
332C-----------
333C GATHER NODAL VARIABLES FOR TOTAL STRAIN CASE.
334C-----------
335 IF ((ismstr >= 10.AND.ismstr <= 12).AND.jlag > 0) THEN
336 CALL sgcoor3(
337 1 tt, 4, x, ixs(1,nf1),
338 2 x0, y0, z0, vx0,
339 3 vy0, vz0, gbuf%SMSTR,d,
340 4 gbuf%OFF, offg0, nel, xdp,
341 5 mtn, ismstr)
342 IF (ismstr == 11) THEN
343 CALL s4derit3(
344 1 off, voln, ngl, deltax,
345 2 mxt, x0(1,1), x0(1,2), x0(1,3),
346 3 x0(1,4), y0(1,1), y0(1,2), y0(1,3),
347 4 y0(1,4), z0(1,1), z0(1,2), z0(1,3),
348 5 z0(1,4), px1, px2, px3,
349 6 px4, py1, py2, py3,
350 7 py4, pz1, pz2, pz3,
351 8 pz4, rx, ry, rz,
352 9 sx, sy, sz, tx,
353 a ty, tz, pm, voldp,
354 b nel, iformdt)
355 ELSE
356 IF (ismstr == 12.AND.idtmin(1)==3.AND.ism12_11==0) THEN
357C!!!!!!calcul local rep for ISMSTR 10 to 11 (offg>un)
358 CALL s4rcoor12(
359 1 gbuf%OFF,nc1, nc2, nc3,
360 2 nc4, x, xdp, d,
361 3 e1x, e2x, e3x, e1y,
362 4 e2y, e3y, e1z, e2z,
363 5 e3z, nel, jcvt)
364 END IF
365 CALL s4derito3(
366 1 off, voln, x0(1,1), x0(1,2),
367 2 x0(1,3), x0(1,4), y0(1,1), y0(1,2),
368 3 y0(1,3), y0(1,4), z0(1,1), z0(1,2),
369 4 z0(1,3), z0(1,4), px1, px2,
370 5 px3, px4, py1, py2,
371 6 py3, py4, pz1, pz2,
372 7 pz3, pz4, rx, ry,
373 8 rz, sx, sy, sz,
374 9 tx, ty, tz, gbuf%JAC_I,
375 a nel)
376 END IF !(ISMSTR == 11) THEN
377 CALL s4defot3(
378 1 px1, px2, px3, px4,
379 2 py1, py2, py3, py4,
380 3 pz1, pz2, pz3, pz4,
381 4 vx0(1,1),vx0(1,2),vx0(1,3),vx0(1,4),
382 5 vy0(1,1),vy0(1,2),vy0(1,3),vy0(1,4),
383 6 vz0(1,1),vz0(1,2),vz0(1,3),vz0(1,4),
384 7 mfxx, mfxy, mfxz, mfyx,
385 8 mfyy, mfyz, mfzx, mfzy,
386 9 mfzz, nel)
387 IF (isorth == 0) THEN
388 DO i=lft,llt
389 gama(i,1) = one
390 gama(i,2) = zero
391 gama(i,3) = zero
392 gama(i,4) = zero
393 gama(i,5) = one
394 gama(i,6) = zero
395 ENDDO
396 ELSE
397 CALL sreploc3(
398 1 rx, ry, rz, sx,
399 2 sy, sz, tx, ty,
400 3 tz, e1x, e2x, e3x,
401 4 e1y, e2y, e3y, e1z,
402 5 e2z, e3z, llt)
403 CALL sorthdir3(
404 1 rx, ry, rz, sx,
405 2 sy, sz, tx, ty,
406 3 tz, e1x, e2x, e3x,
407 4 e1y, e2y, e3y, e1z,
408 5 e2z, e3z, gbuf%GAMA,gama,
409 6 nel, irep)
410 ENDIF
411 IF (ismstr == 12.AND.ism12_11==0.AND.idtmin(1)==3) THEN
412 CALL sordeft12(lft,llt,mfxx, mfxy, mfxz,
413 . mfyx, mfyy, mfyz,
414 . mfzx, mfzy, mfzz,
415 . e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z,gbuf%OFF)
416 ENDIF
417 ENDIF
418C
419 IF(jale+jlag /= 0)THEN
420 IF (ismstr /= 11) THEN
421C!!!!!!also add exception for ismstr10+offg>un
422 CALL s4deri3(
423 1 off, voln, ngl, deltax,
424 2 mxt, xd1, xd2, xd3,
425 3 xd4, yd1, yd2, yd3,
426 4 yd4, zd1, zd2, zd3,
427 5 zd4, px1, px2, px3,
428 6 px4, py1, py2, py3,
429 7 py4, pz1, pz2, pz3,
430 8 pz4, rx, ry, rz,
431 9 sx, sy, sz, tx,
432 a ty, tz, gbuf%SMSTR,gbuf%OFF,
433 b nel, pm, voldp, ismstr,
434 c iformdt, jlag)
435C
436 IF (isorth == 0) THEN
437 DO i=lft,llt
438 gama(i,1) = one
439 gama(i,2) = zero
440 gama(i,3) = zero
441 gama(i,4) = zero
442 gama(i,5) = one
443 gama(i,6) = zero
444 ENDDO
445 ELSE
446 CALL sreploc3(
447 1 rx, ry, rz, sx,
448 2 sy, sz, tx, ty,
449 3 tz, e1x, e2x, e3x,
450 4 e1y, e2y, e3y, e1z,
451 5 e2z, e3z, llt)
452 CALL sorthdir3(
453 1 rx, ry, rz, sx,
454 2 sy, sz, tx, ty,
455 3 tz, e1x, e2x, e3x,
456 4 e1y, e2y, e3y, e1z,
457 5 e2z, e3z, gbuf%GAMA,gama,
458 6 nel, irep)
459 ENDIF
460 END IF !(ISMSTR /= 11) THEN
461C
462 ELSEIF(jeul/=0)THEN
463 CALL e4pxle3(
464 1 gbuf%VOL,veul, x1, x2,
465 2 x3, x4, y1, y2,
466 3 y3, y4, z1, z2,
467 4 z3, z4, px1, px2,
468 5 px3, px4, py1, py2,
469 6 py3, py4, pz1, pz2,
470 7 pz3, pz4, voln, deltax,
471 8 nel, nft)
472 ENDIF
473C
474 CALL s4defo3(
475 1 px1, px2, px3, px4,
476 2 py1, py2, py3, py4,
477 3 pz1, pz2, pz3, pz4,
478 4 vx1, vx2, vx3, vx4,
479 5 vy1, vy2, vy3, vy4,
480 6 vz1, vz2, vz3, vz4,
481 7 dxx, dxy, dxz, dyx,
482 8 dyy, dyz, dzx, dzy,
483 9 dzz, d4, d5, d6,
484 a wxx, wyy, wzz, nel,
485 b ismstr)
486 IF (idtmin(1)==3.AND.ismstr == 12.AND.ism12_11==0) THEN
487 CALL sordef12(lft,llt,dxx, dyy, dzz,
488 . d4, d5, d6,
489 . e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z,offg0)
490 ENDIF
491C
492 CALL s11fx3(
493 1 pm, flux(1,nf1),ale_connect,ixs,
494 2 ipm, bufmat, nel, nft,
495 3 jale, jeul , nummat, numels+nsvois)
496C-----------------------------------------------------
497C SFEM
498C-----------------------------------------------------
499 IF(jlag > 0.AND.isrot == 3) THEN
500 rho_0 = pm(1,mxt(1))
501 CALL s4voln_m(
502 1 sfem_nodvar, nc1, nc2, nc3,
503 2 nc4, mxt, gbuf%OFF, gbuf%RHO,
504 3 rho_0, mfxx, mfxy, mfxz,
505 4 mfyx, mfyy, mfyz, mfzx,
506 5 mfzy, mfzz, gbuf%VOL, voln,
507 6 lbuf%VOL0DP, voldp, gbuf%AMU, dxx,
508 7 dyy, dzz, mat_elem%MAT_PARAM,nel,
509 8 ismstr, s_sfem_nodvar)
510 ENDIF
511 IF(jale > 0 .AND. isrot == 3 .AND. mtn /= 37 .AND. mtn /= 51 .AND.
512 . mtn /= 18 .AND. mtn /= 11) THEN
513 DO i=lft,llt
514 IF(off(i) /= 0) THEN
515 sum=sfem_nodvar(nc1(i))+sfem_nodvar(nc2(i))+sfem_nodvar(nc3(i))+sfem_nodvar(nc4(i))
516 voln(i)=fourth*sum*gbuf%RHO(i)/pm(1,mxt(i))
517 ENDIF
518 ENDDO
519 ENDIF
520C-------------------------------------------------------------
521C DENSITY
522C-------------------------------------------------------------
523 divde(1:nel) = dt1*(dxx(1:nel)+ dyy(1:nel)+ dzz(1:nel))
524 CALL srho3(
525 1 pm, gbuf%VOL, gbuf%RHO, gbuf%EINT,
526 2 divde, flux(1,nf1), flu1(nf1), voln,
527 3 dvol, ngl, mxt, off,
528 4 iparg(64,ng),gbuf%TAG22, voldp, lbuf%VOL0DP,
529 5 amu, gbuf%OFF, nel, mtn,
530 6 jale, ismstr, jeul, jlag)
531C----for Eint compute----
532 IF (ismstr == 12.AND.idtmin(1)==3.AND.ism12_11==0) THEN
533 CALL sroto12_sig(lft,llt,lbuf%SIG,nel,
534 . e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z,offg0)
535!! temporary replaced by (the same) SROTO12_SIG() in order to do not affect
536!! the other multidimensional buffer ARRAYS which are still not modified
537!! CALL SROTO12(LFT,LLT,LBUF%SIG,
538!! . E1X, E1Y, E1Z, E2X, E2Y, E2Z, E3X, E3Y, E3Z,GBUF%OFF)
539 ENDIF
540
541
542 CALL srota3(
543 1 gbuf%SIG,s1, s2, s3,
544 2 s4, s5, s6, wxx,
545 3 wyy, wzz, nel, mtn,
546 4 ismstr)
547C-----------------------------
548C SMALL STRAIN
549C-----------------------------
550 CALL s4malla3(
551 1 gbuf%SMSTR,gbuf%OFF, off, wxx,
552 2 wyy, wzz, nel, ismstr,
553 3 jlag)
554C--------------------------
555C- UPDATE REF CONFIGURATION (possible future change to small strain option)
556C- Total strain option doesn't change the Ref CONF.
557C--------------------------
558 IF (ismstr <= 3.OR.(ismstr==4.AND.jlag>0)) THEN
559 CALL s4sav3(gbuf%OFF,gbuf%SMSTR,
560 . xd1, xd2, xd3, xd4, yd1, yd2, yd3, yd4,
561 . zd1, zd2, zd3, zd4,nel)
562 END IF !(ISMSTR <= 3) THEN
563C-----------------------------
564C heat transfert
565C
566 IF (jthe < 0 ) THEN
567 DO i = lft, llt
568 tempel(i) = fourth*( temp(nc1(i)) + temp(nc2(i)) +
569 . temp(nc3(i)) + temp(nc4(i)) )
570 ENDDO
571 ENDIF
572c-------------------------------------------
573c COMPUTE Regularized non local variable in Gauss point
574c-------------------------------------------
575 IF (inloc > 0) THEN
576 l_nloc = nloc_dmg%L_NLOC
577 dnl => nloc_dmg%DNL(1:l_nloc) ! DNL = non local variable increment
578 DO i=lft,llt
579 inod(1) = nloc_dmg%IDXI(nc1(i))
580 inod(2) = nloc_dmg%IDXI(nc2(i))
581 inod(3) = nloc_dmg%IDXI(nc3(i))
582 inod(4) = nloc_dmg%IDXI(nc4(i))
583 ipos(1) = nloc_dmg%POSI(inod(1))
584 ipos(2) = nloc_dmg%POSI(inod(2))
585 ipos(3) = nloc_dmg%POSI(inod(3))
586 ipos(4) = nloc_dmg%POSI(inod(4))
587 var_reg(i) = fourth*(dnl(ipos(1)) + dnl(ipos(2)) + dnl(ipos(3)) + dnl(ipos(4)))
588 ENDDO
589 ENDIF
590C------------------------------------------------------
591C CALCUL DES CONTRAINTES SUIVANT LOIS CONSTITUTIVES
592C------------------------------------------------------
593Cfor fluid in spmd: bufvois loi11
594 IF ((itask==0).AND.(imon_mat==1)) CALL startime(timers,35)
595
596 ilay = 1
597 iptr = 1
598 ipts = 1
599 iptt = 1
600 ip = 1
601
602 IF(iboltp /= 0) THEN
603 CALL srota3(
604 1 bpreld(3*nel+1),b1, b2, b3,
605 2 b4, b5, b6, wxx,
606 3 wyy, wzz, nel, mtn,
607 4 ismstr)
608
609 CALL boltst(
610 1 ip, bpreld, lbuf%SIG,tt, nel,
611 2 npt, sensors%NSENSOR,sensors%SENSOR_TAB,
612 3 iparg(67,ng),iparg(68,ng))
613 END IF
614
615 CALL mmain(timers,output,
616 1 elbuf_tab, ng, pm, geo,
617 2 ale_connect, ixs, iparg,
618 3 v, tf, npf, bufmat,
619 4 sti, x, dt2t, neltst,
620 5 ityptst, offset, nel, w,
621 6 off, ngeo, mxt, ngl,
622 7 voln, vd2, dvol, deltax,
623 8 vis, qvis, cxx, s1,
624 9 s2, s3, s4, s5,
625 a s6, dxx, dyy, dzz,
626 b d4, d5, d6, wxx,
627 c wyy, wzz, rx, ry,
628 d rz, sx, sy, sz,
629 e vdx, vdy, vdz, muvoid,
630 f ssp_eq, aire, sigy, et,
631 g bufvois, lbuf%PLA, r3_free, amu,
632 h mfxx, mfxy, mfxz, mfyx,
633 i mfyy, mfyz, mfzx, mfzy,
634 j mfzz, ipm, gama, bid,
635 k bid, bid, bid, bid,
636 l bid, bid, istrain, tempel,
637 m die, iexpan, ilay, mssa,
638 n dmels, iptr, ipts, iptt,
639 o table, bid, bid, bid,
640 p bid, iparg(1,ng), igeo, conde,
641 q itask, nloc_dmg, var_reg, mat_elem,
642 r h3d_strain, jplasol, jsph, sz_bufvois,
643 t snpc, stf, sbufmat, glob_therm,
644 u svis, sz_ix, iresp,
645 v n2d, th_strain, ngroup, tt,
646 . dt1, ntable, numelq, nummat,
647 . numgeo, numnod, numels,
648 . idel7nok, idtmin, maxfunc,
649 . imon_mat, userl_avail, impl_s,
650 . idyna, dt, fheat ,sensors)
651c
652 IF ((itask==0).AND.(imon_mat==1)) CALL stoptime(timers,35)
653C----------- return to global system
654 IF (ismstr == 12.AND.idtmin(1)==3) THEN
655 IF (ism12_11==0) THEN
656 CALL sroto12_sig(lft,llt,lbuf%SIG,nel,
657 . e1x,e2x,e3x,e1y,e2y,e3y,e1z,e2z,e3z,offg0)
658!! temporary replaced by (the same) SROTO12_SIG() in order to do not affect
659!! the other multidimensional buffer ARRAYS which are still not modified
660!! CALL SROTO12(LFT,LLT,LBUF%SIG,
661!! . E1X,E2X,E3X,E1Y,E2Y,E3Y,E1Z,E2Z,E3Z,GBUF%OFF)
662 IF (istrain == 1) THEN
663 CALL sordef12(lft,llt,dxx, dxy, dxz,
664 . d4, d5, d6,
665 . e1x,e2x,e3x,e1y,e2y,e3y,e1z,e2z,e3z,offg0)
666 ENDIF
667 END IF
668 ENDIF
669 IF (istrain == 1) THEN
670 CALL sstra3(
671 1 dxx, dyy, dzz, d4,
672 2 d5, d6, lbuf%STRA,wxx,
673 3 wyy, wzz, off, nel,
674 4 jcvt)
675 ENDIF
676C--------------------------
677C balances by material
678C--------------------------
679 iflag=mod(ncycle,ncpri)
680 IF(ioutprt>0)THEN
681 IF (mtn == 11) THEN
682 eint => elbuf_tab(ng)%GBUF%EINS(1:nel)
683 ELSE
684 eint => elbuf_tab(ng)%GBUF%EINT(1:nel)
685 ENDIF
686 CALL s4bilan(partsav,eint,gbuf%RHO,gbuf%RK,gbuf%VOL,
687 . vx1, vx2, vx3, vx4, vy1, vy2, vy3, vy4,
688 . vz1, vz2, vz3, vz4, voln,iparts,gresav,
689 . grth,igrth,iexpan,gbuf%EINTTH,gbuf%FILL,
690 . x1, x2, x3, x4, y1, y2, y3, y4,
691 . z1, z2, z3, z4,itask,iparg(1,ng),off,sensors,
692 . nel,gbuf%G_WPLA,gbuf%WPLA)
693 ENDIF
694C
695 IF(jlag+jale+jeul == 0)RETURN
696 itet = 1
697 bid(lft:llt)=zero
698 CALL sgeodel3(ngl,gbuf%OFF,voln,deltax,gbuf%VOL,geo(1,ngeo(1)),bid,dt,nel,idel7nok )
699C-----------------------------
700C SMALL STRAIN
701C-----------------------------
702 CALL smallb3(gbuf%OFF,off, nel, ismstr)
703 CALL smallgeo3(ngl, gbuf%OFF ,voln ,deltax, gbuf%VOL ,itet, nel,ismstr,dt)
704 IF (ismstr == 12.AND.idtmin(1)==3) THEN
705 CALL s4sav12(
706 1 gbuf%OFF, offg0, gbuf%SMSTR,xd1,
707 2 xd2, xd3, xd4, yd1,
708 3 yd2, yd3, yd4, zd1,
709 4 zd2, zd3, zd4, nel)
710 IF (ism12_11>0 .AND. isorth == 0) THEN
711 CALL s4upd11t12(
712 1 gbuf%OFF, offg0, xd1, xd2,
713 2 xd3, xd4, yd1, yd2,
714 3 yd3, yd4, zd1, zd2,
715 4 zd3, zd4, gbuf%JAC_I,gbuf%SIG,
716 5 lbuf%SIGL, nel, jcvt)
717 END IF
718 END IF
719C--------------------------
720C update of masses
721C TRANSPORT FORCES
722C----------------------------
723 IF (jale+jeul > 0 .AND. ale%GLOBAL%INCOMP == 0)THEN
724 IF(iparit == 0)THEN
725 CALL a4mass3(
726 1 ms, gbuf%RHO,voln, nc1,
727 2 nc2, nc3, nc4, msnf,
728 3 off, nel)
729 ELSE
730 CALL a4mass3p(
731 1 fskym, gbuf%RHO,voln, iads,
732 2 off, nel, nft)
733 ENDIF
734 ENDIF
735C-----------------------------
736C TRANSPORT FORCES
737C-----------------------------
738 IF(jale == 1 .OR. jeul == 1)THEN
739C
740 CALL a4momt3(
741 1 pm, gbuf%RHO,voln, x1,
742 2 x2, x3, x4, y1,
743 3 y2, y3, y4, z1,
744 4 z2, z3, z4, vx1,
745 5 vx2, vx3, vx4, vy1,
746 6 vy2, vy3, vy4, vz1,
747 7 vz2, vz3, vz4, f11,
748 8 f21, f31, f12, f22,
749 9 f32, f13, f23, f33,
750 a f14, f24, f34, px1,
751 b px2, px3, px4, py1,
752 c py2, py3, py4, pz1,
753 d pz2, pz3, pz4, dxx,
754 e dxy, dxz, dyx, dyy,
755 f dyz, dzx, dzy, dzz,
756 g vdx1, vdx2, vdx3, vdx4,
757 h vdy1, vdy2, vdy3, vdy4,
758 i vdz1, vdz2, vdz3, vdz4,
759 j vdx, vdy, vdz, deltax,
760 k vis, mxt, rx, ry,
761 l rz, sx, sy, sz,
762 m tx, ty, tz, nel,
763 n mtn)
764 ENDIF
765
766 IF(jeul+jale /= 0) CALL check_off_ale(f11,f21,f31,f12,f22,
767 1 f32,f13,f23,f33,f14,
768 2 f24,f34,bid,bid,bid,
769 3 bid,bid,bid,bid,bid,
770 4 bid,bid,bid,bid,gbuf%OFF,
771 5 lft,llt,nel)
772
773C----------------------------
774C INTERNAL FORCES
775C----------------------------
776 CALL s4fint3(gbuf%SIG,
777 . px1, px2, px3, px4,
778 . py1, py2, py3, py4,
779 . pz1, pz2, pz3, pz4,
780 . f11,f21,f31,f12,f22,f32,f13,f23,f33,f14,f24,f34,
781 . voln,qvis,nel,svis)
782C --------------------------
783C --- heat transfert
784C --------------------------
785 IF (jthe < 0 ) THEN
786 imat = mxt(1)
787 IF (mat_elem%MAT_PARAM(imat)%HEAT_FLAG == 1) THEN
788 CALL s4therm (
789 1 pm, imat, voln, nc1,
790 2 nc2, nc3, nc4, px1,
791 3 px2, px3, px4, py1,
792 4 py2, py3, py4, pz1,
793 5 pz2, pz3, pz4, dt1,
794 6 temp, tempel, fheat, them,
795 7 gbuf%OFF,lbuf%OFF,nel,glob_therm%THEACCFACT)
796 ELSE
797 CALL s4therm (
798 1 pm, imat, voln, nc1,
799 2 nc2, nc3, nc4, px1,
800 3 px2, px3, px4, py1,
801 4 py2, py3, py4, pz1,
802 5 pz2, pz3, pz4, dt1,
803 6 temp, tempel, die, them,
804 7 gbuf%OFF,lbuf%OFF,nel,glob_therm%THEACCFACT)
805 END IF
806
807 ENDIF
808c-------------------------
809c Virtual internal forces of regularized non local ddl
810c--------------------------
811 IF (inloc > 0) THEN
812 CALL s4fint_reg(
813 1 nloc_dmg,var_reg, nel, lbuf%OFF,
814 2 voln, nc1, nc2, nc3,
815 3 nc4, px1, px2, px3,
816 4 px4, py1, py2, py3,
817 5 py4, pz1, pz2, pz3,
818 6 pz4, mxt(lft),itask, dt2t,
819 7 gbuf%VOL,nft)
820 ENDIF
821C----------------------------
822C distortion control
823C----------------------------
824 IF (isctl > 0) THEN
825 CALL sdistor_ini(
826 1 nel ,sti_c ,npropm ,nummat ,
827 2 ismstr ,mxt ,istab ,pm ,
828 3 gbuf%SIG ,gbuf%RHO ,cxx ,off ,
829 4 gbuf%OFF ,ll ,voln ,fld ,
830 5 mu ,fqmax )
831 IF (ismstr==1.OR.ismstr>=11) THEN ! XX,YY,ZZ now current conf.
832#include "vectorize.inc"
833 DO i=1,nel
834 x1(i) =x(1,nc1(i))
835 y1(i) =x(2,nc1(i))
836 z1(i) =x(3,nc1(i))
837 x2(i) =x(1,nc2(i))
838 y2(i) =x(2,nc2(i))
839 z2(i) =x(3,nc2(i))
840 x3(i) =x(1,nc3(i))
841 y3(i) =x(2,nc3(i))
842 z3(i) =x(3,nc3(i))
843 x4(i) =x(1,nc4(i))
844 y4(i) =x(2,nc4(i))
845 z4(i) =x(3,nc4(i))
846 ENDDO
847 ENDIF
848 CALL s4for_distor(
849 . x1, x2, x3, x4,
850 . y1, y2, y3, y4,
851 . z1, z2, z3, z4,
852 . vx1, vx2, vx3, vx4,
853 . vy1, vy2, vy3, vy4,
854 . vz1, vz2, vz3, vz4,
855 . f11, f12, f13, f14,
856 . f21, f22, f23, f24,
857 . f31, f32, f33, f34,
858 . sti, fld, sti_c, ll,
859 . mu, fqmax, istab, nel ,
860 . gbuf%EINT_DISTOR,dt1 )
861 ENDIF
862C --------------------------
863 IF(nfilsol/=0) CALL s4fillopt(
864 1 gbuf%FILL,sti, f11, f21,
865 2 f31, f12, f22, f32,
866 3 f13, f23, f33, f14,
867 4 f24, f34, nel)
868C----------------------------
869 IF (iparit == 0) THEN
870 CALL s4cumu3(
871 1 gbuf%OFF,a, nc1, nc2,
872 2 nc3, nc4, stifn, sti,
873 3 f11, f21, f31, f12,
874 4 f22, f32, f13, f23,
875 5 f33, f14, f24, f34,
876 6 them, fthe, condn, conde,
877 7 nel, jthe, glob_therm%NODADT_THERM)
878 ELSE
879 CALL s4cumu3p(
880 1 gbuf%OFF,sti, fsky, fsky,
881 2 iads, f11, f21, f31,
882 3 f12, f22, f32, f13,
883 4 f23, f33, f14, f24,
884 5 f34, them, fthesky, condnsky,
885 6 conde, nel, nft, jthe,
886 7 glob_therm%NODADT_THERM)
887 ENDIF
888 IF (ALLOCATED(var_reg)) DEALLOCATE(var_reg)
889 RETURN
890 END
a4mass3
subroutine a4mass3(ms, rho, volu, nc1, nc2, nc3, nc4, msnf, off, nel)
Definition a4mass3.F:33
a4mass3p
subroutine a4mass3p(fskym, rho, volu, iads, off, nel, nft)
Definition a4mass3p.F:32
a4momt3
subroutine a4momt3(pm, rho, vol, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, vx1, vx2, vx3, vx4, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, vdx1, vdx2, vdx3, vdx4, vdy1, vdy2, vdy3, vdy4, vdz1, vdz2, vdz3, vdz4, vdx, vdy, vdz, deltax, vis, mat, rx, ry, rz, sx, sy, sz, tx, ty, tz, nel, mtn)
Definition a4momt3.F:57
boltst
subroutine boltst(ip, bpreld, sig, tt, nel, npt, nsensor, sensor_tab, fun_id, sens_id)
Definition boltst.F:39
check_off_ale
subroutine check_off_ale(f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, off, lft, llt, nel)
Definition check_off_ale.F:38
my_real
#define my_real
Definition cppsort.cpp:32
e4pxle3
subroutine e4pxle3(vol, veul, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, det, deltax, nel, nft)
Definition e4pxle3.F:37
mmain
subroutine mmain(pm, elbuf_str, ix, nix, x, geo, iparg, nel, skew, bufmat, ipart, ipartel, nummat, matparam, imat, ipm, ngl, pid, npf, tf, mfxx, mfxy, mfxz, mfyx, mfyy, mfyz, mfzx, mfzy, mfzz, rx, ry, rz, sx, sy, sz, gama, voln, dvol, s1, s2, s3, s4, s5, s6, dxx, dyy, dzz, d4, d5, d6, wxx, wyy, wzz)
Definition mmain.F:43
ale_connectivity_mod
Definition ale_connectivity_mod.F:225
ale_mod
Definition ale_mod.F:150
ale_mod::ale
type(ale_) ale
Definition ale_mod.F:253
dt_mod
Definition dt_mod.F:80
nlocal_reg_mod
Definition nlocal_reg_mod.F:119
table_mod
Definition table_mod.F:113
s11fx3
subroutine s11fx3(pm, flux, ale_connect, ixs, ipm, bufmat, nel, nft, jale, jeul, nummat, s_ixs)
Definition s11fx3.F:36
s4bilan
subroutine s4bilan(partsav, eint, rho, rk, vol, vx1, vx2, vx3, vx4, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, vnew, iparts, gresav, grth, igrth, iexpan, eintth, fill, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, itask, iparg, offg, sensors, nel, g_wpla, wpla)
Definition s4bilan.F:40
s4cumu3
subroutine s4cumu3(offg, e, nc1, nc2, nc3, nc4, stifn, sti, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, them, fthe, condn, conde, nel, jthe, nodadt_therm)
Definition s4cumu3.F:37
s4cumu3p
subroutine s4cumu3p(offg, sti, fsky, fskyv, iads, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, them, fthesky, condnsky, conde, nel, nft, jthe, nodadt_therm)
Definition s4cumu3p.F:37
s4derit3
subroutine s4derit3(off, det, ngl, deltax, mxt, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, rx, ry, rz, sx, sy, sz, tx, ty, tz, pm, voldp, nel, iformdt)
Definition s4derit3.F:43
s4derito3
subroutine s4derito3(off, det, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, rx, ry, rz, sx, sy, sz, tx, ty, tz, jac_i, nel)
Definition s4derito3.F:39
s4fillopt
subroutine s4fillopt(fill, sti, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, nel)
Definition s4fillopt.F:33
s4fint3
subroutine s4fint3(sig, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, vol, qvis, nel, svis)
Definition s4fint3.F:34
s4fint_reg
subroutine s4fint_reg(nloc_dmg, var_reg, nel, off, vol, nc1, nc2, nc3, nc4, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, imat, itask, dt2t, vol0, nft)
Definition s4fint_reg.F:38
s4for_distor
subroutine s4for_distor(x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, vx1, vx2, vx3, vx4, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, f11, f12, f13, f14, f21, f22, f23, f24, f31, f32, f33, f34, sti, fld, sti_c, ll, mu, fqmax, istab, nel, e_distor, dt1)
Definition s4for_distor.F:44
s4forc3
subroutine s4forc3(timers, output, elbuf_tab, ng, pm, geo, ixs, x, a, v, ms, w, flux, flu1, veul, fv, ale_connect, iparg, tf, npf, bufmat, partsav, nloc_dmg, dt2t, neltst, ityptst, stifn, fsky, iads, offset, eani, iparts, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, nel, fskym, msnf, ipm, igeo, bufvois, istrain, itask, temp, fthe, fthesky, iexpan, gresav, grth, igrth, mssa, dmels, table, xdp, sfem_nodvar, voln, condn, condnsky, d, sensors, ioutprt, mat_elem, h3d_strain, dt, idel7nok, nsvois, sz_bufvois, snpc, stf, sbufmat, svis, idtmins, iresp, idel7ng, maxfunc, userl_avail, glob_therm, impl_s, idyna, s_sfem_nodvar)
Definition s4forc3.F:110
s4malla3
subroutine s4malla3(sav, offg, off, wxx, wyy, wzz, nel, ismstr, jlag)
Definition s4mall3.F:32
s4rcoor12
subroutine s4rcoor12(off, nc1, nc2, nc3, nc4, x, xdp, d, r11, r12, r13, r21, r22, r23, r31, r32, r33, nel, jcvt)
Definition s4rcoor12.F:36
s4defot3
subroutine s4defot3(px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, vx1, vx2, vx3, vx4, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz)
Definition s4refsta3.F:562
s4sav12
subroutine s4sav12(offg, offg0, sav, xd1, xd2, xd3, xd4, yd1, yd2, yd3, yd4, zd1, zd2, zd3, zd4, nel)
Definition s4sav12.F:33
s4sav3
subroutine s4sav3(offg, sav, xd1, xd2, xd3, xd4, yd1, yd2, yd3, yd4, zd1, zd2, zd3, zd4, nel)
Definition s4sav3.F:31
s4therm
subroutine s4therm(pm, imat, vol, nc1, nc2, nc3, nc4, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, dt1, tempnc, tel, heat, fphi, offg, off, nel, theaccfact)
Definition s4therm.F:36
s4upd11t12
subroutine s4upd11t12(offg, offg0, xd1, xd2, xd3, xd4, yd1, yd2, yd3, yd4, zd1, zd2, zd3, zd4, jac_1, sig, sigl, nel, jcvt)
Definition s4upd11t12.F:36
s4voln_m
subroutine s4voln_m(sfem_nodvar, nc1, nc2, nc3, nc4, mat, offg, rho, rho0, fxx, fxy, fxz, fyx, fyy, fyz, fzx, fzy, fzz, vol0, voln, vol0dp, voldp, amu0, dxx, dyy, dzz, matparam, nel, ismstr, s_sfem_nodvar)
Definition s4voln_m.F:41
sgcoor3
subroutine sgcoor3(time, npe, x, ixs, x0, y0, z0, vx0, vy0, vz0, sav, d, off, off0, nel, xdp, mtn, ismstr)
Definition sgcoor3.F:40
sgeodel3
subroutine sgeodel3(ngl, offg, volg, deltax, volg0, geo, l_max, dt, nel, idel7nok)
Definition sgeodel3.F:42
smallb3
subroutine smallb3(offg, off, nel, ismstr)
Definition smallb3.F:45
smallgeo3
subroutine smallgeo3(ngl, offg, volg, deltax, volg0, itet, nel, ismstr, dt)
Definition smallgeo3.F:41
sordef12
subroutine sordef12(jft, jlt, dxx, dyy, dzz, d4, d5, d6, g1x, g1y, g1z, g2x, g2y, g2z, g3x, g3y, g3z, off)
Definition sordef12.F:35
sordeft12
subroutine sordeft12(jft, jlt, mxx, mxy, mxz, myx, myy, myz, mzx, mzy, mzz, g1x, g1y, g1z, g2x, g2y, g2z, g3x, g3y, g3z, off)
Definition sordeft12.F:35
sorthdir3
subroutine sorthdir3(rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, gama0, gama, nel, irep)
Definition sorthdir3.F:43
srota3
subroutine srota3(sig, s1, s2, s3, s4, s5, s6, wxx, wyy, wzz, nel, mtn, ismstr)
Definition srota3.F:43
sroto12_sig
subroutine sroto12_sig(jft, jlt, sig, nel, g1x, g1y, g1z, g2x, g2y, g2z, g3x, g3y, g3z, off)
Definition sroto12_sig.F:33
sstra3
subroutine sstra3(dxx, dyy, dzz, d4, d5, d6, strain, wxx, wyy, wzz, off, nel, jcvt)
Definition sstra3.F:47
s4coor3
subroutine s4coor3(x, xrefs, ixs, ngl, mxt, ngeo, ix1, ix2, ix3, ix4, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4)
Definition s4coor3.F:40
s4defo3
subroutine s4defo3(px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, vx1, vx2, vx3, vx4, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, dxx, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, d4, d5, d6, wxx, wyy, wzz)
Definition s4defo3.F:34
s4deri3
subroutine s4deri3(vol, veul, geo, igeo, rx, ry, rz, sx, sy, sz, tx, ty, tz, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, jac_i, deltax, det, ngl, ngeo, mxt, pm, voldp)
Definition s4deri3.F:47
sreploc3
subroutine sreploc3(rx, ry, rz, sx, sy, sz, tx, ty, tz, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z)
Definition sreploc3.F:32
srho3
subroutine srho3(pm, volo, rhon, eint, dxx, dyy, dzz, voln, dvol, mat)
Definition srho3.F:31
startime
subroutine startime(event, itask)
Definition timer.F:93
stoptime
subroutine stoptime(event, itask)
Definition timer.F:135
table_mod::ttable
Definition table_mod.F:126