OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
int2rupt.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!|| int2rupt ../engine/source/interfaces/interf/int2rupt.F
25!||--- called by ------------------------------------------------------
26!|| intti1 ../engine/source/interfaces/interf/intti1.F
27!||--- calls -----------------------------------------------------
28!|| i2rupt ../engine/source/interfaces/interf/int2rupt.F
29!||--- uses -----------------------------------------------------
30!|| h3d_mod ../engine/share/modules/h3d_mod.F
31!|| intbufdef_mod ../common_source/modules/interfaces/intbufdef_mod.F90
32!||====================================================================
33 SUBROUTINE int2rupt(IPARI ,MS ,IN ,
34 . X ,V ,A ,STIFN ,IGEO ,
35 . WEIGHT ,FSAV ,ILEV ,NPF ,TF ,
36 . ITAB ,FNCONT ,PDAMA2 ,INTBUF_TAB,H3D_DATA,
37 . FNCONTP ,FTCONTP )
38C-----------------------------------------------
39C M o d u l e s
40C-----------------------------------------------
41 USE intbufdef_mod
42 USE h3d_mod
43C-----------------------------------------------
44C I m p l i c i t T y p e s
45C-----------------------------------------------
46#include "implicit_f.inc"
47C-----------------------------------------------
48C D u m m y A r g u m e n t s
49C-----------------------------------------------
50 INTEGER ILEV
51 INTEGER IPARI(*),WEIGHT(*),IGEO(NPROPGI,*),NPF(*),ITAB(*)
52C REAL
53 my_real
54 . ms(*),in(*), x(*), v(*), a(*),stifn(*),fsav(*),tf(*),
55 . fncont(*),pdama2(*), fncontp(*) ,ftcontp(* )
56
57 TYPE(intbuf_struct_) INTBUF_TAB
58 TYPE (H3D_DATABASE) :: H3D_DATA
59C-----------------------------------------------
60C C o m m o n B l o c k s
61C-----------------------------------------------
62#include "param_c.inc"
63C-----------------------------------------------
64C L o c a l V a r i a b l e s
65C-----------------------------------------------
66 integer
67 . nuvar,nsn,nmn,nrts,nrtm,igtyp,pid,isym
68C=======================================================================
69 nrts = ipari(3)
70 nrtm = ipari(4)
71 nsn = ipari(5)
72 nmn = ipari(6)
73 nuvar = ipari(35)
74 isym = ipari(44)
75c
76 IF (ilev == 10 .OR. ilev == 11 .OR. ilev == 12) THEN
77 pid = ipari(43)
78 IF(pid > 0) THEN ! on test test_i2_10_0001.rad , PID = 0
79 igtyp = igeo(11,pid)
80 ELSE
81 igtyp = -1
82 ENDIF
83 ELSE
84 pid = 0
85 igtyp =-1
86 ENDIF
87C
88 CALL i2rupt(
89 1 x ,v ,a ,ms ,in ,
90 2 stifn ,fsav ,weight ,intbuf_tab%IRECTM,intbuf_tab%NSV,
91 3 intbuf_tab%MSR ,intbuf_tab%IRTLM,intbuf_tab%IRUPT ,intbuf_tab%CSTS ,intbuf_tab%NMAS(1),
92 4 intbuf_tab%NMAS(1+nmn),intbuf_tab%SMAS ,intbuf_tab%SINER ,intbuf_tab%AREAS2,intbuf_tab%UVAR ,
93 5 intbuf_tab%XM0 ,intbuf_tab%DSM ,intbuf_tab%FSM ,intbuf_tab%RUPT ,ipari ,
94 6 nsn ,nmn ,nuvar ,igtyp ,pid ,
95 7 npf ,tf ,itab ,fncont ,pdama2 ,
96 8 isym ,intbuf_tab%INORM,h3d_data,fncontp ,ftcontp )
97C-----------
98 RETURN
99 END SUBROUTINE int2rupt
100C
101!||====================================================================
102!|| i2rupt ../engine/source/interfaces/interf/int2rupt.F
103!||--- called by ------------------------------------------------------
104!|| int2rupt ../engine/source/interfaces/interf/int2rupt.F
105!||--- calls -----------------------------------------------------
106!|| ancmsg ../engine/source/output/message/message.F
107!|| arret ../engine/source/system/arret.f
108!|| ruptint2 ../engine/source/interfaces/interf/ruptint2.F
109!||--- uses -----------------------------------------------------
110!|| h3d_mod ../engine/share/modules/h3d_mod.F
111!|| message_mod ../engine/share/message_module/message_mod.F
112!|| uintbuf_mod ../engine/source/user_interface/uintbuf_mod.F
113!||====================================================================
114 SUBROUTINE i2rupt(X ,V ,A ,MS ,IN ,
115 2 STIFN ,FSAV ,WEIGHT ,IRECT ,NSV ,
116 3 MSR ,IRTL ,IRUPT ,CRST ,MMASS ,
117 4 MINER ,SMASS ,SINER ,AREA ,UVAR ,
118 5 XSM0 ,DSM ,FSM ,PROP ,IPARI ,
119 6 NSN ,NMN ,NUVAR ,IGTYP ,PID ,
120 7 NPF ,TF ,ITAB ,FNCONT ,PDAMA2 ,
121 8 ISYM ,INORM ,H3D_DATA,FNCONTP ,FTCONTP )
122C-----------------------------------------------
123C M o d u l e s
124C-----------------------------------------------
125 USE uintbuf_mod
126 USE h3d_mod
127 USE message_mod
128C-----------------------------------------------
129C I m p l i c i t T y p e s
130C-----------------------------------------------
131#include "implicit_f.inc"
132C-----------------------------------------------
133C D u m m y A r g u m e n t s
134C-----------------------------------------------
135 INTEGER NSN, NMN,NUVAR,PID,IGTYP,ISYM
136 INTEGER IRECT(4,*),MSR(*),NSV(*),IRTL(*),IRUPT(*),IPARI(*),
137 . weight(*),npf(*),itab(*),pdama2(*),inorm(*)
138C REAL
139 my_real
140 . x(3,*),v(3,*),a(3,*),xsm0(3,*),dsm(3,*),fsm(3,*),prop(*),
141 . area(*),stifn(*),siner(*),smass(*),fncont(3,*),
142 . ms(*),in(*),mmass(*),miner(*),crst(2,*),fsav(*),tf(*),
143 . fncontp(3,*) ,ftcontp(3,* )
144 my_real, DIMENSION(NUVAR,NSN) ::
145 . uvar
146 TYPE (H3D_DATABASE) :: H3D_DATA
147C-----------------------------------------------
148C C o m m o n B l o c k s
149C-----------------------------------------------
150#include "com01_c.inc"
151#include "com08_c.inc"
152#include "scr14_c.inc"
153#include "comlock.inc"
154#include "userlib.inc"
155C-----------------------------------------------
156C L o c a l V a r i a b l e s
157C-----------------------------------------------
158 INTEGER NIR,I,J,II,JJ,L,N1,N2,N3,N4,RFLAG,W,IONE,
159 . IMOD,IFILTR,IFUNS,IFUNN,IFUNT,NOINT,IDBG
160C REAL
161 my_real
162 . S,T,SP,SM,TP,TM,AA,INS,MXI,MYI,MZI,XSM,YSM,ZSM,XC,YC,ZC,
163 . X0,X1,X2,X3,X4,Y0,Y1,Y2,Y3,Y4,Z0,Z1,Z2,Z3,Z4,SUM,DX,DY,DZ,
164 . FX,FY,FZ,FXI,FYI,FZI,SX,SY,SZ,TX,TY,TZ,DT12M,DTIME,
165 . DXNORM,DXTANG,DXTAN2,FNORM,FTANG,FTAN2,MCDG,
166 . VX,VY,VZ,VNX,VNY,VNZ,VTX,VTY,VTZ,VUX,VUY,VUZ,FACN,FACT,
167 . FSAV1,FSAV2,FSAV3,FSAV4,FSAV5,FSAV6,FSAV7,FSAV8,FSAV9,FSAV10,
168 . FSAV11,FSAV12,FSAV13,FSAV14,FSAV15,IMPX,IMPY,IMPZ
169 my_real
170 . h(4),dxn(3),dxt(3),dxu(3),fn(3),ft(3),fu(3),fnarea,ftarea,adxtang
171C----
172 type(uintbuf) :: USERBUF
173!
174 CHARACTER OPTION*256
175 INTEGER SIZE
176C=======================================================================
177 NIR =4
178 ione=1
179C
180 IF (tt == zero) THEN
181 DO ii=1,nsn
182 i=nsv(ii)
183 l=irtl(ii)
184 irupt(ii) = 0
185 inorm(ii) = 0
186C initial distance main-secnd
187 s = crst(1,ii)
188 t = crst(2,ii)
189 sp=one + s
190 sm=one - s
191 tp=fourth*(one + t)
192 tm=fourth*(one - t)
193 h(1)=tm*sm
194 h(2)=tm*sp
195 h(3)=tp*sp
196 h(4)=tp*sm
197 n1 = irect(1,l)
198 n2 = irect(2,l)
199 n3 = irect(3,l)
200 n4 = irect(4,l)
201 x0 = x(1,i)
202 y0 = x(2,i)
203 z0 = x(3,i)
204 x1 = x(1,n1)
205 y1 = x(2,n1)
206 z1 = x(3,n1)
207 x2 = x(1,n2)
208 y2 = x(2,n2)
209 z2 = x(3,n2)
210 x3 = x(1,n3)
211 y3 = x(2,n3)
212 z3 = x(3,n3)
213 x4 = x(1,n4)
214 y4 = x(2,n4)
215 z4 = x(3,n4)
216 xc = x1 * h(1) + x2 * h(2) + x3 * h(3) + x4 * h(4)
217 yc = y1 * h(1) + y2 * h(2) + y3 * h(3) + y4 * h(4)
218 zc = z1 * h(1) + z2 * h(2) + z3 * h(3) + z4 * h(4)
219 xsm = x0 - xc
220 ysm = y0 - yc
221 zsm = z0 - zc
222C
223 xsm0(1,ii) = xsm
224 xsm0(2,ii) = ysm
225 xsm0(3,ii) = zsm
226 dsm(1,ii) = zero
227 dsm(2,ii) = zero
228 dsm(3,ii) = zero
229 IF (isym == 1) THEN
230 sx = -(x2 + x3 - x1 - x4)
231 sy = -(y2 + y3 - y1 - y4)
232 sz = -(z2 + z3 - z1 - z4)
233 tx = -(x3 + x4 - x1 - x2)
234 ty = -(y3 + y4 - y1 - y2)
235 tz = -(z3 + z4 - z1 - z2)
236 vnx = sy * tz - sz * ty
237 vny = sz * tx - sx * tz
238 vnz = sx * ty - sy * tx
239 sum = one / sqrt(vnx*vnx + vny*vny + vnz*vnz)
240 vnx = vnx * sum
241 vny = vny * sum
242 vnz = vnz * sum
243 sum = vnx*xsm + vny*ysm + vnz*zsm
244 inorm(ii) = sign(ione, nint(sum))
245 ELSE
246 inorm(ii) = ione
247 ENDIF
248 ENDDO
249 dt12m = one/dt2
250 ELSE
251 dt12m = one/dt12
252 ENDIF
253C---------------
254C Reset initial main & secnd masse /inertia
255C---------------
256 DO ii=1,nmn
257 i=msr(ii)
258 ms(i) = mmass(ii)
259 IF (iroddl /= 0) in(i) = miner(ii)
260 ENDDO
261 DO ii=1,nsn
262 i=nsv(ii)
263 ms(i) = smass(ii)
264 IF (iroddl /= 0) in(i)= siner(ii)
265 ENDDO
266C
267C------------------------------------------------
268 DO ii=1,nsn
269 i=nsv(ii)
270 l=irtl(ii)
271 rflag = irupt(ii)
272C
273 IF (i > 0) THEN
274 IF (irupt(ii) == 1) THEN
275c rupture totale
276 fn(1) = zero
277 fn(2) = zero
278 fn(3) = zero
279 ft(1) = zero
280 ft(2) = zero
281 ft(3) = zero
282 fu(1) = zero
283 fu(2) = zero
284 fu(3) = zero
285 fnorm = zero
286 ftang = zero
287 ELSE
288 s = crst(1,ii)
289 t = crst(2,ii)
290 sp=one + s
291 sm=one - s
292 tp=fourth*(one + t)
293 tm=fourth*(one - t)
294 h(1)=tm*sm
295 h(2)=tm*sp
296 h(3)=tp*sp
297 h(4)=tp*sm
298C
299 n1 = irect(1,l)
300 n2 = irect(2,l)
301 n3 = irect(3,l)
302 n4 = irect(4,l)
303C------------------------------------------------
304C relative deplacement main-secnd
305C------------------------------------------------
306 x0 = x(1,i)
307 y0 = x(2,i)
308 z0 = x(3,i)
309 x1 = x(1,n1)
310 y1 = x(2,n1)
311 z1 = x(3,n1)
312 x2 = x(1,n2)
313 y2 = x(2,n2)
314 z2 = x(3,n2)
315 x3 = x(1,n3)
316 y3 = x(2,n3)
317 z3 = x(3,n3)
318 x4 = x(1,n4)
319 y4 = x(2,n4)
320 z4 = x(3,n4)
321 xc = x1 * h(1) + x2 * h(2) + x3 * h(3) + x4 * h(4)
322 yc = y1 * h(1) + y2 * h(2) + y3 * h(3) + y4 * h(4)
323 zc = z1 * h(1) + z2 * h(2) + z3 * h(3) + z4 * h(4)
324 xsm = x0 - xc
325 ysm = y0 - yc
326 zsm = z0 - zc
327C--------- depl relatif total
328 dx = xsm - xsm0(1,ii)
329 dy = ysm - xsm0(2,ii)
330 dz = zsm - xsm0(3,ii)
331 dsm(1,ii) = dx
332 dsm(2,ii) = dy
333 dsm(3,ii) = dz
334C------------------------------------------------
335C direction normale facette main
336C------------------------------------------------
337 sx = -(x2 + x3 - x1 - x4)
338 sy = -(y2 + y3 - y1 - y4)
339 sz = -(z2 + z3 - z1 - z4)
340 tx = -(x3 + x4 - x1 - x2)
341 ty = -(y3 + y4 - y1 - y2)
342 tz = -(z3 + z4 - z1 - z2)
343 vnx = sy * tz - sz * ty
344 vny = sz * tx - sx * tz
345 vnz = sx * ty - sy * tx
346 sum = inorm(ii) / sqrt(vnx*vnx + vny*vny + vnz*vnz)
347 vnx = vnx * sum
348 vny = vny * sum
349 vnz = vnz * sum
350 sum = one / sqrt(sx*sx + sy*sy + sz*sz)
351 vtx = sx * sum
352 vty = sy * sum
353 vtz = sz * sum
354C------------------------------------------------
355C composantes N/T du deplacement
356C------------------------------------------------
357 dxnorm = vnx*dx + vny*dy + vnz*dz
358 dxn(1) = vnx*dxnorm
359 dxn(2) = vny*dxnorm
360 dxn(3) = vnz*dxnorm
361 dxt(1) = dx - dxn(1)
362 dxt(2) = dy - dxn(2)
363 dxt(3) = dz - dxn(3)
364 dxtang = sqrt(dxt(1)**2 + dxt(2)**2 + dxt(3)**2)
365C repere local
366 IF (dxtang > zero) THEN
367 sum = one / dxtang
368 vtx = dxt(1) * sum
369 vty = dxt(2) * sum
370 vtz = dxt(3) * sum
371 ENDIF
372 vux = vny * vtz - vnz * vty
373 vuy = vnz * vtx - vnx * vtz
374 vuz = vnx * vty - vny * vtx
375C------------------------------------------------
376C force main/secnd
377C------------------------------------------------
378 IF (irupt(ii) == 0) THEN
379C pas de rupture - interface cinematique
380 mcdg = one/(ms(i)+ms(n1)+ms(n2)+ms(n3)+ms(n4))
381 vx = v(1,i)*ms(i)+v(1,n1)*ms(n1)
382 . +v(1,n2)*ms(n2)+v(1,n3)*ms(n3)+v(1,n4)*ms(n4)
383 vy = v(2,i)*ms(i)+v(2,n1)*ms(n1)
384 . +v(2,n2)*ms(n2)+v(2,n3)*ms(n3)+v(2,n4)*ms(n4)
385 vz = v(3,i)*ms(i)+v(3,n1)*ms(n1)
386 . +v(3,n2)*ms(n2)+v(3,n3)*ms(n3)+v(3,n4)*ms(n4)
387 vx = vx * mcdg
388 vy = vy * mcdg
389 vz = vz * mcdg
390 fx = a(1,i) + (v(1,i) - vx)*ms(i)*dt12m
391 fy = a(2,i) + (v(2,i) - vy)*ms(i)*dt12m
392 fz = a(3,i) + (v(3,i) - vz)*ms(i)*dt12m
393 ELSEIF (irupt(ii) == -1) THEN
394C rupture partielle = spring main/secnd
395 fx = stifn(i)*dx
396 fy = stifn(i)*dy
397 fz = stifn(i)*dz
398 ENDIF
399C composantes N/T de la forces nodale
400 fnorm = vnx*fx + vny*fy + vnz*fz
401 ftang = vtx*fx + vty*fy + vtz*fz
402 ftan2 = vux*fx + vuy*fy + vuz*fz
403 fn(1) = vnx*fnorm
404 fn(2) = vny*fnorm
405 fn(3) = vnz*fnorm
406 ft(1) = vtx*ftang
407 ft(2) = vty*ftang
408 ft(3) = vtz*ftang
409 fu(1) = vux*ftan2
410 fu(2) = vuy*ftan2
411 fu(3) = vuz*ftan2
412C------------------------------------------------
413C Call user subroutine - test rupture
414C------------------------------------------------
415 dtime = max(dt1,em20)
416 IF(userl_avail==1.AND.igtyp /= -1)THEN
417 adxtang = abs(dxtang)
418 fnarea=fnorm / area(ii)
419 ftarea=max(abs(ftang),abs(ftan2)) / area(ii)
420 CALL eng_userlib_uintbuf_var(i,area(ii),dxnorm,adxtang,fnarea,ftarea,dtime,rflag)
421
422 ELSE
423 userbuf%ISECND = i
424 userbuf%AREA = area(ii)
425 userbuf%DXN = dxnorm
426 userbuf%DXT = abs(dxtang)
427 userbuf%SIGN = fnorm / area(ii)
428 userbuf%SIGT = max(abs(ftang),abs(ftan2)) / area(ii)
429 userbuf%DT = dtime
430 userbuf%RUPT = rflag
431 ENDIF
432C
433 IF (igtyp == -1) THEN
434 noint = ipari(15)
435 imod = ipari(43)
436 ifiltr= ipari(59)
437 ifuns = ipari(48)
438 ifunn = ipari(49)
439 ifunt = ipari(50)
440 idbg = 0
441 CALL ruptint2(
442 . nsn ,ii ,nuvar ,uvar(1,ii),userbuf ,
443 . prop ,ifuns ,ifunn ,ifunt ,imod ,
444 . ifiltr ,idbg ,npf ,tf ,noint ,
445 . itab ,pdama2 ,isym ,h3d_data )
446 ELSEIF (igtyp == 29) THEN
447 IF(userl_avail==1)THEN
448 CALL eng_userlib_userint(igtyp,
449 . nsn ,ii ,pid ,nuvar ,
450 . uvar(1,ii),userbuf )
451 ELSE
452 ! ----------------
453 ! ERROR to be printed & exit
454 option='INTERFACE TYPE2 RUPTURE MODEL'
455 size=len_trim(option)
456 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
457 CALL arret(2)
458 ! ----------------
459 ENDIF
460
461 ELSEIF (igtyp == 30) THEN
462 IF(userl_avail==1)THEN
463 CALL eng_userlib_userint(igtyp,
464 . nsn ,ii ,pid ,nuvar ,
465 . uvar(1,ii),userbuf )
466 ELSE
467 ! ----------------
468 ! ERROR to be printed & exit
469 option='INTERFACE TYPE2 RUPTURE MODEL'
470 size=len_trim(option)
471 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
472 CALL arret(2)
473 ! ----------------
474 ENDIF
475 ELSEIF (igtyp == 31) THEN
476 IF(userl_avail==1)THEN
477 CALL eng_userlib_userint(igtyp,
478 . nsn ,ii ,pid ,nuvar ,
479 . uvar(1,ii),userbuf )
480 ELSE
481 ! ----------------
482 ! ERROR to be printed & exit
483 option='INTERFACE TYPE2 RUPTURE MODEL'
484 size=len_trim(option)
485 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
486 CALL arret(2)
487 ! ----------------
488 ENDIF
489 ELSEIF (igtyp == 37) THEN
490 IF(userl_avail==1)THEN
491 CALL eng_userlib_userint(igtyp,
492 . nsn ,ii ,pid ,nuvar ,
493 . uvar(1,ii),userbuf )
494 ELSE
495 ! ----------------
496 ! ERROR to be printed & exit
497 option='INTERFACE TYPE2 RUPTURE MODEL'
498 size=len_trim(option)
499 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
500 CALL arret(2)
501 ! ----------------
502 ENDIF
503 ELSEIF (igtyp == 38) THEN
504 IF(userl_avail==1)THEN
505 CALL eng_userlib_userint(igtyp,
506 . nsn ,ii ,pid ,nuvar ,
507 . uvar(1,ii),userbuf )
508 ELSE
509 ! ----------------
510 ! ERROR to be printed & exit
511 option='INTERFACE TYPE2 RUPTURE MODEL'
512 size=len_trim(option)
513 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
514 CALL arret(2)
515 ! ----------------
516 ENDIF
517 ELSEIF (igtyp == 39) THEN
518 IF(userl_avail==1)THEN
519 CALL eng_userlib_userint(igtyp,
520 . nsn ,ii ,pid ,nuvar ,
521 . uvar(1,ii),userbuf )
522 ELSE
523 ! ----------------
524 ! ERROR to be printed & exit
525 option='INTERFACE TYPE2 RUPTURE MODEL'
526 size=len_trim(option)
527 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
528 CALL arret(2)
529 ! ----------------
530 ENDIF
531 ELSEIF (igtyp == 40) THEN
532 IF(userl_avail==1)THEN
533 CALL eng_userlib_userint(igtyp,
534 . nsn ,ii ,pid ,nuvar ,
535 . uvar(1,ii),userbuf )
536 ELSE
537 ! ----------------
538 ! ERROR to be printed & exit
539 option='INTERFACE TYPE2 RUPTURE MODEL'
540 size=len_trim(option)
541 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
542 CALL arret(2)
543 ! ----------------
544 ENDIF
545 ELSEIF (igtyp == 41) THEN
546 IF(userl_avail==1)THEN
547 CALL eng_userlib_userint(igtyp,
548 . nsn ,ii ,pid ,nuvar ,
549 . uvar(1,ii),userbuf )
550 ELSE
551 ! ----------------
552 ! ERROR to be printed & exit
553 option='INTERFACE TYPE2 RUPTURE MODEL'
554 size=len_trim(option)
555 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
556 CALL arret(2)
557 ! ----------------
558 ENDIF
559 ELSEIF (igtyp == 42) THEN
560 IF(userl_avail==1)THEN
561 CALL eng_userlib_userint(igtyp,
562 . nsn ,ii ,pid ,nuvar ,
563 . uvar(1,ii),userbuf )
564 ELSE
565 ! ----------------
566 ! ERROR to be printed & exit
567 option='INTERFACE TYPE2 RUPTURE MODEL'
568 size=len_trim(option)
569 CALL ancmsg(msgid=257,c1=option(1:size),anmode=aninfo)
570 CALL arret(2)
571 ! ----------------
572 ENDIF
573 ENDIF
574C
575 facn = userbuf%FACN
576 fact = userbuf%FACT
577 irupt(ii) = nint(userbuf%RUPT)
578 fn(1) = fn(1) * facn
579 fn(2) = fn(2) * facn
580 fn(3) = fn(3) * facn
581 ft(1) = ft(1) * fact
582 ft(2) = ft(2) * fact
583 ft(3) = ft(3) * fact
584 fu(1) = fu(1) * fact
585 fu(2) = fu(2) * fact
586 fu(3) = fu(3) * fact
587 fnorm = fnorm*facn
588 ftang = ftang*fact
589 ftan2 = ftan2*fact
590 fsm(1,ii) = vnx*fnorm + vtx*ftang + vux*ftan2
591 fsm(2,ii) = vny*fnorm + vty*ftang + vuy*ftan2
592 fsm(3,ii) = vnz*fnorm + vtz*ftang + vuz*ftan2
593C
594 ENDIF
595C------------------------------------------------
596C Save TH
597C------------------------------------------------
598 w = weight(i)
599 impx=fn(1)
600 impy=fn(2)
601 impz=fn(3)
602 fsav1 = impx*dt12
603 fsav2 = impy*dt12
604 fsav3 = impz*dt12
605 fsav8 = abs(impx)
606 fsav9 = abs(impy)
607 fsav10= abs(impz)
608 fsav11= fnorm
609 impx=ft(1)
610 impy=ft(2)
611 impz=ft(3)
612 fsav4 = impx*dt12
613 fsav5 = impy*dt12
614 fsav6 = impz*dt12
615 fsav12= abs(impx)
616 fsav13= abs(impy)
617 fsav14= abs(impz)
618 fsav15= ftang
619#include "lockon.inc"
620 fsav(1) = fsav(1) + fsav1*w
621 fsav(2) = fsav(2) + fsav2*w
622 fsav(3) = fsav(3) + fsav3*w
623 fsav(4) = fsav(4) + fsav4*w
624 fsav(5) = fsav(5) + fsav5*w
625 fsav(6) = fsav(6) + fsav6*w
626 fsav(8) = fsav(8) + fsav8*w
627 fsav(9) = fsav(9) + fsav9*w
628 fsav(10)= fsav(10)+ fsav10*w
629 fsav(11)= fsav(11)+ fsav11*w
630 fsav(12)= fsav(12)+ fsav12*w
631 fsav(13)= fsav(13)+ fsav13*w
632 fsav(14)= fsav(14)+ fsav14*w
633 fsav(15)= fsav(15)+ fsav15*w
634#include "lockoff.inc"
635C---
636 IF(anim_v(13)+h3d_data%N_VECT_CONT2>0) THEN
637 fncont(1,i) = - (fn(1)+ft(1)) * w
638 fncont(2,i) = - (fn(2)+ft(2)) * w
639 fncont(3,i) = - (fn(3)+ft(3)) * w
640 DO jj=1,nir
641 j=irect(jj,l)
642 fncont(1,j) = fncont(1,j) + w *(fn(1)+ft(1))*h(jj)
643 fncont(2,j) = fncont(2,j) + w *(fn(2)+ft(2))*h(jj)
644 fncont(3,j) = fncont(3,j) + w *(fn(3)+ft(3))*h(jj)
645 ENDDO
646 ENDIF
647
648 IF(anim_v(27)+h3d_data%N_VECT_PCONT2>0) THEN ! Normal/Tangential forces output
649 fncontp(1,i) = - (fn(1)+ft(1)) * w
650 fncontp(2,i) = - (fn(2)+ft(2)) * w
651 fncontp(3,i) = - (fn(3)+ft(3)) * w
652 DO jj=1,nir
653 j=irect(jj,l)
654 fncontp(1,j) = fncontp(1,j) - fncontp(1,i)*h(jj)
655 fncontp(2,j) = fncontp(2,j) - fncontp(2,i)*h(jj)
656 fncontp(3,j) = fncontp(3,j) - fncontp(3,i)*h(jj)
657 ENDDO
658
659 ftcontp(1,i) = vnx * w
660 ftcontp(2,i) = vny * w
661 ftcontp(3,i) = vnz * w
662 DO jj=1,nir
663 j=irect(jj,l)
664 ftcontp(1,j) = ftcontp(1,j) - ftcontp(1,i)*h(jj)
665 ftcontp(2,j) = ftcontp(2,j) - ftcontp(2,i)*h(jj)
666 ftcontp(3,j) = ftcontp(3,j) - ftcontp(3,i)*h(jj)
667 ENDDO
668 ENDIF
669 ENDIF
670 ENDDO
671C-----------
672 RETURN
673 END SUBROUTINE i2rupt
my_real
#define my_real
Definition cppsort.cpp:32
area
subroutine area(d1, x, x2, y, y2, eint, stif0)
Definition hm_read_prop32.F:567
int2rupt
subroutine int2rupt(ipari, ms, in, x, v, a, stifn, igeo, weight, fsav, ilev, npf, tf, itab, fncont, pdama2, intbuf_tab, h3d_data, fncontp, ftcontp)
Definition int2rupt.F:38
i2rupt
subroutine i2rupt(x, v, a, ms, in, stifn, fsav, weight, irect, nsv, msr, irtl, irupt, crst, mmass, miner, smass, siner, area, uvar, xsm0, dsm, fsm, prop, ipari, nsn, nmn, nuvar, igtyp, pid, npf, tf, itab, fncont, pdama2, isym, inorm, h3d_data, fncontp, ftcontp)
Definition int2rupt.F:122
max
#define max(a, b)
Definition macros.h:21
h3d_mod
Definition h3d_mod.F:308
message_mod
Definition message_mod.F:1249
uintbuf_mod
Definition uintbuf_mod.F:29
ruptint2
subroutine ruptint2(nsl, isl, nuvar, uvar, userbuf, prop, ifuns, ifunn, ifunt, imod, ifiltr, idbg, npf, tf, noint, itab, pdama2, isym, h3d_data)
Definition ruptint2.F:38
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889
arret
subroutine arret(nn)
Definition arret.F:87
uintbuf_mod::uintbuf
Definition uintbuf_mod.F:35