Manual QuickField
Manual QuickField
QuickField
Finite Element Analysis System
Acknowledgments:
TetGen (www.tetgen.org)
CGAL (www.cgal.org)
VTK (www.vtk.org)
Open CASCADE Technology (www.opencascade.org)
All other brand and product names are trademarks or registered trademarks of their
respective owners
iii
Contents
Stress Analysis 24
Index 251
9
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.
11
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:
Professional QuickField in a ZIP-archive - unpack the archive
preserving the directory structure and run Autorun.exe from the root
of the unpacked directory tree.
Student QuickField in a ZIP-archive - unpack the archive in the same
way and run Setup.exe.
QuickField on a CD - insert the CD and, if not started automatically,
run Autorun.exe from the CD root.
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.
12 Chapter 1 Getting Started
Configuration Notes
To solve very large problems on a computer with insufficient memory it is
essential that virtual memory is configured optimally.
To manage virtual memory settings:
1. Bring up Control Panel and double-click System.
2. Switch to Performance tab.
3. See Windows Help for details.
15
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 diamonds, QuickField shows the corresponding rectangle where the
window can be docked if you release the mouse button.
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
18 Chapter 1 Getting Started
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 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 based on a vector potential
formulation. Following options are available for magnetic analysis:
Material properties: air, orthotropic materials with constant
permeability, ferromagnets, current carrying conductors, and permanent
magnets. B-H curves for ferromagnets can easily be defined through an
interactive curve editor, see the "Editing Curves" section in Chapter 5.
QuickField Installation 19
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 are electric
current (and its source and induced component), voltage, generated Joule
heat, magnetic flux density, field intensity, forces, torques, impedance and
inductance. The AC magnetic field simulation can be coupled with
electric circuit. The circuit can contain arbitrarily connected resistors,
capacitors, inductances, and solid conductors located in the magnetic field
region.
A special type of AC magnetic is nonlinear analysis. It allows estimating
with certain precision the behavior of a system with ferromagnets, which
otherwise would require much lengthier transient analysis.
Following options are available for AC magnetic analysis:
Material properties: air, orthotropic materials with constant permeability
or isotropic ferromagnets, current carrying conductors with known current
or voltage.
The electrical conductivity of material can depend on temperature. . The
dependence of conductivity on temperature is given in tabular form using
the Curve Editor. The temperature value can be defined for each block
QuickField Installation 21
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.
Loading sources: voltages, and electric charge density.
Boundary conditions: prescribed potential values (voltages), prescribed
values for normal derivatives (surface charges), and prescribed constraints
for constant potential boundaries with given total charges.
Postprocessing results: voltages, electric fields, gradients of electric
field, flux densities (electric displacements), surface charges, self and
mutual capacitances, forces, torques, and electric energy.
Special features: An integral calculator can evaluate user-defined
integrals along specified contours and surfaces. Floating conductors with
22 Chapter 1 Getting Started
unknown voltages and given charges can be modeled. Electric forces can
be imported into stress analysis (electro-structural coupling). A
Capacitance Wizard is available for calculation of the self- and mutual
capacitance of the conductors.
3D analysis limitations: current QuickField version does not support the
anisotropic materials and coupled problems in 3D Electrostatic analysis.
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:
Material properties: orthotropic materials with constant conductivity.
The electrical conductivity of material can depend on temperature. The
dependence of conductivity on temperature is given in tabular form using
the Curve Editor. The temperature value can be defined for each block
by a number or a formula of coordinates.
Loading sources: voltages, electric current density.
Boundary conditions: prescribed potential values (voltages), prescribed
values for normal derivatives (surface current densities), and prescribed
constraints for constant potential boundaries.
Postprocessing results: voltages, current densities, electric fields, electric
current through a surface, and power losses.
Special features: An integral calculator can evaluate user-defined
integrals along specified contours and surfaces. The electric power losses
can be used as heat sources for thermal analysis (electro-thermal
coupling).
3D analysis limitations: current QuickField version does not support
anisotropic materials, temperature-dependent conductivity and coupled
problems in 3D DC Conduction analysis.
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 quantities of interest are dielectric losses, voltage, electric
field components, forces, and torques.
The following options are available for AC conduction analysis:
QuickField Installation 23
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.
QuickField can perform linear and nonlinear thermal analysis for two-
dimensional and axisymmetric models; and linear steady-state thermal
analysis for three-dimensional models. The program is based on heat
conduction equation with convection and radiation boundary conditions.
Following options are available for thermal analysis:
Material properties: orthotropic materials with constant thermal
conductivity, isotropic temperature dependent conductivities, temperature
dependent specific heat.
Loading sources: constant and temperature dependent volume heat
densities, convective and radiative sources, Joule heat sources imported
from DC or AC conduction or AC or transient magnetic analysis.
Boundary conditions: prescribed temperatures, boundary heat flows,
convection, radiation, and prescribed constraints for constant temperature
boundaries.
Postprocessing results: temperatures, thermal gradients, heat flux
densities, and total heat losses or gains on a given part; with transient
analysis: graphs and tables of time dependency of any quantity in any
given point of a region.
Special features: A postprocessing calculator is available for evaluating
user-defined integrals on given curves and surfaces. Plate models with
varying thickness can be used for thermal analysis. The temperatures can
be used for thermal stress analysis (thermo-structural coupling). Special
type of inter-problem link is provided to import temperature distribution
from another problem as initial state for transient thermal analysis.
3D analysis limitations: current QuickField version does not support
anisotropic materials, non-linear dependencies (dependence of the thermal
conductivity or volume power of the heat source from the temperature),
radiation boundary condition, transient analysis and coupled problems.
Stress Analysis
Stress analysis plays an important role in design of many different
mechanical and electrical components. Generally, the quantities of interest
QuickField Installation 25
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.
To analyze the results, click View Results in the Problem menu or
context menu.
Problem type: Select the type of analysis, which your problem belongs
to.
Editing Problems 29
Precision: Select the precision you need. Note that higher precision leads
to longer solution time.
Formulation: Select the formulation of planar stress analysis problem.
Frequency: Type the value of frequency for the time-harmonic problem.
Note the difference between frequency f and angular frequency :
= 2f.
Files: Edit the file names of your model, data files, and circuit file (if
applicable). You may use long file names. If the name is given without the
full path, it is assumed with respect to the problem description file. You
can also click Browse to select file in any folder on your hard disk or the
network.
Edit: Instantly loads selected file into the new QuickField window.
To establish a link between the problem that imports data and the problem
that originates them, click Links tab in problem description dialog box.
the coupled problem provides the initial temperature distribution then the
imported temperature has a priority.
1
for magnetic problems.
g
The ratio 2/4 is evaluated in all the mesh elements in the model, and the
smallest value is used as an initial time step size.
As the solution progresses, the time steps are adjusted automatically by an
adaptive time stepping scheme.
The next time step is adjusted by
tn +1 = ktn,
where k is a scaling factor varying from 0.25 to 4.0 (with discrete values
of 0.25; 0.5; 1.0; 2.0; 4.0) and dependent on behavior of potential and its
time derivative, as well as all the time- and coordinate-dependent sources
and boundary conditions in the model.
The two factors are taken into account when choosing the value of k:
The norm of time derivative variation on previous time step in all
mesh nodes:
u n u n 1
un 2
u n u n 1
u n Fn F n 1
T
n ,
un K T un
T
n
Chosen units are associated with each particular problem, which gives you
freedom to use different units for different problems. Usually units of
length are chosen before creating the model geometry. It is possible to
change units of length later, but it does not affect physical dimensions of
the model. So, if you create your geometry as a square with 1 m side and
then switch to centimeters, you will get a square measured 100 cm by
100 cm, which is the same as it was before. To actually change size of the
model you should rather use Scaling option of the Move Selection
command of the Model Editor (see page 42 for details).
The choice of length units does not affect units for other physical
parameters, which always use standard SI units. E.g., the current density is
always measured in A/m2 and never in A/mm2. The only physical quantity
that is measured in chosen units of length, is the displacement vector in
stress analysis problems.
places in QuickField where you can make choice between Cartesian and
polar coordinate systems. Using Coordinate System section in problem
description dialog box you can define the default coordinate system
associated with a problem. The same option is also available in the Model
Editor and in the Postprocessor. Definition of orthotropic material
properties, some loads and boundary conditions depends on the choice of
the coordinate system. You can choose Cartesian or polar coordinate
system for each element of data individually and independently from the
default coordinate system associated with the problem. This choice is
available in the dialog boxes of the Data Editor.
existing link, simply clear the source problem filename from the
corresponding property.
36 Chapter 3 Problem Description
37
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 its label. The mesh spacing value defines the approximate
distance between mesh nodes in the neighborhood of the vertex. Define
vertex label to link a vertex with, for example, a line source or load.
Each Edge represents a linear segment or a circular arc connecting two
vertices. Model edges do not intersect each other. Creating new model
edge QuickField splits it as many times as needed at intersection points
with existing model edges and at the points represented by existing model
vertices. QuickField also automatically creates new model vertices
representing intersection points of the new edge and splits the old model
edges at these points. Define edge label to link an edge with, for example,
related boundary conditions.
Each Block represents a continuous subregion of the model plane.
External block boundary is a sequence of edges. Blocks might contain
38 Chapter 4 Model Geometry Definition
Geometry Description
Model development consists of three stages:
Geometry description and manipulation;
Definition of properties, field sources and boundary conditions;
Mesh generation.
Note. Keep in mind that when you click inside a block QuickField select
neither boundary edges nor vertices. Similarly, when you click in the
middle of an edge QuickField does not select either of its ending vertices.
This might be important for correct understanding of such model
operations as Delete, Duplicate, and Move.
If you want to select a block and its boundary edges or an edge and its
ending vertices, drag the mouse to select the required objects with a
rubberband rectangle.
You can also use Select All and Unselect All commands in the Edit or
context menu. Note that you can select objects of different types - blocks,
edges or vertices - at once.
The set of selected model objects is shared between the windows
displaying the model. If several windows display the same model, selected
objects are highlighted in all of them.
Keyboard shortcuts:
Select All CTRL+A
Unselect All CTRL+D
To select all model objects having the same label, click this label in the
Problem Tree View.
Geometry Description 41
The first copy of a model object is always the result of the specified
transformation applied to the object itself. When the transformation allows
to create several copies of every involved object simultaneously, the
second and the following copies of any object are the results of the
transformation applied to the preceding copies
You can also move the selected objects to another location. The only
limitation is that QuickField will not perform moves that change the
model topology. You cannot move vertices or edges into any block or out
of the containing block. To move selected objects, choose Move Selection
in the Edit or context menu. The displayed Move Selection dialog is
similar to the Duplicate Selection dialog described above.
Successful Move preserves all labels and spacing values. Mesh is
preserved in the blocks that are not reshaped.
QuickField always removes the mesh from the reshaped blocks before
checking that the topology remains unchanged. So, if you try a move that
changes the model topology QuickField will block it displaying the
42 Chapter 4 Model Geometry Definition
corresponding message, and in result of the operation you might find that
some of the blocks are no longer meshed.
If you do not like the results of your operation, use Undo to restore the
previous state of the model
Geometric transformations available with move and copy operations are:
Displacement — parallel displacement is applied to selected objects
for specified displacement vector. With copy operation, several
copies can be asked for, it means that copying operation will be
performed several times, each time being applied to the previous
result. Parameters needed are displacement vector components.
Rotation — selected objects are rotated around the specified point for
the specified angle. With copy operation, several copies can be asked
for, it means that copying operation will be performed several times,
each time being applied to the previous result. Parameters needed are
center of rotation coordinates and angle measured in degrees.
Symmetry — selected objects are mirrored; symmetry line is
specified by coordinates of any point on it and the angle between the
horizontal axis and the symmetry line. Positive value of an angle
means counter-clockwise direction. This transformation is available
for copy operation only.
Scaling — selected objects are dilated (constricted) by means of
homothetic transformation. Parameters needed are center of
homothety and scaling factor. This transformation is available for
move operation only.
There is also a more simple method of copying and moving of the
geometric objects – mouse dragging (see Drag and Drop and Clipboard
Editing). Drag-and drop is possible within the same or different model
editor windows.
Deleting Objects
To delete geometric objects:
1. Select the objects you want to delete.
2. In the Edit or context menu, click Delete Selection.
If the selection contains the vertex (vertices) adjacent to exactly two
remaining edges that could be merged together, QuickField, having
deleted the separating vertex, automatically performs the merge.
Otherwise, when one of vertices being deleted is adjacent to one or several
of remaining edges, QuickField adds the adjacent edges to the list of
objects to be deleted and requests the user to confirm the action.
This feature is frequently used for "clipping" of the obsolete parts of
model edges.
Example:
Geometry Description 43
Consider the model shown in Pic.1 below with the semicircles having the
radiuses of 2 and 3 and the common center at (0, 0). Suppose that you
need to create several horizontal edges inside the block with the distance
between consecutive edges equal to 0.5.
The fastest way to create them would be the following:
Set focus to the model window clicking inside it.
Choose Grid Settings from View menu and set Spacing to 0.5.
Press INS to enter the Insert Mode.
Drag mouse from (0, 3) to (4, 3) to create the new edge connecting
these points.
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).
44 Chapter 4 Model Geometry Definition
Pic.3 Pic.4
Keyboard shortcuts:
Delete DEL
Note. Placing the pointer over selection might be difficult when the
selected set of objects does not contain blocks and snap-to-grid option is
on. In such case we suggest you to place the pointer over one of the
vertices you are going to drag.
Keep in mind that if you press a mouse button with the pointer outside of
selection, QuickField, instead of dragging, initiates a rubberband
selection. In such case the shape of the cursor and the color of selected
objects do not change when you press the mouse button down.
The difference between dragging with left and right mouse buttons is
described in Actions Performed on Drop.
Geometry Description 45
Note: the rubberband does not contain any of selected isolated vertices.
This does not mean that these vertices will not be moved or copied. When
isolated vertices constitute the whole selection, the only things that move
during drag are the cursor and the anchor.
When you change the drag mode the rubberband feedback and the shape
of the cursor are changed appropriately.
Drag Modes and Drop Effects
Dragging of model objects can be performed in different modes. The drag
mode used immediately before the drop defines the actions performed by
QuickField.
Drag mode is defined by:
the mouse button you keep pressed while dragging; and
the state of CTRL and ALT keyboard keys before the drop.
There is no way to change the mouse button in the middle of the drag -
you press it at the beginning and release to perform the drop. On the other
hand, you can change the state of CTRL and ALT keyboard keys at any
moment.
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.
To get the specific drop effect choose the drag mode according to the
following rules:
To choose the required action from the displayed context menu drag
with right mouse button and keep control keys released before the
drop.
To move the objects inside the same model preserving connections
between the moved and the stationary parts drag with left mouse
button and keep control keys released before the drop.
Geometry Description 47
The Move operation never moves objects from one model to another.
Drag-moving can move objects to another model breaking the
connections with the stationary part of the source model.
The Move operation does not allow changing the model topology - it
does not allow moving vertices to other blocks, or creation of
intersections. Drag-moving inside the same model does not have any
limitations.
The Move operation always preserves the labels of the related
objects. Drag-moving might cause labels to change.
See also Dragging to another Model.
Undo after Drag and Drop
QuickField Model Editor performs Undo/Redo operations on per-model
basis. It maintains separate stack of model states for every model and,
when you request Undo or Redo, restores the state of the model before or
after the corresponding operation.
If several model windows are opened simultaneously, Model Editor
performs Undo/Redo for the model displayed in the active window. To
make another window active, click anywhere inside it. Click scrollbar if
you do not want to change the current selection set of the model.
Unlike other Model Editor operations, Drag and Drop might affect two
different models at once. When you drag a group of items moving them
from one model to another, QuickField changes both the source and the
target models. In case you want to Undo the effects of the whole operation
you need to do it for each of these models separately. If you decide to
Undo the effects only on one of the models, you should be careful and
prior to performing Undo make sure that proper model window is active at
the moment.
Dragging to another View
In some cases moving model objects could be quite inconvenient. For
example, this would be the case when you need to move or copy relatively
small objects across relatively large spaces. Inconveniency would be
caused by the fact that it is impossible simultaneously select the small
source objects and fit the target place inside the window.
This inconveniency could be easily eliminated with Drag and Drop
between different views of the same model. Try the following:
Open the second window for the same model choosing New Window
from Window menu.
Arrange the windows so that both are visible.
Zoom in the first window on the source objects. Select those objects
you want to move or copy.
Zoom in the second window on the target place.
Drag the selected objects from the first to the second window.
Geometry Description 49
The pasted objects remain selected after the operation. All other objects in
the target model become unselected.
Keyboard shortcuts:
Copy CTRL+C
Cut CTRL+X
Paste CTRL+V
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 geometrical model database associated with
these operations. The topmost increment on stack corresponds to the latest
editing operation performed. QuickField also keeps track of its current
stack position that steps in top-to-bottom direction with every Undo action
and steps in the opposite direction with every Redo action.
The depth of this internal stack defines maximum number of operations
you can undo. When the total number of editing operations on the model
exceeds the stack depth, the database increments corresponding to the
eldest operations are destroyed to free stack positions for new increments.
When you start QuickField the stack depth has its default value of 25 that
allows to undo last 25 editing operations done on the model.
Maintaining geometrical model database increments impacts QuickField
memory requirements. To make you able to influence this impact
QuickField provides the possibility to reduce the depth of internal undo
stack. You can do it at any time setting the depth value to any integer
between 0 and 100 (both limits included). Ultimately, setting the stack
Geometry Description 51
Undoable Operations
You can undo the following types of geometric operations (we use menu
item labels whenever appropriate):
Add Edge
Add Vertices
Build Mesh
Cut
Delete Edges
Delete Vertices
Delete Selection
Drag
Drop
Duplicate
Import DXF
Move Selection
Paste
Properties
Refine Mesh
Remove Mesh
Select
Select All
Unselect All
Meshing Technology
Having described a reasonable part of the model geometry you can start
building the finite element mesh. It is important to remember that the
mesh you build even for a highly complex geometry can be rather non-
uniform. Namely, you can set the size of mesh elements in some of the
model blocks to be much less than in the other blocks. In such cases the
meshing technology we call geometric decomposition will automatically
produce smooth transition from large to small elements of the mesh. As a
rule, the places where the mesh has to be the finest are those with the
highest field gradient and those where you need higher precision.
In case of simple model geometry we suggest to leave calculation of
element sizes to QuickField. The same applies to the case of preliminary
design analysis when draft precision should perfectly suit your needs.
Click the Build Mesh toolbar button and QuickField will automatically
generate suitable mesh.
It might happen, however, that the quality of the automatically generated
mesh would not satisfy you. For such cases, QuickField provides you with
the way to set the mesh density manually. You control the mesh density
defining mesh spacing values for particular model vertices. The mesh
spacing value defined for a vertex specifies the approximate distance
between adjacent mesh nodes in the neighborhood of the vertex.
You never need to define the spacing for all model vertices. To obtain a
uniform mesh it would be enough to set the spacing for a single vertex. If
you need non-uniform mesh start with defining the spacing values only
for the vertices where you need the finest and the roughest mesh. The
spacing values will be automatically interpolated to other model vertices
smoothing the mesh density distribution across the meshed blocks. Use
54 Chapter 4 Model Geometry Definition
The figure shows the statistics related to the geometry model as well as
the properties of the selected geometric objects (blocks, edges, and nodes).
Some properties like the object label or the discretization step may be
changed by typing it in or selecting the new value from the drop-down
list.
The label value is set separately for every type of geometric object.
To set the step of discretization, switch the mode to Manual and enter the
new Spacing value or select one from the list.
Zooming the Model - gives you the ability to see more or less of
your model to deal with small or large objects.
Setting Model Discretization Visibility - makes the picture more
suitable to specific tasks of model creation process.
Changing Background Grid Settings - simplifies creation of model
vertices and edges.
Consider also opening several windows for the same model and tuning
them differently. To do so, choose New Window from the Window
menu.
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 -
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.
The Show Grid option allows to switch grid visibility on and off.
The Snap to Grid option allows to switch grid attraction on and off.
When attraction is on new vertices can appear either at intersection points
of model edges or at grid nodes. This makes model description faster and
results in more consistent models.
The Grid Spacing field defines the sizes of the grid cell. To define
different horizontal and vertical sizes, first check the Anisotropic box.
The Anisotropic option allows to define the grid cell with different
horizontal and vertical sizes.
The Scale with Zoom option allows to switch between different policies
related to changing grid cell sizes with zoom.
When the option is on (default) the grid cell sizes retain approximately the
same values in screen coordinates.
When the option is off the grid cell sizes retain the same values in length
units defined for the model.
The Grid Origin field defines the coordinates of one of grid nodes. This
allows to create vertices at even distances from the point with such
coordinates.
58 Chapter 4 Model Geometry Definition
Ignore extents? is a flag for ignoring the size of the model fragment
recorded in the DXF file. In some cases this information is incorrect and
may prevent the successful model import.
We provide the special add-in to QuickField that allows you to import
sketches from running SolidWorks® application to QuickField models.
Choose Import from SolidWorks... command from the Edit menu to
invoke this add-in.
For more details refer to online help library.
It is also possible to export the model mesh to an ASCII text file. This
feature allows to interact with other finite element analysis programs. In
particular, you can use QuickField as a mesh generator for your own FEM
solver.
The detailed description of the resulting file format can be found in the
QuickField Help. In brief – the file contains the data defining the model
geometry, the finite element mesh and the labels assigned to model
objects. We provide the utility allowing to import files in this format to
MATLAB® in the form compatible with Partial Differential Equations
Toolbox (PDE Toolbox).
C H A P T E R 5
DC conduction .dcf
AC conduction .dec
Transient electric analysis .dtv
Steady state and Transient heat transfer .dht
Stress analysis .dsa
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 labels besides tensor components reflect the coordinate
system (Cartesian or Polar) selected for the property.
For transient problems, QuickField allows specifying non-zero electric
conductivity in any block. In such case, it calculates the distribution of
eddy currents in this block. The electrical conductivity of the material may
depend on temperature. Dependency is given in tabular form and
automatically be approximated by a spline. To specify the electrical
conductivity, which depends on temperature, check the
Function of Temperature check box to get into the curve editor.
Enter a number or a formula into the Temperature field. The formula
describes the temperature dependency on time and coordinates. Please
note that all coordinates in a formula are given in meters no regards to
length units you have choose for the problem
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.
For non-transient problems, as well as for zero electric conductivity
blocks in transient problems, you can define the field source either by
current density value or by the total number of ampere-turns (total
current). For total current, QuickField will, depending on your choice,
either consider the blocks labeled with the label as single conductor or as
several conductors connected in series. Serially connected conductors
always carry the same current with calculated current densities inversely
proportional to their squares.
Providing total number of ampere-turns for axisymmetric problems, you
can additionally specify that current density in your coil varies inversely
66 Chapter 5 Problem Parameters Description
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 that blocks with this 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 labels besides tensor components reflect the coordinate
system (Cartesian or Polar) selected for the property.
When dealing with nonlinear materials, instead of magnetic permeability
you need to define the magnetization curve. To do so, check the
Nonlinear box and QuickField will display the Edit B-H Curve dialog
that allows to define the curve. To reopen the Edit Curve dialog later,
click the Edit B-H Curve button.
Defining a magnetization curve with the Curve Editor you specify the DC-
based B(H) dependency and QuickField automatically recalculates the
curve for the problem-defined frequency. The graph shows the original
DC-based curve in green color while the curve recalculated for the
problem-defined frequency is dashed red.
Enter the value of electrical conductivity to the Conductivity field or
leave zero for non-conductive material. The electrical conductivity of the
material may depend on temperature. Dependency is given in tabular form
and automatically be approximated by a spline. To specify the electrical
conductivity, which depends on temperature, check the
Function of Temperature check box, to get into the curve editor.
Enter a number or a formula into the Temperature field. The formula
describes the temperature dependency on coordinates. Please note that all
coordinates in a formula are given in meters no regards to length units you
have choose for the problem.
You can import the temperature field from a coupled heat transfer
problem. To do that create and solve a heat transfer problem based on the
same geometry model, and link it to the problem. When the imported
temperature field and temperature given by a constant or a formula are
both available, the imported temperature will be used.
Editing Label Data 69
If the field analysis is performed with the electric circuit connected, then
applied voltage or the full current for the solid conductor (with non-zero
conductivity) cannot be defined in the label properties window. Instead of
it, include all the conductive blocks into electric circuit with voltage or
current sources attached to them. Only the parallel or serial connection of
the separate conductors with the same label assigned should be defined in
the label properties window.
The ways you define field sources are different for conductors and
non-conductive blocks. For solid conductors, you specify either applied
voltage or total current. For non-conductive blocks and stranded
conductors you always specify zero conductivity value and the field
source can only be specified by total current or current density values.
You can specify coordinate-dependent current density phase and
magnitude. To do it, enter the required formula in place of numerical
value. Formula syntax and other details are discussed in detail in "Using
Formulas" section later in this chapter.
For total current and applied voltage, QuickField will, depending on your
choice, either consider the blocks labeled with the label as single
conductor or as several conductors connected in series. Serially connected
conductors always carry the same current, whereas current densities will
be calculated when QuickField finishes solution of the problem.
Note. Unless the conductors are connected in series, the value of total
current associated with a block label specifies the gross current in all
blocks labeled with that label.
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.
First term of the formula above corresponds to the hysteresis losses,
second- to the Eddy current losses, and third approximates the additional
magnetic losses, not covered by first two loss types.
Specific loss coefficients for the given material are calculated outside
QuickField by data fitting using the known or measured tables of the
volume losses per flux density and frequency pст = f(B, f) . More details
about the loss coefficient calculations may be found at
www.quickfield.com/glossary/core_loss_coefficients.htm
Editing Label Data 71
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.
2. Loads
The values of allowable stresses do not affect problem solution. They are
used only in postprocessing to calculate the Mohr-Coulomb,
Drucker-Prager, and Hill criteria. You don't need to define allowable
stresses, if these criteria are of no interest to you.
Unlike other finite element packages, QuickField does not require that the
mesh be equal on both boundaries, and is capable of merging potential
values on edges with mismatching edge mesh.
To apply this type of condition, simply check the correspondent box
(Even or Odd Periodical) for the label assigned to edges on the two
boundaries – QuickField will analyze the geometry and detect the
periodicity automatically.
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.
To correct a table value, simply type the new value in the table cell.
To copy table rows to the Clipboard, either select these rows in the table
or select the corresponding points on the graph and press ( CTRL+C). The
copied rows could be later pasted into the table.
To remove the point, select it in the table or on the graph and click Delete
button or press the DEL key.
You may control the scaling of the graph with use of zoom buttons. To
make the whole graph visible, click Zoom to fit button.
Dragging a dialog window boundary allows to resize the window.
QuickField remembers the resulting window size and position for future
use.
The graph image could be copied to the Clipboard and/or saved to a file.
To do that, invoke the corresponding commands in the context (right-
click) menu of the graph.
To finish editing, click Close or press ESC. Note that subsequent canceling
of label data editing with ESC or Cancel will discard all changes including
those made during curve editing.
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.
Physical properties allowing definitions with formulas
The following properties can be defines with formulas (f(x,y) denotes
function of coordinates, f(t) denotes function of time).
Magnetostatics
Editing Label Data 87
Transient magnetics
Block labels f(x,y) f(t)
Current density in stranded conductors + +
Total current +
Voltage applied to solid conductors +
Edge labels f(x,y) f(t)
Magnetic potential + +
Surface current density + +
Vertex labels f(x,y) f(t)
Concentrated current + +
Magnetic potential + +
Electrostatics
Block labels f(x,y) f(t)
Electric charge density +
Edge labels f(x,y) f(t)
88 Chapter 5 Problem Parameters Description
Voltage +
Normal component of electric induction +
Vertex Labels f(x,y) f(t)
Concentrated (linear) charge +
Voltage +
DC Conduction
Block Labels f(x,y) f(t)
Not applicable
Edge Labels f(x,y) f(t)
Voltage +
Normal current density +
Vertex Labels f(x,y) f(t)
Concentrated current +
Voltage +
Transient Electric
Block labels f(x,y) f(t)
Not applicable
Edge labels f(x,y) f(t)
Voltage + +
Normal current density + +
Vertex Labels f(x,y) f(t)
Concentrated current + +
Voltage + +
Editing Label Data 89
Heat Transfer
Block labels f(x,y) f(t)*
Heat source volume power + +
Edge labels f(x,y) f(t)*
Temperature + +
Heat flux + +
Film coefficient and temperature of contacting fluid
+ +
medium
Emissivity coefficient and ambient radiation
+ +
temperature
Vertex labels f(x,y) f(t)*
Concentrated (linear) heat source + +
Temperature + +
Syntax
QuickField formula is an expression composed of the following elements:
Numerical constants
Integral (Example: 123)
Fixed-point (Examples: 123.45 123. 0.123 .123)
Floating-point (Examples: 1e12 5.39e+8 0.1E-12 .2E+2)
Arithmetic operators
+ addition (Ex: 2+2)
- subtraction (Ex: 3-5)
* multiplication (Ex: 1.23*0.12)
/ division (Ex: 1E5/0.01)
90 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 (',').
Functions
Name Formula Comments
abs 𝑣, 𝑖𝑓 𝑣 ≥ 0 The function value is equal
𝑎𝑏𝑠(𝑣) = |𝑣| = {
−𝑣, 𝑖𝑓 𝑣 < 0 to absolute value of the
argument.
impulse 𝑖𝑚𝑝𝑢𝑙𝑠𝑒(𝑣, 𝑙, 𝑟) =
0, 𝑖𝑓 𝑣 < 𝑙
= {1, 𝑖𝑓 𝑙 ≤ 𝑣 ≤ 𝑟
0, 𝑖𝑓 𝑣 > 𝑟
saw saw(,p) =
𝑣/𝑝, 𝑖𝑓 0 ≤ 𝑣 < 𝑝
{𝑠𝑎𝑤(𝑣 + 𝑝), 𝑖𝑓 𝑣 < 0
𝑠𝑎𝑤(𝑣 − 𝑝), 𝑖𝑓 𝑣 ≥ 𝑝
saw(,p,p0) =
𝑣
, 𝑖𝑓 0 ≤ 𝑣 < 𝑝
𝑝
0, 𝑖𝑓 𝑝 ≤ 𝑣 < 𝑝 + 𝑝0
𝑠𝑎𝑤(𝑣 + (𝑝 + 𝑝0)), 𝑖𝑓 𝑣 < 0
{𝑠𝑎𝑤(𝑣 − (𝑝 + 𝑝0)), 𝑖𝑓 𝑣 ≥ 𝑝 + 𝑝0
Editing Label Data 95
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.
Mathematical notation Formula syntax
t (1-t) (2-t) t*(1-t)*(2-t)
2t2 - t – 3 2*t^2 - t - 3
e-t²/2 exp(-t^2 / 2)
2t 2^t
arcsin 2 asin(sqrt(2))
𝑡
tan tan(t / 2.4e-8)
2.4 ∙ 10−8
|2t| abs(2*pi*t)
𝑡, 𝑖𝑓 𝑡 < 0.5
{ t*step(0.5-t) + (1-t)*step(t-0.5)
1 − 𝑡, 𝑖𝑓 𝑡 ≥ 0.5
𝑡, 𝑖𝑓 0 ≤ 𝑡 < 0.5
t*impulse(t,0,0.5) + (1-
{1 − 𝑡, 𝑖𝑓 0.5 ≤ 𝑡 < 1 t)*impulse(t,0.5,1)
0, 𝑜𝑡ℎ𝑒𝑟𝑤𝑖𝑠𝑒
96 Chapter 5 Problem Parameters Description
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,
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.
Editing Circuit
Moving, Copying and Resizing Circuit Elements
To move circuit element to another place:
1. Place cursor over element you want to move. Cursor should have the
shape of four-pointed arrow.
2. Click left mouse button and drag the selected element keeping the
button pressed.
To resize a wire:
Editing Circuit 101
1. Place cursor over the wire end point. Cursor should have the shape of
two-pointed arrow.
2. Click mouse button and drag the end point to a new location.
You can not only resize the wire, but also move it sideway during this
operation.
To move several elements at once:
1. Select circuit elements you want to move.
2. Place cursor over one of selected elements.
3. Click mouse button and drag the selected elements.
Dragging with Attached Elements
When you drag the circuit elements, another element attached to them
could be dragged or resized to preserve connections between elements and
wires. For example, when you drag some electric component the wires
attached to it could be also moved or resized.
Dragging without Attached Elements
Sometimes it is more convenient to drag circuit elements so that attached
wires and elements are not dragged. For example, you may want to
change the circuit topology by placing the electric component to another
position. For this, 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
This chapter describes how to solve the prepared problem, and methods
QuickField uses to solve.
Several conditions have to be met to solve a problem. The problem type,
plane, required precision and other parameters have to be specified in the
problem description file. The model geometry file must contain complete
model with mesh and labels. Each label referred by the model file is to be
defined in the problem's private or library data file.
To obtain the problem solution, click Solve Problem in the Problem
menu or context (right mouse button) menu of the Problem editor. You
may skip this action and directly proceed to the analysis results by
clicking Analyze Results in the Problem or context menu. If the problem
has not been solved yet, or its results are out of date, the solver will be
invoked automatically.
Each solver runs in its separate thread, so you can solve several problems
at once or edit or analyze other problems while the problem is being
solved. There is of course no since in editing any document related to the
problem being solved.
Special bar indicator lets you see the progress of the solution process. To
interrupt it, click Cancel on the indicator's panel. When solving a transient
problem, you have an option to keep the results for already stored time
steps.
Linear problems are solved by using a powerful preconditioned conjugate
gradient method. The preconditioning based on the geometric
decomposition technique guaranties a very high speed and close to linear
dependence between number of nodes and the resulting solution time.
Nonlinear problems are solved using the Newton-Raphson method. The
Jacobian matrix arising at each step of the Newton-Raphson method is
inverted the same way as it is done for linear problems.
We use the Euler’s method (constant time step size) for solving transient
problems, with initial value set to zero or taken from another field
calculation. This method is extremely fast and stable, however we
104 Chapter 7 Solving the Problem
recommend having at least 15-20 time steps for the whole transitional
process to achieve accurate and smooth results.
C H A P T E R 8
Analyzing Solution
1 2
𝑄= 𝑗
𝜎
Building the Field Picture on the Screen 107
1
𝐇 = 𝐁,
𝜇
where is the magnetic permeability tensor.
Complex vectors may be shown in form of momentary, RMS or peak
magnitude.
Time average and peak Joule heat density
108 Chapter 8 Analyzing Solution
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 were coupled with the connected
electric circuit, then the following parameters would 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;
1
σ𝑒 = √ [(𝜎1 − 𝜎2 )2 + (𝜎2 − 𝜎3 )2 + (𝜎3 − 𝜎1 )2 ]
2
where = [+]/[-], and [+] and [-] denote tensile and compressive
allowable stress.
Drucker-Prager criterion:
2
√𝜒 − 𝜒 1 1 − √𝜒
𝜎𝑒 = (1 + √𝜒)𝜎𝑖 − 𝜎̅ + ( ∙ 𝜎̅)
1 + √𝜒 [𝜎− ] 1 + √𝜒
where
1 𝜎1 + 𝜎2 + 𝜎3
σ𝑖 = √ [(𝜎1 − 𝜎2 )2 + (𝜎2 − 𝜎3 )2 + (𝜎3 − 𝜎1 )2 ]; 𝜎̅ =
2 3
With these methods, you can change the grid cell size and the scaling
factors in order to manipulate the appearance.
It is possible to combine several visualization methods in the same picture
to obtain the most expressive result.
QuickField can display several different field pictures for the same
problem. To open a new window, click New Window in the Window
menu.
The Field View dialog box for the stress analysis problem additionally
allows to select tensor quantity visualization.
Sizes of the vector symbols for all vector quantities except the
displacement vector are determined by the corresponding physical value
multiplied by the scaling factor and by the cell size. Similar method is
used for stress tensor components. Unlike other vector quantities, the size
of the displacement vector on the screen does not depend on the cell size.
It is determined by the dimensionless scaling factor, the unit value of
which means that the displacement is shown in its natural scale.
114 Chapter 8 Analyzing Solution
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.
When you change the time value QuickField adjusts both the field
pictures, XY-Plots and tables. The scaling factors will remain unchanged.
For time-harmonic problems the displayed field picture corresponds to
the selected phase. QuickField uses the same combobox to display the
selected phase value in degrees. To change it, specify another value in
Calculator Window 115
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.
Note. Even though the problem results can be stored with varying time
steps, the animation shows the frames in equal time intervals.
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;
116 Chapter 8 Analyzing Solution
Effective value of the current in the circuit component and voltage drop
along it is shown in the tool tip window displayed with small delay after
the mouse cursor points to the circuit component. More details about the
circuit component are available in the component list pane.
Every list item is presented by its definition (shown in semi-bold font) and
two or three groups of data. Electric current and voltage drop are
displayed for every item; the nominal value is also shown for passive
components and sources.
Components selection in the list and schema is synchronized
automatically.
Every circuit component corresponds to several lines in the list. By default
all lines but the first one are hidden. To expand them you should click on
the small + (plus) sign to the left from the component.
In the time-harmonic problems for the current and voltage drop in every
circuit branch component effective value, instantaneous value for the
selected phase, and complex value components in the algebraic and
trigonometric forms are shown.
It is possible to select any number of lines in the circuit components list
and put them into the clipboard (Copy command in the Edit menu).
118 Chapter 8 Analyzing Solution
Selected lines may also be dragged into other application, for example text
editor or spreadsheet.
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.
Parameter Calculation Wizards 119
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 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 equation above becomes extremely simple and you can
get inductance value as:
L = /i,
where is the flux linkage with the coil excited by current i, or
L =2W/i2,
where W is stored magnetic energy and i is the only current.
The first approach gives the self-inductance, if you get the flux linkage
and the current in the same coil and mutual inductance if the coils are
different. The second approach gives only the self-inductance.
Initial page of the inductance wizards invites you to choose between two
approaches described above. After choosing one of them click the Next
button.
The second page of inductance wizard allows you to define, which blocks
represent the cross section of your coil. In general, two blocks represent
each coil in the model plane: forward and return wires. If there is only one
side of the coil in your model, the second one is assumed as being
symmetrical to the first one or as being infinitely distant of the model and
not affecting the field distribution.
Parameter Calculation Wizards 121
To define each side of your coil, simply point the corresponding item in
the Block Labels list and drag it to one of the side list. You can also use
the Add buttons. No matter, which side of your coil you call Right Side
and which Left Side. If only one side of the coil is represented in the
model, drag item Symmetry to the opposite list if return wire of the coil is
symmetrical to the direct one, or leave the list empty if return wire does
not affect the electromagnetic state of your model.
You can select and drag more than one item at once if the cross section of
your coil is split to several blocks.
Enter the Number of Turns for your coil if it is more than one.
As result of any action on the lists or number of turns the Flux Linkage
value will change automatically being calculated as
∫𝐿 𝐴 𝑑𝑠 ∫𝑅 𝐴 𝑑𝑠
𝜙 =𝑁∙( − ) for planar case;
∫𝐿 𝑑𝑠 ∫𝑅 𝑑𝑠
∫𝑅 𝑟𝐴 𝑑𝑠 ∫𝐿 𝑟𝐴 𝑑𝑠
𝜙 = 2π𝑁 ∙ ( − ) for axisymmetric case;
∫𝑅 𝑑𝑠 ∫𝐿 𝑑𝑠
where A is the vector magnetic potential; R and L denote the right and the
left side of the coil accordingly, r is the radius of the point.
For planar problems flux linkage and the inductance are calculated per one
meter of axial depth no matter what length unit you have chosen.
When you finish with flux linkage calculation, click on the Next button.
In the Current page you can select the current exciting the field and
provide a number of turns in your coil.
122 Chapter 8 Analyzing Solution
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 produce or vice versa. Measuring the
charge is a bit more complex than the voltage. It requires you to build the
closed contour surrounding your conductor (but not coinciding with its
surface) before you start the capacitance wizard. The easiest way to
formulate the problem for capacitance calculating is to put constant
Parameter Calculation Wizards 123
In the right side of the page electrodes are listed which charge you have
specified. If you have put voltage boundary condition rather then charge
124 Chapter 8 Analyzing Solution
on your electrodes, the only way to calculate the charge is to build the
contour surrounding it but not coincident with it boundary. If so, you have
to do that before you start the capacitance wizard.
When selecting one or more items in the list, you get the resulting charge
in the Charge box.
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,
where Z is absolute value of the impedance and f is the frequency.
If you select more than one conductor at once, the impedance wizard
considers it as being connected in parallel if the voltage applied to each of
them is equal and as being connected in series otherwise.
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 125
For volume integrals, the domain of integration lies to the left of the
contour.
For surface integrals, the positive normal vector points to the right
relative to the contour direction.
The starting point of the contour corresponds to zero point at the x-axis
of the X-Y plot.
If the plotted or the integrated function has different values to the left
and to the right of the contour, the right-hand value is used.
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:
Select the set of shown quantities. Click X-Y Plot Curves in the View
or context menu.
Zoom the plot in or out.
View the correspondence between quantities and curves (legend).
Copy the picture to clipboard.
Open new X-Y plot window for the same contour.
Calculating Integrals 127
Few quantities having the same unit of measurement can be shown at the
same X-Y plot. According to this, all quantities are combined into groups.
Full list of quantities includes all those available for the color map
representation (see “Interpreted Quantities”), and also normal and
tangential components of vector and scalar quantities.
When you select the appropriate group of quantities, the Curves to Show
list contains the quantities selected for display, and the Available
Quantities list contains available but not selected quantities. You can use
buttons located between the lists, or simply double-click in the lists, to
move some quantity from one list to another.
In the dialog box, you can also modify the range of y coordinate. By
default, it fits all the currently selected curves. You can get the suggested
value of lower or upper limit by selecting the corresponding text box
(Minimum or Maximum) and choosing Suggest button.
In time-harmonic analysis, you can also switch between momentary (at
given phase), time average and peak values of time dependent quantities.
You can turn on or off the switches for displaying coordinate grid and
markers on the curves. The last mode allows you to distinguish between
the coinciding curves.
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
128 Chapter 8 Analyzing Solution
Name,
Formula and Description
ActiveFied constant
Mechanical force
1
qfInt_MaxwellForce 𝐅= ∮ (𝐇(𝐧𝐁) + 𝐁(𝐧𝐇)– 𝐧(𝐇𝐁))𝑑𝑠,
2
𝑆
2𝜋 ∮(𝑟𝐴)𝑑𝑠 for
Ψ =
𝑆𝐶 axisymmetric
case.
Magnetomotive force
𝐹 = ∫ (𝐇 ∙ 𝐭)𝑑𝑙
qfInt_KGrad_t_dl 𝐿
Magnetic flux
𝛷 = ∫ (𝐁 ∙ 𝐧)𝑑𝑠
qfInt_Grad_n_ds 𝑆
𝑊
𝐵
nonlinear
= ∫ (∫ 𝐻(𝐵′)𝑑𝐵′) 𝑑𝑣 case.
0
Surface energy
1
qfInt_GradKGrad_n_ds 𝑊𝑠 = ∫ (𝐁 · 𝐇)𝑑𝑠
2
𝑆
Average volume 1
potential 𝐴𝑉 = ∫ 𝐴𝑑𝑣
𝑉 𝑉
qfInt_Potential_dv
Surface integral of
strength 𝑥 = ∫ (𝐇 ∙ 𝐧)𝑑𝑠
𝑆
qfInt_KGrad_n_ds
Total current
𝐼 = ∫ 𝑗𝑑𝑠
qfInt_Jtotal 𝑆𝑐
External current
𝐼𝑒𝑥𝑡 = ∫ 𝑗𝑒𝑥𝑡 𝑑𝑠
qfInt_Jextern 𝑆𝑐
Eddy current
𝐼𝑒𝑑𝑑𝑦 = ∫ 𝑗𝑒𝑑𝑑𝑦 𝑑𝑠
qfInt_Jeddies 𝑆𝑐
Joule heat 1
𝑃 = ∫ 𝑗 2 𝑑𝑣
qfInt_Power 𝜎
Joule heat power in a volume.
σ – electric conductivity,
j – total current density.
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
134 Chapter 8 Analyzing Solution
Joule heat 1 2
𝑃 = ∫ ∙𝑗 𝑑𝑣
qfInt_Power 𝑉 𝜎 𝑡𝑜𝑡𝑎𝑙
Quadratic value.
Joule heat power in a particular volume.
σ – electric conductivity of the media.
Core loss 𝑘ℎ ∙ 𝑓 ∙ 𝐵 2 + 𝑘𝑐 ∙ 𝑓 2 ∙ 𝐵 2 +
𝑃= ∫ [ 3 ] 𝑑𝑣
qfInt_Steinmetz 𝑉 +𝑘𝑒 ∙ (𝑓 ∙ 𝐵)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.
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
Calculating Integrals 135
point.
Lorentz force
qfInt_LorentzForce 𝐅 = ∫ [𝐣 × 𝐁]𝑑𝑣.
𝑉
Quadratic vector.
The Lorentz force acting on conductors
contained in a particular volume.
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.
Quadratic value.
This formula is used for both linear and
nonlinear cases.
Complex value.
The integral has to be evaluated over a cross
section of the coil, and SC is the area of the
cross section.
Magnetomotive force
𝐹 = ∫ (𝐇 ∙ 𝐭)𝑑𝑙
qfInt_KGrad_t_dl 𝐿
Complex value.
Magnetomotive force.
Magnetic flux
𝛷 = ∫ (𝐁 ∙ 𝐧)𝑑𝑠
qfInt_Grad_n_ds 𝑆
Complex value.
Magnetic flux through a particular surface.
136 Chapter 8 Analyzing Solution
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.
Average volume 1
strength 𝐇𝑎 = ∫ 𝐇𝑑𝑣
𝑉 𝑉
qfInt_KGrad_dv Complex vector.
Surface integral of
strength 𝑥 = ∫ (𝐇 ∙ 𝐧)𝑑𝑠
𝑆
qfInt_KGrad_n_ds Complex value.
Calculating Integrals 137
Name, ActiveFied
Formula and Description
constant
Electrical charge
𝑞 = ∮ 𝐃 ∙ 𝐧 𝑑𝑠
qfInt_KGrad_n_ds 𝑆
Mechanical force
1
qfInt_MaxwellForce 𝐅= ∮ (𝐄(𝐧𝐃) + 𝐃(𝐧𝐄)– 𝐧(𝐄𝐃))𝑑𝑠
2
𝑆
Mechanical torque 𝐓=
qfInt_MaxwellTorque 1 [𝐫 × 𝐄](𝐧𝐃) +
= ∮( ) 𝑑𝑠
2 +[𝐫 × 𝐃](𝐧𝐄)– [𝐫 × 𝐧](𝐄𝐃)
𝑆
Stored energy
1
qfInt_ElectrostaticEnergy 𝑊= ∫ (𝐄 · 𝐃)𝑑𝑣
2
𝑉
Surface energy
qfInt_GradKGrad_n_ds 𝑊𝑆 = ∫ (𝐄 · 𝐃)𝑑𝑠
𝑆
Potential difference
qfInt_Grad_t_dl ∆𝑈 = ∫(𝐄 · 𝐭)𝑑𝑙
𝐿
Average volume 1
potential 𝑈𝑉 = ∫ 𝑈𝑑𝑣
𝑉 𝑉
qfInt_Potential_dv
Average volume 1
displacement 𝐃𝑎 = ∫ 𝐃𝑑𝑣
𝑉 𝑉
qfInt_KGrad_dv Average electric displacement vector in a
particular volume.
Mean square 1
displacement 𝐷𝑎2 = ∫ 𝐷2 𝑑𝑣
𝑉 𝑉
qfInt_KGrad_dv
Line integral of
displacement 𝑥 = ∫ (𝐃 ∙ 𝐭)𝑑𝑙
𝐿
qfInt_KGrad_t_dl
Surface integral of
strength 𝑥 = ∫ (𝐄 ∙ 𝐧)𝑑𝑠
𝑆
qfInt_Grad_n_ds
Name,
Formula and Description
ActiveFied constant
Current through a
surface 𝐼 = ∫ (𝐣 ∙ 𝐧)𝑑𝑠
𝑆
qfInt_KGrad_n_ds Electric current through a particular surface.
Potential difference
∆𝑈 = ∫ (𝐄 ∙ 𝐭)𝑑𝑙
qfInt_Grad_t_dl 𝐿
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
Surface integral of
strength 𝑥 = ∫ (𝐄 ∙ 𝐧)𝑑𝑠
𝑆
qfInt_Grad_n_ds
Name,
Formula and Description
ActiveFied constant
Apparent power
produced in a volume 𝑃𝐴𝑃𝑃 = ∫ (𝐄 ∙ 𝐣𝑨𝑷𝑷 )𝑑𝑣
𝑉
qfInt_PowerApparent Quadratic value
Apparent power produced in a particular
volume.
Mechanical force
1
qfInt_MaxwellForce 𝐅= ∮ (𝐄(𝐧𝐃) + 𝐃(𝐧𝐄)– 𝐧(𝐄𝐃))𝑑𝑠
2
𝑆
Quadratic vector
Electric force acting on bodies contained in
a particular volume. Evaluated by
calculating of Maxwell stress tensor over
volume’s bounding surface.
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.
Surface energy
𝑊𝑆 = ∫ (𝐄 ∙ 𝐃)𝑑𝑠
qfInt_GradKGrad_n_ds 𝑆
Quadratic value.
Calculating Integrals 143
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 volume 1
potential 𝑈𝑉 = ∫ 𝑈𝑑𝑣
𝑉 𝑉
qfInt_Potential_dv Complex value.
Average volume 1
displacement 𝐃𝑎 = ∫ 𝐃𝑑𝑣
𝑉 𝑉
qfInt_KGrad_dv Complex vector.
Average electric displacement vector in a
particular volume.
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.
144 Chapter 8 Analyzing Solution
Line integral of
displacement 𝑥 = ∫ (𝐃 ∙ 𝐭)𝑑𝑙
𝐿
qfInt_KGrad_t_dl Complex value.
Name, ActiveFied
Formula and Description
constant
Electrical charge
𝑞 = ∮ 𝐃 ∙ 𝐧 𝑑𝑠
qfInt_KGrad_n_ds
According to the Gauss theorem, total
electric charge in a particular volume can
be calculated as a flux of electric
displacement over its closed boundary.
Mechanical force
1
qfInt_MaxwellForce 𝐅= ∮ (𝐄(𝐧𝐃) + 𝐃(𝐧𝐄)– 𝐧(𝐄𝐃))𝑑𝑠
2
𝑆
Stored energy 1
𝑊= ∫ (𝐄 ∙ 𝐃)𝑑𝑣
qfInt_ElectrostaticEnergy 2 𝑉
Electric field energy in a particular volume.
Surface energy
𝑊𝑆 = ∫ (𝐄 ∙ 𝐃)𝑑𝑠
qfInt_GradKGrad_n_ds 𝑆
Potential difference
∆𝑈 = ∫ (𝐄 ∙ 𝐭)𝑑𝑙
qfInt_Grad_t_dl 𝐿
Average volume 1
potential 𝑈𝑉 = ∫ 𝑈𝑑𝑣
𝑉 𝑉
qfInt_Potential_dv
Average volume 1
displacement 𝐃𝑎 = ∫ 𝐃𝑑𝑣
𝑉 𝑉
146 Chapter 8 Analyzing Solution
Mean square 1
displacement 𝐷𝑎2 = ∫ 𝐷2 𝑑𝑣
𝑉 𝑉
qfInt_KGrad_dv
Line integral of
displacement 𝑥 = ∫ (𝐃 ∙ 𝐭)𝑑𝑙
𝐿
qfInt_KGrad_t_dl
Surface integral of
strength 𝑥 = ∫ (𝐄 ∙ 𝐧)𝑑𝑠
𝑆
qfInt_Grad_n_ds
Name,
Formula and Description
ActiveFied constant
Heat flux
𝛷 = ∫ (𝐅 ∙ 𝐧)𝑑𝑠
qfInt_KGrad_n_ds 𝑆
Temperature difference
Δ𝑇 = ∫ (𝐆 ∙ 𝐭)𝑑𝑙
qfInt_Grad_t_dl 𝐿
Average surface 1
temperature 𝑇𝑆 = ∫ 𝑇𝑑𝑠
𝑆 𝑆
qfInt_Potential_ds
Calculating Integrals 147
Average volume 1
temperature 𝑇𝑉 = ∫ 𝑇𝑑𝑣
𝑉 𝑉
qfInt_Potential_dv
Average volume 1
temperature gradient 𝐆𝑎 = ∫ 𝐆𝑑𝑣
𝑉 𝑉
qfInt_Grad_dv Mean vector of temperature gradient in a
volume.
Average volume 1
temperature gradient 𝐺𝑎2 = ∫ 𝐺 2 𝑑𝑣
𝑉 𝑉
qfInt_Grad2_dv
Surface integral of
grad(T) 𝑥 = ∫ (𝐆 ∙ 𝐧)𝑑𝑠
𝑆
qfInt_Grad_n_ds
Name,
Formula and Description
ActiveFied constant
Force
𝐅 = ∮ (𝜎 ∙ 𝐧) 𝑑𝑠
qfInt_Force 𝑆
Torque 1
𝐓= ∮ [𝐫 × (𝜎 ∙ 𝐧)]𝑑𝑠
qfInt_Torque 2 𝑆
Total torque of the forces acting on 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.
Lengthen
Δ𝐿 = ∫ (𝜎 ∙ 𝐭)𝑑𝑙
qfInt_lengthen 𝐿
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.
Data Tables 149
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 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.
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.
QuickField applies these settings to all columns at once. To change the
width of a single column drag the right border of the column in the table
header.
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:
150 Chapter 8 Analyzing Solution
When you enter the number of rows in Rows to Fill, QuickField interprets
this value as either the required number of rows in the whole table or the
required number of rows for every contour segment depending on the state
of the radio button (Total / Per Segment) below this field.
When QuickField tabulates the contour, it uses the data in the dialog items
you fill last.
You can combine automatic contour tabulation with manual addition and
removal of rows via the window’s context menu.
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 151
Plots of local field values can simultaneously show the curves for several
field points. At the same time, such plots can show the curves for several
physical quantities with the same measurement units. In cases when the
measurement units differ (as with temperature, its gradient, and the heat
flux), QuickField separates physical quantities into groups with same
measurement units and displays one group of quantities at a time.
The plots of integrated values, on the other hand, always display only one
physical quantity, with graphs of vector quantities represented as several
curves displaying the quantity's coordinate values and its absolute value.
To create a new time plot, choose Time Plot in the View or context menu.
The field value QuickField plots corresponds to the last point you clicked.
The clicks that count are those inside the field picture window and those
inside the Calculator Window.
If the last click is in the Integral Calculator the time plot displays the
clicked integrated quantity. Such plots display only one quantity at a time.
In other cases, the plot displays the local value at the clicked point against
time. If you did not click any point, or if the last point you clicked cannot
be associated with a field value, QuickField displays an empty time plot
window.
To display the curves related to different points on the same time plot,
click the points one by one invoking the Time Plot command after every
click. Or, to specify the exact point coordinates, invoke the
View / Time Plot Curves command available from the Time Plot
window, enter the coordinates, and click Add.
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.
152 Chapter 8 Analyzing Solution
The Choose Quantity Family drop-down list at the bottom allows you to
switch the displayed curve family. Alternatively you can change the
family via the context menu of the time plot or the legend window.
The time plot window immediately reflects all changes in the dialog.
To add a new point, click the very first row in the table, enter the
coordinates in the boxes above, and click Add. To change the coordinates
of the existing point select it in the list and enter the new values. This way
you can also switch the curves associated with points on and off.
The drop-down listbox in the upper left corner of the table window allows
one to toggle between local and integrated quantities. The default mode of
the window depends on whether its contour is defined or not. If no
contour exists, the window displays local quantities (Local values mode).
Otherwise it displays the quantities integrated over the current contour
(Integrals mode).
In Local values you can change the coordinates of the point. By default,
both coordinates are 0. Having changed the coordinates click OK to apply
the changes.
In Time Table window you can:
Copy the entire table or the selected rows to Windows Clipboard. To
do that, select the desired rows and choose Copy (CTRL+C) in the Edit
or context menu.
Drag the selected rows to another application, like Microsoft Excel.
When all rows in a table are selected, the table header is also included
into copying and drag-n-drop operations. To select all rows in the table
choose Select All (CTRL+A) in the Edit or context menu.
Save the entire table in a text file. Choose File / Save As.
We reorganize this system of three second degree equations into six first-
degree equations and append the following additional equation:
𝑑𝑙 𝑑𝑥 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.
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.
156 Chapter 8 Analyzing Solution
The page shows the trajectory length and the total flight time along with
the particle’s velocity and acceleration at the current trajectory point
(marked with cross on screen). The current point follows cursor position.
Export of Field Calculation Results 157
To print the picture, click Print in the File menu. You will have an
option to choose the printer and set up the picture, such as paper size
and orientation, before printing will occur.
To preview the output before printing, click Print Preview in the File
menu. To see how the picture will appear on a printer of your choice,
click Print Setup before.
2. Select needed File Type from the list, and set the name and location of
the file in the File Name field.
1. Click OK.
Export of Field Calculation Results 159
In the problem list on the left hand side choose the problem of interest.
The conductor list (right) shows labels of edges and vertices that can be
considered as electrodes. Icons close to each label show the geometrical
nature of the electrode (an edge or a vertex).
You could set or remove grounding from any of the conductors by a
mouse click in the Info column. Partial capacitances are not calculated for
grounded conductors with zero potential.
Conductor may be excluded from the calculations by a mouse click over
the symbol to the left from the corresponding label. Excluded conductors
are not involved in the capacitance matrix calculation, but unlike
grounded conductors, zero potential is not assigned to them.
When the conductor list is ready, click the Calculate button and wait for
results.
List of conductors shows the conductors with the numbers assigned in the
results window.
The energy calculation results for particular problems are then displayed,
first in single column, and then in the upper-diagonal matrix form.
The matrix of electrostatic induction coefficient is shown next, and then
the matrix of Self and mutual partial capacitances. Use the list of
conductors (above) to find the correspondence between the conductor
number and its label.
Any text from the results window can be copied to the clipboard ( CTRL+C
or the Copy command from the context menu).
163
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.
For more than 20 years QuickField is known as an efficient tool for 2D
simulation of electromagnetic, heat transfer and stress analysis problems.
Engineers and researchers use 2D analysis for fast evaluation and finding
the critical points in the model. It is actual now and will be used in the
foreseen future. However, sometimes 3D analysis is essential even for
very rough approximation of the field distribution and model parameters.
QuickField 6.0 is a first milestone on the way to 3D FEA simulation,
adding 3D extrusion to familiar process of 2D model creation, and
featuring full-scale 3D field analysis. In QuickField 6.1 working with 3D
geometry extended with importing geometry from most CAD program via
STEP format. For the present these 3D features are available for
electrostatic, DC conduction and steady state heat transfer analyses only.
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:
164 Chapter 9 Working with 3D Problem
For the problem with model class 3D Import import the geometry
information from the STEP file (ISO 10303), exported from almost
any 3D CAD software.
Build automatic or manually adjusted mesh of the tetrahedral finite
elements.
Assign the text labels to all geometrical objects – bodies, faces, edges
and vertices.
Define the physical properties, boundary conditions and field sources
for all labeled objects.
Start automatic solution procedure.
When solution is completed, 3D Postprocessor allows selecting of the
most impressive 3D field presentation, getting local results and
calculating integral parameters.
3D model comprises bodies with various physical properties and
surrounding space. Same as in 2D case, 3D model should always be
spatially confined: for open problems (where no physical boundary exists)
the artificial boundary should be created at sufficient distance from the
model objects.
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 approximation means that the geometry, boundary
conditions, field sources, material properties and thus the calculated field
depend on X and Y coordinates, but do not depend on the coordinate Z.
166 Chapter 9 Working with 3D Problem
Axisymmetric problems
Axisymmetric 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
Starting from the version 6.0 QuickField is able to solve 3D problems.
Right now it is applicable to the electrostatic, DC conduction and 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 specified height several times,
which allows making a complicated multi layered 3D body.
how to add the vertices and edges, can assign the labels to geometrical
objects, and able to control the finite-element mesh density using Spacing.
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.
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.
If the object (body, face, edge, vertex) label is selected in the problem
properties window, then all the respectively labeled objects will be
selected in the 3D view window. And vice versa – selection of the object
in the 3D view window causes its highlighting in the problem properties
window.
Majority of operations with 3D view is performed by the mouse cursor
movement while left or right mouse button is pressed. For some actions
one or more modifier keys should be pressed on the keyboard before the
mouse is moved.
Making the random transformations of the image in 3D view may easily
get you lost and lose control over the model view. But it is always easy to
return to one of predefined views using the camera control buttons on the
toolbar.
System of the coordinate axes shown at the bottom left corner of the
model window rotates together with the model and helps to understand the
current model orientation:
Image rotation, panning and zooming 175
Mouse movements listed above allow full control over the zoom level and
camera position relating the model. Besides that, there are several fixed
image settings which may be switched on using the toolbar buttons.
In conclusion here is the list of all methods of the image zoom control:
1. Mouse wheel rotation increases the zoom level in front direction
or decreases the level with backward direction. Screen position
of the model point under the mouse cursor stays intact.
2. After pressing the toolbar button Zoom in you may click by the
mouse in the model to increase the zoom level twice, or by
moving the mouse with left button pressed you may draw the
rectangle on the screen, which will be stretched to occupy the full
window when the mouse button is released.
3. Pressing the button Zoom out undoes the last zoom increase. If
there is nothing to undo the zoom level will be decreased twice.
4. Button Zoom to fit sets the zoom level necessary to fit the full
model into the window.
5. Mouse movement with the left button and keys SHIFT+CTRL
pressed increases the zoom level if directed up, or decreases it if
directed down.
6. Keyboard commands:
CTRL + (keys Ctrl and plus) increases the zoom level twice;
Image rotation, panning and zooming 177
CTRL - (keys Ctrl and minus) decreases the zoom level twice;
CTRL 0 (keys Ctrl and zero) scales the window to fit the full
model.
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.
Group selection using the rubber rectangle. This action requires the
right mouse slide.
In 3D view the body obscuring other bodies can be not only made hidden
but also rendered transparently. To do that use the Transparent command
in context menu. To cancel transparency of selected bodies use the
Opaque command in context menu.
Transparency does not provide an opportunity to select by mouse click the
body located behind a transparent body. To do this, you need to hide the
obscuring body or use a cutting plane, making it invisible.
Image rotation, panning and zooming 179
To turn the cutting plane on or off click on the button Cut by Plane in the
model editor 3D view toolbar, or use the command Cut by Plane in the
View menu.
It is possible to control the position of the cutting plane by the following
methods:
Dragging the cutting plane by the mouse along the plane normal;
Dragging normal vector to the cutting plane by the mouse;
Selecting one of the predefined plane normal orientations along the
coordinate axes in the context menu;
Inside out command of the context menu controls what part of the cut
model should be shown on the screen.
If some object cannot be viewed or selected behind the cutting plane – the
plane may be hidden temporarily using the command Hide cutting plane
in the context menu or in the View menu.
properties window with the current label of the selected object shown in
the corresponding subgroup (Bodies, Faces, Edges or Vertices).
If the view control button for the mesh of some dimension turns on
visibility of the mesh which is not yet built, it will be created
automatically. Volume mesh is easier to view using the model cut by
plane.
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 183
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:
Different ways of the field picture display:
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,
Vector field plotting using scaled directed vectors,
Scalar field plotting using small scattered color balls,
Equipotential surfaces,
Family of parallel slices, colored according to the value of
the chosen physical parameter,
Cut plane, which may be interactively controlled,
2D picture of the field view in the cut plane,
Plot of the chosen physical parameter along the selected
edge or edges.
Local field values in the chosen point. Point may be selected by
mouse or by entering the coordinates.
Integral parameters calculated across the volume of selected bodies or
along their external surface.
Besides this, the field view picture includes the 3D coordinate system
icon, which turns along with the camera movement relating to the model
and informs about the current orientation, and color map which shows the
field parameter levels in physical units.
Show geometry
Show Mesh
Colored surfaces
Vector plot
Scatter plot
186 Chapter 9 Working with 3D Problem
Isosurface plot
Slice plot
Cut pane
2D section
Pic. Field view and settings pane: model geometry, colored surfaces,
vector plot, scatter plot and plane cut. Selected presentations are shown by
red circles in the settings pane at the left.
Finite Element Mesh presentation shows the finite element mesh on the
visible surfaces. Finite Element Mesh command in the View menu or the
field view toolbar button turns this presentation on and off. It does not
have adjustment settings.
Colored surfaces
Colored surface presentation is turned on by default and displays the
distribution of the chosen physical parameter on the visible surfaces as a
color map, where every level of the physical parameter value corresponds
to some color. Field view toolbar switches the colored surface
presentation on and off.
Colored surface settings include:
1. Selection of the scalar physical quantity from the list. As a
default parameter for electrostatic problems the electric potential
V is used, for thermal problems – temperature T is used etc. List
of alternative parameters includes the field potential, modules
and coordinate components of vectors, and media properties
(dielectric permeability etc.).
2. Minimum and maximum values of the parameter for the legend
adjustment. Default range corresponds to the maximum and
minimum parameter values in the whole calculation area.
3. Number of color grades, with 256 colors as a default.
4. Whether the legend should be displayed or not.
module in chosen scale. The Vectors command in the View menu or the
button in the field view toolbar turn the vector plot on and off.
Pic. 2D field picture in the cut plane: geometry model, colored surface and
isolines.
Mouse gestures which control the field view are the same as in 3D case
with obvious exceptions: model cannot be rotated, cut plane is always
positioned in the screen plane.
Cut plane allows browsing the local field values similarly to how it could
be done in 3D window.
The tab Integral calculator in the 2D window provides calculation of the
integral values over the whole planar section, e.g. vector flux across the
section. Surface integrals are calculated over all the visible bodies cut by
the section plane.
Field distribution in the section plane may be also displayed as a plot or
table. Table and plot present the field along the broken line contour
consisting from one or more segments.
When the contour is ready, press one of the buttons in the toolbar to
activate the plot or table window. If the table or plot window is already
open, pressing the button 1D Plot or 1D Table changes the window
content accordingly to the current contour.
Contour for plots and tables
Contour for making a plot or table of the field distribution is a flat broken
line, built from one or more segments.
To start creating a contour click the button on the toolbox. The cursor
will change its shape to a cross-hair. Click mouse button to create the
initial node and then drag the mouse till the next node. Each left mouse
196 Chapter 9 Working with 3D Problem
click generates the new node. Right mouse click makes the last node. If
the last node coincides with the first node, the contour will be closed
polygon.
After creating the last contour node, it turns to the mode which allows
dragging existent nodes by the mouse or creation new nodes on the edges.
To pull the node you need to move the cursor to the node to change its
color, and drag the cursor with the right mouse button pressed. Doing the
same with the middle point of the contour will create a new intermediate
node.
When the mouse cursor is over the contour, right mouse button click
activates the context menu with contour editing commands:
Two upper commands select one of the alternative contour modes: adding
new segments to the end of the broken line, or dragging the existent
nodes. Next two commands remove the selected node or the whole
contour.
The contour is always directed, i.e. it is known which nodes are first
and last. Contour direction may be changed using the Change
Direction command.
Translate Contour command allows dragging the contour as a whole
by the mouse.
Context menu commands for the contour manipulation are also included
into the Contour menu.
Multiple parallel cuts – slice plot
It is often useful to use not one, but multiple cut planes for complete field
analysis. Presentation “Slice plot” is provided for such situations.
Command Parallel Slices from the View menu or the button in the
field view toolbar turns this presentation on and off.
Result analysis 197
The plot shows distribution of the selected physical quantity along the
selected edge or the sequence of edges. The plot is automatically redrawn
with the change of edge selection. Sequence of edges for plotting should
form the spatial contiguous not self-crossing. The edge selection suitable
for plotting should form a contiguous line without self-crossing segments.
Beginning of this line may coincide with the end – no other crossings are
allowed.
Toolbar above the plot area includes:
Zoom buttons. To the left from the plot the legend pane is
shown, which may be dragged to any side of the plot or left floating.
The legend pane allows selection of the physical parameter for
plotting. Radio buttons provide the choice between the parameter
families, expandable set of the flags turns on/off the specific
parameter within the family.
This button returns to the field view which is associated with this
plot. If the plot is drawn along the selected edges, the button activates
the 3D field view picture. If the plot is drawn along the flat contour,
the button activates the 2D section window, which includes the
contour.
This button activates the table window, which includes the
distribution of the physical parameters along the edges or 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.
To specify the number of the table rows you need to press the button in
the toolbar, which displays the dialogue window Table rows.
Number of the rows in the table may be defined directly, or by setting the
desired space between the nodes in the model length units measured along
the contour. If the contour consists from several segments, or there are
several edges selected, dialogue allows defining different steps for
different segments. All the segment ends are mandatory included into the
table.
In the upper part of the pane three fields X, Y and Z are prepared for
entering the probing point coordinates. These fields are automatically
refreshed when the mouse is clicked within the model.
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 are shown in the Integral Calculator tab of the 2D
section window.
In the electrostatic problems the following integrals may be calculated:
Result analysis 201
𝑞 = ∮ 𝐃 ∙ 𝐝𝐬
𝑆
w = 1/2 ∫ 𝐄𝐃 ⋅ 𝐝𝐯
𝑉
I = ∮ 𝐉 ∙ 𝐝𝐬
𝑆
Ф = ∮ 𝜆 ∙ grad 𝑇 𝐝𝐬
𝑆
C H A P T E R 1 0
Add-ins
For more information about writing your own add-ins, see Creating Add-
ins in ActiveField help.
Menu items attached to this add-in. This box contains the list of menu
items attached to the add-in. Typically, one menu item is attached to every
add-in. However, for more advanced add-ins, you can associate several
menu items with one add-in.
To add a menu item, click Add Item.
To delete a menu item, click Delete.
To change menu item properties, click Edit.
Description Properties
Friendly name. It is recommended to specify a readable name here for
add-in.
Description. It is recommended to provide detailed description of the add-
in here.
206 Chapter 10 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, Bmagnetic 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 𝜕𝐴 𝜕𝐻𝑐𝑦 𝜕𝐻𝑐𝑥
( )+ ( ) = −𝑗 + ( − );
𝜕𝑥 𝜇𝑦 𝜕𝑥 𝜕𝑦 𝜇𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑦
and for axisymmetric case is
𝜕 1 𝜕(𝑟𝐴) 𝜕 1 𝜕𝐴 𝜕𝐻𝑐𝑟 𝜕𝐻𝑐𝑧
( )+ ( ) = −𝑗 + ( − ).
𝜕𝑟 𝑟𝜇𝑧 𝜕𝑟 𝜕𝑧 𝜇𝑟 𝜕𝑧 𝜕𝑧 𝜕𝑟
where components of magnetic permeability tensor x and y (z and r),
components of coercive force vector Hcx and Hcy (Hcz and Hcr), and current
density j are constants within each block of the model.
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 density or total number
of ampere-turns associated with the block density associated with the
block. Current density in a coil can be obtained from the equation
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.
Magnetostatics 209
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.
Note. Zero Dirichlet condition is defaulted at the axis of rotation for the
axisymmetric problems.
210 Chapter 11 Theoretical Description
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:
𝐵
𝐵 = (𝐵) · (𝐻 + 𝐻𝑐 ); (𝐵) =
𝐻 + 𝐻𝑐
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.
Magnetostatics 211
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
𝑆
1
𝐓= ∮ ([𝐫 × 𝐇](𝐧𝐁) + [𝐫 × 𝐁](𝐧𝐇) – [𝐫 × 𝐧](𝐇𝐁))𝑑𝑠,
2
𝑆
1
𝑊= ∫ (𝐇 · 𝐁)𝑑𝑣 linear case;
2
𝑉
212 Chapter 11 Theoretical Description
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.
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
Transient Magnetics 213
each conductor. Note that the current should be turned on in all conductors
corresponding to one coil.
In plane-parallel case the inductance is calculated per unit length in
z-direction.
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 = − σ ∙ 𝐠𝐫𝐚𝐝𝑈,
𝜕𝐀
𝐣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
214 Chapter 11 Theoretical Description
𝜕 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 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.
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 either by the electric current density or total number
of ampere-turns associated with the block density associated with the
block. Current density in a coil can be obtained from the equation
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.
Several blocks with the same number of ampere-turns specified can be
considered as connected in series. In that case current density in each
Transient Magnetics 215
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.
.
Note. Zero Dirichlet condition is defaulted at the axis of rotation for the
axisymmetric problems.
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.
Transient Magnetics 217
1 𝜕(𝑟𝐴) 𝜕𝐴
𝐵𝑧 = ∙ , 𝐵𝑟 = − for axisymmetric case;
𝑟 𝜕𝑟 𝜕𝑧
Vector of magnetic field intensity
1
𝐇 = 𝐁,
𝜇
where is the magnetic permeability tensor.
Joule heat density
1 2
𝑄 = 𝑗
𝜎
Magnetic field energy density w = (BH)/2;
218 Chapter 11 Theoretical Description
𝐼 = ∫ 𝑗 𝑑𝑠
1
𝐅= ∮ (𝐇(𝐧𝐁) + 𝐁(𝐧𝐇)– 𝐧(𝐇𝐁))𝑑𝑠,
2
𝑆
1
𝐓= ∮ ([𝐫 × 𝐇](𝐧𝐁) + [𝐫 × 𝐁](𝐧𝐇) – [𝐫 × 𝐧](𝐇𝐁))𝑑𝑠,
2
𝑆
1
𝑊= ∫ (𝐇 · 𝐁)𝑑𝑣 linear case;
2
𝑉
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.
ellipse. The semi-major axis of the ellipse corresponds to the peak value
of the vector. The ratio between minor and major axes of the ellipse
defines the coefficient of polarization. The coefficient of polarization is
assumed to be positive for the counterclockwise and negative for the
clockwise rotation. Zero coefficient corresponds to the linear polarization.
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 oscillating magnetic field
𝐣 = 𝐣0 + 𝐣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 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.
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 represents a surface current in three-dimensional
world. It is specified by the Neumann boundary condition for the edge.
There are several ways to specify space-distributed electric current. In a
massive conductor, you can specify either a total current or a voltage
applied to the conductor. In planar problems, voltage drop is specified per
unit depth of the model, and in axisymmetric case voltage is assumed per
one turn around the axis of symmetry. Nonzero voltage applied to a
conductor in axisymmetric problem means that the conductor has a radial
cut, and the voltage is applied to sides of the cut. In practice this option
could be used to describe known voltage applied to massive spiral wiring,
in which case the total voltage drop for the coil should be divided by
number of turns in the coil.
Several blocks with the same value of total current or voltage applied can
be considered as connected in series. In that case each conductor carries
the same total current, and voltage (if any) is applied to the terminals of
the whole group of conductors connected in series.
Note. The meanings of zero total current and zero voltage applied to a
conductor are very different. Zero voltage means that the conductor’s ends
are short circuit, and zero value of the total current means open ends of the
conductor.
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.
Note. In order to properly model thin wire coils, the source current density
j0 in non-conductive areas is assumed to be uniform in both plane and
axisymmetric cases. Its behavior is different for massive conductors,
where source current density varies as ~(1/r) in axisymmetric case.
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,
Ht+ - Ht- = σ 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.
Note. Zero Dirichlet condition is defaulted at the axis of rotation for the
axisymmetric problems.
1 𝜕(𝑟𝐴) 𝜕𝐴
𝐵𝑧 = ∙ , 𝐵𝑟 = − for axisymmetric case;
𝑟 𝜕𝑟 𝜕𝑧
Complex vector of magnetic field intensity
1
𝐇 = 𝐁,
𝜇
where is the magnetic permeability tensor;
Time average and peak Joule heat density
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
𝐼 = ∫ 𝑗 𝑑𝑠
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
𝑆
1
𝐓= ∮ ([𝐫 × 𝐇](𝐧𝐁) + [𝐫 × 𝐁](𝐧𝐇) – [𝐫 × 𝐧](𝐇𝐁))𝑑𝑠,
2
𝑆
𝐅 = ∫ [𝐣 × 𝐁]𝑑𝑣.
𝑉
𝐓 = ∫[𝐫 × [𝐣 × 𝐁]]𝑑𝑣.
𝑉
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
226 Chapter 11 Theoretical Description
𝑉
𝐿= sin(φ𝑉 − φ𝐼 ),
𝐼 ∙ 2𝜋𝑓
To get mutual inductance between two conductors, you can specify non-
zero total current in one of them, make the ends of the other open
(applying zero total current), and measure the voltage induced in the
second conductor by the current in the first one.
Note. As in planar case, voltage is applied and measured per unit length,
the impedance is also calculated per unit length in z-direction.
Electrostatics
Electrostatic problems are described by the Poisson's equation for scalar
electric potential U (E = gradU, Eelectric 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.
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.
Electrostatics 227
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 from one entity to another, but have to be adjusted to
avoid discontinuities at entities' junction points.
𝜕𝑈 𝜕𝑈 for axisymmetric
𝐸𝑧 = − , 𝐸𝑟 = −
𝜕𝑧 𝜕𝑟 case;
Tensor of the gradient of electric field intensity G = gradE
𝜕𝐸𝑥 𝐸𝑦 𝜕𝐸𝑥 𝜕𝐸𝑦 —for planar
G𝑥𝑥 = , G𝑦𝑦 = , 𝐺𝑥𝑦 = ½ ( + )
𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑥 case;
𝑞 = ∮ 𝐃 ∙ 𝐧 𝑑𝑠,
1
𝐅= ∮ (𝐄(𝐧𝐃) + 𝐃(𝐧𝐄)– 𝐧(𝐄𝐃))𝑑𝑠,
2
𝑆
1
𝐓= ∮ ([𝐫 × 𝐄](𝐧𝐃) + [𝐫 × 𝐃](𝐧𝐄) – [𝐫 × 𝐧](𝐄𝐃))𝑑𝑠,
2
𝑆
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 capacitance of a conductor, put constant
potential boundary condition on its surface, specify an arbitrary non zero
electric charge in one of the vertices on the surface of the conductor (in
fact, the charge will be distributed over the conductor’s surface), and turn
off all other field sources in the model. Once the problem is solved, go to
the Postprocessor and take the value of electric potential somewhere on
the surface of the conductor. The capacitance of the conductor can be
obtained from the equation
𝑞
𝐶= ,
𝑈
where q is the electric charge and U is the potential of the conductor.
To calculate mutual capacitance between two conductors put a charge on
one conductor and measure electric potential on another. Constant
potential boundary condition has to be applied to the surfaces of both
conductors.
𝑞1
𝐶12 = .
𝑈2
Self- and mutual partial capacitance matrix calculation in the multi-
conductor system is discussed in the “Partial Capacitance Matrix
Calculation for the System of Conductors” chapter.
230 Chapter 11 Theoretical Description
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.
The equation for planar case is
𝜕 𝜕𝑈 𝜕 𝜕𝑈
(𝜎 )+ (𝜎 ) = 0,
𝜕𝑥 𝑥 𝜕𝑥 𝜕𝑦 𝑦 𝜕𝑦
for axisymmetric case:
1 𝜕 𝜕𝑈 𝜕 𝜕𝑈
(𝜎 𝑟 ) + (𝜎𝑧 ) = 0,
𝑟 𝜕𝑟 𝑟 𝜕𝑟 𝜕𝑧 𝜕𝑧
and for 3D formulation:
𝜕 𝜕𝑈 𝜕 𝜕𝑈 𝜕 𝜕𝑈
(𝜎 ) + (𝜎 ) + (𝜎 ) = 0,
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧
where components of electric conductivity tensor σx, σy or σz, σr are
constant within each model block. Anisotropy is not accounted for in 3D
formulation.
The electric current density j can be obtained from the equation
j σgrad U.
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 density by the relationship I = 2r.
Edge-bound current density in the plane of model represents a
surface-bound external current density in three-dimensional world. It is
specified by the Neumann boundary condition for the edge.
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
DC Conduction Analysis 231
Note. The edge described as possessing constant potential should not have
common points with any Dirichlet edge. In that case the constant potential
edge has to be described by the Dirichlet condition with an appropriate
potential value.
𝜕𝑈 𝜕𝑈
𝐸𝑧 = − , 𝐸𝑟 = − —for axisymmetric case;
𝜕𝑧 𝜕𝑟
Vector of current density j σE, where σ - is electric conductivity
tensor.
Integral quantities:
Electric current through a given surface
𝐼 = ∫ 𝐣 ∙ 𝐧 𝑑𝑠,
𝑆
𝑊 = ∫(𝐄 · 𝐣)𝑑𝑣.
𝑉
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.
Similar to problems of AC magnetics, variation of the field with respect to
time is assumed to be sinusoidal. All field components and electric current
vary with time like
z = z0 cos(t + z),
where z0 is a peak value of z, z — its phase angle, and — the angular
frequency.
Complex arithmetic is used to represent harmonic time-dependency.
The problem formulation combines equations for electrostatics (E = )
and dc conduction (j = -iσ), taking into account the Ohm's law, j = σE.
Final equation with respect to electric potential U is
𝑖𝜎
∇ ([𝜀 − ] ∇𝑈) = 0;
𝜔
where electric conductivity σ and components of electric permittivity
tensor z and y (z and r) are constants within each block of the model.
AC Conduction Analysis 233
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 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 density by the relationship
I = 2r. Edge-bound current density in the plane of model represents a
surface-bound external current density in three-dimensional world. It is
specified by the Neumann boundary condition for the edge. Every
condition is defined by its magnitude and phase.
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.
Note. The edge described as possessing constant potential should not have
common points with any Dirichlet edge. In that case the constant potential
edge has to be described by the Dirichlet condition with an appropriate
potential value.
𝜕𝑈 𝜕𝑈
𝐸𝑧 = − , 𝐸𝑟 = − 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
𝐼 = ∫ (𝐣 ∙ 𝐧) 𝑑𝑠
𝑆
𝑃 = ∫ (𝐄 ∙ 𝐣) 𝑑𝑣
𝑉
Transient Electric Analysis 235
1
𝐅= ∮ (𝐄(𝐧𝐃) + 𝐃(𝐧𝐄)– 𝐧(𝐄𝐃))𝑑𝑠,
2
𝑆
1
𝐓= ∮ ([𝐫 × 𝐄](𝐧𝐃) + [𝐫 × 𝐃](𝐧𝐄) – [𝐫 × 𝐧](𝐄𝐃))𝑑𝑠,
2
𝑆
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
236 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.
Note. The edge described as possessing constant potential should not have
common points with any Dirichlet edge. In that case the constant potential
edge has to be described by the Dirichlet condition with an appropriate
potential value.
𝐼 = ∫ 𝐣 ∙ 𝐧 𝑑𝑠,
𝑆
𝑃 = ∫(𝐄 · 𝐣)𝑑𝑣.
𝑉
1
𝐅= ∮ (𝐄(𝐧𝐃) + 𝐃(𝐧𝐄)– 𝐧(𝐄𝐃))𝑑𝑠,
2
𝑆
1
𝐓= ∮ ([𝐫 × 𝐄](𝐧𝐃) + [𝐫 × 𝐃](𝐧𝐄) – [𝐫 × 𝐧](𝐄𝐃))𝑑𝑠,
2
𝑆
238 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;
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧
𝜕 𝜕𝑇 𝜕 𝜕𝑇
(𝜆(𝑇) ) + (𝜆(𝑇) ) =
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦
𝜕𝑇 planar case;
= −𝑞(𝑇) − 𝑐(𝑇)𝜌
𝜕𝑡
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);
Heat Transfer 239
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 the xy-plane corresponds to a linear
string-like heater which is perpendicular to the plane of the model, and is
described by the generated power per unit length. In axisymmetric case
the vertex heat source represents a heating circle around the axis of
symmetry or a point heater located on the axis. To incorporate both these
cases a total generated power value is associated with the vertex. For the
heating circle the total power is connected with its linear density by the
relationship q 2rql. Edge-bound heat source in the plane of model
represents a surface heat source in three-dimensional world. It is described
by power per unit area and is specified by the Neumann boundary
condition for the edge. The volume power density associated with a block
corresponds to the volume heat source.
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 flux boundary condition (known also as boundary condition of the
second kind) is defined by the following equations:
Fn = -qs at outward boundaries,
+
Fn - Fn- = -qs at inner boundaries,
where Fn is a normal component of heat flux density, "+" and ""
superscripts denote quantities to the left and to the right side of the
boundary. For inner boundary qs denotes the generated power per unit
240 Chapter 11 Theoretical Description
area, for outward boundary it specifies the known value of the heat flux
density across the boundary. If qs value is zero, the boundary condition is
called homogeneous. The homogeneous condition at the outward
boundary indicates vanishing of the heat flux across the surface. This type
of boundary condition is the natural one, it is defaulted at all outward
boundary parts where no explicit boundary condition is specified. This
kind of boundary condition is used at an outward boundary of the region
which is formed by the symmetry plane of the problem.
If the surface heat source is to be specified at the plane of problem
symmetry and this plane constitutes the outward boundary of the region,
the surface power has to be halved.
This boundary condition sometimes is called the boundary condition of
the second kind.
Convection boundary condition can be specified at outward boundary of
the region. It describes convective heat transfer and is defined by the
following equation:
Fn = (T - T0),
where is a film coefficient, and T0temperature of contacting fluid
medium. Parameters and T0 may differ from part to part of the
boundary.
This boundary condition sometimes is called the boundary condition of
the third kind.
Radiation boundary condition can be specified at outward boundary of
the region. It describes radiative heat transfer and is defined by the
following equation:
Fn = kSB(T4 - T04),
where kSB is a Stephan-Boltsman constant, is an emissivity coefficient,
and T0—ambient radiation temperature. Parameters and T0 may differ
from part to part of the boundary.
𝜕𝑇 𝜕𝑇 for axisymmetric
𝐹𝑧 = −𝜆𝑧 , 𝐹𝑟 = −𝜆𝑟
𝜕𝑧 𝜕𝑟 case;
𝜕𝑇 𝜕𝑇
𝐹𝑥 = −𝜆 , 𝐹𝑦 = −𝜆 ,
𝜕𝑥 𝜕𝑦
𝜕𝑇 for 3D case;
𝐹𝑧 = −𝜆
𝜕𝑧
The postprocessor can calculate the heat flux through an arbitrary closed
or unclosed surface
𝛷 = ∫ (𝐅 ∙ 𝐧) 𝑑𝑠
𝑆
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.
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.
242 Chapter 11 Theoretical Description
𝛿𝑥
{𝛿} = { } for plane problems;
𝛿𝑦
𝛿
{𝛿} = { 𝑧 } for axisymmetric problems.
𝛿𝑟
Only three components of strain and stress tensors are independent in both
plane stress and plane strain cases. The strain-displacement relationship is
defined as:
𝜕𝛿𝑥
𝜕𝑥
𝜀𝑥
𝜕𝛿𝑦
{𝜀} = { 𝜀𝑦 } =
𝛾𝑥𝑦 𝜕𝑦
𝜕𝛿𝑥 𝜕𝛿𝑦
+
{ 𝜕𝑦 𝜕𝑥 }
The corresponding stress components:
𝜎𝑥
{𝜎} = { 𝜎𝑦 }
𝜏𝑥𝑦
The axisymmetric problem formulation also includes the out-of-plane
direct strain , caused by the radial deformation. The strain-displacement
relationship is defined as:
𝜕𝛿𝑧
𝜕𝑧
𝜀𝑧 𝜕𝛿𝑟
𝜀 𝜕𝑟
{𝜀} = { 𝜀𝑟 } =
𝜃 𝛿𝑟
𝛾𝑟𝑧 𝑟
𝜕𝛿𝑧 𝜕𝛿𝑟
{ 𝜕𝑟 + 𝜕𝑧 }
The corresponding stress components:
𝜎𝑧
𝜎
{𝜎} = { 𝜎𝑟 }
𝜃
𝜏𝑟𝑧
The equilibrium equations for the plane problems are:
𝜕𝜎𝑥 𝜕𝜏𝑥𝑦
+ = −𝑓𝑥
𝜕𝑥 𝜕𝑦
𝜕𝜏𝑥𝑦 𝜕𝜎𝑦
+ = −𝑓𝑦
{ 𝜕𝑥 𝜕𝑦
Stress Analysis 243
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
Stress Analysis 245
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.
Note. The concentrated load is specified by the force per thickness unit
for plane problems and by the total force value for axisymmetric ones. In
the last case the force can be applied to the point at the axis of symmetry
or distributed along the circle around the axis.
Any surface force which is directed along the normal to the surface can be
described as a pressure. The pressure is considered positive if it is directed
inside region at its outward boundary or from right to left at the inner
boundary. Left and right are referred relative to the edge intrinsic
direction, which is always counterclockwise for arcs and is determined for
line segments by the order of picking vertices when the edge is created.
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.
246 Chapter 11 Theoretical Description
1
σ𝑒 = √ [(𝜎1 − 𝜎2 )2 + (𝜎2 − 𝜎3 )2 + (𝜎3 − 𝜎1 )2 ]
2
where
1 𝜎1 + 𝜎2 + 𝜎3
σ𝑖 = √ [(𝜎1 − 𝜎2 )2 + (𝜎2 − 𝜎3 )2 + (𝜎3 − 𝜎1 )2 ]; 𝜎̅ =
2 3
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.
Integral quantities:
Total force acting on a particular volume
𝐅 = ∮ (𝜎 ∙ 𝐧)𝑑𝑠
𝑆
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:
248 Chapter 11 Theoretical Description
Note. In spite of the requirement that coupled problems must use the same
model file, the geometrical region for the problems need not coincide, i.e.
some subregions those are in use in one problem, could be excluded from
consideration in the other one.
Index
.
finite element mesh 187
integral values 183, 200
local values 183, 199
.m3d (model) files 27
objects hiding 186
result analysis 183
2 results 184
results window 183
2D and 3D problems 165 results, colored surfaces 189
results, isosurface plot 191
3
results, scatter plot 190
results, vector plot 189
slice plot 196
3D
solving 182
2D sketch 164
XY plot 197
extrusion 164, 166, 168
geometry model 167
import 166
level 168, 169
A
level, add 170 ac conduction 232
level, change 170 AC conduction 22
level, remove 170 data for 76
problem 163 ac magnetic 113, 219
problem definition 166, 167 AC magnetic 20
view manipulation 175, 186 data for 67
workflow 163 adaptive meshing 104
3D model add-in 160
boundary conditions 182 analysis capabilities 18
camera 176 animation 115
cutting by plane 179 attraction distance 39
import 171 axisymmetric problem 207, 209, 215, 216,
labels 179 222
mesh spacing 180
object hiding 178
object selection 177 B
view 171, 172
background grid 57
zooming 176
beam of charged particles 154
3D problem
B-H curve 18, 19, 64, 208, 214
cut plane 192, 194
block 37, 207
cut plane normal direction 193
and dc conduction data 74
electric charge 201
and electric charge 226
electric force 201
and electric current 208, 214
energy of electric field 201
and electric data 78
field calculator 183
and electrostatic data 72
field presentation modes 184
252 Index
C
creating
circuit 98
edge 39
CAD systems 58
label 63
calculated quantities 106
model 38
calculator 115, 119
problem 28
capacitance 21, 119, 122
vertex 39
calculation 229
criterion
wizard 122
Drucker-Prager 83, 110, 246
capacitance matrix 160
Hill 83
charge
Mohr-Coulomb 83, 110, 246
electric 21, 122, 128, 226, 228, 229
Tresca 110, 246
circuit 97
Tsai-Hill 110, 247
element 98
von Mises 110, 246
coefficient
current 128, 208, 221
of convection 240
active 109, 234, 237
of emissivity 240
alternating 20, 22, 219, 232
of polarization 220
apparent 109, 234, 237
of thermal expansion 244, 245
density 19, 21, 22, 63, 65, 67, 108, 208,
coercive force See force: coercive
220, 230
coil
density, units of 33
spiral 209, 215, 221
density, volume 220
color legend 153
eddy 20, 106, 107, 213, 217, 220, 223
color map 111
electric 22, 213, 220
Index 253
D
deleting 42
direction of 245
labeling 52
data
editing
ac conduction 76
circuit 100
ac magnetics 67
contour 124
dc conduction 74
curve 85
electrostatics 72
formula 86
heat transfer 79
label data 63
magnetostatics 63
model 37–60
stress analysis 81
problem 28
transient electric 77
electric
transient magnetics 63
boundary condition 227, 230
dc conduction 230
electric conductivity 74, 235
DC conduction 22
electric motor 18, 19, 20
data for 74
electric permittivity 21, 72, 76, 78, 108, 109,
dcf(dc conduction data) files 27
226, 228, 232, 234, 235, 237
dec(ac conduction data) files 27
electrostatic analysis 21, 226
deformed shape 111
data for 72
des (electrostatic data) files 27
emissivity coefficient 240
describing
emitter 154
geometry 38
energy
dhe (ac magnetic data) files 27
of deformation 110
dht (steady-state and transient heat transfer
of electric field 21, 108, 122, 229
data) files 27
of magnetic field 19, 20, 21, 106, 108,
dielectric constant See electric permittivity
120, 212, 217, 218, 223, 224
displacement 25, 110, 113, 242, 246
exporting
boundary condition 83, 84
DXF file 58
components 242
pictures 59
prescribed 25
units of 33
dms (magnetostatic and transient magnetic
data) files 27
F
dsa (stress analysis data) files 27 field 52
dtv (transient electric data) files 27 electric 213, 235
DXF file electric, transient 23
export 58 lines 111
import 58 local data 116
magnetic 208, 210, 213, 216, 217, 222
E
magnetic, uniform external 19, 21, 209,
215, 222
picture 111, 114
edge 37, 38, 207, 210, 223
presentation methods 111
and ac conduction data 77
source 38, 52, 61, 128, 207, 208, 214,
and ac magnetic data 71
221, 226, 230, 233, 235
and boundary condition 209, 215, 222,
film coefficient See coefficient: of convection
227, 230, 233, 236, 239, 245
floating conductor 21
and dc conduction data 75
254 Index
J
editing 86
frequency 219, 221
angular 219, 232
Joule heat 20, 106, 132, 134, 217, 218, 223,
224, 232, See power: losses
G
geometric objects L
copying 41
label 37, 38, 103
moving 41
assigning data to 63
selecting 40
assigning to geometric objects 52
geometry
copying data 96
describing 38
deleting data 96
editing data 63
H renaming 96
legend 153
hardware requirements 11 loading See also field source
heat AC magnetic analysis 21
specific 24 DC conduction analysis 22
heat flow 110 electrostatic analysis 21
heat transfer 20, 21, 24, 238 force 25, 245
coupling 29 heat transfer analysis 24
data for 79 magnetostatic analysis 19
transient 24 pressure 25, 245
thermal 25, 82, 114, 244
I
local field data 116
copying to clipboard 116
impedance 20, 21, 119, 124, 225, See also
inductance
wizard 124
M
importing magnetic permeability 19, 20, 63, 64, 65, 67,
data 29 68, 106, 107, 108, 208, 210, 211, 214, 217,
DXF file 58 218, 220, 223
electric forces to stress analysis 249 magnetostatic analysis 18, 208
Joule heat to heat transfer analysis 248 data for 63
magnetic forces to stress analysis 249 manipulating
Index 255
O
properties window
problem 18, 34
objects
manipulating 40
ohmic losses See power: losses
R
radiation 81, 240
P
res (results) files 27
resistivity
electric 22, 23, 108, 232
particle 154
particles trajectory 154
password 13
pbm (problem description) files 27
S
peak value 219, 232 selecting
periodic boundary conditions 84 circuit elements 99
permanent magnet 18, 19, 63, 208, 210, 213, geometric objects 40
214, 216, 217, 221 vertices 54
equivalent current 210, 216 serial conductor 69
phase angle 219, 232 shear modulus 244
physical quantities 211, 217, 223, 231, 234, solenoid 18, 20
237 source
heat transfer 241 thermal 21, 22, 23, 239, 240
stress analysis 246 specific heat 239
picture spiral coil 215
copying 59, 158 Stephan-Boltsman constant 240
exporting 59 strain 110
printing 157 components 242
picture, of the field 112 plane 25, 207
plane strain See strain: plane thermal 25, 243, 244
plane stress 246, See stress: plane stress 110, 113
Poisson's ratio 244 allowable 110, 246
polarization axisymmetric problem 241
256 Index
T
deleting 42
labeling 52
voltage 20, 21, 22, 67, 214, 221, See
tables
potential: electric
in postprocessor 114, 148
volume power 79
versus time 152
temperature 24, 29, 80, 85, 110, 238
ambient, for convection 240
ambient, for radiation 240
W
curve 80 window
difference 114 document 17
terminology 17, 37 problem 17
three-dimensional problem See 3D:problem tools 18
torque 128 wizard 119
of electric forces 21, 22, 228, 235, 237 capacitance 22, 122
of magnetic forces 18, 19, 20, 211, 218, impedance 21, 124
224, 225 inductance 19, 21, 120
transformer 19
X
transient analysis
electric 235
transient analysis
X-Y plot 114, 126
heat transfer 24
copying 158
data for 79
printing 157
initial condition 29
versus time 150
magnetic 19, 29, 213
data for 63
plot versus time 150
time layer selecting 114
Y
timing 31, 32 Young's modulus 244
Young's modulus 82
U
units of length 33
Z
zooming 56, 114
V
vertex 37, 38, 207, 208, 221
and ac conduction data 77