def test_geometric(self):
self.section.calculate_geometric_properties()
val_list = []
val_list.append({"prop": "area", "val": 100 * 50, "tol": None})
val_list.append({"prop": "ea", "val": 1 * 100 * 50, "tol": None})
val_list.append({"prop": "qx", "val": 100 * 50 * 50, "tol": None})
val_list.append({"prop": "qy", "val": 100 * 50 * 25, "tol": None})
val_list.append({"prop": "ixx_g", "val": 50 * 100 ** 3 / 3,
"tol": None})
val_list.append({"prop": "iyy_g", "val": 100 * 50 ** 3 / 3,
"tol": None})
val_list.append({"prop": "ixy_g", "val": 100 * 50 * 50 * 25,
"tol": None})
val_list.append({"prop": "cx", "val": 50 / 2, "tol": None})
val_list.append({"prop": "cy", "val": 100 / 2, "tol": None})
val_list.append({"prop": "ixx_c", "val": 50 * 100 ** 3 / 12,
"tol": None})
val_list.append({"prop": "iyy_c", "val": 100 * 50 ** 3 / 12,
"tol": None})
val_list.append({"prop": "ixy_c", "val": 0, "tol": None})
val_list.append({"prop": "zxx_plus", "val": 50 * 100 ** 2 / 6,
"tol": None})
val_list.append({"prop": "zxx_minus", "val": 50 * 100 ** 2 / 6,
"tol": None})
val_list.append({"prop": "zyy_plus", "val": 100 * 50 ** 2 / 6,
"tol": None})
val_list.append({"prop": "zyy_minus", "val": 100 * 50 ** 2 / 6,
"tol": None})
val_list.append({"prop": "rx",
"val": (50 * 100 ** 3 / 12 / 100 / 50) ** 0.5,
"tol": None})
val_list.append({"prop": "ry",
"val": (100 * 50 ** 3 / 12 / 100 / 50) ** 0.5,
"tol": None})
val_list.append({"prop": "i11_c", "val": 50 * 100 ** 3 / 12,
"tol": None})
val_list.append({"prop": "i22_c", "val": 100 * 50 ** 3 / 12,
"tol": None})
val_list.append({"prop": "phi", "val": 0, "tol": None})
val_list.append({"prop": "z11_plus", "val": 50 * 100 ** 2 / 6,
"tol": None})
val_list.append({"prop": "z11_minus", "val": 50 * 100 ** 2 / 6,
"tol": None})
val_list.append({"prop": "z22_plus", "val": 100 * 50 ** 2 / 6,
"tol": None})
val_list.append({"prop": "z22_minus", "val": 100 * 50 ** 2 / 6,
"tol": None})
val_list.append({"prop": "r11",
"val": (50 * 100 ** 3 / 12 / 100 / 50) ** 0.5,
"tol": None})
val_list.append({"prop": "r22",
"val": (100 * 50 ** 3 / 12 / 100 / 50) ** 0.5,
"tol": None})
validate_properties(self, val_list, self.section)
def test_plastic(self):
self.section.calculate_geometric_properties()
self.section.calculate_plastic_properties()
val_list = []
val_list.append({"prop": "x_pc", "val": 50 / 2, "tol": None})
val_list.append({"prop": "y_pc", "val": 100 / 2, "tol": None})
val_list.append({"prop": "x11_pc", "val": 50 / 2, "tol": None})
val_list.append({"prop": "y22_pc", "val": 100 / 2, "tol": None})
val_list.append({"prop": "sxx", "val": 50 * 100**2 / 4, "tol": None})
val_list.append({"prop": "syy", "val": 100 * 50**2 / 4, "tol": None})
val_list.append({"prop": "s11", "val": 50 * 100**2 / 4, "tol": None})
val_list.append({"prop": "s22", "val": 100 * 50**2 / 4, "tol": None})
val_list.append({"prop": "sf_xx_plus", "val": 1.5, "tol": None})
val_list.append({"prop": "sf_xx_minus", "val": 1.5, "tol": None})
val_list.append({"prop": "sf_yy_plus", "val": 1.5, "tol": None})
val_list.append({"prop": "sf_yy_minus", "val": 1.5, "tol": None})
val_list.append({"prop": "sf_11_plus", "val": 1.5, "tol": None})
val_list.append({"prop": "sf_11_minus", "val": 1.5, "tol": None})
val_list.append({"prop": "sf_22_plus", "val": 1.5, "tol": None})
val_list.append({"prop": "sf_22_minus", "val": 1.5, "tol": None})
validate_properties(self, val_list, self.section)
def test_angle(self):
"""Section properties are validated against results from the Strand7
beam section utility."""
geometry = sections.AngleSection(d=150,
b=90,
t=12,
r_r=10,
r_t=5,
n_r=8)
mesh = geometry.create_mesh(mesh_sizes=[2.5])
section = CrossSection(geometry, mesh)
val_list = []
val_list.append({"prop": "area", "val": 2746.73, "tol": 2e-4})
val_list.append({"prop": "perimeter", "val": 471, "tol": 1e-3})
val_list.append({"prop": "cx", "val": 21.2255, "tol": 2e-4})
val_list.append({"prop": "cy", "val": 50.9893, "tol": 2e-4})
val_list.append({"prop": "ixx_g", "val": 1.3428e7, "tol": 2e-4})
val_list.append({"prop": "iyy_g", "val": 2.95629e6, "tol": 2e-4})
val_list.append({"prop": "ixy_g", "val": 1.08669e6, "tol": 2e-4})
val_list.append({"prop": "ixx_c", "val": 6.28678e6, "tol": 2e-4})
val_list.append({"prop": "iyy_c", "val": 1.71882e6, "tol": 3e-4})
val_list.append({"prop": "ixy_c", "val": -1.88603e6, "tol": 3e-4})
val_list.append({"prop": "zxx_plus", "val": 63496.0, "tol": 2e-4})
val_list.append({"prop": "zxx_minus", "val": 123296, "tol": 2e-4})
val_list.append({"prop": "zyy_plus", "val": 24992.1, "tol": 3e-4})
val_list.append({"prop": "zyy_minus", "val": 80979.0, "tol": 2e-4})
val_list.append({"prop": "rx", "val": 47.8416, "tol": 2e-4})
val_list.append({"prop": "ry", "val": 25.0154, "tol": 2e-4})
val_list.append({"prop": "i11_c", "val": 6.96484e6, "tol": 2e-4})
val_list.append({"prop": "i22_c", "val": 1.04076e6, "tol": 2e-4})
val_list.append({"prop": "phi", "val": 19.7744 - 180, "tol": 2e-4})
val_list.append({"prop": "z11_plus", "val": 97751.9, "tol": 2e-4})
val_list.append({"prop": "z11_minus", "val": 69403.3, "tol": 2e-4})
val_list.append({"prop": "z22_plus", "val": 27959.0, "tol": 2e-4})
val_list.append({"prop": "z22_minus", "val": 20761.6, "tol": 3e-4})
val_list.append({"prop": "r11", "val": 50.3556, "tol": 2e-4})
val_list.append({"prop": "r22", "val": 19.4656, "tol": 2e-4})
val_list.append({"prop": "sxx", "val": 113541, "tol": 2e-4})
val_list.append({"prop": "syy", "val": 45724.6, "tol": 2e-4})
val_list.append({
"prop": "sf_xx_plus",
"val": 113541 / 63496.0,
"tol": 2e-4
})
val_list.append({
"prop": "sf_xx_minus",
"val": 113541 / 123296,
"tol": 2e-4
})
val_list.append({
"prop": "sf_yy_plus",
"val": 45724.6 / 24992.1,
"tol": 3e-4
})
val_list.append({
"prop": "sf_yy_minus",
"val": 45724.6 / 80979.0,
"tol": 2e-4
})
val_list.append({"prop": "s11", "val": 121030, "tol": 2e-4})
val_list.append({"prop": "s22", "val": 43760.6, "tol": 2e-4})
val_list.append({
"prop": "sf_11_plus",
"val": 121030 / 97751.9,
"tol": 2e-4
})
val_list.append({
"prop": "sf_11_minus",
"val": 121030 / 69403.3,
"tol": 2e-4
})
val_list.append({
"prop": "sf_22_plus",
"val": 43760.6 / 27959.0,
"tol": 2e-4
})
val_list.append({
"prop": "sf_22_minus",
"val": 43760.6 / 20761.6,
"tol": 3e-4
})
val_list.append({"prop": "j", "val": 135333, "tol": 1e-3})
val_list.append({"prop": "gamma", "val": 1.62288e8, "tol": 5e-4})
val_list.append({"prop": "A_s11", "val": 885.444, "tol": 2e-4})
val_list.append({"prop": "A_s22", "val": 1459.72, "tol": 4e-4})
val_list.append({"prop": "x11_se", "val": 28.719, "tol": 1e-3})
val_list.append({"prop": "y22_se", "val": 35.2348, "tol": 5e-4})
section.calculate_geometric_properties()
section.calculate_plastic_properties()
section.calculate_warping_properties()
validate_properties(self, val_list, section)
def test_custom(self):
"""Section properties are validated against results from the Strand7
beam section utility."""
points = [[-10, 0], [110, 0], [100, 10], [55, 10], [55, 90], [100, 90],
[110, 100], [110, 110], [-10, 110], [-10, 100], [0, 90],
[45, 90], [45, 10], [-10, 10]]
facets = [[0, 1], [1, 2], [2, 3], [3, 4], [4, 5], [5, 6], [6,
7], [7, 8],
[8, 9], [9, 10], [10, 11], [11, 12], [12, 13], [13, 0]]
holes = []
control_points = [[0, 5]]
geometry = sections.CustomSection(points, facets, holes,
control_points)
mesh = geometry.create_mesh(mesh_sizes=[5])
section = CrossSection(geometry, mesh)
val_list = []
val_list.append({"prop": "area", "val": 4250, "tol": None})
val_list.append({"prop": "cx", "val": 49.3333, "tol": None})
val_list.append({"prop": "cy", "val": 65.0196, "tol": None})
val_list.append({"prop": "ixx_g", "val": 2.56725e7, "tol": None})
val_list.append({"prop": "iyy_g", "val": 1.41858e7, "tol": None})
val_list.append({"prop": "ixy_g", "val": 1.37979e7, "tol": None})
val_list.append({"prop": "ixx_c", "val": 7.70542e6, "tol": None})
val_list.append({"prop": "iyy_c", "val": 3.84228e6, "tol": None})
val_list.append({"prop": "ixy_c", "val": 165472, "tol": None})
val_list.append({"prop": "zxx_plus", "val": 171306, "tol": None})
val_list.append({"prop": "zxx_minus", "val": 118509, "tol": None})
val_list.append({"prop": "zyy_plus", "val": 63334.2, "tol": None})
val_list.append({"prop": "zyy_minus", "val": 64757.5, "tol": None})
val_list.append({"prop": "rx", "val": 42.5798, "tol": None})
val_list.append({"prop": "ry", "val": 30.0677, "tol": None})
val_list.append({"prop": "phi", "val": 177.552 - 180, "tol": 1e-4})
val_list.append({"prop": "i11_c", "val": 7.71249e6, "tol": None})
val_list.append({"prop": "i22_c", "val": 3.8352e6, "tol": None})
val_list.append({"prop": "z11_plus", "val": 162263, "tol": None})
val_list.append({"prop": "z11_minus", "val": 114268, "tol": None})
val_list.append({"prop": "z22_plus", "val": 60503.0, "tol": None})
val_list.append({"prop": "z22_minus", "val": 62666.1, "tol": None})
val_list.append({"prop": "r11", "val": 42.5993, "tol": None})
val_list.append({"prop": "r22", "val": 30.04, "tol": None})
val_list.append({"prop": "sxx", "val": 153196, "tol": None})
val_list.append({"prop": "syy", "val": 101494, "tol": None})
val_list.append({
"prop": "sf_xx_plus",
"val": 153196 / 171306,
"tol": None
})
val_list.append({
"prop": "sf_xx_minus",
"val": 153196 / 118509,
"tol": None
})
val_list.append({
"prop": "sf_yy_plus",
"val": 101494 / 63334.2,
"tol": None
})
val_list.append({
"prop": "sf_yy_minus",
"val": 101494 / 64757.5,
"tol": None
})
val_list.append({"prop": "s11", "val": 153347, "tol": None})
val_list.append({"prop": "s22", "val": 101501, "tol": None})
val_list.append({
"prop": "sf_11_plus",
"val": 153347 / 162263,
"tol": None
})
val_list.append({
"prop": "sf_11_minus",
"val": 153347 / 114268,
"tol": None
})
val_list.append({
"prop": "sf_22_plus",
"val": 101501 / 60503.0,
"tol": None
})
val_list.append({
"prop": "sf_22_minus",
"val": 101501 / 62666.1,
"tol": None
})
val_list.append({"prop": "j", "val": 347040, "tol": 5e-3})
val_list.append({"prop": "gamma", "val": 7.53539e9, "tol": 1e-3})
val_list.append({"prop": "A_s11", "val": 2945.53, "tol": 5e-4})
val_list.append({"prop": "A_s22", "val": 956.014, "tol": 5e-4})
val_list.append({"prop": "x11_se", "val": 1.9134, "tol": 5e-3})
val_list.append({"prop": "y22_se", "val": 3.02028, "tol": 5e-3})
section.calculate_geometric_properties()
section.calculate_plastic_properties()
section.calculate_warping_properties()
validate_properties(self, val_list, section)