9 - 1 1 R e s e a r c h papers

Analysis of the Mass and Potential Energy

of World Trade Center Tower 1

 

Gregory H. Urich

B.S. Electrical and Computer Engineering

1                     Abstract

The mass and potential energy of one of the Twin Towers is calculated based on available data. The mass for each floor is established based on floor types, documented design loads, and estimated in-service live loads. The calculated mass of 288,100 metric tons (317,500 short tons) is found to correspond with two other comparable structures in terms of mass per unit floor area, NIST’s SAP2000 model, and the reported amount of recovered debris. The calculated mass refutes the popular notion that the building weighed 500,000 tons.


2                     Introduction

 

In the aftermath of the World Trade Center Disaster, the Federal Emergency Management Agency (FEMA) and the National Institute of Standards and Technology (NIST) conducted investigations of considerable scope regarding building performance and the collapse of the World Trade Center Towers. The FEMA report described a pancake-type progressive floor collapse scenario causing the removal of lateral support on several floors leading to buckling in unshored columns which were weakened by fire and partially damaged by aircraft impact.3 The FEMA report was not rigorous and the conclusions regarding collapse initiation and progressive collapse can only be considered to be an educated guess by the investigators. The more rigorous NIST reports described aircraft impact damage and collapse initiation based on forensic evidence and computer simulations. Unlike scientific articles, the NIST investigation reports do not provide enough information to be able to reproduce the models or any results derived from the models.

 

Analyses from independent researchers regarding aircraft impact damage and collapse scenarios have appeared during and after the official investigations. Earlier analyses were severely limited by a lack of information and were overly simplified. Later analyses have been more substantial, but as seen in Bazant et al. (2007)10, the mass and potential energy are probably grossly overestimated.

 

Due to these limitations, many have questioned both the government’s specific account of collapse initiation and the general theory of gravity-driven, progressive collapse. These questions can only be answered by better modeling and a truly scientific approach. To be valid, further analyses and models must be based on the correct mass and mass distribution throughout the building.

 

The purpose of this paper is to establish a substantiated mass, mass distribution and potential energy in World Trade Center Tower 1 (north tower) within a reasonable margin of error. Here, the NIST “Federal Building and Fire Safety Investigation of the World Trade Center Disaster” (called NCSTAR) documents provide a wealth of information regarding the structural design, dimensions, building materials, contemporaneous building codes and an approach to modeling.

2.1                  Popular numbers

Many references can be found with different values for the mass of and the amount of potential energy stored in the WTC twin towers. A number of references are shown in Table 1 below. None of these references provide any data or calculation method on which the mass and potential energy are based.

 

Table 1: Different values for mass and potential energy

 Source

Mass

Potential Energy

Ashley 5

500,000 tons

 

Bazant et al. 10

576,000 tons *

 

Eagar and Musso 12

500,000 tons

 

Hamburger, et al. (FEMA) 4

 

> 4 E+11 J

Sunder et al. (NIST NCSTAR1) 8

> 250,000 tons

 

Tyson 2

500,000 tons

 

 

* The Bazant et al. number is calculated here based on the following:

 

“Near the top, the specific mass (mass per unit height) µ = 1.02 × 106 kg/m. In view of proportionality to the cross section area of columns, µ = 1.05 × 106 kg/m at the impact level (81st floor) of South Tower. Generally, we assume that µ(z) = k0ek2z + k1 (where k0, k1, k2 = constants), with a smooth transition at the 81st floor to a linear variation all the way down (precise data on µ(z) are unavailable). The condition that ∫0 H µ(z)dz be equal to the total mass of tower (known to be almost 500,000 tons) gives µ = 1.46 × 106 kg/m at the base.” 10

 

Since µ(z) is unknown we can approximate the value for floors 82-110 using a linear variation from the value at floor 81 to the value at floor 110 (29 floors) and the proportion of the height for those floors. The height of WTC1 from the base to the roof is 437.69 m. The total number of floors is 116. µ(z)avg81-110  = 1.035 × 106 kg/m.  µ(z)avgB6-81 = 1.2475 × 106 kg/m.

 

Mass82-110 = µ(z)avg81-110 x (29/116) x h = 113.3 × 106 kg

 

