summaryrefslogtreecommitdiff
diff options
context:
space:
mode:
-rw-r--r--coil_parasitics.py8
-rw-r--r--twisted_coil_gen_twolayer.py6
2 files changed, 10 insertions, 4 deletions
diff --git a/coil_parasitics.py b/coil_parasitics.py
index 5bb7b02..47d5cd2 100644
--- a/coil_parasitics.py
+++ b/coil_parasitics.py
@@ -73,14 +73,14 @@ def elmer_grid(infile, outfile=None, intype=None, outtype=None, cwd=None, **kwar
if outfile is not None:
kwargs['out'] = str(outfile)
- args = ['ElmerGrid', intype, outtype, infile]
+ args = ['ElmerGrid', intype, outtype, str(infile)]
for key, value in kwargs.items():
args.append(f'-{key}')
if isinstance(value, (tuple, list)):
args.extend(str(v) for v in value)
else:
args.append(str(value))
- subprocess.run(args, cwd=cwd)
+ subprocess.run(args, cwd=cwd, check=True)
def elmer_solver(cwd):
subprocess.run(['ElmerSolver'], cwd=cwd)
@@ -138,7 +138,7 @@ def run_capacitance_simulation(mesh_file, sim_dir):
sim.write_startinfo(tmpdir)
sim.write_sif(tmpdir)
# Convert mesh from gmsh to elemer formats. Also scale it from 1 unit = 1 mm to 1 unit = 1 m (SI units)
- elmer_grid(mesh_file.name, 'mesh', cwd=tmpdir, scale=[1e-3, 1e-3, 1e-3])
+ elmer_grid(mesh_file.absolute(), 'mesh', cwd=tmpdir, scale=[1e-3, 1e-3, 1e-3])
elmer_solver(tmpdir)
capacitance_matrix = np.loadtxt(tmpdir / 'capacitance.txt')
@@ -212,7 +212,7 @@ def run_inductance_simulation(mesh_file, sim_dir):
sim.write_startinfo(tmpdir)
sim.write_sif(tmpdir)
# Convert mesh from gmsh to elemer formats. Also scale it from 1 unit = 1 mm to 1 unit = 1 m (SI units)
- elmer_grid(mesh_file.name, 'mesh', cwd=tmpdir, scale=[1e-3, 1e-3, 1e-3])
+ elmer_grid(mesh_file.absolute(), 'mesh', cwd=tmpdir, scale=[1e-3, 1e-3, 1e-3])
elmer_solver(tmpdir)
diff --git a/twisted_coil_gen_twolayer.py b/twisted_coil_gen_twolayer.py
index f603cd0..1ecde20 100644
--- a/twisted_coil_gen_twolayer.py
+++ b/twisted_coil_gen_twolayer.py
@@ -253,6 +253,7 @@ def traces_to_gmsh_mag(traces, mesh_out, bbox, model_name='gerbonara_board', log
print('Cutting substrate')
occ.cut([(3, substrate)], [(3, toplevel_tag)], removeObject=True, removeTool=False)
+ occ.fragment([(3, substrate)], [(2, interface_tag_top), (2, interface_tag_bottom)])
occ.fragment([(3, airbox)], [(3, substrate), (3, toplevel_tag)])
print('Synchronizing')
@@ -263,6 +264,11 @@ def traces_to_gmsh_mag(traces, mesh_out, bbox, model_name='gerbonara_board', log
trace_physical = gmsh.model.add_physical_group(3, [toplevel_tag], name='trace')
interface_top_physical = gmsh.model.add_physical_group(2, [interface_tag_top], name='interface_top')
interface_bottom_physical = gmsh.model.add_physical_group(2, [interface_tag_bottom], name='interface_bottom')
+
+ print('first disk', first_disk)
+ print('bbox', occ.getBoundingBox(2, interface_tag_top))
+ print('last disk', last_disk)
+ print('bbox', occ.getBoundingBox(2, interface_tag_bottom))
airbox_adjacent = set(gmsh.model.getAdjacencies(3, airbox)[1])
in_bbox = {tag for _dim, tag in gmsh.model.getEntitiesInBoundingBox(x1+eps, y1+eps, z0+eps, x2-eps, y2-eps, z0+ab_h-eps, dim=2)}