Abaqus 6.11: Abaqus Scripting User's Manual
Abaqus 6.11: Abaqus Scripting User's Manual
Abaqus 6.11
Scripting User’s Manual
Abaqus ID:
Printed on:
Abaqus Scripting User’s Manual
Abaqus ID:
Printed on:
Legal Notices
CAUTION: This documentation is intended for qualified users who will exercise sound engineering judgment and expertise in the use of the Abaqus
Software. The Abaqus Software is inherently complex, and the examples and procedures in this documentation are not intended to be exhaustive or to apply
to any particular situation. Users are cautioned to satisfy themselves as to the accuracy and results of their analyses.
Dassault Systèmes and its subsidiaries, including Dassault Systèmes Simulia Corp., shall not be responsible for the accuracy or usefulness of any analysis
performed using the Abaqus Software or the procedures, examples, or explanations in this documentation. Dassault Systèmes and its subsidiaries shall not
be responsible for the consequences of any errors or omissions that may appear in this documentation.
The Abaqus Software is available only under license from Dassault Systèmes or its subsidiary and may be used or reproduced only in accordance with the
terms of such license. This documentation is subject to the terms and conditions of either the software license agreement signed by the parties, or, absent
such an agreement, the then current software license agreement to which the documentation relates.
This documentation and the software described in this documentation are subject to change without prior notice.
No part of this documentation may be reproduced or distributed in any form without prior written permission of Dassault Systèmes or its subsidiary.
The Abaqus Software is a product of Dassault Systèmes Simulia Corp., Providence, RI, USA.
© Dassault Systèmes, 2011
Abaqus, the 3DS logo, SIMULIA, CATIA, and Unified FEA are trademarks or registered trademarks of Dassault Systèmes or its subsidiaries in the United
States and/or other countries.
Other company, product, and service names may be trademarks or service marks of their respective owners. For additional information concerning
trademarks, copyrights, and licenses, see the Legal Notices in the Abaqus 6.11 Release Notes.
Abaqus ID:
Printed on:
Locations
SIMULIA Worldwide Headquarters Rising Sun Mills, 166 Valley Street, Providence, RI 02909–2499, Tel: +1 401 276 4400,
Fax: +1 401 276 4408, simulia.support@3ds.com, http://www.simulia.com
SIMULIA European Headquarters Stationsplein 8-K, 6221 BT Maastricht, The Netherlands, Tel: +31 43 7999 084,
Fax: +31 43 7999 306, simulia.europe.info@3ds.com
Abaqus ID:
Printed on:
Preface
This section lists various resources that are available for help with using Abaqus Unified FEA software.
Support
Both technical engineering support (for problems with creating a model or performing an analysis) and
systems support (for installation, licensing, and hardware-related problems) for Abaqus are offered through
a network of local support offices. Regional contact information is listed in the front of each Abaqus manual
and is accessible from the Locations page at www.simulia.com.
Training
All offices and representatives offer regularly scheduled public training classes. The courses are offered in
a traditional classroom form and via the Web. We also provide training seminars at customer sites. All
training classes and seminars include workshops to provide as much practical experience with Abaqus as
possible. For a schedule and descriptions of available classes, see www.simulia.com or call your local office
or representative.
Feedback
We welcome any suggestions for improvements to Abaqus software, the support program, or documentation.
We will ensure that any enhancement requests you make are considered for future releases. If you wish to
make a suggestion about the service or products, refer to www.simulia.com. Complaints should be addressed
by contacting your local office or through www.simulia.com by visiting the Quality Assurance section of
the Support page.
Abaqus ID:
Printed on:
CONTENTS
Contents
3. Simple examples
Creating a part 3.1
Reading from an output database 3.2
Summary 3.3
4. Introduction to Python
Python and Abaqus 4.1
Python resources 4.2
Using the Python interpreter 4.3
Object-oriented basics 4.4
The basics of Python 4.5
Programming techniques 4.6
Further reading 4.7
Abaqus ID:cmd-toc
Printed on: Mon February 28 -- 12:29:52 2011
CONTENTS
ii
Abaqus ID:cmd-toc
Printed on: Mon February 28 -- 12:29:52 2011
CONTENTS
iii
Abaqus ID:cmd-toc
Printed on: Mon February 28 -- 12:29:52 2011
Part I: An introduction to the Abaqus Scripting
Interface
The Abaqus Scripting Interface is an application programming interface (API) to the models and data used by
Abaqus. The Abaqus Scripting Interface is an extension of the Python object-oriented programming language;
Abaqus Scripting Interface scripts are Python scripts. You can use the Abaqus Scripting Interface to do the
following:
• Create and modify the components of an Abaqus model, such as parts, materials, loads, and steps.
• Create, modify, and submit Abaqus analysis jobs.
• Read from and write to an Abaqus output database.
• View the results of an analysis.
You use the Abaqus Scripting Interface to access the functionality of Abaqus/CAE from scripts (or programs).
(The Visualization module of Abaqus/CAE is also licensed separately as Abaqus/Viewer; therefore, the
Abaqus Scripting Interface can also be used to access the functionality of Abaqus/Viewer.) Because the
Abaqus Scripting Interface is a customized extension of standard Python, further extension of Abaqus base
types to create user-defined classes is not allowed.
This section provides an introduction to the Abaqus Scripting Interface. The following topics are
covered:
• Chapter 1, “An overview of the Abaqus Scripting User’s Manual”
• Chapter 2, “Introduction to the Abaqus Scripting Interface”
• Chapter 3, “Simple examples”
Abaqus ID:
Printed on:
AN OVERVIEW OF THE Abaqus SCRIPTING USER’S MANUAL
Simple examples
Two simple examples are provided to introduce you to programming with the Abaqus Scripting
Interface.
• Creating a part.
• Reading from an output database.
An introduction to Python
This section is intended as a basic introduction to the Python programming language and is not an
exhaustive description of the language. There are several books on the market that describe Python,
and these books are listed as references. Additional resources, such as Python-related sites, are also
listed.
Example scripts
This section provides a set of example scripts that lead you through the cantilever beam tutorial
found in Appendix B, “Creating and Analyzing a Simple Model in Abaqus/CAE,” of Getting Started
with Abaqus: Interactive Edition. Additional examples are provided that read from an output
database, display a contour plot, and print a contour plot from each step of the analysis. The final
example illustrates how you can read from a model database created by Abaqus/CAE, parameterize
the model, submit a set of analysis jobs, and generate results from the resulting output databases.
1–1
Abaqus ID:
Printed on:
AN OVERVIEW OF THE Abaqus SCRIPTING USER’S MANUAL
1–2
Abaqus ID:
Printed on:
Abaqus/CAE AND THE Abaqus SCRIPTING INTERFACE
2–1
Abaqus ID:
Printed on:
HOW DOES THE Abaqus SCRIPTING INTERFACE INTERACT WITH Abaqus/CAE?
Abaqus/CAE
command
line
GUI interface script
(CLI)
commands
Python replay
interpreter files
Abaqus/CAE
kernel
input file
Abaqus/Standard
Abaqus/Explicit
Abaqus/CFD
Abaqus/Design
output database
Abaqus Scripting Interface commands can be issued to the Abaqus/CAE kernel from one of the
following:
• The graphical user interface (GUI). For example, when you click OK or Apply in a dialog box,
the GUI generates a command based on your options and settings in the dialog box. You can use
the Macro Manager to record a sequence of the generated Abaqus Scripting Interface commands
2–2
Abaqus ID:
Printed on:
HOW DOES THE Abaqus SCRIPTING INTERFACE INTERACT WITH Abaqus/CAE?
in a macro file. For more information, see “Creating and running a macro,” Section 9.5.5 of the
Abaqus/CAE User’s Manual.
• Click in the lower left corner of the main window to display the command line interface (CLI).
You can type a single command or paste in a sequence of commands from another window; the
command is executed when you press [Enter]. You can type any Python command into the command
line; for example, you can use the command line as a simple calculator.
Note: When you are using Abaqus/CAE, errors and messages are posted into the message area.
Click in the lower left corner of the main window to display the message area.
• If you have more than a few commands to execute or if you are repeatedly executing the same
commands, it may be more convenient to store the set of statements in a file called a script. A
script contains a sequence of Python statements stored in plain ASCII format. For example, you
might create a script that opens an output database, displays a contour plot of a selected variable,
customizes the legend of the contour plot, and prints the resulting image on a local PostScript printer.
In addition, scripts are useful for starting Abaqus/CAE in a predetermined state. For example, you
can define a standard configuration for printing, create remote queues, and define a set of standard
materials and their properties.
You can use one of the following methods to run a script:
Running a script when you start Abaqus/CAE
You can run a script when you start an Abaqus/CAE session by typing the following command:
abaqus cae script=myscript.py
where myscript.py is the name of the file containing the script. The equivalent command
for Abaqus/Viewer is
abaqus viewer script=myscript.py
Arguments can be passed into the script by entering -- on the command line, followed by
the arguments separated by one or more spaces. These arguments will be ignored by the
Abaqus/CAE execution procedure, but they will be accessible within the script. For more
information, see “Abaqus/CAE execution,” Section 3.2.4 of the Abaqus Analysis User’s
Manual, and “Abaqus/Viewer execution,” Section 3.2.5 of the Abaqus Analysis User’s
Manual.
2–3
Abaqus ID:
Printed on:
HOW DOES THE Abaqus SCRIPTING INTERFACE INTERACT WITH Abaqus/CAE?
The Abaqus/CAE kernel is started without the GUI. Running a script without the
Abaqus/CAE GUI is useful for automating pre- or postanalysis processing tasks without the
added expense of running a display. When the script finishes running, the Abaqus/CAE kernel
terminates. If you execute a script without the GUI, the script cannot interact with the user,
monitor jobs, or generate animations. When running a script without the user interface, jobs
are always run interactively. If a job queue is specified, it will be ignored.
execfile('myscript.py')
where myscript.py is the name of the file containing the script and the file in this example
is in the current directory. Figure 2–2 shows an example script being run from the command
line interface.
Figure 2–2 Scripts can be run from the command line interface.
Click in the lower left corner of the main window to switch from the message area to the
command line interface.
2–4
Abaqus ID:
Printed on:
CREATING A PART
3. Simple examples
Programming with the Abaqus Scripting Interface is straightforward and logical. To illustrate how easy
it is to write your own programs, the following sections describe two simple Abaqus Scripting Interface
scripts.
• “Creating a part,” Section 3.1
• “Reading from an output database,” Section 3.2
You are not expected to understand every line of the examples; the terminology and the syntax will
become clearer as you read the detailed explanations in the following chapters. “Summary,” Section 3.3,
describes some of the principles behind programming with Python and the Abaqus Scripting Interface.
The first example shows how you can use an Abaqus/CAE script to replicate the functionality of
Abaqus/CAE. The script does the following:
• Creates a new model in the model database.
• Creates a two-dimensional sketch.
• Creates a three-dimensional, deformable part.
• Extrudes the two-dimensional sketch to create the first geometric feature of the part.
• Creates a new viewport.
• Displays a shaded image of the new part in the new viewport.
The new viewport and the shaded part are shown in Figure 3–1.
The example scripts from this manual can be copied to the user’s working directory by using the
Abaqus fetch utility:
where scriptName.py is the name of the script (see “Fetching sample input files,” Section 3.2.13 of the
Abaqus Analysis User’s Manual). Use the following command to retrieve the script for this example:
Note: Abaqus does not install the sample scripts by default during the installation procedure. As a result,
if the Abaqus fetch utility fails to find the sample script, the script may be missing from your Abaqus
installation. You must rerun the installation procedure and request Abaqus sample problems from
the list of items to install.
3–1
Abaqus ID:
Printed on:
CREATING A PART
3–2
Abaqus ID:
Printed on:
CREATING A PART
import sketch
import part
for i in range(len(xyCoordsInner)-1):
mySketch.Line(point1=xyCoordsInner[i],
point2=xyCoordsInner[i+1])
for i in range(len(xyCoordsOuter)-1):
mySketch.Line(point1=xyCoordsOuter[i],
point2=xyCoordsOuter[i+1])
myPart.BaseSolidExtrude(sketch=mySketch, depth=20.0)
myViewport.setValues(displayedObject=myPart)
myViewport.partDisplay.setValues(renderStyle=SHADED)
3–3
Abaqus ID:
Printed on:
CREATING A PART
This statement makes the basic Abaqus objects accessible to the script. It also provides access to a default
model database using the variable named mdb. The statement, from abaqusConstants import
*, makes the Symbolic Constants defined by the Abaqus Scripting Interface available to the script.
import sketch
import part
These statements provide access to the objects related to sketches and parts. sketch and part are
called Python modules.
The next statement in the script is shown in Figure 3–2. You can read this statement from right to
left as follows:
1. Create a new model named Model A.
2. Store the new model in the model database mdb.
3. Assign the new model to a variable called myModel.
3–4
Abaqus ID:
Printed on:
CREATING A PART
This loop creates the inner profile of the letter “A” in mySketch. Four lines are created using the X-
and Y-coordinates of the vertices in xyCoordsInner to define the beginning point and the end point
of each line. Note the following:
• Python uses only indentation to signify the start and the end of a loop. Python does not use the
brackets {} of C and C++.
• The len() function returns the number of coordinate pairs in xyCoordsInner—five in our
example.
• The range() function returns a sequence of integers. In Python, as in C and C++, the index of an
array starts at zero. In our example range(4) returns 0,1,2,3. For each iteration of the loop
in the example the variable i is assigned to the next value in the sequence of integers.
Similarly, a second loop creates the outer profile of the “A” character.
myPart.BaseSolidExtrude(sketch=mySketch, depth=20.0)
This statement creates a base solid extrude feature in myPart by extruding mySketch through a depth
of 20.0.
3–5
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
myViewport.setValues(displayedObject=myPart)
This statement tells Abaqus to display the new part, myPart, in the new viewport, myViewport.
myViewport.partDisplayOptions.setValues(renderStyle=SHADED)
This statement sets the render style of the part display options in myViewport to shaded. As a result,
myPart appears in the shaded render style.
The second example shows how you can use the Abaqus Scripting Interface to read an output database,
manipulate the data, and display the results using the Visualization module in Abaqus/CAE. The Abaqus
Scripting Interface allows you to display the data even though you have not written it back to an output
database. The script does the following:
• Opens the tutorial output database.
• Creates variables that refer to the first and second steps in the output database.
• Creates variables that refer to the last frame of the first and second steps.
• Creates variables that refer to the displacement field output in the last frame of the first and second
steps.
• Creates variables that refer to the stress field output in the last frame of the first and second steps.
• Subtracts the displacement field output from the two steps and puts the result in a variable called
deltaDisplacement.
• Subtracts the stress field output from the two steps and puts the result in a variable called
deltaStress.
• Selects deltaDisplacement as the default deformed variable.
• Selects the von Mises invariant of deltaStress as the default field output variable.
• Plots a contour plot of the result.
The resulting contour plot is shown in Figure 3–3.
Use the following commands to retrieve the script and the output database that is read by the script:
abaqus fetch job=odbExample
abaqus fetch job=viewer_tutorial
3–6
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
S − S, Mises
(Avg: 75%)
+4.012e−02
+3.678e−02
+3.344e−02
+3.011e−02
+2.677e−02
+2.343e−02
+2.010e−02
+1.676e−02
+1.342e−02
+1.009e−02
+6.750e−03
+3.414e−03
+7.711e−05
myOdb = visualization.openOdb(path='viewer_tutorial.odb')
3–7
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
myViewport.setValues(displayedObject=myOdb)
firstStep = myOdb.steps['Step-1']
secondStep = myOdb.steps['Step-2']
frame1 = firstStep.frames[-1]
frame2 = secondStep.frames[-1]
displacement1 = frame1.fieldOutputs['U']
displacement2 = frame2.fieldOutputs['U']
stress1 = frame1.fieldOutputs['S']
stress2 = frame2.fieldOutputs['S']
myViewport.odbDisplay.setDeformedVariable(deltaDisplacement)
myViewport.odbDisplay.setPrimaryVariable(field=deltaStress,
outputPosition=INTEGRATION_POINT,
refinement=(INVARIANT, 'Mises'))
myViewport.odbDisplay.display.setValues(plotState=(CONTOURS_ON_DEF,))
3–8
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
import visualization
This statement provides access to the commands that replicate the functionality of the Visualization
module in Abaqus/CAE (Abaqus/Viewer).
odbPath = 'viewer_tutorial.odb'
This statement creates a path to the tutorial output database.
myOdb = session.openOdb(path=odbPath)
This statement uses the path variable odbPath to open the output database and to assign it to the variable
myOdb.
myViewport.setValues(displayedObject=myOdb)
This statement displays the default plot of the output database in the viewport.
firstStep = myOdb.steps['Step-1']
secondStep = myOdb.steps['Step-2']
These statements assign the first and second steps in the output database to the variables firstStep
and secondStep.
frame1 = firstStep.frames[-1]
frame2 = secondStep.frames[-1]
These statements assign the last frame of the first and second steps to the variables frame1 and frame2.
In Python an index of 0 refers to the first item in a sequence. An index of −1 refers to the last item in a
sequence.
3–9
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
displacement1 = frame1.fieldOutputs['U']
displacement2 = frame2.fieldOutputs['U']
These statements assign the displacement field output in the last frame of the first and second steps to
the variables displacement1 and displacement2.
stress1 = frame1.fieldOutputs['S']
stress2 = frame2.fieldOutputs['S']
Similarly, these statements assign the stress field output in the last frame of the first and second steps to
the variables stress1 and stress2.
This statement subtracts the displacement field output from the last frame of the two steps and puts the
resulting field output into a new variable deltaDisplacement.
Similarly, this statement subtracts the stress field output and puts the result in the variable
deltaStress.
myViewport.odbDisplay.setDeformedVariable(deltaDisplacement)
This statement uses deltaDisplacement, the displacement field output variable that we created
earlier, to set the deformed variable. This is the variable that Abaqus will use to display the shape of the
deformed model.
myViewport.odbDisplay.setPrimaryVariable(field=deltaStress,
outputPosition=INTEGRATION_POINT,
refinement=(INVARIANT, 'Mises'))
This statement uses deltaStress, our stress field output variable, to set the primary variable. This is
the variable that Abaqus will display in a contour or symbol plot.
myViewport.odbDisplay.display.setValues(plotState=(CONTOURS_ON_DEF,))
The final statement sets the plot state to display a contour plot on the deformed model shape.
3–10
Abaqus ID:
Printed on:
SUMMARY
3.3 Summary
The examples illustrate how a script can operate on a model in a model database or on the data stored in
an output database. The details of the commands in the examples are described in later sections; however,
you should note the following:
• You can run a script from the Abaqus/CAE startup screen when you start a session. After a session
has started, you can run a script from the File→Run Script menu or from the command line
interface.
• A script is a sequence of commands stored in ASCII format and can be edited with a standard text
editor.
• A set of example scripts are provided with Abaqus. Use the abaqus fetch command to retrieve
a script and any associated files.
• You must use the import statement to make the required set of Abaqus Scripting Interface
commands available. For example, the statement import part provides the commands that
create and operate on parts.
• A command that creates something (an “object” in object-oriented programming terms) is called a
constructor and starts with an uppercase character. For example, the following statement uses the
Model constructor to create a model object.
• You can use a variable to refer to an object. Variables make your scripts easier to read and
understand. myModel refers to a model object in the previous example.
• A Python script can include a loop. The start and end of a loop is controlled by indentation in the
script.
• Python includes a set of built-in functions. For example, the len() function returns the length of
a sequence.
• You can use commands to replicate any operation that can be performed interactively when you are
working with Abaqus/CAE; for example, creating a viewport, displaying a contour plot, and setting
the step and the frame to display.
3–11
Abaqus ID:
Printed on:
Part II: Using the Abaqus Scripting Interface
This section provides an introduction to the Python programming language and a discussion of how you can
combine Python statements and the Abaqus Scripting Interface to create your own scripts. The following
topics are covered:
• Chapter 4, “Introduction to Python”
• Chapter 5, “Using Python and the Abaqus Scripting Interface”
• Chapter 6, “Using the Abaqus Scripting Interface with Abaqus/CAE”
Abaqus ID:
Printed on:
PYTHON AND Abaqus
4. Introduction to Python
This section provides a basic introduction to the Python programming language. You are encouraged to
try the examples and to experiment with Python statements. The Python language is used throughout
Abaqus, not only in the Abaqus Scripting Interface. Python is also used by Abaqus/Design to perform
parametric studies and in the Abaqus/Standard, Abaqus/Explicit, Abaqus/CFD, and Abaqus/CAE
environment file (abaqus_v6.env). For more information, see Chapter 19, “Parametric Studies,”
of the Abaqus Analysis User’s Manual, and “Using the Abaqus environment settings,” Section 3.3.1
of the Abaqus Analysis User’s Manual.
The following topics are covered:
4–1
Abaqus ID:
Printed on:
PYTHON RESOURCES
4–2
Abaqus ID:
Printed on:
USING THE PYTHON INTERPRETER
Python newsgroups
Discussions of Python programming can be found at:
• comp.lang.python
• comp.lang.python.announce
• If you have Abaqus installed on your UNIX or Windows workstation, type abaqus python at
the system prompt. Python enters its interpretive mode and displays the >>> prompt.
You can enter Python statements at the >>> prompt. To see the value of a variable or expression,
type the variable name or expression at the Python prompt. To exit the Python interpreter, type
[Ctrl]+D on UNIX systems or [Ctrl]+Z[Enter] on Windows systems.
You can also use Python to run a script directly by typing abaqus python scriptname.py
at the system prompt. Abaqus will run the script through the Python interpreter and return you to the
system prompt. For an example of running a Python script using Abaqus, see “Creating functions,”
Section 4.6.1.
• You can also use the Python interpreter provided in the command line interface by Abaqus/CAE.
The command line is at the bottom of the Abaqus/CAE window and is shared with the message
area. Abaqus/CAE displays the Python >>> prompt in the command line interface.
Click in the lower left corner of the main window to display the command line interface.
You may want to drag the handle at the top of the command line interface to increase the number
of lines displayed.
4–3
Abaqus ID:
Printed on:
OBJECT-ORIENTED BASICS
If Abaqus/CAE displays new messages while you are using the command line interface, the message
area tab turns red.
4–4
Abaqus ID:
Printed on:
THE BASICS OF PYTHON
The following sections introduce you to the basics of the Python language. The following topics are
covered:
• “Variable names and assignment,” Section 4.5.1
• “Python data types,” Section 4.5.2
• “Determining the type of a variable,” Section 4.5.3
• “Sequences,” Section 4.5.4
• “Sequence operations,” Section 4.5.5
• “Python None,” Section 4.5.6
• “Continuation lines and comments,” Section 4.5.7
• “Printing variables using formatted output,” Section 4.5.8
• “Control blocks,” Section 4.5.9
4–5
Abaqus ID:
Printed on:
THE BASICS OF PYTHON
The second line assigns the string ’Einstein’ to the variable myName and assigns the string
’Newton’ to the variable yourName. The third line assigns the string ’Einstein’ to both
myName and yourName.
The following naming rules apply:
• Variable names must start with a letter or an underscore character and can contain any number of
letters, digits, or underscores. load_3 and _frictionStep are legal variable names; 3load,
load_3$, and $frictionStep are not legal names. There is no limit on the length of a variable
name.
• Some words are reserved and cannot be used to name a variable; for example, print, while,
return, and class.
• Python is case sensitive. A variable named Load is different from a variable named load.
When you assign a variable in a Python program, the variable refers to a Python object, but the
variable is not an object itself. For example, the expression numSpokes=3 creates a variable that refers
to an integer object; however, numSpokes is not an object. You can change the object to which a
variable refers. numSpokes can refer to a real number on one line, an integer on the next line, and a
viewport on the next line.
The first example script in “Creating a part,” Section 3.1, created a model using the following
statement:
myModel = mdb.Model(name='Model A')
The constructor mdb.Model(name=’Model A’) creates an instance of a model, and this instance is
a Python object. The object created is mdb.models[’Model A’], and the variable myModel refers
to this object.
An object always has a type. In our example the type of mdb.models[’Model A’] is
Model. An object’s type cannot be changed. The type defines the data encapsulated by an object—its
members—and the functions that can manipulate those data—its methods. Unlike most programming
languages, you do not need to declare the type of a variable before you use it. Python determines the
type when the assignment statement is executed. The Abaqus Scripting Interface uses the term “object”
to refer to a specific Abaqus type as well as to an instance of that type; for example, a Model object
refers to a Model type and to an instance of a Model type.
4–6
Abaqus ID:
Printed on:
THE BASICS OF PYTHON
An integer is based on a C long and can be compared to a FORTRAN integer*4 or *8. For extremely
large integer values, you should declare a long integer. The size of a long integer is essentially
unlimited. The “L” at the end of the number indicates that it is a long integer.
>>> nodes = 2000000L
>>> bigNumber = 120L**21
Use int(n) to convert a variable to an integer; use long(n) to convert a variable to a long integer.
>>> load = 279.86
>>> iLoad = int(load)
>>> iLoad
279
>>> a = 2
>>> b = 64
>>> bigNumber = long(a)**b
>>> print 'bigNumber = ', bigNumber
bigNumber = 18446744073709551616
Note: All Abaqus Scripting Interface object types begin with an uppercase character; for example,
a Part or a Viewport. An integer is another kind of object and follows the same convention. The
Abaqus Scripting Interface refers to an integer object as an “Int.” Similarly, the Abaqus Scripting
Interface refers to a floating-point object as a “Float.”
Float
Floats represent floating-point numbers or real numbers. You can use exponential notation for floats.
>>> pi = 22.0/7.0
>>> r = 2.345e-6
>>> area = pi * r * r
>>> print 'Area = ', area
Area = 1.728265e-11
A float is based on a C double and can be compared to a FORTRAN real*8. Use float(n) to
convert a variable to a float.
Complex
Complex numbers use the “j” notation to indicate the imaginary part of the number. Python provides
methods to manipulate complex numbers. The conjugate method calculates the conjugate of a
complex number.
>>> a = 2 + 4j
>>> a.conjugate()
(2-4j)
4–7
Abaqus ID:
Printed on:
THE BASICS OF PYTHON
A complex number has two members, the real member and the imaginary member.
>>> a = 2 + 4j
>>> a.real
2.0
>>> a.imag
4.0
Python provides complex math functions to operate on complex variables. You need to import the
cmath module to use the complex square root function.
>>> import cmath
>>> y = 3 + 4j
>>> print cmath.sqrt(y)
(2+1j)
Remember, functions of a type are called methods; data of a type are called members. In our example
conjugate is a method of a complex type; a.real refers to the real member of a complex type.
Sequences
Sequences include strings, lists, tuples, and arrays. Sequences are described in “Sequences,”
Section 4.5.4, and “Sequence operations,” Section 4.5.5.
4.5.4 Sequences
Sequences are important and powerful data types in Python. A sequence is an object containing a
series of objects. There are three types of built-in sequences in Python—list, tuple, and string. In
addition, imported modules allow you to use arrays in your scripts. The following table describes the
characteristics of list, tuple, string, and array sequences.
4–8
Abaqus ID:
Printed on:
THE BASICS OF PYTHON
List
Lists are mutable heterogeneous sequences (anything that can be modified is called mutable). A list
can be a sequence of strings, integers, floats, or any combination of these. In fact, a list can contain
any type of object; for example,
>>> myIntegerList = [7,6,5,4]
>>> myFloatList = [7.1,6.5,5.3,4.8]
You can refer to individual items from a sequence using the index of the item. Indices start at
zero. Negative indices count backward from the end of a sequence.
>>> myList = [1,2,3]
>>> myList[0]
1
>>> myList[1] = 9
>>> myList
[1, 9, 3]
>>> myNewList = [1.0,2.0,myList]
>>> myNewList
[1.0, 2.0, [1, 9, 3]]
>>> myNewList[-1]
[1, 9, 3]
Lists are heterogeneous, which means they can contain objects of different type.
>>> myList=[1,2.5,'steel']
A list can contain other lists.
>>> myList=[[0,1,2],[3,4,5],[6,7,8]]
4–9
Abaqus ID:
Printed on:
THE BASICS OF PYTHON
>>> myList[0]
[0, 1, 2]
>>> myList[2]
[6,7,8]
myList[1][2] refers to the third item in the second list. Remember, indices start at zero.
>>> myList[1][2]
5
Python has built-in methods that allow you to operate on the items in a sequence.
>>> myList
[1, 9, 3]
>>> myList.append(33)
>>> myList
[1, 9, 3, 33]
>>> myList.remove(9)
>>> myList
[1, 3, 33]
The following are some additional built-in methods that operate on lists:
count()
Return the number of times a value appears in the list.
index()
Return the index indicating the first time an item appears in the list.
>>> myList.index(5)
11
>>> myList.index(4)
8
insert()
Insert a new element into a list at a specified location.
>>> myList.insert(2,22)
>>> myList
[0, 1, 22, 2, 1, 2, 3, 2, 3, 4, 3, 4, 5]
4–10
Abaqus ID:
Printed on:
THE BASICS OF PYTHON
reverse()
Reverse the elements in a list.
>>> myList.reverse()
>>> myList
[5, 4, 3, 4, 3, 2, 3, 2, 1, 2, 22, 1, 0]
sort()
Sort the elements in a list.
>>> myList.sort()
>>> myList
[0, 1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 5, 22]
Tuple
Tuples are very similar to lists; however, they are immutable heterogeneous sequences, which means
that you cannot change them after you create them. You can think of a tuple as a list that cannot be
modified. Tuples have no methods; you cannot append items to a tuple, and you cannot modify or
delete the items in a tuple. The following statement creates an empty tuple:
myTuple = ()
The following statement creates a tuple with one element:
myTuple = (5.675,)
You can use the tuple() function to convert a list or a string to a tuple.
>>> myList = [1, 2, "stress", 4.67]
>>> myTuple = tuple(myList)
>>> print myTuple
(1, 2, 'stress', 4.67)
>>> myString = 'Failure mode'
>>> myTuple = tuple(myString)
>>> print myTuple
('F', 'a', 'i', 'l', 'u', 'r', 'e', ' ', 'm', 'o', 'd', 'e')
The following statements create a tuple and then try to change the value of an item in the tuple.
An AttributeError error message is generated because a tuple is immutable.
>>> myTuple = (1,2,3,4,5)
>>> type(myTuple)
<type 'tuple'>
>>> myTuple[2]=3
Traceback (innermost last):
4–11
Abaqus ID:
Printed on:
THE BASICS OF PYTHON
>>> annotationString
'Symbol plot from load1.odb'
>>> annotationString.upper()
'SYMBOL PLOT FROM LOAD1.ODB'
>>> annotationString.split()
['Symbol', 'plot', 'from', 'load1.odb']
As with all sequences, you use negative indices to index backward from the end of a string.
Array
Arrays are mutable homogeneous sequences. The numpy module allows you to create and operate
on multidimensional arrays. Python determines the type of elements in the array; you do not have
to declare the type when you create the array. For more information about the numpy module, see
numpy.scipy.org.
4–12
Abaqus ID:
Printed on:
THE BASICS OF PYTHON
Repeat a sequence
>>> x=(1,2)
>>> x*2
(1, 2, 1, 2)
>>> s = 'Hoop Stress'
>>> s*2
>>> 'Hoop StressHoop Stress'
4–13
Abaqus ID:
Printed on:
THE BASICS OF PYTHON
4–14
Abaqus ID:
Printed on:
THE BASICS OF PYTHON
>>> a = [1, 3, 7, 5]
>>> print a.sort()
None
>>> import sys
>>> x = sys.path.append('.')
>>> print x
None
If a statement breaks at any other place, you must include a “\” character at the end of the line to indicate
that it is continued on the next line. For example,
distance = mdb.models['Model-1'].parts['housing'].\
getDistance(entity1=node1, entity2=node2)
When you are running Python from a local UNIX or Windows window, the prompt changes to the “. . .”
characters to indicate that you are on a continuation line.
Comments in a Python script begin with the “#” character and continue to the end of the line.
4–15
Abaqus ID:
Printed on:
THE BASICS OF PYTHON
4–16
Abaqus ID:
Printed on:
THE BASICS OF PYTHON
indentation that governs a block. When your script returns to the original indentation, the block ends.
For example,
max = 5
i = 0
while i <= max:
square = i**2
cube = i**3
print i, square, cube
i = i + 1
print 'Loop completed'
When you are using the Python interpreter from the Abaqus/CAE command line interface or if you
are running Python from a local UNIX or Windows window, the prompt changes to the “. . .” characters
to indicate that you are in a block controlled by indentation.
if, elif, and else
>>> load = 10
>>> if load > 6.75:
... print 'Reached critical load'
... elif load < 2.75:
... print 'Minimal load'
... else:
... print 'Typical load'
while
>>> load = 10
>>> length = 3
>>> while load < 1E4:
... load = load * length
... print load
>>> while 1:
... x = raw_input(Enter a number or 0 to quit:')
... if x == '0':
... break
... else:
... print x
Use continue to skip the rest of the loop and to go to the next iteration.
4–17
Abaqus ID:
Printed on:
PROGRAMMING TECHNIQUES
>>> load = 10
>>> length = -3
>>> while load < 1E6: #Continue jumps up here
... load = load * length
... if load < 0:
... continue #Do not print if negative
... print load
for
Use a sequence to control the start and the end of for loops. The range() function is an easy
way to create a sequence.
>>> for i in range(5):
... print i
...
0
1
2
3
4
The following sections introduce you to some of the techniques you will need to program with Python.
The following topics are covered:
• “Creating functions,” Section 4.6.1
• “Using dictionaries,” Section 4.6.2
• “Reading and writing from files,” Section 4.6.3
• “Error handling,” Section 4.6.4
• “Functions and modules,” Section 4.6.5
• “Writing your own modules,” Section 4.6.6
4–18
Abaqus ID:
Printed on:
PROGRAMMING TECHNIQUES
4–19
Abaqus ID:
Printed on:
PROGRAMMING TECHNIQUES
You can now call the function without providing all the arguments. The function assigns the default
value to any missing arguments.
It is good programming practice to use a documentation string that indicates the purpose of a
function and the arguments expected. A documentation string appears at the top of a function and is
delimited by triple quotes """. You can use the __doc__ method to obtain the documentation string
from a function while running the Python interpreter. For example,
>>>def calculateCylinderVolume(radius=0.5,height=1.0):
... """
... Calculates the volume of a cylinder.
...
... Takes two optional arguments, radius (default=0.5)
... and height (default=1.0).
... """
... from math import *
... volume = pi * radius**2 * height
... return volume
...
>>> print calculateCylinderVolume.__doc__
You can retrieve the documentation string for the methods in the Abaqus Scripting Interface. For
example,
>>> mdb.Model.__doc__
'Mdb.Model(name <, description, stefanBoltzmann, absoluteZero>) ->
This method creates a Model object.'
>>> session.Viewport.__doc__
'Session.Viewport(name <, origin, width, height, border, titleBar,
titleStyle, customTitleString>)
-> This method creates a Viewport object with the specified
origin and dimensions.'
4–20
Abaqus ID:
Printed on:
PROGRAMMING TECHNIQUES
The documentation string shows the name of each argument name and whether the argument is required
or optional. The string also shows a brief description of the method.
You can use the sys module to retrieve command line arguments and pass them to a function. For
example, the following script takes two arguments—the X- and Y-coordinates of a point—and calculates
the distance from the point to the origin. The script uses the following modules:
• The sys module to retrieve the command line arguments.
• The math module to calculate the square root.
import sys, math
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
def distance(x, y):
"""
Prints distance from origin to (x, y).
a = x**2 + y**2
return math.sqrt(a)
x = float(sys.argv[1])
y = float(sys.argv[2])
d = distance(x, y)
4–21
Abaqus ID:
Printed on:
PROGRAMMING TECHNIQUES
4–22
Abaqus ID:
Printed on:
PROGRAMMING TECHNIQUES
>>> myPart.__methods__
['clear', 'copy', 'get', 'has_key', 'items', 'keys',
'update', 'values']
The keys() method returns a list of the dictionary keys.
>>> myPart.keys()
['size', 'weight', 'number', 'material', 'cost', 'color']
The values() method returns a list of the values of each entry in the dictionary.
>>> myPart.values()
[3.0, 376.0, 667, 'Steel', 10.34, 'Red']
The items() method returns a list of tuples. Each tuple contains the key and its value.
>>> myPart.items()
[('size', 3.0), ('number', 667), ('material', 'Steel'),
('color', 'Red'), ('weight', 376.0), ('cost', 10.34),]
You use the has_key() method to see if a key exists. A return value of 1 indicates the key exists in
the dictionary. A return value of 0 indicates the key does not exist.
>>> myPart.has_key('color')
1
Python’s del statement allows you to delete a variable.
>>> del myPart
You can also use del to delete an item from a dictionary.
>>> del myPart['color']
>>> myPart.has_key('color')
0
You can use the keys(), values(), or items() methods to loop through a dictionary. In the
following example, items() returns two values; the first is assigned to property, and the second is
assigned to setting.
>>> for property, setting in myPart.items():
... print property, setting
...
size 3.0
weight 376.0
number 667
material Steel
cost 10.34
4–23
Abaqus ID:
Printed on:
PROGRAMMING TECHNIQUES
The first line opens an existing file in the crash_test directory called fender.txt. The file is
opened in read-only mode; myInputFile is a variable that refers to a file object. The second line
creates and opens a new file object in the local directory called peak_deflection.txt. This file is
opened in read and write mode.
Use the “__methods__” technique that we saw earlier to see the methods of a file object.
The readline() method reads a single line from a file into a string, including the new line character
that terminates the string. The readlines() method reads all the lines in a file into a list. The
write() function writes a string to a file. Look at the standard Python documentation on the official
Python web site (www.python.org) for a description of functions that operate on files. File objects are
described in the Built-in Types section of the Python Library Reference.
The following example reads each line of a text file and changes the line to uppercase characters:
4–24
Abaqus ID:
Printed on:
PROGRAMMING TECHNIQUES
The first line opens the input file; you don’t need the ’r’ because read-only is the default access mode.
The next line opens a new file to which you will write. You read the lines in the input file into a list.
Finally, you enter a loop that converts each line to uppercase characters and writes the result to the output
file. The final two lines close the files.
if y == 0:
raise ValueError, 'y argument cannot be zero'
else:
4–25
Abaqus ID:
Printed on:
PROGRAMMING TECHNIQUES
return x/y
try:
print myFunction(temperature, velocity)
except ValueError, error:
print error
Exception handling is discussed in more detail in “Error handling in the Abaqus Scripting Interface,”
Section 5.5.
4–26
Abaqus ID:
Printed on:
PROGRAMMING TECHNIQUES
>>> sin(x)
0.707106781187
The first line imports all of the names from the math module. The second line uses pi, a float number
defined in the math module. The third line refers to a sin() function. Python can use the sin()
function because you imported it from the math module.
To import only the sin() function, you could have typed
>>> from math import sin
You need to import a module only once during a session. Once a module is imported, its functions,
methods, and attributes are always available to you. You cannot “unload” a module after you import it.
To see a list of all the functions that come with the math module, look at the Miscellaneous
Services section of the Python Library Reference. You can download public-domain
modules, and you can create your own modules.
Python provides a second approach to importing modules. For example,
>>> import math
>>> x = 22.0/(7.0 * 4.0)
>>> math.sin(x)
0.707330278085
The “import” approach shown above imports the module as a unit, and you must qualify the name of an
object from the module. To access a function from the math module in our example, you must prepend
the function with math.; the math. statement is said to “qualify” the sin() function.
What is the difference between the two approaches to importing modules? If two modules contain
an object with the same name, Python cannot distinguish between the objects if you use the “from
modulename import *” approach. If two objects have the same name, Python uses the object most
recently imported. However, if you use the “import modulename” approach, modulename qualifies the
name of the object and makes it unique.
import math
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
def distance(x, y):
4–27
Abaqus ID:
Printed on:
PROGRAMMING TECHNIQUES
"""
Prints distance from origin to (x, y).
a = x**2 + y**2
return math.sqrt(a)
You must import the module to make use of the functions and constants that it contains.
import myUtilities
4–28
Abaqus ID:
Printed on:
FURTHER READING
reload(myModule)
maxStress = myModule.calculateStress(odb)
4–29
Abaqus ID:
Printed on:
EXECUTING SCRIPTS
If your script accesses and manipulates data in an output database, you can execute the script using
either of the methods already described:
• Type abaqus python scriptname .py at the system prompt. The script must contain the
following statement:
5–1
Abaqus ID:
Printed on:
Abaqus SCRIPTING INTERFACE DOCUMENTATION STYLE
• Select File→Run Script from the Abaqus/CAE main menu bar, and select the file to execute. The
script must contain the following statement:
from visualization import *
This section describes the style that is used to describe a command in the Abaqus Scripting Reference
Manual. You may want to refer to the Abaqus Scripting Reference Manual while you read this section
and compare the style of a documented command with the descriptions provided here. The following
topics are covered:
• “How the commands are ordered,” Section 5.2.1
• “Access,” Section 5.2.2
• “Path,” Section 5.2.3
• “Arguments,” Section 5.2.4
• “Return value,” Section 5.2.5
5–2
Abaqus ID:
Printed on:
Abaqus SCRIPTING INTERFACE DOCUMENTATION STYLE
• Within each object description, the commands are listed in the following order:
– Constructors (in alphabetical order)
– Methods (in alphabetical order)
– Members
• Some methods are not associated with an object and appear at the end of a chapter; for example,
the evaluateMaterial() method appears at the end of Chapter 28, “Material commands,” of
the Abaqus Scripting Reference Manual.
5.2.2 Access
The description of each object in the Abaqus Scripting Reference Manual begins with a section that
describes how you access an instance of the object. The import statements are provided for completeness.
Abaqus/CAE imports all modules when you start a session, and you do not need to include the import
module name statement in your scripts. However, you must import the Abaqus Scripting Interface
Symbolic Constants with the following statement:
These should be the first statement in all your Abaqus Scripting Interface scripts.
The following is the access description for the Material object:
import material
mdb.models[name].materials[name]
The first line of the access description indicates the module that Abaqus/CAE imported to make this
object, and its methods and members, available to your script.
The access description also specifies where instances of the object are located in the data model.
In the previous example the second line indicates how your script can access Material objects from a
particular model. You must qualify a material object, command, or member with the variable mdb, as
described in “Functions and modules,” Section 4.6.5. For example,
mdb.models[crash].Material[steel]
mdb.models[crash].materials[steel].Elastic(
table=((30000000.0, 0.3), ))
elasticityType = mdb.models[crash].materials[steel].elastic.type
Similarly, if you are reading from an output database, the following is the access description for the
HistoryRegion object:
import odbAccess
session.odbs[name].steps[name].historyRegions[name]
5–3
Abaqus ID:
Printed on:
Abaqus SCRIPTING INTERFACE DOCUMENTATION STYLE
The first line indicates that Abaqus/CAE imported the odbAccess module to make the Odb objects,
methods, and members available to your Abaqus Scripting Interface script. The second line indicates
how your script can access HistoryRegion objects from a particular step.
The Access description for the FieldOutput object is
session.odbs[name].steps[name].frames[i].fieldOutputs[name]
The following statements show how you use the object described by this Access description:
sideLoadStep = session.odbs['Forming loads'].steps['Side load']
lastFrame = sideLoadStep.frames[-1]
stressData = lastFrame.fieldOutputs['S']
integrationPointData = stressData.getSubset(
position=INTEGRATION_POINT)
invariantsData = stressData.validInvariants
• The next to last line shows the getSubset method of the FieldOutput object.
• The last line shows the validInvariants member of the FieldOutput object.
5.2.3 Path
A method that creates an object is called a “constructor.” The Abaqus Scripting Interface uses the
convention that constructors begin with an uppercase character. In contrast, methods that operate on an
object begin with a lowercase character. The description of each constructor in the Abaqus Scripting
Reference Manual includes a path to the command. For example, the following describes the path to
the Viewport constructor:
session.Viewport
Some constructors include more than one path. For example, you can create a datum that is associated
with either a Part object or the RootAssembly object, and each path is listed.
mdb.models[name].parts[name].DatumAxisByCylFace
mdb.models[name].rootAssembly.DatumAxisByCylFace
The path is not listed if the method is not a constructor.
If you are using the Abaqus Scripting Interface to read data from an output database, the objects exist
when you open the output database, and you do not have to use constructors to create them. However, if
you are creating or writing to an output database, you may need to use constructors to create new objects,
such as part instances and steps. The documentation describes the path to the constructors that create
objects in an output database.
For example, the Path description for the FieldOutput constructor is
session.odbs[name].steps[name].frames[i].FieldOutput
5–4
Abaqus ID:
Printed on:
Abaqus SCRIPTING INTERFACE DOCUMENTATION STYLE
5.2.4 Arguments
The ellipsis (...) in the command description indicates that the command takes one or more
arguments. For example, the Viewport constructor takes arguments.
Viewport(...)
In contrast, the makeCurrent method takes no arguments.
makeCurrent()
Some arguments of a command are required, and some arguments are optional. In the Abaqus
Scripting Reference Manual the required arguments are listed first, followed by the optional arguments.
If the argument is optional, the default value is provided. The default value is the value of an optional
argument when you call a method and omit the argument.
The setValues method is a special case. All of the arguments to the setValues method are
optional, but any argument that you omit retains its current value; Abaqus does not assign a default value
to the argument.
Some objects have no constructors; Abaqus creates the objects for you. For such objects the
documentation describes the “initial value” of an optional argument. The initial value given for the
argument is the initial value assigned to the corresponding member when Abaqus creates the object.
For example, the defaultViewportAnnotationOptions object has no constructor; Abaqus creates the
defaultViewportAnnotationOptions object when you start a session. When you create a new viewport,
the settings are copied from the current viewport.
You can use the setValues method to modify the value of a member; for example, to modify
the value of the triad member of the defaultViewportAnnotationsOptions object. When you call
session.defaultViewportAnnotationOptions.setValues(triad=OFF), the value
of the triad member is set to off. The other member values remain unchanged; this behavior is called
“as is” behavior because the values remain “as is.” The setValuesInStep method displays similar
“as is” behavior.
Keyword and positional arguments are described in “Creating functions,” Section 4.6.1. We
recommend that you use keyword arguments since they can be supplied in any order and they make
your scripts easier to read and debug; for example,
newViewport = session.Viewport(name='myViewport',
origin=(10, 10), width=100, height=50)
5–5
Abaqus ID:
Printed on:
Abaqus SCRIPTING INTERFACE DOCUMENTATION STYLE
If you choose not to use keywords, the arguments must be provided in the order in which they are
documented.
newViewport = session.Viewport('myViewport',
(10, 10), 100, 50)
You can use a combination of keyword and positional arguments. Keyword arguments can be supplied
after positional arguments; however, positional arguments cannot be entered after keyword arguments.
For example, you can use the following statement:
newViewport = session.Viewport('myViewport',
(10, 10), width=100, height=50)
newViewport = session.Viewport(name='myViewport',
(10, 10), 100, 50)
You will find it easier to use keyword arguments so that you do not have to concern yourself with the
positional requirements.
newViewport = session.Viewport(name='myViewport',
origin=(10, 10), width=100, height=50)
You can use the object returned by a command in subsequent statements. For example, the titlebar
member of a Viewport object is a Boolean specifying whether the viewport title bar is displayed and can
have a value of either ON or OFF. The following statement tests the titlebar member of the new viewport
created by the previous statement:
if newViewport.titleBar:
print 'The title bar will be displayed.'
5–6
Abaqus ID:
Printed on:
Abaqus SCRIPTING INTERFACE DATA TYPES
The standard Python data types described in “Python data types,” Section 4.5.2, include integers, floats,
strings, and sequences. The Abaqus Scripting Interface adds over 500 additional data types. The
following sections describe the most common Abaqus Scripting Interface data types:
• “SymbolicConstants,” Section 5.3.1
• “Booleans,” Section 5.3.2
• “Repositories,” Section 5.3.3
5.3.1 SymbolicConstants
Some arguments require that you provide a SymbolicConstant. SymbolicConstants are defined by the
Abaqus Scripting Interface and are written in all capital letters. If your script uses a SymbolicConstant
defined by the Abaqus Scripting Interface, you must import the SymbolicConstant with the following
statement before you refer to it:
from abaqusConstants import *
When an argument to a command is a SymbolicConstant, the description in the Abaqus Scripting
Reference Manual lists all its possible values. For example, when you are printing an image, the image
can be rendered in one of the following formats: BLACK_AND_WHITE, GREYSCALE, or COLOR.
Similarly, a data member can be a SymbolicConstant. For example, the type member
of the Elastic object can be one of the following SymbolicConstants: ISOTROPIC,
ORTHOTROPIC, ANISOTROPIC, ENGINEERING_CONSTANTS, LAMINA, TRACTION, or
COUPLED_TRACTION.
If the SymbolicConstants provided by the Abaqus Scripting Interface do not meet your needs,
you can create your own SymbolicConstants using the SymbolicConstant constructor. For more
information, see “SymbolicConstant object,” Section 52.1 of the Abaqus Scripting Reference Manual.
5.3.2 Booleans
Python defines two Boolean values, True and False. The type of a Python Boolean is <type ’bool’>.
myPythonBoolean = True
type(myPythonBoolean)
<type 'bool'>
In addition, the Abaqus Scripting Interface defines a Boolean object, derived from the SymbolicConstant
object, which can take the values ON and OFF. For example, noPartsInputFile is a member of a Model
5–7
Abaqus ID:
Printed on:
Abaqus SCRIPTING INTERFACE DATA TYPES
object that indicates whether the input file will be written with parts and assemblies. The type of the
noPartsInputFile member is <type ’AbaqusBoolean’>.
Abaqus recommends that you use the Python Boolean in your scripts and that you convert existing
scripts to use the Python Boolean.
The value of a Boolean argument can appear to be ambiguous; for example,
newModel = mdb.ModelFromInputFile(name='beamTutorial',
inputFileName='Deform')
newModel.setValues(noPartsInputFile=False)
print newModel.noPartsInputFile
OFF
Because of this ambiguity, you should test a Boolean for a positive or negative value, as opposed to
comparing it to a specific value like 0, OFF, or False. For example, the following statements show how
you should test the value of a Boolean member:
if (newModel.noPartsInputFile):
print 'Input file will be written without parts \
and assemblies. '
else:
print 'Input file will be written with parts \
and assemblies.'
5.3.3 Repositories
Repositories are containers that store a particular type of object; for example, the steps repository
contains all the steps defined in the model. A repository maps to a set of information and is similar to
a Python dictionary; for more information, see “Using dictionaries,” Section 4.6.2. However, only a
constructor can add an object to a repository. In addition, all the objects in a repository are of the same
type. For example, the following repository contains all the models in the model database:
mdb.models
In turn, the following repository contains all the parts in the model Model-1:
mdb.models['Model-1'].parts
As with dictionaries, you can refer to an object in a repository using its key. The key is typically
the name you provided in the constructor command when the object was created. For example, the
Viewport constructor creates a new Viewport object in the viewports repository.
session.Viewport(name='Side view',
origin = (10,10), width=50, height=50)
5–8
Abaqus ID:
Printed on:
OBJECT-ORIENTED PROGRAMMING AND THE Abaqus SCRIPTING INTERFACE
The key to this new Viewport object in the viewports repository is Side view. You use this key
to access this particular Viewport object. For example,
session.viewports['Side view'].viewportAnnotationOptions.\
setValues(legend=OFF, title=OFF)
You can make your scripts more readable by assigning a variable to an object in a repository. For
example, you could rewrite the previous statement after assigning the Viewport object to the variable
myViewport:
In general, if the user can create the object, its repository key is a string. In some cases Abaqus/CAE
creates an object, and the key can be a string, an integer, or a SymbolicConstant.
As with dictionaries, you can use the keys() method to access the repository keys.
>>> session.Viewport(name='Side view')
>>> session.Viewport(name='Top view')
>>> session.Viewport(name='Front view')
>>> for key in session.viewports.keys():
...
print key
Front view
Top view
Side view
You can use the keys()[ i] method to access an individual key; however, most repositories are not
ordered, and this is not recommended.
You can use the changeKey() method to change the name of a key in a repository. For example,
myPart = mdb.models['Model-1'].Part(name='housing',
dimensionality=THREE_D, type=DEFORMABLE_BODY)
mdb.models['Model-1'].parts.changeKey(fromName='housing',
toName='form')
You should now be familiar with some of the concepts behind object-oriented programming, such as
objects, constructors, methods, and members. This section describes how object-oriented programming
5–9
Abaqus ID:
Printed on:
OBJECT-ORIENTED PROGRAMMING AND THE Abaqus SCRIPTING INTERFACE
relates to the Abaqus Scripting Interface and summarizes some of the terminology. The following topics
are covered:
• “The Abaqus Scripting Interface and methods,” Section 5.4.1
• “The Abaqus Scripting Interface and members,” Section 5.4.2
• “Object-oriented programming and the Abaqus Scripting Interface—a summary,” Section 5.4.3
5–10
Abaqus ID:
Printed on:
OBJECT-ORIENTED PROGRAMMING AND THE Abaqus SCRIPTING INTERFACE
The standard Python statement object.__members__ lists all the members of an object. For example,
the following statement lists all the members of a Viewport object:
session.viewports['myViewport'].__members__
The values of members are specific to each instance of the object. For example, the value of the width
member of a Viewport object is specific to each viewport.
Members of an Abaqus object are read-only; consequently, you cannot change their value with a
simple assignment statement. You use the setValues() method to change the value of a member. For
example, the setValues() statement in the following script changes the thickness of a shell section:
>>> import section
>>> shellSection = mdb.models['Model-1'].HomogeneousShellSection(
name='Steel Shell', thickness=1.0, material='Steel')
>>> print 'Original shell section thickness = ' \
, shellSection.thickness
Original shell section thickness = 1.0
>>> shellSection.setValues(thickness=2.0)
>>> print 'Final shell section thickness = ' \
, shellSection.thickness
Final shell section thickness = 2.0
You cannot use assignment to change the value of the Shell object.
>>> myShell.thickness = 2.0
TypeError: readonly Attribute
The following statements illustrate the use of constructors, methods, and members:
>>> # Create a Section object
>>> mySection = mdb.models['Model-1'].HomogeneousSolidSection(
name='solidSteel', material='Steel', thickness=1.0)
>>> # Display the type of the object
>>> print 'Section type = ', type(mySection)
Section type = <type 'HomogeneousSolidSection'>
>>> # List the members of the object
>>> print 'Members of the section are:' , mySection.__members__
Members of the section are: ['category', 'dimension',
'layout', 'material', 'name',
'thickness']
>>> # List the methods of the object
>>> print 'Methods of the section are: ', mySection.__methods__
Methods of the section are: ['setValues']
>>> # Print the value of each member in a nice format
>>> for member in mySection.__members__:
5–11
Abaqus ID:
Printed on:
OBJECT-ORIENTED PROGRAMMING AND THE Abaqus SCRIPTING INTERFACE
...
print 'mySection.%s = %s' % (member,
getattr(mySection, member))
mySection.category = SOLID
mySection.dimension = THREE_DIM
mySection.layout = HOMOGENEOUS
mySection.material = Steel
mySection.name = solidSteel
mySection.thickness = 1.0
You use the Access description provided with each object in the Abaqus Scripting Reference
Manual to determine how you access the object. You append a method or member to this description
when you are writing a script. Similarly, you use the Path description provided with each constructor
in the Abaqus Scripting Reference Manual to determine the path to the constructor.
5–12
Abaqus ID:
Printed on:
ERROR HANDLING IN THE Abaqus SCRIPTING INTERFACE
• You use the Path description provided with each constructor in the Abaqus Scripting Reference
Manual to determine the path to the constructor.
• You use the setValues() method to modify the members of an Abaqus Scripting Interface
object.
session.viewports['Side view'].setValues(origin=(20,20))
The basics of Python’s exception handling are described in “Error handling,” Section 4.6.4; and the
same techniques apply to the Abaqus Scripting Interface. If certain circumstances arise while a script
is running, Python allows you to take the necessary action and still allows the script to continue.
Alternatively, when Abaqus/CAE issues (or “throws”) an exception and the exception is not handled
by the script, Abaqus/CAE displays the exception message in the message area and the script stops
executing.
The following topics are covered:
• “Standard Python exceptions,” Section 5.5.1
• “Standard Abaqus Scripting Interface exceptions,” Section 5.5.2
• “Additional Abaqus Scripting Interface exceptions,” Section 5.5.3
• “Exception handling,” Section 5.5.4
5–13
Abaqus ID:
Printed on:
ERROR HANDLING IN THE Abaqus SCRIPTING INTERFACE
InvalidNameError
You specified an invalid name. Abaqus/CAE enforces a naming convention for objects that you
create. Names must adhere to the following rules:
• The name can have up to 38 characters.
• The name can include spaces and most punctuation marks and special characters; however,
only 7-bit ASCII characters are supported.
WARNING: While Python allows most punctuation marks and special characters, some of the
strings you provide will be used in an Abaqus input file; therefore, you cannot use the following
characters: $&*~!()[]{}|;’‘",.?/\ when you are naming an object, such as a part, a model or
a job.
RangeError
A numeric value is out of range.
AbaqusError
Context-dependent message.
AbaqusException
Context-dependent message.
Note: The command descriptions in the Abaqus Scripting Reference Manual list the type of standard
Abaqus Scripting Interface exceptions that can occur; however, the exception messages are not included
with the command description.
5–14
Abaqus ID:
Printed on:
ERROR HANDLING IN THE Abaqus SCRIPTING INTERFACE
Standard Abaqus
Scripting Interface exception
Additional
exception
Exception message
Exception type
Figure 5–1 The layout of a typical exception description in the online documentation.
try:
session.Viewport(name='tiny',width=1, height=1)
except RangeError, message:
print 'Viewport too small:', message
print 'Script continues running and prints this line'
5–15
Abaqus ID:
Printed on:
EXTENDING THE Abaqus SCRIPTING INTERFACE
You can extend the functionality of the Abaqus Scripting Interface by writing your own modules that
contain classes and functions to accomplish tasks that are not directly available in Abaqus. For example,
you can write a function to print the names of all materials that have a density specified, or you can write
a function that creates a contour plot using a custom set of contour plot options. Creating functions and
modules in Python is described in “Creating functions,” Section 4.6.1, and “Functions and modules,”
Section 4.6.5.
This section describes how you can extend the functionality of the Abaqus Scripting Interface. The
following topics are covered:
• “Storing custom data in the model database or in other objects,” Section 5.6.1
• “Interaction with the GUI,” Section 5.6.2
• “CommandRegister class,” Section 5.6.3
• “Repositories,” Section 5.6.4
• “Repository methods,” Section 5.6.5
• “RepositorySupport,” Section 5.6.6
• “Registered dictionaries,” Section 5.6.7
• “Registered lists,” Section 5.6.8
• “Registered tuples,” Section 5.6.9
• “Session data,” Section 5.6.10
• “Saving application data in a model database,” Section 5.6.11
• “Checking a model database when it is opened,” Section 5.6.12
5–16
Abaqus ID:
Printed on:
EXTENDING THE Abaqus SCRIPTING INTERFACE
mdb.saveAs('custom-test.cae')
mdb.close()
If you start a new session and open the model database, custom-test.cae, you can refer to the
variables that you saved. For example,
>>> import customKernel
mdb = openMdb('custom-test.cae')
>>> print mdb.customData.myString, mdb.customData.myNumber
The width is 58
You can store almost any type of Python object under mdb.customData; for example, strings,
numbers, and Python classes. However, there are some restrictions; for example, you cannot store file
objects. These restrictions are due to the fact that the Abaqus/CAE infrastructure uses Python’s pickle
module to store the customData object in the model database. The pickle module allows the Python
programmer to write a data structure to a file and then recreate that data structure when reading from
the file. For details on the restrictions imposed by the pickle module, see the official Python web site
(www.python.org).
If your code creates a custom class and stores an instance of the class in the model database, the
custom module that defined that custom class must be available for Python to unpickle the data when
the database is subsequently opened. Consequently, if a user saves custom data to a model database
and then passes that model database to another user, the other user must also have access to the custom
modules that produced the custom data. Otherwise, they will not be able to load the custom data into
their Abaqus/CAE session.
Abaqus/CAE does not keep track of changes made to the customData object. As a result, when
the user quits a session, Abaqus/CAE will not prompt them to save their changes if they changed only
objects under customData.
• Notification to the GUI when custom kernel data change. For example, you can have a manager
dialog box that lists the objects in a repository. When the contents of the repository change, you can
be notified and take the appropriate action to update the list of objects in the manager dialog box.
To make use of these features, you must derive your custom kernel objects from the classes listed in
the following sections. For more details on GUI customization, see the Abaqus GUI Toolkit Reference
Manual.
5–17
Abaqus ID:
Printed on:
EXTENDING THE Abaqus SCRIPTING INTERFACE
5.6.4 Repositories
Repositories are containers that hold objects that are keyed by strings. It may be convenient to store your
custom kernel objects in repositories, in the same way that Abaqus/CAE part objects are stored in the
Parts repository.
The customData object is an instance of a RepositorySupport class, which provides a
Repository method that allows you to create a repository as an attribute of the instance. For more
information, see “RepositorySupport,” Section 5.6.6. The arguments to the Repository method are
the name of the repository and a constructor or a sequence of constructors. Those constructors must
have name as their first argument, and the infrastructure will automatically assign that value to a member
called name. Instances of these constructors will be stored in the repository. For more information, see
“Repository object,” Section 52.3 of the Abaqus Scripting Reference Manual.
Since repositories are designed to notify the GUI when their contents change, the objects placed
inside them should be derived from either CommandRegister or RepositorySupport to extend
this capability to its fullest.
The Abaqus Scripting Interface uses the following conventions:
• The name of a repository is a plural noun with all lowercase letters.
• A constructor is a capitalized noun (or a combination of capitalized nouns and adjectives).
• The first argument to the constructor must be name.
For example, the Part constructor creates a part object and stores it in the parts repository.
You can access the part object from the repository using the same name argument that you passed
in with the Part constructor. In some cases, more than one constructor can create instances that
are stored in the same repository. For example, the HomogeneousSolidSection and the
5–18
Abaqus ID:
Printed on:
EXTENDING THE Abaqus SCRIPTING INTERFACE
HomogeneousShellSection constructors both create section objects that are stored in the
sections repository. For more information, see “Abstract base type,” Section 6.1.5. For example,
the following script creates a blocks repository, and the Block constructor creates a block object
in the blocks repository:
from customKernel import CommandRegister
class Block(CommandRegister):
def __init__(self, name):
CommandRegister.__init__(self)
mdb.customData.Repository('blocks', Block)
block = mdb.customData.Block(name='Block-1')
print mdb.customData.blocks['Block-1'].name Block-1
mdb.customData.Repository('blocks', Block)
mdb.customData.Block(name='Block-1')
mdb.customData.Block(name='Block-2')
print 'The original repository keys are: ',
mdb.customData.blocks.keys()
5–19
Abaqus ID:
Printed on:
EXTENDING THE Abaqus SCRIPTING INTERFACE
print mdb.customData.blocks.has_key('Block-2')
print mdb.customData.blocks.has_key('Block-3')
mdb.customData.blocks.changeKey('Block-1', 'Block-11')
print 'The modified repository keys are: ',
mdb.customData.blocks.keys()
print 'The name member is ',
mdb.customData.blocks['Block-11'].name
print 'The repository size is', len(mdb.customData.blocks)
The resulting output is
The original repository keys are ['Block-1', 'Block-2']
1
0
The modified repository keys are ['Block-11', 'Block-2']
The name member is Block-11
The repository size is 2
5.6.6 RepositorySupport
You can use the RepositorySupport class to derive a class that can contain one or more repositories.
However, if you do not intend to create a repository as an attribute of your class, you should derive your
class from CommandRegister, not from RepositorySupport.
Using the RepositorySupport class allows you to create a hierarchy of repositories; for
example, in the Abaqus Scripting Interface the parts repository is a child of the models repository.
The first argument passed into your constructor is stored as name; it is created automatically by the
infrastructure. To create a hierarchy of repositories, derive your class from RepositorySupport
and use its Repository method to create child repositories as shown below. The Repository
method is described in “Repositories,” Section 5.6.4.
from abaqus import *
from customKernel import CommandRegister, RepositorySupport
class Block(CommandRegister):
def __init__(self, name):
CommandRegister.__init__(self)
class Model(RepositorySupport):
def __init__(self, name):
RepositorySupport.__init__(self)
self.Repository('blocks', Block)
mdb.customData.Repository('models', Model)
5–20
Abaqus ID:
Printed on:
EXTENDING THE Abaqus SCRIPTING INTERFACE
mdb.customData.Model('Model-1')
mdb.customData.models['Model-1'].Block('Block-1')
5–21
Abaqus ID:
Printed on:
EXTENDING THE Abaqus SCRIPTING INTERFACE
mdb.customData.Repository('blocks', Block)
block1 = mdb.customData.Block(name='Block-1')
block2 = mdb.customData.Block(name='Block-2')
tuple = (block1, block2)
mdb.customData.myTuple = RegisteredTuple(tuple)
5–22
Abaqus ID:
Printed on:
EXTENDING THE Abaqus SCRIPTING INTERFACE
You can call the setAppData method only once per application name, which prevents unauthorized
changes to the method. However, the setAppData method may be called multiple times using different
application names to allow more than one application to register with the same model database.
When the user saves a model database, Abaqus copies the session.customData.appData object to
the mdb.customData.appData object.
Kernel initialization scripts specified by the startup command line option are executed by
Abaqus/CAE after it has finished its initialization process. By that time, a new model database or
a database specified on the command line using the database option has already been opened. A
utility method called processInitialMdb has been created to automatically process the initial
5–23
Abaqus ID:
Printed on:
EXTENDING THE Abaqus SCRIPTING INTERFACE
model database for you. If the initial model database does not have any customData or does not have
customData for your particular application, your initializeMdb method will be called. If the initial
model database has customData for your application, your verifyMdb method will be called.
The following example shows how you can use the verifyMdb, intializeMdb, and
processInitialMdb methods. You should execute the example using the startup command line
option when you start Abaqus/CAE. For more information, see “Abaqus/CAE execution,” Section 3.2.4
of the Abaqus Analysis User’s Manual.
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
def initializeMdb():
# Initialize the MDB here
5–24
Abaqus ID:
Printed on:
EXTENDING THE Abaqus SCRIPTING INTERFACE
5–25
Abaqus ID:
Printed on:
THE Abaqus OBJECT MODEL
We have already discussed how Python provides built-in objects like integers, lists, dictionaries, and
strings. When you are writing Abaqus Scripting Interface scripts, you need to access these built-in
objects together with the objects used by Abaqus/CAE. These Abaqus Scripting Interface objects extend
Python with new types of objects. The hierarchy and the relationship between these objects is called
the Abaqus object model. The following sections describe the Abaqus object model in more detail. The
following topics are covered:
• “An overview of the Abaqus object model,” Section 6.1.1
• “Using tab completion to explore the object model,” Section 6.1.2
• “The Model object model,” Section 6.1.3
• “Using the object model,” Section 6.1.4
• “Abstract base type,” Section 6.1.5
• “Importing modules to extend the object model,” Section 6.1.6
6–1
Abaqus ID:
Printed on:
THE Abaqus OBJECT MODEL
• Definitions of the relationships between the objects. These relationships form the structure or
the hierarchy of the object model. The relationships between the objects are summarized in the
following list:
Ownership
The ownership hierarchy defines the access path to the objects in the Abaqus model.
Associations
Associations describe the relationships between the objects; for example, whether one object
refers to another and whether an object is an instance of another.
Abaqus extends Python with approximately 500 additional objects, and there are many relationships
between these objects. As a result, the complete Abaqus object model is too complex to illustrate in a
single figure.
In general terms the Abaqus object model is divided into the Session, the Mdb, and the Odb objects,
as shown in Figure 6–1.
session mdb
viewports
jobs
fieldReportOptions
models
Model
sketches
parts
odb
rootAssembly = Container
parts = Singular object
sectionCategories
steps
6–2
Abaqus ID:
Printed on:
THE Abaqus OBJECT MODEL
Singular object
Objects that are not containers are shown as a Singular object. A singular object contains no
other objects of a similar type; for example, the Session object and the Mdb object. There is only
one Session object and only one Mdb object in the Abaqus object model.
The “...” at the end of the object models shown in this section indicates that there are additional objects
in the model that are not included in the figure. For clarity, the figures show only the most commonly
used objects in the object model.
The statement from abaqus import * imports the Session object (named session) and the
Mdb object (named mdb) and makes them available to your scripts. The statement from odbAccess
import * allows you to access Abaqus output results from your script. The Session, Mdb, and Odb
objects are described as follows:
Session
Session objects are objects that are not saved between Abaqus/CAE sessions; for example, the
objects that define viewports, remote queues, and user-defined views, as shown in Figure 6–2. The
viewports container is owned by the Session object, as shown in Figure 6–3.
Mdb
The statement from abaqus import * creates an instance of the Mdb object called mdb. Mdb
objects are objects that are saved in a model database and can be recovered between Abaqus/CAE
sessions. Mdb objects include the Model object and the Job object. The Model object, in turn, is
comprised of Part objects, Section objects, Material objects, Step objects, etc. Figure 6–4 shows
the basic structure of the objects under the Model object. For more information, see “The Model
object model,” Section 6.1.3.
Odb
Odb objects are saved in an output database and contain both model and results data, as shown in
Figure 6–5.
Most of the commands in the Abaqus Scripting Interface begin with either the Session, the Mdb, or the
Odb object. For example,
session.viewports['Viewport-1'].bringToFront()
mdb.models['wheel'].rootAssembly.regenerate()
stress = odb.steps['Step-1'].frames[3].fieldOutputs['S']
6–3
Abaqus ID:
Printed on:
THE Abaqus OBJECT MODEL
session
odbs
defaultOdbDisplay
displayGroups
colors
printOptions
psOptions
epsOptions
= Container pngOptions
= Singular object tiffOptions
aviOptions
animationController
views
viewports
defaultViewportAnnotationOptions
sketcherOptions
fieldReportOptions
probeOptions
probeReport
selectedProbeValues
paths
xyDataObjects
xyPlots
xyReportOptions
queues
6–4
Abaqus ID:
Printed on:
THE Abaqus OBJECT MODEL
session
viewports
Viewport
assemblyDisplay
displayedObject
displayMode
= Container
odbDisplay
= Singular object
partDisplay
view
viewportAnnotationOptions
= Container materials
= Singular object parts
profiles
rootAssembly
sections
sketches
steps
Figure 6–4 The structure of the objects under the Model object.
line interface. When you press the [Tab] key, the command line cycles through the parts in the model.
When you press [Shift]+[Tab], the command line cycles backwards through the parts in the model.
6–5
Abaqus ID:
Printed on:
THE Abaqus OBJECT MODEL
odb
rootAssembly
parts
sectionCategories
steps
= Container Step
= Singular object frames
Frame
historyRegions fieldOutputs
HistoryRegion
historyPoints
historyOutputs
Tab completion also searches the file system when it detects an incomplete string. For example,
from part import THR[Tab]
from part import THREE_D
openMdb('hinge_t[Tab]
openMdb('hinge_tutorial.mdb')
from odbAccess import *
myOdb=openOdb('vi[Tab]
myOdb=openOdb('viewer_tutorial.odb')
In most cases when you type in a constructor or a method and include the opening parenthesis, tab
completion prompts you to provide a value for a keyword argument. For example,
mdb.models['Model-1'].Part([Tab]
mdb.models['Model-1'].Part(name=
When you press the [Tab] key, the command line cycles through the arguments to the method.
You can use tab completion when you are accessing an output database. For example,
p=myOdb.parts[[Tab]
p=myOdb.parts['Part-1']
You can also use tab completion when you are accessing an output database from the Abaqus Python
prompt. For example,
6–6
Abaqus ID:
Printed on:
THE Abaqus OBJECT MODEL
abaqus python
>>>from odbAccess import *
>>>myOdb=openOdb('viewer_tutorial.odb')
>>>p=myOdb.parts[[Tab]
>>>p=myOdb.parts['Part-1']
models
parts
Part
features
datums
cells
faces
edges
vertices
= Container elements
= Singular object
nodes
referencePoints
sets
reinforcements
The Job object is separate from the Model object. The object model for the Job object is
straightforward; the Job object owns no other objects. The Job object refers to a Model object but is
not owned by the Model object.
6–7
Abaqus ID:
Printed on:
THE Abaqus OBJECT MODEL
models
rootAssembly
features
datums
sets
= Container
surfaces
= Singular object
instances
PartInstance
datums
cells
faces
edges
vertices
elements
nodes
referencePoints
sets
• The object model describes the relationships between objects. For example, in object-oriented
programming terms a geometry object, such as a Cell, Face, Edge, or Vertex object, is said to be
“owned” by the Part object. The Part object, in turn, is owned by the Model object. This ownership
relationship between objects is referred to as the ownership hierarchy of the object model.
Ownership implies that if an object is copied, everything owned by that object is also copied.
Similarly, if an object is deleted, everything owned by the object is deleted. This concept is similar
to parent-child relationships in Abaqus/CAE. If you delete a Part, all the children of the part—such
as geometry, datums, and regions—are also deleted.
• The relationships between objects are described in the Path and Access descriptions in the
command reference. For example, the following statement uses the path to a Cell object:
cell4 = mdb.models['block'].parts['crankcase'].cells[4]
The statement mirrors the structure of the object model. The Cell object is owned by a Part object,
the Part object is owned by a Model object, and the Model object is owned by the Mdb object.
• The associations between the objects are captured by the object model. Objects can refer to other
objects; for example, the section objects refer to a material, and the interaction objects refer to a
6–8
Abaqus ID:
Printed on:
THE Abaqus OBJECT MODEL
region, to steps, and possibly to amplitudes. An object that refers to another object usually has a
data member that indicates the name of the object to which it is referring. For example, material is
a member of the section objects, and createStepName is a member of the interaction objects.
Step
InitialStep AnalysisStep
StaticStep BuckleStep
In contrast the object model figures described at the beginning of this section show what object-oriented
programmers call “has a” relationships between objects. For example, a session has a viewport
repository, and a model has a root assembly.
Abaqus uses the name of the abstract base type as the name of the repository that contains
objects of similar types. For example, the StaticStep, BuckleStep, and FrequencyStep
constructors all create objects in the steps repository. Other abstract base types include Amplitude,
BoundaryCondition, Datum, Field, Interaction, and Section.
The term “abstract” implies that the Abaqus object model does not contain an object that has the
type of an abstract base type. For example, there are no objects of type Load or Step in the Abaqus
object model. In contrast, the Feature object is a base type, but it is not abstract. The Abaqus object
model includes Feature objects.
6–9
Abaqus ID:
Printed on:
THE Abaqus OBJECT MODEL
6–10
Abaqus ID:
Printed on:
COPYING AND DELETING Abaqus SCRIPTING INTERFACE OBJECTS
The following section describes how you copy and delete Abaqus Scripting Interface objects. The
following topics are covered:
• “Creating a copy of an object,” Section 6.2.1
• “More on copying objects,” Section 6.2.2
• “Deleting objects,” Section 6.2.3
6–11
Abaqus ID:
Printed on:
COPYING AND DELETING Abaqus SCRIPTING INTERFACE OBJECTS
import material
import section
impactModel = mdb.Model(name='Model A')
mySteel = impactModel.Material(name='Steel')
# Create a section
firstSection = impactModel.HomogeneousSolidSection(
name='steelSection 1', material='Steel',
thickness=1.0)
secondSection = impactModel.Section(
name='steelSection 2', objectToCopy=firstSection)
myMaterial = mdb.models['Model-1'].Material(name='aluminum')
You can use the Python del statement to delete an object, but you must provide the full path to the
object.
del mdb.models['Model-1'].materials['aluminum']
The variable myMaterial that referred to the object still exists; however, the variable no longer refers
to the object. You can use the del statement to delete the variable.
del myMaterial
Conversely, if you create the object as before but delete the variable that referred to the object, only
the variable is deleted; the object still exists. You can assign a new variable to the object.
6–12
Abaqus ID:
Printed on:
Abaqus/CAE SEQUENCES
myMaterial = mdb.models['Model-1'].Material(name='aluminum')
del myMaterial
myNewMaterial = mdb.models['Model-1'].materials['aluminum']
The previous explanation does not apply to the few Abaqus/CAE objects that are not members of
either an Mdb object or a Session object; for example, XYData and Leaf objects. These objects are
sometimes referred to as “temporary,” and the delete semantics for these objects are the same as for
standard Python objects. The object exists as long as there is a reference to it. If you delete the reference,
the object is also deleted.
6–13
Abaqus ID:
Printed on:
Abaqus/CAE SEQUENCES
GeomSequence
A GeomSequence is a sequence of geometry objects, such as Vertices or Edges. An Edge sequence
is derived from the GeomSequence object. Use the len() function to determine the number of
objects in a GeomSequence. A GeomSequence has methods and members too.
For example, the following creates a three-dimensional part by extruding a 70 × 70 square
through a distance of 20. The members of the resulting Part object are listed along with some
information about the sequence of Edge objects.
mdb.Model('Body')
mySketch = mdb.models['Body'].ConstrainedSketch(
name='__profile__', sheetSize=200.0)
mySketch.rectangle(point1=(0.0, 0.0),
point2=(70.0, 70.0))
switch = mdb.models['Body'].Part(name='Switch',
dimensionality=THREE_D, type=DEFORMABLE_BODY)
switch.BaseSolidExtrude(sketch=mySketch, depth=20.0)
The following statement displays the members of the resulting three-dimensional part:
>>> print mdb.models['Body'].parts['Switch'].__members__
['allInternalSets', 'allInternalSurfaces', 'allSets',
'allSurfaces', 'cell', 'cells', 'datum', 'datums', 'edge',
'edges', 'elemEdge', 'elemEdges', 'elemFace', 'elemFaces',
'element', 'elements', 'engineeringFeatures', 'face',
'faces', 'feature', 'featureById', 'features',
'featuresById', 'geometryPrecision', 'geometryRefinement',
'geometryValidity', 'ip', 'ips', 'isOutOfDate', 'modelName',
'name', 'node', 'nodes', 'referencePoint', 'referencePoints',
'reinforcement', 'reinforcements', 'sectionAssignments',
'sets', 'space', 'surfaces', 'twist', 'type',
'vertex', 'vertices']
The edges, faces, vertices, cells, and ips members are all derived from the GeomSequence object.
The following statements display some information about the edges sequence:
>>> print 'Single edge type = ', type(switch.edges[0])
Single edge type = <type 'Edge'>
6–14
Abaqus ID:
Printed on:
NAMESPACE
MeshSequence
A sequence of Nodes or Elements.
SurfSequence
A sequence of Surfaces.
6.4 Namespace
Namespace is an important concept for the Abaqus Scripting Interface programmer. A namespace can
be thought of as a program execution environment, and namespaces are independent of each other.
Namespaces prevent conflict between variable names. You can use the same variable name to refer
to different objects in different name spaces. Figure 6–9 illustrates how commands interact with the
Abaqus/CAE kernel. Abaqus Scripting Interface commands are issued to the Python interpreter from
either the GUI, the command line interface, or a script. Abaqus/CAE executes these commands in one
of two namespaces.
Script namespace
Abaqus Scripting Interface commands issued from scripts and from the command line interface are
executed in the script namespace. Commands issued from File→Run Script are also executed in
the script namespace. For example, if you enter the following statement from the command line
interface to create a viewport:
myViewport = session.Viewport(name='newViewport',
width=100, height=100)
the variable myViewport exists only in the script namespace. The name of the script namespace
is main.
Journal namespace
Abaqus Scripting Interface commands issued by the GUI are executed in the journal namespace.
For example, if you use the GUI to partition an edge, Abaqus/CAE writes the following statements
to the replay file, abaqus.rpy:
p1 = mdb.models['Model A'].parts['Part 3D A']
e = p1.edges
edges =(e[23], )
6–15
Abaqus ID:
Printed on:
NAMESPACE
Abaqus/CAE
command
line
GUI interface script
(CLI)
commands
Python replay
interpreter files
Abaqus/CAE
kernel
input file
Abaqus/Standard
Abaqus/Explicit
Abaqus/CFD
Abaqus/Design
output database
Figure 6–9 The Abaqus Scripting Interface and the Abaqus/CAE kernel.
p1.PartitionEdgeByParam(edges=edges, parameter=0.5)
The variables defined in the replay file (p1, e, and edges, in the above example) exist only in
the journal namespace. Abaqus/CAE issues an exception if you attempt to refer to one of these
variables from the script namespace. For example, the following statement was issued from the
command line interface and tries to partition the same edge:
p1.PartitionEdgeByParam(edges=edges, parameter=0.75)
NameError: p1
The name of the journal namespace is journaling.
The statement from abaqus import * described in “Executing scripts,” Section 5.1, imports
the mdb variable into the script namespace. You can then use the mdb variable in your scripts to access the
6–16
Abaqus ID:
Printed on:
SPECIFYING WHAT IS DISPLAYED IN THE VIEWPORT
objects in the object model. Although variables in one namespace are not visible to the other namespace,
the object repositories are now available in both. As a result, an object created in one namespace can still
be referred to in another namespace if you use its full path (mdb.models[’Model A’]...) and its
repository key.
For example, although the variable p1 in the above statement cannot be accessed from the script
namespace, you can still use the command line interface to access the part to which p1 referred.
The model and part repositories are available in both the journal and script namespaces. You can also
create your own variable p1 from the command line interface or from a script.
p1 = myPart
The variable p1 in the script namespace is independent of the variable p1 in the journal namespace.
• A part
• The assembly
• A sketch
• Data from an output database
• An X–Y plot
• Empty
In some cases you will want to update the contents of the viewport as the model changes; for example, to
illustrate how the assembly was partitioned prior to meshing. However, frequent updates to a viewport
will slow down your script, and you may want to leave the viewport empty until the script has completed.
Alternatively, you can display an object that the script is not operating on; for example, you can display
a part while the script operates on the assembly.
You use the following command to change the contents of a specified viewport:
session.viewports[name].setvalues(displayedObject=object)
The displayedObject argument can be a Part, Assembly, Sketch, Odb, or XYPlot object or None.
If displayedObject=None, Abaqus/CAE displays an empty viewport. For more information, see
“setValues,” Section 11.4.33 of the Abaqus Scripting Reference Manual.
6–17
Abaqus ID:
Printed on:
SPECIFYING A REGION
centerNSet = odb.rootAssembly.nodeSets['CENTER']
The following script illustrates how you can create a region. Regions are created from each of the
following:
• A sequence of tuples indicating the vertices, edges, faces, or cells in the region. The sequence can
include multiple tuples of the same type.
6–18
Abaqus ID:
Printed on:
SPECIFYING A REGION
• A sequence of tuples indicating a combination of the vertices, edges, faces, and cells in the region.
The tuples can appear in any order in the sequence. In addition, you can include multiple tuples of
the same type, and you can omit any type from the sequence.
• A Surface object specifying an entity and the side of the entity.
Use the following command to retrieve the script:
abaqus fetch job=createRegions
"""
createRegions.py
myModel = mdb.models['Model-1']
myViewport=session.Viewport(name='Region syntax',
origin=(20, 20), width=200, height=100)
mySketch.rectangle(point1=(-40.0, 30.0),
point2=(-10.0, 0.0))
mySketch.rectangle(point1=(10.0, 30.0),
point2=(40.0, 0.0))
6–19
Abaqus ID:
Printed on:
SPECIFYING A REGION
door = myModel.Part(name='Door',
dimensionality=THREE_D, type=DEFORMABLE_BODY)
door.BaseSolidExtrude(sketch=mySketch, depth=20.0)
myAssembly = myModel.rootAssembly
doorInstance = myAssembly.Instance(name='Door-1',
part=door)
pillarVertices = doorInstance.vertices.findAt(
((-40,30,0),), ((40,0,0),) )
myModel.StaticStep(name='impact',
previous='Initial', initialInc=1, timePeriod=1)
myPillarLoad = myModel.ConcentratedForce(
name='pillarForce', createStepName='impact',
region=(pillarVertices,), cf1=12.50E4)
topFace = doorInstance.faces.findAt(((-25,30,10),))
bottomFace = doorInstance.faces.findAt(((-25,0,10),))
myFenderLoad = myModel.Pressure(
6–20
Abaqus ID:
Printed on:
SPECIFYING A REGION
name='pillarPressure', createStepName='impact',
region=((topFace+bottomFace, SIDE1),),
magnitude=10E4)
myEdge1 = doorInstance.edges.findAt(((10,15,20),))
myEdge2 = doorInstance.edges.findAt(((10,15,0),))
myDisplacementBc= myModel.DisplacementBC(
name='xBC', createStepName='impact',
region=(pillarVertices, myEdge1+myEdge2,
topFace), u1=5.0)
faceRegion = doorInstance.faces.findAt(
((-30,15,20), ), ((30,15,20),))
mySurface = myModel.rootAssembly.Surface(
name='exterior', side1Faces=faceRegion)
myFoundation = myModel.ElasticFoundation(
name='elasticFloor', createStepName='Initial',
surface=mySurface, stiffness=1500)
myViewport.setValues(displayedObject=myAssembly)
myViewport.assemblyDisplay.setValues(step='impact',
loads=ON, bcs=ON, fields=ON)
6–21
Abaqus ID:
Printed on:
PROMPTING THE USER FOR INPUT
You may want to request input from a user while an Abaqus Scripting Interface script is executing. There
are many reasons for requesting input; for example, to specify design parameters, to enable a macro to
take action based on the input received, or to force parts of the script to be repeated. The Abaqus Scripting
Interface provides three functions that request input from the user and return the data entered by the user.
• The getInput function requests a single input from the user from a text field in a dialog box. For
more information, see “Requesting a single input from the user,” Section 6.7.1.
• The getInputs function requests multiple inputs from the user from text fields in a dialog box.
For more information, see “Requesting multiple inputs from the user,” Section 6.7.2.
• The getWarningReply function requests a reply to a warning message from the user from
a warning dialog box. For more information, see “Requesting a warning reply from the user,”
Section 6.7.3.
Note: You cannot use a script that contains getInput, getInputs or getWarningReply if you
are running the script from the command line and passing the script name to the command line options
-start,-replay or -noGUI.
The float function on the third line converts the string returned by getInput into a floating point
number. The following figure shows the dialog box that appears when this script is executed:
6–22
Abaqus ID:
Printed on:
PROMPTING THE USER FOR INPUT
The next example shows how to modify a macro recorded by the Macro Manager in Abaqus/CAE
to use the getInput function. The following text shows a macro named createViewport
that was recorded by Abaqus/CAE while the user created a viewport. Macros are stored in the file
abaqusMacros.py in your local or home directory.
from abaqus import *
def createViewport():
session.Viewport(name='Viewport: 2',
origin=(15.0,15.0), width=145.0,
height=90.0)
session.viewports['Viewport: 2'].makeCurrent()
The following shows how you can modify the macro to accept input from the user. Default values
for the viewport width and height have been added to the input request.
from abaqus import *
def createViewport():
name = getInput('Enter viewport name:')
prompt = 'Enter viewport width, height (mm):'
w, h = eval(getInput(prompt, '100,50'))
vp = session.Viewport(name=name, width=w, height=h)
vp.restore()
vp.makeCurrent()
The eval function in the third line of the macro converts the string returned by the getInput
function into two integers. When you supply the default values shown in this example to the getInput
function, the prompt and the text field in the dialog box that appears are shown in the following figure.
If the user clicks OK or presses [Enter], the default values are accepted and returned to the getInput
function. If the user clicks Cancel, None is returned.
6–23
Abaqus ID:
Printed on:
PROMPTING THE USER FOR INPUT
The fields argument to the getInputs method is a sequence of sequences of Strings. The inner
sequence is a pair of Strings that specifies the description of the text field and the default value of the
field. If the text field does not have a default value, you must specify an empty string; for example,
fields = (('Width',''), ('Length', ''), ('Height', ''))
length, width, height =
getInputs(fields=fields, label='Specify block dimensions:')
The label argument to the getInputs method is an optional label that appears across the top of the
dialog box. The dialogTitle argument is an optional string that appears in the title bar of the dialog box.
If the user clicks Cancel, the getInputs method returns a sequence of None objects. You can
check the first value in the sequence to determine if the user clicked Cancel; for example:
6–24
Abaqus ID:
Printed on:
INTERACTING WITH Abaqus/Standard, Abaqus/Explicit, AND Abaqus/CFD
The following figure shows the dialog box that appears when this script is executed:
The Job commands include methods that allow you to submit jobs to Abaqus/Standard, Abaqus/Explicit,
and Abaqus/CFD. This section describes how you can interact with Abaqus/Standard, Abaqus/Explicit,
and Abaqus/CFD and synchronize your scripts with the analysis job. The following topics are covered:
• “Processing messages from Abaqus/Standard, Abaqus/Explicit, and Abaqus/CFD,” Section 6.8.1
• “Waiting for a job to complete,” Section 6.8.2
• “An example of a callback function,” Section 6.8.3
6–25
Abaqus ID:
Printed on:
INTERACTING WITH Abaqus/Standard, Abaqus/Explicit, AND Abaqus/CFD
6–26
Abaqus ID:
Printed on:
INTERACTING WITH Abaqus/Standard, Abaqus/Explicit, AND Abaqus/CFD
– HEALER_TYPE
– INTERRUPTED
– ITERATION
– JOB_ABORTED
– JOB_COMPLETED
– JOB_INTERRUPTED
– JOB_SUBMITTED
– MONITOR_DATA
– ODB_FILE
– ODB_FRAME
– SIMULATION_ABORTED
– SIMULATION_COMPLETED
– SIMULATION_INTERRUPTED
– SIMULATION_SUBMITTED
– STARTED
– STATUS
– STEP
– WARNING
• data: A DataObject object containing the message data. The following list describes the members
of the DataObject object:
– clientHost: A String specifying the host name of the machine that is running the analysis.
– clientName: A String specifying the name of the client that sent the message. Possible values
are
– “BatchPre” (the input file preprocessor)
– “Packager” (the Abaqus/Explicit preprocessor packager)
– “Standard” (the Abaqus/Standard analysis)
– “Explicit” (the Abaqus/Explicit analysis)
– “Calculator” (the postprocessing calculator)
– phase: A SymbolicConstant specifying the phase of the analysis. Possible values are
• BATCHPRE_PHASE
• PACKAGER_PHASE
• STANDARD_PHASE
• EXPLICIT_PHASE
• CFD_PHASE
6–27
Abaqus ID:
Printed on:
INTERACTING WITH Abaqus/Standard, Abaqus/Explicit, AND Abaqus/CFD
• CALCULATOR_PHASE
• HEALER_PHASE
– processId: An Int specifying the process ID of the analysis program.
– threadId: An Int specifying the thread ID of the analysis program. Threads are used for parallel
or multiprocessing; in most cases threadId is set to zero.
– timeStamp: An Int specifying the time the message was sent in seconds since 00:00:00 UTC,
January 1, 1970.
• userData: Any Python object or None. This object is passed as the userData argument to
addMessageCallback.
The following script is an example of how you can use the messaging capability of the Abaqus
Scripting Interface. The callback function will intercept all messages from Abaqus/Standard,
Abaqus/Explicit, or Abaqus/CFD and print the messages in the Abaqus/CAE command line interface.
Use the following command to retrieve the example script:
abaqus fetch job=simpleMonitor
To execute the script, do the following:
• From the Abaqus/CAE command line interface type from simpleMonitor import
printMessages
• Submit an analysis job as usual.
• To start printing the messages, type printMessages(ON) from the Abaqus/CAE command line
interface.
• To stop printing the messages, type printMessages(OFF) from the Abaqus/CAE command
line interface.
"""
simpleMonitor.py
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
def simpleCB(jobName, messageType, data, userData):
"""
This callback prints out all the
members of the data objects
6–28
Abaqus ID:
Printed on:
INTERACTING WITH Abaqus/Standard, Abaqus/Explicit, AND Abaqus/CFD
"""
members = dir(data)
for member in members:
memberValue = getattr(data, member)
memberType = type(memberValue).__name__
print format%(member, memberType, memberValue)
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
def printMessages(start=ON):
"""
Switch message printing ON or OFF
"""
if start:
monitorManager.addMessageCallback(ANY_JOB,
ANY_MESSAGE_TYPE, simpleCB, None)
else:
monitorManager.removeMessageCallback(ANY_JOB,
ANY_MESSAGE_TYPE, simpleCB, None)
6–29
Abaqus ID:
Printed on:
INTERACTING WITH Abaqus/Standard, Abaqus/Explicit, AND Abaqus/CFD
myJob1.submit()
try:
myJob1.waitForCompletion(6000)
except AbaqusException, message:
print "Job timed out", message
If you do not provide a timeOut argument, the script will wait indefinitely for the job to complete.
In the following example, the script submits myJob1 and waits for it to complete before submitting
myJob2.
myJob1 = mdb.Job(name='Job-1')
myJob2 = mdb.Job(name='Job-2')
myJob1.submit()
myJob1.waitForCompletion()
myJob2.submit()
myJob2.waitForCompletion()
If you submit more than one job and then issue a waitForCompletion statement, Abaqus waits until
the job associated with the waitForCompletion statement is complete before checking the status of
the second job. If the second job has already completed, the waitForCompletion method returns
immediately. In the following example the script will not check the status of myJob2 until myJob1 has
completed.
myJob1 = mdb.Job(name='Job-1')
myJob2 = mdb.Job(name='Job-2')
myJob1.submit()
myJob2.submit()
myJob1.waitForCompletion()
myJob2.waitForCompletion()
If the job submission is delayed using either the at or wait arguments, the interval specified by the
timeOut argument includes both the time required to wait and the time required to execute. It is possible
for the waitForCompletion method to time-out before the job is submitted. If this occurs, the script
continues, and the waitForCompletion method throws an exception to indicate that the job did not
run to completion because of a time-out; however, the computer system will still submit the job. As a
result, if you use the at argument, you should add statements that check the status of the job and proceed
accordingly.
6–30
Abaqus ID:
Printed on:
INTERACTING WITH Abaqus/Standard, Abaqus/Explicit, AND Abaqus/CFD
import visualization
viewport.odbDisplay.commonOptions.setValues(renderStyle=FILLED)
The following statements show how the example script can be modified to use the callback function.
After the first statement is executed, the callback function responds to all messages from the job named
Deform. The final two statements create the job and submit it for analysis; the example script has now
finished executing. When the job is complete, the callback function opens the resulting output database
and displays a contour plot.
...
myJobName = 'Deform'
monitorManager.addMessageCallback(jobName='myJobName',
messageType=ANY_MESSAGE_TYPE, callback=onMessage,
userData=myViewport)
6–31
Abaqus ID:
Printed on:
USING Abaqus SCRIPTING INTERFACE COMMANDS IN YOUR ENVIRONMENT FILE
def onCaeGraphicsStartup():
# Graphics preferences
#
session.defaultGraphicsOptions.setValues(
displayLists=OFF, dragMode=AS_IS)
def onCaeStartup():
# Print preferences
#
session.printOptions.setValues(vpDecorations=OFF,
vpBackground=OFF, rendition=COLOR,
printCommand='lpr')
session.psOptions.setValues(date=OFF)
# Job preferences
#
def setJobPreferences(module, userData):
import job
6–32
Abaqus ID:
Printed on:
USING Abaqus SCRIPTING INTERFACE COMMANDS IN YOUR ENVIRONMENT FILE
session.Queue(name='long', hostName='server',
queueName='large', directory='/tmp')
addImportCallback('job', setJobPreferences)
# Visualization preferences
#
def setVisPreferences(module, userData):
import visualization
session.defaultOdbDisplay.contourOptions.setValues(
renderStyle=SHADED, visibleEdges=EXTERIOR,
contourStyle=CONTINUOUS)
addImportCallback('visualization', setVisPreferences)
The addImportCallback statement instructs Abaqus to call a function when the user first
imports a module. In this example Abaqus calls the setJobPreferences function when the user
first enters the Job module, and Abaqus calls the setVisPreferences function when the user first
enters the Visualization module. The setJobPreferences function creates a queue on a remote
host. The setVisPreferences function sets default options for contour plots.
The example environment file uses the onCaeStartup() function to control a set of Python
statements that are executed when Abaqus/CAE first starts. The environment file can also contain the
following:
• The onJobStartup() function controls a set of statements that execute when an analysis job
starts. For example,
def onJobStartup():
import os, shutil
restartDir = savedir + id + '_restart'
if (os.path.exists(restartDir)):
shutil.rmtree(restartDir)
• The onJobCompletion() function controls a set of statements that execute when an analysis
job completes. For example,
def onJobCompletion():
import os
extensions = ('res','stt','mdl','prt','abq','pac')
restartDir = savedir + os.sep + id + '_restart'
if (not os.path.exists(restartDir)):
os.mkdir(restartDir)
for extension in extensions:
fileName = id + '.' + extension
if (os.path.exists(savedir + os.sep + fileName)):
6–33
Abaqus ID:
Printed on:
USING Abaqus SCRIPTING INTERFACE COMMANDS IN YOUR ENVIRONMENT FILE
savedir
The path to the directory from which the job was submitted.
scrdir
The path to the scratch directory.
analysisType
The type of analysis to be executed. Possible values are STANDARD and EXPLICIT.
For a list of the variables that are available outside of the onJobStartup() and
onJobCompletion() functions, see “Job variables,” Section 4.1.9 of the Abaqus Installation and
Licensing Guide.
For more information on the environment file, see “Using the Abaqus environment settings,”
Section 3.3.1 of the Abaqus Analysis User’s Manual, and Chapter 4, “Customizing the Abaqus
environment,” of the Abaqus Installation and Licensing Guide.
6–34
Abaqus ID:
Printed on:
Part III: The Abaqus Python development
environment
This section describes the Abaqus Python development environment (PDE). The Abaqus PDE provides a
simple interface that you can use to develop—create, edit, test, and debug—Python scripts. The Abaqus PDE
is primarily intended for use with Abaqus/CAE user interface (GUI) and kernel scripts, including plug-ins,
but you can also use it to work on scripts that function independently from Abaqus/CAE. The following topic
is covered:
• Chapter 7, “Using the Abaqus Python development environment”
Abaqus ID:
Printed on:
AN OVERVIEW OF THE Abaqus PYTHON DEVELOPMENT ENVIRONMENT
The Abaqus PDE is a separate application that you can access from within Abaqus/CAE or launch
independently to work on Python scripts. It is intended primarily for use with scripts that use the
Abaqus/CAE graphical user interface (GUI) or kernel commands, including plug-ins, but you can
also use it to work on scripts that are unrelated to Abaqus. The Abaqus PDE also enables you to set
breakpoints to pause script execution at a particular line in any Python script, including an Abaqus
plug-in.
Figure 7–1 shows a .guiLog file in the Abaqus PDE. The script creates an extruded solid
rectangular part named “box1” and was recorded by logging the actions to complete the task in the
Abaqus/CAE user interface.
The PDE controls allow you to complete the following tasks:
7–1
Abaqus ID:
Printed on:
Abaqus PDE BASICS
7–2
Abaqus ID:
Printed on:
Abaqus PDE BASICS
Note: Using this method starts Abaqus/CAE without any local user preference settings. Ignoring
user preferences allows you to record and run .guiLog tests using the consistent default startup
settings.
where abaqus is the command used to start Abaqus, and filenames are the names, including the
directory paths, of scripts to be opened at startup.
The -script option allows you to enter the name, including the directory path, of a main file
to be opened at startup. The Abaqus PDE will create a new blank script if the named file does not
exist in the specified directory. If the directory does not exist, the Abaqus PDE generates an error
message.
Note: File names and paths specified without the -script option are opened for editing—not as
the main file.
The -pde option is used to specify options for use with Abaqus/CAE if you run a script in the
Abaqus PDE that requires the Abaqus/CAE kernel or user interface. You can also add command
line options for Abaqus/CAE using the Settings menu. For more information, see “Selecting the
settings for use with a file,” Section 7.2.4.
The sections that follow describe how to use the menus and tools within the Abaqus PDE.
7–3
Abaqus ID:
Printed on:
Abaqus PDE BASICS
7–4
Abaqus ID:
Printed on:
Abaqus PDE BASICS
New
Create a new file. The Abaqus PDE creates a new main file and displays it in the main window.
The file is created using the default naming convention _abaqus#_.guiLog, where # starts at
1 and is incremented as you create more files in the current directory. You can also click the New
Open
Open a script. You can also click the Open file icon to open a script.
If you have not yet opened or created another script, Abaqus automatically makes the first
opened file the main file for testing. Otherwise, the file opened becomes the current file viewed in
the main window, but it is not the main file used for testing.
Tip: You can drag and drop script files from the desktop or from Windows Explorer into
the Abaqus PDE for editing.
Open a script as the main file for testing. You can also click the Open main file icon to open
a script as the main file.
Save
Save changes to the current file. You can also click Save to save the current file.
Save As
Save the current file with a new name.
Reload Modules
Reload user interface modules to capture any changes that you made since they were first loaded.
You can also click Reload Modules to reload the user interface modules. The Abaqus PDE
reloads user interface modules in the Abaqus/CAE GUI and Abaqus/CAE kernel processes unless
the current setting for the Run Script In option is “local,” in which case any changed modules are
reloaded in the local PDE process.
Close File
Close the current file.
7–5
Abaqus ID:
Printed on:
Abaqus PDE BASICS
Filename.py
The name and file extension of the current main file, if one is selected.
Clicking here shows a list of dependent files that were found when the main file was run. If
the current main file has not been run in the Abaqus PDE, this list will be empty.
Recent Files
A list of all files that you have opened, regardless of whether you opened them to view and edit
them or opened them as the main file for testing. Recent Files from previous sessions will be read
from the abaqus_pde.deps file, if it exists in the current directory.
The recently used files lists are stored in the abaqus_pde.deps file in the directory from which you
start the current Abaqus PDE session. If you start an Abaqus PDE session from another location, the
lists contain only the files that you used the last time you opened a session in that directory. If you have
not previously used the Abaqus PDE in the current directory, a new set of recently used files is recorded
as you work.
7–6
Abaqus ID:
Printed on:
Abaqus PDE BASICS
7–7
Abaqus ID:
Printed on:
Abaqus PDE BASICS
Run Script In
Select whether the main file is run in the Abaqus/CAE GUI, the Abaqus/CAE kernel, or run locally.
By default, .guiLog files are run in the GUI, and .py and other file types are run in the kernel.
You can also set this option using the GUI, Kernel, and Local radio buttons located above the main
window.
If the Abaqus PDE was opened without Abaqus/CAE and you run a script with the GUI or
Kernel process, the Abaqus PDE will start Abaqus/CAE to run the script.
Line Animation
Highlight the line currently being executed in the main window. The following animation settings
are available:
• No animation.
• Animate main file (default). Highlights only the statements in the main function or method.
Functions called from the main script are not highlighted.
• Animate main file functions. Highlights the main script statements and the statements in
functions that are defined within the main file.
• Animate all files. Highlights the main script statements and statements within all functions for
which the source code is available.
Python Code
Control the appearance and editing behavior of Python scripts in the Abaqus PDE main window.
Syntax Coloring
Display the code using various font colors according to its purpose. This option is selected by
default.
You can view or change the color selections with the Choose Syntax Colors option.
7–8
Abaqus ID:
Printed on:
Abaqus PDE BASICS
Python Editing
Edit scripts with Python formatting, such as indentation, included automatically. This option
is selected by default.
message area. Click in the bottom left corner of the Abaqus PDE main window to switch from the
message area to the GUI command line interface. The GUI and kernel processes in Abaqus/CAE run
separately, each using its own Python interpreter. You can use the GUI command line interface to type
Python commands and to evaluate mathematical expressions using the Python interpreter that is built into
the Abaqus/CAE GUI. You can use the kernel command line interface in Abaqus/CAE for similar tasks.
Each command line interface includes primary (>>>) and secondary (...) prompts to indicate when you
must indent commands to comply with Python syntax. After you use the GUI command line interface,
7–9
Abaqus ID:
Printed on:
USING THE Abaqus PDE
The following sections contain detailed information that you can use to create and work with files in the
Abaqus PDE:
• “Creating .guiLog files,” Section 7.3.1
• “Running a script,” Section 7.3.2
• “Using the debugger,” Section 7.3.3
• “Using breakpoints,” Section 7.3.4
• “Using the Abaqus PDE with plug-ins,” Section 7.3.5
• “Using the Abaqus PDE with custom applications,” Section 7.3.6
Tip: You can also click the New guiLog icon to create a new .guiLog file.
2. Click the Start Recording icon to begin recording actions from Abaqus/CAE.
Abaqus writes the following two lines to begin the file:
from abaqusTester import *
import abaqusGui
3. Complete all the desired actions in the Abaqus/CAE session to record them in the .guiLog file.
Note: When you record .guiLog files, do not use mouse button 2 to close the dialog box for a
procedure. Instead, use the buttons in the dialog box to close it. Using mouse button 2 adds multiple
7–10
Abaqus ID:
Printed on:
USING THE Abaqus PDE
dialog box closing commands to the recorded .guiLog file. Since only one command is needed
to close the dialog, the extra commands will result in an error when the recorded script is played.
of Main File to position the cursor at the end of the file, then repeat Step 2 through Step 4.
7. Select File→Save to save the file or File→Save As to save the file with a new name; new files
automatically use Save As.
Use the other buttons—Next Line , Stop , Go to Start , and Go to End —to
execute the main file one line at a time, stop running the file, or reposition the cursor at the beginning or
end of the file, respectively.
As you run a script, you might want to specify a breakpoint to pause script execution at a particular
line. For more information about breakpoints, see “Using breakpoints,” Section 7.3.4
Window→Debugger or click Start debugger . If you have a script paused in the main window,
7–11
Abaqus ID:
Printed on:
USING THE Abaqus PDE
the debugger opens at the current position of the test. If you do not have a paused script, the debugger
automatically begins running the main file and positions the cursor at the start of the script.
The debugger consists of a call stack area, action buttons, and the debugger command line interface
(CLI) window, as shown in Figure 7–4. The debugger is positioned between the Abaqus PDE main
window and the message area.
Debugger Command
Call Stack area Debug buttons Line Interface (CLI)
The debugger uses a custom Python module named atxPdb, based on the Python pdb module. You can
enter Pdb commands in the debugger CLI; you can also enter any Python statements in the debugger
CLI. Python statements are executed in the same process—GUI, kernel, or local—that is running the
current script (for more information, see “Running a script,” Section 7.3.2).
Note: If you enter a command in the debugger CLI and it does not seem to work, it may be conflicting
with a Pdb command.
The call stack area shows the commands that are currently being executed.
7–12
Abaqus ID:
Printed on:
USING THE Abaqus PDE
The debugger contains a tool to watch variables as scripts are executed. To add a variable to the
watch list, click mouse button 3 over the variable name in the Abaqus PDE main window and select Add
Watch: variable name from the menu that appears. The Abaqus PDE adds the variable to the watch
list, indicating the namespace that the variable is defined within, the variable name, the type of data the
variable can store, the current value, and the file and line where the variable is located. When you execute
or step through the script, the variable information is updated as needed. You can also activate the watch
list tool after starting the debugger by selecting Show Watch from the action buttons below the call
stack area or by selecting Window→Debug Windows→Watch List. Abaqus displays the Watch List
area below the debugger or below the main window if the debugger is not open.
You can also customize the following aspects of variable display in the Watch List:
• You can switch a variable’s display format between repr and str formats. The repr() function
returns a string that is better suited to being read by the Python interpreter, while the str() function
returns a string better suited to printing. For more information about the built-in Python functions
repr() and str(), refer to the official Python web site (www.python.org).
To toggle between these settings, click mouse button 3 on a watch variable row and select
Display repr (not str) value of variable name or Display str (not repr) value of variable
name from the list that appears. If the variable is a local variable and the program is not accessing
that section of code, the variable value will be set to “not evaluated.” Variable values are also set to
“not evaluated” if the program is not running.
• You can prompt the Abaqus PDE to pause when the program reaches a line in which the value of
a selected watch variable has changed. To toggle on the “stop on change” option for a particular
watch variable, click mouse button 3 on the variable’s line and select Stop on change to variable
name from the list that appears. When this option is selected, Abaqus PDE stops at the line after
the change.
• You can remove any watch variables from the debugger by clicking mouse button 3 on the variable’s
line and selecting Delete watch var variable name from the list that appears.
You can also add and remove breakpoints using the breakpoint tool located above the main window
or the [F9] key.
Breakpoints are indicated by an asterisk to the right of the line number in the Abaqus PDE. If syntax
colors are active, the line number, asterisk, and the line of code are colored using the current breakpoint
color selection (for more information, see “Selecting the settings for use with a file,” Section 7.2.4).
7–13
Abaqus ID:
Printed on:
USING THE Abaqus PDE
7–14
Abaqus ID:
Printed on:
USING THE Abaqus PDE
Note: You cannot start the custom application and launch the Abaqus PDE from within
Abaqus/CAE since the initial startup processes would already be complete.
2. Open the file that you want to debug.
3. Position the cursor where you want to add a breakpoint. Click mouse button 3, and select Add
Breakpoint. (For more information, see “Using breakpoints,” Section 7.3.4.)
4. Click Start CAE at the top right of the Abaqus PDE to start Abaqus/CAE with the custom startup
commands.
5. The application code appears in the Abaqus PDE window, stopped at a breakpoint.
6. Use the Abaqus PDE controls and options described in the previous sections to step through the
execution of the custom application.
7–15
Abaqus ID:
Printed on:
Part IV: Putting it all together: examples
The section provides examples that illustrate how you can combine Abaqus Scripting Interface commands
and Python statements to create your own scripts. You can use the scripts to create Abaqus/CAE models,
submit jobs for analysis, and view the results. The following topic is covered:
• Chapter 8, “Abaqus Scripting Interface examples”
For examples of scripts that read and write from an output database, see “Example scripts that access
data from an output database,” Section 9.10.
Abaqus ID:
Printed on:
Abaqus SCRIPTING INTERFACE EXAMPLES
• How you can use commands from the Abaqus Scripting Interface to create a simple model, submit
it for analysis, and view the results. “Reproducing the cantilever beam tutorial,” Section 8.1,
uses Abaqus Scripting Interface commands to reproduce the cantilever beam tutorial described in
Appendix B, “Creating and Analyzing a Simple Model in Abaqus/CAE,” of Getting Started with
Abaqus: Interactive Edition.
• How you can use the Abaqus Scripting Interface to control the output from the Visualization module
in Abaqus/CAE (Abaqus/Viewer).
– “Opening the tutorial output database,” Section 8.2.1, explains how to use abaqus fetch to
retrieve the Abaqus/CAE tutorial output database.
– “Opening an output database and displaying a contour plot,” Section 8.2.2, explains how to
open the tutorial output database, display a contour plot, and print the resulting viewport to a
file.
– “Printing a contour plot at the end of each step,” Section 8.2.3, explains how to open the tutorial
output database, customize the legend, display a contour plot at the end of each step, and print
the resulting viewports to a file.
• How you can introduce more complex programming techniques into your Abaqus Scripting
Interface scripts. “Investigating the skew sensitivity of shell elements,” Section 8.3, reproduces the
problem found in “Skew sensitivity of shell elements,” Section 2.3.4 of the Abaqus Benchmarks
Manual. You use Abaqus/CAE to create the model, and you use Abaqus Scripting Interface
commands to parameterize an evaluation of the model by changing its geometry and element type.
The example investigates the sensitivity of the shell elements in Abaqus to skew distortion when
they are used as thin plates.
• How you can use functions available in the caePrefsAccess module to edit the display
preferences and GUI settings in the abaqus_v6.11.gpr file. “Editing display preferences and
GUI settings,” Section 8.4, describes how to query for and set several default display and GUI
behaviors in Abaqus/CAE.
The example scripts from this manual can be copied to the user’s working directory by using the
Abaqus fetch utility:
where scriptName.py is the name of the script (see “Fetching sample input files,” Section 3.2.13 of the
Abaqus Analysis User’s Manual).
8–1
Abaqus ID:
Printed on:
REPRODUCING THE CANTILEVER BEAM TUTORIAL
This example uses Abaqus Scripting Interface commands to reproduce the cantilever beam tutorial
described in Appendix B, “Creating and Analyzing a Simple Model in Abaqus/CAE,” of Getting Started
with Abaqus: Interactive Edition. Figure 8–1 illustrates the model that you will create and analyze.
0.5 MPa
20 mm
25 mm
200 mm
8–2
Abaqus ID:
Printed on:
REPRODUCING THE CANTILEVER BEAM TUTORIAL
"""
beamExample.py
# Create a model.
myModel = mdb.Model(name='Beam')
#-----------------------------------------------------
8–3
Abaqus ID:
Printed on:
REPRODUCING THE CANTILEVER BEAM TUTORIAL
import part
mySketch = myModel.ConstrainedSketch(name='beamProfile',
sheetSize=250.)
mySketch.rectangle(point1=(-100,10), point2=(100,-10))
myBeam.BaseSolidExtrude(sketch=mySketch, depth=25.0)
#-----------------------------------------------------
import material
# Create a material.
mySteel = myModel.Material(name='Steel')
#-------------------------------------------------------
import section
8–4
Abaqus ID:
Printed on:
REPRODUCING THE CANTILEVER BEAM TUTORIAL
mySection = myModel.HomogeneousSolidSection(name='beamSection',
material='Steel', thickness=1.0)
region = (myBeam.cells,)
myBeam.SectionAssignment(region=region,
sectionName='beamSection')
#-------------------------------------------------------
import assembly
myAssembly = myModel.rootAssembly
myInstance = myAssembly.Instance(name='beamInstance',
part=myBeam, dependent=OFF)
#-------------------------------------------------------
import step
myModel.StaticStep(name='beamLoad', previous='Initial',
timePeriod=1.0, initialInc=0.1,
description='Load the top of the beam.')
#-------------------------------------------------------
import load
endFaceCenter = (-100,0,12.5)
endFace = myInstance.faces.findAt((endFaceCenter,) )
8–5
Abaqus ID:
Printed on:
REPRODUCING THE CANTILEVER BEAM TUTORIAL
endRegion = (endFace,)
myModel.EncastreBC(name='Fixed',createStepName='beamLoad',
region=endRegion)
topFaceCenter = (0,10,12.5)
topFace = myInstance.faces.findAt((topFaceCenter,) )
#-------------------------------------------------------
import mesh
region = (myInstance.cells,)
elemType = mesh.ElemType(elemCode=C3D8I, elemLibrary=STANDARD)
myAssembly.setElementType(regions=region, elemTypes=(elemType,))
myAssembly.seedPartInstance(regions=(myInstance,), size=10.0)
myAssembly.generateMesh(regions=(myInstance,))
myViewport.assemblyDisplay.setValues(mesh=ON)
myViewport.assemblyDisplay.meshOptions.setValues(meshTechnique=ON)
myViewport.setValues(displayedObject=myAssembly)
8–6
Abaqus ID:
Printed on:
GENERATING A CUSTOMIZED PLOT
#-------------------------------------------------------
import job
jobName = 'beam_tutorial'
myJob = mdb.Job(name=jobName, model='Beam',
description='Cantilever beam tutorial')
myJob.submit()
myJob.waitForCompletion()
#-------------------------------------------------------
import visualization
myViewport.odbDisplay.commonOptions.setValues(renderStyle=FILLED)
The following section provides examples of Abaqus Scripting Interface scripts that open an output
database and generate a customized plot. In effect, these scripts reproduce the functionality of the
Visualization module in Abaqus/CAE. The following examples are provided:
8–7
Abaqus ID:
Printed on:
GENERATING A CUSTOMIZED PLOT
8–8
Abaqus ID:
Printed on:
GENERATING A CUSTOMIZED PLOT
odbPath = "viewer_tutorial.odb"
myOdb = visualization.openOdb(path=odbPath)
myViewport.setValues(displayedObject=myOdb)
myViewport.odbDisplay.display.setValues(plotState=(CONTOURS_ON_DEF,))
myViewport.odbDisplay.setFrame(step=2, frame=0)
myViewport.odbDisplay.contourOptions.setValues(numIntervals=10,
maxAutoCompute=OFF, maxValue=0.10,
minAutoCompute=OFF, minValue=0.0,)
session.printOptions.setValues(rendition=COLOR,
vpDecorations=OFF, vpBackground=OFF)
session.printToFile(fileName='contourPlot', format=PNG,
canvasObjects=(myViewport,))
8–9
Abaqus ID:
Printed on:
GENERATING A CUSTOMIZED PLOT
myViewport=session.Viewport(name=
'Print contour plot after each step', origin=(10, 10),
width=150, height=100)
try:
myOdb = visualization.openOdb(path='viewer_tutorial.odb')
myViewport.setValues(displayedObject=myOdb)
myViewport.odbDisplay.display.setValues(plotState=(CONTOURS_ON_DEF,))
8–10
Abaqus ID:
Printed on:
GENERATING A CUSTOMIZED PLOT
mySteps = myOdb.steps
numSteps = len(mySteps)
myViewport.odbDisplay.contourOptions.setValues(numIntervals=10,
maxAutoCompute=OFF, maxValue=0.1,
minAutoCompute=OFF, minValue=0.0)
session.printOptions.setValues(vpBackground=OFF)
session.psOptions.setValues(orientation=LANDSCAPE)
myViewport.viewportAnnotationOptions.setValues(
triad=OFF,title=OFF,state=OFF)
myViewport.odbDisplay.basicOptions.setValues(
coordSystemDisplay=OFF, )
for i in range(numSteps):
stepKey = mySteps.keys()[i]
step = mySteps[stepKey]
numFrames = len(step.frames)
myViewport.odbDisplay.setFrame(step=i, frame=numFrames-1)
myViewport.odbDisplay.display.setValues(plotState=(CONTOURS_ON_DEF,))
# Remove white space from the step key and use the result
# to name the file.
fileName=stepKey.replace(' ','')
8–11
Abaqus ID:
Printed on:
INVESTIGATING THE SKEW SENSITIVITY OF SHELL ELEMENTS
In this example you will use Abaqus/CAE to create the model and store the model in a model database.
The script opens the model database and performs a parametric study on the model. The example
illustrates how you can use a combination of Abaqus/CAE and the Abaqus Scripting Interface to analyze
a problem.
This example uses Abaqus Scripting Interface commands to evaluate the sensitivity of the shell
elements in Abaqus to skew distortion when they are used as thin plates. Further details can be found
in “Skew sensitivity of shell elements,” Section 2.3.4 of the Abaqus Benchmarks Manual. The problem
investigates the effects on the accuracy of the bending moment computed at the center of a shell using:
• different shell formulations and
• at different angles.
Figure 8–2 illustrates the basic geometry of the simply supported skew plate with a uniform
distributed load.
y
x, y, z displacements
constrained on boundary
m
0
1.
δ
x
1.0 m
The plate is loaded by a uniform pressure of 1.0 × 10−6 MPa applied over the entire surface. The edges of
the plate are all simply supported. The analysis is performed for five different values of the skew angle,
: 90°, 80°, 60°, 40°, and 30°. The analysis is performed for two different quadrilateral elements: S4
and S8R.
The example is divided into two scripts. The controlling script, skewExample.py, imports
skewExampleUtils.py. Use the fetch utility to retrieve the scripts:
abaqus fetch job=skewExample
abaqus fetch job=skewExampleUtils
8–12
Abaqus ID:
Printed on:
INVESTIGATING THE SKEW SENSITIVITY OF SHELL ELEMENTS
Material
Create a material, and name it Steel. The Young’s modulus is 30 MPa, and the Poisson’s ratio is
0.3.
Section
Create a homogeneous shell section that refers to the material called Steel. Name the section
Shell. The plate thickness is 0.01 m. The length/thickness ratio is, thus, 100/1 so that the plate is
thin in the sense that transverse shear deformation should not be significant. Assign the section to
the plate.
Assembly
Create the assembly using a single instance of the part. Abaqus/CAE names the part instance
Plate-1.
Step
Create a static step and name it Step-1. Enter Apply pressure for the step Description.
Accept the default time period of 1.0 and the default initial increment of 1.0.
8–13
Abaqus ID:
Printed on:
INVESTIGATING THE SKEW SENSITIVITY OF SHELL ELEMENTS
Boundary condition
Create a displacement boundary condition, and name it Pinned. The boundary condition pins the
exterior edges of the plate.
Load
Create a pressure load, and name it Pressure. Apply the load to the face of the plate. Accept the
default side of the plate and use a magnitude of 1.0. This positive pressure will result in a negative
displacement in the 3-direction.
Set
Partition the plate into quarters. Create a set that contains the vertex at the center of the plate, and
name the set CENTER.
Mesh
Create a 4 × 4 swept mesh of quadrilateral elements on the plate.
Keyword editor
You must use the Keyword Editor to request output of section forces and moments for the node
at the center of the plate. When you edited the output database request to select the variable, SF,
Abaqus/CAE requested output of section forces and moments at element integration points. To
modify this output request, you must add position=NODES to the OUTPUT REQUESTS block
as follows:
Job
Create a job, and name it skew. The job must refer to the model Model-1.
8–14
Abaqus ID:
Printed on:
INVESTIGATING THE SKEW SENSITIVITY OF SHELL ELEMENTS
4. From the Sketcher toolbox, select the dimension tool and dimension the angle at the lower left
corner of the plate as shown in Figure 8–3.
Figure 8–3 Dimension the angle at the lower left corner of the plate.
To determine the indices of the dimension to modify and the vertices to move:
8–15
Abaqus ID:
Printed on:
INVESTIGATING THE SKEW SENSITIVITY OF SHELL ELEMENTS
7. The example script, skewExample.py, contains a similar statement that modifies the angular
dimension of the plate. The index of the angular dimension in your model must be the same as the
index in the example script. If the indices are not the same, you must edit the example script and
enter the correct indices.
d[0].setValues(value=angle, )
Save the model database, and name it skew. Abaqus/CAE saves the model database in a file called
skew.cae. The example script opens this model database and parameterizes the model it contains.
8–16
Abaqus ID:
Printed on:
INVESTIGATING THE SKEW SENSITIVITY OF SHELL ELEMENTS
import part
import mesh
from mesh import S4, S8R, STANDARD, STRUCTURED
import job
from skewExampleUtils import getResults, createXYPlot
model = mdb.models['Model-1']
part = model.parts['Plate']
feature = part.features['Shell planar-1']
assembly = model.rootAssembly
instance = assembly.instances['Plate-1']
job = mdb.jobs['skew']
allFaces = instance.faces
regions =(allFaces[0], allFaces[1], allFaces[2], allFaces[3])
assembly.setMeshControls(regions=regions,
technique=STRUCTURED)
face1 = allFaces.findAt((0.,0.,0.), )
face2 = allFaces.findAt((0.,1.,0.), )
face3 = allFaces.findAt((1.,1.,0.), )
face4 = allFaces.findAt((1.,0.,0.), )
allVertices = instance.vertices
v1 = allVertices.findAt((0.,0.,0.), )
v2 = allVertices.findAt((0.,.5,0.), )
v3 = allVertices.findAt((0.,1.,0.), )
v4 = allVertices.findAt((.5,1.,0.), )
v5 = allVertices.findAt((1.,1.,0.), )
v6 = allVertices.findAt((1.,.5,0.), )
v7 = allVertices.findAt((1.,0.,0.), )
v8 = allVertices.findAt((.5,0.,0.), )
v9 = allVertices.findAt((.5,.5,0.), )
8–17
Abaqus ID:
Printed on:
INVESTIGATING THE SKEW SENSITIVITY OF SHELL ELEMENTS
elemName = repr(elemCode)
dispData[elemName], maxMomentData[elemName], \
minMomentData[elemName] = [], [], []
elemType = mesh.ElemType(elemCode=elemCode,
elemLibrary=STANDARD)
assembly.setElementType(regions=(instance.faces,),
elemTypes=(elemType,))
d[0].setValues(value=angle, )
feature.setValues(sketch=tmpSketch)
part.regenerate()
8–18
Abaqus ID:
Printed on:
INVESTIGATING THE SKEW SENSITIVITY OF SHELL ELEMENTS
assembly.regenerate()
assembly.setLogicalCorners(
region=face1, corners=(v1,v2,v9,v8))
assembly.setLogicalCorners(
region=face2, corners=(v2,v3,v4,v9))
assembly.setLogicalCorners(
region=face3, corners=(v9,v4,v5,v6))
assembly.setLogicalCorners(
region=face4, corners=(v8,v9,v6,v7))
assembly.generateMesh(regions=(instance,))
job.submit()
job.waitForCompletion()
print 'Completed job for %s at %s degrees' % (elemName,
angle)
disp, maxMoment, minMoment = getResults()
dispData[elemName].append((angle, disp))
maxMomentData[elemName].append((angle, maxMoment))
minMomentData[elemName].append((angle, minMoment))
The script imports two functions from skewExampleUtils. The functions do the following:
• Retrieve the displacement and calculate the maximum bending moment at the center of the plate.
• Display curves of theoretical and computed results in a new viewport.
"""
skewExampleUtils.py
8–19
Abaqus ID:
Printed on:
INVESTIGATING THE SKEW SENSITIVITY OF SHELL ELEMENTS
import visualization
from miscUtils import sorted
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
def getResults():
"""
Retrieve the displacement and calculate the minimum
and maximum bending moment at the center of plate.
"""
odb = visualization.openOdb('skew.odb')
centerNSet = odb.rootAssembly.nodeSets['CENTER']
frame = odb.steps['Step-1'].frames[-1]
dispField = frame.fieldOutputs['U']
dispSubField = dispField.getSubset(region=centerNSet)
disp = dispSubField.values[0].data[2]
momentField = frame.fieldOutputs['SM']
momentSubField = momentField.getSubset(region=centerNSet,
position=ELEMENT_NODAL)
m1, m2, m3 = 0, 0, 0
for value in momentSubField.values:
m1 = m1 + value.data[0]
m2 = m2 + value.data[1]
m3 = m3 + value.data[2]
numElements = len(momentSubField.values)
m1 = m1 / numElements
m2 = m2 / numElements
m3 = m3 / numElements
momentA = 0.5 * (abs(m1) + abs(m2))
momentB = sqrt(0.25 * (m1 - m2)**2 + m3**2)
maxMoment = momentA + momentB
minMoment = momentA - momentB
odb.close()
8–20
Abaqus ID:
Printed on:
EDITING DISPLAY PREFERENCES AND GUI SETTINGS
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
def createXYPlot(vpOrigin, vpName, plotName, data):
"""
Display curves of theoretical and computed results in
a new viewport.
"""
vp = session.Viewport(name=vpName, origin=vpOrigin,
width=150, height=100)
xyPlot = session.XYPlot(plotName)
chart = xyPlot.charts.values()[0]
curveList = []
for elemName, xyValues in sorted(data.items()):
xyData = session.XYData(elemName, xyValues)
curve = session.Curve(xyData)
curveList.append(curve)
chart.setValues(curvesToPlot=curveList)
chart.axes1[0].axisData.setValues(useSystemTitle=False,title='Skew Angle')
chart.axes2[0].axisData.setValues(useSystemTitle=False,title=plotName)
vp.setValues(displayedObject=xyPlot)
WARNING: Editing the abaqus_v6.11.gpr file is for experienced users only. Do not use the
functions in the caePrefsAccess module unless you are comfortable with the Abaqus Scripting
Interface and understand the structure of the abaqus_v6.11.gpr file. In addition, you should not
have Abaqus/CAE running when you make changes to the graphical preferences file.
You can retrieve the location of your abaqus_v6.11.gpr file using the
getGuiPrefsFileName function. The file records default settings in two sections: display options
reside in the sessionOptions section, and GUI settings reside in the guiPreferences section.
Editing the options in one section does not have any effect on the options in the other section.
8–21
Abaqus ID:
Printed on:
EDITING DISPLAY PREFERENCES AND GUI SETTINGS
sessionOptions
The session options consist of the settings that you can save using the File→Save Display Options
menu option. In Abaqus/CAE you can save these options in the current directory or in your home
directory.
You can display and edit session options using the openSessionOptions function.
> abaqus Python
...
>>> import caePrefsAccess
>>> sessionOptions = caePrefsAccess.openSessionOptions()
>>> printValuesList(sessionOptions)
...
sessionOptions['session.animationController.animationOptions']\
['frameCounter']:[type:bool] True
sessionOptions['session.animationController.animationOptions']\
['frameRate']:[type:int] 100
sessionOptions['session.aviOptions']['compressionMethod']:\
[type:SymbolicConstant] CODEC
sessionOptions['session.aviOptions']['compressionQuality']:[type:int] 75
...
The following statement changes the frame rate to 50. You should confirm that the data type
you specify matches the type of the existing value.
>>> sessionOptions['session.animationController.animationOptions']\
['frameRate'] = 50
You can save the options you change to the original file by issuing the following command:
>>> sessionOptions.save()
guiPreferences
The GUI preferences control many default behaviors in the Abaqus/CAE graphical interface,
including size and location of the main window, size and location of the dialog boxes within
Abaqus/CAE, and the number of recent files listed in the Start Session dialog box and in the
File menu.
Abaqus/CAE saves guiPreferences settings to your home directory when you exit
the application. A separate guiPreferences record is stored in the preferences file for each
display you use, so you must specify the displayName you want to modify when you open the
guiPreferences settings. You can obtain a list of the available displayName settings by
calling the getDisplayNamesInGuiPreferences function, and you can edit these settings
by using the openGuiPreferences function and specifying the displayName of the settings
that you want to modify.
In the following example, the openGuiPreferences function is used to examine the X-
and Y-location and the height and width of the following components of Abaqus/CAE:
8–22
Abaqus ID:
Printed on:
EDITING DISPLAY PREFERENCES AND GUI SETTINGS
You can change the geometry of the Abaqus/CAE main window by issuing a command like the
following:
>>> guiPreferences['Abaqus/CAE']['Geometry']['AFXMainWindow'] = '193,67,800,600'
You can save the options you change to the original file by issuing the following command:
>>> sessionOptions.save()
8–23
Abaqus ID:
Printed on:
Part V: Accessing an output database
This section describes how you access the data in an output database using either the Abaqus Scripting
Interface or the C++ Application Programming Interface (API). You can read model data and field and
history data from an output database. You can also write field and history data to an output database. The
following topics are covered:
• Chapter 9, “Using the Abaqus Scripting Interface to access an output database”
• Chapter 10, “Using C++ to access an output database”
The Abaqus Scripting Interface commands that read and write data from an output database are described
in Chapter 33, “Odb commands,” of the Abaqus Scripting Reference Manual.
The C++ commands that read and write data from an output database are described in Chapter 60, “Odb
commands,” of the Abaqus Scripting Reference Manual.
Abaqus ID:
Printed on:
HOW THE OBJECT MODEL FOR THE OUTPUT DATABASE RELATES TO COMMANDS
9.2 How the object model for the output database relates to
commands
You need to understand the object model for the output database both to read data from it and to write data
to it. An object model describes the relationship between objects. The object model for the Abaqus/CAE
model is described in “The Abaqus object model,” Section 6.1.
9–1
Abaqus ID:
Printed on:
OBJECT MODEL FOR THE OUTPUT DATABASE
For example, consider the object model for field output data shown in Figure 9–1. The Odb object
at the top of the figure is created when you issue the command to open or create an output database. As
you move down the object model, an OdbStep object is a member of the Odb object; similarly, a Frame
object is a member of the OdbStep object. The FieldOutput object has two members—fieldValue and
fieldLocation.
The object model translates directly to the structure of an Abaqus Scripting Interface command.
For example, the following command refers to a Frame object in the sequence of frames contained in an
OdbStep object:
odb.steps['10 hz vibration'].frames[3]
Similarly, the following command refers to the sequence of field data contained in a FieldOutput object.
odb.steps['10 hz vibration'].frames[3].\
fieldOutputs['U'].values[47]
You use commands to access objects by stepping through the hierarchy of objects in the object model. The
Access and Path descriptions in Chapter 33, “Odb commands,” of the Abaqus Scripting Reference
Manual describe the interface definition of the command. The interface definition of the command
reflects the hierarchy of objects in the object model. If you are unsure of the structure of the output
database, you can issue the objectname.__members__ command from the command line interface to
view the members of an object.
An output database generated from an Abaqus analysis contains both model and results data as shown
in Figure 9–1.
Model data
Model data describe the parts and part instances that make up the root assembly; for example, nodal
coordinates, set definitions, and element types. Model data are explained in more detail in “Model
data,” Section 9.3.1.
Results data
Results data describe the results of your analysis; for example, stresses, strains, and displacements.
You use output requests to configure the contents of the results data. Results data can be either field
output data or history output data; for a more detailed explanation, see “Results data,” Section 9.3.2.
Note: For a description of object models, see “An overview of the Abaqus object model,” Section 6.1.1.
You can find more information on the format of the output database in “Output to the output database,”
Section 4.1.3 of the Abaqus Analysis User’s Manual.
9–2
Abaqus ID:
Printed on:
OBJECT MODEL FOR THE OUTPUT DATABASE
odb
rootAssembly
parts
sections
= Container
= Singular object materials
Model Data
steps
frames
fieldOutputs
historyRegions
point
historyOutputs
Results Data
9–3
Abaqus ID:
Printed on:
OBJECT MODEL FOR THE OUTPUT DATABASE
Part instances
A part instance is a usage of a part within the assembly. All characteristics (such as mesh and
section definitions) defined for a part become characteristics for each instance of that part—they
are inherited by the part instances. Each part instance is positioned independently within the root
assembly.
Materials
Materials contain material models comprised of one or more material property definitions. The same
material models may be used repeatedly within a model; each component that uses the same material
model shares identical material properties. Many materials may exist within a model database, but
only the materials that are used in the assembly are copied to the output database.
Sections
Sections add the properties that are necessary to define completely the geometric and material
properties of an element. Various element types require different section types to complete their
definitions. For example, shell elements in a composite part require a section that provides a
thickness, multiple material models, and an orientation for each material model; all these pieces
combine to complete the composite shell element definition. Like materials, only those sections
that are used in the assembly are copied to the output database.
Section assignments
Section assignments link section definitions to the regions of part instances. Section assignments
in the output database maintain this association. Sections are assigned to each part in a model, and
the section assignments are propagated to each instance of that part.
Section categories
You use section categories to group the regions of the model that use the same section definitions;
for example, the regions that use a shell section with five section points. Within a section category,
you use the section points to identify the location of results; for example, you can associate section
point 1 with the top surface of a shell and section point 5 with the bottom surface.
9–4
Abaqus ID:
Printed on:
OBJECT MODEL FOR THE OUTPUT DATABASE
Rigid bodies
You use rigid bodies to define a collection of nodes, elements, and/or surfaces whose motion is
governed by the motion of a single node, called the rigid body reference node.
Pretension Sections
Pretension sections are used to associate a pre-tension node with a pre-tension section. The pre-
tension section can be defined using a surface for continuum elements or using an element for truss
or beam elements.
Interactions
Interactions are used to define contact between surfaces in an analysis. Only contact interactions
defined using contact pairs are written to the output database.
Interaction properties
Interaction properties define the physical behavior of surfaces involved in an interaction. Only
tangential friction behavior is written to the output database.
odb
rootAssembly
parts elementSets
sectionCategories nodeSets
materials surfaces
= Container instances
= Singular object elementSets
description nodeSets
name surfaces
elements
nodes
Model Data
steps
Step
frames
Frame
fieldOutputs
historyRegions
HistoryRegion
point
historyOutputs
Results Data
9–5
Abaqus ID:
Printed on:
OBJECT MODEL FOR THE OUTPUT DATABASE
The objects stored as model data in an output database are similar to the objects stored in an
Abaqus/CAE model database. However, the output database does not require a model name because
an analysis job always refers to a single model and the resulting output database can contain only one
model. For example, the following Abaqus Scripting Interface statements refer to an Instance object
in the model database:
mdb = openMdb(pathName='/users/smith/mdb/hybridVehicle')
myModel = mdb.models['Transmission']
myPart = myModel.rootAssembly.instances['housing']
Similar statements refer to an Instance object in the output database.
odb = openOdb(path='/users/smith/odb/transmission.odb')
myPart = odb.rootAssembly.instances['housing']
You can use the prettyPrint method to display a text representation of an output database
and to view the structure of the model data in the object model. For example, the following shows the
output from prettyPrint applied to the output database created by the Abaqus/CAE cantilever beam
tutorial:
from odbAccess import *
from textRepr import *
odb=openOdb('Deform.odb')
prettyPrint(odb,2)
({'analysisTitle': 'Cantilever beam tutorial',
'closed': False,
'description': 'DDB object',
'diagnosticData': ({'analysisErrors': 'OdbSequenceAnalysisError object',
'analysisWarnings': 'OdbSequenceAnalysisWarning object',
'jobStatus': JOB_STATUS_COMPLETED_SUCCESSFULLY,
'jobTime': 'OdbJobTime object',
'numberOfAnalysisErrors': 0,
'numberOfAnalysisWarnings': 0,
'numberOfSteps': 1,
'numericalProblemSummary': 'OdbNumericalProblemSummary object',
'steps': 'OdbSequenceDiagnosticStep object'}),
'isReadOnly': False,
'jobData': ({'analysisCode': ABAQUS_STANDARD,
'creationTime': 'date time year',
'machineName': '',
'modificationTime': 'date time year',
'name': 'Deform.odb',
'precision': SINGLE_PRECISION,
'productAddOns': 'tuple object',
'version': 'Abaqus/Standard release'}),
'name': 'Deform.odb',
'parts': {'BEAM': 'Part object'},
'path': 'C:/Deform.odb',
'rootAssembly': ({'connectorOrientations': 'ConnectorOrientationArray object',
'datumCsyses': 'Repository object',
'elementSet': 'Repository object',
'elementSets': 'Repository object',
'elements': 'OdbMeshElementArray object',
'instance': 'Repository object',
'instances': 'Repository object',
'name': 'ASSEMBLY',
'nodeSet': 'Repository object',
'nodeSets': 'Repository object',
'nodes': 'OdbMeshNodeArray object',
'sectionAssignments': 'Sequence object',
'surface': 'Repository object',
'surfaces': 'Repository object'}),
'sectionCategories': {'solid < STEEL >': 'SectionCategory object'},
9–6
Abaqus ID:
Printed on:
OBJECT MODEL FOR THE OUTPUT DATABASE
'sectorDefinition': None,
'steps': {'Beamload': 'OdbStep object'},
'userData': ({'annotations': 'Repository object',
'xyData': 'Repository object',
'xyDataObjects': 'Repository object'})})
For more information, see “prettyPrint,” Section 51.1.4 of the Abaqus Scripting Reference Manual.
Frames
Each step contains a sequence of frames, where each increment of the analysis that resulted in output
to the output database is called a frame. In a frequency or buckling analysis each eigenmode is stored
as a separate frame. Similarly, in a steady-state harmonic response analysis each frequency is stored
as a separate frame.
Field output
Field output is intended for infrequent requests for a large portion of the model and can be used
to generate contour plots, animations, symbol plots, and displaced shape plots in the Visualization
module of Abaqus/CAE. You can also use field output to generate an X–Y data plot. Only complete
sets of basic variables (for example, all the stress or strain components) can be requested as field
output. Field output is composed of a “cloud of data values” (e.g., stress tensors at each integration
point for all elements). Each data value has a location, type, and value. You use the regions defined
in the model data, such as an element set, to access subsets of the field output data. Figure 9–3
shows the field output data object model within an output database.
History output
History output is output defined for a single point or for values calculated for a portion of the model
as a whole, such as energy. History output is intended for relatively frequent output requests for
small portions of the model and can be displayed in the form of X–Y data plots in the Visualization
module of Abaqus/CAE. Individual variables (such as a particular stress component) can be
requested.
Depending on the type of output expected, a HistoryRegion object can be defined for one of
the following:
• a node
• an integration point
9–7
Abaqus ID:
Printed on:
OBJECT MODEL FOR THE OUTPUT DATABASE
odb
rootAssembly
parts
sections
materials
Model Data
= Container
= Singular object
steps
frames
fieldOutputs
fieldValues
historyRegions
Results Data
• a region
• the whole model
The output from all history requests that relate to a particular point or region is then collected in
one HistoryRegion object. Figure 9–4 shows the history output data object model within an output
database.
odb
rootAssembly
parts
= Container sections
= Singular object
materials
Model Data
steps
frames
fieldOutputs
historyRegions
point
historyOutputs
data
Results Data
9–8
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
The following sections describe how you use Abaqus Scripting Interface commands to read data from
an output database. The following topics are covered:
• “The Abaqus/CAE Visualization module tutorial output database,” Section 9.5.1
• “Making the Odb commands available,” Section 9.5.2
• “Opening an output database,” Section 9.5.3
• “Reading model data,” Section 9.5.4
• “Reading results data,” Section 9.5.5
• “Reading field output data,” Section 9.5.6
• “Using regions to read a subset of field output data,” Section 9.5.7
• “Reading history output data,” Section 9.5.8
• “An example of reading node and element information from an output database,” Section 9.5.9
• “An example of reading field data from an output database,” Section 9.5.10
9–9
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
9–10
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
After you open the output database, you can access its contents using the methods and members of the
Odb object returned by the openOdb method. In the above example the Odb object is referred to by the
variable odb. For a full description of the openOdb command, see “openOdb,” Section 33.32.5 of the
Abaqus Scripting Reference Manual.
Part instances
Part instances are stored in the instances repository under the OdbAssembly object. The
following statements display the repository keys of the part instances in the tutorial output database:
for instanceName in odb.rootAssembly.instances.keys():
print instanceName
The output database contains only one part instance, and the resulting output is
PART-1-1
Regions
Regions in the output database are OdbSet objects. Regions refer to the part and assembly sets stored
in the output database. A part set refers to elements or nodes in an individual part and appears in
each instance of the part in the assembly. An assembly set refers to the elements or nodes in part
instances in the assembly. A region can be one of the following:
• A node set
• An element set
• A surface
For example, the following statement displays the node sets in the OdbAssembly object:
print 'Node sets = ',odb.rootAssembly.nodeSets.keys()
9–11
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
Materials
Materials are stored in the materials repository under the Odb object.
Access the materials repository using the command:
allMaterials = odb.materials
for materialName in allMaterials.keys():
print 'Material Name : ',materialName
To print isotropic elastic material properties in a material object:
for material in allMaterials.values():
if hasattr(material,'elastic'):
elastic = material.elastic
if elastic.type == ISOTROPIC:
print 'isotropic elastic behavior, type = %s' \
% elastic.moduli
title1 = 'Young modulus Poisson\'s ratio '
title2 = ''
if elastic.temperatureDependency == ON:
title2 = 'Temperature '
9–12
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
dep = elastic.dependencies
title3 = ''
for x in range(dep):
title3 += ' field # %d' % x
print '%s %s %s' % (title1,title2,title3)
for dataline in elastic.table:
print dataline
Some Material definitions have suboptions. For example, to access the smoothing type used
for biaxial test data specified for a hyperelastic material:
if hasattr(material,'hyperelastic'):
hyperelastic = material.hyperelastic
testData = hyperelastic.testData
if testData == ON:
if hasattr(hyperelastic,'biaxialTestData'):
biaxialTestData = hyperelastic.biaxialTestData
print 'smoothing type : ',biaxialTestData.smoothing
Chapter 28, “Material commands,” of the Abaqus Scripting Reference Manual describes the
Material object commands in more detail.
Sections
Sections are stored in the sections repository under the Odb object.
The following statements display the repository keys of the sections in an output database:
allSections = odb.sections
for sectionName in allSections.keys():
print 'Section Name : ',sectionName
The Section object can be one of the various section types. The type command provides
information on the section type. For example, to determine whether a section is of type
“homogeneous solid section” and to print it’s thickness and associated material name:
for mySection in allSections.values():
if type(mySection) == HomogeneousSolidSectionType:
print 'material name = ', mySection.material
print 'thickness = ', mySection.thickness
Similarily, to access the beam profile repository:
allProfiles = odb.profiles.values()
numProfiles = len(allProfiles)
print 'Total Number of profiles in the ODB : %d' \
% numProfiles
9–13
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
The Profile object can be one of the various profile types. The type command provides
information on the profile type. For example, to output the radius of all circular profiles in the odb:
Section assignments
Section assignments are stored in the odbSectionAssignmentArray repository under the
OdbAssembly object.
All elements in an Abaqus analysis need to be associated with section and material properties.
Section assignments provide the relationship between elements in a part instance and their
section properties. The section properties include the associated material name. To access the
sectionAssignments repository from the PartInstance object:
instances = odb.rootAssembly.instances
for instance in instances.values():
assignments = instance.sectionAssignments
print 'Instance : ',instance.name
for sa in assignments:
region = sa.region
elements = region.elements
print ' Section : ',sa.sectionName
print ' Elements associated with this section : '
for e in elements:
print ' label : ',e.label
Rigid bodies
Rigid bodies are stored in the odbRigidBodyArray. The OdbPart object, OdbInstance object,
and OdbAssembly object each have an odbRigidBodyArray.
Pretension sections
Pretension sections are stored in odbPretensionSectionArray under the OdbAssembly
object.
9–14
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
Frames
Each step contains a sequence of frames, where each increment of the analysis (or each mode in an
eigenvalue analysis) that resulted in output to the output database is called a frame. The following
statement assigns a variable to the last frame in the first step:
lastFrame = odb.steps['Step-1'].frames[-1]
9–15
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
COPEN TARGET/IMPACTOR
CPRESS TARGET/IMPACTOR
CSHEAR1 TARGET/IMPACTOR
CSLIP1 TARGET/IMPACTOR
LE
RF
RM3
S
U
UR3
Different variables can be written to the output database at different frequencies. As a result, not all
frames will contain all the field output variables.
You can use the following to view all the available field data in a frame:
for f in lastFrame.fieldOutputs.values():
print f.name, ':', f.description
print 'Type: ', f.type
The resulting print output lists all the field output variables in a particular frame, along with their type
and position.
9–16
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
RF : Reaction force
Type: VECTOR
Position: NODAL
S : Stress components
Type: TENSOR_2D_PLANAR
Position: INTEGRATION_POINT
U : Spatial displacement
Type: VECTOR
Position: NODAL
In turn, a FieldOutput object has a member values that is a sequence of FieldValue objects that
contain data. Each data value in the sequence has a particular location in the model. You can query the
FieldValue object to determine the location of a data value; for example,
displacement=lastFrame.fieldOutputs['U']
fieldValues=displacement.values
for v in fieldValues:
print 'Node = %d U[x] = %6.4f, U[y] = %6.4f' % (v.nodeLabel,
v.data[0], v.data[1])
9–17
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
9–18
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
The arguments to getSubset are a region, an element type, a position, or section point data. The
following is a second example that uses an element set to define the region and generates formatted output.
For more information on tuples, the len() function, and the range() function, see “Sequences,”
Section 4.5.4, and “Sequence operations,” Section 4.5.5.
topCenter = \
odb.rootAssembly.instances['PART-1-1'].elementSets['CENT']
stressField = odb.steps['Step-2'].frames[3].fieldOutputs['S']
field = stressField.getSubset(region=topCenter,
position=INTEGRATION_POINT, elementType = 'CAX4')
fieldValues = field.values
for v in fieldValues:
print 'Element label = ', v.elementLabel,
if v.integrationPoint:
print 'Integration Point = ', v.integrationPoint
else:
print
# For each tensor component.
9–19
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
9–20
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
odb
rootAssembly
parts
= Container sections
= Singular object
materials
Model Data
steps
frames
fieldOutputs
historyRegions
point
historyOutputs
data
Results Data
sequence is a tuple of (frequency, value). In a modal domain analysis (domain=MODAL) the sequence
is a tuple of (mode, value).
In the analysis that generated the Abaqus/CAE Visualization module tutorial output database, the
user asked for the following history output:
At the rigid body reference point (Node 1000)
• U
• V
• A
The history output data can be retrieved from the HistoryRegion objects in the output database. The
tutorial output database contains HistoryRegion objects that relate to the rigid body reference point and
the integration points of the corner element as follows:
• ’Node PART-1-1.1000’
• ’Element PART-1-1.1 Int Point 1’
• ’Element PART-1-1.1 Int Point 2’
9–21
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
odb = openOdb(path='viewer_tutorial.odb')
step2 = odb.steps['Step-2']
region = step2.historyRegions['Node PART-1-1.1000']
u2Data = region.historyOutputs['U2'].data
dispFile = open('disp.dat','w')
for time, u2Disp in u2Data:
dispFile.write('%10.4E %10.4E\n' % (time, u2Disp))
dispFile.close()
The output in this example is a sequence of tuples containing the frame time and the displacement value.
The example uses nodal history data output. If the analysis requested history output from an element, the
output database would contain one HistoryRegion object and one HistoryPoint object for each integration
point.
9–22
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
if len(sys.argv) != 2:
print 'Error: you must supply the name \
of an odb on the command line'
sys.exit(1)
odbPath = sys.argv[1]
odb = openOdb(path=odbPath)
assembly = odb.rootAssembly
numNodes = numElements = 0
n = len(instance.nodes)
print 'Number of nodes of instance %s: %d' % (name, n)
numNodes = numNodes + n
print
print 'NODAL COORDINATES'
9–23
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
if instance.embeddedSpace == THREE_D:
print ' X Y Z'
for node in instance.nodes:
print node.coordinates
else:
print ' X Y'
for node in instance.nodes:
print node.coordinates
n = len(instance.elements)
print 'Number of elements of instance ', name, ': ', n
numElements = numElements + n
print
print 'Number of instances: ', len(assembly.instances)
print 'Total number of elements: ', numElements
print 'Total number of nodes: ', numNodes
9–24
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
# odbRead.py
# A script to read the Abaqus/CAE Visualization module tutorial
# output database and read displacement data from the node at
# the center of the hemispherical punch.
odb = openOdb(path='viewer_tutorial.odb')
lastFrame = odb.steps['Step-1'].frames[-1]
displacement = lastFrame.fieldOutputs['U']
center = odb.rootAssembly.instances['PART-1-1'].\
nodeSets['PUNCH']
centerDisplacement = displacement.getSubset(region=center)
for v in centerDisplacement.values:
print 'Position = ', v.position,'Type = ',v.type
print 'Node label = ', v.nodeLabel
print 'X displacement = ', v.data[0]
print 'Y displacement = ', v.data[1]
print 'Displacement magnitude =', v.magnitude
9–25
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
odb.close()
You can write your own data to an output database, and you can use Abaqus/CAE to view the data.
Writing to an output database is very similar to reading from an output database. When you open an
existing database, the Odb object contains all the objects found in the output database, such as instances,
steps, and field output data. In contrast, when you are writing to a new output database, these objects do
not exist. As a result you must use a constructor to create the objects. For example, you use the Part
constructor to create a Part object, the Instance constructor to create an OdbInstance object, and the
Step constructor to create an OdbStep object.
After you create an object, you use methods of the objects to enter or modify the data associated
with the object. For example, if you are creating an output database, you first create an Odb object.
You then use the Part constructor to create a part. After creating the part, you use the addNodes and
addElements methods of the Part object to add nodes and elements, respectively. Similarly, you use
the addData method of the FieldOutput object to add field output data to the output database. After
creating an output database, you should use the save method on the Odb object to save the output
database.
The example script in “Creating an output database,” Section 9.10.2, also illustrates how you can
write to an output database.
The following topics are covered:
9–26
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
9–27
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
part1.addNodes(nodeData=nodeData, nodeSetName='nset-1')
For a full description of the addNodes command, see “addNodes,” Section 33.20.4 of the Abaqus
Scripting Reference Manual.
After you have created nodes, you can use the NodeSetFromNodeLabels constructor to
create a node set from the node labels. For more information, see “NodeSetFromNodeLabels,”
Section 33.24.2 of the Abaqus Scripting Reference Manual.
Similarly, you use the addElements method to add elements to the part using a sequence
of element labels, element connectivity, and element type. You can also define an optional element
set and an optional section category. For example,
sCat = odb.SectionCategory(name='S5',
description='Five-Layered Shell')
spBot = sCat.SectionPoint(number=1,
description='Bottom')
spMid = sCat.SectionPoint(number=3,
description='Middle')
spTop = sCat.SectionPoint(number=5,
description='Top')
For a full description of the addElements command, see “addElements,” Section 33.20.2 of the
Abaqus Scripting Reference Manual.
odb.rootAssembly
You can create both instances and regions on the RootAssembly object.
Part instances
You use the Instance constructor to create part instances of the parts you have already defined
using the Part constructor. For example,
9–28
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
a = odb.rootAssembly
instance1 = a.Instance(name='part-1-1', object=part1)
You can also supply an optional local coordinate system that specifies the rotation and translation
of the part instance. You can add nodes and elements only to a part; you cannot add elements
and nodes to a part instance. As a result, you should create the nodes and elements that define the
geometry of a part before you instance the part. For a full description of the Instance command,
see “OdbInstance object,” Section 33.16 of the Abaqus Scripting Reference Manual.
Regions
Region commands are used to create sets from element labels, node labels, and element faces. You
can create a set on a part, part instance, or the root assembly. Node and element labels are unique
within an instance but not within the assembly. As a result, a set on the root assembly requires
the names of the part instances associated with the nodes and elements. You can also use region
commands to create surfaces. For example,
# An element set on an instance
eLabels = [9,99]
elementSet = instance1.ElementSetFromElementLabels(
name='elsetA',elementLabels=eLabels)
# A node set on the rootAssembly
nodeLabels = (5,11)
instanceName = 'part-1-1'
nodeSet = assembly.NodeSetFromNodeLabels(
name='nodesetRA',((instanceName,nodeLabels),))
The region commands are described in Chapter 44, “Region commands,” of the Abaqus
Scripting Reference Manual.
Materials
You use the Material object to list material properties.
Materials are stored in the materials repository under the Odb object.
To create an isotropic elastic material, with a Young’s modulus of 12000.0 and an effective
Poisson’s ratio of 0.3 in the output database:
materialName = "Elastic Material"
material_1 = odb.Material(name=materialName)
material_1.Elastic(type=ISOTROPIC,table=((12000,0.3),))
For more information, see Chapter 28, “Material commands,” of the Abaqus Scripting
Reference Manual.
Sections
You use the Section object to create sections and profiles.
9–29
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
Sections are stored in the sections repository under the Odb object.
The following code creates a homogeneous solid section object. A Material object must be
present before creating a Section object. An exception is thrown if the material does not exist.
Section assignments
You use the SectionAssignment object to assign sections and their associated material properties
to regions of the model. SectionAssignment objects are members of the Odb object. For a full
description of the assignSection method, see “assignSection,” Section 33.16.7 of the Abaqus
Scripting Reference Manual.
All Elements in an Abaqus analysis need to be associated with section and material properties.
Section assignments provide the relationship between elements in an Instance object and their
section properties. The section properties include the associated material name. To create an
element set and assign a section:
elLabels = (1,2)
elset = instance.ElementSetFromElementLabels(
name=materialName, elementLabels=elLabels)
instance.assignSection(region=elset,section=section)
step1 = odb.Step(name='step-1',
description='', domain=TIME, timePeriod=1.0)
9–30
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
The Step constructor has an optional previousStepName argument that specifies the step after
which this step must be inserted in the steps repository. For a full description of the Step
command, see “Step,” Section 33.25.1 of the Abaqus Scripting Reference Manual.
Frames
You use the Frame constructor to create a frame for field output. For example,
frame1 = step1.Frame(incrementNumber=1,
frameValue=0.1, description='')
For a full description of the Frame command, see “Frame,” Section 33.15.3 of the Abaqus Scripting
Reference Manual.
part1 = odb.Part(name='part-1',
embeddedSpace=THREE_D, type=DEFORMABLE_BODY)
a = odb.rootAssembly
instance1 = a.Instance(name='part-1-1', object=part1)
uField = frame1.FieldOutput(name='U',
description='Displacements', type=VECTOR)
nodeLabelData = (1, 2, 3, 4, 5, 6)
dispData = ((1,2,3),
(4,5,6),
(7,8,9),
9–31
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
(10,11,12),
(13, 14, 15),
(16,17,18))
uField.addData(position=NODAL, instance=instance1,
labels=nodeLabelData, data=dispData)
step1.setDefaultDeformedField(uField)
For a full description of the FieldOutput constructor, see “FieldOutput,” Section 33.6.1 of the Abaqus
Scripting Reference Manual.
The type argument to the FieldOutput constructor describes the type of the data—tensor, vector,
or scalar. The properties of the different tensor types are:
Full tensor
A tensor that has six components and three principal values. Full three-dimensional rotation of the
tensor is possible.
The valid components and invariants for the different data types are given in Table 9–1.
9–32
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
Table 9–1 Valid components and invariants for Abaqus data types.
For example, the following statements add element data to the FieldOutput object:
elementLabelData = (1, 2)
9–33
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
sField = frame1.FieldOutput(name='S',
description='Stress', type=TENSOR_3D_PLANAR,
componentLabels=('S11', 'S22', 'S33',
'S12'), validInvariants=(MISES,))
sField.addData(position=INTEGRATION_POINT,
sectionPoint=spTop, instance=instance1,
labels=elementLabelData, data=topData,
localCoordSystem=transform)
sField.addData(position=INTEGRATION_POINT,
sectionPoint=spBot, instance=instance1,
labels=elementLabelData, data=bottomData,
localCoordSystem=transform)
step1.setDefaultField(sField)
For a full description of the addData command, see “addData,” Section 33.6.3 of the Abaqus Scripting
Reference Manual.
As a convenience, localCoordSystem can be a single transform or a list of transforms. If
localCoordSystem is a single transform, it applies to all values. If localCoordSystem is a list of
transforms, the number of items in the list must match the number of data values.
9–34
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
a contour plot; for example, stress. Similarly, the default deformed field variable determines the variable
that distinguishes a deformed plot from an undeformed plot. Typically, you will use displacement for the
default deformed field variable; you cannot specify an invariant or a component. The default variable
settings apply for each frame in the step. For example, the following statements use the deformation ’U’
as the default setting for both field variable and deformed field variable settings during a particular step:
field=odb.steps['impact'].frames[1].fieldOutputs['U']
odb.steps['impact'].setDefaultField(field)
odb.steps['impact'].setDefaultDeformedField(field)
You can set a different default field variable and deformed field variable for different steps. You
will need to use a loop to set the defaults for each step. For example,
for step in odb.steps.values():
step.setDefaultField(field)
9–35
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
h1_u1 = h1.HistoryOutput(name='U1',
description='Displacement', type=SCALAR)
h1_rf1 = h1.HistoryOutput(name='RF1',
description='Reaction Force', type=SCALAR)
step2 = odb.Step(name='step-2',
description='', domain=TIME, timePeriod=1.0)
h2 = step2.HistoryRegion(name='my history',
description='my stuff', point=point1)
h2_u1 = h2.HistoryOutput(name='U1',
description='Displacement', type=SCALAR)
h2_rf1 = h2.HistoryOutput(name='RF1',
description='Reaction Force', type=SCALAR)
h1_u1.addData(frameValue=timeData, value=u1Data)
h1_rf1.addData(frameValue=timeData, value=rf1Data)
h2_u1.addData(frameValue=timeData, value=u1Data)
h2_rf1.addData(frameValue=timeData, value=rf1Data)
9–36
Abaqus ID:
Printed on:
COMPUTATIONS WITH Abaqus RESULTS
9–37
Abaqus ID:
Printed on:
COMPUTATIONS WITH Abaqus RESULTS
9–38
Abaqus ID:
Printed on:
COMPUTATIONS WITH Abaqus RESULTS
9–39
Abaqus ID:
Printed on:
COMPUTATIONS WITH Abaqus RESULTS
9–40
Abaqus ID:
Printed on:
COMPUTATIONS WITH Abaqus RESULTS
Note: Stress results for three-dimensional continuum elements transformed into a cylindrical
system would have the hoop stress in S22, which is consistent with the coordinate system axis
but inconsistent with the stress state for a three-dimensional axisymmetric elements having hoop
stress in S33.
9–41
Abaqus ID:
Printed on:
IMPROVING THE EFFICIENCY OF YOUR SCRIPTS
• When you are transforming a tensor, the location or integration point always takes into account the
deformation. The location of the coordinate system depends on the model, as follows:
– If the system is fixed, the coordinate system is fixed.
– If the system is model based, you must supply a displacement field that determines the
instantaneous location and orientation of the coordinate system.
• Abaqus will perform transformations of tensor results for shells, membranes, and planar elements
as rotations of results about the element normal at the element result location. The element normal
is the normal computed for the frame associated with the field by Abaqus, and you cannot redefine
the normal. Abaqus defines the location of the results location from the nodal locations. You
specify optional arguments if you want to use the deformed nodal locations to transform results.
For rectangular, cylindrical, and spherical coordinate systems the second component direction for
the transformed results will be determined by one of the following:
– The Y-axis in a rectangular coordinate system.
– The -axis in a cylindrical coordinate system.
– The -axis in a spherical coordinate system.
– A user-specified datum axis projected onto the element plane.
If the coordinate system used for projection and the element normal have an angle less than the
specified tolerance (the default is 30°), Abaqus will use the next axis and generate a warning.
If you are accessing large amounts of data from an output database, you should be aware of potential
inefficiencies in your script and techniques that will help to speed up your scripts.
• “Creating objects to hold temporary variables,” Section 9.9.1
9–42
Abaqus ID:
Printed on:
IMPROVING THE EFFICIENCY OF YOUR SCRIPTS
count = 0
for v in stressField.values:
if v.mises > stressCap:
if v.integrationPoint:
print 'Element label = ', v.elementLabel, \
'Integration Point = ', v.integrationPoint
else:
print 'Element label = ', v.elementLabel
for component in strainField.values[count].data:
print '%-10.5f' % component,
print
count = count + 1
In this example every time the script accesses a strain component from strainField.value, Abaqus
must reconstruct the sequence of FieldValue objects. This reconstruction could result in a significant
performance degradation, particularly for a large model.
A slight change in the script greatly improves its performance, as shown in the following example:
stressField = frame.fieldOutputs['MISES']
strainFieldValues = frame.fieldOutputs['LE'].values
count = 0
for v in stressField.values:
if v.mises > stressCap:
if v.integrationPoint:
print 'Element label = ', v.elementLabel, \
'Integration Point = ', v.integrationPoint
else:
print 'Element label = ', v.elementLabel
for component in strainFieldValues[count].data:
print '%-10.5f' % component,
print
count = count + 1
for i in range(len(odb.steps[name].frames)-1):
9–43
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
frame[i] = odb.steps[name].frames[i]
Creating a temporary variable to hold the frame repository provides the same functionality and speeds
up the process:
frameRepository = odb.steps[name].frames
for i in range(len(frameRepository)-1):
frame[i] = frameRepository[i]
Such a potential loss of performance will not be a problem when accessing a load case frame. Accessing
a load case frame does not result in the creation of a frame repository and, thus, does not suffer from a
corresponding loss of performance.
The following examples illustrate how you use the output database commands to access data from an
output database:
9–44
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
9–45
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
def rightTrim(input,suffix):
if (input.find(suffix) == -1):
input = input + suffix
return input
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
def getMaxMises(odbName,elsetName):
""" Print max mises location and value given odbName
and elset(optional)
"""
elset = elemset = None
region = "over the entire model"
""" Open the output database """
odb = openOdb(odbName)
assembly = odb.rootAssembly
9–46
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
""" Close the output database before exiting the program """
odb.close()
#==================================================================
# S T A R T
#
if __name__ == '__main__':
odbName = None
elsetName = None
argList = argv
argc = len(argList)
i=0
while (i < argc):
if (argList[i][:2] == "-o"):
i += 1
name = argList[i]
9–47
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
odbName = rightTrim(name,".odb")
elif (argList[i][:2] == "-e"):
i += 1
elsetName = argList[i]
elif (argList[i][:2] == "-h"):
print __doc__
exit(0)
i += 1
if not (odbName):
print ' **ERROR** output database name is not provided'
print __doc__
exit(1)
getMaxMises(odbName,elsetName)
"""odbWrite.py
Script to create an output database and add model,
field, and history data. The script also reads
history data, performs an operation on the data, and writes
the result back to the output database.
usage: abaqus python odbWrite.py
"""
from odbAccess import *
from odbMaterial import *
9–48
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
def createODB():
9–49
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
(1, 1,0,0),
(2, 2,0,0),
(3, 2,1,0.1),
(4, 1,1,0.1),
(5, 2,-1,-0.1),
(6, 1,-1,-0.1),
)
part1.addNodes(nodeData=nodeData,
nodeSetName='nset-1')
elementData = (
(1, 1,2,3,4),
(2, 6,5,2,1),
)
part1.addElements(elementData=elementData, type='S4',
elementSetName='eset-1', sectionCategory=sCat)
# Field data:
step1 = odb.Step(name='step-1',
description='first analysis step',
domain=TIME, timePeriod=1.0)
analysisTime=0.1
frame1 = step1.Frame(incrementNumber=1,
frameValue=analysisTime,
description=\
'results frame for time '+str(analysisTime))
9–50
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
uField = frame1.FieldOutput(name='U',
description='Displacements', type=VECTOR)
nodeLabelData = (1, 2, 3, 4, 5, 6)
dispData = (
(1,2,3),
(4,5,6),
(7,8,9),
(10,11,12),
(13, 14, 15),
(16,17,18)
)
uField.addData(position=NODAL, instance=instance1,
labels=nodeLabelData,
data=dispData)
step1.setDefaultDeformedField(uField)
elementLabelData = (1, 2)
topData = (
(1.,2.,3.,4.),
(1.,2.,3.,4.),
(1.,2.,3.,4.),
(1.,2.,3.,4.),
(1.,2.,3.,4.),
(1.,2.,3.,4.),
(1.,2.,3.,4.),
9–51
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
(1.,2.,3.,4.),
)
bottomData = (
(1.,2.,3.,4.),
(1.,2.,3.,4.),
(1.,2.,3.,4.),
(1.,2.,3.,4.),
(1.,2.,3.,4.),
(1.,2.,3.,4.),
(1.,2.,3.,4.),
(1.,2.,3.,4.),
)
transform = (
(1.,0.,0.),
(0.,1.,0.),
(0.,0.,1.)
)
sField = frame1.FieldOutput(name='S',
description='Stress', type=TENSOR_3D_PLANAR,
componentLabels=('S11', 'S22', 'S33','S12'),
validInvariants=(MISES,))
sField.addData(position=INTEGRATION_POINT,
sectionPoint=spTop, instance=instance1,
labels=elementLabelData, data=topData,
localCoordSystem=transform)
sField.addData(position=INTEGRATION_POINT,
sectionPoint=spBot, instance=instance1,
labels=elementLabelData, data=bottomData,
localCoordSystem=transform)
step1.setDefaultField(sField)
# History data:
9–52
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
hRegionStep1 = step1.HistoryRegion(name='historyNode0',
description='Displacement and reaction force',
point=instance1.nodes[0])
hOutputStep1U1 = hRegionStep1.HistoryOutput(name='U1',
description='Displacement', type=SCALAR)
hOutputStep1Rf1 = hRegionStep1.HistoryOutput(name='RF1',
description='Reaction Force', type=SCALAR)
hOutputStep1U1.addData(frameValue=timeData1,
value=u1Data)
hOutputStep1Rf1.addData(frameValue=timeData1,
value=rf1Data)
9–53
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
hOutputStep2U1.addData(frameValue=timeData2,
value=u1Data)
hOutputStep2Rf1.addData(frameValue=timeData2,
value=rf1Data)
hOutputStep2sumU1 = hRegionStep2.HistoryOutput(
name='squareU1',
description='Square of displacements',
type=SCALAR)
hOutputStep2sumU1.addData(data=u1SquaredFromStep2.data)
odb.save()
odb.close()
if __name__ == "__main__":
createODB()
9–54
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
import sys
from odbAccess import *
from types import IntType
step = odb.steps.values()[0]
9–55
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
try:
coords = step.frames[0].fieldOutputs['COORD']
except:
err = "The analysis must include a field output request \
for variable COORD."
print err
sys.exit(1)
factors = []
for mode in modes:
try:
frame = step.frames[mode]
except IndexError:
print 'Input error: mode %s does not exist' % mode
sys.exit(1)
factors.append(float(raw_input(
'Enter imperfection factor for mode %s: '% mode)))
coords = coords + factors[-1] * frame.fieldOutputs['U']
9–56
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
field1 = odb.steps['LC1'].frames[1].fieldOutputs['U']
field2 = odb.steps['LC2'].frames[1].fieldOutputs['U']
newStep = odb.Step(name='user',
description='user defined results', domain= TIME, timePeriod=0)
9–57
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
odb.save()
#
# fieldValue operators example problem:
#
# sum and average stress field values in a region
#
#
# get field
#
odb = openOdb(path='sumRegionFieldValue.odb')
endSet = odb.rootAssembly.elementSets['END1']
field = odb.steps.values()[-1].frames[-1].fieldOutputs['S']
subField = field.getSubset(region=endSet)
9–58
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
# sum values
#
sum = 0
for val in subField.values:
sum = sum + val
ave = sum / len(subField.values)
#
# print results
#
9–59
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
# allowable displacement.
#
#
# get historyRegion for the node in nodeSet TIP
#
odb = openOdb(path='compDispMagHist.odb')
endSet = odb.rootAssembly.instances['BEAM-1-1'].nodeSets['TIP']
histPoint = HistoryPoint(node=endSet.nodes[0])
tipHistories = odb.steps['Step-2'].getHistoryRegion(
point=histPoint)
#
# Compute and scale magnitude.
#
maxAllowableDisp = 5.0
sum = 0
componentLabels = ('U1', 'U2', 'U3')
for name in componentLabels:
sum = sum + power(tipHistories.historyOutputs[name], 2.0)
sum = sqrt(sum) / maxAllowableDisp
#
# Print magnitude.
#
9–60
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
• Retrieves the information describing the new load combination from the command line.
• Retrieves the frames for each load case.
• Computes the new stresses and displacements.
• Saves data computed to the output database as a new load combination.
The command line arguments provide the following:
• odbName: The output database file name.
• stepName: The name of the step containing the load cases.
• loadCaseNames: The load case names.
• scaling: The scale factors to apply to each load case.
Use the following command to retrieve the example script:
abaqus fetch job=createLoadComb
The fetch command also retrieves an input file that you can use to generate an output database that can
be read by the example script.
import types
from odbAccess import *
loadCaseNames = eval(raw_input( \
'Enter new load case as: \
[\'loadCase1Name\', ..., \'loadCaseNName\']'))
if type(loadCaseNames) == types.TupleType:
loadCaseNames = list(loadCaseNames)
lcName = raw_input('Enter new load case name')
scaling = eval(raw_input( \
'Enter new load case as:(scaleFactor1, .., scaleFactorN)'))
odb = openOdb(odbName)
step = odb.steps[stepName]
9–61
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
odb.save()
odb.close()
9–62
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
odb.save()
odb.close()
9–63
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
odb = openOdb(odbName)
step = odb.steps[stepName]
frame = step.frames[frameNo]
displacement = frame.fieldOutputs['U']
coordSys = odb.rootAssembly.DatumCsysByThreePoints(name='cylC',
coordSysType=CYLINDRICAL, origin=(0,0,0),
point1=(1.0, 0.0, 0), point2=(0.0, 0.0, 1.0) )
cylindricalDisp = displacement.getTransformedField(
datumCsys=coordSys)
radialDisp = cylindricalDisp.getScalarField(componentLabel='U1')
sum = 0.0
for val in radialDisp.values:
sum = sum + val.data
aveDisp = sum / len(radialDisp.values)
# Compute distortion.
9–64
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
frame.FieldOutput(field=radialDisp)
fieldDescription = 'Distortion ( \
average radial displacement = ' + str(aveDisp) + ')'
frame.FieldOutput(name='Distortion',
description=fieldDescription, field=distortion)
odb.save()
odb.close()
• An output database generated by the two-dimensional cross-section analysis. The script reads cross-
section data, including the out-of-plane warping function, from this output database.
• An output database generated by the beam analysis. The script reads generalized section strains
(SE) from this output database.
You must run the script from Abaqus/CAE by selecting File→Run Script from the main menu bar. The
script uses getInputs to display a dialog box that prompts you for the name of the output databases
generated by the two-dimensional cross-section analysis and by the beam analysis. The names are case-
insensitive, and you can omit the .odb file suffix. The files must be in the local directory. The dialog
box also prompts you for the following:
9–65
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
9–66
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
where is the acoustic pressure at a distance from the reference point, is the wave number, and
is the acoustic far-field pressure. The acoustic pressure decibel value is defined as
where is the magnitude of the acoustic pressure at a point, is the root mean square acoustic
pressure, and is the decibel reference value given as user input. The far-field pressure decibel
value is defined in the same manner as , using the same reference value ( ).
Note: If (in SI units), corresponds to
The script also calculates the far-field acoustic intensity, which is defined as
where is the far-field rms pressure, is the fluid density, and c is the speed of sound in the
medium.
Before you execute the script, you must run a direct-solution, steady-state dynamics acoustics
analysis that includes three-dimensional acoustic infinite elements (ACIN3D3, ACIN3D4, ACIN3D6,
and ACIN3D8). In addition, the output database must contain results for the following output variables:
• INFN, the acoustic infinite element normal vector.
• INFR, the acoustic infinite element “radius,” used in the coordinate map for these elements.
9–67
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
9–68
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
abaqus job=singleTriangularElementModel
The results from running the script twice using the single triangular-element model, changing the radius
of the sphere, and writing the data back to the original output database are shown in Figure 9–6.
9–69
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
POR
+1.000e−03
+9.294e−04
+8.587e−04
+7.881e−04
+7.174e−04
+6.468e−04
+5.761e−04
+5.055e−04
+4.348e−04
+3.642e−04
+2.935e−04
+2.229e−04
+1.523e−04
This model simulates the response of a sphere in “breathing” mode (a uniform radial
expansion/compression mode). The model consists of one triangular ACIN3D3 element. Each node of
the element is placed on a coordinate axis at a distance of 1.0 from the origin that serves as the reference
point for the infinite element. The acoustic material properties do not have physical significance; the
values used are for convenience only. The loading consists of applying an in-phase pressure boundary
condition to all the nodes. Under this loading and geometry, the model behaves as a spherical source
(an acoustic monopole) radiating in the radial direction only. The acoustic pressure, , and the acoustic
far-field pressure, , at a distance from the center of the sphere are
and
where is the known acoustic pressure at some reference distance and is the wave number.
For this single-element example, you should enter a value of 1.0 for the speed of sound; thus,
, where is the frequency in Hz. in this model is 1, and is 0.001. The equations for the
acoustic pressure, , and the acoustic far-field pressure, , reduce to
and
9–70
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
9–71
Abaqus ID:
Printed on:
EXAMPLE SCRIPTS THAT ACCESS DATA FROM AN OUTPUT DATABASE
Before executing the script, run an analysis that creates an output database file containing the
appropriate output. This analysis includes, for example, output for the elements and the integration
point coordinates of the elements. Execute the script using the following command:
The script prompts for other information, such as the desired postprocessing option, part name, etc. The
script processes the data and produces a text file or a new output database that contains the information
required to visualize the elbow element results.
“Elastic-plastic collapse of a thin-walled elbow under in-plane bending and internal pressure,”
Section 1.1.2 of the Abaqus Example Problems Manual, contains several figures that can be created
with the aid of this program.
9–72
Abaqus ID:
Printed on:
WHAT DO YOU NEED TO ACCESS THE OUTPUT DATABASE?
10.1 Overview
The C++ interface to an output database is related closely to the Abaqus Scripting Interface. Disparities
between the two interfaces are due to fundamental differences in the programming languages. The C++
interface is intended for users with high-performance requirements; others are encouraged to use the
Abaqus Scripting Interface.
A working knowledge of the C++ programming language is assumed.
10–1
Abaqus ID:
Printed on:
Abaqus SCRIPTING INTERFACE DOCUMENTATION STYLE
This section describes the style that is used to describe a command in the Abaqus Scripting Reference
Manual. You may want to refer to the Abaqus Scripting Reference Manual while you read this section
and compare the style of a documented command with the descriptions provided here. The following
topics are covered:
• “How the commands are ordered,” Section 10.3.1
• “Access,” Section 10.3.2
• “Path,” Section 10.3.3
• “Prototype,” Section 10.3.4
• “Return value,” Section 10.3.5
10–2
Abaqus ID:
Printed on:
Abaqus SCRIPTING INTERFACE DOCUMENTATION STYLE
10.3.2 Access
The description of each object in the Abaqus Scripting Reference Manual begins with a section that
describes how you access an instance of the object.
The following is the access description for the Part object:
odb.parts()[name]
The access description specifies where instances of the object are located in the data model. The Part
object can accordingly be accessed as:
odb.steps()[name].frames(i).fieldOutputs()[name]
The following statements show how you use the object described by this Access description:
10–3
Abaqus ID:
Printed on:
Abaqus SCRIPTING INTERFACE DOCUMENTATION STYLE
10.3.3 Path
A method that creates an object is called a “constructor.” The Abaqus C++ API uses the convention that
constructors begin with an uppercase character. In contrast, methods that operate on an object begin with
a lowercase character. The description of each constructor in the Abaqus Scripting Reference Manual
includes a path to the command. For example, the following describes the path to the Part constructor:
odb.Part
Some constructors include more than one path. For example, you can create a nodeSet that is associated
with either a Part object or the RootAssembly object, and each path is listed.
odb.parts()[name].NodeSet
odb.rootAssembly().NodeSet
The path is not listed if the method is not a constructor.
If you are using the Abaqus C++ API to read data from an output database, the objects exist when
you open the output database, and you do not have to use constructors to create them. However, if you
are creating or writing to an output database, you may need to use constructors to create new objects,
such as part instances and steps. The documentation describes the path to the constructors that create
objects in an output database.
For example, the Path description for the FieldOutput constructor is
odb.steps()[name].frames(i).FieldOutput
The following statement creates a FieldOutput object:
odb_StepContainer stepCon = odb.steps();
odb_Step step = stepCon["Side load"];
odb_SequenceFrame frameSeq = step.frames();
odb_Frame frame = frameSeq.Get( frameSeq.Size() -1 );
odb_FieldOutput& myFieldOutput = frame.FieldOutput("S",
"stress", odb_Enum::TENSOR_3D_FULL);
10.3.4 Prototype
Chapter 60, “Odb commands,” of the Abaqus Scripting Reference Manual, contains a prototype section
for each C++ command. The prototype provides the type returned by the command, the name of the
command, and a list of all its arguments along with the type of each argument. Required arguments
appear first in the list followed by default arguments along with their default value. For example, the
Frame constructor is given as
odb_Frame Frame(int incrementNumber, float frameValue,
10–4
Abaqus ID:
Printed on:
HOW THE OBJECT MODEL FOR THE OUTPUT DATABASE RELATES TO COMMANDS
10.4 How the object model for the output database relates to
commands
You need to understand the object model for the output database both to read data from it and to write data
to it. An object model describes the relationship between objects. The object model for the Abaqus/CAE
model is described in “The Abaqus object model,” Section 6.1.
For example, consider the object model for field output data shown in Figure 10–1. The Odb object
at the top of the figure is created when you issue the command to open or create an output database. As
you move down the object model, an OdbStep object is a member of the Odb object; similarly, a Frame
object is a member of the OdbStep object. The FieldOutput object has two members—fieldValue and
fieldLocation.
The object model translates directly to the structure of an Abaqus C++ API command. For example,
the following command refers to a Frame object in the sequence of frames contained in an OdbStep
object:
odb.steps()["10 hz"].frames(3);
Similarly, the following command refers to the sequence of field data contained in a FieldOutput object.
odb.steps()["10 hz"].frames.get(3).
10–5
Abaqus ID:
Printed on:
OBJECT MODEL FOR THE OUTPUT DATABASE
odb
rootAssembly
parts
sections
= Container
= Singular object materials
Model Data
steps
frames
fieldOutputs
historyRegions
point
historyOutputs
Results Data
fieldOutputs()["U"].values(47);
You use commands to access objects by stepping through the hierarchy of objects in the object model.
The Access, Path, and Prototype descriptions in Chapter 60, “Odb commands,” of the Abaqus
Scripting Reference Manual describe the interface definition of the command. The interface definition
of the command reflects the hierarchy of objects in the object model.
An output database generated from an Abaqus analysis contains both model and results data as shown
in Figure 10–1.
Model data
Model data describe the parts and part instances that make up the root assembly; for example, nodal
coordinates, set definitions, and element types. Model data are explained in more detail in “Model
data,” Section 10.5.1.
10–6
Abaqus ID:
Printed on:
OBJECT MODEL FOR THE OUTPUT DATABASE
Results data
Results data describe the results of your analysis; for example, stresses, strains, and displacements.
You use output requests to configure the contents of the results data. Results data can be either
field output data or history output data; for a more detailed explanation, see “Results data,”
Section 10.5.2.
Note: For a description of object models, see “An overview of the Abaqus object model,” Section 6.1.1.
You can find more information on the format of the output database in “Output to the output database,”
Section 4.1.3 of the Abaqus Analysis User’s Manual.
Part instances
A part instance is a usage of a part within the assembly. All characteristics (such as mesh and
section definitions) defined for a part become characteristics for each instance of that part—they
are inherited by the part instances. Each part instance is positioned independently within the root
assembly.
10–7
Abaqus ID:
Printed on:
OBJECT MODEL FOR THE OUTPUT DATABASE
Materials
Materials contain material models comprised of one or more material property definitions. The same
material models may be used repeatedly within a model; each component that uses the same material
model shares identical material properties. Many materials may exist within a model database, but
only the materials that are used in the assembly are copied to the output database.
Sections
Sections add the properties that are necessary to define completely the geometric and material
properties of an element. Various element types require different section types to complete their
definitions. For example, shell elements in a composite part require a section that provides a
thickness, multiple material models, and an orientation for each material model; all these pieces
combine to complete the composite shell element definition. Like materials, only those sections
that are used in the assembly are copied to the output database.
Section assignments
Section assignments link section definitions to the regions of part instances. Section assignments
in the output database maintain this association. Sections are assigned to each part in a model, and
the section assignments are propagated to each instance of that part.
Section categories
You use section categories to group the regions of the model that use the same section definitions;
for example, the regions that use a shell section with five section points. Within a section category,
you use the section points to identify the location of results; for example, you can associate section
point 1 with the top surface of a shell and section point 5 with the bottom surface.
Rigid bodies
You use rigid bodies to define a collection of nodes, elements, and/or surfaces whose motion is
governed by the motion of a single node, called the rigid body reference node.
Pretension Sections
Pretension sections are used to associate a pre-tension node with a pre-tension section. The pre-
tension section can be defined using a surface for continuum elements or using an element for truss
or beam elements.
10–8
Abaqus ID:
Printed on:
OBJECT MODEL FOR THE OUTPUT DATABASE
Interactions
Interactions are used to define contact between surfaces in an analysis. Only contact interactions
defined using contact pairs are written to the output database.
Interaction properties
Interaction properties define the physical behavior of surfaces involved in an interaction. Only
tangential friction behavior is written to the output database.
steps
Step
frames
Frame
fieldOutputs
historyRegions
HistoryRegion
point
historyOutputs
Results Data
10–9
Abaqus ID:
Printed on:
OBJECT MODEL FOR THE OUTPUT DATABASE
Frames
Each step contains a sequence of frames, where each increment of the analysis that resulted in output
to the output database is called a frame. In a frequency or buckling analysis each eigenmode is stored
as a separate frame. Similarly, in a steady-state harmonic response analysis each frequency is stored
as a separate frame.
Field output
Field output is intended for infrequent requests for a large portion of the model and can be used
to generate contour plots, animations, symbol plots, and displaced shape plots in the Visualization
module of Abaqus/CAE. You can also use field output to generate an X–Y data plot. Only complete
sets of basic variables (for example, all the stress or strain components) can be requested as field
output. Field output is composed of a “cloud of data values” (e.g., stress tensors at each integration
point for all elements). Each data value has a location, type, and value. You use the regions defined
in the model data, such as an element set, to access subsets of the field output data. Figure 10–3
shows the field output data object model within an output database.
History output
History output is output defined for a single point or for values calculated for a portion of the model
as a whole, such as energy. History output is intended for relatively frequent output requests for
small portions of the model and can be displayed in the form of X–Y data plots in the Visualization
module of Abaqus/CAE. Individual variables (such as a particular stress component) can be
requested.
Depending on the type of output expected, a HistoryRegion object can be defined for one of
the following:
• a node
• an integration point
• a region
• the whole model
The output from all history requests that relate to a particular point or region is then collected in
one HistoryRegion object. Figure 10–4 shows the history output data object model within an output
database.
10–10
Abaqus ID:
Printed on:
COMPILING AND LINKING YOUR C++ SOURCE CODE
odb
rootAssembly
parts
sections
materials
Model Data
= Container
= Singular object
steps
frames
fieldOutputs
fieldValues
historyRegions
Results Data
odb
rootAssembly
parts
= Container sections
= Singular object
materials
Model Data
steps
frames
fieldOutputs
historyRegions
point
historyOutputs
data
Results Data
10–11
Abaqus ID:
Printed on:
ACCESSING THE C++ INTERFACE FROM AN EXISTING APPLICATION
the programs cannot contain a C++ main routine. Instead, the programs must begin with a C++ function
called ABQmain.
#include <odb_API.h>
This section provides information that may be helpful to users who need to access results in an output
database from an existing application. Most users should find that the abaqus make utility is sufficient
for their postprocessing needs. Since linking and executing with dynamically linked runtime libraries
is highly system dependent, the information in this section is intended for users who have an advanced
working knowledge of compilation and linking with runtime libraries.
It is important to ensure that the compiler used to compile and link the existing application
is consistent with the compilers used to generate the Abaqus release. The “System Requirements”
document lists the name and version of the compiler used for the Abaqus release on each supported
platform. You can access this document through the System Information section of the Support
page at www.simulia.com. You can also find information on compiling and linking with the
C++ interface to an output database in the Dassault Systèmes DSX.ECO Knowledge Base at
www.3ds.com/support/knowledge-base or the SIMULIA Online Support System, which is accessible
through the My Support page at www.simulia.com.
The following topics are covered in this section:
• “Initializing the C++ interface,” Section 10.7.1
• “Link library location,” Section 10.7.2
• “Runtime library location,” Section 10.7.3
• “Header file location,” Section 10.7.4
10–12
Abaqus ID:
Printed on:
ACCESSING THE C++ INTERFACE FROM AN EXISTING APPLICATION
odb_initializeAPI();
This call is generated automatically when the abaqus make utility is run but must be included in any
application that is not compiled and linked using the abaqus make utility. After all calls to the C++
interface have been completed, the interface may be deactivated by including a call to
odb_finalizeAPI();
If the finalization call is not made explicitly, the finalize routine will be called automatically when
the application exits.
Windows
abaqus_dir\exec\lbr
where abaqus_dir is the name of the directory in which Abaqus is installed. To determine the location
of abaqus_dir at your site, type abaqus whereami at an operating system prompt.
During linking, the ABQodb library and several other libraries shipped with the Abaqus release are
used to resolve all the functions available in the interface to the output database. The command used
by Abaqus to link runtime libraries (for example, for user subroutines) is available through the Abaqus
environment variable link_sl. Additional information about linking with the Abaqus libraries, including
the names of all libraries which must be specified as part of the link command, may be obtained by
running the abaqus make utility in verbose mode with a verbosity level of 3.
Windows
abaqus_dir\exec\lbr
where abaqus_dir is the name of the directory in which Abaqus is installed. To determine the location
of abaqus_dir at your site, type abaqus whereami at an operating system prompt.
10–13
Abaqus ID:
Printed on:
THE Abaqus C++ API ARCHITECTURE
The correct path to the Abaqus runtime libraries must be specified prior to starting the user
application. The runtime library path is typically set using the system environment variable
LD_LIBRARY_PATH, but the method used to set the path may vary depending on your operating
system configuration. The ABQodb library and several utility libraries resolve all the functions
available in the interface to the output database, as described in “Link library location,” Section 10.7.2.
At runtime these libraries depend on many of the underlying Abaqus libraries. As a result, if you do not
define the correct runtime library path, your application will not run.
UNIX
abaqus_dir/include
Windows
abaqus_dir\include
where abaqus_dir is the name of the directory in which Abaqus is installed. To determine the location
of abaqus_dir at your site, type abaqus whereami at an operating system prompt.
Only odb_API.h must be included to access the C++ interface, but the path to the header files
must be provided during compilation.
This section describes the architecture of the Abaqus C++ interface to an output database. The output
database is an object-oriented database, which means that the data are held by “objects” (C++ classes)
that have certain behavior (C++ methods). The methods of an object in the database allow access to
and manipulation of the data held by the object. The data members of an object can be either primitives
(integer, floating point, string) or other objects.
The following topics are covered:
10–14
Abaqus ID:
Printed on:
THE Abaqus C++ API ARCHITECTURE
10.8.2 Constructors
A constructor is a method that creates an object. The Abaqus C++ API uses the following three types of
constructors:
Constructors for nonpersistent objects
Constructors for nonpersistent objects are the standard C++ constructors. For example,
odb_String partName("New_Part");
10–15
Abaqus ID:
Printed on:
UTILITY INTERFACE
The Abaqus C++ API provides a set of utilities that allow a user to access certain commonly used
functionality (such as strings, sequences (lists), and repositories) quickly and easily using a set of
supported and maintained interfaces.
The following topics are covered:
• “Utility interface classes,” Section 10.9.1
• “Utility interface examples,” Section 10.9.2
Sequences
An odb_Sequence class is a container used to hold an ordered list of objects of a specific type. Data
can be appended and retrieved from the sequence.
The following odb_Sequence objects are provided to store integer, float, and enumeration data:
• odb_SequenceInt
• odb_SequenceFloat
• odb_SequenceString
• odb_SequenceInvariant
• odb_SequenceElementFace
10–16
Abaqus ID:
Printed on:
UTILITY INTERFACE
Repositories
Repositories are provided to store objects retrieved by name. Both the repositories and the content of
the repositories are created by the API; the user can only retrieve objects from repositories. Iterators
are provided to navigate the repositories.
The following Abaqus repositories are provided:
• odb_PartRepository
• odb_FieldOutputRepository
• odb_SectionCategoryRepository
• odb_HistoryRegionRepository
• odb_SetRepository
• odb_HistoryOutputRepository
• odb_StepRepository
• odb_InstanceRepository
More detail on these interface utility objects can be found in Chapter 60, “Odb commands,” of the Abaqus
Scripting Reference Manual.
10–17
Abaqus ID:
Printed on:
UTILITY INTERFACE
Strings
odb_String type = stressField.baseElementTypes()[0];
odb_String elementType =
odb_String("Element type is ") + type;
cout << elementType.CStr() << endl;
Sequences
odb_Set& mySurface = rootAssy.surfaces()["TARGET"];
const odb_String instanceName = "PART-1-1";
const odb_SequenceElementFace allFaces =
mySurface.faces(instanceName);
odb_SequenceSequenceElementFace newFaces;
int allFaces_size = allFaces.size();
for (int i=0; i<allFaces_size; i++) {
const odb_SequenceElementFace fList = allFaces[i];
odb_SequenceElementFace newList;
int fList_size = fList.size();
for (int j=0; j<fList_size; j++) {
const odb_Enum::odb_ElementFaceEnum face = fList[j];
newList.append(face);
}
newFaces.append(newList);
}
Repositories
odb_StepRepository stepCon = odb.steps();
odb_StepRepositoryIT iter (stepCon);
for (iter.first(); !iter.isDone(); iter.next()) {
cout << "step name : " << iter.currentKey().CStr() << endl;
const odb_Step& step = iter.currentValue();
cout << "step description : " << step.description().CStr();
cout << endl;
}
10–18
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
The following sections describe how you use Abaqus C++ API commands to read data from an output
database. The following topics are covered:
• “The Abaqus/CAE Visualization module tutorial output database,” Section 10.10.1
• “Making the Odb commands available,” Section 10.10.2
• “Opening an output database,” Section 10.10.3
• “Reading model data,” Section 10.10.4
• “Reading results data,” Section 10.10.5
• “Reading field output data,” Section 10.10.6
• “Using bulk data access to an output database,” Section 10.10.7
• “Using regions to read a subset of field output data,” Section 10.10.8
• “Reading history output data,” Section 10.10.9
• “An example of reading field data from an output database,” Section 10.10.10
10–19
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
#include <odb_API.h>
To make the material and section Odb commands available to your program, you also need to include
their output database classes:
#include <odb_MaterialTypes.h>
#include <odb_SectionTypes.h>
After you open the output database, you can access its contents using the methods and members of the
Odb object returned by the openOdb method. In the above example the Odb object is referred to by the
variable odb. For a full description of the openOdb command, see “openOdb,” Section 60.33.5 of the
Abaqus Scripting Reference Manual.
10–20
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
Part instances
Part instances are stored in the instance repository under the OdbAssembly object. The
following statements display the repository keys of the part instances in the tutorial output
database:
odb_InstanceRepositoryIT instIter(rootAssy.instances());
for (instIter.first(); !instIter.isDone(); instIter.next())
cout << instIter.currentKey().CStr() << endl;
The output database contains only one part instance, and the resulting output is
PART-1-1
From a part instance or part you can retrieve the node and element information as follows:
{
odb_Instance& instance1 =
rootAssy.instances()["PART-1-1"];
odb_Enum::odb_DimensionEnum instanceType =
instance1.embeddedSpace();
const odb_SequenceNode& nodeList = instance1.nodes();
int nodeListSize = nodeList.size();
if (instanceType == odb_Enum::THREE_D) {
for (int n=0; n<nodeListSize; n++) {
const odb_Node node = nodeList[n];
int nodeLabel = node.label();
const float* const coord = node.coordinates();
cout << "Xcoord: " << coord[0] << " , Ycoord: "
<< coord[1] << " , Zcoord: " << coord[2] << endl;
}
}
else if((instanceType == odb_Enum::TWO_D_PLANAR) ||
(instanceType == odb_Enum::AXISYMMETRIC)) {
for (int n=0; n<nodeListSize; n++) {
const odb_Node node = nodeList[n];
int nodeLabel = node.label();
const float* const coord = node.coordinates();
cout << "Xcoord: " << coord[0] << " , Ycoord: "
<< coord[1] << endl;
}
}
10–21
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
Regions
Regions in the output database are OdbSet objects. Regions refer to the part and assembly sets stored
in the output database. A part set refers to elements or nodes in an individual part and appears in
each instance of the part in the assembly. An assembly set refers to the elements or nodes in part
instances in the assembly. A region can be one of the following:
• A node set
• An element set
• A surface
For example, the following statement displays the node sets in the OdbAssembly object:
cout << "Node set keys:" << endl;
odb_SetRepositoryIT setIter( rootAssy.nodeSets() );
for (setIter.first(); !setIter.isDone(); setIter.next())
cout << setIter.currentKey().CStr() << endl;
10–22
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
The following statements display the node sets and the element sets in the PART-1-1 part instance:
{
odb_InstanceRepository& iCon =
odb.rootAssembly().instances();
odb_Instance& instance = iCon["PART-1-1"];
The type of the object to which topNodeSet refers is OdbSet. After you create a variable that
refers to a region, you can use the variable to refer to a subset of field output data, as described in
“Using regions to read a subset of field output data,” Section 10.10.8.
10–23
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
// surface information
odb_Set& surface = instance.surfaces()["IMPACTOR"];
const odb_SequenceElement& elementList =
surface.elements();
const odb_SequenceElementFace& faces =
surface.faces();
The set information in an assembly set is keyed by instance name and can be accessed using
the following:
// assembly surface information
odb_Set& aSurface = rootAssy.surfaces()["TARGET"];
odb_SequenceString instanceNames =
aSurface.instanceNames();
int totalNames = instanceNames.size();
for (int name=0; name<totalNames; name++) {
const odb_String& iName = instanceNames[name];
const odb_SequenceElement& els =
aSurface.elements(iName);
const odb_SequenceElementFace& face =
aSurface.faces(iName);
10–24
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
Materials
Materials are stored in the materials repository under the Odb object.
Extend the Material commands available to the Odb object using the following statement:
odb_MaterialApi materialApi;
odb.extendApi(odb_Enum::odb_MATERIAL,materialApi);
Access the materials repository using the command:
odb_MaterialContainer& materialContainer = materialApi.materials();
odb_MaterialContainerIT matIT(materialContainer);
for (matIT.first(); !matIT.isDone(); matIT.next()) {
cout << "Material Name : " << matIT.currentKey().CStr() << endl;
const odb_Material& myMaterial = matIT.currentValue();
Some Material definitions have suboptions. For example, to access the smoothing type used
for biaxial test data specified for a hyperelastic material:
odb_Hyperelastic hyperelastic = myMaterial.hyperelastic();
if (hyperelastic.hasValue()) {
bool testData = hyperelastic.testData();
odb_BiaxialTestData biaxialTestData =
10–25
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
hyperelastic.biaxialTestData();
odb_String smoothingType("smoothing type: ");
if (biaxialTestData.hasValue()) {
odb_Union smoothing = biaxialTestData.smoothing();
switch(smoothing.type()) {
case (odb_UNION_STRING):
smoothingType.append(smoothing.getString());
break;
case (odb_UNION_INT):
smoothingType.append(smoothing.getInt());
break;
case (odb_UNION_FLOAT):
smoothingType.append(smoothing.getFloat());
break;
case (odb_UNION_DOUBLE):
smoothingType.append(smoothing.getDouble());
break;
case (odb_UNION_BOOL):
smoothingType.append(smoothing.getBool());
break;
}
cout << smoothingType.CStr() << endl;
}
}
Chapter 59, “Material commands,” of the Abaqus Scripting Reference Manual, describes
the Material object commands in more detail; the odb_Union object is defined in “Union object,”
Section 63.6 of the Abaqus Scripting Reference Manual.
Sections
Sections are stored in the sections repository under the Odb object.
Extend the Section commands available to the Odb object using the following statement:
odb_SectionApi sectionApi;
odb.extendApi(odb_Enum::odb_SECTION,sectionApi);
The following statements display the repository keys of the sections in an output database:
odb_SectionContainer& sectionContainer =
sectionApi.sections();
odb_SectionContainerIT scIT(sectionContainer);
for (scIT.first(); !scIT.isDone(); scIT.next()) {
cout << "Section Name : " << scIT.currentKey().CStr() << endl;
}
10–26
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
The Section object can be one of the various section types. The odb_isA method can be used
to determine the section type. For example, to determine whether a section is of type “homogeneous
solid section” and to print it’s thickness and associated material name:
odb_ProfileContainer profileContainer =
sectionApi.profiles();
int numProfiles = sectionApi.numProfiles();
cout << "Total Number of profiles in the ODB: "
<< numProfiles << endl;
The Profile object can be one of the various profile types. The odb_isA method can be used
to determine the profile type. For example, to output the radius of all circular profiles in the odb:
odb_ProfileContainerIT pcIT(profileContainer);
for (pcIT.first(); !pcIT.isDone(); pcIT.next()) {
const odb_Profile& myProfile = pcIT.currentValue();
if (odb_isA(odb_CircularProfile,myProfile)) {
odb_CircularProfile circularProfile =
odb_dynamicCast( odb_CircularProfile, myProfile );
cout << "profile name = " << myProfile.name().CStr()
<< " radius = " << circularProfile.r();
}
}
Section assignments
Section assignments are stored in the sectionAssignments repository under the OdbAssembly
object.
10–27
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
All elements in an Abaqus analysis need to be associated with section and material properties.
Section assignments provide the relationship between elements in a part instance and their
section properties. The section properties include the associated material name. To access the
sectionAssignments repository from the PartInstance object:
odb_InstanceRepository& instanceRepository =
odb.rootAssembly().instances();
odb_InstanceRepositoryIT instIT(instanceRepository);
for (instIT.first(); !instIT.isDone(); instIT.next()) {
const odb_Instance& instance = instIT.currentValue();
odb_SequenceSectionAssignment sectionAssignmentSeq =
instance.sectionAssignments();
int sects = sectionAssignmentSeq.size();
cout << "Instance : " << instance.name().CStr() << endl;
for (int s = 0; s < sects; ++s) {
odb_SectionAssignment sa = sectionAssignmentSeq[s];
odb_String sectionName = sa.sectionName();
cout << " Section : " << sectionName.CStr() << endl;
odb_Set set = sa.region();
const odb_SequenceElement& elements = set.elements();
int size = elements.size();
cout << " Elements associated with this section : "
<< endl;
for (int e = 0; e< size; ++e)
cout << elements[e].label() << endl;
}
}
10–28
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
Step-1
Step-2
Step-3
Frames
Each step contains a sequence of frames, where each increment of the analysis (or each mode in an
eigenvalue analysis) that resulted in output to the output database is called a frame. The following
statement assigns a variable to the last frame in the first step:
odb_FieldOutputRepository& fieldOutputRep =
lastFrame.fieldOutputs();
odb_FieldOutputRepositoryIT fieldIter( fieldOutputRep );
for (fieldIter.first(); !fieldIter.isDone(); fieldIter.next())
cout << fieldIter.currentKey().CStr() << endl;
S
U
LE
CSHEAR1 ASURF/BSURF
CSLIP1 ASURF/BSURF
CPRESS ASURF/BSURF
COPEN ASURF/BSURF
UR3
10–29
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
Different variables can be written to the output database at different frequencies. As a result, not all
frames will contain all the field output variables.
You can use the following to view all the available field data in a frame:
The resulting print output lists all the field output variables in a particular frame, along with their type
and position.
S : Stress components
Type: 7
Number of fieldValues : 135
Number of locations : 1
U : Spatial displacement
Type: 3
Number of fieldValues : 161
Number of locations : 1
In turn, a FieldOutput object has a method values that returns a reference to a sequence of
FieldValue objects that contain data. Each FieldValue object in the sequence corresponds to a particular
location in the model. You can obtain the data corresponding to each FieldValue object using the data
method, which returns a pointer to an array that contains the results at the current location. For example,
10–30
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
int numComp = 0;
for (int i=0; i<numValues; i++) {
const odb_FieldValue val = displacements[i];
cout << "Node = " << val.nodeLabel();
const float* const U = val.data(numComp);
cout << ", U = ";
for (int comp=0;comp<numComp;comp++)
cout << U[comp] << " ";
cout << endl;
}
10–31
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
The numberOfElements method returns the number of elements in a block. When you
are accessing the results for elements, the numberOfElements method is useful in determining
the number of output locations per element. For example, when you are accessing element data at
integration points, you may need to determine the number of integration points per element. You can
determine the number of integration points per element by dividing the length of the block, which is the
total number of output locations, by the number of elements in the block. For example,
10–32
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
For more information, see “FieldBulkData object,” Section 60.5 of the Abaqus Scripting Reference
Manual.
The bulkDataBlocks method is an alternative to the values method of a FieldOutput object,
described in “Reading field output data,” Section 10.10.6. The values method of a FieldOutput object
returns a reference to a sequence of FieldValue objects that contain data. Each FieldValue object in the
sequence provides data for a unique location in the model.
Performance can be increased with the bulk data interface because the field data stored in a bulk
data block are made available in a single array of floating point numbers. If you access the same data
in nonbulk form, you must loop over a sequence of FieldValue objects and then access the data for each
location separately. Traversing an array can prove to be significantly faster than traversing a sequence of
objects and extracting data stored within the objects. As a result, accessing the data in an output database
using the bulk data interface can be significantly faster than the nonbulk form.
If you do not need to access large amounts of data, you may get better performance with the nonbulk
access method. This is especially true if the number of output locations you are accessing is smaller than
the number of elements in a class. Similarly, the nonbulk access method may be faster if the number of
nodes you are accessing is smaller than the number of nodes in an instance. The nonbulk access method
10–33
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
is also better suited for random access to an output database, where successive output locations to be
accessed may lie in completely different blocks.
10–34
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
odb_FieldOutput fieldIP =
fieldCAX4.getSubset(odb_Enum::INTEGRATION_POINT);
odb_FieldOutput fieldTopCenter = fieldIP.getSubset(topCenter);
const odb_SequenceFieldValue& vals = fieldTopCenter.values();
int valSize = vals.size();
int dSize = 0;
for (int l=0; l<valSize; l++) {
const odb_FieldValue val = vals[l];
cout << "Element label = " << val.elementLabel();
cout << " Integration Point = " << val.integrationPoint();
cout << endl;
const float* const data = val.data(dSize);
cout << " S : ";
for (int k=0; k < dSize; k++) {
cout << data[k] << " ";
}
cout << endl;
}
10–35
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
If the requested field values are not found in the output database at the specified
odb_Enum::ELEMENT_NODAL or odb_Enum::CENTROID positions, they are extrapolated from
the field data at the odb_Enum::INTEGRATION_POINT position.
10–36
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
odb
rootAssembly
parts
= Container sections
= Singular object
materials
Model Data
steps
frames
fieldOutputs
historyRegions
point
historyOutputs
data
Results Data
• S22
The history output data can be retrieved from the HistoryRegion objects in the output database. The
tutorial output database contains HistoryRegion objects that relate to the rigid body reference point and
the integration points of the corner element as follows:
• ’Node PART-1-1.1000’
• ’Element PART-1-1.1 Int Point 1’
• ’Element PART-1-1.1 Int Point 2’
• ’Element PART-1-1.1 Int Point 3’
• ’Element PART-1-1.1 Int Point 4’
The following statements read the tutorial output database and write the U2 history data from the second
step to an ASCII file that can be plotted by Abaqus/CAE:
10–37
Abaqus ID:
Printed on:
READING FROM AN OUTPUT DATABASE
step.getHistoryRegion(hPoint);
odb_HistoryOutputRepository& hoCon =
histRegion.historyOutputs();
odb_HistoryOutput& histOutU2 = hoCon["U2"];
odb_SequenceSequenceFloat data = histOutU2.data();
int numHDat = data.size();
for (int i=0; i<numHDat; i++) {
odb_SequenceFloat pair = data[i];
cout << pair.constGet(0) << " "
<< pair.constGet(1) << endl;
}
The output in this example is a sequence of tuples containing the frame time and the displacement value.
The example uses nodal history data output. If the analysis requested history output from an element, the
output database would contain one HistoryRegion object and one HistoryPoint object for each integration
point.
Create a variable that refers to the last frame of the first step.
Create a variable that refers to the displacement ’U’ in the last frame of the first step.
odb_FieldOutput& displacements =
lastFrame.fieldOutputs().get("U");
10–38
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
Create a variable that refers to the node set ’PUNCH’ in the part instance’PART-1–1’ :
odb_Instance& instance =
odb.rootAssembly().instances()["PART-1-1"];
odb_Set& nodeSet =
instance.nodeSets()["PUNCH"];
Create a variable that refers to the displacement of the node set in the last frame of the first step:
odb_FieldOutput myDisplacements =
displacements.getSubset(nodeSet);
Finally, print some field output data from each node in the node set (a single node in this example).
const odb_FieldValue val = myDisplacements.values()[0];
const float* const data = val.data(numComp);
cout << " Node: " << val.nodeLabel() << endl;
cout << " U = ";
for (int comp=0;comp<numComp;comp++)
cout << data[comp] << " ";
cout << endl;
cout << " Magnitude = " << val.magnitude();
The resulting output is
Node : 1000
U = 0.0000 , -76.4554
Magnitude = 76.4554
You can write your own data to an output database, and you can use Abaqus/CAE to view the data.
Writing to an output database is very similar to reading from an output database. When you open an
existing database, the Odb object contains all the objects found in the output database, such as instances,
steps, and field output data. In contrast, when you are writing to a new output database, these objects do
not exist. As a result you must use a constructor to create the objects. For example, you use the Part
constructor to create a Part object, the Instance constructor to create an OdbInstance object, and the
Step constructor to create an OdbStep object.
After you create an object, you use methods of the objects to enter or modify the data associated
with the object. For example, if you are creating an output database, you first create an Odb object.
You then use the Part constructor to create a part. After creating the part, you use the addNodes and
addElements methods of the Part object to add nodes and elements, respectively. Similarly, you use
10–39
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
the addData method of the FieldOutput object to add field output data to the output database. After
creating an output database, you should use the save method on the Odb object to save the output
database.
The example program in “Creating an output database,” Section 10.15.2, also illustrates how you
can write to an output database.
The following topics are covered:
• “Creating a new output database,” Section 10.11.1
• “Writing model data,” Section 10.11.2
• “Writing results data,” Section 10.11.3
• “Writing field output data,” Section 10.11.4
• “Default display properties,” Section 10.11.5
• “Writing history output data,” Section 10.11.6
For a full description of the Odb command, see “Odb object,” Section 60.1 of the Abaqus Scripting
Reference Manual. Abaqus creates the RootAssembly object when you create or open an output database.
You use the save method to save the output database.
odb.save();
For a full description of the save command, see “save,” Section 60.1.4 of the Abaqus Scripting
Reference Manual.
10–40
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
Part
If the part was created by Abaqus/CAE, the description of the native Abaqus/CAE geometry
is stored in the model database, but it is not stored in the output database. A part is stored in
an output database as a collection of nodes, elements, surfaces, and sets. You use the Part
constructor to add a part to the Odb object. You can specify the type of the part; however, only
DEFORMABLE_BODY is currently supported. For example,
For a full description of the Part constructor, see “OdbPart object,” Section 60.21 of the Abaqus
Scripting Reference Manual. The new Part object is empty and does not contain geometry. After
you create the Part object, you then add nodes and elements.
You use the addNodes method to add nodes by defining node labels and coordinates. You
can also define an optional node set. For example,
odb_SequenceInt nodeLabels;
nodeLabels.append(1);
nodeLabels.append(2);
nodeLabels.append(3);
nodeLabels.append(5);
nodeLabels.append(7);
nodeLabels.append(11);
double c[6][3] = { {2.0, 1.0, 0.0},
{1.0, 1.0, 0.0},
{1.0, 0.0, 0.0},
{2.0, 0.0, 0.0},
{1.0, 0.0, 1.0},
{2.0, 0.0, 1.0} };
odb_SequenceSequenceFloat nodeCoor;
for (int n=0; n<nodeLabels.size(); n++) {
odb_SequenceFloat loc;
for (int i=0; i<3; i++)
loc.append(c[n][i]);
nodeCoor.append(loc);
}
part1.addNodes(nodeLabels, nodeCoor, "nodes_1");
10–41
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
For a full description of the addNodes command, see “addNodes,” Section 60.21.3 of the Abaqus
Scripting Reference Manual.
Similarly, you use the addElements method to add elements to the part using a sequence
of element labels, element connectivity, and element type. You can also define an optional element
set and an optional section category. For example,
odb_SequenceInt elLabels;
elLabels.append(9);
elLabels.append(99);
odb_SequenceSequenceInt connect;
const int numNodePerEl = 4;
int conn[2][numNodePerEl] = {{1, 2, 3, 5},
{5, 3, 7, 11}};
for (int e=0; e<elLabels.size(); e++) {
odb_SequenceInt l;
for (int i=0; i<numNodePerEl; i++)
l.append(conn[e][i]);
connect.append(l);
}
part1.addElements(elLabels, connect, "S4R",
"s4_els", shellCat);
For a full description of the addElements command, see “addElements,” Section 60.21.2 of the
Abaqus Scripting Reference Manual.
You can create both instances and regions on the RootAssembly object.
Part instances
You use the Instance constructor to create part instances of the parts you have already defined
using the Part constructor. For example,
odb_Instance& instanceA =
odb.rootAssembly().Instance("part-1-1", part1);
10–42
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
You can also supply an optional local coordinate system that specifies the rotation and translation
of the part instance. You can add nodes and elements only to a part; you cannot add elements
and nodes to a part instance. As a result, you should create the nodes and elements that define the
geometry of a part before you instance the part. For a full description of the Instance command,
see “OdbInstance object,” Section 60.17 of the Abaqus Scripting Reference Manual.
Regions
Region commands are used to create sets from element labels, node labels, and element faces. You
can create a set on a part, part instance, or the root assembly. Node and element labels are unique
within an instance but not within the assembly. As a result, a set on the root assembly requires
the names of the part instances associated with the nodes and elements. You can also use region
commands to create surfaces.
For example,
// An ElementSet on an instance
odb_SequenceInt eLabelsA(2);
eLabelsA.append(9);
eLabelsA.append(99);
instanceA.ElementSet("elSetA", eLabelsA);
odb_SequenceSequenceInt nodeLabelsRA;
odb_SequenceString namesRA;
namesRA.append("part-1-1");
odb_SequenceInt nodeLabelsRA_A;
nodeLabelsRA_A.append(5);
nodeLabelsRA_A.append(11);
nodeLabelsRA.append(nodeLabelsRA_A);
const odb_Set& nSetRA = rootAssy.NodeSet("nodeSetRA",
namesRA, nodeLabelsRA);
Materials
You use the Material object to list material properties. Material objects are members of the Odb
object.
Materials are stored in the materials repository under the Odb object.
Extend the Material commands available to the Odb object using the following statement:
odb_MaterialApi materialApi;
odb.extendApi(odb_Enum::odb_MATERIAL,materialApi);
10–43
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
To create an isotropic elastic material, with a Young’s modulus of 12000.0 and an effective
Poisson’s ratio of 0.3 in the output database:
odb_String materialName("Elastic Material");
odb_Material& material = materialApi.Material(materialName);
odb_SequenceSequenceFloat myTable;
odb_SequenceFloat myData;
myData.append(12000.0); myData.append(0.3);
myTable.append(myData);
odb_String type("ISOTROPIC");
material.Elastic(myTable,type);
For more information, see Chapter 59, “Material commands,” of the Abaqus Scripting
Reference Manual.
Sections
You use the Section object to create sections and profiles. Section objects are members of the Odb
object.
Sections are stored in the sections repository under the Odb object.
Extend the API commands available to the Odb object using the following statement:
odb_SectionApi sectionApi;
odb.extendApi(odb_Enum::odb_SECTION,
sectionApi);
The following code creates a homogeneous solid section object. A Material object must be
present before creating a Section object. An exception is thrown if the material does not exist.
Section assignments
You use the SectionAssignment object to assign sections and their associated material properties
to regions of the model. SectionAssignment objects are members of the Odb object. For a full
10–44
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
description of the assignSection method, see “assignSection,” Section 60.17.7 of the Abaqus
Scripting Reference Manual.
All Elements in an Abaqus analysis need to be associated with section and material properties.
Section assignments provide the relationship between elements in an Instance object and their
section properties. The section properties include the associated material name. To create an
element set and assign a section:
odb_SequenceInt setLabels;
setLabels.append(1);
setLabels.append(2);
elsetName = "Material 1";
odb_Set& elset = instance.ElementSet(elsetName,setLabels);
// section assignment on instance
instance.assignSection(elset,section);
The Step constructor has an optional previousStepName argument that specifies the step after
which this step must be inserted in the steps repository. For a full description of the Step
command, see “Step,” Section 60.26.1 of the Abaqus Scripting Reference Manual.
Frames
You use the Frame constructor to create a frame for field output. For example,
10–45
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
For a full description of the Frame command, see “Frame,” Section 60.16.3 of the Abaqus Scripting
Reference Manual.
// vector
odb_SequenceString vectorCompLabels;
vectorCompLabels.append("1");
vectorCompLabels.append("2");
vectorCompLabels.append("3");
odb_SequenceInvariant vectorInvar;
vectorInvar.append(odb_Enum::MAGNITUDE);
odb_FieldOutput& vectorField = frameOne.FieldOutput("U",
"displacement vector",
odb_Enum::VECTOR,
vectorCompLabels, vectorInvar);
odb_SequenceInt labels2;
labels2.append(3);
labels2.append(5);
odb_SequenceSequenceFloat vecDat;
odb_SequenceFloat v1;
v1.append(1.1); v1.append(1.2); v1.append(1.3);
vecDat.append(v1);
odb_SequenceFloat v2;
v2.append(2.1); v2.append(2.2); v2.append(2.3);
vecDat.append(v2);
vectorField.addData(odb_Enum::NODAL, instanceA,
labels2, vecDat);
10–46
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
For a full description of the FieldOutput constructor, see “FieldOutput,” Section 60.7.1 of the Abaqus
Scripting Reference Manual.
The type argument to the FieldOutput constructor describes the type of the data—tensor, vector,
or scalar. The properties of the different tensor types are:
Full tensor
A tensor that has six components and three principal values. Full three-dimensional rotation of the
tensor is possible.
The valid components and invariants for the different data types are given in Table 10–1.
Table 10–1 Valid components and invariants for Abaqus data types.
10–47
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
For example, the following statements add element data to the FieldOutput object:
odb_SequenceString tensorCompLabels;
tensorCompLabels.append("s11");
tensorCompLabels.append("s22");
tensorCompLabels.append("s33");
tensorCompLabels.append("s12");
tensorCompLabels.append("s13");
tensorCompLabels.append("s23");
odb_SequenceInvariant tensorInvar;
tensorInvar.append(odb_Enum::MISES);
tensorInvar.append(odb_Enum::TRESCA);
tensorInvar.append(odb_Enum::MAX_PRINCIPAL);
tensorInvar.append(odb_Enum::MID_PRINCIPAL);
tensorInvar.append(odb_Enum::MIN_PRINCIPAL);
odb_SequenceInt tensorLabels;
10–48
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
tensorLabels.append(9);
tensorLabels.append(99);
odb_SequenceSequenceFloat tensorDat;
odb_SequenceFloat t1;
t1.append(1.0); t1.append(2.0); t1.append(3.0);
t1.append(0.0); t1.append(0.0); t1.append(0.0);
odb_SequenceFloat t2;
t2.append(120.0); t2.append(-55.0); t2.append(-85.0);
t2.append(-55.0); t2.append(-75.0); t2.append(33.0);
tensorDat.append(t1);
tensorDat.append(t2);
For a full description of the addData command, see “addData,” Section 60.7.6 of the Abaqus Scripting
Reference Manual.
As a convenience, localCoordSystem can be a single transform or a list of transforms. If
localCoordSystem is a single transform, it applies to all values. If localCoordSystem is a list of
transforms, the number of items in the list must match the number of data values.
step1.setDefaultField(tensorField);
step1.setDefaultDeformedField(vectorField);
You can set a different default field variable and deformed field variable for different steps.
10–49
Abaqus ID:
Printed on:
WRITING TO AN OUTPUT DATABASE
odb_HistoryPoint hPoint1(instanceA.elements(0));
For a full description of the HistoryPoint command, see “HistoryPoint,” Section 60.10.1 of the
Abaqus Scripting Reference Manual.
You then use the HistoryRegion constructor to create a HistoryRegion object:
For a full description of the HistoryRegion command, see “HistoryRegion,” Section 60.11.1 of the
Abaqus Scripting Reference Manual.
You use the HistoryOutput constructor to add variables to the HistoryRegion object.
10–50
Abaqus ID:
Printed on:
EXCEPTION HANDLING IN AN OUTPUT DATABASE
ho1.addData(0.001, 0.1);
odb_SequenceFloat timeData;
odb_SequenceFloat values;
timeData.append(0.001);
values.append(0.1);
ho1.addData(timeData, values);
10–51
Abaqus ID:
Printed on:
COMPUTATIONS WITH Abaqus RESULTS
For more information, see “BaseException object,” Section 63.1 of the Abaqus Scripting Reference
Manual.
10–52
Abaqus ID:
Printed on:
COMPUTATIONS WITH Abaqus RESULTS
not the same, the local coordinate systems of both fieldValues will be disregarded and the fieldValue
computed will have no local coordinate system.
• The operations will not be performed on the conjugate data (the imaginary portion of a complex
result).
The FieldOutput operations are significantly more efficient than the FieldValue operators. You can
save the computed FieldOutput objects with the following procedure:
• Create a new FieldOutput object in the output database.
• Use the addData method to add the new computed field objects to the new FieldOutput object.
For example,
10–53
Abaqus ID:
Printed on:
COMPUTATIONS WITH Abaqus RESULTS
10–54
Abaqus ID:
Printed on:
COMPUTATIONS WITH Abaqus RESULTS
odb_SequenceFieldOutput flds =
10–55
Abaqus ID:
Printed on:
COMPUTATIONS WITH Abaqus RESULTS
maxEnvelope(odb_SequenceFieldOutput& fields,
odb_Enum::odb_InvariantEnum invariant);
odb_SequenceFieldOutput flds =
minEnvelope(odb_SequenceFieldOutput& fields,
odb_Enum::odb_InvariantEnum invariant);
odb_SequenceFieldOutput flds =
maxEnvelope(odb_SequenceFieldOutput& fields,
const odb_String& componentLabel);
odb_SequenceFieldOutput flds =
minEnvelope(odb_SequenceFieldOutput& fields,
const odb_String& componentLabel);
The envelope commands return two FieldOutput objects.
• The first object contains the requested extreme values.
• The second object contains the indices of the fields for which the extreme values were found. The
indices derive from the order in which you supplied the fields to the command.
The optional invariant argument is a odb_Enum::odb_DataTypeEnum specifying the invariant to be used
when comparing vectors or tensors. The optional componentLabel argument is an odb_String specifying
the component of the vector or tensor to be used for selecting the extreme value.
The following rules apply to envelope calculations:
• Abaqus compares the values using scalar data. If you are looking for the extreme value of a vector
or a tensor, you must supply an invariant or a component label for the selection of the extreme value.
For example, for vectors you can supply the MAGNITUDE invariant and for tensors you can supply
the MISES invariant.
• The fields being compared must be similar. For example,
– VECTOR and TENSOR_3D_FULL fields cannot appear in the same list.
– The output region of all the fields must be the same. All the fields must apply to the whole
model, or all the fields must apply to the same set.
10–56
Abaqus ID:
Printed on:
COMPUTATIONS WITH Abaqus RESULTS
• Beams, truss, and axisymmetric shell element results will not be transformed.
• The component directions 1, 2, and 3 of the transformed results will correspond to the system
directions X, Y, and Z for rectangular coordinate systems; R, , and Z for cylindrical coordinate
systems; and R, , and for spherical coordinate systems.
Note: Stress results for three-dimensional continuum elements transformed into a cylindrical
system would have the hoop stress in S22, which is consistent with the coordinate system axis
but inconsistent with the stress state for a three-dimensional axisymmetric elements having hoop
stress in S33.
• When you are transforming a tensor, the location or integration point always takes into account the
deformation. The location of the coordinate system depends on the model, as follows:
– If the system is fixed, the coordinate system is fixed.
– If the system is model based, you must supply a displacement field that determines the
instantaneous location and orientation of the coordinate system.
• Abaqus will perform transformations of tensor results for shells, membranes, and planar elements
as rotations of results about the element normal at the element result location. The element normal
is the normal computed for the frame associated with the field by Abaqus, and you cannot redefine
the normal. Abaqus defines the location of the results location from the nodal locations. You
specify optional arguments if you want to use the deformed nodal locations to transform results.
For rectangular, cylindrical, and spherical coordinate systems the second component direction for
the transformed results will be determined by one of the following:
– The Y-axis in a rectangular coordinate system.
– The -axis in a cylindrical coordinate system.
– The -axis in a spherical coordinate system.
– A user-specified datum axis projected onto the element plane.
If the coordinate system used for projection and the element normal have an angle less than the
specified tolerance (the default is 30°), Abaqus will use the next axis and generate a warning.
10–57
Abaqus ID:
Printed on:
IMPROVING THE EFFICIENCY OF YOUR SCRIPTS
If you are accessing large amounts of data from an output database, you should be aware of potential
inefficiencies in your program and techniques that will help to speed up your scripts.
• “Creating objects to hold loop counters,” Section 10.14.1
• “Creating objects to hold temporary variables,” Section 10.14.2
• “Using references to objects,” Section 10.14.3
You can make the program more efficient if you create an object to hold the value of the number of nodes.
You can use this technique only if the maximum value of the loop counter remains fixed for the duration
of the loop.
10–58
Abaqus ID:
Printed on:
IMPROVING THE EFFICIENCY OF YOUR SCRIPTS
In this example every time the script calls the strainField.values method, Abaqus must
reconstruct the sequence of FieldValue objects. This reconstruction could result in a significant
performance degradation, particularly for a large model.
10–59
Abaqus ID:
Printed on:
IMPROVING THE EFFICIENCY OF YOUR SCRIPTS
A slight change in the program greatly improves its performance, as shown in the following
example:
odb_FieldOutputRepository& fieldRep =
frame1.fieldOutputs();
odb_FieldOutput& stressField = fieldRep.get("S");
odb_FieldOutput& strainField = fieldRep.get("LE");
const odb_SequenceFieldValue& seqStressVal =
stressField.values();
const odb_SequenceFieldValue& seqStrainVal =
strainField.values();
int numFV = seqStressVal.size();
int strainComp = 0;
for (int loc=0; loc < numFV; loc++) {
const odb_FieldValue stressVal = seqStressVal[loc];
if (stressVal.mises() > stressCap) {
cout << "Element label = " << stressVal.elementLabel()
<< endl;
cout << "Integration Point = "
<< stressVal.integrationPoint() << endl;
const odb_FieldValue strainVal = seqStrainVal[loc];
const float* data = strainVal.data(strainComp);
cout << " LE : ";
for (int comp = 0; comp < strainComp; comp++)
cout << data[comp];
cout << endl;
}
}
}
Similarly, if you expect to retrieve more than one frame from an output database, you should create
a temporary variable that holds the entire frame repository. You can then provide the logic to retrieve the
desired frames from the repository and avoid recreating the repository each time. For example, executing
the following statements could be very slow:
10–60
Abaqus ID:
Printed on:
IMPROVING THE EFFICIENCY OF YOUR SCRIPTS
Creating a temporary variable to hold the frame repository provides the same functionality and speeds
up the process:
Such a potential loss of performance will not be a problem when accessing a load case frame. Accessing
a load case frame does not result in the creation of a frame repository and, thus, does not suffer from a
corresponding loss of performance.
In the above case a copy of the nodeSequence object has to be created in memory.
Many of the methods in the Abaqus Scripting Interface that provide access to an output database
return a reference to an object rather than the object itself. It is much more efficient to modify the previous
example to specify the returned type to be a reference:
In this case no new object is created and no copy constructors are called.
10–61
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
The following examples illustrate how you use the output database commands to access data from an
output database:
• “Finding the maximum value of von Mises stress,” Section 10.15.1
• “Creating an output database,” Section 10.15.2
• “Reading data from an output database,” Section 10.15.3
• “Decreasing the amount of data in an output database by retaining data at specific frames,”
Section 10.15.4
• “Stress range for multiple load cases,” Section 10.15.5
• “A C++ version of FELBOW,” Section 10.15.6
10–62
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
/***************************************************************
odbMaxMises.C
Code to determine the location and value of the maximum
von-mises stress in an output database.
Usage: abaqus odbMaxMises -odb odbName -elset(optional)
elsetName
Requirements:
1. -odb : Name of the output database.
2. -elset : Name of the assembly level element set.
Search will be done only for element belonging
to this set. If this parameter is not provided,
search will be performed over the entire model.
3. -help : Print usage
****************************************************************/
#if (defined(HP) && (! defined(HKS_HPUXI)))
#include <iostream.h>
#else
#include <iostream>
using namespace std;
#endif
#include <odb_API.h>
#include <sys/stat.h>
/*
***************
utility functions
***************
*/
bool fileExists(const odb_String &string);
void rightTrim(odb_String &string,const char* char_set);
void printExecutionSummary();
/***************************************************************/
10–63
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
10–64
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
{
if (myAssembly.elementSets().isMember(elsetName))
{
myElset = myAssembly.elementSets()[elsetName];
region = " in the element set : " + elsetName;
}
else
{
cerr<<"An assembly level elset " << elsetName.CStr()
<< " does not exist in the output database :"
<< myOdb.name().CStr() << endl;
myOdb.close();
exit(0);
}
}
// Initialize maximum values.
float maxMises = -0.1;
int numFV = 0;
int maxElem = 0;
odb_String maxStep = "__None__";
int maxFrame = -1;
static const odb_String Stress = "S";
bool isStressPresent = false;
int numBD = 0,numElems = 0, numIP = 0, numComp = 0, position = 0;
// Iterate over all available steps
odb_StepRepository& sRep1 = myOdb.steps();
odb_StepRepositoryIT sIter1 (sRep1);
for (sIter1.first(); !sIter1.isDone(); sIter1.next())
{
odb_Step& step = sRep1[sIter1.currentKey()];
cout<<"Processing Step: "<<step.name().CStr()<<endl;
odb_SequenceFrame& frameSequence = step.frames();
int numFrames = frameSequence.size();
for (int f = 0; f<numFrames; f++)
{
odb_Frame& frame = frameSequence[f];
odb_FieldOutputRepository& fieldRep = frame.fieldOutputs();
if (fieldRep.isMember(Stress))
{
isStressPresent = true;
odb_FieldOutput field = fieldRep.get(Stress);
10–65
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
if (ifElset)
field = field.getSubset(myElset);
const odb_SequenceFieldBulkData& seqVal =
field.bulkDataBlocks();
int numBlocks = seqVal.size();
for ( int iblock=0; iblock<numBlocks; iblock++)
{
const odb_FieldBulkData& bulkData = seqVal[iblock];
numBD = bulkData.length();
numElems = bulkData.numberOfElements();
numIP = numBD/numElems;
numComp = bulkData.width();
float* mises = bulkData.mises();
int* elementLabels = bulkData.elementLabels();
int* integrationPoints = bulkData.integrationPoints();
for (int elem=0; elem<numElems; elem++)
{
for (int ip=0; ip<numIP; ip++)
{
position = elem*numIP+ip;
float misesData = mises[position];
if (misesData > maxMises)
{
maxMises = misesData;
maxElem = elementLabels[elem];
maxStep = step.name();
maxFrame = frame.incrementNumber();
}
}
}
}
}
}
}
if (isStressPresent)
{
cout << "Maximum von Mises stress " << region.CStr()
<< " is " << maxMises << " in element "
<< maxElem << endl;
cout << "Location: frame # " << maxFrame << " step: "
10–66
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
void printExecutionSummary()
{
cout << " Code to determine the location and value of the\n"
<< " maximum von-mises stress in an output database.\n"
<< " Usage: abaqus odbMaxMises -odb odbName \n"
<< " -elset(optional), -elsetName\n"
<< " Requirements:\n"
<< " 1. -odb : Name of the output database.\n"
<< " 2. -elset : Name of the assembly level element set.\n"
<< " Search will be done only for element \n"
<< " belonging to this set.\n"
<< " If this parameter is not provided, search \n"
<< " will be performed over the entire model.\n"
10–67
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
////////////////////////////////////////////////////
// Code to create an output database and add model,
// field, and history data. The code also reads
// history data, performs an operation on the data, and writes
// the result back to the output database.
//
// SECTION: System includes
//
#include <math.h>
//
// Begin local includes
//
#include <odb_API.h>
#include <odb_MaterialTypes.h>
#include <odb_SectionTypes.h>
//
// End local includes
//
10–68
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
odb_String path("odbWriteC.odb");
odb_Odb& odb = Odb(name,
analysisTitle,
description,
path);
// Model data:
// Set up the section categories.
odb_String sectionCategoryName("S5");
odb_String sectionCategoryDescription("Five-Layered Shell");
odb_SectionCategory& sCat =
odb.SectionCategory(sectionCategoryName,
sectionCategoryDescription);
int sectionPointNumber = 1;
odb_String sectionPointDescription("Bottom");
odb_SectionPoint spBot =
sCat.SectionPoint(sectionPointNumber,
sectionPointDescription);
sectionPointNumber = 3;
sectionPointDescription = "Middle";
odb_SectionPoint spMid =
sCat.SectionPoint(sectionPointNumber,
sectionPointDescription);
sectionPointNumber = 5;
sectionPointDescription = "Top";
odb_SectionPoint spTop =
sCat.SectionPoint(sectionPointNumber,
sectionPointDescription);
10–69
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
noTension,
temperatureDependency,
dependencies,
moduli);
odb_SequenceInt nodeLabels;
for(n=1; n<7; n++)
nodeLabels.append(n);
odb_SequenceInt elLabels;
elLabels.append(1);
elLabels.append(2);
odb_SequenceSequenceInt connect;
const int numNodePerEl = 4;
int conn[2][numNodePerEl] = { {1, 2, 3, 4},
10–70
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
{6, 5, 2, 1} };
for(int e=0; e<elLabels.size(); e++) {
odb_SequenceInt l;
for(int i=0; i<numNodePerEl; i++)
l.append(conn[e][i]);
connect.append(l);
}
odb_String elType("S4");
odb_String elsetName("eset-1");
part1.addElements(elLabels,
connect,
elType,
elsetName,
sCat);
10–71
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
odb_SequenceSequenceFloat dispData;
odb_SequenceFloat dispData1[6];
// create some displacement values
for(n=0; n<6; n++) {
for(int m=1; m<4; m++)
dispData1[n].append(n*3+m);
dispData.append(dispData1[n]);
}
uField.addData(odb_Enum::NODAL,
instance1,
nodeLabels,
dispData);
step1.setDefaultDeformedField(uField);
odb_SequenceSequenceFloat topData;
odb_SequenceFloat topData1;
for(n=1; n<5; n++)
topData1.append(n);
odb_SequenceSequenceFloat bottomData;
odb_SequenceFloat bottomData1;
for(n=1; n<5; n++)
bottomData1.append(n);
odb_SequenceSequenceFloat transform;
10–72
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
}
transform.append(transform1);
}
odb_SequenceString componentLabels;
componentLabels.append("S11");
componentLabels.append("S22");
componentLabels.append("S33");
componentLabels.append("S12");
odb_SequenceInvariant validInvariants;
validInvariants.append(odb_Enum::MISES);
fieldName = "S";
fieldDescription = "Stress";
odb_FieldOutput& sField =
frame1.FieldOutput(fieldName,
fieldDescription,
odb_Enum::TENSOR_3D_PLANAR,
componentLabels,
validInvariants);
sField.addData(odb_Enum::INTEGRATION_POINT,
instance1,
elLabels,
topData,
spTop,
transform);
sField.addData(odb_Enum::INTEGRATION_POINT,
instance1,
elLabels,
bottomData,
spBot,
transform);
step1.setDefaultField(sField);
// History data:
// Create a HistoryRegion for a specific point.
odb_HistoryPoint hPoint1(instance1.getNodeFromLabel(1));
odb_String historyRegionName("historyNode0");
odb_String historyRegionDescription(
"Displacement and reaction force");
odb_HistoryRegion& hRegionStep1 =
step1.HistoryRegion(historyRegionName,
historyRegionDescription,
10–73
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
hPoint1);
odb_String historyOutputName("U1");
odb_String historyOutputDescription("Displacement");
odb_HistoryOutput& hOutputStep1U1 =
hRegionStep1.HistoryOutput(historyOutputName,
historyOutputDescription,
odb_Enum::SCALAR);
historyOutputName = "RF1";
historyOutputDescription = "Reaction Force";
odb_HistoryOutput& hOutputStep1Rf1 =
hRegionStep1.HistoryOutput(historyOutputName,
historyOutputDescription,
odb_Enum::SCALAR);
// Add history data for step1.
hOutputStep1U1.addData(0.0, 0.0);
hOutputStep1Rf1.addData(0.0,0.0);
hOutputStep1U1.addData(0.1, 0.1);
hOutputStep1Rf1.addData(0.1,0.1);
hOutputStep1U1.addData(0.3, 0.3);
hOutputStep1Rf1.addData(0.3,0.3);
hOutputStep1U1.addData(1.0, 0.5);
hOutputStep1Rf1.addData(1.0,0.5);
odb_HistoryPoint hPoint2(instance1.getNodeFromLabel(1));
odb_HistoryRegion& hRegionStep2 =
step2.HistoryRegion(historyRegionName,
historyRegionDescription,
hPoint2);
10–74
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
hRegionStep2.HistoryOutput(historyOutputName,
historyOutputDescription,
odb_Enum::SCALAR);
historyOutputName = "RF1";
historyOutputDescription = "Reaction Force";
odb_HistoryOutput& hOutputStep2Rf1 =
hRegionStep2.HistoryOutput(historyOutputName,
historyOutputDescription,
odb_Enum::SCALAR);
hOutputStep2U1.addData(1.2, 0.8);
hOutputStep2Rf1.addData(1.2,0.9);
hOutputStep2U1.addData(1.9, 0.9);
hOutputStep2Rf1.addData(1.9,1.1);
hOutputStep2U1.addData(3.0, 1.3);
hOutputStep2Rf1.addData(3.0,1.3);
hOutputStep2U1.addData(4.0, 1.5);
hOutputStep2Rf1.addData(4.0,1.5);
odb_HistoryOutputRepository& historyOutputs1 =
hRegionStep1.historyOutputs();
historyOutputName = "U1";
odb_HistoryOutput& u1FromStep1 =
historyOutputs1[historyOutputName];
odb_HistoryOutputRepository& historyOutputs2 =
hRegionStep2.historyOutputs();
odb_HistoryOutput& u1FromStep2 =
historyOutputs2[historyOutputName];
10–75
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
odb.save();
odb.close();
return 0;
}
10–76
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
• Sections
• Steps
• Frames
• Fields
• Field values
• Field bulk data
• Field locations
• History regions
• History output
• History points
Use the following command to retrieve the example program:
abaqus fetch job=odbDump
-largeOdb odbName
The name of the large output database generated by the original problem. The program copies
selected frames from this output database.
10–77
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
-history
Copy all history output from all available steps in the large output database. By default, history
output is not copied.
WARNING: Copying large amounts of history data can result in the program creating a very large
output database.
-debug
Print a detailed report of all the operations performed during the running of the program. By default,
no debug information is generated.
WARNING: If you are extracting data from a large output database, the debug option can generate
large amounts of information.
You can also run the example with only the -help parameter for a summary of the usage.
The following is an example of how you can use this program in conjunction with the output
database generated by the problem described in “Free ring under initial velocity: comparison of rate-
independent and rate-dependent plasticity,” Section 1.3.4 of the Abaqus Benchmarks Manual. Use the
following commands to retrieve the example program and the benchmark input file:
abaqus fetch job=odbFilter.C
abaqus fetch job=ringshell.inp
1. Run an analysis using the benchmark input file:
abaqus job=ringshell
This creates an output database called ringshell.odb that contains 100 frames of data.
2. Run a datacheck analysis to obtain a new output database called ringshell_datacheck.odb
that contains the same model data as ringshell.odb:
abaqus job=ringshell_datacheck -input ringshell datacheck
3. Create the executable program:
abaqus make job=odbFilter.C
The program displays the number of frames available in each step. For each step you must specify
the number of increments between frames, which is the frequency at which the data will be copied to the
new output database. Data for the first and last increment in each step are always copied. For example,
if a step has 100 frames, and you enter a frame interval of 37, the program will copy data for frames 0,
37, 74, and 100.
The following statement will run the executable program and read data from the small output
database containing only model data and the large output database created by the benchmark example:
abaqus odbFilter -smallOdb ringshell_datacheck -largeOdb ringshell
10–78
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
////////////////////////////////////////////////////
// Code to compute a stress range from
// all the load cases in a step.
//
10–79
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
// System includes
#if (defined(HP) && (! defined(HKS_HPUXI)))
#include <iostream.h>
#else
#include <iostream>
using namespace std;
#endif
odb_String odbName(argv[1]);
odb_String stepName(argv[2]);
cout << "Computing for odb \"" << odbName.CStr() << "\"";
cout << " and step \"" << stepName.CStr() << "\"." << endl;
10–80
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
return 0;
}
odb_FieldOutput
computeStressRange(odb_Step& step)
{
// collect stress fields for all load cases
odb_SequenceFieldOutput sFields;
odb_LoadCaseRepositoryIT iter(step.loadCases());
for( iter.first(); !iter.isDone(); iter.next() ) {
odb_Frame frame = step.getFrame( iter.currentValue() );
odb_FieldOutput& stressField = frame.fieldOutputs()["S"];
sFields.append(stressField.getScalarField("S11"));
};
10–81
Abaqus ID:
Printed on:
EXAMPLE PROGRAMS THAT ACCESS DATA FROM AN OUTPUT DATABASE
To use option 2, you must ensure that the integration point coordinates (COORD) are written to
the output database. For option 1 the X-data are data for the distance around the circumference of the
elbow element, measured along the middle surface, and the Y-data are data for the output variable. For
option 2 the X–Y data are the current coordinates of the middle-surface integration points around the
circumference of the elbow element, projected to a local coordinate system in the plane of the deformed
cross-section. The origin of the local system coincides with the center of the cross-section; the plane of
the deformed cross-section is defined as the plane that contains the center of the cross-section.
You should specify the name of the output database during program execution. The program
prompts for additional information, depending on the option that was chosen; this information includes
the following:
• Your choice for storing results (ASCII file or a new output database)
• File name based on the above choice
• The postprocessing option (1 or 2)
• The part name
• The step name
• The frame number
• The element output variable (option 1 only)
• The component of the variable (option 1 only)
• The section point number (option 1 only)
• The element number or element set name
Before program execution, compile and link the C++ program using the abaqus make utility:
The program prompts for other information, such as the desired postprocessing option, part name, etc.
The program processes the data and produces a text file or a new output database file that contains the
information required to visualize the elbow element results.
“Elastic-plastic collapse of a thin-walled elbow under in-plane bending and internal pressure,”
Section 1.1.2 of the Abaqus Example Problems Manual, contains several figures that can be created
with the aid of this program.
10–82
Abaqus ID:
Printed on:
Abaqus, the 3DS logo, SIMULIA, CATIA, SolidWorks, DELMIA, ENOVIA, 3DVIA, and Unified FEA are trademarks or registered trademarks of Dassault Systèmes or its
subsidiaries in the US and/or other countries. Other company, product, and service names may be trademarks or service marks of their respective owners. © Dassault Systèmes, 2011
About SIMULIA
SIMULIA is the Dassault Systèmes brand that delivers a scalable portfolio of
Realistic Simulation solutions including the Abaqus product suite for Unified Finite
Element Analysis; multiphysics solutions for insight into challenging engineering
problems; and lifecycle management solutions for managing simulation data,
processes, and intellectual property. By building on established technology,
respected quality, and superior customer service, SIMULIA makes realistic
simulation an integral business practice that improves product performance,
reduces physical prototypes, and drives innovation. Headquartered in Providence,
RI, USA, with R&D centers in Providence and in Vélizy, France, SIMULIA provides
sales, services, and support through a global network of regional offices and
distributors. For more information, visit www.simulia.com.
www.3ds.com