33!||
sforc3 ../engine/source/elements/solid/solide/
sforc3.f
37 1 OFFG, STI, FSKY, FSKYV,
38 2 IADS, F11, F21, F31,
46 A NC8, AR, FR_WAVE, FR_WAV,
53 H THEM, FTHESKY, CONDNSKY,CONDE,
54 I NEL, NFT, JTHE, ISROT,
55 J IPARTSPH ,NODADT_THERM)
59#include "implicit_f.inc"
74 INTEGER,
INTENT(IN) :: NEL
75 INTEGER,
INTENT(IN) :: NFT
76 INTEGER,
INTENT(IN) :: JTHE
77 INTEGER,
INTENT(IN) :: ISROT
78 INTEGER,
INTENT(IN) :: IPARTSPH
79 INTEGER,
INTENT(IN) :: NODADT_THERM
81 INTEGER NC1(*), NC2(*), NC3(*), NC4(*), NC5(*), NC6(*), NC7(*),
84 . OFFG(*),FSKYV(LSKY,8),FSKY(8,LSKY),STI(*),
85 . F11(*),F21(*),F31(*),F12(*),F22(*),F32(*),
86 . F13(*),F23(*),F33(*),F14(*),F24(*),F34(*),
87 . F15(*),F25(*),F35(*),F16(*),F26(*),F36(*),
88 . F17(*),F27(*),F37(*),(*),F28(*),F38(*),
89 . AR(3,*),FR_WAVE(*),FR_WAV(*),
90 . MX1(*),MY1(*),MZ1(*),MX2(*),MY2(*),MZ2(*),
91 . MX3(*),MY3(*),MZ3(*),MX4(*),MY4(*),MZ4(*),
92 . MX5(*),MY5(*),MZ5(*),MX6(*),MY6(*),MZ6(*),
93 . mx7(*),my7(*),mz7(*),mx8(*),my8(*),mz8(*),
94 . them(mvsiz,8),fthesky(*),condnsky(*),conde(*)
105 off_l =
min(off_l,offg(i))
141 IF(offg(i)<=zero)
THEN
177 IF(offg(i)<=zero)
THEN
183 IF(nodadt_therm == 1)
THEN
199 IF(nodadt_therm == 1 )
THEN
201 conde(i)=one_over_8*conde(i)
207#include "vectorize.inc"
315#include "vectorize.inc"
323 fthesky(k) = them(i,1)
330 fthesky(k) = them(i,7)
337 fthesky(k) = them(i,2)
344 fthesky(k) = them(i,8)
351 fthesky(k) = them(i,3)
358 fthesky(k) = them(i,5)
365 fthesky(k) = them(i,4)
372 fthesky(k) = them(i,6)
375 IF(nodadt_therm == 1)
THEN
383 fthesky(k) = them(i,1)
384 condnsky(k) = conde(i)
391 fthesky(k) = them(i,7)
392 condnsky(k) = conde(i)
399 fthesky(k) = them(i,2)
400 condnsky(k) = conde(i)
407 fthesky(k) = them(i,8)
408 condnsky(k) = conde(i)
415 condnsky(k) = conde(i)
422 fthesky(k) = them(i,5)
423 condnsky(k) = conde(i)
430 fthesky(k) = them(i,4)
431 condnsky(k) = conde(i)
438 fthesky(k) = them(i,6)
439 condnsky(k) = conde(i)
449 fthesky(k) = them(i,1)
456 fthesky(k) = them(i,7)
463 fthesky(k) = them(i,2)
477 fthesky(k) = them(i,3)
484 fthesky(k) = them(i,5)
491 fthesky(k) = them(i,4)
498 fthesky(k) = them(i,6)
537#include "vectorize.inc"
590#include "vectorize.inc"
649 IF(fr_wave(nc1(i))==0.0)fr_wave(nc1(i))=-fr_wav(i)
650 IF(fr_wave(nc2(i))==0.0)fr_wave(nc2(i))=-fr_wav(i)
651 IF(fr_wave(nc3(i))==0.0)fr_wave(nc3(i))=-fr_wav(i)
652 IF(fr_wave(nc4(i))==0.0)fr_wave(nc4(i))=-fr_wav(i)
653 IF(fr_wave(nc5(i))==0.0)fr_wave(nc5(i))=-fr_wav(i)
654 IF(fr_wave(nc6(i))==0.0)fr_wave(nc6(i))=-fr_wav(i)
655 IF(fr_wave(nc7(i))==0.0)fr_wave(nc7(i))=-fr_wav(i)
656 IF(fr_wave(nc8(i))==0.0)fr_wave(nc8(i))=-fr_wav(i)
658#include "lockoff.inc"
subroutine s8cforc3(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, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nel, icp, icsig, smr, sms, smt, mfxx, mfxy, mfxz, mfyx, mfyy, mfyz, mfzx, mfzy, mfzz, nvc, ipm, itask, istrain, temp, fthe, fthesky, iexpan, igeo, nnpt, gresav, grth, igrth, mssa, dmels, table, xdp, voln, condn, condnsky, sensors, ioutprt, mat_elem, h3d_strain, dt, nodadt, dtfac1, dtmin1, idtmin, snpc, stf, sbufmat, svis, nsvois, idtmins, iresp, maxfunc, userl_avail, glob_therm, impl_s, idyna)
subroutine s8sforc3(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, icp, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nel, nvc, ipm, itask, istrain, temp, fthe, fthesky, iexpan, gresav, grth, igrth, mssa, dmels, table, igeo, xdp, voln, condn, condnsky, d, ioutprt, mat_elem, h3d_strain, snpc, stf, sbufmat, svis, nsvois, idtmins, iresp, idtmin, maxfunc, userl_avail, dt, glob_therm, sensors)
subroutine scumu3p(offg, sti, fsky, fskyv, iads, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, ar, fr_wave, fr_wav, mx1, my1, mz1, mx2, my2, mz2, mx3, my3, mz3, mx4, my4, mz4, mx5, my5, mz5, mx6, my6, mz6, mx7, my7, mz7, mx8, my8, mz8, them, fthesky, condnsky, conde, nel, nft, jthe, isrot, ipartsph, nodadt_therm)
subroutine sforc3(timers, output, elbuf_tab, ng, pm, geo, ixs, x, nv46, a, v, ms, w, flux, flu1, veul, fv, ale_connect, iparg, tf, npf, bufmat, partsav, itab, dt2t, neltst, ityptst, stifn, fsky, iads, offset, eani, iparts, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nel, fskym, msnf, isky, fskyi, nvc, ipm, igeo, ar, vr, in, fr_wave, dr, bufvois, itask, qmv, istrain, temp, fthe, fthesky, iexpan, gresav, grth, igrth, mssa, dmels, table, phi1, phi2, vf, af, df, wf, ffsky, afglob, msf, iparg1, xdp, por, icontact, ifoam, voln, condn, condnsky, agrav, igrv, lgrav, sensors, skew, nale, d, ioutprt, nloc_dmg, mat_elem, h3d_strain, dt, idel7nok, nsvois, sz_bufvois, snpc, stf, sbufmat, svis, idtmins, iresp, idel7ng, maxfunc, userl_avail, glob_therm, impl_s, idyna, wfext)