OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s10coork.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!|| s10coork ../engine/source/elements/solid/solide10/s10coork.F
25!||--- called by ------------------------------------------------------
26!|| s10ke3 ../engine/source/elements/solid/solide10/s10ke3.F
27!||--- uses -----------------------------------------------------
28!|| element_mod ../common_source/modules/elements/element_mod.F90
29!||====================================================================
30 SUBROUTINE s10coork(
31 1 X, IXS, IXS10, XX,
32 2 YY, ZZ, OFFG, OFF,
33 3 SAV, NC, NGL, MXT,
34 4 NGEO, K11, K12, K13,
35 5 K14, K15, K16, K17,
36 6 K18, K19, K10, K22,
37 7 K23, K24, K25, K26,
38 8 K27, K28, K29, K20,
39 9 K33, K34, K35, K36,
40 A K37, K38, K39, K30,
41 B K44, K45, K46, K47,
42 C K48, K49, K40, K55,
43 D K56, K57, K58, K59,
44 E K50, K66, K67, K68,
45 F K69, K60, K77, K78,
46 G K79, K70, K88, K89,
47 H K80, K99, K90, K00,
48 I NEL, ISMSTR)
49 use element_mod , only : nixs
50C-----------------------------------------------
51C I m p l i c i t T y p e s
52C-----------------------------------------------
53#include "implicit_f.inc"
54C-----------------------------------------------
55C G l o b a l P a r a m e t e r s
56C-----------------------------------------------
57#include "mvsiz_p.inc"
58C-----------------------------------------------
59C C o m m o n B l o c k s
60C-----------------------------------------------
61#include "scr18_c.inc"
62C-----------------------------------------------
63C D u m m y A r g u m e n t s
64C-----------------------------------------------
65 INTEGER, INTENT(IN) :: ISMSTR
66 INTEGER NEL
67C REAL
68 my_real
69 . x(3,*),
70 . xx(mvsiz,10), yy(mvsiz,10), zz(mvsiz,10),
71 . offg(*), off(*)
72 DOUBLE PRECISION
73 . SAV(NEL,30)
74 my_real
75 . K11(9,*) ,K12(9,*) ,K13(9,*) ,K14(9,*) ,K15(9,*) ,
76 . K16(9,*) ,K17(9,*) ,K18(9,*) ,K19(9,*) ,K10(9,*) ,
77 . K22(9,*) ,K23(9,*) ,K24(9,*) ,K25(9,*) ,K26(9,*) ,
78 . K27(9,*) ,K28(9,*) ,K29(9,*) ,K20(9,*) ,K33(9,*) ,
79 . K34(9,*) ,K35(9,*) ,K36(9,*) ,K37(9,*) ,K38(9,*) ,
80 . K39(9,*) ,K30(9,*) ,K44(9,*) ,K45(9,*) ,K46(9,*) ,
81 . K47(9,*) ,K48(9,*) ,K49(9,*) ,K40(9,*) ,K55(9,*) ,
82 . K56(9,*) ,K57(9,*) ,K58(9,*) ,K59(9,*) ,K50(9,*) ,
83 . k66(9,*) ,k67(9,*) ,k68(9,*) ,k69(9,*) ,k60(9,*) ,
84 . k77(9,*) ,k78(9,*) ,k79(9,*) ,k70(9,*) ,k88(9,*) ,
85 . k89(9,*) ,k80(9,*) ,k99(9,*) ,k90(9,*) ,k00(9,*)
86 INTEGER NC(MVSIZ,10),MXT(*), NGL(*),NGEO(*)
87 INTEGER IXS(NIXS,*), IXS10(6,*)
88C-----------------------------------------------
89C L o c a l V a r i a b l e s
90C-----------------------------------------------
91 INTEGER I, IPERM1(10),IPERM2(10),N,N1,N2,NN,IUN,J,MXT_1
92C REAL
93 DATA iperm1/0,0,0,0,1,2,3,1,2,3/
94 DATA iperm2/0,0,0,0,2,3,1,4,4,4/
95C-----------------------------------------------
96 iun=1
97C
98 mxt_1 = ixs(1,1)
99 DO i=1,nel
100 ngeo(i)=ixs(10,i)
101 ngl(i)=ixs(11,i)
102 mxt(i)=mxt_1
103 nc(i,1)=ixs(2,i)
104 nc(i,2)=ixs(4,i)
105 nc(i,3)=ixs(7,i)
106 nc(i,4)=ixs(6,i)
107 nc(i,5) =ixs10(1,i)
108 nc(i,6) =ixs10(2,i)
109 nc(i,7) =ixs10(3,i)
110 nc(i,8) =ixs10(4,i)
111 nc(i,9) =ixs10(5,i)
112 nc(i,10)=ixs10(6,i)
113 ENDDO
114C----------------------------
115C NODAL COORDINATES |
116C----------------------------
117 DO n=1,10
118 DO i=1,nel
119 nn = max(iun,nc(i,n))
120 xx(i,n)=x(1,nn)
121 yy(i,n)=x(2,nn)
122 zz(i,n)=x(3,nn)
123 off(i) = min(one,abs(offg(i)))
124 ENDDO
125 ENDDO
126C-----------
127 IF(ismstr==1.OR.(ismstr==2.AND.idtmin(1)==3))THEN
128 DO n=1,10
129 DO i=1,nel
130 IF(abs(offg(i))>one)THEN
131 xx(i,n)=sav(i,n)
132 yy(i,n)=sav(i,n+10)
133 zz(i,n)=sav(i,n+20)
134 off(i) = abs(offg(i))-one
135 ELSE
136 off(i) = offg(i)
137 ENDIF
138 ENDDO
139 ENDDO
140C
141 ENDIF
142C
143 DO n=5,10
144 n1=iperm1(n)
145 n2=iperm2(n)
146 DO i=1,nel
147 IF(nc(i,n)==0)THEN
148 xx(i,n) = half*(xx(i,n1)+xx(i,n2))
149 yy(i,n) = half*(yy(i,n1)+yy(i,n2))
150 zz(i,n) = half*(zz(i,n1)+zz(i,n2))
151 ENDIF
152 ENDDO
153 ENDDO
154C-----------
155 DO j=1,9
156 DO i=1,nel
157 k11(j,i)=zero
158 k12(j,i)=zero
159 k13(j,i)=zero
160 k14(j,i)=zero
161 k15(j,i)=zero
162 k16(j,i)=zero
163 k17(j,i)=zero
164 k18(j,i)=zero
165 k19(j,i)=zero
166 k10(j,i)=zero
167 k22(j,i)=zero
168 k23(j,i)=zero
169 k24(j,i)=zero
170 k25(j,i)=zero
171 k26(j,i)=zero
172 k27(j,i)=zero
173 k28(j,i)=zero
174 k29(j,i)=zero
175 k20(j,i)=zero
176 k33(j,i)=zero
177 k34(j,i)=zero
178 k35(j,i)=zero
179 k36(j,i)=zero
180 k37(j,i)=zero
181 k38(j,i)=zero
182 k39(j,i)=zero
183 k30(j,i)=zero
184 k44(j,i)=zero
185 k45(j,i)=zero
186 k46(j,i)=zero
187 k47(j,i)=zero
188 k48(j,i)=zero
189 k49(j,i)=zero
190 k40(j,i)=zero
191 k55(j,i)=zero
192 k56(j,i)=zero
193 k57(j,i)=zero
194 k58(j,i)=zero
195 k59(j,i)=zero
196 k50(j,i)=zero
197 k66(j,i)=zero
198 k67(j,i)=zero
199 k68(j,i)=zero
200 k69(j,i)=zero
201 k60(j,i)=zero
202 k77(j,i)=zero
203 k78(j,i)=zero
204 k79(j,i)=zero
205 k70(j,i)=zero
206 k88(j,i)=zero
207 k89(j,i)=zero
208 k80(j,i)=zero
209 k99(j,i)=zero
210 k90(j,i)=zero
211 k00(j,i)=zero
212 ENDDO
213 ENDDO
214C-----------
215 RETURN
216 END
my_real
#define my_real
Definition cppsort.cpp:32
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
s10coork
subroutine s10coork(x, ixs, ixs10, xx, yy, zz, offg, off, sav, nc, ngl, mxt, ngeo, k11, k12, k13, k14, k15, k16, k17, k18, k19, k10, k22, k23, k24, k25, k26, k27, k28, k29, k20, k33, k34, k35, k36, k37, k38, k39, k30, k44, k45, k46, k47, k48, k49, k40, k55, k56, k57, k58, k59, k50, k66, k67, k68, k69, k60, k77, k78, k79, k70, k88, k89, k80, k99, k90, k00, nel, ismstr)
Definition s10coork.F:49