<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 TRANSITIONAL//EN">
<HTML>
<HEAD>
<META HTTP-EQUIV="Content-Type" CONTENT="text/html; CHARSET=UTF-8">
<META NAME="GENERATOR" CONTENT="GtkHTML/3.18.3">
</HEAD>
<BODY>
In my last .pro, I (hopefully) used the integral over the whole _area_ of E=int dArea*mu^(-1)*0.5*B(r)^2.<BR>
I neglected the factor 0.5.<BR>
Honestly, I recommend to forget about batch files, install Cygwin or MinGW and use bash scripting in Windows.<BR>
You will have so much more possibilities.<BR>
<BR>
Also, I think that the h field inside is not correct in your outputted pos file.<BR>
greetings,<BR>
bernhard<BR>
<BR>
On Wed, 2008-07-23 at 11:07 +0200, Sebastian wrote:
<BLOCKQUOTE TYPE=CITE>
<PRE>
Dear Bernhard:
Today I am working in a BATCH file ( MSDOS script
file, I am using Gmsh in Windows ) for
moving the magnets ( Geometry ) and obtain data for each step.
Do you known the FEMM package ?
<A HREF="http://femm.foster-miller.net/wiki/HomePage">http://femm.foster-miller.net/wiki/HomePage</A>
This is a very interesting Freeware software, for 2D ElectroMagnetic problems.
In this software the force is calculated by integration
F=integration((H(B*n)+B(H*n).(H*B)n)/2)
where n denotes the direction normal to the surface at the point of interest.
This integration is calculated using the data in
a perimeter around the object that you want to
find the force. See the FEMM manual for more
details, about integration of Maxwell's Stress
Tensor, see page 44. I am sure that may interest you.
I want to find how to obtain this integration in GetDP.
In my webpage you can see how I have created a
Lua script to move magnets ( geometry ), and
rotor generator, and obtaining images and data in each step.
<A HREF="http://www.tinet.org/~sje/femm/acgen.htm">http://www.tinet.org/~sje/femm/acgen.htm</A>
Thanks.
At 10:15 23/07/2008, you wrote:
>Dear Sebastian,
>I just played around a bit with your files, and
>maybe I have a method that somehow can calculate the force.
>You create different meshes, where the magnets are varied in position.
>As an output for each mesh, you calculate the total field energy.
>the force should then be the variation of the
>field energy divided by the dislocation distance.
>but I have not verified things, only thing is I
>wrote an script that does the variation.
>after the script (that should work with linux
>after setting the paths correctly, or with cywin
>in windows) you can do the manual calculations to get the force.
>I _think_ that you just need to add the field
>energies of all domains. but what confuses me is
>that the magnets have negative energy?? I quite
>in a hurry now, and cannot investigate much further right now.
>nice greetings,
>bernhard
>
>On Tue, 2008-07-22 at 13:45 +0200, Sebastian wrote:
>>
>>
>>Dear friends:
>>
>>Here do you have a 2D magnetostatic simulation.
>>
>>The question is. How can I obtain the mechanical force in N, produced
>>by one magnet over the other ?
>>
>>Thanks.
>>
>>Sebastian Jardi Estadella
>><<A HREF="mailto:sjardi@telefonica.net">mailto:sjardi@telefonica.net</A>><A HREF="mailto:sjardi@telefonica.net">sjardi@telefonica.net</A>
>><A HREF="http://www.tinet.org/~sje/index_en.htm">http://www.tinet.org/~sje/index_en.htm</A>
>>_______________________________________________
>>getdp mailing list
>><<A HREF="mailto:getdp@geuz.org">mailto:getdp@geuz.org</A>><A HREF="mailto:getdp@geuz.org">getdp@geuz.org</A>
>><A HREF="http://www.geuz.org/mailman/listinfo/getdp">http://www.geuz.org/mailman/listinfo/getdp</A>
>________________________________________
>
>arsenal research
>
>Geschäftsfeld MEA /
>
>Business MEA
>
>DI Bernhard Kubicek
>
>Österreichisches Forschungs- und Prüfzentrum Arsenal Ges.m.b.H.
>FN 165088b; UID-Nr.: ATU 46577208
>Sitz der Gesellschaft: Wien; Gerichtsstand: Wien
>
>Austria, 1210 Wien, Giefinggasse 2
>
>ph: +43 (0) 50550-DW, f: +43 (0) 50550-DW
>
>mobile: +43 (0)664/ 12 34 567
>
><<A HREF="mailto:bernhard.kubicek@arsenal.ac.at">mailto:bernhard.kubicek@arsenal.ac.at</A>><A HREF="mailto:bernhard.kubicek@arsenal.ac.at">bernhard.kubicek@arsenal.ac.at</A>
>
>www.arsenal.ac.at
>
>
>
>
>
>
>
>
>
>
>
>
</PRE>
</BLOCKQUOTE>
</BODY>
</HTML>