发明名称 一种基于压电传感网络的工程结构冲击位置确定方法
摘要 本发明公开了一种基于压电传感网络的工程结构冲击位置确定方法,包括布置分布式传感器网络的步骤、记录受载后的结构信息的步骤、进行载荷危险判别的步骤、进行特征信息提取的步骤、进行传感子网络的选择的步骤、冲击位置概率成像的步骤,本发明的基于压电传感网络的工程结构冲击位置确定方法结合了信号处理与概率计算的方式进行冲击位置图像的计算,并合理的筛选了定位子网的形式,不仅计算时间、计算量少,而且冲击定位的精度高,稳定性好,尤其能够用于加筋壁板这类复杂的结构中,为航天结构健康监测技术提供了一个可行的方案。
申请公布号 CN103438794A 申请公布日期 2013.12.11
申请号 CN201310352408.0 申请日期 2013.08.13
申请人 南京航空航天大学 发明人 裘进浩;张超;王晓宇;吴郁程;孙维;单胜博;张立新
分类号 G01B7/00(2006.01)I;G01M7/08(2006.01)I 主分类号 G01B7/00(2006.01)I
代理机构 南京经纬专利商标代理有限公司 32200 代理人 叶连生
主权项 一种基于压电传感网络的工程结构冲击位置确定方法,其特征在于,包括下列步骤:步骤一:布置分布式传感器网络:在结构中均匀地布置压电传感器构成传感网络,并保证传感器覆盖待测工程结构的全部区域;步骤二:通过步骤一布置的分布式传感器网络监测工程结构中不同位置处的信号,当结构受到冲击载荷后,记录一组受载后的结构信息。步骤三:进行载荷危险判别:选择步骤二所得的一组传感信号中幅值最大的信号进行分析,利用傅里叶变换求解信号中不同频率的能量,若对结构安全产生影响的频带中含有超过危险判别阈值的成分,则预判此次冲击是危险的,并进入冲击位置确定的过程,否则忽略此次冲击,进入下一次冲击信号的采集过程,重复进行步骤二;步骤四:进行特征信息提取:对冲击信号进行固定中心频率和带宽的小波变换,提取出带宽分量,如下: <mrow> <mi>CWT</mi> <mrow> <mo>(</mo> <mi>a</mi> <mo>,</mo> <mi>b</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mi>a</mi> </msqrt> </mfrac> <msubsup> <mo>&Integral;</mo> <mrow> <mo>-</mo> <mo>&infin;</mo> </mrow> <mrow> <mo>+</mo> <mo>&infin;</mo> </mrow> </msubsup> <mi>v</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <msup> <mi>&psi;</mi> <mo>*</mo> </msup> <mrow> <mo>(</mo> <mfrac> <mrow> <mi>t</mi> <mo>-</mo> <mi>b</mi> </mrow> <mi>a</mi> </mfrac> <mo>)</mo> </mrow> <mi>dt</mi> </mrow>其中,CWT(a,b)为小波系数,其中a为尺度因子,b为时间因子,v(t)为冲击信号,ψ(t)为母小波函数,具体表达式如下: <mrow> <mi>&psi;</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <msqrt> <mi>&pi;&gamma;</mi> </msqrt> </mfrac> <msup> <mi>e</mi> <mrow> <mi>j</mi> <msub> <mi>&omega;</mi> <mn>0</mn> </msub> <mi>t</mi> </mrow> </msup> <msup> <mi>e</mi> <mrow> <msup> <mrow> <mo>-</mo> <mi>t</mi> </mrow> <mn>2</mn> </msup> <mo>/</mo> <mi>&gamma;</mi> </mrow> </msup> </mrow>其中ω为中心频率,γ为频带宽度;步骤五:进行传感子网络的选择:对冲击信号完成特征提取后,选出小波系数值较大的传感器所处的子网络对冲击进行定位;步骤六:冲击位置概率成像:根据所选的传感器子网络计算冲击位置时间差概率函数PoTD(t),计算次数由子网络包含的传感器个数N决定,其计算的公式如下所示:PoTDi,j(t)=Ri,j(t)/Max(Ri,j(t))其中i,j为传感器的序号,R(t)为传感器小波系数的互相关函数;为了综合所有求的PoTD(t)函数并求得强度图表示冲击位置概率,强度图像素矩阵由如下公式确定: <mrow> <mi>I</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mrow> <mi>&Sigma;</mi> <msub> <mi>PoTD</mi> <mrow> <mi>i</mi> <mo>,</mo> <mi>j</mi> </mrow> </msub> <mrow> <mo>(</mo> <mi>t</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>)</mo> </mrow> </mrow> <msub> <mi>N</mi> <mi>p</mi> </msub> </mfrac> </mrow>其中,I(x,y)表示冲击位置概率图,(x,y)表示像素的位置,t(x,y)为像素与传感器之间的距离差,具体计算如下:t(x,y)=(di‑dj)/cgNp=N*(N‑1)/2di表示当前计算的像素值位置与第i号传感器位置的距离,cg为所选频率下弹性波在待测结构中传播的群速度。
地址 210016 江苏省南京市御道街29号