Test problem 12#
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, 4, 0, 10**4, 10**4, 10**4)
s.add_member(4, 0, 8, 3, E, I, A)
s.add_member(8, 3, 11, -1, E, I, A)
s.apply_load(0, 0, 15, 0, -1)
s.apply_load(8, 3, 16, 270, -1)
s.apply_load(0, 0, 6, 270, 0, 4, 0)
s.apply_load(4, 0, 6, 270, 0, 8, 3)
s.apply_load(8, 3, 6, 270, 0, 11, -1)
s.apply_support(11,-1,"fixed")
s.draw(show_load_values=True);
s.solve_for_reaction_loads()
{R_v (x=11,y=-1): -100, R_h (x=11,y=-1): -15, T (x=11,y=-1): -444}
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=11.00,y=-1.00) (unwrapped x=14.00) = -100
R_h (x=11.00,y=-1.00) (unwrapped x=14.00) = -15
T (x=11.00,y=-1.00) (unwrapped x=14.00) = -444
Points of Interest - Bending Moment:
bending_moment at (x=0.00,y=0.00) (unwrapped x=0.00) = 0
bending_moment at (x=2.00,y=0.00) (unwrapped x=2.00) = -12
bending_moment at (x=4.00,y=0.00) (unwrapped x=4.00) = -48
bending_moment at (x=4.00,y=0.00) (unwrapped x=4.00) = -48
bending_moment at (x=6.00,y=1.50) (unwrapped x=6.50) = -267/2
bending_moment at (x=8.00,y=3.00) (unwrapped x=9.00) = -249
bending_moment at (x=8.00,y=3.00) (unwrapped x=9.00) = -249
bending_moment at (x=9.50,y=1.00) (unwrapped x=11.50) = -1341/4
bending_moment at (x=11.00,y=-1.00) (unwrapped x=14.00-) = -443.999952000002
Points of Interest - Shear Force:
shear_force at (x=0.00,y=0.00) (unwrapped x=0.00+) = 0
shear_force at (x=4.00,y=0.00) (unwrapped x=4.00-) = -23.9999940000000
shear_force at (x=4.00,y=0.00) (unwrapped x=4.00+) = -141/5
shear_force at (x=8.00,y=3.00) (unwrapped x=9.00-) = -52.1999952000000
shear_force at (x=8.00,y=3.00) (unwrapped x=9.00+) = -30
shear_force at (x=11.00,y=-1.00) (unwrapped x=14.00-) = -47.9999964000000
Points of Interest - Axial Force:
axial_force at (x=0.00,y=0.00) (unwrapped x=0.00+) = -15
axial_force at (x=4.00,y=0.00) (unwrapped x=4.00-) = -15
axial_force at (x=4.00,y=0.00) (unwrapped x=4.00+) = 12/5
axial_force at (x=8.00,y=3.00) (unwrapped x=9.00-) = 20.3999964000000
axial_force at (x=8.00,y=3.00) (unwrapped x=9.00+) = -65
axial_force at (x=11.00,y=-1.00) (unwrapped x=14.00-) = -88.9999952000000