OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
fdefo3.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!|| fdefo3 ../engine/source/elements/solid/solide/fdefo3.F
25!||--- called by ------------------------------------------------------
26!|| sforc3 ../engine/source/elements/solid/solide/sforc3.F
27!||====================================================================
28 SUBROUTINE fdefo3(
29 1 PX1, PX2, PX3, PX4,
30 2 PY1, PY2, PY3, PY4,
31 3 PZ1, PZ2, PZ3, PZ4,
32 4 VX1, VX2, VX3, VX4,
33 5 VX5, VX6, VX7, VX8,
34 6 VY1, VY2, VY3, VY4,
35 7 VY5, VY6, VY7, VY8,
36 8 VZ1, VZ2, VZ3, VZ4,
37 9 VZ5, VZ6, VZ7, VZ8,
38 A DXX, DXY, DXZ, DYX,
39 B DYY, DYZ, DZX, DZY,
40 C DZZ, D4, D5, D6,
41 D WXX, WYY, WZZ, VXR,
42 E VYR, VZR, DD, NEL,
43 F JHBE, JCVT, ISROT)
44C-----------------------------------------------
45C I m p l i c i t T y p e s
46C-----------------------------------------------
47#include "implicit_f.inc"
48C-----------------------------------------------
49C G l o b a l P a r a m e t e r s
50C-----------------------------------------------
51#include "mvsiz_p.inc"
52C-----------------------------------------------
53C C o m m o n B l o c k s
54C-----------------------------------------------
55#include "com08_c.inc"
56#include "impl1_c.inc"
57C-----------------------------------------------
58C D u m m y A r g u m e n t s
59C-----------------------------------------------
60 INTEGER, INTENT(IN) :: NEL
61 INTEGER, INTENT(IN) :: JHBE
62 INTEGER, INTENT(IN) :: JCVT
63 INTEGER, INTENT(IN) :: ISROT
64 my_real
65 . VX1(*), VX2(*), VX3(*), VX4(*),
66 . VX5(*), VX6(*), VX7(*), VX8(*),
67 . VY1(*), VY2(*), VY3(*), VY4(*),
68 . VY5(*), VY6(*), VY7(*), VY8(*),
69 . VZ1(*), VZ2(*), VZ3(*), VZ4(*),
70 . VZ5(*), VZ6(*), VZ7(*), VZ8(*),
71 . PX1(*), PX2(*), PX3(*), PX4(*),
72 . PY1(*), PY2(*), PY3(*), PY4(*),
73 . PZ1(*), PZ2(*), PZ3(*), PZ4(*),
74 . DXX(*), DXY(*), DXZ(*),
75 . dyx(*), dyy(*), dyz(*),
76 . dzx(*), dzy(*), dzz(*), d4(*), d5(*), d6(*),
77 . wxx(*), wyy(*), wzz(*), dd(*)
78 my_real
79 . vxr(*),vyr(*),vzr(*)
80C-----------------------------------------------
81C L o c a l V a r i a b l e s
82C-----------------------------------------------
83 INTEGER I, J
84 my_real
85 . VX17(MVSIZ), VY17(MVSIZ), VZ17(MVSIZ),
86 . VX28(MVSIZ), VY28(MVSIZ), VZ28(MVSIZ),
87 . VX35(MVSIZ), VY35(MVSIZ), VZ35(MVSIZ),
88 . VX46(MVSIZ), VY46(MVSIZ), VZ46(MVSIZ),
89 . DT1D2, DT1D
90 my_real
91 . PXX2,PYY2,PZZ2,PXX2P,PYY2P,PZZ2P
92 my_real
93 . aaa,bbb
94C-----------------------------------------------
95 DO i=1,nel
96 vx17(i)=vx1(i)-vx7(i)
97 vx28(i)=vx2(i)-vx8(i)
98 vx35(i)=vx3(i)-vx5(i)
99 vx46(i)=vx4(i)-vx6(i)
100 vy17(i)=vy1(i)-vy7(i)
101 vy28(i)=vy2(i)-vy8(i)
102 vy35(i)=vy3(i)-vy5(i)
103 vy46(i)=vy4(i)-vy6(i)
104 vz17(i)=vz1(i)-vz7(i)
105 vz28(i)=vz2(i)-vz8(i)
106 vz35(i)=vz3(i)-vz5(i)
107 vz46(i)=vz4(i)-vz6(i)
108 ENDDO
109C
110 DO i=1,nel
111 dxx(i)=px1(i)*vx17(i)+px2(i)*vx28(i)+
112 . px3(i)*vx35(i)+px4(i)*vx46(i)
113 dyy(i)=py1(i)*vy17(i)+py2(i)*vy28(i)+
114 . py3(i)*vy35(i)+py4(i)*vy46(i)
115 dzz(i)=pz1(i)*vz17(i)+pz2(i)*vz28(i)+
116 . pz3(i)*vz35(i)+pz4(i)*vz46(i)
117 dxy(i)=py1(i)*vx17(i)+py2(i)*vx28(i)+
118 . py3(i)*vx35(i)+py4(i)*vx46(i)
119 dxz(i)=pz1(i)*vx17(i)+pz2(i)*vx28(i)+
120 . pz3(i)*vx35(i)+pz4(i)*vx46(i)
121 dyx(i)=px1(i)*vy17(i)+px2(i)*vy28(i)+
122 . px3(i)*vy35(i)+px4(i)*vy46(i)
123 dyz(i)=pz1(i)*vy17(i)+pz2(i)*vy28(i)+
124 . pz3(i)*vy35(i)+pz4(i)*vy46(i)
125 dzx(i)=px1(i)*vz17(i)+px2(i)*vz28(i)+
126 . px3(i)*vz35(i)+px4(i)*vz46(i)
127 dzy(i)=py1(i)*vz17(i)+py2(i)*vz28(i)+
128 . py3(i)*vz35(i)+py4(i)*vz46(i)
129 ENDDO
130C
131C
132 dt1d2=half*dt1
133 IF (iscau>0)dt1d2=dt1
134C
135 IF (jcvt > 0) THEN
136 IF (impl_s==0.OR.idyna > 0) THEN
137 DO i=1,nel
138 dxx(i) = dxx(i)
139 . -dt1d2*(dxx(i)*dxx(i)+dyx(i)*dyx(i)+dzx(i)*dzx(i))
140 dyy(i) = dyy(i)
141 . -dt1d2*(dyy(i)*dyy(i)+dzy(i)*dzy(i)+dxy(i)*dxy(i))
142 dzz(i) = dzz(i)
143 . -dt1d2*(dzz(i)*dzz(i)+dxz(i)*dxz(i)+dyz(i)*dyz(i))
144 aaa = dt1d2*(dxx(i)*dxy(i)+dyx(i)*dyy(i)+dzx(i)*dzy(i))
145 dxy(i) = dxy(i) - aaa
146 dyx(i) = dyx(i) - aaa
147 d4(i) = dxy(i)+dyx(i)
148 aaa = dt1d2*(dyy(i)*dyz(i)+dzy(i)*dzz(i)+dxy(i)*dxz(i))
149 dyz(i) = dyz(i) - aaa
150 dzy(i) = dzy(i) - aaa
151 d5(i) = dyz(i)+dzy(i)
152 aaa = dt1d2*(dzz(i)*dzx(i)+dxz(i)*dxx(i)+dyz(i)*dyx(i))
153 dxz(i) = dxz(i) - aaa
154 dzx(i) = dzx(i) - aaa
155 d6(i) = dxz(i)+dzx(i)
156 wxx(i)=zero
157 wyy(i)=zero
158 wzz(i)=zero
159 ENDDO
160C-------implicit static---------
161 ELSEIF (iline>0.AND.iscau==0) THEN
162 DO i=1,nel
163 d4(i) = dxy(i)+dyx(i)
164 d5(i) = dyz(i)+dzy(i)
165 d6(i) = dxz(i)+dzx(i)
166 wxx(i)=zero
167 wyy(i)=zero
168 wzz(i)=zero
169 ENDDO
170 ELSE
171 dt1d=two*dt1d2
172 DO i=1,nel
173 d4(i) = dxy(i)+dyx(i)
174 . -dt1d*(dxx(i)*dxy(i)+dyx(i)*dyy(i)+dzx(i)*dzy(i))
175 d5(i) = dyz(i)+dzy(i)
176 . -dt1d*(dyy(i)*dyz(i)+dzy(i)*dzz(i)+dxy(i)*dxz(i))
177 d6(i) = dxz(i)+dzx(i)
178 . -dt1d*(dzz(i)*dzx(i)+dxz(i)*dxx(i)+dyz(i)*dyx(i))
179 dxx(i) = dxx(i)
180 . -dt1d2*(dxx(i)*dxx(i)+dyx(i)*dyx(i)+dzx(i)*dzx(i))
181 dyy(i) = dyy(i)
182 . -dt1d2*(dyy(i)*dyy(i)+dzy(i)*dzy(i)+dxy(i)*dxy(i))
183 dzz(i) = dzz(i)
184 . -dt1d2*(dzz(i)*dzz(i)+dxz(i)*dxz(i)+dyz(i)*dyz(i))
185 wxx(i)=zero
186 wyy(i)=zero
187 wzz(i)=zero
188 ENDDO
189 ENDIF ! IF (IMPL_S==0.OR.IDYNA>0)
190C
191 ELSEIF (jhbe >= 2) THEN
192 DO i=1,nel
193 dxx(i) = dxx(i)
194 . -dt1d2*(dxx(i)*dxx(i)+dyx(i)*dyx(i)+dzx(i)*dzx(i))
195 dyy(i) = dyy(i)
196 . -dt1d2*(dyy(i)*dyy(i)+dzy(i)*dzy(i)+dxy(i)*dxy(i))
197 dzz(i) = dzz(i)
198 . -dt1d2*(dzz(i)*dzz(i)+dxz(i)*dxz(i)+dyz(i)*dyz(i))
199 aaa = dt1d2*(dxx(i)*dxy(i)+dyx(i)*dyy(i)+dzx(i)*dzy(i))
200 dxy(i) = dxy(i) - aaa
201 dyx(i) = dyx(i) - aaa
202 d4(i) = dxy(i)+dyx(i)
203 aaa = dt1d2*(dyy(i)*dyz(i)+dzy(i)*dzz(i)+dxy(i)*dxz(i))
204 dyz(i) = dyz(i) - aaa
205 dzy(i) = dzy(i) - aaa
206 d5(i) = dyz(i)+dzy(i)
207 aaa = dt1d2*(dzz(i)*dzx(i)+dxz(i)*dxx(i)+dyz(i)*dyx(i))
208 dxz(i) = dxz(i) - aaa
209 dzx(i) = dzx(i) - aaa
210 d6(i) = dxz(i)+dzx(i)
211 pxx2 = px1(i)*px1(i)+px2(i)*px2(i)
212 . +px3(i)*px3(i)+px4(i)*px4(i)
213 pyy2 = py1(i)*py1(i)+py2(i)*py2(i)
214 . +py3(i)*py3(i)+py4(i)*py4(i)
215 pzz2 = pz1(i)*pz1(i)+pz2(i)*pz2(i)
216 . +pz3(i)*pz3(i)+pz4(i)*pz4(i)
217 wzz(i)=dt1*(pyy2*dyx(i)-pxx2*dxy(i))/(pxx2+pyy2)
218 wxx(i)=dt1*(pzz2*dzy(i)-pyy2*dyz(i))/(pyy2+pzz2)
219 wyy(i)=dt1*(pxx2*dxz(i)-pzz2*dzx(i))/(pzz2+pxx2)
220 ENDDO
221C
222 ELSE
223 DO i=1,nel
224 d4(i) = dxy(i)+dyx(i)
225 d5(i) = dyz(i)+dzy(i)
226 d6(i) = dxz(i)+dzx(i)
227 dxx(i) = dxx(i)
228 dyy(i) = dyy(i)
229 dzz(i) = dzz(i)
230 wzz(i)=dt1d2*(dyx(i)-dxy(i))
231 wyy(i)=dt1d2*(dxz(i)-dzx(i))
232 wxx(i)=dt1d2*(dzy(i)-dyz(i))
233 ENDDO
234 ENDIF
235C
236
237 IF(isrot > 0)THEN
238 DO i=1,nel
239 dxy(i)=dxy(i) + vzr(i)
240 dyx(i)=dyx(i) - vzr(i)
241 dyz(i)=dyz(i) + vxr(i)
242 dzy(i)=dzy(i) - vxr(i)
243 dzx(i)=dzx(i) + vyr(i)
244 dxz(i)=dxz(i) - vyr(i)
245 ENDDO
246 aaa = 0.
247 bbb = dt1 - aaa*dt1
248 DO i=1,nel
249 wxx(i)=aaa*wxx(i) + bbb*vxr(i)
250 wyy(i)=aaa*wyy(i) + bbb*vyr(i)
251 wzz(i)=aaa*wzz(i) + bbb*vzr(i)
252 ENDDO
253 ENDIF
254C pour le calcul du pas du temps
255 DO i=1, nel
256 dd(i) = - dxx(i)- dyy(i) - dzz(i)
257 ENDDO
258C
259 RETURN
260 END SUBROUTINE fdefo3
fdefo3
subroutine fdefo3(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, dxy, dxz, dyx, dyy, dyz, dzx, dzy, dzz, d4, d5, d6, wxx, wyy, wzz, vxr, vyr, vzr, dd, nel, jhbe, jcvt, isrot)
Definition fdefo3.F:44