MassB6-81 = µ(z)avgB6-81 x (87/116) x h = 409.5 × 106 kg

 

The total mass is then 522.8 × 106 kg or, converting to short tons, 576,000 tons. Bazant et al. most likely assumed metric tons for the popular 500,000 ton number but that doesn’t explain why we get 522.8 × 106 kg. The maximum error of using the linear approximation instead of the exponential equation is less than 2 × 106 kg. If Bazant et al. used the nominal height of the building (414.63 m from the concourse level to the roof) the result would be 493.9 × 106 kg which corresponds better to the statement “known to be almost 500,000 tons” assuming metric tons.

2.2                  Original Design

 

A number of original design documents are provided in NCSTAR1-1 and NCSTAR1-1A. NIST NCSTAR1-1A (p. 5)9 presents definitions from the original design as follows:

 

1.      “Floor inside of core”. That part of the floor bounded by the outside faces of columns 501, 508, 1001 and 1008.

2.      “Floor outside of core”. That part of the floor between the outside walls and the “Floor inside of core”.

3.      “Code live load”. The load specified in the New York Building Code for a given occupancy.

4.      “Live load for floor design”. The actual live load used for the design of the parts of the floor which load may not be less than the “Code live load”, and may be reduced for tributary areas as defined in “Live load reduction”.

5.      “Live load for column design”. The code live load, reduced as defined in “Live load reduction” for columns.

6.      “Construction dead load”. The weight of the bare structure (i.e. the slab and beam) used in design of unshored composite beams.

7.      “Construction live load”. The allowance for the weight of any equipment and/or forms which is not permanent and does not form part of the total load summation.

8.      “Superimposed dead load”. The weight of ceilings, floor finish, walls or partitions of known location, mechanical and electrical equipment and similar items not included in the “Superimposed live load” or “Construction dead load”.

9.      “Dead load”. The sum of items 6 and 8 above.

10.  “Superimposed live load”. The weight of the design live load, based on occupancy, plus the weight of partitions if their location is subject to change.

 

Values for construction dead load (CDL), superimposed dead load (SDL) and superimposed live load (SLL) are also given in the design documents presented in NCSTAR1-1A for some of the different types of floors within the building, inside and outside the core. CDLs include steel used in floors such and beams, trusses, deck and concrete reinforcement.

2.3                  Amount of Steel

NIST gives the total the weight of structural steel in the two WTC towers as 200,000 tons.11 NIST describes steel contracts in NCSTAR1-3 (p.16), and the values are shown in Table 2 below.3 These contracts do not include trusses outside the core, steel deck, concrete reinforcements or grillages.

Table 2: Weight of steel from supplier contracts

Structural component

Weight (short tons)

Weight per tower (short tons)

external columns w/ spandrels

55 800

27 900

rolled core columns and beams

25 900

12 950

bifurcation columns

6 800

3 400

external box columns

13 600

6 800

core box below floor 9

13 000

6 500

core box above floor 9

31 000

15 500

slab supports below grade

12 000

6 000

total

158 100

79 050

 

2.4                  NIST Reference Models

 

In the NCSTAR1-2 series, NIST presents the methods used for developing the reference structural models the WTC towers. These models were used to assess the towers’ ability to withstand gravity and wind loads and to establish the reserve capacity in the structures to withstand unanticipated events. According to NIST:

 

“The reference models included the following: Two global models of the primary structural components and systems for each of the two towers (and) two models, one of a typical truss-framed floor (tenant floor) and one of a typical beam-framed floor (mechanical floor), within the impact and fire regions. All reference models were linearly elastic and three-dimensional, and were developed using the Computers and Structures, Inc. SAP2000 software. SAP2000 is a commercial finite-element software package that is customarily used for the analysis and design of structures.” 7

 

The databases for the reference models were developed based on original structural drawings. The databases were reviewed and checked against the original drawing books. According to NIST:

 

“The original structural drawings of the WTC Towers were issued in two main formats: (1) Large-size drawing sheets containing plan and elevation information, and (2) Smaller book-sized drawings containing details and tabulated information of cross-sectional dimensions and material properties. The larger-sized drawings referred to the structural drawing books in their notes, section and details. The structural databases, developed in Microsoft Excel file format, were generated from these drawing books and included modifications made after construction. The databases were generated for use in the development of reference global models of the towers.” 7

 

