The 9/11 Forum

Intelligent and evidence-based discussion of 9/11 issues

Skip to content

v
Welcome
Welcome!

Our vision is to provide a home to sincere 9/11 researchers free from biased moderation and abusive tirades from other members.

You are currently viewing our boards as a guest, which only gives you access to view the discussions. Feel free to register to request membership. By joining our free community, you will have access to post topics, communicate privately with other members (PM), respond to polls, upload content, and access many other special features. All potential members will be subject to an interview via email and only sincere and responsible researchers will be approved. See the forum guidelines for more information.

WTCCS Spreadsheet Crush-Down Model

Analysis of airplane impacts, fires and collapse theories and examination of related evidence.

Re: WTCCS Spreadsheet Crush-Down Model

Postby OneWhiteEye on Thu Jun 11, 2009 6:18 pm

May I suggest doing whatever it takes to make your simulation match. It's not dry labbing, it's exploratory data analysis. All you have to play with is mass and stretch, and you have, to some extent. Take it all the way. Use your default mass distribution and increase stretch until you get a match. Use a best guess stretch and add mass above until it's a fit. Allow for +/- at least 0.5 seconds start time sync error. Put a box around this problem. If the solution space which does fit is not entirely unreasonable, then there's not much that can be said, except perhaps that the model works better than it should.

Personally, I think there's a lot of mileage yet to be obtained with 1D simplified models. That doesn't mean there's any expectation it can apply to any and all specific aspects of the real collapse. I get where you're going with the resistance of zero and all that, but there are modes of failure that allow for faster failure propagation than slab impact, and that's a fact. Whether they apply is another discussion, but the fact that we're examining ejecta from a limited spatial region most certainly admits it, and runs somewhat contrary to the slab model.

I would be interested in how much mass you have to add to the run that came 3.5 stories short. There may be problems with my code, and the assumption of uniform mass of impacted floors means they'll never quite agree, but we should be getting closer results if everything's in order. I'll check my code.

Due to schedule uncertainties, my participation may be spotty over the next few days, or not. If so, I'm not igging anyone.
OneWhiteEye
 
Posts: 1220
Joined: Sat Jul 05, 2008 9:40 pm

Re: WTCCS Spreadsheet Crush-Down Model

Postby Hambone on Thu Jun 11, 2009 9:32 pm

femr2 wrote:Has anyone produced an estimation of the hat truss mass ?


My paper and spreadsheet include it but the information was rather sketchy in the NIST reports.
Hambone
 
Posts: 236
Joined: Wed Jul 02, 2008 9:09 pm

Re: WTCCS Spreadsheet Crush-Down Model

Postby femr2 on Thu Jun 11, 2009 10:46 pm

OneWhiteEye wrote:Remember this only applies to the southern end, it's application of the model to a more specific vertical slice in an attempt to match a later ejection in that horizontal region.

Of course. Going through the process is leading me towards dropping the initial 'failure' floor down anyway, so am just trying to ensure that shift is backed up with visual and descriptive clarification.

It would represent the extreme lowest point, but that's precisely the case. IF the model can work under these circumstances (...), it can only be expected to work in a 1D slice representative of the conditions at that location.

I've been having a think about this, and may consider setting up a 'grid' consisting of multiple instances of the spreadsheet model to try and express more defined regional behaviour. This would entail at minimum a 3x3 grid of spreadsheets, the core, and the 8 surrounding regions. The grid could be made finer, but the work involved may be OTT. Each spreadsheet base parameters could then be made more specific to each region, but the results combined within the visual module...
Of course, there are fairly serious issues with such a method, not least the probable initial lack of interactivity, but it's worth thinking about IMO.

What if the problem were reversed?

I assume you mean matching to slowest rate of descent. If so, then some resistance/energy sink would be included. The discussion is still for a zero-resistance, 100% mass, no energy conversion scenario after all. (I started the "What Does Collapse Time Mean" thread for defining the various different 'types' of collapse time. It may become useful here a bit later.)

