Note

Go to the end to download the full example code.

Thermal-Stress Analysis of a Cooled Turbine Blade#

This example demonstrate 3D thermal-stress analysis of a cooled turbine Blade. Convection, mass flow rate, temperatures are applied to Thermal fluid bodies. Adiabatic surfaces are fixed. Thermal stress analysis is solved and post-process temperature distribution and equivalent von-mises stress for solid region.

Import necessary libraries#

import os
import os.path

from ansys.mechanical.core import App
from ansys.mechanical.core.embedding import AddinConfiguration
from ansys.mechanical.core.examples import delete_downloads, download_file
from matplotlib import image as mpimg
from matplotlib import pyplot as plt

Embed mechanical and set global variables

config = AddinConfiguration("Mechanical")
# config = AddinConfiguration("WorkBench")
config.no_act_addins = True
app = App(config=config)
app.update_globals(globals())
print(app)

cwd = os.path.join(os.getcwd(), "out")


def display_image(image_name):
    plt.figure(figsize=(16, 9))
    plt.imshow(mpimg.imread(os.path.join(cwd, image_name)))
    plt.xticks([])
    plt.yticks([])
    plt.axis("off")
    plt.show()

Ansys Mechanical [Ansys Mechanical Enterprise]
Product Version:252
Software build date: 06/13/2025 11:25:56

Configure graphics for image export#

Graphics.Camera.SetSpecificViewOrientation(ViewOrientationType.Iso)
Graphics.Camera.SetFit()
image_export_format = GraphicsImageExportFormat.PNG
settings_720p = Ansys.Mechanical.Graphics.GraphicsImageExportSettings()
settings_720p.Resolution = GraphicsResolutionType.EnhancedResolution
settings_720p.Background = GraphicsBackgroundType.White
settings_720p.Width = 1280
settings_720p.Height = 720
settings_720p.CurrentGraphicsDisplay = False
Graphics.Camera.Rotate(180, CameraAxisType.ScreenY)

Download and open MECHDAT file#

Open the mechdat file

mechdat_filepath = download_file(
    "Cooled_Turbine_Blade_Meshed.mechdat", "pymechanical", "embedding"
)
app.open(mechdat_filepath)

app.plot()

Static Scene

Interactive Scene

Define Unit System#

Specify unit system

ExtAPI.Application.ActiveUnitSystem = MechanicalUnitSystem.StandardMKS
ExtAPI.Application.ActiveMetricTemperatureUnit = MetricTemperatureUnitType.Kelvin

Define python variables#

Store python variables

GEOM = Model.Geometry
CS_GRP = Model.CoordinateSystems
MSH = Model.Mesh

Define Named Selections#

Specify variables for named selection objects

NS_Passage1 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Passage 1"
][0]
NS_Passage2 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Passage 2"
][0]
NS_Passage3 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Passage 3"
][0]
NS_Passage4 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Passage 4"
][0]
NS_Passage5 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Passage 5"
][0]
NS_Passage6 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Passage 6"
][0]
NS_Passage7 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Passage 7"
][0]
NS_Passage8 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Passage 8"
][0]
NS_Passage9 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Passage 9"
][0]
NS_Passage10 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Passage 10"
][0]
NS_Hole1 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Hole 1"
][0]
NS_Hole2 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Hole 2"
][0]
NS_Hole3 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Hole 3"
][0]
NS_Hole4 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Hole 4"
][0]
NS_Hole5 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Hole 5"
][0]
NS_Hole6 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Hole 6"
][0]
NS_Hole7 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Hole 7"
][0]
NS_Hole8 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Hole 8"
][0]
NS_Hole9 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Hole 9"
][0]
NS_Hole10 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Hole 10"
][0]
NS_Inlet1 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Inlet 1"
][0]
NS_Inlet2 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Inlet 2"
][0]
NS_Inlet3 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Inlet 3"
][0]
NS_Inlet4 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Inlet 4"
][0]
NS_Inlet5 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Inlet 5"
][0]
NS_Inlet6 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Inlet 6"
][0]
NS_Inlet7 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Inlet 7"
][0]
NS_Inlet8 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Inlet 8"
][0]
NS_Inlet9 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Inlet 9"
][0]
NS_Inlet10 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Inlet 10"
][0]
NS_Path1 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Path1"
][0]
NS_Path2 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Path2"
][0]
NS_Faces4 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Faces4"
][0]
NS_Face1 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Face1"
][0]
NS_Face2 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Face2"
][0]
NS_Body1 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Body1"
][0]
NS_Bodies10 = [
    x
    for x in Model.GetChildren[Ansys.ACT.Automation.Mechanical.NamedSelection](True)
    if x.Name == "Bodies10"
][0]

Material Assignment#

Assign materials to blade and fluid bodies and model type for line bodies

GEOM.Children[0].Material = "Blade"
for i in range(1, 11):
    GEOM.Children[i].Material = "Fluid"
    GEOM.Children[i].Children[0].ModelType = PrototypeModelType.ModelPhysicsTypeFluid
    GEOM.Children[i].Children[
        0
    ].FluidDiscretization = FluidDiscretizationType.FluidUpwindExponential

Define coordinate systems and paths#

LCS01 = CS_GRP.AddCoordinateSystem()
LCS01.OriginLocation = NS_Path1

LCS02 = CS_GRP.AddCoordinateSystem()
LCS02.OriginLocation = NS_Path2

CONST_GEOM = Model.AddConstructionGeometry()
PATH01 = CONST_GEOM.AddPath()
PATH01.PathType = PathScopingType.Edge
PATH01.Location = NS_Passage1

PATH02 = CONST_GEOM.AddPath()
PATH02.PathType = PathScopingType.Points
PATH02.StartCoordinateSystem = LCS01
PATH02.EndCoordinateSystem = LCS02

Define Analysis#

Add Steady-state thermal analysis

STAT_THERM = Model.AddSteadyStateThermalAnalysis()
STAT_THERM_SOLN = STAT_THERM.Solution

Define Loads and Boundary Conditions#

Setup loads and supports in steady-state thermal analysis

CONV01 = STAT_THERM.AddConvection()
CONV01.Location = NS_Hole1
CONV01.FilmCoefficient.Output.DiscreteValues = [Quantity("295430 [W m^-1 m^-1 K^-1]")]
CONV01.HasFluidFlow = True
CONV01.FluidFlowSelection = NS_Passage1

CONV02 = STAT_THERM.AddConvection()
CONV02.Location = NS_Hole2
CONV02.FilmCoefficient.Output.DiscreteValues = [Quantity("296290 [W m^-1 m^-1 K^-1]")]
CONV02.HasFluidFlow = True
CONV02.FluidFlowSelection = NS_Passage2

CONV03 = STAT_THERM.AddConvection()
CONV03.Location = NS_Hole3
CONV03.FilmCoefficient.Output.DiscreteValues = [Quantity("300760 [W m^-1 m^-1 K^-1]")]
CONV03.HasFluidFlow = True
CONV03.FluidFlowSelection = NS_Passage3

CONV04 = STAT_THERM.AddConvection()
CONV04.Location = NS_Hole4
CONV04.FilmCoefficient.Output.DiscreteValues = [Quantity("314160 [W m^-1 m^-1 K^-1]")]
CONV04.HasFluidFlow = True
CONV04.FluidFlowSelection = NS_Passage4

CONV05 = STAT_THERM.AddConvection()
CONV05.Location = NS_Hole5
CONV05.FilmCoefficient.Output.DiscreteValues = [Quantity("314950 [W m^-1 m^-1 K^-1]")]
CONV05.HasFluidFlow = True
CONV05.FluidFlowSelection = NS_Passage5

CONV06 = STAT_THERM.AddConvection()
CONV06.Location = NS_Hole6
CONV06.FilmCoefficient.Output.DiscreteValues = [Quantity("301990 [W m^-1 m^-1 K^-1]")]
CONV06.HasFluidFlow = True
CONV06.FluidFlowSelection = NS_Passage6

CONV07 = STAT_THERM.AddConvection()
CONV07.Location = NS_Hole7
CONV07.FilmCoefficient.Output.DiscreteValues = [Quantity("302470 [W m^-1 m^-1 K^-1]")]
CONV07.HasFluidFlow = True
CONV07.FluidFlowSelection = NS_Passage7

CONV08 = STAT_THERM.AddConvection()
CONV08.Location = NS_Hole8
CONV08.FilmCoefficient.Output.DiscreteValues = [Quantity("443430 [W m^-1 m^-1 K^-1]")]
CONV08.HasFluidFlow = True
CONV08.FluidFlowSelection = NS_Passage8

CONV09 = STAT_THERM.AddConvection()
CONV09.Location = NS_Hole9
CONV09.FilmCoefficient.Output.DiscreteValues = [Quantity("285270 [W m^-1 m^-1 K^-1]")]
CONV09.HasFluidFlow = True
CONV09.FluidFlowSelection = NS_Passage9

CONV010 = STAT_THERM.AddConvection()
CONV010.Location = NS_Hole10
CONV010.FilmCoefficient.Output.DiscreteValues = [Quantity("895860 [W m^-1 m^-1 K^-1]")]
CONV010.HasFluidFlow = True
CONV010.FluidFlowSelection = NS_Passage10

