Download full text - European Laboratory for Structural Assessment

More documents

Recommendations

Info

The present document contains the presentation slides and the traces ofproposed exercises/examples related to the Course: Numerical Simulationof Fast Transient Dynamic Phenomena in Fluid-Structure Systems,given by F. Casadei at the Universitat Politècnica de Catalunya, Barcelona,14-17 May 2007.This is the fourth, competely revised edition of the document:(1st Edition): “Proceedings of a Course on: Numerical Simulation of Fast Transient Dynamic Phenomenain Fluid-Structure Systems, Presentation Slides and Trace of Proposed Exercises/Examples, JRC —Ispra, 11-15 April 2005”, JRC Special Publication N. S.P.I.05.93, 2005.(2nd Edition): “Proceedings of a Course on “Numerical Simulation of Fast Transient Dynamic Phenomenain Fluid-Structure Systems. Presentation slides and trace of proposed exercises/examples”, Coursegiven at Universitat Politècnica de Catalunya, Barcelona, on May 17-20, 2005, JRC Special PublicationN. I.05.139, June 2005.(3rd Edition): “Proceedings of a Course on “Numerical Simulation of Fast Transient Dynamic Phenomenain Fluid-Structure Systems. Presentation slides and trace of proposed exercises/examples”, Coursegiven at Universitat Politècnica de Catalunya, Barcelona, on May 15-18, 2006, and at Royal MilitaryAcademy, Brussels, on June 26-30, 2006, JRC Special Publication N. PB/2006/IPSC/1956, June 2006.Numerous corrections and ameliorations have been carried out in the presentation slides and a numberof new proposed practical exercises and examples have been added.4 of 556
Programme of the Course1. Introductiona. Introductory example of a FSI problemi.Application spectrum and goalsb. Modeling the structural domaini.Equilibrium equationsii.Explicit time integration schemec. Treatment of essential boundary conditionsi.The Lagrange multipliers methodii.Solving the linear system2. ALE formulationa. Modeling the fluid domaini.Euler equationsii.Finite Element discretizationiii.Finite Volume discretizationb. Mesh rezoning algorithmsi.Motivationii.Mean-based algorithmsiii.Giuliani’s algorithmc. Free surface modeling3. ALE Fluid-Structure Interactiona. Motivationb. Classification of FSI algorithmsc. Equilibrium-based methodsi.The Uniform Pressure (UP) methodii.Shortcomingsd. Geometrical methodsi.The FSA/FSR methodii.Shortcomingse. A combined methodi.The FSCR method4. Advanced topics and applicationsa. ALE description of structuresb. Non-conforming FSIc. Lagrangian contacti.Classical methodsii.Pinballsiii.SPHd. Spectral elementse. Spatial time step partitioningf. Domain decomposition5 of 556
Page 2 and 3: The mission of the Institute for th
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:
Page 214 and 215:
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 432 and 433:
The resulting final deformed mesh w
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:
Page 466 and 467:
The reaction force is:Approximatean
Page 468 and 469:
The final velocities:8468 of 556
Page 470 and 471:
Page 472 and 473:
Page 474 and 475:
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
show all

Download full text - European Laboratory for Structural Assessment

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

Delete template?

Save as template?