2222# Define the Imports
2323# ==================
2424# Here we bring in the primative section shapes and also the more generic
25- # Shapely `Polygon` object. We also use the nifty `tqdm` package to provide
26- # a progress bar on computations.
25+ # Shapely `Polygon` object.
2726import numpy as np
2827import matplotlib .pyplot as plt
2928from shapely .geometry import Polygon
3837# It's better to collect the relevant section property calculation in a
3938# single function. We are only interested in the torsion constant, so this
4039# is straightforward enough. `geom` is the Section Property Gemoetry object;
41- # `ms` is the mesh size. The "cgs" solver helps avoid singular matrices due
42- # to small angles.
40+ # `ms` is the mesh size.
4341def get_section_j (geom , ms , plot_geom = False ):
4442 geom .create_mesh (mesh_sizes = [ms ])
4543 section = Section (geom )
@@ -59,7 +57,7 @@ def get_section_j(geom, ms, plot_geom=False):
5957n = 100 # mesh density
6058
6159#%%
62- # Create and Analyze the Section
60+ # Create and Analyse the Section
6361# ==============================
6462# This function accepts the width `b` and a slope `S` to create the trapezoid.
6563# Since we are only interested in relative results, the nominal dimensions are
@@ -101,7 +99,7 @@ def do_section(b, S, d_mid=1, plot_geom=False):
10199#%%
102100# Create Loop Variables
103101# =====================
104- # # The slope `S` is 0 for a rectangle, and 1 for a triangle and is
102+ # The slope `S` is 0 for a rectangle, and 1 for a triangle and is
105103# defined per the plot below. A range of `S`, between 0.0 and 1.0 and
106104# a range of `b` are considered, between 1 and 10 here (but can be extended)
107105b_list = np .logspace (0 , np .log10 (10 ), 10 )
0 commit comments