This article provides a step-by-step tutorial for a thermomechanical analysis of an engine piston under maximum pressure and temperatures.
This tutorial teaches how to:
We are following the typical SimScale workflow:
Firstly, you can click the button below and it will copy the tutorial project containing the geometry into your Workbench.
The following picture shows what should be visible after importing the tutorial project.
The first step of a simulation setup is to create a new simulation.
Hitting the ‘Create Simulation’ button leads to the following options:
Choose ‘Thermomechanical’ as the analysis type and click the ‘Create Simulation’ button to confirm. After that, a simulation tree showing you all the necessary steps to set up your simulation will appear in your workbench.
Before setting up the simulation, we recommend creating topological entity sets. You can do this at any point of the simulation setup, however, it is easier to do it in the beginning.
The sets for the top, the rings, and the ring grooves have already been created for you according to the following picture:
Now, you will only need to create the set for the interior and the skirt of the piston, which is everything other than what is already assigned to sets. Here are the steps:
Did you know?
You can find out more quick selection tips in the following article: Viewer Tips & Tricks
In this section, we will define the physics of our model.
Moving from the top to the bottom of the simulation tree, we have the Model tab, where we define the magnitude and direction of the gravity:
The magnitude of the gravity that we use is ‘9.81’ \(m^2/s\) in the direction ‘-1’ in \(e_y\).
Afterward, we will define the material of our piston, which will be Aluminium. After clicking on the ‘+’ button next to Materials, the following list appears:
Choose ‘Aluminium’ from the list and confirm by pressing ‘Apply’. Since the CAD model from this tutorial contains a single volume, it will automatically receive the aluminium definition.
A thermomechanical analysis will need two sets of boundary conditions: thermal and mechanical. In this section, we will define both sets of boundary conditions.
a. Thermal Boundary Conditions
We will use several Convective heat flux definitions as our thermal boundary conditions. You can therefore add a boundary condition by hitting the ‘+’ and selecting a ‘Convective heat flux’ boundary condition:
A convective heat flux boundary condition generates a heat flux based on the surface area, heat transfer coefficient, and temperature difference between the Reference temperature and the surface temperature.
The first convective heat flux boundary condition will have the following settings:
The \(T_0\) Reference temperature is 741 \(°C\), with a Heat transfer coefficient of 450 \(W/(K.m^2)\). This boundary condition is assigned to the ‘top’ topological entity set. Alternatively, you can click directly on the top face in the viewer.
Follow the same procedure for the rows within the table below:
|Entity Set||Reference temperature [°C]||Heat transfer coefficient [W/(K*m^2)]|
|Groove – Ring 1||180||1000|
|Groove – Ring 2||160||400|
|Groove – Ring 3||140||400|
|Interior and skirt||120||650|
b. Mechanical Boundary Conditions
Now it is time for the mechanical conditions, which will be three Pressure definitions (top and the first two rings) and two Remote displacements. At first, we will create a new boundary condition, this time choosing ‘Pressure’ from the list, and give it the following definition:
The first pressure boundary condition is defined to the ‘top’ face of the piston, with a value of ‘2e7’ \(Pa\). The pressure boundary conditions for the other parts can be seen in the table below:
Now, you can create a ‘Remote displacement’ boundary condition and give it the following definition:
The first remote displacement is assigned to the top pin of the piston with the following coordinates for the External point:
Following the same workflow, create a second ‘Remote displacement’ boundary condition, this time for the other side of the piston:
For the second remote displacement at the opposite side of the piston, please use the following coordinates for the External point:
Now all boundary conditions are assigned and we can proceed to create the mesh.
We will use the standard algorithm for our mesh, which is a good choice in general as it is quite automated and delivers good results for most geometries.
You will only need to change the Sizing to Manual and set our maximum edge length as 1.8e-3 \(m\). Make sure your settings look like the picture below:
When the mesh has been generated, it will look like this:
With default settings, a second-order mesh is created. Such meshes have additional nodes in-between two connecting nodes, which helps to calculate deformations more accurately.
At this point, you can run your thermomechanical analysis for the engine piston by clicking on the ‘+’ button next to Simulation runs.
Moreover, you can change the name of your simulation to your liking and start the simulation by clicking ‘Start’
After the simulation has finished, you can access the simulation results by either clicking on ‘Solution Fields’ or by pressing ‘Post-process results’:
You can choose what results to visualize by going to the Filters panel:
For example, we will start by plotting the Temperature field:
To achieve this visualization, we performed a couple of changes from the default state:
We can also visualize the von Mises stress distribution with the same method:
By lowering the maximum value of the color scale, we can create a threshold and find all the regions in the part where the stress exceeds this value. Those regions are colored in red in Figure 23 and are found in the support holes and the top face, where applied and reaction forces appear.
Finally, we want to visualize the deformation of the piston head. For this, make sure that the ‘Displacement’ option from the toolbar is active:
You might notice that the deformation is not apparent. This happens because the magnitude of the deformations is relatively small. To fix this, we changed the Scaling factor to 20:
Also, to be able to find the magnitude of the deformations, you can switch the Coloring to ‘Displacement > Magnitude’:
Finally, we come up with the following plot:
We can see how the piston head deforms due to the applied loads, boundary, and thermal conditions. The top face sinks, while the ring portion bulks out. We can also see the expansion at the bottom portion, which can be important for design purposes.
Congratulations! You finished the thermomechanical analysis of an engine piston tutorial!
Last updated: August 22nd, 2022
We appreciate and value your feedback.