CN103559332A - 一种一阶螺旋杆频散特性的提取方法 - Google Patents
一种一阶螺旋杆频散特性的提取方法 Download PDFInfo
- Publication number
- CN103559332A CN103559332A CN201310473916.4A CN201310473916A CN103559332A CN 103559332 A CN103559332 A CN 103559332A CN 201310473916 A CN201310473916 A CN 201310473916A CN 103559332 A CN103559332 A CN 103559332A
- Authority
- CN
- China
- Prior art keywords
- screw rod
- wave number
- frequency
- order screw
- matrix
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明涉及一种一阶螺旋杆频散特性的提取方法,首先建立一阶螺旋杆模型,并在该模型母线上沿螺旋中心轴Z轴方向取等间距的M个点作为信号采集点,通过有限元计算各采集点的加速度瞬态信号,即形成一个行列矩阵X;将X中的每个行向量分别进行时域傅里叶变换;接下来对上步中得到的矩阵的每个列向量进行空间域的傅里叶变换,将导波信号的位置信息变换为波数信息,得到一阶螺旋杆结构波数—频率分布图;根据波数—频率分布图中能量的分布特点,绘制出一阶螺旋杆结构的频率—波数频散曲线,进而绘制出一阶螺旋杆结构的频率—相速度频散曲线。
Description
技术领域
本发明涉及一种基于二维傅里叶变换提取一阶螺旋杆结构频散特性的方法,属于超声导波无损检测领域。
背景技术
钢绞线被广泛应用于工程上的各个领域,如斜拉桥索、高层建筑电梯绳与索道等生命线工程。一阶螺旋杆结构是钢绞线中的基础结构,作为系统中重要的支撑构件,其机械性能状况的好坏直接影响到整体结构的安全和使用寿命。因此需要应用无损检测手段对钢绞线结构进行定期检测,以及时发现缺陷、评估其健康状态,以避免构件失效带来损失和灾害。超声导波技术作为一种新的无损检测方法,具有检测传播距离长、覆盖范围全面、效率高等特点,已被广泛应用于管、杆、板等结构的缺陷检测中。
频散特性是超声导波检测过程中确定检测方案的重要参考和依据,现有的波导频散特性提取方法普遍是通过理论计算来获得特定结构频散特性的数值描述。然而如钢绞线和圆柱螺旋压缩弹簧等含有一阶螺旋杆结构的复杂构件,由于其结构的复杂性无法通过理论计算直接得到其频散特性的数值描述,因此研究一阶螺旋杆结构的频散特性提取方法具有非常重要的现实意义。
发明内容
本发明的目的在于提供一种直接通过检测信号来提取一阶螺旋杆结构频散特性曲线的方法。一阶螺旋杆结构示意及坐标系如附图2所示,图中:(X,Y,Z)为固定参考坐标;Z轴为螺旋中心;(x,y,s)为沿螺旋线的曲线坐标系;θp为空间轨迹在XY平面投影圆弧所经过的角度,称作旋转角;Rh为螺旋线中心在XY平面投影的半径,称作螺旋半径;,Lp为螺旋线旋转一周沿Z向行进的距离,称作捻距。首先通过有限元仿真计算获得一阶螺旋杆模型沿螺旋中心轴Z轴方向上等间距的多组加速度瞬态信号,然后对得到的加速度瞬态信号进行二维傅里叶变换获得波数—频率分布图,最后分别对多模态中每个模态的极大值进行追踪和运算,提取出一阶螺旋杆结构的相速度频散曲线。因此,该方法可以用于一阶螺旋杆结构频散特性的提取,能够对钢绞线结构的频散特性研究提供有力支持。
本发明的具体实现步骤如下:
1)建立一阶螺旋杆模型,并在该模型母线上沿螺旋中心轴Z轴方向取等间距的M个点作为信号采集点,通过有限元计算各采集点的加速度瞬态信号,即形成一个行列矩阵X,其中行向量表示螺旋杆某一位置上采集到的不同时刻的加速度信息,列向量代表某一时刻所有信号采集点上的加速度信息;
2)将步骤1得到的矩阵X中的每个行向量分别进行时域傅里叶变换,得到经过第一次傅立叶变换后的矩阵,其中第i行第k列元素的变化公式如下:
其中:i=0,1,2…M-1,k=0,1,2…Ns-1,xi[n]代表矩阵X中第i个行向量,j代表虚部,M为信号采集点总个数,Ns为信号采样时刻;
3)对步骤2中得到的矩阵的每个列向量进行空间域的傅里叶变换,将导波信号的位置信息变换为波数信息,得到一阶螺旋杆结构波数—频率分布图,该分布图横坐标为频率f,纵坐标为波数k,其中第i行第k列元素的变化公式如下:
其中:Am[k]代表矩阵X经过步骤2变换后得到的第m行第k列元素,i=0,1,2…M-1,k=0,1,2…Ns-1,j代表虚部,M为信号采集点总个数;
4)根据步骤3中得到的波数—频率分布图中能量的分布特点,求取每个频率下能量的极大值点,从而获得各极大值点所对应的波数值,根据各个极大值点对应的频率和波数值,即可绘制出一阶螺旋杆结构的频率—波数频散曲线,其中频率f为横坐标,波数k为纵坐标;
5)相速度Vp和波数k、频率f有如下对应关系:
Vp=f/k
将步骤4中每个极大值点所对应的频率f与波数k代入上述公式,即可得到相速度Vp,以频率f为横坐标,相速度Vp为纵坐标绘制出一阶螺旋杆结构的频率—相速度频散曲线。
有益效果:
(1)操作简便,易于实施。(2)能够准确提取出一阶螺旋杆这一非规则结构波导的特定模态频散曲线,且与理论频散曲线具有很高的拟合度,适应性好。
附图说明:
图1、是本发明一种基于二维傅里叶变换一阶螺旋杆频散特性的提取方法流程图;
图2、一阶螺旋杆结构及坐标系示意图;
图3、是有限元仿真计算的一阶螺旋杆模型的多组瞬态加速度信号;
图4、是对瞬态信号进行二维傅里叶变换后得到的波数—频率分布图;
图5、1M频率下波数-频率曲线图;
图6、1M频率下相速度频散曲线图;
具体实施方式:
下面结合附图说明,说明本发明的具体实施方式。
1)在有限元仿真软件中建立一阶螺旋杆模型,设定材料参数(密度、杨氏模量和泊松比)。选择上升沿0.1μs、下降沿0.2μs的三角波脉冲信号,作为激励信号加载到螺旋杆一端的中心位置,并在螺旋杆母线Z方向上距离加载面50mm处每隔1mm设定400个接收节点,用于采集瞬态加速度信号。设定采样频率、采样点数以及采样时间。通过仿真计算,得到该一阶螺旋杆模型在径向x方向上瞬态加速度信号x(t),如附图3所示,400个接收节点上得到的加速度信号构成一个行列矩阵X;
2)将步骤1得到的矩阵中的每个行向量分别进行时域傅里叶变换,得到经过第一次傅立叶变换后的矩阵,其中第i行第k列元素的变化公式如下:
其中:i=0,1,2…M-1,k=0,1,2…Ns-1,xi[n]代表矩阵X中第i个行向量,即代表第i个采样点上的加速度数据,j代表虚部,信号采集点总个数M=400,采样时刻点数Ns=10000。
3)对步骤2中得到的矩阵中每个列向量进行空间傅里叶变换,可将导波信号的位置信息变换为波数信息,得到一阶螺旋杆结构波数—频率分布图,如图4所示。其中第i行第k列元素的变化公式如下:
其中:Am[k]代表矩阵X经过步骤2变换后得到的第m行第k列元素,即代表沿螺旋杆中心轴线Z轴方向的时域傅里叶变换的频谱值,i=0,1,2…M-1,k=0,1,2…Ns-1,j代表虚部,信号采集点总个数M=400。
4)根据步骤3中得到的波数—频率分布图中能量的分布特点,图中出现能量较大的部分(即图中颜色较亮位置)即为所求频散曲线。求取每个频率下能量的极大值点,从20kHz频率至930kHz频率,逐1kHz进行计算每个频率下能量的极大值,将每个频率下的波数-能量数据进行求导,导数为零的点即为该频率下的极大值点,从而获得各极大值点所对应的波数值。根据各个极大值点对应的频率和波数值,即可绘制出一阶螺旋杆结构的频率—波数频散曲线,其中频率f为横坐标,波数k为纵坐标。如附图5所示该频率段内的三条频率—波数频散曲线,分别为F(1,1),F(1,2)和L(2,1)模态的频率—波数频散曲线;
5)相速度Vp和波数k、频率f有如下对应关系:
Vp=f/k
将步骤4中得到的每一个极大值对应的频率f与波数k代入上述公式,即可得到对应导波模态连续的相速度Vp,以频率f为横坐标,相速度Vp为纵坐标绘制出一阶螺旋杆结构的频率—相速度频散曲线,如附图6所示该频率段内的三条频率—相速度曲线,分别为F(1,1),F(1,2)和L(2,1)模态的频率—相速度频散曲线;
最后应说明的是:以上实施例仅用以说明本发明而并非限制本发明所描述的技术方案;因此,尽管本说明书参照上述的各个实施例对本发明已进行了详细的说明,但是,本领域的普通技术人员应当理解,仍然可以对本发明进行修改或等同替换;而一切不脱离发明的精神和范围的技术方案及其改进,其均应涵盖在本发明的权利要求范围当中。
Claims (1)
1.一种一阶螺旋杆频散特性的提取方法,其特征在于包含以下步骤:
1)建立一阶螺旋杆模型,并在该模型母线上沿螺旋中心轴Z轴方向取等间距的M个点作为信号采集点,通过有限元计算各采集点的加速度瞬态信号,即形成一个行列矩阵X,其中行向量表示螺旋杆某一位置上采集到的不同时刻的加速度信息,列向量代表某一时刻所有信号采集点上的加速度信息;
2)将步骤1得到的矩阵X中的每个行向量分别进行时域傅里叶变换,得到经过第一次傅立叶变换后的矩阵,其中第i行第k列元素的变化公式如下:
其中:i=0,1,2…M-1,k=0,1,2…Ns-1,xi[n]代表矩阵X中第i个行向量,j代表虚部,M为信号采集点总个数,Ns为信号采样时刻;
3)对步骤2中得到的矩阵的每个列向量进行空间域的傅里叶变换,将导波信号的位置信息变换为波数信息,得到一阶螺旋杆结构波数—频率分布图,该分布图横坐标为频率f,纵坐标为波数k,其中第i行第k列元素的变化公式如下:
其中:Am[k]代表矩阵X经过步骤2变换后得到的第m行第k列元素,i=0,1,2…M-1,k=0,1,2…Ns-1,j代表虚部,M为信号采集点总个数;
4)根据步骤3中得到的波数—频率分布图中能量的分布特点,求取每个频率下能量的极大值点,从而获得各极大值点所对应的波数值,根据各个极大值点对应的频率和波数值,即可绘制出一阶螺旋杆结构的频率—波数频散曲线,其中频率f为横坐标,波数k为纵坐标;
5)相速度Vp和波数k、频率f有如下对应关系:Vp=f/k
将步骤4中每个极大值点所对应的频率f与波数k代入上述公式,即可得到相速度Vp,以频率f为横坐标,相速度Vp为纵坐标绘制出一阶螺旋杆结构的频率—相速度频散曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310473916.4A CN103559332B (zh) | 2013-10-11 | 2013-10-11 | 一种一阶螺旋杆频散特性的提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310473916.4A CN103559332B (zh) | 2013-10-11 | 2013-10-11 | 一种一阶螺旋杆频散特性的提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103559332A true CN103559332A (zh) | 2014-02-05 |
CN103559332B CN103559332B (zh) | 2016-04-06 |
Family
ID=50013578
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310473916.4A Active CN103559332B (zh) | 2013-10-11 | 2013-10-11 | 一种一阶螺旋杆频散特性的提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103559332B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107368644A (zh) * | 2017-07-14 | 2017-11-21 | 重庆交通大学 | 基于导波模态分叉角的钢绞线预应力检测方法 |
CN112464524A (zh) * | 2020-11-07 | 2021-03-09 | 西南交通大学 | 一种道岔变截面钢轨导波传播特性确定方法 |
CN112949121A (zh) * | 2021-02-05 | 2021-06-11 | 南京工业大学 | 一种导波传播特性的求解方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110280102A1 (en) * | 2010-05-17 | 2011-11-17 | Smith International, Inc. | Methods for making acoustic anisotropy logging while drilling measurements |
CN103336054A (zh) * | 2013-06-03 | 2013-10-02 | 北京工业大学 | 基于超声Lamb波的对接焊缝无损检测方法 |
-
2013
- 2013-10-11 CN CN201310473916.4A patent/CN103559332B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110280102A1 (en) * | 2010-05-17 | 2011-11-17 | Smith International, Inc. | Methods for making acoustic anisotropy logging while drilling measurements |
CN103336054A (zh) * | 2013-06-03 | 2013-10-02 | 北京工业大学 | 基于超声Lamb波的对接焊缝无损检测方法 |
Non-Patent Citations (3)
Title |
---|
FABIEN TREYSSEDE等: "A semi-analytical finite element method for elastic guided waves propagating in helical structures", 《ACOUSTICS 08 PARIS》, 4 July 2008 (2008-07-04) * |
吴斌等: "波导结构频散分析的特征频率法及在板条结构中的应用", 《计算力学学报》, vol. 30, no. 4, 15 August 2013 (2013-08-15) * |
周方俊等: "螺旋杆结构中弹性导波频散特性研究", 《华中科技大学学报(自然科学版)》, vol. 41, no. 1, 23 January 2013 (2013-01-23) * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107368644A (zh) * | 2017-07-14 | 2017-11-21 | 重庆交通大学 | 基于导波模态分叉角的钢绞线预应力检测方法 |
CN107368644B (zh) * | 2017-07-14 | 2019-10-25 | 重庆交通大学 | 基于导波模态分叉角的钢绞线预应力检测方法 |
CN112464524A (zh) * | 2020-11-07 | 2021-03-09 | 西南交通大学 | 一种道岔变截面钢轨导波传播特性确定方法 |
CN112949121A (zh) * | 2021-02-05 | 2021-06-11 | 南京工业大学 | 一种导波传播特性的求解方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN103559332B (zh) | 2016-04-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101561379B (zh) | 一种用于结构损伤检测的敲击扫描方法 | |
CN101762347B (zh) | 一种用半波法测量多跨钢拉索索力的方法 | |
CN112985671B (zh) | 基于全索系索力测试和误差自适应分析的损伤判定方法 | |
CN109839440B (zh) | 一种基于静置车辆测试的桥梁损伤定位方法 | |
CN202101909U (zh) | 锚杆拉拔力无损动力检测装置 | |
WO2014101832A1 (zh) | 敲击扫描式桥梁损伤检测系统 | |
CN106990169B (zh) | 基于前向散射波和c均值聚类算法的板类缺陷定位方法 | |
CN111256993A (zh) | 一种风电机组主轴承故障类型诊断方法及系统 | |
CN108519444A (zh) | 一种接触线缺陷位置的精确测量方法 | |
CN103559332B (zh) | 一种一阶螺旋杆频散特性的提取方法 | |
CN104314801B (zh) | 基于非消声环境的喷水推进泵噪声测试系统及其测试方法 | |
CN103592365A (zh) | 一种转子裂纹快速检测方法 | |
CN104406761A (zh) | 霍尔推力器低频振荡时间尺度内羽流发散角测量方法 | |
CN106078508B (zh) | 一种砂轮磨损在线监测方法 | |
Lei et al. | An automatic extraction algorithm for measurement of installed rock bolt length based on stress wave reflection | |
CN102681009B (zh) | 多通道瞬态瑞雷波测试系统的测试方法 | |
CN104101542A (zh) | 一种梁结构损伤监测装置及其监测方法 | |
CN101701882B (zh) | 塔式结构刚度的快速识别方法 | |
CN204359518U (zh) | 旋转电机力矩波动测试装置 | |
CN112985672B (zh) | 基于非接触空间振动测试的预应力索系索力分析方法 | |
CN103308386B (zh) | 一种测试土动本构关系与动强度参数的方法 | |
CN205313975U (zh) | 一种基于四翼无人机的桥梁检测系统 | |
CN103630604B (zh) | 离心式压缩机半开式叶轮裂纹故障识别方法 | |
CN204662519U (zh) | 一种基桩完整性的检测装置及其重锤 | |
CN108414624B (zh) | 基于全波形反演法的起重机异形截面结构起重伸缩臂的检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |