摘要

筆者在技術鄰最初的技術帖子中，《包絡譜識別沖擊振動在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)

查看濾波后時域特點：

總結：

本例表明，只要頻譜上存在頻率相近的幅值（固有頻率，調制頻率），則可以復現拍現象；