Test problem 7

Test problem 7#

import micropip
await micropip.install('../../packages/sympy-1.15.0.dev0-py3-none-any.whl')
#takes a while

%pip install git+https://github.com/saiudayagiri/sympy.git@borek-27130

from sympy.physics.continuum_mechanics.structure2d import Structure2d
from sympy import symbols
s = Structure2d()
E, I, A = symbols('E I A')
s.add_member(0, 0, 3, 4, 10**4, 10**4, 10**4)
s.add_member(3, 4, 6, 0, E, I, A)
s.apply_load(3, 4, 60, 270, -1)
s.apply_support(0, 0,"pin")
s.apply_support(6, 0,"pin")

s.draw(show_load_values=True);

../../_images/0be838a44148607b92622d9f99c3f3a1b6e7de18a4f391c4731e7eb020666f06.png

s.solve_for_reaction_loads()

{R_v (x=0,y=0): -30,
 R_h (x=0,y=0): 7920/427,
 R_v (x=6,y=0): -30,
 R_h (x=6,y=0): -7920/427}

s.plot_axial_force_on_structure();

../../_images/e47407e79a172f670a364c84f57080587d58e24f9773ddd9443e0aa13f9c23b0.png

s.plot_shear_force_on_structure();

../../_images/5d208643f72485684fcb6b172aaed7602fdc4b8b4210067189c6fff8a5fc10e3.png

s.plot_bending_moment_on_structure();

../../_images/bb332dfc3c1b138421a3b77d3e443c3ccd3e76723f1d66702faf73ebc4af4e41.png

s.plot_deformation_on_structure();

../../_images/670ec43d6e4a5813117c4681cf49501d5bf0d3ee1839a59960db51b05f45bd52.png

s.summary();

===================== Structure Summary =====================

Reaction Loads:
R_v   (x=0.00,y=0.00)  (unwrapped x=0.00)          = -30
R_h   (x=0.00,y=0.00)  (unwrapped x=0.00)          = 7920/427
R_v   (x=6.00,y=0.00)  (unwrapped x=10.00)         = -30
R_h   (x=6.00,y=0.00)  (unwrapped x=10.00)         = -7920/427

Points of Interest - Bending Moment:
bending_moment at (x=0.00,y=0.00)  (unwrapped x=0.00) = 0
bending_moment at (x=1.50,y=2.00)  (unwrapped x=2.50) = 3375/427
bending_moment at (x=3.00,y=4.00)  (unwrapped x=5.00) = 6750/427
bending_moment at (x=3.00,y=4.00)  (unwrapped x=5.00) = 6750/427
bending_moment at (x=4.50,y=2.00)  (unwrapped x=7.50) = 3375/427
bending_moment at (x=6.00,y=0.00)  (unwrapped x=10.00-) = 0.00000316159250246528

Points of Interest - Shear Force:
shear_force at (x=0.00,y=0.00)  (unwrapped x=0.00+) = 1350/427
shear_force at (x=3.00,y=4.00)  (unwrapped x=5.00-) = 1350/427
shear_force at (x=3.00,y=4.00)  (unwrapped x=5.00+) = -1350/427
shear_force at (x=6.00,y=0.00)  (unwrapped x=10.00-) = -1350/427

Points of Interest - Axial Force:
axial_force at (x=0.00,y=0.00)  (unwrapped x=0.00+) = -15000/427
axial_force at (x=3.00,y=4.00)  (unwrapped x=5.00-) = -15000/427
axial_force at (x=3.00,y=4.00)  (unwrapped x=5.00+) = -15000/427
axial_force at (x=6.00,y=0.00)  (unwrapped x=10.00-) = -15000/427