MFLOW_RT01 = STAT_THERM.AddMassFlowRate()
MFLOW_RT01.Location = NS_Passage1
MFLOW_RT01.Magnitude.Output.DiscreteValues = [
    Quantity("0[kg sec^-1]"),
    Quantity("-0.0228[kg sec^-1]"),
]

MFLOW_RT02 = STAT_THERM.AddMassFlowRate()
MFLOW_RT02.Location = NS_Passage2
MFLOW_RT02.Magnitude.Output.DiscreteValues = [
    Quantity("0[kg sec^-1]"),
    Quantity("-0.0239[kg sec^-1]"),
]

MFLOW_RT03 = STAT_THERM.AddMassFlowRate()
MFLOW_RT03.Location = NS_Passage3
MFLOW_RT03.Magnitude.Output.DiscreteValues = [
    Quantity("0[kg sec^-1]"),
    Quantity("-0.0228[kg sec^-1]"),
]

MFLOW_RT04 = STAT_THERM.AddMassFlowRate()
MFLOW_RT04.Location = NS_Passage4
MFLOW_RT04.Magnitude.Output.DiscreteValues = [
    Quantity("0[kg sec^-1]"),
    Quantity("-0.0243[kg sec^-1]"),
]

MFLOW_RT05 = STAT_THERM.AddMassFlowRate()
MFLOW_RT05.Location = NS_Passage5
MFLOW_RT05.Magnitude.Output.DiscreteValues = [
    Quantity("0[kg sec^-1]"),
    Quantity("-0.0239[kg sec^-1]"),
]

MFLOW_RT06 = STAT_THERM.AddMassFlowRate()
MFLOW_RT06.Location = NS_Passage6
MFLOW_RT06.Magnitude.Output.DiscreteValues = [
    Quantity("0[kg sec^-1]"),
    Quantity("-0.0242[kg sec^-1]"),
]

MFLOW_RT07 = STAT_THERM.AddMassFlowRate()
MFLOW_RT07.Location = NS_Passage7
MFLOW_RT07.Magnitude.Output.DiscreteValues = [
    Quantity("0[kg sec^-1]"),
    Quantity("-0.0232[kg sec^-1]"),
]

MFLOW_RT08 = STAT_THERM.AddMassFlowRate()
MFLOW_RT08.Location = NS_Passage8
MFLOW_RT08.Magnitude.Output.DiscreteValues = [
    Quantity("0[kg sec^-1]"),
    Quantity("-0.00799[kg sec^-1]"),
]

MFLOW_RT09 = STAT_THERM.AddMassFlowRate()
MFLOW_RT09.Location = NS_Passage9
MFLOW_RT09.Magnitude.Output.DiscreteValues = [
    Quantity("0[kg sec^-1]"),
    Quantity("-0.00499[kg sec^-1]"),
]

MFLOW_RT10 = STAT_THERM.AddMassFlowRate()
MFLOW_RT10.Location = NS_Passage10
MFLOW_RT10.Magnitude.Output.DiscreteValues = [
    Quantity("0[kg sec^-1]"),
    Quantity("-0.00253[kg sec^-1]"),
]

TEMP01 = STAT_THERM.AddTemperature()
TEMP01.Location = NS_Inlet1
TEMP01.Magnitude.Output.DiscreteValues = [Quantity("0[K]"), Quantity("348.83[K]")]

TEMP02 = STAT_THERM.AddTemperature()
TEMP02.Location = NS_Inlet2
TEMP02.Magnitude.Output.DiscreteValues = [Quantity("0[K]"), Quantity("349.32[K]")]

TEMP03 = STAT_THERM.AddTemperature()
TEMP03.Location = NS_Inlet3
TEMP03.Magnitude.Output.DiscreteValues = [Quantity("0[K]"), Quantity("339.49[K]")]

TEMP04 = STAT_THERM.AddTemperature()
TEMP04.Location = NS_Inlet4
TEMP04.Magnitude.Output.DiscreteValues = [Quantity("0[K]"), Quantity("342.3[K]")]

TEMP05 = STAT_THERM.AddTemperature()
TEMP05.Location = NS_Inlet5
TEMP05.Magnitude.Output.DiscreteValues = [Quantity("0[K]"), Quantity("333.99[K]")]

TEMP06 = STAT_THERM.AddTemperature()
TEMP06.Location = NS_Inlet6
TEMP06.Magnitude.Output.DiscreteValues = [Quantity("0[K]"), Quantity("364.95[K]")]

