Wednesday, 1 December 2010

Modal Analysis of a Cantilever Beam


Introduction

This tutorial was created using ANSYS 7.0 The purpose of this tutorial is to outline the steps required to do a simple modal analysis of the cantilever beam shown below.

Preprocessing: Defining the Problem

The simple cantilever beam is used in all of the Dynamic Analysis Tutorials. If you haven't created the model in ANSYS, please use the links below. Both the command line codes and the GUI commands are shown in the respective links.

Solution: Assigning Loads and Solving


  1. Define Analysis Type
  2. Solution > Analysis Type > New Analysis > Modal
    ANTYPE,2

  3. Set options for analysis type:

    • Select: Solution > Analysis Type > Analysis Options..The following window will appear

    • As shown, select the Subspace method and enter 5 in the 'No. of modes to extract'
    • Check the box beside 'Expand mode shapes' and enter 5 in the 'No. of modes to expand'
    • Click 'OK'Note that the default mode extraction method chosen is the Reduced Method. This is the fastest method as it reduces the system matrices to only consider the Master Degrees of Freedom (see below). The Subspace Method extracts modes for all DOF's. It is therefore more exact but, it also takes longer to compute (especially when the complex geometries).
    • The following window will then appear
      For a better understanding of these options see the Commands manual.
    • For this problem, we will use the default options so click on OK.
  4. Apply Constraints
  5. Solution > Define Loads > Apply > Structural > Displacement > On KeypointsFix Keypoint 1 (ie all DOFs constrained).
  6. Solve the System
  7. Solution > Solve > Current LS
    SOLVE

Postprocessing: Viewing the Results


  1. Verify extracted modes against theoretical predictions

    • Select: General Postproc > Results Summary...The following window will appear

      The following table compares the mode frequencies in Hz predicted by theory and ANSYS.

      ModeTheoryANSYSPercent Error
      18.3118.3000.1
      251.9452.010.2
      3145.68145.640.0
      4285.69285.510.0
      5472.22472.540.1

      Note: To obtain accurate higher mode frequencies, this mesh would have to be refined even more (i.e. instead of 10 elements, we would have to model the cantilever using 15 or more elements depending upon the highest mode frequency of interest).
  2. View Mode Shapes

    • Select: General Postproc > Read Results > First SetThis selects the results for the first mode shape
    • Select General Postproc > Plot Results > Deformed shape . Select 'Def + undef edge'The first mode shape will now appear in the graphics window.
    • To view the next mode shape, select General Postproc > Read Results > Next Set . As above choose General Postproc > Plot Results > Deformed shape . Select 'Def + undef edge'.
    • The first four mode shapes should look like the following:
  3. Animate Mode Shapes

    • Select Utility Menu (Menu at the top) > Plot Ctrls > Animate > Mode ShapeThe following window will appear

    • Keep the default setting and click 'OK'
    • The animated mode shapes are shown below.

      • Mode 1
      • Mode 2
      • Mode 3
      • Mode 4

Using the Reduced Method for Modal Analysis

This method employs the use of Master Degrees of Freedom. These are degrees of freedom that govern the dynamic characteristics of a structure. For example, the Master Degrees of Freedom for the bending modes of cantilever beam are
For this option, a detailed understanding of the dynamic behavior of a structure is required. However, going this route means a smaller (reduced) stiffness matrix, and thus faster calculations.
The steps for using this option are quite simple.

  • Instead of specifying the Subspace method, select the Reduced method and specify 5 modes for extraction.
  • Complete the window as shown below
    Note:For this example both the number of modes and frequency range was specified. ANSYS then extracts the minimum number of modes between the two.
  • Select Solution > Master DOF > User Selected > Define
  • When prompted, select all nodes except the left most node (fixed).The following window will appear:

  • Select UY as the 1st degree of freedom (shown above).The same constraints are used as above.
The following table compares the mode frequencies in Hz predicted by theory and ANSYS (Reduced).

    ModeTheoryANSYSPercent Error
    18.3118.3000.1
    251.9452.010.1
    3145.68145.660.0
    4285.69285.710.0
    5472.22473.660.3
As you can see, the error does not change significantly. However, for more complex structures, larger errors would be expected using the reduced method.

ANSYS Command Listing

FINISH
/CLEAR

/TITLE, Dynamic Analysis
/PREP7

K,1,0,0 ! Enter keypoints
K,2,1,0

L,1,2 ! Create line

ET,1,BEAM3 ! Element type

R,1,0.0001,8.33e-10,0.01 ! Real Const: area,I,height

MP,EX,1,2.068e11 ! Young's modulus
MP,PRXY,1,0.33 ! Poisson's ratio
MP,DENS,1,7830 ! Density

LESIZE,ALL,,,10 ! Element size
LMESH,1 ! Mesh line

FINISH
/SOLU

ANTYPE,2 ! Modal analysis
MODOPT,SUBSP,5 ! Subspace, 5 modes
EQSLV,FRONT ! Frontal solver
MXPAND,5 ! Expand 5 modes

DK,1,ALL ! Constrain keypoint one

SOLVE
FINISH

/POST1 ! List solutions
SET,LIST

SET,FIRST
PLDISP ! Display first mode shape

ANMODE,10,0.5, ,0 ! Animate mode shape

No comments:

Post a Comment