ANALYSIS FOR THE STABILITY OF AN ARCH(PROGRAM FRAMEX.BAS)
The structure chosen is one that could have stability problems: a two hinged arch with a span of 100 feet and a rise of 15 ft. and a relatively light load. It is intended to show the importance of including nonlinear deflections and moments in the design of an arch. The analysis was made with the program FRAMEX.BAS previously described. The arches are on 20 ft. centers, and have a dead load of 20 psf and a live load of 15 psf on the slope (not on the horizontal projection). The loads at the nodes will the all be the same. Note that if the arch is curved, the loads should be concentrated at the nodes, otherwise there will be an erroneous moment at the nodes. The arch is divided into 12 sections as shown in the following sketch.
First we will make a preliminary design. The maximum moment in the arch will occur at the quarter point for the live load on that half of the span, or at the opposite quarter point. It can be demonstrated that the moment due to live load is WL2/64, where W is the live load, and L is the span of the arch. On a parabolic arch with a uniform load, there will be no moment due to the dead load, so we will assume for a circular arch, the moment is low. The horizontal thrust is, approximately, H=(D + L/2)xS²/8R where D is the dead load, L is the live load, S is the span, and R is the rise of the arch: H = (20x20+20x15/2)x100²/8x15 = 45,833 pounds = 44.83 kips The thrust at the quarter point is approximately equal to: H/cosA, where A is the angle with the horizontal. In this case it is 16.7 degrees and H = 50.46 kips The moment at the quarter point is: M = 20x15x100²/64 = 46,875 fp, =46.88 fk The section chosen for investigation is a W16x26, with A=7.68 si, I=301, and S= 38.4. Using these values, and the equation f=P/A +/ M/(S), the maximum stress is 21.22 ksi. Input for the program is shown below: 5000 DATA ARCH FOR WEB SITE  DL + 1/2 LL
5010 ' Parameters
The results from running this program are as follows:
Note the great increase in deflections and moments, and that the maximum moment is on the right half of the arch. This is an extreme case but illustrates the stability calculations for an arch. The final deflection, D(final), can be predicted from the initial deflection by assuming that the shape of the moment diagram is similar to the shape of the deflection diagram: D(final)=D(initial) + (D(initial)/M)D(final)P where P is the thrust at the section, and M is the moment. Solve for D(final): D(final) = D(initial)/(1  PD(initial)/M) Using the valuses obtained from the first pass: D(final) = 0.373/(141.42*0.373/76.71) = 0.467 ft This compares to 0.474 for the program. If you would like to play around with this file (It will require a Zip program:)
