Download full text - European Laboratory for Structural Assessment

More documents

Recommendations

Info

The system is initially at rest, and starting at the initial time a step-wise external pressure isapplied onto a square region near the top of the tower, simulating the impact zone. Thepressure stays constant for a duration of 500 ms, then goes linearly to 0 in 100 ms.CALCULATION:The calculation is performed up to 2000 ms. At the final time, the tower has undergone largeplastifications in the impact zone and to some extent in the surrounding area.RESULTS:Although no comparison with experiments can be done on this problem, the results appearphysically correct.POST-TREATMENTSeveral animations of the computed results from this calculation are available on theEUROPLEXUS Consortium Web site.REFERENCES:This calculation is detailed in:1) J.J. Lopez Cela, F. Casadei and P. Pegon: "Fast Transient Analysis of Thin ShellReinforced Concrete Structures with Drucker-Prager Model." SMiRT-14 Conference, Lyon(F), August 17-22, 1997.For details on the reinforced concrete material model, see:2) J.J. Lopez Cela, P. Pegon, F. Casadei: "Brittle Material Law with Drucker Prager YieldSurface and Softening Behaviour." Technical Note N. I.96.34, February 1996.The layered shell model is described in:3) F. Casadei, J.J. Lopez Cela: "A Multilayer Formulation for Shell Elements in PLEXIS-3C." Technical Note N. I.96.14, February 1996. (all these documents are available on theEUROPLEXUS Consortium Web site).Shell modelLayers282 of 556
Numerical SolutionsTOWE02The mesh generation file (K2000):*%siz 50*opti echo 1;*opti trac psc ftra 'towe02_mesh.ps';*opti titr 'TOWE - 02';opti dime 3;*ncir=10;nlon=20;*p11=47.90 0.0 0.0;p12=34.43 0.0 123.90;p13=28. 0.0 50.;p21=0.0 47.90 0.0;*opti elem seg2;line1 = para nlon p11 p13 p12;opti donn 5;*p4=0.0 0.0 0.0;p5=0.0 0.0 123.90;*opti elem tri3;surf1= line1 ncir ROTA 90. p4 p5;*oeilz=0.0 0.0 1000.0;oeilx=1000.0 0.0 0.0;oeily=0.0 1000.0 0.0;*p21=0.0 47.90 0.0;p22=0.0 34.43 123.90;p23=0.0 28.0 50.0;*opti elem seg2;line2 = para nlon p21 p23 p22;*opti elem tri3;surf2= line2 ncir ROTA 90. p4 p5;*p31=-47.90 0.0 0.0;p32=-34.43 0.0 123.90;p33=-28. 0.0 50.;opti elem seg2;line3 = para nlon p31 p33 p32;*surf=surf1 et surf2;tol = 0.0001;elim tol surf;**opti elem seg2;*bloc1 = cerc ncir p11 p4 p21;bloc2 =cerc ncir p21 p4 p31;bloc = bloc1 et bloc2;**mesh=surf et bloc et line3;elim tol mesh;*opti sauv form 'towe02.msh';sauv form mesh;*opti trac psc;trac cach mesh;fin;The EUROPLEXUS input file is:TOWE - 02$ECHO!CONV win383 of 556
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 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?