Finite Element Method Magnetics: Octavefemm: User'S Manual

Finite Element Method Magnetics: OctaveFEMM
Version 1.2 Users Manual

October 2, 2006
David Meeker [email protected] http://femm.foster-miller.com c 2006
1 Introduction
OctaveFEMM is a Matlab toolbox that allows for the operation of Finite Element Method Magnetics (FEMM) via a set of Matlab functions. The toolbox works with Octave, a Matlab clone. Since Octave does not support ActiveX, the toolbox use an alternative interprocess communication scheme based on temporary les for passing messages between Octave and FEMM. This toolbox can also run in Matlab using Matlabs built-in support of ActiveX. The use and functionality of the toolbox is essentially identical under either Matlab or Octave. The only difference is speed Matlabs ActiveX interface is noticeably faster than the le-based interprocess communications method used when running with Matlab. When OctaveFEMM starts up a FEMM process, the usual FEMM user interface is displayed and is fully functional. The user then has the choice of accomplishing modeling and analysis tasks either exclusively through functions implemented by the toolbox, or by a combination of manual and programmatic operations whichever is easiest for the task at hand. The syntax of the OctaveFEMM toolbox closely mirrors that of FEMMs existing Lua scripting language interface associated with FEMM v4.2. However, there are some differences between the Lua functions and the analogous Octave/Matlab implementations: All strings are enclosed in single quotes, rather than double quotes as in Lua. Functions in Lua that have no arguments require a set of empty parentheses after the function name (e.g. mi analyze()). In Octave or Matlab, no parentheses should be used needed (e.g. mi analyze with the OctaveFEMM toolbox). Several commands have also been added to OctaveFEMM that have no analog in Lua. These commands streamline the drawing of new geometries with the OctaveFEMM toolbox, as well is the collection of data from solutions. Perhaps the most remarkable difference between Lua and OctaveFEMM, however, is due to the matrix-oriented nature of Octave/Matlab. In just about any OctaveFEMM function in which it would be desiable to enter an array of points such that multiple copies of an operation are performed, OctaveFEMM will correctly interpret the input perform the requested operation on every element in the array. In addition, for any function in which the coordinates of a point are required, that point can be specied as an array with two elements instead of specifying each element separately. In functions that require the specication of multiple points, those points can be entered as an array of two-element arrays.
2 Installation and Startup

The OctaveFEMM distribution is automatically installed with FEMM 4.2 in the mfiles subdirectory. The absolute directory is typically c:\Program Files\femm42\mfiles. To use OctaveFEMM with Matlab, c:\Program Files\femm42\mfiles must be added to the Matlab search path. This can be accomplished via Matlabs interactive pathtool command. For use with Octave, you need to add the OctaveFEMM Toolbox to Octaves search path. This can be done permanently by either adding to or creating the .octaverc le that octave runs when it starts up. For example, in the Octave 2.1.50 release, Octave looks for .octaverc in the directory: 2
C:\Program Files\GNU Octave 2.1.50\octave_files To insert OctaveFEMM into the search path, one can create the .octaverc file and add the line: LOADPATH = [ /cygdrive/c/progra1/femm42/mfiles//, LOADPATH ]; to add the directory containing the OctaveFEMM commands to the Octave search path. Next, you need to may need to edit the command used to invoke FEMM, especially if you have installed FEMM in a directory that is not the default one. Edit c:\Program Files\femm42\mfiles\openfemm.m to make sure that OctaveFEMM is pointing to the right FEMM executable. The temporary les dened in this directory should also be located in the same directory as the FEMM executable, so you may need to edit those denitions as well. To start an OctaveFEMM session, use the openfemm function. This function starts up a FEMM process to which OctaveFEMM calls are sent. When done with OctaveFEMM, the FEMM process can be shut down via the closefemm function. A number of examples that use OctaveFEMM to analyze various problems are included in the directory cd c:\Program Files\femm42\examples
3 Common Command Set

There are a number of FEMM-specic Octave that are not associated with any particular problem type. These functions manipulate the appearance of the main window and other top-level components like the Lua console and Point Properties output window. clearconsole Clears the output window of the Lua console. newdocument(doctype) Creates a new preprocessor document and opens up a new preprocessor window. Specify doctype to be 0 for a magnetics problem, 1 for an electrostatics problem, 2 for a heat ow problem, or 3 for a current ow problem. Alternative syntax for this function is create(doctype) hideconsole Hides the oating Lua console window. hidepointprops Hides the oating FEMM Properties display window. messagebox(message) displays the message string to the screen in a pop-up message box. opendocument(filename) Opens a document specied by filename. print(item1,item2,...) This is standard Lua print command directed to the output of the Lua console window. Any number of comma-separated items can be printed at once via the print command. prompt(message) This function allows a FEMM script to prompt a user for input. When this command is used, a dialog box pops up with the message string on the title bar of the dialog box. The user can enter in a single line of input via the dialog box. Prompt returns the users input to Octave and parses it using Octaves eval command. 3
showconsole Displays the oating Lua console window. showpointprops Displays the oating FEMM Properties display window. main minimize minimizes the main FEMM window. main maximize maximizes the main FEMM window. main restore restores the main FEMM window from a minimized or maximized state. main resize(width,height) resizes the main FEMM window client area to width height.
4 Magnetics Preprocessor Command Set

A number of different commands are available in the preprocessor.
4.1 Object Add/Remove Commands

mi addnode(x,y) Add a new node at x,y mi addsegment(x1,y1,x2,y2) Add a new line segment from node closest to (x1,y1) to node closest to (x2,y2) mi addblocklabel(x,y) Add a new block label at (x,y) mi addarc(x1,y1,x2,y2,angle,maxseg) Add a new arc segment from the nearest node to (x1,y1) to the nearest node to (x2,y2) with angle angle divided into maxseg segments. mi drawline(x1,y1,x2,y2) Adds nodes at (x1,y1) and (x2,y2) and adds a line between the nodes. mi drawpolyline([x1,y1;x2,y2...]) Adds nodes at each of the specied points and connects them with segments. mi drawpolygon([x1,y1;x2,y2...]) Adds nodes at each of the specied points and connects them with segments to form a closed contour. mi drawarc(x1,y1,x2,y2,angle,maxseg) Adds nodes at (x1,y1) and (x2,y2) and adds an arc of the specied angle and discretization connecting the nodes. mi drawrectangle(x1,y1,x2,y2) Adds nodes at the corners of a rectangle dened by the points (x1,y1) and (x2,y2), then adds segments connecting the corners of the rectangle. mi deleteselected Delete all selected objects. mi deleteselectednodes Delete selected nodes. mi deleteselectedlabels Delete selected block labels. 4
mi deleteselectedsegments Delete selected segments. mi deleteselectedarcsegments Delete selects arcs.
4.2 Geometry Selection Commands

mi clearselected Clear all selected nodes, blocks, segments and arc segments. mi selectsegment(x,y) Select the line segment closest to (x,y) mi selectnode(x,y) Select the node closest to (x,y). Returns the coordinates of the selected node. mi selectlabel(x,y) Select the label closet to (x,y). Returns the coordinates of the selected label. mi selectarcsegment(x,y) Select the arc segment closest to (x,y) mi selectgroup(n) Select the nth group of nodes, segments, arc segments and blocklabels. This function will clear all previously selected elements and leave the editmode in 4 (group)
4.3 Object Labeling Commands

mi setnodeprop(propname,groupno) Set the selected nodes to have the nodal property propname and group number groupno. mi setblockprop(blockname, automesh, meshsize, incircuit, magdir, group, turns) Set the selected block labels to have the properties: Block property blockname. automesh: 0 = mesher defers to mesh size constraint dened in meshsize, 1 = mesher automatically chooses the mesh density. meshsize: size constraint on the mesh in the block marked by this label. Block is a member of the circuit named incircuit The magnetization is directed along an angle in measured in degrees denoted by the parameter magdir A member of group number group The number of turns associated with this label is denoted by turns. mi setsegmentprop(propname, elementsize, automesh, hide, group) Set the selected segments to have: Boundary property propname Local element size along segment no greater than elementsize automesh: 0 = mesher defers to the element constraint dened by elementsize, 1 = mesher automatically chooses mesh size along the selected segments 5
hide: 0 = not hidden in post-processor, 1 == hidden in post processor A member of group number group mi setarcsegmentprop(maxsegdeg, propname, hide, group) Set the selected arc segments to: Meshed with elements that span at most maxsegdeg degrees per element Boundary property propname hide: 0 = not hidden in post-processor, 1 == hidden in post processor A member of group number group
4.4 Problem Commands

mi probdef(freq,units,type,precision,depth,minangle) changes the problem definition. Set freq to the desired frequency in Hertz. The units parameter species the units used for measuring length in the problem domain. Valid units entries are inches, millimeters, centimeters, mils, meters, and micrometers. Set the parameter problemtype to planar for a 2-D planar problem, or to axi for an axisymmetric problem. The precision parameter dictates the precision required by the solver. For example, entering 1E-8 requires the RMS of the residual to be less than 108 . A fth parameter, representing the depth of the problem in the into-the-page direction for 2-D planar problems. Specify the depth to be zero for axisymmetric problems. The sixth parameter represents the minimum angle constraint sent to the mesh generator 30 degrees is the usual choice for this parameter. mi analyze(flag) runs the magnetics solver. The flag parameter controls whether the solver window is visible or minimized. For a visible window, either specify no value for flag or specify 0. For a minimized window, flag should be set to 1. mi loadsolution loads and displays the solution corresponding to the current geometry. mi setfocus(documentname) Switches the magnetics input le upon which commands are to act. If more than one magnetics input le is being edited at a time, this command can be used to switch between les so that the mutiple les can be operated upon programmatically. documentname should contain the name of the desired document as it appears on the windows title bar. mi saveas(filename) saves the le with name filename.
4.5 Mesh Commands

mi createmesh runs triangle to create a mesh. Note that this is not a necessary precursor of performing an analysis, as mi analyze will make sure the mesh is up to date before running an analysis. The number of elements in the mesh is pushed back onto the lua stack. mi showmesh shows the mesh. 6
mi purgemesh clears the mesh out of both the screen and memory.
4.6 Editing Commands

mi copyrotate(bx, by, angle, copies ) bx, by base point for rotation angle angle by which the selected objects are incrementally shifted to make each copy. angle is measured in degrees. copies number of copies to be produced from the selected objects. mi copyrotate2(bx, by, angle, copies, editaction ) bx, by base point for rotation angle angle by which the selected objects are incrementally shifted to make each copy. angle is measured in degrees. copies number of copies to be produced from the selected objects. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group mi copytranslate(dx, dy, copies) dx,dy distance by which the selected objects are incrementally shifted. copies number of copies to be produced from the selected objects. mi copytranslate2(dx, dy, copies, editaction) dx,dy distance by which the selected objects are incrementally shifted. copies number of copies to be produced from the selected objects. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group mi createradius(x,y,r)turnsacornerlocatedat(x,y)intoacurveofradiusr. mi moverotate(bx,by,shiftangle) bx, by base point for rotation shiftangle angle in degrees by which the selected objects are rotated. mi moverotate2(bx,by,shiftangle, editaction) bx, by base point for rotation shiftangle angle in degrees by which the selected objects are rotated. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group mi movetranslate(dx,dy) dx,dy distance by which the selected objects are shifted. 7
mi movetranslate2(dx,dy,editaction) dx,dy distance by which the selected objects are shifted. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group mi scale(bx,by,scalefactor) bx, by base point for scaling scalefactor a multiplier that determines how much the selected objects are scaled mi scale2(bx,by,scalefactor,editaction) bx, by base point for scaling scalefactor a multiplier that determines how much the selected objects are scaled editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group mi mirror(x1,y1,x2,y2) mirror the selected objects about a line passing through the points (x1,y1) and (x2,y2). mi mirror2(x1,y1,x2,y2,editaction) mirror the selected objects about a line passing through the points (x1,y1) and (x2,y2). Valid editaction entries are 0 for nodes, 1 for lines (segments), 2 for block labels, 3 for arc segments, and 4 for groups. mi seteditmode(editmode) Sets the current editmode to: nodes - nodes segments - line segments arcsegments - arc segments blocks - block labels group - selected group This command will affect all subsequent uses of the other editing commands, if they are used WITHOUT the editaction parameter.
4.7 Zoom Commands

mi zoomnatural zooms to a natural view with sensible extents. mi zoomout zooms out by a factor of 50%. mi zoomin zoom in by a factor of 200%. mi zoom(x1,y1,x2,y2) Set the display area to be from the bottom left corner specied by (x1,y1) to the top right corner specied by (x2,y2).
4.8 View Commands

mi_showgrid Show the grid points. mi_hidegrid Hide the grid points points. mi_grid_snap(flag) Setting flag to on turns on snap to grid, setting flag to off turns off snap to grid. mi_setgrid(density,type) Change the grid spacing. The density parameter species the space between grid points, and the type parameter is set to cart for cartesian coordinates or polar for polar coordinates. mi refreshview Redraws the current view. mi minimize minimizes the active magnetics input view. mi maximize maximizes the active magnetics input view. mi restore restores the active magnetics input view from a minimized or maximized state. mi resize(width,height) resizes the active magnetics input window client area to width height.
4.9 Object Properties

mi addmaterial(matname, mu x, mu y, H c, J, Cduct, Lam d, Phi hmax, lam fill, LamType, Phi hx, Phi hy, nstr, dwire) adds a new material with called matname with the material properties: mu x Relative permeability in the x- or r-direction. mu y Relative permeability in the y- or z-direction. H c Permanent magnet coercivity in Amps/Meter. J Applied source current density in Amps/mm2. Cduct Electrical conductivity of the material in MS/m. Lam d Lamination thickness in millimeters. Phi hmax Hysteresis lag angle in degrees, used for nonlinear BH curves. Lam fill Fraction of the volume occupied per lamination that is actually lled with iron (Note that this parameter defaults to 1 in the femm preprocessor dialog box because, by default, iron completely lls the volume) Lamtype Set to 0 Not laminated or laminated in plane 1 laminated x or r 2 laminated y or z 3 magnet wire 9
4 plain stranded wire 5 Litz wire 6 square wire Phi hx Hysteresis lag in degrees in the x-direction for linear problems. Phi hy Hysteresis lag in degrees in the y-direction for linear problems. nstr Number of strands in the wire build. Should be 1 for Magnet or Square wire. dwire Diameter of each of the wires constituent strand in millimeters. Note that not all properties need be dened properties that arent dened are assigned default values. mi addbhpoint(blockname,b,h) Adds a B-H data point the the material specied by the string blockname. The point to be added has a ux density of b in units of Teslas and a eld intensity of h in units of Amps/Meter. mi clearbhpoints(blockname) Clears all B-H data points associatied with the material specied by blockname. mi addpointprop(pointpropname,a,j) adds a new point property of name pointpropname with either a specied potential a in units Webers/Meter or a point current j in units of Amps. Set the unused parameter pairs to 0. mi addboundprop(propname, A0, A1, A2, Phi, Mu, Sig, c0, c1, BdryFormat) adds a new boundary property with name propname For a Prescribed A type boundary condition, set the A0, A1, A2 and Phi parameters as required. Set all other parameters to zero. For a Small Skin Depth type boundary condtion, set the Mu to the desired relative permeability and Sig to the desired conductivity in MS/m. Set BdryFormat to 1 and all other parameters to zero. To obtain a Mixed type boundary condition, set C1 and C0 as required and BdryFormat to 2. Set all other parameters to zero. For a Strategic dual image boundary, set BdryFormat to 3 and set all other parameters to zero. For a Periodic boundary condition, set BdryFormat to 4 and set all other parameters to zero. For an Anti-Perodic boundary condition, set BdryFormat to 5 set all other parameters to zero. mi addcircprop(circuitname, i, circuittype) adds a new circuit property with name circuitname with a prescribed current. The circuittype parameter is 0 for a parallel-connected circuit and 1 for a series-connected circuit.
10
mi deletematerial(materialname) deletes the material named materialname. mi deleteboundprop(propname) deletes the boundary property named propname. mi deletecircuit(circuitname) deletes the circuit named circuitname. mi deletepointprop(pointpropname) deletes the point property named pointpropname mi_modifymaterial(BlockName,propnum,value) This function allows for modication of a materials properties without redening the entire material (e.g. so that current can be modied from run to run). The material to be modied is specied by BlockName. The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below: propnum 0 1 2 3 4 5 6 7 8 9 10 11 Symbol BlockName x y Hc J dlam hmax LamFill LamType hx hy Description Name of the material x (or r) direction relative permeability y (or z) direction relative permeability Coercivity, Amps/Meter Source current density, MA/m2 Electrical conductivity, MS/m Lamination thickness, mm Hysteresis lag angle for nonlinear problems, degrees Iron ll fraction 0 = None/In plane, 1 = parallel to x, 2=parallel to y Hysteresis lag in x-direction for linear problems, degrees Hysteresis lag in y-direction for linear problems, degrees
mi_modifyboundprop(BdryName,propnum,value) This function allows for modication of a boundary property. The BC to be modied is specied by BdryName. The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below:
11
propnum 0 1 2 3 4 5 6 7 8 9
Symbol BdryName A0 A1 A2 c0 c1 BdryFormat
Description Name of boundary property Prescribed A parameter Prescribed A parameter Prescribed A parameter Prescribed A phase Small skin depth relative permeability Small skin depth conductivity, MS/m Mixed BC parameter Mixed BC parameter Type of boundary condition: 0 = Prescribed A 1 = Small skin depth 2 = Mixed 3 = Strategic Dual Image 4 = Periodic 5 = Antiperiodic
mi_modifypointprop(PointName,propnum,value) This function allows for modication of a point property. The point property to be modied is specied by PointName. The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below: propnum 0 1 2 Symbol PointName A J Description Name of the point property Nodal potential, Weber/Meter Nodal current, Amps
mi_modifycircprop(CircName,propnum,value) This function allows for modication of a circuit property. The circuit property to be modied is specied by CircName. The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below: propnum 0 1 2 Symbol CircName i CircType Description Name of the circuit property Total current 0 = Parallel, 1 = Series
mi setcurrent(CircName,i) sets the current in the circuit specied by CircName to i.
4.10 Miscellaneous
mi savebitmap(filename) saves a bitmapped screenshot of the current view to the le specied by filename. 12
mi savemetafile(filename) saves a metale screenshot of the current view to the le specied by filename. mi refreshview Redraws the current view. mi close Closes current magnetics preprocessor document and destroys magnetics preprocessor window. mi shownames(flag) This function allow the user to display or hide the block label names on screen. To hide the block label names, flag should be 0. To display the names, the parameter should be set to 1. mi readdxf(filename) This function imports a dxf le specied by filename. mi defineouterspace(Zo,Ro,Ri) denes an axisymmetric external region to be used in conjuction with the Kelvin Transformation method of modeling unbounded problems. The Zo parameter is the z-location of the origin of the outer region, the Ro parameter is the radius of the outer region, and the Ri parameter is the radius of the inner region (i.e. the region of interest). In the exterior region, the permeability varies as a function of distance from the origin of the external region. These parameters are necessary to dene the permeability variation in the external region. mi attachouterspace marks all selected block labels as members of the external region used for modeling unbounded axisymmetric problems via the Kelvin Transformation. mi detachouterspace undenes all selected block labels as members of the external region used for modeling unbounded axisymmetric problems via the Kelvin Transformation.
5 Magnetics Post Processor Command Set

There are a number of scripting commands designed to operate in the postprocessing environment.
5.1 Data Extraction Commands

mo lineintegral(type) Calculate the line integral for the dened contour
type 0 1 2 3 4 5 name B.n H.t Contour length Stress Tensor Force Stress Tensor Torque (B.n)2 values 1 total B.n total H.t surface area DC r/x force DC torque total (B.n)2 values 2 avg B.n avg H.t DC y/z force 2 torque avg (B.n)2 values 3 2 r/x force values 4 2 y/z force -
Returns typically two values. The rst value is the result of the integral calculation, and the second value is the average of the quantity of interest over the contour. The only exception is integral 3, which evaluates Maxwells stress tensor. This integral can return up to four 13
results. For force and torque results, the 2 results are only relevant for problems where = 0. mo blockintegral(type) Calculate a block integral for the selected blocks Type 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 Denition AJ A Magnetic eld energy Hysteresis and/or lamination losses Resistive losses Block cross-section area Total losses Total current Integral of Bx (or Br ) over block Integral of By (or Bz ) over block Block volume x (or r) part of steady-state Lorentz force y (or z) part of steady-state Lorentz force x (or r) part of 2 Lorentz force y (or z) part of 2 Lorentz force Steady-state Lorentz torque 2 component of Lorentz torque Magnetic eld coenergy x (or r) part of steady-state weighted stress tensor force y (or z) part of steady-state weighted stress tensor force x (or r) part of 2 weighted stress tensor force y (or z) part of 2 weighted stress tensor force Steady-state weighted stress tensor torque 2 component of weighted stress tensor torque R2 (i.e. moment of inertia / density)
mo getpointvalues(x,y) Get the values associated with the point at (x,y). The function returns an array whose contents are, in order:
14
Symbol A B1 B2 Sig E H1 H2 Je Js Mu1 Mu2 Pe Ph ff
Denition Potential A or ux Bx if planar, Br if axisymmetric By if planar, Bz if axisymmetric conductivity stored energy density Hx if planar, Hr if axisymmetric Hy if planar, Hz if axisymmetric eddy current density source current density x if planar, r if axisymmetric y if planar, z if axisymmetric Power density dissipated through ohmic losses Power density dissipated by hysteresis Winding ll factor
The following series of functions retrieves smaller subsets of these results. mo geta(x,y) Get the potential associated with the point at (x,y). For planar problems, the reported potential is vector potential A. For axisymmetric problems, 2rA is reported. mo getb(x,y) Get the magnetic ux density associated with the point at (x,y). The return value is a list with two element representing Bx and By for planar problems and Br and Bz for axisymmetric problems. mo getconductivity(x,y) Gets the conductivity associated with the point at (x,y). mo getenergydensity(x,y) Gets the magnetic eld energy density associated with the point at (x,y). mo geth(x,y) Get the magnetic eld intensity associated with the point at (x,y). The return value is a list with two element representing Hx and Hy for planar problems and Hr and Hz for axisymmetric problems. mo getj(x,y) Get the electric current density associated with the point at (x,y). mo getmu(x,y) Get the relative magnetic permeability associated with the point at (x,y). The return value is a list with two element representing x and y for planar problems and r and z for axisymmetric problems. mo getpe(x,y) Get the ohmic loss density associated with the point at (x,y). mo getph(x,y) Get the hysteresis/laminated eddy current loss density associated with the point at (x,y). mo getfill(x,y) Get the winding factor (i.e. the average fraction of the volume lled with conductor) associated with the point at (x,y).
15
mo_makeplot(PlotType,NumPoints,Filename,FileFormat) Allows Octave access to FEMMs X-Y plot routines. If only PlotType or only PlotType and NumPoints are specied, the command is interpreted as a request to plot the requested plot type to the screen. If, in addition, the Filename parameter is specied, the plot is instead written to disk to the specied le name as an extended metale. If the FileFormat parameter is also, the command is instead interpreted as a command to write the data to disk to the speced le name, rather than display it to make a graphical plot. Valid entries for PlotType are: PlotType 0 1 2 3 4 5 6 7 8 Valid le formats are FileFormat 0 1 2 Denition Multi-column text with legend Multi-column text with no legend Mathematica-style formatting Denition Potential |B| Bn Bt |H| H n H t Jeddy Jsource + Jeddy
For example, if one wanted to plot B n to the screen with 200 points evaluated to make the graph, the command would be: mo makeplot(2,200) If this plot were to be written to disk as a metale, the command would be: mo_makeplot(2,200,c:\temp\myfile.emf) To write data instead of a plot to disk, the command would be of the form: mo_makeplot(2,200,c:\temp\myfile.txt,0) mo_getprobleminfo Returns info on problem description. Returns two values: Return value Denition 1 problem type 2 frequency in Hz mo_getcircuitproperties(circuit) Used primarily to obtain impedance information associated with circuit properties. Properties are returned for the circuit property named circuit. Six values are returned by the function. In order, these parameters are: 16
current Current carried by the circuit. volts Voltage drop across the circuit in the circuit. flux Circuits ux linkage
5.2 Selection Commands

mo seteditmode(mode) Sets the mode of the postprocessor to point, contour, or area mode. Valid entries for mode are point, contour, and area. mo selectblock(x,y) Select the block that contains point (x,y). mo groupselectblock(n) Selects all of the blocks that are labeled by block labels that are members of group n. If no number is specied (i.e. mo groupselectblock ), all blocks are selected. mo addcontour(x,y) Adds a contour point at (x,y). If this is the rst point then it starts a contour, if there are existing points the contour runs from the previous point to this point. The mo addcontour command has the same functionality as a right-button-click contour point addition when the program is running in interactive mode. mo bendcontour(angle,anglestep) Replaces the straight line formed by the last two points in the contour by an arc that spans angle degrees. The arc is actually composed of many straight lines, each of which is constrained to span no more than anglestep degrees. The angle parameter can take on values from -180 to 180 degrees. The anglestep parameter must be greater than zero. If there are less than two points dened in the contour, this command is ignored. mo selectpoint(x,y) Adds a contour point at the closest input point to (x,y). If the selected point and a previous selected points lie at the ends of an arcsegment, a contour is added that traces along the arcsegment. The mo selectpoint command has the same functionality as the left-button-click contour point selection when the program is running in interactive mode. mo clearcontour Clear a prevously dened contour mo clearblock Clear block selection
5.3 Zoom Commands

mo_zoomnatural Zoom to the natural boundaries of the geometry. mo_zoomin Zoom in one level. mo_zoomout Zoom out one level. mo zoom(x1,y1,x2,y2) Zoom to the window dened by lower left corner (x1,y1) and upper right corner (x2,y2). 17
5.4 View Commands

mo_showmesh Show the mesh. mo_hidemesh Hide the mesh. mo_showpoints Show the node points from the input geometry. mo_hidepoints Hide the node points from the input geometry. mo smooth(flag) This function controls whether or not smoothing is applied to the B and H elds, which are naturally piece-wise constant over each element. Setting flag equal to on turns on smoothing, and setting flag to off turns off smoothing. mo_showgrid Show the grid points. mo_hidegrid Hide the grid points points. mo_grid_snap(flag) Setting flag to on turns on snap to grid, setting flag to off turns off snap to grid. mo_setgrid(density,type) Change the grid spacing. The density parameter species the space between grid points, and the type parameter is set to cart for cartesian coordinates or polar for polar coordinates. mo_hidedensityplot hides the ux density plot. mo_showdensityplot(legend,gscale,upper_B,lower_B,type) Shows the ux density plot with options: legend Set to 0 to hide the plot legend or 1 to show the plot legend. gscale Set to 0 for a colour density plot or 1 for a grey scale density plot. upper_B Sets the upper display limit for the density plot. lower_B Sets the lower display limit for the density plot. type Type of density plot to display. Valid entries are mag, real, and imag for magnitude, real component, and imaginary component of B, respectively. Alternatively, current density can be displayed by specifying jmag, jreal, and jimag for magnitude, real component, and imaginary component of J, respectively. mo_hidecontourplot Hides the contour plot. mo_showcontourplot(numcontours,lower_A,upper_A,type) shows the A contour plot with options: numcontours Number of A equipotential lines to be plotted. upper_A Upper limit for A contours. lower_A Lower limit for A contours.
18
type Choice of real, imag, or both to show either the real, imaginary of both real and imaginary components of A. mo minimize minimizes the active magnetics output view. mo maximize maximizes the active magnetics output view. mo restore restores the active magnetics output view from a minimized or maximized state. mo resize(width,height) resizes the active magnetics output window client area to width height.
5.5 Miscellaneous
mo close Closes the current post-processor instance. mo refreshview Redraws the current view. mo reload Reloads the solution from disk. mo savebitmap(filename) saves a bitmapped screen shot of the current view to the le specied by filename. mo savemetafile(filename) saves a metale screenshot of the current view to the le specied by filename. mo shownames(flag) This function allow the user to display or hide the block label names on screen. To hide the block label names, flag should be 0. To display the names, the parameter should be set to 1.
6 Electrostatics Preprocessor Command Set

A number of different commands are available in the preprocessor. Two naming conventions can be used: one which separates words in the command names by underscores, and one that eliminates the underscores.

ei addnode(x,y) Add a new node at x,y ei addsegment(x1,y1,x2,y2) Add a new line segment from node closest to (x1,y1) to node closest to (x2,y2) ei addblocklabel(x,y) Add a new block label at (x,y) ei addarc(x1,y1,x2,y2,angle,maxseg) Add a new arc segment from the nearest node to (x1,y1) to the nearest node to (x2,y2) with angle angle divided into maxseg segments.
19
ei drawline(x1,y1,x2,y2) Adds nodes at (x1,y1) and (x2,y2) and adds a line between the nodes. ei drawpolyline([x1,y1;x2,y2...]) Adds nodes at each of the specied points and connects them with segments. ei drawpolygon([x1,y1;x2,y2...]) Adds nodes at each of the specied points and connects them with segments to form a closed contour. ei drawarc(x1,y1,x2,y2,angle,maxseg) Adds nodes at (x1,y1) and (x2,y2) and adds an arc of the specied angle and discretization connecting the nodes. ei drawrectangle(x1,y1,x2,y2) Adds nodes at the corners of a rectangle dened by the points (x1,y1) and (x2,y2), then adds segments connecting the corners of the rectangle. ei deleteselected Delete all selected objects. ei deleteselectednodes Delete selected nodes. ei deleteselectedlabels Delete selected block labels. ei deleteselectedsegments Delete selected segments. ei deleteselectedarcsegments Delete selects arcs.

ei clearselected Clear all selected nodes, blocks, segments and arc segments. ei selectsegment(x,y) Select the line segment closest to (x,y) ei selectnode(x,y) Select the node closest to (x,y). Returns the coordinates of the selected node. ei selectlabel(x,y) Select the label closet to (x,y). Returns the coordinates of the selected label. ei selectarcsegment(x,y) Select the arc segment closest to (x,y) ei selectgroup(n) Select the nth group of nodes, segments, arc segments and block labels. This function will clear all previously selected elements and leave the edit mode in 4 (group)

ei setnodeprop(propname,groupno, inconductor) Set the selected nodes to have the nodal property propname and group number groupno. The inconductor string species which conductor the node belongs to. If the node doesnt belong to a named conductor, this parameter can be set to <None>.
20
ei setblockprop(blockname, automesh, meshsize, group) Set the selected block labels to have the properties: Block property blockname. automesh: 0 = mesher defers to mesh size constraint dened in meshsize, 1 = mesher automatically chooses the mesh density. meshsize: size constraint on the mesh in the block marked by this label. A member of group number group ei setsegmentprop(name, elmsize, automesh, hide, group, inconductor) Set the select segments to have: Boundary property name Local element size along segment no greater than elmsize automesh: 0 = mesher defers to the element constraint dened by elementsize, 1 = mesher automatically chooses mesh size along the selected segments hide: 0 = not hidden in post-processor, 1 == hidden in post processor A member of group number group A member of the conductor specied by the string inconductor. If the segment is not part of a conductor, this parameter can be specied as <None>. ei setarcsegmentprop(maxsegdeg, propname, hide, group, inconductor) Set the selected arc segments to: Meshed with elements that span at most maxsegdeg degrees per element Boundary property propname hide: 0 = not hidden in post-processor, 1 == hidden in post processor A member of group number group A member of the conductor specied by the string inconductor. If the segment is not part of a conductor, this parameter can be specied as <None>.

ei probdef(units,type,precision,depth,minangle) changes the problem denition. The units parameter species the units used for measuring length in the problem domain. Valid units entries are inches, millimeters, centimeters, mils, meters, and micrometers. Set problemtype to planar for a 2-D planar problem, or to axi for an axisymmetric problem. The precision parameter dictates the precision required by the solver. For example, entering 1.E-8 requires the RMS of the residual to be less than 108 . The fourth parameter represents the depth of the problem in the into-the-page direction for 2-D planar problemsjust enter 0 as a placeholder for axisymmetric problems. The sixth parameter represents the minimum angle constraint sent to the mesh generator (usually 30 degrees). 21
ei analyze(flag) runs the electrostatics solver. The flag parameter controls whether the solver window is visible or minimized. For a visible window, either specify no value for flag or specify 0. For a minimized window, flag should be set to 1. ei loadsolution loads and displays the solution corresponding to the current geometry. ei setfocus(documentname) Switches the electrostatics input le upon which commands are to act. If more than one electrostatics input le is being edited at a time, this command can be used to switch between les so that the mutiple les can be operated upon programmatically. documentname should contain the name of the desired document as it appears on the windows title bar. ei saveas(filename) saves the le with name filename.
6.5 Mesh Commands

ei createmesh runs triangle to create a mesh. Note that this is not a necessary precursor of performing an analysis, as ei analyze will make sure the mesh is up to date before running an analysis. The number of elements in the mesh is pushed back onto the lua stack. ei showmesh toggles the ag that shows or hides the mesh. ei purgemesh clears the mesh out of both the screen and memory.

ei copyrotate(bx, by, angle, copies ) bx, by base point for rotation angle angle by which the selected objects are incrementally shifted to make each copy. angle is measured in degrees. copies number of copies to be produced from the selected objects. ei copyrotate2(bx, by, angle, copies, editaction ) bx, by base point for rotation angle angle by which the selected objects are incrementally shifted to make each copy. angle is measured in degrees. copies number of copies to be produced from the selected objects. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group ei copytranslate(dx, dy, copies) dx,dy distance by which the selected objects are incrementally shifted. copies number of copies to be produced from the selected objects.
22
ei copytranslate2(dx, dy, copies, editaction) dx,dy distance by which the selected objects are incrementally shifted. copies number of copies to be produced from the selected objects. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group ei createradius(x,y,r)turnsacornerlocatedat(x,y)intoacurveofradiusr. ei moverotate(bx,by,shiftangle) bx, by base point for rotation shiftangle angle in degrees by which the selected objects are rotated. ei moverotate2(bx,by,shiftangle, editaction) bx, by base point for rotation shiftangle angle in degrees by which the selected objects are rotated. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group ei movetranslate(dx,dy) dx,dy distance by which the selected objects are shifted. ei movetranslate2(dx,dy,editaction) dx,dy distance by which the selected objects are shifted. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group ei scale(bx,by,scalefactor) bx, by base point for scaling scalefactor a multiplier that determines how much the selected objects are scaled ei scale2(bx,by,scalefactor,editaction) bx, by base point for scaling scalefactor a multiplier that determines how much the selected objects are scaled editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group ei mirror(x1,y1,x2,y2) mirror the selected objects about a line passing through the points (x1,y1) and (x2,y2). ei mirror2(x1,y1,x2,y2,editaction) mirror the selected objects about a line passing through the points (x1,y1) and (x2,y2). Valid editaction entries are 0 for nodes, 1 for lines (segments), 2 for block labels, 3 for arc segments, and 4 for groups. ei seteditmode(editmode) Sets the current editmode to: 23
nodes - nodes segments - line segments arcsegments - arc segments blocks - block labels group - selected group This command will affect all subsequent uses of the other editing commands, if they are used WITHOUT the editaction parameter.
6.7 Zoom Commands

ei zoomnatural zooms to a natural view with sensible extents. ei zoomout zooms out by a factor of 50%. ei zoomin zoom in by a factor of 200%. ei zoom(x1,y1,x2,y2) Set the display area to be from the bottom left corner specied by (x1,y1) to the top right corner specied by (x2,y2).
6.8 View Commands

ei showgrid Show the grid points. ei hidegrid Hide the grid points points. ei gridsnap(flag) Setting ag to on turns on snap to grid, setting ag to offturns off snap to grid. ei setgrid(density,type) Change the grid spacing. The density parameter species the space between grid points, and the type parameter is set to cart for Cartesian coordinates or polar for polar coordinates. ei refreshview Redraws the current view. ei minimize minimizes the active electrostatics input view. ei maximize maximizes the active electrostatics input view. ei restore restores the active electrostatics input view from a minimized or maximized state. ei resize(width,height) resizes the active electrostatics input window client area to width height.
24

ei addmaterial(matname, ex, ey, qv) adds a new material with called matname with the material properties: ex Relative permittivity in the x- or r-direction. ey Relative permittivity in the y- or z-direction. qv Volume charge density in units of C/m3 ei addpointprop(pointname,Vp,qp) adds a new point property of name pointname with either a specied potential Vp a point charge density qp in units of C/m. . ei addboundprop(boundname, Vs, qs, c0, c1, BdryFormat) adds a new boundary property with name boundname For a Fixed Voltage type boundary condition, set the Vs parameter to the desired voltage and all other parameters to zero. To obtain a Mixed type boundary condition, set C1 and C0 as required and BdryFormat to 1. Set all other parameters to zero. To obtain a prescribes surface charge density, set qs to the desired charge density in C/m2 and set BdryFormat to 2. For a Periodic boundary condition, set BdryFormat to 3 and set all other parameters to zero. For an Anti-Perodic boundary condition, set BdryFormat to 4 set all other parameters to zero. ei addconductorprop(conductorname, Vc, qc, conductortype) adds a new conductor property with name conductorname with either a prescribed voltage or a prescribed total charge. Set the unused property to zero. The conductortype parameter is 0 for prescribed charge and 1 for prescribed voltage. ei deletematerial(materialname) deletes the material named materialname. ei deleteboundprop(boundname) deletes the boundary property named boundname. ei deleteconductor(conductorname) deletes the conductor named conductorname. ei deletepointprop(pointname) deletes the point property named pointname ei modifymaterial(BlockName,propnum,value) This function allows for modication of a materials properties without redening the entire material (e.g. so that current can be modied from run to run). The material to be modied is specied by BlockName. The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below:
25
propnum 0 1 2 3
Symbol BlockName ex ey qs
Description Name of the material x (or r) direction relative permittivity y (or z) direction relative permittivity Volume charge
ei modifyboundprop(BdryName,propnum,value) This function allows for modication of a boundary property. The BC to be modied is specied by BdryName. The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below: propnum 0 1 2 3 4 5 Symbol BdryName Vs qs c0 c1 BdryFormat Description Name of boundary property Fixed Voltage Prescribed charge density Mixed BC parameter Mixed BC parameter Type of boundary condition: 0 = Prescribed V 1 = Mixed 2 = Surface charge density 3 = Periodic 4 = Antiperiodic
ei modifypointprop(PointName,propnum,value) This function allows for modication of a point property. The point property to be modied is specied by PointName. The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below: propnum Symbol Description 0 PointName Name of the point property Vp Prescribed nodal voltage 1 2 qp Point charge density in C/m ei modifyconductorprop(ConductorName,propnum,value) This function allows for modication of a conductor property. The conductor property to be modied is specied by ConductorName. The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below: propnum Symbol Description 0 ConductorName Name of the conductor property 1 Vc Conductor voltage 2 qc Total conductor charge 3 ConductorType 0 = Prescribed charge, 1 = Prescribed voltage
26
6.10 Miscellaneous
ei savebitmap(filename) saves a bitmapped screenshot of the current view to the le specied by filename. ei savemetafile(filename) saves a metale screenshot of the current view to the le specied by filename. ei refreshview Redraws the current view. ei close closes the preprocessor window and destroys the current document. ei shownames(flag) This function allow the user to display or hide the block label names on screen. To hide the block label names, flag should be 0. To display the names, the parameter should be set to 1. ei readdxf(filename) This function imports a dxf le specied by filename. ei defineouterspace(Zo,Ro,Ri) denes an axisymmetric external region to be used in conjuction with the Kelvin Transformation method of modeling unbounded problems. The Zo parameter is the z-location of the origin of the outer region, the Ro parameter is the radius of the outer region, and the Ri parameter is the radius of the inner region (i.e. the region of interest). In the exterior region, the permeability varies as a function of distance from the origin of the external region. These parameters are necessary to dene the permeability variation in the external region. ei attachouterspace marks all selected block labels as members of the external region used for modeling unbounded axisymmetric problems via the Kelvin Transformation. ei detachouterspace undenes all selected block labels as members of the external region used for modeling unbounded axisymmetric problems via the Kelvin Transformation.
7 Electrostatics Post Processor Command Set

There are a number of scripting commands designed to operate in the postprocessor. As with the preprocessor commands, these commands can be used with either the underscore naming or with the no-underscore naming convention.

eo lineintegral(type) Calculates the line integral for the dened contour type 0 1 2 3 4 Integral E t Dn Contour length Force from stress tensor Torque from stress tensor 27
This integral a single result for eld intensity, ux density, contour length, and torque. For force, an array with two elements is returned representing force in the x and y or r and z directions. eo blockintegral(type) Calculates a block integral for the selected blocks type 0 1 2 3 4 5 6 Integral Stored Energy Block Cross-section Block Volume Average D over the block Average E over the block Weighted Stress Tensor Force Weighted Stress Tensor Torque
This integral a single result for energy, volume and torque. For force, D and E, the function returns an array with two elements is returned representing components in the x and y or r and z directions. eo getpointvalues(x,y) Gets the values associated with the point at (x, y). The function returns an array of results whose elements represent: Symbol V Dx Dy Ex Ey ex ey nrg Denition Voltage x- or r- direction component of displacement y- or z- direction component of displacement x- or r- direction component of electric eld intensity y- or z- direction component of electric eld intensity x- or r- direction component of permittivity y- or z- direction component of permittivity electric eld energy density
eo getv(x,y) Gets the voltage associated with the point at (x,y). eo getd(x,y) Gets the electric ux density associated with the point at (x,y). The return value is a list with two element representing Dx and Dy for planar problems and Dr and Dz for axisymmetric problems. eo gete(x,y) Gets the electric eld intensity associated with the point at (x,y). The return value is a list with two element representing Ex and Ey for planar problems and Er and Ez for axisymmetric problems. eo getenergydensity(x,y) Gets the electric eld energy density associated with the point at (x,y). eo getperm(x,y) Gets the electric permittivity associated with the point at (x,y). The return value is a list with two element representing x and y for planar problems and r and z for axisymmetric problems.
28
eo makeplot(PlotType,NumPoints,Filename,FileFormat) Allows Octave to access to the X-Y plot routines. If only PlotType or only PlotType and NumPoints are specied, the command is interpreted as a request to plot the requested plot type to the screen. If, in addition, the Filename parameter is specied, the plot is instead written to disk to the specied le name as an extended metale. If the FileFormat parameter is also, the command is instead interpreted as a command to write the data to disk to the speced le name, rather than display it to make a graphical plot. Valid entries for PlotType are: PlotType 0 1 2 3 4 5 6 Denition V (Voltage) |D| (Magnitude of flux density) D . n (Normal flux density) D . t (Tangential flux density) |E| (Magnitude of field intensity) E . n (Normal field intensity) E . t (Tangential field intensity)
Valid le formats are: FileFormat Denition 0 Multi-column text with legend 1 Multi-column text with no legend 2 Mathematica-style formatting For example, if one wanted to plot V to the screen with 200 points evaluated to make the graph, the command would be: eo makeplot(0,200) If this plot were to be written to disk as a metale, the command would be: eo makeplot(0,200,c:temp.emf) To write data instead of a plot to disk, the command would be of the form: eo makeplot(0,200,c:temp.txt,0) eo getprobleminfo Returns info on problem description. Returns two values: the Problem type (0 for planar and 1 for axisymmetric) and the depth assumed for planar problems in units of meters. eo getconductorproperties(conductor) Properties are returned for the conductor property named conductor. Two values are returned: The voltage of the specied conductor, and the charge carried on the specied conductor.

eo seteditmode(mode) Sets the mode of the postprocessor to point, contour, or area mode. Valid entries for mode are point, contour, and area. eo selectblock(x,y) Select the block that contains point (x,y).
29
eo groupselectblock(n) Selects all of the blocks that are labeled by block labels that are members of group n. If no number is specied (i.e. eo groupselectblock), all blocks are selected. eo selectconductor(name) Selects all nodes, segments, and arc segments that are part of the conductor specied by the string (name). This command is used to select conductors for the purposes of the weighted stress tensor force and torque integrals, where the conductors are points or surfaces, rather than regions (i.e. cant be selected with eo selectblock). eo addcontour(x,y) Adds a contour point at (x,y). If this is the rst point then it starts a contour, if there are existing points the contour runs from the previous point to this point. The eo addcontour command has the same functionality as a right-button-click contour point addition when the program is running in interactive mode. eo bendcontour(angle,anglestep) Replaces the straight line formed by the last two points in the contour by an arc that spans angle degrees. The arc is actually composed of many straight lines, each of which is constrained to span no more than anglestep degrees. The angle parameter can take on values from -180 to 180 degrees. The anglestep parameter must be greater than zero. If there are less than two points dened in the contour, this command is ignored. eo selectpoint(x,y) Adds a contour point at the closest input point to (x,y). If the selected point and a previous selected points lie at the ends of an arcsegment, a contour is added that traces along the arcsegment. The selectpoint command has the same functionality as the left-button-click contour point selection when the program is running in interactive mode. eo clearcontour Clear a prevously dened contour eo clearblock Clear block selection
7.3 Zoom Commands

eo zoomnatural Zoom to the natural boundaries of the geometry. eo zoomin Zoom in one level. eo zoomout Zoom out one level. eo zoom(x1,y1,x2,y2) Zoom to the window dened by lower left corner (x1,y1) and upper right corner (x2,y2).
7.4 View Commands

eo showmesh Show the mesh. eo hidemesh Hide the mesh. eo showpoints Show the node points from the input geometry. 30
eo hidepoints Hide the node points from the input geometry. eo smooth(flag) This function controls whether or not smoothing is applied to the D and E elds which are naturally piece-wise constant over each element. Setting ag equal to on turns on smoothing, and setting ag to off turns off smoothing. eo showgrid Show the grid points. eo hidegrid Hide the grid points points. eo gridsnap(flag) Setting ag to on turns on snap to grid, setting ag to off turns off snap to grid. eo setgrid(density,type) Change the grid spacing. The density parameter species the space between grid points, and the type parameter is set to cart for Cartesian coordinates or polar for polar coordinates. eo hidedensityplot hides the ux density plot. eo showdensityplot(legend,gscale,upper D,lower D,type) Shows the ux density plot with options: legend Set to 0 to hide the plot legend or 1 to show the plot legend. gscale Set to 0 for a colour density plot or 1 for a grey scale density plot. upper D Sets the upper display limit for the density plot. lower D Sets the lower display limit for the density plot. type Sets the type of density plot. A value of 0 plots voltage, 1 plots the magnitude of D, and 2 plots the magnitude of E eo hidecontourplot Hides the contour plot. eo showcontourplot(numcontours,lower V,upper V) shows the V contour plot with options: numcontours Number of equipotential lines to be plotted. upper V Upper limit for contours. lower V Lower limit for contours. eo showvectorplot(type,scalefactor) controls the display of vectors denoting the eld strength and direction. The parameters taken are the type of plot, which should be set to 0 for no vector plot, 1 for ux density D, and 2 for eld intensity E. The scalefactor determines the relative length of the vectors. If the scale is set to 1, the length of the vectors are chosen so that the highest ux density corresponds to a vector that is the same length as the current grid size setting. eo minimize minimizes the active electrostatics output view. eo maximize maximizes the active electrostatics output view. 31
eo restore restores the active electrostatics output view from a minimized or maximized state. eo resize(width,height) resizes the active electrostatics output window client area to width height.
7.5 Miscellaneous
eo close close the current postprocessor window. eo refreshview Redraws the current view. eo reload Reloads the solution from disk. eo savebitmap(filename) saves a bitmapped screen shot of the current view to the le specied by filename. eo savemetafile(filename) saves a metale screenshot of the current view to the le specied by filename. eo shownames(flag) This function allow the user to display or hide the block label names on screen. To hide the block label names, flag should be 0. To display the names, the parameter should be set to 1.
8 Heat Flow Preprocessor Lua Command Set


hi addnode(x,y) Add a new node at x,y hi addsegment(x1,y1,x2,y2) Add a new line segment from node closest to (x1,y1) to node closest to (x2,y2) hi addblocklabel(x,y) Add a new block label at (x,y) hi addarc(x1,y1,x2,y2,angle,maxseg) Add a new arc segment from the nearest node to (x1,y1) to the nearest node to (x2,y2) with angle angle divided into maxseg segments. hi deleteselected Delete all selected objects. hi deleteselectednodes Delete selected nodes. hi deleteselectedlabels Delete selected block labels.
32
hi deleteselectedsegments Delete selected segments. hi deleteselectedarcsegments Delete selects arcs.

hi clearselected() Clear all selected nodes, blocks, segments and arc segments. hi selectsegment(x,y) Select the line segment closest to (x,y) hi selectnode(x,y) Select the node closest to (x,y). Returns the coordinates of the selected node. hi selectlabel(x,y) Select the label closet to (x,y). Returns the coordinates of the selected label. hi selectarcsegment(x,y) Select the arc segment closest to (x,y) hi selectgroup(n) Select the nth group of nodes, segments, arc segments and block labels. This function will clear all previously selected elements and leave the edit mode in 4 (group)

hi setnodeprop("propname",groupno, "inconductor") Set the selected nodes to have the nodal property "propname" and group number groupno. The "inconductor" string species which conductor the node belongs to. If the node doesnt belong to a named conductor, this parameter can be set to "<None>". hi setblockprop("blockname", automesh, meshsize, group) Set the selected block labels to have the properties: Block property "blockname". automesh: 0 = mesher defers to mesh size constraint dened in meshsize, 1 = mesher automatically chooses the mesh density. meshsize: size constraint on the mesh in the block marked by this label. A member of group number group hi setsegmentprop("propname", elementsize, automesh, hide, group, "inconductor") Set the select segments to have: Boundary property "propname" Local element size along segment no greater than elementsize automesh: 0 = mesher defers to the element constraint dened by elementsize, 1 = mesher automatically chooses mesh size along the selected segments hide: 0 = not hidden in post-processor, 1 == hidden in post processor A member of group number group 33
A member of the conductor specied by the string "inconductor". If the segment is not part of a conductor, this parameter can be specied as "<None>". hi setarcsegmentprop(maxsegdeg, "propname", hide, group, "inconductor") Set the selected arc segments to: Meshed with elements that span at most maxsegdeg degrees per element Boundary property "propname" hide: 0 = not hidden in post-processor, 1 == hidden in post processor A member of group number group A member of the conductor specied by the string "inconductor". If the segment is not part of a conductor, this parameter can be specied as "<None>".

hi probdef(units,type,precision,(depth),(minangle)) changes the problem denition. The units parameter species the units used for measuring length in the problem domain. Valid "units" entries are "inches", "millimeters", "centimeters", "mils", "meters, and "micrometers". Set problemtype to "planar" for a 2-D planar problem, or to "axi" for an axisymmetric problem. The precision parameter dictates the precision required by the solver. For example, entering 1.E-8 requires the RMS of the residual to be less than 108 . A fourth parameter, representing the depth of the problem in the into-thepage direction for 2-D planar problems, can also be specied for planar problems. A sixth parameter represents the minimum angle constraint sent to the mesh generator. hi analyze(flag) runs hsolv to solve the problem. The flag parameter controls whether the hsolve window is visible or minimized. For a visible window, either specify no value for flag or specify 0. For a minimized window, flag should be set to 1. hi loadsolution() loads and displays the solution corresponding to the current geometry. hi setfocus(documentname) Switches the heat ow input le upon which Lua commands are to act. If more than one heat ow input le is being edited at a time, this command can be used to switch between les so that the mutiple les can be operated upon programmatically via Lua. documentname should contain the name of the desired document as it appears on the windows title bar. hi saveas("filename") saves the le with name "filename". Note if you use a path you must use two backslashes e.g. c:\\temp\\myfile.feh
8.5 Mesh Commands

hi createmesh() runs triangle to create a mesh. Note that this is not a necessary precursor of performing an analysis, as hi analyze() will make sure the mesh is up to date before running an analysis. The number of elements in the mesh is pushed back onto the lua stack.
34
hi showmesh() toggles the ag that shows or hides the mesh. hi purgemesh() clears the mesh out of both the screen and memory.

hi copyrotate(bx, by, angle, copies, (editaction) ) bx, by base point for rotation angle angle by which the selected objects are incrementally shifted to make each copy. angle is measured in degrees. copies number of copies to be produced from the selected objects. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group hi copytranslate(dx, dy, copies, (editaction)) dx,dy distance by which the selected objects are incrementally shifted. copies number of copies to be produced from the selected objects. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group mi createradius(x,y,r)turnsacornerlocatedat(x,y)intoacurveofradiusr. hi moverotate(bx,by,shiftangle (editaction)) bx, by base point for rotation shiftangle angle in degrees by which the selected objects are rotated. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group hi movetranslate(dx,dy,(editaction)) dx,dy distance by which the selected objects are shifted. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group hi scale(bx,by,scalefactor,(editaction)) bx, by base point for scaling scalefactor a multiplier that determines how much the selected objects are scaled editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group hi mirror(x1,y1,x2,y2,(editaction)) mirror the selected objects about a line passing through the points (x1,y1) and (x2,y2). Valid editaction entries are 0 for nodes, 1 for lines (segments), 2 for block labels, 3 for arc segments, and 4 for groups. hi seteditmode(editmode) Sets the current editmode to: "nodes" nodes "segments" - line segments 35
"arcsegments" - arc segments "blocks" - block labels "group" - selected group This command will affect all subsequent uses of the other editing commands, if they are used WITHOUT the editaction parameter.
8.7 Zoom Commands

hi zoomnatural() zooms to a natural view with sensible extents. hi zoomout() zooms out by a factor of 50%. hi zoomin() zoom in by a factor of 200%. hi zoom(x1,y1,x2,y2) Set the display area to be from the bottom left corner specied by (x1,y1) to the top right corner specied by (x2,y2).
8.8 View Commands

hi showgrid() Show the grid points. hi hidegrid() Hide the grid points points. hi gridsnap("flag") Setting ag to on turns on snap to grid, setting ag to offturns off snap to grid. hi setgrid(density,"type") Change the grid spacing. The density parameter species the space between grid points, and the type parameter is set to "cart" for Cartesian coordinates or "polar" for polar coordinates. hi refreshview() Redraws the current view. hi minimize() minimizes the active heat ow input view. hi maximize() maximizes the active heat ow input view. hi restore() restores the active heat ow input view from a minimized or maximized state. hi resize(width,height) resizes the active heat ow input window client area to width height.
36

hi addmaterial("materialname", kx, ky, qv) adds a new material with called "materialname" with the material properties: kx Thermal conductivity in the x- or r-direction. ky Thermal conductivity in the y- or z-direction. qv Volume heat generation density in units of W/m3 hi addpointprop("pointpropname",Tp,qp) adds a new point property of name "pointpropname" with either a specied temperature Tp or a point heat generation density qp in units of W/m. . hi addboundprop("boundpropname", BdryFormat, Tset, qs, Tinf, h, beta) adds a new boundary property with name "boundpropname". For a Fixed Temperature type boundary condition, set the Tset parameter to the desired temperature and all other parameters to zero. To obtain a Heat Flux type boundary condition, set qs to be the heat ux density and BdryFormat to 1. Set all other parameters to zero. To obtain a convection boundary condition, set h to the desired heat transfer coefcient and Tinf to the desired external temperature and set BdryFormat to 2. For a Radiation boundary condition, set beta equal to the desired emissivity and Tinf to the desired external temperature and set BdryFormat to 3. For a Periodic boundary condition, set BdryFormat to 4 and set all other parameters to zero. For an Anti-Perodic boundary condition, set BdryFormat to 5 set all other parameters to zero. hi addconductorprop("conductorname", Tc, qc, conductortype) adds a new conductor property with name "conductorname" with either a prescribed temperature or a prescribed total heat ux. Set the unused property to zero. The conductortype parameter is 0 for prescribed heat ux and 1 for prescribed temperature. hi deletematerial("materialname") deletes the material named "materialname". hi deleteboundprop("boundpropname") deletes the boundary property named "boundpropname". hi deleteconductor("conductorname") deletes the conductor named conductorname. hi deletepointprop("pointpropname") deletes the point property named "pointpropname" hi modifymaterial("BlockName",propnum,value) This function allows for modication of a materials properties without redening the entire material (e.g. so that current can be modied from run to run). The material to be modied is specied by "BlockName". The next parameter is the number of the property to be set. The last number is the value to 37
be applied below: propnum 0 1 2 3
to the specied property. The various properties that can be modied are listed Symbol BlockName kx ky qs Description Name of the material x (or r) direction thermal conductivity y (or z) direction thermal conductivity Volume heat generation
hi modifyboundprop("BdryName",propnum,value) This function allows for modication of a boundary property. The BC to be modied is specied by "BdryName". The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below: propnum Symbol Description 0 BdryName Name of boundary property 1 BdryFormat Type of boundary condition: 0 = Prescribed temperature 1 = Heat Flux 2 = Convection 3 = Radiation 4 = Periodic 5 = Antiperiodic 2 Tset Fixed Temperature 3 qs Prescribed heat ux density Tinf External temperature 4 5 h Heat transfer coefcient 6 beta Emissivity hi modifypointprop("PointName",propnum,value) This function allows for modication of a point property. The point property to be modied is specied by "PointName". The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below: propnum Symbol Description 0 PointName Name of the point property 1 Tp Prescribed nodal temperature 2 qp Point heat generation in W/m hi modifyconductorprop("ConductorName",propnum,value) This function allows for modication of a conductor property. The conductor property to be modied is specied by "ConductorName". The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below:
38
propnum 0 1 2 3
Symbol ConductorName Tc qc ConductorType
Description Name of the conductor property Conductor temperature Total conductor heat ux 0 = Prescribed heat ow, 1 = Prescribed temperature
hi addtkpoint("materialname",T,k) adds the point (T,k) to the thermal conductivity vs. temperature curve for the material specied by "materialname". hi cleartkpoints("materialname") erases all of the thermal conductivity points that have been dened for the material named "materialname".
8.10 Miscellaneous
hi savebitmap("filename") saves a bitmapped screenshot of the current view to the le specied by "filename", subject to the printf-type formatting explained previously for the hi saveas command. hi savemetafile("filename") saves a metale screenshot of the current view to the le specied by "filename", subject to the printf-type formatting explained previously for the hi saveas command. hi refreshview() Redraws the current view. hi close() closes the preprocessor window and destroys the current document. hi shownames(flag) This function allow the user to display or hide the block label names on screen. To hide the block label names, flag should be 0. To display the names, the parameter should be set to 1. hi readdxf("filename") This function imports a dxf le specied by "filename". hi defineouterspace(Zo,Ro,Ri) denes an axisymmetric external region to be used in conjuction with the Kelvin Transformation method of modeling unbounded problems. The Zo parameter is the z-location of the origin of the outer region, the Ro parameter is the radius of the outer region, and the Ri parameter is the radius of the inner region (i.e. the region of interest). In the exterior region, the permeability varies as a function of distance from the origin of the external region. These parameters are necessary to dene the permeability variation in the external region. hi attachouterspace() marks all selected block labels as members of the external region used for modeling unbounded axisymmetric problems via the Kelvin Transformation. hi detachouterspace() undenes all selected block labels as members of the external region used for modeling unbounded axisymmetric problems via the Kelvin Transformation.
39
9 Heat Flow Post Processor Command Set

There are a number of Lua scripting commands designed to operate in the postprocessor. As with the preprocessor commands, these commands can be used with either the underscore naming or with the no-underscore naming convention.

ho lineintegral(type) Calculate the line integral for the dened contour type 0 1 2 3 Integral Temperature difference (G t) Heat ux through the contour (F n) Contour length Average Temperature
This integral returns either 1 or 2 values, depending on the integral type, e.g. : Ftot, Favg = ho lineintegral(2) ho blockintegral(type) Calculate a block integral for the selected blocks type 0 1 2 3 4 Integral Average T over the block Block Cross-section Block Volume Average F over the block Average G over the block
Returns one or two oating point values as results, e.g.: Gx, Gy = ho blockintegral(4) ho getpointvalues(X,Y) Get the values associated with the point at x,y The return values, in order, are: Symbol Denition V Temperature Fx x- or r- direction component of heat ux density Fy y- or z- direction component of heat ux density Gx x- or r- direction component of temperature gradient Gy y- or z- direction component of temperature gradient kx x- or r- direction component of thermal conductivity ky y- or z- direction component of thermal conductivity Example: To catch all values at (0.01,0) use T,Fx,Fy,Gx,Gy,kx,ky= ho getpointvalues(0.01,0) ho makeplot(PlotType,NumPoints,Filename,FileFormat) Allows Lua access to the X-Y plot routines. If only PlotType or only PlotType and NumPoints are specied, the command is interpreted as a request to plot the requested plot type to the screen. If, in addition, 40
the Filename parameter is specied, the plot is instead written to disk to the specied le name as an extended metale. If the FileFormat parameter is also, the command is instead interpreted as a command to write the data to disk to the speced le name, rather than display it to make a graphical plot. Valid entries for PlotType are: PlotType 0 1 2 3 4 5 6 Denition V (Temperature) |D| (Magnitude of heat flux density) D . n (Normal heat flux density) D . t (Tangential heat flux density) |E| (Magnitude of field intensity) E . n (Normal field intensity) E . t (Tangential field intensity)
Valid le formats are: FileFormat Denition 0 Multi-column text with legend 1 Multi-column text with no legend 2 Mathematica-style formatting For example, if one wanted to plot V to the screen with 200 points evaluated to make the graph, the command would be: ho makeplot(0,200) If this plot were to be written to disk as a metale, the command would be: ho makeplot(0,200,"c:temp.emf") To write data instead of a plot to disk, the command would be of the form: ho makeplot(0,200,"c:temp.txt",0) ho getprobleminfo() Returns info on problem description. Returns two values: the Problem type (0 for planar and 1 for axisymmetric) and the depth assumed for planar problems in units of meters. ho getconductorproperties("conductor") Properties are returned for the conductor property named conductor. Two values are returned: The temperature of the specied conductor, and the total heat ux through the specied conductor.

ho seteditmode(mode) Sets the mode of the postprocessor to point, contour, or area mode. Valid entries for mode are "point", "contour", and "area". ho selectblock(x,y) Select the block that contains point (x,y). ho groupselectblock(n) Selects all of the blocks that are labeled by block labels that are members of group n. If no number is specied (i.e. ho groupselectblock()), all blocks are selected. 41
ho selectconductor("name") Selects all nodes, segments, and arc segments that are part of the conductor specied by the string ("name"). This command is used to select conductors for the purposes of the weighted stress tensor force and torque integrals, where the conductors are points or surfaces, rather than regions (i.e. cant be selected with ho selectblock). ho addcontour(x,y) Adds a contour point at (x,y). If this is the rst point then it starts a contour, if there are existing points the contour runs from the previous point to this point. The ho addcontour command has the same functionality as a right-button-click contour point addition when the program is running in interactive mode. ho bendcontour(angle,anglestep) Replaces the straight line formed by the last two points in the contour by an arc that spans angle degrees. The arc is actually composed of many straight lines, each of which is constrained to span no more than anglestep degrees. The angle parameter can take on values from -180 to 180 degrees. The anglestep parameter must be greater than zero. If there are less than two points dened in the contour, this command is ignored. ho selectpoint(x,y) Adds a contour point at the closest input point to (x,y). If the selected point and a previous selected points lie at the ends of an arcsegment, a contour is added that traces along the arcsegment. The selectpoint command has the same functionality as the left-button-click contour point selection when the program is running in interactive mode. ho clearcontour() Clear a prevously dened contour ho clearblock() Clear block selection
9.3 Zoom Commands

ho zoomnatural() Zoom to the natural boundaries of the geometry. ho zoomin() Zoom in one level. ho zoomout() Zoom out one level. ho zoom(x1,y1,x2,y2) Zoom to the window dened by lower left corner (x1,y1) and upper right corner (x2,y2).
9.4 View Commands

ho showmesh() Show the mesh. ho hidemesh() Hide the mesh. ho showpoints() Show the node points from the input geometry. ho hidepoints() Hide the node points from the input geometry.
42
ho smooth("flag") This function controls whether or not smoothing is applied to the F and G elds which are naturally piece-wise constant over each element. Setting ag equal to "on" turns on smoothing, and setting ag to "off" turns off smoothing. ho showgrid() Show the grid points. ho hidegrid() Hide the grid points points. ho gridsnap("flag") Setting ag to on turns on snap to grid, setting ag to off turns off snap to grid. ho setgrid(density,"type") Change the grid spacing. The density parameter species the space between grid points, and the type parameter is set to "cart" for Cartesian coordinates or "polar" for polar coordinates. ho hidedensityplot() hides the heat ux density plot. ho showdensityplot(legend,gscale,upper,lower,type) Shows the heat ux density plot with options: legend Set to 0 to hide the plot legend or 1 to show the plot legend. gscale Set to 0 for a colour density plot or 1 for a grey scale density plot. upper Sets the upper display limit for the density plot. lower Sets the lower display limit for the density plot. type Sets the type of density plot. A value of 0 plots temperature, 1 plots the magnitude of F, and 2 plots the magnitude of G ho hidecontourplot() Hides the contour plot. ho showcontourplot(numcontours,lower V,upper V) shows the V contour plot with options: numcontours Number of equipotential lines to be plotted. upper V Upper limit for contours. lower V Lower limit for contours. If ho numcontours is -1 all parameters are ignored and default values are used, e.g. show contour plot(-1) ho showvectorplot(type,scalefactor) controls the display of vectors denoting the eld strength and direction. The parameters taken are the type of plot, which should be set to 0 for no vector plot, 1 for heat ux density F, and 2 for temperature gradient G. The scalefactor determines the relative length of the vectors. If the scale is set to 1, the length of the vectors are chosen so that the highest ux density corresponds to a vector that is the same length as the current grid size setting. ho minimize() minimizes the active heat ow input view. ho maximize() maximizes the active heat ow input view. 43
ho restore() restores the active heat ow input view from a minimized or maximized state. ho resize(width,height) resizes the active heat ow input window client area to width height.
9.5 Miscellaneous
ho close() close the current postprocessor window. ho refreshview() Redraws the current view. ho reload() Reloads the solution from disk. ho savebitmap("filename") saves a bitmapped screen shot of the current view to the le specied by "filename". Note that if you use a path you must use two backslashes (e.g. "c:\\temp\\myfile.bmp"). If the le name contains a space (e.g. le names like c:\program les\stuff) you must enclose the le name in (extra) quotes by using a \" sequence. For example: ho savebitmap("\"c:\\temp\\screenshot.bmp\"") ho savemetafile("filename") saves a metale screenshot of the current view to the le specied by "filename", subject to the printf-type formatting explained previously for the savebitmap command. ho shownames(flag) This function allow the user to display or hide the block label names on screen. To hide the block label names, flag should be 0. To display the names, the parameter should be set to 1.
10 Current Flow Preprocessor Lua Command Set


ci addnode(x,y) Add a new node at x,y ci addsegment(x1,y1,x2,y2) Add a new line segment from node closest to (x1,y1) to node closest to (x2,y2) ci addblocklabel(x,y) Add a new block label at (x,y) ci addarc(x1,y1,x2,y2,angle,maxseg) Add a new arc segment from the nearest node to (x1,y1) to the nearest node to (x2,y2) with angle angle divided into maxseg segments. ci deleteselected Delete all selected objects. 44
ci deleteselectednodes Delete selected nodes. ci deleteselectedlabels Delete selected block labels. ci deleteselectedsegments Delete selected segments. ci deleteselectedarcsegments Delete selects arcs.

ci clearselected() Clear all selected nodes, blocks, segments and arc segments. ci selectsegment(x,y) Select the line segment closest to (x,y) ci selectnode(x,y) Select the node closest to (x,y). Returns the coordinates of the selected node. ci selectlabel(x,y) Select the label closet to (x,y). Returns the coordinates of the selected label. ci selectarcsegment(x,y) Select the arc segment closest to (x,y) ci selectgroup(n) Select the nth group of nodes, segments, arc segments and block labels. This function will clear all previously selected elements and leave the edit mode in 4 (group)

ci setnodeprop("propname",groupno, "inconductor") Set the selected nodes to have the nodal property "propname" and group number groupno. The "inconductor" string species which conductor the node belongs to. If the node doesnt belong to a named conductor, this parameter can be set to "<None>". ci setblockprop("blockname", automesh, meshsize, group) Set the selected block labels to have the properties: Block property "blockname". automesh: 0 = mesher defers to mesh size constraint dened in meshsize, 1 = mesher automatically chooses the mesh density. meshsize: size constraint on the mesh in the block marked by this label. A member of group number group ci setsegmentprop("propname", elementsize, automesh, hide, group, "inconductor",) Set the select segments to have: Boundary property "propname" Local element size along segment no greater than elementsize automesh: 0 = mesher defers to the element constraint dened by elementsize, 1 = mesher automatically chooses mesh size along the selected segments 45
hide: 0 = not hidden in post-processor, 1 == hidden in post processor A member of group number group A member of the conductor specied by the string "inconductor". If the segment is not part of a conductor, this parameter can be specied as "<None>". ci setarcsegmentprop(maxsegdeg, "propname", hide, group, "inconductor") Set the selected arc segments to: Meshed with elements that span at most maxsegdeg degrees per element Boundary property "propname" hide: 0 = not hidden in post-processor, 1 == hidden in post processor A member of group number group A member of the conductor specied by the string "inconductor". If the segment is not part of a conductor, this parameter can be specied as "<None>".

ci probdef(units,type,frequency,precision,(depth),(minangle)) changes the problem denition. The units parameter species the units used for measuring length in the problem domain. Valid "units" entries are "inches", "millimeters", "centimeters", "mils", "meters, and "micrometers". Set problemtype to "planar" for a 2-D planar problem, or to "axi" for an axisymmetric problem. The frequency parameter species the frequency in Hz at which the analysis ois to be performed. The precision parameter dictates the precision required by the solver. For example, entering 1.E-8 requires the RMS of the residual to be less than 108 . A fourth parameter, representing the depth of the problem in the into-the-page direction for 2-D planar problems, can also be specied for planar problems. A sixth parameter represents the minimum angle constraint sent to the mesh generator. ci analyze(flag) runs belasolv to solve the problem. The flag parameter controls whether the Belasolve window is visible or minimized. For a visible window, either specify no value for flag or specify 0. For a minimized window, flag should be set to 1. ci loadsolution() loads and displays the solution corresponding to the current geometry. ci setfocus("documentname") Switches the electrostatics input le upon which Lua commands are to act. If more than one electrostatics input le is being edited at a time, this command can be used to switch between les so that the mutiple les can be operated upon programmatically via Lua. documentname should contain the name of the desired document as it appears on the windows title bar. ci saveas("filename") saves the le with name "filename". Note if you use a path you must use two backslashes e.g. c:\\temp\\myfemmfile.fee
46
10.5 Mesh Commands

ci createmesh() runs triangle to create a mesh. Note that this is not a necessary precursor of performing an analysis, as ci analyze() will make sure the mesh is up to date before running an analysis. The number of elements in the mesh is pushed back onto the lua stack. ci showmesh() toggles the ag that shows or hides the mesh. ci purgemesh() clears the mesh out of both the screen and memory.

ci copyrotate(bx, by, angle, copies, (editaction) ) bx, by base point for rotation angle angle by which the selected objects are incrementally shifted to make each copy. angle is measured in degrees. copies number of copies to be produced from the selected objects. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group ci copytranslate(dx, dy, copies, (editaction)) dx,dy distance by which the selected objects are incrementally shifted. copies number of copies to be produced from the selected objects. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group mi createradius(x,y,r)turnsacornerlocatedat(x,y)intoacurveofradiusr. ci moverotate(bx,by,shiftangle (editaction)) bx, by base point for rotation shiftangle angle in degrees by which the selected objects are rotated. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group ci movetranslate(dx,dy,(editaction)) dx,dy distance by which the selected objects are shifted. editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group ci scale(bx,by,scalefactor,(editaction)) bx, by base point for scaling scalefactor a multiplier that determines how much the selected objects are scaled editaction 0 nodes, 1 lines (segments), 2 block labels, 3 arc segments, 4- group ci mirror(x1,y1,x2,y2,(editaction)) mirror the selected objects about a line passing through the points (x1,y1) and (x2,y2). Valid editaction entries are 0 for nodes, 1 for lines (segments), 2 for block labels, 3 for arc segments, and 4 for groups. 47
ci seteditmode(editmode) Sets the current editmode to: "nodes" nodes "segments" - line segments "arcsegments" - arc segments "blocks" - block labels "group" - selected group This command will affect all subsequent uses of the other editing commands, if they are used WITHOUT the editaction parameter.
10.7 Zoom Commands

