[Getdp] Nonlinear BH function

janvrot at infolink.com.br janvrot at infolink.com.br
Fri May 10 05:21:03 CEST 2002


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,
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