Portal Rigid Plane Frame Analysis Design Excel Sheet

Description

Name of File :- Portal Rigid Plane Frame Analysis Design Excel Sheet

“FRAME” is a spreadsheet program written in MS-Excel for the purpose of plane frame analysis of portal and gable rigid plane frames subjected to various types of loading. Specifically, the “stiffness matrix” method of analysis is used to determine the unknown joint displacements, support reactions, and member end forces. Individual frame members are also analyzed to determine the shears and intermediate moments. Plots of both the shear and moment diagrams are also produced. Also, the frame is drawn for visual confirmation of geometry/configuration.

This program is a workbook consisting of three (3) worksheets, described as follows:

Worksheet Name           Description
Doc                                This documentation sheet
Portal Frame                  Portal rigid plane frame analysis
Gable Frame                 Gable rigid plane frame analysis

Program Assumptions and Limitations:

1. This program uses the “stiffness matrix” method of analysis and four (4) following basic analysis assumptions:
   a. Members must be of constant cross section (E and I are constant for entire length).
   b. Deflections must not significantly alter the geometry of the problem.
   c. Stress must remain within the “elastic” region.
   d. Since this is a “first-order”, linear analysis, the effects of “P-delta” and shear deformation are not included.
   (However, significant effects due to shear deformation are limited to very short and deep members.)
2. Additional assumptions and features are as follows:
   a. Frame support joints may each be either fixed or pinned.
   b. Frame support joints may be at different levels (elevations).
   c. Columns must be vertical (cannot be sloped).
   d. For a portal frame, the top (roof) member may be flat or sloped in either direction.
3. A vertical load, horizontal load, and externally moment may be applied to any of the joints of the frame. These joint loads are to be applied in “global” axes directions. Note: Joint loads applied directly at supports are merely added directly to support reactions and are not reflected in member end force values.
4. On any individual member, this program will handle up to five (5) full uniform, partial uniform, triangular, or trapezoidal loads, up to ten (10) point loads, and up to four (4) externally applied moments. For vertical members, distributed loads and point loads are input in a “X-Global” sence of direction. For flat or sloped top (roof) members, distributed loads may be applied global over actual member length or applied global over the “projected” member length. Program designations are “Y-Global”, “Y-Projected”, “X-Global”, and “X-Projected”.
   For a flat top (roof) member of a portal frame, “Y-Global” and “Y-Projected” loads produce the same results.
   Uniformly distributed gravity (dead or live) load would be an example of a “Y-Global” distributed load on a sloped top (roof) member, while lateral uniformly distributed wind load on sloped top (roof) member would be an example of an “X-Projected” distributed load. A uniformly distributed load such as wind suction perpendicular (normal) to a sloped top (roof) member must be resolved into Y-Global and X-Global component values by user.
5. This program will calculate the member end reactions, the member end forces (axial, shear, and moment),
   the member maximum positive and negative moments (if applicable), and the joint displacements.
   The calculated values for the maximum moments are determined from dividing the member into fifty (50) equal segments with fifty-one (51) points, and including all of the point load and applied moment locations as well. (Note: the actual point of maximum moment occurs where the shear = 0, or passes through zero.)
6. The user is also given the ability to select an AISC W, S, C, MC, or HSS (rectangular tube) shape to aide in obtaining the required moment of inertia for input. (This facility is located off to the right of the main page.)
7. This program contains numerous “comment boxes” which contain a wide variety of information including explanations of input or output items, equations used, data tables, etc. (Note: presence of a “comment box” is denoted by a “red triangle” in the upper right-hand corner of a cell. Merely move the mouse pointer to the desired cell to view the contents of that particular “comment box”.)

Procedure for Stiffness Method of Frame Analysis:

