Discovery geometry scripting - Pipe geometry creation from scratch

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 - Current post
3. Mechanical python script to pre and post-process the static structural model: https://discuss.ansys.com/discussion/10617/mechanical-scripting-pipe-model-pre-post-processing?new=1

Discovery script: create a simple pipe geometry from scratch

In Discovery, the script editor can be found in the Design tab --> "Script Editor", and you can record yourself to see what commands correspond to any user action.
Note: you can change the way object selection is written to script in the dropdown of the "Insert selection" button. It is usually more useful to switch to "Index" so the selection logic is based on Tree items:

The following python code has been written for version 25R2:
see comment below

gdebennetot

The code:

ClearAll()


#USER INPUTS

#Right pipe

Z1_1=400

Z1_2=150

rint1=90

rout1=100


#T pipe

#	right section

Z2_1=150

Z2_2=80

rint2=90

rout2=100


#	middle section

Z3_1=80

Z3_2=-80

rint3=90

rout3=110

#	Edge fillet at section change

R_EF1=12


#Branch pipe

X4_1=200

X4_2=100

rout4=25

rint4=20

R_EF_4=15

#END OF USER INPUT


def create_circular_face(Z,rint,rout,name):

    origin_point=Point.Create(MM(0),MM(0),MM(Z))

    frame = Frame.Create(origin_point, Direction.DirX,Direction.DirY)

    circle_int = Circle.Create(frame,MM(rint))

    circle_ext = Circle.Create(frame,MM(rout))

    SketchCircle.Create(circle_int)

    SketchCircle.Create(circle_ext)

    result = ViewHelper.SetViewMode(InteractionMode.Solid)

    face=result.CreatedBodies[0].Faces[0]


    face_selection=Selection.Create(face)

    if not name is None:

        NamedSelection.Create(face_selection,Selection.Empty(),name)

    return face


def create_Zpipe(Z1,Z2,rint,rout,fname1,fname2,pipe_name):

    face1=create_circular_face(Z1,rint,rout,fname1)

    face2=create_circular_face(Z2,rint,rout,fname2)

    # Create Blend

    selection = FaceSelection.Create(face1, face2)

    options = LoftOptions()

    options.ExtrudeType = ExtrudeType.ForceAdd

    result = Loft.Create(selection, None, options)

    result.CreatedBodies[0].Name=pipe_name

    return result.CreatedBodies[0]


# Extrude main pipe

main_pipe_R=create_Zpipe(Z1_1,Z1_2,rint1,rout1,"header",None,"Main pipe right")

ViewHelper.SetObjectVisibility(Selection.Create(main_pipe_R), VisibilityType.Hide, inSelectedView=False, faceLevel=False)


main_pipe_L=create_Zpipe(-Z1_1,-Z1_2,rint1,rout1,"outlet",None,"Main pipe left")

ViewHelper.SetObjectVisibility(Selection.Create(main_pipe_L), VisibilityType.Hide, inSelectedView=False, faceLevel=False)


#T pipe creation

#right section

T_pipe=create_Zpipe(Z2_1,Z2_2,rint2,rout2,None,None,"T_pipe")

#pipe=Selection.Create(T_pipe)

ViewHelper.SetObjectVisibility(Selection.Create(T_pipe), VisibilityType.Hide, inSelectedView=False, faceLevel=False)


#middle section

T_pipe1=create_Zpipe(Z3_1,Z3_2,rint3,rout3,None,None,"T_pipe1")

pipe1=Selection.Create(T_pipe1)

faces=pipe1.ConvertToFaces()

max_face=faces.FilterMaxArea()

NamedSelection.Create(max_face,Selection.Empty(),"OD_main_pipe")

ViewHelper.SetObjectVisibility(pipe1, VisibilityType.Hide, inSelectedView=False, faceLevel=False)


#left section

T_pipe2=create_Zpipe(-Z2_1,-Z2_2,rint2,rout2,None,None,"T_pipe2")

pipe2=Selection.Create(T_pipe2)


#Combine all 3 sections

ViewHelper.SetObjectVisibility(Selection.Create(T_pipe), VisibilityType.Show, inSelectedView=False, faceLevel=False)

ViewHelper.SetObjectVisibility(Selection.Create(T_pipe1), VisibilityType.Show, inSelectedView=False, faceLevel=False)

targets = BodySelection.Create(T_pipe, T_pipe1)

result = Combine.Merge(targets)


#Edge fillet right

edge=Selection.Create(result.CreatedEdges)

options = ConstantRoundOptions()

result = ConstantRound.Execute(edge, MM(R_EF1), options)


targets = BodySelection.Create(T_pipe, T_pipe2)

result = Combine.Merge(targets)

T_pipe=Selection.Create(result.CreatedBodies)


#Edge fillet left

edge=Selection.Create(result.CreatedEdges)

options = ConstantRoundOptions()

result = ConstantRound.Execute(edge, MM(R_EF1), options)


