平板振动拍现象实例分析

摘要

笔者在技术邻最初的技术帖子中，《包络谱识别冲击振动在Ansys软件中如何仿真（一）》，使用APDL脚本语言，对一个结构进行连续冲击激励，然后提取结构上某点的振动响应。

提取振动：

/POST26

NSOL,2,7,ACC,Z !提取变量

*DIM,ACCE,ARRAY,NI!定义参数

VGET,ACCE(1),2 !将变量写入参数

*CFOPEN,ACCE,txt

*VWRITE,ACCE(1)!将参数写入文件

(F12.9)

*CFCLOS

振动时域：

t=np.arange(0,4,1/200)

  with open('ACCE.txt','r') as f:

      accy=f.read()

  accy=accy.split()

  accy=[float(x) for x in accy]

  plt.plot(t,accy)

  plt.show()

没有明显拍现象

振动频谱：

xf=np.fft.fft(accy)

  af=np.abs(xf)/len(xf)

  afn=np.short(af.size/2)

  f=np.arange(0,200/2,200/len(accy))

  plt.plot(f,af[:afn])

  plt.show()

滤波，保留红圈内频率：

N,Wn=sig.buttord([28,36],[25,40],3,30,fs=200)

  b,a=sig.butter(N,Wn,btype='bandpass',fs=200)

  acc2=sig.lfilter(b,a,accy)

查看滤波后时域特点：

总结：

本例表明，只要频谱上存在频率相近的幅值（固有频率，调制频率），则可以复现拍现象；

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删