Parametric Shape Variation

Shape parametrization can be accomplished in an Abaqus input file by:

  • parametrizing nodal coordinates; or

  • relating nodal coordinates to shape parameters using shape variations.

The different approaches to shape parametrization are described in this section.

This page discusses:

Parametrization of Nodal Coordinates

Any individual nodal coordinates can be parametrized directly. This is usually of limited value because it often leads to designs with irregular shape that cannot be manufactured easily. In addition, parametrization of individual nodal coordinates generally requires an excessive number of parameters to define the parametrized shape.

Parametrization of nodal coordinates used in conjunction with node generation in Abaqus provides a more practical method of shape parametrization. However, this method is still of somewhat limited practical use because the simple node generation capabilities available in Abaqus cannot describe complex shapes.

Direct Parametrization of Individual Nodal Coordinates

The simplest form of parametrization of nodal coordinates is to define individual parameters and use them in place of the nodal coordinates to be parametrized, as described in Parametric Modeling. For example,

PARAMETER
x_coord_node_1 = 10.
y_coord_node_1 = 20.
NODE
1, <x_coord_node_1>, <y_coord_node_1>

Parametrization of Nodal Coordinates Using Node Generation

Shape parametrization can be accomplished by parametrizing the coordinates of some nodes, then using these nodes to generate other nodes and their coordinates. For example:

PARAMETER
x_coord_node_1 = 10.
x_coord_node_11 = 20.
NODE
1, <x_coord_node_1>, 50.
11, <x_coord_node_11>, 50.
NGEN
1, 11, 1

This method of shape parametrization reduces the number of user-defined parameters necessary for shape parametrization by implicitly making the nodal coordinates of the generated nodes dependent on the shape parameters.

Shape Change by Linear Combination of Shape Variations

The definition of shape in Abaqus includes a basic shape plus any number of additional shape variations that are added to the basic shape using a linear combination. Mathematically, we can express the nodal coordinates, x, as

x=x0+i=1Nsipi,

where x0 is the basic shape, si is the ith shape variation, and pi is the value of the ith shape parameter. This calculation is always done in the global rectangular Cartesian coordinate system. Although it is not necessarily so, it is frequently the case that the input to define a shape variation is simply the gradient of the basic shape x0 taken with respect to the corresponding shape parameter.

You specify the basic shape of a model in the Abaqus input file by providing nodal definitions either directly or through node generation; see Node Definition.

You can specify shape variations and associated shape parameters, as described here.*PARAMETER SHAPE VARIATION

In addition, you can specify perturbations of the shape as a linear combination of other shapes (for example, buckling mode shapes); see Introducing a Geometric Imperfection into a Model.

The definition of the nodal coordinates for a model in the Abaqus input file is then possible using a combination of four types of methods:

  • You can directly define individual nodes and their respective coordinates; these coordinates are part of the definition of the basic shape, x0, and can be parametrized.

  • Node generation can be used to create nodes and their coordinates according to geometrically simple mappings that rely on existing node definitions; these generated coordinates are also part of the definition of the basic shape, x0. If necessary, the node generation input can be parametrized.

  • Parameter shape variations can be used to vary the coordinates of nodes defined using the above methods.

  • Geometric imperfections can be used to perturb nodal coordinates previously defined using any combination of the above three types of methods.

Shape Parametrization Using Shape Variations

Instead of parametrizing nodal coordinates directly, you can specify shape variations. Each shape variation must be associated with a single shape parameter. The names of the parameters associated with the shape variations must be chosen such that the names remain unique when interpreted in a case-insensitive manner. The values of the shape parameters are assigned using parameter definitions.

A parameter shape variation can be defined more than once for the same parameter so that different parts of a shape variation can be specified separately. In these cases if the same node is specified in multiple parameter shape variation definitions, the last definition for the node prevails.

A node that is specified under a parameter shape variation definition that has not also been defined directly or through node generation will be ignored.

You can specify shape variations using a combination of three possibilities: directly specifying them, reading them from an alternate input file, and reading them from the results files of auxiliary analyses. These methods are described in the following sections.

Defining Shape Variations Directly or Reading Them from an Alternate Input File