#Create branch pipe profile

origin_point=Point.Create(MM(X4_1), MM(0), MM(0))

x_Direction = Direction.DirZ

y_Direction = Direction.DirY

frame=Frame.Create(origin_point,Direction.DirX)    

circle_int = Circle.Create(frame,MM(rint4))

circle_ext = Circle.Create(frame,MM(rout4))

SketchCircle.Create(circle_int)

SketchCircle.Create(circle_ext)

result = ViewHelper.SetViewMode(InteractionMode.Solid)

face=result.CreatedBodies[0].Faces[0]

NamedSelection.Create(Selection.Create(face),Selection.Empty(),"branch")


#extrude branch pipe

selection = Selection.Create(face)

upToSelection = max_face

options = ExtrudeFaceOptions()

options.ExtrudeType = ExtrudeType.ForceIndependent

result = ExtrudeFaces.UpTo(selection, Direction.DirX, upToSelection,Point.Create(MM(X4_2), MM(0), MM(0)), options)

targets = BodySelection.Create(GetRootPart().GetBodies("T_pipe")[0], result.CreatedBodies[0])

result = Combine.Merge(targets)


edge=Selection.Create(result.CreatedEdges[0])

options = ConstantRoundOptions()

result = ConstantRound.Execute(edge, MM(R_EF_4), options)


#cut opening into main pipe

origin_point=Point.Create(MM(X4_1), MM(0), MM(0))

x_Direction = Direction.DirZ

y_Direction = Direction.DirY

frame=Frame.Create(origin_point,Direction.DirX)    

circle_int = Circle.Create(frame,MM(rint4))

SketchCircle.Create(circle_int)

result = ViewHelper.SetViewMode(InteractionMode.Solid)

face=result.CreatedBodies[0].Faces[0]


selection = Selection.Create(face)

options = ExtrudeFaceOptions()

options.ExtrudeType = ExtrudeType.ForceCut

result = ExtrudeFaces.Execute(selection, MM(-X4_1), options)


part=GetRootPart()

all_bodies=Selection.Create(part.GetAllBodies())

ViewHelper.SetObjectVisibility(all_bodies, VisibilityType.Show, inSelectedView=False, faceLevel=False)



plane=DatumPlaneCreator.Create(Point.Create(MM(Z2_1), MM(0), MM(0)),-Direction.DirX)

datum=Selection.Create(plane.CreatedPlanes[0])

# Slice Bodies by Plane

selection = Selection.Create(part.GetBodies("T_pipe")[0])

result = SplitBody.ByCutter(selection,datum)

result.CreatedBodies[0].Name='Branch pipe'

body=Selection.Create(result.CreatedBodies[0])

NamedSelection.Create(body,Selection.Empty(),"Branch _pipe")


#Create remaining Named selections

body=Selection.Create(part.GetBodies("T_pipe")[0])

NamedSelection.Create(Selection.Create(part.GetBodies("T_pipe")[0]),Selection.Empty(),"main_t_section")

NamedSelection.Create(Selection.Create(part.GetAllBodies()),Selection.Empty(),"all_bodies")


#Main pipe id od

selection=Selection.CreateByNames("Main pipe right","Main pipe left")

selection=selection.ConvertToFaces()     

cylindrical_faces = []

for face in selection.Items:

    if isinstance(face.Shape.Geometry,Cylinder) :

        cylindrical_faces.append(face)

NamedSelection.Create(Selection.Create(cylindrical_faces),Selection.Empty(),"main_pipe_id_od")


#branch pipe id od

selection=Selection.CreateByNames("Branch pipe")

selection=selection.ConvertToFaces()     

cylindrical_faces = []

for face in selection.Items:

    if isinstance(face.Shape.Geometry,Cylinder) :

        cylindrical_faces.append(face)

NamedSelection.Create(Selection.Create(cylindrical_faces),Selection.Empty(),"branch_pipe_id_od")


#differential_pressure

selection=Selection.Create(GetRootPart().GetAllBodies())

selection=selection.ConvertToFaces()     

cylindrical_innerfaces = []

for face in selection.Items:

    if isinstance(face.Shape.Geometry,Cylinder) :

        if face.Shape.Geometry.Radius==MM(rint2) or face.Shape.Geometry.Radius==MM(rint4):

            cylindrical_innerfaces.append(face)

NamedSelection.Create(Selection.Create(cylindrical_innerfaces),Selection.Empty(),"differentialpressure")


#jointfaces

selection=Selection.Create(GetRootPart().GetBodies("Main pipe right")[0].Faces[0])

jointfaces = PowerSelection.Faces.ByArea(selection, 

    PowerSelectOptions(True), 

    SearchCriteria.SizeComparison.Equal)


result = NamedSelection.Create(jointfaces, Selection.Empty(), "join_faces")


# Share Topology

options = ShareTopologyOptions()

options.Tolerance = MM(0.5)

result = ShareTopology.Create(options)

# EndBlock



```