Carnegie Mellon

Mechanical Engineering

Self-paced learning on the Web
FEM/ANSYS

 

V1 Turbine Blade
Home Course Info Problems Test Problems Students Reference

Vibration #1: Modal Analysis of a Turbine

 

Introduction: In the area of dynamics and vibrations the natural frequencies of a structure is of great importance to determine whether a structure can withstand excitation from the surroundings. In this example, we will learn to model a turbine and then determine its first few natural frequencies.

Physical Problem: To determine the natural frequencies of the turbine shown in the figure. Modal analysis means the calculation of the natural frequencies of a mechanical system. It also involves the calculation of the mode shapes.

Problem Description:

 

bullet

We will model the turbine as a disk with blades fixed on it ('blisk'=bladed disk). The inner radius of the hub is 10 cm, outer radius is 40 cm, blade length is 20 cm, blade width is 5 cm, and the thickness is 2.5 mm

bullet

Material: Assume the structure is made of steel with modulus of elasticity E=210 GPa and has a Poisson ratio of 0.3 and density of 7.21e3 kg/cubic meter.

bullet

Boundary conditions: The blisk is fixed around the inner diameter of the disk.

bullet

Loading: The blisk is not loaded.

bullet

Objective:
bullet

To determine first three family of modes.

bullet

To animate the mode shape of the first 3 modes.

bullet

You are required to hand in print outs for the above. You don't have to hand in the animation files but you will have to give at least 3 captured frames of the animation for each of the three modes.

bullet

Figure:

 

STARTING ANSYS:

 

bullet

Click on ANSYS 6.1 in the programs menu.

bullet

Select Interactive.

bullet

The following menu that comes up. Enter the working directory. All your files will be stored in this directory. Also enter 64 for Total Workspace and 32 for Database.

bullet

Click on Run.

 


   
MODELING THE STRUCTURE:

 

bullet

We will model one quarter of the blisk and then reflect it to create the complete blisk.3

bullet

Go to the ANSYS Utility Menu.
bullet

Click Workplane>WP Settings.

bullet

The following window comes up:

 

 

bullet

Check the Cartesian and Grid Only buttons

bullet

Enter the values shown in the figure above.

 

bullet

The following is the quarter blisk we will model first:

 

 

bullet

Go to the ANSYS Main Menu Preprocessor>Modeling>Create>Areas>Rectangle>By 2 corners.

bullet

Select the two corners for the horizontal rectangle and click Apply. Remember the rectangles (blades) have a thickness half the actual thickness since we are modeling only a quarter of the blisk.

 

 

bullet

Now similarly create the vertical rectangle.

 

 

bullet

Now we will create the quarter disk. Go to Preprocessor>Modeling>Create>Areas>Circle>Partial Annulus. The following window comes up:

 

 

bullet

Enter the values as shown and click OK. The model looks like the one below:

 

 

bullet

Now we will reflect the areas we have created about the YZ plane and then all the areas about the XZ plane. Go to Preprocessor>Modeling>Reflect>Areas. Click on "Pick All". the following window comes up:

 

 

bullet

Select the YZ plane and say OK. The figure will look like the following:

 

 

bullet

Now repeat the same process and reflect the whole figure about the XZ plane. The figure will look like this now.

 

 

bullet

Now we will add the areas up. Go to Preprocessor>Modeling>Operate>Booleans>Add>Areas.

bullet

In the window that comes up click "Pick all".            

   

MATERIAL PROPERTIES:

 

bullet

Go to the ANSYS Main Menu

bullet

Click Preprocessor>Material Props>Material Models.  In the window that comes up, select Structural>Linear>Elastic>Isotropic.

 

 

bullet

Enter 1 for the Material Property Number and click OK. The following window comes up.

 

 

bullet

Fill in 2.1e11 for the Young's modulus and 0.3 for minor Poisson's Ratio.  From the Material Model window, select Structural>Density and enter 7.21e3 for the density. Click OK.

 

 

bullet

Now the material 1 has the properties defined in the above table. We will use this material for the structure.

 

ELEMENT PROPERTIES:

 

bullet

SELECTING ELEMENT TYPE:

bullet

Click Preprocessor>Element Type>Add/Edit/Delete... In the 'Element Types' window that opens click on Add... The following window opens:

 

 

bullet

Type 1 in the Element type reference number.

bullet

Click on Structural Shell and select Elastic 4node 63. Click OK. Close the 'Element types' window.

bullet

Now we need to define the thickness for this element.

bullet

Go to Preprocessor>Real Constants>Add/Edit/Delete...

bullet

In the "Real Constants" dialog box that comes up click on Add

bullet

In the "Element Type for Real Constants" that comes up click OK. The following window comes up.

 

 

bullet

Fill in the relevant values and click on OK.

bullet

We have now defined the thickness of the element.

 

MESHING:

 

bullet

Go to Preprocessor>Meshing>Size Controls>Manual Size>Lines>Picked Lines. Pick all the lines on the outer boundary of the figure and click OK.

bullet

In the menu that comes up type 0.025 in the field for 'Element edge length'.

 

 

bullet

Click on OK.

bullet

Now go to Preprocessor>Meshing>Mesh>Areas>Free.

bullet

Click "pick all" in the "Mesh Areas" dialog box. The meshed model looks like this.

 

 

bullet

Now the blisk is divided into Shell elements.

 

BOUNDARY CONDITIONS AND CONSTRAINTS:

 

bullet

APPLYING BOUNDARY CONDITIONS

bullet

The blisk is fixed around the inner diameter of the disk.

bullet

Go to Main Menu Preprocessor>Loads>Define Loads>Apply>Structural>Displacement>On Lines.

bullet

Select the lines on the inner circumference of the disk and click OK. The following window comes up:

 

 

bullet

Select All DOF and click OK. 

bullet

The model now looks like this:

 

 

bullet

Now the Modeling of the problem is done.

 

SOLUTION:

 

bullet

Go to ANSYS Main Menu>Solution>Analysis Type>New Analysis. The following window comes up:

 

 

bullet

Select Modal and click on OK.

bullet

Now go to Main Menu>Solution>Analysis Type>Analysis Options. The following window comes up:

 

 

bullet

Enter the values shown in the window above and click OK. The following window comes up.

 

 

bullet

Enter 100000000 for the End Frequency.  Then Click OK.

bullet

Go to Solution>Solve>Current LS.

bullet

Wait for ANSYS to solve the problem.

bullet

Click on OK and close the 'Information' window.

 

POST-PROCESSING:

 

bullet

To list the first three frequencies, go to Main Menu>General Postprocessing>Results Summary. The following window will be displayed:

 

 

bullet

To animate the mode shapes, go to Main Menu>General Postprocessing>Read Results>First Set.

bullet

Go to Utility Menu>Plot Controls>Animate>Mode Shape. The following window will come up:

 

 

bullet

Select the required animation: in this case Deformed Shape and click OK.

bullet

The animation will be similar to the ones below. (Don't capture images from these files, they are not the solutions. Just similar to solutions.)
bullet

ANIMATION 1

bullet

ANIMATION 2

bullet

ANIMATION 3

 

MODIFICATIONS:

 

bullet

To plot the deformed shape, go to Main Menu>General Postprocessing>Read Results>First Set.

bullet

Now in the same window go to Plot Results>Contour Plot>Nodal Solution. The following window comes up:

 

 

bullet

Select DOF solution, and select USUM.

bullet

Check the Def + undeformed button.

bullet

Click on OK. The contour plot will look similar to the figure below.

 

 

 

Home Course Info Problems Test Problems Students Reference
Send mail to the Teaching Staff with questions or comments about this web site.