Differences

This shows you the differences between two versions of the page.

--- users_manual [2020/10/20 21:22]
stanzurek [A.3.1 Truncation of Outer Boundaries]
+++ users_manual [2021/11/22 21:25] (current)
stanzurek [A.1 Modeling Permanent Magnets]
@@ Line 3227: / Line 3227: @@
 (A.1)  $  Hc = \frac{5·10^5· \sqrt{E}}{π}  $
-where //E// is the energy product in MGOe and the resulting //Hc// is in units of A/m (e.g. 40 MGOe ≈ 106 A/m).
+where //E// is the energy product in MGOe and the resulting //Hc// is in units of A/m (e.g. 40 MGOe ≈ 10<sup>6</sup> A/m).
 {{fig_a-3.png}}
@@ Line 3314: / Line 3314: @@
 === A.3.2 Asymptotic Boundary Conditions ===
-A thorough review of open boundary techniques is contained in [17]. Perhaps the simple way to
-approximate an “open” boundary (other than truncation) described in [17] is to use asymptotic
+A thorough review of open boundary techniques is contained in [17]. Perhaps the simple way to approximate an “open” boundary (other than truncation) described in [17] is to use asymptotic boundary conditions. The result is that by carefully specifying the parameters for the “mixed” boundary condition, and then applying this boundary condition to a circular outer boundary, the unbounded solution can be closely approximated. An example that employs an asymptotic boundary condition to obtain an unbounded field solution is the ''axi1.fem'' example included in the distribution.
-boundary conditions. The result is that by carefully specifying the parameters for the “mixed”
-boundary condition, and then applying this boundary condition to a circular outer boundary, the
+Consider a 2-D planar problem in polar coordinates. The domain is a circular shell of radius //ro// in an unbounded region. As //r → ∞//, vector potential //A// goes to zero. On the surface of the circle, the vector is a prescribed function of //θ//. This problem has an analytical solution, which is:
-unbounded solution can be closely approximated. An example that employs an asymptotic boundary
-condition to obtain an unbounded field solution is the axi1.fem example included in the distribution.
+(A.12)  $ A(r,θ) = \sum_{m=1}^{∞} \frac{a_m}{r^m} cos(mθ + α_m) $
-Consider a 2-D planar problem in polar coordinates. The domain is a circular shell of radius ro
-in an unbounded region. As r →¥, vector potential A goes to zero. On the surface of the circle,
+where the $a_m$ and $α_m$ parameters are chosen so that the solution matches the prescribed potential on the surface of the circle.
-the vector is a prescribed function of q. This problem has an analytical solution, which is:
-A(r,q) =
+One could think of this solution as describing the solution exterior to a finite element problem with a circular outer boundary. The solution is described inside the circle via a finite element solution. The trick is to knit together the analytical solution outside the circle to the finite element solution inside the circle.
-¥å
-m=1
+From inspecting (A.12), one can see that the higher-numbered harmonic, the faster the magnitude of the harmonic decays with respect to increasing r. After only a short distance, the higher numbered harmonics decay to the extent that almost all of the open-space solution is described by only the leading harmonic. If n is the number of the leading harmonic, the open-field solution for large, but not infinite, r is closely described by:
-am
-rm cos(mq+am) (A.12)
+(A.13) $ A(r,θ) ≈ \frac{a_n}{r^n} cos(nθ + α_n) $
-where the am and am parameters are chosen so that the solution matches the prescribed potential
-on the surface of the circle.
-One could think of this solution as describing the solution exterior to a finite element problem
-with a circular outer boundary. The solution is described inside the circle via a finite element
-solution. The trick is to knit together the analytical solution outside the circle to the finite element
-solution inside the circle.
-From inspecting (A.12), one can see that the higher-numbered harmonic, the faster the magnitude
-of the harmonic decays with respect to increasing r. After only a short distance, the highernumbered
-harmonics decay to the extent that almost all of the open-space solution is described by
-only the leading harmonic. If n is the number of the leading harmonic, the open-field solution for
-large, but not infinite, r is closely described by:
-A(r,q) ≈
-an
-rn cos(nq+an) (A.13)
 Differentiating with respect to r yields:
