7.5. Example: Eigenvalue Buckling Analysis (GUI Method)

In this sample problem, you will analyze the buckling of a bar with hinged ends.

7.5.1. Problem Description

Determine the critical buckling load of an axially loaded long slender bar of length with hinged ends. The bar has a cross-sectional height h, and area A. Only the upper half of the bar is modeled because of symmetry. The boundary conditions become free-fixed for the half-symmetry model. The moment of inertia of the bar is calculated as I = Ah2/12 = 0.0052083 in4.

7.5.2. Problem Specifications

The model is assumed to act only in the X-Y plane.

The following material properties are used:

E = 30 x 106 psi

The following geometric properties are used:

= 200 in
A = 0.25 in2
h = 0.5 in

Loading is a follows:

F = 1 lb.

7.5.3. Problem Sketch

Figure 7.3: Bar with Hinged Ends

Bar with Hinged Ends

7.5.3.1. Set the Analysis Title

After you enter the Mechanical APDL program, follow these steps to set the title.

  1. Select menu path Utility Menu> File> Change Title.

  2. Enter the text "Buckling of a Bar with Hinged Ends" and click on OK.

7.5.3.2. Define the Element Type

In this step, you define BEAM188 as the element type.

  1. Select menu path Main Menu> Preprocessor> Element Type> Add/Edit/Delete. The Element Types dialog box appears.

  2. Click Add. The Library of Element Types dialog box appears.

  3. In the scroll box on the left, click on "Structural Beam" to select it.

  4. In the scroll box on the right, click on "2 Node 188" to select it.

  5. Click OK. The Library of Element Types dialog box closes.

  6. Click Options in the Element Types dialog box.

  7. Select Element Behavior K3 : Cubic Form. Click OK.

  8. Click Close in the Element Types dialog box.

