[Getdp] vector normal to a magnet

Joseph Frippiat frippiat.j at kst.be
Thu Oct 17 13:34:19 CEST 2013


The development version of GetDP that I downloaded yesterday evening
corrects this problem.

 

Joseph

 

--

Joseph Frippiat
Senior Developer - UPS Division

KST S.A. - Site of Battice
Cour Lemaire, 16
B-4651 BATTICE
BELGIUM
Tel : +.32.4.247.01.67
Fax : +.32.87.31.42.21
E-mail : frippiat.j at kst.be
Web : www.kst-ups.com

Head office : Cour Lemaire, 16 - 4651 Battice (Belgium)
RPM Liège (Belgium) : 0423.472.603

 

From: Joseph Frippiat [mailto:frippiat.j at kst.be] 
Sent: mardi 15 octobre 2013 13:23
To: 'getdp at geuz.org'
Subject: RE: [Getdp] vector normal to a magnet

 

This doesn’t work (b = 0 everywhere):

  For k In {1:nbMagnets}

    br[ Rotor_Magnet~{k} ] = b_remanent * Vector[ 0, 1, 0 ];

  EndFor

 

This works (with the wrong amplitude) :

  For k In {1:nbMagnets}

    br[ Rotor_Magnet~{k} ] = b_remanent * Vector[ 0, Y[], 0 ];

  EndFor

 

Rotate doesn’t work for me (b = 0 everywhere).

 

Joseph

 

--

Joseph Frippiat
Senior Developer - UPS Division

KST S.A. - Site of Battice
Cour Lemaire, 16
B-4651 BATTICE
BELGIUM
Tel : +.32.4.247.01.67
Fax : +.32.87.31.42.21
E-mail : frippiat.j at kst.be
Web : www.kst-ups.com

Head office : Cour Lemaire, 16 - 4651 Battice (Belgium)
RPM Liège (Belgium) : 0423.472.603

 

From: Patrick Dular [mailto:Patrick.Dular at ulg.ac.be] 
Sent: mardi 15 octobre 2013 08:51
To: Joseph Frippiat
Cc: getdp at geuz.org
Subject: Re: [Getdp] vector normal to a magnet

 

Joseph,

Because Atan2[Y[],X[]]] gives you the angle made by the position vector (
X[] , Y [] ) (of each current point in the XY-plane, i.e. each point where
this function has to be evaluated) (see http://en.wikipedia.org/wiki/Atan2),
the Vector[] function

  Vector[ Cos[Atan2[Y[],X[]]], Sin[Atan2[Y[],X[]]], 0 ]

gives a radial vector at this point (which can be checked on your field map:
the direction of the field in the magnets is mainly the radial one.

Cos[Atan2[Y[],X[]]] is the X-projection of the unit radial vector.
Sin[Atan2[Y[],X[]]] is the Y-projection of the unit radial vector.

To define the direction you wish for the remanent induction in each magnet,
you should use the function I sent you yesterday. If you wish, send me your
.geo and .pro files and I will implement it for you.

Patrick


On 14/10/13 20:42, Joseph Frippiat wrote: 

I don't understand the Vector[] usage..

With this :

  For k In {1:nbMagnets}
    br[ Rotor_Magnet~{k} ] = b_remanent * Vector[ Cos[Atan2[Y[],X[]]],
Sin[Atan2[Y[],X[]]], 0 ];
  EndFor

I have a b field (not the wanted one) :



With this :

  For k In {1:nbMagnets}
    br[ Rotor_Magnet~{k} ] = b_remanent * Vector[ Cos[Atan2[Y[],X[]]+Pi/4],
Sin[Atan2[Y[],X[]]+Pi/4], 0 ];
  EndFor

I have nothing :



Why ?


Joseph

-- 

Joseph Frippiat
Senior Developer - UPS Division

KST S.A. - Site of Battice
Cour Lemaire, 16
B-4651 BATTICE
BELGIUM
Tel : +.32.4.247.01.67
Fax : +.32.87.31.42.21
E-mail : frippiat.j at kst.be
Web : www.kst-ups.com

Head office : Cour Lemaire, 16 - 4651 Battice (Belgium)
RPM Liège (Belgium) : 0423.472.603





_______________________________________________
getdp mailing list
getdp at geuz.org
http://www.geuz.org/mailman/listinfo/getdp

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20131017/320eafbe/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.png
Type: image/png
Size: 69929 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20131017/320eafbe/attachment.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image002.png
Type: image/png
Size: 15711 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20131017/320eafbe/attachment-0001.png>


More information about the getdp mailing list