OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s6ccumg3.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!|| s6ccumg3 ../engine/source/elements/thickshell/solide6c/s6ccumg3.F
25!||--- called by ------------------------------------------------------
26!|| s6cke3 ../engine/source/elements/thickshell/solide6c/s6cke3.F
27!||--- calls -----------------------------------------------------
28!|| set_ri33 ../engine/source/elements/shell/coqueba/cbasumg3.f
29!|| setprojks6 ../engine/source/elements/thickshell/solide6c/setprojks6.F
30!|| systran3 ../engine/source/elements/solid/solide8z/systran3.F
31!||====================================================================
32 SUBROUTINE s6ccumg3(
33 1 R11, R21, R31, R12,
34 2 R22, R32, R13, R23,
35 3 R33, K11, K12, K13,
36 4 K14, K15, K16, K22,
37 5 K23, K24, K25, K26,
38 6 K33, K34, K35, K36,
39 7 K44, K45, K46, K55,
40 8 K56, K66, X1, X2,
41 9 X3, X4, X5, X6,
42 A Y1, Y2, Y3, Y4,
43 B Y5, Y6, Z1, Z2,
44 C Z3, Z4, Z5, Z6,
45 D NEL)
46C-----------------------------------------------
47C I m p l i c i t T y p e s
48C-----------------------------------------------
49#include "implicit_f.inc"
50C-----------------------------------------------
51C G l o b a l P a r a m e t e r s
52C-----------------------------------------------
53#include "mvsiz_p.inc"
54C-----------------------------------------------
55C D u m m y A r g u m e n t s
56C-----------------------------------------------
57 INTEGER, INTENT(IN) :: NEL
58C REAL
59 my_real
60 . R11(*), R12(*), R13(*),
61 . R21(*), R22(*), R23(*),
62 . R31(*), R32(*), R33(*),
63 . K11(9,*) ,K12(9,*) ,K13(9,*) ,K14(9,*) ,K15(9,*) ,
64 . K16(9,*) ,K22(9,*) ,K23(9,*) ,K24(9,*) ,K25(9,*) ,
65 . K26(9,*) ,K33(9,*) ,K34(9,*) ,K35(9,*) ,K36(9,*) ,
66 . K44(9,*) ,K45(9,*) ,K46(9,*) ,K55(9,*) ,K56(9,*) ,
67 . K66(9,*)
68 my_real
69 . X1(*), X2(*), X3(*), X4(*),X5(*), X6(*),
70 . y1(*), y2(*), y3(*), y4(*),y5(*), y6(*),
71 . z1(*), z2(*), z3(*), z4(*),z5(*), z6(*)
72C-----------------------------------------------
73C C o m m o n B l o c k s
74C-----------------------------------------------
75C-----------------------------------------------
76C L o c a l V a r i a b l e s
77C-----------------------------------------------
78 INTEGER I,ISYM,L,J
79C REAL
80 my_real
81 . q(3,3,mvsiz)
82 my_real
83 . dr(3,3,mvsiz),
84 . r1(3,3,mvsiz),r2(3,3,mvsiz),r3(3,3,mvsiz),r4(3,3,mvsiz),
85 . r5(3,3,mvsiz),r6(3,3,mvsiz),di(6),xx,yy,zz,xy,xz,yz,rtr(6),
86 . abc,xxyz2,zzxy2,yyxz2,deta
87C-----------------------------------------------
88C-----PROJECTION----
89 DO i=1,nel
90 xx = x1(i)*x1(i)+x2(i)*x2(i)+x3(i)*x3(i)
91 1 +x4(i)*x4(i)+x5(i)*x5(i)+x6(i)*x6(i)
92 yy = y1(i)*y1(i)+y2(i)*y2(i)+y3(i)*y3(i)
93 1 +y4(i)*y4(i)+y5(i)*y5(i)+y6(i)*y6(i)
94 xy = x1(i)*y1(i)+x2(i)*y2(i)+x3(i)*y3(i)
95 1 +x4(i)*y4(i)+x5(i)*y5(i)+x6(i)*y6(i)
96 xz = x1(i)*z1(i)+x2(i)*z2(i)+x3(i)*z3(i)
97 1 +x4(i)*z4(i)+x5(i)*z5(i)+x6(i)*z6(i)
98 yz = y1(i)*z1(i)+y2(i)*z2(i)+y3(i)*z3(i)
99 1 +y4(i)*z4(i)+y5(i)*z5(i)+y6(i)*z6(i)
100 zz = z1(i)*z1(i)+z2(i)*z2(i)+z3(i)*z3(i)
101 1 +z4(i)*z4(i)+z5(i)*z5(i)+z6(i)*z6(i)
102 rtr(1)= yy+zz
103 rtr(2)= xx+zz
104 rtr(3)= xx+yy
105 rtr(4)= -xy
106 rtr(5)= -xz
107 rtr(6)= -yz
108C
109 abc = rtr(1)*rtr(2)*rtr(3)
110 xxyz2 = rtr(1)*rtr(6)*rtr(6)
111 yyxz2 = rtr(2)*rtr(5)*rtr(5)
112 zzxy2 = rtr(3)*rtr(4)*rtr(4)
113 deta = abc + two*rtr(4)*rtr(5)*rtr(6)-xxyz2-yyxz2-zzxy2
114 IF (deta<em20) THEN
115 deta=one
116 ELSE
117 deta=one/deta
118 ENDIF
119 di(1) = (abc-xxyz2)*deta/rtr(1)
120 di(2) = (abc-yyxz2)*deta/rtr(2)
121 di(3) = (abc-zzxy2)*deta/rtr(3)
122 di(4) = (rtr(5)*rtr(6)-rtr(4)*rtr(3))*deta
123 di(5) = (rtr(6)*rtr(4)-rtr(5)*rtr(2))*deta
124 di(6) = (rtr(4)*rtr(5)-rtr(6)*rtr(1))*deta
125 dr(1,1,i)= di(1)
126 dr(2,2,i)= di(2)
127 dr(3,3,i)= di(3)
128 dr(1,2,i)= di(4)
129 dr(1,3,i)= di(5)
130 dr(2,3,i)= di(6)
131 dr(2,1,i)= dr(1,2,i)
132 dr(3,1,i)= dr(1,3,i)
133 dr(3,2,i)= dr(2,3,i)
134 ENDDO
135C--------ini-RJ---------
136 DO i=1,nel
137 DO l=1,3
138 DO j=1,3
139 r1(l,j,i)=zero
140 r2(l,j,i)=zero
141 r3(l,j,i)=zero
142 r4(l,j,i)=zero
143 r5(l,j,i)=zero
144 r6(l,j,i)=zero
145 ENDDO
146 ENDDO
147 END DO
148C-------------Q<-Q^t---------------
149 DO i=1,nel
150 q(1,1,i)=r11(i)
151 q(1,2,i)=r21(i)
152 q(1,3,i)=r31(i)
153 q(2,1,i)=r12(i)
154 q(2,2,i)=r22(i)
155 q(2,3,i)=r32(i)
156 q(3,1,i)=r13(i)
157 q(3,2,i)=r23(i)
158 q(3,3,i)=r33(i)
159 ENDDO
160C
161 CALL set_ri33(x1 ,y1, z1 ,r1 ,1,nel)
162 CALL set_ri33(x2 ,y2, z2 ,r2 ,1,nel)
163 CALL set_ri33(x3 ,y3, z3 ,r3 ,1,nel)
164 CALL set_ri33(x4 ,y4, z4 ,r4 ,1,nel)
165 CALL set_ri33(x5 ,y5, z5 ,r5 ,1,nel)
166 CALL set_ri33(x6 ,y6, z6 ,r6 ,1,nel)
167 CALL setprojks6(dr ,r1 ,r2 ,r3 ,r4 ,
168 1 r5 ,r6 ,
169 2 k11,k12,k13,k14,k15,k16,k22,k23,
170 3 k24,k25,k26,k33,k34,k35,k36,k44,
171 4 k45,k46,k55,k56,k66,
172 5 q ,1,nel )
173C
174 RETURN
175C----------no projection ---------
176 DO i=1,nel
177 q(1,1,i)=r11(i)
178 q(2,1,i)=r21(i)
179 q(3,1,i)=r31(i)
180 q(1,2,i)=r12(i)
181 q(2,2,i)=r22(i)
182 q(3,2,i)=r32(i)
183 q(1,3,i)=r13(i)
184 q(2,3,i)=r23(i)
185 q(3,3,i)=r33(i)
186 ENDDO
187 isym=1
188C------------QKQ^t
189 CALL systran3(1,nel,q,k11,isym)
190 CALL systran3(1,nel,q,k22,isym)
191 CALL systran3(1,nel,q,k33,isym)
192 CALL systran3(1,nel,q,k44,isym)
193 CALL systran3(1,nel,q,k55,isym)
194 CALL systran3(1,nel,q,k66,isym)
195 isym=0
196 CALL systran3(1,nel,q,k12,isym)
197 CALL systran3(1,nel,q,k13,isym)
198 CALL systran3(1,nel,q,k14,isym)
199 CALL systran3(1,nel,q,k15,isym)
200 CALL systran3(1,nel,q,k16,isym)
201 CALL systran3(1,nel,q,k23,isym)
202 CALL systran3(1,nel,q,k24,isym)
203 CALL systran3(1,nel,q,k25,isym)
204 CALL systran3(1,nel,q,k26,isym)
205 CALL systran3(1,nel,q,k34,isym)
206 CALL systran3(1,nel,q,k35,isym)
207 CALL systran3(1,nel,q,k36,isym)
208 CALL systran3(1,nel,q,k45,isym)
209 CALL systran3(1,nel,q,k46,isym)
210 CALL systran3(1,nel,q,k56,isym)
211C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
212 RETURN
213 END
cbasumg3
subroutine cbasumg3(jft, jlt, vqn, vq, nplat, iplat, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, m11, m12, m13, m14, m22, m23, m24, m33, m34, m44, mf11, mf12, mf13, mf14, mf22, mf23, mf24, mf33, mf34, mf44, fm12, fm13, fm14, fm23, fm24, fm34, ke11, ke22, ke33, ke44, ke12, ke13, ke14, ke23, ke24, ke34, vcore, idril, iorth)
Definition cbasumg3.F:43
set_ri33
subroutine set_ri33(xi, yi, zi, ri, jft, jlt)
Definition cbasumg3.F:903
s6ccumg3
subroutine s6ccumg3(r11, r21, r31, r12, r22, r32, r13, r23, r33, k11, k12, k13, k14, k15, k16, k22, k23, k24, k25, k26, k33, k34, k35, k36, k44, k45, k46, k55, k56, k66, x1, x2, x3, x4, x5, x6, y1, y2, y3, y4, y5, y6, z1, z2, z3, z4, z5, z6, nel)
Definition s6ccumg3.F:46
setprojks6
subroutine setprojks6(dr, r1, r2, r3, r4, r5, r6, k11, k12, k13, k14, k15, k16, k22, k23, k24, k25, k26, k33, k34, k35, k36, k44, k45, k46, k55, k56, k66, vq, jft, jlt)
Definition setprojks6.F:39
systran3
subroutine systran3(jft, jlt, vq, kk, isym)
Definition systran3.F:30