7.5.3.3. Define the Real Constants and Material Properties

  1. Select menu path Main Menu> Preprocessor> Sections> Beam> Common Sections. The BeamTool dialog box appears.

  2. Enter .5 for B and .5 for H.

  3. Click OK.

  4. Select menu path Main Menu> Preprocessor> Material Props> Material Models. The Define Material Model Behavior dialog box appears.

  5. In the Material Models Available window, double-click on the following options: Structural, Linear, Elastic, Isotropic. A dialog box appears.

  6. Enter 30e6 for EX (Young's modulus), and click on OK. Material Model Number 1 appears in the Material Models Defined window on the left.

  7. Select menu path Material> Exit to remove the Define Material Model Behavior dialog box.

7.5.3.4. Define Nodes and Elements

  1. Select menu path Main Menu> Preprocessor> Modeling> Create> Nodes> In Active CS. The Create Nodes in Active Coordinate System dialog box appears.

  2. Enter 1 for node number.

  3. Click Apply. Node location defaults to 0,0,0.

  4. Enter 11 for node number.

  5. Enter 0,100,0 for the X, Y, Z coordinates.

  6. Click OK. The two nodes appear in the Graphics window.


    Note:  The triad, by default, hides the node number for node 1. To turn the triad off, select menu path Utility Menu> PlotCtrls> Window Controls> Window Options and select the "Not Shown" option for Location of triad. Then click OK to close the dialog box.


  7. Select menu path Main Menu> Preprocessor> Modeling> Create> Nodes> Fill between Nds. The Fill between Nds picking menu appears.

  8. Click on node 1, then 11, and click on OK. The Create Nodes Between 2 Nodes dialog box appears.

  9. Click OK to accept the settings (fill between nodes 1 and 11, and number of nodes to fill 9).

  10. Select menu path Main Menu> Preprocessor> Modeling> Create> Elements> Auto Numbered> Thru Nodes. The Elements from Nodes picking menu appears.

  11. Click on nodes 1 and 2, then click on OK.

  12. Select menu path Main Menu> Preprocessor> Modeling> Copy> Elements> Auto Numbered. The Copy Elems Auto-Num picking menu appears.

  13. Click Pick All. The Copy Elements (Automatically-Numbered) dialog box appears.

  14. Enter 10 for total number of copies and enter 1 for node number increment.

  15. Click OK. The remaining elements appear in the Graphics window.

7.5.3.5. Define the Boundary Conditions

  1. Select menu path Main Menu> Solution> Unabridged Menu> Analysis Type> New Analysis. The New Analysis dialog box appears.

  2. Click OK to accept the default of "Static."

  3. Select menu path Main Menu> Solution> Analysis Type> Analysis Options. The Static or Steady-State Analysis dialog box appears.

  4. To include prestress effects, click the check box to "YES" for "[PSTRES] Incl prestess effects?"

  5. Click OK.

  6. Select menu path Main Menu> Solution> Define Loads> Apply> Structural> Displacement> On Nodes. The Apply U,ROT on Nodes picking menu appears.

  7. Click on node 1 in the Graphics window, then click on OK in the picking menu. The Apply U,ROT on Nodes dialog box appears.

  8. Click on "UY" and “ROTZ” to select them, and click on OK.

  9. Select menu path Main Menu> Solution> Define Loads> Apply> Structural> Displacements> On Nodes. The Apply U,ROT on Nodes picking menu appears.

  10. Click on node 11 in the Graphics window, then click on OK in the picking menu. The Apply U,ROT on Nodes dialog box appears.

  11. Click on "UX" to select it, and click on OK.

  12. Select menu path Main Menu> Solution> Define Loads> Apply> Structural> Force/Moment> On Nodes. The Apply F/M on Nodes picking menu appears.

  13. Click on node 11, then click OK. The Apply F/M on Nodes dialog box appears.

  14. In the scroll box for Direction of force/mom, scroll to "FY" to select it.

  15. Enter -1 for the force/moment value, and click on OK. The force symbol appears in the Graphics window.

  16. Select menu path Main Menu> Solution> Define Loads> Apply> Structural> Displacement> Symmetry B.C.> On Nodes. The Apply SYMM on Nodes box dialog appears.

  17. In the scroll box for Norml symm surface is normal to, scroll to “z-axis” and click on OK.

7.5.3.6. Solve the Static Analysis

  1. Select menu path Main Menu> Solution> Solve> Current LS.

  2. Carefully review the information in the status window, and click on Close.

  3. Click OK in the Solve Current Load Step dialog box to begin the solution.

  4. Click Close in the Information window when the solution is finished.

7.5.3.7. Solve the Buckling Analysis

  1. Select menu path Main Menu> Solution> Analysis Type> New Analysis.

  2. In the New Analysis dialog box, click the "Eigen Buckling" option on, then click on OK.


    Note:  Click Close in the Warning window if the following warning appears: Changing the analysis type is only valid within the first load step. Pressing OK will cause you to exit and reenter SOLUTION. This will reset the load step count to 1.


  3. Select menu path Main Menu> Solution> Analysis Type> Analysis Options. The Eigenvalue Buckling Options dialog box appears.

  4. Enter 1 for number of modes to extract.

  5. Click OK.

  6. Select menu path Main Menu> Solution> Load Step Opts> ExpansionPass> Single Expand> Expand Modes.

  7. Enter 1 for number of modes to expand, and click on OK.

  8. Select menu path Main Menu> Solution> Solve> Current LS.

  9. Carefully review the information in the status window, and click on Close.

  10. Click OK in the Solve Current Load Step dialog box to begin the solution.

  11. Click Close in the Information window when the solution is finished.

7.5.3.8. Review the Results

  1. Select menu path Main Menu> General Postproc> Read Results> First Set.

  2. Select menu path Main Menu> General Postproc> Plot Results> Deformed Shape. The Plot Deformed Shape dialog box appears.

  3. Click the "Def + undeformed" option on. Click OK. The deformed and undeformed shapes appear in the Graphics window as seen in Figure 7.4: Buckling Analysis: Deformed and Undeformed Shapes.

Figure 7.4: Buckling Analysis: Deformed and Undeformed Shapes

Buckling Analysis: Deformed and Undeformed Shapes

7.5.3.9. Exit Mechanical APDL

  1. In the Toolbar, click Quit.

  2. Select the save option you want and click OK.