[Getdp] Nonlinear BH function
Christophe Geuzaine
geuzaine at acm.caltech.edu
Thu May 9 19:02:39 CEST 2002
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