Quickfield Manual Pocket
Quickfield Manual Pocket
Acknowledgments:
TetGen (www.tetgen.org)
CGAL (www.cgal.org)
VTK (www.vtk.org)
C3D Toolkit (www.c3dlabs.com)
VisTools/Mesh (www.vki.com)
All other brand and product names are trademarks or registered trademarks of their
respective owners
iii
Contents
Index 390
13
What Is QuickField?
Welcome to QuickField Finite Elements Analysis System.
QuickField is a PC-oriented interactive environment for
electromagnetic, thermal and stress analysis. Standard
analysis types include:
• Electrostatics.
• DC and AC conduction analysis.
• Linear and nonlinear DC and transient magnetics.
• AC magnetics (involving eddy current analysis).
• Linear and nonlinear, steady state and transient heat
transfer and diffusion.
• Linear stress analysis.
• Coupled problems.
During a 15-minute session, you can describe the problem
(geometry, material properties, sources and other
conditions), obtain solution with high accuracy and analyze
field details looking through full color picture. With
QuickField, complicated field problems can be solved on
your PC instead of large mainframes or workstations.
Conventions
In this manual we use SMALL CAPITAL LETTERS to specify
the names of keys on your keyboard. For example, ENTER,
ESC, or ALT. Four arrows on the keyboard, collectively
named the DIRECTION keys, are named for the direction the
key points: UP ARROW, DOWN ARROW, RIGHT ARROW, and
LEFT ARROW.
A plus sign (+) between key names means to hold down the
first key while you press the second key. A comma (,)
between key names means to press the keys one after the
other.
Bold type is used for QuickField menu and dialog options.
16
C H A P T E R 1
Getting Started
QuickField Installation
QuickField can be supplied on a CD, or packed in ZIP-
archive. Depending on the format, do the following:
QuickField Installation 17
Autorun Applet
On the left side of the Autorun screen you can see several
menu topics organized in a scrollable tree. When you
highlight a topic, additional topic-related information
appears in the bottom pane. To execute the command
associated with this topic double-click it or click Run in the
right-bottom corner of the window.
Menu topics allow you to:
• See the complete QuickField User's Guide in Adobe PDF
format (Read User’s Guide command);
• Learn QuickField interactively (Virtual classroom
command);
• Find technical support and sales contact information
(Contact Us);
• Install QuickField (Install QuickField command
group).
With the Student version of QuickField the last command
starts QuickField installer.
With QuickField Professional the installation steps depend
on the type of license you purchased. For single-user
license, choose the Single-user QuickField option below.
18 Chapter 1 Getting Started
Application-Wide Options
QuickField Options panel may be displayed using the
command Tools>Options… in the main menu. Left pane of
the window shows the list of the options, and the pages for
viewing/editing the chosen option are shown in the right
pane. Pressing the OK button saves the options. Options are
stored separately for each user of the computer.
22 Chapter 1 Getting Started
C H A P T E R 2
Introductory Guide
Problem Window
The problem window is normally docked on the left side of
the main QuickField window. When several problems are
open at once, their windows can be docked side by side, or
in a column, or they can be tabbed in a single pane, leaving
maximum space for graphics. This window can also be left
floating on top of the other windows.
To move the problem window, simply drag it to the new
position holding by the window title. While dragging, the
possible docking positions are shown by the diamond
shaped arrows. When you move the pointer over the
28 Chapter 2 Introductory Guide
Document Windows
QuickField document windows, such as the model editor,
the field plot, or the electric circuit window, occupy the
main are of the QuickField workspace; they cannot be
docked. For fast switching between these windows, there is
a tab bar near the bottom edge of the main window, similar
to the Windows task bar.
A document window can be minimized to an icon, or
maximized, or arranged with its regular size and position,
which you can change by dragging any corner or edge of the
window. This is useful to display several document
windows at once. QuickField can also tile all document
windows automatically, when you choose Tile Vertically
or Tile Horizontally in the Window menu.
Some document windows can be split into two or four
panes. To split the window, drag the small gray rectangle on
top of the vertical scrollbar or on the left of the horizontal
scrollbar. You can also choose Split in the Window menu.
To switch between panes, click it with the mouse or use F6.
To restore the single view, double-click the splitter or drag
it to the window border until it disappears.
Tool Windows
Finally, the tool windows–the field calculator, color legend,
circuit elements list, etc. –are usually docked within the
corresponding document window. Like the problem
Overview of Analysis Capabilities 29
window, you can drag and dock tools within their parent’s
boundaries. When floating, tools can be dragged anywhere
on the screen, even to another monitor.
Properties Window
The Properties window can be opened using the Properties
command in the View menu. This window is docked to the
problem window (as shown in the picture) by default or can
be switched to floating. The Properties window displays
different editing fields relevant to the current object (the
problem, geometry model, etc.) Some properties are for
information only (shown in grey), the others can be changed
by typing in the new value or selecting from the dropdown
list. The changed property value comes into effect
immediately.
Magnetostatic Analysis
Magnetic analysis is used to design or analyze variety of
devices such as solenoids, electric motors, magnetic shields,
permanent magnets, magnetic disk drives, and so forth.
Generally, the quantities of interest in magnetostatic
analysis are magnetic flux density, field intensity, forces,
torques, inductance, and flux linkage.
QuickField can perform linear and nonlinear magnetostatic
analysis for 2-D and axisymmetric models. The program is
30 Chapter 2 Introductory Guide
AC Magnetic Analysis
AC magnetic analysis is used to analyze magnetic field
caused by alternating currents and, vise versa, electric
currents induced by alternating magnetic field (eddy
currents). This kind of analysis is useful with different
inductor devices, solenoids, electric motors, and so forth.
Generally the quantities of interest in AC magnetic analysis
Overview of Analysis Capabilities 33
Electrostatic Analysis
Electrostatic analysis is used to design or analyze variety of
capacitive systems such as fuses, transmission lines and so
forth. Generally the quantities of interest in electrostatic
analysis are voltages, electric fields, capacitances, and
electric forces.
QuickField can perform linear electrostatic analysis for two-
dimensional, axisymmetric and three-dimensional models.
The program is based on Poisson's equation. Following
options are available for electrostatic analysis:
Material properties: air, orthotropic materials with
constant permittivity.
Overview of Analysis Capabilities 35
DC Conduction Analysis
DC conduction analysis is used to analyze variety of
conductive systems. Generally, the quantities of interest in
DC conduction analysis are voltages, current densities,
electric power losses (Joule heat).
QuickField can perform linear DC conduction analysis for
two-dimensional, axisymmetric and three-dimensional
models. The program is based on Poisson's equation.
Following options are available for DC conduction analysis:
36 Chapter 2 Introductory Guide
AC Conduction Analysis
AC conduction analysis is used to analyze electric field
caused by alternating currents and voltages in imperfect
dielectric media. This kind of analysis is mostly used with
complex insulator systems and capacitors. Generally, the
Overview of Analysis Capabilities 37
Thermal Analysis
Thermal analysis plays an important role in design of many
different mechanical and electrical systems. Generally the
quantities of interest in thermal analysis are temperature
distribution, thermal gradients, and heat losses. Transient
analysis allows you to simulate transition of heat
distribution between two heating states of a system.
Overview of Analysis Capabilities 39
Stress Analysis
Stress analysis plays an important role in design of many
different mechanical and electrical components. Generally,
the quantities of interest in stress analysis are displacements,
strains and different components of stresses.
QuickField can perform linear stress analysis for 2-D plane
stress, plane strain, and axisymmetric models. The program
is based on Navier equations of elasticity. Following options
are available for stress analysis:
Material properties: isotropic and orthotropic materials.
Loading sources: concentrated loads, body forces,
pressure, thermal strains, and imported electric or magnetic
forces from electric or magnetic analysis.
Boundary conditions: prescribed displacements, elastic
spring supports.
Postprocessing results: displacements, stress components,
principal stresses, von Mises stress, Tresca, Mohr-Coulomb,
Drucker-Prager, and Hill criteria.
41
C H A P T E R 3
Problem Description
Editing Problems
• To create a new, empty problem description, click New
in the File menu and then select QuickField problem in
the list that appears. Then enter the name and path of the
new problem. You can also create a new problem as a
copy of another problem being currently opened. In that
case new problem inherits all the properties of the sample
one and the referenced model and data documents are
copied if necessary.
• To open an existing document, click Open in the File
menu, or use drag and drop features of Windows.
Open problem documents are shown in a special view to the
left of main QuickField window. In problem view, you can
edit problem description options and references to files. The
tree shows the names of files, which the problem currently
references.
• To change problem settings or file names, click
Problem Properties in the Problem menu or context
(right mouse button) menu.
• To start editing a referenced document (model, data,
secondary data or other problem referenced as coupling
link), double-click its name in the tree, or click Edit File
in the context menu, or click correspondent item in Edit
menu.
• To solve the problem, click Solve Problem in the
Problem menu or context (right mouse button) menu.
Editing Problems 43
C H A P T E R 4
Terminology
Geometric Model, or simply Model, is the name we use for
the collection containing all geometric shapes of a problem.
Besides being an object container, the model helps to link
the contained objects with related material properties, field
sources, and boundary conditions.
Vertex, edge and block are three basic types of geometric
objects contained by QuickField models.
Each Vertex represents a point. Point coordinates could be
either explicitly specified by user or automatically
calculated by QuickField at the intersection of two edges.
For each vertex you can define its mesh spacing value and
Terminology 55
Geometry Description
Model development consists of three stages:
• Geometry description and manipulation;
• Definition of properties, field sources and boundary
conditions;
• Mesh generation.
Creating Vertices
To create new vertices:
• Choose Insert Mode in the Edit menu, or click the
Insert Vertices/Edges toolbar button or context menu
item, or press INS, to switch model view into insert mode.
• Make sure that current coordinate grid settings fit
coordinates of the vertices you want to create.
• Use mouse or DIRECTION keys to move the cursor to the
vertex insertion point and double-click the left mouse
button or press ENTER.
Or:
• Choose Add Vertices from the Edit menu.
• Enter new vertex coordinates and click Add. Repeat if
you need more vertices.
• Click Close.
Attraction Distance
To avoid small unrecognizable inaccuracies in geometry
definition, new vertices or edges cannot be created very
close to the existing objects. Creation of new geometric
objects is controlled by the value we denote by and call the
attraction distance.
The following rules apply to creation of new vertices and
edges.
• New vertices cannot be created within 2-neighborhood
of the existing vertex.
• New edge cannot connect the ends of the existing edge
and lie inside its -neighborhood.
Geometry Description 59
• Press INS to leave the Insert Mode. You will get the
model shown in Pic.2.
Pic.1 Pic.2
• Select the new edge dragging the left mouse button from
(-0.25, 3.25) to (4.25, 2.75).
• Choose Duplicate Selection from the Edit menu, set
displacement ordinate to -0.5, set Copies to 12, and click
OK. You will get the model shown in Pic.3.
• Select the right ends of horizontal edges dragging the left
mouse button from (3.75, 3.25) to (4.25, -3.25).
• Choose Delete Selection from the Edit menu and click
Yes to confirm deletion.
• Select the left ends of the edges dragging the left mouse
button from (-0.25, 1.75) to (0.25, -1.75) and delete them
similarly. You will get the required model (see Pic.4).
Geometry Description 65
Pic.3 Pic.4
Keyboard shortcuts:
Delete DEL
When you drag the objects the anchor point is also dragged.
QuickField keeps displaying it as a bright red dot. The
dragged anchor always coincides with one of the model
vertices or, unless the snap to grid option is off, with one of
the background grid nodes. You can see the coordinates of
the dragged anchor point in the status bar.
After the drop QuickField calculates the difference between
initial and final anchor positions and shifts all dragged
objects exactly for the length of that vector.
Example: Suppose that you want to move a group of model
objects containing the point with coordinates (a, b). After
the move the new coordinates of the point should be (c, d).
Here is the sequence of required actions:
• If there is no vertex at (a, b), add it choosing Add Nodes
from the Edit menu and entering the coordinates in the
dialog.
• If there is no vertex at (c, d), add it in the same way.
• Select the objects to move including the vertex at (a, b).
• Place mouse pointer over this vertex and press the left
button. You will see the anchor at (a, b).
• Drag the objects until the anchor coincides with the
vertex at (c, d) and release the mouse. The first vertex
will be moved exactly to (c, d).
• Delete one or both of the created vertices (in most cases,
the first vertex will not exist after move) if you no longer
need them.
Visual Drag Effects
To help you drag and drop objects correctly QuickField
provides visual feedback consisting of:
68 Chapter 4 Model Geometry Definition
Note. If you drag with right mouse button make sure that the
ALT key is released before the drop. If you release the right
mouse button with the ALT key pressed QuickField will do
nothing.
Undo/Redo Operations
To undo the latest operation, make sure that the active
window shows the geometric model you are editing and
choose Undo <your last operation> from the Edit menu.
To redo the last operation undone, make sure that the active
window shows the geometric model you are editing and
choose Redo <last operation undone> from the Edit
menu. QuickField modifies the corresponding menu items
to show you which operations would be undone and redone.
By default QuickField allows you to undo 25 latest
operations for every model. You cannot increase the number
of undoable operations above 100 but you can make it any
number between 0 and 100 at any time. You can find the
detailed description of this feature in Undo Settings.
The operations that could be undone are listed in Undoable
Operations.
Keyboard shortcuts:
Undo CTRL+Z
Redo CTRL+Y
Undo Settings
To be able to undo and redo your editing operations
QuickField maintains internal stack of increments to
76 Chapter 4 Model Geometry Definition
If the new depth value exceeds the old one, only the depth
of the stack is affected. Stack contents remains unchanged
allowing you to undo and redo the operations performed
before this change. The same is true if the new depth value
is less than the old one but still exceeds the number of
accumulated stack positions.
Suppose, for example, that you started QuickField,
performed 10 editing operations and invoked Undo 5 times.
At this moment the stack accumulates 10 stack positions. If
you decide to reduce stack depth to 10 no database
increments will be lost. If, however, you decide to perform
one more editing operation before reducing the depth, you
lose the possibility to redo anything and the number of
accumulated positions becomes equal to 6. After that you
can set the depth to 6 and lose nothing.
The only case you lose some of accumulated data changing
undo stack depth is when you reduce the depth to the value
that is less than the number of accumulated stack positions.
In such case, QuickField performs the following actions:
• if the new depth value exceeds the number of positions
below current stack position, QuickField discards
several topmost database increments to make the number
of remaining increments equal to new stack depth;
• if the new depth value is less than the number of
positions below current stack position, QuickField
retains only the database increments at the current stack
position and immediately below it making the number of
retained database increments equal to new stack depth.
Suppose once again that you started QuickField, performed
10 editing operations and invoked Undo 5 times. If at this
78 Chapter 4 Model Geometry Definition
• Unselect All
Meshing Technology
Having described a reasonable part of the model geometry
you can start building the finite element mesh. It is
80 Chapter 4 Model Geometry Definition
Zooming
To match window size with the size of the model:
• Click Zoom to Fit toolbar button.
To magnify the picture:
• Click Zoom In toolbar button.
• With left button pressed drag the mouse diagonally
drawing a rubberband rectangle around that part of the
model you want to fill the window, and release the
button.
Or
• Click inside the window. QuickField will use
magnification factor 2 relative to the clicked point.
Shortcut: CTRL + emulates click at the
central point of the window.
To see more of the model:
• Click Zoom Out toolbar button.
Shortcut: CTRL -
Tuning the View of the Model 85
Background Grid
Grid makes creation of model vertices and edges easier and
helps to check correctness of the model. To change the grid,
choose Grid Settings in the Edit or context menu and
change the dialog fields described below according to your
needs.
86 Chapter 4 Model Geometry Definition
C H A P T E R 5
Cartesian Polar
α
α
α α α α
α
When you associate data with a new label, the text boxes for
magnetic permeability components contain None. The word
None in these boxes or absence of the values mean blocks
with the label will be excluded from calculations. To define
the material’s properties (thereby including the blocks into
calculation), type in the required value of magnetic
permeability.
As long as Anisotropic is unchecked, QuickField changes
the components of magnetic permeability tensor
synchronously. To specify different components, check
Anisotropic before entering the required values. The dialog
Editing Label Data 97
The ways you define field sources for transient and non-
transient problems are slightly different. Also transient
problems allow using two types of conductors, stranded and
solid. QuickField distinguishes between these types by the
specified electric conductivity. Zero electric conductivity
value implies stranded conductor and no eddy currents in
the block. Non-zero value implies solid conductor and
forces QuickField to calculate eddy current distribution for
the block.
98 Chapter 5 Problem Parameters Description
There are two main causes for the magnetic material losses:
ohmic losses generated by eddy currents, and losses caused
by the cyclic reversal of the magnetization and proportional
to the area of the hysteresis loop. If the non-zero
conductivity is specified, then for this material eddy current
losses are calculated automatically and there is no need to
specify the loss coefficients separately.
Situation is different for the laminated cores, which are
made of thin isolated sheets. In the laminated cores the
electrical conductivity should be set as 0, otherwise the eddy
currents simulated with QuickField will be too large.
However, even the small eddy currents generate the losses.
They are usually taken into account using some empirical
formula.
QuickField uses the empirical Bertotti formula for core loss
calculation:
p = kh·f·B2 + kc·f 2· B2 + ke· (f·B)1.5
Here B – magnitude of the module of the flux density vector
per period, f – problem frequency, kh, kc, ke – volume power
loss coefficients for specific magnetic material. Default zero
values of the loss coefficient will exclude the corresponding
loss component from the calculations.
106 Chapter 5 Problem Parameters Description
1. Elasticity
When you associate data with a new label, the text boxes for
Young's moduli contain None. The word None in these
boxes or absence of the values means that blocks with the
label will be excluded from calculations. To define the
material’s properties (thereby including the blocks into
calculation), type in the required values of Young's moduli.
As long as Anisotropic is unchecked, only two of the values
entered on this dialog page remain independent. As you type
in a new value, QuickField automatically updates the rest of
them. To describe orthotropic materials with seven
independent values, check Anisotropic before entering the
required values.
122 Chapter 5 Problem Parameters Description
2. Loads
3. Allowable Stresses
Editing Curves
Curves in QuickField represent dependencies between two
physical properties, e.g., magnetic field intensity and flux
density or temperature and thermal conductivity. To define
a curve in QuickField, you open this curve’s editor and enter
coordinates of several points on its graph. QuickField
accumulates the points in a table and provides you with
graphical representation of the curve interpolating between
the table points with cubic splines. QuickField solver always
uses the interpolated values displayed on screen during the
edit session.
Temperature Units
In QuickField Options it is possible to choose the
temperature units (Kelvins, Celsius or Fahrenheit degrees)
which then will be used everywhere in QuickField for this
user for all QuickField problems. These units will be used
for display of numerical data or plots, but they do not affect
128 Chapter 5 Problem Parameters Description
Using Formulas
Solving a problem you might need to specify a boundary
condition or a field source as a function of time or
coordinates. To do that, QuickField allows entering
formula-defined field values instead of constant numerical
values. The fields accepting formula-defined values are
specified above in this chapter. Besides that, the fields
accepting formula-defined values can be distinguished by
adjacent button and by the appropriate tooltip text.
Formula in QuickField is a mathematical expression
constructed of numbers, arithmetical operators, parentheses,
built-in constants and functions and predefined variables.
Formula syntax is typical for most algorithmic languages
and standard mathematical notation.
When you specify a formula-defined value QuickField
checks the formula syntax and reports syntax errors to you.
If the syntax is correct, calculator tries to calculate the result
using current values of predefined variables. This
calculation might also result in error if, for example, the
value of a function argument does not belong to the
function’s domain of definition. As with syntax errors,
QuickField reports such errors to you.
130 Chapter 5 Problem Parameters Description
• Unary operators
+ sign retaining (Ex: +180)
- sign change (Ex: -180)
• Embedded functions
abs - absolute value
sign - sign
max - maximum
min - minimum
step - step-function by 1
impulse - impulse segment
sin - sine
cos - cosine
tan - tangent
asin - arc sine
acos - arc cosine
atan - arc tangent
atan2 - angle of the vector defined by two
arguments
exp - exponent
log - natural logarithm
sqrt - square root
pow - raising to a power
saw - saw-tooth periodic function
• Embedded constants
pi - pi
e -e
• Predefined variables
t - current time
x - Cartesian coordinate x
y - Cartesian coordinate y
r - polar coordinate r
136 Chapter 5 Problem Parameters Description
Notes:
1. Numerical values should not contain group separators.
Use dots ('.') as decimal separators regardless of regional
settings.
2. QuickField allows using both e and E to separate
mantissa and magnitude in floating point values.
3. The names of embedded functions, constants and
predefined variables are case insensitive.
4. QuickField allows enclosing of names in double quotes.
For example, sin (t) is equivalent to "sin"(t).
5. Operation precedence (highest to lowest): ^, then * and
/, then + and -. Use parentheses, if you need to change this
order.
6. QuickField allows to insert any number of spaces inside
a formula without impact, provided the inserted spaces are
not inside names.
7. Place function arguments inside parentheses after the
name of the function and separate them with commas (',').
Constants and Predefined Variables
pi 3.141592653589793238462643 The constant equal
to the ratio of the
length of any circle
to its diameter.
Editing Label Data 137
= {1, 𝑖𝑓 𝑙 ≤ 𝑣 ≤ 𝑟 v
0, 𝑖𝑓 𝑣 > 𝑟 l r
saw(,p,p0) = saw(v,p,p0)
𝑣 1
, 𝑖𝑓 0 ≤ 𝑣 < 𝑝
𝑝
p p0 p p0 v
0, 𝑖𝑓 𝑝 ≤ 𝑣 < 𝑝 + 𝑝0
𝑠𝑎𝑤(𝑣 + (𝑝 + 𝑝0)),
𝑖𝑓 𝑣 < 0
𝑠𝑎𝑤(𝑣 − (𝑝 + 𝑝0)),
{ 𝑖𝑓 𝑣 ≥ 𝑝 + 𝑝0
142 Chapter 5 Problem Parameters Description
Examples
The table below contains examples you can use to learn
writing your own QuickField formulas. The left column
contains mathematical formulas with the corresponding
QuickField expressions contained in its right counterpart.
You can evaluate expressions and view a plot using our
online utility
www.quickfield.com/qf_formula_plotter.htm
Mathematical notation Formula syntax
2t2 - t - 3 2*t^2 - t - 3
e-t²/2 exp(-t^2 / 2)
{
𝑡, 𝑖𝑓 𝑡 < 0.5 t*step(0.5-t) + (1-t)*step(t-
1 − 𝑡, 𝑖𝑓 𝑡 ≥ 0.5 0.5)
0, 𝑖𝑓 𝑡 < 0
{
𝑡, 𝑖𝑓 0 ≤ 𝑡 < 0.5 t*impulse(t,0,0.5) + (1-
1 − 𝑡, 𝑖𝑓 0.5 ≤ 𝑡 < 1 t)*impulse(t,0.5,1)
0, 𝑖𝑓 𝑡 ≥ 1
sin 𝑡 , 𝑖𝑓 sin 𝑡 > cos 𝑡
{
cos 𝑡 , 𝑖𝑓 sin 𝑡 ≤ cos 𝑡 max(sin(t), cos(t))
𝑡
, 𝑖𝑓 0 ≤ 𝑡 < 2
{ 2 saw(t, 2)
periodic with period 2
10 ∙ 𝑒 5𝑡 , 𝑖𝑓 0 ≤ 𝑡 < 2
{ 10, 𝑖𝑓 2 ≤ 𝑡 < 3 10 * exp(5 * saw(t,2,1))
periodic with period 3
10 ∙ 𝑒 5𝑡 , 𝑖𝑓 0 ≤ 𝑡 < 2
{ 0, 𝑖𝑓 2 ≤ 𝑡 < 3
10 * exp(5 * saw(t,3)) *
periodic with period 3 impulse(saw(t,3), 0, 2/3)
𝑒 𝑡−1 , 𝑖𝑓 0 ≤ 𝑡 < 1
{ 𝑒 1−𝑡 , 𝑖𝑓 1 ≤ 𝑡 < 2
exp(saw(t,1,1)-1) +
periodic with period 2 exp(saw(2-t,1,1)-1) - exp(-1)
To copy a label:
1. In the list of labels, choose Copy from the
label’s context menu.
2. Switch to destination window and click Paste
in the Edit menu or context menu.
Or,
1. Drag the label to destination position with left
mouse button.
To delete a label:
• In the list of labels, choose Delete from the label’s
context menu, or
• Select the label and click Delete in the Edit menu.
To move (cut and paste) a label:
1. In the list of labels, choose Cut from the label’s
context menu.
2. Switch to destination window and click Paste
in the Edit menu or context menu.
Or,
1. Holding SHIFT pressed, drag the label to
destination position with left mouse button.
145
C H A P T E R 6
What is a Circuit?
Electric circuit consists of circuit components connected
with wires. Circuit components in QuickField can be of two
kinds:
1. First group includes usual electric circuitry
components, such as:
• resistors,
• capacitors,
• inductors,
• voltage sources,
146 Chapter 6 Electric Circuit Definition
• current sources.
2. Second group is specific for QuickField and
represents blocks of the geometric model. These
elements are used to provide interaction between
circuit and other parts of QuickField problem.
Notes:
• The point where you click the mouse button will be the
left pin of the new device.
• Circuit components are always aligned to a nearest grid
point. It means that component will be placed so that its
pins are at grid points.
• To insert component in the middle of some wire you
can just click on this wire. Wire segment will be split in
two and component will be inserted between them.
Notes:
• This way you may add either the vertical wire segment,
or horizontal wire segment, or two wire segments, vertical
and horizontal, making the right angle. To create wire of
more complex shape you should repeat this operation
several times.
• Wires are always aligned to a nearest grid node. It
means that the wire will be placed so that its end points are
at grid nodes.
Adding Junction Points
To add a junction point, you should just place a wire so that
one of end points belongs to some existing wire. Junction
point will be added in this point automatically.
Placing a wire so that it intersects another wire in a middle
point for both wires is different. In this case the junction
point will not be added, and the wires will be considered as
not connected.
Editing Circuit
Moving, Copying and Resizing Circuit Elements
To move circuit element to another place:
Editing Circuit 151
press alt and keep it pressed until you release the mouse
button.
Copying Elements
Instead of moving elements, you can make a copy of
selected elements. For this, press ctrl and keep it pressed
until you release the mouse button.
C H A P T E R 7
C H A P T E R 8
Analyzing Solution
be shown as:
• momentary value at a given phase 0 = –t0
• z = Re[z0ei(z – 0)] = z0 cos(z – 0);
0
1 2 3
𝑄= 𝑗 + 𝑘ℎ ∙ 𝑓 ∙ 𝐵2 + 𝑘𝑐 ∙ 𝑓 2 ∙ 𝐵2 + 𝑘𝑒 ∙ (𝑓 ∙ 𝐵)2
𝜎
• Time average and peak magnetic field energy density
w = (BH)/2;
• Time average Poynting vector (local power flow)
S = [ExH];
• Time average Lorentz force density vector F = [jxB];
• Magnetic permeability (its largest component in
anisotropic media);
• Electric conductivity σ.
Besides that, if the solved problem is coupled with the
electric circuit, then the following parameters can be
displayed in the circuit window:
• Effective I, amplitude Iabs, momentary (for the chosen
phase) and complex (Ire, Iim) values of the current in the
circuit branches;
• Effective U, amplitude Uabs, momentary (for the chosen
phase) and complex (Ure, Uim) values of the voltage drop
across the circuit component;
• Defined parameter value (resistance R, inductance L,
capacitance C) for the passive circuit components.
For the electrostatic problem:
• Scalar electric potential (voltage) U;
• Vector of electric field intensity E = −grad U;
• Tensor of gradient of electric field G = grad E;
• Vector of electrostatic induction D = E;
• Electric permittivity (or its largest component in
anisotropic media);
• Electric field energy density w = (ED)/2.
Building the Field Picture on the Screen 163
be shown as:
• momentary value at a given phase 0 = –t0
• z = Re[z0ei(z – 0)] = z0 cos(z – 0);
0
1
σ𝑒 = √ [(𝜎1 − 𝜎2 )2 + (𝜎2 − 𝜎3 )2 + (𝜎3 − 𝜎1 )2 ]
2
where 1, 2 and 3 denote the principal stresses in
descending order.
• Tresca criterion (maximum shear): e = 1 - 3;
• Mohr-Coulomb criterion: e = 1 - 3,
where = [+]/[-], and [+] and [-] denote tensile and
compressive allowable stress.
• Drucker-Prager criterion:
2
√𝜒 − 𝜒 1 1 − √𝜒
𝜎𝑒 = (1 + √𝜒)𝜎𝑖 − 𝜎̅ + ( ∙ 𝜎̅)
1 + √𝜒 [𝜎− ] 1 + √𝜒
where
166 Chapter 8 Analyzing Solution
1
σ𝑖 = √ [(𝜎1 − 𝜎2 )2 + (𝜎2 − 𝜎3 )2 + (𝜎3 − 𝜎1 )2 ]; 𝜎̅
2
𝜎1 + 𝜎2 + 𝜎3
=
3
• Tsai-Hill failure index for orthotropic materials:
𝜎12 𝜎1 𝜎2 𝜎22 𝜏12
2
𝐶𝑡ℎ = 2− 2 + 2+ 2
𝑋1 𝑋1 𝑋2 𝑆12
where 1, 2 and 12 are computed stresses in the material
directions and,
X1 = X1T if 1>0; X1 = X1C if 1<0
X2 = X2T if 2>0; X2 = X2C if 2<0
S12 = S12+ if 12>0; S12 = S12- if 12<0,
where X1T, X2T, X1C, X2C, S12+ and S12- are tensile,
compressive and shear allowable stresses.
The Tsai-Hill failure index is calculated only for those
materials, where allowable stresses were defined (while
editing the block data, see “Problem Parameters
Description”). If any pair of allowable stresses is not
given, the corresponding term is dropped while
calculating the Tsai-Hill Index.
The Field View dialog box for the stress analysis problem
additionally allows to select tensor quantity visualization.
Building the Field Picture on the Screen 171
Zooming
Zooming in postprocessor view is very similar to Model
Editor.
To magnify the picture:
1. Click Zoom In on the toolbar
2. Select (click and drag diagonally) the
rectangular part of the picture to fill the whole
window with.
To see more of the model:
• Click Zoom Out on the toolbar; or
• Click Zoom to Fit to see the whole model.
Animation
When the analysis results depend on time (transient and
time-harmonic problems), QuickField can present animated
field pictures. To start animation, choose View / Animation
or click the corresponding toolbar button .
During the animation, two speed control buttons appear on
the toolbar allowing you to change the animation speed. The
moment of time or the phase on the toolbar changes
synchronously with the picture.
During the first loop, QuickField accumulates the animation
frames in memory. Depending on the size of the problem,
this can take considerable time, and the speed control has no
effect.
Animation stops automatically upon any user action that
changes the contents or the scaling of the field picture.
174 Chapter 8 Analyzing Solution
Calculator Window
Calculator Window is a window normally docked to the left
side of the field view.
To open the calculator window, choose
Calculator Window command in the View menu or
corresponded button on the postprocessor toolbar. The
calculator window also opens when choosing Local Values,
Integral Values or one of the Wizard commands in View
menu.
The calculator window is organized in several trees, which
root items correspond to several kinds of numerical data.
These are:
• Local Values shows several field quantities at a point of
interest;
• Integral calculator lists available quantities calculated
by integration over given line, surface or volume;
• Inductance Wizard opens wizard, which helps you
calculate self or mutual inductance of the coils and
conductors;
• Capacitance Wizard opens wizard guiding you through
steps needed to calculate self or mutual capacitance of
your conductors in electrostatics problems;
• Impedance Wizard opens wizard, which helps you
calculate the impedance of the conductors in AC
magnetic problems.
Examining Local Field Data 175
You can resize the panes. You can also hide panes if there
is not enough of window space for all of them. Twin
boundary frames indicate the existence of hidden panes in
between.
The plots are displayed curves are plotted for each selected
circuit component. The correspondence between the curves
and circuit elements is shown by using the same color of a
curve and related line in the circuit element list. Curves are
redrawn when the circuit elements are selected or deselected
in the list or in the electric schema drawing.
The scale of XY-plot could be adjusted by corresponding
toolbar buttons: . Pressing the zoom-in button (with
180 Chapter 8 Analyzing Solution
plus sign) changes the cursor shape to a cross, after that you
can drag over the rectangle of interest.
The graph image could be printed, copied to the Clipboard,
and/or saved to a file in any of the supported raster or vector
formats. To do that, invoke the corresponding commands in
the context (right-click) menu of the graph.
The table of circuit currents and voltages could be printed,
copied to the Clipboard, and/or saved to a text file. To do
that, invoke the corresponding commands in the context
(right-click) menu of the table.
Inductance Wizard
Inductance wizard helps you to calculate self and mutual
inductance of your coils in the problem of DC or AC
magnetics.
When your model contains several coils that carry different
currents, the flux linkage with one of them can be calculated
as
𝜙𝑘 = 𝐿𝑘𝑘 ∙ 𝑖𝑘 + ∑ 𝑀𝑛𝑘 ∙ 𝑖𝑛
𝑛
where Lkk is the self inductance of the coil k, Mnk is the
mutual inductance between the coils n and k, ik is the current
in the coil k.
On the other hand, stored magnetic energy also derives from
current and inductance:
1
𝑊= (∑ 𝐿𝑘𝑘 ∙ 𝑖𝑘2 + ∑ 𝑀𝑛𝑘 ∙ 𝑖𝑛 ∙ 𝑖𝑘 )
2
𝑘 𝑛≠𝑘
Before using the inductance wizard, you have to formulate
your problem in such a way that all the currents (space,
surface or linearly distributed) but one are set to zero. There
must be no permanent magnets in your model. In that case
182 Chapter 8 Analyzing Solution
∫𝐿 𝐴 𝑑𝑠 ∫𝑅 𝐴 𝑑𝑠
𝜙 =𝑁∙( − ) ⎯ for planar case;
∫𝐿 𝑑𝑠 ∫𝑅 𝑑𝑠
∫𝑅 𝑟𝐴 𝑑𝑠 ∫𝐿 𝑟𝐴 𝑑𝑠 ⎯ for axisymmetric
𝜙 = 2π𝑁 ∙ ( − )
∫𝑅 𝑑𝑠 ∫𝐿 𝑑𝑠 case;
Capacitance Wizard
Capacitance wizard helps you to calculate self and mutual
capacitance of your conductors.
When your model contains several conductors, the charge
of one of them can be calculated as:
𝑞𝑘 = 𝐶𝑘𝑘 ∙ 𝑈𝑘 + ∑ 𝐶𝑛𝑘 ∙ 𝑈𝑛
𝑛
where Ckk is the self-capacitance of the conductor k, Cnk is
the mutual capacitance between the conductors n and k, Uk
is the voltage drop on the conductor k.
On the other hand, stored energy also derives from charge
and capacitance as:
1 𝑞𝑘2 𝑞𝑛 ∙ 𝑞𝑘
𝑊 = (∑ +∑ )
2 𝐶𝑘𝑘 𝐶𝑛𝑘
𝑘 𝑛≠𝑘
and from the voltage and capacitance as:
1
𝑊= (∑ 𝐶𝑘𝑘 ∙ 𝑈𝑘2 + ∑ 𝐶𝑛𝑘 ∙ 𝑈𝑛 ∙ 𝑈𝑘 )
2
𝑘 𝑛≠𝑘
Before using the capacitance wizard, you have to formulate
your problem in such a way that all field sources (space,
surface or linear distributed charge or voltage) but one are
set to zero. In that case equation above becomes extremely
simple and you can get capacitance value if you know any
two of these three quantities: charge, voltage, stored energy.
When formulating your problem, you can apply known
voltage to the conductor and measure the charge it produces
or vice versa. Measuring the charge is a bit more complex
186 Chapter 8 Analyzing Solution
Impedance Wizard
Impedance wizard helps you to calculate the impedance of
your conductors. It is simple and contains only one page. To
get the impedance value and its real and imaginary parts
(resistance and reactance accordingly) the impedance
wizard simply divides complex values of voltage by current:
Z = U/I, R = Re[Z], XL = Im[Z], L = XL / 2f,
188 Chapter 8 Analyzing Solution
Editing Contours
The contour is a directed curved line consisting of line
segments and arcs (including the edges of the model). Some
rules are applied to the contours:
• The contour may not intersect itself.
• Open and closed contours are discerned.
• Multiply connected contours have sense only for
calculating integral quantities.
Editing Contours 189
X-Y Plots
QuickField postprocessor can display field distribution
along contours. To open new X-Y plot window, choose X-Y
plot in View menu or context (right mouse button) menu in
field picture window, in which the contour is already
defined.
In X-Y plot view, you can:
X-Y Plots 191
Calculating Integrals
QuickField calculates line, surface and volume integrals. In
plane-parallel problem, a contour defines cylindrical (in
generalized sense) surface of infinite depth, or volume of
that cylinder for volume integral. Therefore, in
plane-parallel formulation surface and volume integrals are
calculated per unit depth. In axisymmetric problem, a
contour defines toroidal surface, or toroid for volume
integral.
Positive direction of a contour is counter-clockwise. The
direction of the contour is accounted as follows:
Calculating Integrals 193
Name,
Formula and Description
ActiveFied constant
Total current
𝐼𝑡𝑜𝑡𝑎𝑙 = ∫ 𝑗𝑡𝑜𝑡𝑎𝑙 𝑑𝑠
qfInt_Jtotal 𝑆𝑐
Complex value.
Electric current through a
particular surface.
External current
𝐼𝑒𝑥𝑡 = ∫ 𝑗𝑒𝑥𝑡 𝑑𝑠
qfInt_Jextern 𝑆𝑐
Complex value.
External current through a
particular surface.
Eddy current
𝐼𝑒𝑑𝑑𝑦 = ∫ 𝑗𝑒𝑑𝑑𝑦 𝑑𝑠
qfInt_Jeddies 𝑆𝑐
Complex value
Eddy current through a
particular surface.
Joule heat 1 2
𝑃 = ∫ ∙𝑗 𝑑𝑣
qfInt_Power 𝑉 𝜎 𝑡𝑜𝑡𝑎𝑙
Quadratic value.
Joule heat power in a particular
volume.
σ – electric conductivity of the
media.
Calculating Integrals 203
Core loss 𝑘ℎ ∙ 𝑓 ∙ 𝐵 2
2 2
qfInt_Steinmetz 𝑃 = ∫ [ +𝑘𝑐 ∙ 𝑓 ∙ 𝐵 ] 𝑑𝑣
𝑉 3
+𝑘𝑒 ∙ (𝑓 ∙ 𝐵)2
where B is the magnetic flux
density peak value,
and kh, kc, ke are user defined
coefficients.
Power flow
𝑃𝑆 = ∫ (𝐒 ∙ 𝐧)𝑑𝑠
qfInt_EnergyFlow 𝑆
Quadratic value.
Power flow through the given
surface (Poynting vector flow)
Here S – is a Pointing vector S =
[EH].
Maxwell force 𝐇(𝐧𝐁)
1
qfInt_MaxwellForce 𝐅 = ∮ (+ 𝐁(𝐧𝐇)) 𝑑𝑠,
2
𝑆 – 𝐧(𝐇𝐁)
Quadratic vector.
Maxwell force acting on bodies
contained in a particular
volume.
The integral is evaluated over
the boundary of the volume, and
n denotes the vector of the
outward unit normal.
204 Chapter 8 Analyzing Solution
Maxwell torque 𝐓
qfInt_MaxwellTorque [𝐫 × 𝐇](𝐧𝐁)
1
= ∮ ( + [𝐫 × 𝐁](𝐧𝐇)) 𝑑𝑠
2
𝑆 – [𝐫 × 𝐧](𝐇𝐁)
Quadratic value.
Maxwell force torque acting on
bodies contained in a particular
volume, where r is a radius
vector of the point of
integration.
The torque vector is parallel to z-
axis in the planar case, and is
identically equal to zero in the
axisymmetric one. The torque is
considered relative to the origin
of the coordinate system. The
torque relative to any other
arbitrary point can be obtained
by adding extra term of [F×r0],
where F is the total force and r0
is the radius vector of the point.
Lorentz force
qfInt_LorentzForce 𝐅 = ∫ [𝐣 × 𝐁]𝑑𝑣.
𝑉
Quadratic vector.
The Lorentz force acting on
conductors contained in a
particular volume.
Calculating Integrals 205
Lorentz torque
qfInt_LorentzTorque 𝐓 = ∫ [𝐫 × [𝐣 × 𝐁]]𝑑𝑣.
𝑉
Quadratic value.
The Lorentz force torque acting
on bodies contained in a
particular volume. The torque is
considered relative to the origin
of the coordinate system.
Magnetic field energy
1
qfInt_MagneticEnergy 𝑊= ∫ (𝐇 · 𝐁)𝑑𝑣
2
𝑉
Quadratic value.
This formula is used for both
linear and nonlinear cases.
Flux linkage per one ∮ 𝐴𝑑𝑠 ⎯ for planar
turn Ψ =
𝑆𝐶 case;
qfInt_FluxLinkage
Ψ ⎯ for
2𝜋 ∮(𝑟𝐴)𝑑𝑠 axisymmetric
=
𝑆𝐶 case.
Complex value.
The integral has to be evaluated
over a cross section of the coil,
and SC is the area of the cross
section.
206 Chapter 8 Analyzing Solution
Magnetomotive force
𝐹 = ∫ (𝐇 ∙ 𝐭)𝑑𝑙
qfInt_KGrad_t_dl 𝐿
Complex value.
Magnetomotive force.
Magnetic flux
𝛷 = ∫ (𝐁 ∙ 𝐧)𝑑𝑠
qfInt_Grad_n_ds 𝑆
Complex value.
Magnetic flux through a
particular surface.
Surface energy
1
qfInt_GradKGrad_n_ds 𝑊𝑠 = ∫ (𝐁 · 𝐇)𝑑𝑠
2
𝑆
Quadratic value.
The integral is evaluated over
the surface swept by the
movement of the contour.
Average surface 1
potential 𝐴𝑆 = ∫ 𝐴𝑑𝑠
𝑆 𝑆
qfInt_Potential_ds Complex value.
Integrated over a particular
volume
Average volume 1
potential 𝐴𝑉 = ∫ 𝐴𝑑𝑣
𝑉 𝑉
qfInt_Potential_dv Complex value.
Calculating Integrals 207
Potential difference
qfInt_Grad_t_dl ∆𝑈 = ∫ (𝐄 · 𝐭)𝑑𝑙
𝐿
The potential difference
between the ending and
started points of a contour can
be calculated as a line integral
over the contour of electric
field strength.
Average surface 1
potential 𝑈𝑆 = ∫ 𝑈𝑑𝑠
𝑆 𝑆
qfInt_Potential_ds
Average volume 1
potential 𝑈𝑉 = ∫ 𝑈𝑑𝑣
𝑉 𝑉
qfInt_Potential_dv
Calculating Integrals 211
Average surface 1
potential 𝑈𝑆 = ∫ 𝑈𝑑𝑠
𝑆 𝑆
qfInt_Potential_ds
Average volume 1
potential 𝑈𝑉 = ∫ 𝑈𝑑𝑣
𝑉 𝑉
qfInt_Potential_dv
Average volume 1
strength 𝐄𝑎 = ∫ 𝐄𝑑𝑣
𝑉 𝑉
qfInt_Grad_dv
Average volume 1
current density 𝐣𝑎 = ∫ 𝐣𝑑𝑣
𝑉 𝑉
qfInt_KGrad_dv
Mean square strength 1
𝐸𝑎2 = ∫ 𝐸 2 𝑑𝑣
qfInt_Grad2_dv 𝑉 𝑉
Mean square current 1
density 𝑗𝑎2 = ∫ 𝑗 2 𝑑𝑣
𝑉 𝑉
qfInt_KGrad2_dv
Surface integral of
strength 𝑥 = ∫ (𝐄 ∙ 𝐧)𝑑𝑠
𝑆
qfInt_Grad_n_ds
214 Chapter 8 Analyzing Solution
Mechanical torque 𝐓
qfInt_MaxwellTorque [𝐫 × 𝐄](𝐧𝐃)
1
= ∮ ( + [𝐫 × 𝐃](𝐧𝐄) ) 𝑑𝑠
2
𝑆 – [𝐫 × 𝐧](𝐄𝐃)
Quadratic value.
Electric force torque acting on
bodies contained in a
particular volume, where r is
a radius vector of the point of
integration.
The torque vector is parallel to
z-axis in the planar case, and is
identically equal to zero in the
axisymmetric one. The torque
is considered relative to the
origin of the coordinate
system. The torque relative to
any other arbitrary point can be
obtained by adding extra term
of [F×r0], where F is the total
force and r0 is the radius vector
of the point.
Electric field energy 1
𝑊= ∫ (𝐄 ∙ 𝐃)𝑑𝑣
qfInt_ElectrostaticEnergy 2 𝑉
Quadratic value.
Electric field energy in a
particular volume.
218 Chapter 8 Analyzing Solution
Surface energy
𝑊𝑆 = ∫ (𝐄 ∙ 𝐃)𝑑𝑠
qfInt_GradKGrad_n_ds 𝑆
Quadratic value.
Potential difference
∆𝑈 = ∫ (𝐄 ∙ 𝐭)𝑑𝑙
qfInt_Grad_t_dl 𝐿
Complex value.
The potential difference
between the ending and
started points of a contour can
be calculated as a line integral
over the contour of electric
field strength.
Average surface potential 1
𝑈𝑆 = ∫ 𝑈𝑑𝑠
qfInt_Potential_ds 𝑆 𝑆
Complex value.
Average volume 1
potential 𝑈𝑉 = ∫ 𝑈𝑑𝑣
𝑉 𝑉
qfInt_Potential_dv Complex value.
Average volume strength 1
𝐄𝑎 = ∫ 𝐄𝑑𝑣
qfInt_Grad_dv 𝑉 𝑉
Complex vector.
Average electric field strength
vector in a particular volume.
Calculating Integrals 219
Average volume 1
displacement 𝐃𝑎 = ∫ 𝐃𝑑𝑣
𝑉 𝑉
qfInt_KGrad_dv Complex vector.
Average electric displacement
vector in a particular volume.
Mean square strength 1
𝐸𝑎2 = ∫ 𝐸 2 𝑑𝑣
qfInt_Grad2_dv 𝑉 𝑉
Quadratic value.
Mean square 1
displacement 𝐷𝑎2 = ∫ 𝐷2 𝑑𝑣
𝑉 𝑉
qfInt_KGrad2_dv Quadratic value.
Electric charge
𝑄𝑆 = ∫ (𝐃 ∙ 𝐧)𝑑𝑠
qfInt_KGrad_n_ds 𝑆
Complex value.
The total electric charge in a
particular volume can be
calculated as a flux of electric
displacement over the
volume’s closed boundary.
Line integral of
displacement 𝑥 = ∫ (𝐃 ∙ 𝐭)𝑑𝑙
𝐿
qfInt_KGrad_t_dl Complex value.
220 Chapter 8 Analyzing Solution
Potential difference
∆𝑈 = ∫ (𝐄 ∙ 𝐭)𝑑𝑙
qfInt_Grad_t_dl 𝐿
The potential difference between
the ending and started points of a
contour can be calculated as a
line integral over the contour of
electric field strength.
Average surface potential 1
𝑈𝑆 = ∫ 𝑈𝑑𝑠
qfInt_Potential_ds 𝑆 𝑆
Average volume 1
potential 𝑈𝑉 = ∫ 𝑈𝑑𝑣
𝑉 𝑉
qfInt_Potential_dv
Calculating Integrals 223
Average volume 1
temperature 𝐆𝑎 = ∫ 𝐆𝑑𝑣
𝑉 𝑉
gradient
Mean vector of temperature gradient
qfInt_Grad_dv in a volume.
Average volume 1
heat flux density 𝐅𝑎 = ∫ 𝐅𝑑𝑣
𝑉 𝑉
qfInt_KGrad_dv Mean vector of heat flux density in
a volume.
Average volume 1
temperature 𝐺𝑎2 = ∫ 𝐺 2 𝑑𝑣
𝑉 𝑉
gradient
qfInt_Grad2_dv
Mean square heat 1
flux density 𝐹𝑎2 = ∫ 𝐹 2 𝑑𝑣
𝑉 𝑉
qfInt_KGrad_dv
Line integral of
heat flux density 𝑥 = ∫ (𝐅 ∙ 𝐭)𝑑𝑙
𝐿
qfInt_KGrad_t_dl
Surface integral of
grad(T) 𝑥 = ∫ (𝐆 ∙ 𝐧)𝑑𝑠
𝑆
qfInt_Grad_n_ds
For problems of stress analysis:
For the stress analysis problems, the most interesting
integral values are: force, torque and lengthen.
226 Chapter 8 Analyzing Solution
Data Tables
QuickField can display the field data at discrete points,
distributed along the currently selected contour, in table
view. To open new table window, choose Table in the View
menu or context (right mouse button) menu in field picture
window, in which the contour is already defined.
In table view, you can:
• Select the list of shown quantities (table columns).
Choose Columns in View or context menu.
• Select how the points are distributed along the contour
(table rows). Choose Rows in View or context menu.
• Insert additional rows at specified distance from the
beginning of the contour. Choose Insert in Edit or
context menu.
• Copy the set of rows or the whole table to Windows
clipboard. In latter case (when all of the rows are
selected), column headers are also copied. To copy the
header only, click the right mouse button within the
header and choose Copy Header from the context menu.
228 Chapter 8 Analyzing Solution
Table Columns
To change the set of visible table columns or the contents of
their headers, choose Columns from the View menu or from
context (right mouse button) menu in the Table window.
The Table Columns dialog appears on the screen:
The left part of the dialog
window contains the list of
all known columns. You
choose the set of visible
columns marking or
clearing the check boxes in
the list. Clicking on
Select All or Clear All
button respectively marks
and clears all check boxes in
the list.
The state of the check boxes
in Header Contains
defines which column
identification details should
be included in its header.
The
Adjust Columns Width
radio button state specifies
whether the width of a
column should be based on
the width of its data or on
the width of its header.
Data Tables 229
Table Rows
To fill the table with rows of values calculated at points on
a contour, choose Rows from the View menu or from
context (right mouse button) menu in the Table window.
The Table Rows dialog appears on the screen:
The dialog provides two alternative
ways to tabulate the contour. You
either use a fixed step between
contour points or fill the required
number of rows.
To use fixed step, enter the step
length in Tabulating Step.
If you mark
Always Include Segment Ends,
QuickField adds ends of all contour
segments to the table starting a new
step at every segment end. When you
enter a value in Tabulating Step,
QuickField shows the total number of
steps in Rows to Fill.
230 Chapter 8 Analyzing Solution
Time Plot
For transient problems QuickField provides the means to
plot time-dependent field quantities versus time. You can
display both the plots of local field values at the given points
and the plots of the values integrated over the current
contour.
Plots and Tables versus Time 231
Alternately, you can invoke the Time Plot command via the
context menu of the Field Picture window and the
Time Plot Curves command via the context menu of the
Time Plot window. The context menu also provides the
way to easily switch the displayed time plot between
different groups of local values or move from local to
integrated values and back.
With time plot view, you can:
• Define the set of curves displayed for various groups of
physical values. To do it, choose Time Plot Curves
from the View or context menu.
• Zoom the plot in or out with toolbar buttons .
• View the legend showing the correspondence between
quantities and curves.
• Copy the picture to the clipboard with
Edit / Copy Picture or save it to a file with
File / Save As.
𝑑𝑙 𝑑𝑥 2 𝑑𝑦 2 𝑑𝑧 2
√
= ( ) +( ) +( ) ,
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡
defining the length l(t) of the trajectory covered by the
particle in time t. We integrate the resulting system using
the Runge-Kutta-Merson method with automatically
defined integration step. Numerical integration stops
immediately before the finite element’s boundary, the step
leading outside of the element being excluded. At the last
point in the element, we extrapolate the trajectory with
cubical segment of its Taylor series relative to time and
solve the resulting equation using Tartaglia-Cardano
formula and taking into account possible decrease of the
equation’s degree in homogeneous or zero fields.
Using Trajectories
Choosing Particle Trajectory from View menu opens the
modeless dialog window that allows changing beam
properties and viewing the calculation results on screen.
238 Chapter 8 Analyzing Solution
The property values entered into the dialog fields come into
effect when you click Apply button.
The Particle dialog page allows to choose particle type from
the list or define its charge and mass manually. It also
provides you with the possibility to define the initial
absolute value of the particle’s velocity or its initial kinetic
energy.
Next tab Emitter to specify parameters of monochromatic
point emitter.
Trajectories of Charged Particles 239
C H A P T E R 9
Introduction
This section includes detailed information about three-
dimensional analysis with QuickField. The workflow and
technique of working with 3D problems mostly relies on
knowledge of 2D analysis. If you feel yourself not quite
familiar with such things as material properties, boundary
conditions, field sources, geometric labels and spacing and
so on, we recommend first browse through Chapters 2, 3, 4
and 5 of this manual as well as videos, samples, papers and
other teaching material available online at
www.quickfield.com.
Workflow for 3D model definition and analysis in general
follows the same steps as in 2D. You should:
• Create the new QuickField problem, and assign the
model type 3D Extrusion or 3D Import.
• For the problem with model class 3D Extrusion describe
the geometry in familiar Model Editor – first in XY
plane, and when add the heights to all planar objects,
or
Introduction 251
2D and 3D problems
The choice between plane-parallel, axisymmetric or 3D
formulation is performed by selecting the proper model
class in the QuickField problem properties.
2D problems
Majority of QuickField analysis types are formulated as 2D,
i.e. plane-parallel or axisymmetric approximations.
Plane-parallel problems
Plane-parallel y
approximation means that
the geometry, boundary model
x
conditions, field sources,
material properties and thus the calculated field depend on
X and Y coordinates, but do not depend on the coordinate Z.
254 Chapter 9 Working with 3D Problem
r
Axisymmetric problems
Axisymmetric model
z
approximation means that
the geometry, boundary
conditions, field sources, material properties and thus the
calculated field depend of Z and R, but do not depend on the
angular coordinate θ.
It means that all the half cross-sections, which include the
rotational axis X and positioned on one side from it, are the
same (i.e. independent from the angular coordinate θ).
In the screen plane we work with upper half of the axial
cross section of the model
3D problems
QuickField is able to solve 3D Electrostatic, 3D DC
conduction and 3D Steady state heat transfer analyses only.
In 3D problems the model geometry, material properties,
field sources and/or boundary conditions may depend upon
all three spatial coordinates X, Y and Z.
In QuickField the 3D geometry is built by one of two ways:
• By extrusion objects from the flat XY geometry into Z
direction. Each planar object may be extruded for
model x
z
Problem definition – adding 3D features 255
Switching to 3D View
After the height levels are set, the extrusion may be
performed by choosing a command 3D in the View menu,
or by pressing the button (Toggle 2D/3D in the model
editor toolbar.
For successful extrusion following conditions should be
met:
a. At least one block of the 2D model should
have the finite-element mesh and two or more
height levels assigned.
b. None of 2D model meshed blocks should have
one height level assigned.
If either of these conditions is not met, error message is
displayed, and if it was caused by “b” then the blocks where
the condition does not work are highlighted. In case of
success the bodies extruded from the blocks which satisfy
“a”, and lower dimension objects at their boundaries are
displayed in 3D.
this case with the extension .m3d. Double click on the file
name in the problem tree will open the model editor window
in its 3D view. The button Import STEP in the 3D model
editor toolbar should be clicked next. The STEP file (usually
with the extension .step or .stp) should be selected and the
button Open clicked then.
Background Region
If required, the imported geometry entities may be
complemented by the background region.
The background region represents the rectangular box,
oriented on the coordinate axes. The background region
dimensions are automatically chosen to include all the
geometrical objects with the predefined gaps around. The
gap is defined by the user in percent from the maximal
model size.
Extra vertices are usually used to fine tune the mesh density
distribution. To achieve this the needed mesh step may be
assigned to the extra vertex. Extra vertex may also have a
label and a boundary condition or field source associated
with this label.
Selected vertex
Selected edge
Selected face
Body “column”
Object selection
3D geometric objects: bodies, faces, edges, vertices are
displayed in the 3D window. Each of objects may be in the
normal or selected state. Object selection is needed to
change its properties (assign the label etc.) or hide the object
to see what is hidden behind it.
272 Chapter 9 Working with 3D Problem
Show spacing
on vertices
Solving 3D problems
To obtain the problem solution, click the Solve button on
the toolbar, or choose Solve Problem in the Problem menu.
You may skip this action and directly proceed to the analysis
results by clicking Analyze Results on the toolbar or in the
Problem menu. If the problem has not been solved yet, or
its results are out of date, the solver will be invoked
automatically.
Result analysis
Problem solution in QuickField starts by choosing the
command Solve in the Problem menu, or by pressing the
button in the toolbar.
3D solution results may be presented in QuickField by
following methods:
• ways of the field picture display:
Different
• Color painting of the surfaces according to the value
of the chosen physical parameter,
• Geometry model presentation using the edges,
• Presenting the finite element mesh on the visible
surfaces,
Result analysis 281
from the View menu or the button from the field view
toolbar switches this presentation on or off.
color surfaces will be shown from one side of the plane and
within the plane. Cut plane does not affect other
presentations (scatter plot, vector plot, isosurface plot etc.).
The command Cut by plane from the View menu or the
button from the field view toolbar turn the cut plane
presentation on and off.
Tables
Field parameter distribution along the selected edges or
contour may be presented not only on a graph, but also as a
table. Table window is activated by pressing the button in
the 3D field view toolbar (table along the selected edges) or
in the 2D section toolbar (table along the contour).
Table rows are related to the geometry model nodes, evenly
distributed along the contour or the selected edges group.
Table columns are the physical parameter values.
Table columns displayed may be modified by pressing the
button in the table toolbar or in the context menu using
the Table columns dialogue.
Integral calculation
In addition to the local field values QuickField can calculate
integral parameters, such as electric charge, mechanical
force and electric field energy.
Integration region for the volume integrals includes the
selected bodies, and for the surface integrals – the selected
faces or the planar section of the 3D model.
Integration area of the surface integrals, calculated over the
planar section includes the cross-sections of all the visible
bodies which were cut by the plane. Turning their visibility
on/off will control the integration area. Integration results
Result analysis 307
𝑞 = ∮ 𝐃 ∙ 𝐝𝐬
𝑆
where integral is evaluated over the boundary of the
volume, and n denotes the vector of the outward unit
normal.
• Total electric force acting on bodies contained in a
particular volume
1
𝐅= ∮(𝐄(𝐧𝐃) + 𝐃(𝐧𝐄)– 𝐧(𝐄𝐃))𝑑𝑠,
2
𝑆
• Energy of electric field
w = 1/2 ∫ (𝐄 ⋅ 𝐃) 𝐝𝐯
𝑉
In the DC conduction problems the following integrals may
be calculated:
• Electric current through a particular surface
I = ∮ 𝐉 ∙ 𝐝𝐬
𝑆
308 Chapter 9 Working with 3D Problem
Ф = ∮ 𝜆 ∙ grad 𝑇 𝐝𝐬
𝑆
where integral is evaluated over the boundary of the
volume.
Integral calculation mode may be turned on by pressing the
button on the toolbar, or by selecting the pane Integral
Calculator from the field calculator dialog, or by using the
command Integral Values from the View menu. This mode
is indicated by the cursor shape.
All available by the moment integral values are calculated
by integration across the selected bodies (or their external
surfaces). All changes in this selection invoke automatic
calculation of the new integral values.
309
C H A P T E R 1 0
Add-ins
C H A P T E R 1 1
Theoretical Description
Magnetostatics
QuickField can solve both linear and nonlinear magnetic
problems. Magnetic field may be induced by the
concentrated or distributed currents, permanent magnets or
external magnetic fields.
The magnetic problem is formulated as the Poisson's
equation for vector magnetic potential A (B = curl A, B⎯
magnetic flux density vector). The flux density is assumed
to lie in the plane of model (xy or zr), while the vector of
electric current density j and the vector potential A are
orthogonal to it. Only jz and Az in planar or j and A in
Magnetostatics 317
Field Sources
The field sources can be specified in blocks, at the edges or
at the individual vertices of the model. Possible field sources
include space, surface and linear electric currents and
permanent magnets. The coercive force is chosen to be
primary characteristic for the permanent magnets.
A point source in the xy-plane describes a linear current in
out-of-plane direction. In axisymmetric case the point
source represents the current in a thin ring around the axis
of symmetry. Edge-bound source in the plane of model
represents a surface current in three-dimensional world. It is
specified by the Neumann boundary condition for the edge.
The space current is described either by the electric current
318 Chapter 11 Theoretical Description
Boundary Conditions
The following boundary conditions can be specified at
outward and inner boundaries of the region.
Dirichlet condition specifies a known value of vector
magnetic potential A0 at the vertex or at the edge of the
model. This boundary condition defines normal component
of the flux density vector. It is often used to specify
vanishing value of this component, for example at the axis
of symmetry or at the distant boundary. QuickField also
supports the Dirichlet condition with a function of
coordinates. This approach allows you to model a uniform
external field by specifying non-zero normal component of
the flux density at arbitrary straight boundary segment.
Magnetostatics 319
Permanent Magnets
Since the coercive force is considered in QuickField to be
the piecewise constant function, its contribution to the
equation is equivalent to surface currents which flow along
the surface of the permanent magnet in direction orthogonal
to the model plane. The density of such effective current is
equal to jump of the tangent component of the coercive force
across the magnet boundary. For example, rectangular
magnet with the coercive force Hc directed along x-axis can
be replaced by two oppositely directed currents at its upper
and lower surfaces. The current density at the upper edge is
numerically equal to Hc, and −Hc at the lower edge.
Therefore, the permanent magnet can be specified by either
coercive force or Neumann boundary conditions at its edges.
You can choose more convenient and obvious way in each
particular case.
Permanent magnet with nonlinear magnetic properties
needs some special consideration. Magnetic permeability is
assumed to be defined by the following equation:
Magnetostatics 321
𝐵
𝐵 = (𝐵) · (𝐻 + 𝐻𝑐 ); (𝐵) =
𝐻 + 𝐻𝑐
It must be pointed out that (B) dependence is different from
the analogous curve for the same material but without
permanent magnetism. If the real characteristic for the
magnet is not available for you, it is possible to use row
material curve as an approximation. If you use such
approximation and magnetic field value inside magnet is
much smaller than its coercive force, it is recommended to
replace the coercive force by the following effective value
1
𝐻𝑐′ = 𝐵
𝜇(𝐵𝑟 ) 𝑟
where Br is remanent induction.
1 𝜕(𝑟𝐴)
𝐵𝑧 =∙ , 𝐵𝑟 ⎯ for
𝑟 𝜕𝑟
𝜕𝐴 axisymmetric
=− case;
𝜕𝑧
• Vector of magnetic field intensity
1
𝐇 = 𝐁,
𝜇
where is the magnetic permeability tensor.
Integral quantities:
• Total magnetic force acting on bodies contained in a
particular volume
1
𝐅= ∮ (𝐇(𝐧𝐁) + 𝐁(𝐧𝐇)– 𝐧(𝐇𝐁))𝑑𝑠,
2
𝑆
where integral is evaluated over the boundary of the
volume, and n denotes the vector of the outward unit
normal.
• Total torque of magnetic forces acting on
bodies contained in a particular volume
1
𝐓= ∮ ([𝐫 × 𝐇](𝐧𝐁)
2
𝑆
+ [𝐫 × 𝐁](𝐧𝐇) – [𝐫 × 𝐧](𝐇𝐁))𝑑𝑠,
where r is a radius vector of the point of integration.
The torque vector is parallel to z-axis in the planar case,
and is identically equal to zero in the axisymmetric one.
The torque is considered relative to the origin of the
Magnetostatics 323
1 ⎯ linear
𝑊= ∫ (𝐇 · 𝐁)𝑑𝑣
2 case;
𝑉
𝐵
⎯ nonlinear
𝑊 = ∫ (∫ 𝐻(𝐵′)𝑑𝐵′) 𝑑𝑣
case.
0
Inductance Calculation
To get self inductance of a coil, leave the current on in this
coil only and make sure that all other currents are turned off.
After solving the problem go to the Postprocessor and obtain
flux linkage for the contour coinciding with the cross section
of the coil. Once you’ve done that, the inductance of the coil
can be obtained from the following equation:
𝑛Ψ
𝐿= ,
𝐼
where n is a number of turns in the coil, is a flux linkage,
j is a current per one turn of the coil.
Mutual inductance between two coils can be obtained in a
similar way. The only difference from the previous case is
that electric current has to be turned on in one coil, and the
flux linkage has to be evaluated over the cross section of
another.
𝑛2 Ψ2
𝐿12 = ,
𝐼1
In plane-parallel case every coil has to be represented by at
least two conductors with equal but opposite currents. In
some cases both conductors are modeled, in other cases only
one of two conductors is included in the model and the rest
is replaced by the boundary condition A = 0 at the plane of
symmetry. If the magnetic system is symmetric, the
inductance can be obtained based on the flux linkage for one
of the conductors only. The result has to be then multiplied
by a factor of two to account for the second conductor. If the
model is not symmetric, then the total inductance can be
obtained by adding up the analogous terms for each
Transient Magnetics 325
Transient Magnetics
Transient magnetic analysis is the generalized form of
computation of electric and magnetic field, induced by
direct or time-varying currents (alternating, impulse, etc.),
permanent magnets, or external magnetic fields, in linear or
nonlinear (ferromagnetic) media, and takes into account
eddy current (skin) effect in conductors of electric current.
The formulation is derived from Maxwell's equations for
vector magnetic potential A (B = curl A) and scalar electric
potential U (E = -grad U):
𝟏
𝐜𝐮𝐫𝐥 𝐜𝐮𝐫𝐥 𝐀 = 𝐣 + 𝐜𝐮𝐫𝐥 𝐇𝒄 ,
𝜇
𝜕𝐀
𝐣 = σ𝐄 = −σ − σ ∙ 𝐠𝐫𝐚𝐝𝑈
𝜕𝑡
where 1/μ is an inverse permeability tensor, and g is electric
conductivity. In accordance with the second equation,
vector j of the total current in a conductor can be considered
as a combination of a source current produced by the
external voltage and an eddy current induced by the time-
varying magnetic field
𝐣 = 𝐣0 + 𝐣eddy ,
where
𝐣0 = − σ ∙ 𝐠𝐫𝐚𝐝𝑈,
326 Chapter 11 Theoretical Description
𝜕𝐀
𝐣eddy = −σ
.
𝜕𝑡
If a field simulation is coupled with an electric circuit, the
branch equation for a conductor is:
𝑈 𝜕𝐀
𝐼= −σ∫ 𝑑𝑠
𝑅 𝜕𝑡
Ω
where U is the voltage difference between the two terminals
of the solid conductor, and R is the DC resistance of the
conductor.
The flux density is assumed to lie in the plane of model (xy
or zr), while the vector of electric current density j and the
vector potential A are orthogonal to it. Only jz and Az in
planar or j and A in axisymmetric case are not equal to
zero. We will denote them simply j and A. Finally, the
equation for planar case is
𝜕 1 𝜕𝐴 𝜕 1 𝜕𝐴 𝜕𝐀
( )+ ( )−σ
𝜕𝑥 𝜇𝑦 𝜕𝑥 𝜕𝑦 𝜇𝑥 𝜕𝑦 𝜕𝑡
𝜕𝐻𝑐𝑦 𝜕𝐻𝑐𝑥
= −𝑗0 + ( − );
𝜕𝑥 𝜕𝑦
and for axisymmetric case is
𝜕 1 𝜕(𝑟𝐴) 𝜕 1 𝜕𝐴 𝜕𝐀
( )+ ( )−σ
𝜕𝑟 𝑟𝜇𝑧 𝜕𝑟 𝜕𝑧 𝜇𝑟 𝜕𝑧 𝜕𝑡
𝜕𝐻𝑐𝑟 𝜕𝐻𝑐𝑧
= −𝑗0 + ( − ),
𝜕𝑧 𝜕𝑟
where components of magnetic permeability tensor x and
y (z and r), components of coercive force vector Hcx and
Hcy (Hcz and Hcr) are constants within each block of the
Transient Magnetics 327
Field Sources
The field sources can be specified in blocks, at the edges or
at the individual vertices of the model. Possible field sources
include space, surface and linear electric currents; voltages
applied to conductive areas; and permanent magnets. The
coercive force is chosen to be primary characteristic for the
permanent magnets.
A point source in the xy-plane describes a linear current in
out-of-plane direction. In axisymmetric case the point
source represents the current in a thin ring around the axis
of symmetry. Edge-bound source in the plane of model
represents a surface current in three-dimensional world. It is
specified by the Neumann boundary condition for the edge.
The space-distributed current is defined differently in areas,
where the eddy current effect is considered (non-zero
conductivity is specified) or not considered (conductivity is
set to zero). In latter case, the space current is described
328 Chapter 11 Theoretical Description
Boundary Conditions
The following boundary conditions can be specified at
outward and inner boundaries of the region.
Dirichlet condition specifies a known value of vector
magnetic potential A0 at the vertex or at the edge of the
model. This boundary condition defines normal component
of the flux density vector. It is often used to specify
vanishing value of this component, for example at the axis
of symmetry or at the distant boundary. QuickField also
supports the Dirichlet condition with a function of
coordinates. This approach allows you to model a uniform
external field by specifying non-zero normal component of
the flux density at arbitrary straight boundary segment.
H t+ - H t- = σ ⎯ at inner boundaries,
where Ht is a tangent component of magnetic field intensity,
"+" and "−" superscripts denote quantities to the left and to
the right side of the boundary and σ is a linear density of the
surface current. If value is zero, the boundary condition is
called homogeneous. This kind of boundary condition is
often used at an outward boundary of the region that is
formed by the plane of magnetic antisymmetry of the
problem (opposite sources in symmetrical geometry). The
homogeneous Neumann condition is the natural one, it is
assumed by default at all outward boundary parts where no
explicit boundary condition is specified.
Permanent Magnets
Since the coercive force is considered in QuickField to be
the piecewise constant function, its contribution to the
equation is equivalent to surface currents which flow along
the surface of the permanent magnet in direction orthogonal
to the model plane. The density of such effective current is
equal to jump of the tangent component of the coercive force
across the magnet boundary. For example, rectangular
magnet with the coercive force Hc directed along x-axis can
be replaced by two oppositely directed currents at its upper
and lower surfaces. The current density at the upper edge is
numerically equal to Hc, and −Hc at the lower edge.
Therefore, the permanent magnet can be specified by either
coercive force or Neumann boundary conditions at its edges.
You can choose more convenient and obvious way in each
particular case.
Permanent magnet with nonlinear magnetic properties
needs some special consideration. Magnetic permeability is
assumed to be defined by the following equation
𝐵
𝐵 = (𝐵) · (𝐻 + 𝐻𝑐 ); (𝐵) =
𝐻 + 𝐻𝑐
Pleae note that (B) dependence is different from the
analogous curve for the same material but without
permanent magnetism. If the real characteristic for the
magnet is not available for you, it is possible to use row
material curve as an approximation. If you use such
approximation and magnetic field value inside magnet is
much smaller than its coercive force, it is recommended to
replace the coercive force by the following effective value
332 Chapter 11 Theoretical Description
1
𝐻𝑐′ = 𝐵
𝜇(𝐵𝑟 ) 𝑟
where Br is remanent induction.
1
𝐇 = 𝐁,
𝜇
where is the magnetic permeability tensor.
• Joule heat density
1 2
𝑄 = 𝑗
𝜎
• Magnetic field energy density w = (BH)/2;
• Magnetic permeability (its largest
component in anisotropic media);
• Electric conductivity σ.
Integral quantities:
• Electric current through a particular surface
𝐼 = ∫ 𝑗 𝑑𝑠
1
𝐅= ∮ (𝐇(𝐧𝐁) + 𝐁(𝐧𝐇)– 𝐧(𝐇𝐁))𝑑𝑠,
2
𝑆
where integral is evaluated over the boundary of the
volume, and n denotes the vector of the outward unit
normal.
334 Chapter 11 Theoretical Description
1
𝐓= ∮ ([𝐫 × 𝐇](𝐧𝐁)
2
𝑆
+ [𝐫 × 𝐁](𝐧𝐇) – [𝐫 × 𝐧](𝐇𝐁))𝑑𝑠,
where r is a radius vector of the point of integration.
The torque vector is parallel to z-axis in the planar case,
and is identically equal to zero in the axisymmetric one.
The torque is considered relative to the origin of the
coordinate system. The torque relative to any other
arbitrary point can be obtained by adding extra term of
F r0, where F is the total force and r0 is the radius vector
of the point.
• Magnetic field energy
1 ⎯ linear
𝑊= ∫ (𝐇 · 𝐁)𝑑𝑣
2 case;
𝑉
𝐵
⎯ nonlinear
𝑊 = ∫ (∫ 𝐻(𝐵′)𝑑𝐵′) 𝑑𝑣
case.
0
2𝜋 ∮(𝑟𝐴)𝑑𝑠 ⎯ for
Ψ = axisymmetric
𝑆
case.
the integral has to be evaluated over the cross section of
the coil, and S is the area of the cross section.
For planar problems, all integral quantities are considered
per unit length in z-direction.
The domain of integration is specified in the plane of the
model as a closed contour consisting of line segments and
circular arcs.
AC Magnetic
AC magnetic analysis is the study of magnetic and electric
fields arising from the application of an alternating (AC)
current source, or an imposed alternating external field.
Variation of the field with respect to time is assumed to be
sinusoidal. All field components and electric currents vary
with time like
𝑧 = 𝑧0 cos(𝑡 + 𝑍 ),
where z0 is a peak value of z, z — its phase angle, and —
the angular frequency.
336 Chapter 11 Theoretical Description
Im
z0 t t+z Re
- 2 z
z
𝑈 𝜕𝐀
𝐼= −σ∫ 𝑑𝑠
𝑅 𝜕𝑡
Ω
where U is the voltage difference between the two terminals
of the solid conductor, and R is the DC resistance of the
conductor.
The problem is formulated as a partial differential equation
for the complex amplitude of vector magnetic potential A
(B = curl A, B⎯magnetic flux density vector). The flux
density is assumed to lie in the plane of model (xy or zr),
while the vector of electric current density j and the vector
potential A are orthogonal to it. Only jz and Az in planar or
j and A in axisymmetric case are not equal to zero. We will
denote them simply j and A. The equation for planar case is
𝜕 1 𝜕𝐴 𝜕 1 𝜕𝐴
( )+ ( ) − 𝑖𝜔𝜎𝐴 = −𝑗0
𝜕𝑥 𝜇𝑦 𝜕𝑥 𝜕𝑦 𝜇𝑥 𝜕𝑦
and for axisymmetric case is
𝜕 1 𝜕(𝑟𝐴) 𝜕 1 𝜕𝐴
( )+ ( ) − 𝑖𝜔𝜎𝐴 = −𝑗0
𝜕𝑟 𝑟𝜇𝑧 𝜕𝑟 𝜕𝑧 𝜇𝑟 𝜕𝑧
where electric conductivity σ and components of magnetic
permeability tensor x and y (z and r) are constants within
each block of the model. Source current density j0 is
assumed to be constant within each model block in planar
case and vary as ~(1/r) in axisymmetric case.
338 Chapter 11 Theoretical Description
Field Sources
The field sources can be specified in the blocks, at the edges
or at the individual vertices of the model. Possible field
sources include space, surface and linear electric currents
and voltages applied to conductive areas.
A point source in the xy-plane corresponds to a linear current
in out-of-plane direction. In axisymmetric case the point
source represents the current in a thin ring around the axis
of symmetry. Edge-bound source in the plane of model
AC Magnetic 339
j = n · I / S,
where n is a number of turns, I is a total current, and S is a
cross-sectional area of the coil.
Boundary Conditions
The following boundary conditions can be specified at
outward and inner boundaries of the region.
Dirichlet condition specifies a known value of vector
magnetic potential A0 at the vertex or at the edge of the
model. This boundary condition defines normal component
of the flux density vector. It is often used to specify
vanishing value of this component, for example at the axis
of symmetry or at the distant boundary. QuickField also
supports the Dirichlet condition with a function of
coordinates. This approach allows you to model a uniform
external field by specifying non-zero normal component of
the flux density at arbitrary straight boundary segment.
Neumann condition has the following form
Ht = σ ⎯ at outward boundaries,
H t+ - H t- = σ ⎯ at inner boundaries,
where Ht is a tangent component of magnetic field intensity,
"+" and "−" superscripts denote quantities to the left and to
the right side of the boundary and σ is a linear density of the
AC Magnetic 341
1 2 3
𝑄 = 𝑗 + 𝑘ℎ ∙ 𝑓 ∙ 𝐵2 + 𝑘𝑐 ∙ 𝑓 2 ∙ 𝐵2 + 𝑘𝑒 ∙ (𝑓 ∙ 𝐵)2
𝜎
• Time average and peak magnetic field energy
density w = B H / 2;
• Time average Poynting vector (local power
flow) S = E×H;
• Time average Lorentz force density vector
F = j×B;
• Magnetic permeability (its largest
component in anisotropic media);
• Electric conductivity σ.
Integral quantities:
• Complex magnitude of electric current through a
particular surface and its source and eddy components
𝐼 = ∫ 𝑗 𝑑𝑠
Oscillating
average component
peak
t
- 2
344 Chapter 11 Theoretical Description
𝑃𝑐𝑜𝑟𝑒 = ∫ [𝑘ℎ · 𝑓 · 𝐵2 + 𝑘𝑐 · 𝑓 2 · 𝐵2 + 𝑘𝑒
3
· (𝑓 · 𝐵)2 ] 𝑑𝑣,
where B – is the flux density magnetide and kh, kc, ke are
user defined coefficients.
• Time average and peak magnetic field energy
1
𝑊= ∫(𝐇 · 𝐁)𝑑𝑣
2
𝑉
• Time average and peak power flow through the
given surface (Poynting vector flow)
𝑆 = ∫ (𝐒 · 𝐧)𝑑𝑠
1
𝐅= ∮ (𝐇(𝐧𝐁) + 𝐁(𝐧𝐇)– 𝐧(𝐇𝐁))𝑑𝑠,
2
𝑆
where integral is evaluated over the boundary of the
volume, and n denotes the vector of the outward unit
normal.
• Time average and peak Maxwell force torque
acting on bodies contained in a particular
volume
AC Magnetic 345
1
𝐓= ∮ ([𝐫 × 𝐇](𝐧𝐁)
2
𝑆
+ [𝐫 × 𝐁](𝐧𝐇) – [𝐫 × 𝐧](𝐇𝐁))𝑑𝑠,
where r is a radius vector of the point of integration.
• Time average and oscillating part of Lorentz force acting
on conductors contained in a particular volume
𝐅 = ∫ [𝐣 × 𝐁]𝑑𝑣.
𝑉
• Time average and peak Lorentz force torque
acting on bodies contained in a particular
volume
𝐓 = ∫ [𝐫 × [𝐣 × 𝐁]]𝑑𝑣.
𝑉
where r is a radius vector of the point of integration.
The torque vector is parallel to z-axis in the planar case,
and is identically equal to zero in the axisymmetric one.
The torque is considered relative to the origin of the
coordinate system. The torque relative to any other
arbitrary point can be obtained by adding extra term of
F r0, where F is the total force and r0 is the radius vector
of the point.
346 Chapter 11 Theoretical Description
Impedance Calculation
Impedance in ac magnetic analysis is a complex coefficient
between complex values of current and voltage, V = ZI. Its
real part represents active resistance of the conductor,
calculated with the skin effect taken into account. The
imaginary part of the impedance is the inductance
multiplied by the angular frequency .
𝑍 = 𝑅 + 𝑖𝜔𝐿
As values of voltage and current in any conductor are easily
accessible in the postprocessor, you can determine the
impedance by dividing voltage by current using complex
arithmetic. Let V and I be peak values of voltage and current,
and V and I be phases of those quantities. Then the active
resistance is calculated as
𝑉
𝑅 = cos(φ𝑉 − φ𝐼 ),
𝐼
and the inductance as
𝑉
𝐿= sin(φ𝑉 − φ𝐼 ),
𝐼 ∙ 2𝜋𝑓
Electrostatics 347
Electrostatics
Electrostatic problems are described by the Poisson's
equation for scalar electric potential U (E = −gradU, E⎯
electric field intensity vector). The equation for planar case
is
𝜕 𝜕𝑈 𝜕 𝜕𝑈
(𝜀𝑥 ) + (𝜀𝑦 ) = −ρ,
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦
for axisymmetric case:
1𝜕 𝜕𝑈 𝜕 𝜕𝑈
(𝜀𝑟 𝑟 ) + (𝜀𝑧 ) = −ρ,
𝑟 𝜕𝑟 𝜕𝑟 𝜕𝑧 𝜕𝑧
And for 3D formulation:
𝜕 𝜕𝑈 𝜕 𝜕𝑈 𝜕 𝜕𝑈
(𝜀 ) + (𝜀 ) + (𝜀 ) = −ρ,
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧
where components of electric permittivity tensor x, y, z or
z, r and electric charge density are constants within each
block of the model. Anisotropy is not accounted for in 3D
formulation.
348 Chapter 11 Theoretical Description
Field Sources
QuickField provides possibility to specify electric charges
located in the bodies, on the faces, at the edges or at the
individual vertices of the model. In planar case, the electric
charge specified at a point of the xy-plane corresponds to a
charged string that is perpendicular to the plane of the
model, and is described by the linear charge density. In
axisymmetric case the vertex charge represents a charged
circle around the axis of symmetry or a point charge located
on the axis. To incorporate both these cases a total charge
value is associated with the vertex. For the charged circle
the total charge is connected with its linear density by the
relationship q = 2r. Edge-bound charge in the plane of
model represents a surface-bound charge in
three-dimensional world. It is described by surface charge
density and is specified by the Neumann boundary condition
for the edge. The charge density associated with a body or
within a block is equivalent to the space charge.
Boundary Conditions
The following boundary conditions can be specified at
outward and inner boundaries of the region.
Dirichlet condition specifies a known value of electric
potential U0 at the vertex, at the edge, or on the face of the
model (for example on a capacitor plate). This kind of
boundary condition is also useful at an outward boundary of
the region that is formed by the plane of electric
antisymmetry of the problem (opposite charges in
symmetrical geometry). U0 value can be specified as a linear
function of coordinates. The function parameters can vary
Electrostatics 349
𝜕𝐸𝑥 𝐸𝑦
G𝑥𝑥 = , G𝑦𝑦 = ,𝐺
𝜕𝑥 𝜕𝑦 𝑥𝑦
𝜕𝐸𝑥 —for planar
=½ (
𝜕𝑦 case;
𝜕𝐸𝑦
+ )
𝜕𝑥
𝜕𝐸𝑧 𝐸𝑟
G𝑧𝑧 = , G𝑟𝑟 = , 𝐺𝑧𝑟 —for
𝜕𝑥 𝜕𝑟
𝜕𝐸𝑧 𝜕𝐸𝑟 axisymmetric
=½ ( + ) case;
𝜕𝑟 𝜕𝑧
and also its principal components G1 and G2.
• Vector of electric induction D = E, where is electric
permittivity tensor.
Integral quantities:
• Total electric charge in a particular volume
𝑞 = ∮ 𝐃 ∙ 𝐧 𝑑𝑠,
1
𝐅= ∮ (𝐄(𝐧𝐃) + 𝐃(𝐧𝐄)– 𝐧(𝐄𝐃))𝑑𝑠,
2
𝑆
352 Chapter 11 Theoretical Description
1
𝐓= ∮ ([𝐫 × 𝐄](𝐧𝐃)
2
𝑆
+ [𝐫 × 𝐃](𝐧𝐄) – [𝐫 × 𝐧](𝐄𝐃))𝑑𝑠,
where r is a radius vector of the point of integration.
The torque vector is parallel to z-axis in planar case, and
is identically equal to zero in axisymmetric one. The
torque is considered relative to the origin of the
coordinate system. The torque relative to any other
arbitrary point can be obtained by adding extra term of
F r0, where F is the total force and r0 is the radius vector
of the point.
• Energy of electric field
1
𝑊= ∫ (𝐄 · 𝐃)𝑑𝑣
2
𝑉
For planar problems all integral quantities are considered
per unit length in z-direction.
The domain of integration is specified in the plane of the
model as a closed contour consisting of line segments and
circular arcs.
Capacitance Calculation
There are several ways to calculate capacitance using
QuickField. The easiest one of them is based on measuring
an electric potential produced by a known charge. To get
DC Conduction Analysis 353
DC Conduction Analysis
QuickField is able to calculate the distribution of electric
current in systems of conductors. The problems of current
distribution are described by the Poisson's equation for
scalar electric potential U.
354 Chapter 11 Theoretical Description
Field Sources
With the problems of dc conduction, the field sources are
external currents supplied to the boundary of a conductor.
QuickField provides possibility to specify external current
density at the edges or at the individual vertices of the
model. The current density specified at a point of the
xy-plane corresponds to a knife-edge current collector that
is perpendicular to the plane of the model, and is described
by the linear current density. In axisymmetric case the
vertex source represents a circular collector around the axis
of symmetry or a point collector located on the axis. To
incorporate both these cases, a total current value is
associated with the vertex. For the circular knife-edge
collector the total current value is connected with its linear
DC Conduction Analysis 355
Boundary Conditions
The following boundary conditions can be specified at
outward and inner boundaries of the region.
Dirichlet condition specifies a known value of electric
potential U0 at the vertex or at the edge of the model. U0
value at the edge can be specified as a linear function of
coordinates. The function parameters can vary from one
edge to another, but have to be adjusted to avoid
discontinuities at edges' junction points.
𝜕𝑈
𝐸𝑧 = − , 𝐸𝑟 —for
𝜕𝑧 axisymmetric
𝜕𝑈
=− case;
𝜕𝑟
• Vector of current density j = σE, where σ - is
electric conductivity tensor.
Integral quantities:
• Electric current through a given surface
𝐼 = ∫ 𝐣 ∙ 𝐧 𝑑𝑠,
𝑆
where n denotes the vector of the unit normal.
• Joule heat produced in a volume
𝑊 = ∫(𝐄 · 𝐣)𝑑𝑣.
𝑉
For planar problems all integral quantities are considered
per unit length in z direction.
The domain of integration is specified in the plane of the
model as a closed contour consisting of line segments and
circular arcs.
AC Conduction Analysis
AC conduction analysis is the study of electric field, current
and losses arising in conductors and imperfect (lossy)
dielectrics from the application of an alternating (AC)
voltage or external current to electrodes.
358 Chapter 11 Theoretical Description
Field Sources
With ac conduction problems, the field sources are external
currents or voltages supplied to the boundary of a conductor.
QuickField provides possibility to specify external current
density at the edges or at the individual vertices of the
model. The current density specified at a point of the
xy-plane corresponds to a knife-edge current collector,
which is perpendicular to the plane of the model, and is
described by the linear current density. In axisymmetric case
the vertex source represents a circular collector around the
axis of symmetry or a point collector located on the axis. To
AC Conduction Analysis 359
Boundary Conditions
The following boundary conditions can be specified at
outward and inner boundaries of the region.
Dirichlet condition specifies a known value of electric
potential U0 at the vertex or at the edge of the model. U0
value at the edge can be specified as a linear function of
coordinates. The function parameters can vary from one
edge to another, but have to be adjusted to avoid
discontinuities at edges' junction points.
and to the right side of the boundary, and j at right hand side
is a density of the external current. If j value is zero, the
boundary condition is called homogeneous. This kind of
boundary condition is used at an outward boundary of the
region that is formed by the symmetry plane of the problem.
The homogeneous Neumann condition is the natural one; it
is defaulted at all outward boundary parts where no explicit
boundary condition is specified.
If the surface-bound current density is to be specified at the
plane of problem symmetry and this plane is the outward
boundary of the region, the surface current density has to be
halved.
Constant potential boundary condition is used to describe
surface of a conductor having much greater conductivity
than the surrounding medium. This conductor is assumed to
have constant but unknown potential value.
𝜕𝑈
𝐸𝑥 = − , 𝐸𝑦
𝜕𝑥 ⎯ for planar case;
𝜕𝑈
= −
𝜕𝑦
𝜕𝑈
𝐸𝑧 = − , 𝐸𝑟 ⎯ for
𝜕𝑧 axisymmetric
𝜕𝑈
= − case;
𝜕𝑟
• Complex vector of active jactive = σE, reactive
jreactive = iE and apparent japparent = jactive +
jreactive current density;
• Time average and peak active power (losses)
density Qactive = jactive · E, reactive Qreactive
= jreactive · E, and apparent Qapparent = japparent · E
power density;
• Electric permittivity (its largest component in
anisotropic media);
• Electric conductivity g (its largest component
in anisotropic media).
Integral quantities:
• Complex magnitude of electric current (active Iactive,
reactive Ireactive and apparent I) through a given surface
𝐼 = ∫ (𝐣 ∙ 𝐧) 𝑑𝑠
𝑆
where n denotes the vector of the unit normal.
• Time average and peak active Pactive, reactive
Preactive, and apparent P power produced in a
volume
362 Chapter 11 Theoretical Description
𝑃 = ∫ (𝐄 ∙ 𝐣) 𝑑𝑣
𝑉
• Time average and oscillating part of electric force acting
on bodies contained in a particular volume
1
𝐅= ∮ (𝐄(𝐧𝐃) + 𝐃(𝐧𝐄)– 𝐧(𝐄𝐃))𝑑𝑠,
2
𝑆
where integral is evaluated over the boundary of the
volume, and n denotes the vector of the outward unit
normal.
• Time average and peak electric force torque
acting on bodies contained in a particular
volume
1
𝐓= ∮ ([𝐫 × 𝐄](𝐧𝐃)
2
𝑆
+ [𝐫 × 𝐃](𝐧𝐄) – [𝐫 × 𝐧](𝐄𝐃))𝑑𝑠,
where r is a radius vector of the point of integration.
The torque vector is parallel to z-axis in the planar case,
and is identically equal to zero in the axisymmetric one.
The torque is considered relative to the origin of the
coordinate system. The torque relative to any other
arbitrary point can be obtained by adding extra term of
F r0, where F is the total force and r0 is the radius vector
of the point.
Transient Electric Analysis 363
Field Sources
The field sources are external currents or voltages supplied
to the boundary of a conductor. QuickField provides
possibility to specify external current density at the edges or
at the individual vertices of the model. The current density
specified at a point of the xy-plane corresponds to a knife-
edge current collector, which is perpendicular to the plane
of the model, and is described by the linear current density.
In axisymmetric case the vertex source represents a circular
collector around the axis of symmetry or a point collector
located on the axis. To incorporate both these cases, a total
current value is associated with the vertex. For the circular
knife-edge collector the total current value is connected with
364 Chapter 11 Theoretical Description
Boundary Conditions
The following boundary conditions can be specified at
outward and inner boundaries of the region.
Dirichlet condition specifies a known value of electric
potential U0 at the vertex or at the edge of the model. U0
value at the edge can be specified as a linear function of
coordinates. The function parameters can vary from one
edge to another, but have to be adjusted to avoid
discontinuities at edges' junction points.
𝜕𝐸𝑥 𝐸𝑦
G𝑥𝑥 = , G𝑦𝑦 = ,𝐺
𝜕𝑥 𝜕𝑦 𝑥𝑦
𝜕𝐸𝑥 —for planar
=½ (
𝜕𝑦 case;
𝜕𝐸𝑦
+ )
𝜕𝑥
𝜕𝐸𝑧 𝐸𝑟
G𝑧𝑧 = , G𝑟𝑟 = , 𝐺𝑧𝑟 —for
𝜕𝑥 𝜕𝑟
𝜕𝐸𝑧 𝜕𝐸𝑟 axisymmetric
=½ ( + ) case;
𝜕𝑟 𝜕𝑧
and also its principal components G1 and G2.
• Vector of active j = σE and displacement
𝜕
jd = (𝐄) current;
𝜕𝑡
• Active power (losses) density Q = j · E;
• Electric permittivity (E);
• Electric conductivity σ(E).
Integral quantities:
• Electric current (active Iactive and reactive Ireactive) through
a given surface
𝐼 = ∫ 𝐣 ∙ 𝐧 𝑑𝑠,
𝑆
where n denotes the vector of the unit normal.
• Active Pactive and reactive Preactive power
produced in a volume
Transient Electric Analysis 367
𝑃 = ∫ (𝐄 · 𝐣)𝑑𝑣.
𝑉
• Electric force acting on bodies contained in a particular
volume
1
𝐅= ∮ (𝐄(𝐧𝐃) + 𝐃(𝐧𝐄)– 𝐧(𝐄𝐃))𝑑𝑠,
2
𝑆
where integral is evaluated over the boundary of the
volume, and n denotes the vector of the outward unit
normal.
• Electric force torque acting on bodies
contained in a particular volume
1
𝐓= ∮ ([𝐫 × 𝐄](𝐧𝐃)
2
𝑆
+ [𝐫 × 𝐃](𝐧𝐄) – [𝐫 × 𝐧](𝐄𝐃))𝑑𝑠,
where r is a radius vector of the point of integration.
The torque vector is parallel to z-axis in the planar case,
and is identically equal to zero in the axisymmetric one.
The torque is considered relative to the origin of the
coordinate system. The torque relative to any other
arbitrary point can be obtained by adding extra term of
F r0, where F is the total force and r0 is the radius vector
of the point. The domain of integration is specified in the
plane of the model as a closed contour consisting of line
segments and circular arcs.
368 Chapter 11 Theoretical Description
Heat Transfer
With QuickField you can analyze linear and nonlinear
temperature fields in one of two formulations: steady state
or transient: heating or cooling of the system.
Heat-transfer equation for linear problems is:
𝜕 𝜕𝑇 𝜕 𝜕𝑇
(𝜆𝑥 ) + (𝜆𝑦 )
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦
⎯ planar case;
𝜕𝑇
= −𝑞 − 𝑐𝜌
𝜕𝑡
1𝜕 𝜕𝑇 𝜕 𝜕𝑇
(𝜆𝑟 𝑟 ) + (𝜆𝑧 ) ⎯
𝑟 𝜕𝑟 𝜕𝑥 𝜕𝑧 𝜕𝑧 axisymmetric
𝜕𝑇
= −𝑞 − 𝑐𝜌 case;
𝜕𝑡
𝜕 𝜕𝑇 𝜕 𝜕𝑇
(𝜆 ) + (𝜆 )
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦
𝜕 𝜕𝑇 ⎯ 3D case;
+ (𝜆 )
𝜕𝑧 𝜕𝑧
= −𝑞
for nonlinear problems:
𝜕 𝜕𝑇
(𝜆(𝑇) )
𝜕𝑥 𝜕𝑥
𝜕 𝜕𝑇
+ (𝜆(𝑇) )
𝜕𝑦 𝜕𝑦 ⎯ planar case;
= −𝑞(𝑇)
𝜕𝑇
− 𝑐(𝑇)𝜌
𝜕𝑡
Heat Transfer 369
1𝜕 𝜕𝑇
(𝜆(𝑇)𝑟 )
𝑟 𝜕𝑟 𝜕𝑥
𝜕 𝜕𝑇 ⎯
+ (𝜆(𝑇) )
𝜕𝑧 𝜕𝑧 axisymmetric
= −𝑞(𝑇) case;
𝜕𝑇
− 𝑐(𝑇)𝜌
𝜕𝑡
where:
T — temperature;
t — time;
x(y,z,r) — components of heat conductivity tensor,
anisotropy is not accounted for in 3D formulation;
(T) — heat conductivity as a function of
temperature approximated by cubic spline (anisotropy is
not supported in nonlinear case);
q(T) — volume power of heat sources, in linear
case—constant, in nonlinear case—function of temperature
approximated by cubic spline.
c(T) — specific heat, in nonlinear case—function
of temperature approximated by cubic spline;
ρ — density of the substance.
In steady state case the last term in these equations equals
zero. In linear case all the parameters are constants within
each block of the model.
Heat Sources
QuickField provides possibility to specify the heat sources
located in the blocks, at the edges or at the individual
vertices of the model. The heat source specified at a point of
370 Chapter 11 Theoretical Description
Boundary Conditions
The following boundary conditions can be specified at
outward and inner boundaries of the region.
Known temperature boundary condition (known also as
boundary condition of the first kind) specifies a known
value of temperature T0 at the vertex or at the edge of the
model (for example on a liquid-cooled surface). T0 value at
the edge can be specified as a linear function of coordinates.
The function parameters can vary from one edge to another,
but have to be adjusted to avoid discontinuities at edges'
junction points.
This boundary condition sometimes is called the boundary
condition of the first kind.
Heat Transfer 371
𝛷 = ∫ (𝐅 ∙ 𝐧) 𝑑𝑠
𝑆
where n denotes the unit vector of normal to the surface. The
surface is specified by a contour consisting of line segments
and circular arcs in the plane of the model.
374 Chapter 11 Theoretical Description
Stress Analysis
Within QuickField package, the plane stress, the plane strain
and the axisymmetric stress models are available with both
isotropic and orthotropic materials. The plane stress model
is suitable for analyzing structures that are thin in the
out-of-plane direction, e.g., thin plates subject to in-plate
loading. Out-of-plane direct stress and shear stresses are
assumed to be negligible. The plane strain model is
formulated by assuming that out-of-plane strains are
negligible. This model is suitable for structures that are thick
in the out-of-plane direction.
Thermal Strain
Temperature strain is determined by the coefficients of
thermal expansion and difference of temperatures between
strained and strainless states. Components of the thermal
strain for plane stress and isotropic material are defined by
the following equation:
𝛼
{𝜀0 } = {𝛼 } Δ𝑇;
0
plane stress, orthotropic material:
𝛼𝑥
{𝜀0 } = {𝛼𝑦 } Δ𝑇;
0
plane strain, isotropic material:
𝛼
{𝜀0 } = (1 + υ) {𝛼 } Δ𝑇;
0
plane strain, orthotropic material:
𝛼𝑥 + υ𝑧𝑥 𝛼𝑧
{𝜀0 } = {𝛼𝑦 + υ𝑧𝑦 𝛼𝑧 } Δ𝑇;
0
380 Chapter 11 Theoretical Description
External Forces
QuickField provides way to specify concentrated loads,
surface and body forces. The concentrated loads are defined
at vertices as two components of the corresponding vector.
The surface forces at the edges of the model are specified by
the vector components or by the normal pressure. The body
forces are defined by their components within blocks of the
model. Each component of the body force vector can be
specified as a function of the coordinates. This feature can
be used, for example, to model centrifugal forces. The
normal pressure also can be a function of the coordinates
that is useful for hydrostatic pressure.
Stress Analysis 381
Restriction Conditions
Rigid constraint condition along one or both axes can be
specified at any vertex or along any edge of the model.
Prescribed displacement at restrained edge can be specified
as a linear function of the coordinates.
Elastic support condition describes a vertex subject to
springy force which is proportional to difference between
actual and predetermined displacement. The elastic support
condition is characterized by the predetermined
displacement and the support elasticity.
382 Chapter 11 Theoretical Description
1
σ𝑒 = √ [(𝜎1 − 𝜎2 )2 + (𝜎2 − 𝜎3 )2 + (𝜎3 − 𝜎1 )2 ]
2
Stress Analysis 383
1
σ𝑖 = √ [(𝜎1 − 𝜎2 )2 + (𝜎2 − 𝜎3 )2 + (𝜎3 − 𝜎1 )2 ]; 𝜎̅
2
𝜎1 + 𝜎2 + 𝜎3
=
3
• Tsai-Hill failure index for orthotropic materials:
𝜎12 𝜎1 𝜎2 𝜎22 𝜏12
2
𝐶𝑡ℎ = 2− 2 + 2+ 2
𝑋1 𝑋1 𝑋2 𝑆12
where 1, 2 and 12 are computed stresses in the material
directions and,
X1 = X1T if 1>0; X1 = X1C if 1<0
X2 = X2T if 2>0; X2 = X2C if 2<0
384 Chapter 11 Theoretical Description
𝐅 = ∮ (𝜎 ∙ 𝐧)𝑑𝑠
𝑆
where σ is the stress tensor.
The integral is evaluated over the boundary of the
volume, and n denotes the vector of the outward unit
normal.
• Total torque of the forces acting on a particular
volume
1
𝐓= ∮ [𝐫 × (𝜎 ∙ 𝐧)]𝑑𝑠
2 𝑆
where r is a radius vector of the point of integration.
The torque vector is parallel to z-axis in the planar case,
and is identically equal to zero in the axisymmetric one.
The torque is considered relative to the origin of the
coordinate system. The torque relative to any other
arbitrary point can be obtained by adding extra term of
F r0, where F is the total force and r0 is the radius vector
of the point.
Coupled Problems 385
Coupled Problems
QuickField is capable of importing loads (distributed
sources) calculated in some problem into the problem of
another type. Following coupling types are supported:
• Heat transfer caused by Joule heat generated in the
transient or AC magnetic problem, or DC or AC
conduction problem.
• Thermal stresses based on a calculated temperature
distribution.
• Stress analysis of the system loaded by magnetic of
electric forces.
• Electrical conductivity depending on temperature.
• A special case of coupling allows for importing of the
field distribution in some steady state or transient
problem into another transient problem as its initial state.
This applies to transient magnetic and transient heat
transfer analysis.
In addition to imported loads, you can define any other loads
and boundary conditions, similar to non-coupled problem.
You can combine several coupling types in one problem.
E.g., after calculating currents distribution, electrostatic and
magnetic fields as separate problems based on the same
model file, you can calculate temperature distribution from
Joule heat and then find stresses caused by temperature and
magnetic and electric forces at once. However, such
problems are rather rare.
Further we will call the problem, from which the data are
being loaded, the source problem, and the problem, which
imports the data, the target problem.
386 Chapter 11 Theoretical Description
Index
object selection 271
. view 262, 265
zooming 270
.m3d (model) files 41 3D problem
cut plane 294, 297
2 cut plane normal direction
296
2D and 3D problems 253 electric charge 307
electric force 307
3
energy of electric field 307
field calculator 281
3D field presentation modes
2D sketch 251 283
extrusion 251, 254, 257 field view presentations
geometry model 256 284
import 255 finite element mesh 288
level 258, 259 integral values 282, 306
level, add 260 local values 282, 305
level, change 260 objects hiding 286
level, remove 261 result analysis 280
problem 250 results 282
problem definition 254, 255 results window 281
view manipulation 268, 285 results, colored surfaces
workflow 250 289
3D model results, isosurface plot 293
boundary conditions 280 results, scatter plot 292
camera 270 results, vector plot 290
cutting by plane 274 slice plot 300
import 262 solving 280
labels 276 XY plot 302
mesh spacing 277
object hiding 273
391
B
318, 319, 320, 329, 330,
340, 341, 348, 349, 355,
background grid 85 359, 364
beam of charged particles 235 displacement 123, 124, 381
B-H curve 30, 31, 95, 317, 327 elastic support 381
block 54, 55, 316 electric potential 348, 355,
and DC conduction data 359, 364
110 heat flux 371
and electric charge 348 known temperature 370
and electric current 318, Neumann 30, 32, 34, 37,
328 38, 317, 319, 320, 327,
and electric data 115 329, 330, 331, 339, 340,
and electrostatic data 108 341, 348, 349, 355, 359,
and heat source 370 364, 370
and heat transfer data 117 periodic 124
and magnetic data 95 radiation 372
and mesh 82 zero flux 319, 330, 341
and stress analysis data 120
copying 60 C
CAD systems 87
392 Index
S
superconductor 30, 32, 34, 320,
330, 341
selecting
circuit elements 149 T
vertices 81
serial conductor 104 tables
shear modulus 379 in postprocessor 173, 227
solenoid 29, 32 versus time 233
source temperature 38, 39, 45, 118,
thermal 34, 36, 37, 369, 126, 165, 368
371 ambient, for convection
specific heat 369 372
spiral coil 328 ambient, for radiation 372
Stephan-Boltsman constant 372 curve 118
strain 165 difference 171
components 375 terminology 27, 54
plane 40, 315 three-dimensional problem See
thermal 40, 376, 379 3D:problem
stress 165, 171 torque 193, 194
allowable 165, 383 of electric forces 35, 37,
axisymmetric problem 374 352, 362, 367
components 376 of magnetic forces 29, 31,
plane 40, 315 33, 322, 334, 344, 345
principal 40, 165, 168, 382, transformer 31
383 transient analysis
tensor 168 electric 363
thermal 39 transient analysis
von Mises 165, 382 heat transfer 38
data for 117
399
U W
units of length 50 window
document 28
V problem 27
tools 28
vertex 54, 56, 316, 317, 338 wizard 180
and AC conduction data capacitance 35, 185
114 impedance 34, 187
and AC magnetic data 107 inductance 30, 34, 181
and boundary condition
318, 329, 340, 355, 359,
364, 370, 381
X
and DC conduction data X-Y plot 173, 190
112 copying 241
and electric charge 348 printing 241
and electric current 354, versus time 230
358, 363
and electric data 116
and electrostatic data 109
Y
and heat source 369 Young's modulus 379
and heat transfer data 120 Young's modulus 121
and magnetic data 100
and mesh spacing 57, 80,
81 Z
and stress analysis data 124 zooming 84, 172
copying 60