zcy
/
openEMS
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

python tutorials: use new/better automesh options for CRLH examples 

Signed-off-by: Thorsten Liebig <Thorsten.Liebig@gmx.de>
pull/107/head
Thorsten Liebig 2023-01-07 20:52:36 +01:00
parent 595c8effbd
commit 0342eefd27
1 changed files with 7 additions and 16 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

23
python/Tutorials/CRLH_Extraction.py
View File

@ -18,6 +18,7 @@ from pylab import *
from CSXCAD import ContinuousStructure
from openEMS import openEMS
from openEMS.physical_constants import *
from openEMS.automesh import mesh_hint_from_box
### Class to represent single CRLH unit cells
class CRLH_Cells:
@ -42,39 +43,29 @@ class CRLH_Cells:
self.edge_resolution = res
def createCell(self, translate = [0,0,0]):
mesh = [[],[],None]
mesh = [None,None,None]
third_res = self.edge_resolution/3
translate = array(translate)
start = [-self.LL/2 , -self.LW/2, self.Top] + translate
stop = [-self.GLT/2, self.LW/2, self.Top] + translate
box = self.props['metal_top'].AddBox(start, stop, priority=10)
# create edge mesh at +x and -+y
mesh[0] += [stop[0]-third_res, stop[0]+2*third_res]
mesh[1] += [stop[1]-third_res, stop[1]+2*third_res]
mesh[1] += [start[1]-2*third_res, start[1]+third_res]
mesh = mesh_hint_from_box(box, 'x', metal_edge_res=self.edge_resolution, down_dir=False, mesh=mesh)
mesh = mesh_hint_from_box(box, 'y', metal_edge_res=self.edge_resolution, mesh=mesh)
start = [+self.LL/2 , -self.LW/2, self.Top] + translate
stop = [+self.GLT/2, self.LW/2, self.Top] + translate
box = self.props['metal_top'].AddBox(start, stop, priority=10)
# create edge mesh at -x
mesh[0] += [stop[0]-2*third_res, stop[0]+third_res]
mesh = mesh_hint_from_box(box, 'x', metal_edge_res=self.edge_resolution, up_dir=False, mesh=mesh)
start = [-(self.LL-self.GLB)/2, -self.LW/2, self.Bot] + translate
stop = [+(self.LL-self.GLB)/2, self.LW/2, self.Bot] + translate
box = self.props['metal_bot'].AddBox(start, stop, priority=10)
# create edge mesh at -+x
mesh[0] += [start[0]+third_res, start[0]-2*third_res]
mesh[0] += [stop[0]-third_res, stop[0]+2*third_res]
mesh = mesh_hint_from_box(box, 'x', metal_edge_res=self.edge_resolution, mesh=mesh)
start = [-self.SW/2, -self.LW/2-self.SL, self.Bot] + translate
stop = [+self.SW/2, self.LW/2+self.SL, self.Bot] + translate
box = self.props['metal_bot'].AddBox(start, stop, priority=10)
# create edge mesh at -+x and -+y
mesh[0] += [start[0]+third_res, start[0]-2*third_res]
mesh[0] += [stop[0]-third_res, stop[0]+2*third_res]
mesh[1] += [start[1]+third_res, start[1]-2*third_res]
mesh[1] += [stop[1]-third_res, stop[1]+2*third_res]
mesh = mesh_hint_from_box(box, 'xy', metal_edge_res=self.edge_resolution, mesh=mesh)
start = [0, -self.LW/2-self.SL+self.SW/2, 0 ] + translate
stop = [0, -self.LW/2-self.SL+self.SW/2, self.Bot] + translate