Get the section properties of build-up sections

A forum dedicated to feature requests and the future direction of OpenSees, i.e. what would you like, what do you need, what should we explore

Moderators: silvia, selimgunay, Moderators

Post Reply
Anran
Posts: 1
Joined: Mon May 01, 2023 8:03 am

Get the section properties of build-up sections

Post by Anran »

Hi Community,
I have created an L-shape section using the patch command. I want to know the section properties of the created section, such as the moment of inertia. One solution is to do a static analysis and compare the obtained displacement with the theoretical one. But, I'm curious if there is a command to get the build-up fiber section properties.
Attached is my code for the build-up L-shape
import matplotlib.pyplot as plt
import numpy as np
import openseespy.opensees as ops
import openseespy.postprocessing.ops_vis as opsv
# Angle Section
y1 = np.array([-0.98, -0.98, -0.73, -0.73])
z1 = np.array([-0.237, 1.263, 1.263, -0.237])
y2 = np.array([-0.98, -0.98, 2.02, 2.02])
z2 = np.array([-0.487, -0.237, -0.237, -0.487])
xy1 = np.c_[y1, z1].ravel()
xy2 = np.c_[y2, z2].ravel()

fib_sec_3 = [['section', 'Fiber', 3, '-GJ', 1.0e6],
['patch', 'quad', 1, 8, 4, *xy1], # noqa: E501
['patch', 'quad', 1, 4, 12, *xy2], # noqa: E501
]
matcolor = ['r', 'lightgrey', 'gold', 'w', 'w', 'w']
opsv.plot_fiber_section(fib_sec_3, matcolor=matcolor)
plt.axis('equal')
plt.show()


Thanks,
Ran
Post Reply