-¶A
-¶r
+(A.14) $ \frac{ \partial A }{\partial r}  = - \frac{n a_n}{r^{n+1}} cos(nθ + α_n) $
-= −
-nan
-rn+1 cos(nq+an) (A.14)
 If (A.14) is solved for an and substituted into (A.13), the result is:
-¶A
-¶r
+(A.15) $ \frac{ \partial A }{\partial r} + \left( \frac{n}{r} \right) A  = 0 $
-+
-n
+Now, (A.15) is a very useful result. This is the same form as the “mixed” boundary condition supported by FEMM. If the outer edge of the solution domain is circular, and the outer finite element boundary is somewhat removed from the area of primary interest, the open domain solution can be closely approximated by applying (A.15) the circular boundary.
-r
+To apply the Asymptotic Boundary Condition, define a new, mixed-type boundary condition. Then, pick the parameters so that:
-A = 0 (A.15)
-Now, (A.15) is a very useful result. This is the same form as the “mixed” boundary condition
+(A.16)  $ c_0 = \frac{n}{μ_0 r_0} $
-supported by FEMM. If the outer edge of the solution domain is circular, and the outer finite element
-boundary is somewhat removed from the area of primary interest, the open domain solution
+(A.17)  $ c_1 = 0 $
-can be closely approximated by applying (A.15) the circular boundary.
-To apply the Asymptotic Boundary Condition, define a new, mixed-type boundary condition.
+where $r_0$ is the outer radius of the region in meters (regardless of the working length units), and $μ_o = 4 π 10^{−7}$.
-Then, pick the parameters so that:
-c0 =
+Although the above derivation was specifically for 2-D problems, it turns out that when the same derivation is done for the axisymmetric case, the definition of the mixed boundary condition coefficients are exactly the same as (A.16).
-n
-μoro
+To apply the Asymptotic Boundary Condition to electrostatics problems, pick the parameters so that:
-(A.16)
-c1 = 0 (A.17)
+(A.18)  $ c_0 = \frac{ε_0 n}{r_0} $
-where ro is the outer radius of the region in meters (regardless of the working length units), and
-μo = 4p(10−7).
+$ c_1 = 0 $
-Although the above derivation was specifically for 2-D problems, it turns out that when the
-same derivation is done for the axisymmetric case, the definition of the mixed boundary condition
+where $r_0$ is the outer radius of the region in meters (regardless of the working length units), and ε0 = 8.85418781762e-012. Note that ε0 is defined in the Lua implementation in both the pre- and post-processors as the global variable ''eo'', which can be used in any script or edit box in the program.
-coefficients are exactly the same as (A.16).
-To apply the Asymptotic Boundary Condition to electrostatics problems, pick the parameters
+{{fig_a-7.png}}
-so that:
-c0 =
-eon
-ro
-(A.18)
-c1 = 0
-where ro is the outer radius of the region in meters (regardless of the working length units), and
-eo = 8.85418781762e-012. Note that eo is defined in the Lua implementation in both the preand
-post-processors as the global variable eo, which can be used in any script or edit box in the
-program.
 Figure A.7: Air-cored coil with “open” boundary condition
-Just like magnetics problems, it turns out that 2-D problems are also described by (A.18). One
-subtle difference, however, is that n =1 in the axisymmetric case corresponds to the case in which
+Just like magnetics problems, it turns out that 2-D problems are also described by (A.18). One subtle difference, however, is that n =1 in the axisymmetric case corresponds to the case in which there is a net charge with (i.e. the geometry looks like a point charge when viewed from a distance), whereas the n= 1 case corresponds to a dipole charge distribution in 2D planar problems. If charge is conserved in the geometry of interest in the axisymmetric case, one needs to use n =2 in Eq. (A.18). It should be noted that this is a departure from the magnetostatic case with the vector potential formulation in which n = 1 corresponds to a dipole arrangement in the axisymmetric case.
-there is a net charge with (i.e. the geometry looks like a point charge when viewed from a distance),
-whereas the n= 1 case corresponds to a dipole charge distribution in 2D planar problems. If charge
+Some care must be used in applying this boundary condition. Most of the time, it is sufficient to take n = 1 (i.e the objects in the solution region look like a dipole when viewed from a large distance). However, there are other cases (e.g. a 4-pole halbach permanent magnet array) in which the leading harmonic is something other than n= 1. You need to use your insight into your specific problem to pick the appropriate n for the leading harmonic. You also must put the objects of interest roughly in the center of the circular finite element domain to minimize the magnitude higher-order field components at the outer boundary.
-is conserved in the geometry of interest in the axisymmetric case, one needs to use n =2 in Eq.
-(A.18). It should be noted that this is a departure from the magnetostatic case with the vector
+Although the application of this boundary condition requires some thought on the part of the user, the results can be quite good. Figure A.7, corresponding to the axi1 example, represents the field produced by an air-cored coil in free space. The asymptotic boundary condition has been applied to the circular outer boundary. Inspecting the solution, flux lines appear to cross the circular boundary as if the solution domain were truly unbounded.
-potential formulation in which n = 1 corresponds to a dipole arrangement in the axisymmetric
-case.
+A quick note on computational efficiency: applying the absorbing boundary condition imposes no additional computing cost on the problem. The ABC is computationally no more time-consuming to apply than enforcing A = 0 at the outer boundary. Solution times for the PCG solver are equivalent in either case. It can also readily be derived that the ABC works exactly the same for harmonics problems. (To see this, just assume that the am in (A.12) can be complex valued, and follow the same derivation).
-Some care must be used in applying this boundary condition. Most of the time, it is sufficient
-to take n = 1 (i.e the objects in the solution region look like a dipole when viewed from a large
-distance). However, there are other cases (e.g. a 4-pole halbach permanent magnet array) in which
+=== A.3.3 Kelvin Transformation ===
-the leading harmonic is something other than n= 1. You need to use your insight into your specific
-problemto pick the appropriate n for the leading harmonic. You alsomust put the objects of interest
+A particularly good approach to “open boundary” problems is the Kelvin Transformation, a technique first discussed in the context of computational magnetics in [18] and [19]. The strengths of this technique are:
-roughly in the center of the circular finite element domain to minimize the magnitude higher-order
+  * the effects of the exterior region are, in theory, exactly modeled by this approach;
-field components at the outer boundary.
+  * a sparse matrix representation of the problem is retained (unlike FEM-BEM methods, which give the same “exact solution” but densely couples together the boundary nodes).
-Although the application of this boundary condition requires some thought on the part of the
+  * requires no “special” features in the finite element solver to implement the technique, other than the ability to apply periodic boundary conditions.
-user, the results can be quite good. Figure A.7, corresponding to the axi1 example, represents
-the field produced by an air-cored coil in free space. The asymptotic boundary condition has been
+The purposes of this note are to explain what the Kelvin transformation is derived and to show how it is implemented in the context of the FEMM finite element program.
-applied to the circular outer boundary. Inspecting the solution, flux lines appear to cross the circular
-boundary as if the solution domain were truly unbounded.
+** Derivation **
-A quick note on computational efficiency: applying the absorbing boundary condition imposes
+In the “far field” region, the material is typically homogeneous (e.g. air and free of sources. In this case, the differential equation that describes vector potential A is the Laplace equation:
-no additional computing cost on the problem. The ABC is computationally no more timeconsuming
-to apply than enforcing A = 0 at the outer boundary. Solution times for the PCG solver
+(A.19) $ \nabla^2 A = 0  $
-are equivalent in either case. It can also readily be derived that the ABC works exactly the same
-for harmonics problems. (To see this, just assume that the am in (A.12) can be complex valued,
+If we write (A.19) in polar notation, //A// is described by:
-and follow the same derivation).
-A.3.3 Kelvin Transformation
+(A.20) $ \frac{1}{r} \left( r \frac{\partial A}{\partial r} \right) + \frac{1}{r^2} \frac{\partial^2 A}{\partial θ^2} = 0  $
-A particularly good approach to “open boundary” problems is the Kelvin Transformation, a technique
-first discussed in the context of computational magnetics in [18] and [19]. The strengths of
+Assume that the “near field” region of the problem can be contained in a circle of radius //ro// centered at the origin. The far-field region is then everything outside the circle.
-this technique are:
-• the effects of the exterior region are, in theory, exactly modeled by this approach;
+One approach to unbounded problems is to attempt to map the unbounded region onto a bounded region, wherein problems can more easilby be solved. Specifically, we desire a way to transform the unbounded region outside the circle into a bounded region. One simple way to make such a mapping is to define another variable, //R//, that is related to //r// by:
-• a sparse matrix representation of the problem is retained (unlike FEM-BEM methods, which
-give the same “exact solution” but densely couples together the boundary nodes).
+(A.21) $ R = \frac{r_0^2}{r} $
-• requires no “special” features in the finite element solver to implement the technique, other
-than the ability to apply periodic boundary conditions.
+By inspecting (A.21), it can be seen that this relationship maps the exterior region onto a circle of radius //ro//.
-The purposes of this note are to explain what the Kelvin transformation is derived and to show how
-it is implemented in the context of the FEMM finite element program.
+The next step is to transform (A.19), the differential equation that the field must satisfy, into the mapped space. That is, (A.19) must be written in terms of //R// and //θ// rather than //r// and //θ//. We can evaluate derivatives in terms of //R// instead of //r// by employing the chain rule:
-Derivation
-In the “far field” region, the material is typically homogeneous (e.g. air and free of sources. In this
+(A.22) $ \frac{\partial}{\partial r} = \frac{\partial}{\partial R} \left( \frac{\partial R}{\partial r} \right) = - \frac{\partial}{\partial R} \left( \frac{R}{r_0} \right)^2  $
-case, the differential equation that describes vector potential A is the Laplace equation:
-Ñ2A = 0 (A.19)
+Now, we can note that at //r// = //R// = //ro//,
-If we write (A.19) in polar notation, A is described by:
+(A.23)  $ \frac{\partial A}{\partial r} = -\frac{\partial A}{\partial R} $
-r
-¶
-¶r
-r
-¶A
-¶r
-+
-r2
-¶2A
-¶q2 = 0 (A.20)
-Assume that the “near field” region of the problem can be contained in a circle of radius ro centered
-at the origin. The far-field region is then everything outside the circle.
-One approach to unbounded problems is to attempt to map the unbounded region onto a
-bounded region, wherein problems can more easilby be solved. Specifically, we desire a way
-to transform the unbounded region outside the circle into a bounded region. One simple way to
-make such a mapping is to define another variable, R, that is related to r by:
-R =
-r2
-o
-r
-(A.21)
-By inspecting (A.21), it can be seen that this relationship maps the exterior region onto a circle of
-radius ro.
-The next step is to transform (A.19), the differential equation that the field must satisfy, into
-the mapped space. That is, (A.19) must be written in terms of R and q rather than r and q. We can
-evaluate derivatives in terms of R instead of r by employing the chain rule:
-¶
-¶r
-=
-¶
-¶R
-dR
-dr
-= −
-¶
-¶R
-R
-ro
-(A.22)
-Now, we can note that at r = R = ro,
-¶A
-¶r
-= −
-¶A
-¶R
-(A.23)
 and we can substitute (A.22) into (A.19) to yield, after some algebraic manipulation:
-R
+(A.24) $  \frac{1}{R} \frac{\partial}{\partial R} \left(R \frac{\partial A}{\partial R} \right) + \frac{1}{R^2}  \frac{\partial^2 A}{\partial θ^2}  = 0 $
-¶
-¶R
+Eq. (A.24), the transformed equation for the outer region, has exactly the same form as inner region, only in terms of //R// rather than //r//. The implication is that for the 2-D planar problem, the exterior can be modeled simply by creating a problem domain consisting of two circular regions: on circular region containing the items of interest, and an additional circular region to represent the “far field.” Then, periodic boundary conditions must be applied to corresponding edges of the circle to enforce the continuity of //A// at the edges of the two regions. The is continuity of //A// at the boundary between the exterior and interior regions. For a finite element formulation consisting of first-order triangles, (A.23) is enforced automatically at the boundaries of the two regions. The second circular region exactly models the infinite space solution, but does it on a bounded domain - one could always back out the field for any point in space by applying the inverse of (A.21).
-R
+** Kelvin Transformation Example ** – ''open1.fem''
-¶A
-¶R
+As an example, consider an E-core lamination stack with a winding around it. Suppose that the objective is to determine the field around the E-core in the absence of any flux return path (i.e. when the magnetic circuit is open). In this case, the flux is not constrained to flow in a path that is a priori well defined, because the laminations that complete the flux path have been removed.
-+
+The geometry was chosen arbitrarily, the purpose here being more the procedure than the actual problem. The E-core was chosen to have a 0.5” thick center leg, 0.25” thick outer legs, and a slot depth of 0.75”. The material for the core is linear with a relative permeability of 2500. The coil carries a bulk current density of 2 MA/m2. The input geometry is picture in Figure A.8.
-R2
+In Figure A.8, the core is placed within a circular region, and a second circular region is drawn next to the region containing the core. Periodic boundary conditions are applied to the arcs that define the boundaries as shown in Figure A.8. The way that periodic boundary conditions are implemented in FEMM, each periodic boundary condition defined for the problem is to be applied to two and only two corresponding entities. In this case, each boundary circle is composed of two arcs, so two periodic boundary conditions must be defined to link together each arc with in the domain with the core to its corresponding arc in the domain representing the exterior region.
-¶2A
-¶q2 = 0 (A.24)
+Also notice that a point has been drawn in the center of the exterior region. A point property has been applied to this point that specifies that //A// = 0 at this reference point. The center of the circle maps to infinity in the analogous open problem, so it makes sense to define, in effect, //A// = 0 at infinity. If no reference point is defined, it is fairly easy to see that the solution is only unique to within a constant. The situation is analogous to a situation where Neumann boundary conditions have been defined on all boundaries, resulting in a non-unique solution for A. Due to the type of solver that FEMM employs, the problem can most likely be solved even if a reference point is not defined. However, defining a reference point eliminates the possibility of numerical difficulties due to uniqueness issues.
-Eq. (A.24), the transformed equation for the outer region, has exactly the same form as inner
-region, only in terms of R rather than r. The implication is that for the 2-D planar problem, the
+{{fig_a-8.png}}
-exterior can be modeled simply by creating a problem domain consisting of two circular regions:
-on circular region containing the items of interest, and an additional circular region to represent
-the “far field.” Then, periodic boundary conditions must be applied to corresponding edges of the
-circle to enforce the continuity of A at the edges of the two regions. The is continuity of A at the
-boundary between the exterior and interior regions. For a finite element formulation consisting
-of first-order triangles, (A.23) is enforced automatically at the boundaries of the two regions. The
-second circular region exactly models the infinite space solution, but does it on a bounded domain–
-one could always back out the field for any point in space by applying the inverse of (A.21).
-Kelvin Transformation Example – open1.fem
-As an example, consider an E-core lamination stack with a winding around it. Suppose that the
-objective is to determine the field around the E-core in the absence of any flux return path (i.e.
-when the magnetic circuit is open). In this case, the flux is not constrained to flow in a path that is
-a priori well defined, because the laminations that complete the flux path have been removed.
-The geometry was chosen arbitrarily, the purpose here being more the procedure than the actual
-problem. The E-core was chosen to have a 0.5” thick center leg, 0.25” thick outer legs, and a slot
-depth of 0.75”. The material for the core is linear with a relative permeability of 2500. The coil
-carries a bulk current density of 2 MA/m2. The input geometry is picture in Figure A.8.
-In Figure A.8, the core is placed within a circular region, and a second circular region is drawn
-next to the region containing the core. Periodic boundary conditions are applied to the arcs that
-define the boundaries as shown in Figure A.8. The way that periodic boundary conditions are
-implemented in FEMM, each periodic boundary condition defined for the problem is to be applied
-to two and only two corresponding entities. In this case, each boundary circle is composed of two
-arcs, so two periodic boundary conditions must be defined to link together each arc with in the
-domain with the core to its corresponding arc in the domain representing the exterior region.
-Also notice that a point has been drawn in the center of the exterior region. A point property
-has been applied to this point that specifies that A = 0 at this reference point. The center of the
-circle maps to infinity in the analogous open problem, so it makes sense to define, in effect, A = 0
 Figure A.8: Example input geometry.
-at infinity. If no reference point is defined, it is fairly easy to see that the solution is only unique to
-within a constant. The situation is analogous to a situation where Neumann boundary conditions
-have been defined on all boundaries, resulting in a non-unique solution for A. Due to the type of
+The resulting solution is shown in Figure A.9. As is the intention, the flux lines appear to cross out of the of the region containing the core as if unaffected by the presence of the boundary. The flux lines reappear in the domain representing the exterior region, completing their flux paths through the exterior region.
-solver that FEMM employs, the problem can most likely be solved even if a reference point is not
-defined. However, defining a reference point eliminates the possibility of numerical difficulties
-due to uniqueness issues.
+==== A.4 Nonlinear Time Harmonic Formulation ====
-The resulting solution is shown in Figure A.9. As is the intention, the flux lines appear to
-cross out of the of the region containing the core as if unaffected by the presence of the boundary.
+Starting with the the 3.3 version of FEMM, the program includes a “nonlinear time harmonic” solver. In general, the notion of a “nonlinear time harmonic” analysis is something of a kludge. To obtain a purely sinusoidal response when a system is driven with a sinusoidal input, the system must, by definition, be linear. The nonlinear time harmonic analysis seeks to include the effects of nonlinearities like saturation and hysteresis on the fundamental of the response, while ignoring higher harmonic content. This is a notion similar to “describing function analysis,” a widely used tool in the analysis of nonlinear control systems. There are several subtly different variations of the formulation that can yield slightly different results, so documentation of what has actually been implement is important to the correct interpretation of the results from this solver.
-The flux lines reappear in the domain representing the exterior region, completing their flux paths
-through the exterior region.
+An excellent description of this formulation is contained in [21]. FEMM formulates the nonlin-ear time harmonic problem as described in this paper. Similar to Jack and Mecrow, FEMM derives an apparent BH curve by taking H to be the sinusoidally varying quantity. The amplitude of B is obtained by taking the first coefficient in a Fourier series representation of the resulting B. For the purposes of this Fourier series computation, FEMM interpolates linearly between the user-defined points on the BH curve to get a set of points with the same H values as the input set, but with an adjusted B level. The rationale for choosing H to be the sinusoidal quantity (rather than B) is that choosing B to be sinusoidal shrinks the defined BH curve–the B values stay fixed while the H values become smaller. It then becomes hard to define a BH curve that does not get interpolated. In contrast, with H sinusoidal, the B points are typically larger than the DC flux density levels, creating a curve with an expanded range.
-A.4 Nonlinear Time Harmonic Formulation
-Starting with the the 3.3 version of FEMM, the program includes a “nonlinear time harmonic”
+{{fig_a-9.png}}
-solver. In general, the notion of a “nonlinear time harmonic” analysis is something of a kludge.
-To obtain a purely sinusoidal response when a system is driven with a sinusoidal input, the system
-must, by definition, be linear. The nonlinear time harmonic analysis seeks to include the effects
-of nonlinearities like saturation and hysteresis on the fundamental of the response, while ignoring
-higher harmonic content. This is a notion similar to “describing function analysis,” a widely used
-tool in the analysis of nonlinear control systems. There are several subtly different variations of
-the formulation that can yield slightly different results, so documentation of what has actually been
-implement is important to the correct interpretation of the results from this solver.
-An excellent description of this formulation is contained in [21]. FEMMformulates the nonlin-
 Figure A.9: Solved problem.
-ear time harmonic problem as described in this paper. Similar to Jack and Mecrow, FEMM derives
-an apparent BH curve by taking H to be the sinusoidally varying quantity. The amplitude of B is
+A “nonlinear hysteresis lag” parameter is then applied to the effective BH curve. The lag is assumed to be proportional to the permeability, which gives a hysteresis loss that is always proportional to |B|<sup>2</sup>. This form was suggested by O’Kelly [22]. It has been suggested that that the Steinmetz equation could be used to specify hysteresis lag, but the Steinmetz equation is badly behaved at low flux levels (i.e. one can’t solve for a hysteresis lag that produces the Steinmetz |B|<sup>1.6</sup> form for the loss as B goes to zero.)
-obtained by taking the first coefficient in a Fourier series representation of the resulting B. For the
-purposes of this Fourier series computation, FEMM interpolates linearly between the user-defined
+For nonlinear in-plane laminations, an additional step is taken to obtain an effective BH curve that also includes eddy current effects. At each H level on the user-defined BH curve, a 1D nonlinear time harmonic finite element problem is solved to obtain the total flux that flows in the lamination as a function of the H applied at the edge of the lamination. Then dividing by the lamination thickness and accounting for fill factor, and effective B that takes into account saturation, hysteresis, and eddy currents in the lamination is obtained for each H.
-points on the BH curve to get a set of points with the same H values as the input set, but with
-an adjusted B level. The rationale for choosing H to be the sinusoidal quantity (rather than B) is
-that choosing B to be sinusoidal shrinks the defined BH curve–the B values stay fixed while the H
-values become smaller. It then becomes hard to define a BH curve that does not get interpolated.
-In contrast, with H sinusoidal, the B points are typically larger than the DC flux density levels,
-creating a curve with an expanded range.
-A “nonlinear hysteresis lag” parameter is then applied to the effective BH curve. The lag
-is assumed to be proportional to the permeability, which gives a hysteresis loss that is always
-proportional to |B|2. This form was suggested by O’Kelly [22]. It has been suggested that that
-the Steinmetz equation could be used to specify hysteresis lag, but the Steinmetz equation is badly
-behaved at low flux levels (i.e. one can’t solve for a hysteresis lag that produces the Steinmetz
-|B|1.6 form for the loss as B goes to zero.)
-For nonlinear in-plane laminations, an additional step is taken to obtain an effective BH curve
-that also includes eddy current effects. At each H level on the user-defined BH curve, a 1D nonlinear
-time harmonic finite element problem is solved to obtain the total flux that flows in the
-lamination as a function of the H applied at the edge of the lamination. Then dividing by the lamination
-thickness and accounting for fill factor, and effective B that takes into account saturation,
-hysteresis, and eddy currents in the lamination is obtained for each H.

Unofficial FEMM documentation by Dr Stan Zurek

User Tools

Site Tools

Differences

Page Tools

Unofficial
FEMM
documentation
by Dr Stan Zurek