OK. There are your mistakes. (I assume you are using floor 97 with 33,000 tonnes total mass above as usual.)
In accord with the correct load distribution calculated by Urich [2007] (which concurs with S. W. Banovic, T. Foecke, W.E. Luecke, et al. “The role of metallurgy in the NIST investigation of the World Trade Center towers collapse”, JOM, vol. 59, no. 11, pp. 22-29, November 2007.)
the values you should be using are:
The mass acting on the core is closer to 19,500 tonnes.
The total core demand is 1,9081E+08 N.
Correct calculation of buckling stress here using the actual properties of the steel as opposed to your assumptions.
The load capacity of the core is 4,3603E+08 N.
Core DCR = 44% (FoS = 2.29) prior to aircraft impact damage.
Core DCR = 51% (FoS = 1.96) after aircraft impact damage.
Greg, as you know, I have been working on a computer program to calculate stresses in a purely elastic, axial collision of a WTC top (floors 93 to 111) with the WTC base (floors 92 and down). In the process, I calculate a spring constant for each floor, using the mass of columns, only (no spandrels, e.g.). I also calculate the total mass for each floor, using the assumption that the mass remaining when you subtract the absolute total mass from the column masses, is to be evenly divided between 111 floors. (So, I ignore differences in spandrel weights, service floors, etc., and I completely ignore the basement floors).
Also, to repeat, I haven't carefully double checked everything, yet. Caveat emptor.
Even so, I read this thread and got curious as to how theoretical maximum load that is within the elastic limit compares to the static weight. I just added a few lines of code,
Output:
Spring Constant at floor 92: 3965840390490.68
Total Mass floors 93 - 111: 43824141.6906478
Therefore, static load from WTC top is 429,476,588.568348
Force to get .002 compression at floor 92 is 30,695,604,622.3978
Do you have any comments on this wide discrepancy? Somehow, taking a normal eccentricity into account doesn't seem like it should make a column over 10X as weak, but then again, I'm not an engineer.
Whoa, after I posted the above, I noticed that my spring constant calcs were using total mass, instead of column mass. Here are the corrected figures:
Static load from WTC top is 429,476,588.568348
Force to get .002 compression at floor 92 is 1,971,244,221.96645
===================================
Spring constant calcs:
// compute spring constants; from F = (YA / L) * Delta_L = (YV / L^2) * Delta_L ,
// k = (Y(m/rho) / L^2 ) where Y is Young's modulus = 200 GPa
Last edited: