[Getdp] Nonlinear BH function

janvrot at infolink.com.br janvrot at infolink.com.br
Fri May 17 00:42:46 CEST 2002


Thanks, Christophe

I was out and also had a computer hardware fault, so I could not test it
before.
It works the way you suggested.

Thank you again

Janvrot

Christophe Geuzaine wrote:
> 
> janvrot at infolink.com.br wrote:
> >
> > Hi, All
> >
> > I'm using a function to express a nonlinear BH curve.
> >
> > So, I have mu as a function of H in a Magnetic Static Scalar Potential
> > formulation. Is that ok, or it shoul be the other way? Here is the code,
> 
> Yes, that's the way to do it. Are you sure that Domain, Steel and Magnet
> do not overlap? It would be better to write:
> 
> Domain = Region[{NonSteel,Steel}];
> mu[NonSteel] = mu0;
> mu[Steel] = 1./(as + bs * Norm[$1]);
> 
> and simply
> 
> Galerkin { [ - mu[{d phi}] * Dof{d phi} , {d phi} ];
>              In Domain; Jacobian MyJac; Integration MyInt; }
> 
> > and the error message.
> >
> > Function {
> >   mu0 = 4.e-7 * Pi ;
> >   murMagnet = 1.1;
> >   as = 417.497;
> >   bs = 38.3881;
> >   mu[Steel] = 1./(as + bs * Norm[$1]);
> >
> >   mu[Magnet]  = murMagnet * mu0;
> >   mu[Air]     = mu0;
> >   hc[ Magnet1 ] = Vector[0., 0.,  10000.];
> >   hc[ Magnet2 ] = Vector[0., 0., -10000.];
> > }
> >
> > Formulation {
> >   { Name MagSta_phi; Type FemEquation;
> >     Quantity {
> >       { Name phi; Type Local; NameOfSpace Hgrad_phi; }
> >     }
> >     Equation {
> >       Galerkin { [ - mu[] * Dof{d phi} , {d phi} ];
> >                  In Domain; Jacobian MyJac; Integration MyInt; }
> >       Galerkin { [ - mu[{d phi}] * Dof{d phi} , {d phi} ];
> >                   In Steel; Jacobian MyJac; Integration MyInt; }
> >       Galerkin { [ - mu[] * hc[] , {d phi} ];
> >                  In Magnet; Jacobian MyJac; Integration MyInt; }
> >     }
> >   }
> > }
> >
> > I get the message:
> > P r o c e s s i n g . . .
> > Operation : Generate[A]
> > Solver    : Loading parameter file 'solver.par'
> > Error     : Wrong type of argument for function 'Norm'
> >
> > Thanks for any help.
> > --
> > Janvrot IVM
> > janvrot at infolink.com.br
> > _______________________________________________
> > GetDP mailing list
> > GetDP at geuz.org
> > http://www.geuz.org/mailman/listinfo/getdp
> 
> --
> Christophe Geuzaine
> 
> Tel: (626) 395-4552    http://www.geuz.org
> Fax: (626) 578-0124    mailto:geuzaine at acm.caltech.edu
> _______________________________________________
> GetDP mailing list
> GetDP at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp

-- 
Janvrot IVM
janvrot at infolink.com.br