# Workshop 2 Session 3: Aircraft jet engine vibration analysis

**URL:** <https://www.simscale.com/forum/t/workshop-2-session-3-aircraft-jet-engine-vibration-analysis/87889>\
**Category:** Aerospace Workshop\
**Created:** [March 5, 2020, 12:33am UTC](https://www.simscale.com/forum/t/workshop-2-session-3-aircraft-jet-engine-vibration-analysis/87889 "2020-03-05T00:33:33Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![Ricardopg](https://www.simscale.com/forum/user_avatar/www.simscale.com/ricardopg/32/26909_2.png) [@Ricardopg](https://www.simscale.com/forum/u/Ricardopg)\
**Post date:** [March 5, 2020, 12:33am UTC](https://www.simscale.com/forum/t/workshop-2-session-3-aircraft-jet-engine-vibration-analysis/87889/1 "2020-03-05T00:33:33Z")

</div>

## Recording

[Aerospace Workshop II feat. EUROAVIA (Session 3) ― Vibration Analysis of a Jet Engine](https://www.youtube.com/watch?v=39YdR_z2Kvc)

## Exercise

Have you ever thought that what exactly will happen to the jet engine if one of the fan blade breaks out? It will rather produce vibration due to overall change in center of gravity of a jet engine. But how much it will vibrate? In this exercise we will be looking in to the vibration of the jet engine due to blade loss.

The cutout of the jet engine geometry used in this exercise is shown in the figure below. Only the outer cowling (with strator and struts) and pylon are taken for this analysis, whereas the internal whole jet engine and wing will be ignored since we are only interested in the impact of its vibration on to the cowling via strator and struts.

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/f/0/f02d77a5867780d2c0440ef6c071289af30960cd.jpeg)

- _The base model for whole assembly except the internals of jet engine is provided by [Aisak](https://grabcad.com/aisak-2) from [GrabCAD](https://grabcad.com/library/airbus-a360-800-prototype)._
- _The base model for the internals of jet engine is provided by [Goutam Das](https://grabcad.com/goutam.das-1) from [GrabCAD](https://grabcad.com/library/turbo_fan_jet_engine-1)._

This exercise focuses on the vibration study of the jet engine cowling due to missing single or multiple fan blades. First, a frequency analysis will be performed in order to get the natural frequencies of the the model. Next, harmonic analysis will be performed on a specific frequency range of interest with the missing blade vibration load in order to see at which frequency the vibration will be highest. Harmonic analysis produces a sinusoidal load through which one can detect the peak displacement value at a specific frequency. The model with missing single and multiple blades are shown in the figure below.

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/4/d/4d3cc42ac39e385d6463a1a9a21f8ba6db647ecd.png)  
_Jet engine assembly with 1 missing blade._

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/b/6/b6947ac2f0d3e5076f742b42474cdf0662099ca3.png)  
_Jet engine assembly with 2 missing blades._

The different combinations for performing the simulations are:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/b/7/b77e882848e41eff23517b0a0cb8f239d6fe2171.png)

## _Please note that separate simulation needs to be set up for each configuration_

## Step-by-Step

**Import the project by clicking the link below.** (Crtl + Click to open in a new tab)

[Click to Import the Homework Project](https://www.simscale.com/workbench?publiclink=ce93bfd3-86f2-4776-b2c0-5e74eedcab78)

We will start off by clicking on the geometry **jet engine model** and then on **Create Simulation**.

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/1/8/189f3405ce09347f6e628e5a40ce177be511a9d7.jpeg)

For the first simulation, we will choose a **Frequency Analysis**.

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/6/0/6008eade3ceb32939cc3d4335286a325a2b87e9c.png)

### Meshing

In the simulation tree, navigate to **Mesh**. The **Algorithm** will be **Standard**. Furthermore, change the **Fineness** to **5**. Afterwards, click on **Generate** to start the meshing operation. The number of processors is automatically chosen by SimScale to optimize the usage of core hours.

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/f/e/fe104a6cb0735351449d471f642ed8837c57ccfa.png)

Now it’s time to create topological entity sets. These will help us to assign boundary conditions later on.  
To create a new set, follow these steps:

- Make sure the **Geometry** is selected in the simulation tree;
- Select a surface or surfaces that will constitute the new set;
- In the right-hand side panel, click on the **+** icon next to **Topological Entity Sets** ;
- Name your set as you seem fit.

The first topological entity set will be named **fixation**. It consists of the fixed supports:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/3/c/3ce62868d3622d55baeadf50da29d22e367bb668.png)

