BeamDeflection剧情

Question

我遇到了脚本没有显示剧情的问题。 绘图必须显示光束的偏转作为整个光束的x坐标的函数。我不知道如果我能做出声明： “×[I]> A [V]” 如果不提​​供X ...

import numpy as np #Imports NumPy 
import matplotlib.pyplot as plt 

def beamPlot(beamLength, loadPositions, loadForces, beamSupport): 

    l=beamLength  #Scalar 
    a=loadPositions #Vector 
    W=loadForces  #Vector 
    x=np.array(range(0,l)) 

    E=200*10**9  #Constant [N/m^2] 
    I=0.001   #Constant [m^4] 


#Makes an empty vector with the same size as x 
    y=np.empty_like(x) 

    for i in range(np.size(x)): #Continues as long as the vector x 
     for v in range(np.size(a)): 

      if a[v]==[ ] and W[v]==[ ]: 
       return np.zeros(np.size(x)) 

      elif beamSupport=="both" and x[i]<a[v]: 
       y[i]=np.sum(((W[v]*(l-a[v])*x[i])/(6*E*I*l))*(l**2-x[i]**2-(l-a[v])**2)) 
      elif beamSupport=="both" and x[i]>=a[v]: 
       y[i]=np.sum(W[v]*a[v]*(l-x[i])/(6*E*I*l)*(l**2-(l-x[i])**2-a[v]**2))       
      elif beamSupport=="cantilever" and x[i]<a[v]: 
       y[i]=np.sum((W[v]*x[i]**2)/(6*E*I)*(3*a[v]-x[i])) 
      elif beamSupport=="cantilever" and x[i]>=a[v]: 
       y[i]=np.sum((W[v]*a[v]**2)/(6*E*I)*(3*x[i]-a[v])) 

     deflection=y 


    plt.ylim([0,10000]) 
    plt.xlim([0,l]) 
    plt.title("Beam deflection") 
    plt.plot(x, deflection) 
plt.show() 

Answer 1

你的阵列x与来自range(0,l)整数列表创建，这意味着阵列中的元素的类型为int。您使用np.epty_like()创建了y阵列，这意味着它也具有int类型的元素。除非您使用巨大的负载值，否则当您将计算创建的浮点值转换为int时，四舍五入为0，因此该图在y = 0时为扁平线。

y=np.empty_like(x, dtype=float) 

你也应该删除plt.ylim(0,10000)，而是让matplotlib自动缩放你的y轴，因为：

您可以通过指定y应该包含浮点值时，它是通过将dtype=float创建解决这个问题对于任何合理的负载值（考虑到您的刚度），位移可能不会很大。