<DIV>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US style="mso-ansi-language: EN-US"><FONT size=3><FONT face="Times New Roman">Dear Christophe:<?xml:namespace prefix = o ns = "urn:schemas-microsoft-com:office:office" /><o:p></o:p></FONT></FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US style="mso-ansi-language: EN-US"><FONT size=3><FONT face="Times New Roman"> <o:p></o:p></FONT></FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US style="mso-ansi-language: EN-US"><FONT size=3><FONT face="Times New Roman">I sent this mail to Johan Gynselink,<SPAN style="COLOR: black"> </SPAN>but I do not have response.<SPAN style="COLOR: black"><o:p></o:p></SPAN></FONT></FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><FONT size=3><FONT face="Times New Roman"><SPAN lang=EN-US style="mso-ansi-language: EN-US">I have reviewed your study "Analysis of a Permanent Magnet Generator with Dual Formulations Using Periodicity Conditions and Moving Band" and I would like reconstructing a quarter motor geometry.</SPAN><SPAN lang=EN-US style="COLOR: black; mso-ansi-language: EN-US; mso-bidi-font-size: 10.0pt"> <o:p></o:p></SPAN></FONT></FONT></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US style="mso-ansi-language: EN-US"><FONT size=3><FONT face="Times New Roman">Until now I have utilized the function MB for all the motor and I have not problems, but requires larger calculation time.<SPAN style="COLOR: black"><o:p></o:p></SPAN></FONT></FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US style="COLOR: black; mso-ansi-language: EN-US"><FONT size=3><FONT face="Times New Roman"> <o:p></o:p></FONT></FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US style="COLOR: black; mso-ansi-language: EN-US; mso-bidi-font-size: 10.0pt"><FONT size=3><FONT face="Times New Roman">I wait for your comments.<o:p></o:p></FONT></FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US style="mso-ansi-language: EN-US"><FONT size=3><FONT face="Times New Roman"> <o:p></o:p></FONT></FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><FONT size=3><FONT face="Times New Roman">Best regard...</FONT></FONT></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><FONT size=3><FONT face="Times New Roman"></FONT></FONT> </P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><FONT size=3><FONT face="Times New Roman">Christian Geikowsky R.</FONT></FONT></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><FONT size=3><FONT face="Times New Roman"><SPAN lang=EN-US style="mso-ansi-language: EN-US"><o:p></o:p></SPAN></FONT></FONT> </P></DIV>
<DIV><B><I>"Christian Geikowsky R." <christiangeikowsky@yahoo.es></I></B> escribió:</DIV>
<BLOCKQUOTE class=replbq style="PADDING-LEFT: 5px; MARGIN-LEFT: 5px; BORDER-LEFT: #1010ff 2px solid">Fecha: Fri, 4 Nov 2005 14:35:58 +0100 (CET)<BR>De: "Christian Geikowsky R." <christiangeikowsky@yahoo.es><BR>Asunto: Moving Band<BR>Para: Johan Gyselinck <johan.gyselinck@ulb.ac.be><BR><BR>
<DIV>
<DIV>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><FONT size=3><FONT face="Times New Roman"><SPAN lang=EN-US style="mso-ansi-language: EN-US">Estimated Johan:</SPAN><SPAN lang=EN-US style="COLOR: black; mso-ansi-language: EN-US; mso-bidi-font-size: 10.0pt"><o:p></o:p></SPAN></FONT></FONT></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US style="COLOR: black; mso-ansi-language: EN-US; mso-bidi-font-size: 10.0pt"><o:p><FONT face="Times New Roman" size=3> </FONT></o:p></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><FONT size=3><FONT face="Times New Roman"><SPAN lang=EN-US style="mso-ansi-language: EN-US">Attachment a simple geometry of a fourth part of a 4 poles motor.</SPAN><SPAN lang=EN-US style="COLOR: black; mso-ansi-language: EN-US; mso-bidi-font-size: 10.0pt"><o:p></o:p></SPAN></FONT></FONT></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US style="COLOR: black; mso-ansi-language: EN-US; mso-bidi-font-size: 10.0pt"><FONT size=3><FONT face="Times New Roman">I would like know<o:p></o:p></FONT></FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US style="COLOR: black; mso-ansi-language: EN-US; mso-bidi-font-size: 10.0pt"><FONT size=3><FONT face="Times New Roman"><SPAN style="mso-spacerun: yes"> </SPAN><o:p></o:p></FONT></FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><FONT size=3><FONT face="Times New Roman"><SPAN lang=EN-US style="mso-ansi-language: EN-US">How must I apply the boundary conditions in the limits S000,</SPAN><SPAN lang=EN-US style="COLOR: black; mso-ansi-language: EN-US; mso-bidi-font-size: 10.0pt"> </SPAN><SPAN lang=EN-US style="mso-ansi-language: EN-US">S090,</SPAN><SPAN lang=EN-US style="COLOR: black; mso-ansi-language: EN-US; mso-bidi-font-size: 10.0pt"> </SPAN><SPAN lang=EN-US style="mso-ansi-language: EN-US">R000 and R090 and in the air gap?</SPAN><SPAN lang=EN-US style="COLOR: black; mso-ansi-language: EN-US; mso-bidi-font-size: 10.0pt"> <o:p></o:p></SPAN></FONT></FONT></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US style="COLOR: black; mso-ansi-language: EN-US; mso-bidi-font-size: 10.0pt"><o:p><FONT face="Times New Roman" size=3> </FONT></o:p></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US style="mso-ansi-language: EN-US"><FONT size=3><FONT face="Times New Roman">How must give form to the function MovingBand for a fourth of the geometry?<o:p></o:p></FONT></FONT></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><SPAN lang=EN-US style="mso-ansi-language: EN-US"><o:p><FONT face="Times New Roman" size=3> </FONT></o:p></SPAN></P>
<P class=MsoNormal style="MARGIN: 0cm 0cm 0pt"><FONT size=3><FONT face="Times New Roman">Best regard...<SPAN style="COLOR: black; mso-bidi-font-size: 10.0pt"><o:p></o:p></SPAN></FONT></FONT></P></DIV></DIV><BR><BR>Saludos...<BR><BR>Christian Geikowsky R.
<P>
<HR SIZE=1>
<BR><FONT face=Verdana size=-2>Correo Yahoo!<BR>Comprueba qué es nuevo, <A href="http://us.rd.yahoo.com/mail/es/whatsnew/*http://es.whatsnew.mail.yahoo.com/">aquí</A><BR>http://correo.yahoo.es</FONT>R1 = 0.5; // Radio rotor<BR>R2 = 0.6; // Radio interno Estator<BR><BR><BR>Point(1) = {0, 0, 0, 1}; <BR><BR>Point(2) = {1, 0, 0, 1}; // Circulo externo estator<BR>Point(3) = {0, 1, 0, 1};<BR><BR>Circle(1) = {2,1,3};<BR><BR>Point(12) = {R1, 0, 0, 1}; // Circulo rotor<BR>Point(13) = {0, R1, 0, 1};<BR><BR>Circle(11) = {12,1,13};<BR><BR>Point(22) = {R2, 0, 0, 1}; // Circulo interno estator<BR>Point(23) = {0, R2, 0, 1};<BR><BR>Circle(21) = {22,1,23};<BR><BR>Line(25) = {1,12}; // 0º Rotor<BR>Line(26) = {22,2}; // 0º Estator<BR>Line(27) = {1,13}; // 90º Rotor<BR>Line(28) = {23,3}; // 90º Estator<BR><BR>Line(30) = {12,22}; // Entrehierro<BR>Line(31) = {13,23};<BR><BR>Physical Line(33) = {1}; // Circulo externo estator<BR><BR>Physical Line(37) = {21}; // Circulo interno estator<BR><BR>Physical
 Line(41) = {11}; // Circulo rotor<BR><BR>Physical Line(45) = {26}; // 0º Estator<BR>Physical Line(46) = {28}; // 90º Estator<BR><BR>Physical Line(49) = {25}; // 0º Rotor<BR>Physical Line(50) = {27}; // 90º Rotor<BR><BR>Line Loop(53) = {28,-1,-26,21}; // 1º Cuadrante estator<BR>Plane Surface(54) = {53};<BR><BR>Line Loop(61) = {11,-27,25}; // 1º Cuadrante rotor<BR>Plane Surface(62) = {61};<BR><BR>Line Loop(69) = {21,-31,-11,30}; // Entrehierro<BR>Plane Surface(71) = {69};<BR><BR>Physical Surface(72) = {54}; // 1º Cuadrante estator<BR><BR>Physical Surface(76) = {62}; // 1º Cuadrante rotor<BR><BR>Physical Surface(80) = {71}; // Moving Band<BR><BR><BR><BR>Mesh.CharacteristicLengthFactor = 0.1;<BR><BR><BR><BR><BR>Group {<BR><BR>Stator1 = #{72};<BR><BR>Rotor1 = #{76};<BR><BR><BR>StatorRotor = #{Stator1, Rotor1} ; // 1 poles<BR><BR>S000 = #45; S090 = #46; <BR>R000 = #49; R090 = #50; <BR><BR>MB_R1 = #41; <BR>MB_R = #{MB_R1}; <BR><BR>MB_S1 = #37; <BR>MB_S = #{MB_S1};<BR><BR><BR>Outer = #{33};
 <BR><BR>MB[5000] = MovingBand2D [ #MB_R , #MB_S , 1] ;<BR><BR>Domain_Mag = #{StatorRotor,MB} ;<BR>DomainS_Mag = #{};<BR><BR>NbrPoles = 1;<BR>Rotor = #{Rotor1,MB_R};<BR><BR>Boundaries = #{S000, S090, R000, R090, MB_S, MB_R, Outer };<BR><BR>}<BR><BR><BR>/* --------------------------------------------------------------------------*/<BR><BR>Function {<BR>mu0 = 4.e-7*Pi ;<BR><BR>nu [StatorRotor] = 1. / mu0*1773.4;<BR>nu [MB] = 1. / mu0;<BR><BR>RotatePZ[] = Rotate[ Vector[$X,$Y,$Z], 0, 0, $1 ] ; <BR>RR[] = Vector[Cos[($1/180.*Pi)#1]#2 *$X- Sin[#1]#3*$Y, #3*$X+#2*$Y, $Z] ; <BR>RRrad[] = Vector[Cos[$1#1]#2 *$X- Sin[#1]#3*$Y, #3*$X+#2*$Y, $Z] ; <BR><BR>Freq = 50;<BR>Om = 2*Pi*Freq; // mech angular speed (rad/s)<BR>Period = 1./Freq;<BR><BR>NbrStep = 10;<BR><BR>theta[] = $1*Om;<BR>}<BR><BR>/* --------------------------------------------------------------------------*/<BR><BR>Constraint {<BR>{ Name MagneticVectorPotential_2D ;<BR>Case {<BR>{ Region Outer ; Value Cos[2*Atan2[$Y,$X]];
 TimeFunction F_Cos_wt_p[]{314,0}; }<BR>}<BR>}<BR>}<BR><BR><BR>Jacobian {<BR>{ Name Vol ;<BR>Case { { Region All ; Jacobian Vol ; } }<BR>}<BR>}<BR><BR><BR>Integration {<BR><BR>{ Name CurlCurl ;<BR>Case { {Type Gauss ;<BR>Case { { GeoElement Triangle ; NumberOfPoints 4 ; }<BR>{ GeoElement Quadrangle ; NumberOfPoints 4 ; }<BR>{ GeoElement Tetrahedron ; NumberOfPoints 4 ; }<BR>{ GeoElement Hexahedron ; NumberOfPoints 6 ; }<BR>{ GeoElement Line ; NumberOfPoints 2 ; } <BR>}<BR>}<BR>}<BR>}<BR><BR>}<BR><BR><BR><BR>FunctionSpace {<BR>{ Name Hcurl_a_Mag_2D ; Type Form1P ;<BR>BasisFunction {<BR>{ Name se ; NameOfCoef ae ; Function BF_PerpendicularEdge ;<BR>Support #{Domain_Mag} ; Entity NodesOf[ All ] ; }<BR>}<BR>Constraint {<BR>{ NameOfCoef ae ;<BR>EntityType NodesOf ; NameOfConstraint MagneticVectorPotential_2D ; }<BR>}<BR>}<BR>}<BR><BR><BR>Formulation {<BR>{ Name Magnetostatics_a_2D ; Type FemEquation ;<BR>Quantity {<BR>{ Name a ; Type Local ; NameOfSpace Hcurl_a_Mag_2D ; }<BR>}<BR>Equation
 {<BR>Galerkin { [ nu[] * Dof{d a} , {d a} ] ; In #Domain_Mag ;<BR>Jacobian Vol ; Integration CurlCurl ; }<BR><BR>/* Galerkin { [ 0* Dof{d a} , {d a} ] ; In #{MB_R,MB_S} ;<BR>Jacobian Vol ; Integration CurlCurl ; }<BR><BR>Galerkin { [ -js[] , {a} ] ; In DomainS_Mag ;<BR>Jacobian Vol ; Integration CurlCurl ; }<BR>*/<BR>}<BR>}<BR>}<BR><BR><BR>Resolution {<BR>{ Name MagSta_a_2D_MB_time ;<BR>System {<BR>{ Name Sys_Mag ; NameOfFormulation Magnetostatics_a_2D ; }<BR>}<BR>Operation {<BR><BR>InitSolution Sys_Mag ;<BR>Init_MovingBand2D{MB} ;<BR>Mesh_MovingBand2D { MB} ; <BR><BR>TimeLoopTheta [Period/NbrStep,Period*0.9999,Period/NbrStep, 1]<BR>{<BR>Generate Sys_Mag ;<BR>Solve Sys_Mag ;<BR>SaveSolution Sys_Mag ;<BR>ChangeOfCoordinates [NodesOf[Rotor], RRrad[$DTime*Om]]; <BR>Mesh_MovingBand2D { MB} ;<BR>PostOperation [ Map_a ];<BR>}<BR>}<BR>}<BR>}<BR><BR><BR>PostProcessing {<BR>{ Name MagSta_a_2D ; NameOfFormulation Magnetostatics_a_2D ;<BR>PostQuantity{<BR>{ Name a ; Value { Term { [ {a} ] ; In
 Domain_Mag ; } } }<BR>{ Name az ; Value { Term { [ CompZ[{a}] ] ; In Domain_Mag ; } } }<BR>{ Name b ; Value { Term { [ {d a} ]; In Domain_Mag; Jacobian Vol; } } }<BR>{ Name h ; Value { Term { [ nu[] * {d a} ]; In Domain_Mag; Jacobian Vol ; } } }<BR>{ Name bm ; Value { Term { [ Norm[{d a}] ] ; In Domain_Mag ; Jacobian Vol ; } } }<BR>{ Name one; Value { Term { [ 1 ] ; In Boundaries ; Jacobian Vol ; } } }<BR>}<BR>}<BR>}<BR><BR><BR><BR>PostOperation Map_a UsingPost MagSta_a_2D {<BR>Print[ az, OnElementsOf Domain_Mag, File "a"] ;<BR>//Print[ one, OnElementsOf Boundaries, File "bound"] ;<BR>}<BR><BR><BR><BR><BR><BR><BR><BR><BR><BR><BR><BR><BR><BR><BR><BR></BLOCKQUOTE><p>
                <hr size=1><br><font face="Verdana" size="-2">Correo Yahoo!<br>Comprueba qué es nuevo, <a href="http://us.rd.yahoo.com/mail/es/whatsnew/*http://es.whatsnew.mail.yahoo.com/">aquí</a><br>http://correo.yahoo.es</font>