OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
mhvis3.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!|| mhvis3 ../engine/source/airbag/mhvis3.F
25!||--- called by ------------------------------------------------------
26!|| cforc3 ../engine/source/elements/shell/coque/cforc3.F
27!|| cforc3_crk ../engine/source/elements/xfem/cforc3_crk.F
28!||====================================================================
29 SUBROUTINE mhvis3(JFT,JLT ,PM ,THK,HOUR,
30 2 OFF,PX1 ,PX2,PY1,PY2 ,
31 3 IXC,DT1C,SSP,RHO,STI ,
32 4 EANI,GEO ,PID,STIR,MAT,
33 5 THK0,VISCMX,ALPE,IPARTC ,PARTSAV,
34 6 IHBE ,NFT ,ISMSTR , RX1,
35 7 RX2,RX3,RX4,RY1,RY2,
36 8 RY3,RY4,VX1,VX2,VX3,
37 9 VX4,VY1,VY2,VY3,VY4,
38 A VZ1,VZ2,VZ3,VZ4,B11,
39 B B12,B13,B14,B21,B22,
40 C B23,B24,AREA,YM,PR,
41 D VHX,VHY,H11,H12,H13,
42 E H14,H21,H22,H23,H24,
43 F H31,H32,H33,H34,H1 ,
44 G H2,IGEO,NEL,MTN,A1 )
45C-----------------------------------------------
46C I m p l i c i t T y p e s
47C-----------------------------------------------
48#include "implicit_f.inc"
49#include "comlock.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 C o m m o n B l o c k s
56C-----------------------------------------------
57#include "com04_c.inc"
58#include "com06_c.inc"
59#include "param_c.inc"
60#include "scr02_c.inc"
61#include "scr06_c.inc"
62#include "scr14_c.inc"
63#include "scr16_c.inc"
64#include "scr17_c.inc"
65#include "sms_c.inc"
66C-----------------------------------------------
67C D u m m y A r g u m e n t s
68C-----------------------------------------------
69 INTEGER IXC(NIXC,*),IPARTC(*), JFT, JLT,PID(*),
70 . IHBE ,NFT ,ISMSTR,IGEO(NPROPGI, *),NEL,MTN
71 INTEGER MAT(MVSIZ)
72C REAL
73 my_real
74 . PM(NPROPM,*), GEO(NPROPG,*), THK(*), HOUR(NEL,5), OFF(*),
75 . PX1(*), PX2(*), PY1(*), PY2(*),DT1C(*),EANI(*),
76 . SSP(MVSIZ), RHO(MVSIZ),STI(MVSIZ),STIR(*),
77 . H1(MVSIZ), H2(MVSIZ),
78 . THK0(MVSIZ),VISCMX(MVSIZ), ALPE(MVSIZ),PARTSAV(NPSAV,*)
79C REAL
80 my_real
81 . B11(MVSIZ), B12(MVSIZ), B13(MVSIZ), B14(MVSIZ), B21(MVSIZ),
82 . B22(MVSIZ), B23(MVSIZ), B24(MVSIZ), H11(MVSIZ), H12(MVSIZ),
83 . H13(MVSIZ), H14(MVSIZ), H21(MVSIZ), H22(MVSIZ), H23(MVSIZ),
84 . h24(mvsiz), h31(mvsiz), h32(mvsiz), h33(mvsiz), h34(mvsiz),
85 . rx1(mvsiz), rx2(mvsiz), rx3(mvsiz), rx4(mvsiz), ry1(mvsiz),
86 . ry2(mvsiz), ry3(mvsiz), ry4(mvsiz), vhx(mvsiz), vhy(mvsiz),
87 . vx1(mvsiz), vx2(mvsiz), vx3(mvsiz),
88 . vx4(mvsiz), vy1(mvsiz), vy2(mvsiz), vy3(mvsiz), vy4(mvsiz),
89 . vz1(mvsiz), vz2(mvsiz), vz3(mvsiz), vz4(mvsiz), area(mvsiz),
90 . a11r(mvsiz)
91 my_real,DIMENSION(MVSIZ), INTENT(IN) :: a1
92C-----------------------------------------------
93C L o c a l V a r i a b l e s
94C-----------------------------------------------
95 INTEGER I, MX,IPID,IGTYP,IGMAT,IPGMAT
96 my_real H1L(MVSIZ), H2L(MVSIZ), H1Q(MVSIZ), H2Q(MVSIZ), HG1(MVSIZ),
97 . HG2(MVSIZ), FAC(MVSIZ), YM(MVSIZ), PR(MVSIZ),
98 . GAMA1(MVSIZ), GAMA2(MVSIZ), GAMA3(MVSIZ), GAMA4(MVSIZ),
99 . H4(MVSIZ), H4L(MVSIZ), H4Q(MVSIZ),THK02(MVSIZ),
100 . G(MVSIZ) , B1(MVSIZ), B2(MVSIZ),A11(MVSIZ),EHOU(MVSIZ),
101 . PX1V, PX2V, PY1V, PY2V, EHOURT,VV ,SCALE(MVSIZ),FAC1
102C-----------------------------------------------
103 IF(ISMSTR /=3 .AND. ihbe >= 1)THEN
104 DO i=jft,jlt
105 px1v = px1(i)*vhx(i)
106 px2v = px2(i)*vhx(i)
107 py1v = py1(i)*vhy(i)
108 py2v = py2(i)*vhy(i)
109 gama1(i)= off(i)*( one- px1v-py1v)
110 gama3(i)= off(i)*( one+ px1v+py1v)
111 gama2(i)= off(i)*(-one- px2v-py2v)
112 gama4(i)= off(i)*(-one+ px2v+py2v)
113 ENDDO
114 ELSE
115 DO i=jft,jlt
116 gama1(i)= off(i)
117 gama3(i)= off(i)
118 gama2(i)= -off(i)
119 gama4(i)= -off(i)
120 ENDDO
121 ENDIF
122C
123 IF(invstr >= 35)THEN
124 mx=mat(jft)
125 DO i=jft,jlt
126 ipid=pid(i)
127 h4(i) =geo(17,ipid)
128 ENDDO
129C
130 ELSE
131 mx=mat(jft)
132 DO i=jft,jlt
133 h4(i) =pm(91,mx)
134 ENDDO
135C
136 ENDIF
137C
138 DO i=jft,jlt
139 fac(i)=fourth*rho(i)*thk0(i)
140 h1l(i)= zero
141 h1q(i)= zero
142 h2l(i)= zero
143 h2q(i)= zero
144 h4l(i)=fac(i)*sqrt(hvisc*h4(i)*area(i))
145 h4q(i)=sqrt(hvisc*h4(i))*h4l(i)*hundred
146 h4l(i)=h4l(i)*ssp(i)
147 ENDDO
148C
149 DO i=jft,jlt
150 thk02(i)= thk0(i)*thk0(i)
151 b1(i) = px1(i)*px1(i)+py1(i)*py1(i)
152 b2(i) = px2(i)*px2(i)+py2(i)*py2(i)
153 fac(i)=fourth*ym(i)*thk0(i)*dt1c(i)*helas
154 h1(i)=h1(i)*fac(i)
155 h2(i)=h2(i)*fac(i)
156 ENDDO
157C-----------------------------------
158C TRIANGLES
159C-----------------------------------
160 DO i=jft,jlt
161 IF(ixc(4,i)/=ixc(5,i))cycle
162 h1q(i)=zero
163 h1l(i)=zero
164 h2q(i)=zero
165 h2l(i)=zero
166 h4l(i)=zero
167 h4q(i)=zero
168 h1(i)=zero
169 h2(i)=zero
170 ENDDO
171C-----------------------------------
172C STIFFNESS - DT
173C-----------------------------------
174 DO i=jft,jlt
175 scale(i)= zero
176 ENDDO
177 ipid=pid(1)
178 igtyp = igeo(11,ipid)
179 igmat = igeo(98,ipid)
180 ipgmat = 700
181 IF(nodadt /= 0 .OR. idt1sh == 1.OR. idtmins == 2)THEN
182C
183 DO i=jft,jlt
184 scale(i)= max(gama1(i)*gama1(i),gama2(i)*gama2(i),gama3(i)*gama3(i),gama4(i)*gama4(i)) *
185 . dt1c(i)*max(h1(i)+h1l(i),h2(i)+h2l(i),h4l(i)) / max(dt1c(i)*dt1c(i),em20)
186 sti(i)=sti(i) + scale(i)
187 ENDDO
188 IF(igtyp == 11 .AND. igmat > 0) THEN
189 DO i=jft,jlt
190 a11(i) = geo(ipgmat +5 ,pid(i))
191 g(i) = geo(ipgmat+4,pid(i))
192 a11r(i) = geo(ipgmat+7,pid(i))
193 IF (off(i)==zero) THEN
194 sti(i) = zero
195 stir(i) = zero
196 ELSE
197 vv = viscmx(i) * viscmx(i) * alpe(i)
198 fac1 = max(b1(i),b2(i)) / (area(i) * vv)
199 sti(i) = sti(i) + fac1* thk0(i) * a11(i)
200 stir(i) = fac1 * a11r(i)*one_over_12*thk0(i)**3 +
201 . fac1 * a11(i)*thk0(i)*area(i)*one_over_9 +
202 . fac1*scale(i)*(one_over_12*thk0(i)**2 + area(i)*one_over_9)
203 ENDIF
204 ENDDO
205 ELSE !
206 mx = mat(jft)
207 DO i=jft,jlt
208 a11(i) =pm(24,mx)
209 g(i) =pm(22,mx)
210 ENDDO
211 IF (mtn==58) a11(jft:jlt)=a1(jft:jlt)
212 DO i=jft,jlt
213 IF (off(i)==zero) THEN
214 sti(i) = zero
215 stir(i) = zero
216 ELSE
217 vv = viscmx(i) * viscmx(i) * alpe(i)
218 sti(i) = sti(i) + max(b1(i),b2(i)) * thk0(i) * a11(i) / (area(i) * vv)
219 stir(i) = sti(i)*(thk02(i) * one_over_12 + area(i) * one_over_9)
220 ENDIF
221 ENDDO
222 ENDIF
223 ENDIF
224C-----------------------------------
225C ANTI-HOURGLASS MEMBRANE FORCES
226C-----------------------------------
227 IF(ismstr == 3 .OR. ihbe < 1)THEN
228 DO i=jft,jlt
229 hg1(i)=vx1(i)-vx2(i)+vx3(i)-vx4(i)
230 hg2(i)=vy1(i)-vy2(i)+vy3(i)-vy4(i)
231 ENDDO
232
233 DO i=jft,jlt
234 hour(i,1)=hour(i,1)+hg1(i)*h1(i)
235 hour(i,2)=hour(i,2)+hg2(i)*h1(i)
236 hg1(i)=hg1(i)*(h1l(i)+h1q(i)*abs(hg1(i)))
237 hg2(i)=hg2(i)*(h1l(i)+h1q(i)*abs(hg2(i)))
238 h11(i)= hour(i,1)+hg1(i)
239 h12(i)=-hour(i,1)-hg1(i)
240 h13(i)= hour(i,1)+hg1(i)
241 h14(i)=-hour(i,1)-hg1(i)
242 h21(i)= hour(i,2)+hg2(i)
243 h22(i)=-hour(i,2)-hg2(i)
244 h23(i)= hour(i,2)+hg2(i)
245 h24(i)=-hour(i,2)-hg2(i)
246 ENDDO
247 ELSE
248 DO i=jft,jlt
249 hg1(i)=vx1(i)*gama1(i)+vx2(i)*gama2(i)+vx3(i)*gama3(i)+vx4(i)*gama4(i)
250 hg2(i)=vy1(i)*gama1(i)+vy2(i)*gama2(i)+vy3(i)*gama3(i)+vy4(i)*gama4(i)
251 ENDDO
252 DO i=jft,jlt
253 hour(i,1)=hour(i,1)+hg1(i)*h1(i)
254 hour(i,2)=hour(i,2)+hg2(i)*h1(i)
255 hg1(i)=hg1(i)*(h1l(i)+h1q(i)*abs(hg1(i)))
256 hg2(i)=hg2(i)*(h1l(i)+h1q(i)*abs(hg2(i)))
257 h11(i)=(hour(i,1)+hg1(i))*gama1(i)
258 h12(i)=(hour(i,1)+hg1(i))*gama2(i)
259 h13(i)=(hour(i,1)+hg1(i))*gama3(i)
260 h14(i)=(hour(i,1)+hg1(i))*gama4(i)
261 h21(i)=(hour(i,2)+hg2(i))*gama1(i)
262 h22(i)=(hour(i,2)+hg2(i))*gama2(i)
263 h23(i)=(hour(i,2)+hg2(i))*gama3(i)
264 h24(i)=(hour(i,2)+hg2(i))*gama4(i)
265 ENDDO
266 ENDIF
267 ehourt = 0.
268 DO i=jft,jlt
269 ehou(i) = vx1(i)*h11(i) + vx2(i)*h12(i) + vx3(i)*h13(i) + vx4(i)*h14(i)
270 . + vy1(i)*h21(i) + vy2(i)*h22(i) + vy3(i)*h23(i) + vy4(i)*h24(i)
271 ENDDO
272C----------------------------------------------------------------
273C HYPERBOLIC SHAPE ANTI-HOURGLASS BENDING FORCES
274C----------------------------------------------------------------
275 IF(ismstr==3.OR.ihbe<1)THEN
276 DO i=jft,jlt
277 hg1(i)=vz1(i)-vz2(i)+vz3(i)-vz4(i)
278 ENDDO
279
280 DO i=jft,jlt
281 hour(i,3)=hour(i,3)+hg1(i)*h2(i)
282 hg1(i)=hg1(i)*(h2l(i)+h2q(i)*abs(hg1(i)))
283 h31(i)= hour(i,3)+hg1(i)
284 h32(i)=-hour(i,3)-hg1(i)
285 h33(i)= hour(i,3)+hg1(i)
286 h34(i)=-hour(i,3)-hg1(i)
287 ENDDO
288 ELSE
289 DO i=jft,jlt
290 hg1(i)=vz1(i)*gama1(i)+vz2(i)*gama2(i)+vz3(i)*gama3(i)+vz4(i)*gama4(i)
291 ENDDO
292 DO i=jft,jlt
293 hour(i,3)=hour(i,3)+hg1(i)*h2(i)
294 hg1(i)=hg1(i)*(h2l(i)+h2q(i)*abs(hg1(i)))
295 h31(i)=(hour(i,3)+hg1(i))*gama1(i)
296 h32(i)=(hour(i,3)+hg1(i))*gama2(i)
297 h33(i)=(hour(i,3)+hg1(i))*gama3(i)
298 h34(i)=(hour(i,3)+hg1(i))*gama4(i)
299 ENDDO
300 ENDIF
301C------------------------------------------------------
302C V-SHAPE ANTI-HOURGLASS BENDING FORCES
303C------------------------------------------------------
304 DO i=jft,jlt
305 hg1(i)=+vz1(i)+vz2(i)-vz3(i)-vz4(i)
306 ENDDO
307C
308 DO i=jft,jlt
309 hg1(i)=hg1(i)*(h4l(i)+h4q(i)*abs(hg1(i)))
310 h31(i)=h31(i) +hg1(i)
311 h32(i)=h32(i) +hg1(i)
312 h33(i)=h33(i) -hg1(i)
313 h34(i)=h34(i) -hg1(i)
314 ENDDO
315C
316 DO i=jft,jlt
317 hg1(i)=vz1(i)-vz2(i)-vz3(i)+vz4(i)
318 ENDDO
319C
320 DO i=jft,jlt
321 hg1(i)=hg1(i)*(h4l(i)+h4q(i)*abs(hg1(i)))
322 h31(i)=h31(i) +hg1(i)
323 h32(i)=h32(i) -hg1(i)
324 h33(i)=h33(i) -hg1(i)
325 h34(i)=h34(i) +hg1(i)
326 ENDDO
327C------------------------------------------------
328 DO i=jft,jlt
329 b11(i)= zero
330 b12(i)= zero
331 b13(i)= zero
332 b14(i)= zero
333 b21(i)= zero
334 b22(i)= zero
335 b23(i)= zero
336 b24(i)= zero
337 ENDDO
338C------------------------------------------------
339 DO i=jft,jlt
340 ehou(i) = ehou(i) + vz1(i)*h31(i) + vz2(i)*h32(i) + vz3(i)*h33(i) + vz4(i)*h34(i)
341 ehou(i) = dt1c(i) * ehou(i)
342 ehourt = ehourt + ehou(i)
343 ENDDO
344
345 DO i=jft,jlt
346 mx = ipartc(i)
347 partsav(8,mx)=partsav(8,mx) + ehou(i)
348 ENDDO
349
350!$OMP ATOMIC
351 ehour = ehour + ehourt
352
353 DO i=jft,jlt
354 eani(nft+numels+i) = eani(nft+numels+i)+ehou(i)
355 ENDDO
356 RETURN
357 END
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
max
#define max(a, b)
Definition macros.h:21
mhvis3
subroutine mhvis3(jft, jlt, pm, thk, hour, off, px1, px2, py1, py2, ixc, dt1c, ssp, rho, sti, eani, geo, pid, stir, mat, thk0, viscmx, alpe, ipartc, partsav, ihbe, nft, ismstr, rx1, rx2, rx3, rx4, ry1, ry2, ry3, ry4, vx1, vx2, vx3, vx4, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, b11, b12, b13, b14, b21, b22, b23, b24, area, ym, pr, vhx, vhy, h11, h12, h13, h14, h21, h22, h23, h24, h31, h32, h33, h34, h1, h2, igeo, nel, mtn, a1)
Definition mhvis3.F:45