OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s6cfint3.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!|| s6cfint3 ../engine/source/elements/thickshell/solide6c/s6cfint3.F
25!||--- called by ------------------------------------------------------
26!|| s6cforc3 ../engine/source/elements/thickshell/solide6c/s6cforc3.F
27!||====================================================================
28 SUBROUTINE s6cfint3(
29 1 SIG, PX1, PX2, PX3,
30 2 PX4, PY1, PY2, PY3,
31 3 PY4, PZ1, PZ2, PZ3,
32 4 PZ4, PX1H, PX2H, PX3H,
33 5 PY1H, PY2H, PY3H, PZ1H,
34 6 PZ2H, PZ3H, JI33, B1X,
35 7 B1Y, B2Y, B2X, B1122,
36 8 B1221, B2212, B1121, B1XH,
37 9 B1YH, B2XH, B2YH, B1122H,
38 A B1221H, B2212H, B1121H, F11,
39 B F21, F31, F12, F22,
40 C F32, F13, F23, F33,
41 D F14, F24, F34, F15,
42 E F25, F35, F16, F26,
43 F F36, VOL, QVIS, EINT,
44 G RHO, Q, EPLA, EPSD,
45 H EPSDM, SIGM, EINTM, RHOM,
46 I QM, EPLASM, ZI, WI,
47 J VOLG, OFF, NU, VOL0,
48 K VOL0G, G_PLA, G_EPSD, NEL,
49 L SVIS, G_WPLA, L_WPLA, G_WPLA_FLAG)
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 D u m m y A r g u m e n t s
60C-----------------------------------------------
61 INTEGER G_PLA,G_EPSD,NEL
62C REAL
63 my_real
64 . SIG(NEL,6),
65 . PX1(*), PX2(*), PX3(*), PX4(*),
66 . PY1(*), PY2(*), PY3(*), PY4(*),
67 . PZ1(*), PZ2(*), PZ3(*), PZ4(*),
68 . PX1H(*), PX2H(*), PX3H(*),
69 . PY1H(*), PY2H(*), PY3H(*),
70 . PZ1H(*), PZ2H(*), PZ3H(*),
71 . JI33(*),B1X(MVSIZ,2),B1Y(MVSIZ,2),B2X(MVSIZ,2),B2Y(MVSIZ,2),
72 . B1XH(MVSIZ,2),B1YH(MVSIZ,2),B2XH(MVSIZ,2),B2YH(MVSIZ,2),
73 . B1122(*),B1221(*),B2212(*),B1121(*),
74 . B1122H(*),B1221H(*),B2212H(*),B1121H(*),
75 . F11(*),F21(*),F31(*),F12(*),F22(*),F32(*),
76 . F13(*),F23(*),F33(*),F14(*),F24(*),F34(*),
77 . F15(*),F25(*),F35(*),F16(*),F26(*),F36(*),
78 . VOL(*),QVIS(*),
79 . EINT(*),RHO(*),Q(*),EPLASM(*),EPLA(*),
80 . SIGM(NEL,6),EINTM(*),RHOM(*),QM(*),EPSD(*),EPSDM(*),
81 . ZI,WI,VOLG(*),OFF(*),NU(*),VOL0(*),VOL0G(*)
82 my_real, DIMENSION(MVSIZ,6), INTENT(INOUT) :: svis
83 INTEGER, INTENT(IN) :: G_WPLA_FLAG
84 my_real,DIMENSION(NEL*G_WPLA_FLAG), INTENT(INOUT) :: g_wpla
85 my_real,DIMENSION(NEL*G_WPLA_FLAG), INTENT(IN) :: l_wpla
86C-----------------------------------------------
87C L o c a l V a r i a b l e s
88C-----------------------------------------------
89 INTEGER I, J
90C REAL
91 my_real
92 . S1(MVSIZ), S2(MVSIZ), S3(MVSIZ),
93 . S4(MVSIZ), S5(MVSIZ), S6(MVSIZ),
94 . RX1(MVSIZ), RY1(MVSIZ), SX1(MVSIZ), SY1(MVSIZ),
95 . FINT,FINTX,FINTY,FINTZ,FAC(MVSIZ),FXC,FYC,FINSX,FINSY
96 my_real
97 . NU1,FINSZ
98C=======================================================================
99 DO I=1,nel
100 s1(i)=(sig(i,1)+svis(i,1)-qvis(i))*vol(i)
101 s2(i)=(sig(i,2)+svis(i,2)-qvis(i))*vol(i)
102 s3(i)=(sig(i,3)+svis(i,3)-qvis(i))*vol(i)
103 s4(i)=(sig(i,4)+svis(i,4))*vol(i)
104 s5(i)=(sig(i,5)+svis(i,5))*vol(i)
105 s6(i)=(sig(i,6)+svis(i,6))*vol(i)
106 ENDDO
107C 12
108C -------constant part---------
109 DO i=1,nel
110 fint=s1(i)*px1(i)+s4(i)*py1(i)
111 fxc = ji33(i)*s6(i)
112 fyc = ji33(i)*s5(i)
113 fintx=s1(i)*px4(i)+s4(i)*py4(i) + three*fxc
114 finsx=s6(i)*(b1x(i,1)-b1x(i,2))-s5(i)*(b2x(i,1)-b2x(i,2))-fxc
115 f11(i)=f11(i)-fint+fintx+finsx
116 f14(i)=f14(i)-fint-fintx-finsx
117 fint=s2(i)*py1(i)+s4(i)*px1(i)
118 finty=s2(i)*py4(i)+s4(i)*px4(i) + three*fyc
119 finsy=s6(i)*(b1y(i,1)-b1y(i,2))-s5(i)*(b2y(i,1)-b2y(i,2))-fyc
120 f21(i)=f21(i)-fint+finty+finsy
121 f24(i)=f24(i)-fint-finty-finsy
122 fint=s3(i)*pz1(i)+half*(s6(i)*px1(i)+s5(i)*py1(i))
123 fintz=s3(i)*pz4(i)
124 f31(i)=f31(i)-fint+fintz
125 f34(i)=f34(i)-fint-fintz
126C
127 fint=s1(i)*px2(i)+s4(i)*py2(i)
128 finsx=s6(i)*(b1221(i)+b1x(i,2))-s5(i)*(b1121(i)+b2x(i,2))
129 f12(i)=f12(i)-fint+fintx+finsx
130 f15(i)=f15(i)-fint-fintx-finsx
131 fint=s2(i)*py2(i)+s4(i)*px2(i)
132 finsy=s6(i)*(b2212(i)+b1y(i,2))-s5(i)*(b1122(i)+b2y(i,2))
133 f22(i)=f22(i)-fint+finty+finsy
134 f25(i)=f25(i)-fint-finty-finsy
135 fint=s3(i)*pz2(i)+(s6(i)*px2(i)+s5(i)*py2(i))*half
136 f32(i)=f32(i)-fint+fintz
137 f35(i)=f35(i)-fint-fintz
138C
139 fint=s1(i)*px3(i)+s4(i)*py3(i)
140 finsx=-s6(i)*(b1122(i)+b1x(i,1))+s5(i)*(b1121(i)+b2x(i,1))
141 f13(i)=f13(i)-fint+fintx+finsx
142 f16(i)=f16(i)-fint-fintx-finsx
143 fint=s2(i)*py3(i)+s4(i)*px3(i)
144 finsy=-s6(i)*(b2212(i)+b1y(i,1))+s5(i)*(b1221(i)+b2y(i,1))
145 f23(i)=f23(i)-fint+finty+finsy
146 f26(i)=f26(i)-fint-finty-finsy
147 fint=s3(i)*pz3(i)+(s6(i)*px3(i)+s5(i)*py3(i))*half
148 f33(i)=f33(i)-fint+fintz
149 f36(i)=f36(i)-fint-fintz
150 ENDDO
151C --- non constante part------------
152C 12
153 DO i=1,nel
154 s1(i) = zi*s1(i)
155 s2(i) = zi*s2(i)
156 s3(i) = zi*s3(i)
157 s4(i) = zi*s4(i)
158 s5(i) = zi*s5(i)
159 s6(i) = zi*s6(i)
160 ENDDO
161 DO i=1,nel
162 nu1 = nu(i)/(one - nu(i))
163 fxc =s1(i)-(s2(i)+s3(i))*nu(i)
164 fyc =s2(i)-(s1(i)+s3(i))*nu(i)
165 finsx = s1(i)*nu1
166 finsy = s2(i)*nu1
167 finsz = s3(i)*nu1
168 fint=fxc*px1h(i)+s4(i)*py1h(i)
169 fintx=(s1(i)-finsz)*px1(i)+s4(i)*py1(i)
170 fintx=fintx+s6(i)*(b1xh(i,1)-b1xh(i,2))
171 . -s5(i)*(b2xh(i,1)-b2xh(i,2))
172 f11(i)=f11(i)-fint+fintx
173 f14(i)=f14(i)-fint-fintx
174 fint=(s2(i)-finsz)*py1h(i)+s4(i)*px1h(i)
175 finty=fyc*py1(i)+s4(i)*px1(i)
176 finty=finty+s6(i)*(b1yh(i,1)-b1yh(i,2))
177 . -s5(i)*(b2yh(i,1)-b2yh(i,2))
178 f21(i)=f21(i)-fint+finty
179 f24(i)=f24(i)-fint-finty
180 fint=(s3(i)-finsy)*pz1h(i)+
181 . (s6(i)*px1h(i)+s5(i)*py1h(i))*half
182 fintz=(s3(i)-finsx)*pz1(i)
183 f31(i)=f31(i)-fint+fintz
184 f34(i)=f34(i)-fint-fintz
185C
186 fint=fxc*px2h(i)+s4(i)*py2h(i)
187 fintx=(s1(i)-finsz)*px2(i)+s4(i)*py2(i)
188 fintx=fintx+
189 . s6(i)*(b1221h(i)+b1xh(i,2))-s5(i)*(b1121h(i)+b2xh(i,2))
190 f12(i)=f12(i)-fint+fintx
191 f15(i)=f15(i)-fint-fintx
192 fint=(s2(i)-finsz)*py2h(i)+s4(i)*px2h(i)
193 finty=fyc*py2(i)+s4(i)*px2(i)
194 finty=finty+
195 . s6(i)*(b2212h(i)+b1yh(i,2))-s5(i)*(b1122h(i)+b2yh(i,2))
196 f22(i)=f22(i)-fint+finty
197 f25(i)=f25(i)-fint-finty
198 fint=(s3(i)-finsy)*pz2h(i)+
199 . (s6(i)*px2h(i)+s5(i)*py2h(i))*half
200 fintz=(s3(i)-finsx)*pz2(i)
201 f32(i)=f32(i)-fint+fintz
202 f35(i)=f35(i)-fint-fintz
203C
204 fint=fxc*px3h(i)+s4(i)*py3h(i)
205 fintx=(s1(i)-finsz)*px3(i)+s4(i)*py3(i)
206 fintx=fintx
207 . -s6(i)*(b1122h(i)+b1xh(i,1))+s5(i)*(b1121h(i)+b2xh(i,1))
208 f13(i)=f13(i)-fint+fintx
209 f16(i)=f16(i)-fint-fintx
210 fint=(s2(i)-finsz)*py3h(i)+s4(i)*px3h(i)
211 finty=fyc*py3(i)+s4(i)*px3(i)
212 finty=finty
213 . -s6(i)*(b2212h(i)+b1yh(i,1))+s5(i)*(b1221h(i)+b2yh(i,1))
214 f23(i)=f23(i)-fint+finty
215 f26(i)=f26(i)-fint-finty
216 fint=(s3(i)-finsy)*pz3h(i)+
217 . (s6(i)*px3h(i)+s5(i)*py3h(i))*half
218 fintz=(s3(i)-finsx)*pz3(i)
219 f33(i)=f33(i)-fint+fintz
220 f36(i)=f36(i)-fint-fintz
221 ENDDO
222C---------------------------------------------- /
223C - post-traitement-valeur moyenne au sens a'=(_/ a dv ) /v
224 DO i=1,nel
225 fac(i) = off(i)*vol(i)/volg(i)
226 sigm(i,1) = sigm(i,1) + fac(i) * sig(i,1)
227 sigm(i,2) = sigm(i,2) + fac(i) * sig(i,2)
228 sigm(i,3) = sigm(i,3) + fac(i) * sig(i,3)
229 sigm(i,4) = sigm(i,4) + fac(i) * sig(i,4)
230 sigm(i,5) = sigm(i,5) + fac(i) * sig(i,5)
231 sigm(i,6) = sigm(i,6) + fac(i) * sig(i,6)
232 rhom(i) = rhom(i) + fac(i) * rho(i)
233 eintm(i) = eintm(i) + eint(i)*vol0(i)/vol0g(i)
234 IF (g_wpla_flag > 0) g_wpla(i) = g_wpla(i) + l_wpla(i)
235 qm(i) = qm(i) + fac(i) * q(i)
236 ENDDO
237 IF (g_pla > 0) THEN
238 DO i=1,nel
239 eplasm(i) = eplasm(i) + fac(i) * epla(i)
240 ENDDO
241 ENDIF
242 IF (g_epsd > 0) THEN
243 DO i=1,nel
244 epsdm(i) = epsdm(i) + fac(i) * epsd(i)
245 ENDDO
246 ENDIF
247C-----------
248 RETURN
249 END
s6cfint3
subroutine s6cfint3(sig, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, px1h, px2h, px3h, py1h, py2h, py3h, pz1h, pz2h, pz3h, ji33, b1x, b1y, b2y, b2x, b1122, b1221, b2212, b1121, b1xh, b1yh, b2xh, b2yh, b1122h, b1221h, b2212h, b1121h, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, vol, qvis, eint, rho, q, epla, epsd, epsdm, sigm, eintm, rhom, qm, eplasm, zi, wi, volg, off, nu, vol0, vol0g, g_pla, g_epsd, nel, svis, g_wpla, l_wpla, g_wpla_flag)
Definition s6cfint3.F:50