.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "examples/fluent-mechanical/wf_fm_02_mechanical.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_examples_fluent-mechanical_wf_fm_02_mechanical.py>`
        to download the full example code.

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_examples_fluent-mechanical_wf_fm_02_mechanical.py:


.. _ref_fluent_mechanical_02-mechanical:

Thermo-mechanical assessment of representative exhaust manifold model
#####################################################################

Background
----------

The turbine housing/exhaust manifold is a critical component as it is subjected to extreme
temperature cycling. This induces fatigue in the housing material and leads to early
failure. Simulation helps to design housing and other parts for extended service life.

Thermo-mechanical workflow
--------------------------

A comprehensive CFD analysis of the exhaust manifold component is executed as a
steady-state examination, focusing on the specified duty cycle. This process involves
two recurring duty cycles, during which the temperature and other crucial operating
parameters, such as mass flow and pressure, maintain a consistent pseudo-steady state.

Three unique operating conditions are chosen from the duty cycle, and steady-state
calculations are performed corresponding to the exhaust gas temperature in each case.
The resulting Heat Transfer Coefficients (HTCs) and temperatures at the solid-fluid
interface are exported to a CSV file, which serves as the input for subsequent thermal
calculations in the Mechanical run.

The Mechanical run encompasses Structural Transient Thermal and Non-linear (NL) Static
analyses, employing temperature-dependent Bilinear Kinematic hardening material
properties for a thorough structural evaluation.

In the Transient Thermal run, imported HTCs and reference temperatures from the
CFD runs are mapped onto the structural mesh of the exhaust manifold assembly. External
convection loads and other boundary conditions are applied to calculate the temperature
distribution across the exhaust manifold assembly.

Subsequently, a non-linear static analysis is conducted to determine the plastic strain,
taking into account the temperature distribution calculated in the Transient Thermal
solve and other relevant structural boundary conditions.

This approach leads to a more impactful and precise understanding of the exhaust
manifold's performance subjected to thermal cycling.

.. GENERATED FROM PYTHON SOURCE LINES 67-77

.. code-block:: Python


    import os
    from pathlib import Path

    from ansys.mechanical.core import launch_mechanical
    from ansys.mechanical.core.examples import download_file
    from matplotlib import image as mpimg
    from matplotlib import pyplot as plt









.. GENERATED FROM PYTHON SOURCE LINES 84-89

Parameters for the script
-------------------------
The following parameters are used to control the script execution. You can
modify these parameters to suit your needs.


.. GENERATED FROM PYTHON SOURCE LINES 89-93

.. code-block:: Python

    GRAPHICS_BOOL = False  # Set to True to display the graphics
    OUTPUT_DIR = Path(Path(__file__).parent, "outputs")  # Output directory









.. GENERATED FROM PYTHON SOURCE LINES 98-101

Start a PyMechanical app
------------------------


.. GENERATED FROM PYTHON SOURCE LINES 101-114

.. code-block:: Python

    mechanical = launch_mechanical(batch=True, cleanup_on_exit=False)
    print(mechanical)


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






.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    Ansys Mechanical [Ansys Mechanical Enterprise]
    Product Version:242
    Software build date: 06/03/2024 09:35:09





.. GENERATED FROM PYTHON SOURCE LINES 115-119

Input files needed for the simulation
---------------- --------------------
Download the input files needed for the simulation.


.. GENERATED FROM PYTHON SOURCE LINES 119-158

.. code-block:: Python

    geometry_path = download_file(
        "Exhaust_Manifold_Geometry.pmdb", "pyansys-workflow", "exhaust-manifold", "pymechanical"
    )
    material_path = download_file(
        "Nonlinear_Material.xml", "pyansys-workflow", "exhaust-manifold", "pymechanical"
    )

    # Files necessary for the thermal simulation from fluent analysis

    temp_htc_data_high_path = os.path.join(OUTPUT_DIR, "htc_temp_mapping_HIGH_TEMP.csv")
    temp_htc_data_med_path = os.path.join(OUTPUT_DIR, "htc_temp_mapping_MEDIUM_TEMP.csv")
    temp_htc_data_low_path = os.path.join(OUTPUT_DIR, "htc_temp_mapping_LOW_TEMP.csv")

    all_input_files = {
        "geometry_path": geometry_path,
        "material_path": material_path,
        "temp_htc_data_high_path": temp_htc_data_high_path,
        "temp_htc_data_med_path": temp_htc_data_med_path,
        "temp_htc_data_low_path": temp_htc_data_low_path,
    }

    # Upload to Mechanical Remote session server and get the file paths

    project_directory = mechanical.project_directory
    print(f"project directory = {project_directory}")
    for input_file_name, input_file_path in all_input_files.items():

        # Upload the file to the project directory.
        mechanical.upload(file_name=input_file_path, file_location_destination=project_directory)

        # Build the path relative to project directory.
        base_name = os.path.basename(input_file_path)
        combined_path = os.path.join(project_directory, base_name)
        server_file_path = combined_path.replace("\\", "\\\\")
        mechanical.run_python_script(f"{input_file_name} = '{server_file_path}'")
        result = mechanical.run_python_script(f"{input_file_name}")
        print(f"path of {input_file_name} on server: {result}")






