def draw_FEMM(self,
nodeprop=None,
maxseg=None,
propname=None,
hide=False,
group=None):
"""Draw the Arc object in FEMM and assign the property
Parameters
----------
nodeprop :
Nodal property
(Default value = None)
maxseg :
Meshed with elements that span at most maxsegdeg degrees per element
(Default value = None)
propname :
Boundary property ’propname’
(Default value = None)
hide :
0 = not hidden in post-processor, 1 == hidden in post processor
(Default value = False)
group :
the group the Arc1 object belongs
(Default value = None)
Returns
-------
None
"""
# Get BC (if any)
if self.label in boundary_prop:
propname = boundary_prop[self.label]
# Create the nodes
begin = self.get_begin()
end = self.get_end()
X1, Y1 = begin.real, begin.imag
X2, Y2 = end.real, end.imag
femm.mi_addnode(X1, Y1)
femm.mi_selectnode(X1, Y1)
femm.mi_setnodeprop(nodeprop, group)
femm.mi_clearselected()
femm.mi_addnode(X2, Y2)
femm.mi_selectnode(X2, Y2)
femm.mi_setnodeprop(nodeprop, group)
femm.mi_clearselected()
# Create the arc
angle = self.get_angle(is_deg=True)
if angle > 0:
femm.mi_addarc(X1, Y1, X2, Y2, angle, 2)
else:
femm.mi_addarc(X2, Y2, X1, Y1, -angle, 2)
# Set Arc properties
Zm = self.get_middle()
femm.mi_selectarcsegment(Zm.real, Zm.imag)
femm.mi_setarcsegmentprop(maxseg, propname, hide, group)
femm.mi_clearselected()
def draw_FEMM(
self,
nodeprop=None,
propname=None,
elementsize=None,
automesh=None,
hide=False,
group=None,
):
"""< Draw the segment in FEMM and assign the property
Parameters
----------
nodeprop :
Nodal property
(Default value = None)
propname :
Boundary property ’propname’
(Default value = None)
elementsize :
Local element size along segment no greater than elementsize
(Default value = None)
automesh :
0 = mesher defers to the element constraint defined by
elementsize, 1 = mesher automatically chooses mesh size along
the selected segments
(Default value = None)
hide :
0 = not hidden in post-processor, 1 == hidden in post processorc
(Default value = False)
group :
group the segment belongs
(Default value = None)
Returns
-------
"""
# Get BC (if any)
if self.label in boundary_prop:
propname = boundary_prop[self.label]
# Add the nodes
X1, Y1 = self.begin.real, self.begin.imag
X2, Y2 = self.end.real, self.end.imag
femm.mi_addnode(X1, Y1)
femm.mi_selectnode(X1, Y1)
femm.mi_setnodeprop(nodeprop, group)
femm.mi_clearselected()
femm.mi_addnode(X2, Y2)
femm.mi_selectnode(X2, Y2)
femm.mi_setnodeprop(nodeprop, group)
femm.mi_clearselected()
# add the segment
femm.mi_addsegment(X1, Y1, X2, Y2)
# Set property
femm.mi_selectsegment((X1 + X2) / 2, (Y1 + Y2) / 2)
femm.mi_setsegmentprop(propname, elementsize, automesh, hide, group)
femm.mi_clearselected()
def draw_FEMM(self, nodeprop=None, maxseg=None, propname=None, hide=False, group=None):
"""Draw the Arc object in FEMM and assign the property
Parameters
----------
nodeprop :
Nodal property
(Default value = None)
maxseg :
Meshed with elements that span at most maxsegdeg degrees per element
(Default value = None)
propname :
Boundary property ’propname’
(Default value = None)
hide :
0 = not hidden in post-processor, 1 == hidden in post processor
(Default value = False)
group :
the group the Arc1 object belongs
(Default value = None)
Returns
-------
None
"""
# Get BC (if any)
for bound_label in boundary_prop:
if bound_label in self.label:
propname = boundary_prop[bound_label]
# split if arc angle > 180
angle = self.get_angle(is_deg=True)
# Create the nodes
begin = self.get_begin()
mid = self.get_middle()
end = self.get_end()
X1, Y1 = begin.real, begin.imag
X2, Y2 = mid.real, mid.imag
X3, Y3 = end.real, end.imag
femm.mi_addnode(X1, Y1)
femm.mi_selectnode(X1, Y1)
femm.mi_setnodeprop(nodeprop, group)
femm.mi_clearselected()
femm.mi_addnode(X3, Y3)
femm.mi_selectnode(X3, Y3)
femm.mi_setnodeprop(nodeprop, group)
femm.mi_clearselected()
if abs(angle) > 180:
femm.mi_addnode(X2, Y2)
femm.mi_selectnode(X2, Y2)
femm.mi_setnodeprop(nodeprop, group)
femm.mi_clearselected()
# invert the nodes
if angle < 0:
angle = -angle
X1, Y1, X3, Y3 = X3, Y3, X1, Y1
if angle > 180:
femm.mi_addarc(X1, Y1, X2, Y2, angle / 2, 2)
femm.mi_addarc(X2, Y2, X3, Y3, angle / 2, 2)
else:
femm.mi_addarc(X1, Y1, X3, Y3, angle, 2)
# Set Arc properties
if angle > 180:
cent = self.get_center()
mid1 = cent + (mid - cent) * exp(1j * angle / 4)
mid2 = cent + (mid - cent) * exp(-1j * angle / 4)
femm.mi_selectarcsegment(mid1.real, mid1.imag)
femm.mi_selectarcsegment(mid2.real, mid1.imag)
else:
femm.mi_selectarcsegment(X2, Y2)
femm.mi_setarcsegmentprop(maxseg, propname, hide, group)
femm.mi_clearselected()