基于损伤力学的复合材料结构损伤识别方法技术

本发明专利技术提供了一种基于损伤力学的复合材料结构损伤识别方法，包括：建立复合材料结构的有限元模型并修正；获取未损伤的复合材料结构的位移模态，并结合修正后的有限元模型求解未损伤的复合材料结构的单元应变能；获取损伤的复合材料结构的位移模态，并结合修正后的有限元模型求解损伤的复合材料结构的单元应变能；根据未损伤的复合材料结构的单元应变能和损伤的复合材料结构的单元应变能对复合材料结构的损伤位置进行定位识别；利用损伤力学选取复合材料结构的多个损伤参量对复合材料结构的损伤位置的损伤程度进行定量化识别。本发明专利技术能够实现对复合材料结构损伤的定位和定量化识别，并能够提高识别的效率，增加对复合材料结构损伤识别的针对性。

Damage identification method of composite structure based on damage mechanics

【技术实现步骤摘要】
基于损伤力学的复合材料结构损伤识别方法
本专利技术涉及结构损伤识别
，具体涉及一种基于损伤力学的复合材料结构损伤识别方法。
技术介绍
鉴于复合材料高比强度、高比刚度的优越性，复合材料在结构中的使用率越来越大，应用结构件也从次承力结构逐步扩展到主承力结构。然而，由于复合材料在服役过程中在受到冲击或压缩载荷作用下容易产生多种损伤形式，如纤维断裂，基体断裂，分层，剪切破坏等损伤类型，因此需要检测复合材料在服役过程中的损伤情况。传统地，通常采用无损检测技术，例如超声波法、射线探伤法、红外热成像法以及声发射法，但是，这些方法是对结构属性的局部检测，需要经验支持，预先估计可能的结构损伤位置，而且不能做到实时监控和在线检测。此外，这种周期性的检测方式在增加了结构维护费用的同时，对于长时间在天空或者太空中运行的飞行器结构，难以实时检测。因此，为了解决传统无损检测技术存在的问题，研究人员研究了动力指纹法来检测复合材料的结构损伤情况。然而，目前通常采用的动力指纹法大多存在损伤识别效率较低、难以有效地对复合材料的损伤情况进行定量识别的问题。
技术实现思路
本专利技术旨在至少在一定程度上解决上述技术中的技术问题之一。为此，本专利技术的目的在于提出一种基于损伤力学的复合材料结构损伤识别方法，能够实现对复合材料结构损伤的定位和定量化识别，并能够提高对复合材料结构损伤识别的效率，增加对复合材料结构损伤识别的针对性。为达到上述目的，本专利技术实施例提出了一种基于损伤力学的复合材料结构损伤识别方法，包括以下步骤：建立复合材料结构的有限元模型并修正；获取未损伤的所述复合材料结构的位移模态，并结合修正后的所述有限元模型求解未损伤的所述复合材料结构的单元应变能；获取损伤的所述复合材料结构的位移模态，并结合修正后的所述有限元模型求解损伤的所述复合材料结构的单元应变能；根据未损伤的所述复合材料结构的单元应变能和损伤的所述复合材料结构的单元应变能对所述复合材料结构的损伤位置进行定位识别；利用损伤力学选取所述复合材料结构的多个损伤参量对所述复合材料结构的所述损伤位置的损伤程度进行定量化识别。根据本专利技术实施例的基于损伤力学的复合材料结构损伤识别方法，通过建立复合材料结构的有限元模型并修正，其次获取未损伤的复合材料结构的位移模态和未损伤的复合材料结构的位移模态，并结合修正后的有限元模型求解未损伤的复合材料结构和损伤的复合材料结构的单元应变能，然后根据未损伤的复合材料结构和损伤的复合材料结构的单元应变能对复合材料结构的损伤位置进行定位识别，最后利用损伤力学选取复合材料结构的多个损伤参量对复合材料结构损伤位置的损伤程度进行定量化识别，由此，能够实现对复合材料结构损伤的定位和定量化识别，并能够提高对复合材料结构损伤识别的效率，增加对复合材料结构损伤识别的针对性。另外，根据本专利技术上述实施例提出的基于损伤力学的复合材料结构损伤识别方法还可以具有如下附加的技术特征：根据本专利技术的一个实施例，建立复合材料结构的有限元模型并修正包括：对所述复合材料结构进行测量，并根据测量数据建立所述复合材料结构的有限元模型；对所述有限元模型进行模态分析；对未损伤的所述复合材料结构进行试验模态测量；根据所述模态分析的结果、所述试验模态测量的数值和对所述复合材料结构测量的数据构建误差函数；根据所述误差函数对所述有限元模型进行修正。进一步地，根据所述误差函数对所述有限元模型进行修正包括：设定修正阈值；若所述误差函数小于所述修正阈值，则停止对所述有限元模型的修正。进一步地，所述误差函数为：Δf＝||λe-λa(Eij,bk)||2其中，‖‖2为2范数，λe和λa分别为所述试验模态测量的各阶固有频率和所述模态分析的各阶固有频率，Eij,bk分别为待修正的所述有限模型的弹性常数和结构尺寸。进一步地，修正所述有限元模型所用的优化模型为：findEij,bks.t.ELBij≤E≤EUBijbLBk≤bk≤bUBij其中，Wε为加权系数矩阵，LB,UB分别为下边界和上边界，T为矩阵的转置符号。根据本专利技术的一个实施例，获取所述复合材料结构的多阶模态的位移模态振型，并根据所述位移模态振型求解所述复合材料结构中各阶模态的单元应变能。进一步地，所述单元应变能为：其中，i为模态阶数，j为单元编号，Uij为第i阶模态j单元的单元应变能，为第i阶位移模态振型，Kj为第j个单元的刚度矩阵。进一步地，根据损伤的所述复合材料结构和未损伤的所述复合材料结构中各阶模态的单元应变能构建所述复合材料结构的损伤位置识别指标，所述损伤位置识别指标为：其中，d为损伤，m为节点位置，n为选用的模态数，ω为加权系数。根据本专利技术的一个实施例，根据所述复合材料结构的多个损伤参量构建所述复合材料结构的所述损伤位置的损伤程度识别指标，所述损伤程度识别指标为：findd1,d2,d3s.t.0≤d1,d2,d3≤1其中，Ff为所述复合材料结构固有频率的函数，FMAC为所述复合材料结构位移模态振型的函数，C为常数，表示所述复合材料结构的固有频率和位移模态振型的权重关系，d1,d2,d3为所述损伤参量，φε为实验模态振型参数，φα为分析模态参数。进一步地，所述复合材料结构固有频率的函数为：所述复合材料结构位移模态振型的函数为：附图说明图1为本专利技术实施例的基于损伤力学的复合材料结构损伤识别方法的流程图；图2为本专利技术一个具体实施例的简支板的单元编号示意图；图3为本专利技术一个具体实施例的简支板的损伤定位指标柱状图；图4为本专利技术一个具体实施例的简支板的损伤程度识别过程中的目标函数收敛图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。图1为本专利技术实施例的基于损伤力学的复合材料结构损伤识别方法的流程图。如图1所示，本专利技术实施例的基于损伤力学的复合材料结构损伤识别方法，包括以下步骤：S1，建立复合材料结构的有限元模型并修正。具体地，可先对复合材料结构进行测量，并根据测量数据建立复合材料结构的有限元模型；其次对有限元模型进行模态分析，并对未损伤的复合材料结构进行试验模态测量；然后根据模态分析的结果和试验模态测量的数值构建误差函数；最后根据误差函数对有限元模型进行修正。其中，根据误差函数对有限元模型进行修正包括设定修正阈值，若误差函数小于修正阈值，则停止对有限元模型的修正。更具体地，可先测量复合材料结构的弹性常数和结构尺寸，并根据复合材料结构的弹性常数和结构尺寸建立复合材料结构的有限元模型；其次求解有限元模型的固有频率和位移模本文档来自技高网...