ci zoomnatural() zooms to a natural view with sensible extents. ci zoomout() zooms out by a factor of 50%. ci zoomin() zoom in by a factor of 200%. ci zoom(x1,y1,x2,y2) Set the display area to be from the bottom left corner specied by (x1,y1) to the top right corner specied by (x2,y2).
10.8 View Commands

ci showgrid() Show the grid points. ci hidegrid() Hide the grid points points. ci gridsnap("flag") Setting ag to on turns on snap to grid, setting ag to offturns off snap to grid. ci setgrid(density,"type") Change the grid spacing. The density parameter species the space between grid points, and the type parameter is set to "cart" for Cartesian coordinates or "polar" for polar coordinates. ci refreshview() Redraws the current view. ci minimize() minimizes the active magnetics input view. ci maximize() maximizes the active magnetics input view. ci restore() restores the active magnetics input view from a minimized or maximized state. ci resize(width,height) resizes the active magnetics input window client area to width height.
48

ci addmaterial("materialname", ox, oy, ex, ey, ltx, lty) adds a new material with called "materialname" with the material properties: ox Electrical conductivity in the x- or r-direction in units of S/m. oy Electrical conductivity in the y- or z-direction in units of S/m. ex Relative permittivity in the x- or r-direction. ey Relative permittivity in the y- or z-direction. ltx Dielectric loss tangent in the x- or r-direction. lty Dielectric loss tangent in the y- or z-direction. ci addpointprop("pointpropname",Vp,jp) adds a new point property of name "pointpropname" with either a specied potential Vp a point current density jp in units of A/m. ci addboundprop("boundpropname", Vs, js, c0, c1, BdryFormat) adds a new boundary property with name "boundpropname" For a Fixed Voltage type boundary condition, set the Vs parameter to the desired voltage and all other parameters to zero. To obtain a Mixed type boundary condition, set C1 and C0 as required and BdryFormat to 1. Set all other parameters to zero. To obtain a prescribes surface current density, set js to the desired current density in A/m2 and set BdryFormat to 2. For a Periodic boundary condition, set BdryFormat to 3 and set all other parameters to zero. For an Anti-Perodic boundary condition, set BdryFormat to 4 set all other parameters to zero. ci addconductorprop("conductorname", Vc, jc, conductortype) adds a new conductor property with name "conductorname" with either a prescribed voltage or a prescribed total current. Set the unused property to zero. The conductortype parameter is 0 for prescribed current and 1 for prescribed voltage. ci deletematerial("materialname") deletes the material named "materialname". ci deleteboundprop("boundpropname") deletes the boundary property named "boundpropname". ci deleteconductor("conductorname") deletes the conductor named conductorname. ci deletepointprop("pointpropname") deletes the point property named "pointpropname" ci modifymaterial("BlockName",propnum,value) This function allows for modication of a materials properties without redening the entire material (e.g. so that current can be modied from run to run). The material to be modied is specied by "BlockName". The next parameter is the number of the property to be set. The last number is the value to 49
be applied below: propnum 0 1 2 3 4 5 6
to the specied property. The various properties that can be modied are listed Symbol BlockName ox oy ex ey ltx lty Description Name of the material x (or r) direction conductivity y (or z) direction conductivity x (or r) direction relative permittivity y (or z) direction relative permittivity x (or r) direction dielectric loss tangent y (or z) direction dielectric loss tangent
ci modifyboundprop("BdryName",propnum,value) This function allows for modication of a boundary property. The BC to be modied is specied by "BdryName". The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below: propnum 0 1 2 3 4 5 Symbol BdryName Vs js c0 c1 BdryFormat Description Name of boundary property Fixed Voltage Prescribed current density Mixed BC parameter Mixed BC parameter Type of boundary condition: 0 = Prescribed V 1 = Mixed 2 = Surface current density 3 = Periodic 4 = Antiperiodic
ci modifypointprop("PointName",propnum,value) This function allows for modication of a point property. The point property to be modied is specied by "PointName". The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below: propnum Symbol Description 0 PointName Name of the point property 1 Vp Prescribed nodal voltage 2 jp Point current density in A/m ci modifyconductorprop("ConductorName",propnum,value) This function allows for modication of a conductor property. The conductor property to be modied is specied by "ConductorName". The next parameter is the number of the property to be set. The last number is the value to be applied to the specied property. The various properties that can be modied are listed below:
50
propnum 0 1 2 3
Symbol ConductorName Vc jc ConductorType
Description Name of the conductor property Conductor voltage Total conductor current 0 = Prescribed current, 1 = Prescribed voltage
10.10 Miscellaneous
ci savebitmap("filename") saves a bitmapped screenshot of the current view to the le specied by "filename", subject to the printf-type formatting explained previously for the ci saveas command. ci savemetafile("filename") saves a metale screenshot of the current view to the le specied by "filename", subject to the printf-type formatting explained previously for the ci saveas command. ci refreshview() Redraws the current view. ci close() closes the preprocessor window and destroys the current document. ci shownames(flag) This function allow the user to display or hide the block label names on screen. To hide the block label names, flag should be 0. To display the names, the parameter should be set to 1. ci readdxf("filename") This function imports a dxf le specied by "filename". ci defineouterspace(Zo,Ro,Ri) denes an axisymmetric external region to be used in conjuction with the Kelvin Transformation method of modeling unbounded problems. The Zo parameter is the z-location of the origin of the outer region, the Ro parameter is the radius of the outer region, and the Ri parameter is the radius of the inner region (i.e. the region of interest). In the exterior region, the permeability varies as a function of distance from the origin of the external region. These parameters are necessary to dene the permeability variation in the external region. ci attachouterspace() marks all selected block labels as members of the external region used for modeling unbounded axisymmetric problems via the Kelvin Transformation. ci detachouterspace() undenes all selected block labels as members of the external region used for modeling unbounded axisymmetric problems via the Kelvin Transformation.
11 Current Flow Post Processor Command Set