None of the original structural drawings were released by NIST. However, the larger drawing sheets for WTC-1 (north tower) were leaked subsequently to the general public and are generally available.17 The smaller drawing books still have not been made public.

2.5                  NIST’s “Tower and Aircraft Impact Models”

 

NIST describes the “Tower and Aircraft Impact Models” in NIST NCSTAR 1-2. These models were developed using the LS-DYNA 2003 software package.

 

“The WTC models for the impact analysis required considerably greater sophistication and detail than was required for the reference models described in Chapter 2. The reference models provided a basis for the more detailed models required for the impact simulations. The impact models of the towers, which utilized the structural databases described in Chapter 2 (see also NIST NCSTAR 1-2A), included the following refinements…” 7 (p. 93)

 

The loading of the structure for the analysis was determined by NIST as follows:

 

“The densities of the tower components (workstations and gypsum walls) were scaled by the appropriate ratios to obtain the desired distribution of live loads in the core and truss floor areas. The densities of all the remaining tower structural components were scaled proportionately to obtain the desired superimposed dead loads. These additional loads were important for obtaining an accurate mass distribution in the towers and inertial effects in the impact response. The in-service live load used was assumed to be 25% of the design live load on the floors inside and outside the core. The in-service live load was selected based on a survey of live loads in office buildings (Culver 1976) and on engineering judgment.” 7 (p. 106)

 

NIST NCSTAR 1-2B (p. 53) gives an SDL (36.2 psf) which is in fact applied to the structural components (columns).13 The SDL mass being applied to columns, is not a problem when calculating the mass. However, the impact analysis must be significantly affected by reducing the probability of debris coming into contact with core contents.  The effect is that impacting debris has a free shot at core structural members and is more likely to pass all the way through the core. It is unclear if the partitions are included in the SDL, SLL or both.

 

NIST NCSTAR 1-2B (p. 53) gives a summary of superimposed dead loads and live loads and floor areas to which they are applied.13 The values are shown in table 3 below.

 

Table 3: Summary of superimposed dead loads and live loads

 

Area (sq ft)

Weighting (psf)

Core Dead Load (SDL)

8,694

36.2

Outer Dead Load (SDL)

31,257

11.5

Core Live Load

8,694

19.7

Outer Live Load

31,257

16.2

 


3                     Method

The mass for the building is calculated on a floor by floor basis based on information in the NIST reports and the architectural drawings. In some cases there are deviations from NIST values and motivations for alternative values are described. In cases where there is not enough information in the NIST reports, dimensions or materials are used from similar areas of the building. As described in the introduction above, the design documentation for WTC1 has the structural loads divided into construction dead loads, superimposed dead loads, and superimposed live loads. These divisions are also used here.

 

3.1                  Floor Areas

 

According to NIST, the floor areas inside the core and outside the core are 8,694 sq ft and 31,257 sq ft respectively (see Table 3). However, the architectural drawings give the distance between the center of the external columns on one side to the center of the external columns on the other side as 207’-8”.17 NIST gives the width of the external column flanges as 13.5” and the spandrel thickness as 5/8”. Together these are roughly 14” contributing approximately 7” on each side to the 207’-8” dimension. Thus the overall floor dimensions must be 206’-6” x 206’-6” with a gross floor area of 42,642 sq ft. The outer dimensions of the core were 137’ x 85’ giving a gross core area of 11,745 sq ft. Thus the floor area outside the core is 30,897 sq ft. It may be that NIST subtracted the areas taken up by core columns, elevator shafts and utility shafts in the core area, which would account for the difference of roughly 25%. Generally in this analysis, the floor areas used inside the core and outside the core are 11,745 sq ft and 30,897 sq ft respectively.

 