【技术保护点】
1.一种基于损伤力学的复合材料结构损伤识别方法，其特征在于，包括以下步骤：/n建立复合材料结构的有限元模型并修正；/n获取未损伤的所述复合材料结构的位移模态，并结合修正后的所述有限元模型求解未损伤的所述复合材料结构的单元应变能；/n获取损伤的所述复合材料结构的位移模态，并结合修正后的所述有限元模型求解损伤的所述复合材料结构的单元应变能；/n根据未损伤的所述复合材料结构的单元应变能和损伤的所述复合材料结构的单元应变能对所述复合材料结构的损伤位置进行定位识别；/n利用损伤力学选取所述复合材料结构的多个损伤参量对所述复合材料结构的所述损伤位置的损伤程度进行定量化识别。/n

【技术特征摘要】
1.一种基于损伤力学的复合材料结构损伤识别方法，其特征在于，包括以下步骤：
建立复合材料结构的有限元模型并修正；
获取未损伤的所述复合材料结构的位移模态，并结合修正后的所述有限元模型求解未损伤的所述复合材料结构的单元应变能；
获取损伤的所述复合材料结构的位移模态，并结合修正后的所述有限元模型求解损伤的所述复合材料结构的单元应变能；
根据未损伤的所述复合材料结构的单元应变能和损伤的所述复合材料结构的单元应变能对所述复合材料结构的损伤位置进行定位识别；
利用损伤力学选取所述复合材料结构的多个损伤参量对所述复合材料结构的所述损伤位置的损伤程度进行定量化识别。


2.根据权利要求1所述的基于损伤力学的复合材料结构损伤识别方法，其特征在于，建立复合材料结构的有限元模型并修正包括：
对所述复合材料结构进行测量，并根据测量数据建立所述复合材料结构的有限元模型；
对所述有限元模型进行模态分析；
对未损伤的所述复合材料结构进行试验模态测量；
根据所述模态分析的结果、所述试验模态测量的数值和对所述复合材料结构测量的数据构建误差函数；
根据所述误差函数对所述有限元模型进行修正。


3.根据权利要求2所述的基于损伤力学的复合材料结构损伤识别方法，其特征在于，根据所述误差函数对所述有限元模型进行修正包括：
设定修正阈值；
若所述误差函数小于所述修正阈值，则停止对所述有限元模型的修正。


4.根据权利要求3所述的基于损伤力学的复合材料结构损伤识别方法，其特征在于，所述误差函数为：
Δf＝||λe-λa(Eij,bk)||2
其中，‖‖2为2范数，λe和λa分别为所述试验模态测量的各阶固有频率和所述模态分析的各阶固有频率，Eij,bk分别为待修正的所述有限模型的弹性常数和结构尺寸。


5.根据权利要求4所述的基于损伤力学的复合材料结构损伤识别方法，其特征在于，其中，修正所述有限元模型所用的优化模型为：
findEij,bk
min...

【专利技术属性】
技术研发人员：石庆贺，胡可军，韩文钦，王璐，
申请(专利权)人：江苏理工学院，
类型：发明
国别省市：江苏;32