Test problem 6#
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, 0, -1)
s.apply_support(0, 0,"pin")
s.apply_support(6, 0,"pin")
s.draw(show_load_values=True);
s.solve_for_reaction_loads()
{R_v (x=0,y=0): 40, R_h (x=0,y=0): -30, R_v (x=6,y=0): -40, R_h (x=6,y=0): -30}
s.plot_axial_force_on_structure();
s.plot_shear_force_on_structure();
s.plot_bending_moment_on_structure();
s.plot_deformation_on_structure();
s.summary();
===================== Structure Summary =====================
Reaction Loads:
R_v (x=0.00,y=0.00) (unwrapped x=0.00) = 40
R_h (x=0.00,y=0.00) (unwrapped x=0.00) = -30
R_v (x=6.00,y=0.00) (unwrapped x=10.00) = -40
R_h (x=6.00,y=0.00) (unwrapped x=10.00) = -30
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) = 0
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) = 0
bending_moment at (x=6.00,y=0.00) (unwrapped x=10.00-) = 0
Points of Interest - Shear Force:
shear_force at (x=0.00,y=0.00) (unwrapped x=0.00+) = 0
shear_force at (x=3.00,y=4.00) (unwrapped x=5.00-) = 0
shear_force at (x=3.00,y=4.00) (unwrapped x=5.00+) = 0
shear_force at (x=6.00,y=0.00) (unwrapped x=10.00-) = 0
Points of Interest - Axial Force:
axial_force at (x=0.00,y=0.00) (unwrapped x=0.00+) = 50
axial_force at (x=3.00,y=4.00) (unwrapped x=5.00-) = 50
axial_force at (x=3.00,y=4.00) (unwrapped x=5.00+) = -50
axial_force at (x=6.00,y=0.00) (unwrapped x=10.00-) = -50