def test_distal_forces(self):
j = FreeJoint('joint')
# Constant loading
j.distal_forces = np.array([2, 3.1, 2.1, 4.3, 2.5, 1.0])
j.calc_external_loading()
assert_array_equal(j.applied_forces[:6], 0)
assert_array_equal(j.applied_forces[6:], j.distal_forces)
assert_array_equal(j.applied_stress, 0)
# Loading function
j.distal_forces = lambda element, t: np.ones(6)
j.calc_external_loading()
assert_array_equal(j.applied_forces[:6], 0)
assert_array_equal(j.applied_forces[6:], 1)
assert_array_equal(j.applied_stress, 0)
def test_interal_forces(self):
# NB minus sign because convention for applied_stress is that
# stiffness loads are positive.
j = FreeJoint('joint')
# Constant loading
j.internal_forces = np.array([2, 3.1, 2.1, 4.3, 2.5, 1.0])
j.calc_external_loading()
assert_array_equal(j.applied_forces, 0)
assert_array_equal(j.applied_stress, -j.internal_forces)
# Loading function
j.internal_forces = lambda element, t: np.ones(6)
j.calc_external_loading()
assert_array_equal(j.applied_forces, 0)
assert_array_equal(j.applied_stress, -1)
def test_interal_forces(self):
# NB minus sign because convention for applied_stress is that
# stiffness loads are positive.
j = FreeJoint('joint')
# Constant loading
j.internal_forces = np.array([2, 3.1, 2.1, 4.3, 2.5, 1.0])
j.calc_external_loading()
assert_array_equal(j.applied_forces, 0)
assert_array_equal(j.applied_stress, -j.internal_forces)
# Loading function
j.internal_forces = lambda element, t: np.ones(6)
j.calc_external_loading()
assert_array_equal(j.applied_forces, 0)
assert_array_equal(j.applied_stress, -1)
