[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