def plotBMD(my_beam): x=np.arange(-0.00099,my_beam.BEAM_LENGTH+0.01,0.01) supp_pos=[] supp_base=[] for i in my_beam.VERTICAL_REACTIONS: supp_pos.append(i) supp_base.append(0) moments=[] for i in x: moments.append(moment.momentX(my_beam,i)) plt.xlim(-my_beam.BEAM_LENGTH/2,my_beam.BEAM_LENGTH*3/2) plt.ylim(-max(moments)*3/2,max(moments)*3/2) plt.plot(supp_pos,supp_base,'^') #supports plt.plot(x,x*0,'b--') #beam plt.plot(x,moments,'g') #bending moment diagram plt.show()
def plotSFD(my_beam): x=np.arange(-0.00099,my_beam.BEAM_LENGTH+0.01,0.01) forces=[] for i in x: forces.append(forceAtX(my_beam,i)) moments=[] for i in x: moments.append(moment.momentX(my_beam,i)) supp_pos=[] supp_base=[] for i in my_beam.VERTICAL_REACTIONS: supp_pos.append(i) supp_base.append(0) plt.xlim(-my_beam.BEAM_LENGTH/2,my_beam.BEAM_LENGTH*3/2) plt.plot(supp_pos,supp_base,'^') #supports plt.plot(x,x*0,'b--') #beam plt.plot(x,forces,'r') #shear force diagram plt.show()