From ddd0641649f0a76b4bb8aa4d7f36dfc3c840e930 Mon Sep 17 00:00:00 2001 From: jaseg Date: Wed, 14 Jun 2023 11:39:44 +0200 Subject: basic coils working --- coil_gen.py | 86 ++++++++++++++++++++++++++++++++++++------------------------- 1 file changed, 51 insertions(+), 35 deletions(-) diff --git a/coil_gen.py b/coil_gen.py index 642a77d..78610d2 100644 --- a/coil_gen.py +++ b/coil_gen.py @@ -73,16 +73,18 @@ def generate(infile, outfile, polygon, start_angle, stop_radius, trace_width, cl segment_heights = [point_line_distance((cx, cy), (x1, y1), (x2, y2)) for (x1, y1), (x2, y2) in segments] segment_foo = list(zip(segment_heights, segments)) - closest_points = [] + midpoints = [] + for h, ((x1, y1), (x2, y2)) in segment_foo: + xb = (x1 + x2) / 2 + yb = (y1 + y2) / 2 + midpoints.append((xb, yb)) + + normals = [] for h, ((x1, y1), (x2, y2)) in segment_foo: - dc1 = dist((x1, y1), (cx, cy)) d12 = dist((x1, y1), (x2, y2)) - db = sqrt(dc1**2 - h**2) - xn = (x2 - x1) / d12 - yn = (y2 - y1) / d12 - xb = x1 + xn * db - yb = y1 + yn * db - closest_points.append((xb, yb)) + dx = x2 - x1 + dy = y2 - y1 + normals.append((-dy/d12, dx/d12)) smallest_radius = min(segment_heights) #trace_radius = smallest_radius - stop_radius @@ -90,24 +92,21 @@ def generate(infile, outfile, polygon, start_angle, stop_radius, trace_width, cl num_windings = floor((trace_radius - trace_width) / (clearance + trace_width)) print(f'Going for {num_windings} windings') - segment_foo = list(zip(segment_heights, segments, segment_angles, closest_points)) + segment_foo = list(zip(segment_heights, segments, segment_angles, midpoints, normals)) - dbg_lines = [] + dbg_lines1, dbg_lines2 = [], [] spiral_points = [] dr_tot = 0 for n in range(num_windings): - for (ha, (pa1, pa2), aa, ma), (hb, (pb1, pb2), ab, mb) in zip(segment_foo[-1:] + segment_foo[:-1], segment_foo): + for (ha, (pa1, pa2), aa, ma, na), (hb, (pb1, pb2), ab, mb, nb) in zip(segment_foo[-1:] + segment_foo[:-1], segment_foo): pitch = clearance + trace_width dr_tot_a = dr_tot dr_tot_b = dr_tot + ab/(2*pi) * pitch xma, yma = ma + xna, yna = na xmb, ymb = mb - - xra = (xma - cx) / ha - yra = (yma - cy) / ha - xrb = (xmb - cx) / hb - yrb = (ymb - cy) / hb + xnb, ynb = nb xa1, ya1 = pa1 xa2, ya2 = pa2 @@ -117,27 +116,39 @@ def generate(infile, outfile, polygon, start_angle, stop_radius, trace_width, cl dma = dist(pa2, ma) dmb = dist(pb1, mb) - qa = dr_tot_a*dma/ha - dra = hypot(dr_tot_a, qa) - xea = xa2 + (cx - xa2) / dist((cx, cy), pa2) * dra - yea = ya2 + (cy - ya2) / dist((cx, cy), pa2) * dra + x_cons_a, y_cons_a = p_cons_a = line_line_intersection((pa2, (cx, cy)), (ma, (xma-xna, yma-yna))) + d_cons_a = dist(p_cons_a, ma) + qa = dma * dr_tot_a / d_cons_a + dra = hypot(qa, dr_tot_a) + + nrax = (xa2 - cx) / dist((cx, cy), pa2) + nray = (ya2 - cy) / dist((cx, cy), pa2) + + xea = xa2 - nrax*dra + yea = ya2 - nray*dra + + x_cons_b, y_cons_b = p_cons_b = line_line_intersection((pb1, (cx, cy)), (mb, (xmb-xnb, ymb-ynb))) + d_cons_b = dist(p_cons_b, mb) + qb = dmb * dr_tot_b / d_cons_b + drb = hypot(qb, dr_tot_b) + + nrbx = (xb1 - cx) / dist((cx, cy), pb1) + nrby = (yb1 - cy) / dist((cx, cy), pb1) - qb = dr_tot_b*dmb/hb - drb = hypot(dr_tot_b, qb) - xeb = xb1 + (cx - xb1) / dist((cx, cy), pb1) * drb - yeb = yb1 + (cy - yb1) / dist((cx, cy), pb1) * drb + xeb = xb1 - nrbx*drb + yeb = yb1 - nrby*drb - xsa = xma - xra*dr_tot_a - ysa = yma - yra*dr_tot_a + xsa = xma - xna*dr_tot_a + ysa = yma - yna*dr_tot_a - xsb = xmb - xrb*dr_tot_b - ysb = ymb - yrb*dr_tot_b + xsb = xmb - xnb*dr_tot_b + ysb = ymb - ynb*dr_tot_b l1 = (xsa, ysa), (xea, yea) l2 = (xsb, ysb), (xeb, yeb) - dbg_lines.append(l1) - dbg_lines.append(l2) + dbg_lines1.append(l1) + dbg_lines2.append(l2) pic = line_line_intersection(l1, l2) spiral_points.append(pic) @@ -172,13 +183,18 @@ def generate(infile, outfile, polygon, start_angle, stop_radius, trace_width, cl f.write(f'\n') f.write(f'\n') f.write(f'\n') - f.write(f'\n') - for x, y in closest_points: - f.write(f'\n') + + for (x1, y1), (x2, y2) in dbg_lines1: + f.write(f'') + + for (x1, y1), (x2, y2) in dbg_lines2: + f.write(f'') + + for x, y in midpoints: f.write(f'') + f.write(f'\n') - for (x1, y1), (x2, y2) in dbg_lines: - f.write(f'') + f.write(f'\n') f.write('\n') if __name__ == '__main__': -- cgit