It is my opinion that the northernmost, uppermost slice of the building is still above the top of the green region. Can you tell me what tweaks need to be applied to the spreadsheet model to match that? (in case the tweaks should result in arrest, you'll be obligated to provide a mechanism for continuation! haha)

As above, some 'resistance' would need to be included. Piece of pie. I can include structural support failure energy, concrete crush energy, mass loss, etc...all of which can result in longer collapse times.
Treatment of collapse arrest has been in the model since day one. The method employed is that if the energetics calcs result in a deficit of energy, the resultant full post floor impact velocity is set to zero, the failure floor is recorded, and progression is allowed by allowing the cap mass (including any coalesced mass) to drop from rest...continues...
http://www.youtube.com/watch?v=gohc5lMO5Q8
Image
Image is post arrest, so the debris is 'cut-off', the rendering of the descending mass is halted, but the cyan bar still progresses the tower after the collapse is allowed to progress.

As example, if I say that this view most accurately represents the maximum split elevation, it's just an impression. If I'm wrong, oh well.

It's a great video for seeing the initial north-south behaviour of WTC 1. I've included time and original frame numbering on the video so it can be used for further study.

The mere presence of the ejecta at the level indicated by the lower line indicates this is a more accurate position for an extreme front, especially towards the south.

I'm inclined to agree there (for the purpose of making use of 1D model more valid).

As you note, the video also depicts the progression south to north. In discriminating conditions along the north-south axis, this is worth - on the outside - a second of head start or offset down one floor.

I'll be rendering a video synchronised around the time of the annotated still frame...to see what happens. (I'll have to use a dynamic camera setup and a video source split, so it'll take a little while to do...)

...Visible evidence from that video must be reconciled with other knowns. One of which is the bowing extent I mentioned earlier. This is strong, independent corroboration for choosing a value no higher than 94. All was I was asking for was two stories beyond that!

Yeah, but I'm English and stubborn :) For the specific purpose of speculative research, no probs.

The other known is the crop of one of your frames I posted above. Going back to those calculations, and judging crush front location by exactly the same standard as the five second mark, it's already 5 stories ahead.

I could really do with a higher resolution copy of the source video. Anyone..... ?

So, either this gets applied to the initial conditions for the model or the model tries to fit this early, huge discrepancy, and I don't think it would succeed...

Will be syncing to early events to try and force a 'fit'.

Can we agree on 93-94 as nominal, and 95 to 92 as extremes, for the purposes of this context? If so, then it's 92 we must consider, and so you have.

In this context, 92. In a wider sense, I'm leaning toward 95, with slightly earlier start point.

The mass distribution may be quite different

Yes, and floor heights are not constant.

Best figures are always best(!) but there has to be caution in applying precision inappropriately.

More to do with adding mass to the cap itself, rather than dropping the start floor down.

Do you have a handy graph of the upper floor masses contained in your sheet?

Image
http://femr2.ucoz.com/photo/1-0-68
And, AHA, looks like I already have the truss mass in place. Have checked and my source is the NIST data. Excellent.

I'll cut to the chase now, since I'm not going to defend to death a 1D simple model in an attempt to match reality. If it comes up short, I will fail the model, simple as that.

We haven't begun to include anything other than conservation of momentum that will retard the descent in any way yet...modelled collapse time increases significantly when such is included.
The model can't be used to 'prove' anything as such, but I do think that it would be very useful, to me at least, to progress through discussion of inclusion of each of the available but currently omitted factors...mass loss, some crush...
I hope the thread can delve into the implications of each of the factors, and also hopefully result in addition of further factors.

without video of the WTC1 S face

There are a couple of useful video's. I'll locate and link.
femr2
 
Posts: 1010
Joined: Thu Apr 23, 2009 12:08 am
Location: UK

Re: WTCCS Spreadsheet Crush-Down Model

Postby femr2 on Thu Jun 11, 2009 10:49 pm

Hambone wrote:
femr2 wrote:Has anyone produced an estimation of the hat truss mass ?


My paper and spreadsheet include it but the information was rather sketchy in the NIST reports.

Hm. I may be double-dipping then. I'll trace back through source data and confirm. (I have NIST data included)
femr2
 
Posts: 1010
Joined: Thu Apr 23, 2009 12:08 am
Location: UK

Re: WTCCS Spreadsheet Crush-Down Model

Postby femr2 on Thu Jun 11, 2009 11:28 pm

OneWhiteEye wrote:May I suggest doing whatever it takes to make your simulation match. It's not dry labbing, it's exploratory data analysis. All you have to play with is mass and stretch, and you have, to some extent. Take it all the way.

...

Use your default mass distribution and increase stretch until you get a match.

Increasing floor thickness to 1.13m generated the desired result. (floor 60 @ 5s)
At the rates of progression, +-0.5s around 5s takes the range down to floor 54, or up to floor 64/65.

Use a best guess stretch and add mass above until it's a fit.

Adding 4e7 kg to the cap generated the desired result, nearly doubling the initial cap mass to 7.6e7 kg. No change in 'stretch' was used, remaining at 0.7492m.
At the rates of progression, +-0.5s around 5s takes the range down to floor 54, or up to floor 65.

Allow for +/- at least 0.5 seconds start time sync error.

With original figures...
At the rates of progression, +-0.5s around 5s takes the range down to floor 58, or up to floor 68.

Personally, I think there's a lot of mileage yet to be obtained with 1D simplified models. That doesn't mean there's any expectation it can apply to any and all specific aspects of the real collapse.

Agreed.

I get where you're going with the resistance of zero and all that, but there are modes of failure that allow for faster failure propagation than slab impact, and that's a fact. Whether they apply is another discussion, but the fact that we're examining ejecta from a limited spatial region most certainly admits it, and runs somewhat contrary to the slab model.

Agreed, but could you outline the modes of failure ? (May be able to address them)
femr2
 
Posts: 1010
Joined: Thu Apr 23, 2009 12:08 am
Location: UK

Re: WTCCS Spreadsheet Crush-Down Model

Postby femr2 on Sat Jun 13, 2009 4:26 am

(Slight derail, but couldn't resist)

LS-Dyna alternative:

Image
http://www.youtube.com/watch?v=LPsW9cVPMgw

It's based on the following technology:
http://extremeloading.com/

What I find very curious is it's specific development since 1995 for 'progressive collapse', the contrasting intimate relationship with the demolition companies and that they have included the Murrah building as one of their simulations.

I don't suggest it's discussed on this thread, but...very interesting.
femr2
 
Posts: 1010
Joined: Thu Apr 23, 2009 12:08 am
Location: UK

Re: WTCCS Spreadsheet Crush-Down Model

Postby OneWhiteEye on Sat Jun 13, 2009 9:58 pm

femr2 wrote:I've been having a think about this, and may consider setting up a 'grid' consisting of multiple instances of the spreadsheet model to try and express more defined regional behaviour.

Lots of work! It would be cool. Don't underestimate the value of running the existing sheet using only an adjustment to the top mass.

I assume you mean matching to slowest rate of descent.

Yes, there's a real min and max displacement value seen in the videos. It's racing way ahead on the south side and actually remains stationary for a while on the NW corner. A model of aggregate behavior, not just uniform in the horizontal dimension, but also requiring full participation of all mass at each level, cannot be compared directly to a heterogeneous front proceeding at different rates. The best to be hoped for is fairness in including the opposite extreme in some meaningful way.

One way is by fitting the slowest extreme, too, to see the difference.

As above, some 'resistance' would need to be included.

OK, so you can jack the energy up until it arrests then....

Treatment of collapse arrest has been in the model since day one. The method employed is that if the energetics calcs result in a deficit of energy, the resultant full post floor impact velocity is set to zero, the failure floor is recorded, and progression is allowed by allowing the cap mass (including any coalesced mass) to drop from rest...continues...

...just let it drop again at the appropriate time, all controlled by parametric settings, piece of pie. What physical correspondence does it have? None. There's nothing driven from a constitutive stance, because it's outside the scope of the model. The NW corner receives a substantial punch from above and survives. Nothing in this model is going to get it going again except the whim of the programmer. It has arrested. The model has 12 stories of building still on top of it! In reality, there's only air above, being undermined from below is what got it going, and it may have initially dropped at rates greater than freefall. How would that be addressed in the model? Arbitrarily reduce static capacity to zero AND temporarily increase gravity?

Why can't we just cheat similarly when trying to model the fastest advance?

Better but still questionable and loose, would be considering the extremes together, making a single weighted composite displacement which the model tries to match. There's the crush front we're discussing at one extreme, and the NW corner standing there while everything else falls around it at the other extreme. I doubt the discrete slab model would have any difficulty fitting that profile with some energy sinks and mass loss included. Looks like your first illustration incorporating stretch is at least a candidate for a match. It's questionable how much detail can be gleaned from such a fitting, but at least it's not arbitrary and it's arguably within the scope of the model.

If you take the spire tip as the slowest extreme, what sort of dissipative energies must be assigned, and are they remotely realistic? This is the only point I was trying to make in suggesting a fit to the slowest extreme. It's a no-go from the predictive modeling perspective, the behavior dictated by decree not by dynamics. If dissipated energy is increased until arrest, I guess you've found the amount required to stop the collapse, that's pretty interesting. That's what 1D models are for.

There is another approach to fitting arrest. Mass can be removed from above to give arrest immediately, that would be closer to the real situation in the NW corner than increasing energy dissipation. Re-initiation being completely outside the scope of this model.

Why resist going the other way and adding mass to simulate the other side, where the mass was last seen tipping towards, and coincidentally is the fastest visible progression? How much mass has to come off one side to arrest under energy sinks you feel comfortable with? Great, now add that amount to the other side, since that's where it went anyway, and run again. If you can't match the front at 5 seconds, starting from 92 and with a realistic stretch, then maybe the fast extreme is also outside the scope of the model?

Incidentally, I do have a 1D physics-based slab simulation, with reasonable parameters including static strength, that can reinitiate crush-down from ongoing crush-up. There's nothing that says there only needs to be one generalized coordinate in 1D. But I still can't model the NW corner in this framework, even if (not necessarily unrealistic) vibrational energy accumulated in the lower structure could re-precipitate global failure, because there's still 12 floors of mass on top and the NW corner had none! The model I'm currently working with doesn't allow any mass bypass, and neither does yours. This is a double-edged sword.

If so, then some resistance/energy sink would be included.

While it's unrealistic to exclude energy sinks per se, this model has other unrealistic aspects which may tend to counter the omission, to what degree is infinitely debatable and some I've mentioned. Fully inelastic collisions make a lot of KE go away, because of a definition or constraint imposed intellectually. It may as well be magic because there are very few materials that display near-perfect inelastic behavior over a wide range of velocities. Remember Silly Putty? It would bounce if thrown at a hard surface, but was like modeling clay when working in the hands.

It is known for this model, implied in the collision, the mass below is swept up and accelerated to the upper block velocity. It's implicit in the dynamics that have been defined. Conversely, and also by definition, KE is forced to go into unrecoverable energy distribution: permanent deformation, fracture, lossy ejections, heat, sound, the works. And it does, obviously. The distribution over time and location of the PE is quite debatable, indeed it's part of the central mechanics debate. The simple model doesn't care whether all of those sinks are truly unrecoverable (like vibration), and to what extent energy can be mediated ahead of Zone B. I don't know how the energy distribution actually went, so there's little I can add to sorting this out. I'm quite confident it does goes beyond simple inelastic collision losses, so adding dissipative sinks would be appropriate for the purposes of comparison.

But comparison to what? Only an aggregate quantity, such as overall collapse time as defined by (say) 90% of the mass reaching rest at pile location. Comparison to the extreme leading front in a limited area? No, not a valid comparison, per above. An interesting one, for sure, and even educational, but not one from which one can determine that there was some sort of violation of physical laws governing a PE-driven collapse. It's reaching too far.

The divvying of energy is one thing, the motion of mass is another. The biggest drawback of this model is not inelastic collision, in fact that seems to be a virtue; the hitch is that mass elements must remain in their original sequence. Nothing can bypass, or race ahead, or lag behind - all mass is entrained at collision time. I can look right at any of the videos and see that's not how it went, therefore am thankful for whatever comes out of a model that has such assumptions.

Truth is, we're not that far off from a match under the apparent conditions of the slightest cheating. Given the model, I'd call that significant. And a piece of pie.

In this context, 92. In a wider sense, I'm leaning toward 95, with slightly earlier start point.

I tend to stick around 96-98, or higher, never used 92 before this discussion.

Image


Looks pretty uniform except for the last 4, and that adding ~3 story-masses at the top would make a simple uniform model more accurate. Did I read it correctly?

We haven't begun to include anything other than conservation of momentum that will retard the descent in any way yet...

Except fully inelastic, head-on collisions that immediately entrain all the mass at a given level...

The model can't be used to 'prove' anything as such, but I do think that it would be very useful, to me at least, to progress through discussion of inclusion of each of the available but currently omitted factors...mass loss, some crush...

No doubt very useful.


---------------------

Agreed, but could you outline the modes of failure ? (May be able to address them)

Some has been touched on already above, including elasticity, but these two words encapsulate a lot: deflection and ablation.

The first is any wedging effect. Axe, ship bow, upper block tilting in that direction before dropping. Material in front of the advancing front is no longer required to be swept up in totality. Look at WTC2's upper block shove the E perimeter out. It was like a paring knife. Much of the lower block on the N side directly below the failure point continued to stand well after the passage of the block. All of the mass from above went down, not all of the mass below was entrained, easily observed directly on video with my own eyes - slab model not generally applicable. QED

I think a similar thing happened with WTC1. Not only does tipping redistribute the mass you've so carefully assembled into a 1D world, it creates a geometry that allows a fairly strong portion of the upper block to do more damage than it receives. The structure is not homogeneous, so it isn't always like against like. The upper perimeter in WTC2 may have been badly damaged in its sliding contact with the lower while shearing it away, but it retained the majority of its shape a long time. Even if it was rubble-ized by explosives!

Why should WTC1 be radically different? One difference may be a greater percentage staying over the perimeter area longer. The antenna mast motion is consistent with gradual continued rotation of the upper block, with the hat truss area staying more or less cohesive through the first half of the fall. This is a dense and more rigid subset of building impinging on the structure below in an orientation favorable for cleaving and or shoving material aside instead of sweeping it up. It's also consistent with the lateral location of the leading front seen in the video.

What about the part of the upper block that's dragging over the core, some which remained standing after the wave passed? This is where ablation comes in. We're working with an exclusive accretion model for the dynamics of collision, so it's not a big deal that mass shedding is left out. In an local ablative model. KE can be retained via the loss of mass, as opposed to continued application of resisting force. It may be more energetically favorable for the core impulses to shear the upper block than to absorb KE from it in a sustained fashion. The rubble falls more slowly, but it is detached.

Think road rash. I saw a guy go off the back of his motorcycle and slide across pavement like he was on ice. I'm not sure when he would have stopped if hadn't started tumbling like his bike. Getting a leg sucked underneath the rest of his body imparted angular momentum, stopped the ablation, which brought things to a more rapid conclusion. His arse certainly wasn't transmitting the drag force to the rest of his body, but it was leaving a ghastly stripe on the road.
OneWhiteEye
 
Posts: 1220
Joined: Sat Jul 05, 2008 9:40 pm

Re: WTCCS Spreadsheet Crush-Down Model

Postby David B. Benson on Sat Jun 13, 2009 11:28 pm

From NCSTAR1--6 and NCSTAR1--6D I determined the extra mass at the top of WTC 1 as the equivalent of 8.8 of the ordinary stories below.
David B. Benson
 
Posts: 861
Joined: Sun Jul 06, 2008 12:29 am

Re: WTCCS Spreadsheet Crush-Down Model

Postby femr2 on Sat Jun 13, 2009 11:39 pm

OneWhiteEye wrote:Yes, there's a real min and max displacement value seen in the videos. It's racing way ahead on the south side and actually remains stationary for a while on the NW corner.

I started this thread (Link) to clarify the various rates of descent. I would view the NW corner remnant as that, a remnant, in the same way that remnants of the core remained standing for a short time after both the lowest 'crush front' and the upper/latest descending masses had already passed.
I think it's important to separate such remnants, specifically external perimeter remnants, from the internal descending masses.
(I could add mass transfer data to the sheet, such that the mass of that segment is applied to a lower floor, but it would make almost nil difference, and the validity of doing so would be a bit suspect also.)

A model of aggregate behavior, not just uniform in the horizontal dimension, but also requiring full participation of all mass at each level, cannot be compared directly to a heterogeneous front proceeding at different rates.

True, but one thing I am very curious about is whether it is possible for such aggregate behaviour to result in a descent time which exceeds that of the 1d alternative. My mindset there is that a series of lower mass impacts occuring over a wider span of time would seem to lean towards a longer 'general' descent time.
One 'simulation' which may be useful is to determine what KE by surface area is required to pass through a floor assembly, and whether there would then be sufficient KE to pass through additional floors (as in punch a hole through)

OK, so you can jack the energy up until it arrests then....

Yes.

...just let it drop again at the appropriate time, all controlled by parametric settings, piece of pie.

Not quite sure what you mean. An 'arrested' descending mass is immediately allowed to fall from an initial velocity of zero. What do you mean by parametric settings ?

The NW corner receives a substantial punch from above and survives. Nothing in this model is going to get it going again except the whim of the programmer.

If it was a floor section I'd agree, but as it's part of the perimeter, I don't. There's no asynchronous collision implied by the late descent of the corner, as it's still attached to the floors below. I suggest transferring it's mass to the floors below, as outlined above.

It has arrested. The model has 12 stories of building still on top of it!

That's part of the reasoning behind immediately starting from zero velocity, to emulate load based failure. It's a freefall calc which doesn't involve and resistances (even if they are in the scenario) so it leans in favour of progression. A 'reasonable' mechanism given the limitations of synchronous impacts I think.

In reality, there's only air above, being undermined from below is what got it going, and it may have initially dropped at rates greater than freefall.

Could you clarify 'dropped faster than freefall' ? Do you mean propelled from a prior impact ? (I assume so, just clarifying)

How would that be addressed in the model? Arbitrarily reduce static capacity to zero AND temporarily increase gravity?

I think it's approaching out of scope. I'm looking to perhaps go to code for an asynchronous impact model which draws some data (mostly mass) from the sheet model. What I'd prefer is to get my hands on the ExtremeLoading software linked-in above :)

Better but still questionable and loose, would be considering the extremes together, making a single weighted composite displacement which the model tries to match.

Perhaps, though I'd prefer to progress mechanics to include cap crush-up, and also perhaps gradual additional compaction of the crush zone, to favour generating the lowest descent times possible, then including some resistances and trying to validate each additional factor.
My view is that a few tons of debris falling above the primary 'pile driving mass' is not really relevant to what I'm trying to achieve, which is basically...with everything set to values which have been discussed and validated in some form...what is the simulated collapse time ?

If you take the spire tip as the slowest extreme, what sort of dissipative energies must be assigned, and are they remotely realistic? This is the only point I was trying to make in suggesting a fit to the slowest extreme. It's a no-go from the predictive modeling perspective, the behavior dictated by decree not by dynamics. If dissipated energy is increased until arrest, I guess you've found the amount required to stop the collapse, that's pretty interesting. That's what 1D models are for.

The amount of energy to arrest is really dependant upon what factors are included (structural resistance, mass loss or crush), and at what point arrest occurs. It's quite possible to set parameters such that arrest occurs at various points for various reasons, and the effect of the combined use of factors can produxe some quite interesting results. That's one reason I'd like to discuss the inclusion of them, so that I can clarify/improve the mechanisms I use...for example, I have switches which select between using either a constant value for the energy required to break the structural supports for a floor, or a ramped value. Crush mechanism which selects between a constant or uses impact KE to determine crush scale...all that kind of thing. Might be worthwhile separating each factor into a separate thread of you are interested in discussing ?

How much mass has to come off one side to arrest under energy sinks you feel comfortable with? Great, now add that amount to the other side, since that's where it went anyway, and run again. If you can't match the front at 5 seconds, starting from 92 and with a realistic stretch, then maybe the fast extreme is also outside the scope of the model?

Will model and post a bit later.

The model I'm currently working with doesn't allow any mass bypass, and neither does yours. This is a double-edged sword.

Could fudge it, but it's going to make very little difference to the subsequent behaviour.

Fully inelastic collisions make a lot of KE go away, because of a definition or constraint imposed intellectually.

I make the full energy loss available to subsequent calcs, so I use it in 'zero KE loss' sense...other than the velocity change.

The distribution over time and location of the PE is quite debatable, indeed it's part of the central mechanics debate.

Indeed.

The divvying of energy is one thing, the motion of mass is another. The biggest drawback of this model is not inelastic collision, in fact that seems to be a virtue; the hitch is that mass elements must remain in their original sequence. Nothing can bypass, or race ahead, or lag behind - all mass is entrained at collision time.

Yes. It would be great to progress to asynchronous collisions, elastic collisions and fragmentation, but it's probably out of scope. ExtremeLoading :)

Image
Looks pretty uniform except for the last 4, and that adding ~3 story-masses at the top would make a simple uniform model more accurate. Did I read it correctly?

Image
http://femr2.ucoz.com/photo/1-0-70-3
Needs to be extrapolated to include the rest of the floors IMO.

I'll get to the later half of the post shortly.
femr2
 
Posts: 1010
Joined: Thu Apr 23, 2009 12:08 am
Location: UK

Re: WTCCS Spreadsheet Crush-Down Model

Postby femr2 on Sat Jun 13, 2009 11:51 pm

David B. Benson wrote:From NCSTAR1--6 and NCSTAR1--6D I determined the extra mass at the top of WTC 1 as the equivalent of 8.8 of the ordinary stories below.

The figures tend to relate to the rated maximum load rather than specifying what that load actually was and what it was comprised of.

I have also used NIST data (NCSTAR1-2a), and specified the roof, aircraft and antenna masses separately.

Mass of Roof and Antenna 4300000 kg
Mass of Aircraft 128639 kg

(These are in addition to the graphed values)

There are many conflicts within the NIST report.

Without wading through all 940 pages, have you looked at how that figure was derived ?
If you could post the breakdown it would be appreciated.
femr2
 
Posts: 1010
Joined: Thu Apr 23, 2009 12:08 am
Location: UK

Re: WTCCS Spreadsheet Crush-Down Model

Postby OneWhiteEye on Sun Jun 14, 2009 2:23 am

To keep central points from being obscured in a sea of words:

The maximum velocity obtained from a basic 1D inelastic slab model is not a cosmic speed limit of some kind. Part of what you say acknowledges this, but the posed question asks why there is such a discrepancy. After lowering the initiation zone and adding stretch, both values within sensible bounds, your simulation came up 3.5 stories short. That might be enough to justify further inquiry, I think it does. On the other hand, I find it fascinating that a slab model can get so close to the actual ejecta front.


femr2 wrote:True, but one thing I am very curious about is whether it is possible for such aggregate behaviour to result in a descent time which exceeds that of the 1d alternative.

Me. too. It's a very worthwhile question. I would be surprised if it could. Therein lies some of the value of the 1D model. Unfortunately, I doubt it's possible to reliably map the actual 1D mass distribution over time in order to form a comparison. What we're look at here is a expulsion front because it can be observed and measured with some accuracy. After 5 seconds, it could be up to two floors ahead of the bulk of mass in a crush front, the actual distribution totally unknown.

My mindset there is that a series of lower mass impacts occuring over a wider span of time would seem to lean towards a longer 'general' descent time.

By definition, if occurring over a wider span of time, they will lead to longer descent times. The question really is, is the reduced mass in this scenario associated with even more reduced resistance and energy dissipation? I don't know if it is but it could be. Then it will fall faster and lead to shorter times.

I would view the NW corner remnant as that, a remnant, in the same way that remnants of the core remained standing for a short time after both the lowest 'crush front' and the upper/latest descending masses had already passed.
I think it's important to separate such remnants, specifically external perimeter remnants, from the internal descending masses.

In other words, it's outside the scope of the model. That's all I'm saying about the extreme visible ejecta at 5 seconds. It's a remnant of another sort. Precursor? Please indulge me in the comparison between extremes; I don't view them that differently when it comes to applicability of model, thus conceptually avoiding any consideration of causative mechanisms behind the discrepancies. It's a binary choice, model sufficient or not? Reasons (relatively) unimportant.

One 'simulation' which may be useful is to determine what KE by surface area is required to pass through a floor assembly, and whether there would then be sufficient KE to pass through additional floors (as in punch a hole through)

Yes. With a big heavy thing behind it, perhaps rigidly attached. Pressure is the keyword here.

Not quite sure what you mean. An 'arrested' descending mass is immediately allowed to fall from an initial velocity of zero. What do you mean by parametric settings ?

What you choose to do, as the programmer, is allow it to fall immediately and, because it's more useful than just letting it sit there, it's a good thing. But does it reflect the real dynamics? If it does, it's only by coincidence, matching the case where arrest occurs followed immediately by continuation. That's not what happened in the NW corner, it remained at rest for a period following the punch, then went down abruptly later. If you seek to match its displacement over time, there is no mechanism to account for the delay between impact and descent. Adding a delay parameter for a local t0 would allow matching initiation, but it would be arbitrary unless based on actual underlying physical principles outside the scope of the model.

That's part of the reasoning behind immediately starting from zero velocity, to emulate load based failure.

It's an excellent approach, I find no fault whatever with that. It is, however, a specialization of the general case of subsequent overload, where the ensuing overload is effective immediately. Subsequent overload from above is not what happened with the NW corner, so that in itself is making a round hole with a square peg. If the NW corner, as a remnant, does not apply, I have some data on a cardhouse collapse that does. Displacement stops, arrest occurs and, at a later time, crush down resumes. Surprisingly 1D, Zone B/C laying on top, motion arrested... I'll post a graph. There is a real reason for the structure re-initiating later; good luck analyzing it in a quantitative fashion.

It's a freefall calc which doesn't involve and resistances (even if they are in the scenario) so it leans in favour of progression.

In that context, it can only progress, as only structural resistance can produce arrest. Thus, to model the NW corner as a minima using this model, resistance (of arresting magnitude!) must be supplied, as you indicate. It is whatever it needs to be to arrest, even if that's not realistic, and the model's applicability stops there (sans enhancement) because it cannot predict the delay before next motion. Conceptually, it's pretty much the same as reducing resistance for the sake of matching the other extreme. Making the structural resistance zero is overly simplified, but a good starting point.

Progression or not is different from progression rate, too.

Could you clarify 'dropped faster than freefall' ? Do you mean propelled from a prior impact ? (I assume so, just clarifying)

Only admitting the possibility that, as a mere remnant, any impact below of previously falling material could impart impulse or even ongoing downard force to the corner, yanking it down. It's not a free body, and there are other things impacting the foundation to which it's connected. Outside the model, for sure.

If it was a floor section I'd agree, but as it's part of the perimeter, I don't. There's no asynchronous collision implied by the late descent of the corner, as it's still attached to the floors below. I suggest transferring it's mass to the floors below, as outlined above.

I may be a step or more behind the current scheme, definitely behind what could be. Transferring mass is fine, it's effectively the same as loading the top more. I already found an initial mass distribution (reasonable according to DBB), using stretch and story values close to optimal, ignored tip and advanced pressurization and satisfied the displacement. This is a surprising result to me for apparently the opposite reasons!

Could fudge it, but it's going to make very little difference to the subsequent behaviour.

For this spreadsheet it would almost have to be arbitrary guesswork but, in general modeling case, I could not disagree more.

I make the full energy loss available to subsequent calcs, so I use it in 'zero KE loss' sense...other than the velocity change.

I can understand keeping it around for accounting, but how do you apply it back in to the dynamics?

Yes. It would be great to progress to asynchronous collisions, elastic collisions and fragmentation, but it's probably out of scope. ExtremeLoading

It's good to use more sophisticated methods and tools to explore the problem domain. Personally, I'm starting simple intentionally, and remaining so with more intention. I believe one of the shortcomings of the NIST simulations is the unprecedented size and complexity of the models without having the basis of prior comparable modeling projects, or at least a long chain of substantial, incrementally scaled prototypes. If their simulations are correct by any measure, it's either unprecedented genius or pure dumb luck.

Image

Righteous!!!! Thank you! Big bumps at the mechanical floors.
(my eyes glaze over at grids of numbers)
OneWhiteEye
 
Posts: 1220
Joined: Sat Jul 05, 2008 9:40 pm

Re: WTCCS Spreadsheet Crush-Down Model

Postby femr2 on Sun Jun 14, 2009 3:33 am

OneWhiteEye wrote:After lowering the initiation zone and adding stretch, both values within sensible bounds, your simulation came up 3.5 stories short.

I should make it clear that the additional 3.5 storeys is achieved in ~0.3s, so it's a very marginal difference. That's why I'm pushing for the inclusion for some forms of resistance :)

My mindset there is that a series of lower mass impacts occuring over a wider span of time would seem to lean towards a longer 'general' descent time.

By definition, if occurring over a wider span of time, they will lead to longer descent times. The question really is, is the reduced mass in this scenario associated with even more reduced resistance and energy dissipation? I don't know if it is but it could be. Then it will fall faster and lead to shorter times.

Could you clarify how resistance would be reduced and less energy dissipated ?
A virtual view would be that for small mass impacts, failure of a floor will not occur until enough mass has built up to cause the floor to fail due to load rather than impulse. An over simplification, but that's part of what I'm saying with the above.

Please indulge me in the comparison between extremes; I don't view them that differently when it comes to applicability of model, thus conceptually avoiding any consideration of causative mechanisms behind the discrepancies.

Of course. Am posing only minimal resistance. (Pun intended :wink: )

With a big heavy thing behind it, perhaps rigidly attached. Pressure is the keyword here.

Agreed, though the behaviour there would surely be that the smaller area mass would, say, punch a hole in the floor causing no other failures (just eroding the floor a bit) and the effective impact would not occur until the big heavy thing actually hit that floor ?

That's part of the reasoning behind immediately starting from zero velocity, to emulate load based failure.

It is, however, a specialization of the general case of subsequent overload, where the ensuing overload is effective immediately.

Agreed. As it didn't appear to happen in reality, it's there simply to force progression. Early versions of the model with very high volumes of concrete crush were prone to arrest at the mechanical floor regions.

This is a surprising result to me for apparently the opposite reasons!

Do you mean the "d"/!"d" word ?

For this spreadsheet it would almost have to be arbitrary guesswork but, in general modeling case, I could not disagree more.

Definitely meant in terms of the spreadsheet.

I make the full energy loss available to subsequent calcs, so I use it in 'zero KE loss' sense...other than the velocity change.

I can understand keeping it around for accounting, but how do you apply it back in to the dynamics?

It's simply part of available KE, so for example, if the KE value 'used' to overcome structural resistance is equal to the energy loss 'gained' from the COM calc, then there is no change in subsequent velocity, and no change in collapse time. (That's never going to happen exactly, but hope it makes sense)
If the 'resistance' KE 'used' is less that that 'gained' through COM, subsequent velocity is not altered and the KE remainder is 'lost'.
If the resistance KE used is more than that gained through COM, subsequent velocity is reduced accordingly.
All based on remaining 'output' KE.

I believe one of the shortcomings of the NIST simulations is the unprecedented size and complexity of the models without having the basis of prior comparable modeling projects, or at least a long chain of substantial, incrementally scaled prototypes. If their simulations are correct by any measure, it's either unprecedented genius or pure dumb luck.

The ExtremeLoading system has apparently been used within the demolition industry for a couple of years now, at least in terms of proof of comparison. It's not FEA as such, but Applied Element Method (AEM). NIST did not use such a system.
Did you watch the video ?
http://www.youtube.com/watch?v=LPsW9cVPMgw

Thank you! Big bumps at the mechanical floors.

Would another which separates core from not core be useful ?
Image
http://femr2.ucoz.com/photo/1-0-71-3
femr2
 
Posts: 1010
Joined: Thu Apr 23, 2009 12:08 am
Location: UK

Re: WTCCS Spreadsheet Crush-Down Model

Postby OneWhiteEye on Sun Jun 14, 2009 4:58 am

femr2 wrote:I should make it clear that the additional 3.5 storeys is achieved in ~0.3s, so it's a very marginal difference. That's why I'm pushing for the inclusion for some forms of resistance

Alright then, bring on the resistance! Best figures you think you have all the way around, but start at 92 with s=0.2, let's see how far off it is.

Could you clarify how resistance would be reduced and less energy dissipated?

The upper block south quarter, as a triangular wedge, versus the floor connections to the south perimeter wall?

Such a wedge represents a small cross-sectional area at the tip, even if accordioned column ends and compressed debris, but has behind it a tremendous amount of inertia, provided that wedge can stay mostly intact for a period of time. This sort of arrangement can come to exist through ablation of the inner portion towards the core in conjunction with tipping, eroding the weak floor areas in the upper block. The tip of the wedge, if inside the lower perimeter, can exert tremendous pressure on a confined, vulnerable region. Floors are pushed downward, yes, but not exclusively, and only a small area near the edge. Perimeter pushed mostly outward.

Speculation.

But not at all implausible. Evidence of a relatively intact hat truss exists on the one hand, by way of antenna motion in the first half of collapse, and a obliterated northern half because of the spire, on the other. The tilt gives core intrusion on the diagonal, the squarish block means better retention by the perimeter than the prior example of WTC2, add some wastage of brittle unsupported floors experiencing maximum engagement and you've got a recipe for a balanced knife slicing off a peel. Three seconds would be about when this starts to take shape and, by 5 seconds, be gathering a serious head of steam. Peak at 7-8 seconds, then the angle becomes too severe. ???

Do you mean the "d"/!"d" word ?

Oh my, no. I expected the discrepancy to be greater, and then speculation like above would have to be more rigorously developed. It's too slow! A little bit. I thought this would be harder. A well-configured slab model tuned for full-floor axial impacts, no tilt, should NOT match this particular localized crush front at 5 seconds. In my opinion.

The ExtremeLoading system has apparently been used within the demolition industry for a couple of years now, at least in terms of proof of comparison. It's not FEA as such, but Applied Element Method (AEM). NIST did not use such a system.

Right tool for the job, and experience using it.

Did you watch the video ?

No, but I will.

Would another which separates core from not core be useful ?

Definitely!

Edit: chart added since I last looked - splendid!! Thank you.
OneWhiteEye
 
Posts: 1220
Joined: Sat Jul 05, 2008 9:40 pm

Re: WTCCS Spreadsheet Crush-Down Model

Postby femr2 on Sun Jun 14, 2009 11:51 am

OneWhiteEye wrote:Alright then, bring on the resistance! Best figures you think you have all the way around, but start at 92 with s=0.2, let's see how far off it is.

I don't really have 'best figures', just unclarified mechanisms. I hope the process of focussed discussion will result in refinement of the mechanisms, and a narrowing of the range of values used.

To kick the process off, I've set up the following scenario, which includes the most 'unverified' mechanisms, but does not have any more specific per-floor values set for elements such as mass loss.

Image
Simulation at 5s
YouTube Video
Image
Velocity Graph
Image
Mass Graph
Image
Kinetic Graph

Spreadsheet

Basic Parameters:
* 0.7492m Floor Thickness
* Initiation Floor 92
* Ramped Structural Support Energy (With core remnant modelling included)
* Impact KE based Concrete Crush (20% crushed. Not the same as ejected)
* Impact KE based Mass Loss (For Scenario, ranges from 1.6% at top to 20% at base. Not linear.)
* External Column Steel Mass Shed during impact
* Mass Loss Average Lateral Ejection Velocity set to 5m/s

First Stage Collapse Time = 30.6s
Second Stage Collapse Time = 32.8s

Initial Collapse Arrest Floor = 51

I'm sure there will be plenty of questions (or outrage :wink: ) so rather than go into detail now, I'll respond in detail.
femr2
 
Posts: 1010
Joined: Thu Apr 23, 2009 12:08 am
Location: UK

Re: WTCCS Spreadsheet Crush-Down Model

Postby OneWhiteEye on Sun Jun 14, 2009 5:38 pm

First Stage Collapse Time = 30.6s
Second Stage Collapse Time = 32.8s

Now I see why you're so anxious to get the energy sinks in!

I'm sure there will be plenty of questions (or outrage )...

Once I change my drawers, everything will be fine.

Must ponder this. I foresee a fairly long back and forth, after which one or both of us will come away changed.

Thanks.
OneWhiteEye
 
Posts: 1220
Joined: Sat Jul 05, 2008 9:40 pm

PreviousNext



Return to WTC1 and WTC2

Who is online

Users browsing this forum: No registered users and 0 guests