nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

Mechanical scripting - Pipe model pre/post-processing

gdebennetot

Objective:

This post is part of a series of post that will detail each step to run a completely automated Workbench Workflow. The scripts are the following:
1. Workbench python script & the batch command to launch this script from a single command line: https://discuss.ansys.com/discussion/10615/scripting-a-complete-workbench-workflow-from-geometry-to-post-processing?new=1
2. Discovery python script to create a geometry from scratch: https://discuss.ansys.com/discussion/10616/discovery-geometry-scripting-pipe-geometry-creation-from-scratch?new=1
3. Mechanical python script to pre and post-process the static structural model - Current post

Mechanical Script: Pre process, solve and post-process the pipe geometry model

In Mechanical, the script editor can be found in the Automation tab --> "Scripting, and you can record yourself to see what commands correspond to any user action.

Credits to Vishnu Venkataraman for this script:
code split into 2 comments below due to size restrictions

Find more posts tagged with

Mechanical

Ansys Workbench

Comments

gdebennetot

PART 1:

import re

def nsbyname(myname):
    ns = ExtAPI.DataModel.Project.Model.GetChildren(DataModelObjectCategory.NamedSelection,True)
    my_ns  = filter(lambda x: re.search(myname,x.Name, re.IGNORECASE), ns)
    return my_ns
    
def rpbyname(myname):
    rps = ExtAPI.DataModel.Project.Model.GetChildren(DataModelObjectCategory.RemotePoint,True)
    my_rp  = filter(lambda x: re.search(myname,x.Name, re.IGNORECASE), rps)
    return my_rp


def csbyname(myname):
    rps = ExtAPI.DataModel.Project.Model.GetChildren(DataModelObjectCategory.CoordinateSystem,True)
    my_cs  = filter(lambda x: re.search(myname,x.Name, re.IGNORECASE), rps)
    return my_cs
    
def midpoint(P1,P2):
    P3 = []
    for x1,x2 in zip(P1,P2):
        P3.Add((x1+x2)/2.0)
    return P3

def GetFacesOfSameSize(nsName):
    nsHeader = nsbyname(nsName)[0]
    headerArea = round(DataModel.GeoData.GeoEntityById(nsHeader.Ids[0]).Area,5)
    
    headerBody = DataModel.GeoData.GeoEntityById(nsHeader.Ids[0]).Bodies[0]

    matchingFaces = []
    for face in headerBody.Faces:
      if round(face.Area,5) == headerArea:
          matchingFaces.append(face.Id)
    return matchingFaces
   

def GetFacesOfSameSizeInAllGeos(nsName):
    nsHeader = nsbyname(nsName)[0]
    headerArea = round(DataModel.GeoData.GeoEntityById(nsHeader.Ids[0]).Area,5)
    
    g = ExtAPI.DataModel.Project.Model.Geometry
    allbodies = g.GetChildren[Ansys.ACT.Automation.Mechanical.Body](True)
    allGeoBodies = [x.GetGeoBody() for x in allbodies]

    matchingFaces = []
    for body in allGeoBodies:
        for face in body.Faces:
          if round(face.Area,5) == headerArea:
              matchingFaces.append(face.Id)
    return matchingFaces
rp = ExtAPI.DataModel.Project.Model.AddRemotePoint()
rp.Location=nsbyname('outlet')[0]
rp.Name = 'Remote Point Outlet'
rp = ExtAPI.DataModel.Project.Model.AddRemotePoint()
rp.Location=nsbyname('header')[0]
rp.Name = 'Remote Point Header'
rp = ExtAPI.DataModel.Project.Model.AddRemotePoint()
rp.Location=nsbyname('branch')[0]
rp.Name = 'Remote Point Branch'

# Mesh sizing
mesh = ExtAPI.DataModel.Project.Model.Mesh
mesh.ElementOrder = ElementOrder.Linear
mesh.ElementSize = Quantity(30,"mm")
mesh.UseAdaptiveSizing = False

sizing = ExtAPI.DataModel.Project.Model.Mesh.AddSizing()
sizing.Location = nsbyname('main_pipe_id_od')[0]
sizing.ElementSize = Quantity('15 [mm]')
sizing.Name = 'Element Size on main_pipe_ID_OD'

sizing = ExtAPI.DataModel.Project.Model.Mesh.AddSizing()
sizing.Location = nsbyname('branch_pipe_id_od')[0]
sizing.ElementSize = Quantity('5.0 [mm]')
sizing.InternalObject.UseProximitySize = True
sizing.Name = 'Element Size on branch_pie_ID_OD'

sizing = ExtAPI.DataModel.Project.Model.Mesh.AddSizing()
sizing.Location = nsbyname('main_t_section')[0]
sizing.Type=SizingType.SphereOfInfluence

# Coordinate system for mesh sizing object
mainSectionNS = nsbyname("main_t_section")[0]
mainSectionBody = DataModel.GeoData.GeoEntityById(mainSectionNS.Ids[0])
centroidMainSectionBody = mainSectionBody.Centroid