And then name this set as **fixation** :

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/d/c/dc40b736c977ecbc247cc3febd93bf185b02d6d8.png)

The second and last topological entity set will be named **engine-vibration**. The 4 inner round surfaces of the cowling will be selected:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/0/3/03aecb473582c84716a940fff2569a4cb2306167.jpeg)

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/1/e/1e5e0cba283ab11e1a3468f6e8ecce33181a72e2.jpeg)

In the simulation tree, navigate to **Materials**. From the list, select **Aluminium** and hit **Apply**.

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/e/5/e5e8b9dd047fae692ce8c42100368cdb202d052c.png)

Assign it to the two volumes ( **1** and **2** ), then click on the blue mark to save.

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/c/8/c84787051616cfcea6ce8e8a4dabe8fd7932a74a.png)

### Boundary conditions

For the frequency analysis, the only boundary condition needed is the fixation of the wing. Click on the **+** sign next to **Boundary conditions**. From the drop down menu, choose **Fixed value**. The displacement in all directions will be 0. Assign it to the fixation topological entity set. You can rename this boundary condition **fixation** :

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/9/5/955420302d5ca0e72c91be5a9e48f6a1d12b5d80.png)

In **Simulation control** , make sure that **Eigenfrequency Scope** is set to **First modes** and that **20** modes will be analyzed:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/b/b/bb4ea491699a838cbb350964fbea192e144f6108.png)

Before starting the simulation, expand **Contacts** in the simulation tree. Click on **Bonded 1** and lower the **Tolerance** to **1e-4**.

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/6/f/6f5cd2bcd7f1d8d43ebd97fd65d7962c2f109ed2.png)

Now we can start our simulation. Go to **Simulation Runs** and click **Start**.

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/e/c/ecf6913cf09c141cdda1935fd0dedfefad5e43ce.png)

After the simulation finishes running, navigate to **Plots** to inspect the Eigenmode _versus_ Eigenfrequency (Hz) plot:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/9/e/9eaf85b5eac538696e7f5bfaeb7e8b6877dc23a1.png)

Jet engines usually run between 3000 and 6000 RPM (50 to 100 Hz). A total of 5 modes have eigenfrequencies under 100 Hz. 5 modes have their frequencies between 100 and 150 Hz.

### Create further configurations

**Configuration 2**

Based on the frequency analysis result, we’re going to perform a series of harmonic analysis. Configuration 2 will simulate one missing blade. To start, create a new simulation, this time a **Harmonic** analysis:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/4/2/42a309a44cfd8b3f778f626a5e6cfdf71f832f8e.png)

To make the simulations easily distinguishable, rename this simulation **Vibration analysis - 1 missing blade**.

In **Contacts** , set **Tolerance** to **1e-4** :

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/c/c/cc176cdb72f3ddd4bca54c783b4a77bfb128bbd5.png)

**Materials** will be **Aluminium** , just like the previous simulation. Assign it to volumes **1** and **2** :

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/a/e/ae0ac608c2d0bdcf39a58a338d68d4466b3ff775.png)

Now, in **Boundary conditions** , the first one will be a **Fixed value** , with **0**  **displacement** in all directions. It will be assigned to the **fixation** topological entity set:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/e/4/e43ff73a8d8edcf5711e4dab23611766e7170cf2.png)

The applied force load due to the missing blades is estimated considering the following formulaton:

F=m\*Δx\*4π^2\*f^2

Where:

- m is the mass of a missing blade. Considering a blade made of titanium (density equivalent to _4500_ kg/m³), the mass of a blade will be roughly 7.65 kg.
- Δx is the distance of the dislocated center of mass, which is found out to be _5e-3_ m;

To apply this condition, we will create a **Force** boundary condition. Rename the newly created BC to **Jet engine load - Fy**. Then follow the steps below:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/7/d/7d7f85eed372abe2446ce78cae225ed27767b6a8.png)

And input the following formula: 7.65\*5\*e-3\*4\*pi^2\*(f^2)

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/d/c/dcd424967ad44bc25c9301c4a7629d8e26fb3384.png)

Lastly, assign this boundary condition to the **engine-vibration** topological entity set:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/b/d/bde9b38cf72e19fc44d15449649341835c6fee10.png)

We will now make a similar boundary condition, but this time for the z direction. To do so, right click on the **Jet engine load - Fy** boundary condition and **Duplicate** it:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/3/c/3c2720102f5c800a7fe70c53b7b620ec385f6564.png)

Rename this boundary condition **Jet engine load - Fz**. Change **fy** to **0** and **fz** to **z**. Also give a value of **90** to **Phase angle** :

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/7/9/7987795472a595de163bd9e20d28a0a13611de92.png)

