Realistic Finite Element Problems Might Consist of Up To Hundreds of
Realistic Finite Element Problems Might Consist of Up To Hundreds of
The data lines basically provide data, if required, that is associated with the
option block used. In the example in Case 1, element identification (i.d.) and
the nodes that make up the element are the necessary data required. Note
that the information provided in the data lines would vary with some of the
parameters defined in the keyword line. For example, if the element type
being used is a 2D plane stress element, then the data lines would require
different data, as shown in Case 2.
Note that, in this case, the element type CPS4 which represents four-nodal,
2D solid elements is being used, and correspondingly, the data lines must
include the element i.d. (as before), and four nodes that make up each
element instead of two for the case of the beam element previously.
Using Sets
In the last section, it was seen in Case 1 that elements can be grouped into
a set for future reference by other option blocks. A set can be a grouping of
nodes or a grouping of elements. The analyst will usually provide a name for
the set that contains between 1-80 characters. For example, in Case 1,
‘BEAM’ is the name of the set containing elements 1 and 2; and in Case 2,
‘PLSTRESS’ is the name of the set containing elements 1 and 2. In both
examples, the sets are defined together with the definition of the elements
themselves in the element block. However, sets can also be defined as a
separate block on their own. In Case 3, the pinned support of the 1D beam is
to be defined. Nodes 1 and 11 (provided in the data line for the ‘NSET’ block)
are first grouped in the node set called ‘SUPPORT’. Then using the
‘BOUNDARY’ option block, the node set, ‘SUPPORT’ is referenced to
constrain the DOFs, 1 and 2 (x translation and y translation) to zero. In other
words, rather than having four data lines for the two nodes 1 and 11 (each
node having to constrain two DOFs), we now have only two lines with the
reference to the node set ‘SUPPORT’. Another thing to note in this example is
the use of comment lines. Comment lines exist in ABAQUS too, just as in
most programming languages. Comment lines begin with two asterisks, **,
and whatever follows in that line after that will not be read by ABAQUS as
input information defining the model.
This is, of course, a very simple case, and the reduction in the number of
data lines from four to two may not seem very significant. However, imagine if
the model is a huge 3D model, and one whole surface containing about 100
nodes is to be prescribed a boundary condition. If the nodes on this surface
were not grouped into node sets, then the user would end up with 100 x (No.
of DOFs to be constrained) data lines just to prescribe a boundary condition.
A more efficient way, of course, is to group the nodes in this particular surface
in a node set, and then like Case 3, write down the data lines referencing the
node set to be constrained. In this way, the number of data lines for the
‘BOUNDARY’ block would be equal to the number of DOFs to be constrained.
Similar use of sets can be applied to elements as well. One common use of
element sets (ELSET) is the referencing of element properties to the elements
in a particular set (see the example in Section 13.5.3). Sets are thus the basic
referencing tool in ABAQUS.
• A line can include no more than 256 characters including blanks.
• Parameter values are not usually case-sensitive. The only exceptions to
this rule are those imposed externally on ABAQUS, such as file names on
case-sensitive operating systems.
• If a parameter has a value, the equal sign (=) is used. The value can be an
integer, a floating point number, or a character string, depending on the
context. For example,
Data Lines
The data lines must immediately follow a keyword line if they are
required. The following rules apply when entering a data line:
• A data line can include no more than 256 characters including blanks.
• All data items are separated by commas. An empty data field is specified
by omitting data between commas. ABAQUS will use values of zeros for any
required numeric data that are omitted, unless there is a default value
allocated. If a data line contains only a single data item, the data item should
be followed by a comma.
• Floating point numbers can occupy a maximum of 20 spaces including the
sign, decimal point, and any exponential notation.
• Floating point numbers can be given with or without an exponent. Any
exponent, if input, must be preceded by E or D and an optional (-) or (+) to
indicate the sign of the exponent.
• Character strings can be up to 80 characters long and are not case-
sensitive.
• Continuation lines are allowed in specific instances, such as when defining
elements with a large number of nodes. If allowed, such lines are indicated by
a comma as the last character of the preceding line.
Labels
Examples of labels are set names, surface names and material names, and
they are case-sensitive, unlike the other entries in the keyword line. Labels
can be up to 80 characters long. All spaces within a label are ignored unless
the label is enclosed in quotation marks, in which case all spaces within the
label are maintained. A label that is not enclosed within quotation marks may
not include a period (.), and should not contain characters such as commas
and equal signs. If a label is defined using quotation marks, the quotation
marks are stored as part of the label. Any subsequent reference or use of the
label should also include the quotation marks. Labels cannot begin and end
with a double underscore (e.g._ALU_).
This label format is reserved for use internally within ABAQUS.
Model Data
Some of the data that must be included in the model data are as follows:
• Geometry of the model: The geometry of the model is described by its
elements and nodes.
• Material definitions, which are usually associated with parts of the
geometry.
History Data
As mentioned, in the history data, the entries are divided into steps. Each
step begins with the *STEP option and ends with the *END STEP option.
There are generally two kinds of step that can be defined in ABAQUS – the
general response analysis steps (can be linear or nonlinear); and the linear
perturbation steps. A general analysis step is one in which the effects of any
nonlinearities present in the model can be included. The starting condition for
each general analysis step is the ending condition from the last general
analysis step. The response of each general analysis step contributes to the
overall history of the response of the model. A linear perturbation analysis
step, on the other hand, is used to calculate the linear perturbation response
from the base state. The base state is the present state of the model at the
end of the last general analysis response. For the perturbation step, the
response does not contribute to the history of the overall response, and hence
can be called for at any time in between general steps. For cases where the
general step or the linear perturbation step is the first step, then the initial
conditions defined will define the starting condition or the base state,
respectively. The following is a list of the analysis types that uses linear
perturbation procedures:
• *BUCKLE (Eigenvalue buckling prediction)
Except for the above analysis types and for the * STATIC (where both general
and perturbation steps can be used), all other analysis types are general
analysis steps.
Some of the data that must be included in the history data or within a
step are:
• Analysis type: an option to define the analysis procedure type which
ABAQUS will perform. This must appear immediately after the *STEP option.
Other optional data include:
• Loading: some form of external loading can be defined. Loadings can be
in the form of concentrated loads, distributed loads, thermal loads, and so on.
Loadings can also be prescribed as a function of time following the amplitude
curve defined in the model data. If an amplitude curve is not defined,
ABAQUS will assume that the loading varies linearly over the step (ramp
loading), or that the loading is applied instantaneously at the beginning of the
step (step loading).
• Boundary conditions: zero-valued or non-zero boundary conditions can
be added, modified or removed. Note that if defined in the model data, only
zero-valued and symmetrical boundary conditions can be included.
• Output control: controls the requested output from the analysis. Output
variables depend upon the type of analysis and the type of elements used.
• Auxiliary controls: options are provided to allow the user to overwrite the
solution controls that are built into ABAQUS.
• Element and surface removal/reactivation: portions of the model can be
removed or reactivated from step to step.
Next will be writing the model data. First, the nodes of the problem must be
defined, since elements must be made up of nodes, and both nodes and
elements make up the geometry of the problem.
Using the *NODE option, the nodes at the end are first defined. We then use
the option *NGEN to generate evenly distributed nodes between the first and
last nodes. *NGEN is one of the several mesh generation capabilities
provided by ABAQUS. We could also define all 11 nodes individually by
specifying their coordinates, but using *NGEN would be more efficient for
large problems. So we have now defined 11 nodes uniformly along the length
of the beam. Next, the elements would be defined:
Here, the *ELEMENT option is used to define the first element that consists of
nodes 1 and 2. The TYPE parameter is included to specify what type of
element is being defined. In this case, B23 refers to a 1D beam elment in a
plane with cubic interpolation. Users can refer to the ABAQUS manual for the
element library to check the codes for other element types. Similar to the
definition of nodes, *ELGEN is used to generate 1-10 elements subsequently.
The elements are then grouped into a node set called RECT_BEAM. This will
make the referencing of element properties much easier later. So we have
now defined 11 nodes and 10 elements as shown in Figure 13.2.
The next step will be to define the element properties:
In the *BEAM SECTION keyword line, the element set RECT_BEAM
defined earlier is now referenced, meaning that the elements grouped under
RECT_BEAM will all have the properties defined in this option block. We also
provide the information that the beam has a rectangular (RECT) cross-
section. There are other cross-sections available in ABAQUS, such as circular
cross-sections (CIRC), trapezoidal cross-sections (TRAPEZOID), closed thin-
walled sections (BOX, HEX and PIPE) and open thin-walled sections (I-
section and L-section). ABAQUS also provides for a ‘general’ cross-section by
specifying geometrical quantitites necessary to define the section. The
material associated with the elements is also defined as ALU, where the
properties will be defined later. It is a good time to note that, unlike most
programming languages, the ABAQUS input file need not follow a top-down
approach when ABAQUS is assessing the file. For example, the material ALU
is already referenced at this point under the *BEAM SECTION option block,
though its material properties are actually defined further down the input file.
There will not be any error stating that the material ALU is invalid regardless
of where the material is defined unless it is not defined at all throughout the
input file. This is true for all other entries into the input file. Let us now look at
the data lines. The first data line in the *BEAM SECTION basically defines the
dimensions of the cross-section (0.025 x 0.04 m). Note that the dimensions
here are converted to metres to be consistent with the coordinates of the
nodes. The second data line basically defines the direction cosines indicating
the local beam axis. What is given is the default values, and this line can
actually be omitted in this case.
The next entry in the model data would be the material properties
definition:
The material for our example is aluminium, and we name it ALU for short.
All the properties option block will follow after the *MATERIAL option block,
which does not require any data lines by itself. The *ELASTIC option defines
elastic properties, and TYPE=ISOTROPIC defines the material as an isotropic
material, i.e. the material properties are the same in all directions. The data
line for the * ELASTIC option includes the values for the Young’s modulus and
Poisson’s ratio. Depending upon the type of analysis carried out, or the type
of material being defined, other properties may need to be defined. For
example, if a dynamic analysis is required, then the *DENSITY option would
also need to be included; or when the material exhibits viscoelastic behaviour,
then the
* VISCOELASTIC option would be required.
At this point, we have almost completed describing the model in the model
data. What is left are the boundary conditions. Note that the boundary
conditions can also be defined in the history data. What can be defined in the
model data is only the zero valued conditions.
From what we learned from the finite element method, we can actually
deduce that certain output variables are direct nodal variables like
displacements, while others like stress and strain are actually determined as a
distribution in the element using the shape functions. In ABAQUS, this
difference is categorized into nodal output variables and element output
variables. *NODF PRINT outputs the results of the required nodal variables in
an ASCII text file (.dat file), while the *NODF FILE ouputs the results in a
binary format (.fil file). The binary format can be read by post-processors in
which the results can be displayed. Similarly, *ELFMENT PRINT outputs
element variables in ASCII format, while *F,LFMF,NT FILE outputs them in
binary format. A list of the different output variables can be obtained in the
ABAQUS manuals. For our case, U in the data lines for *NODF PRINT and
*NODF FILF will output all the components of the nodal displacements. S and
F represent all the components of stress and strain, respectively. So if the
analysis is run, there will be altogether three tables: one showing the nodal
displacements, one showing the stresses in the elements, and the last one
showing the strains in the elements. The last thing to do now is end the step
by including *FNDSTFP. If multiple steps are present, this would separate the
different steps in the history data.