Test problem 10

Test problem 10#

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, 7, 1, E, I, A)
s.apply_load(3, 4, 60, 270, -1)
s.apply_load(3, 4, 18, 0, 0, 7, 1)
s.apply_support(0, 0,"fixed")
s.apply_support(7, 1,"roller")

s.draw(show_load_values=True);

../../_images/d2c265f522a30d701e154a3228db3ea7c5d0d72ed8f2aa40417724beb904b49a.png

s.solve_for_reaction_loads()

{R_v (x=0,y=0): -44421/1120,
 R_h (x=0,y=0): -90,
 T (x=0,y=0): 42021/160,
 R_v (x=7,y=1): -22779/1120}

s.plot_axial_force_on_structure();

../../_images/25b15364fd585c9f6ed9bc9a0e1b18e9db924f12c8c0e10748d218d7fb80a819.png

s.plot_shear_force_on_structure();

../../_images/0072c25387ad7ca7bfebde40040065da01d7a9546fb78dcfa9413ccda53fedd4.png

s.plot_bending_moment_on_structure();

../../_images/0af9c6e7ec67b87770371f0b30f6b1fd8c1c753f8e64cfc8d5d7846452981534.png

s.plot_deformation_on_structure();

../../_images/22dedff1d40e6c893451686a0e7fd3207a75af49a8086c8418cdead33b517774.png

s.summary();

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

Reaction Loads:
R_v   (x=0.00,y=0.00)  (unwrapped x=0.00)          = -44421/1120
R_h   (x=0.00,y=0.00)  (unwrapped x=0.00)          = -90
T     (x=0.00,y=0.00)  (unwrapped x=0.00)          = 42021/160
R_v   (x=7.00,y=1.00)  (unwrapped x=10.00)         = -22779/1120

Points of Interest - Bending Moment:
bending_moment at (x=0.00,y=0.00)  (unwrapped x=0.00) = -42021/160
bending_moment at (x=1.50,y=2.00)  (unwrapped x=2.50) = -51831/2240
bending_moment at (x=3.00,y=4.00)  (unwrapped x=5.00) = 60579/280
bending_moment at (x=3.00,y=4.00)  (unwrapped x=5.00) = 60579/280
bending_moment at (x=5.00,y=2.50)  (unwrapped x=7.50) = 41679/560
bending_moment at (x=7.00,y=1.00)  (unwrapped x=10.00-) = 0.0000162707196977863

Points of Interest - Shear Force:
shear_force at (x=0.00,y=0.00)  (unwrapped x=0.00+) = -oo
shear_force at (x=3.00,y=4.00)  (unwrapped x=5.00-) = 536463/5600
shear_force at (x=3.00,y=4.00)  (unwrapped x=5.00+) = -98379/1400
shear_force at (x=7.00,y=1.00)  (unwrapped x=10.00-) = -16.2707250857143

Points of Interest - Axial Force:
axial_force at (x=0.00,y=0.00)  (unwrapped x=0.00+) = 31179/1400
axial_force at (x=3.00,y=4.00)  (unwrapped x=5.00-) = 31179/1400
axial_force at (x=3.00,y=4.00)  (unwrapped x=5.00+) = 334863/5600
axial_force at (x=7.00,y=1.00)  (unwrapped x=10.00-) = -12.2030213142857