<html>

<head>

<meta http-equiv="Content-Type" content="text/html; charset=utf-8">

</head>

<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">

<br class="">

<div>

<blockquote type="cite" class="">

<div class="">On 13 Apr 2017, at 09:38, Pablo Xavier Jara Palacios <<a href="mailto:pabloxjara@gmail.com" class="">pabloxjara@gmail.com</a>> wrote:</div>

<br class="Apple-interchange-newline">

<div class="">

<div dir="ltr" class="">

<p class="MsoNormal" style="text-align:start">Thanks Ruth, I really appreciate your help. Do you mean that un is used for the purpose of compute the normal direction of the force on each magnet?.

</p>

</div>

</div>

</blockquote>

<div>Not on each magnet but on each element that touches the boundary of the magnet.</div>

<blockquote type="cite" class="">

<div class="">

<div dir="ltr" class="">

<p class="MsoNormal" style="text-align:start">On the other hand, in the weak formulation this part  "Galerkin { [ 0 * Dof{un~{i}} , {un~{i}} ]" is a separetely equation or is part of the same  equation of the weak formulation either scalar o vector potential?</p>

</div>

</div>

</blockquote>

<div>You can see it as a dummy equation, as it has no influence on the system you are actually solving. </div>

<div>It is an additional equation with the only objective of generating the degrees of freedom (Dofs), they are set to zero unless indicated otherwise via a constraint.</div>

<div><br class="">

</div>

<blockquote type="cite" class="">

<div class="">

<div dir="ltr" class="">

<p class="MsoNormal" style="text-align:start">The thing is that I do not undestand how un is computed if there is a 0 that is multipliying inside the integranl, so I would suppose that everything is canceled out in the integral.

</p>

</div>

</div>

</blockquote>

<div>As I said, it is not computed. All the coefficients associated with that basis function are set to zero except the ones fixed via the constraint => the boundary is set to one.</div>

<br class="">

<blockquote type="cite" class="">

<div class="">

<div dir="ltr" class="">

<p class="MsoNormal" style="text-align:start">In this point is where I do not know how un is calculated at the end. Please maybe could you recommend me any book or paper where I could find a similar example?</p>

</div>

</div>

</blockquote>

I do not think there is an example of this as is. The 0 multiplying an equation, it is just a GetDP trick for getting the local normal.</div>

<div><br class="">

</div>

<div>The papers you can have a look at are those dealing with the Maxwell Stress Tensor, also called Arkkio method in the electrical-machine modelling community.</div>

<span class="">

<div class=""><span class=""><br class="">

</span></div>

<div class=""><span class="">Have a look at:</span></div>

A. Arkkio, Analysis of induction motors based on the numerical solution of the magnetic field and circuit equations , Acta Polytechnica Scandinavica, 1987, p. 56<br class="">

Z. Ren, A. Razek, Local force computation in deformable bodies using edge elements , IEEE Trans. Magnetics, 28(2):1212–1215, 1992.<br class="">

F. Henrotte, G. Deliége, K. Hameyer, The eggshell approach for the computations of electromagnetic forces in 2D and 3D , COMPEL, 23(4), 996–1005, 2004.</span>

<div class=""><br class="">

</div>

<div class=""><span class="">Regards,</span></div>

<div class=""><span class="">Ruth<br class="">

</span><span class=""><br class="">

</span><span class=""><br class="">

</span><span class=""><br class="">

<blockquote type="cite" class=""><br class="">

regards,<br class="">

<br class="">

<br class="">

<br class="">

<br class="">

2017-04-12 7:43 GMT-07:00 Ruth Vazquez Sabariego <<a href="mailto:ruth.sabariego@kuleuven.be" class="">ruth.sabariego@kuleuven.be</a>>:<br class="">

Dear Pablo Xavier, <br class="">

<br class="">

The force is calculated by means of the Maxwell Stress Tensor with the local normal computed with the help of an auxiliary basis function. <br class="">

The dummy term added in the formulation is then only used in the post-processing stage for the force computation.<br class="">

<br class="">

The unknowns linked to this basis function un are indeed fully fixed and linked to a layer of elements touching the magnet. <br class="">

We impose 1 at the boundary of the magnet and the rest of the values are  0, as the function varies linearly per element, the gradient gives the normal modulus the size of the element. <br class="">

<br class="">

HTH,<br class="">

Ruth<br class="">

<br class="">

<br class="">

—<br class="">

Prof. Ruth V. Sabariego<br class="">

KU Leuven  <br class="">

Dept. Electrical Engineering ESAT/Electa, EnergyVille<br class="">

http://www.esat.kuleuven.be/electa<br class="">

http://www.energyville.be<br class="">

<br class="">

Free software: http://gmsh.info | http://getdp.info | http://onelab.info<br class="">

<br class="">

<br class="">

<br class="">

<br class="">

<br class="">

<br class="">

<br class="">

<blockquote type="cite" class="">On 12 Apr 2017, at 08:52, Pablo Xavier Jara Palacios <pabloxjara@gmail.com> wrote:<br class="">

<br class="">

Hi,<br class="">

I was reviewing the .pro file in example of magnet which is posted in onelab's page, but I do not understand how the force is calculated in this problem. There is a dummy term in the formulation that I do not know what it is. Please help me with this doubt.<br class="">

<br class="">

Regards,<br class="">

<br class="">

<br class="">

_______________________________________________<br class="">

getdp mailing list<br class="">

getdp@onelab.info<br class="">

http://onelab.info/mailman/listinfo/getdp<br class="">

</blockquote>

<br class="">

<br class="">

</blockquote>

</span><br class="">

</div>

</body>

</html>