Excavation Supported by Cantilevered Sheet Pile Wall

This example is from file Excavation Supported by Cantilevered Sheet Pile Wall on the OpenSees website and has been converted using opst.pre.tcl2py The python model script can be found here excavation.py

import openseespy.opensees as ops

Load the FEM model function FEMmodel from file excavation.py. [excavation.py](https://github.com/yexiang1992/opstool/blob/master/docs/examples/post/excavation/excavation.py)

from utils.excavation import FEMmodel

import opstool as opst
import opstool.vis.pyvista as opsvis

FEMmodel()

opsvis.set_plot_props(point_size=1, font_size=9)
opsvis.plot_model(show_node_numbering=True, show_ele_numbering=True).show()

Static Scene

Interactive Scene

Create output database (ODB) file. Since some elements and nodes will be removed in subsequent analyses, ensure that model_update=True.

ODB = opst.post.CreateODB(
    odb_tag=1,
    model_update=True,
    compute_mechanical_measures=True,
    project_gauss_to_nodes="copy",
)

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:179: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:179: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:179: UserWarning:

eleTag 3 have no fiber sec!

GRAVITY ANALYSIS (w/ INITIAL STATE ANALYSIS TO RESET DISPLACEMENTS)

define analysis parameters for gravity phase

ops.constraints("Transformation")
ops.test("NormDispIncr", 1e-05, 50, 0)
ops.algorithm("Newton")
ops.numberer("RCM")
ops.system("BandGeneral")
ops.integrator("LoadControl", 1)
ops.analysis("Static")

Perform an initial state analysis, where elements with tags 1001–1042 are [BeamContact2D](https://opensees.berkeley.edu/wiki/index.php?title=BeamContact2D).

# turn on initial state analysis feature
ops.InitialStateAnalysis("on")
# ensure soil material intially considers linear elastic behavior
ops.updateMaterialStage("-material", 1, "-stage", 0)
# set contact elements to be frictionless for gravity analysis
ops.setParameter("-val", 0, "-eleRange", 1001, 1042, "friction")
# analysis 4 steps, and fetch response
for _ in range(4):
    ops.analyze(1)
    ODB.fetch_response_step()

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

Update soil material to consider elastoplastic behavior and analyze a few more steps:

# update soil material to consider elastoplastic behavior and analyze a few more steps
ops.updateMaterialStage("-material", 1, "-stage", 1)
# analysis 4 steps, and fetch response
for _ in range(4):
    ops.analyze(1)
    ODB.fetch_response_step()
# designate end of initial state analysis (zeros displacements, keeps state variables)
ops.InitialStateAnalysis("off")
# turn on frictional behavior for beam contact elements
ops.setParameter("-val", 1, "-eleRange", 1001, 1042, "friction")

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

REMOVE ELEMENTS TO SIMULATE EXCAVATION

# define analysis parameters for excavation phase
ops.wipeAnalysis()
ops.constraints("Transformation")
ops.test("NormDispIncr", 0.0001, 60)
ops.algorithm("KrylovNewton")
ops.numberer("RCM")
ops.system("BandGeneral")
ops.integrator("LoadControl", 1)
ops.analysis("Static")

We first define a function to avoid repetitive removal of elements and nodes, and then proceed with several steps of analysis.

def remove_components(ele_tags, node_tags, nsteps=4):
    for etag in ele_tags:
        ops.remove("element", etag)
    for ntag in node_tags:
        ops.remove("node", ntag)
    # run analysis after object removal
    for _ in range(nsteps):
        ops.analyze(1)
        ODB.fetch_response_step()

Remove objects associated with lift 1:

# soil elements
ele_tags = [191, 192, 193, 194, 195, 196, 197, 198, 199, 200]
ele_tags += [1042]  # contact element
# soil nodes
node_tags = [430, 437, 446, 455, 461, 468, 473, 476, 480, 482, 484]
node_tags += [1042]  # lagrange multiplier node

remove_components(ele_tags, node_tags, nsteps=4)

print("Lift 1 removed")

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

Lift 1 removed

We can then remove the remaining 9 lifts：

# remove objects associated with lift 2
# soil elements
ele_tags = [181, 182, 183, 184, 185, 186, 187, 188, 189, 190]
ele_tags += [1040]  # contact element
# soil nodes
node_tags = [412, 424, 433, 444, 453, 460, 466, 471, 475, 479, 483]
node_tags += [1040]  # lagrange multiplier node

remove_components(ele_tags, node_tags, nsteps=4)

print("Lift 2 removed")

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

Lift 2 removed

remove objects associated with lift 3

# soil elements
ele_tags = [171, 172, 173, 174, 175, 176, 177, 178, 179, 180]
ele_tags += [1038]  # contact element
# soil nodes
node_tags = [387, 405, 418, 429, 442, 450, 458, 464, 470, 477, 481]
node_tags += [1038]  # lagrange multiplier node

remove_components(ele_tags, node_tags, nsteps=4)

print("Lift 3 removed")

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

Lift 3 removed

remove objects associated with lift 4

# soil elements
ele_tags = [161, 162, 163, 164, 165, 166, 167, 168, 169, 170]
ele_tags += [1036]  # contact element
# soil nodes
node_tags = [363, 380, 398, 414, 427, 439, 448, 457, 465, 472, 478]
node_tags += [1036]  # lagrange multiplier node

remove_components(ele_tags, node_tags, nsteps=4)

print("Lift 4 removed")

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

Lift 4 removed

remove objects associated with lift 5

# soil elements
ele_tags = [151, 152, 153, 154, 155, 156, 157, 158, 159, 160]
ele_tags += [1034]  # contact element
# soil nodes
node_tags = [336, 353, 378, 395, 411, 425, 440, 449, 459, 467, 474]
node_tags += [1034]  # lagrange multiplier node

remove_components(ele_tags, node_tags, nsteps=4)

print("Lift 5 removed")

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

Lift 5 removed

remove objects associated with lift 6

# soil elements
ele_tags = [141, 142, 143, 144, 145, 146, 147, 148, 149, 150]
ele_tags += [1032]  # contact element
# soil nodes
node_tags = [308, 326, 347, 369, 392, 408, 426, 441, 452, 462, 469]
node_tags += [1032]  # lagrange multiplier node

remove_components(ele_tags, node_tags, nsteps=4)

print("Lift 6 removed")

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

Lift 6 removed

remove objects associated with lift 7

# soil elements
ele_tags = [131, 132, 133, 134, 135, 136, 137, 138, 139, 140]
ele_tags += [1030]  # contact element
# soil nodes
node_tags = [281, 304, 322, 345, 370, 394, 415, 428, 443, 454, 463]
node_tags += [1030]  # lagrange multiplier node

remove_components(ele_tags, node_tags, nsteps=4)

print("Lift 7 removed")

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

Lift 7 removed

remove objects associated with lift 8

# soil elements
ele_tags = [121, 122, 123, 124, 125, 126, 127, 128, 129, 130]
ele_tags += [1028]  # contact element
# soil nodes
node_tags = [260, 278, 302, 325, 346, 374, 399, 419, 434, 447, 456]
node_tags += [1028]  # lagrange multiplier node

remove_components(ele_tags, node_tags, nsteps=4)

print("Lift 8 removed")

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

Lift 8 removed

remove objects associated with lift 9

# soil elements
ele_tags = [111, 112, 113, 114, 115, 116, 117, 118, 119, 120]
ele_tags += [1026]  # contact element
# soil nodes
node_tags = [241, 253, 277, 306, 327, 350, 379, 406, 422, 438, 451]
node_tags += [1026]  # lagrange multiplier node

remove_components(ele_tags, node_tags, nsteps=4)

print("Lift 9 removed")

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

Lift 9 removed

remove objects associated with lift 10

# soil elements
ele_tags = [101, 102, 103, 104, 105, 106, 107, 108, 109, 110]
ele_tags += [1024]  # contact element
# soil nodes
node_tags = [218, 239, 259, 282, 307, 335, 364, 389, 413, 431, 445]
node_tags += [1024]  # lagrange multiplier node

remove_components(ele_tags, node_tags, nsteps=4)

print("Lift 10 removed")

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 1 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 4 have no fiber sec!

G:\opstool\opstool\post\_get_response\_get_fiber_sec_resp.py:314: UserWarning:

eleTag 3 have no fiber sec!

Lift 10 removed

We can save all previous responses to a file: zlib compression is used to reduce file size.

ODB.save_response(zlib=True)

OPSTOOL ::  All responses data with _odb_tag = 1 saved in G:\opstool\docs\.opstool.output/RespStepData-1.nc!

Post-processing

import matplotlib.pyplot as plt

import opstool as opst
import opstool.vis.pyvista as opsvis

Since the result data has already been saved, we can read it at any time for post-processing:

opsvis.set_plot_props(point_size=0, line_width=5, cmap="turbo")
opsvis.set_plot_props(
    scalar_bar_kargs={
        "label_font_size": 12,
        "title_font_size": 13,
        "position_x": 0.85,  # 0--1
    }
)

Nodal responses

opsvis.plot_nodal_responses(
    odb_tag=1,
    slides=True,
    defo_scale=20,
    resp_type="disp",
    resp_dof=["UX", "UY"],
    unit_symbol="m",
).show()

Static Scene

Interactive Scene

OPSTOOL ::  Loading response data from G:\opstool\docs\.opstool.output/RespStepData-1.nc ...

We can create animations:

fig = opsvis.plot_nodal_responses_animation(
    odb_tag=1,
    framerate=20,
    defo_scale=25,
    savefig="images/NodalRespAnimation-excavation.gif",
    resp_type="disp",
    resp_dof=["UX", "UY"],
    unit_symbol="m",
)
fig.close()

OPSTOOL ::  Loading response data from G:\opstool\docs\.opstool.output/RespStepData-1.nc ...
Animation has been saved to images/NodalRespAnimation-excavation.gif!

Frame elements responses

plotter = opsvis.plot_frame_responses(
    odb_tag=1,
    resp_type="sectionForces",
    resp_dof="MZ",
    unit_symbol="kN·m",
    show_values="eleMaxMin",
    scale=3,
    slides=True,
    style="surface",
    show_model=False,  # plot all model
    opacity=1.0,
    show_bc=False,
)
plotter.show()

Static Scene

Interactive Scene

OPSTOOL ::  Loading response data from G:\opstool\docs\.opstool.output/RespStepData-1.nc ...

opsvis.plot_frame_responses_animation(
    odb_tag=1,
    resp_type="sectionForces",
    resp_dof="MZ",
    unit_symbol="kN·m",
    show_values=False,
    framerate=20,
    scale=3,
    style="surface",
    opacity=1.0,
    show_model=True,  # plot all model
    show_bc=False,
    savefig="images/FrameForcesMZ-excavation.gif",
).close()

OPSTOOL ::  Loading response data from G:\opstool\docs\.opstool.output/RespStepData-1.nc ...
Animation has been saved as images/FrameForcesMZ-excavation.gif!

Plane elements response

pl = opsvis.plot_unstruct_responses(
    odb_tag=1, slides=True, ele_type="Plane", resp_type="StressesAtNodes", resp_dof="sigma22", unit_symbol="kPa"
)
pl.show()

Static Scene

Interactive Scene

OPSTOOL ::  Loading response data from G:\opstool\docs\.opstool.output/RespStepData-1.nc ...

opsvis.set_plot_props(show_mesh_edges=False)
opsvis.plot_unstruct_responses(
    odb_tag=1,
    slides=True,
    ele_type="Plane",
    resp_type="StressesAtNodes",
    resp_dof="sigma12",
    show_defo=True,
    defo_scale=30,
    unit_symbol="MPa",
).show()

Static Scene

Interactive Scene

OPSTOOL ::  Loading response data from G:\opstool\docs\.opstool.output/RespStepData-1.nc ...

opsvis.plot_unstruct_responses(
    odb_tag=1, slides=True, ele_type="Plane", resp_type="StressesAtNodes", resp_dof="sigma_vm"
).show()

Static Scene

Interactive Scene

OPSTOOL ::  Loading response data from G:\opstool\docs\.opstool.output/RespStepData-1.nc ...

Read data from ODB

Reading the response of the contact element

data = opst.post.get_element_responses(odb_tag=1, ele_type="Contact")
data

OPSTOOL ::  Loading Contact response data from G:\opstool\docs\.opstool.output/RespStepData-1.nc ...

<xarray.Dataset> Size: 91kB
Dimensions:       (time: 49, eleTags: 42, globalDOFs: 3, localDOFs: 3,
                   slipDOFs: 2)
Coordinates:
  * eleTags       (eleTags) int32 168B 1001 1002 1003 1004 ... 1040 1041 1042
  * globalDOFs    (globalDOFs) <U2 24B 'Px' 'Py' 'Pz'
  * localDOFs     (localDOFs) <U2 24B 'N' 'Tx' 'Ty'
  * slipDOFs      (slipDOFs) <U2 16B 'Tx' 'Ty'
  * time          (time) float32 196B 0.0 1.0 2.0 3.0 ... 37.0 38.0 39.0 40.0
Data variables:
    globalForces  (time, eleTags, globalDOFs) float32 25kB -0.0 -0.0 ... nan nan
    localForces   (time, eleTags, localDOFs) float32 25kB 0.0 0.0 ... nan nan
    localDisp     (time, eleTags, localDOFs) float32 25kB 0.0 0.0 ... nan nan
    slips         (time, eleTags, slipDOFs) float32 16kB 0.0 0.0 0.0 ... nan nan
Attributes:
    Px:       Global force in the x-direction on the constrained node
    Py:       Global force in the y-direction on the constrained node
    Pz:       Global force in the z-direction on the constrained node
    N:        Normal force or deformation
    Tx:       Tangential force or deformation in the x-direction
    Ty:       Tangential force or deformation in the y-direction

xarray.Dataset

• Dimensions:
  • time: 49
  • eleTags: 42
  • globalDOFs: 3
  • localDOFs: 3
  • slipDOFs: 2
• Coordinates: (5)
  • eleTags
    (eleTags)
    int32
    1001 1002 1003 ... 1040 1041 1042

    array([1001, 1002, 1003, 1004, 1005, 1006, 1007, 1008, 1009, 1010, 1011, 1012,
           1013, 1014, 1015, 1016, 1017, 1018, 1019, 1020, 1021, 1022, 1023, 1024,
           1025, 1026, 1027, 1028, 1029, 1030, 1031, 1032, 1033, 1034, 1035, 1036,
           1037, 1038, 1039, 1040, 1041, 1042])

  • globalDOFs
    (globalDOFs)
    <U2
    'Px' 'Py' 'Pz'

    array(['Px', 'Py', 'Pz'], dtype='<U2')

  • localDOFs
    (localDOFs)
    <U2
    'N' 'Tx' 'Ty'

    array(['N', 'Tx', 'Ty'], dtype='<U2')

  • slipDOFs
    (slipDOFs)
    <U2
    'Tx' 'Ty'

    array(['Tx', 'Ty'], dtype='<U2')

  • time
    (time)
    float32
    0.0 1.0 2.0 3.0 ... 38.0 39.0 40.0

    array([ 0.,  1.,  2.,  3.,  4.,  5.,  6.,  7.,  8.,  1.,  2.,  3.,  4.,  5.,
            6.,  7.,  8.,  9., 10., 11., 12., 13., 14., 15., 16., 17., 18., 19.,
           20., 21., 22., 23., 24., 25., 26., 27., 28., 29., 30., 31., 32., 33.,
           34., 35., 36., 37., 38., 39., 40.], dtype=float32)

• Data variables: (4)
  • globalForces
    (time, eleTags, globalDOFs)
    float32
    -0.0 -0.0 0.0 -0.0 ... nan nan nan

    array([[[-0.0000000e+00, -0.0000000e+00,  0.0000000e+00],
            [-0.0000000e+00, -0.0000000e+00,  0.0000000e+00],
            [-0.0000000e+00, -0.0000000e+00,  0.0000000e+00],
            ...,
            [-0.0000000e+00, -0.0000000e+00,  0.0000000e+00],
            [-0.0000000e+00, -0.0000000e+00,  0.0000000e+00],
            [-0.0000000e+00, -0.0000000e+00,  0.0000000e+00]],
    
           [[-2.2277697e+01,  1.2898215e-16,  0.0000000e+00],
            [ 2.2277697e+01, -1.2898215e-16,  0.0000000e+00],
            [-4.3412949e+01, -0.0000000e+00,  0.0000000e+00],
            ...,
            [ 2.2848918e+00, -1.6235133e-14,  0.0000000e+00],
            [-5.7122296e-01, -4.0587833e-15,  0.0000000e+00],
            [ 5.7122296e-01, -4.0587833e-15,  0.0000000e+00]],
    
           [[-2.2277697e+01,  3.4214551e-16,  0.0000000e+00],
            [ 2.2277697e+01, -3.4214551e-16,  0.0000000e+00],
            [-4.3412949e+01,  2.1286442e-16,  0.0000000e+00],
            ...,
    ...
            ...,
            [           nan,            nan,            nan],
            [-4.1343015e-02,  9.7618526e-04,  0.0000000e+00],
            [           nan,            nan,            nan]],
    
           [[-2.6858183e+01,  1.3562502e-01,  0.0000000e+00],
            [ 8.8400278e+00, -1.0248453e-01,  0.0000000e+00],
            [-3.9127411e+01,  3.6174440e-01,  0.0000000e+00],
            ...,
            [           nan,            nan,            nan],
            [-4.1335627e-02,  9.7796472e-04,  0.0000000e+00],
            [           nan,            nan,            nan]],
    
           [[-2.6858181e+01,  1.3562506e-01,  0.0000000e+00],
            [ 8.8400249e+00, -1.0248456e-01,  0.0000000e+00],
            [-3.9127403e+01,  3.6174452e-01,  0.0000000e+00],
            ...,
            [           nan,            nan,            nan],
            [-4.1335266e-02,  9.7783026e-04,  0.0000000e+00],
            [           nan,            nan,            nan]]], dtype=float32)

  • localForces
    (time, eleTags, localDOFs)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 nan nan nan

    array([[[ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00],
            [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00],
            [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00],
            ...,
            [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00],
            [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00],
            [ 0.0000000e+00,  0.0000000e+00,  0.0000000e+00]],
    
           [[ 2.2277697e+01, -2.5796431e-16,  0.0000000e+00],
            [ 2.2277697e+01,  2.5796431e-16,  0.0000000e+00],
            [ 4.3412949e+01,  0.0000000e+00,  0.0000000e+00],
            ...,
            [ 2.2848918e+00,  0.0000000e+00,  0.0000000e+00],
            [ 5.7122296e-01,  0.0000000e+00,  0.0000000e+00],
            [ 5.7122296e-01,  0.0000000e+00,  0.0000000e+00]],
    
           [[ 2.2277697e+01, -6.8429102e-16,  0.0000000e+00],
            [ 2.2277697e+01,  6.8429102e-16,  0.0000000e+00],
            [ 4.3412949e+01, -4.2572884e-16,  0.0000000e+00],
            ...,
    ...
            ...,
            [           nan,            nan,            nan],
            [ 4.1345987e-02, -1.6815730e-03,  0.0000000e+00],
            [           nan,            nan,            nan]],
    
           [[ 2.6858387e+01, -1.7246141e-01,  0.0000000e+00],
            [ 8.8402014e+00,  1.7245381e-01,  0.0000000e+00],
            [ 3.9128014e+01, -5.7866025e-01,  0.0000000e+00],
            ...,
            [           nan,            nan,            nan],
            [ 4.1338608e-02, -1.6851805e-03,  0.0000000e+00],
            [           nan,            nan,            nan]],
    
           [[ 2.6858385e+01, -1.7246147e-01,  0.0000000e+00],
            [ 8.8401985e+00,  1.7245387e-01,  0.0000000e+00],
            [ 3.9128006e+01, -5.7866043e-01,  0.0000000e+00],
            ...,
            [           nan,            nan,            nan],
            [ 4.1338246e-02, -1.6849140e-03,  0.0000000e+00],
            [           nan,            nan,            nan]]], dtype=float32)

  • localDisp
    (time, eleTags, localDOFs)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 nan nan nan

    array([[[ 0.,  0.,  0.],
            [ 0.,  0.,  0.],
            [ 0.,  0.,  0.],
            ...,
            [ 0.,  0.,  0.],
            [ 0.,  0.,  0.],
            [ 0.,  0.,  0.]],
    
           [[ 0.,  0.,  0.],
            [ 0.,  0.,  0.],
            [ 0.,  0.,  0.],
            ...,
            [ 0.,  0.,  0.],
            [ 0.,  0.,  0.],
            [ 0.,  0.,  0.]],
    
           [[ 0.,  0.,  0.],
            [ 0.,  0.,  0.],
            [ 0.,  0.,  0.],
            ...,
    ...
            ...,
            [nan, nan, nan],
            [ 0.,  0.,  0.],
            [nan, nan, nan]],
    
           [[ 0.,  0.,  0.],
            [ 0.,  0.,  0.],
            [ 0.,  0.,  0.],
            ...,
            [nan, nan, nan],
            [ 0.,  0.,  0.],
            [nan, nan, nan]],
    
           [[ 0.,  0.,  0.],
            [ 0.,  0.,  0.],
            [ 0.,  0.,  0.],
            ...,
            [nan, nan, nan],
            [ 0.,  0.,  0.],
            [nan, nan, nan]]], dtype=float32)

  • slips
    (time, eleTags, slipDOFs)
    float32
    0.0 0.0 0.0 0.0 ... 0.0 0.0 nan nan

    array([[[ 0.,  0.],
            [ 0.,  0.],
            [ 0.,  0.],
            ...,
            [ 0.,  0.],
            [ 0.,  0.],
            [ 0.,  0.]],
    
           [[ 0.,  0.],
            [ 0.,  0.],
            [ 0.,  0.],
            ...,
            [ 0.,  0.],
            [ 0.,  0.],
            [ 0.,  0.]],
    
           [[ 0.,  0.],
            [ 0.,  0.],
            [ 0.,  0.],
            ...,
    ...
            ...,
            [nan, nan],
            [ 0.,  0.],
            [nan, nan]],
    
           [[ 0.,  0.],
            [ 0.,  0.],
            [ 0.,  0.],
            ...,
            [nan, nan],
            [ 0.,  0.],
            [nan, nan]],
    
           [[ 0.,  0.],
            [ 0.,  0.],
            [ 0.,  0.],
            ...,
            [nan, nan],
            [ 0.,  0.],
            [nan, nan]]], dtype=float32)

• Indexes: (5)
  • eleTags
    PandasIndex

    PandasIndex(Index([1001, 1002, 1003, 1004, 1005, 1006, 1007, 1008, 1009, 1010, 1011, 1012,
           1013, 1014, 1015, 1016, 1017, 1018, 1019, 1020, 1021, 1022, 1023, 1024,
           1025, 1026, 1027, 1028, 1029, 1030, 1031, 1032, 1033, 1034, 1035, 1036,
           1037, 1038, 1039, 1040, 1041, 1042],
          dtype='int32', name='eleTags'))

  • globalDOFs
    PandasIndex

    PandasIndex(Index(['Px', 'Py', 'Pz'], dtype='object', name='globalDOFs'))

  • localDOFs
    PandasIndex

    PandasIndex(Index(['N', 'Tx', 'Ty'], dtype='object', name='localDOFs'))

  • slipDOFs
    PandasIndex

    PandasIndex(Index(['Tx', 'Ty'], dtype='object', name='slipDOFs'))

  • time
    PandasIndex

    PandasIndex(Index([ 0.0,  1.0,  2.0,  3.0,  4.0,  5.0,  6.0,  7.0,  8.0,  1.0,  2.0,  3.0,
            4.0,  5.0,  6.0,  7.0,  8.0,  9.0, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0,
           16.0, 17.0, 18.0, 19.0, 20.0, 21.0, 22.0, 23.0, 24.0, 25.0, 26.0, 27.0,
           28.0, 29.0, 30.0, 31.0, 32.0, 33.0, 34.0, 35.0, 36.0, 37.0, 38.0, 39.0,
           40.0],
          dtype='float32', name='time'))

• Attributes: (6)

  Px :

  Global force in the x-direction on the constrained node

  Py :

  Global force in the y-direction on the constrained node

  Pz :

  Global force in the z-direction on the constrained node

  N :

  Normal force or deformation

  Tx :

  Tangential force or deformation in the x-direction

  Ty :

  Tangential force or deformation in the y-direction

data["localForces"].sel(eleTags=1001).plot.line(x="time")
plt.show()

Let’s examine the response of contact element #1034. Since it is removed during the fifth lift, its response is truncated at time=16, and subsequent data will be filled with numpy.nan.

data["localForces"].sel(eleTags=1034).plot.line(x="time")
plt.show()

data["localForces"].sel(eleTags=1034).data

Reading the response of the beam element

data = opst.post.get_element_responses(odb_tag=1, ele_type="Frame")
print(data)

OPSTOOL ::  Loading Frame response data from G:\opstool\docs\.opstool.output/RespStepData-1.nc ...
<xarray.Dataset> Size: 309kB
Dimensions:              (time: 49, eleTags: 21, localDofs: 12, basicDofs: 6,
                          secPoints: 3, secDofs: 6, locs: 3)
Coordinates:
  * eleTags              (eleTags) int32 84B 401 402 403 404 ... 418 419 420 421
  * localDofs            (localDofs) <U3 144B 'FX1' 'FY1' 'FZ1' ... 'MY2' 'MZ2'
  * basicDofs            (basicDofs) <U3 72B 'N' 'MZ1' 'MZ2' 'MY1' 'MY2' 'T'
  * secPoints            (secPoints) int32 12B 1 2 3
  * secDofs              (secDofs) <U2 48B 'N' 'MZ' 'VY' 'MY' 'VZ' 'T'
  * locs                 (locs) <U5 60B 'alpha' 'X' 'Y'
  * time                 (time) float32 196B 0.0 1.0 2.0 3.0 ... 38.0 39.0 40.0
Data variables:
    localForces          (time, eleTags, localDofs) float32 49kB -0.0 ... 0.0...
    basicForces          (time, eleTags, basicDofs) float32 25kB 0.0 ... 0.0
    basicDeformations    (time, eleTags, basicDofs) float32 25kB 0.0 ... 0.0
    plasticDeformation   (time, eleTags, basicDofs) float32 25kB 0.0 ... 0.0
    sectionForces        (time, eleTags, secPoints, secDofs) float32 74kB 0.0...
    sectionDeformations  (time, eleTags, secPoints, secDofs) float32 74kB 0.0...
    sectionLocs          (time, eleTags, secPoints, locs) float32 37kB 0.1127...
Attributes:
    localDofs:  local coord system dofs at end 1 and end 2
    basicDofs:  basic coord system dofs at end 1 and end 2
    secPoints:  section points No.
    secDofs:    section forces and deformations Dofs. Note that the section D...
    Notes:      Note that the deformations are displacements and rotations in...

data = opst.post.get_element_responses(odb_tag=1, ele_type="Plane")
print(data)

OPSTOOL ::  Loading Plane response data from G:\opstool\docs\.opstool.output/RespStepData-1.nc ...
<xarray.Dataset> Size: 7MB
Dimensions:                (time: 49, eleTags: 400, GaussPoints: 4,
                            stressDOFs: 5, strainDOFs: 3, nodeTags: 462,
                            measures: 4)
Coordinates:
  * eleTags                (eleTags) int32 2kB 1 2 3 4 5 ... 396 397 398 399 400
  * GaussPoints            (GaussPoints) int32 16B 1 2 3 4
  * stressDOFs             (stressDOFs) <U7 140B 'sigma11' 'sigma22' ... 'eta_r'
  * strainDOFs             (strainDOFs) <U5 60B 'eps11' 'eps22' 'eps12'
  * nodeTags               (nodeTags) int32 2kB 1 2 3 4 5 ... 481 482 483 484
  * time                   (time) float32 196B 0.0 1.0 2.0 ... 38.0 39.0 40.0
  * measures               (measures) <U8 128B 'p1' 'p2' 'sigma_vm' 'tau_max'
Data variables:
    Stresses               (time, eleTags, GaussPoints, stressDOFs) float32 2MB ...
    Strains                (time, eleTags, GaussPoints, strainDOFs) float32 941kB ...
    StressesAtNodes        (time, nodeTags, stressDOFs) float32 453kB 0.0 ......
    StrainsAtNodes         (time, nodeTags, strainDOFs) float32 272kB 0.0 ......
    StressAtNodesErr       (time, nodeTags, stressDOFs) float32 453kB 0.0 ......
    StrainsAtNodesErr      (time, nodeTags, strainDOFs) float32 272kB 0.0 ......
    StressMeasures         (time, eleTags, GaussPoints, measures) float32 1MB ...
    StrainMeasures         (time, eleTags, GaussPoints, measures) float32 1MB ...
    StressMeasuresAtNodes  (time, nodeTags, measures) float32 362kB 0.0 ... nan
    StrainMeasuresAtNodes  (time, nodeTags, measures) float32 362kB 0.0 ... nan
Attributes:
    p1, p2:                     Principal stresses (strains).
    sigma11, sigma22, sigma12:  Normal stress and shear stress (strain) in th...
    eta_r:                      Ratio between the shear (deviatoric) stress a...
    sigma_vm:                   Von Mises stress.
    tau_max:                    Maximum shear stress (strains).