1. Identifiy members and joints in frame
2. Specify near (start) joint and far (end) joint for each member in frame
3. Establish global coordinate system
4. Calculate fixed-end moments (FEM’s) and shears for each member due to applied member loads
5. Specify x, y, and z coding components (3 in all) at each joint as follows:
   a. Use lowest numbers to identify unknown joint displacements (for partioning overall matrix)
   b. Use remaining numbers to indentify known displacements
6. From the problem, establish the known displacements, Dk, and known external forces and reactions, Qk
7. Determine 6×6 stiffness matrix, k’, for each of the member expressed in global coordinates
   As an example, a member numbered with 1,2, and 3 at start and 4,5, and 6 at end would be:

where: A = A*E/L*lx^2+12*E*I/L^3*ly^2 Note: lx and ly are the “direction cosines”
B = (A*E/L-12*E*I/L^3)*lx*ly lx = (xj-xi)/L and ly = (yj-yi)/L
C = 6*E*I/L^2*ly (xi,yi) = joint start coordinates
D = A*E/L*ly^2+12*E*I/L^3*lx^2 (xj,yj) = joint end coordinates
E = 6*E*I/L^2*lx L = SQRT((xj-xi)^2+(yj-yi)^2)
F = 4*E*I/L
G = 2*E*I/L

• Merge individual member stiffness matrices into stiffness matrix, K, for entire frame
• Partition the structure stiffness matrix, K, as follows:
  End Forces Vector Partitioned Stiff. Matrix, K Displacements Vector
  Qk = K11 K12 * Du
  Qu K21 K22 Dk
  Expansion then leads to: Qk = K11 * Du + K12 * Dk (Eqn. 1)
  Qu = K21 * Du + K22 * Dk (Eqn. 2)
  where: Qk and Dk = known external loads and displacements (typ. Dk = 0, all)
  Qu and Du = unknown support reactions and displacements
  K11, K12, K21, K22 = submatrices of partitioned stiffness matrix, K
• Solve for unknown displacements, Du, from Equation 1:
  Du = (Qk – K12 * Dk) / K11 = Qk / K11 where: Dk = 0, all
• With the solved displacements, solve for unknown support reactions, Qu, from Equation 2
  Qu = K21 * Du + K22 * Dk = K21 * Du where: Dk = 0, all
• Solve for internal member end forces from: q = k’ * T * D
  where: T = 6×6 displacement transformation matrix as follows for each member:

and k’ is the member stiffness matrix (from above), and D is the displacements vector