.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    project directory = /tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/

    Uploading Exhaust_Manifold_Geometry.pmdb to dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/.:   0%|          | 0.00/1.51M [00:00<?, ?B/s]
    Uploading Exhaust_Manifold_Geometry.pmdb to dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/.: 100%|██████████| 1.51M/1.51M [00:00<00:00, 366MB/s]
    path of geometry_path on server: /tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/Exhaust_Manifold_Geometry.pmdb

    Uploading Nonlinear_Material.xml to dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/.:   0%|          | 0.00/470k [00:00<?, ?B/s]
    Uploading Nonlinear_Material.xml to dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/.: 100%|██████████| 470k/470k [00:00<00:00, 312MB/s]
    path of material_path on server: /tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/Nonlinear_Material.xml

    Uploading htc_temp_mapping_HIGH_TEMP.csv to dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/.:   0%|          | 0.00/1.61M [00:00<?, ?B/s]
    Uploading htc_temp_mapping_HIGH_TEMP.csv to dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/.: 100%|██████████| 1.61M/1.61M [00:00<00:00, 1.19GB/s]
    path of temp_htc_data_high_path on server: /tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/htc_temp_mapping_HIGH_TEMP.csv

    Uploading htc_temp_mapping_MEDIUM_TEMP.csv to dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/.:   0%|          | 0.00/1.61M [00:00<?, ?B/s]
    Uploading htc_temp_mapping_MEDIUM_TEMP.csv to dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/.: 100%|██████████| 1.61M/1.61M [00:00<00:00, 1.36GB/s]
    path of temp_htc_data_med_path on server: /tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/htc_temp_mapping_MEDIUM_TEMP.csv

    Uploading htc_temp_mapping_LOW_TEMP.csv to dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/.:   0%|          | 0.00/1.61M [00:00<?, ?B/s]
    Uploading htc_temp_mapping_LOW_TEMP.csv to dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/.: 100%|██████████| 1.61M/1.61M [00:00<00:00, 577MB/s]
    path of temp_htc_data_low_path on server: /tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/htc_temp_mapping_LOW_TEMP.csv




.. GENERATED FROM PYTHON SOURCE LINES 159-162

Configure graphics for image export
-----------------------------------


.. GENERATED FROM PYTHON SOURCE LINES 162-181

.. code-block:: Python


    mechanical.run_python_script(
        """
    ExtAPI.Graphics.Camera.SetSpecificViewOrientation(
        Ansys.Mechanical.DataModel.Enums.ViewOrientationType.Iso
    )
    ExtAPI.Graphics.Camera.SetFit()
    image_export_format = Ansys.Mechanical.DataModel.Enums.GraphicsImageExportFormat.PNG
    settings_720p = Ansys.Mechanical.Graphics.GraphicsImageExportSettings()
    settings_720p.Resolution = (
        Ansys.Mechanical.DataModel.Enums.GraphicsResolutionType.EnhancedResolution
    )
    settings_720p.Background = Ansys.Mechanical.DataModel.Enums.GraphicsBackgroundType.White
    settings_720p.Width = 1280
    settings_720p.Height = 720
    settings_720p.CurrentGraphicsDisplay = False
    """
    )





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    ''



.. GENERATED FROM PYTHON SOURCE LINES 182-185

Import geometry
---------------


.. GENERATED FROM PYTHON SOURCE LINES 185-216

.. code-block:: Python


    mechanical.run_python_script(
        """
    import os
    geometry_import_group = Model.GeometryImportGroup
    geometry_import = geometry_import_group.AddGeometryImport()
    geometry_import_format = (
        Ansys.Mechanical.DataModel.Enums.GeometryImportPreference.Format.Automatic
    )
    geometry_import_preferences = Ansys.ACT.Mechanical.Utilities.GeometryImportPreferences()
    geometry_import_preferences.ProcessNamedSelections = True
    geometry_import_preferences.NamedSelectionKey = ""
    geometry_import_preferences.ProcessMaterialProperties = True
    geometry_import_preferences.ProcessCoordinateSystems = True
    geometry_import.Import(
        geometry_path, geometry_import_format, geometry_import_preferences
    )
    project_directory = ExtAPI.DataModel.Project.ProjectDirectory
    ExtAPI.Graphics.Camera.SetFit()
    ExtAPI.Graphics.ExportImage(
        os.path.join(project_directory, "geometry.png"), image_export_format, settings_720p
    )
    """
    )

    # Download the geometry image and display it
    mechanical.download(files=os.path.join(project_directory, "geometry.png"), target_dir=OUTPUT_DIR)
    if GRAPHICS_BOOL:
        display_image("geometry.png")





.. image-sg:: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_001.png
   :alt: wf fm 02 mechanical
   :srcset: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_001.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/geometry.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/geometry.png:   0%|          | 0.00/86.2k [00:00<?, ?B/s]
    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/geometry.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/geometry.png: 100%|██████████| 86.2k/86.2k [00:00<00:00, 338MB/s]




.. GENERATED FROM PYTHON SOURCE LINES 217-220

Import material, assign it to the bodies and create Named Selections
--------------------------------------------------------------------


