Three defects, all of them the same mistake in three places: the cut face was built from something other than the geometry it claims to show. 1. A prismatic cut face removed a feature by its BOUNDING BOX. `Hole.removal_box` is a square, so the gap punched in the cut face was the full diameter wherever the plane passed. On the plan's 80x50x8 plate with a Ho10 hole at (40, 25): a plane at y=29.9 grazes the hole (true chord 1.990) and the section still showed an 8 mm gap through solid metal - 560.0 mm2 hatched against a true 624.08, a 64 mm2 error in the hatched area AND a wrong drawn outline. `Slot` had the same obround-as-box. Features now answer `cut_spans(part, p1, unit)` with their own exact chord: `_circle_chord` for a hole, the union of two circles and a rotated rectangle for a slot, the polygon's own crossings for a pocket, and an exact slab clip of `removal_box` as the default. `_box_span` (the corner projection, which is the box's shadow and not a chord at all) is gone. 2. A transverse revolved section ignored every feature and the honesty rule never fired. A cross section through a DIN 6885-1 keyway came out as a plain hatched disc with `omitted == ()` - the commonest mechanical cross section there is, hatched straight over the 12x5 keyway. The face is now built by `_transverse_face`: its radii come from `radial_profile` (so a groove, an undercut, a bore or an ISO 6410 internal thread is in the face already, which the raw segment diameters could not report), a non-axisymmetric feature is asked for its own `cross_section` notch record, and one the plane really passes through with no such record is listed in `omitted` with its reason rather than hatched over. `Keyway._notch` is the one calculation the end view and the cut face both come out of. 3. A section drew its own cut boundary as a HIDDEN (dashed) line, because the longitudinal branch reused `silhouette_prims` and the silhouette draws the bore dashed - correct for an un-cut view, wrong for a section, where that line IS the hatch loop's boundary. ISO 128-50 bounds a surface lying in the cutting plane with a continuous wide line and ISO 128-3 leaves hidden detail out of a sectional view; via ROLE_LAYER['hidden'] this would also have put the cut boundary on the HIDDEN layer in Task 7. New `section_profile_prims` emits the cut outline all-continuous, and `half=True` keeps the dashed bore on the -Y half only, which is the half a half section does not cut. And the claim the last commit made without earning it: the hatched area and the drawn outline are now literally one calculation. The longitudinal loop is the *emitted* outline flattened (`_flatten_path`), so a chamfer or fillet that substitutes a silhouette corner is cut out of the hatch exactly where it is drawn - a 2 mm 45 deg chamfer takes its 2 mm2 off the 1000 mm2 half face instead of being hatched over. Eleven new tests, every one of them red on the parent commit: the hole chord at 0/4/4.9 mm off centre, a tangent plane leaving the face whole, the 45 deg slot, the keyway cross section against its own analytic area, the feature a cross section cannot project, the feature the plane misses, the groove read off the profile, the continuous cut boundary, the half section's dashed half, and the chamfered cut face. No new ActiveX member is called: all three fixes are pure geometry in the backend-free `engineering/mech` package. Co-Authored-By: Claude Opus 5 (1M context) <noreply@anthropic.com> Claude-Session: https://claude.ai/code/session_01FKFStjEFzpGyLq5bMpWxTd
474 lines
21 KiB
Python
474 lines
21 KiB
Python
"""Sections, cut faces, material hatch and detail views.
|
|
|
|
The cut face and the silhouette come out of the same radial profile, so the
|
|
tests assert that the hatched loop closes on the real geometry - including the
|
|
ISO 6410 rule that an internal thread's hatch runs to the major diameter, and
|
|
the ISO 128-3 rule that a rib cut longitudinally is not hatched.
|
|
"""
|
|
|
|
from __future__ import annotations
|
|
|
|
import math
|
|
|
|
import pytest
|
|
|
|
from engineering.mech.features import (
|
|
AxialBore,
|
|
Chamfer,
|
|
CrossHole,
|
|
Hole,
|
|
Keyway,
|
|
Pocket,
|
|
RetainingGroove,
|
|
Slot,
|
|
Thread,
|
|
)
|
|
from engineering.mech.part import PrismaticPart, RevolvedPart, Segment
|
|
from engineering.mech.primitives import Arc, Circle, HatchArea, Line, Text
|
|
from engineering.mech.views import (
|
|
CIRCLE_SEGMENTS,
|
|
SECTION_PLANE,
|
|
SECTION_STYLES,
|
|
VIEW_KINDS,
|
|
build_view,
|
|
cut_loops,
|
|
detail_marker_prims,
|
|
normalise_plane,
|
|
radial_profile,
|
|
)
|
|
|
|
EPS = 1e-9
|
|
|
|
|
|
def shaft(*features) -> RevolvedPart:
|
|
return RevolvedPart(
|
|
name="shaft",
|
|
segments=(Segment(20.0, 25.0), Segment(40.0, 30.0), Segment(15.0, 20.0)),
|
|
features=tuple(features),
|
|
)
|
|
|
|
|
|
def plate(*features) -> PrismaticPart:
|
|
return PrismaticPart(
|
|
name="plate",
|
|
outline=((0.0, 0.0), (80.0, 0.0), (80.0, 50.0), (0.0, 50.0)),
|
|
thickness=8.0,
|
|
features=tuple(features),
|
|
)
|
|
|
|
|
|
def axis_plane(**kw) -> dict:
|
|
return normalise_plane({"p1": (0.0, 0.0), "p2": (100.0, 0.0), **kw})
|
|
|
|
|
|
def polygon_area(points) -> float:
|
|
total = 0.0
|
|
count = len(points)
|
|
for i in range(count):
|
|
x1, y1 = points[i]
|
|
x2, y2 = points[(i + 1) % count]
|
|
total += x1 * y2 - x2 * y1
|
|
return abs(total) / 2.0
|
|
|
|
|
|
def test_the_engine_now_declares_five_kinds_and_four_section_styles():
|
|
assert VIEW_KINDS == ("front", "side", "top", "section", "detail")
|
|
assert SECTION_STYLES == ("full", "half", "offset", "revolved")
|
|
assert set(SECTION_PLANE) == {"p1", "p2", "label", "direction", "style"}
|
|
|
|
|
|
def test_normalise_plane_fills_the_pinned_shape_and_refuses_a_zero_length_plane():
|
|
plane = normalise_plane({"p1": (0.0, 0.0), "p2": (10.0, 0.0)})
|
|
assert plane["label"] == "A"
|
|
assert plane["direction"] == (0.0, -1.0)
|
|
assert plane["style"] == "full"
|
|
with pytest.raises(ValueError, match="zero length"):
|
|
normalise_plane({"p1": (1.0, 1.0), "p2": (1.0, 1.0)})
|
|
|
|
|
|
# -- revolved longitudinal sections ------------------------------------------
|
|
|
|
|
|
def test_a_solid_shaft_full_section_hatches_both_halves_to_the_axis():
|
|
loops = cut_loops(shaft(), axis_plane(), style="full")
|
|
assert len(loops) == 2
|
|
assert all(hatched for _, hatched in loops)
|
|
upper, lower = (points for points, _ in loops)
|
|
assert max(y for _, y in upper) == 15.0
|
|
assert min(y for _, y in upper) == 0.0
|
|
assert min(y for _, y in lower) == -15.0
|
|
# the cut face is the real cross-section: 20x12.5 + 40x15 + 15x10 per half
|
|
assert polygon_area(upper) == pytest.approx(20 * 12.5 + 40 * 15.0 + 15 * 10.0)
|
|
|
|
|
|
def test_a_bore_is_not_hatched_because_the_loop_closes_on_it():
|
|
part = RevolvedPart(name="sleeve", segments=(Segment(30.0, 40.0, 20.0),))
|
|
(upper, _), (lower, _) = cut_loops(part, axis_plane(), style="full")
|
|
assert polygon_area(upper) == pytest.approx(30.0 * (20.0 - 10.0))
|
|
assert min(y for _, y in upper) == 10.0
|
|
assert max(y for _, y in lower) == -10.0
|
|
|
|
|
|
def test_a_half_section_hatches_one_half_only():
|
|
loops = cut_loops(shaft(), axis_plane(style="half"), style="half")
|
|
assert len(loops) == 1
|
|
assert min(y for _, y in loops[0][0]) == 0.0
|
|
|
|
|
|
def test_iso6410_an_internal_thread_hatch_runs_to_the_major_diameter():
|
|
thread = Thread(id="th", segment=0, designation="M20", length=25.0, internal=True)
|
|
# drilled to the ISO 6410 minor (0.8 x 20 = 16 mm); the hatch must stop at
|
|
# the 20 mm major, not at the drilled 16.
|
|
part = RevolvedPart(name="boss", segments=(Segment(30.0, 50.0, 16.0),), features=(thread,))
|
|
bands = radial_profile(part)
|
|
threaded = [b for b in bands if b.x1 <= 25.0 + EPS]
|
|
assert {round(b.r_in, 6) for b in threaded} == {10.0}
|
|
unthreaded = [b for b in bands if b.x0 >= 25.0 - EPS]
|
|
assert {round(b.r_in, 6) for b in unthreaded} == {8.0}
|
|
(upper, _), _ = cut_loops(part, axis_plane(), style="full")
|
|
# The cut face closes on the major over the threaded span - and on the real
|
|
# drilled bore over the 5 mm that is drilled but not threaded. Asserting the
|
|
# minimum over the *whole* loop would be asserting that the drill does not
|
|
# exist beyond the thread, which is a different (and wrong) drawing.
|
|
assert min(y for x, y in upper if x < 25.0 - EPS) == pytest.approx(10.0)
|
|
assert min(y for x, y in upper if x > 25.0 + EPS) == pytest.approx(8.0)
|
|
|
|
|
|
def test_a_section_draws_its_own_cut_boundary_continuous_and_never_dashed():
|
|
"""ISO 128-50: a surface lying in the cutting plane is bounded by a
|
|
continuous wide line, and ISO 128-3 leaves hidden detail out of a sectional
|
|
view. The bore wall is the case that catches a section built out of the
|
|
un-cut silhouette: there the material hides the bore, so the silhouette
|
|
draws it dashed - but y = +-10 is the hatch loop's own boundary, and a
|
|
dashed line drawn on the hatch boundary means the outline and the hatched
|
|
area disagree at identical coordinates.
|
|
"""
|
|
part = RevolvedPart(name="sleeve", segments=(Segment(30.0, 40.0, 20.0),))
|
|
view = build_view(part, "section", plane=axis_plane(), style="full")
|
|
assert "hidden" not in {getattr(prim, "role", None) for prim in view.prims}
|
|
bore = [
|
|
line
|
|
for line in view.prims
|
|
if isinstance(line, Line)
|
|
and abs(abs(line.p1[1]) - 10.0) < EPS
|
|
and abs(abs(line.p2[1]) - 10.0) < EPS
|
|
]
|
|
assert len(bore) == 2
|
|
assert {line.role for line in bore} == {"visible"}
|
|
(upper, _), (lower, _) = cut_loops(part, axis_plane(), style="full")
|
|
assert min(y for _, y in upper) == pytest.approx(10.0)
|
|
assert max(y for _, y in lower) == pytest.approx(-10.0)
|
|
|
|
|
|
def test_a_half_section_keeps_its_dashed_bore_only_on_the_half_it_does_not_cut():
|
|
part = RevolvedPart(name="sleeve", segments=(Segment(30.0, 40.0, 20.0),))
|
|
view = build_view(part, "section", plane=axis_plane(style="half"), style="half")
|
|
hidden = [prim for prim in view.prims if getattr(prim, "role", None) == "hidden"]
|
|
assert hidden, "the un-sectioned half still shows its bore, dashed"
|
|
assert all(line.p1[1] < 0.0 and line.p2[1] < 0.0 for line in hidden)
|
|
assert all(abs(abs(line.p1[1]) - 10.0) < EPS for line in hidden)
|
|
# the cut half carries the hatch, and nothing dashed on it
|
|
((upper, _),) = cut_loops(part, axis_plane(style="half"), style="half")
|
|
assert min(y for _, y in upper) == pytest.approx(10.0)
|
|
|
|
|
|
def test_a_chamfer_is_cut_out_of_the_hatched_area_exactly_where_it_is_drawn():
|
|
"""The loop is the emitted outline flattened, so a corner substitution
|
|
reaches the cut face too: a 2 mm 45 deg chamfer takes a 2 mm2 triangle off
|
|
the 1000 mm2 half face, and a loop built from the raw profile would hatch
|
|
over the chamfer it just drew."""
|
|
(upper, _), _lower = cut_loops(
|
|
shaft(Chamfer(id="c1", at="start", size=2.0)), axis_plane(), style="full"
|
|
)
|
|
assert polygon_area(upper) == pytest.approx(20 * 12.5 + 40 * 15.0 + 15 * 10.0 - 2.0)
|
|
|
|
|
|
def test_the_section_view_carries_a_hatch_area_per_loop_with_the_parts_material():
|
|
view = build_view(shaft(), "section", plane=axis_plane(), style="full")
|
|
hatches = [p for p in view.prims if isinstance(p, HatchArea)]
|
|
assert len(hatches) == 2
|
|
assert {h.material for h in hatches} == {"steel"}
|
|
assert view.label == "A-A"
|
|
|
|
|
|
def test_an_unknown_material_is_refused_before_a_section_is_built():
|
|
part = RevolvedPart(name="x", segments=(Segment(10.0, 20.0),), material="unobtainium")
|
|
with pytest.raises(ValueError, match="unobtainium"):
|
|
build_view(part, "section", plane=axis_plane())
|
|
|
|
|
|
def test_an_offset_plane_on_a_revolved_part_is_refused_with_the_alternative():
|
|
with pytest.raises(ValueError) as excinfo:
|
|
build_view(shaft(), "section", plane=axis_plane(style="offset"), style="offset")
|
|
assert "style='full'" in str(excinfo.value)
|
|
|
|
|
|
# -- revolved transverse (cross) sections ------------------------------------
|
|
|
|
|
|
def test_a_transverse_section_hatches_the_annulus_at_that_station():
|
|
part = RevolvedPart(name="sleeve", segments=(Segment(30.0, 40.0, 20.0),))
|
|
plane = normalise_plane({"p1": (15.0, -30.0), "p2": (15.0, 30.0), "label": "B"})
|
|
view = build_view(part, "section", plane=plane, style="full")
|
|
hatches = [p for p in view.prims if isinstance(p, HatchArea)]
|
|
assert len(hatches) == 1
|
|
outer, island = hatches[0].loops
|
|
assert len(outer) == CIRCLE_SEGMENTS and len(island) == CIRCLE_SEGMENTS
|
|
# An inscribed regular CIRCLE_SEGMENTS-gon is exactly
|
|
# 1 - (n / 2pi) sin(2pi / n) = 2.0307e-4 under the circle, so the tolerance
|
|
# is the flattening the constant declares and nothing looser.
|
|
flattening = 1.0 - (CIRCLE_SEGMENTS / (2.0 * math.pi)) * math.sin(
|
|
2.0 * math.pi / CIRCLE_SEGMENTS
|
|
)
|
|
assert flattening == pytest.approx(2.0307e-4, rel=1e-3)
|
|
assert polygon_area(outer) == pytest.approx(math.pi * 20.0**2, rel=2.1e-4)
|
|
assert polygon_area(island) == pytest.approx(math.pi * 10.0**2, rel=2.1e-4)
|
|
assert view.label == "B-B"
|
|
circles = [p for p in view.prims if isinstance(p, Circle)]
|
|
assert {round(c.radius, 6) for c in circles} == {20.0, 10.0}
|
|
|
|
|
|
def test_a_revolved_style_cross_section_draws_its_boundary_as_phantom():
|
|
part = RevolvedPart(name="sleeve", segments=(Segment(30.0, 40.0, 20.0),))
|
|
plane = normalise_plane({"p1": (15.0, -30.0), "p2": (15.0, 30.0)})
|
|
view = build_view(part, "section", plane=plane, style="revolved")
|
|
circles = [p for p in view.prims if isinstance(p, Circle)]
|
|
assert {c.role for c in circles} == {"phantom"}
|
|
|
|
|
|
def keyed_shaft(*features) -> RevolvedPart:
|
|
"""A plain 40 mm shaft; DIN 6885-1 gives its keyway 12 mm x 5 mm deep."""
|
|
return RevolvedPart(name="keyed", segments=(Segment(60.0, 40.0),), features=tuple(features))
|
|
|
|
|
|
def test_a_cross_section_through_a_keyway_cuts_the_notch_out_of_the_disc():
|
|
"""The commonest mechanical cross section there is. The plane at x=20 is
|
|
squarely inside the slot, so the face is a notched disc - not the plain
|
|
hatched circle a face built from the segment diameter alone produces."""
|
|
keyway = Keyway(id="kw", segment=0, length=30.0, offset=5.0)
|
|
plane = normalise_plane({"p1": (20.0, -40.0), "p2": (20.0, 40.0)})
|
|
view = build_view(keyed_shaft(keyway), "section", plane=plane, style="full")
|
|
assert view.omitted == ()
|
|
arcs = [prim for prim in view.prims if isinstance(prim, Arc)]
|
|
assert len(arcs) == 1 and arcs[0].radius == pytest.approx(20.0)
|
|
assert not [prim for prim in view.prims if isinstance(prim, Circle)]
|
|
|
|
radius, half, floor = 20.0, 6.0, 15.0
|
|
top = math.sqrt(radius * radius - half * half)
|
|
# the disc, less the notch rectangle, less the cap the slot mouth removes
|
|
cap = radius * radius * math.acos(top / radius) - top * half
|
|
expected = math.pi * radius * radius - 2.0 * half * (top - floor) - cap
|
|
(hatch,) = [prim for prim in view.prims if isinstance(prim, HatchArea)]
|
|
assert polygon_area(hatch.loops[0]) == pytest.approx(expected, rel=2.1e-4)
|
|
# and it really is a notch: the whole disc is 56 mm2 larger
|
|
assert polygon_area(hatch.loops[0]) < math.pi * radius * radius - 50.0
|
|
# the arc the face is bounded by is the arc the end view draws
|
|
assert arcs[0].start_deg == pytest.approx(math.degrees(math.atan2(top, -half)) % 360.0)
|
|
assert arcs[0].end_deg == pytest.approx(math.degrees(math.atan2(top, half)) % 360.0)
|
|
|
|
|
|
def test_a_cross_section_reports_the_feature_it_cannot_project():
|
|
"""The honesty rule, on the one view where silence is indistinguishable
|
|
from a correct drawing: the face is hatched either way."""
|
|
cross = CrossHole(id="xh", x=20.0, diameter=6.0)
|
|
plane = normalise_plane({"p1": (20.0, -40.0), "p2": (20.0, 40.0)})
|
|
view = build_view(keyed_shaft(cross), "section", plane=plane, style="full")
|
|
assert [entry["feature"] for entry in view.omitted] == ["xh"]
|
|
assert "transverse" in view.omitted[0]["reason"]
|
|
|
|
|
|
def test_a_feature_the_transverse_plane_misses_is_reported_as_off_the_plane():
|
|
keyway = Keyway(id="kw", segment=0, length=30.0, offset=5.0)
|
|
plane = normalise_plane({"p1": (50.0, -40.0), "p2": (50.0, 40.0)})
|
|
view = build_view(keyed_shaft(keyway), "section", plane=plane, style="full")
|
|
assert [entry["feature"] for entry in view.omitted] == ["kw"]
|
|
assert "not on the cutting plane" in view.omitted[0]["reason"]
|
|
|
|
|
|
def test_a_transverse_face_is_read_off_the_radial_profile_not_the_segment():
|
|
"""A groove is in the profile, so it is in the cross section - at its own
|
|
reduced diameter, and without being reported missing."""
|
|
groove = RetainingGroove(id="rg", x=30.0, d=40.0, m=1.85, d2=37.5)
|
|
plane = normalise_plane({"p1": (30.0, -40.0), "p2": (30.0, 40.0)})
|
|
view = build_view(keyed_shaft(groove), "section", plane=plane, style="full")
|
|
assert view.omitted == ()
|
|
circles = [prim for prim in view.prims if isinstance(prim, Circle)]
|
|
assert {round(circle.radius, 6) for circle in circles} == {37.5 / 2.0}
|
|
|
|
|
|
def test_a_half_transverse_section_is_refused_with_its_reason():
|
|
part = RevolvedPart(name="sleeve", segments=(Segment(30.0, 40.0, 20.0),))
|
|
plane = normalise_plane({"p1": (15.0, -30.0), "p2": (15.0, 30.0)})
|
|
with pytest.raises(ValueError, match="half section"):
|
|
build_view(part, "section", plane=plane, style="half")
|
|
|
|
|
|
# -- prismatic sections ------------------------------------------------------
|
|
|
|
|
|
def test_a_prismatic_full_section_is_the_chord_times_the_thickness():
|
|
plane = normalise_plane({"p1": (-10.0, 25.0), "p2": (100.0, 25.0), "label": "C"})
|
|
view = build_view(plate(), "section", plane=plane, style="full")
|
|
(hatch,) = [p for p in view.prims if isinstance(p, HatchArea)]
|
|
assert polygon_area(hatch.loops[0]) == pytest.approx(80.0 * 8.0)
|
|
assert view.bbox == (0.0, 0.0, 80.0, 8.0)
|
|
|
|
|
|
def test_a_through_hole_splits_the_cut_face_in_two():
|
|
hole = Hole(id="h1", x=40.0, y=25.0, diameter=10.0)
|
|
plane = normalise_plane({"p1": (-10.0, 25.0), "p2": (100.0, 25.0)})
|
|
view = build_view(plate(hole), "section", plane=plane, style="full")
|
|
hatches = [p for p in view.prims if isinstance(p, HatchArea)]
|
|
assert len(hatches) == 2
|
|
total = sum(polygon_area(h.loops[0]) for h in hatches)
|
|
assert total == pytest.approx((80.0 - 10.0) * 8.0)
|
|
|
|
|
|
@pytest.mark.parametrize("offset", [0.0, 4.0, 4.9])
|
|
def test_a_hole_is_cut_on_its_chord_and_not_on_the_box_around_it(offset):
|
|
"""The gap a hole punches in a cut face is the chord the plane actually
|
|
cuts, and it shrinks to nothing as the plane leaves the centre. Removing
|
|
the hole by its bounding box reports the full 10 mm diameter at every
|
|
offset: 32 mm2 of missing cut face at 4 mm off centre and 64 mm2 at 4.9,
|
|
with an 8 mm gap drawn through solid metal where the plane merely grazes
|
|
the hole."""
|
|
hole = Hole(id="h1", x=40.0, y=25.0, diameter=10.0)
|
|
plane = normalise_plane({"p1": (-10.0, 25.0 + offset), "p2": (100.0, 25.0 + offset)})
|
|
view = build_view(plate(hole), "section", plane=plane, style="full")
|
|
chord = 2.0 * math.sqrt(5.0**2 - offset**2)
|
|
total = sum(polygon_area(prim.loops[0]) for prim in view.prims if isinstance(prim, HatchArea))
|
|
assert total == pytest.approx((80.0 - chord) * 8.0)
|
|
if offset > EPS:
|
|
assert total != pytest.approx((80.0 - 10.0) * 8.0)
|
|
|
|
|
|
def test_a_plane_tangent_to_a_hole_leaves_the_cut_face_whole():
|
|
hole = Hole(id="h1", x=40.0, y=25.0, diameter=10.0)
|
|
plane = normalise_plane({"p1": (-10.0, 30.0), "p2": (100.0, 30.0)})
|
|
view = build_view(plate(hole), "section", plane=plane, style="full")
|
|
hatches = [prim for prim in view.prims if isinstance(prim, HatchArea)]
|
|
assert len(hatches) == 1
|
|
assert polygon_area(hatches[0].loops[0]) == pytest.approx(80.0 * 8.0)
|
|
|
|
|
|
def test_an_angled_slot_is_cut_on_its_obround_and_not_on_the_box_around_it():
|
|
"""A 45 degree slot's bounding box is 24.14 mm across the cutting line; the
|
|
slot itself is 14.14. The box would erase 10 mm of material the slot never
|
|
touches."""
|
|
slot = Slot(id="s1", x=40.0, y=25.0, length=30.0, width=10.0, angle=45.0)
|
|
plane = normalise_plane({"p1": (-10.0, 25.0), "p2": (100.0, 25.0)})
|
|
view = build_view(plate(slot), "section", plane=plane, style="full")
|
|
total = sum(polygon_area(prim.loops[0]) for prim in view.prims if isinstance(prim, HatchArea))
|
|
assert total == pytest.approx((80.0 - 10.0 * math.sqrt(2.0)) * 8.0)
|
|
|
|
|
|
def test_a_pocket_reduces_the_cut_face_to_the_remaining_thickness():
|
|
pocket = Pocket(
|
|
id="p1",
|
|
outline=((20.0, 15.0), (60.0, 15.0), (60.0, 35.0), (20.0, 35.0)),
|
|
depth=3.0,
|
|
)
|
|
plane = normalise_plane({"p1": (-10.0, 25.0), "p2": (100.0, 25.0)})
|
|
view = build_view(plate(pocket), "section", plane=plane, style="full")
|
|
hatches = [p for p in view.prims if isinstance(p, HatchArea)]
|
|
total = sum(polygon_area(h.loops[0]) for h in hatches)
|
|
assert total == pytest.approx(40.0 * 8.0 + 40.0 * 5.0)
|
|
|
|
|
|
def test_iso128_3_a_rib_cut_longitudinally_is_drawn_but_not_hatched():
|
|
rib = Pocket(
|
|
id="rib",
|
|
outline=((20.0, 15.0), (60.0, 15.0), (60.0, 35.0), (20.0, 35.0)),
|
|
depth=3.0,
|
|
no_section_hatch=True,
|
|
)
|
|
plane = normalise_plane({"p1": (-10.0, 25.0), "p2": (100.0, 25.0)})
|
|
loops = cut_loops(plate(rib), plane, style="full")
|
|
hatched = [points for points, flag in loops if flag]
|
|
unhatched = [points for points, flag in loops if not flag]
|
|
assert len(unhatched) == 1
|
|
assert polygon_area(unhatched[0]) == pytest.approx(40.0 * 5.0)
|
|
assert sum(polygon_area(p) for p in hatched) == pytest.approx(40.0 * 8.0)
|
|
view = build_view(plate(rib), "section", plane=plane, style="full")
|
|
assert len([p for p in view.prims if isinstance(p, HatchArea)]) == 2
|
|
|
|
|
|
def test_an_offset_prismatic_section_accumulates_along_the_dog_leg():
|
|
plane = normalise_plane(
|
|
{
|
|
"p1": (-10.0, 15.0),
|
|
"p2": (100.0, 35.0),
|
|
"via": [[30.0, 15.0], [30.0, 35.0]],
|
|
"style": "offset",
|
|
}
|
|
)
|
|
view = build_view(plate(), "section", plane=plane, style="offset")
|
|
hatches = [p for p in view.prims if isinstance(p, HatchArea)]
|
|
total = sum(polygon_area(h.loops[0]) for h in hatches)
|
|
assert total == pytest.approx(80.0 * 8.0)
|
|
|
|
|
|
def test_a_revolved_style_section_of_a_prismatic_part_is_refused():
|
|
plane = normalise_plane({"p1": (-10.0, 25.0), "p2": (100.0, 25.0)})
|
|
with pytest.raises(ValueError, match="revolved section"):
|
|
build_view(plate(), "section", plane=plane, style="revolved")
|
|
|
|
|
|
def test_a_plane_that_misses_the_part_is_refused_rather_than_hatched_empty():
|
|
plane = normalise_plane({"p1": (-10.0, 200.0), "p2": (100.0, 200.0)})
|
|
with pytest.raises(ValueError, match="does not cut"):
|
|
build_view(plate(), "section", plane=plane, style="full")
|
|
|
|
|
|
# -- detail views ------------------------------------------------------------
|
|
|
|
|
|
def test_a_detail_scales_the_region_and_labels_both_places():
|
|
detail = {"center": (20.0, 12.5), "radius": 8.0, "scale": 5.0, "label": "D"}
|
|
view = build_view(shaft(), "detail", detail=detail)
|
|
assert view.kind == "detail"
|
|
assert view.scale == 5.0
|
|
assert view.label == "D (5:1)"
|
|
width = view.bbox[2] - view.bbox[0]
|
|
# straight geometry is trimmed to the circle, so the enlargement can never
|
|
# be wider than the circle itself at the detail scale
|
|
assert width <= 2.0 * 8.0 * 5.0 + 1e-6
|
|
assert width > 8.0 * 5.0
|
|
marker = detail_marker_prims(detail)
|
|
circle = next(p for p in marker if isinstance(p, Circle))
|
|
text = next(p for p in marker if isinstance(p, Text))
|
|
assert circle.center == (20.0, 12.5) and circle.radius == 8.0
|
|
assert text.text == "D"
|
|
|
|
|
|
def test_a_detail_reports_the_features_that_fall_outside_its_circle():
|
|
from engineering.mech.features import RetainingGroove
|
|
|
|
groove = RetainingGroove(id="rg", x=70.0, d=20.0, m=1.6, d2=18.6)
|
|
detail = {"center": (5.0, 12.5), "radius": 6.0, "scale": 5.0, "label": "D"}
|
|
view = build_view(shaft(groove), "detail", detail=detail)
|
|
assert any(entry["feature"] == "rg" for entry in view.omitted)
|
|
assert "detail circle D" in next(
|
|
entry["reason"] for entry in view.omitted if entry["feature"] == "rg"
|
|
)
|
|
|
|
|
|
def test_a_detail_without_a_radius_or_with_a_bad_scale_is_refused():
|
|
with pytest.raises(ValueError, match="radius"):
|
|
build_view(shaft(), "detail", detail={"center": (20.0, 12.5), "scale": 5.0})
|
|
with pytest.raises(ValueError, match="scale"):
|
|
build_view(shaft(), "detail", detail={"center": (20.0, 12.5), "radius": 8.0, "scale": 0.0})
|
|
|
|
|
|
def test_a_section_without_a_plane_is_refused_by_name():
|
|
with pytest.raises(ValueError, match="plane"):
|
|
build_view(shaft(), "section")
|
|
|
|
|
|
def test_an_axial_bore_shows_its_flat_bottom_in_the_section():
|
|
bore = AxialBore(id="ab", at="start", diameter=8.0, depth=30.0)
|
|
view = build_view(shaft(bore), "section", plane=axis_plane(), style="full")
|
|
(hatch_upper, _hatch_lower) = [p for p in view.prims if isinstance(p, HatchArea)]
|
|
assert min(y for _, y in hatch_upper.loops[0]) == 0.0
|
|
assert 4.0 in {round(y, 6) for _, y in hatch_upper.loops[0]}
|