The scalar potential method as implemented in SOLID5, SOLID96, and SOLID98 for 3D magnetostatic fields is discussed in this section. Magnetostatics means that time varying effects are ignored. This reduces Maxwell's equations for magnetic fields to:

(5–14)

(5–15)

In the domain Ω₀ and Ω₁ of a magnetostatic field problem (Ω₂ is not considered for magnetostatics) a solution is sought which satisfies the relevant Maxwell's Equation 5–14 and Equation 5–15 and the constitutive relation Equation 5–10 in the following form (Gyimesi([142]) and Gyimesi([150])):

(5–16)

(5–17)

where:

The development of {H_g} varies depending on the problem and the formulation. Basically, {H_g} must satisfy Ampere's law (Equation 5–14) so that the remaining part of the field can be derived as the gradient of the generalized scalar potential φ_g. This ensures that φ_g is singly valued. Additionally, the absolute value of {H_g} must be greater than that of Δφ_g. In other words, {Hg} should be a good approximation of the total field. This avoids difficulties with cancellation errors (Gyimesi([150])).

This framework allows for a variety of scalar potential formulation to be used. The appropriate formulation depends on the characteristics of the problem to be solved. The process of obtaining a final solution may involve several steps (controlled by the MAGOPT solution option).

As mentioned above, the selection of {H_g} is essential to the development of any of the following scalar potential strategies. The development of {H_g} always involves the Biot-Savart field {H_s} which satisfies Ampere's law and is a function of source current {J_s}. {H_s} is obtained by evaluating the integral:

(5–18)

where:

The above volume integral can be reduced to the following surface integral (Gyimesi et al.([175]))

(5–19)

where:

Evaluation of this integral is automatically performed upon initial solution execution or explicitly (controlled by the BIOT command). The values of {J_s} are obtained either directly as input by:

or indirectly calculated by electric field calculation using:

Depending upon the current configuration, the integral given in Equation 5–19 is evaluated in a closed form and/or a numerical fashion (Smythe([151])).

Three different solution strategies emerge from the general framework discussed above:

 

5.1.2.1. RSP Strategy

Applicability

If there are no current sources ({J_s} = 0), the RSP strategy is applicable. The RSP strategy is also applicable, in general, if there are current sources but there is no iron ([μ] = [μ_o]) within the problem domain. This formulation is developed by Zienkiewicz([76]).

Procedure

The RSP strategy uses a one-step procedure (MAGOPT,0). Equation 5–16 and Equation 5–17 are solved making the following substitution:

(5–20)

Saturation is considered if the magnetic material is nonlinear. Permanent magnets are also considered.

 

5.1.2.2. DSP Strategy

Applicability

The DSP strategy is applicable when current sources and singly connected iron regions exist within the problem domain ({J_s} ≠ {0}) and ([μ] ≠ [μ_o]). A singly connected iron region does not enclose a current. In other words, a contour integral of {H} through the iron must approach zero as u → .

(5–21)

This formulation is developed by Mayergoyz([120]).

Procedure

The DSP strategy uses a two-step solution procedure. The first step (MAGOPT,2) makes the following substitution into Equation 5–16 and Equation 5–17:

(5–22)

subject to:

(5–23)

This boundary condition is satisfied by using a very large value of permeability in the iron (internally set by the program). Saturation and permanent magnets are not considered. This step produces a near zero field in the iron region which is subsequently taken to be zero according to:

(5–24)

and in the air region:

(5–25)

The second step (MAGOPT,3) uses the fields calculated on the first step as the preliminary field for Equation 5–16 and Equation 5–17:

(5–26)

(5–27)

Here saturation and permanent magnets are considered. This step produces the following fields:

(5–28)

and

(5–29)

which are the final results to the applicable problems.

 

5.1.2.3. GSP Strategy

Applicability

