section_force_diagram_3d

opsvis.section_force_diagram_3d(sf_type, Ew, sfac=1.0, nep=17, fmt_secforce='b-', dir_plt=0)[source]

Display section forces diagram of a 3d beam column model.

This function plots section forces diagrams for 3d beam column elements with or without element loads. For now only ‘-beamUniform’ constant transverse or axial element loads are supported.

Parameters

  • sf_type (str) – type of section force: ‘N’ - normal force, ‘Vy’ or ‘Vz’ - shear force, ‘My’ or ‘Mz’ - bending moments, ‘T’ - torsional moment.

  • Ew (dict) – Ew Python dictionary contains information on non-zero element loads, therfore each item of the Python dictionary is in the form: ‘ele_tag: [‘-beamUniform’, Wy, Wz, Wx]’.

  • sfac (float) – scale factor by wich the values of section forces are multiplied.

  • nep (int) – number of evaluation points including both end nodes (default: 17)

  • fmt_secforce (str) – format line string for section force distribution curve. The format contains information on line color, style and marks as in the standard matplotlib plot function. (default: fmt_secforce = ‘b-’ # blue solid line)

  • {0 (dir_plt) – direction in which to plot the load effects: 0 (default) - as defined in the code for each load effect type 1 - in the y-axis (default for N, Vy, T, Mz) 2 - in the z-axis (default for Vz, My)

  • 1 – direction in which to plot the load effects: 0 (default) - as defined in the code for each load effect type 1 - in the y-axis (default for N, Vy, T, Mz) 2 - in the z-axis (default for Vz, My)

  • 2} – direction in which to plot the load effects: 0 (default) - as defined in the code for each load effect type 1 - in the y-axis (default for N, Vy, T, Mz) 2 - in the z-axis (default for Vz, My)

Usage:

Wy, Wz, Wx = -5., 0., 0.
Ew = {3: ['-beamUniform', Wy, Wz, Wx]}
sfacMz = 1.e-1
plt.figure()
minY, maxY = opsv.section_force_diagram_3d('Mz', Ew, sfacMz)
plt.title(f'Bending moments Mz, max = {maxY:.2f}, min = {minY:.2f}')

Todo:

Add support for other element loads available in OpenSees: partial (trapezoidal) uniform element load, and ‘beamPoint’ element load.

See example Statics of a 3d 3-element cantilever beam