- Page 2 and 3:
The mission of the Institute for th
- Page 4 and 5:
The present document contains the p
- Page 6 and 7:
6 of 556
- Page 8:
8 LOCALWWW.CONTEXTODEDURANGO.COM.MX
- Page 11 and 12:
Introductory Example (2)5Introducto
- Page 13 and 14:
Goals/Characteristics• Simulate f
- Page 15 and 16:
Mu f extB Tσ dVe eV= −∑ ∫Com
- Page 17:
Direct Time Integration (4)• The
- Page 21 and 22:
Geometric non-linearitiesA large-di
- Page 23 and 24:
Advantages of the Method• Transie
- Page 25 and 26:
Finding the Lagrange Multipliers (2
- Page 28:
Exercise 0 - Ideal ballistics (5)
- Page 32 and 33:
Exercise 1 - Suspended mass (7)•
- Page 34 and 35:
Exercise 2 - Wave propagation (4)
- Page 36 and 37:
Exercise 3 - Impact on Cooling Towe
- Page 38 and 39:
38 of 556
- Page 40 and 41:
40 of 556
- Page 42 and 43:
RESULTS:Results are in perfect agre
- Page 44 and 45:
VIEW 0.00000E+00 0.00000E+00 -1.000
- Page 46 and 47:
Trajectories:These results are in p
- Page 48 and 49:
848 of 556
- Page 50 and 51:
Hence:11 −8ES 2× 10 ⋅ 2.5×10
- Page 52 and 53:
The stress history is:Note that the
- Page 54 and 55:
AnalyticalTEST11Same as TEST01 but
- Page 56 and 57:
Comparison of natural vs. engineeri
- Page 58 and 59:
The total external forces at the tw
- Page 60 and 61:
1260 of 556
- Page 62 and 63:
TEST14Same as TEST13 but uses samee
- Page 64 and 65:
The resulting displacement is shown
- Page 66 and 67:
The load is reversed when the cable
- Page 68 and 69:
The (dynamic) equilibrium is expres
- Page 70 and 71:
Simple 1-D wave propagation problem
- Page 72 and 73:
Numerical SolutionsBARI01We discret
- Page 74 and 75:
BARI02We introduce some damping (10
- Page 76 and 77:
BARI09Use a 2-D geometry (elements
- Page 78 and 79:
or the velocities:or the Von Mises
- Page 80 and 81:
1480 of 556
- Page 82 and 83:
The system is initially at rest, an
- Page 84 and 85:
CAST MESHTRID NONL LAGR$$ Dimension
- Page 86 and 87:
Final deformation (with superposed
- Page 88 and 89:
88 of 556
- Page 90 and 91:
Detailed Contents• Modeling the f
- Page 92 and 93:
Modeling the fluid domain• The fl
- Page 94 and 95:
Euler equations (3)• For a compre
- Page 96 and 97:
Time integrationEach time increment
- Page 98 and 99:
Time integration (5)10. Account for
- Page 100 and 101:
Euler Equations (FV)• They are wr
- Page 102 and 103:
Mesh rezoning - motivations• Rezo
- Page 104 and 105:
Giuliani’s automatic rezoning•
- Page 106 and 107:
• Analytical solution (self-simil
- Page 108 and 109:
Shock tube (6)• Influence of pseu
- Page 110 and 111:
Exercise 2 - Explosion inair-filled
- Page 112 and 113:
Exercise 3 - Bubble expansionin a T
- Page 114 and 115:
Exercise 4 - External blast on two
- Page 116 and 117:
116 of 556
- Page 118 and 119:
118 of 556
- Page 120 and 121:
Analytical solution (see e.g. Harlo
- Page 122 and 123:
3.05000E+01 5.00000E-01 3.10000E+01
- Page 124 and 125:
AnalyticalConclusion: the pseudo-vi
- Page 126 and 127:
Finally, we investigate the effect
- Page 128 and 129:
“Optimal” solution (for the pre
- Page 130 and 131:
Initial conditions: FE vs. FVThe tw
- Page 132 and 133:
Example of shock tube problem solve
- Page 134 and 135:
The pressure distribution at t = 0.
- Page 136 and 137:
TANK01Lagrangian solution: the whol
- Page 138 and 139:
The computed fluid pressures at the
- Page 140 and 141:
The computed pressure histories are
- Page 142 and 143:
8142 of 556
- Page 144 and 145:
TUBE06Lagrangian solution: the whol
- Page 146 and 147:
And the final velocity distribution
- Page 148 and 149:
The computed displacements are:To s
- Page 150 and 151:
43 30 30 43 44 31 31 44 45 3232 45
- Page 152 and 153:
10152 of 556
- Page 154 and 155:
p6 = 70 20 0;p7 = 70 40 0;p8 = 40 4
- Page 156 and 157:
$ high-pressure perfect gas (explos
- Page 158 and 159:
*r_cais = 965.00e-3 ;h_cais = 1560.
- Page 160 and 161:
flut ro 1.3 eint 0.21978e6 gamm 1.3
- Page 162 and 163:
The fluid pressures and velocities
- Page 164 and 165:
mess 'NB_POIN =' (NBNO tout) ;mess
- Page 166 and 167:
The fluid pressures and velocities
- Page 168 and 169:
168 of 556
- Page 170 and 171:
Further FSI ExampleElectric arc in
- Page 172 and 173:
FSI Motivation (2)Fully coupled ana
- Page 174 and 175:
Permanent FSI treatment (2)• Upon
- Page 176 and 177:
Geometrically complex cases• Bila
- Page 178 and 179:
The FSA algorithm (2)• Effect of
- Page 180 and 181:
The UP algorithmThe method simply r
- Page 182 and 183:
Application to Finite Volumes• Th
- Page 184 and 185:
Exercise 2 - Explosions insimple de
- Page 186 and 187:
Exercise/Example 3 - Wavepropagatio
- Page 188 and 189:
Exercise/Example 4 - CONT problemTh
- Page 190 and 191:
Exercise/Example 4 - CONT problem (
- Page 192 and 193:
Exercise/Example 6 : Woodward-Colel
- Page 194 and 195:
Exercise/Example 7 : Exfs Test(Comp
- Page 196 and 197:
Exercise/Example 8 : Steam Explosio
- Page 198 and 199:
Exercise/Example 9Explosion in Seco
- Page 200 and 201:
Boundary Condition Elements: CLxx
- Page 202 and 203:
Boundary Condition Elements: CLxx (
- Page 204 and 205:
Exercise/Example 12 - Sloshing (Cou
- Page 206 and 207:
Spectrum at point P5Exercise/Exampl
- Page 208 and 209:
Exercise/Example 12 - Sloshing (9)
- Page 210 and 211:
Exercise/Example 13 - Building Vuln
- Page 212 and 213:
Exercise/Example 13 - Building Vuln
- Page 214 and 215:
Exercise/Example 13 - Building Vuln
- Page 216 and 217:
Exercise 14 - Improving Initial Con
- Page 218 and 219:
218 of 556
- Page 220 and 221:
DIMEPT3L 23 PT2L 143 FL24 145 ED01
- Page 222 and 223:
The final pressure field in the flu
- Page 224 and 225:
ALE solution with FSA: it is not po
- Page 226 and 227:
120 119 138 139121 120 139 140122 1
- Page 228 and 229:
and the final velocity field:The su
- Page 230 and 231:
The final velocities:The final liqu
- Page 232 and 233:
*----------------------------------
- Page 234 and 235:
TWIS07We study the phenomenon by a
- Page 236 and 237:
6236 of 556
- Page 238 and 239:
DIMEPT2L 293 PT3L 70 ZONE 3ED41 33
- Page 240 and 241:
2.10000E+01 8.33330E+00 2.10000E+01
- Page 242 and 243:
FLUT RO 2.4278E3 EINT 0. GAMM 0.75D
- Page 244 and 245:
The final mesh (colors indicate mat
- Page 246 and 247:
GRIL LAGR LECT stru TERMALE LECT fl
- Page 248 and 249:
The final bubble material mass frac
- Page 250 and 251:
Geometric data:MaterialsThe explosi
- Page 252 and 253:
Some results: global deformed mesh
- Page 254 and 255:
Detail of first diaphragm:Detail of
- Page 256 and 257:
BOUNDARY CONDITIONS:The step is ent
- Page 258 and 259:
The EUROPLEXUS input file reads:WOC
- Page 260 and 261:
WOCO3DThe mesh generation file (K20
- Page 262 and 263:
*RESU ALIC GARD PSCR*SORT GRAP*AXTE
- Page 264 and 265:
PT3L 16389PT6L 504 NBLE 16389NALE 1
- Page 266 and 267:
Final pressure distribution:4266 of
- Page 268 and 269:
BOUNDARY CONDITIONS:The vessel is e
- Page 270 and 271:
as01=daller c1 c2 c3 c4 plan;*c1=p1
- Page 272 and 273:
*trac oeil cach fluidstr;*opti donn
- Page 274 and 275:
cam3=cam31 et cam32;camb=cam1 et ca
- Page 276 and 277:
log 1dtmlREZO GAM0 0.5FLMP EPS1 5.0
- Page 278 and 279:
Final structure deformation:Final v
- Page 280 and 281:
concrete, and by the VMSF (von Mise
- Page 282 and 283:
flut=flui elem 'TRI3';fluq=flui ele
- Page 284 and 285:
COMPEPAI 1.15 LECT STR1 TERM0.61 LE
- Page 286 and 287:
Intermediate structure velocities:F
- Page 288 and 289:
10288 of 556
- Page 290 and 291:
LOADING:The event is initiated by t
- Page 292 and 293:
FLUT RO .242777373 EINT 6.865E5 GAM
- Page 294 and 295:
Final structure velocities:6294 of
- Page 296 and 297:
RESULTS:A paper by Bermudez and Rod
- Page 298 and 299:
The applied displacement and the co
- Page 300 and 301:
PROBLEM:A rigid tank with a deforma
- Page 302 and 303:
The EUROPLEXUS input file reads:GRA
- Page 304 and 305:
The velocity at 200 ms is presented
- Page 306 and 307:
compute a reasonable initial distri
- Page 308 and 309:
FSA LECT 2 3 4 5 6 7 8 9 10 1112 13
- Page 310 and 311:
$$$91 92 93 94 95 96 97 98 99 100 T
- Page 312 and 313:
The computed displacements of the t
- Page 314 and 315:
TITLE:OILCUP2: uniformly accelerate
- Page 316 and 317:
2526 2527 2528 2529 2530 2531 25322
- Page 318 and 319:
Intermediate and final oil mass fra
- Page 320 and 321:
The QUAS STAT option is used to int
- Page 322 and 323:
t8 = (r4*c45) (r4*c45) h2;*Group a:
- Page 324 and 325:
t7t8 = cerc t7 (0. 0. h2) t8 dini l
- Page 326 and 327:
fill0 = fil7 et fil8 et fil9 et fil
- Page 328 and 329:
d13d3 = u13u3 proj coni a0 plan d0
- Page 330 and 331:
icolo0 et ibwal1 et iswal0 et ichoi
- Page 332 and 333:
sph0 = sph0 et (sph0 syme poin ec0)
- Page 334 and 335:
Some (qualitative) results are show
- Page 336 and 337:
Some selected snapshots:t = 2 msT =
- Page 338 and 339:
CALCULATIONS:Preliminary calculatio
- Page 340 and 341:
QUAS02Further preliminary calculati
- Page 342 and 343:
s1 = p2p d 5 p3p;s2 = p3p d 10 p4p;
- Page 344 and 345:
The velocities and pressures during
- Page 346 and 347:
346 of 556
- Page 348 and 349:
Detailed contents• ALE descriptio
- Page 350 and 351:
Example 1 - Taylor bar impact7Examp
- Page 352 and 353:
Example 3 - CoiningBilateral,Plane
- Page 354 and 355:
Example 3 - Coining (5)Bilateral, 2
- Page 356 and 357:
Example 4 - Explosion in a 2D box19
- Page 358 and 359:
Example 5 - Explosion in a corridor
- Page 360 and 361:
Lagrangian contact• Non-permanent
- Page 362 and 363:
Example 6a - Deep drawing A (2)Velo
- Page 364 and 365:
Tube Crash (2)35Conventional contac
- Page 366 and 367:
The Basic Pinball Method[Belytschko
- Page 368 and 369:
Hierarchic Pinball Method (2)• Hi
- Page 370 and 371:
Example 7 - Cable impact (2)lLv 0MW
- Page 372 and 373:
Example 7b - Indentation (4)• 3D
- Page 374 and 375:
Example 7b - Indentation (8)• Com
- Page 376 and 377:
Example 7d - Falling Rock Catcher
- Page 378 and 379:
Smoothed Particles Hydrodynamics(Co
- Page 380 and 381:
SPH Formulation (2)• Basic idea:
- Page 382 and 383: SPH impact [Courtesy of Samtech/Son
- Page 384 and 385: PRGL01:Example 8 - SPH impacts(Cour
- Page 386 and 387: Spectral Elements• Motivation: lo
- Page 388 and 389: Spectral Elements (5)Convergence pr
- Page 390 and 391: Example 9 - Closed FE/SE interface
- Page 392 and 393: Example 11 - Cylindrical valley91Cy
- Page 394 and 395: Performance Optimization• All exa
- Page 396 and 397: Time Step Partitioning (4)• Intri
- Page 398 and 399: Time Step Partitioning (8)• Activ
- Page 400 and 401: Treatment of links in partitioning
- Page 402 and 403: Example 12a - Taylor bar impact rev
- Page 404 and 405: Coupling at the Interfaces• Two a
- Page 406 and 407: Coupling at theInterfaces (5)MUT−
- Page 408 and 409: Coupling at the Interfaces (9)Time
- Page 410 and 411: Multiple scales in time (3)• Expr
- Page 412 and 413: Multiple scales in space• Further
- Page 414 and 415: Multiple scales in frequency• Som
- Page 416 and 417: Example: power plant139Example: pow
- Page 418 and 419: Example: aircraft (3)Thanks toCPU g
- Page 420 and 421: Example 14 - Domains in 3D• Thick
- Page 422 and 423: Example 15 - Bar Impact Revisited (
- Page 424 and 425: 424 of 556
- Page 426 and 427: *mesh=stru et symax et viti et lili
- Page 428 and 429: The displacements are:4428 of 556
- Page 430 and 431: c2=p1 d 8 p2;c3=p2 d 20 p3;c4=p3 d
- Page 434 and 435: TERM*------------------------------
- Page 436 and 437: 8436 of 556
- Page 438 and 439: p4s=p4 plus (0 0);tol=0.001;c1=p1 d
- Page 440 and 441: The deformed final mesh at 4 ms wit
- Page 442 and 443: The input file:INFS - 02*----------
- Page 444 and 445: The final mesh and pressures are:Th
- Page 446 and 447: ZONE 2PMAT 2Q4GS 3904ECRO 858884BLO
- Page 448 and 449: The final deformed mesh is:The fina
- Page 450 and 451: 6450 of 556
- Page 452 and 453: *----------------------------------
- Page 454 and 455: The final deformed mesh is:The fina
- Page 456 and 457: pc = 4.5 -1 -1.5;pd = 5.5 -1 -1.5;a
- Page 458 and 459: The initial configuration (with par
- Page 460 and 461: 6460 of 556
- Page 462 and 463: LOADING:The indenter is pushed into
- Page 464 and 465: The initial configuration (with par
- Page 466 and 467: The reaction force is:Approximatean
- Page 468 and 469: The final velocities:8468 of 556
- Page 470 and 471: The reaction force is:Approximatean
- Page 472 and 473: The initial configuration (with par
- Page 474 and 475: The reaction force is:Approximatean
- Page 476 and 477: is implicitly contained within the
- Page 478 and 479: The input file is:plaque poinçon p
- Page 480 and 481: The configurations at 0.5 ms (half-
- Page 482 and 483:
The first ruptured elements appear
- Page 484 and 485:
The same results, more detailed, ar
- Page 486 and 487:
The final view is:12486 of 556
- Page 488 and 489:
The two supports are clamped at the
- Page 490 and 491:
sauv form mesh;*fin;The input file
- Page 492 and 493:
UP -2.75643E-01 3.29707E-01 9.02948
- Page 494 and 495:
poin sphp!pinb cdescolo papelima on
- Page 496 and 497:
The final configuration (geometry o
- Page 498 and 499:
NETCF8Solution with a cables ruptur
- Page 500 and 501:
The strain in the cables at 53 ms i
- Page 502 and 503:
-5.10000E-01 1.33333E+00 2.00000E+0
- Page 504 and 505:
Some hardening quantity in the targ
- Page 506 and 507:
Final plastic streen in the target
- Page 508 and 509:
Final plastic strain in the target:
- Page 510 and 511:
10510 of 556
- Page 512 and 513:
opti echo 0;opti donn 'D:\Users\Fol
- Page 514 and 515:
*opti dime 2 elem qua4;p0 = 0 0;p1
- Page 516 and 517:
The final X-displacement in the hyb
- Page 518 and 519:
*opti dime 2 elem qua4;opti titr 'V
- Page 520 and 521:
VALLPSThis calculation assumes a pu
- Page 522 and 523:
The receiver signal in the coupled
- Page 524 and 525:
8524 of 556
- Page 526 and 527:
VEC1=2. 0.;*P1=P11;R1=P13;S1=P13 PL
- Page 528 and 529:
S3=P33 PLUS VEC1;*N1=2;N2=3;N3=4;*P
- Page 530 and 531:
CHPO2=CHPO2 / SCAL1;TSTAT2 . &BCL9=
- Page 532 and 533:
AXTE 1E3 'Temps (ms)'*-------------
- Page 534 and 535:
The final displacement field in the
- Page 536 and 537:
P9=40. 5. 5.;P14=20. 5. 5.;*BASE =
- Page 538 and 539:
CLIM1=(P5 ET P6 ET P7 ET P8);CLIM1=
- Page 540 and 541:
SUITPost-treatmentECHORESU ALIC TEM
- Page 542 and 543:
The final displacement field in the
- Page 544 and 545:
tol=0.01E-3;*base=p0 d 5 p1;stru=ba
- Page 546 and 547:
eli=stru elem i;si (ega j 0);sq42=e
- Page 548 and 549:
h3=8.1e-3;h4=16.2e-3;h5=24.3e-3;*n1
- Page 550 and 551:
The characteristic displacements in
- Page 552 and 553:
The final yield stress field in the
- Page 554 and 555:
554 of 556
- Page 556:
The mission of the Joint Research C