The GSP strategy is applicable when current sources ({J_s ≠ {0}) exist within the problem domain in conjunction with a multiply connected iron ([μ] ≠ [μ_o]) region. A multiply connected iron region encloses some current source. This means that a contour integral of {H} through the iron region is not zero:

(5–30)

where:

= refers to the dot product

This formulation is developed by Gyimesi([142], [150], [203]).

Procedure

The GSP strategy uses a three-step solution procedure. The first step (MAGOPT,1) performs a solution only in the iron with the following substitution into Equation 5–16 and Equation 5–17:

(5–31)

subject to:

(5–32)

Here S₁ is the surface of the iron air interface. Saturation can optimally be considered for an improved approximation of the generalized field but permanent magnets are not. The resulting field is:

(5–33)

The second step (MAGOPT,2) performs a solution only in the air with the following substitution into Equation 5–16 and Equation 5–17:

(5–34)

subject to:

(5–35)

This boundary condition is satisfied by automatically constraining the potential solution φ_g at the surface of the iron to be what it was on the first step (MAGOPT,1). This step produces the following field:

(5–36)

Saturation or permanent magnets are of no consequence since this step obtains a solution only in air.

The third step (MAGOPT,3) uses the fields calculated on the first two steps as the preliminary field for Equation 5–16 and Equation 5–17:

(5–37)

(5–38)

Here saturation and permanent magnets are considered. The final step allows for the total field to be computed throughout the domain as:

(5–39)

The vector potential method, implemented in PLANE13 and PLANE233 for 2D electromagnetic fields, is discussed in this section. Considering static and dynamic fields and neglecting displacement currents (quasi-stationary limit), the following subset of Maxwell's equations apply:

(5–40)

(5–41)

(5–42)

The usual constitutive equations for magnetic and electric fields apply as described by Equation 5–11 and Equation 5–12. Although some restrictions on anisotropy and nonlinearity do occur in the formulations mentioned below.

In the entire domain, Ω, of an electromagnetic field problem a solution is sought which satisfies the relevant Maxwell's Equation 5–40 thru Equation 5–41. See Figure 5.1: Electromagnetic Field Regions for a representation of the problem domain Ω.

A solution can be obtained by introducing potentials which allow the magnetic field {B} and the electric field {E} to be expressed as (Biro([121])):

(5–43)

(5–44)

where:

These specifications ensure the satisfaction of two of Maxwell's equations, Equation 5–41 and Equation 5–42. What remains to be solved is Ampere's law, Equation 5–40 in conjunction with the constitutive relations, Equation 5–11, and the divergence free property of current density. The resulting differential equations are:

(5–45)

(5–46)

(5–47)

These equations are subject to the appropriate boundary conditions.

This system of simplified Maxwell's equations with the introduction of potential functions is used for the solutions of 2D static and dynamic fields. Silvester([73]) presents a 2D static formulation. Chari([70]), Brauer([71]) and Tandon([72]) discuss the 2D eddy current problem and Weiss([95]) and Garg([96]) discuss 2D eddy current problems which allow for skin effects (eddy currents present in the source conductor).

An edge-based method, using edge-based elements with a discontinuous normal component of magnetic vector potential, is implemented in the 3D electromagnetic SOLID236 and SOLID237 elements.

The differential equations governing SOLID236 and SOLID237 elements are the following:

(5–48)

(5–49)

(5–50)

These equations are subject to the appropriate magnetic and electrical boundary conditions.

The uniqueness of edge-based magnetic vector potential is ensured by the tree gauging procedure (GAUGE command) that sets the edge-flux degrees of freedom corresponding to the spanning tree of the finite element mesh to zero.

In a general dynamic problem, any field quantity, q(r,t) depends on the space, r, and time, t, variables. In a harmonic analysis, the time dependence can be described by periodic functions:

(5–51)

or

(5–52)

where:

In an electromagnetic analysis, q(r,t) can be the flux density, {B}, the magnetic field, {H}, the electric field, {E}, the current density, J, the vector potential, {A}, or the scalar potential, V. Note, however, that q(r,t) can not be the Joule heat, Q^j, the magnetic energy, W, or the force, F^jb, because they include a time-constant term.

The quantities in Equation 5–51 and Equation 5–52 are related by

(5–53)

(5–54)

(5–55)

(5–56)

In Equation 5–51) a(r), φ(r), c(r) and s(r) depend on space coordinates but not on time. This separation of space and time is taken advantage of to minimize the computational cost. The originally 4 (3 space + 1 time) dimensional real problem can be reduced to a 3 (space) dimensional complex problem. This can be achieved by the complex formalism.

The measurable quantity, q(r,t), is described as the real part of a complex function:

(5–57)

Q(r) is defined as:

(5–58)

where:

