CN101915733A - 一种基于频率变化的结构损伤程度评估方法 - Google Patents
一种基于频率变化的结构损伤程度评估方法 Download PDFInfo
- Publication number
- CN101915733A CN101915733A CN 201010244573 CN201010244573A CN101915733A CN 101915733 A CN101915733 A CN 101915733A CN 201010244573 CN201010244573 CN 201010244573 CN 201010244573 A CN201010244573 A CN 201010244573A CN 101915733 A CN101915733 A CN 101915733A
- Authority
- CN
- China
- Prior art keywords
- damage
- degree
- frequency
- injury
- finite element
- 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
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明提出一种基于频率变化的结构损伤程度评估方法,在已知损伤位置时可以利用结构实测频率的变化来评估损伤的严重程度。该方法对损伤单元设立损伤程度系数,基于基准有限元模型和损伤结构有限的低阶实测模态频率,利用结构特征值建立了方程组,通过迭代求解,可以准确求解该损伤程度系数。本发明利用损伤前后的频率信息的差异进行损伤程度评估,不需要结构的振型信息和高阶的模态信息。本发明仅利用结构损伤前后的有限的低阶频率信息,就可以进行结构的损伤程度评估。该方法适用于单损伤、多种损伤工况,能够准确评估结构的损伤程度,具有一定的实际应用价值。
Description
技术领域
本发明涉及一种基于频率变化的结构损伤程度评估新方法,在已知损伤位置时可以利用结构实测频率的变化来评估损伤的严重程度。
背景技术
海洋平台等大型工程结构在其服役期间损伤是不可避免的。为确保人员的生命安全、减少财产损失,唯一方法是尽早诊断出结构物的结构损伤,并对损伤的程度进行评估,以便为及时的修复加固提供科学依据和技术支持。对于如构件断裂的损伤,目测能够判定损伤。然而,对于材料老化造成的损伤、损伤在不可靠近的区域或者被油漆、锈以及海洋生物覆盖着是很难利用目测检测损伤的,所以利用目测进行无损检测是不可靠的。其它无损检测技术,如X光检测、超声检测、工业CT等方法,都属于结构局部损伤检测方法。这类技术要求损伤的大致区域必须已知的情况下才能使用,而且必须位于测试仪器可及的地方,需要专业的技术人员对检测结果分析,费时费力。目前基于结构动力特性变化的结构损伤检测方法正受到重视。众所周知,任何结构都可以看作由刚度、质量、阻尼矩阵组成的动力学系统,结构一旦出现损伤,结构的特性参数随之发生改变,从而导致系统的模态参数(固有频率、阻尼和振型等)发生变化,所以结构的模态参数的改变可用于判定结构是否有损伤发生,并进而评估损伤的严重程度。利用损伤发生前后结构模态参数的变化来诊断结构损伤的方法,其优点是可将环境荷载作为结构的振动激励源,损伤检测的过程不影响结构的正常使用,能方便地完成结构损伤的在线检测和评估。
结构损伤检测主要包括四个层次:(1)判断结构是否发生损伤(损伤识别);(2)确定结构的损伤位置(损伤定位);(3)评估结构的损伤程度(损伤程度评估);(4)结构剩余寿命的预测。基于结构动力特性变化的损伤诊断方法经过多年发展,研究人员提出并发展了许多方法,但大都集中在一、二层次,结构损伤程度评估的研究相对较少。而且这些方法在实际运用于大型结构存在一定的局限性,主要表现在以下几个方面:
(1)一些方法进行损伤检测时需要完备或高阶的模态信息:实际测试得到的结构模态参数是不完备的,模态不完备包含量测自由度不完备和量测模态阶数不完备。前者是由于传感器数量的限制和某些自由度信息(如转动自由度)难以量测造成的。而后者是由于激振荷载的频率范围和识别技术等方面的原因造成的,无法得到高阶模态信息。
(2)一些方法需要质量归一化的阵型信息:对实际的大型结构进行振动测试时,人工激励是很困难的。因此常利用环境荷载(如风、波浪)进行激励来进行振动测试,此时激励信息是不知道的;在仅仅知道输出(振动响应)条件下,模态参数识别得到的振型无法质量归一化。
如何克服这些缺点,仅利用识别得到的有限的低阶模态信息进行损伤检测,是目前研究的重点和难点。
发明内容
本发明为解决现有结构损伤程度评估技术需要振型信息、需要高阶模态的的问题,提出一种基于频率变化的结构损伤程度评估方法,该方法不需要振型信息,仅利用结构损伤前后的有限的低阶频率信息即可进行损伤程度的准确评估。
为解决上述技术问题,本发明采用以下技术方案予以实现:
1、一种基于频率变化的结构损伤程度评估方法,其特征在于,包括如下步骤:
A、建立基准有限元模型;
C、损伤程度评估:
其中K为基准有限元模型的总体刚度矩阵;和αn分别为第n个损伤单元的单元刚度矩阵和损伤程度;Nd为损伤单元的个数;Φi和分别为基准有限元模型和损伤结构的第i阶振型;ωi和分别为基准有限元模型和损伤结构的第i阶频率;上标T代表矩阵或向量的转置;
第1步:迭代初始赋值。k=0,假设损伤程度为0,即α(0)=0,Φ*(0)=Φ,上标“0”代表迭代初值;
第2步:迭代开始,k=k+1。由α=(CTC)-1CTb计算α(k);
第4步:设定迭代终止条件。如果max{|α(k)-α(k-1)|}≤tol成立,则迭代停止,其中tol为预先设定的容许误差(如0.001);否则返回第2步继续;
C3、通过最小二乘法解出
α=(CTC)-1CTb
αi(i=1,…,Nd)即代表了相应单元的损伤程度。
进一步地,所述B步骤中包括如下步骤:
B1、利用传感器测量结构物损伤后的结构动力响应数据,并将数据存储入存储器中;
本发明在已知损伤位置时可以利用结构实测频率的变化来评估损伤的严重程度。该方法对损伤单元设立损伤程度系数,基于基准有限元模型和损伤结构有限的低阶实测模态频率,利用结构特征值建立了方程组,通过迭代求解,可以准确求解该损伤程度系数。与现有技术相比,本发明的优点和积极效果在于:(1)本发明不需要结构的振型信息,避开了实测振型完备性、质量归一化以及识别精度较低的缺点;(2)本发明仅利用损伤前后有限的低阶频率信息进行损伤程度评估,这是非常方便的。一方面,目前振型的识别精度低于频率的识别精度,这样可以把模态参数识别误差对损伤程度评估的影响降至最低;另一方面,频率反映了结构的整体状况,可以在结构的任意部位(当然以方便为准)布置少量(甚至一个)传感器进行振动测试,然后识别其频率,可以有效地避开水下布置传感器或结构内部不宜达到的位置布置传感器造成的不便。本发明适用于单损伤、多种损伤工况,能够准确评估结构的损伤程度,具有一定的实际应用价值。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为:本发明的实施海洋平台结构有限元模型图;
图2为:本发明的实施海洋平台结构的损伤工况示意图;
图3为:本发明的实施海洋平台结构的损伤程度评估效果(工况1);
图4为:本发明的实施海洋平台结构的损伤程度评估效果(工况2);
图5为:本发明的实施海洋平台结构的损伤程度评估效果(工况3);
图6为:本发明的实施海洋平台结构的损伤程度评估效果(工况4);
图7为:本发明的实施海洋平台结构的损伤程度评估效果(工况5);
图8为:本发明的实施海洋平台结构的损伤程度评估效果(工况6);
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
下面结合附图和具体实施方式对本发明作进一步详细的说明。
一、具体算法推导
以K和M表示基准有限元模型(即未损伤结构)的刚度矩阵和质量矩阵,其特征值问题可以表示为
其中ωi和Φi分别为基准有限元模型的第i阶频率和振型。
以K*和M*表示损伤结构的整体刚度矩阵和整体质量矩阵。结构的损伤一般只引起结构刚度的变化,而对结构的质量影响很小,假设M*=M。则其特征值问题为
其中上标T代表转置,即行列互换。
假设结构有Nd个单元发生了损伤,且损伤位置已知,则损伤结构的整体刚度矩阵可以表示为
其中,αn和ln分别第n个损伤单元的损伤程度和第n个损伤单元的单元号。将公式(4)代入公式(3),可得
写成矩阵形式,即
Cα=b (7)
假设已经测得了损伤结构的前m阶频率信息,则C是维数为m×Nd的矩阵,α和b分别是维数为Nd×1和m×1的列向量。通过最小二乘法求解公式(7),得到
α=(CTC)-1CTb (8)
α即为结构单元的损伤程度。
需要说明的是,利用公式(8)求解结构的损伤程度需要用到损伤结构完备的振型信息,这在实际应用中是不可能的。为了克服该限制条件,本发明提出了仅利用实测的低阶频率信息,采用迭代技术来计算结构的损伤程度,具体步骤如下:
第1步:迭代初始赋值。k=0,假设损伤程度为0,即α(0)=0,Φ*(0)=Φ,上标“0”代表迭代初值。
第2步:迭代开始,k=k+1。由α=(CTC)-1CTb计算α(k);
第4步:设定迭代终止条件。如果max{|α(k)-α(k-1)|}≤tol成立,则迭代停止,其中tol为预先设定的容许误差(如0.001);否则返回第2步继续。
二、建立海洋平台有限元模型
本算例模拟研究的海洋平台结构由柱、梁、水平斜撑、竖向斜撑等构件组成,共68个单元,如图1所示。利用MATLAB编写有限元程序建立有限元模型,作为未损伤结构的基准有限元模型。然后再模拟不同的损伤工况,得出模拟实测的模态参数。本算例模拟了六种损伤工况,包括各种构件单个损伤,多损伤以及不同程度的损伤。损伤位置示意图如图2所示。
三、损伤程度评估分析
基于前三阶实测频率,利用本发明的方法进行损伤程度评估,各工况说明如下。
损伤工况一,61号单元为斜撑杆件,发生刚度损失30%,用本发明方法可以准确地评估损伤程度,迭代次数为4次,见图3所示。
损伤工况二,35号单元为横撑杆件,发生刚度损失50%,用本发明方法可以准确地评估损伤程度,迭代次数为5次,见图4所示。
损伤工况三,17号单元为立柱,发生刚度损失40%,用本发明方法可以准确地评估损伤程度,迭代次数为10次,见图5所示。
损伤工况四,35号横撑单元和61号斜撑单元分别发生刚度损失30%和40%,用本发明方法可以准确地评估损伤程度,迭代次数为4次,见图6所示。验证了该方法对发生两处损伤(横撑和斜撑)可以很好的作出评估。
损伤工况五,17号立柱单元与61号斜撑单元同时发生损伤,刚度损失分别为20%和40%,用本发明方法可以准确地评估损伤程度,迭代次数为5次,见图7所示。验证了该方法对发生两处损伤(立柱和斜撑)可以很好的作出评估。
损伤工况六,17号立柱单元、35号横撑单元与61号斜撑单元同时发生损伤,刚度损失分别为30%、20%和40%,用本发明方法可以准确地评估各个单元的损伤程度,迭代次数为5次,见图8所示。验证了该方法对多处损伤工况的损伤程度评估能力。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (2)
1.一种基于频率变化的结构损伤程度评估方法,其特征在于,包括如下步骤:
A、建立基准有限元模型;
B、识别损伤结构模态频率
C、损伤程度评估:
其中K为基准有限元模型的总体刚度矩阵;和αn分别为第n个损伤单元的单元刚度矩阵和损伤程度;Nd为损伤单元的个数;Φi和分别为基准有限元模型和损伤结构的第i阶振型;ωi和分别为基准有限元模型和损伤结构的第i阶频率;上标T代表矩阵或向量的转置;
第1步:迭代初始赋值,k=0,假设损伤程度为0,即α(0)=0,Φ*(0)=Φ,上标“0”代表迭代初值;
第2步:迭代开始,k=k+1,由α=(CTC)-1CTb计算α(k);
第4步:设定迭代终止条件,如果max{|α(k)-α(k-1)|}≤tol成立,则迭代停止,其中tol为预先设定的容许误差(如0.001);否则返回第2步继续;
C3、通过最小二乘法解出
α=(CTC)-1CTb
αi(i=1,…,Nd)即代表了相应单元的损伤程度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102445730A CN101915733B (zh) | 2010-07-30 | 2010-07-30 | 一种基于频率变化的结构损伤程度评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102445730A CN101915733B (zh) | 2010-07-30 | 2010-07-30 | 一种基于频率变化的结构损伤程度评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101915733A true CN101915733A (zh) | 2010-12-15 |
CN101915733B CN101915733B (zh) | 2012-07-25 |
Family
ID=43323295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102445730A Active CN101915733B (zh) | 2010-07-30 | 2010-07-30 | 一种基于频率变化的结构损伤程度评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101915733B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102495134A (zh) * | 2011-12-08 | 2012-06-13 | 大连理工大学 | 基于附加质量的桁架结构损伤定位方法 |
CN102759573A (zh) * | 2012-07-25 | 2012-10-31 | 中国海洋石油总公司 | 基于频率变化的结构损伤定位与损伤程度的评估方法 |
CN103076394A (zh) * | 2013-01-06 | 2013-05-01 | 中国海洋石油总公司 | 基于振动识别频率和振型综合的海洋平台安全评定的方法 |
CN103592220A (zh) * | 2013-11-14 | 2014-02-19 | 西安交通大学 | 一种鼓风机叶轮在位裂纹检测方法 |
CN103852277A (zh) * | 2013-09-23 | 2014-06-11 | 中国海洋大学 | 老龄海洋平台阶段性结构损伤评估方法 |
CN105005695A (zh) * | 2015-07-14 | 2015-10-28 | 中国海洋大学 | 一种用于时域疲劳分析的波浪散布图组块等效方法 |
CN106017906A (zh) * | 2016-07-26 | 2016-10-12 | 南京航空航天大学 | 一种冷却塔模态测试方法 |
CN107687872A (zh) * | 2017-08-14 | 2018-02-13 | 深圳市智能机器人研究院 | 基于动态模型更新的桥梁结构健康状态监测方法及系统 |
CN108226399A (zh) * | 2018-01-23 | 2018-06-29 | 中冶建筑研究总院有限公司 | 一种基于模态参数的张弦结构损伤组合识别方法 |
CN113987873A (zh) * | 2021-10-25 | 2022-01-28 | 王永亮 | 板体平直裂纹损伤识别的有限元自适应求解方法和系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5327358A (en) * | 1991-08-07 | 1994-07-05 | The Texas A&M University System | Apparatus and method for damage detection |
CN1740770A (zh) * | 2005-09-18 | 2006-03-01 | 中国海洋大学 | 大型结构物无损伤的健康检测方法 |
CN101013449A (zh) * | 2006-12-31 | 2007-08-08 | 中国海洋大学 | 利用交叉模型交叉模态的结构模型修正方法 |
-
2010
- 2010-07-30 CN CN2010102445730A patent/CN101915733B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5327358A (en) * | 1991-08-07 | 1994-07-05 | The Texas A&M University System | Apparatus and method for damage detection |
CN1740770A (zh) * | 2005-09-18 | 2006-03-01 | 中国海洋大学 | 大型结构物无损伤的健康检测方法 |
CN101013449A (zh) * | 2006-12-31 | 2007-08-08 | 中国海洋大学 | 利用交叉模型交叉模态的结构模型修正方法 |
Non-Patent Citations (6)
Title |
---|
《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑》 20060115 方辉 海洋平台结构损伤诊断方法研究 A010-22 1-2 , 第1期 2 * |
《中国优秀博硕士学位论文全文数据库(硕士) 工程科技II辑》 20021215 王爱文 建筑结构动力特性参数反演方法研究 C038-54 1-2 , 第2期 2 * |
《中国博士学位论文全文数据库 工程科技II辑》 20091115 王俊荣 海洋平台结构物损伤检测与模型修正方法研究 C034-73 1-2 , 第11期 2 * |
《振动与冲击》 20071231 刘爱霞等 一种海洋平台广义模型修正方法的研究 112-116 1-2 第26卷, 第1期 2 * |
《海洋工程》 20060531 刁延松等 基于不完备模态信息的海洋平台损伤诊断研究 14-20 1-2 第24卷, 第2期 2 * |
《青岛海洋大学学报》 20030731 王树青等 基于有限元模型和模态参数的损伤识别 650-656 1-2 第33卷, 第4期 2 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102495134B (zh) * | 2011-12-08 | 2013-08-07 | 大连理工大学 | 基于附加质量的桁架结构损伤定位方法 |
CN102495134A (zh) * | 2011-12-08 | 2012-06-13 | 大连理工大学 | 基于附加质量的桁架结构损伤定位方法 |
CN102759573B (zh) * | 2012-07-25 | 2015-09-09 | 中国海洋石油总公司 | 基于频率变化的结构损伤定位与损伤程度的评估方法 |
CN102759573A (zh) * | 2012-07-25 | 2012-10-31 | 中国海洋石油总公司 | 基于频率变化的结构损伤定位与损伤程度的评估方法 |
CN103076394A (zh) * | 2013-01-06 | 2013-05-01 | 中国海洋石油总公司 | 基于振动识别频率和振型综合的海洋平台安全评定的方法 |
CN103076394B (zh) * | 2013-01-06 | 2015-06-03 | 中国海洋石油总公司 | 基于振动识别频率和振型综合的海洋平台安全评定的方法 |
CN103852277A (zh) * | 2013-09-23 | 2014-06-11 | 中国海洋大学 | 老龄海洋平台阶段性结构损伤评估方法 |
CN103852277B (zh) * | 2013-09-23 | 2016-04-27 | 中国海洋大学 | 老龄海洋平台阶段性结构损伤评估方法 |
CN103592220A (zh) * | 2013-11-14 | 2014-02-19 | 西安交通大学 | 一种鼓风机叶轮在位裂纹检测方法 |
CN103592220B (zh) * | 2013-11-14 | 2015-05-20 | 西安交通大学 | 一种鼓风机叶轮在位裂纹检测方法 |
CN105005695A (zh) * | 2015-07-14 | 2015-10-28 | 中国海洋大学 | 一种用于时域疲劳分析的波浪散布图组块等效方法 |
CN105005695B (zh) * | 2015-07-14 | 2018-03-09 | 中国海洋大学 | 一种用于时域疲劳分析的波浪散布图组块等效方法 |
CN106017906A (zh) * | 2016-07-26 | 2016-10-12 | 南京航空航天大学 | 一种冷却塔模态测试方法 |
CN107687872A (zh) * | 2017-08-14 | 2018-02-13 | 深圳市智能机器人研究院 | 基于动态模型更新的桥梁结构健康状态监测方法及系统 |
CN107687872B (zh) * | 2017-08-14 | 2020-10-27 | 深圳市智能机器人研究院 | 基于动态模型更新的桥梁结构健康状态监测方法及系统 |
CN108226399A (zh) * | 2018-01-23 | 2018-06-29 | 中冶建筑研究总院有限公司 | 一种基于模态参数的张弦结构损伤组合识别方法 |
CN113987873A (zh) * | 2021-10-25 | 2022-01-28 | 王永亮 | 板体平直裂纹损伤识别的有限元自适应求解方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101915733B (zh) | 2012-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101915733B (zh) | 一种基于频率变化的结构损伤程度评估方法 | |
Wickramasinghe et al. | Vibration characteristics and damage detection in a suspension bridge | |
CN102759573B (zh) | 基于频率变化的结构损伤定位与损伤程度的评估方法 | |
Esfandiari et al. | Structural model updating using frequency response function and quasi-linear sensitivity equation | |
Bernal | Load vectors for damage localization | |
Jensen et al. | The use of updated robust reliability measures in stochastic dynamical systems | |
CN106896156A (zh) | 逐跨均匀荷载面曲率差的梁结构损伤识别方法 | |
CN103076394B (zh) | 基于振动识别频率和振型综合的海洋平台安全评定的方法 | |
CN115983062B (zh) | 一种基于有限元模型修正的高拱坝地震损伤评估方法及系统 | |
CN100575924C (zh) | 大型结构物无损伤的健康检测方法 | |
CN101034053A (zh) | 一种基于模态信息的结构损伤诊断新方法 | |
CN105064420A (zh) | 一种基于结构残余模态力的高桩码头基桩损伤诊断方法 | |
Tabrizian et al. | A New Damage Detection Method: Big Bang‐Big Crunch (BB‐BC) Algorithm | |
CN103884593B (zh) | 基于少量静态测量数据的结构损伤定位与程度识别方法 | |
He-zhen et al. | Damage Localization of Offshore Platforms Under Ambient Excitation" YANG He-zhen (ffili)", LI Hua-jun (% 43) and WANG Shu-qing (Eff) | |
CN106529049A (zh) | 一种用于模型修正的应变模态相关性评价方法 | |
CN104133960A (zh) | 一种改进的静力传感器的优化布设方法 | |
CN110008520B (zh) | 基于位移响应协方差参数和贝叶斯融合的结构损伤识别方法 | |
CN106338372A (zh) | 基于残差应变能的海洋平台损伤定位方法及系统 | |
Izadi et al. | Structural damage identification by a sensitivity-based finite element model updating method using transmissibility function data | |
Jiang et al. | Structural stiffness identification of traditional mortise-tenon joints based on statistical process control chart | |
Abrahamsen et al. | Hydroelastic response of concrete shells during impact on calm water | |
CN102998367B (zh) | 一种基于虚拟派生结构的损伤识别方法 | |
Dincal et al. | Nondestructive damage detection in Euler–Bernoulli beams using nodal curvatures—Part I: Theory and numerical verification | |
Hu et al. | Mechanistically-informed damage detection using dynamic measurements: Extended constitutive relation error |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |