def CreateMultiPatch():
## create arc 1
r1 = 1.0
arc1_ptr = geometry_factory.CreateSmallArc(
[0.0, 0.0, 0.0], 'z', r1, 0.0,
90.0) # put 'y' will give interesting ring
arc1 = arc1_ptr.GetReference()
arc1.Id = 1
## create arc 2
r2 = 2.0
arc2_ptr = geometry_factory.CreateSmallArc([0.0, 0.0, 0.0], 'z', r2, 0.0,
90.0)
arc2 = arc2_ptr.GetReference()
arc2.Id = 2
## create ring patch by connect the two arcs
ring_patch_ptr = bsplines_patch_util.CreateLoftPatch(arc2, arc1)
ring_patch = ring_patch_ptr.GetReference()
ring_patch.Id = 1
######create multipatch
mpatch = MultiPatch2D()
mpatch.AddPatch(ring_patch_ptr)
return mpatch
def CreateMultiPatch():
####### create arc 1
r1 = 1.0
arc1_ptr = geometry_factory.CreateSmallArc([0.0, 0.0, 0.0], 'z', r1, 0.0,
45.0)
arc1 = arc1_ptr.GetReference()
# create line 1
b1 = 4.0
line1_ptr = geometry_factory.CreateLine([b1, 0.0, 0.0], [b1, b1, 0.0],
arc1.Order(0))
line1 = line1_ptr.GetReference()
# create patch 1
patch1_ptr = bsplines_patch_util.CreateConnectedPatch(line1, arc1)
patch1 = patch1_ptr.GetReference()
patch1.Id = 1
####### create arc 2
arc2_ptr = geometry_factory.CreateSmallArc([0.0, 0.0, 0.0], 'z', r1, 45.0,
90.0)
arc2 = arc2_ptr.GetReference()
# create line 2
line2_ptr = geometry_factory.CreateLine([b1, b1, 0.0], [0.0, b1, 0.0],
arc2.Order(0))
line2 = line2_ptr.GetReference()
# create patch 2
patch2_ptr = bsplines_patch_util.CreateConnectedPatch(line2, arc2)
patch2 = patch2_ptr.GetReference()
patch2.Id = 2
######create multipatch
mpatch = MultiPatch2D()
mpatch.AddPatch(patch1_ptr)
mpatch.AddPatch(patch2_ptr)
bsplines_patch_util.MakeInterface(patch1, BoundarySide.Right, patch2,
BoundarySide.Left,
BoundaryDirection.Forward)
multipatch_refine_util.DegreeElevate(patch1_ptr, [1, 1])
multipatch_refine_util.InsertKnots(patch1_ptr, [[0.5], [0.5]])
return mpatch
def CreateMultiPatch():
#############################################
############### patch 1 #####################
####### create arc 1
R = 10.5
L = 50
scale_value = 0.98
arc1_ptr = geometry_factory.CreateSmallArc([0.0, 0.0, 0.0], 'y', R, 180.0,
90.0)
arc1 = arc1_ptr.GetReference()
####### create arc 2
arc2_ptr = geometry_factory.CreateSmallArc([0.0, L / 2.0, 0.0], 'y',
R * scale_value, 180.0, 90.0)
arc2 = arc2_ptr.GetReference()
# create patch 1
patch1_ptr = bsplines_patch_util.CreateLoftPatch(arc1, arc2)
patch1 = patch1_ptr.GetReference()
patch1.Id = 1
############################################
############### patch 2 ####################
####### create arc 2_1
arc3_ptr = geometry_factory.CreateSmallArc([0.0, L / 2.0, 0.0], 'y',
R * scale_value, 180.0, 90.0)
arc3 = arc3_ptr.GetReference()
####### create arc 2_2
arc4_ptr = geometry_factory.CreateSmallArc([0.0, L, 0.0], 'y', R, 180.0,
90.0)
arc4 = arc4_ptr.GetReference()
# create patch 2
patch2_ptr = bsplines_patch_util.CreateLoftPatch(arc3, arc4)
patch2 = patch2_ptr.GetReference()
patch2.Id = 2
############################################
############### patch 3 ####################
####### create arc 3_1
arc5_ptr = geometry_factory.CreateSmallArc([0.0, 0.0, 0.0], 'y', R, 90.0,
0.0)
arc5 = arc5_ptr.GetReference()
####### create arc 3_2
arc6_ptr = geometry_factory.CreateSmallArc([0.0, L / 2.0, 0.0], 'y',
R * scale_value, 90.0, 0.0)
arc6 = arc6_ptr.GetReference()
# create patch 2
patch3_ptr = bsplines_patch_util.CreateLoftPatch(arc5, arc6)
patch3 = patch3_ptr.GetReference()
patch3.Id = 3
############################################
############### patch 4 ####################
####### create arc 4_1
arc7_ptr = geometry_factory.CreateSmallArc([0.0, L / 2.0, 0.0], 'y',
R * scale_value, 90.0, 0.0)
arc7 = arc7_ptr.GetReference()
####### create arc 4_2
arc8_ptr = geometry_factory.CreateSmallArc([0.0, L, 0.0], 'y', R, 90.0,
0.0)
arc8 = arc8_ptr.GetReference()
# create patch 2
patch4_ptr = bsplines_patch_util.CreateLoftPatch(arc7, arc8)
patch4 = patch4_ptr.GetReference()
patch4.Id = 4
#############################################
############### patch 5 #####################
####### create arc 1
arc9_ptr = geometry_factory.CreateSmallArc([0.0, 0.0, 0.0], 'y', R, 0.0,
-90.0)
arc9 = arc9_ptr.GetReference()
####### create arc 2
arc10_ptr = geometry_factory.CreateSmallArc([0.0, L / 2.0, 0.0], 'y',
R * scale_value, 0.0, -90.0)
arc10 = arc10_ptr.GetReference()
# create patch 1
patch5_ptr = bsplines_patch_util.CreateLoftPatch(arc9, arc10)
patch5 = patch5_ptr.GetReference()
patch5.Id = 5
############################################
############### patch 6 ####################
####### create arc 2_1
arc11_ptr = geometry_factory.CreateSmallArc([0.0, L / 2.0, 0.0], 'y',
R * scale_value, 0.0, -90.0)
arc11 = arc11_ptr.GetReference()
####### create arc 2_2
arc12_ptr = geometry_factory.CreateSmallArc([0.0, L, 0.0], 'y', R, 0.0,
-90.0)
arc12 = arc12_ptr.GetReference()
# create patch 2
patch6_ptr = bsplines_patch_util.CreateLoftPatch(arc11, arc12)
patch6 = patch6_ptr.GetReference()
patch6.Id = 6
############################################
############### patch 7 ####################
####### create arc 3_1
arc13_ptr = geometry_factory.CreateSmallArc([0.0, 0.0, 0.0], 'y', R, -90.0,
-180.0)
arc13 = arc13_ptr.GetReference()
####### create arc 3_2
arc14_ptr = geometry_factory.CreateSmallArc([0.0, L / 2.0, 0.0], 'y',
R * scale_value, -90.0, -180.0)
arc14 = arc14_ptr.GetReference()
# create patch 2
patch7_ptr = bsplines_patch_util.CreateLoftPatch(arc13, arc14)
patch7 = patch7_ptr.GetReference()
patch7.Id = 7
############################################
############### patch 8 ####################
####### create arc 4_1
arc15_ptr = geometry_factory.CreateSmallArc([0.0, L / 2.0, 0.0], 'y',
R * scale_value, -90.0, -180.0)
arc15 = arc15_ptr.GetReference()
####### create arc 4_2
arc16_ptr = geometry_factory.CreateSmallArc([0.0, L, 0.0], 'y', R, -90.0,
-180.0)
arc16 = arc16_ptr.GetReference()
# create patch 2
patch8_ptr = bsplines_patch_util.CreateLoftPatch(arc15, arc16)
patch8 = patch8_ptr.GetReference()
patch8.Id = 8
######create multipatch
mpatch = MultiPatch2D()
mpatch.AddPatch(patch1_ptr)
mpatch.AddPatch(patch2_ptr)
mpatch.AddPatch(patch3_ptr)
mpatch.AddPatch(patch4_ptr)
mpatch.AddPatch(patch5_ptr)
mpatch.AddPatch(patch6_ptr)
mpatch.AddPatch(patch7_ptr)
mpatch.AddPatch(patch8_ptr)
bsplines_patch_util.MakeInterface(patch1, BoundarySide.Top, patch2,
BoundarySide.Bottom,
BoundaryDirection.Forward)
bsplines_patch_util.MakeInterface(patch3, BoundarySide.Top, patch4,
BoundarySide.Bottom,
BoundaryDirection.Forward)
bsplines_patch_util.MakeInterface(patch5, BoundarySide.Top, patch6,
BoundarySide.Bottom,
BoundaryDirection.Forward)
bsplines_patch_util.MakeInterface(patch7, BoundarySide.Top, patch8,
BoundarySide.Bottom,
BoundaryDirection.Forward)
bsplines_patch_util.MakeInterface(patch7, BoundarySide.Right, patch1,
BoundarySide.Left,
BoundaryDirection.Forward)
bsplines_patch_util.MakeInterface(patch8, BoundarySide.Right, patch2,
BoundarySide.Left,
BoundaryDirection.Forward)
bsplines_patch_util.MakeInterface(patch1, BoundarySide.Right, patch3,
BoundarySide.Left,
BoundaryDirection.Forward)
bsplines_patch_util.MakeInterface(patch2, BoundarySide.Right, patch4,
BoundarySide.Left,
BoundaryDirection.Forward)
bsplines_patch_util.MakeInterface(patch3, BoundarySide.Right, patch5,
BoundarySide.Left,
BoundaryDirection.Forward)
bsplines_patch_util.MakeInterface(patch4, BoundarySide.Right, patch6,
BoundarySide.Left,
BoundaryDirection.Forward)
bsplines_patch_util.MakeInterface(patch5, BoundarySide.Right, patch7,
BoundarySide.Left,
BoundaryDirection.Forward)
bsplines_patch_util.MakeInterface(patch6, BoundarySide.Right, patch8,
BoundarySide.Left,
BoundaryDirection.Forward)
return mpatch
#importing Kratos modules
from KratosMultiphysics import *
from KratosMultiphysics.IsogeometricApplication import *
kernel = Kernel() #defining kernel
nurbs_fespace_library = BSplinesFESpaceLibrary()
grid_lib = ControlGridLibrary()
multipatch_util = MultiPatchUtility()
bsplines_patch_util = BSplinesPatchUtility()
mpatch_export = MultiNURBSPatchMatlabExporter()
import geometry_factory
## create arc 1
r1 = 1.0
arc1_ptr = geometry_factory.CreateSmallArc(
[0.0, 0.0, 0.0], 'z', r1, 0.0, 90.0) # put 'y' will give interesting ring
arc1 = arc1_ptr.GetReference()
arc1.Id = 1
## create arc 2
r2 = 2.0
arc2_ptr = geometry_factory.CreateSmallArc([0.0, 0.0, 0.0], 'z', r2, 0.0, 90.0)
arc2 = arc2_ptr.GetReference()
arc2.Id = 2
## create ring patch by connect the two arcs
ring_patch_ptr = bsplines_patch_util.CreateConnectedPatch(arc1, arc2)
ring_patch = ring_patch_ptr.GetReference()
ring_patch.Id = 3
print(ring_patch)
mpatch_export.Export(ring_patch, "ring.m")
from KratosMultiphysics import * from KratosMultiphysics.IsogeometricApplication import * kernel = Kernel() #defining kernel nurbs_fespace_library = BSplinesFESpaceLibrary() grid_lib = ControlGridLibrary() multipatch_util = MultiPatchUtility() multipatch_refine_util = MultiPatchRefinementUtility() bsplines_patch_util = BSplinesPatchUtility() mpatch_export = MultiNURBSPatchMatlabExporter() import geometry_factory ####### create arc 1 r1 = 1.0 arc1_ptr = geometry_factory.CreateSmallArc([0.0, 0.0, 0.0], 'z', r1, 0.0, 90.0) arc1 = arc1_ptr.GetReference() ####### create arc 2 r2 = 2.0 arc2_ptr = geometry_factory.CreateSmallArc([0.0, 0.0, 0.0], 'z', r2, 0.0, 90.0) arc2 = arc2_ptr.GetReference() # create patch 2 patch_ptr = bsplines_patch_util.CreateConnectedPatch(arc1, arc2) patch = patch_ptr.GetReference() patch.Id = 1 ######create multipatch mpatch = MultiPatch2D() mpatch.AddPatch(patch)
##### all rights reserved #####
##################################################################
##################################################################
##################################################################
##################################################################
import sys
import os
kratos_root_path=os.environ['KRATOS_ROOT_PATH']
##################################################################
##################################################################
#importing Kratos modules
from KratosMultiphysics import *
from KratosMultiphysics.BRepApplication import *
from KratosMultiphysics.IsogeometricApplication import *
kernel = Kernel() #defining kernel
import geometry_factory
mpatch_util = MultiPatchUtility()
mpatch_export = MultiNURBSPatchMatlabExporter()
arc_ptr = geometry_factory.CreateSmallArc([0.0, 0.0, 0.0], 'y', 1.0, 0.0, 60.0)
arc = arc_ptr.GetReference()
print(arc)
mpatch_export.Export(arc, "arc.m")
## use this command to check in matlab
# for x: arc;nrbplot(nurbs,10);xlabel('x');ylabel('y');zlabel('z');view([1 0 0]);
# for y: arc;nrbplot(nurbs,10);xlabel('x');ylabel('y');zlabel('z');view([0 1 0]);
# for z: arc;nrbplot(nurbs,10);xlabel('x');ylabel('y');zlabel('z');view([0 0 1]);