The complex exponential in Equation 5–57 can be expressed by sine and cosine as

(5–59)

Substituting Equation 5–59 into Equation 5–57 provides Equation 5–58

(5–60)

Comparing Equation 5–51 with Equation 5–60 reveals:

(5–61)

(5–62)

In words, the complex real, Q_r(r), and imaginary, Q_i(r), parts are the same as the measurable cosine, c(r), and sine, s(r), amplitudes.

A harmonic analysis provides two solution sets: the real and imaginary components of a complex solution. According to Equation 5–51, and Equation 5–61 the magnitude of the real and imaginary sets describe the measurable field at t = 0 and at ωt = -90 degrees, respectively. Comparing Equation 5–52 and Equation 5–61 provides:

(5–63)

(5–64)

Equation 5–63 expresses the amplitude (peak) and phase angle of the measurable harmonic field quantities by the complex real and imaginary parts.

The time average of harmonic fields such as A, E, B, H, J, or V is zero at point r. This is not the case for P, W, or F because they are quadratic functions of B, H, or J. To derive the time dependence of a quadratic function - for the sake of simplicity - we deal only with a Lorentz force, F, which is product of J and B. (This is a cross product but, for simplicity, the components are not shown. The space dependence is also omitted.)

(5–65)

where:

The time average of cos² and sin² terms is 1/2 whereas that of the sin cos term is zero. Therefore, the time average force density is:

(5–66)

Thus, the force can be obtained as the sum of "real" and "imaginary" forces. In a similar manner the time averaged Joule power density, Q^j, and magnetic energy density, W, can be obtained as:

(5–67)

(5–68)

Neglecting the time-derivative of magnetic flux density (the quasistatic approximation), the system of Maxwell's equations (Equation 5–1 through Equation 5–4) reduces to:

(5–69)

(5–70)

(5–71)

(5–72)

As follows from Equation 5–70, the electric field {E} is irrotational, and can be derived from:

(5–73)

where:

In the time-varying electromagnetic field governed by Equation 5–69 through Equation 5–72, the electric and magnetic fields are uncoupled. If only the electric solution is of interest, replacing Equation 5–69 by the continuity Equation 5–5 and eliminating Equation 5–71 produces the system of differential equations governing the quasistatic electric field.

Repeating Equation 5–12 and Equation 5–13 without velocity effects, the constitutive equations for the electric fields become:

(5–74)

(5–75)

where:

ρxx = resistivity in the x-direction (input as RSVX on MP command)

εxx = permittivity in the x-direction (input as PERX on MP command)

The conditions for {E}, {J}, and {D} on an electric material interface are:

(5–76)

(5–77)

(5–78)

where:

Two cases of the electric scalar potential approximation are considered below.

 

5.1.6.1. Quasistatic Electric Analysis

In this analysis, the relevant governing equations are Equation 5–73 and this continuity equation:

(5–79)

Substituting the constitutive Equation 5–74 and Equation 5–75 into Equation 5–79, and taking into account Equation 5–73, one obtains the differential equation for electric scalar potential:

(5–80)

Equation 5–80 is used to approximate a time-varying electric field in elements PLANE230, SOLID231, and SOLID232. It takes into account both the conductive and dielectric effects in electric materials. Neglecting time-variation of the electric potential, Equation 5–80 reduces to the governing equation for steady-state electric conduction:

(5–81)

In the case of a time-harmonic electric field analysis, the complex formalism allows Equation 5–80 to be re-written as:

(5–82)

where:

ω = angular frequency

Equation 5–82 is the governing equation for a time-harmonic electric analysis using elements PLANE121, SOLID122, and SOLID123.

In a time-harmonic analysis, the loss tangent tan δ can be used instead of or in addition to the electrical conductivity [σ] to characterize losses in dielectric materials. In this case, the conductivity matrix [σ] is replaced by the effective conductivity [σ^eff] defined as:

(5–83)

where:

tan δ = loss tangent (input as LSST on MP command)

 

5.1.6.2. Electrostatic Analysis

Electric scalar potential equation for electrostatic analysis is derived from governing Equation 5–72 and Equation 5–73, and constitutive Equation 5–75:

(5–84)

Equation 5–84, subject to appropriate boundary conditions, is solved in an electrostatic field analysis of dielectrics using elements PLANE121, SOLID122, and SOLID123.