Influence line support reaction A

Click –> Live Code on the top right corner of this screen to investigate some influence lines!

import micropip
await micropip.install("ipympl")

import numpy as np
import sympy as sym
import matplotlib.pylab as plt
from ipywidgets import widgets, interact
from matplotlib.patches import Arc,StepPatch

%matplotlib widget

EI, x = sym.symbols('EI, x',real=True)

Bv, Cv, Cm = sym.symbols('Bv, Cv, Cm',real=True)
a, b, L1, L2 = sym.symbols('a, b, L1, L2',positive=True,real=True)
C1, C2, C3, C4 = sym.symbols('C1, C2, C3, C4')

q = - Bv * sym.SingularityFunction(x,0,-1) - Cv * sym.SingularityFunction(x,L1,-1) + Cm * sym.SingularityFunction(x,L1,-2) + 1 * sym.SingularityFunction(x,a,-1)
V = -sym.integrate(q,x)+C1
M = sym.integrate(V,x)+C2
kappa = M / EI
phi = sym.integrate(kappa,x)+C3
w = sym.integrate(phi,x)+C4
eq1 = sym.Eq(V.subs(x,-1),0)
eq2 = sym.Eq(M.subs(x,-1),0)
eq3 = sym.Eq(w.subs(x,0),0)
eq4 = sym.Eq(w.subs(x,L1),0)
eq5 = sym.Eq(phi.subs(x,L1),0)
eq6 = sym.Eq(M.subs(x,L1+1),0)
eq7 = sym.Eq(V.subs(x,L1+1),0)
sol = sym.solve([eq1,eq2,eq3,eq4,eq5,eq6,eq7],[C1,C2,C3,C4,Bv, Cv, Cm])
Bv_sol = sol[Bv].rewrite(sym.Piecewise)
Bv_numpy = sym.lambdify([L1,a,EI],Bv_sol)
Cv_sol = sol[Cv].rewrite(sym.Piecewise)
Cv_numpy = sym.lambdify([L1,a,EI],Cv_sol)
Cm_sol = sol[Cm].rewrite(sym.Piecewise)
Cm_numpy = sym.lambdify([L1,a,EI],Cm_sol)
V_sol = V.subs(sol)
V_subs = V_sol.subs([(L1,10),(a,2.5),(EI,5000)])
V_numpy = sym.lambdify([L1,a,EI,x],V_sol.rewrite(sym.Piecewise))
M_sol = M.subs(sol)
M_subs = M_sol.subs([(L1,10),(a,2.5),(EI,5000)])
M_numpy = sym.lambdify([L1,a,EI,x],M_sol.rewrite(sym.Piecewise))

fig, axs = plt.subplots(2, 1, figsize=(7, 6))
x = np.linspace(0,10,100)
def func(b):
    axs[0].clear()  # Clear the existing plot
    axs[1].clear()
    
    axs[0].plot([0,10],[0,0],color='black',linewidth=2)
    axs[0].annotate(text='', xy=(0,0), xytext=(0,Bv_numpy(10,b,5000)), arrowprops=dict(arrowstyle='simple'))
    axs[0].annotate('%.2f kN' % Bv_numpy(10,b,5000),xy = [0.2,0.2])
    axs[0].annotate(text='', xy=(10,0), xytext=(10,Cv_numpy(10,b,5000)), arrowprops=dict(arrowstyle='simple'))
    axs[0].annotate(text='', xy=(10,0.5*Cm_numpy(10,b,5000)/-2), xytext=(10,-0.5*Cm_numpy(10,b,5000)/-2), arrowprops=dict(arrowstyle='simple',connectionstyle="arc3,rad=-.5"))
    axs[0].annotate(text='', xy=(b,0), xytext=(b,-1), arrowprops=dict(arrowstyle='simple'))
    axs[0].set_ylim([-1,1])
    axs[0].invert_yaxis()
    axs[0].axis('off')
    title0 = 'Free body diagram for unit load at $x_F = '+str(b)+'$'
    axs[0].set_title(title0)
    axs[1].plot([0,10],[0,0],color='black',linewidth=2)
    axs[1].grid()
    axs[1].axis('equal')
    axs[1].plot(x,Bv_numpy(10,x,5000),color='blue')
    axs[1].plot(b,Bv_numpy(10,b,5000),marker='o')
    axs[1].annotate('%.2f kN' % Bv_numpy(10,b,5000),xy = [b,Bv_numpy(10,b,5000)+0.4])
    axs[1].invert_yaxis()
    title1 = 'Influence line for support reaction at B'
    axs[1].set_title(title1)

    axs[1].spines['right'].set_color('none')
    axs[1].spines['top'].set_color('none')
    axs[1].spines['bottom'].set_position('zero')
    axs[1].spines['left'].set_position('zero')
    plt.draw()  

    
interact(func, b = widgets.FloatSlider(min=0, max=10, value=2, step=0.1, description="Location unit load x_F =  ... (m)",readout_format='.1f',style= {'description_width': '180px'},layout = {'width': '400px'}) );