[Getdp] RV: Moving Band

Wed Nov 9 22:16:45 CET 2005

Dear Christophe:

I sent this mail to Johan Gynselink, but I do not have response.

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. 

Until now I have utilized the function MB for all the motor and I have not problems, but requires larger calculation time.

I wait for your comments.

Best regard...

Christian Geikowsky R.

"Christian Geikowsky R." <christiangeikowsky at yahoo.es> escribió:
Fecha: Fri, 4 Nov 2005 14:35:58 +0100 (CET)
De: "Christian Geikowsky R." <christiangeikowsky at yahoo.es>
Asunto: Moving Band
Para: Johan Gyselinck <johan.gyselinck at ulb.ac.be>

Estimated Johan:

Attachment a simple geometry of a fourth part of a 4 poles motor.

I would like know

How must I apply the boundary conditions in the limits S000, S090, R000 and R090 and in the air gap? 

How must give form to the function MovingBand for a fourth of the geometry?

Best regard...

Saludos...

Christian Geikowsky R.

---------------------------------

Correo Yahoo!
Comprueba qué es nuevo, aquí
http://correo.yahoo.esR1 = 0.5; // Radio rotor
R2 = 0.6; // Radio interno Estator

Point(1) = {0, 0, 0, 1}; 

Point(2) = {1, 0, 0, 1}; // Circulo externo estator
Point(3) = {0, 1, 0, 1};

Circle(1) = {2,1,3};

Point(12) = {R1, 0, 0, 1}; // Circulo rotor
Point(13) = {0, R1, 0, 1};

Circle(11) = {12,1,13};

Point(22) = {R2, 0, 0, 1}; // Circulo interno estator
Point(23) = {0, R2, 0, 1};

Circle(21) = {22,1,23};

Line(25) = {1,12}; // 0º Rotor
Line(26) = {22,2}; // 0º Estator
Line(27) = {1,13}; // 90º Rotor
Line(28) = {23,3}; // 90º Estator

Line(30) = {12,22}; // Entrehierro
Line(31) = {13,23};

Physical Line(33) = {1}; // Circulo externo estator

Physical Line(37) = {21}; // Circulo interno estator

Physical Line(41) = {11}; // Circulo rotor

Physical Line(45) = {26}; // 0º Estator
Physical Line(46) = {28}; // 90º Estator

Physical Line(49) = {25}; // 0º Rotor
Physical Line(50) = {27}; // 90º Rotor

Line Loop(53) = {28,-1,-26,21}; // 1º Cuadrante estator
Plane Surface(54) = {53};

Line Loop(61) = {11,-27,25}; // 1º Cuadrante rotor
Plane Surface(62) = {61};

Line Loop(69) = {21,-31,-11,30}; // Entrehierro
Plane Surface(71) = {69};

Physical Surface(72) = {54}; // 1º Cuadrante estator

Physical Surface(76) = {62}; // 1º Cuadrante rotor

Physical Surface(80) = {71}; // Moving Band

Mesh.CharacteristicLengthFactor = 0.1;

