Test problem 5

Test problem 5#

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(0, 0, 60, 0, 0, 3, 4)
s.apply_load(3, 4, 60, 0, 0, 6, 0)
s.apply_support(0, 0,"pin")
s.apply_support(6, 0,"pin")

s.draw(show_load_values=True);

../../_images/7ffb814ceddb1dfbf31624fa054291d6ed10e1deab8b96dce5d988cc934d88a4.png

s.solve_for_reaction_loads()

{R_v (x=0,y=0): 200,
 R_h (x=0,y=0): -300,
 R_v (x=6,y=0): -200,
 R_h (x=6,y=0): -300}

s.plot_axial_force_on_structure();

../../_images/e1e03855ac17fc180a52cdabd7167ea01e33f8971d1b0554f2008920d1d2173a.png

s.plot_shear_force_on_structure();

../../_images/b5634f296873aeeb96cefd119012d1716fefd50a10d8ed1911cff92059a6559c.png

s.plot_bending_moment_on_structure();

../../_images/a5f3c1f03ba5ab18c19dfc0ff14851cba1ac8c4794aa3f8645965990d6632e1d.png

s.plot_deformation_on_structure();

../../_images/7146b11302fb362481363c523283e2fe83cf7c9ac8d7b16a1d9e38ce26b997c4.png

s.summary();

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

Reaction Loads:
R_v   (x=0.00,y=0.00)  (unwrapped x=0.00)          = 200
R_h   (x=0.00,y=0.00)  (unwrapped x=0.00)          = -300
R_v   (x=6.00,y=0.00)  (unwrapped x=10.00)         = -200
R_h   (x=6.00,y=0.00)  (unwrapped x=10.00)         = -300

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) = 150
bending_moment at (x=3.00,y=4.00)  (unwrapped x=5.00) = 0
bending_moment at (x=3.00,y=4.00)  (unwrapped x=5.00) = 0
bending_moment at (x=4.50,y=2.00)  (unwrapped x=7.50) = -150
bending_moment at (x=6.00,y=0.00)  (unwrapped x=10.00-) = -0.000119999976050167

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

Points of Interest - Axial Force:
axial_force at (x=0.00,y=0.00)  (unwrapped x=0.00+) = 340
axial_force at (x=3.00,y=4.00)  (unwrapped x=5.00-) = 160.000036000000
axial_force at (x=3.00,y=4.00)  (unwrapped x=5.00+) = -160
axial_force at (x=6.00,y=0.00)  (unwrapped x=10.00-) = -339.999964000000