OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
c3coori.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!|| c3coori ../starter/source/elements/sh3n/coque3n/c3coori.F
25!||--- called by ------------------------------------------------------
26!|| c3init3 ../starter/source/elements/sh3n/coque3n/c3init3.f
27!|| cdkinit3 ../starter/source/elements/sh3n/coquedk/cdkinit3.F
28!|| inirig_mat ../starter/source/elements/initia/inirig_mat.F
29!|| inivoid ../starter/source/elements/initia/inivoid.F
30!|| lslocal ../starter/source/elements/xfem/lslocal.F
31!||--- calls -----------------------------------------------------
32!|| ancmsg ../starter/source/output/message/message.F
33!||--- uses -----------------------------------------------------
34!|| message_mod ../starter/share/message_module/message_mod.F
35!||====================================================================
36 SUBROUTINE c3coori(X ,XREFTG,IXP ,NGL ,
37 . X1 ,X2 ,X3 ,Y1 ,Y2 ,Y3 ,
38 . Z1 ,Z2 ,Z3 ,IX1 ,IX2 ,IX3 )
39C-----------------------------------------------
40 USE message_mod
41C-----------------------------------------------
42C I m p l i c i t T y p e s
43C-----------------------------------------------
44#include "implicit_f.inc"
45C-----------------------------------------------
46C G l o b a l P a r a m e t e r s
47C-----------------------------------------------
48#include "mvsiz_p.inc"
49C-----------------------------------------------
50C C o m m o n B l o c k s
51C-----------------------------------------------
52#include "scr03_c.inc"
53#include "vect01_c.inc"
54C-----------------------------------------------
55C V a r i a b l e s
56C-----------------------------------------------
57 INTEGER NGL(MVSIZ),IXP(NIXTG,*),IX1(MVSIZ),IX2(MVSIZ),IX3(MVSIZ)
58 my_real
59 . x(3,*),xreftg(3,3,*),
60 . x1(mvsiz),x2(mvsiz),x3(mvsiz),
61 . y1(mvsiz),y2(mvsiz),y3(mvsiz),
62 . z1(mvsiz),z2(mvsiz),z3(mvsiz)
63C-----------------------------------------------
64C L o c a l V a r i a b l e s
65C-----------------------------------------------
66 INTEGER I, TEST0, TEST1, TEST2
67 my_real E3X, E3Y, E3Z, SURF,
68 . xn13, xn23, xn12, x13, y13, z13, x12, y12, z12
69C=======================================================================
70C CONNECTIVITES ET MATERIEL
71C--------------------------------------------------
72 DO i=lft,llt
73 ngl(i)=ixp(6,i)
74 ix1(i)=ixp(2,i)
75 ix2(i)=ixp(3,i)
76 ix3(i)=ixp(4,i)
77 ENDDO
78C COORDONNEES
79C----------------------------
80 IF (nxref == 0) THEN
81 DO i=lft,llt
82 x1(i)=x(1,ix1(i))
83 y1(i)=x(2,ix1(i))
84 z1(i)=x(3,ix1(i))
85 x2(i)=x(1,ix2(i))
86 y2(i)=x(2,ix2(i))
87 z2(i)=x(3,ix2(i))
88 x3(i)=x(1,ix3(i))
89 y3(i)=x(2,ix3(i))
90 z3(i)=x(3,ix3(i))
91 ENDDO
92 ELSE
93 DO i=lft,llt
94 x1(i)=xreftg(1,1,i)
95 y1(i)=xreftg(1,2,i)
96 z1(i)=xreftg(1,3,i)
97 x2(i)=xreftg(2,1,i)
98 y2(i)=xreftg(2,2,i)
99 z2(i)=xreftg(2,3,i)
100 x3(i)=xreftg(3,1,i)
101 y3(i)=xreftg(3,2,i)
102 z3(i)=xreftg(3,3,i)
103 ENDDO
104 ENDIF
105c
106 DO i=lft,llt
107 xn13=(x1(i)-x3(i))**2+(y1(i)-y3(i))**2+(z1(i)-z3(i))**2
108 xn23=(x2(i)-x3(i))**2+(y2(i)-y3(i))**2+(z2(i)-z3(i))**2
109 xn12=(x1(i)-x2(i))**2+(y1(i)-y2(i))**2+(z1(i)-z2(i))**2
110 x13=x3(i)-x1(i)
111 y13=y3(i)-y1(i)
112 z13=z3(i)-z1(i)
113 x12=x2(i)-x1(i)
114 y12=y2(i)-y1(i)
115 z12=z2(i)-z1(i)
116 e3x=y12*z13-z12*y13
117 e3y=z12*x13-x12*z13
118 e3z=x12*y13-y12*x13
119 surf=half*sqrt(e3x*e3x+e3y*e3y+e3z*e3z)
120C
121 test0=0
122 test1=0
123 test2=0
124 IF(xn13+xn12<em20) THEN
125 CALL ancmsg(msgid=35,
126 . msgtype=msgerror,
127 . anmode=aninfo_blind_1,
128 . i1=ixp(nixtg,i))
129 test0=1
130 ELSE IF(surf<=em20) THEN
131 CALL ancmsg(msgid=35,
132 . msgtype=msgerror,
133 . anmode=aninfo_blind_1,
134 . i1=ixp(nixtg,i))
135 test1=1
136 ELSE IF(xn12<=em20.OR.xn23<=em20.OR.xn13<=em20) THEN
137 CALL ancmsg(msgid=36,
138 . msgtype=msgerror,
139 . anmode=aninfo_blind_1,
140 . i1=ixp(nixtg,i))
141 test2=1
142 ENDIF
143 IF(test0==1.OR.test1==1.OR.test2==1) THEN
144 x1(i)=zero
145 y1(i)=zero
146 z1(i)=zero
147 x2(i)=one
148 y2(i)=zero
149 z2(i)=zero
150 x3(i)=one
151 y3(i)=one
152 z3(i)=zero
153 ENDIF
154 ENDDO
155C-----------
156 RETURN
157 END
c3coori
subroutine c3coori(x, xreftg, ixp, ngl, x1, x2, x3, y1, y2, y3, z1, z2, z3, ix1, ix2, ix3)
Definition c3coori.F:39
c3init3
subroutine c3init3(elbuf_str, ixtg, pm, x, geo, igrsh4n, xmas, in, nvc, dtelem, igrsh3n, xreftg, offset, nel, ithk, thk, isigsh, sigsh, stifn, stifr, partsav, v, ipart, mstg, intg, ptg, skew, iparg, nsigsh, igeo, iuser, etnod, nshnod, sttg, ptsh3n, ipm, bufmat, sh3tree, mcp, mcptg, temp, cpt_eltens, part_area, itage, itagn, ixfem, npf, tf, sh3trim, xfem_str, isubstack, stack, rnoise, drape, sh3ang, iddlevel, geo_stack, igeo_stack, strtg, perturb, ish3n, iyldini, ele_area, nloc_dmg, ng, group_param, idrape, drapeg, mat_param, fail_fractal, fail_brokmann, glob_therm)
Definition c3init3.F:79
message_mod
Definition message_mod.F:1249
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889
starter
program starter
Definition starter.F:39