Group {

Stator1 = #{72};

Rotor1 = #{76};

StatorRotor = #{Stator1, Rotor1} ; // 1 poles

S000 = #45; S090 = #46; 
R000 = #49; R090 = #50; 

MB_R1 = #41; 
MB_R = #{MB_R1}; 

MB_S1 = #37; 
MB_S = #{MB_S1};

Outer = #{33}; 

MB[5000] = MovingBand2D [ #MB_R , #MB_S , 1] ;

Domain_Mag = #{StatorRotor,MB} ;
DomainS_Mag = #{};

NbrPoles = 1;
Rotor = #{Rotor1,MB_R};

Boundaries = #{S000, S090, R000, R090, MB_S, MB_R, Outer };

}

/* --------------------------------------------------------------------------*/

Function {
mu0 = 4.e-7*Pi ;

nu [StatorRotor] = 1. / mu0*1773.4;
nu [MB] = 1. / mu0;

RotatePZ[] = Rotate[ Vector[$X,$Y,$Z], 0, 0, $1 ] ; 
RR[] = Vector[Cos[($1/180.*Pi)#1]#2 *$X- Sin[#1]#3*$Y, #3*$X+#2*$Y, $Z] ; 
RRrad[] = Vector[Cos[$1#1]#2 *$X- Sin[#1]#3*$Y, #3*$X+#2*$Y, $Z] ; 

Freq = 50;
Om = 2*Pi*Freq; // mech angular speed (rad/s)
Period = 1./Freq;

NbrStep = 10;

theta[] = $1*Om;
}

/* --------------------------------------------------------------------------*/

Constraint {
{ Name MagneticVectorPotential_2D ;
Case {
{ Region Outer ; Value Cos[2*Atan2[$Y,$X]]; TimeFunction F_Cos_wt_p[]{314,0}; }
}
}
}

Jacobian {
{ Name Vol ;
Case { { Region All ; Jacobian Vol ; } }
}
}

Integration {

{ Name CurlCurl ;
Case { {Type Gauss ;
Case { { GeoElement Triangle ; NumberOfPoints 4 ; }
{ GeoElement Quadrangle ; NumberOfPoints 4 ; }
{ GeoElement Tetrahedron ; NumberOfPoints 4 ; }
{ GeoElement Hexahedron ; NumberOfPoints 6 ; }
{ GeoElement Line ; NumberOfPoints 2 ; } 
}
}
}
}

}

FunctionSpace {
{ Name Hcurl_a_Mag_2D ; Type Form1P ;
BasisFunction {
{ Name se ; NameOfCoef ae ; Function BF_PerpendicularEdge ;
Support #{Domain_Mag} ; Entity NodesOf[ All ] ; }
}
Constraint {
{ NameOfCoef ae ;
EntityType NodesOf ; NameOfConstraint MagneticVectorPotential_2D ; }
}
}
}

Formulation {
{ Name Magnetostatics_a_2D ; Type FemEquation ;
Quantity {
{ Name a ; Type Local ; NameOfSpace Hcurl_a_Mag_2D ; }
}
Equation {
Galerkin { [ nu[] * Dof{d a} , {d a} ] ; In #Domain_Mag ;
Jacobian Vol ; Integration CurlCurl ; }

/* Galerkin { [ 0* Dof{d a} , {d a} ] ; In #{MB_R,MB_S} ;
Jacobian Vol ; Integration CurlCurl ; }

Galerkin { [ -js[] , {a} ] ; In DomainS_Mag ;
Jacobian Vol ; Integration CurlCurl ; }
*/
}
}
}

Resolution {
{ Name MagSta_a_2D_MB_time ;
System {
{ Name Sys_Mag ; NameOfFormulation Magnetostatics_a_2D ; }
}
Operation {

InitSolution Sys_Mag ;
Init_MovingBand2D{MB} ;
Mesh_MovingBand2D { MB} ; 

TimeLoopTheta [Period/NbrStep,Period*0.9999,Period/NbrStep, 1]
{
Generate Sys_Mag ;
Solve Sys_Mag ;
SaveSolution Sys_Mag ;
ChangeOfCoordinates [NodesOf[Rotor], RRrad[$DTime*Om]]; 
Mesh_MovingBand2D { MB} ;
PostOperation [ Map_a ];
}
}
}
}

PostProcessing {
{ Name MagSta_a_2D ; NameOfFormulation Magnetostatics_a_2D ;
PostQuantity{
{ Name a ; Value { Term { [ {a} ] ; In Domain_Mag ; } } }
{ Name az ; Value { Term { [ CompZ[{a}] ] ; In Domain_Mag ; } } }
{ Name b ; Value { Term { [ {d a} ]; In Domain_Mag; Jacobian Vol; } } }
{ Name h ; Value { Term { [ nu[] * {d a} ]; In Domain_Mag; Jacobian Vol ; } } }
{ Name bm ; Value { Term { [ Norm[{d a}] ] ; In Domain_Mag ; Jacobian Vol ; } } }
{ Name one; Value { Term { [ 1 ] ; In Boundaries ; Jacobian Vol ; } } }
}
}
}

PostOperation Map_a UsingPost MagSta_a_2D {
Print[ az, OnElementsOf Domain_Mag, File "a"] ;
//Print[ one, OnElementsOf Boundaries, File "bound"] ;
}

---------------------------------

Correo Yahoo!
Comprueba qué es nuevo, aquí
http://correo.yahoo.es
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20051109/3ee07727/attachment.html>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: circle1.geo
URL: <http://www.geuz.org/pipermail/getdp/attachments/20051109/3ee07727/attachment.geo>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: circle1.pro
URL: <http://www.geuz.org/pipermail/getdp/attachments/20051109/3ee07727/attachment.pro>