OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
coqini.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!|| coqini ../engine/source/elements/shell/coqini.F
25!||--- called by ------------------------------------------------------
26!|| radioss2 ../engine/source/engine/radioss2.F
27!||====================================================================
28 SUBROUTINE coqini
29C-----------------------------------------------
30C I m p l i c i t T y p e s
31C-----------------------------------------------
32#include "implicit_f.inc"
33C-----------------------------------------------
34C C o m m o n B l o c k s
35C-----------------------------------------------
36#include "com20_c.inc"
37C-----------------------------------------------
38C L o c a l V a r i a b l e s
39C-----------------------------------------------
40 INTEGER I, J
41C REAL
42 my_real
43 . z01(11,11), wf1(11,11), wm1(11,11), zn1(12,11)
44C-----------------------------------------------
45 DATA z01/
46 1 0. ,0. ,0. ,0. ,0. ,
47 1 0. ,0. ,0. ,0. ,0. ,0. ,
48 2 -.5 ,0.5 ,0. ,0. ,0. ,
49 2 0. ,0. ,0. ,0. ,0. ,0. ,
50 3 -.5 ,0. ,0.5 ,0. ,0. ,
51 3 0. ,0. ,0. ,0. ,0. ,0. ,
52 4 -.5 ,-.1666667,0.1666667,0.5 ,0. ,
53 4 0. ,0. ,0. ,0. ,0. ,0. ,
54 5 -.5 ,-.25 ,0. ,0.25 ,0.5 ,
55 5 0. ,0. ,0. ,0. ,0. ,0. ,
56 6 -.5 ,-.3 ,-.1 ,0.1 ,0.3 ,
57 6 0.5 ,0. ,0. ,0. ,0. ,0. ,
58 7 -.5 ,-.3333333,-.1666667,0.0 ,0.1666667,
59 7 0.3333333,0.5 ,0. ,0. ,0. ,0. ,
60 8 -.5 ,-.3571429,-.2142857,-.0714286,0.0714286,
61 8 0.2142857,0.3571429,0.5 ,0. ,0. ,0. ,
62 9 -.5 ,-.375 ,-.25 ,-.125 ,0.0 ,
63 9 0.125 ,0.25 ,0.375 ,0.5 ,0. ,0. ,
64 a -.5 ,-.3888889,-.2777778,-.1666667,-.0555555,
65 a 0.0555555,0.1666667,0.2777778,0.3888889,0.5 ,0. ,
66 b -.5 ,-.4 ,-.3 ,-.2 ,-.1 ,
67 b 0. ,0.1 ,0.2 ,0.3 ,0.4 ,0.5 /
68C
69 DATA wf1/
70 1 1. ,0. ,0. ,0. ,0. ,
71 1 0. ,0. ,0. ,0. ,0. ,0. ,
72 2 0.5 ,0.5 ,0. ,0. ,0. ,
73 2 0. ,0. ,0. ,0. ,0. ,0. ,
74 3 0.25 ,0.5 ,0.25 ,0. ,0. ,
75 3 0. ,0. ,0. ,0. ,0. ,0. ,
76 4 0.1666667,0.3333333,0.3333333,0.1666667,0. ,
77 4 0. ,0. ,0. ,0. ,0. ,0. ,
78 5 0.125 ,0.25 ,0.25 ,0.25 ,0.125 ,
79 5 0. ,0. ,0. ,0. ,0. ,0. ,
80 6 0.1 ,0.2 ,0.2 ,0.2 ,0.2 ,
81 6 0.1 ,0. ,0. ,0. ,0. ,0. ,
82 7 0.0833333,0.1666667,0.1666667,0.1666667,0.1666667,
83 7 0.1666667,0.0833333,0. ,0. ,0. ,0. ,
84 8 0.0714286,0.1428571,0.1428571,0.1428571,0.1428571,
85 8 0.1428571,0.1428571,0.0714286,0. ,0. ,0. ,
86 9 0.0625 ,0.125 ,0.125 ,0.125 ,0.125 ,
87 9 0.125 ,0.125 ,0.125 ,0.0625 ,0. ,0. ,
88 a 0.0555556,0.1111111,0.1111111,0.1111111,0.1111111,
89 a 0.1111111,0.1111111,0.1111111,0.1111111,0.0555556,0. ,
90 b 0.05 ,0.1 ,0.1 ,0.1 ,0.1 ,
91 b 0.1 ,0.1 ,0.1 ,0.1 ,0.1 ,0.05 /
92C
93C 022
94C
95C WM1i = Z01i / (n-1) if |Z01i| # 0.5
96C WM1i = 1/2(n-1) [1/4 - 1/3(n-1)] if Z01i = 0.5
97C WM1i = -1/2(n-1) [1/4 - 1/3(n-1)] if Z01i = -0.5
98C
99 DATA wm1/
100 1 0. ,0. ,0. ,0. ,0. ,
101 1 0. ,0. ,0. ,0. ,0. ,0. ,
102 2 -.0833333,0.0833333,0. ,0. ,0. ,
103 2 0. ,0. ,0. ,0. ,0. ,0. ,
104 3 -.0833333,0. ,0.0833333,0. ,0. ,
105 3 0. ,0. ,0. ,0. ,0. ,0. ,
106 4 -.0648148,-.0555556,0.0555556,0.0648148,0. ,
107 4 0. ,0. ,0. ,0. ,0. ,0. ,
108 5 -.0520833,-.0625 ,0. ,0.0625 ,0.0520833,
109 5 0. ,0. ,0. ,0. ,0. ,0. ,
110 6 -.0433333,-.0600000,-.0200000,0.0200000,0.0600000,
111 6 0.0433333,0.0000000,0.0000000,0.0000000,0.0000000,0.0000000,
112 7 -.0370370,-.0555556,-.0277778,0.0000000,0.0277778,
113 7 0.0555556,0.0370370,0.0000000,0.0000000,0.0000000,0.0000000,
114 8 -.0323129,-.0510204,-.0306122,-.0102041,0.0102041,
115 8 0.0306122,0.0510204,0.0323129,0.0000000,0.0000000,0.0000000,
116 9 -.0286458,-.0468750,-.0312500,-.0156250,0.0000000,
117 9 0.0156250,0.0312500,0.0468750,0.0286458,0.0000000,0.0000000,
118 a -.0257202,-.0432099,-.0308642,-.0185185,-.0061728,
119 a 0.0061728,0.0185185,0.0308642,0.0432099,0.0257202,0.0000000,
120 b -.0233333,-.0400000,-.0300000,-.0200000,-.0100000,
121 b 0.0000000,0.0100000,0.0200000,0.0300000,0.0400000,0.0233333/
122C
123 DATA zn1/
124 1 0. ,0. ,0. ,0. ,0. ,0. ,
125 1 0. ,0. ,0. ,0. ,0. ,0. ,
126 2 -.5 ,0.5 ,0. ,0. ,0. ,0. ,
127 2 0. ,0. ,0. ,0. ,0. ,0. ,
128 3 -.5 ,-.25 ,0.25 ,0.5 ,0. ,0. ,
129 3 0. ,0. ,0. ,0. ,0. ,0. ,
130 4 -.5 ,-.3333333,0. ,0.3333333,0.5 ,0. ,
131 4 0. ,0. ,0. ,0. ,0. ,0. ,
132 5 -.5 ,-.375 ,-0.125 ,0.125 ,0.375 ,0.5 ,
133 5 0. ,0. ,0. ,0. ,0. ,0. ,
134 6 -.5 ,-.4 ,-.2 ,0.0 ,0.2 ,0.4 ,
135 6 0.5 ,0. ,0. ,0. ,0. ,0. ,
136 7 -.5 ,-.4166667,-.25 ,-.0833333,0.0833333,0.25 ,
137 7 0.4166667,0.5 ,0. ,0. ,0. ,0. ,
138 8 -.5 ,-.4285715,-.2857143,-.1428572,0.0 ,0.1428572,
139 8 0.2857143,0.4285715,0.5 ,0. ,0. ,0. ,
140 9 -.5 ,-.4375 ,-.3125 ,-.1875 ,-.0625 ,0.0625 ,
141 9 0.1875 ,0.3125 ,0.4375 ,0.5 ,0. ,0. ,
142 a -.5 ,-.4444444,-.3333333,-.2222222,-.1111111,0. ,
143 a 0.1111111,0.2222222,0.3333333,0.4444444,0.5 ,0. ,
144 b -.5 ,-.45 ,-.35 ,-.25 ,-.15 ,-.05 ,
145 b 0.05 ,0.15 ,0.25 ,0.35 ,0.45 ,0.5 /
146
147 DO i=1,11
148 DO j=1,11
149 z0(i,j)=z01(i,j)
150 wf(i,j)=wf1(i,j)
151 wm(i,j)=wm1(i,j)
152 ENDDO
153 ENDDO
154 DO i = 1,12
155 DO j = 1, 11
156 zth(i,j)=zn1(i,j)
157 ENDDO
158 ENDDO
159 RETURN
160 END
161
162!||====================================================================
163!|| coqini_wm ../engine/source/elements/shell/coqini.F
164!||--- called by ------------------------------------------------------
165!|| mulawc ../engine/source/materials/mat_share/mulawc.F90
166!||====================================================================
167 SUBROUTINE coqini_wm(WM)
168C-----------------------------------------------
169C I m p l i c i t T y p e s
170C-----------------------------------------------
171#include "implicit_f.inc"
172C-----------------------------------------------
173C D u m m y A r g u m e n t s
174C-----------------------------------------------
175 my_real, intent(out):: wm(11,11)
176C-----------------------------------------------
177C L o c a l V a r i a b l e s
178C-----------------------------------------------
179 INTEGER I, J
180C REAL
181 my_real
182 . wm1(11,11)
183C-----------------------------------------------
184C 022
185C
186C WM1i = Z01i / (n-1) if |Z01i| # 0.5
187C WM1i = 1/2(n-1) [1/4 - 1/3(n-1)] if Z01i = 0.5
188C WM1i = -1/2(n-1) [1/4 - 1/3(n-1)] if Z01i = -0.5
189C
190 DATA wm1/
191 1 0. ,0. ,0. ,0. ,0. ,
192 1 0. ,0. ,0. ,0. ,0. ,0. ,
193 2 -.0833333,0.0833333,0. ,0. ,0. ,
194 2 0. ,0. ,0. ,0. ,0. ,0. ,
195 3 -.0833333,0. ,0.0833333,0. ,0. ,
196 3 0. ,0. ,0. ,0. ,0. ,0. ,
197 4 -.0648148,-.0555556,0.0555556,0.0648148,0. ,
198 4 0. ,0. ,0. ,0. ,0. ,0. ,
199 5 -.0520833,-.0625 ,0. ,0.0625 ,0.0520833,
200 5 0. ,0. ,0. ,0. ,0. ,0. ,
201 6 -.0433333,-.0600000,-.0200000,0.0200000,0.0600000,
202 6 0.0433333,0.0000000,0.0000000,0.0000000,0.0000000,0.0000000,
203 7 -.0370370,-.0555556,-.0277778,0.0000000,0.0277778,
204 7 0.0555556,0.0370370,0.0000000,0.0000000,0.0000000,0.0000000,
205 8 -.0323129,-.0510204,-.0306122,-.0102041,0.0102041,
206 8 0.0306122,0.0510204,0.0323129,0.0000000,0.0000000,0.0000000,
207 9 -.0286458,-.0468750,-.0312500,-.0156250,0.0000000,
208 9 0.0156250,0.0312500,0.0468750,0.0286458,0.0000000,0.0000000,
209 a -.0257202,-.0432099,-.0308642,-.0185185,-.0061728,
210 a 0.0061728,0.0185185,0.0308642,0.0432099,0.0257202,0.0000000,
211 b -.0233333,-.0400000,-.0300000,-.0200000,-.0100000,
212 b 0.0000000,0.0100000,0.0200000,0.0300000,0.0400000,0.0233333/
213
214 DO i=1,11
215 DO j=1,11
216 wm(i,j)=wm1(i,j)
217 ENDDO
218 ENDDO
219
220 RETURN
221 END
222
coqini_wm
subroutine coqini_wm(wm)
Definition coqini.F:168
coqini
subroutine coqini
Definition coqini.F:29
my_real
#define my_real
Definition cppsort.cpp:32