robbievanleeuwen / section-properties

Analysis of an arbitrary cross-section in python using the finite element method.
https://sectionproperties.rtfd.io
MIT License
412 stars 92 forks source link

MatrixRankWarning for angle section with thickness equal to toe radius #421

Closed Ernst79 closed 2 months ago

Ernst79 commented 4 months ago

Issue

First of all, thanks for this wonderful piece of work.

I have a small issue, when creating and analyzing an angle section.

from sectionproperties.pre.library import angle_section
from sectionproperties.analysis import Section

geom = angle_section(d=20, b=20, t=3, r_r=3.5, r_t=3, n_r=16)
geom.create_mesh(mesh_sizes=[20 / 20])

# create mesh
sec = Section(geometry=geom)

# analyse mesh
sec.calculate_geometric_properties()
sec.calculate_warping_properties()
sec.calculate_plastic_properties()

ixx_c, iyy_c, ixy_c = sec.get_ic()
zxx_plus, zxx_minus, zyy_plus, zyy_minus = sec.get_z()

print("area:", int(sec.get_area()))
print("Ixx:", int(ixx_c))
print("Iyy:", int(iyy_c))
print("Ixy:", int(ixy_c))

The above code return the following warning

image

The problem only occurs if the toe radius is the same as the thickness of the angle, probably due to the vertical part of the bottom flange is getting a length of t - r_t = 0.00000.

github-actions[bot] commented 4 months ago

Thanks for opening your first issue in sectionproperties :raised_hands: Pull requests are always welcome :wink:

robbievanleeuwen commented 2 months ago

Hi @Ernst79, thanks for submitting the issue. This has been fixed with #429. There will be a new release shortly that incorporates this fix!

Ernst79 commented 2 months ago

Thanks!