Electromagnetism

In CalculiX, certain types of electromagnetic calculations are possible. These include:

- electrostatic calculations: these can be performed as a special case of thermal calculations, cf. Section 6.9.13.
- magnetostatic calculations. Due to the absence of time derivatives the interaction between electric and magnetic fields drops out and the magnetic equations can be considered on their own.
- magnetic induction calculations. These are calculations without fixed electric charges and in which the displacement current can be neglected. The major industrial application for this type of calculation is inductive heating.

In this section only the last two applications are treated. The governing Maxwell equations run like (the displacement current term was dropped in Equation (510)):

where is the electric field, is the electric displacement field, is the magnetic field, is the magnetic intensity, is the electric current density and is the electric charge density. These fields are connected by the following constitutive equations:

and

Here, is the permittivity, is the magnetic permeability and is the electrical conductivity. For the present applications and are not needed and Equation (507) can be discarded. It will be assumed that these relationships are linear and isotropic, the material parameters, however, can be temperature dependent. So no hysteresis is considered, which basically means that only paramagnetic and diamagnetic materials are considered. So far, no ferromagnetic materials are allowed.

Due to electromagnetism, an additional basic unit is needed, the Ampère (A). All other quantities can be written using the SI-units A, m, s, kg and K, however, frequently derived units are used. An overview of these units is given in Table 17 (V=Volt, C=Coulomb, T=Tesla, F=Farad, S=Siemens).

symbol | meaning | unit |

I | current | A |

electric field | ||

electric displacement field | ||

magnetic field | ||

magnetic intensity | ||

current density | ||

permittivity | ||

magnetic permeability | ||

electrical conductivity | ||

P | magnetic scalar potential | A |

V | electric scalar potential |
V |

magnetic vector potential | ||

In what follows the references [77] and [38] have been used. In inductive heating applications the domain of interest consists of the objects to be heated (= workpiece), the surrounding air and the coils providing the current leading to the induction. It will be assumed that the coils can be considered seperately as a driving force without feedback from the system. This requires the coils to be equiped with a regulating system counteracting any external influence trying to modify the current as intended by the user.

Let us first try to understand what happens physically. In the simplest case the volume to be analyzed consists of a simply connected body surrounded by air, Figure 139. A body is simply connected if any fictitious closed loop within the body can be reduced to a point without leaving the body. For instance, a sphere is simply connected, a ring is not. The coil providing the current is located within the air. Turning on the current leads to a magnetic intensity field through Equation (510) and a magnetic field through Equation (512) everywhere, in the air and in the body. If the current is not changing in time, this constitutes the solution to the problem.

If the current is changing in time, so is the magnetic field, and through Equation (508) one obtains an electric field everywhere. This electric field generates a current by Equation (513) (called Eddy current) in any part which is electrically conductive, i.e. generally in the body, but not in the air. This current generates a magnetic intensity field by virtue of Equation (510), in a direction which is opposite to the original magnetic intensity field. Thus, the Eddy currents oppose the generation of the magnetic field in the body. Practically, this means that the magnetic field in the body is not built up at once. Rather, it is built up gradually, in the same way in which the temperature in a body due to heat transfer can only change gradually. As a matter of fact, both phenomena are described by first order differential equations in time. The Ohm-losses of the Eddy currents are the source of the heat generation used in industrial heat induction applications.

From these considerations one realizes that in the body (domain 2, cf. Figure 139; notice that domain 1 and 2 are interchanged compared to [77]) both the electric and the magnetic field have to be calculated, while in the air it is sufficient to consider the magnetic field only (domain 1). Therefore, in the air it is sufficient to use a scalar magnetic potential P satisfying:

Here, is the magnetic intensity due to the coil current in infinite free space. can be calculated using the Biot-Savart relationships [23]. The body fields can be described using a vector magnetic potential and a scalar electric potential V satisfying:

In practice, it is convenient to set , leading to

(517) |

This guarantees that the resulting matrices will be symmetric.

If the body is multiply connected, the calculational domain consists of three domains. The body (or bodies) still consist of domain 2 governed by the unknowns and . The air, however, has to split into two parts: one part which is such that, if added to the bodies, makes them simply connected. This is domain 3 and it is described by the vector magnetic potential . It is assumed that there are no current conducting coils in domain 3. The remaining air is domain 1 described by the scalar magnetic potential P.

In the different domains, different equations have to be solved. In domain 1 the electric field is not important, since there is no conductance. Therefore, it is sufficient to calculate the magnetic field, and only Equations (509) and (510) have to be satisfied. Using the ansatz in Equation (514), Equation (510) is automatically satiesfied, since it is satisfied by and the curl of the gradient vanishes. The only equation left is (509). One arrives at the equation

In domain 2, Equations (508), (509) and (510) have to be satisfied, using the approach of Equations (515) and (516). Taking the curl of Equation (516) yields Equation (508). Taking the divergence of Equation (515) yields Equation (509). Substituting Equations (515) and (516) into Equation (510) leads to:

The magnetic vector potential is not uniquely defined by Equation (515). The divergence of can still be freely defined. Here, we take the Coulomb gauge, which amounts to setting

Notice that the fulfillment of Equation (510) automatically satisfies the conservation of charge, which runs in domain 2 as

since there is no concentrated charge. Thus, for a simply connected body we arrive at the Equations (518) (domain 1), (519) (domain 2) and (520) (domain 2). In practice, Equations (519) and (520) are frequently combined to yield

This, however, is not any more equivalent to the solution of Equation (510) and consequently the satisfaction of Equation (521) has now to be requested explicitly:

Consequently, the equations to be solved are now Equations (522) (domain 2), (523) (domain 2), and (518) (domain 1).

