```python
%matplotlib widget
```


```python
from sympy import *
x, y, z = symbols('x y z')
init_printing(use_unicode=True, wrap_line=False)
```


```python
from sympy.physics.continuum_mechanics.beam import Beam
E, I = symbols('E, I')
b = Beam(9, E, I)
```


```python
b.apply_load(-12, 9, -1)
b.apply_load(50, 5, -2)
b.apply_load(-8, 0, 0, end=5)
```


```python
b.bc_deflection.append((0, 0))
b.bc_slope.append((0, 0))
```


```python
R, M = symbols('R, M')
b.apply_load(R, 0, -1)
b.apply_load(M, 0, -2)
b.load
```


```python
b.solve_for_reaction_loads(R, M)
b.reaction_loads
```


```python
b.load
```


```python
b.shear_force()
```


```python
b.bending_moment()
```


```python
b.slope()
```


```python
b.deflection()
```


```python
b.plot_shear_force()
```


    Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …


    C:\Users\spart\miniconda3\envs\bmcs_env\lib\site-packages\sympy\plotting\plot.py:1486: RuntimeWarning: invalid value encountered in double_scalars
      cos_theta = dot_product / (vector_a_norm * vector_b_norm)





    <sympy.plotting.plot.Plot at 0x1fd4c670940>




```python
b.plot_bending_moment()
```


    Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …





    <sympy.plotting.plot.Plot at 0x1fd4c5be7c0>




```python
b.plot_slope(subs={E: 20E9, I: 3.25E-6})
```


    Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …





    <sympy.plotting.plot.Plot at 0x1fd4d991cd0>




```python
b.plot_deflection(subs={E: 20E9, I: 3.25E-6})
```


    Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …





    <sympy.plotting.plot.Plot at 0x1fd4de00ac0>




```python
b.plot_loading_results(subs={E: 20E9, I: 3.25E-6})
```


    Canvas(toolbar=Toolbar(toolitems=[('Home', 'Reset original view', 'home', 'home'), ('Back', 'Back to previous …





    <sympy.plotting.plot.PlotGrid at 0x1fd4e4f4af0>



e102009dff2bd254ba91387d49e06ba14969c122
ipynb
Jupyter Notebook
notebooks/deflection/external_libraries/sympy_beam/sympy_beam_example.ipynb
bmcs-group/bmcs_beam
b53967d0d0461657ec914a3256ec40f9dcff80d5
MIT
2021-05-07T11:10:27.000Z
2021-05-07T11:10:27.000Z
notebooks/deflection/external_libraries/sympy_beam/sympy_beam_example.ipynb
bmcs-group/bmcs_beam
b53967d0d0461657ec914a3256ec40f9dcff80d5
MIT
notebooks/deflection/external_libraries/sympy_beam/sympy_beam_example.ipynb
bmcs-group/bmcs_beam
b53967d0d0461657ec914a3256ec40f9dcff80d5
MIT
Qwen/Qwen-72B
1. YES
2. YES