<!doctype html public "-//w3c//dtd html 4.0 transitional//en">
<html>
janvrot@infolink.com.br wrote:
<blockquote TYPE=CITE>> You just have to replace dhdb[] and h[] by dbdh[]
and b[]
<br>> respectively. This was already partly prepared in the file with the
<br>> tables I sent you a while ago.
<br>>
<p>Sorry, but I didn't get it.</blockquote>
I was referring to this:
<p> Jan_mu_h2 = ListAlt[Jan_h2, Jan_mu] ;
<br> mu_Jan[] = InterpolationAkima[$1]{List[Jan_mu_h2]} ;
<br> dmudh2_Jan[] = dInterpolationAkima[$1]{List[Jan_mu_h2]} ;
<br> b_Jan[] = mu_Jan[(SquNorm[$1#3])] * #3 ;
<br> dbdh_Jan[] = TensorDiag[1,1,1] * mu_Jan[SquNorm[$1#1]#2] + 2*dmudh2_Jan[#2]
* SquDyadicProduct[#1] ;
<p>which was in the file I sent you, wasn't it?
<blockquote TYPE=CITE>
<br>So, is this correct for the scalar magnetic potential in a case where
<br>there are no current loops?
<p>Formulation {
<br> { Name MagSta_phi_nl; Type FemEquation;
<br> Quantity {
<br> { Name phi; Type Local; NameOfSpace
Hgrad_phi; }
<br> }
<br> Equation {
<br> Galerkin { [ - mu[] * Dof{d phi} , {d
phi} ];
<br>
In NonSteel; Jacobian MyJac; Integration MyInt; }
<p> Galerkin { JacNL[ dbdh[{d phi}] * Dof{d
phi} , {d phi} ] ;
<br>
In Steel ; Jacobian MyJac ; Integration MyInt ; }
<p> Galerkin { [ b[{d phi}] , {d phi} ]
;
<br>
In Steel ; Jacobian MyJac ; Integration MyInt ; }
<p> Galerkin { [ - mu[] * hc[] , {d phi}
];
<br>
In Magnet; Jacobian MyJac; Integration MyInt; }
<br> }
<br> }
<br>}</blockquote>
You should add a minus sign to the two terms in the middle! (Remember
that in the linear case, dbdh = mu !)
<br>Besides that it looks ok.
<br>
<p>Johan
<p>PS For those trying to solve a 3D magnetostatic problem with an imposed
current density using GetDP, maybe you should consider the possibility
of modelling the excitation (in an approximate way) by means of one or
several permanent magnets. In that case, the formulation is very
simple.
<pre>--
Dr. ir. Johan Gyselinck
University of Liege
Dept. of Electrical Engineering (ELAP), Institut Montefiore
Campus du Sart Tilman, Parking P32, Bâtiment B28, B-4000 LIEGE, Belgique
phone: +32-4-3663732 -- fax: +32-4-3662910 -- <A HREF="http://elap.montefiore.ulg.ac.be/">http://elap.montefiore.ulg.ac.be/</A></pre>
</html>