Test problem 1

Test problem 1#

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

s.draw(show_load_values=True);

../../_images/9837d1ffb08351fc1b7e33e38c41abc62f1e862b850f2deea008bcbda68b9a50.png

print(s.solve_for_reaction_loads())

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

s.plot_axial_force_on_structure();

../../_images/6e0f8e77707ce85030d360e03ba2fb7e0e2771975e1fe43c0b8c6a09204d7c77.png

s.plot_shear_force_on_structure();

../../_images/eb47423fe6c6b0439a747f420cad8072fa65b9967fb59a03df0690ead5324446.png

s.plot_bending_moment_on_structure();

../../_images/8ddc91cdaebcab10bc22efe0a41630b90b245ddf13214b239ca61018ef54a573.png

s.plot_deformation_on_structure();

../../_images/fc8bf29c37620eb4d9c483839c15a21c53dd4717ad3782cacc1dbf017728a7f2.png

s.summary();

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

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

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) = -80
bending_moment at (x=4.00,y=0.00)  (unwrapped x=4.00) = -80
bending_moment at (x=6.00,y=1.50)  (unwrapped x=6.50) = -395/2
bending_moment at (x=8.00,y=3.00)  (unwrapped x=9.00) = -330
bending_moment at (x=8.00,y=3.00)  (unwrapped x=9.00) = -330
bending_moment at (x=9.50,y=1.00)  (unwrapped x=11.50) = -765/2
bending_moment at (x=11.00,y=-1.00)  (unwrapped x=14.00-) = -434.999979000000

Points of Interest - Shear Force:
shear_force at (x=0.00,y=0.00)  (unwrapped x=0.00+) = 0
shear_force at (x=2.00,y=0.00)  (unwrapped x=2.00-) = -11.9999940000000
shear_force at (x=2.00,y=0.00)  (unwrapped x=2.00+) = -28
shear_force at (x=4.00,y=0.00)  (unwrapped x=4.00-) = -39.9999940000000
shear_force at (x=4.00,y=0.00)  (unwrapped x=4.00+) = -41
shear_force at (x=8.00,y=3.00)  (unwrapped x=9.00-) = -53.0000000000000
shear_force at (x=8.00,y=3.00)  (unwrapped x=9.00+) = -21
shear_force at (x=11.00,y=-1.00)  (unwrapped x=14.00-) = -21.0000000000000

Points of Interest - Axial Force:
axial_force at (x=0.00,y=0.00)  (unwrapped x=0.00+) = -15
axial_force at (x=2.00,y=0.00)  (unwrapped x=2.00-) = -15
axial_force at (x=2.00,y=0.00)  (unwrapped x=2.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
axial_force at (x=8.00,y=3.00)  (unwrapped x=9.00-) = 21.0000000000000
axial_force at (x=8.00,y=3.00)  (unwrapped x=9.00+) = -53
axial_force at (x=11.00,y=-1.00)  (unwrapped x=14.00-) = -53.0000000000000