OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i24iply_pxfem.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!|| i24iply_pxfem ../engine/source/interfaces/int24/i24iply_pxfem.F
25!||--- called by ------------------------------------------------------
26!|| i24mainf ../engine/source/interfaces/int24/i24main.F
27!||--- uses -----------------------------------------------------
28!|| plyxfem_mod ../engine/share/modules/plyxfem_mod.F
29!|| tri7box ../engine/share/modules/tri7box.F
30!||====================================================================
31 SUBROUTINE i24iply_pxfem(
32 1 JLT ,CAND_E ,MSEGTYP , IX1,IX2 ,
33 2 IX3 ,IX4 ,PENE,MS_PLY,INOD_PXFEM ,
34 3 IPLY ,ITAB)
35C-----------------------------------------------
36C M o d u l e s
37C-----------------------------------------------
38 USE tri7box
39 USE plyxfem_mod
40C-------------------------------------------f----
41C I m p l i c i t T y p e s
42C-----------------------------------------------
43#include "implicit_f.inc"
44C-----------------------------------------------
45C G l o b a l P a r a m e t e r s
46C-----------------------------------------------
47#include "mvsiz_p.inc"
48#include "comlock.inc"
49#include "param_c.inc"
50C-----------------------------------------------
51C D u m m y A r g u m e n t s
52C-----------------------------------------------
53 INTEGER JLT, CAND_E(*),IX1(*), IX2(*), IX3(*),
54 . IX4(*), INOD_PXFEM(*),
55 . IPLY(4,*),MSEGTYP(*),ITAB(*)
56 my_real
57 . pene(*),ms_ply(nplyxfe,*)
58C-----------------------------------------------
59C L o c a l V a r i a b l e s
60C-----------------------------------------------
61 INTEGER I,J1,J2,J3,J4,NN,ILY,ISEGTYP(MVSIZ),ITAG(MVSIZ),
62 . ISEG(MVSIZ),J5,J6,J7,J8,J9,J10,J11,J12
63 my_real
64 . xo16, xo17, xo15, xo14, x163, x153, x152, x142, x141,
65 . yo16, yo17, yo15, yo14, y163, y153, y152, y142, y141,
66 . zo16, zo17, zo15, zo14, z163, z153, z152, z142, z141,
67 . x51 , y51,z51,x52,y52,z52,x53,y53,z53,x54,y54,z54 ,
68 . x164,y164,z164,x174,y174,z174,x171,y171,z171,nx1,
69 . ny1,nz1,aaa,gapn,nx2,ny2,nz2,nx3,ny3,nz3,dx,dy,dz,
70 . nx4,ny4,nz4,
71 . sx125(mvsiz),sx235(mvsiz),sx345(mvsiz),sx415(mvsiz),
72 . sy125(mvsiz),sy235(mvsiz),sy345(mvsiz),sy415(mvsiz),
73 . sz125(mvsiz),sz235(mvsiz),sz345(mvsiz),sz415(mvsiz),
74 . sx1250(mvsiz),sx2350(mvsiz),sx3450(mvsiz),sx4150(mvsiz),
75 . sy1250(mvsiz),sy2350(mvsiz),sy3450(mvsiz),sy4150(mvsiz),
76 . sz1250(mvsiz),sz2350(mvsiz),sz3450(mvsiz),sz4150(mvsiz),
77 . sx2114(mvsiz),sx3215(mvsiz),sx4316(mvsiz),sx1417(mvsiz),
78 . sy2114(mvsiz),sy3215(mvsiz),sy4316(mvsiz),sy1417(mvsiz),
79 . sz2114(mvsiz),sz3215(mvsiz),sz4316(mvsiz),sz1417(mvsiz),
80 . sx21140(mvsiz),sx32150(mvsiz),sx43160(mvsiz),sx14170(mvsiz),
81 . sy21140(mvsiz),sy32150(mvsiz),sy43160(mvsiz),sy14170(mvsiz),
82 . sz21140(mvsiz),sz32150(mvsiz),sz43160(mvsiz),sz14170(mvsiz),
83 . nxx(mvsiz,17),nyy(mvsiz,17),nzz(mvsiz,17),bbb,
84 . iplyv(8,mvsiz)
85
86C ---------------------------------------------------------------------
87 DO i=1,jlt
88 isegtyp(i) = msegtyp(cand_e(i))
89 ENDDO
90C--------------------------------------------------------
91C
92C for computing contact ply forces.
93C
94 DO i=1,jlt
95C or by node segment if it needed
96 iply(1,i) = 0
97 iply(2,i) = 0
98 iply(3,i) = 0
99 iply(4,i) = 0
100
101 IF(pene(i) == zero) cycle
102 IF(isegtyp(i) == zero) cycle
103
104 IF(isegtyp(i) > 0) THEN
105 j1 = ix1(i)
106 IF(j1 > 0) THEN
107 nn = inod_pxfem(j1)
108 IF(nn > 0 ) THEN
109 ily = nplymax
110 DO WHILE(ms_ply(nn,ily) == zero)
111 ily = ily - 1
112 ENDDO
113 iply(1,i) = ily
114 ENDIF
115 ENDIF
116 j1 = ix2(i)
117 IF(j1 > 0) THEN
118 nn = inod_pxfem(j1)
119 IF(nn > 0) THEN
120 ily = nplymax
121 DO WHILE(ms_ply(nn,ily) == zero)
122 ily = ily - 1
123 ENDDO
124 iply(2,i) = ily
125 ENDIF
126 ENDIF
127 j1 = ix3(i)
128 IF(j1 > 0)THEN
129 nn = inod_pxfem(j1)
130 IF( nn > 0) THEN
131 ily = nplymax
132 DO WHILE(ms_ply(nn,ily) == zero)
133 ily = ily - 1
134 ENDDO
135 iply(3,i) = ily
136 ENDIF
137 ENDIF
138 j1 = ix4(i)
139 IF(j1 > 0) THEN
140 nn = inod_pxfem(j1)
141 IF( nn > 0) THEN
142 ily = nplymax
143 DO WHILE(ms_ply(nn,ily) == zero)
144 ily = ily - 1
145 ENDDO
146 iply(4,i) = ily
147 ENDIF
148 ENDIF
149 ELSE !
150 j1= ix1(i)
151 IF(j1 > 0) THEN
152 nn = inod_pxfem(j1)
153 IF( nn > 0) THEN
154 ily = 1
155 DO WHILE(ms_ply(nn,ily) == zero)
156 ily = ily + 1
157 ENDDO
158 iply(1,i) = ily
159 ENDIF
160 ENDIF
161 j1= ix2(i)
162 IF(j1 > 0)THEN
163 nn = inod_pxfem(j1)
164 IF(nn > 0) THEN
165 ily = 1
166 DO WHILE(ms_ply(nn,ily) == zero)
167 ily = ily + 1
168 ENDDO
169 iply(2,i) = ily
170 ENDIF
171 ENDIF
172C
173 j1= ix3(i)
174 IF(j1 > 0)THEN
175 nn = inod_pxfem(j1)
176 IF( nn > 0 ) THEN
177 ily = 1
178 DO WHILE(ms_ply(nn,ily) == zero)
179 ily = ily + 1
180 ENDDO
181 iply(3,i) = ily
182 ENDIF
183 ENDIF
184C
185 j1= ix4(i)
186 IF(j1> 0)THEN
187 nn = inod_pxfem(j1)
188 IF( nn > 0) THEN
189 ily = 1
190 DO WHILE(ms_ply(nn,ily) == zero)
191 ily = ily + 1
192 ENDDO
193 iply(4,i) = ily
194 ENDIF
195 ENDIF
196 ENDIF
197 ENDDO
198C
199 RETURN
200 END
my_real
#define my_real
Definition cppsort.cpp:32
i24iply_pxfem
subroutine i24iply_pxfem(jlt, cand_e, msegtyp, ix1, ix2, ix3, ix4, pene, ms_ply, inod_pxfem, iply, itab)
Definition i24iply_pxfem.F:35
plyxfem_mod
Definition plyxfem_mod.F:31