OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s8fint3.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!|| s8fint3 ../engine/source/elements/solid/solide8/s8fint3.F
25!||--- called by ------------------------------------------------------
26!|| s8forc3 ../engine/source/elements/solid/solide8/s8forc3.F
27!||--- uses -----------------------------------------------------
28!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
29!||====================================================================
30 SUBROUTINE s8fint3(BUFLY,VOLGP,QVIS,
31 . PX1,PX2,PX3,PX4,PX5,PX6,PX7,PX8,
32 . PY1,PY2,PY3,PY4,PY5,PY6,PY7,PY8,
33 . PZ1,PZ2,PZ3,PZ4,PZ5,PZ6,PZ7,PZ8,NEL,
34 . F11, F12, F13, F14,
35 . F15, F16, F17, F18, F21,
36 . F22, F23, F24, F25, F26,
37 . F27, F28, F31, F32, F33,
38 . F34, F35, F36, F37, F38,SVIS)
39C-----------------------------------------------
40C M o d u l e s
41C-----------------------------------------------
42 USE elbufdef_mod
43C-----------------------------------------------
44C I m p l i c i t T y p e s
45C-----------------------------------------------
46#include "implicit_f.inc"
47C-----------------------------------------------
48C G l o b a l P a r a m e t e r s
49C-----------------------------------------------
50#include "mvsiz_p.inc"
51C-----------------------------------------------
52C C o m m o n B l o c k s
53C-----------------------------------------------
54C-----------------------------------------------
55C D u m m y A r g u m e n t s
56C-----------------------------------------------
57 INTEGER NEL
58C REAL
59 my_real VOLGP(MVSIZ,8),QVIS(NEL),
60 . PX1(MVSIZ,8),PX2(MVSIZ,8),PX3(MVSIZ,8),PX4(MVSIZ,8),
61 . PX5(MVSIZ,8),PX6(MVSIZ,8),PX7(MVSIZ,8),PX8(MVSIZ,8),
62 . PY1(MVSIZ,8),PY2(MVSIZ,8),PY3(MVSIZ,8),PY4(MVSIZ,8),
63 . PY5(MVSIZ,8),PY6(MVSIZ,8),PY7(MVSIZ,8),PY8(MVSIZ,8),
64 . PZ1(MVSIZ,8),PZ2(MVSIZ,8),PZ3(MVSIZ,8),PZ4(MVSIZ,8),
65 . pz5(mvsiz,8),pz6(mvsiz,8),pz7(mvsiz,8),pz8(mvsiz,8)
66 my_real :: f11(mvsiz),f12(mvsiz),f13(mvsiz),f14(mvsiz),
67 . f15(mvsiz), f16(mvsiz), f17(mvsiz), f18(mvsiz), f21(mvsiz),
68 . f22(mvsiz), f23(mvsiz), f24(mvsiz), f25(mvsiz), f26(mvsiz),
69 . f27(mvsiz), f28(mvsiz), f31(mvsiz), f32(mvsiz), f33(mvsiz),
70 . f34(mvsiz), f35(mvsiz), f36(mvsiz), f37(mvsiz), f38(mvsiz)
71 my_real, DIMENSION(MVSIZ,6), INTENT(INOUT) :: svis
72 TYPE (BUF_LAY_) ,TARGET :: BUFLY
73C-----------------------------------------------
74C L o c a l V a r i a b l e s
75C-----------------------------------------------
76 INTEGER I,J,IPT,II(6)
77C REAL
78 my_real
79 . s1, s2, s3,
80 . s4, s5, s6
81 TYPE(l_bufel_) ,POINTER :: LBUF
82C-----------------------------------------------
83!
84 DO i=1,6
85 ii(i) = nel*(i-1)
86 ENDDO
87!
88 DO i=1,nel
89 f11(i)=zero
90 f21(i)=zero
91 f31(i)=zero
92 f12(i)=zero
93 f22(i)=zero
94 f32(i)=zero
95 f13(i)=zero
96 f23(i)=zero
97 f33(i)=zero
98 f14(i)=zero
99 f24(i)=zero
100 f34(i)=zero
101 f15(i)=zero
102 f25(i)=zero
103 f35(i)=zero
104 f16(i)=zero
105 f26(i)=zero
106 f36(i)=zero
107 f17(i)=zero
108 f27(i)=zero
109 f37(i)=zero
110 f18(i)=zero
111 f28(i)=zero
112 f38(i)=zero
113 ENDDO
114C
115 DO ipt=1,8
116 lbuf => bufly%LBUF(1,1,ipt)
117 DO i=1,nel
118 j = 6*(i-1)
119 s1=(lbuf%SIG(ii(1)+i)+svis(i,1)-qvis(i))*volgp(i,ipt)
120 s2=(lbuf%SIG(ii(2)+i)+svis(i,2)-qvis(i))*volgp(i,ipt)
121 s3=(lbuf%SIG(ii(3)+i)+svis(i,3)-qvis(i))*volgp(i,ipt)
122 s4=(lbuf%SIG(ii(4)+i)+svis(i,4))*volgp(i,ipt)
123 s5=(lbuf%SIG(ii(5)+i)+svis(i,5))*volgp(i,ipt)
124 s6=(lbuf%SIG(ii(6)+i)+svis(i,6))*volgp(i,ipt)
125C
126 f11(i)=f11(i)
127 - -(s1*px1(i,ipt)+s4*py1(i,ipt)+s6*pz1(i,ipt))
128 f21(i)=f21(i)
129 - -(s2*py1(i,ipt)+s4*px1(i,ipt)+s5*pz1(i,ipt))
130 f31(i)=f31(i)
131 - -(s3*pz1(i,ipt)+s6*px1(i,ipt)+s5*py1(i,ipt))
132 f12(i)=f12(i)
133 - -(s1*px2(i,ipt)+s4*py2(i,ipt)+s6*pz2(i,ipt))
134 f22(i)=f22(i)
135 - -(s2*py2(i,ipt)+s4*px2(i,ipt)+s5*pz2(i,ipt))
136 f32(i)=f32(i)
137 - -(s3*pz2(i,ipt)+s6*px2(i,ipt)+s5*py2(i,ipt))
138 f13(i)=f13(i)
139 - -(s1*px3(i,ipt)+s4*py3(i,ipt)+s6*pz3(i,ipt))
140 f23(i)=f23(i)
141 - -(s2*py3(i,ipt)+s4*px3(i,ipt)+s5*pz3(i,ipt))
142 f33(i)=f33(i)
143 - -(s3*pz3(i,ipt)+s6*px3(i,ipt)+s5*py3(i,ipt))
144 f14(i)=f14(i)
145 - -(s1*px4(i,ipt)+s4*py4(i,ipt)+s6*pz4(i,ipt))
146 f24(i)=f24(i)
147 - -(s2*py4(i,ipt)+s4*px4(i,ipt)+s5*pz4(i,ipt))
148 f34(i)=f34(i)
149 - -(s3*pz4(i,ipt)+s6*px4(i,ipt)+s5*py4(i,ipt))
150 f15(i)=f15(i)
151 - -(s1*px5(i,ipt)+s4*py5(i,ipt)+s6*pz5(i,ipt))
152 f25(i)=f25(i)
153 - -(s2*py5(i,ipt)+s4*px5(i,ipt)+s5*pz5(i,ipt))
154 f35(i)=f35(i)
155 - -(s3*pz5(i,ipt)+s6*px5(i,ipt)+s5*py5(i,ipt))
156 f16(i)=f16(i)
157 - -(s1*px6(i,ipt)+s4*py6(i,ipt)+s6*pz6(i,ipt))
158 f26(i)=f26(i)
159 - -(s2*py6(i,ipt)+s4*px6(i,ipt)+s5*pz6(i,ipt))
160 f36(i)=f36(i)
161 - -(s3*pz6(i,ipt)+s6*px6(i,ipt)+s5*py6(i,ipt))
162 f17(i)=f17(i)
163 - -(s1*px7(i,ipt)+s4*py7(i,ipt)+s6*pz7(i,ipt))
164 f27(i)=f27(i)
165 - -(s2*py7(i,ipt)+s4*px7(i,ipt)+s5*pz7(i,ipt))
166 f37(i)=f37(i)
167 - -(s3*pz7(i,ipt)+s6*px7(i,ipt)+s5*py7(i,ipt))
168 f18(i)=f18(i)
169 - -(s1*px8(i,ipt)+s4*py8(i,ipt)+s6*pz8(i,ipt))
170 f28(i)=f28(i)
171 - -(s2*py8(i,ipt)+s4*px8(i,ipt)+s5*pz8(i,ipt))
172 f38(i)=f38(i)
173 - -(s3*pz8(i,ipt)+s6*px8(i,ipt)+s5*py8(i,ipt))
174 ENDDO
175C
176 ENDDO
177c-----------
178 RETURN
179 END
s8fint3
subroutine s8fint3(bufly, volgp, qvis, px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, nel, f11, f12, f13, f14, f15, f16, f17, f18, f21, f22, f23, f24, f25, f26, f27, f28, f31, f32, f33, f34, f35, f36, f37, f38, svis)
Definition s8fint3.F:39