For the purposes of establishing CDLs in the core, the floor areas inside the core were adjusted to account for empty space due to elevator and utility shafts. The actual floor areas were approximated by sampling a number of representative floors using the architectural drawings.17 Two sizes of elevators predominated and the other shafts were split into three groups: small, medium and large. The areas for the shafts in each group were established by taking the dimensions of all shafts on floors 11-16 from the architectural drawings (core plans), grouping them, and taking the average size for each group.17 Elevators and shafts were then counted on the representative floors and grouped by size. Elevators and shafts on average take up 41% of the core floor area. The sampled floors, number of elevators and shafts, area with no floor, and the percentage of empty space in the core are shown in Table 4. See Diagram 1 for examples of elevators and shafts.

 


Table 4: Elevators and shafts on representative floors

floor

count elev. 1 89.5 sq ft

count elev. 2  203.8 sq ft

count shaft 1 10.5  sq ft

count shaft 2 27.3  sq ft

count shaft 3 48.1  sq ft

area w/ no floor sq ft

% core

6

27

23

12

16

0

7 231

0,62

14

27

23

9

8

10

7 462

0,64

26

15

23

10

10

10

6 456

0,55

33

15

23

10

10

10

6 456

0,55

41

3

23

10

10

10

5 385

0,46

50

26

14

11

7

6

5 508

0,47

66

14

12

12

8

7

4 153

0,35

77

3

14

11

5

6

3 401

0,29

83

24

2

8

5

5

2 974

0,25

94

19

2

11

5

5

2 559

0,22

105

7

2

8

4

5

1 429

0,12

average

 

 

 

 

 

4 819

0,41

 

 

 

Diagram 1: Example of architectural drawings - core plan floors 11-16. (Colored areas with number designations are examples of the groups described above.)


3.2                  Floor Types

A table of floors and diagrams of 15 different floor framing types are given in NIST NCSTAR 1-2A, Appendix G (p192-196).14 The table shows which type of floor framing was used for each particular floor. The diagrams show how the different types of framing (i.e. truss or beam) were used in different floor types. The approximate percentages of truss versus beam areas can easily be deduced from the diagrams. Unfortunately, there is no indication of concrete types or thicknesses. For the purposes of this analysis, floors are divided up into normal, mechanical, special and sublevel floors.

3.2.1                   Normal Floors

All floors are considered normal unless they are mechanical floors, sublevels, or special floors as described below. The floor numbers for normal floors are 10-40, 44-74 and 78-106. These floors, which are predominantly offices and related areas, comprise eleven floor types (type 1-11) that are predominantly truss framed. Some of these types have sections of beam framed floor and two types have heavier angles or have reinforced trusses. All of these floor types are treated as type 1 (100% truss framed) to simplify the calculation of mass. The total error induced by this simplification is less than 1/1000 and can be calculated as follows:

 

Err% = AAvg% x S% x (B-T)/B% = 0.073%

 

Where AAvg is the average proportion of beam area (1.45%), S = the floor frame steel proportion of the total mass of the building (approx. 10%), T is the truss design construction dead load (10 psf) 9 and B is the design construction dead load (20 psf) 9 for beam framed floors. See table 5 for calculation of AAvg.

 

Table 5: Floor types, count and beam framed area for calculation of average beam area.

Floor type

Beam area/floor area

count

1

0

74

2

0

3

3

0,0338

4

4

0,2772

1

5

0,2017

1

6

0

1

7

0,1689

1

8

0,1014

1

9

0,0338

4

*10

0,5

1

11

0

0

total

1,3168

91

avg %

0,0145

 

 

* Note: Floor 106 is type 10 which has reinforced trusses. The reinforced trusses are assumed to be heavier than normal trusses and lighter than beam frames. Thus the floor is given as being 50% beam framed to account for the extra weight for the purpose of calculating the error due to simplification.

 


3.2.2                   Mechanical floors

 

The mechanical floors are 7-9, 41-43, 75-77, and 108-110, which are all beam framed floors. In each group of three floors, the upper and lower floors are type 12 and the middle floor is type 13 (mechanical mezzanine). The mechanical mezzanines were 50% open (no floor) outside the core so the floor area is 15,448.5 sq ft.

 

3.2.3                   Sublevels

 

Sublevel floors were beam framed floors, designated B1-B6, and are type 14. As seen in Table 2, 6000 tons of steel were used for slab support below grade. There is a minor discrepancy between the NIST documentation and the architectural drawings. In the architectural drawings, the floor below floor 1 is called the “Service Level” and the five floors below are named B1-B5.

 