TEMP07 = STAT_THERM.AddTemperature()
TEMP07.Location = NS_Inlet7
TEMP07.Magnitude.Output.DiscreteValues = [Quantity("0[K]"), Quantity("343.37[K]")]

TEMP08 = STAT_THERM.AddTemperature()
TEMP08.Location = NS_Inlet8
TEMP08.Magnitude.Output.DiscreteValues = [Quantity("0[K]"), Quantity("365.41[K]")]

TEMP09 = STAT_THERM.AddTemperature()
TEMP09.Location = NS_Inlet9
TEMP09.Magnitude.Output.DiscreteValues = [Quantity("0[K]"), Quantity("408.78[K]")]

TEMP10 = STAT_THERM.AddTemperature()
TEMP10.Location = NS_Inlet10
TEMP10.Magnitude.Output.DiscreteValues = [Quantity("0[K]"), Quantity("453.18[K]")]

TEMP11 = STAT_THERM.AddTemperature()
TEMP11.Location = NS_Faces4
TEMP11.Magnitude.Output.DiscreteValues = [Quantity("0[K]"), Quantity("568[K]")]

Insert results#

TEMP_RST01 = STAT_THERM_SOLN.AddTemperature()
TEMP_RST01.Location = NS_Body1
TEMP_RST02 = STAT_THERM_SOLN.AddTemperature()
TEMP_RST02.Location = NS_Bodies10
TEMP_RST03 = STAT_THERM_SOLN.AddTemperature()
TEMP_RST03.ScopingMethod = GeometryDefineByType.ResultFileItem
TEMP_RST03.ItemType = ResultFileItemType.ElementNameIDs
TEMP_RST03.SolverComponentIDs = "SURF152"
TEMP_RST04 = STAT_THERM_SOLN.AddTemperature()
TEMP_RST04.Location = PATH01

Solve#

STAT_THERM_SOLN.Solve(True)
STAT_THERM_SS = STAT_THERM_SOLN.Status

Messages#

Messages = ExtAPI.Application.Messages
if Messages:
    for message in Messages:
        print(f"[{message.Severity}] {message.DisplayString}")
else:
    print("No [Info]/[Warning]/[Error] Messages")

[Info] The requested license was received from the License Manager after 22 seconds.

Results#

Temperature Plots

Tree.Activate([TEMP_RST01])
Graphics.Camera.SetFit()
Graphics.ExportImage(
    os.path.join(cwd, "temperature_turbine.png"), image_export_format, settings_720p
)
display_image("temperature_turbine.png")

Define Analysis#

Add static structural analysis

STAT_STRUC = Model.AddStaticStructuralAnalysis()
STAT_STRUC.ImportLoad(Model.Analyses[0])
STAT_STRUC_SOLN = STAT_STRUC.Solution

Define Loads and Boundary Conditions#

Setup loads and supports in static analysis

FIX_SUP01 = STAT_STRUC.AddFixedSupport()
FIX_SUP01.Location = NS_Face1

FIX_SUP02 = STAT_STRUC.AddFixedSupport()
FIX_SUP02.Location = NS_Face2

Insert results#

EQV_STRS01 = STAT_STRUC_SOLN.AddEquivalentStress()
EQV_STRS01.Location = NS_Body1

Solve#

STAT_STRUC_SOLN.Solve(True)
STAT_STRUC_SS = STAT_STRUC_SOLN.Status

Messages#

Messages = ExtAPI.Application.Messages
if Messages:
    for message in Messages:
        print(f"[{message.Severity}] {message.DisplayString}")
else:
    print("No [Info]/[Warning]/[Error] Messages")

[Warning] Thermal Fluid Line bodies are not included in non-thermal solution.
[Warning] Not all Named Selections were successfully written because one or more names are not valid for the solver
[Info] The Imported Body Temperature load will set the Program Controlled Tabular Loading option to Off when Mapped Data property is set to Input file, and the specified analysis time values match the step end times.
[Info] The requested license was received from the License Manager after 22 seconds.

Results#

Equivalent Stress Plots

Tree.Activate([EQV_STRS01])
Graphics.Camera.SetFit()
Graphics.ExportImage(
    os.path.join(cwd, "equivalent_stess_turbine.png"),
    image_export_format,
    settings_720p,
)
display_image("equivalent_stess_turbine.png")

Cleanup#

Save project

app.save(os.path.join(cwd, "cooled_turbine.mechdat"))
app.new()

Delete example file

delete_downloads()

True

Total running time of the script: (0 minutes 27.782 seconds)

Gallery generated by Sphinx-Gallery