CN106529049A - 一种用于模型修正的应变模态相关性评价方法 - Google Patents
一种用于模型修正的应变模态相关性评价方法 Download PDFInfo
- Publication number
- CN106529049A CN106529049A CN201611019190.7A CN201611019190A CN106529049A CN 106529049 A CN106529049 A CN 106529049A CN 201611019190 A CN201611019190 A CN 201611019190A CN 106529049 A CN106529049 A CN 106529049A
- Authority
- CN
- China
- Prior art keywords
- strain
- finite element
- mode
- test
- 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]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种用于模型修正的应变模态相关性评价方法,具体步骤为:1)建立结构的有限元模型并分析:在有限元软件中,采用合适的单元对结构进行建模,并基于有限元求解器计算结构中的应变模态;2)试验测试方案的设计;3)试验应变模态识别:通过模态识别软件或是基于计算机语言平台,编写模态识别程序,从而识别出试验测试的应变模态参数;4)有限元计算结果保存读取;5)采用模态置信准则进行相关性分析;该方法将应变模态引入相关性评价,具有描述两模型局部差异的能力;该方法可以评价两模型描述应变的能力,并为后续的结构动力学优化设计、结构健康监测,有限元模型修正或结构响应预测等提供了有力的保证。
Description
技术领域
本发明属于结构建模及评估及动力学反问题领域,特别是一种用于模型修正的应变模态相关性评价方法。
背景技术
随着有限元技术的发展以及计算机应用水平的提高,建模及分析技术在结构动力学优化设计、结构健康监测、有限元模型修正及结构响应预测等领域得到广泛的应用。在有限元建模的初期阶段,模型中的参数大多是根据相关标准的名义值或是分析者的经验设定。由于制造、装配等多方面因素的原因,这些赋予模型的初始参数值均与其真实值之间存在偏差,导致由此模型计算出的响应与真实结构响应之间存在一定的偏差。相比于仿真分析而言,试验测试是针对实际工程结构而进行的,在工程实践中认为所测得的响应是准确的。为了能够在结构动力学优化设计、结构健康监测或是结构响应预测中使用有限元模型,就需要对所使用的有限元模型进行评价。在上述的各领域中,现使用较多的方法是基于振动信号的评价,即基于加速度、速度、位移等测试信号识别的结构位移模态信息对有限元模型和试验模型之间的相关性进行评价。
应变作为一种反应结构局部特征的物理量,在试验测试和仿真分析中均易获得,因此在近年也在上述的各大领域中被广泛采用。关于应变响应,较多的学者采用相对应变振型差来描述应变振型的相关性,但是对于航空航天、高铁列车等模态较为密集的结构模型,如何判断两个模型中的模态属于同一阶次,并没有相关的研究和报道。找到试验模型和仿真模型应变模态的相互对应关系是评价两者模型相关性的关键步骤,也是基于模型进行后续的优化设计和响应预测的必要条件。然而,现阶段现有的技术中关于评价应变模态相关性的研究尚未见报道。
发明内容
本发明针对现有技术的不足,公开了一种用于模型修正的应变模态相关性评价方法,该方法能够准确评价仿真分析的应变振型反映实测应变振型的能力,并为基于模型的结构动力学优化设计、结构健康监测、有限元模型修正和结构响应预测给予准确的模型支持;尤其对于模态较为密集的模型,能够为后续的应变响应计算提供一个合理的有限元模型。
本发明是这样实现的:一种用于模型修正的应变模态相关性评价方法,其特征在于,具体步骤如下:
1)建立结构的有限元模型并分析:在有限元软件中,采用合适的单元对结构进行建模,并基于有限元求解器计算结构中的应变模态;
2)试验测试方案的设计:根据步骤1)中有限元计算的结构应变模态振型选取结构上应变敏感的位置和方向,选择应变片贴于各测点处组成应变测试桥路;在结构上施加随机激励,并用信号采集分析仪记录各测点的应变信号;
3)试验应变模态识别:通过模态识别软件或是基于计算机语言平台,编写模态识别程序,从而识别出对应的应变模态参数;
4)有限元计算结果保存读取:基于matlab平台编写读写程序,读取有限元模型计算结果中各单元在不同方向的应变值,并按照列格式存储,从而获得整个结构仿真分析的应变振型;
5)相关性分析:采用模态置信准则,即Modal Assurance Criterion(MAC),MAC计算仿真分析与试验测试得到的应变模态之间的相关性,公式如下:
式中,分别为第i阶试验应变振型以及第j阶仿真应变振型,均为列向量;
进一步,所述的步骤1)中的结构包括薄壁结构、实体结构等,在划分有限元网格时,需要针对结构的具体类型选择合适的单元类型,例如薄壁结构采用二维板单元、实体结构采用三维实体单元等。对结构进行网格划分时应该采用四边形网格或六面体网格。此外,对于网格尺度较大的模型,在对应试验测试位置的周围的网格要进行加密处理,即该区域的网格大小应与测试传感器尺寸近似或更小。
进一步,所述的步骤2)中的试验设计因素还包括:试验测点位置的分布、各测点的测试方向设定、激励的类型及施加方式、信号采集仪器的选取及测试桥路组成方式。
进一步,所述的步骤3)中应变模态参数包括应变模态频率、应变模态振型和模态阻尼比。
进一步,所述的步骤4)中读取的是与试验测试对应的有限元模型中节点位置处的应变值,方向也与试验测试一致。
进一步,所述的步骤5)计算的是仿真分析和试验测试的应变模态振型的相关性,可以是单个方向应变振型的相关性分析,也可是多方向同时进行相关性分析。
本发明相对于现有技术的有益效果在于:在模型相关性的评价中,将应变模态作为响应特征,提出了能够评价应变模态相关性的方法;应变是结构的局部特征量,将应变模态引入相关性评价,具有描述两模型局部差异的能力;通过应变模态相关性分析,可以评价两模型描述应变的能力,并为后续的结构动力学优化设计、结构健康监测,有限元模型修正或结构响应预测等提供了有力的保证。
附图说明
图1是本发明一种用于模型修正的应变模态相关性评价方法的流程图;
图2是本发明一种用于模型修正的应变模态相关性评价方法的实例结构的有限元模型;
图3是本发明一种用于模型修正的应变模态相关性评价方法的实例结构的测试方案;
图4是本发明一种用于模型修正的应变模态相关性评价方法的实例中结构前10阶应变振型相关性结果,即MAC矩阵的三维柱状图。
具体实施方式
本发明提供一种评价应变模态相关性方法,为使本发明的目的、技术方案及效果更加清楚,明确,以及参照附图并举实例对本发明进一步详细说明。应当指出此处所描述的具体实施仅用以解释本发明,并不用于限定本发明。
如图1所示,本发明的一种评价应变模态相关性方法的具体步骤如下:
1)建立结构的有限元模型并分析:在有限元软件中,采用合适的单元对结构进行建模,并基于有限元求解器计算结构中的应变模态;
2)试验测试方案的设计:根据步骤1)中有限元计算结果,在结构上选取响应敏感的位置和方向,选择应变片贴于各测点处组成应变测试桥路;在结构上施加随机激励,并用信号采集分析仪记录各测点的应变信号;
3)试验应变模态识别:通过模态识别软件或是基于计算机语言平台,编写模态识别程序,从而识别出对应的模态参数;
4)有限元计算结果保存读取:基于matlab平台编写读写程序,读取有限元模型计算结果中各单元在不同方向的应变值,并按照列的格式存储,从而获得整个结构仿真分析的应变振型;
5)相关性分析:采用模态置信准则,即Modal Assurance Criterion(MAC),MAC计算仿真分析与试验测试得到的应变模态之间的相关性。
具体的实施例如下:
本发明采用一块长550mm,宽120mm,厚3mm的矩形钢板作为验证对象,具体实施步骤如下:
1)建立矩形钢板的有限元模型并分析:采用二维四边形板单元建立矩形板的有限元模型。单元尺寸取为20mm,模型共有203个节点,168个单元。矩形板在一端固支约束,另一端自由,构造成悬臂板结构,悬臂板的有限元模型示意图如图2所示。基于MSC.Nastran求解获得该悬臂板的应变模态;
2)试验方案设计及分析:根据有限元分析结果选取结构响应敏感的位置及方向,并据此安排测试方案。在实施过程中,在悬臂板的长度方向分布5个测点,宽度方向分布3个测点,即共有15个测点。对于所有测点,均测试其X向和Y向的应变值。测点分布及测试方向如图3所示;
3)应变模态识别:在实施过程中,采用仿真模型代替试验模型,即在悬臂板有限元模型中加10%的噪声的方式构造试验模型,并通过输出与测点相同位置处的应变以获得悬臂板试验应变模态;
4)仿真结果读取:基于matlab平台读取出矩形板的应变值,并按照列模式存储成应变模态矩阵,输出保存;
5)相关性分析:将试验识别结果和有限元仿真分析的读取结果分别代入到(1)式中,即可计算得到仿真和试验应变模态的MAC值,前10阶的应变模态MAC矩阵的对角线元素如表1所示,对应的MAC矩阵三维柱状图如图4所示。
表1悬臂板仿真模型和试验模型MAC匹配
工程中关于MAC的评价指标为,MAC矩阵对角线元素越接近于1、非对角线元素越小说明两者之间的相关性较好。此外,工程中认为MAC矩阵中元素值大于0.7时所对应的两阶模态所表达的振型是一致的。从表1的结果看,仿真分析和试验测试的频率误差绝对平均值为4.71%,MAC矩阵的对角线元素均在0.94以上,MAC矩阵对角线元素的均值为0.968。以上数据表明,悬臂板仿真模型和所构造的试验模型之间相关性较好,但仍存在一定的误差。以上所述仅为本发明的较佳可行实施例,并非因此局限本发明的专利范围,故凡是运用本发明说明书及附图内容所作的等效变化,均包含与本发明的保护范围。
Claims (6)
1.一种用于模型修正的应变模态相关性评价方法,其特征在于,具体步骤如下:
1)建立结构的有限元模型并分析:在有限元软件中,采用合适的单元对结构进行建模,并基于有限元求解器计算结构中的应变模态;
2)试验测试方案的设计:根据步骤1)中有限元计算结果,根据结构应变模态振型选取结构上应变敏感的位置和方向,选择应变片贴于各测点处组成应变测试桥路;在结构上施加随机激励,并用信号采集分析仪记录各测点的应变信号;
3)试验应变模态识别:通过模态识别软件或是基于计算机语言平台,编写模态识别程序,从而识别出试验测试的应变模态参数;
4)有限元计算结果保存读取:基于matlab平台编写读写程序,读取有限元模型计算结果中的应变值,并按照每阶模态振型一列的格式存储,从而获得整个结构仿真分析的应变振型;
5)相关性分析:采用模态置信准则,即Modal Assurance Criterion,MAC,计算仿真分析与试验测试得到的应变模态之间的相关性,公式如下:
式中,分别为第i阶试验应变振型以及第j阶仿真应变振型,均为列向量。
2.根据权利要求1所述的一种用于模型修正的应变模态相关性评价方法,其特征在于,所述的步骤1)中的结构包括薄壁结构、实体结构,对结构进行网格划分时采用四边形网格或六面体网格。
3.根据权利要求2所述的一种用于模型修正的应变模态相关性评价方法,其特征在于,所述的步骤2)中的试验设计因素还包括:试验测点位置的分布、各测点的测试方向设定、激励的类型及施加方式、信号采集仪器的选取及测试桥路组成方式。
4.根据权利要求1所述的一种用于模型修正的应变模态相关性评价方法,其特征在于,所述的步骤3)中应变模态参数包括应变模态频率、应变模态振型和模态阻尼比。
5.根据权利要求1所述的一种用于模型修正的应变模态相关性评价方法,其特征在于,所述的步骤4)中读取的是与试验测试对应的有限元模型中节点位置处的应变值,方向与试验测试一致。
6.根据权利要求1所述的一种用于模型修正的应变模态相关性评价方法,其特征在于,所述的步骤5)中计算的是仿真分析和试验测试的应变模态振型的相关性;所述的相关性分析为单个方向应变振型的相关性分析或者多方向同时进行相关性分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611019190.7A CN106529049A (zh) | 2016-11-18 | 2016-11-18 | 一种用于模型修正的应变模态相关性评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611019190.7A CN106529049A (zh) | 2016-11-18 | 2016-11-18 | 一种用于模型修正的应变模态相关性评价方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106529049A true CN106529049A (zh) | 2017-03-22 |
Family
ID=58352401
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611019190.7A Pending CN106529049A (zh) | 2016-11-18 | 2016-11-18 | 一种用于模型修正的应变模态相关性评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106529049A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107644134A (zh) * | 2017-09-20 | 2018-01-30 | 北京安怀信科技股份有限公司 | 一种动力学验模工具 |
CN107844669A (zh) * | 2017-12-01 | 2018-03-27 | 中国直升机设计研究所 | 一种直升机全机静力有限元模型修正及验证方法 |
CN109883389A (zh) * | 2019-03-22 | 2019-06-14 | 西安交通大学 | 一种旋转叶片动应变场测量方法及其系统 |
CN109902350A (zh) * | 2019-01-26 | 2019-06-18 | 北京工业大学 | 对变截面梁的截面惯性矩进行模型修正中克服模态交换的方法 |
CN110321634A (zh) * | 2019-07-03 | 2019-10-11 | 中汽研(天津)汽车工程研究院有限公司 | 一种用于多学科协同仿真的模态自动优化方法 |
CN111079320A (zh) * | 2018-10-19 | 2020-04-28 | 上海天纳克排气系统有限公司 | 排气系统的虚拟路谱仿真分析方法 |
CN111177848A (zh) * | 2019-12-26 | 2020-05-19 | 中国航空工业集团公司西安飞机设计研究所 | 一种基于有限元模型的应变理论值的获取方法和装置 |
CN113673050A (zh) * | 2021-08-11 | 2021-11-19 | 中国科学院力学研究所 | 一种基于采掘头动臂应变响应的采掘物识别方法 |
-
2016
- 2016-11-18 CN CN201611019190.7A patent/CN106529049A/zh active Pending
Non-Patent Citations (1)
Title |
---|
展铭: "基于代理模型的L型双加筋板结构动力学优化设计", 《中国优秀硕士学位论文全文数据库工程科技Ⅱ辑》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107644134A (zh) * | 2017-09-20 | 2018-01-30 | 北京安怀信科技股份有限公司 | 一种动力学验模工具 |
CN107844669A (zh) * | 2017-12-01 | 2018-03-27 | 中国直升机设计研究所 | 一种直升机全机静力有限元模型修正及验证方法 |
CN111079320A (zh) * | 2018-10-19 | 2020-04-28 | 上海天纳克排气系统有限公司 | 排气系统的虚拟路谱仿真分析方法 |
CN111079320B (zh) * | 2018-10-19 | 2023-02-21 | 上海天纳克排气系统有限公司 | 排气系统的虚拟路谱仿真分析方法 |
CN109902350A (zh) * | 2019-01-26 | 2019-06-18 | 北京工业大学 | 对变截面梁的截面惯性矩进行模型修正中克服模态交换的方法 |
CN109883389A (zh) * | 2019-03-22 | 2019-06-14 | 西安交通大学 | 一种旋转叶片动应变场测量方法及其系统 |
CN109883389B (zh) * | 2019-03-22 | 2020-08-25 | 西安交通大学 | 一种旋转叶片动应变场测量方法及其系统 |
CN110321634A (zh) * | 2019-07-03 | 2019-10-11 | 中汽研(天津)汽车工程研究院有限公司 | 一种用于多学科协同仿真的模态自动优化方法 |
CN111177848A (zh) * | 2019-12-26 | 2020-05-19 | 中国航空工业集团公司西安飞机设计研究所 | 一种基于有限元模型的应变理论值的获取方法和装置 |
CN111177848B (zh) * | 2019-12-26 | 2023-05-23 | 中国航空工业集团公司西安飞机设计研究所 | 一种基于有限元模型的应变理论值的获取方法和装置 |
CN113673050A (zh) * | 2021-08-11 | 2021-11-19 | 中国科学院力学研究所 | 一种基于采掘头动臂应变响应的采掘物识别方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106529049A (zh) | 一种用于模型修正的应变模态相关性评价方法 | |
CN106529055A (zh) | 一种基于应变模态振型相关性的模型修正方法 | |
CN101915733B (zh) | 一种基于频率变化的结构损伤程度评估方法 | |
CN108959686A (zh) | 一种基于灵敏度分析的有限元模型修正方法 | |
US9964468B1 (en) | Optimizing sensor placement for structural health monitoring | |
CN103776416B (zh) | 一种基于应变传感器的天线结构变形的间接测量方法 | |
CN106156501A (zh) | 一种基于挠度影响线的桥梁局部损伤量化方法 | |
CN108595301A (zh) | 一种基于机器学习的服务器能耗预测方法和系统 | |
CN107885908A (zh) | 一种基于振动模态的层合板等效动力学模型的建立方法 | |
JP3569899B2 (ja) | 橋梁の損傷推定システムおよびプログラム | |
CN107862170B (zh) | 一种基于动态缩聚的有限元模型修正方法 | |
CN102254068A (zh) | 大跨度桥梁抖振响应的多尺度分析方法 | |
CN103049608B (zh) | 基于约束边应变极值坐标的载荷辨识系统及方法 | |
Stutz et al. | A flexibility-based continuum damage identification approach | |
CN107016166A (zh) | 一种基于模态应变能的新型结构刚度损伤定位方法 | |
CN109839441A (zh) | 一种桥梁模态参数识别方法 | |
CN106997410A (zh) | 一种基于模态应变能的损伤发生的判断方法 | |
CN109781863A (zh) | 基于快速振动测试的结构腐蚀二阶段探测方法及其系统 | |
CN113465734A (zh) | 一种结构振动的实时估计方法 | |
CN116227045B (zh) | 一种结构试件的局部应力应变场构造方法及系统 | |
CN104199744A (zh) | 一种超级计算机应用性能稳定性判断方法及装置 | |
CN103399974B (zh) | 量化比较随机振动仿真数据与实验数据的方法 | |
CN108896625B (zh) | 一种应变场融合机敏网结构裂缝监测方法 | |
CN102998367A (zh) | 一种基于虚拟派生结构的损伤识别方法 | |
CN113627048B (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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170322 |