Test problem 2

Test problem 2#

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, 270, -1)
s.apply_load(2, 0, 16, 0, -1)
s.apply_load(0, 0, 6, 0, 0, 4, 0)
s.apply_load(4, 0, 6, 0, 0, 6, 1.5)
s.apply_support(11,-1,"fixed")

s.draw(show_load_values=True);

../../_images/b0f1d8d461b724b2527af1121f20cc5190851a95cd54f022f9ea6038941a032d.png

s.solve_for_reaction_loads()

{R_v (x=11,y=-1): -15, R_h (x=11,y=-1): -55, T (x=11,y=-1): -395/4}

s.plot_axial_force_on_structure();

../../_images/23dcff44f9819a07bb8ab3f359f5a328fda5c7c5caf4c880e09ffdd0864c9010.png

s.plot_shear_force_on_structure();

../../_images/b671a3c6d1cead120d59205d746d43d9ea0064d7525b5a121c4e020d52fd79a7.png

s.plot_bending_moment_on_structure();

../../_images/449635b251483c7076b72eac71e559e7dd4ad543effda5bea5b9f3c057a8a095.png

s.plot_deformation_on_structure();

../../_images/678f1eb7fc46e06579d562c0c684a6987056cf6f31320665b9f97e52fd499eb5.png

s.summary();

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

Reaction Loads:
R_v   (x=11.00,y=-1.00)  (unwrapped x=14.00)       = -15
R_h   (x=11.00,y=-1.00)  (unwrapped x=14.00)       = -55
T     (x=11.00,y=-1.00)  (unwrapped x=14.00)       = -395/4

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) = -30
bending_moment at (x=4.00,y=0.00)  (unwrapped x=4.00) = -60
bending_moment at (x=4.00,y=0.00)  (unwrapped x=4.00) = -60
bending_moment at (x=6.00,y=1.50)  (unwrapped x=6.50) = -645/4
bending_moment at (x=8.00,y=3.00)  (unwrapped x=9.00) = -1095/4
bending_moment at (x=8.00,y=3.00)  (unwrapped x=9.00) = -1095/4
bending_moment at (x=9.50,y=1.00)  (unwrapped x=11.50) = -745/4
bending_moment at (x=11.00,y=-1.00)  (unwrapped x=14.00-) = -98.7500350000000

Points of Interest - Shear Force:
shear_force at (x=0.00,y=0.00)  (unwrapped x=0.00+) = -15
shear_force at (x=2.00,y=0.00)  (unwrapped x=2.00-) = -15
shear_force at (x=2.00,y=0.00)  (unwrapped x=2.00+) = -15
shear_force at (x=4.00,y=0.00)  (unwrapped x=4.00-) = -15
shear_force at (x=4.00,y=0.00)  (unwrapped x=4.00+) = -36
shear_force at (x=8.00,y=3.00)  (unwrapped x=9.00-) = -45.0000000000000
shear_force at (x=8.00,y=3.00)  (unwrapped x=9.00+) = 35
shear_force at (x=11.00,y=-1.00)  (unwrapped x=14.00-) = 35.0000000000000

Points of Interest - Axial Force:
axial_force at (x=0.00,y=0.00)  (unwrapped x=0.00+) = 0
axial_force at (x=2.00,y=0.00)  (unwrapped x=2.00-) = -11.9999940000000
axial_force at (x=2.00,y=0.00)  (unwrapped x=2.00+) = -28
axial_force at (x=4.00,y=0.00)  (unwrapped x=4.00-) = -39.9999940000000
axial_force at (x=4.00,y=0.00)  (unwrapped x=4.00+) = -23
axial_force at (x=8.00,y=3.00)  (unwrapped x=9.00-) = -35.0000000000000
axial_force at (x=8.00,y=3.00)  (unwrapped x=9.00+) = -45
axial_force at (x=11.00,y=-1.00)  (unwrapped x=14.00-) = -45.0000000000000