There are a number of Lua scripting commands designed to operate in the postprocessor. As with the preprocessor commands, these commands can be used with either the underscore naming or with the no-underscore naming convention.
51

co lineintegral(type) Calculate the line integral for the dened contour type 0 1 2 3 4 5 Integral E t Jn Contour length Average voltage over contour Force from stress tensor Torque from stress tensor
This integral returns either 1 or 2 values, depending on the integral type, e.g. : Fx, Fy = co lineintegral(4) co blockintegral(type) Calculate a block integral for the selected blocks type 0 1 2 3 4 5 6 7 8 9 10 11 Integral Real Power Reactive Power Apparent Power Time-Average Stored Energy Block cross-section area Block volume x (or r) direction Weighted Stress Tensor force, DC component y (or z) direction Weighted Stress Tensor force, DC component x (or r) direction Weighted Stress Tensor force, 2x frequency component y (or z) direction Weighted Stress Tensor force, 2x frequency component Weighted Stress Tensor torque, DC component Weighted Stress Tensor torque, 2x frequency component
Returns a value that can be complex, as necessary. co getpointvalues(X,Y) Get the values associated with the point at x,y The return values, in order, are:
52
Symbol V Jx Jy Kx Ky Ex Ey ex ey Jdx Jdy ox oy Jcx Jcy
Denition Voltage x- or r- direction component of current density y- or z- direction component of current density x- or r- direction component of AC conductivity y- or z- direction component of AC conductivity x- or r- direction component of electric eld intensity y- or z- direction component of electric eld intensity x- or r- direction component of permittivity y- or z- direction component of permittivity x- or r- direction component of displacement current density y- or z- direction component of displacement current density x- or r- direction component of permittivity y- or z- direction component of permittivity x- or r- direction component of conduction current density y- or z- direction component of conduction current density
co makeplot(PlotType,NumPoints,Filename,FileFormat) Allows Lua access to the X-Y plot routines. If only PlotType or only PlotType and NumPoints are specied, the command is interpreted as a request to plot the requested plot type to the screen. If, in addition, the Filename parameter is specied, the plot is instead written to disk to the specied le name as an extended metale. If the FileFormat parameter is also, the command is instead interpreted as a command to write the data to disk to the speced le name, rather than display it to make a graphical plot. Valid entries for PlotType are: PlotType 0 1 2 3 4 5 6 7 8 9 10 11 12 Denition V (Voltage) |J| (Magnitude of current density) J . n (Normal current density) J . t (Tangential current density) |E| (Magnitude of field intensity) E . n (Normal field intensity) E . t (Tangential field intensity) |Jc| (Magnitude of conduction current density) Jc . n (Normal conduction current density) Jc . t (Tangential conduction current density) |Jd| (Magnitude of displacement current density) Jd . n (Normal displacement current density) Jd . t (Tangential displacement current density)
Valid le formats are: FileFormat Denition 0 Multi-column text with legend 1 Multi-column text with no legend 2 Mathematica-style formatting For example, if one wanted to plot V to the screen with 200 points evaluated to make the 53
graph, the command would be: co makeplot(0,200) If this plot were to be written to disk as a metale, the command would be: co makeplot(0,200,"c:temp.emf") To write data instead of a plot to disk, the command would be of the form: co makeplot(0,200,"c:temp.txt",0) co getprobleminfo() Returns info on problem description. Returns three values: Return value 1 2 3 Denition problem type frequency in Hz depth assumed for planar problems in meters.
co getconductorproperties("conductor") Properties are returned for the conductor property named conductor. Two values are returned: The voltage of the specied conductor, and the current on the specied conductor.

