Validation#
sectionproperties
has a growing suite of tests. The testing suite serves to
verify functionality and find exceptions from edge cases, but also validate the results
against known cases. Most importantly, it ensures that whenever any changes are made to
the code base, the expected output is preserved. Further, since this library performs
engineering calculations, it should have some manner of proving its accuracy. Each
analyst who uses it is responsible for their own projects and calculations, but having a
high level of confidence that the software can produce correct results,
given the correct inputs, is beneficial to all users. Some test results and
explanations from the latter category will be outlined on this page.
Textbook Examples#
An obvious starting location is replicating examples from academic texts. The main
reference text for sectionproperties
is “Analysis and Design of Elastic Beams” by
Walter D. Pilkey [Pilkey]. This text contains several examples and validation tests
which are covered in the testing suite, see
source.
Selected validation tests from [Pilkey] are found in the examples section of the
documentation, notably the Pilkey - Symmetric Channel Section,
Pilkey - Circular Arc and Pilkey - Composite Rectangular Strip
examples.
“Aircraft Structures” by David J. Peery is also a highly regarded text in the field of aerospace structures [Peery]. Several examples from [Peery] are described below, and can also be found in the Peery - Beams in Complex Bending example.
Peery - Symmetric Sections#
The simplest example of a realistic section problem is the symmetric I-Beam, with a free-fixed boundary condition and a transverse tip load. The free-body-diagram and shear and moment diagrams are shown in the problem statement, referenced below. This problem is Example 1 in Section 6.2.
Check #1#
In sectionproperties
, there are multiple ways to set this problem up. We could
different shapely geometries and merge together, or a set of custom points, or a
built-in constructor. For the sake of simplicity, this simpler I-section is identical to
the Nastran I-section definition, so it makes sense to utilise the built-in constructor
from nastran_i()
.
Using an arbitrarily coarse mesh, the properties can then be directly calculated from
the class method
calculate_geometric_properties()
.
Peery lists the second moment of area about the primary bending axis as a value of 43.3 in4. For the automated tests in this library, we check against this hardcoded value, with a tolerance of \(\pm 0.1\)%.
Check #2#
As a final check against this example, we can calculate the maximum bending stress on the I-beam. From simple statics, the maximum moment from the FBD will be 800,000 in-lbs at the fixed end. Applying this moment to our section from before will allow computation of stress over the FEM.
Peery quotes the peak value at 55.5 ksi, which is rounded to the nearest decimal place. From the equatio listed in the text, the theoretical value is actually 55,427.3 psi.
Again, the automated test against this checks the hardcoded value with a tolerance of \(\pm 0.1\)%. For accuracy, 55,427.3 is used instead of the rounded value.
For full details and the most updated code of this example, see the Peery - Beams in Complex Bending example. For the exact test code execution, check the source.
Peery - Unsymmetric Sections#
For a more complex example, we can turn to Example 1 in Section 7.2 of Peery. Here, we
have a still-simplified Z-section, but bending about two axes. Note axes definitions in
the problem statement. Beam axial direction in sectionproperties
is always
referenced as the z-axis, and loads must be applied in this coordinate system.
The construction of this geometry takes a similar approach to Ex 6.2.1, and utilises a
built-in factory, nastran_zed()
.
The only difference you may notice in the test code is usage of a custom class for ease
of initialisation. This is not necessary.
Using an arbitrarily coarse mesh, the properties can then be directly calculated from
the class method
calculate_geometric_properties()
.
Each property listed directly by Peery is taken as a hardcoded value and checked
against, within the testing suite.
Property |
Peery Value |
---|---|
I_x |
693.3 [in4] |
I_y |
173.3 [in4] |
I_xy |
-240 [in4] |
I_p |
787.1 [in4] |
I_q |
79.5 [in4] |
theta |
21.35 [deg] |
For stress results, the theoretical values follow the biaxial bending equation. These
values are checked against automatically in the testing suite. Note that again Peery
rounds the values quoted directly, for simplicity. The testing suite also verifies that
the theoretical value as per the equation matches the theoretical value quoted in the
text, which also matches the computed value from the sectionproperties
FEM.
Point |
x |
y |
-494x |
-315y |
\(f_b\), [psi] |
---|---|---|---|---|---|
A |
-5 |
4 |
2470 |
-1260 |
1210.0 = 1210 |
B |
-5 |
6 |
2470 |
-1890 |
580.0 = 580 |
C |
1 |
6 |
-494 |
-1890 |
-2384.0 = -2380 |
For full details and the most updated code of this example, see the Peery - Beams in Complex Bending example. For the exact test code execution, check the source.