SHARE
TWEET

Untitled

a guest Jan 12th, 2017 74 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. v_aluminium = 0.33*1e-6
  2.  
  3. def mobility_plate(E=Youngs_modulus, p=density, v=v_aluminium, h=thickness):
  4.     return 1./(8.*(h**2.) * (np.sqrt((E*p)/(12.*(1-v**2.)))))
  5.  
  6. Y_flr_section = mobility_plate()
  7. PWL_ref = 10*1e-12
  8. def w_in(vel_in=accel_to_vel(), Y=Y_flr_section, Yi=Y_Zmounts()[1], Ys=Y_Zengine()[1], Yb=Y_Zbeam()[1], dB=False):
  9.     """input power from the engine to the vehicle, velocity calculated from the
  10.    accleration spectrum, and going through vib isolation"""
  11.     top = [np.real(Y)*((vel_in1)**2) for vel_in1 in vel_in]
  12.     bottom = [pow(np.abs(Y+Yi1+Ys1+Yb1), 2) for Yi1,Ys1,Yb1 in zip(Yi,Ys,Yb)]
  13.     w = [t/b for t,b in zip(top,bottom)]
  14.     if dB == False:
  15.         return w
  16.     if dB == True:
  17.         Lp = [10*np.log10(w1/PWL_ref) for w1 in w]
  18.         return Lp
RAW Paste Data
Pastebin PRO Summer Special!
Get 60% OFF on Pastebin PRO accounts!
Top