co seteditmode(mode) Sets the mode of the postprocessor to point, contour, or area mode. Valid entries for mode are "point", "contour", and "area". co selectblock(x,y) Select the block that contains point (x,y). co groupselectblock(n) Selects all of the blocks that are labeled by block labels that are members of group n. If no number is specied (i.e. co groupselectblock()), all blocks are selected. co selectconductor("name") Selects all nodes, segments, and arc segments that are part of the conductor specied by the string ("name"). This command is used to select conductors for the purposes of the weighted stress tensor force and torque integrals, where the conductors are points or surfaces, rather than regions (i.e. cant be selected with co selectblock). co addcontour(x,y) Adds a contour point at (x,y). If this is the rst point then it starts a contour, if there are existing points the contour runs from the previous point to this point. The co addcontour command has the same functionality as a right-button-click contour point addition when the program is running in interactive mode. co bendcontour(angle,anglestep) Replaces the straight line formed by the last two points in the contour by an arc that spans angle degrees. The arc is actually composed of many straight lines, each of which is constrained to span no more than anglestep degrees. The angle parameter can take on values from -180 to 180 degrees. The anglestep parameter must be greater than zero. If there are less than two points dened in the contour, this command is ignored.
54
co selectpoint(x,y) Adds a contour point at the closest input point to (x,y). If the selected point and a previous selected points lie at the ends of an arcsegment, a contour is added that traces along the arcsegment. The selectpoint command has the same functionality as the left-button-click contour point selection when the program is running in interactive mode. co clearcontour() Clear a prevously dened contour co clearblock() Clear block selection
11.3 Zoom Commands

co zoomnatural() Zoom to the natural boundaries of the geometry. co zoomin() Zoom in one level. co zoomout() Zoom out one level. co zoom(x1,y1,x2,y2) Zoom to the window dened by lower left corner (x1,y1) and upper right corner (x2,y2).
11.4 View Commands

co showmesh() Show the mesh. co hidemesh() Hide the mesh. co showpoints() Show the node points from the input geometry. co hidepoints() Hide the node points from the input geometry. co smooth("flag") This function controls whether or not smoothing is applied to the D and E elds which are naturally piece-wise constant over each element. Setting ag equal to "on" turns on smoothing, and setting ag to "off" turns off smoothing. co showgrid() Show the grid points. co hidegrid() Hide the grid points points. co gridsnap("flag") Setting ag to on turns on snap to grid, setting ag to off turns off snap to grid. co setgrid(density,"type") Change the grid spacing. The density parameter species the space between grid points, and the type parameter is set to "cart" for Cartesian coordinates or "polar" for polar coordinates. co hidedensityplot() hides the current density plot.
55
co showdensityplot(legend,gscale,upper,lower,type) Shows the current density plot with options: legend Set to 0 to hide the plot legend or 1 to show the plot legend. gscale Set to 0 for a colour density plot or 1 for a grey scale density plot. upper Sets the upper display limit for the density plot. lower Sets the lower display limit for the density plot. type Sets the type of density plot. Specic choices for the type of density plot include: type 0 1 2 3 4 5 6 7 8 Description |V | |Re(V )| |Im(V )| |J| |Re(J)| |Im(J)| |E| |Re(E)| |Im(E)|
co hidecontourplot() Hides the contour plot. co showcontourplot(numcontours,lower V,upper V),type shows the V contour plot with options: numcontours Number of equipotential lines to be plotted; upper V Upper limit for contours; lower V Lower limit for contours; type the type of contour plot to be rendered. If co numcontours is -1 all parameters are ignored and default values are used, e.g. show contour plot(-1) The type can take on the values of "real", "imag", or "both", denoting the real part of voltage, the imaginary part of voltage, or both components of voltage. co showvectorplot(type,scalefactor) controls the display of vectors denoting the eld strength and direction. The type parameter can take on the following values: type 0 1 2 3 4 5 6 Description No vector plot Re(J) Re(E) The scalefactor determines the relative length of the vectors. Im(J) Im(E) Re(J) and Im(J) Re(E) and Im(E) 56
If the scale is set to 1, the length of the vectors are chosen so that the highest eld magnitude corresponds to a vector that is the same length as the current grid size setting. co minimize() minimizes the active magnetics input view. co maximize() maximizes the active magnetics input view. co restore() restores the active magnetics input view from a minimized or maximized state. co resize(width,height) resizes the active magnetics input window client area to width height.
11.5 Miscellaneous
co close() close the current postprocessor window. co refreshview() Redraws the current view. co reload() Reloads the solution from disk. co savebitmap("filename") saves a bitmapped screen shot of the current view to the le specied by "filename". Note that if you use a path you must use two backslashes (e.g. "c:\\temp\\myfile.bmp"). If the le name contains a space (e.g. le names like c:\program les\stuff) you must enclose the le name in (extra) quotes by using a \" sequence. For example: co savebitmap("\"c:\\temp\\screenshot.bmp\"") co savemetafile("filename") saves a metale screenshot of the current view to the le specied by "filename", subject to the printf-type formatting explained previously for the savebitmap command. co shownames(flag) This function allow the user to display or hide the block label names on screen. To hide the block label names, flag should be 0. To display the names, the parameter should be set to 1.
57

Finite Element Method Magnetics: Octavefemm: User'S Manual

Uploaded by

Copyright:

Available Formats

Finite Element Method Magnetics: Octavefemm: User'S Manual

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Finite Element Method Magnetics: Octavefemm: User'S Manual

Uploaded by

Copyright:

Available Formats

Finite Element Method Magnetics: OctaveFEMM

Version 1.2 Users Manual

David Meeker [email protected] http://femm.foster-miller.com c 2006

2 Installation and Startup

3 Common Command Set

4 Magnetics Preprocessor Command Set

4.1 Object Add/Remove Commands

mi deleteselectedsegments Delete selected segments. mi deleteselectedarcsegments Delete selects arcs.

4.2 Geometry Selection Commands

4.3 Object Labeling Commands

4.4 Problem Commands

4.5 Mesh Commands

4.6 Editing Commands

4.7 Zoom Commands

4.8 View Commands

4.9 Object Properties

Symbol BdryName A0 A1 A2 c0 c1 BdryFormat

mi setcurrent(CircName,i) sets the current in the circuit specied by CircName to i.

5 Magnetics Post Processor Command Set

5.1 Data Extraction Commands

Symbol A B1 B2 Sig E H1 H2 Je Js Mu1 Mu2 Pe Ph ff

5.2 Selection Commands

5.3 Zoom Commands

5.4 View Commands

6 Electrostatics Preprocessor Command Set

6.1 Object Add/Remove Commands

6.2 Geometry Selection Commands

6.3 Object Labeling Commands

6.4 Problem Commands

6.5 Mesh Commands

6.6 Editing Commands

6.7 Zoom Commands

6.8 View Commands

6.9 Object Properties

7 Electrostatics Post Processor Command Set

7.1 Data Extraction Commands

7.2 Selection Commands

7.3 Zoom Commands

7.4 View Commands

8 Heat Flow Preprocessor Lua Command Set

8.1 Object Add/Remove Commands

hi deleteselectedsegments Delete selected segments. hi deleteselectedarcsegments Delete selects arcs.

8.2 Geometry Selection Commands

8.3 Object Labeling Commands

8.4 Problem Commands

8.5 Mesh Commands

8.6 Editing Commands

8.7 Zoom Commands

8.8 View Commands

8.9 Object Properties

be applied below: propnum 0 1 2 3

Symbol ConductorName Tc qc ConductorType

9 Heat Flow Post Processor Command Set

9.1 Data Extraction Commands

9.2 Selection Commands

9.3 Zoom Commands

9.4 View Commands

10 Current Flow Preprocessor Lua Command Set

10.1 Object Add/Remove Commands