OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
cdk6fint3.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!|| cdk6fint3 ../engine/source/elements/sh3n/coquedk6/cdk6fint3.F
25!||--- called by ------------------------------------------------------
26!|| cdk6forc3 ../engine/source/elements/sh3n/coquedk6/cdk6forc3.F
27!||====================================================================
28 SUBROUTINE cdk6fint3(JFT,JLT,VOL,THK0,FOR,MOM,
29 1 PX2,PY2,PX3,PY3,PB1,PB2,PB3,
30 2 F11,F12,F13,F21,F22,F23,F31,F32,F33,
31 3 F14,F15,F16,F24,F25,F26,F34,F35,F36,
32 4 N4X,N4Y,N4Z,N5X,N5Y,N5Z,N6X,N6Y,N6Z,
33 5 NEL)
34C--------------------------------------------------------------------------------------------------
35C-----------------------------------------------
36C I M P L I C I T T Y P E S
37C-----------------------------------------------
38#include "implicit_f.inc"
39#include "mvsiz_p.inc"
40C-----------------------------------------------
41C D U M M Y A R G U M E N T S
42C-----------------------------------------------
43 INTEGER JFT,JLT,NEL
44 MY_REAL
45 . FOR(NEL,5),MOM(NEL,3),
46 . PX2(*), PY2(*), PX3(*), PY3(*),
47 . PB1(MVSIZ,3,3),PB2(MVSIZ,3,3),PB3(MVSIZ,3,6),
48 . f11(*), f12(*), f13(*),
49 . f21(*), f22(*), f23(*), f31(*), f32(*), f33(*),
50 . n4x(*), n4y(*), n4z(*),n5x(*), n5y(*), n5z(*),
51 . n6x(*), n6y(*), n6z(*),
52 . f14(*), f15(*), f16(*),f24(*), f25(*), f26(*),
53 . f34(*), f35(*), f36(*),thk0(*),vol(*)
54C-----------------------------------------------
55C L O C A L V A R I A B L E S
56C-----------------------------------------------
57 INTEGER I
58 MY_REAL
59 . C2,FXX,FYY,FXY,FZ4,FZ5,FZ6,MXX,MYY,MXY,FXN,FYN,FXYN
60C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
61 DO I=jft,jlt
62 f12(i) =vol(i)*(px2(i)*for(i,1)+py2(i)*for(i,3))
63 f22(i) =vol(i)*(py2(i)*for(i,2)+px2(i)*for(i,3))
64 f13(i) =vol(i)*(px3(i)*for(i,1)+py3(i)*for(i,3))
65 f23(i) =vol(i)*(py3(i)*for(i,2)+px3(i)*for(i,3))
66 f11(i)=-f12(i)-f13(i)
67 f21(i)=-f22(i)-f23(i)
68 ENDDO
69 DO i=jft,jlt
70 c2=vol(i)*thk0(i)
71 mxx=c2*mom(i,1)
72 myy=c2*mom(i,2)
73 mxy=c2*mom(i,3)
74 f31(i) = pb1(i,1,1)*mxx+pb1(i,2,1)*myy+pb1(i,3,1)*mxy
75 f32(i) = pb1(i,1,2)*mxx+pb1(i,2,2)*myy+pb1(i,3,2)*mxy
76 f33(i) = pb1(i,1,3)*mxx+pb1(i,2,3)*myy+pb1(i,3,3)*mxy
77C
78 fz4 = pb2(i,1,1)*mxx+pb2(i,2,1)*myy+pb2(i,3,1)*mxy
79 f14(i) = fz4*n4x(i)
80 f24(i) = fz4*n4y(i)
81 f34(i) = fz4*n4z(i)
82 fz5 = pb2(i,1,2)*mxx+pb2(i,2,2)*myy+pb2(i,3,2)*mxy
83 f15(i) = fz5*n5x(i)
84 f25(i) = fz5*n5y(i)
85 f35(i) = fz5*n5z(i)
86 fz6 = pb2(i,1,3)*mxx+pb2(i,2,3)*myy+pb2(i,3,3)*mxy
87 f16(i) = fz6*n6x(i)
88 f26(i) = fz6*n6y(i)
89 f36(i) = fz6*n6z(i)
90C
91 fxx = pb3(i,1,1)*mxx + pb3(i,2,1)*myy + pb3(i,3,1)*mxy
92 fyy = pb3(i,1,2)*mxx + pb3(i,2,2)*myy + pb3(i,3,2)*mxy
93 fxy = pb3(i,1,3)*mxx + pb3(i,2,3)*myy + pb3(i,3,3)*mxy
94 fxn = n5x(i)*fxx
95 fyn = n6x(i)*fyy
96 fxyn = n4x(i)*fxy
97 f11(i) = f11(i) + pb3(i,1,4)*fxn+pb3(i,2,4)*fyn+pb3(i,3,4)*fxyn
98 f12(i) = f12(i) + pb3(i,1,5)*fxn+pb3(i,2,5)*fyn+pb3(i,3,5)*fxyn
99 f13(i) = f13(i) + pb3(i,1,6)*fxn+pb3(i,2,6)*fyn+pb3(i,3,6)*fxyn
100 fxn = n5y(i)*fxx
101 fyn = n6y(i)*fyy
102 fxyn = n4y(i)*fxy
103 f21(i) = f21(i) + pb3(i,1,4)*fxn+pb3(i,2,4)*fyn+pb3(i,3,4)*fxyn
104 f22(i) = f22(i) + pb3(i,1,5)*fxn+pb3(i,2,5)*fyn+pb3(i,3,5)*fxyn
105 f23(i) = f23(i) + pb3(i,1,6)*fxn+pb3(i,2,6)*fyn+pb3(i,3,6)*fxyn
106 fxn = n5z(i)*fxx
107 fyn = n6z(i)*fyy
108 fxyn = n4z(i)*fxy
109 f31(i) = f31(i) + pb3(i,1,4)*fxn+pb3(i,2,4)*fyn+pb3(i,3,4)*fxyn
110 f32(i) = f32(i) + pb3(i,1,5)*fxn+pb3(i,2,5)*fyn+pb3(i,3,5)*fxyn
111 f33(i) = f33(i) + pb3(i,1,6)*fxn+pb3(i,2,6)*fyn+pb3(i,3,6)*fxyn
112 ENDDO
113C
114 RETURN
115 END
cdk6fint3
subroutine cdk6fint3(jft, jlt, vol, thk0, for, mom, px2, py2, px3, py3, pb1, pb2, pb3, f11, f12, f13, f21, f22, f23, f31, f32, f33, f14, f15, f16, f24, f25, f26, f34, f35, f36, n4x, n4y, n4z, n5x, n5y, n5z, n6x, n6y, n6z, nel)
Definition cdk6fint3.F:34