Test problem 13

Test problem 13#

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(8, 3, 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, 8, 3)

s.apply_support(11,-1,"fixed")

s.draw(show_load_values=True);

../../_images/9c2f6ea39e82972ef8d80d483c6e0eeafc0520410761d3f56f425b218b25c366.png

s.solve_for_reaction_loads()

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

s.plot_axial_force();

../../_images/a84febe61ce279c8628db9b1e9faabd1656383eff471feaf585ea4cad8821124.png

s.plot_shear_force_on_structure();

../../_images/51c6a7f97f3370db2e29a7e5fd7b6444dcb2c53ec2dd864896e24d55a469cb95.png

s.plot_bending_moment_on_structure();

../../_images/ecb8ddc19964aac12dca5ea41795b4a571c9b653e4157784ad49c9501706bc99.png

s.plot_deformation_on_structure();

../../_images/06b2a4bdc2bf2f19926501c55718d484584c9285699f23db55d799c63ce86446.png

s.summary();

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

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

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) = -175
bending_moment at (x=8.00,y=3.00)  (unwrapped x=9.00) = -300
bending_moment at (x=8.00,y=3.00)  (unwrapped x=9.00) = -300
bending_moment at (x=9.50,y=1.00)  (unwrapped x=11.50) = -375
bending_moment at (x=11.00,y=-1.00)  (unwrapped x=14.00-) = -449.999970000000

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+) = -32
shear_force at (x=8.00,y=3.00)  (unwrapped x=9.00-) = -55.9999952000000
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-) = -30.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-) = 0
axial_force at (x=2.00,y=0.00)  (unwrapped x=2.00+) = 0
axial_force at (x=4.00,y=0.00)  (unwrapped x=4.00-) = 0
axial_force at (x=4.00,y=0.00)  (unwrapped x=4.00+) = 24
axial_force at (x=8.00,y=3.00)  (unwrapped x=9.00-) = 41.9999964000000
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-) = -65.0000000000000