You can define the shape variation data directly by specifying the node number and corresponding variations of coordinate components. Alternatively, the data can be given in an ASCII file.

Defining Shape Variations in Alternative Coordinate Systems

By default, the shape variation data are interpreted in the global rectangular Cartesian coordinate system. You can specify the shape variation data (either directly or in an alternate input file) in cylindrical or spherical coordinate systems. In such cases the computation of the shape variation is done as follows. The nodal coordinate components that define the basic shape are first transformed from the global rectangular Cartesian coordinate system in which they are stored to the specified coordinate system. The shape variation coordinate components are then added to give updated coordinate components, which are transformed back to the global rectangular Cartesian coordinate system. Finally, the shape variation is taken as the difference between the updated coordinate components and the original coordinate components, using the components expressed in the global rectangular Cartesian coordinate system. The value of the shape parameter associated with the shape variation is not used at any point in the calculation of the shape variation.

Using Auxiliary Analyses to Generate Shape Variations

Auxiliary models are additional finite element models that are used to generate shape variations for a primary model. Rather than defining shape variations directly on a node-by-node basis, auxiliary models can be used to simplify this process. Auxiliary analyses are finite element analyses of these auxiliary models.

An auxiliary model usually has the same geometry, element connectivity, and material type as the primary model. However, the boundary conditions are usually different. Applying loading to an auxiliary model results in sets of displacements that we may interpret as shape variations. For example, we may be interested in studying the sensitivity of the nonlinear buckling behavior of a structure with respect to imperfections in the structure. In this case we could perform an auxiliary eigenvalue linear buckling analysis and then use the resulting mode shapes as shape variations to be added to the basic geometry of the primary model. (This particular problem could also be addressed by using a geometric imperfection.)

Abaqus reads the shape variation data from auxiliary analyses through the user node labels. Abaqus does not check model compatibility between both analysis runs. Shape variation data cannot be read from the results file for models defined in terms of an assembly of part instances (Assembly Definition).

Reading Shape Variations from a Static Analysis Results File

To define a shape variation based on the deformed geometry of a previous static analysis, specify the results file and step from a previous static analysis. Optionally, you can specify the increment number from which displacement data are read. (By default, Abaqus will read data from the last increment available for the specified step on the results file.) In addition, you can read shape variation data for a specified node set.

Reading Shape Variations from an Eigenvalue Analysis Results File

To define a shape variation based on a mode shape from a previous eigenvalue analysis, specify the results file and step from a previous eigenfrequency extraction or eigenvalue buckling prediction analysis. Optionally, you can specify the mode number from which eigenvector data are read. (By default, Abaqus will read data from the first eigenvector available for the specified step on the results file.) In addition, you can read eigenmode data for a specified node set.

Shape Parametrization and Design Sensitivity Analysis

For the purpose of design sensitivity analysis with Abaqus/Design (Direct Design Sensitivity Analysis) if the parameter specified for a parameter shape variation is also specified as a design parameter, the shape variation is used to define the design gradient of the nodal coordinates and nodal normals with respect to the design parameter. If you wish to perform design sensitivity analysis for the basic shape, all shape parameters must be given a value of zero. In addition, if any parameter specified in a parameter shape variation definition is also specified as a design parameter, the parameters of all parameter shape variations must be specified as design parameters.

In DSA calculations for shell and beam elements Abaqus always computes the design gradients of nodal normals using the design gradients of nodal coordinates. To overwrite the gradients computed by Abaqus, you must provide the nodal normal as part of the node definition and design gradients of the normals using a parameter shape variation. To prescribe a design-independent normal, you must provide a zero design gradient explicitly. For shape variations read from the results file, Abaqus computes the gradients of the normals based on the displacements and ignores the nodal rotations.

For beam elements Abaqus computes the design gradients for the n1-direction of the beam cross-section using the gradients of the node coordinates and the gradients for the n2-direction specified using a parameter shape variation. You cannot provide the shape variation for the n1-direction. Abaqus ignores any such design gradients implicitly provided in either the beam section definition or as an extra node in the beam element connectivity.

In cases where the data defining a shape variation are given in a cylindrical or spherical coordinate system it is important that you understand how the shape variation is calculated from the data. This calculation is described in the previous section.