Hi, I've written before about making a simulation of a levitation device
with diamagnetic materials in FEMM. I was interested in determing the
equilibrium position of the levitating magnet by calculating the force for a
number of positions between the diamagnetic plates.
Following the advice of the FEMM manual, I used the coenergy method for
calculating the field energy, and then looking for the force from that. But
what I found was that the Coenergy displayed no observable trend with height
of the magnet. On a lark, I tried the same simulation and calculated the
Force from the stress tensor, using a fixed contour. Surprisingly, I got
very believable results that showed the stable and unstable equilibrium
points for the levitating magnet quite nicely. I was content with that, but
I'm still intrigued why the coenergy method doesn't seem to be working in
this case. We tried moving the outer boundary to increase the area for the
problem. I tried increasing the mesh density in the region where the force
on the magnet is being calculated. But none of those things seems to be able
to get the coenergy method to provide meaningful results. I might be tempted
to say that my stress tensor calculations are just misleading, but they
mirror so nicely the behavior of the actual system that I find that hard to
believe.
Would anybody be willing to take a look at the simulation I've put together
and make some suggestions as to where the coenergy calculation might be
messing up, or where I could improve the model to get the "right" behavior?
I've attached the femm file, the associated pre and post lua files, and an
excel file with the coenergy numbers, and graphs of both the coenergy and
stress tensor calculations of the force with height (gravity is included
manually).
If you look at the Excel file, you can't possibly discern a pattern of
behavior from the coenergy graph, while the stress tensor graph is pretty
good, though not as smooth as you might like.
Any help would be appreciated.
Thanks,
Jim
------=_NextPart_001_001A_01C2C870.9B82AD00
Content-Type: text/html;
charset="US-ASCII"
Content-Transfer-Encoding: quoted-printable
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML xmlns=3D"http://www.w3.org/TR/REC-html40"; xmlns:o =3D=20
"urn:schemas-microsoft-com:office:office" xmlns:w =3D=20
"urn:schemas-microsoft-com:office:word"><HEAD>
<META http-equiv=3DContent-Type content=3D"text/html; charset=3Dus-ascii">
<META content=3DWord.Document name=3DProgId>
<META content=3D"MSHTML 6.00.2722.900" name=3DGENERATOR>
<META content=3D"Microsoft Word 10" name=3DOriginator><LINK=20
href=3D"cid:filelist.xml@xxxxxxxxxxxxxxxxx"; rel=3DFile-List><!--[if gte mso=
9]><xml>
<o:OfficeDocumentSettings>
<o:DoNotRelyOnCSS/>
</o:OfficeDocumentSettings>
</xml><![endif]--><!--[if gte mso 9]><xml>
<w:WordDocument>
<w:SpellingState>Clean</w:SpellingState>
<w:GrammarState>Clean</w:GrammarState>
<w:DocumentKind>DocumentEmail</w:DocumentKind>
<w:EnvelopeVis/>
<w:DrawingGridHorizontalSpacing>2.85 pt</w:DrawingGridHorizontalSpacing>
<w:DisplayVerticalDrawingGridEvery>2</w:DisplayVerticalDrawingGridEvery>
<w:Compatibility>
<w:BreakWrappedTables/>
<w:SnapToGridInCell/>
<w:WrapTextWithPunct/>
<w:UseAsianBreakRules/>
</w:Compatibility>
<w:BrowserLevel>MicrosoftInternetExplorer4</w:BrowserLevel>
</w:WordDocument>
</xml><![endif]-->
<STYLE>@page Section1 {size: 8.5in 11.0in; margin: 1.0in 1.0in .5in 1.0in; =
mso-header-margin: .5in; mso-footer-margin: .5in; mso-paper-source: 0; }
P.MsoNormal {
FONT-SIZE: 12pt; MARGIN: 0in 0in 0pt; FONT-FAMILY: "Times New Roman"; mso-=
style-parent: ""; mso-pagination: widow-orphan; mso-fareast-font-family: "T=
imes New Roman"
}
LI.MsoNormal {
FONT-SIZE: 12pt; MARGIN: 0in 0in 0pt; FONT-FAMILY: "Times New Roman"; mso-=
style-parent: ""; mso-pagination: widow-orphan; mso-fareast-font-family: "T=
imes New Roman"
}
DIV.MsoNormal {
FONT-SIZE: 12pt; MARGIN: 0in 0in 0pt; FONT-FAMILY: "Times New Roman"; mso-=
style-parent: ""; mso-pagination: widow-orphan; mso-fareast-font-family: "T=
imes New Roman"
}
A:link {
COLOR: blue; TEXT-DECORATION: underline; text-underline: single
}
SPAN.MsoHyperlink {
COLOR: blue; TEXT-DECORATION: underline; text-underline: single
}
A:visited {
COLOR: purple; TEXT-DECORATION: underline; text-underline: single
}
SPAN.MsoHyperlinkFollowed {
COLOR: purple; TEXT-DECORATION: underline; text-underline: single
}
SPAN.EmailStyle17 {
COLOR: windowtext; FONT-FAMILY: Arial; mso-style-type: personal-compose; m=
so-style-noshow: yes; mso-ansi-font-size: 10.0pt; mso-bidi-font-size: 10.0p=
t; mso-ascii-font-family: Arial; mso-hansi-font-family: Arial; mso-bidi-fon=
t-family: Arial
}
DIV.Section1 {
page: Section1
}
</STYLE>
<!--[if gte mso 10]>
<style>
/* Style Definitions */=20
table.MsoNormalTable
{mso-style-name:"Table Normal";
mso-tstyle-rowband-size:0;
mso-tstyle-colband-size:0;
mso-style-noshow:yes;
mso-style-parent:"";
mso-padding-alt:0in 5.4pt 0in 5.4pt;
mso-para-margin:0in;
mso-para-margin-bottom:.0001pt;
mso-pagination:widow-orphan;
font-size:10.0pt;
font-family:"Times New Roman";}
</style>
<![endif]--></HEAD>
<BODY lang=3DEN-US style=3D"tab-interval: .5in" vLink=3Dpurple link=3Dblue>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN class=3D833462320-30=
012003>Hi,=20
I've written before about making a simulation of a levitation device with=20
diamagnetic materials in FEMM. I was interested in determing the=
=20
equilibrium position of the levitating magnet by calculating the force for =
a=20
number of positions between the diamagnetic=20
plates. </SPAN></FONT></DIV>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN=20
class=3D833462320-30012003></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN=20
class=3D833462320-30012003>Following the advice of the FEMM manual, I =
used=20
the coenergy method for calculating the field energy, and then looking for =
the=20
force from that. But what I found was that the Coenergy displayed no observ=
able=20
trend with height of the magnet. On a lark, I tried the same simulation and=
=20
calculated the Force from the stress tensor, using a fixed contour.=20
Surprisingly, I got very believable results that showed the stable and unst=
able=20
equilibrium points for the levitating magnet quite nicely. I was content wi=
th=20
that, but I'm still intrigued why the coenergy method doesn't seem to =
be=20
working in this case. We tried moving the outer boundary to increase the ar=
ea=20
for the problem. I tried increasing the mesh density in the region where th=
e=20
force on the magnet is being calculated. But none of those things seems to =
be=20
able to get the coenergy method to provide meaningful results. I might be=20
tempted to say that my stress tensor calculations are just misleading, but =
they=20
mirror so nicely the behavior of the actual system that I find that hard to=
=20
believe. </SPAN></FONT></DIV>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN=20
class=3D833462320-30012003></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN class=3D833462320-30=
012003>Would=20
anybody be willing to take a look at the simulation I've put together and m=
ake=20
some suggestions as to where the coenergy calculation might be messing up, =
or=20
where I could improve the model to get the "right" behavior?</SPAN></FONT><=
/DIV>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN=20
class=3D833462320-30012003></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN class=3D833462320-30=
012003>I've=20
attached the femm file, the associated pre and post lua files, and an excel=
file=20
with the coenergy numbers, and graphs of both the coenergy and stress tenso=
r=20
calculations of the force with height (gravity is included=20
manually).</SPAN></FONT></DIV>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN=20
class=3D833462320-30012003></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN class=3D833462320-30=
012003>If you=20
look at the Excel file, you can't possibly discern a pattern of behavior fr=
om=20
the coenergy graph, while the stress tensor graph is pretty good, though no=
t as=20
smooth as you might like. </SPAN></FONT></DIV>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN=20
class=3D833462320-30012003></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN class=3D833462320-30=
012003>Any=20
help would be appreciated.</SPAN></FONT></DIV>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN=20
class=3D833462320-30012003></SPAN></FONT> </DIV>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN=20
class=3D833462320-30012003>Thanks,</SPAN></FONT></DIV>
<DIV><FONT face=3DArial color=3D#0000ff size=3D2><SPAN=20
class=3D833462320-30012003>Jim</SPAN></FONT></DIV></BODY></HTML>
------=_NextPart_001_001A_01C2C870.9B82AD00--
Attachment:
bin00069.bin
Description: Binary data
Attachment:
bin00070.bin
Description: Binary data
Attachment:
bin00071.bin
Description: Binary data
Attachment:
xls00002.xls
Description: MS-Excel spreadsheet
Attachment:
xls00003.xls
Description: MS-Excel spreadsheet