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