.. GENERATED FROM PYTHON SOURCE LINES 220-261

.. code-block:: Python

    mechanical.run_python_script(
        """
    materials = ExtAPI.DataModel.Project.Model.Materials
    materials.Import(material_path)
    materials.RefreshMaterials()

    PRT1 = [x for x in ExtAPI.DataModel.Tree.AllObjects if x.Name == r"Geom-2\Geom-1\solid"][0]

    # Assign it to the bodies

    nmat = "1_HiSi_Model3_Exhaust Manifold updated"
    PRT1.Material = nmat


    # Select MKS units
    ExtAPI.Application.ActiveUnitSystem = MechanicalUnitSystem.StandardNMM

    # Store all main tree nodes as variables
    GEOM = Model.Geometry
    MAT_GRP = Model.Materials

    # Create NS for Named Selection.

    NS_GRP = ExtAPI.DataModel.Project.Model.NamedSelections
    BRACKET_FIX_NS = [x for x in ExtAPI.DataModel.Tree.AllObjects if x.Name == "bracket_fix"][0]
    INTERFACE_SURFACE_NS = [
        x for x in ExtAPI.DataModel.Tree.AllObjects if x.Name == "interface_surface"
    ][0]
    EXHAUST_MANIFOLD_NS = [x for x in ExtAPI.DataModel.Tree.AllObjects if x.Name == "exhaust_manifold"][
        0
    ]
    TOP_BRACKET_SURFACE_NS = [
        x for x in ExtAPI.DataModel.Tree.AllObjects if x.Name == "top_bracket_surface"
    ][0]
    SPACERS_NS = [x for x in ExtAPI.DataModel.Tree.AllObjects if x.Name == "spacers"][0]
    EM_OUTER_SURFACE_NS = [x for x in ExtAPI.DataModel.Tree.AllObjects if x.Name == "em_outer_surface"][
        0
    ]
    """
    )





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/wf_fm_02_mechanical.py:221: SyntaxWarning: invalid escape sequence '\G'
      """

    ''



.. GENERATED FROM PYTHON SOURCE LINES 262-265

Set up the mesh and generate
----------------------------


.. GENERATED FROM PYTHON SOURCE LINES 265-289

.. code-block:: Python

    mechanical.run_python_script(
        """
    MESH = Model.Mesh

    MESH.UseAdaptiveSizing = True
    MESH.TransitionOption = 1

    Tree.Activate([MESH])
    MESH.GenerateMesh()

    # Export mesh image

    ExtAPI.Graphics.Camera.SetFit()
    ExtAPI.Graphics.ExportImage(
        os.path.join(project_directory, "mesh.png"), image_export_format, settings_720p
    )
    """
    )

    # Download the mesh image and display it
    mechanical.download(files=os.path.join(project_directory, "mesh.png"), target_dir=OUTPUT_DIR)
    if GRAPHICS_BOOL:
        display_image("mesh.png")




.. image-sg:: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_002.png
   :alt: wf fm 02 mechanical
   :srcset: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_002.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/mesh.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/mesh.png:   0%|          | 0.00/136k [00:00<?, ?B/s]
    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/mesh.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/mesh.png: 100%|██████████| 136k/136k [00:00<00:00, 480MB/s]




.. GENERATED FROM PYTHON SOURCE LINES 290-293

Add Transient Thermal Analysis and set up the analysis settings
---------------------------------------------------------------


.. GENERATED FROM PYTHON SOURCE LINES 293-481

.. code-block:: Python


    mechanical.run_python_script(
        """
    Model.AddTransientThermalAnalysis()

    # Store all main tree nodes as variables
    TRANS_THERM = Model.Analyses[0]
    TRANS_THERM_SOLN = TRANS_THERM.Solution
    ANA_SETTINGS = TRANS_THERM.Children[1]
    # ANA_SETTINGS = TRANS_THERM.AnalysisSettings

    # Setup transient thermal analysis settings
    ANA_SETTINGS.SolverType = SolverType.Direct
    ANA_SETTINGS.NonLinearFormulation = NonLinearFormulationType.Full

    ANA_SETTINGS.NumberOfSteps = 1
    ANA_SETTINGS.SetStepEndTime(1, Quantity('720[s]'))
    ANA_SETTINGS.NumberOfSteps = 14
    analysis_step = (
        (1, Quantity('1e-3[s]')),
        (2, Quantity('2e-3[s]')),
        (3, Quantity('20[s]')),
        (4, Quantity('30[s]')),
        (5, Quantity('320[s]')),
        (6, Quantity('330[s]')),
        (7, Quantity('350[s]')),
        (8, Quantity('360[s]')),
        (9, Quantity('380[s]')),
        (10, Quantity('390[s]')),
        (11, Quantity('680[s]')),
        (12, Quantity('690[s]')),
        (13, Quantity('710[s]')),
        (14, Quantity('720[s]'))
    )
    for i, q in analysis_step:
        ANA_SETTINGS.SetStepEndTime(i, q)
    ANA_SETTINGS.Activate()

    External_Convection_Load_1 = TRANS_THERM.AddConvection()
    selection = NS_GRP.Children[8]
    External_Convection_Load_1.Location = selection

    External_Convection_Load_1.FilmCoefficient.Inputs[0].DiscreteValues = [
        Quantity('0[s]'), Quantity('1e-3[s]'),
        Quantity('2e-3[s]'), Quantity('20[s]'),
        Quantity('30[s]'), Quantity('320[s]'),
        Quantity('330[s]'), Quantity('350[s]'),
        Quantity('360[s]'), Quantity('380[s]'),
        Quantity('390[s]'), Quantity('680[s]'),
        Quantity('690[s]'), Quantity('710[s]'),
        Quantity('720[s]')
    ]

    External_Convection_Load_1.FilmCoefficient.Output.DiscreteValues = [
        Quantity('60[W m^-1 m^-1 K^-1]'), Quantity('60[W m^-1 m^-1 K^-1]'),
        Quantity('60[W m^-1 m^-1 K^-1]'), Quantity('60[W m^-1 m^-1 K^-1]'),
        Quantity('60[W m^-1 m^-1 K^-1]'), Quantity('60[W m^-1 m^-1 K^-1]'),
        Quantity('60[W m^-1 m^-1 K^-1]'), Quantity('60[W m^-1 m^-1 K^-1]'),
        Quantity('60[W m^-1 m^-1 K^-1]'), Quantity('60[W m^-1 m^-1 K^-1]'),
        Quantity('60[W m^-1 m^-1 K^-1]'), Quantity('60[W m^-1 m^-1 K^-1]'),
        Quantity('60[W m^-1 m^-1 K^-1]'), Quantity('60[W m^-1 m^-1 K^-1]'),
        Quantity('60[W m^-1 m^-1 K^-1]')
    ]

    External_Convection_Load_1.AmbientTemperature.Inputs[0].DiscreteValues = [
        Quantity('0[s]'), Quantity('1e-3[s]'),
        Quantity('2e-3[s]'), Quantity('20[s]'),
        Quantity('30[s]'), Quantity('320[s]'),
        Quantity('330[s]'), Quantity('350[s]'),
        Quantity('360[s]'), Quantity('380[s]'),
        Quantity('390[s]'), Quantity('680[s]'),
        Quantity('690[s]'), Quantity('710[s]'),
        Quantity('720[s]')
    ]

    External_Convection_Load_1.AmbientTemperature.Output.DiscreteValues = [
        Quantity('473.15[K]'), Quantity('473.15[K]'),
        Quantity('473.15[K]'), Quantity('473.15[K]'),
        Quantity('473.15[K]'), Quantity('473.15[K]'),
        Quantity('473.15[K]'), Quantity('473.15[K]'),
        Quantity('473.15[K]'), Quantity('473.15[K]'),
        Quantity('473.15[K]'), Quantity('473.15[K]'),
        Quantity('473.15[K]'), Quantity('473.15[K]'),
        Quantity('473.15[K]')
    ]

    External_Convection_Load_2 = TRANS_THERM.AddConvection()
    selection = NS_GRP.Children[7]
    External_Convection_Load_2.Location = selection

    External_Convection_Load_2.FilmCoefficient.Inputs[0].DiscreteValues = [
        Quantity('0[s]'), Quantity('1e-3[s]'),
        Quantity('2e-3[s]'), Quantity('20[s]'),
        Quantity('30[s]'), Quantity('320[s]'),
        Quantity('330[s]'), Quantity('350[s]'),
        Quantity('360[s]'), Quantity('380[s]'),
        Quantity('390[s]'), Quantity('680[s]'),
        Quantity('690[s]'), Quantity('710[s]'),
        Quantity('720[s]')
    ]

    External_Convection_Load_2.FilmCoefficient.Output.DiscreteValues = [
        Quantity('20[W m^-1 m^-1 K^-1]'), Quantity('20[W m^-1 m^-1 K^-1]'),
        Quantity('20[W m^-1 m^-1 K^-1]'), Quantity('20[W m^-1 m^-1 K^-1]'),
        Quantity('20[W m^-1 m^-1 K^-1]'), Quantity('20[W m^-1 m^-1 K^-1]'),
        Quantity('20[W m^-1 m^-1 K^-1]'), Quantity('20[W m^-1 m^-1 K^-1]'),
        Quantity('20[W m^-1 m^-1 K^-1]'), Quantity('20[W m^-1 m^-1 K^-1]'),
        Quantity('20[W m^-1 m^-1 K^-1]'), Quantity('20[W m^-1 m^-1 K^-1]'),
        Quantity('20[W m^-1 m^-1 K^-1]'), Quantity('20[W m^-1 m^-1 K^-1]'),
        Quantity('20[W m^-1 m^-1 K^-1]')
    ]

    External_Convection_Load_2.AmbientTemperature.Inputs[0].DiscreteValues = [
        Quantity('0[s]'), Quantity('1e-3[s]'),
        Quantity('2e-3[s]'),Quantity('20[s]'),
        Quantity('30[s]'), Quantity('320[s]'),
        Quantity('330[s]'), Quantity('350[s]'),
        Quantity('360[s]'), Quantity('380[s]'),
        Quantity('390[s]'), Quantity('680[s]'),
        Quantity('690[s]'), Quantity('710[s]'),
        Quantity('720[s]')
    ]

    External_Convection_Load_2.AmbientTemperature.Output.DiscreteValues = [
        Quantity('498.15[K]'), Quantity('498.15[K]'),
        Quantity('498.15[K]'), Quantity('498.15[K]'),
        Quantity('498.15[K]'), Quantity('498.15[K]'),
        Quantity('498.15[K]'), Quantity('498.15[K]'),
        Quantity('498.15[K]'), Quantity('498.15[K]'),
        Quantity('498.15[K]'), Quantity('498.15[K]'),
        Quantity('498.15[K]'), Quantity('498.15[K]'),
        Quantity('498.15[K]')
    ]

    External_Convection_Load_3 = TRANS_THERM.AddConvection()
    selection = NS_GRP.Children[6]
    External_Convection_Load_3.Location = selection

    External_Convection_Load_3.FilmCoefficient.Inputs[0].DiscreteValues = [
        Quantity('0[s]'), Quantity('1e-3[s]'),
        Quantity('2e-3[s]'), Quantity('20[s]'),
        Quantity('30[s]'), Quantity('320[s]'),
        Quantity('330[s]'), Quantity('350[s]'),
        Quantity('360[s]'), Quantity('380[s]'),
        Quantity('390[s]'), Quantity('680[s]'),
        Quantity('690[s]'), Quantity('710[s]'),
        Quantity('720[s]')
    ]

    External_Convection_Load_3.FilmCoefficient.Output.DiscreteValues = [
        Quantity('500[W m^-1 m^-1 K^-1]'), Quantity('500[W m^-1 m^-1 K^-1]'),
        Quantity('500[W m^-1 m^-1 K^-1]'), Quantity('500[W m^-1 m^-1 K^-1]'),
        Quantity('500[W m^-1 m^-1 K^-1]'), Quantity('500[W m^-1 m^-1 K^-1]'),
        Quantity('500[W m^-1 m^-1 K^-1]'), Quantity('500[W m^-1 m^-1 K^-1]'),
        Quantity('500[W m^-1 m^-1 K^-1]'), Quantity('500[W m^-1 m^-1 K^-1]'),
        Quantity('500[W m^-1 m^-1 K^-1]'), Quantity('500[W m^-1 m^-1 K^-1]'),
        Quantity('500[W m^-1 m^-1 K^-1]'), Quantity('500[W m^-1 m^-1 K^-1]'),
        Quantity('500[W m^-1 m^-1 K^-1]')
    ]

    External_Convection_Load_3.AmbientTemperature.Inputs[0].DiscreteValues = [
        Quantity('0[s]'), Quantity('1e-3[s]'),
        Quantity('2e-3[s]'), Quantity('20[s]'),
        Quantity('30[s]'), Quantity('320[s]'),
        Quantity('330[s]'), Quantity('350[s]'),
        Quantity('360[s]'), Quantity('380[s]'),
        Quantity('390[s]'), Quantity('680[s]'),
        Quantity('690[s]'), Quantity('710[s]'),
        Quantity('720[s]')
    ]

    External_Convection_Load_3.AmbientTemperature.Output.DiscreteValues = [
        Quantity('373.15[K]'), Quantity('373.15[K]'),
        Quantity('373.15[K]'), Quantity('373.15[K]'),
        Quantity('373.15[K]'), Quantity('373.15[K]'),
        Quantity('373.15[K]'), Quantity('373.15[K]'),
        Quantity('373.15[K]'), Quantity('373.15[K]'),
        Quantity('373.15[K]'), Quantity('373.15[K]'),
        Quantity('373.15[K]'), Quantity('373.15[K]'),
        Quantity('373.15[K]')
    ]

    group_list = [External_Convection_Load_1, External_Convection_Load_2, External_Convection_Load_3]
    grouping_folder = Tree.Group(group_list)
    tree_grouping_folder_70 = DataModel.GetObjectsByName("New Folder")
    """
    )





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    ''



.. GENERATED FROM PYTHON SOURCE LINES 482-486

Use the output from Fluent to import the temperature and HTC data
-----------------------------------------------------------------

Add imported convection

.. GENERATED FROM PYTHON SOURCE LINES 486-661

.. code-block:: Python

    result = mechanical.run_python_script(
        """
    Imported_Load_Group = TRANS_THERM.AddImportedLoadExternalData()
    imported_load_group_61=Imported_Load_Group
    imported_convection_62 = Imported_Load_Group.AddImportedConvection()
    """
    )
    result = mechanical.run_python_script(
        """
    external_data_files = Ansys.Mechanical.ExternalData.ExternalDataFileCollection()
    external_data_files.SaveFilesWithProject = False

    external_data_file_1 = Ansys.Mechanical.ExternalData.ExternalDataFile()
    external_data_files.Add(external_data_file_1)
    external_data_file_1.Identifier = "File1"
    external_data_file_1.Description = "High"
    external_data_file_1.IsMainFile = True
    external_data_file_1.FilePath = temp_htc_data_high_path
    external_data_file_1.ImportSettings = (
        Ansys.Mechanical.ExternalData.ImportSettingsFactory.GetSettingsForFormat(
            MechanicalEnums.ExternalData.ImportFormat.Delimited
        )
    )
    import_settings = external_data_file_1.ImportSettings
    import_settings.SkipRows = 1
    import_settings.SkipFooter = 0
    import_settings.Delimiter = ","
    import_settings.AverageCornerNodesToMidsideNodes = False
    import_settings.UseColumn(
        0, MechanicalEnums.ExternalData.VariableType.NodeId, "", "Node ID@A"
    )
    import_settings.UseColumn(
        1, MechanicalEnums.ExternalData.VariableType.XCoordinate, "m", "X Coordinate@B"
    )
    import_settings.UseColumn(
        2, MechanicalEnums.ExternalData.VariableType.YCoordinate, "m", "Y Coordinate@C"
    )
    import_settings.UseColumn(
        3, MechanicalEnums.ExternalData.VariableType.ZCoordinate, "m", "Z Coordinate@D"
    )
    import_settings.UseColumn(
        4, MechanicalEnums.ExternalData.VariableType.Temperature, "K", "Temperature@E"
    )
    import_settings.UseColumn(
        5, MechanicalEnums.ExternalData.VariableType.HeatTransferCoefficient,
        "W m^-2 K^-1", "Heat Transfer Coefficient@F"
    )

    external_data_file_2 = Ansys.Mechanical.ExternalData.ExternalDataFile()
    external_data_files.Add(external_data_file_2)
    external_data_file_2.Identifier = "File2"
    external_data_file_2.Description = "Med"
    external_data_file_2.IsMainFile = False
    external_data_file_2.FilePath = temp_htc_data_med_path
    external_data_file_2.ImportSettings = (
        Ansys.Mechanical.ExternalData.ImportSettingsFactory.GetSettingsForFormat(
            MechanicalEnums.ExternalData.ImportFormat.Delimited
        )
    )
    import_settings = external_data_file_2.ImportSettings
    import_settings.SkipRows = 1
    import_settings.SkipFooter = 0
    import_settings.Delimiter = ","
    import_settings.AverageCornerNodesToMidsideNodes = False
    import_settings.UseColumn(
        0, MechanicalEnums.ExternalData.VariableType.NodeId, "", "Node ID@A"
    )
    import_settings.UseColumn(
        1, MechanicalEnums.ExternalData.VariableType.XCoordinate, "m", "X Coordinate@B"
    )
    import_settings.UseColumn(
        2, MechanicalEnums.ExternalData.VariableType.YCoordinate, "m", "Y Coordinate@C"
    )
    import_settings.UseColumn(
        3, MechanicalEnums.ExternalData.VariableType.ZCoordinate, "m", "Z Coordinate@D"
    )
    import_settings.UseColumn(
        4, MechanicalEnums.ExternalData.VariableType.Temperature, "K", "Temperature@E"
    )
    import_settings.UseColumn(
        5, MechanicalEnums.ExternalData.VariableType.HeatTransferCoefficient,
        "W m^-2 K^-1", "Heat Transfer Coefficient@F"
    )

    external_data_file_3 = Ansys.Mechanical.ExternalData.ExternalDataFile()
    external_data_files.Add(external_data_file_3)
    external_data_file_3.Identifier = "File3"
    external_data_file_3.Description = "Low"
    external_data_file_3.IsMainFile = False
    external_data_file_3.FilePath = temp_htc_data_low_path
    external_data_file_3.ImportSettings = (
        Ansys.Mechanical.ExternalData.ImportSettingsFactory.GetSettingsForFormat(
            MechanicalEnums.ExternalData.ImportFormat.Delimited
        )
    )
    import_settings = external_data_file_3.ImportSettings
    import_settings.SkipRows = 1
    import_settings.SkipFooter = 0
    import_settings.Delimiter = ","
    import_settings.AverageCornerNodesToMidsideNodes = False
    import_settings.UseColumn(
        0, MechanicalEnums.ExternalData.VariableType.NodeId, "", "Node ID@A"
    )
    import_settings.UseColumn(
        1, MechanicalEnums.ExternalData.VariableType.XCoordinate, "m", "X Coordinate@B"
    )
    import_settings.UseColumn(
        2, MechanicalEnums.ExternalData.VariableType.YCoordinate, "m", "Y Coordinate@C"
    )
    import_settings.UseColumn(
        3, MechanicalEnums.ExternalData.VariableType.ZCoordinate, "m", "Z Coordinate@D"
    )
    import_settings.UseColumn(
        4, MechanicalEnums.ExternalData.VariableType.Temperature, "K", "Temperature@E"
    )
    import_settings.UseColumn(
        5, MechanicalEnums.ExternalData.VariableType.HeatTransferCoefficient,
        "W m^-2 K^-1", "Heat Transfer Coefficient@F"
    )

    imported_load_group_61.ImportExternalDataFiles(external_data_files)
    """
    )

    result = mechanical.run_python_script(
        """
    table = imported_load_group_61.Children[0].GetTableByName("Film Coefficient")
    numofsteps = 15
    Film_Coeff = [
        "File1:Heat Transfer Coefficient@F",
        "File2:Heat Transfer Coefficient@F",
        "File3:Heat Transfer Coefficient@F"
    ]
    Amb_Temp = [
        "File1:Temperature@E",
        "File2:Temperature@E",
        "File3:Temperature@E"
    ]
    Ana_time = [
        "0", "1e-3", "2e-3", "20", "30", "320", "330", "350", "360", "380", "390",
        "680", "690", "710", "720"
    ]

    for i in range(numofsteps - 1):
        table.Add(None)

    for i in range(numofsteps):
        table[i][0] = Film_Coeff[i % 3]
        table[i][1] = Amb_Temp[i % 3]
        table[i][2] = Ana_time[i]

    selection = NS_GRP.Children[4]
    imported_convection_62.Location = selection
    imported_load_id = imported_convection_62.ObjectId
    imported_load = DataModel.GetObjectById(imported_load_id)
    """
    )

    mechanical.run_python_script(
        """
    imported_load.ImportLoad()

    Tree.Activate([imported_load])
    ExtAPI.Graphics.Camera.SetFit()
    ExtAPI.Graphics.ExportImage(
        os.path.join(project_directory, "imported_temperature.png"), image_export_format, settings_720p
    )
    """
    )
    mechanical.download(
        files=os.path.join(project_directory, "imported_temperature.png"), target_dir=OUTPUT_DIR
    )
    if GRAPHICS_BOOL:
        display_image("imported_temperature.png")




.. image-sg:: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_003.png
   :alt: wf fm 02 mechanical
   :srcset: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_003.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/imported_temperature.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/imported_temperature.png:   0%|          | 0.00/113k [00:00<?, ?B/s]
    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/imported_temperature.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/imported_temperature.png: 100%|██████████| 113k/113k [00:00<00:00, 453MB/s]




.. GENERATED FROM PYTHON SOURCE LINES 662-665

Solve and post-process the results
----------------------------------


.. GENERATED FROM PYTHON SOURCE LINES 665-695

.. code-block:: Python

    mechanical.run_python_script(
        """
    # Insert results objects

    Temp = TRANS_THERM_SOLN.AddTemperature()
    Temp.DisplayTime = Quantity("680 [s]")

    # Run Solution: Transient Thermal Simulation

    TRANS_THERM_SOLN.Solve(True)
    TRANS_THERM_SS = TRANS_THERM_SOLN.Status

    # Export temperature image

    Tree.Activate([Temp])
    ExtAPI.Graphics.ViewOptions.ResultPreference.ExtraModelDisplay = (
        Ansys.Mechanical.DataModel.MechanicalEnums.Graphics.ExtraModelDisplay.NoWireframe
    )
    ExtAPI.Graphics.ExportImage(
        os.path.join(project_directory, "temperature.png"), image_export_format, settings_720p
    )
    """
    )

    # Download the temperature image and display it
    mechanical.download(files=os.path.join(project_directory, "temperature.png"), target_dir=OUTPUT_DIR)
    if GRAPHICS_BOOL:
        display_image("temperature.png")





.. image-sg:: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_004.png
   :alt: wf fm 02 mechanical
   :srcset: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_004.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/temperature.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/temperature.png:   0%|          | 0.00/103k [00:00<?, ?B/s]
    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/temperature.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/temperature.png: 100%|██████████| 103k/103k [00:00<00:00, 430MB/s]




.. GENERATED FROM PYTHON SOURCE LINES 696-699

Setup Structural Analysis
-------------------------


.. GENERATED FROM PYTHON SOURCE LINES 699-762

.. code-block:: Python

    mechanical.run_python_script(
        """
    Model.AddStaticStructuralAnalysis()

    # Define analysis settings

    # Setup static structural analysis settings
    STAT_STRUC = Model.Analyses[1]
    STAT_STRUC_SOLN = STAT_STRUC.Solution
    STAT_STRUC_ANA_SETTING = STAT_STRUC.Children[0]

    STAT_STRUC_ANA_SETTING.NumberOfSteps = 1
    STAT_STRUC_ANA_SETTING.SetStepEndTime(1, Quantity('720[s]'))
    STAT_STRUC_ANA_SETTING.NumberOfSteps = 14

    analysis_step = (
        (1, Quantity('1e-3[s]')),
        (2, Quantity('2e-3[s]')),
        (3, Quantity('20[s]')),
        (4, Quantity('30[s]')),
        (5, Quantity('320[s]')),
        (6, Quantity('330[s]')),
        (7, Quantity('350[s]')),
        (8, Quantity('360[s]')),
        (9, Quantity('380[s]')),
        (10, Quantity('390[s]')),
        (11, Quantity('680[s]')),
        (12, Quantity('690[s]')),
        (13, Quantity('710[s]')),
        (14, Quantity('720[s]'))
    )
    for i, q in analysis_step:
        STAT_STRUC_ANA_SETTING.SetStepEndTime(i,q)
    STAT_STRUC_ANA_SETTING.Activate()


    # Add Imported Body Temperature load from Transient Thermal Run

    STAT_STRUC.ImportLoad(Model.Analyses[0])
    imported_load = DataModel.GetObjectsByName("Imported Body Temperature")[0]

    table = imported_load.GetTableByName("Source Time")
    numofsteps = 14
    nCol = 2
    Ana_time = ["1e-3","2e-3","20","30","320","330","350","360","380","390","680","690","710","720"]

    for i in range(numofsteps-1):
        table.Add(None)

    for i in range(numofsteps):
        for j in range(nCol):
            table[i][j] = Ana_time[i]

    imported_load.ImportLoad()

    # Apply Fixed Support Condition

    Fixed_Support = STAT_STRUC.AddFixedSupport()
    selection = NS_GRP.Children[3]
    Fixed_Support.Location = selection
    """
    )





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    ''



.. GENERATED FROM PYTHON SOURCE LINES 763-766

Solve and post-process the results
----------------------------------


.. GENERATED FROM PYTHON SOURCE LINES 766-836

.. code-block:: Python

    mechanical.run_python_script(
        """
    SOLN = STAT_STRUC.Solution

    TOT_DEF1 = SOLN.AddTotalDeformation()
    TOT_DEF1.DisplayTime = Quantity("680 [s]")

    EQV_STRS1 = SOLN.AddEquivalentStress()
    EQV_STRS1.DisplayTime = Quantity("680 [s]")

    EQV_PLAS_STRN1 = SOLN.AddEquivalentPlasticStrain()
    EQV_PLAS_STRN1.DisplayTime = Quantity("680 [s]")

    THERM_STRN1 = SOLN.AddThermalStrain()
    THERM_STRN1.DisplayTime = Quantity("680 [s]")

    # Solve Nonlinear Static Simulation

    SOLN.Solve(True)
    STAT_STRUC_SS = SOLN.Status

    # Export results images

    Tree.Activate([TOT_DEF1])
    ExtAPI.Graphics.ViewOptions.ResultPreference.ExtraModelDisplay = (
        Ansys.Mechanical.DataModel.MechanicalEnums.Graphics.ExtraModelDisplay.NoWireframe
    )
    ExtAPI.Graphics.ExportImage(
        os.path.join(project_directory, "deformation.png"), image_export_format, settings_720p
    )

    Tree.Activate([EQV_STRS1])
    ExtAPI.Graphics.ExportImage(
        os.path.join(project_directory, "stress.png"), image_export_format, settings_720p
    )

    Tree.Activate([EQV_PLAS_STRN1])
    ExtAPI.Graphics.ExportImage(
        os.path.join(project_directory, "plastic_strain.png"), image_export_format, settings_720p
    )
    """
    )

    # Download the results images to local directory
    mechanical.download(files=os.path.join(project_directory, "deformation.png"), target_dir=OUTPUT_DIR)
    mechanical.download(files=os.path.join(project_directory, "stress.png"), target_dir=OUTPUT_DIR)
    mechanical.download(
        files=os.path.join(project_directory, "plastic_strain.png"), target_dir=OUTPUT_DIR
    )

    # Deformation
    if GRAPHICS_BOOL:
        display_image("deformation.png")


    # Stress
    if GRAPHICS_BOOL:
        display_image("stress.png")


    # Plastic strain
    if GRAPHICS_BOOL:
        display_image("plastic_strain.png")


    # ###############################################################################
    # Close the Mechanical
    # --------------------
    #
    mechanical.exit()



.. rst-class:: sphx-glr-horizontal


    *

      .. image-sg:: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_005.png
         :alt: wf fm 02 mechanical
         :srcset: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_005.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_006.png
         :alt: wf fm 02 mechanical
         :srcset: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_006.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_007.png
         :alt: wf fm 02 mechanical
         :srcset: /examples/fluent-mechanical/images/sphx_glr_wf_fm_02_mechanical_007.png
         :class: sphx-glr-multi-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/deformation.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/deformation.png:   0%|          | 0.00/126k [00:00<?, ?B/s]
    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/deformation.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/deformation.png: 100%|██████████| 126k/126k [00:00<00:00, 503MB/s]

    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/stress.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/stress.png:   0%|          | 0.00/146k [00:00<?, ?B/s]
    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/stress.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/stress.png: 100%|██████████| 146k/146k [00:00<00:00, 626MB/s]

    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/plastic_strain.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/plastic_strain.png:   0%|          | 0.00/132k [00:00<?, ?B/s]
    Downloading dns:///127.0.0.1:10000:/tmp/ANSYS.root.1/AnsysMech259F/Project_Mech_Files/plastic_strain.png to /home/runner/work/pyansys-workflows/pyansys-workflows/fluent-mechanical/outputs/plastic_strain.png: 100%|██████████| 132k/132k [00:00<00:00, 619MB/s]





.. rst-class:: sphx-glr-timing

   **Total running time of the script:** (7 minutes 23.728 seconds)


.. _sphx_glr_download_examples_fluent-mechanical_wf_fm_02_mechanical.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: wf_fm_02_mechanical.ipynb <wf_fm_02_mechanical.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: wf_fm_02_mechanical.py <wf_fm_02_mechanical.py>`

    .. container:: sphx-glr-download sphx-glr-download-zip

      :download:`Download zipped: wf_fm_02_mechanical.zip <wf_fm_02_mechanical.zip>`