OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
c3evec3.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!|| c3evec3 ../starter/source/elements/sh3n/coque3n/c3evec3.F
25!||--- called by ------------------------------------------------------
26!|| c3epsini ../starter/source/elements/sh3n/coque3n/c3epsini.f
27!|| c3fint_reg_ini ../starter/source/elements/sh3n/coque3n/c3fint_reg_ini.F
28!|| c3init3 ../starter/source/elements/sh3n/coque3n/c3init3.F
29!|| cdkinit3 ../starter/source/elements/sh3n/coquedk/cdkinit3.F
30!|| fsigtini ../starter/source/constraints/fxbody/fsigtini.F
31!|| inirig_mat ../starter/source/elements/initia/inirig_mat.f
32!|| inivoid ../starter/source/elements/initia/inivoid.f
33!||====================================================================
34 SUBROUTINE c3evec3(JFT ,JLT ,AREA,
35 . X1 ,X2 ,X3 ,Y1 ,Y2 ,Y3 ,
36 . Z1 ,Z2 ,Z3 ,E1X ,E2X ,E3X ,
37 . E1Y ,E2Y ,E3Y ,E1Z ,E2Z ,E3Z ,
38 . X31 ,Y31 ,Z31 ,X2L ,X3L ,Y3L )
39C-----------------------------------------------
40C I m p l i c i t T y p e s
41C-----------------------------------------------
42#include "implicit_f.inc"
43C-----------------------------------------------
44C G l o b a l P a r a m e t e r s
45C-----------------------------------------------
46#include "mvsiz_p.inc"
47C-----------------------------------------------
48C D u m m y A r g u m e n t s
49C-----------------------------------------------
50 INTEGER JFT, JLT
51 my_real X1(MVSIZ), X2(MVSIZ), X3(MVSIZ),
52 . Y1(MVSIZ), Y2(MVSIZ), Y3(MVSIZ),
53 . Z1(MVSIZ), Z2(MVSIZ), Z3(MVSIZ),
54 . e1x(mvsiz), e1y(mvsiz), e1z(mvsiz),
55 . e2x(mvsiz), e2y(mvsiz), e2z(mvsiz),
56 . e3x(mvsiz), e3y(mvsiz), e3z(mvsiz),
57 . x31(mvsiz), y31(mvsiz), z31(mvsiz),
58 . x32(mvsiz), y32(mvsiz), z32(mvsiz),
59 . x21(mvsiz), y21(mvsiz), z21(mvsiz), area(mvsiz),
60 . x2l(mvsiz), x3l(mvsiz), y3l(mvsiz)
61C-----------------------------------------------
62 INTEGER I
63 my_real SUM
64C-----------------------------------------------
65 DO I=jft,jlt
66 x21(i)=x2(i)-x1(i)
67 y21(i)=y2(i)-y1(i)
68 z21(i)=z2(i)-z1(i)
69 x31(i)=x3(i)-x1(i)
70 y31(i)=y3(i)-y1(i)
71 z31(i)=z3(i)-z1(i)
72 x32(i)=x3(i)-x2(i)
73 y32(i)=y3(i)-y2(i)
74 z32(i)=z3(i)-z2(i)
75 ENDDO
76C
77 DO i=jft,jlt
78 e1x(i)= x21(i)
79 e1y(i)= y21(i)
80 e1z(i)= z21(i)
81 x2l(i) = sqrt(e1x(i)*e1x(i)+e1y(i)*e1y(i)+e1z(i)*e1z(i))
82 e1x(i)=e1x(i)/x2l(i)
83 e1y(i)=e1y(i)/x2l(i)
84 e1z(i)=e1z(i)/x2l(i)
85 ENDDO
86C
87 DO i=jft,jlt
88 e3x(i)=y31(i)*z32(i)-z31(i)*y32(i)
89 e3y(i)=z31(i)*x32(i)-x31(i)*z32(i)
90 e3z(i)=x31(i)*y32(i)-y31(i)*x32(i)
91 sum = sqrt(e3x(i)*e3x(i)+e3y(i)*e3y(i)+e3z(i)*e3z(i))
92 e3x(i)=e3x(i)/sum
93 e3y(i)=e3y(i)/sum
94 e3z(i)=e3z(i)/sum
95 area(i) = half * sum
96 ENDDO
97C
98 DO i=jft,jlt
99 e2x(i)=e3y(i)*e1z(i)-e3z(i)*e1y(i)
100 e2y(i)=e3z(i)*e1x(i)-e3x(i)*e1z(i)
101 e2z(i)=e3x(i)*e1y(i)-e3y(i)*e1x(i)
102 sum = sqrt(e2x(i)*e2x(i)+e2y(i)*e2y(i)+e2z(i)*e2z(i))
103 e2x(i)=e2x(i)/sum
104 e2y(i)=e2y(i)/sum
105 e2z(i)=e2z(i)/sum
106 y3l(i)=e2x(i)*x31(i)+e2y(i)*y31(i)+e2z(i)*z31(i)
107 x3l(i)=e1x(i)*x31(i)+e1y(i)*y31(i)+e1z(i)*z31(i)
108 ENDDO
109C-----------
110 RETURN
111 END
subroutine c3epsini(elbuf_str, mat_param, jft, jlt, ismstr, mtn, ithk, pm, geo, ixtg, x, xreftg, for, thk, eint, gstr, nlay, px1g, py1g, py2g, x2s, x3s, y3s, uvar, ipm, nel, dir_a, dir_b, sigi, npf, tf, irep, ish3nfr, imat)
Definition c3epsini.F:45
subroutine area(d1, x, x2, y, y2, eint, stif0)
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 inivoid(elbuf_str, ixc, ixs, ixtg, 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, xrefc, xreftg, xrefs, mssa, volnod, bvolnod, vns, bns, sh3trim, isubstack, stack, rnoise, perturb, ele_area, part_area, iparttr, ixt, ipartp, ixp, mst, msp, stt, stp, strp, inp, stifint, mcpp, inr, msr, msrt, str, ipartr, itab, ixr, imerge2, iadmerge2, nel, defaults, glob_therm, ibeam_vector, rbeam_vector)
Definition inivoid.F:76
subroutine c3evec3(jft, jlt, area, x1, x2, x3, y1, y2, y3, z1, z2, z3, e1x, e2x, e3x, e1y, e2y, e3y, e1z, e2z, e3z, x31, y31, z31, x2l, x3l, y3l)
Definition c3evec3.F:39
program starter
Definition starter.F:39