OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
c_pon.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!|| c_pon ../starter/source/restart/ddsplit/c_pon.F
26!||--- called by ------------------------------------------------------
27!|| ddsplit ../starter/source/restart/ddsplit/ddsplit.F
28!||--- uses -----------------------------------------------------
29!||====================================================================
30 SUBROUTINE c_pon(ADDCNE ,NUMNOD_L,NODGLOB ,LCNE_L,I2NSNT ,
31 2 ADDCNI2,LCNI2_L ,MONVOL ,
32 3 CEP ,PROC ,NNMV_L ,NNMVC_L,ADDCNE_PXFEM,
33 4 LCNEPXFEM_L,INODPXFEM,ADDCNE_CRKXFEM,LCNECRKXFEM_L,
34 5 INODCRKXFEM,ADDCNCND,LCNCND_L,IGRSURF)
35C-----------------------------------------------
36C M o d u l e s
37C-----------------------------------------------
38 USE groupdef_mod
39C-----------------------------------------------
40C I m p l i c i t T y p e s
41C-----------------------------------------------
42#include "implicit_f.inc"
43C-----------------------------------------------
44C C o m m o n B l o c k s
45C-----------------------------------------------
46#include "com01_c.inc"
47#include "com04_c.inc"
48#include "param_c.inc"
49#include "com_xfem1.inc"
50C-----------------------------------------------
51C D u m m y A r g u m e n t s
52C-----------------------------------------------
53 INTEGER NUMNOD_L, LCNE_L, LCNI2_L, I2NSNT, NNMV_L, PROC,NNMVC_L,
54 . ADDCNE(0:*), NODGLOB(*), ADDCNI2(0:*),
55 . MONVOL(*), CEP(*),
56 . ADDCNE_PXFEM(0:*),LCNEPXFEM_L,INODPXFEM(*),
57 . addcne_crkxfem(0:*),lcnecrkxfem_l,inodcrkxfem(*),
58 . addcncnd(0:*),lcncnd_l
59 TYPE (SURF_) , DIMENSION(NSURF) :: IGRSURF
60C-----------------------------------------------
61C L o c a l V a r i a b l e s
62C-----------------------------------------------
63 INTEGER N,I,K0,K1,K2,K6,IS,NN,IAD,J,ITY,OFFC,OFFTG,NN_L,N1,
64 . ITYP,NVENT,IV,IADHOL,KIBHOL,KIBJET,NNC,KAD,NAV,IPVENT
65C-----------------------------------------------
66C
67C taille Matrice connectivite locale
68C
69 lcne_l = 0
70 DO i = 1, numnod_l
71 n = nodglob(i)
72 lcne_l = lcne_l + addcne(n+1)-addcne(n)
73 ENDDO
74C
75C taille Matrice connectivite locale (interface type 2)
76C
77 lcni2_l = 0
78 IF(i2nsnt>0) THEN
79 DO i = 1, numnod_l
80 n = nodglob(i)
81 lcni2_l = lcni2_l + addcni2(n+1)-addcni2(n)
82 END DO
83 END IF
84C
85 lcnepxfem_l = 0
86 IF(iplyxfem > 0) THEN
87 DO i = 1, numnod_l
88 n = nodglob(i)
89 n1 = inodpxfem(n)
90 IF(n1 > 0) THEN
91 lcnepxfem_l = lcnepxfem_l
92 . + addcne_pxfem(n1+1)-addcne_pxfem(n1)
93 ENDIF
94 ENDDO
95 ENDIF
96C-----Itet2 S10
97 lcncnd_l = 0
98 IF(ns10e>0) THEN
99 DO i = 1, numnod_l
100 n = nodglob(i)
101 lcncnd_l = lcncnd_l + addcncnd(n+1)-addcncnd(n)
102 END DO
103 END IF
104C
105C crack xfem layered shell
106C
107 lcnecrkxfem_l = 0
108 IF(icrack3d > 0) THEN
109 DO i = 1, numnod_l
110 n = nodglob(i)
111 n1 = inodcrkxfem(n) ! Num noeud systeme xfem global
112 IF (n1 > 0) THEN
113 lcnecrkxfem_l = lcnecrkxfem_l
114 . + addcne_crkxfem(n1+1) - addcne_crkxfem(n1)
115 ENDIF
116 ENDDO
117 ENDIF
118C
119C Taille connectivite airbag
120C
121 offc = numels+numelq
122 offtg =numels+numelq+ numelc+numelt+numelp+numelr
123 k0 = 0
124 k1 = 1
125 k2 = 1 + nimv*nvolu
126 kibjet = k2 + licbag
127 kibhol = kibjet + libagjet
128 k6 = 0
129 nn_l = 0
130 DO n = 1, nvolu
131 ityp = monvol(k1+1)
132 nav = monvol(k1+2)
133 is = monvol(k1+3)
134 nvent = monvol(k1+10)
135 nn = igrsurf(is)%NSEG
136 iadhol= kibhol+monvol(k1+11)
137 DO j = 1, nn
138 ity = igrsurf(is)%ELTYP(j)
139 i = igrsurf(is)%ELEM(j)
140 IF (ity==3) THEN
141 IF(cep(offc+i)==proc) nn_l = nn_l + 1
142 ELSEIF (ity==7) THEN
143 IF(cep(offtg+i)==proc) nn_l = nn_l + 1
144 ELSE
145 ENDIF
146 ENDDO
147C
148 IF(ityp==3.OR.ityp==4.OR.ityp==5.OR.ityp==7.OR.ityp==9) THEN
149 DO iv = 1, nvent
150 ipvent = monvol(iadhol+nibhol*(iv-1)+2-1)
151 IF(ipvent/=0) THEN
152 nnc=igrsurf(ipvent)%NSEG
153 DO j = 1, nnc
154 ity = igrsurf(ipvent)%ELTYP(j)
155 i = igrsurf(ipvent)%ELEM(j)
156 IF (ity==3) THEN
157 IF(cep(i+offc)==proc) THEN
158 k0 = k0 + 1
159 END IF
160 ELSEIF (ity==7) THEN
161 IF(cep(i+offtg)==proc) THEN
162 k0 = k0 + 1
163 END IF
164 END IF
165 END DO
166 END IF
167 END DO
168 END IF
169 IF(ityp==4.OR.ityp==5.OR.ityp==7.OR.ityp==9)THEN
170 DO iv = 1, nav
171 ipvent = monvol(k2+nicbag*(iv-1)+2-1)
172 IF(ipvent/=0) THEN
173 nnc=igrsurf(ipvent)%NSEG
174 DO j = 1, nnc
175 ity = igrsurf(ipvent)%ELTYP(j)
176 i = igrsurf(ipvent)%ELEM(j)
177 IF (ity==3) THEN
178 IF(cep(i+offc)==proc) THEN
179 k0 = k0 + 1
180 END IF
181 ELSEIF (ity==7) THEN
182 IF(cep(i+offtg)==proc) THEN
183 k0 = k0 + 1
184 END IF
185 END IF
186 END DO
187 END IF
188 END DO
189 END IF
190 k1 = k1 + nimv
191 k2 = k2 + nicbag * nav
192 k6 = k6 + nn
193 ENDDO
194 nnmv_l = nn_l
195 nnmvc_l = k0
196C
197 RETURN
198 END
c_pon
subroutine c_pon(addcne, numnod_l, nodglob, lcne_l, i2nsnt, addcni2, lcni2_l, monvol, cep, proc, nnmv_l, nnmvc_l, addcne_pxfem, lcnepxfem_l, inodpxfem, addcne_crkxfem, lcnecrkxfem_l, inodcrkxfem, addcncnd, lcncnd_l, igrsurf)
Definition c_pon.F:35
groupdef_mod
Definition groupdef_mod.F:662