3.2.4                   Special floors

 

Special floors include the Concourse level (floor 1), Plaza level or mezzanine (floor 2), and the roof, which are beam framed floors. Floors 3-6 had no floors outside the core. The Concourse level which was a high pedestrian traffic area is type 14 and probably had stronger than normal floors. The Plaza level was type 15 and was partially open. Floor 107 was the restaurant “Windows on the World” and had beam-framed floors.

3.3                  Gravity Loads

 

3.3.1                   Foundation

 

The mass of the foundation provides no load on structural components other than itself and contributes a negligible amount to potential energy. The mass of the foundation is nonetheless approximated based on the film footage from the Port Authority of New York and New Jersey.1 Dimensions are established by comparison to objects of known size, i.e. humans. The total mass of the foundation is shown in Table 7.

 

The foundation for the core columns was comprised of steel reinforced concrete footers and steel grillages built up out of I-beams. One steel grillage is made up of 17 I-beams, each with approximate dimensions 0.75m x 0.2m x 2m and a plate thickness of around 0.03m. Each grillage also had a base plate for the core column with average approximate dimensions 1m x 1m x 0.1m. It is assumed that there is one grillage per core column. Using a density of 7.784 metric tons per cubic meter for the density of A36 steel, the total mass for the grillages is approximately 484 metric tons. Each grillage was placed on a concrete footer with approximate dimensions 2.5m x 2.5m x 2m. Using a density of 2.4 metric tons per cubic meter, the total mass for the concrete footers was approximately 1410 metric tons.

 

The foundation for the external columns was comprised of a continuous, steel reinforced, concrete footer and base plates ranging from 7 to 9 sq ft (approx. 0.74 m2).6 The thickness of the base plate is unknown but a thickness of 3 cm is assumed. Using a total number of 80 exterior columns (transition to 238 columns at 7th floor), the total mass of the base plates is approximately 14 metric tons. The concrete footer for the external columns had a perimeter of 252 m. The other dimensions of the footer are unknown but are approximated using 2 m for depth and 2 m for width. The total mass for the concrete footer was approximately 2420 metric tons.

 

Table 6: Mass of the foundation

Component

Mass

(short tons)

Mass

(metric tons)

Core steel grillage w/ base plate

513

466

Core concrete footer

1555

1410

External column steel base plates

15

14

External column concrete footer

2670

2420

Total mass foundation

4753

4310

 

 

3.3.2                   Amount of Core Column Steel

 

As described in the introduction, the steel contracts included 6,500 tons for core box columns below the 9th floor, 15,500 tons for core box columns above the 9th floor and 12,950 short tons for rolled columns and beams. The amount of steel attributed to rolled columns (wide flange shapes) is calculated in Appendix 2 as 3,268 short tons. Thus the total core column steel is 25,268 short tons.

3.3.3                   Variation of Core Column Steel

 

Core columns dimensions have been extracted from NIST’s SAP2000 model, which was released based on a Freedom of Information Act (FOIA) request. These dimensions are currently available on the internet.15 It can be seen in this data that the variation of core columns steel is non-linear in the areas from floor B6 to floor 7 and from floor 107 to the roof. There are also non-linear variations at the mechanical floors where the columns were somewhat heavier, but these are ignored. The variation of core column steel mass is shown in Table 7, which is based on calculations of core column steel per floor for selected floors (see Table 19 in Appendix 3).

 

Table 7: Variation of Core Column Steel

Floor range

Variation

Varies from (tons)

Varies to (tons)

B6-001

Linear

380.63

427.14

002-007

Individual floor

n.a.

427.14

008-053

Linear

427.14

181.57

054-106

Linear

181.57

30.91

107

Individual floor

n.a.

35.81

108

Individual floor

n.a.

41.42

109

Individual floor

n.a.

35.81

110

Individual floor

n.a.

35.81

111(roof)

Individual floor

n.a.

31.20

 

When the core column steel mass is varied in this manner, the total core column steel becomes 24,576 tons with 5,801 tons below floor 9. This amount of core column steel below floor 9 should be 6,500 tons according to the steel contra