OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
sdlen3.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/.
23C
24!||====================================================================
25!|| sdlen3 ../starter/source/elements/solid/solide/sdlen3.F
26!||--- called by ------------------------------------------------------
27!|| inirig_mat ../starter/source/elements/initia/inirig_mat.f
28!|| inivoid ../starter/source/elements/initia/inivoid.F
29!|| multifluid_init3 ../starter/source/multifluid/multifluid_init3.F
30!|| scinit3 ../starter/source/elements/thickshell/solidec/scinit3.F
31!|| sinit3 ../starter/source/elements/solid/solide/sinit3.F
32!|| suinit3 ../starter/source/elements/elbuf_init/suinit3.F
33!||--- calls -----------------------------------------------------
34!|| slen ../starter/source/elements/solid/solide/slen.F
35!||====================================================================
36 SUBROUTINE sdlen3(
37 . X1 ,X2 ,X3 ,X4 ,X5 ,X6 ,X7 ,X8 ,
38 . Y1 ,Y2 ,Y3 ,Y4 ,Y5 ,Y6 ,Y7 ,Y8 ,
39 . Z1 ,Z2 ,Z3 ,Z4 ,Z5 ,Z6 ,Z7 ,Z8 ,
40 . DELTAX, VOLN)
41C
42C CHARACTERISTIC LENGTH : DX = V / Smax
43C
44C-----------------------------------------------
45C I m p l i c i t T y p e s
46C-----------------------------------------------
47#include "implicit_f.inc"
48C-----------------------------------------------
49C G l o b a l P a r a m e t e r s
50C-----------------------------------------------
51#include "mvsiz_p.inc"
52C-----------------------------------------------
53C C o m m o n B l o c k s
54C-----------------------------------------------
55#include "vect01_c.inc"
56C-----------------------------------------------
57C L o c a l V a r i a b l e s
58C-----------------------------------------------
59 INTEGER I, NFAC
60 my_real
61 . VOLN(*), XIOFF(MVSIZ), AREAM(MVSIZ),
62 . ATEST(MVSIZ), AREA(6,MVSIZ),
63 . x1(*), x2(*), x3(*), x4(*), x5(*), x6(*), x7(*), x8(*), y1(*),
64 . y2(*), y3(*), y4(*), y5(*), y6(*), y7(*), y8(*), z1(*), z2(*),
65 . z3(*), z4(*), z5(*), z6(*), z7(*), z8(*), deltax(*)
66C
67C---------------------------------------------------------------------
68 DO i=lft,llt
69 xioff(i) = one
70 aream(i) = zero
71 ENDDO
72
73 CALL slen(x1,x2,x3,x4,y1,y2,y3,y4,z1,z2,z3,z4,1, area, aream)
74 CALL slen(x5,x6,x7,x8,y5,y6,y7,y8,z5,z6,z7,z8,2, area, aream)
75 CALL slen(x1,x2,x6,x5,y1,y2,y6,y5,z1,z2,z6,z5,3, area, aream)
76 CALL slen(x2,x3,x7,x6,y2,y3,y7,y6,z2,z3,z7,z6,4, area, aream)
77 CALL slen(x3,x4,x8,x7,y3,y4,y8,y7,z3,z4,z8,z7,5, area, aream)
78 CALL slen(x4,x1,x5,x8,y4,y1,y5,y8,z4,z1,z5,z8,6, area, aream)
79
80 DO i=lft,llt
81 atest(i) = em4*aream(i)
82 ENDDO
83
84 DO i=lft,llt
85 nfac=0
86 IF(area(1,i)<atest(i)) nfac=nfac+1
87 IF(area(2,i)<atest(i)) nfac=nfac+1
88 IF(area(3,i)<atest(i)) nfac=nfac+1
89 IF(area(4,i)<atest(i)) nfac=nfac+1
90 IF(area(5,i)<atest(i)) nfac=nfac+1
91 IF(area(6,i)<atest(i)) nfac=nfac+1
92 IF(nfac>=3) xioff(i)=ep03
93 ENDDO
94
95 DO i=lft,llt
96 deltax(i) = four*voln(i)*xioff(i)/sqrt(aream(i))
97 ENDDO
98
99 RETURN
100 END
subroutine inirig_mat(ixc, ixs, ixtg, ixs10, x, v, pm, geo, ms, in, ptg, msc, mss, mstg, inc, intg, thkc, thkt, partsav, iparts, ipartc, ipartt, veul, dtelem, ihbe, isolnod, nvc, i8mi, msnf, mssf, igeo, etnod, nshnod, stc, sttg, wma, sh4tree, sh3tree, mcp, mcpc, temp, mcps, mssx, mcpsx, ins, stifn, stifr, connec, irig_node, numel, nindx, xrefc, xreftg, xrefs, mssa, sh3trim, isubstack, bufmat, ipm, stack, rnoise, strc, strtg, perturb, nel, group_param, igtyp, defaults, glob_therm)
Definition inirig_mat.F:69
subroutine initia(iparg, elbuf, ms, in, v, x, ixs, ixq, ixc, ixt, ixp, ixr, detonators, geo, pm, rby, npby, lpby, npc, npts, pld, veul, ale_connectivity, skew, fill, ipart, itab, sensors, skvol, ixtg, thk, nloc_dmg, group_param_tab, glob_therm, igrnod, igrsurf, bufsf, vr, bufmat, xlas, las, dtelem, mss, msq, msc, mst, msp, msr, mstg, ptg, inc, nod2eltg, knod2eltg, inp, inr, intg, index, itri, kxx, ixx, xelemwa, iwa, nod2elq, knod2elq, nod2els, knod2els, kxsp, ixsp, nod2sp, ispcond, icode, iskew, iskn, ispsym, xframe, isptag, spbuf, mssx, nsigi, npbyl, lpbyl, rbyl, msnf, mssf, nsigsh, igeo, ipm, nsigs, nsigsph, vns, vnsx, stc, stt, stp, str, sttg, stur, bns, bnsx, volnod, bvolnod, etnod, nshnod, stifint, fxbdep, fxbvit, fxbacc, fxbipm, fxbrpm, fxbelm, fxbsig, fxbmod, ins, ptshel, ptsh3n, ptsol, ptquad, wma, ptsph, fxbnod, mbufel, mdepl, fxani, numel, nsigrs, sh4tree, sh3tree, mcp, temp, imerge2, iadmerge2, slnrbm, nslnrbm, rmstifn, rmstifr, ms_layer, zi_layer, itag, itagel, mcpc, mcptg, xrefc, xreftg, xrefs, mssa, msrt, irbe2, lrbe2, inivol, kvol, nbsubmat, ixs10, ixs16, ixs20, totaddmas, ipmas, stifn, msz2, itagn, sitage, itage, ixr_kj, elbuf_tab, nom_opt, ptr_nopt_rbe2, ptr_nopt_adm, ptr_nopt_fun, sol2sph, irst, sh3trim, xfem_tab, kxig3d, ixig3d, msig3d, knot, nctrlmax, wige, stack, rnoise, drape, sh4ang, sh3ang, geo_stack, igeo_stack, stifintr, strc, strp, strr, strtg, perturb, itagnd, nativ_sms, iloadp, facload, ptspri, nsigbeam, ptbeam, nsigtruss, pttruss, multi_fvm, sigi, sigsh, sigsp, sigsph, sigrs, sigbeam, sigtruss, strsglob, straglob, orthoglob, isigsh, iyldini, ksigsh3, fail_ini, iusolyld, iuser, iddlevel, inimap1d, inimap2d, func2d, fvm_inivel, tagprt_sms, igrbric, igrquad, igrsh4n, igrsh3n, igrpart, totmas, knotlocpc, knotlocel, vnige, bnige, fxbglm, fxbcpm, fxbcps, fxblm, fxbfls, fxbdls, fxb_matrix, fxb_matrix_add, fxb_last_adress, ptr_nopt_fxb, r_skew, knod2el1d, nod2el1d, ebcs_tab, rby_iniaxis, alea, knod2elc, nod2elc, dr, slrbody, drapeg, ipari, intbuf_tab, interfaces, mat_param, npreload_a, preload_a, fail_fractal, fail_brokmann, defaults, ndamp_freq_range, dampr, ibeam_vector, rbeam_vector, ikine)
Definition initia.F:188
subroutine sdlen3(x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, deltax, voln)
Definition sdlen3.F:41
subroutine slen(x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, j, area, aream)
Definition slen.F:31
program starter
Definition starter.F:39