def iges_144_000(): import pyiges.IGESGeomLib as IGES from pyiges.IGESGeomLib import IGESPoint from pyiges.IGESCore import IGEStorage system = IGEStorage() examples.benchmarks.standard_iges_setup(system, "144-000-benchmark.igs") para_spline_surface = IGES.IGESTestSplineSurf() system.Commit(para_spline_surface) circle = IGES.IGESGeomCircle(IGESPoint(6, 7.25, 0), IGESPoint(6.25, 7.25)) system.Commit(circle) para_spline_curve = IGES.IGESTestSplineCurve() system.Commit(para_spline_curve) curve_on_parametric_surface = IGES.IGESCurveOnParametricSurface( para_spline_surface, circle, para_spline_curve, 2) system.Commit(curve_on_parametric_surface) trimmed_parametric_surface = IGES.IGESTrimmedParaSurface( para_spline_surface, curve_on_parametric_surface) system.Commit(trimmed_parametric_surface) system.save("144-000-benchmark.igs") if not os.environ.get('READTHEDOCS', None): print(system) os.startfile("144-000-benchmark.igs")
def testrun(filename="IGESFile.igs"): system = IGEStorage() system.StartSection.Prolog = " " system.GlobalSection.IntegerBits = int(32) system.GlobalSection.SPMagnitude = int(38) system.GlobalSection.SPSignificance = int(6) system.GlobalSection.DPMagnitude = int(38) system.GlobalSection.DPSignificance = int(15) system.GlobalSection.MaxNumberLineWeightGrads = int(8) system.GlobalSection.WidthMaxLineWeightUnits = float(0.016) system.GlobalSection.MaxCoordValue = float(71) P = [[0, 0], [numpy.divide(1, 3), numpy.divide(numpy.pi, 6)], [numpy.divide(2, 3), 1], [1, 1], [1 + numpy.divide(1, 3), 1], [1 + numpy.divide(2, 3), numpy.divide(numpy.pi, 6)], [2, 0], [2 + numpy.divide(1, 3), 0 - numpy.divide(numpy.pi, 6)], [2 + numpy.divide(2, 3), 0 - 1], [3, 0 - 1], [3 + numpy.divide(1, 3), 0 - 1], [3 + numpy.divide(2, 3), 0 - numpy.divide(numpy.pi, 6)], [4, 0]] for i in range(0, 1): P.extend P = numpy.transpose(P) bezi = bezier_curve(P, nTimes=50) polyln = IGES.IGESGeomPolyline() for n in range(0, 5): for i in range(0, len(bezi[0])): polyln.AddPoint(IGESPoint(bezi[0][i], bezi[1][i], 10)) bezi[0] = bezi[0] + 4 system.Commit(polyln) line = IGES.IGESGeomLine(IGESPoint(-2, -5, 0), IGESPoint(22, -5, 0)) system.Commit(line) #system.Commit(IGES.IGESRevolve(polyln, line)) #system.Commit(IGES.IGESExtrude(polyln.DirectoryDataPointer.data, IGESPoint(0,0,10))) print(system) system.save(filename) if not environ.get('READTHEDOCS', None) == 'True': startfile(filename)
def testrun_random_surface(filename="IGESFile_random_surface.igs"): """draw a randomiced surface consisting of b spline surface segments""" system = IGEStorage() system.StartSection.Prolog = " " system.GlobalSection.IntegerBits = int(32) system.GlobalSection.SPMagnitude = int(38) system.GlobalSection.SPSignificance = int(6) system.GlobalSection.DPMagnitude = int(38) system.GlobalSection.DPSignificance = int(15) system.GlobalSection.MaxNumberLineWeightGrads = int(8) system.GlobalSection.WidthMaxLineWeightUnits = float(0.016) system.GlobalSection.MaxCoordValue = float(71) scale = 2.5 # data = numpy.zeros(shape=(3,5), dtype=int) data = numpy.random.rand(100, 100) # print("random data generated:", data.shape) for line_nr in range(0, data.shape[0]-1): # print("line number:", line_nr) for row_nr in range(0, data.shape[1]-1): x1 = line_nr * scale x2 = (line_nr + 1) * scale y1 = row_nr * scale y2 = (row_nr + 1) * scale p1 = IGESPoint(x1, y1, float(data[line_nr][row_nr])) p2 = IGESPoint(x2, y1, float(data[line_nr+1][row_nr])) p3 = IGESPoint(x1, y2, float(data[line_nr][row_nr+1])) p4 = IGESPoint(x2, y2, float(data[line_nr+1][row_nr+1])) square = IGES.IGESRationalBSplineSurface(p1, p2, p3, p4) system.Commit(square) # print("save file:", line_nr) system.save(filename)
def testrun_rotate_polyline(filename="IGESFile_rotate_polyline.igs"): """draw a poly line and rotate it around a line to for a cylindrical object""" system = IGEStorage() system.StartSection.Prolog = " " system.GlobalSection.IntegerBits = int(32) system.GlobalSection.SPMagnitude = int(38) system.GlobalSection.SPSignificance = int(6) system.GlobalSection.DPMagnitude = int(38) system.GlobalSection.DPSignificance = int(15) system.GlobalSection.MaxNumberLineWeightGrads = int(8) system.GlobalSection.WidthMaxLineWeightUnits = float(0.016) system.GlobalSection.MaxCoordValue = float(71) # add poly line from numpy array P = [[0, 0], [numpy.divide(1, 3), numpy.divide(numpy.pi, 6)], [numpy.divide(2, 3), 1], [1, 1], [1 + numpy.divide(1, 3), 1], [1 + numpy.divide(2, 3), numpy.divide(numpy.pi, 6)], [2, 0], [2 + numpy.divide(1, 3), 0 - numpy.divide(numpy.pi, 6)], [2 + numpy.divide(2, 3), 0 - 1], [3, 0 - 1], [3 + numpy.divide(1, 3), 0 - 1], [3 + numpy.divide(2, 3), 0 - numpy.divide(numpy.pi, 6)], [4, 0]] for i in range(0, 1): P.extend P = numpy.transpose(P) bezi = bezier_curve(P, nTimes = 50) polyln = IGES.IGESGeomPolyline() for n in range(0, 5): for i in range(0, len(bezi[0])): polyln.AddPoint(IGESPoint(bezi[0][i], bezi[1][i], 10)) bezi[0] = bezi[0] + 4 system.Commit(polyln) # add straight line line = IGES.IGESGeomLine(IGESPoint(-2, -5, 0), IGESPoint(22, -5, 0)) system.Commit(line) # revolve polyline around line system.Commit(IGES.IGESRevolve(polyln, line)) system.save(filename)
def testrun_extrude_polyline(filename="IGESFile_extrude_polyline.igs"): """draw a poly line and extrude it to a target points""" system = IGEStorage() system.StartSection.Prolog = " " system.GlobalSection.IntegerBits = int(32) system.GlobalSection.SPMagnitude = int(38) system.GlobalSection.SPSignificance = int(6) system.GlobalSection.DPMagnitude = int(38) system.GlobalSection.DPSignificance = int(15) system.GlobalSection.MaxNumberLineWeightGrads = int(8) system.GlobalSection.WidthMaxLineWeightUnits = float(0.016) system.GlobalSection.MaxCoordValue = float(71) # add poly line from numpy array P = [[0, 0], [numpy.divide(1, 3), numpy.divide(numpy.pi, 6)], [numpy.divide(2, 3), 1], [1, 1], [1 + numpy.divide(1, 3), 1], [1 + numpy.divide(2, 3), numpy.divide(numpy.pi, 6)], [2, 0], [2 + numpy.divide(1, 3), 0 - numpy.divide(numpy.pi, 6)], [2 + numpy.divide(2, 3), 0 - 1], [3, 0 - 1], [3 + numpy.divide(1, 3), 0 - 1], [3 + numpy.divide(2, 3), 0 - numpy.divide(numpy.pi, 6)], [4, 0]] for i in range(0, 1): P.extend P = numpy.transpose(P) bezi = bezier_curve(P, nTimes = 50) polyln = IGES.IGESGeomPolyline() for n in range(0, 5): for i in range(0, len(bezi[0])): polyln.AddPoint(IGESPoint(bezi[0][i], bezi[1][i], 10)) bezi[0] = bezi[0] + 4 system.Commit(polyln) # extrude polyline to endpoint system.Commit(IGES.IGESExtrude(polyln, IGESPoint(0, 0, 15))) system.save(filename)
def testrun_arc_circle(filename="IGESFile_arc_circle.igs"): """draw an arc and a circle, both in a plane paralell to the X-Y-plane""" system = IGEStorage() system.StartSection.Prolog = " " system.GlobalSection.IntegerBits = int(32) system.GlobalSection.SPMagnitude = int(38) system.GlobalSection.SPSignificance = int(6) system.GlobalSection.DPMagnitude = int(38) system.GlobalSection.DPSignificance = int(15) system.GlobalSection.MaxNumberLineWeightGrads = int(8) system.GlobalSection.WidthMaxLineWeightUnits = float(0.016) system.GlobalSection.MaxCoordValue = float(71) arc = IGES.IGESGeomArc(-3, IGESPoint(2, 5, 0), IGESPoint(-3, 0, 0), IGESPoint(2, 10, 0)) system.Commit(arc) circ = IGES.IGESGeomCircle(IGESPoint(5, 5, 5), 5) system.Commit(circ) system.save(filename)
def testrun_extrude_arc(filename="IGESFile_extrude_arc.igs"): """draw arc and extrude it to a target points""" system = IGEStorage() system.StartSection.Prolog = " " system.GlobalSection.IntegerBits = int(32) system.GlobalSection.SPMagnitude = int(38) system.GlobalSection.SPSignificance = int(6) system.GlobalSection.DPMagnitude = int(38) system.GlobalSection.DPSignificance = int(15) system.GlobalSection.MaxNumberLineWeightGrads = int(8) system.GlobalSection.WidthMaxLineWeightUnits = float(0.016) system.GlobalSection.MaxCoordValue = float(71) arc = IGES.IGESGeomArc(0, IGESPoint(20, 11, 0), IGESPoint(20.25, 11, 0), IGESPoint(19.75, 11, 0)) system.Commit(arc) ext = IGES.IGESExtrude(arc, IGESPoint(20.4857, 11.2357, -0.9428)) system.Commit(ext) system.save(filename)
def iges_144_000(): import pyiges.IGESGeomLib as IGES from pyiges.IGESCore import IGEStorage from pyiges.IGESGeomLib import IGESPoint filename = "144-000-example-1.igs" system = IGEStorage() examples.benchmarks.standard_iges_setup(system, filename) # Circular flat surface profile circular_flat_surface_profile = IGES.IGESGeomCircle( IGESPoint(0, 0, 0), IGESPoint(10, 10, 0)) system.Commit(circular_flat_surface_profile) # Make the surface from out outer profile, this is to be trimmed circular_flat_surface = IGES.IGESGeomPlane(circular_flat_surface_profile) system.Commit(circular_flat_surface) # We're going to put a bunch of holes into the surface hole = [] # Hole profiles hole_copse = [] # Curve on Parametric Surface entries (required!) for i in range(0, 6): hole.append( IGES.IGESGeomCircle(IGESPoint(-5 + i * 3, 0, 0), IGESPoint(-5 + i * 3, 1, 0))) system.Commit(hole[-1]) # Make those tubes, given we are not using them as a reference, # Lets directly add them to the system. system.Commit(IGES.IGESExtrude(hole[-1], IGESPoint(0, 0, i + 1))) system.Commit(IGES.IGESExtrude(hole[-1], IGESPoint(0, 0, -i - 1))) hole_copse.append( IGES.IGESCurveOnParametricSurface(circular_flat_surface, hole[-1], hole[-1], 1)) system.Commit(hole_copse[-1]) # Setup the surface we are to put holes into, we're using the whole surface, # as indicated by the 0 in the outer profile parameter trimmed_surface = IGES.IGESTrimmedParaSurface(circular_flat_surface, 0) trimmed_surface.N1 = 0 # Lets put those holes in! for _hole in hole_copse: trimmed_surface.add_bounding_profile(_hole) system.Commit(trimmed_surface) system.save(filename) if not os.environ.get('READTHEDOCS', None): print(system) os.startfile(filename)
def iges_414_000(): import pyiges.IGESGeomLib as IGES from pyiges.IGESCore import IGEStorage from pyiges.IGESGeomLib import IGESPoint system = IGEStorage() examples.benchmarks.standard_iges_setup(system, "414-000-benchmark.igs") line_1 = IGES.IGESGeomLine(IGESPoint(-0.13, 0, 0), IGESPoint(0, 0.25, 0)) line_2 = IGES.IGESGeomLine(IGESPoint(0, 0.25, 0), IGESPoint(0.13, 0, 0)) line_3 = IGES.IGESGeomLine(IGESPoint(0, 0.25, 0), IGESPoint(0, 0.44, 0)) line_4 = IGES.IGESGeomLine(IGESPoint(0, 0.38, 0), IGESPoint(0.09, 0.28, 0)) line_5 = IGES.IGESGeomLine(IGESPoint(0, 0.38, 0), IGESPoint(-0.09, 0.28, 0)) circle_1 = IGES.IGESGeomArc(0, IGESPoint(0, 0.5), IGESPoint(0.06, 0.5), IGESPoint(0.06, 0.5)) system.Commit(line_1) system.Commit(line_2) system.Commit(line_3) system.Commit(line_4) system.Commit(line_5) system.Commit(circle_1) person_group = IGES.IGESGroup("PERSON", line_1, line_2, line_3, line_4, line_5, circle_1) system.Commit(person_group) person_array = IGES.IGESCircularArray(person_group, 3, IGESPoint(5.5, 1.5, 0), 0.5, 0.52, 2.09) system.Commit(person_array) system.save("414-000-benchmark.igs") if not os.environ.get('READTHEDOCS', None): print(system) os.startfile("414-000-benchmark.igs")
def testrun_torus(filename="IGESFile_torus.igs"): """Dont know if this works, neither FreeCAD not IGES Viewer will display anyting for a torus element""" system = IGEStorage() system.StartSection.Prolog = " " system.GlobalSection.IntegerBits = int(32) system.GlobalSection.SPMagnitude = int(38) system.GlobalSection.SPSignificance = int(6) system.GlobalSection.DPMagnitude = int(38) system.GlobalSection.DPSignificance = int(15) system.GlobalSection.MaxNumberLineWeightGrads = int(8) system.GlobalSection.WidthMaxLineWeightUnits = float(0.016) system.GlobalSection.MaxCoordValue = float(71) torus = IGES.IGESGeomTorus(50, 10, IGESPoint(0,0,0), IGESVector(0,0,1)) system.Commit(torus) system.save(filename)
def testrun_spline_curve(filename="IGESFile_spline_curve.igs"): """draw a spline curve""" system = IGEStorage() system.StartSection.Prolog = " " system.GlobalSection.IntegerBits = int(32) system.GlobalSection.SPMagnitude = int(38) system.GlobalSection.SPSignificance = int(6) system.GlobalSection.DPMagnitude = int(38) system.GlobalSection.DPSignificance = int(15) system.GlobalSection.MaxNumberLineWeightGrads = int(8) system.GlobalSection.WidthMaxLineWeightUnits = float(0.016) system.GlobalSection.MaxCoordValue = float(71) spline_curve = IGES.IGESSplineCurve(stype=3, h=1, ndim=3) spline_curve.addSegment(0., [5.817, 0.183686, 7.152560000000000E-007, 0.00391344, 7.42046, 0.116096, -1.430510000000000E-006, -0.0365138, 0.582714, 0.26764, -2.682210000000000E-007, -0.0384366]) spline_curve.addSegment(1., [6.0046, 0.195428, 0.011741, -0.0231148, 7.50004, 0.00655174, -0.109543, 0.0169944, 0.811917, 0.15233, -0.11531, -0.00673003]) spline_curve.addSegment(2., [6.18865, 0.149566, -0.0576034, -0.034095, 7.41404, -0.161551, -0.0585597, 0.0144959, 0.842207, - 0.0984805, -0.1355, 0.000767251]) spline_curve.addSegment(3., [6.24652, -0.0679259, -0.159888, 0.0609804, 7.20843, -0.235182, -0.0150719, 0.0548643, 0.608993, - 0.367179, -0.133198, 0.11326]) spline_curve.addSegment(4., [6.07969, -0.204761, 0.0230529, -0.00538516, 7.01304, -0.100733, 0.149521, -0.0375805, 0.221876, - 0.293796, 0.206582, -0.0375974]) spline_curve.addSegment(5., [5.89259, -0.174811, 0.00689745, 0.041456, 7.02425, 0.0855677, 0.0367796, 0.0150364, 0.0970648, 0.00657602, 0.0937898, 0.0214335]) spline_curve.addSegment(6., [5.76614, -0.0366479, 0.131266, -0.0437549, 7.16163, 0.204236, 0.0818889, -0.0272963, 0.218864, 0.258456, 0.15809, -0.0526968]) spline_curve.addSegment(7., [5.817, 0.0946183, 1.430510000000000E-006, -0.26253, 7.42046, 0.286125, 0., -0.163778, 0.582714, 0.146546, -1.713630000000000E-007, -0.316181]) system.Commit(spline_curve) system.save(filename)
def iges_144_000(): import pyiges.IGESGeomLib as IGES from pyiges.IGESCore import IGEStorage from pyiges.IGESGeomLib import IGESPoint filename = "144-000-example-2.igs" system = IGEStorage() examples.benchmarks.standard_iges_setup(system, filename) # Setup a cone like surface, the cone is boring but whatever! cone_surface_profile = IGES.IGESGeomPolyline(IGESPoint(-10, 0, 0), IGESPoint(0, 0, 10)) system.Commit(cone_surface_profile) # The centre of the revolve center_line = IGES.IGESGeomLine(IGESPoint(0, 0, 0), IGESPoint(0, 0, 10)) system.Commit(center_line) # Make the surface from the profile, this is to be trimmed cone_surface = IGES.IGESRevolve(cone_surface_profile, center_line) system.Commit(cone_surface) #Setup the surface we are to put holes into, we're using the whole surface, # as indicated by the 0 in the outer profile parameter trimmed_surface = IGES.IGESTrimmedParaSurface(cone_surface, 0) trimmed_surface.N1 = 0 hole = [] hole_copse = [] # We're going to put a bunch of holes into the surface, because we don't # care about the index of the hole we can get away with specifying -8 to 8 for holenum in range(-8, 8): hole.append(IGES.IGESGeomPolyline()) hole_data = [[], [], []] # Top line following the cone profile for i in range(21 * holenum, 20 + 21 * holenum): hole_data[0].append(3) hole_data[1].append(i) hole_data[2].append(7) # Line down the cone profile hole_data[0].append(7) hole_data[1].append(hole_data[1][-1]) hole_data[2].append(3) # Bottom line following the cone profile for i in range(20 + 21 * holenum, 21 * holenum, -1): hole_data[0].append(7) hole_data[1].append(i) hole_data[2].append(3) # Line back to start of the cone profile hole_data[0].append(hole_data[0][0]) hole_data[1].append(hole_data[1][0]) hole_data[2].append(hole_data[2][0]) # Convert to rectangular coordinate system hole_rect_cs = cyl_to_rect(hole_data) # Push the points into the common format for i in range(0, len(hole_rect_cs[0])): hole[-1].AddPoint( IGESPoint(hole_rect_cs[0][i], hole_rect_cs[1][i], hole_rect_cs[2][i])) # Commit this hole system.Commit(hole[-1]) # Lets make a really very crazy looking thing! # And revolve the cutting profile by a bit # Setup Revolve Line data random_line_data = cyl_to_rect([[-10, 10], [(20 * holenum) + 90, (20 * holenum) + 90], [10, 10]]) # Make revolve line revolve_line = IGES.IGESGeomLine( IGESPoint(random_line_data[0][0], random_line_data[1][0], random_line_data[2][0]), IGESPoint(random_line_data[0][1], random_line_data[1][1], random_line_data[2][1])) # Commit revolve line system.Commit(revolve_line) # Make the revolve system.Commit(IGES.IGESRevolve(hole[-1], revolve_line, -1, 0)) # Create the Curve on Parametric Surface hole_copse.append( IGES.IGESCurveOnParametricSurface(cone_surface, hole[-1], hole[-1], 2)) system.Commit(hole_copse[-1]) # Lets put that holes in! trimmed_surface.add_bounding_profile(hole_copse[-1]) # commit that surface! system.Commit(trimmed_surface) system.save(filename) if not os.environ.get('READTHEDOCS', None): print(system) os.startfile(filename)
system.GlobalSection.SPSignificance = 15 system.GlobalSection.DPMagnitude = 13 system.GlobalSection.DPSignificance = 15 system.GlobalSection.ProductIdentificationForReceiver = "122-000" system.GlobalSection.ModelSpaceScale = 1.0 system.GlobalSection.Units.setInches() system.GlobalSection.MaxNumberLineWeightGrads = 8 system.GlobalSection.WidthMaxLineWeightUnits = 0.016 system.GlobalSection.MaxUserResolution = 0.0001 system.GlobalSection.MaxCoordValue = 20.4857 system.GlobalSection.NameOfAuthor = "github.com" system.GlobalSection.AuthorOrg = "pyIGES" system.GlobalSection.VersionFlag = 11 system.GlobalSection.DraftStandardFlag = 3 system.GlobalSection.AppProtocol = "pyIGES" arc = IGES.IGESGeomArc(0, IGES.IGESPoint(20, 11), IGES.IGESPoint(20.25, 11), IGES.IGESPoint(19.75, 11)) arc.Color.setYellow() arc.StatusNumber.Hierachy.setGlobalDefer() arc.StatusNumber.EntityUseFlag.setGeometry() arc.StatusNumber.Subordinate.setPhysicallyDependent() arc.StatusNumber.Visablilty.setVisible() system.Commit(arc) extr = IGES.IGESExtrude(arc, IGES.IGESPoint(20.4857, 11.2357, -0.9428)) arc.Color.setYellow() system.Commit(extr) system.save(filename=filename)
numpy.divide(numpy.pi, 6)], [numpy.divide(2, 3), 1], [1, 1], [1 + numpy.divide(1, 3), 1], [1 + numpy.divide(2, 3), numpy.divide(numpy.pi, 6)], [2, 0], [2 + numpy.divide(1, 3), 0 - numpy.divide(numpy.pi, 6)], [2 + numpy.divide(2, 3), 0 - 1], [3, 0 - 1], [3 + numpy.divide(1, 3), 0 - 1], [3 + numpy.divide(2, 3), 0 - numpy.divide(numpy.pi, 6)], [4, 0]] for i in range(0, 1): P.extend P = numpy.transpose(P) bezi = bezier_curve(P, nTimes=50) polyln = IGES.IGESGeomPolyline() system.Commit(polyln) line = IGES.IGESGeomLine(IGESPoint(-2, -5, 0), IGESPoint(22, -5, 0)) print(line) system.Commit(line) line = IGES.IGESGeomLine(IGESPoint(-2, -1, 0), IGESPoint(22, -3, 3)) print(line) system.Commit(line) system.Commit(IGES.IGESRevolve(polyln, line)) # system.Commit(IGES.IGESExtrude(polyln.DirectoryDataPointer.data, IGESPoint(0,0,10))) # system.Commit(IGES.IGESRevolve(line0, line))
system.GlobalSection.DPMagnitude = 13 system.GlobalSection.DPSignificance = 15 system.GlobalSection.ProductIdentificationForReceiver = filename.split(".")[0] system.GlobalSection.ModelSpaceScale = 3.0 system.GlobalSection.Units.setInches() system.GlobalSection.MaxNumberLineWeightGrads = 8 system.GlobalSection.WidthMaxLineWeightUnits = 0.016 system.GlobalSection.MaxUserResolution = 0.0001 system.GlobalSection.MaxCoordValue = 15.5 system.GlobalSection.NameOfAuthor = "github.com" system.GlobalSection.AuthorOrg = "pyIGES" system.GlobalSection.VersionFlag = 11 system.GlobalSection.DraftStandardFlag = 3 system.GlobalSection.AppProtocol = "pyIGES" prop = IGES.IGESPropertyEntity([2, 6., 8.], formNumber=16) prop.StatusNumber.Subordinate.setLogicallyDependent() system.Commit(prop) view = IGES.IGESViewEntity([1, 3., 0, 0, 0, 0, 0, 0]) view.StatusNumber.EntityUseFlag.setOther() view.StatusNumber.Subordinate.setLogicallyDependent() system.Commit(view) drawing = IGES.IGESDrawingEntity([1, 3, 0., 0., 0, 0, 1, 1], formNumber=1) drawing.StatusNumber.EntityUseFlag.setOther() system.Commit(drawing) trans = IGES.IGESGeomTransform( [-1., 0., 0., 3.5, 0., 1., 0., 15., 0., 0., -1., 0.], formNumber=1) trans.StatusNumber.Hierachy.setGlobalDefer()