CN107016166A - 一种基于模态应变能的新型结构刚度损伤定位方法 - Google Patents
一种基于模态应变能的新型结构刚度损伤定位方法 Download PDFInfo
- Publication number
- CN107016166A CN107016166A CN201710137256.0A CN201710137256A CN107016166A CN 107016166 A CN107016166 A CN 107016166A CN 201710137256 A CN201710137256 A CN 201710137256A CN 107016166 A CN107016166 A CN 107016166A
- Authority
- CN
- China
- Prior art keywords
- strain energy
- unit
- modal strain
- stiffness
- modal
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本公开涉及一种基于模态应变能的新型结构刚度损伤定位方法,所述方法包括对分析对象建立有限元模型;对分析对象进行模态测试,测得其前五阶模态振型;基于有限元求解器得到模型的模态计算结果,并与实验测得结果进行匹配;将各单元刚度按照节点位置填入总刚度,得到各单元的刚度贡献值;分别求解出实验与有限元模型中各单元的模态应变能;按照阶次的顺序,组装每个单元的各阶模态应变能为其新的结构刚度损伤定位方法。
Description
技术领域
本发明属于结构建模及评估领域,特别是一种基于模态应变能的新型结构刚度损伤定位方法。
背景技术
结构损伤前后模态应变能变化是较常用的损伤识别算法之一,该方法是基于如下的假设:将结构分解为一系列单元,计算结构损伤前后每个单元的模态应变能变化,而部分模态振型在结构损伤附近发生局部突变,故模态应变能在结构中的分布将发生变化,所以可以通过比较每一单元模态应变能的变化来进行结构损伤识别。
但在实际运用中发现,即使在测量精度足够高,并正确计算的情况下,每阶模态应变能对损伤进行敏感程度依然相差很大。如果将各阶模态应变能变化结果做平均,对损伤敏感的模态识别精确度将会被不敏感模态所降低;如果对各阶模态应变能分别依次判断,则会得到很多错误的损伤单元。
如何通过各阶模态应变能的变化来确定损伤的单元,是该方法从理论走向实际应用中必须解决的问题。
发明内容
本发明针对现有技术的不足,公开了一种基于模态应变能的新型结构刚度损伤定位方法,能够准确定位模型中存在损伤的单元,并为后续基于模型的结构状态评估、剩余寿命计算等工作提供了有力保证。
本发明公开了一种基于模态应变能的新型结构刚度损伤定位方法,具体步骤如下:
1)建立分析对象的有限元模型,在有限元软件中将结构划分为单元结构;
2)数据采集:分析结构中的易受损部分,设计传感器排布方案,保证测点与有限元模型的节点相重合,后测量并记录所述分析对象各点的振动,并汇总求得所述分析对象的各阶模态振型;
3)结构有限元分析及模型修正:基于有限元求解器计算模型的各阶模态振型的频率与振型,以测得的前5阶模态频率为目标值,应用基于模态频率的灵敏度修正方法,对模型的初始物理参数进行修正,并计算各单元的局部刚度矩阵;
4)各单元刚度贡献值的计算:将相邻单元间的连接点记为节点,将所述各单元的局部刚度矩阵按照节点的位置,装填入整体模型的刚度矩阵,从而获得各单元的刚度贡献值;
5)模态应变能的计算:通过所述整体模型的阵型计算得到有限元模型以及实验模型的模态应变能(Modal Strain Energy,MSE);
6)根据所述单元各阶模态应变能向量得到结构刚度损伤位置。
所述有限元模型以及实验模型的模态应变能经由下式计算得出:
其中,{Φr a}、{Φr e}分别表示有限元、实验模型第r阶的模态振型,[K]j表示第j个单元的刚度贡献值,MSEjr a、MSEjr e分别为有限元、实验模型第j个单元的第r阶模态应变能。
进一步的,步骤6)中,为模型中每个单元定义其模态应变能向量,所述模态应变能向量由此单元各阶模态应变能按照模态阶次排列而成,即:
MSEvecj a={MSEj1 a,MSEj2 a,...,MSEjr a} (3)
MSEvecj e={MSEj1 e,MSEj2 e,...,MSEjr e} (4)
其中,MSEvecj a,MSEvecj e分别为第j个单元的有限元、实验模态应变能向量,因此定义新的结构刚度损伤定位方法如下:
MSESFj即为第j个单元的新型结构刚度损伤定位指标,体现分析对象的刚度损伤位置,称为第j个单元的模态应变能比例因子。
进一步的,所述的步骤1)中的单元结构为板单元或实体单元。
进一步的,所述的步骤2)中所选的振型测点数量根据实际情况调整,在被测物体的主振方向上测点数量为单数,且左右对称分布,使其能以较少的传感器数量完备描述被测物体的结构振型。
进一步的,所述的步骤4)中,所述单元中各个节点按照其几何位置,将刚度填充至整体刚度矩阵中,以得到该单元的刚度贡献值。
进一步的,所述的步骤5)中,计算有限元模型、实验模型的模态应变能时,使用的刚度矩阵均为无损情况下的刚度矩阵。
所述的步骤6)中,MSESFj的值越趋近于1,则代表第j个单元存在损伤的概率越低。
本发明的相对于现有技术的有益效果在于:
(1)通过选择合适的单元类型,获得的结构有限元模型为损伤识别计算提供了基准模型。
(2)通过对不同单元之间,新型结构刚度损伤定位方法的对比,能够完整利用损伤敏感模态的损伤信息,从而准确识别出含有损伤的单元。
(3)得到损伤单元的编号及位置,有利于后续基于模型的结构状态评估、剩余寿命计算等工作的开展。
附图说明
图1是本发明一种基于模态应变能的新型结构刚度损伤定位方法的实例结构有限元模型;
图2是本发明一种基于模态应变能的新型结构刚度损伤定位方法的损伤模型与无损模型的前10阶MAC匹配图;
图3是本发明一种基于模态应变能的新型结构刚度损伤定位方法的损伤单元序号识别结果图;
具体实施方式
本发明提供一种基于模态应变能的新型结构刚度损伤定位方法,为使本发明的目的,技术方案及效果更加清楚,明确,以及参照附图并举实例对本发明进一步详细说明。应当理解,此处所描述的具体实施仅用以解释本发明,并不用于限定本发明。
本实施例采用一块长50mm,宽50mm,厚2mm的矩形钢板作为研究对象,以说明所提方法的合理性及有效性。验证实例为矩形钢板,但不仅限于此类结构和如此尺寸。实例的具体实施步骤如下:
1)建立矩形板的有限元模型并计算位移模态:在MSC.Patran软件中,采用二维四边形板单元建立矩形板的有限元模型。将相邻单元间的连接点记为节点,模型中共有36个节点,25个单元,单元编号为1~25,如图1所示。将图中阴影部分所示的19号单元弹性模量降低5%作为损伤单元,并将包含此损伤单元的模型作为损伤模型。基于MSC.Nastran软件求解模型各阶模态,获得矩形板的前10阶频率以及对应的位移振型等信息。
2)试验方案设计及分析:根据实验测量的可行性,通常在测量板状结构时最为注重的是垂直于板平面方向的振动位移,因此为模拟实验测量,只选用损伤模型计算结果的Z向结果作为测量结果,并在各节点的测量结果中加入大小为理论值10%的随机误差,以模拟实验中的测量误差。
3)提取有限元仿真应变模态:基于matlab平台,编写有限元计算结果的读写程序,读出矩形板前10阶模态振型以及各单元的刚度矩阵,并按前文所述方法得到各单元的刚度贡献值。
4)相关性分析:采用模态置信准则计算仿真分析和试验测试所获取的模态振型的相关性,MAC计算公式如下:
MAC矩阵对角线值越接近于1,非对角线元素越接近于0,说明损伤模型与无损模型的相关性越好,其计算结果如图2所示。
本专利实例中损伤模型和无损模型的频率误差及MAC矩阵对角线元素值如表1所示。
表1修正前应变模态匹配
从表1以及图2数据可看出,有限元计算得到的各阶固有频率与实验测得的值误差在1%以内,表明有限元模型的初始弹性模量与实际值已较为接近。
5)有限元模型初始弹性模量的修正:以实验所测得的前10阶固有频率为目标值,对有限元模型的初始弹性模量进行灵敏度修正,此时可将次修正看作一个优化问题:
其中分别为实验测得、有限元计算所得的第r阶固有频率值,E为弹性模量,VLB和VUB分别为其下限和上限。将公式(12)进行Taylor一阶展开,可得
R(E)=S{Δp} (8)
{ΔE}=S+R(E) (9)
其中S为固有频率的一阶灵敏度矩阵,{ΔE}即为弹性模量的修正值。使用matlab平台编写灵敏度修正的程序,进行迭代修正,完成对有限元模型初始弹性模量的修正。
6)选取计算结果的前10阶模态振型,计算各个单元在各模态下的模态应变能,其中第j个单元的第r阶理论、实际模态应变能值分别表示为MSEjr a、MSEjr e,则按照步骤(6),第j个单元的模态应变能向量为:
MSEvecj a={MSEj1 a,MSEj2 a,...,MSEj10 a} (10)
MSEvecj e={MSEj1 e,MSEj2 e,...,MSEj10 e} (11)
7)根据步骤(6)中所述,第j个单元的模态应变能比例因子MSESFj计算步骤如下:
图3为各单元损伤前后模态应变能比例因子值计算结果。如图所示,第19号单元的模态应变能比例因子明显低于其它单元,证明此结构刚度损伤定位方法是有效的。
综上,本发明具有以下优点:
对稳定结构统计可变性小,抗噪性强,能够快速识别结构损伤的位置,克服了原有模态应变能损伤因子的不足;
本发明运算过程简单,可实施性好,适合检验现场操作。
以上所述仅为本发明的实施方式,并非因此限制本发明的专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,或直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。
Claims (8)
1.一种基于模态应变能的新型结构刚度损伤定位方法,其特征在于,具体步骤如下:
1)建立分析对象的有限元模型,在有限元软件中将结构划分为单元结构;
2)数据采集:分析结构中的易受损部分,设计传感器排布方案,保证测点与有限元模型的节点相重合,后测量并记录所述分析对象各点的振动,并汇总求得所述分析对象的各阶模态振型;
3)结构有限元分析及模型修正:基于有限元求解器计算模型的各阶模态振型的频率与振型,以测得的前5阶模态频率为目标值,应用基于模态频率的灵敏度修正方法,对模型的初始物理参数进行修正,并计算各单元的局部刚度矩阵;
4)各单元刚度贡献值的计算:将相邻单元间的连接点记为节点,将所述各单元的局部刚度矩阵按照节点的位置,装填入整体模型的刚度矩阵,从而获得各单元的刚度贡献值;
5)模态应变能的计算:通过所述整体模型的阵型计算得到有限元模型以及实验模型的模态应变能(Modal Strain Energy,MSE);
6)根据所述单元各阶模态应变能向量得到结构刚度损伤位置。
2.根据权利要求1所述的一种基于模态应变能的新型结构刚度损伤定位方法,其特征在于,所述有限元模型以及实验模型的模态应变能经由下式计算得出:
其中,{Φr a}、{Φr e}分别表示有限元、实验模型第r阶的模态振型,[K]j表示第j个单元的刚度贡献值,MSEjr a、MSEjr e分别为有限元、实验模型第j个单元的第r阶模态应变能。
3.根据权利要求1所述的一种基于模态应变能的新型结构刚度损伤定位方法,其特征在于,步骤(6)中,为模型中每个单元定义其模态应变能向量,所述模态应变能向量由此单元各阶模态应变能按照模态阶次排列而成,即:
MSEvecj a={MSEj1 a,MSEj2 a,...,MSEjr a} (3)
MSEvecj a={MSEj1 a,MSEj2 a,...,MSEjr a} (4)
其中,MSEvecj a,MSEvecj e分别为第j个单元的有限元、实验模态应变能向量,因此定义新的结构刚度损伤定位方法如下:
MSESFj即为第j个单元的新型结构刚度损伤定位方法,体现分析对象的刚度损伤位置,称为第j个单元的模态应变能比例因子。
4.根据权利要求1所述的一种基于模态应变能的新型结构刚度损伤定位方法,其特征在于,所述的步骤(1)中的单元结构为板单元或实体单元。
5.根据权利要求1所述的一种基于模态应变能的新型结构刚度损伤定位方法,其特征在于,所述的步骤2)中所选的振型测点数量根据实际情况调整,在被测物体的主振方向上测点数量为单数,且左右对称分布,使其能以较少的传感器数量完备描述被测物体的结构振型。
6.根据权利要求1所述的一种基于模态应变能的新型结构刚度损伤定位方法,其特征在于,所述的步骤4)中,所述单元中各个节点按照其几何位置,将刚度填充至整体刚度矩阵中,以得到该单元的刚度贡献值。
7.根据权利要求1所述的一种基于模态应变能的新型结构刚度损伤定位方法,其特征在于,所述的步骤5)中,计算有限元模型、实验模型的模态应变能时,使用的刚度矩阵均为无损情况下的刚度矩阵。
8.根据权利要求1所述的一种基于模态应变能的新型结构刚度损伤定位方法,其特征在于,所述的步骤6)中,MSESFj的值越趋近于1,则代表第j个单元存在损伤的概率越低。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710137256.0A CN107016166A (zh) | 2017-03-09 | 2017-03-09 | 一种基于模态应变能的新型结构刚度损伤定位方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710137256.0A CN107016166A (zh) | 2017-03-09 | 2017-03-09 | 一种基于模态应变能的新型结构刚度损伤定位方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107016166A true CN107016166A (zh) | 2017-08-04 |
Family
ID=59440215
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710137256.0A Pending CN107016166A (zh) | 2017-03-09 | 2017-03-09 | 一种基于模态应变能的新型结构刚度损伤定位方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107016166A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108287970A (zh) * | 2018-01-31 | 2018-07-17 | 东南大学 | 基于二维正交各向异性复合材料板的热模态对结构参数的灵敏度分析方法 |
CN108981783A (zh) * | 2018-05-16 | 2018-12-11 | 大连理工大学 | 基于灵敏度体积最大化的虚拟质量优化方法 |
CN110108430A (zh) * | 2019-05-15 | 2019-08-09 | 东北大学 | 一种纤维增强复合材料梁的高精度损伤定位方法 |
CN110362855A (zh) * | 2019-05-30 | 2019-10-22 | 中国海洋大学 | 一种基于交叉模态能敏感度的模态优选及模型修正方法 |
CN111261239A (zh) * | 2019-12-11 | 2020-06-09 | 江苏理工学院 | 基于损伤力学的复合材料结构损伤识别方法 |
CN111400940A (zh) * | 2018-12-30 | 2020-07-10 | 长城汽车股份有限公司 | 车架后端结构连接点确定方法 |
CN114969971A (zh) * | 2022-05-24 | 2022-08-30 | 南方天合底盘系统有限公司 | 用于浮动式卡钳盘式制动器低频制动噪音的解决方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101650758A (zh) * | 2009-09-14 | 2010-02-17 | 中南大学 | 车辆结构分析中构造车体板梁结构的方法 |
CN102928247A (zh) * | 2012-10-26 | 2013-02-13 | 重庆大学 | 一种等效模态应变能损伤识别方法 |
CN103344448A (zh) * | 2013-06-26 | 2013-10-09 | 中国路桥工程有限责任公司 | 一种桥梁结构损伤识别方法和系统 |
CN103852277A (zh) * | 2013-09-23 | 2014-06-11 | 中国海洋大学 | 老龄海洋平台阶段性结构损伤评估方法 |
CN103902834A (zh) * | 2014-04-14 | 2014-07-02 | 重庆大学 | 一种基于岭估计和l曲线法的结构损伤识别方法 |
CN104462785A (zh) * | 2014-11-12 | 2015-03-25 | 重庆大学 | 一种两阶段式建筑框架结构损伤检测方法 |
-
2017
- 2017-03-09 CN CN201710137256.0A patent/CN107016166A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101650758A (zh) * | 2009-09-14 | 2010-02-17 | 中南大学 | 车辆结构分析中构造车体板梁结构的方法 |
CN102928247A (zh) * | 2012-10-26 | 2013-02-13 | 重庆大学 | 一种等效模态应变能损伤识别方法 |
CN103344448A (zh) * | 2013-06-26 | 2013-10-09 | 中国路桥工程有限责任公司 | 一种桥梁结构损伤识别方法和系统 |
CN103852277A (zh) * | 2013-09-23 | 2014-06-11 | 中国海洋大学 | 老龄海洋平台阶段性结构损伤评估方法 |
CN103902834A (zh) * | 2014-04-14 | 2014-07-02 | 重庆大学 | 一种基于岭估计和l曲线法的结构损伤识别方法 |
CN104462785A (zh) * | 2014-11-12 | 2015-03-25 | 重庆大学 | 一种两阶段式建筑框架结构损伤检测方法 |
Non-Patent Citations (7)
Title |
---|
V. B. DAWARI,ETC: "Structural Damage Identification Using Modal Strain Energy Method", 《ADVANCES IN STRUCTURAL ENGINEERING》 * |
YINGCHAO LI,ETC: "An improved modal strain energy method for damage detection in offshore platform structures", 《JOURNAL OF MARINE SCIENCE AND APPLICATION》 * |
YOUNG‐JIN CHA,ETC: "Structural Damage Detection Using Modal Strain Energy and Hybrid Multiobjective Optimization", 《SPECIAL ISSUE: COMPUTATIONAL INTELLIGENCE IN STRUCTURAL ENGINEERING AND MECHANICS》 * |
张保强,等: "基于模态频率和有效模态质量的有限元模型修正", 《振动与冲击》 * |
徐典: "结构损伤识别方法与传感器优化布置研究", 《中国博士论文全文数据库 工程科技II辑》 * |
方圣恩: "基于有限元模型修正的结构损伤识别方法研究", 《中国博士论文全文数据库 工程科技II辑》 * |
滕亮: "桥梁结构损伤识别指标比选及损伤程度识别方法研究", 《中国优秀硕士论文全文数据库 工程科技II辑》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108287970A (zh) * | 2018-01-31 | 2018-07-17 | 东南大学 | 基于二维正交各向异性复合材料板的热模态对结构参数的灵敏度分析方法 |
CN108287970B (zh) * | 2018-01-31 | 2019-01-29 | 东南大学 | 基于二维正交各向异性复合材料板的热模态对结构参数的灵敏度分析方法 |
CN108981783A (zh) * | 2018-05-16 | 2018-12-11 | 大连理工大学 | 基于灵敏度体积最大化的虚拟质量优化方法 |
CN111400940A (zh) * | 2018-12-30 | 2020-07-10 | 长城汽车股份有限公司 | 车架后端结构连接点确定方法 |
CN111400940B (zh) * | 2018-12-30 | 2023-04-07 | 长城汽车股份有限公司 | 车架后端结构连接点确定方法 |
CN110108430A (zh) * | 2019-05-15 | 2019-08-09 | 东北大学 | 一种纤维增强复合材料梁的高精度损伤定位方法 |
CN110108430B (zh) * | 2019-05-15 | 2020-10-27 | 东北大学 | 一种纤维增强复合材料梁的高精度损伤定位方法 |
CN110362855A (zh) * | 2019-05-30 | 2019-10-22 | 中国海洋大学 | 一种基于交叉模态能敏感度的模态优选及模型修正方法 |
CN111261239A (zh) * | 2019-12-11 | 2020-06-09 | 江苏理工学院 | 基于损伤力学的复合材料结构损伤识别方法 |
CN111261239B (zh) * | 2019-12-11 | 2023-07-11 | 江苏理工学院 | 基于损伤力学的复合材料结构损伤识别方法 |
CN114969971A (zh) * | 2022-05-24 | 2022-08-30 | 南方天合底盘系统有限公司 | 用于浮动式卡钳盘式制动器低频制动噪音的解决方法 |
CN114969971B (zh) * | 2022-05-24 | 2023-04-18 | 南方天合底盘系统有限公司 | 用于浮动式卡钳盘式制动器低频制动噪音的解决方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107016166A (zh) | 一种基于模态应变能的新型结构刚度损伤定位方法 | |
CN106529055A (zh) | 一种基于应变模态振型相关性的模型修正方法 | |
CN107330176B (zh) | 一种基于结构模态估计的应变计与加速度计联合布设方法 | |
CN101673316B (zh) | 基于混合监测来识别需调整索力的支承索的方法 | |
CN103792143B (zh) | 一种单轴拉伸全程真应力应变曲线的快速获取方法 | |
CN103776416B (zh) | 一种基于应变传感器的天线结构变形的间接测量方法 | |
CN108875178B (zh) | 用于减小结构模态识别不确定性的传感器布置方法 | |
CN105825014A (zh) | 一种基于车桥耦合分析的车桥安全健康评估系统及方法 | |
CN101915733B (zh) | 一种基于频率变化的结构损伤程度评估方法 | |
CN104182585B (zh) | 硅胶弹性体手感分析及疲劳寿命预测的有限元模拟方法 | |
CN106529049A (zh) | 一种用于模型修正的应变模态相关性评价方法 | |
CN102759573B (zh) | 基于频率变化的结构损伤定位与损伤程度的评估方法 | |
CN106777549B (zh) | 一种面向荷载试验的桥梁多尺度有限元模拟方法 | |
CN106354955A (zh) | 一种基于磨机振动模态参数的滑动轴承刚度识别方法 | |
CN102567578B (zh) | 航天器振动试验夹具评价系统 | |
CN104992002A (zh) | 一种面向智能蒙皮天线的应变传感器布局方法 | |
CN106997410A (zh) | 一种基于模态应变能的损伤发生的判断方法 | |
CN105205262A (zh) | 一种基于二次响应面反演的随机模型修正方法 | |
CN115455793A (zh) | 基于多尺度模型修正的高层结构复杂构件受力分析方法 | |
CN110110406B (zh) | 一种基于Excel计算平台实现LS-SVM模型的边坡稳定性预测方法 | |
CN110069822B (zh) | 一种叶片动应变测量的传感器布置方法 | |
CN109726497B (zh) | 一种空间钢结构温度场的获取方法 | |
CN106354944A (zh) | 一种基于振动模态的磨机基础边界支承刚度的识别方法 | |
CN110516297A (zh) | 基于长标距时频域信息的多尺度有限元模型修正方法 | |
CN115687854B (zh) | 高精度土样参数测量方法及其系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170804 |
|
RJ01 | Rejection of invention patent application after publication |