mainBranchNS = nsbyname("Branch")[0]
mainBranchBody = DataModel.GeoData.GeoEntityById(mainBranchNS.Ids[0]).Bodies[0]
centroidBranchBody = mainBranchBody.Centroid

csMidPoint = midpoint(centroidMainSectionBody,centroidBranchBody)
csCoordinates = [Quantity(x, "m") for x in csMidPoint]

cs = DataModel.Project.Model.CoordinateSystems.AddCoordinateSystem()
cs.OriginDefineBy =  CoordinateSystemAlignmentType.Free
cs.OriginX = csCoordinates[0]
cs.OriginY = csCoordinates[1]
cs.OriginZ = csCoordinates[2]
cs.Name = "Coordinate system for mesh sizing"

sizing.SphereCenter = cs
sizing.SphereRadius = Quantity('60.0 [mm]')
sizing.ElementSize = Quantity('5 [mm]')

sizing.Name = 'Sphere of Influence on main_t_section'  
fm = ExtAPI.DataModel.Project.Model.Mesh.AddFaceMeshing()
fm.Location = nsbyname('main_pipe_id_od')[0]
fm.Name = 'Face Meshing on main_pipe_ID_OD'
fm = ExtAPI.DataModel.Project.Model.Mesh.AddFaceMeshing()
fm.Location = nsbyname('branch_pipe_id_od')[0]
fm.Name = 'Face Meshing on main_pipe_ID_OD 2'

fm = ExtAPI.DataModel.Project.Model.Mesh.AddFaceMeshing()
fm.InternalNumberOfDivisions = 2

facesForMappedMesh = []
facesForMappedMesh += GetFacesOfSameSizeInAllGeos("header")
facesForMappedMesh += GetFacesOfSameSizeInAllGeos("Branch")
sel = ExtAPI.SelectionManager.CreateSelectionInfo(SelectionTypeEnum.GeometryEntities)
sel.Ids = facesForMappedMesh
if len(facesForMappedMesh) == 0:
    fm.Suppressed = True
else:
    fm.Location = sel
fm.Name = 'Face Meshing on join_faces'

cs = ExtAPI.DataModel.Project.Model.CoordinateSystems.AddCoordinateSystem()
cs.CoordinateSystemType = CoordinateSystemTypeEnum.Cylindrical
cs.OriginDefineBy = CoordinateSystemAlignmentType.Component
cs.OriginLocation = nsbyname('header')[0]
cs.PrimaryAxis = CoordinateSystemAxisType.PositiveZAxis
cs.Name = 'Pipe Cylindrical CS'

cs = ExtAPI.DataModel.Project.Model.CoordinateSystems.AddCoordinateSystem()
cs.CoordinateSystemType = CoordinateSystemTypeEnum.Cylindrical
cs.OriginDefineBy = CoordinateSystemAlignmentType.Component
cs.OriginLocation = nsbyname('branch')[0]
cs.PrimaryAxis = CoordinateSystemAxisType.PositiveZAxis
cs.PrimaryAxisDefineBy = CoordinateSystemAlignmentType.GlobalY
cs.Name = 'Branch Cylindrical CS'

static_ana = DataModel.Project.Model.Analyses[0]
rf = static_ana.AddRemoteForce()
rf.InternalObject.GeometryDefineBy =  GeometryDefineByType.RemotePoint
rf.Location = rpbyname('Remote Point Branch')[0]
rf.DefineBy=LoadDefineBy.Components
rf.XComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.
rf.XComponent.Output.DiscreteValues = [Quantity("0 [N]"), Quantity("50 [N]")]  # Set the magnitudes.
rf.YComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.
rf.YComponent.Output.DiscreteValues = [Quantity("0 [N]"), Quantity("-20 [N]")]  # Set the magnitudes.
rf.ZComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.
rf.ZComponent.Output.DiscreteValues = [Quantity("0 [N]"), Quantity("35 [N]")]  # Set the magnitudes.
rf.Name = 'Branch Force'

rf = static_ana.AddRemoteForce()
rf.InternalObject.GeometryDefineBy =  GeometryDefineByType.RemotePoint
rf.Location = rpbyname('Remote Point Header')[0]
rf.DefineBy=LoadDefineBy.Components
rf.XComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.
rf.XComponent.Output.DiscreteValues = [Quantity("0 [N]"), Quantity("200 [N]")]  # Set the magnitudes.
rf.YComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.
rf.YComponent.Output.DiscreteValues = [Quantity("0 [N]"), Quantity("100 [N]")]  # Set the magnitudes.
rf.ZComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.
rf.ZComponent.Output.DiscreteValues = [Quantity("0 [N]"), Quantity("500 [N]")]  # Set the magnitudes.
rf.Name = 'Header Force'

