def test_free_crimped(): material = Material(E='27000 MPa', poisson=1, alpha=1) section = Rectangle(25, 44.6) f1 = Load.PontualLoad(-5, x=29) f2 = Load.UniformDistributedLoad(-0.1622, x_begin=0, x_end=113) n1 = Node.Free(x=0) n2 = Node.Crimp(x=113) bar1 = BeamElement([n1, n2], section, material) beam = Beam(loads=[f1, f2], beam_elements=[bar1]) compare(beam=beam, name="free_crimped", ftool_displacement_in_mm=True)
def test_structural__crimped_simple_supported_beam(): material = Material(E=1, poisson=0.3, alpha=1) section = Rectangle(12, 1) f1 = Load.PontualLoad(-200, x=700) n1 = Node.Crimp(x=0) n2 = Node.SimpleSupport(x=1000) bar1 = BeamElement([n1, n2], section, material) beam = Beam(loads=[f1], beam_elements=[bar1]) support_reaction = beam.nodal_efforts assert support_reaction[0] == approx(87.3, abs=0.1) assert support_reaction[1] == approx(27300, abs=10) assert support_reaction[-2] == approx(112.7, abs=0.1) assert support_reaction[-1] == approx(0, abs=10) assert beam.getInternalShearStrength(100) == approx(87.3, abs=0.1) assert beam.getInternalShearStrength(800) == approx(-112.7, abs=0.1) assert beam.getInternalMomentumStrength(e) == approx(-27300, abs=10) assert beam.getInternalMomentumStrength(100) == approx(-18570, abs=10) assert beam.getInternalMomentumStrength(800) == approx(+22540, abs=10) assert beam.getInternalMomentumStrength(1000 - e) == approx(0, abs=10)
def test_v47(): material = Material(E='27000 MPa', poisson=1, alpha=1) section = Rectangle(25, 44.6) f1 = Load.UniformDistributedLoad(-0.1622, x_begin=0, x_end=113) f2 = Load.UniformDistributedLoad(-0.4994, x_begin=113, x_end=583) f3 = Load.UniformDistributedLoad(-0.4196, x_begin=583, x_end=1188) n1 = Node.SimpleSupport(x=0) n2 = Node.SimpleSupport(x=113) n3 = Node.SimpleSupport(x=583) n4 = Node.SimpleSupport(x=1188) bar1 = BeamElement([n1, n2], section, material) bar2 = BeamElement([n2, n3], section, material) bar3 = BeamElement([n3, n4], section, material) beam = Beam(loads=[f1, f2, f3], beam_elements=[bar1, bar2, bar3]) compare(beam=beam, name="v47", ftool_displacement_in_mm=True)
def test_structural_double_crimped_beam(): material = Material(E=1, poisson=0.3, alpha=1) section = Rectangle(12, 1) f1 = Load.PontualLoad(-200, x=200) n1 = Node.Crimp(x=0) n2 = Node.Crimp(x=1000) bar1 = BeamElement([n1, n2], section, material) beam = Beam(loads=[f1], beam_elements=[bar1]) support_reaction = beam.nodal_efforts assert support_reaction[0] == approx(179.2, abs=0.1) assert support_reaction[1] == approx(25600, abs=10) assert support_reaction[-2] == approx(20.8, abs=0.1) assert support_reaction[-1] == approx(-6400, abs=10) assert beam.getInternalShearStrength(150) == approx(179.2, abs=0.1) assert beam.getInternalShearStrength(300) == approx(-20.8, abs=0.1) assert beam.getInternalMomentumStrength(e) == approx(-25600, abs=10) assert beam.getInternalMomentumStrength(200) == approx(+10240, abs=10) assert beam.getInternalMomentumStrength(170) == approx(4860, abs=10) assert beam.getInternalMomentumStrength(210) == approx(10030, abs=10) assert beam.getInternalMomentumStrength(1000 - e) == approx(-6400, abs=10)