Download full text - ELSA - Europa

More documents

Recommendations

Info

Universitat Politècnica de Catalunya, Barcelona, 15-18 May 2006 Numerical Simulation of Fast Transient Phenomena in Fluid-Structure Systems A Short Course by F. Casadei European Commission, Joint Research Centre Institute for the Protection and Security of the Citizen T.P. 480, I-21020 Ispra (VA), Italy. E-mail: Folco.Casadei@jrc.it 1 Contents • Introduction • ALE formulation • ALE Fluid-Structure Interaction • Advanced topics and applications 2 1
Page 1 and 2:
European Laboratory for Structural
Page 3 and 4:
Proceedings of a Course on: Numeric
Page 5:
Programme of the Course 1. Introduc
Page 10 and 11:
Credits & Acknowledgments • Struc
Page 12 and 13:
Introductory Example (4) 7 Applicat
Page 14 and 15:
x Computational Framework • Gover
Page 16 and 17:
n+ 1 n ∆t n n+ 1 u = u + ( u + u
Page 18 and 19:
Scheme start-up and marching (2)
Page 20 and 21:
Stress Update • To solve the equi
Page 22 and 23:
Geometric non-linearities (3) Set u
Page 24 and 25:
Essential Boundary Conditions Essen
Page 26 and 27:
Exercise 0 - Ideal ballistics • M
Page 28 and 29:
Exercise 0 - Ideal ballistics (5)
Page 30 and 31:
Exercise 1 - Suspended mass (3) •
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 41 and 42: TITLE: PMAT04: motion of projectile
Page 43 and 44: sinon; tra2 = tra2 et ele; finsi; f
Page 45 and 46: Some results: horizontal and vertic
Page 47: PMAT01E This test is similar to PMA
Page 50 and 51: Hence: 11 −8 ES 2× 10 ⋅ 2.5×
Page 52 and 53: The stress history is: Note that th
Page 54 and 55: Analytical TEST11 Same as TEST01 bu
Page 56 and 57: Comparison of natural vs. engineeri
Page 58 and 59: The total external forces at the tw
Page 61 and 62: TEST13 Same as TEST01 but uses a 10
Page 63 and 64: The mass returns at the initial pos
Page 65 and 66: TEST22 To verify the independence o
Page 67 and 68: Linear dynamic analysis Consider th
Page 69 and 70: • The two waves f and g propagate
Page 71 and 72: Analytical solution For the bar imp
Page 73 and 74: Note that we have also put the Pois
Page 75 and 76: BARI08 We study the effect of the t
Page 77 and 78: ENDPLAY *==========================
Page 79: BARI10 Same as BARI09 but compares
Page 82 and 83: The system is initially at rest, an
Page 84 and 85: CAST MESH TRID NONL LAGR $ $ Dimens
Page 86: Final deformation (with superposed
Page 91 and 92: Building Vulnerability (2) P max co
Page 93 and 94: Euler equations • These equations
Page 95 and 96: Treatment of momentum equation •
Page 97 and 98: Time integration (3) Some pseudo-vi
Page 99 and 100: ALE description The ALE description
Page 101 and 102: Initial Conditions FE / FV • Equi
Page 103 and 104: Mean-based rezoning • Simple: opt
Page 105 and 106: Free surface modeling • Lagrangia
Page 107 and 108: Shock tube (4) • Assumed paramete
Page 109 and 110: Exercise 2 - Explosion in air-fille
Page 111 and 112: Exercise 3 - Bubble expansion in a
Page 113 and 114: Exercise 4 - External blast on two
Page 115: Exercise 5 - Confined detonation of
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: Analytical Conclusion: the pseudo-v
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. FV The t
Page 132 and 133: Example of shock tube problem solve
Page 134 and 135: The pressure distribution at t = 0.
Page 136 and 137: TANK01 Lagrangian solution: the who
Page 138 and 139:
The computed fluid pressures at the
Page 140 and 141:
The computed pressure histories are
Page 143 and 144:
Geometric data: The calculation is
Page 145 and 146:
The computed vertical displacements
Page 147 and 148:
TUBE04 ALE solution with a multi-ph
Page 149 and 150:
The input file for the test TUBE04
Page 151:
gotr loop 249 offs fich avi cont no
Page 154 and 155:
p6 = 70 20 0; p7 = 70 40 0; p8 = 40
Page 156 and 157:
$ high-pressure perfect gas (explos
Page 158 and 159:
* r_cais = 965.00e-3 ; h_cais = 156
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:
The fluid pressures and velocities
Page 170 and 171:
Further FSI Example Electric arc in
Page 172 and 173:
FSI Motivation (2) Fully coupled an
Page 174 and 175:
Permanent FSI treatment (2) • Upo
Page 176 and 177:
Geometrically complex cases • Bil
Page 178 and 179:
The FSA algorithm (2) • Effect of
Page 180 and 181:
The UP algorithm The method simply
Page 182 and 183:
Application to Finite Volumes • T
Page 184 and 185:
Exercise 2 - Explosions in simple d
Page 186 and 187:
Exercise/Example 3 - Wave propagati
Page 188 and 189:
Exercise/Example 4 - CONT problem T
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 (Com
Page 196 and 197:
Exercise/Example 8 : Steam Explosio
Page 198 and 199:
Exercise/Example 9 Explosion in Sec
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 P5 Exercise/Examp
Page 208 and 209:
Exercise/Example 12 - Sloshing (9)
Page 213 and 214:
Geometric data: For the fluid domai
Page 215 and 216:
ECHO RESU ALIC TEMP GARD PSCR SORT
Page 217 and 218:
The displacements of the horizontal
Page 219 and 220:
0.00000E+00 7.11050E+00 1.06860E+00
Page 221 and 222:
scen geom navi free !vect scco scal
Page 223 and 224:
Alternative 2 Use multi-phase multi
Page 225 and 226:
Geometric data: Complex shell-like
Page 227 and 228:
At 0.5 ms, spurious (non-physical)
Page 229:
Here is an example of the computed
Page 232 and 233:
DIME PT2L 293 PT3L 70 ZONE 3 ED41 3
Page 234 and 235:
2.10000E+01 8.33330E+00 2.10000E+01
Page 236 and 237:
FLUT RO 2.4278E3 EINT 0. GAMM 0.75D
Page 238 and 239:
The final mesh (colors indicate mat
Page 240 and 241:
GRIL LAGR LECT stru TERM ALE LECT f
Page 242 and 243:
The final bubble material mass frac
Page 244 and 245:
Geometric data: Materials The explo
Page 246:
Some results: global deformed mesh
Page 249 and 250:
The well-known Woodward-Colella tes
Page 251 and 252:
Numerical Solutions WOCO2D The mesh
Page 253 and 254:
ECHO * RESU ALIC GARD PSCR * SORT G
Page 255 and 256:
The EUROPLEXUS input file reads: WO
Page 257 and 258:
This example is a patch of four sho
Page 259 and 260:
CONT SPLA NX 1 NY 0 NZ 0 LECT symx
Page 261 and 262:
TITLE: Cavi51: steam explosion in a
Page 263 and 264:
p10p=p10 plus p0; p11=0.75 0 1.7; p
Page 265 and 266:
*opti donn 5; * vol3=vol2 syme plan
Page 267 and 268:
si (nn2 ega 1); str2=ei; sinon; str
Page 269 and 270:
The EUROPLEXUS input file reads: CA
Page 271 and 272:
Some results: final pressure distri
Page 273 and 274:
TITLE: PPLA04: unsteady flow throug
Page 275 and 276:
liq2=daller c1 c2 c3 c4 plan; lag=l
Page 277 and 278:
Some results: final fluid pressure:
Page 279 and 280:
PROBLEM: A deformable tank complete
Page 281 and 282:
s1 = p5 d 4 p9; s2 = p9 d 16 p12; s
Page 283 and 284:
The final pressure distribution in
Page 285 and 286:
RESULTS: Linear theory predicts a 1
Page 288 and 289:
The velocity at 200 ms is presented
Page 290 and 291:
compute a reasonable initial distri
Page 292 and 293:
FSA LECT 2 3 4 5 6 7 8 9 10 11 12 1
Page 294 and 295:
$ $ $ 91 92 93 94 95 96 97 98 99 10
Page 296 and 297:
The computed displacements of the t
Page 298 and 299:
TITLE: OILCUP2: uniformly accelerat
Page 300 and 301:
2526 2527 2528 2529 2530 2531 2532
Page 302:
Intermediate and final oil mass fra
Page 306 and 307:
Detailed contents • ALE descripti
Page 308 and 309:
Example 1 - Taylor bar impact 7 Exa
Page 310 and 311:
Example 3 - Coining Bilateral, Plan
Page 312 and 313:
Example 3 - Coining (5) Bilateral,
Page 314 and 315:
Example 4 - Explosion in a 2D box 1
Page 316 and 317:
Example 5 - Explosion in a corridor
Page 318 and 319:
Lagrangian contact • Non-permanen
Page 320 and 321:
Example 6a - Deep drawing A (2) Vel
Page 322 and 323:
Tube Crash (2) 35 Conventional cont
Page 324 and 325:
The Basic Pinball Method [Belytschk
Page 326 and 327:
Hierarchic Pinball Method (2) • H
Page 328 and 329:
Example 7 - Cable impact (2) l L v
Page 330 and 331:
Example 7b - Indentation (4) • 3D
Page 332 and 333:
Example 7b - Indentation (8) • Co
Page 334 and 335:
Why a Particle-Based Method? (2)
Page 336 and 337:
SPH Formulation (4) • To write th
Page 338 and 339:
Bird Strike [Courtesy of Snecma/CEA
Page 340 and 341:
Example 8 - SPH impacts (2) SONA01:
Page 342 and 343:
Spectral Elements (3) • To obtain
Page 344 and 345:
Spectral Elements (7) Coupling betw
Page 346 and 347:
Example 10 - Sediment valley • Ve
Page 348 and 349:
Example 12 - Matsuzaki valley • M
Page 350 and 351:
Time Step Partitioning (2) • Buil
Page 352 and 353:
Time Step Partitioning (6) • This
Page 354 and 355:
Treatment of links in partitioning
Page 356 and 357:
Example 12a - Taylor bar impact rev
Page 358 and 359:
Example 12a - Taylor bar impact rev
Page 360 and 361:
⎧ M1U 1 = F1+ R1 ⎪ ⎨M U = F +
Page 362 and 363:
Coupling at the Interfaces (7) ⎧
Page 364 and 365:
Multiple scales in time • Use app
Page 366 and 367:
Multiple scales in time (5) • The
Page 368 and 369:
Multiple scales in space (3) •
Page 370 and 371:
Example: simplified engine S 1 S 2
Page 372 and 373:
Example: aircraft Material is linea
Page 374 and 375:
Example 13 - Domains in 2D (2) •
Page 376 and 377:
Example 14 - Domains in 3D (3) •
Page 381 and 382:
Geometric data and materials: See s
Page 383 and 384:
sler cam1 1 nfra 20 trac offs fich
Page 385 and 386:
Page 387 and 388:
POIN TOUS CHAMELEM FICH TPLO FREQ 1
Page 389 and 390:
* p6=1.2E-2 0; p7=1.2E-2 1.0E-2; *
Page 391:
COUR 1 'dx_P1' DEPL COMP 1 NOEU LEC
Page 394 and 395:
p4s=p4 plus (0 0); tol=0.001; c1=p1
Page 396 and 397:
The deformed final mesh at 4 ms wit
Page 398 and 399:
The input file: INFS - 02 *--------
Page 400 and 401:
The final mesh and pressures are: T
Page 402 and 403:
ZONE 2 PMAT 2 Q4GS 3904 ECRO 858884
Page 404 and 405:
The final deformed mesh is: The fin
Page 407 and 408:
Problem description: This example i
Page 409 and 410:
LECT die punch holder TERM MASS 0.1
Page 411 and 412:
Page 413 and 414:
geom !navi free line heou poin sphp
Page 415:
An example of intermediate velociti
Page 418 and 419:
LOADING: The indenter is pushed int
Page 420 and 421:
The initial configuration (with par
Page 422 and 423:
The reaction force is: Approximate
Page 424 and 425:
The final velocities: 8
Page 426 and 427:
Page 428 and 429:
The initial configuration (with par
Page 430 and 431:
Page 432 and 433:
-5.10000E-01 1.33333E+00 2.00000E+0
Page 434 and 435:
Some hardening quantity in the targ
Page 436 and 437:
Final plastic streen in the target
Page 438 and 439:
Final plastic strain in the target:
Page 441 and 442:
Problem description: This example r
Page 443 and 444:
The input file is: HOLE - 06 $ ECHO
Page 445 and 446:
trac 2 1 AXES 1.0 'DISPL. [M]' yzer
Page 447 and 448:
Problem description: This example r
Page 449 and 450:
* fin lab3; nodf=c1p; * tpln=0.0 40
Page 451 and 452:
The input file is: VALL - PS $ ECHO
Page 453:
The final displacement norm in the
Page 456 and 457:
VEC1=2. 0.; * P1=P11; R1=P13; S1=P1
Page 458 and 459:
S3=P33 PLUS VEC1; * N1=2; N2=3; N3=
Page 460 and 461:
CHPO2=CHPO2 / SCAL1; TSTAT2 . &BCL9
Page 462 and 463:
AXTE 1E3 'Temps (ms)' *------------
Page 464 and 465:
The final displacement field in the
Page 466 and 467:
P9=40. 5. 5.; P14=20. 5. 5.; * BASE
Page 468 and 469:
CLIM1=(P5 ET P6 ET P7 ET P8); CLIM1
Page 470 and 471:
SUIT Post-treatment ECHO RESU ALIC
Page 472 and 473:
The final displacement field in the
Page 474 and 475:
tol=0.01E-3; * base=p0 d 5 p1; stru
Page 476 and 477:
eli=stru elem i; si (ega j 0); sq42
Page 478 and 479:
h3=8.1e-3; h4=16.2e-3; h5=24.3e-3;
Page 480 and 481:
The characteristic displacements in
Page 482:
The final yield stress field in the
Page 486:
The mission of the Joint Research C
show all

Download full text - ELSA - Europa

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?