rm = static_ana.AddMoment()
rm.InternalObject.GeometryDefineBy =  GeometryDefineByType.RemotePoint
rm.Location = rpbyname('Remote Point Branch')[0]
rm.DefineBy=LoadDefineBy.Components
rm.XComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.
rm.XComponent.Output.DiscreteValues = [Quantity("0 [N mm]"), Quantity("10 [N mm]")]  # Set the magnitudes.
rm.YComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.
rm.YComponent.Output.DiscreteValues = [Quantity("0 [N mm]"), Quantity("10 [N mm]")]  # Set the magnitudes.
rm.ZComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.
rm.ZComponent.Output.DiscreteValues = [Quantity("0 [N mm]"), Quantity("10 [N mm]")]  # Set the magnitudes.
rm.Name = 'Branch Moment'

gdebennetot

PART 2:

rm = static_ana.AddMoment()

rm.InternalObject.GeometryDefineBy =  GeometryDefineByType.RemotePoint

rm.Location = rpbyname('Remote Point Header')[0]

rm.DefineBy=LoadDefineBy.Components

rm.XComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

rm.XComponent.Output.DiscreteValues = [Quantity("0 [N mm]"), Quantity("-10 [N mm]")]  # Set the magnitudes.

rm.YComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

rm.YComponent.Output.DiscreteValues = [Quantity("0 [N mm]"), Quantity("-10 [N mm]")]  # Set the magnitudes.

rm.ZComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

rm.ZComponent.Output.DiscreteValues = [Quantity("0 [N mm]"), Quantity("-10 [N mm]")]  # Set the magnitudes.

rm.Name = 'Header Moment'


pressure = static_ana.AddPressure() 

pressure.Location = nsbyname('branch')[0] 

pressure.Magnitude.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

pressure.Magnitude.Output.DiscreteValues = [Quantity("0 [MPa]"), Quantity("2 [MPa]")]  # Set the magnitudes.

pressure.Name = 'Branch End Cap Pressure'


pressure = static_ana.AddPressure() 

pressure.Location = nsbyname('header')[0]

pressure.Magnitude.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

pressure.Magnitude.Output.DiscreteValues = [Quantity("0 [MPa]"), Quantity("1 [MPa]")]  # Set the magnitudes.

pressure.Name = 'Header End Cap Pressure'


pressure = static_ana.AddPressure() 

pressure.Location = nsbyname('differentialpressure')[0]

pressure.Magnitude.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

pressure.Magnitude.Output.DiscreteValues = [Quantity("0 [MPa]"), Quantity("0.5 [MPa]")]  # Set the magnitudes.

pressure.Name = 'Differential Pressure'


tc = static_ana.AddThermalCondition() 

tc.Location = nsbyname('all_bodies')[0]

tc.Magnitude.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

tc.Magnitude.Output.DiscreteValues = [Quantity("0 [C]"), Quantity("200 [C]")]  # Set the magnitudes.

tc.Name = 'Thermal Condition'


disp = static_ana.AddDisplacement() 

disp.Location = nsbyname('header')[0]

disp.YComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

disp.YComponent.Output.DiscreteValues = [Quantity("0 [mm]"), Quantity("0 [mm]")]  # Set the magnitudes.

disp.ZComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

disp.ZComponent.Output.DiscreteValues = [Quantity("0 [mm]"), Quantity("0 [mm]")]  # Set the magnitudes.

disp.Name = 'header bc face'

disp.CoordinateSystem = csbyname('Pipe Cylindrical CS')[0]


disp = static_ana.AddDisplacement() 

disp.Location = nsbyname('outlet')[0]

disp.XComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

disp.XComponent.Output.DiscreteValues = [Quantity("0 [mm]"), Quantity("0 [mm]")]  # Set the magnitudes.

disp.YComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

disp.YComponent.Output.DiscreteValues = [Quantity("0 [mm]"), Quantity("0 [mm]")]  # Set the magnitudes.

disp.ZComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

disp.ZComponent.Output.DiscreteValues = [Quantity("0 [mm]"), Quantity("0 [mm]")]  # Set the magnitudes.

disp.Name = 'outlet bc face'

disp.CoordinateSystem = csbyname('Pipe Cylindrical CS')[0]


disp = static_ana.AddDisplacement() 

disp.Location = nsbyname('branch')[0]

disp.YComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

disp.YComponent.Output.DiscreteValues = [Quantity("0 [mm]"), Quantity("0 [mm]")]  # Set the magnitudes.

disp.ZComponent.Inputs[0].DiscreteValues = [Quantity("0 [s]"), Quantity("1 [s]")]  # Set the time values.

disp.ZComponent.Output.DiscreteValues = [Quantity("0 [mm]"), Quantity("0 [mm]")]  # Set the magnitudes.

disp.Name = 'branch bc face'

disp.CoordinateSystem = csbyname('Branch Cylindrical CS')[0]


import context_menu

context_menu.DoGraphicsZoom2Selection(ExtAPI)


static_ana.Solve()


static_ana.Solution.AddEquivalentStress()

static_ana.Solution.AddTotalDeformation()

ud = static_ana.Solution.AddUserDefinedResult()

ud.Location = nsbyname('main_t_section')[0]

ud.Expression = 'S1+S2+S3'

ud.Name = 'Triaxial Stress'


pr = static_ana.Solution.AddPythonResult()

pr.EvaluateAllResults()