• Superimpose member fixed-end moments (FEM’s) and shears with the frame analysis end forces to get final member end forcesReference: “Structural Analysis” – by Russel C. Hibbeler, Macmillan Publishing Company (1985), pages 441 to 497Formulas Used to Determine FEM’s and Shear and Moment in Individual Frame Members

  For Uniform or Distributed Loads:

  Loading functions for each uniform or distributed load evaluated at distance x = L from left end of member:
  FvL = -wb*(L-b-(L-e)) + -1/2*(we-wb)/(e-b)*((L-b)^2-(L-e)^2)+(we-wb)*(L-e)
  FmL = -wb/2*((L-b)^2-(L-e)^2) + -1/6*(we-wb)/(e-b)*((L-b)^3-(L-e)^3)+(we-wb)/2*(L-e)^2
  FqL = -wb/(6*E*I)*((L-b)^3-(L-e)^3) + -1/(24*E*I)*(we-wb)/(e-b)*((L-b)^4-(L-e)^4)+(we-wb)/(6*E*I)*(L-e)^3
  FDL = -wb/(24*E*I)*((L-b)^4-(L-e)^4) + -1/(120*E*I)*(we-wb)/(e-b)*((L-b)^5-(L-e)^5)+(we-wb)/(24*E*I)*(L-e)^4

  Loading functions for each uniform or distributed load evaluated at distance = x from left end of member:
  If x >= e:
  Fvx = -wb*(x-b-(x-e)) + -1/2*(we-wb)/(e-b)*((x-b)^2-(x-e)^2)+(we-wb)*(x-e)
  Fmx = -wb/2*((x-b)^2-(x-e)^2) + -1/6*(we-wb)/(e-b)*((x-b)^3-(x-e)^3)+(we-wb)/2*(x-e)^2
  Fqx = -wb/(6*E*I)*((x-b)^3-(x-e)^3) + -1/(24*E*I)*(we-wb)/(e-b)*((x-b)^4-(x-e)^4)+(we-wb)/(6*E*I)*(x-e)^3
  FDx = -wb/(24*E*I)*((x-b)^4-(x-e)^4) + -1/(120*E*I)*(we-wb)/(e-b)*((x-b)^5-(x-e)^5)+(we-wb)/(24*E*I)*(x-e)^4
  else if x >= b:
  Fvx = -wb*(x-b) + -1/2*(we-wb)/(e-b)*(x-b)^2 else: Fvx = 0
  Fmx = -wb/2*(x-b)^2 + -1/6*(we-wb)/(e-b)*(x-b)^3-(x-e)^3 else: Fmx = 0
  Fqx = -wb/(6*E*I)*(x-b)^3 + -1/(24*E*I)*(we-wb)/(e-b)*(x-b)^4 else: Fqx = 0
  FDx = -wb/(24*E*I)*(x-b)^4 + -1/(120*E*I)*(we-wb)/(e-b)*(x-b)^5 else: FDx = 0

  For Point Loads:

  Loading functions for each point load evaluated at distance x = L from left end of member:
  FvL = -P
  FmL = -P*(L-a)
  FqL = -P*(L-a)^2/(2*E*I)
  FDL = P*(L-a)^3/(6*E*I)

  Loading functions for each point load evaluated at distance = x from left end of member:
  If x > a:
  Fvx = -P else: Fvx = 0
  Fmx = -P*(x-a) else: Fmx = 0
  Fqx = -P*(x-a)^2/(2*E*I) else: Fqx = 0
  FDx = P*(x-a)^3/(6*E*I) else: FDx = 0

  For Applied Moments:

  Loading functions for each applied moment evaluated at distance x = L from left end of member:
  FvL = 0
  FmL = -M
  FqL = -M*(L-c)/(E*I)
  FDL = M*(L-c)^2/(2*E*I)

  Loading functions for each applied moment evaluated at distance = x from left end of member:
  If x >= c:
  Fvx = 0 else: Fvx = 0
  Fmx = -M else: Fmx = 0
  Fqx = -M*(x-c)/(E*I) else: Fqx = 0
  FDx = M*(x-c)^2/(2*E*I) else: FDx = 0

  Formulas Used to Determine FEM’s and Shear and Moment in Individual Frame Members (continued)

  Initial summation values at left end (x = 0) for shear, moment, slope, and deflection:

  Fixed beam: (for determining FEM’s in frame members)
  Vo = -12*E*I/L^3*S(FDL)-6*E*I/L^2*S(FqL)
  Mo = 6*E*I/L^2*S(FDL)+2*E*I/L*S(FqL)
  qo = 0
  Do = 0

  Simple beam:
  Vo = -1/L*S(FmL)
  Mo = 0
  qo = 1/L*S(FDL)+L/(6*E*I)*S(FmL)
  Do = 0

  Summations of shear, moment, slope, and deflection at distance = x from left end of member:

  Shear: Vx = Vo+S(Fvx)
  Moment: Mx = Mo+Vo*x+S(Fmx)
  Slope: qx = qo+Mo*x/(E*I)+Vo*x^2/(2*E*I)+S(Fqx)
  Deflection: Dx = -(Do-qo*x-Mo*x^2/(2*E*I)-Vo*x^3/(6*E*I)+S(FDx)

  Reference: “Modern Formulas for Statics and Dynamics, A Stress-and-Strain Approach”
  by Walter D. Pilkey and Pin Yu Chang, McGraw-Hill Book Company (1978)