In domain 3, only Equations (509) and (510) with have to be satisfied (the coils are supposed to be in domain 1). Using the ansatz from Equation (515), Equation (509) is automatically satisfied and Equation (510) now amounts to

(524) |

The boundary conditions on the interface amount to:

- continuity of the normal component of
- continuity of the tangential component of , and
- no current flow orthogonal to the boundary, or ( is the normal on domain i, pointing away from the domain):

(525) |

(526) |

and

(527) |

all of which have to be satisfied on . In terms of the magnetic vector potential , electric scalar potential V and magnetic scalar potential P this amounts to:

and

on . For uniqueness, the electric potential has to be fixed in one node and the normal component of has to vanish along [77]:

(531) |

To obtain the weak formulation of the above equations they are multiplied with trial functions and integrated. The trial functions will be denoted by and . Starting with Equation(522) one obtains after multiplication with and taking the vector identies

(532) |

(533) |

into account (set in the first vector identity):

(534) |

Integrating one obtains, using Gauss' theorem (it is assumed that has no free boundary, i.e. no boundary not connected to ):

(535) |

The trial functions also have to satisfy the kinematic constraints. Therefore, and the second surface integral is zero.

Applying the vector identity

(536) |

and the boundary condition from Equation (529), the integrand of the first surface integral can be written as:

(537) |

Consequently, the integral now amounts to:

(538) |

Applying the same vector identity from above one further arrives at:

(539) |

Finally, using the vector identity:

(540) |

one obtains

(541) |

The last integral vanishes if the surface is closed due to Stokes' Theorem.

Now the second equation, Equation (523), is being looked at. After multiplication with it can be rewritten as:

(542) |

After integration and application of Gauss' theorem one ends up with the last term only, due to the boundary condition from Equation (530).

Analogously, the third equation, Equation (518) leads to:

(543) |

After integration this leads to (on external faces of , i.e. faces not connected to or the condition is applied) :

(544) |

Applying the boundary condition from Equation (528) leads to:

(545) |

So one finally obtains for the governing equations :

(546) |

(547) |

(548) |

Using the standard shape functions one arrives at (cf. Chapter 2 in [19]):

(549) |

(550) |

(551) | ||

(552) | ||

(553) |

Notice that the first two equations apply to domain 2, the last one applies to domain 1. In domain 3 only the first equation applies, in which the time dependent terms are dropped.

This leads to the following matrices:

(554) |

(555) |

(556) |

(557) |

(558) |

(559) |

(560) |

(561) |

(563) |

Repeated indices imply implicit summation. The matrices are analogous to the conductivity matrix in heat transfer analyses, the matrices are the counterpart of the capacity matrix. represents the force. The resulting system consists of first order ordinary differential equations in time and the corresponding matrices look like in Figure 141. Solution of this system yields the solution for and from which the magnetic field and the electric field can be determined using Equations (515,516).

The internal electromagnetic forces amount to:

(564) |

and

(565) |

They have to be in equilibrium with the external forces.

What does the above theory imply for the practical modeling? The conductor containing the driving current is supposed to be modeled using shell elements. The thickness of the shell elements can vary. The current usually flows near the surface (skin effect), so the modeling with shell elements is not really a restriction. The current and potential in the conductor is calculated using the heat transfer analogy. This means that potential boundary conditions have to be defined as temperature, current boundary conditions as heat flow conditions. The driving current containing conductor is completely separate from the mesh used to calculate the magnetic and electric fields. Notice that the current in the driving electromagnetic coils is not supposed to be changed by the electromagnetic field it generates.

The volumetric domains of interest are and . These three domains represent the air, the conducting workpiece and that part of the air which, if filled with workpiece material, makes the workpiece simply connected. These three domains are to be meshed with volumetric elements. The meshes should not be connected, i.e., one can mesh these domains in a completely independent way. This also applies that one can choose the appropriate mesh density for each domain separately.

Based on the driving current the field is determined in domain 1 with the Biot-Savart law. This part of the code is parallellized, since the Biot-Savart integration is calculationally quite expensive. Because of Equation (562) is also determined on the external faces of domain 2 and 3 which are in contact with domain 1.

The following boundary conditions are imposed (through MPC's):

On the external faces of domain 1 which are in contact with domain 2 or domain 3:

- Calculation of based on the external facial values in domain 2 and 3 (cf. the area integrals in )

On the external faces of domain 2 and 3 which are in contact with domain 1:

- Calculation of based on the external facial values in domain 1.
- Imposition of .

On the faces between domain 2 and 3:

- Continuity of

These MPC's are generated automatically within CalculiX and have not to be taken care of by the user. Finally, the value of V has to be fixed in at least one node of domain 2. This has to be done by the user with a *BOUNDARY condition on degree of freedom 8.

The material data to be defined include:

- the electrical conductivity in the driving coils
- the magnetic permeability in the air
- the density, the thermal conductivity, the specific heat, the electrical conductivity and the magnetic permeability in the workpiece.

To this end the cards *DENSITY, *CONDUCTIVITY, *SPECIFIC HEAT, *ELECTRICAL CONDUCTIVITY, MAGNETIC PERMEABILITY can be used. In the presence of thermal radiation the *PHYSICAL CONSTANTS card is also needed.

The procedure card is *ELECTROMAGNETICS. For magnetostatic calculations the parameter MAGNETOSTATICS is to be used, for athermal electromagnetic calculations the parameter NO HEAT TRANSFER. Default is an electromagnetic calculation with heat transfer.

Available output variables are POT (the electric potential in the driving current coil) on the *NODE FILE card and ECD (electric current density in the driving current coil), EMFE (electric field in the workpiece) and EMFB (magnetic field in the air and the workpiece) on the *EL FILE card. Examples are induction.inp and induction2.inp.