def test_multiple_loads(self):
"""
Case 2: A truss with multiple loads on different joints.
"""
joints = ((0, 0), (5, 0), (0, 4), (2, 6))
bars = (('AC', strong), ('BC', strong), ('CD', strong), ('BD', strong))
loads = [('C', 1000, 0), ('D', 0, -750)]
supports = (('A', {'support_type': 'pin'}), ('B', {'support_type': 'pin'}))
tc.init_truss('Multiple loads')
results = build_from_lists(joints, bars, loads, supports,
sig_figs=3, solution_method=tc.SolveMethod.NUMPY_STD)
tc.plot_diagram(tc.active_truss, results, show_reactions=True)
def test_multiple_loads_on_same_joint(self):
"""
Case 5: A truss with mutliple loads on the same joint which do not cancel out.
"""
joints = ((0, 1), (1, 0), (1, 1))
bars = (('AB', strong), ('BC', strong), ('AC', strong))
loads = [('C', 1, 2), ('C', 0, -1)]
supports = (('A', {'support_type': 'pin', 'pin_rotation': 90}),
('B', {'support_type': 'roller', 'roller_normal': (-1, 2)})
)
tc.init_truss('Multiple loads on the same joint')
results = build_from_lists(joints, bars, loads, supports,
sig_figs=3, solution_method=tc.SolveMethod.NUMPY_STD)
tc.plot_diagram(tc.active_truss, results, show_reactions=True)
def test_unloaded_truss(self):
"""
Case 4: A truss without any applied external loads.
"""
joints = ((0, 1), (1, 0), (1, 1))
bars = (('AB', strong), ('BC', strong), ('AC', strong))
loads = []
supports = (('A', {'support_type': 'pin', 'pin_rotation': 90}),
('B', {'support_type': 'roller', 'roller_normal': (-1, 2)})
)
tc.init_truss('Completely unloaded truss')
results = build_from_lists(joints, bars, loads, supports,
sig_figs=3, solution_method=tc.SolveMethod.NUMPY_STD)
tc.plot_diagram(tc.active_truss, results, show_reactions=True)
def test_with_angled_roller(self):
"""
Case 3: A truss with a roller support at an inclined angle.
"""
joints = ((0, 1), (1, 0), (1, 1))
bars = (('AB', strong), ('BC', strong), ('AC', strong))
loads = [('C', 1, 2)]
supports = (('A', {'support_type': 'pin', 'pin_rotation': 90}),
('B', {'support_type': 'roller', 'roller_normal': (-1, 2)})
)
tc.init_truss('Angled roller support')
results = build_from_lists(joints, bars, loads, supports,
sig_figs=3, solution_method=tc.SolveMethod.NUMPY_STD)
tc.plot_diagram(tc.active_truss, results, show_reactions=True)
def test_with_fully_cancelling_loads(self):
"""
Case 6: A truss with multiple loads on the same joint which do cancel out
giving an effectively unloaded truss.
"""
joints = ((0, 1), (1, 0), (1, 1))
bars = (('AB', strong), ('BC', strong), ('AC', strong))
loads = [('C', 1, 2), ('C', -1, -2)]
supports = (('A', {'support_type': 'pin', 'pin_rotation': 90}),
('B', {'support_type': 'roller', 'roller_normal': (-1, 2)})
)
tc.init_truss('All external loads cancel out')
results = build_from_lists(joints, bars, loads, supports,
sig_figs=3, solution_method=tc.SolveMethod.NUMPY_STD)
tc.plot_diagram(tc.active_truss, results, show_reactions=True)
def test_bridge_structure(self):
"""
Case 1b: A truss bridge.
"""
joints = ((0, 0), (100, 0), (200, 0), (300, 0), (400, 0), (100, 100), (200, 100), (300, 100))
bars = (('AB', strong), ('BC', strong), ('CD', strong), ('DE', strong), ('AF', strong),
('BF', strong), ('CF', strong), ('CG', strong), ('CH', strong), ('DH', strong),
('EH', strong), ('FG', strong), ('GH', strong))
loads = [('A', 0, -100), ('B', 0, -200), ('C', 0, -200), ('D', 0, -200), ('E', 0, -100)]
supports = (('A', {'support_type': 'pin'}), ('E', {'support_type': 'roller', 'roller_normal': (0, 1)}))
tc.init_truss('Bridge')
results = build_from_lists(joints, bars, loads, supports,
sig_figs=3, solution_method=tc.SolveMethod.NUMPY_STD)
tc.active_truss.dump_truss_to_json(filedir='../Saved Trusses')
tc.plot_diagram(tc.active_truss, results, show_reactions=True)
def build_from_lists(joints: tuple[tuple[float]], bars: tuple[tuple[str, dict]],
loads: list[tuple[str, float, float]], supports: tuple[tuple[str, dict]], **res_kwargs):
"""
Allows quick construction of full trusses given lists in an appropriate format:
joints: `((x1, y1), (x2, y2), ...)` named in order by default A, B, C, ...
bars: `(('AB', strong), ('BC', weak), ...)` bar_name is a two-char string,
using the letters from the joints to indicate which ones it goes between
loads: `[('A', x1, y1), ('C', x2, y2), ...]` joint_name is a one-char string, indicating the loaded joint
supports: `(('A', kwargs1), ('B', kwargs2), ...)`
joint_name is a one-char string, indicating the supported joint kwargs can be a dict which fills any of
the following:
`{'support_type': 'pin'/'roller'/'encastre', 'pin_rotation': angle_degrees, 'roller_normal': (x, y)}`
"""
import string
_alpha = string.ascii_uppercase
_nums = string.digits
for i, (x, y) in enumerate(joints):
tc.create_joint('Joint ' + _alpha[i], x, y)
for i, (c, bar_type) in enumerate(bars):
tc.create_bar('Bar ' + c, 'Joint ' + c[0], 'Joint ' + c[1], bar_type)
for i, (j, x, y) in enumerate(loads):
tc.create_load(r'$W_{' + _nums[i] + r'}$', 'Joint ' + j, x, y)
for i, (j, kwargs) in enumerate(supports):
tc.create_support('Support ' + j, 'Joint ' + j, **kwargs)
try: # Get the results of the truss calculation and display graphic
my_results = tc.active_truss.Result(tc.active_truss, **res_kwargs)
print(my_results)
return my_results
except tc.np.linalg.LinAlgError as e: # The truss was badly made, so could not be solved
tc.active_truss.classify_error_in_truss(e)
def test_SDC_truss(self):
"""
Case 1a: standard well-built truss. Uses the module factory functions.
Represents the truss built in the SDC project.
"""
tc.init_truss('SDC: Steel Cantilever')
tc.create_joint('Joint A', 0, 0)
tc.create_joint('Joint B', 290, -90)
tc.create_joint('Joint C', 815, 127.5)
tc.create_joint('Joint D', 290, 345)
tc.create_joint('Joint E', 0, 255)
tc.create_joint('Joint F', 220.836, 127.5)
tc.create_bar('Bar AB', 'Joint A', 'Joint B', medium_2)
tc.create_bar('Bar BC', 'Joint B', 'Joint C', strong)
tc.create_bar('Bar CD', 'Joint C', 'Joint D', medium_1)
tc.create_bar('Bar DE', 'Joint D', 'Joint E', medium_1)
tc.create_bar('Bar EF', 'Joint E', 'Joint F', medium_1)
tc.create_bar('Bar AF', 'Joint F', 'Joint A', medium_2)
tc.create_bar('Bar DF', 'Joint F', 'Joint D', medium_1)
tc.create_bar('Bar BF', 'Joint F', 'Joint B', weak)
tc.create_load('W', 'Joint C', 0, -0.675 * 1)
tc.create_support('Support A', 'Joint A', support_type='encastre')
tc.create_support('Support E', 'Joint E', support_type='pin', pin_rotation=90)
try: # Get the results of the truss calculation and display graphic
my_results = tc.active_truss.Result(tc.active_truss, sig_figs=3,
solution_method=tc.SolveMethod.NUMPY_STD)
print(my_results)
except tc.np.linalg.LinAlgError as e: # The truss was badly made, so could not be solved
tc.active_truss.classify_error_in_truss(e)
tc.active_truss.dump_truss_to_json(filedir='../Saved Trusses')
tc.plot_diagram(tc.active_truss, my_results, show_reactions=True)