In **Simulation control** , change **Excitation frequencies** to **Frequency list**. Keeping in mind that engines usually run between 50 and 100 Hz, we will go from **0** to **120** Hz with steps of **10** Hz:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/2/9/290feac98216099787d9dfb7890eb9e181c7ef8d.png)

Expand **Result control** in the simulation tree. To save some computing time, we will delete some controls that are not of our interest. Please delete cauchy and von Mises stresses, aswell as total strain.

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/7/8/784fd5c3a602c8678731c6f04e6e2f6a64f6d315.png)

Next, we will create extra result controls under **Point data**. Rename the first one to **Dx**. It will be a **Displacement** control in the **X** direction. Change **Complex number** to **Magnitude and phase** :

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/9/b/9b4eb00ca208adff7ab72d107bf348956e5162f5.png)

After clicking on the **+** button next to **Geometry primitives** , we will create one point where the displacement will be assessed. Give the following coordinates (x, y, z) for your point: -2.333779, -0.115045, -0.407097.

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/a/7/a72ceee127faebf9d78749821a6df233ec07dca7.png)

Similarly, create two extra controls for the **Y** and **Z** directions. The same point can be used.

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/a/9/a90271160e804da95f26f6939e281ad8c74724b5.png)

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/2/7/279da8444e44de464952b605861c813886acc8b5.png)

That’s it for configuration 2. You can proceed to running the simulation.

**Configuration 3**

The first step is to duplicate the vibration analysis that was just set up. Rename is to \*\* _Vibration analysis - 2 missing blades_\*\*

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/4/c/4c360116f8ca5268522b83f0170c5d86e4af1fc5.png)

In **Boundary conditions** , change both load formulas to 15.3\*8.42\*e-3\*4\*pi^2\*f^2.

You can proceed to creating a new **Simulation Run** for configuration 3.

**Configuration 4**

For configuration 4, duplicate the simulation named _Vibration analysis - 1 missing blade_.  
Based on the results from ‘Configuration 2’, it is evident that the peak of displacement is between 40 to 100 Hz. Therefore, we will now redefine our frequency range and rerun the case in order to get the exact frequency where the peak displacement is taking place.

Go to **Simulation Control** and change **Start frequency** to **40** Hz, **End frequency** to **100** Hz and **Frequency stepping** to **5**.

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/5/2/52206981ef79803ecb9f45adefe865feba4ffc06.png)

After changing these values, please create a simulation run for configuration 4.

**Configuration 5**

Duplicate the simulation named _Vibration analysis - 2 missing blades_. Change the start and end frequency values, aswell as frequency stepping to those from configuration 4. Then please run the last simulation.

### Post-processing

The vibration was studied for the missing blade vibration load within the frequency range of 0 to 120 Hz (0 to 7200 rpm). The most important result in this case is the displacement plot in y and z direction. Below are the displacement plots of _Configuration 2_ and _Configuration 3_:

Configuration 2 - Dy:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/9/5/95f45a12a98cfc1eee9011103cf8bc423dc09501.png)

Configuration 2- Dz:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/3/c/3c7d29f043ab4f0e104394ed0bfb7d898ad489d0.png)

Configuration 3 - Dy:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/9/c/9c687e951a8be1f3b589f104ae88086e89f18df5.png)

Configuration 3 - Dz:

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/2/3/23b8521f89e111e997e875e629b904a052a1cbc8.png)

Despite the extra blade that was missing, the displacements observed in configurations 2 and 3 are almost the same.

Now taking a closer look at Configurations 4 and 5, in case engine blades are missing, the maximum vibration in the cowling will occur at the frequencies of 45 and 90 Hz:

Configuration 5 - Dy

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/9/f/9fcbadeb3c128ffc11d30f1f82ae6e968b81fafb.png)

Configuration 5 - Dz

 ![image](https://www.simscale.com/forum/uploads/default/original/3X/6/3/63a70e42aa0fb11f5f997a316beeaf6591aa0b3b.png)

---

<div class="post-metadata">

**Author:** ![jousefm](https://www.simscale.com/forum/user_avatar/www.simscale.com/jousefm/32/34908_2.png) [@jousefm](https://www.simscale.com/forum/u/jousefm)\
**Post date:** [March 20, 2020, 2:00pm UTC](https://www.simscale.com/forum/t/workshop-2-session-3-aircraft-jet-engine-vibration-analysis/87889/2 "2020-03-20T14:00:08Z")

</div>