OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s20coor3.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!|| s20coor3 ../starter/source/elements/solid/solide20/s20coor3.F
25!||--- called by ------------------------------------------------------
26!|| s20init3 ../starter/source/elements/solid/solide20/s20init3.F
27!||--- calls -----------------------------------------------------
28!|| checkvolume_8n ../starter/source/elements/solid/solide/checksvolume.F
29!||--- uses -----------------------------------------------------
30!|| message_mod ../starter/share/message_module/message_mod.F
31!||====================================================================
32 SUBROUTINE s20coor3(
33 1 X ,V ,IXS ,IXS20 ,XX ,
34 2 YY ,ZZ ,VX ,VY ,VZ ,
35 3 NC ,NGL ,MXT ,NGEO ,MASS,
36 4 DTELEM ,STI ,SIGG ,EINTG ,RHOG,
37 5 QG ,TEMP0 ,TEMP ,NEL ,NINTEMP)
38C-----------------------------------------------
39 USE message_mod
40 use element_mod , only : nixs
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 "vect01_c.inc"
53#include "com04_c.inc"
54C-----------------------------------------------
55C D u m m y A r g u m e n t s
56C-----------------------------------------------
57 INTEGER NEL
58 INTEGER ,INTENT(IN) :: NINTEMP
59 INTEGER NC(MVSIZ,20),MXT(*),NGL(*),NGEO(*),IXS(NIXS,*),IXS20(12,*)
60C REAL
61 my_real
62 . x(3,*), v(3,*), xx(mvsiz,20), yy(mvsiz,20), zz(mvsiz,20),
63 . vx(mvsiz,20), vy(mvsiz,20), vz(mvsiz,20),
64 . mass(mvsiz),dtelem(*),sti(*),
65 . sigg(nel,6),eintg(*),rhog(*),qg(*),temp(*), temp0(mvsiz)
66C-----------------------------------------------
67C L o c a l V a r i a b l e s
68C-----------------------------------------------
69 INTEGER I, IPERM1(20),IPERM2(20),N,N1,N2,NN,IUN
70C REAL
71 DATA IPERM1/0,0,0,0,0,0,0,0,1,2,3,4,1,2,3,4,5,6,7,8/
72 DATA IPERM2/0,0,0,0,0,0,0,0,2,3,4,1,5,6,7,8,6,7,8,5/
73C-----------------------------------------------
74C E x t e r n a l F u n c t i o n s
75C-----------------------------------------------
76 my_real
77 . checkvolume_8n
78C=======================================================================
79 iun=1
80C
81 DO i=lft,llt
82 ngeo(i)=ixs(10,i)
83 ngl(i) =ixs(11,i)
84 mxt(i) =ixs(1,i)
85 nc(i,1)=ixs(2,i)
86 nc(i,2)=ixs(3,i)
87 nc(i,3)=ixs(4,i)
88 nc(i,4)=ixs(5,i)
89 nc(i,5)=ixs(6,i)
90 nc(i,6)=ixs(7,i)
91 nc(i,7)=ixs(8,i)
92 nc(i,8)=ixs(9,i)
93 nc(i,9) =ixs20(1,i)
94 nc(i,10) =ixs20(2,i)
95 nc(i,11) =ixs20(3,i)
96 nc(i,12) =ixs20(4,i)
97 nc(i,13) =ixs20(5,i)
98 nc(i,14) =ixs20(6,i)
99 nc(i,15) =ixs20(7,i)
100 nc(i,16) =ixs20(8,i)
101 nc(i,17) =ixs20(9,i)
102 nc(i,18) =ixs20(10,i)
103 nc(i,19) =ixs20(11,i)
104 nc(i,20) =ixs20(12,i)
105 dtelem(i)=ep30
106 sti(i)=zero
107 eintg(i)=zero
108 rhog(i)=zero
109 qg(i)=zero
110 sigg(i,1)=zero
111 sigg(i,2)=zero
112 sigg(i,3)=zero
113 sigg(i,4)=zero
114 sigg(i,5)=zero
115 sigg(i,6)=zero
116 IF (checkvolume_8n(x ,ixs(1,i)) < zero) THEN
117C renumber connectivity
118 nc(i,1) =ixs(6,i)
119 nc(i,2) =ixs(7,i)
120 nc(i,3) =ixs(8,i)
121 nc(i,4) =ixs(9,i)
122 nc(i,5) =ixs(2,i)
123 nc(i,6) =ixs(3,i)
124 nc(i,7) =ixs(4,i)
125 nc(i,8) =ixs(5,i)
126 nc(i,9) =ixs20(9,i)
127 nc(i,10)=ixs20(10,i)
128 nc(i,11)=ixs20(11,i)
129 nc(i,12)=ixs20(12,i)
130 nc(i,17)=ixs20(1,i)
131 nc(i,18)=ixs20(2,i)
132 nc(i,19)=ixs20(3,i)
133 nc(i,20)=ixs20(4,i)
134 ixs(2,i) = nc(i,1)
135 ixs(3,i) = nc(i,2)
136 ixs(4,i) = nc(i,3)
137 ixs(5,i) = nc(i,4)
138 ixs(6,i) = nc(i,5)
139 ixs(7,i) = nc(i,6)
140 ixs(8,i) = nc(i,7)
141 ixs(9,i) = nc(i,8)
142 ixs20(1,i) = nc(i,9)
143 ixs20(2,i) = nc(i,10)
144 ixs20(3,i) = nc(i,11)
145 ixs20(4,i) = nc(i,12)
146 ixs20(5,i) = nc(i,13)
147 ixs20(6,i) = nc(i,14)
148 ixs20(7,i) = nc(i,15)
149 ixs20(8,i) = nc(i,16)
150 ixs20(9,i) = nc(i,17)
151 ixs20(10,i) = nc(i,18)
152 ixs20(11,i) = nc(i,19)
153 ixs20(12,i) = nc(i,20)
154 ENDIF
155 ENDDO
156C
157C----------------------------
158C NODAL COORDINATES |
159C----------------------------
160 DO n=1,20
161 DO i=lft,llt
162 nn = max(iun,nc(i,n))
163 xx(i,n)=x(1,nn)
164 yy(i,n)=x(2,nn)
165 zz(i,n)=x(3,nn)
166 vx(i,n)=v(1,nn)
167 vy(i,n)=v(2,nn)
168 vz(i,n)=v(3,nn)
169 ENDDO
170 ENDDO
171C
172 DO i=lft,llt
173 mass(i)=zero
174 ENDDO
175C
176 DO n=9,20
177 n1=iperm1(n)
178 n2=iperm2(n)
179 DO i=lft,llt
180 IF(nc(i,n)==0)THEN
181 xx(i,n) = half*(xx(i,n1)+xx(i,n2))
182 yy(i,n) = half*(yy(i,n1)+yy(i,n2))
183 zz(i,n) = half*(zz(i,n1)+zz(i,n2))
184 vx(i,n) = half*(vx(i,n1)+vx(i,n2))
185 vy(i,n) = half*(vy(i,n1)+vy(i,n2))
186 vz(i,n) = half*(vz(i,n1)+vz(i,n2))
187 ENDIF
188 ENDDO
189 ENDDO
190C
191 IF (jthe < 0 .or. nintemp > 0) THEN
192 IF(nintemp > 0 ) THEN
193 DO n =1,20
194 DO i=lft,llt
195 nn = max(iun,nc(i,n))
196 IF(temp(nn)== zero) temp(nn) = temp0(i)
197 ENDDO
198 ENDDO
199 ELSE
200 DO n =1,20
201 DO i=lft,llt
202 nn = max(iun,nc(i,n))
203 temp(nn) = temp0(i)
204 ENDDO
205 ENDDO
206 ENDIF
207 ENDIF
208C------------
209 RETURN
210 END
checkvolume_8n
function checkvolume_8n(x, ixs)
Definition checksvolume.F:39
max
#define max(a, b)
Definition macros.h:21
message_mod
Definition message_mod.F:1257
s20coor3
subroutine s20coor3(x, v, ixs, ixs20, xx, yy, zz, vx, vy, vz, nc, ngl, mxt, ngeo, mass, dtelem, sti, sigg, eintg, rhog, qg, temp0, temp, nel, nintemp)
Definition s20coor3.F:38