OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
fsigsini.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!|| fsigsini ../starter/source/constraints/fxbody/fsigsini.F
25!||--- called by ------------------------------------------------------
26!|| fxbsini ../starter/source/constraints/fxbody/fxbsini.F
27!||--- calls -----------------------------------------------------
28!|| schrep ../starter/source/constraints/fxbody/fsigsini.F
29!|| sdefoi ../starter/source/constraints/fxbody/fsigsini.F
30!|| sderii ../starter/source/constraints/fxbody/fsigsini.F
31!|| sm1inif ../starter/source/constraints/fxbody/fsigsini.F
32!|| srhoi ../starter/source/constraints/fxbody/fsigsini.F
33!||====================================================================
34 SUBROUTINE fsigsini(FXBELM, IPARG , X , PM, IXS,
35 . GEO , FXBMOD, FXBSIG, R , NELS)
36C-----------------------------------------------
37C I m p l i c i t T y p e s
38C-----------------------------------------------
39#include "implicit_f.inc"
40C-----------------------------------------------
41C G l o b a l P a r a m e t e r s
42C-----------------------------------------------
43#include "mvsiz_p.inc"
44C-----------------------------------------------
45C C o m m o n B l o c k s
46C-----------------------------------------------
47#include "param_c.inc"
48C-----------------------------------------------
49C D u m m y A r g u m e n t s
50C-----------------------------------------------
51 INTEGER FXBELM(*), IPARG(NPARG,*), IXS(NIXS,*), NELS
52 my_real
53 . fxbsig(*), x(3,*), pm(npropm,*), fxbmod(*),
54 . geo(npropg,*), r(3,*)
55C-----------------------------------------------
56C L o c a l V a r i a b l e s
57C-----------------------------------------------
58 INTEGER I, IG, OFFSET, LAST, NFT, NFS, NG, IEL,
59 . n1, n2, n3, n4, n5, n6, n7, n8
60 INTEGER MAT(MVSIZ)
61 my_real
62 . x1(mvsiz), y1(mvsiz), z1(mvsiz),
63 . x2(mvsiz), y2(mvsiz), z2(mvsiz),
64 . x3(mvsiz), y3(mvsiz), z3(mvsiz),
65 . x4(mvsiz), y4(mvsiz), z4(mvsiz),
66 . x5(mvsiz), y5(mvsiz), z5(mvsiz),
67 . x6(mvsiz), y6(mvsiz), z6(mvsiz),
68 . x7(mvsiz), y7(mvsiz), z7(mvsiz),
69 . x8(mvsiz), y8(mvsiz), z8(mvsiz),
70 . d11, d12, d13, d21, d22, d23, d31, d32, d33, d41, d42, d43,
71 . d51, d52, d53, d61, d62, d63, d71, d72, d73, d81, d82, d83,
72 . vx1(mvsiz), vy1(mvsiz), vz1(mvsiz),
73 . vx2(mvsiz), vy2(mvsiz), vz2(mvsiz),
74 . vx3(mvsiz), vy3(mvsiz), vz3(mvsiz),
75 . vx4(mvsiz), vy4(mvsiz), vz4(mvsiz),
76 . vx5(mvsiz), vy5(mvsiz), vz5(mvsiz),
77 . vx6(mvsiz), vy6(mvsiz), vz6(mvsiz),
78 . vx7(mvsiz), vy7(mvsiz), vz7(mvsiz),
79 . vx8(mvsiz), vy8(mvsiz), vz8(mvsiz),
80 . x1p(mvsiz), y1p(mvsiz), z1p(mvsiz),
81 . x2p(mvsiz), y2p(mvsiz), z2p(mvsiz),
82 . x3p(mvsiz), y3p(mvsiz), z3p(mvsiz),
83 . x4p(mvsiz), y4p(mvsiz), z4p(mvsiz),
84 . x5p(mvsiz), y5p(mvsiz), z5p(mvsiz),
85 . x6p(mvsiz), y6p(mvsiz), z6p(mvsiz),
86 . x7p(mvsiz), y7p(mvsiz), z7p(mvsiz),
87 . x8p(mvsiz), y8p(mvsiz), z8p(mvsiz),
88 . px1(mvsiz), px2(mvsiz), px3(mvsiz), px4(mvsiz),
89 . py1(mvsiz), py2(mvsiz), py3(mvsiz), py4(mvsiz),
90 . pz1(mvsiz), pz2(mvsiz), pz3(mvsiz), pz4(mvsiz),
91 . vol(mvsiz), volp(mvsiz), d1(mvsiz), d2(mvsiz), d3(mvsiz),
92 . d4(mvsiz), d5(mvsiz), d6(mvsiz), rhop(mvsiz), sig(6,mvsiz),
93 . eint(mvsiz), sigl(6,mvsiz)
94C
95 DO ig=1,nels,mvsiz
96 offset=ig-1
97 last=min(mvsiz,nels-offset)
98 nft=offset*13
99 nfs=offset*7
100 DO i=1,last
101 ng=fxbelm(nft+13*(i-1)+1)
102 iel=iparg(3,ng)+fxbelm(nft+13*(i-1)+2)
103 mat(i)=ixs(1,iel)
104 x1(i)=x(1,ixs(2,iel))
105 y1(i)=x(2,ixs(2,iel))
106 z1(i)=x(3,ixs(2,iel))
107 x2(i)=x(1,ixs(3,iel))
108 y2(i)=x(2,ixs(3,iel))
109 z2(i)=x(3,ixs(3,iel))
110 x3(i)=x(1,ixs(4,iel))
111 y3(i)=x(2,ixs(4,iel))
112 z3(i)=x(3,ixs(4,iel))
113 x4(i)=x(1,ixs(5,iel))
114 y4(i)=x(2,ixs(5,iel))
115 z4(i)=x(3,ixs(5,iel))
116 x5(i)=x(1,ixs(6,iel))
117 y5(i)=x(2,ixs(6,iel))
118 z5(i)=x(3,ixs(6,iel))
119 x6(i)=x(1,ixs(7,iel))
120 y6(i)=x(2,ixs(7,iel))
121 z6(i)=x(3,ixs(7,iel))
122 x7(i)=x(1,ixs(8,iel))
123 y7(i)=x(2,ixs(8,iel))
124 z7(i)=x(3,ixs(8,iel))
125 x8(i)=x(1,ixs(9,iel))
126 y8(i)=x(2,ixs(9,iel))
127 z8(i)=x(3,ixs(9,iel))
128 n1=fxbelm(nft+13*(i-1)+3)
129 n2=fxbelm(nft+13*(i-1)+4)
130 n3=fxbelm(nft+13*(i-1)+5)
131 n4=fxbelm(nft+13*(i-1)+6)
132 n5=fxbelm(nft+13*(i-1)+7)
133 n6=fxbelm(nft+13*(i-1)+8)
134 n7=fxbelm(nft+13*(i-1)+9)
135 n8=fxbelm(nft+13*(i-1)+10)
136 d11=fxbmod(6*(n1-1)+1)
137 d12=fxbmod(6*(n1-1)+2)
138 d13=fxbmod(6*(n1-1)+3)
139 d21=fxbmod(6*(n2-1)+1)
140 d22=fxbmod(6*(n2-1)+2)
141 d23=fxbmod(6*(n2-1)+3)
142 d31=fxbmod(6*(n3-1)+1)
143 d32=fxbmod(6*(n3-1)+2)
144 d33=fxbmod(6*(n3-1)+3)
145 d41=fxbmod(6*(n4-1)+1)
146 d42=fxbmod(6*(n4-1)+2)
147 d43=fxbmod(6*(n4-1)+3)
148 d51=fxbmod(6*(n5-1)+1)
149 d52=fxbmod(6*(n5-1)+2)
150 d53=fxbmod(6*(n5-1)+3)
151 d61=fxbmod(6*(n6-1)+1)
152 d62=fxbmod(6*(n6-1)+2)
153 d63=fxbmod(6*(n6-1)+3)
154 d71=fxbmod(6*(n7-1)+1)
155 d72=fxbmod(6*(n7-1)+2)
156 d73=fxbmod(6*(n7-1)+3)
157 d81=fxbmod(6*(n8-1)+1)
158 d82=fxbmod(6*(n8-1)+2)
159 d83=fxbmod(6*(n8-1)+3)
160 vx1(i)=r(1,1)*d11+r(1,2)*d12+r(1,3)*d13
161 vy1(i)=r(2,1)*d11+r(2,2)*d12+r(2,3)*d13
162 vz1(i)=r(3,1)*d11+r(3,2)*d12+r(3,3)*d13
163 vx2(i)=r(1,1)*d21+r(1,2)*d22+r(1,3)*d23
164 vy2(i)=r(2,1)*d21+r(2,2)*d22+r(2,3)*d23
165 vz2(i)=r(3,1)*d21+r(3,2)*d22+r(3,3)*d23
166 vx3(i)=r(1,1)*d31+r(1,2)*d32+r(1,3)*d33
167 vy3(i)=r(2,1)*d31+r(2,2)*d32+r(2,3)*d33
168 vz3(i)=r(3,1)*d31+r(3,2)*d32+r(3,3)*d33
169 vx4(i)=r(1,1)*d41+r(1,2)*d42+r(1,3)*d43
170 vy4(i)=r(2,1)*d41+r(2,2)*d42+r(2,3)*d43
171 vz4(i)=r(3,1)*d41+r(3,2)*d42+r(3,3)*d43
172 vx5(i)=r(1,1)*d51+r(1,2)*d52+r(1,3)*d53
173 vy5(i)=r(2,1)*d51+r(2,2)*d52+r(2,3)*d53
174 vz5(i)=r(3,1)*d51+r(3,2)*d52+r(3,3)*d53
175 vx6(i)=r(1,1)*d61+r(1,2)*d62+r(1,3)*d63
176 vy6(i)=r(2,1)*d61+r(2,2)*d62+r(2,3)*d63
177 vz6(i)=r(3,1)*d61+r(3,2)*d62+r(3,3)*d63
178 vx7(i)=r(1,1)*d71+r(1,2)*d72+r(1,3)*d73
179 vy7(i)=r(2,1)*d71+r(2,2)*d72+r(2,3)*d73
180 vz7(i)=r(3,1)*d71+r(3,2)*d72+r(3,3)*d73
181 vx8(i)=r(1,1)*d81+r(1,2)*d82+r(1,3)*d83
182 vy8(i)=r(2,1)*d81+r(2,2)*d82+r(2,3)*d83
183 vz8(i)=r(3,1)*d81+r(3,2)*d82+r(3,3)*d83
184 x1p(i)=x1(i)+vx1(i)
185 y1p(i)=y1(i)+vy1(i)
186 z1p(i)=z1(i)+vz1(i)
187 x2p(i)=x2(i)+vx2(i)
188 y2p(i)=y2(i)+vy2(i)
189 z2p(i)=z2(i)+vz2(i)
190 x3p(i)=x3(i)+vx3(i)
191 y3p(i)=y3(i)+vy3(i)
192 z3p(i)=z3(i)+vz3(i)
193 x4p(i)=x4(i)+vx4(i)
194 y4p(i)=y4(i)+vy4(i)
195 z4p(i)=z4(i)+vz4(i)
196 x5p(i)=x5(i)+vx5(i)
197 y5p(i)=y5(i)+vy5(i)
198 z5p(i)=z5(i)+vz5(i)
199 x6p(i)=x6(i)+vx6(i)
200 y6p(i)=y6(i)+vy6(i)
201 z6p(i)=z6(i)+vz6(i)
202 x7p(i)=x7(i)+vx7(i)
203 y7p(i)=y7(i)+vy7(i)
204 z7p(i)=z7(i)+vz7(i)
205 x8p(i)=x8(i)+vx8(i)
206 y8p(i)=y8(i)+vy8(i)
207 z8p(i)=z8(i)+vz8(i)
208 ENDDO
209C
210 CALL sderii(last, volp,
211 . x1p, x2p, x3p, x4p, x5p, x6p, x7p, x8p,
212 . y1p, y2p, y3p, y4p, y5p, y6p, y7p, y8p,
213 . z1p, z2p, z3p, z4p, z5p, z6p, z7p, z8p,
214 . px1, px2, px3, px4,
215 . py1, py2, py3, py4,
216 . pz1, pz2, pz3, pz4)
217 CALL sderii(last, vol,
218 . x1, x2, x3, x4, x5, x6, x7, x8,
219 . y1, y2, y3, y4, y5, y6, y7, y8,
220 . z1, z2, z3, z4, z5, z6, z7, z8,
221 . px1, px2, px3, px4,
222 . py1, py2, py3, py4,
223 . pz1, pz2, pz3, pz4)
224 CALL sdefoi(last,
225 . px1, px2, px3, px4,
226 . py1, py2, py3, py4,
227 . pz1, pz2, pz3, pz4,
228 . vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8,
229 . vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8,
230 . vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8,
231 . d1, d2 , d3, d4, d5, d6 )
232 CALL srhoi(last, pm, vol, rhop, volp, mat)
233 CALL sm1inif(last, pm, sig, eint,
234 . rhop, d1, d2, d3, d4,
235 . d5, d6, mat)
236 CALL schrep(last, sig, sigl, r)
237C
238 DO i=1,last
239 iel=iparg(3,ng)+fxbelm(nft+13*(i-1)+2)
240 fxbsig(nfs+7*(i-1)+1)=sigl(1,i)
241 fxbsig(nfs+7*(i-1)+2)=sigl(2,i)
242 fxbsig(nfs+7*(i-1)+3)=sigl(3,i)
243 fxbsig(nfs+7*(i-1)+4)=sigl(4,i)
244 fxbsig(nfs+7*(i-1)+5)=sigl(5,i)
245 fxbsig(nfs+7*(i-1)+6)=sigl(6,i)
246 fxbsig(nfs+7*(i-1)+7)=eint(i)
247 ENDDO
248 ENDDO
249C
250 RETURN
251 END SUBROUTINE fsigsini
252!||====================================================================
253!|| sderii ../starter/source/constraints/fxbody/fsigsini.F
254!||--- called by ------------------------------------------------------
255!|| fsigsini ../starter/source/constraints/fxbody/fsigsini.F
256!||====================================================================
257 SUBROUTINE sderii(NEL, DET,
258 . X1, X2, X3, X4, X5, X6, X7, X8,
259 . Y1, Y2, Y3, Y4, Y5, Y6, Y7, Y8,
260 . Z1, Z2, Z3, Z4, Z5, Z6, Z7, Z8,
261 . PX1, PX2, PX3, PX4,
262 . PY1, PY2, PY3, PY4,
263 . PZ1, PZ2, PZ3, PZ4)
264C-----------------------------------------------
265C I m p l i c i t T y p e s
266C-----------------------------------------------
267#include "implicit_f.inc"
268C-----------------------------------------------
269C G l o b a l P a r a m e t e r s
270C-----------------------------------------------
271#include "mvsiz_p.inc"
272C-----------------------------------------------
273C D u m m y A r g u m e n t s
274C-----------------------------------------------
275 INTEGER :: NEL
276 my_real
277 . DET(*),
278 . X1(*), X2(*), X3(*), X4(*), X5(*), X6(*), X7(*), X8(*),
279 . Y1(*), Y2(*), Y3(*), Y4(*), Y5(*), Y6(*), Y7(*), Y8(*),
280 . Z1(*), Z2(*), Z3(*), Z4(*), Z5(*), Z6(*), Z7(*), Z8(*),
281 . PX1(*), PX2(*), PX3(*), PX4(*),
282 . py1(*), py2(*), py3(*), py4(*),
283 . pz1(*), pz2(*), pz3(*), pz4(*)
284C-----------------------------------------------
285C L o c a l V a r i a b l e s
286C-----------------------------------------------
287 INTEGER I
288 my_real
289 . JAC1(MVSIZ),JAC2(MVSIZ),JAC3(MVSIZ),
290 . JAC4(MVSIZ),JAC5(MVSIZ),JAC6(MVSIZ),
291 . DETT(MVSIZ) , JAC7(MVSIZ), JAC8(MVSIZ) , JAC9(MVSIZ),
292 . JACI1(MVSIZ), JACI2(MVSIZ), JACI3(MVSIZ),
293 . JACI4(MVSIZ), JACI5(MVSIZ), JACI6(MVSIZ),
294 . jaci7(mvsiz), jaci8(mvsiz), jaci9(mvsiz),
295 . x17(mvsiz) , x28(mvsiz) , x35(mvsiz) , x46(mvsiz),
296 . y17(mvsiz) , y28(mvsiz) , y35(mvsiz) , y46(mvsiz),
297 . z17(mvsiz) , z28(mvsiz) , z35(mvsiz) , z46(mvsiz),
298 . jac_59_68(mvsiz), jac_67_49(mvsiz), jac_48_57(mvsiz),
299 . jaci12(mvsiz), jaci45(mvsiz), jaci78(mvsiz),
300 . x_17_46(mvsiz) , x_28_35(mvsiz) ,
301 . y_17_46(mvsiz) , y_28_35(mvsiz) ,
302 . z_17_46(mvsiz) , z_28_35(mvsiz)
303C
304 DO i=1,nel
305 x17(i)=x7(i)-x1(i)
306 x28(i)=x8(i)-x2(i)
307 x35(i)=x5(i)-x3(i)
308 x46(i)=x6(i)-x4(i)
309 y17(i)=y7(i)-y1(i)
310 y28(i)=y8(i)-y2(i)
311 y35(i)=y5(i)-y3(i)
312 y46(i)=y6(i)-y4(i)
313 z17(i)=z7(i)-z1(i)
314 z28(i)=z8(i)-z2(i)
315 z35(i)=z5(i)-z3(i)
316 z46(i)=z6(i)-z4(i)
317 ENDDO
318C
319 DO i=1,nel
320 jac1(i)=x17(i)+x28(i)-x35(i)-x46(i)
321 jac2(i)=y17(i)+y28(i)-y35(i)-y46(i)
322 jac3(i)=z17(i)+z28(i)-z35(i)-z46(i)
323 x_17_46(i)=x17(i)+x46(i)
324 x_28_35(i)=x28(i)+x35(i)
325 y_17_46(i)=y17(i)+y46(i)
326 y_28_35(i)=y28(i)+y35(i)
327 z_17_46(i)=z17(i)+z46(i)
328 z_28_35(i)=z28(i)+z35(i)
329 ENDDO
330C
331C Jacobian matrix
332 DO i=1,nel
333 jac4(i)=x_17_46(i)+x_28_35(i)
334 jac5(i)=y_17_46(i)+y_28_35(i)
335 jac6(i)=z_17_46(i)+z_28_35(i)
336 jac7(i)=x_17_46(i)-x_28_35(i)
337 jac8(i)=y_17_46(i)-y_28_35(i)
338 jac9(i)=z_17_46(i)-z_28_35(i)
339 ENDDO
340C
341 DO i=1,nel
342 jac_59_68(i)=jac5(i)*jac9(i)-jac6(i)*jac8(i)
343 jac_67_49(i)=jac6(i)*jac7(i)-jac4(i)*jac9(i)
344 jac_48_57(i)=jac4(i)*jac8(i)-jac5(i)*jac7(i)
345 ENDDO
346C
347 DO i=1,nel
348 det(i)=one_over_64*(jac1(i)*jac_59_68(i)+jac2(i)*jac_67_49(i)
349 . +jac3(i)*jac_48_57(i))
350 ENDDO
351C
352 DO i=1,nel
353 dett(i)=one_over_64/det(i)
354 ENDDO
355C
356C Jacobian matrix inverse
357 DO i=1,nel
358 jaci1(i)=dett(i)*jac_59_68(i)
359 jaci4(i)=dett(i)*jac_67_49(i)
360 jaci7(i)=dett(i)*jac_48_57(i)
361 jaci2(i)=dett(i)*(-jac2(i)*jac9(i)+jac3(i)*jac8(i))
362 jaci5(i)=dett(i)*( jac1(i)*jac9(i)-jac3(i)*jac7(i))
363 jaci8(i)=dett(i)*(-jac1(i)*jac8(i)+jac2(i)*jac7(i))
364 jaci3(i)=dett(i)*( jac2(i)*jac6(i)-jac3(i)*jac5(i))
365 jaci6(i)=dett(i)*(-jac1(i)*jac6(i)+jac3(i)*jac4(i))
366 jaci9(i)=dett(i)*( jac1(i)*jac5(i)-jac2(i)*jac4(i))
367 ENDDO
368C
369 DO i=1,nel
370 jaci12(i)=jaci1(i)-jaci2(i)
371 jaci45(i)=jaci4(i)-jaci5(i)
372 jaci78(i)=jaci7(i)-jaci8(i)
373 ENDDO
374 DO i=1,nel
375 px3(i)= jaci12(i)+jaci3(i)
376 py3(i)= jaci45(i)+jaci6(i)
377 pz3(i)= jaci78(i)+jaci9(i)
378 px4(i)= jaci12(i)-jaci3(i)
379 py4(i)= jaci45(i)-jaci6(i)
380 pz4(i)= jaci78(i)-jaci9(i)
381 ENDDO
382 DO i=1,nel
383 jaci12(i)=jaci1(i)+jaci2(i)
384 jaci45(i)=jaci4(i)+jaci5(i)
385 jaci78(i)=jaci7(i)+jaci8(i)
386 ENDDO
387 DO i=1,nel
388 px1(i)=-jaci12(i)-jaci3(i)
389 py1(i)=-jaci45(i)-jaci6(i)
390 pz1(i)=-jaci78(i)-jaci9(i)
391 px2(i)=-jaci12(i)+jaci3(i)
392 py2(i)=-jaci45(i)+jaci6(i)
393 pz2(i)=-jaci78(i)+jaci9(i)
394 ENDDO
395C
396 RETURN
397 END SUBROUTINE sderii
398!||====================================================================
399!|| sdefoi ../starter/source/constraints/fxbody/fsigsini.F
400!||--- called by ------------------------------------------------------
401!|| fsigsini ../starter/source/constraints/fxbody/fsigsini.F
402!||====================================================================
403 SUBROUTINE sdefoi(NEL,
404 . PX1, PX2, PX3, PX4,
405 . PY1, PY2, PY3, PY4,
406 . PZ1, PZ2, PZ3, PZ4,
407 . VX1, VX2, VX3, VX4, VX5, VX6, VX7, VX8,
408 . VY1, VY2, VY3, VY4, VY5, VY6, VY7, VY8,
409 . VZ1, VZ2, VZ3, VZ4, VZ5, VZ6, VZ7, VZ8,
410 . DXX, DYY, DZZ, D4 , D5 , D6 )
411C-----------------------------------------------
412C I m p l i c i t T y p e s
413C-----------------------------------------------
414#include "implicit_f.inc"
415C-----------------------------------------------
416C G l o b a l P a r a m e t e r s
417C-----------------------------------------------
418#include "mvsiz_p.inc"
419C-----------------------------------------------
420C D u m m y A r g u m e n t s
421C-----------------------------------------------
422 INTEGER :: NEL
423 my_real
424 . VX1(*), VX2(*), VX3(*), VX4(*), VX5(*), VX6(*), VX7(*), VX8(*),
425 . VY1(*), VY2(*), VY3(*), VY4(*), VY5(*), VY6(*), VY7(*), VY8(*),
426 . VZ1(*), VZ2(*), VZ3(*), VZ4(*), VZ5(*), VZ6(*), VZ7(*), VZ8(*),
427 . PX1(*), PX2(*), PX3(*), PX4(*),
428 . PY1(*), PY2(*), PY3(*), PY4(*),
429 . PZ1(*), PZ2(*), PZ3(*), PZ4(*),
430 . DXX(*), DYY(*), DZZ(*), D4(*) , D5(*), D6(*)
431C-----------------------------------------------
432C L o c a l V a r i a b l e s
433C-----------------------------------------------
434 INTEGER I
435 my_real
436 . vx17(mvsiz), vy17(mvsiz), vz17(mvsiz),
437 . vx28(mvsiz), vy28(mvsiz), vz28(mvsiz),
438 . vx35(mvsiz), vy35(mvsiz), vz35(mvsiz),
439 . vx46(mvsiz), vy46(mvsiz), vz46(mvsiz)
440 my_real
441 . dxy(mvsiz), dxz(mvsiz), dyx(mvsiz), dyz(mvsiz), dzx(mvsiz),
442 . dzy(mvsiz)
443C
444 DO i=1,nel
445 vx17(i)=vx1(i)-vx7(i)
446 vx28(i)=vx2(i)-vx8(i)
447 vx35(i)=vx3(i)-vx5(i)
448 vx46(i)=vx4(i)-vx6(i)
449 vy17(i)=vy1(i)-vy7(i)
450 vy28(i)=vy2(i)-vy8(i)
451 vy35(i)=vy3(i)-vy5(i)
452 vy46(i)=vy4(i)-vy6(i)
453 vz17(i)=vz1(i)-vz7(i)
454 vz28(i)=vz2(i)-vz8(i)
455 vz35(i)=vz3(i)-vz5(i)
456 vz46(i)=vz4(i)-vz6(i)
457 ENDDO
458C
459 DO i=1,nel
460 dxx(i)=px1(i)*vx17(i)+px2(i)*vx28(i)+
461 . px3(i)*vx35(i)+px4(i)*vx46(i)
462 dyy(i)=py1(i)*vy17(i)+py2(i)*vy28(i)+
463 . py3(i)*vy35(i)+py4(i)*vy46(i)
464 dzz(i)=pz1(i)*vz17(i)+pz2(i)*vz28(i)+
465 . pz3(i)*vz35(i)+pz4(i)*vz46(i)
466 dxy(i)=py1(i)*vx17(i)+py2(i)*vx28(i)+
467 . py3(i)*vx35(i)+py4(i)*vx46(i)
468 dxz(i)=pz1(i)*vx17(i)+pz2(i)*vx28(i)+
469 . pz3(i)*vx35(i)+pz4(i)*vx46(i)
470 dyx(i)=px1(i)*vy17(i)+px2(i)*vy28(i)+
471 . px3(i)*vy35(i)+px4(i)*vy46(i)
472 dyz(i)=pz1(i)*vy17(i)+pz2(i)*vy28(i)+
473 . pz3(i)*vy35(i)+pz4(i)*vy46(i)
474 dzx(i)=px1(i)*vz17(i)+px2(i)*vz28(i)+
475 . px3(i)*vz35(i)+px4(i)*vz46(i)
476 dzy(i)=py1(i)*vz17(i)+py2(i)*vz28(i)+
477 . py3(i)*vz35(i)+py4(i)*vz46(i)
478 ENDDO
479C
480 DO i=1,nel
481 d4(i) = dxy(i)+dyx(i)
482 d5(i) = dyz(i)+dzy(i)
483 d6(i) = dxz(i)+dzx(i)
484 dxx(i) = dxx(i)
485 dyy(i) = dyy(i)
486 dzz(i) = dzz(i)
487 ENDDO
488C
489 RETURN
490 END SUBROUTINE sdefoi
491!||====================================================================
492!|| srhoi ../starter/source/constraints/fxbody/fsigsini.F
493!||--- called by ------------------------------------------------------
494!|| fsigsini ../starter/source/constraints/fxbody/fsigsini.F
495!||====================================================================
496 SUBROUTINE srhoi(NEL, PM, VOLO, RHON, VOLN, MAT)
497C-----------------------------------------------
498C I m p l i c i t T y p e s
499C-----------------------------------------------
500#include "implicit_f.inc"
501C-----------------------------------------------
502C G l o b a l P a r a m e t e r s
503C-----------------------------------------------
504#include "mvsiz_p.inc"
505C-----------------------------------------------
506C C o m m o n B l o c k s
507C-----------------------------------------------
508#include "param_c.inc"
509C-----------------------------------------------
510C D u m m y A r g u m e n t s
511C-----------------------------------------------
512 INTEGER NEL, MAT(*)
513 my_real
514 . PM(NPROPM,*), VOLO(*), RHON(*), VOLN(*)
515C-----------------------------------------------
516C L o c a l V a r i a b l e s
517C-----------------------------------------------
518 INTEGER I
519 my_real
520 . RHO0(MVSIZ)
521C
522 DO I=1,nel
523 rho0(i)=pm(1,mat(i))
524 rhon(i)=rho0(i)*(volo(i)/voln(i))
525 ENDDO
526C
527 RETURN
528 END SUBROUTINE srhoi
529!||====================================================================
530!|| sm1inif ../starter/source/constraints/fxbody/fsigsini.F
531!||--- called by ------------------------------------------------------
532!|| fsigsini ../starter/source/constraints/fxbody/fsigsini.F
533!||====================================================================
534 SUBROUTINE sm1inif(NEL, PM , SIG , EINT,
535 . RHO, D1 , D2 , D3 , D4 ,
536 . D5 ,D6 , MAT)
537C-----------------------------------------------
538C I m p l i c i t T y p e s
539C-----------------------------------------------
540#include "implicit_f.inc"
541C-----------------------------------------------
542C G l o b a l P a r a m e t e r s
543C-----------------------------------------------
544#include "mvsiz_p.inc"
545C-----------------------------------------------
546C C o m m o n B l o c k s
547C-----------------------------------------------
548#include "param_c.inc"
549C-----------------------------------------------
550C D u m m y A r g u m e n t s
551C-----------------------------------------------
552 INTEGER :: NEL, MAT(*)
553 my_real
554 . PM(NPROPM,*), SIG(6,*), EINT(*), RHO(*),D1(*), D2(*),
555 . D3(*), D4(*), D5(*), D6(*)
556C-----------------------------------------------
557C L o c a l V a r i a b l e s
558C-----------------------------------------------
559 INTEGER I,MX
560 my_real
561 . RHO0(MVSIZ), G1(MVSIZ), G2(MVSIZ), C1(MVSIZ), AMU(MVSIZ),
562 . DAV, E1, E2, E3, E4, E5, E6, NU1(MVSIZ),
563 . CNU1(MVSIZ)
564C
565 DO I=1,nel
566 mx =mat(i)
567 rho0(i)=pm( 1,mx)
568 nu1(i) =pm(21,mx)
569 cnu1(i)=nu1(i)/(one-two*nu1(i))
570 g1(i) =pm(22,mx)
571 g2(i) =two*g1(i)
572 c1(i) =pm(32,mx)
573 amu(i) =rho(i)/rho0(i)-one
574 ENDDO
575C
576 DO i=1,nel
577 dav=cnu1(i)*(d1(i)+d2(i)+d3(i))
578 sig(1,i)=g2(i)*(d1(i)+dav)
579 sig(2,i)=g2(i)*(d2(i)+dav)
580 sig(3,i)=g2(i)*(d3(i)+dav)
581 sig(4,i)=g1(i)*d4(i)
582 sig(5,i)=g1(i)*d5(i)
583 sig(6,i)=g1(i)*d6(i)
584 ENDDO
585C
586 DO i=1,nel
587 e1=d1(i)*sig(1,i)
588 e2=d2(i)*sig(2,i)
589 e3=d3(i)*sig(3,i)
590 e4=d4(i)*sig(4,i)
591 e5=d5(i)*sig(5,i)
592 e6=d6(i)*sig(6,i)
593 eint(i)=e1+e2+e3+e4+e5+e6
594 ENDDO
595C
596 RETURN
597 END SUBROUTINE sm1inif
598!||====================================================================
599!|| schrep ../starter/source/constraints/fxbody/fsigsini.F
600!||--- called by ------------------------------------------------------
601!|| fsigsini ../starter/source/constraints/fxbody/fsigsini.F
602!|| modbufel ../starter/source/constraints/fxbody/modbufel.F
603!||====================================================================
604 SUBROUTINE schrep(NEL, SIG, SIGL, R)
605C-----------------------------------------------
606C I m p l i c i t T y p e s
607C-----------------------------------------------
608#include "implicit_f.inc"
609C-----------------------------------------------
610C D u m m y A r g u m e n t s
611C-----------------------------------------------
612 INTEGER :: NEL
613 my_real
614 . sig(6,*), sigl(6,*), r(3,*)
615C-----------------------------------------------
616C L o c a l V a r i a b l e s
617C-----------------------------------------------
618 INTEGER I
619C
620 DO i=1,nel
621 sigl(1,i)=
622 . r(1,1)*(r(1,1)*sig(1,i)+r(2,1)*sig(4,i)+r(3,1)*sig(6,i))
623 . +r(2,1)*(r(1,1)*sig(4,i)+r(2,1)*sig(2,i)+r(3,1)*sig(5,i))
624 . +r(3,1)*(r(1,1)*sig(6,i)+r(2,1)*sig(5,i)+r(3,1)*sig(3,i))
625 sigl(2,i)=
626 . r(1,2)*(r(1,2)*sig(1,i)+r(2,2)*sig(4,i)+r(3,2)*sig(6,i))
627 . +r(2,2)*(r(1,2)*sig(4,i)+r(2,2)*sig(2,i)+r(3,2)*sig(5,i))
628 . +r(3,2)*(r(1,2)*sig(6,i)+r(2,2)*sig(5,i)+r(3,2)*sig(3,i))
629 sigl(3,i)=
630 . r(1,3)*(r(1,3)*sig(1,i)+r(2,3)*sig(4,i)+r(3,3)*sig(6,i))
631 . +r(2,3)*(r(1,3)*sig(4,i)+r(2,3)*sig(2,i)+r(3,3)*sig(5,i))
632 . +r(3,3)*(r(1,3)*sig(6,i)+r(2,3)*sig(5,i)+r(3,3)*sig(3,i))
633 sigl(4,i)=
634 . r(1,1)*(r(1,2)*sig(1,i)+r(2,2)*sig(4,i)+r(3,2)*sig(6,i))
635 . +r(2,1)*(r(1,2)*sig(4,i)+r(2,2)*sig(2,i)+r(3,2)*sig(5,i))
636 . +r(3,1)*(r(1,2)*sig(6,i)+r(2,2)*sig(5,i)+r(3,2)*sig(3,i))
637 sigl(5,i)=
638 . r(1,2)*(r(1,3)*sig(1,i)+r(2,3)*sig(4,i)+r(3,3)*sig(6,i))
639 . +r(2,2)*(r(1,3)*sig(4,i)+r(2,3)*sig(2,i)+r(3,3)*sig(5,i))
640 . +r(3,2)*(r(1,3)*sig(6,i)+r(2,3)*sig(5,i)+r(3,3)*sig(3,i))
641 sigl(6,i)=
642 . r(1,1)*(r(1,3)*sig(1,i)+r(2,3)*sig(4,i)+r(3,3)*sig(6,i))
643 . +r(2,1)*(r(1,3)*sig(4,i)+r(2,3)*sig(2,i)+r(3,3)*sig(5,i))
644 . +r(3,1)*(r(1,3)*sig(6,i)+r(2,3)*sig(5,i)+r(3,3)*sig(3,i))
645 ENDDO
646C
647 RETURN
648 END SUBROUTINE schrep
649
my_real
#define my_real
Definition cppsort.cpp:32
srhoi
subroutine srhoi(nel, pm, volo, rhon, voln, mat)
Definition fsigsini.F:497
fsigsini
subroutine fsigsini(fxbelm, iparg, x, pm, ixs, geo, fxbmod, fxbsig, r, nels)
Definition fsigsini.F:36
sdefoi
subroutine sdefoi(nel, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, vx1, vx2, vx3, vx4, vx5, vx6, vx7, vx8, vy1, vy2, vy3, vy4, vy5, vy6, vy7, vy8, vz1, vz2, vz3, vz4, vz5, vz6, vz7, vz8, dxx, dyy, dzz, d4, d5, d6)
Definition fsigsini.F:411
sm1inif
subroutine sm1inif(nel, pm, sig, eint, rho, d1, d2, d3, d4, d5, d6, mat)
Definition fsigsini.F:537
sderii
subroutine sderii(nel, det, x1, x2, x3, x4, x5, x6, x7, x8, y1, y2, y3, y4, y5, y6, y7, y8, z1, z2, z3, z4, z5, z6, z7, z8, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4)
Definition fsigsini.F:264
schrep
subroutine schrep(nel, sig, sigl, r)
Definition fsigsini.F:605
min
#define min(a, b)
Definition macros.h:20