CN103049608A - 基于约束边应变极值坐标的载荷辨识系统及方法 - Google Patents
基于约束边应变极值坐标的载荷辨识系统及方法 Download PDFInfo
- Publication number
- CN103049608A CN103049608A CN2012105516436A CN201210551643A CN103049608A CN 103049608 A CN103049608 A CN 103049608A CN 2012105516436 A CN2012105516436 A CN 2012105516436A CN 201210551643 A CN201210551643 A CN 201210551643A CN 103049608 A CN103049608 A CN 103049608A
- Authority
- CN
- China
- Prior art keywords
- load
- coordinate
- extreme value
- strain
- ing
- 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 Strength Of Materials By Application Of Mechanical Stress (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
本发明公开了基于约束边应变极值坐标的载荷辨识系统及方法,属于板壳结构监测技术领域。所述系统包括:用于测量板壳结构应变值的传感器、数据处理器、计算机。所述方法首先建立加载点坐标与对应应变极值坐标的映射关系;接着沿板壳结构边缘均匀布局传感器,测量待判位加载点对应应变极值坐标;然后需找距离待判位加载点对应应变极值坐标最小的样本加载节点;最后根据反距离加权差值法计算得到待判位加载点坐标。本发明所述系统具有结构简单,传感器布局简单的优点。本发明所述方法通过采集应变极值坐标,简化了加载点判位方法;对不同材料参数的结构具有通用性。
Description
技术领域
本发明公开了基于约束边应变极值坐标的载荷辨识系统及方法,属于板壳结构监测技术领域。
背景技术
板壳类结构是飞机中一种重要的典型结构,南京航空航天大学陆观等人的研究成果显示,传感网络信号的有效性和传感网络覆盖位置的正确性依赖于传感网络敏感区的分布,首先需要对整体的结构状态特征进行分析,采用与结构相匹配的算法,优化得到传感器的分布。此外,针对传感器分布常用的构建方法包括有效独立、Guvan模型缩减法。有效独立法从所有可能测点出发,利用复模态矩阵的幂等型,计算有效独立向量,按照目标模态矩阵独立性排序,删除对其秩贡献最小的自由度,从而优化Fisher信息阵而使感兴趣的模态向量尽可能保持线性无关。Guyan模型缩减法也是一种常用的测点选择方法,通过刚度(静力缩减)或质量(动力缩减)子矩阵构成的转换矩阵,可以把那些对模态反应起主要作用的自由度保留下来作为测点的位置。以上传感器分布方法都要对结构整体先进行分析研究,在此基础上优化得到传感器的分布,确定实验系统。
飞机在工作环境中由于承受多各种形式载荷共同作用,板壳类结构会产生不可以预见的损伤。由于飞行任务模式中强过载和加载形式日趋复杂化,容易造成结构强度失效,如复合材料板会造成板结构内部纤维断裂,金属结构会造成局部屈服变形,使得飞机典型承载部件如机体、机翼力学性能大幅下降,直接危及飞行安全。因此必须及时获取载荷加载信息,从而为结构健康状态准确评估和视情维修提供可靠依据。
结构损伤识别的方法主要有:模式匹配法、敏感因子法、模型修正法和参数识别法。由于基于结构静力响应的损伤识别具有稳定性高、数据处理简单的优点,使得目前针对结构静力响应特性的研究逐渐增加。
传统板壳结构载荷加载点判位方法对传感器布置缺少统一的布置规则,因此采用不同布置方法时,传感器测量到的应变/应力数值也不相同。,因此需要以从实际目标结构获取基于样本加载节点的先验结构响应数据作为静载判位的参考基准,根据先验数据布局传感器,这在一定程度上使得判位方法只适用于特定结构,对不同形状的结构需要重新设计传感器测点位置,并再次获取先验数据。随着实际监测场合和承载形式日趋复杂多变,不仅难以直接获取这些实测数据,而且还会使得相关传感器网络配置形式难以确定。
发明内容
本发明所要解决的技术问题是针对上述背景技术的不足,提供了基于约束边应变极值坐标的载荷辨识系统及方法。
本发明为实现上述发明目的采用如下技术方案:
基于约束边应变极值坐标的载荷辨识系统,包括:用于测量板壳结构应变值的传感器、数据处理器、计算机;其中,所述数据处理器输入端与用于测量板桥结构应变值的传感器,输出端与计算机连接。
所述基于约束边应变极值坐标的载荷辨识系统中,用于测量板壳结构应变值的传感器为光纤布拉格光栅传感器,所述数据处理器为光栅解调仪。
所述基于约束边应变极值坐标的载荷辨识系统中,计算机中安装有MATLAB软件,通过MATLAB的定位程序对数据处理器的输出值进行处理。
基于约束边应变极值坐标的载荷辨识方法,使用权利要求1或2或3所述的利用约束边应变极值坐标判定载荷加载点的系统实现,具体包括如下步骤:
步骤1,以板壳结构上的一点为原点建立二维坐标系,建立加载点坐标与对应应变极值坐标的映射关系;
步骤2,在所述板壳结构上布局用于测量板壳结构应变值的传感器:沿板壳结构约束边均匀布置传感器,所述传感器粘贴方向与约束边垂直;
步骤3,对加载点加载载荷,根据按照步骤2所述方式布局的传感器的响应数据,对步骤1确定的映射关系反演得到加载点目标。
所述基于约束边应变极值坐标的载荷辨识方法中,步骤1中建立载荷加载点坐标与应变极值坐标的映射关系的具体方法如下:
步骤A,在板壳结构上取n个样本加载节点,在第k个样本加载节点上加载载荷,利用有限元仿真方法得到第k个样本加载节点对应的应变极值坐标;其中,n为自然数,k为小于n的自然数;
步骤B,构建由每个样本加载节点坐标及其相应应变极值坐标组成的矩阵,并建立样本加载节点坐标及其相应应变极值坐标间映射关系。
所述基于约束边应变极值坐标的载荷辨识方法中,步骤3的具体实施方式如下:
步骤a,对待判位加载点加载载荷,记录传感器的响应数据,根据传感器响应数据得到待判位加载点的应变极值坐标;
步骤b,寻找距离待判位加载点对应应变极值坐标最近的i个样本加载节点对应应变极值坐标,计算待判位加载点对应应变极值坐标与i个样本加载节点对应应变极值坐标的距离di,i为小于n的自然数;
步骤c,利用步骤a所述的i个样本加载节点对应应变极值坐标,反演步骤1得到的影射关系,得到i个样本加载节点加载坐标(xi,yi);利用反距离加权插值法确定待判位加载点坐标(x0,y0)。
本发明采用上述技术方案,具有以下有益效果:本发明所述系统具有结构简单,传感器布局简单的优点。本发明所述方法通过采集应变极值坐标,简化了加载点判位方法;对不同材料参数的结构具有通用性。
附图说明
图1是板结构应变坐标与加载坐标示意图。
图2是本发明判位流程图。
图3是本发明方法静载识别系统构成示意图。
图4是反距离加权法示意图。
图5是各向同性铝合金四边固支板加载后约束边y方向应变图。
图6是各向异性复合材料碳纤维四边固支板加载后约束边y方向应变图。
图7是各向同性铝合金四边固支板不同泊松比对应的x方向应变坐标。
图8是各向同性铝合金四边固支板不同弹性模量对应的x方向应变坐标。
图9是各向异性复合材料碳纤维四边固支板不同泊松比对应的x方向应变坐标。
图10是各向异性复合材料碳纤维四边固支板不同弹性模量对应的x方向应变坐标。
具体实施方式
下面结合附图对发明的技术方案进行详细说明:
基于约束边应变极值坐标的载荷辨识系统及方法如图3所示,包括FBG(Fiber Bragg Grating,光纤布拉格光栅)传感器、Si425光栅解调仪、计算机。Si425光栅解调仪输入端与用于测量板桥结构应变值的传感器,输出端与计算机连接。Si425光栅解调仪将FBG(Fiber Bragg Grating,光纤布拉格光栅)传感器测得的应变值转换为板壳结构应变极值坐标。根据光栅解调仪分析得到的应力极值坐标,在MATLAB编写计算加载点位置坐标的程序,利用程序计算得到加载点位置坐标。
如图1所示,载荷实际加载点为图中(x,y)处,称为加载坐标。a为加载后,板AB边负应变极值处对应的x方向坐标,但实际测量时,传感器无法布置在约束边上,因此以图中EF线上应变数值代替AB边应变数值作为实际测量数据。同理,b表示AD边负应变极值处的y方向坐标。由此可组合构成新的坐标(a,b),该坐标表示加载点(x,y)所对应的应变极值坐标。
如图2所示,对各向同性铝合金四边固支板(试件长×宽×厚为1200×1200×2mm)以及各向异性复合材料碳纤维四边固支板(试件长×宽×厚为540×540×2.54mm)辨识载荷点,包括如下步骤:
步骤1,以板壳结构上的A点为原点建立二维坐标系,建立加载点坐标与对应应变极值坐标的映射关系f;建立载荷加载点坐标与应变极值坐标的映射关系f的具体方法如下:
步骤A,在板壳结构上取n个样本加载节点,在第k个样本加载节点上加载载荷,利用有限元仿真方法得到第k个样本加载节点对应的应变极值坐标(ak,bk);其中,n为自然数,k为小于n的自然数;
步骤B,构建由每个样本加载节点坐标及其相应应变极值坐标组成的矩阵,并通过矩阵运算建立样本加载节点坐标及其相应应变极值坐标间映射关系f:
(ak,bk)=f(xk,yk)。
步骤2,如图3所示,在所述板壳结构上布局FBG(Fiber Bragg Grating,光纤布拉格光栅)传感器:在x轴对应板壳结构约束边(AB边)上、y轴对应板桥结构约束边(AD边)上均匀布置4个FBG传感器,FBG传感器粘贴方向与约束边垂直(由于FBG传感器是有方向性的,且在仿真分析中与实际测量中均得出FBG传感器粘贴方向应与约束边垂直)。x轴对应板壳结构约束边(AB边)上的4个FBG传感器串联后将数据传输至Si425光栅解调仪,y轴对应板桥结构约束边(AD边)上的4个FBG传感器串联后将数据传输至Si425光栅解调仪,计算机根据Si425光栅解调仪的处理结果得到加载点对应应变力极值坐标。此外,其它用于测量板壳结构应变值的传感器都适用于本发明。图5所示各向同性铝合金四边固支板AB边应变分布与x轴坐标的关系,图6是碳纤维板AB应变分布与x轴坐标的关系。边可得应变分布存在单峰性。
步骤3,对加载点加载载荷,根据按照步骤2所述方式布局的FBG传感器的响应数据,对步骤1确定的映射关系反演得到加载点目标,具体实施方式如下:
步骤a,对待判位加载点加载载荷,记录传感器的响应数据,根据传感器响应数据得到待判位加载点的应变极值坐标(a0,b0):根据x轴对应板壳结构约束边(AB边)上、y轴对应板桥结构约束边(AD边)上的四个传感器数据使用三次多项式插值,通过拟合的方法计算出待判位加载点的应变极值坐标(a0,b0)。
步骤b,寻找距离待判位加载点对应应变极值坐标(a0,b0)最近的4个样本加载节点对应应变极值坐标(ai,bi),计算待判位加载点对应应变极值坐标(a0,b0)与4个样本加载节点对应应变极值坐标(ai,bi)的距离di,i=1,2,3,4;
式(1)中,di为第i个样本加载节点的应变坐标(ai,bi)与待判位加载点的应变极值坐标(a0,b0)之间的距离。
步骤c,利用步骤a所述的4个样本加载节点对应应变极值坐标(ai,bi),反演步骤1得到的影射关系,得到4个样本加载节点加载坐标(xi,yi);利用反距离加权插值法确定待判位加载点坐标(x0,y0);其中,i为小于n的自然数:
反距离加权插值法公式是:
距离越小表示待判位加载点与已知样本加载节点距离越近,选取距离最近的4个节点,按距离的反比分配权重,根据公式(2)计算实际加载坐标(x0,y0)。
图7是各向同性铝合金四边固支板不同泊松比μ对应的x方向应变坐标。图8是各向同性铝合金四边固支板不同弹性模量E对应的x方向应变坐标。图9是各向异性复合材料碳纤维四边固支板不同泊松比PRYZ对应的x方向应变坐标。图10是各向异性复合材料碳纤维四边固支板不同弹性模量MOEX对应的x方向应变坐标。从图中可以看出在结构参数(泊松比、弹性模量)改变较大的时候,测试的板壳结构各点应变值大小会发生改变,但是应变极值对应的坐标基本没有变化。可见,利用本发明所述方法判位同形状不同材料的板壳结构载荷加载点时,应变坐标与加载坐标之间的映射f保持不变。
综上所述,本发明通过采集应变极值坐标,简化了加载点判位方法;对不同材料参数的结构具有通用性。
Claims (6)
1.基于约束边应变极值坐标的载荷辨识系统,其特征在于包括:用于测量板壳结构应变值的传感器、数据处理器、计算机;其中,所述数据处理器输入端与用于测量板桥结构应变值的传感器,输出端与计算机连接。
2.根据权利要求1所述的基于约束边应变极值坐标的载荷辨识系统,其特征在于所述用于测量板壳结构应变值的传感器为光纤布拉格光栅传感器,所述数据处理器为光栅解调仪。
3.根据权利要求1或2所述的基于约束边应变极值坐标的载荷辨识系统,其特征在于计算机中安装有MATLAB软件,通过MATLAB的定位程序对数据处理器的输出值进行处理。
4.基于约束边应变极值坐标的载荷辨识方法,其特征在于使用权利要求1或2或3所述的利用约束边应变极值坐标判定载荷加载点的系统实现,具体包括如下步骤:
步骤1,以板壳结构上的一点为原点建立二维坐标系,建立加载点坐标与对应应变极值坐标的映射关系;
步骤2,在所述板壳结构上布局用于测量板壳结构应变值的传感器:沿板壳结构约束边均匀布置传感器,所述传感器粘贴方向与约束边垂直;
步骤3,对加载点加载载荷,根据按照步骤2所述方式布局的传感器的响应数据,对步骤1确定的映射关系反演得到加载点目标。
5.根据权利要求4所述的基于约束边应变极值坐标的载荷辨识方法,其特征在于所述步骤1中建立载荷加载点坐标与应变极值坐标的映射关系的具体方法如下:
步骤A,在板壳结构上取n个样本加载节点,在第k个样本加载节点上加载载荷,利用有限元仿真方法得到第k个样本加载节点对应的应变极值坐标;其中,n为自然数,k为小于n的自然数;
步骤B,构建由每个样本加载节点坐标及其相应应变极值坐标组成的矩阵,并建立样本加载节点坐标及其相应应变极值坐标间映射关系。
6.根据权利要求4或5所述的基于约束边应变极值坐标的载荷辨识方法,其特征在于所述步骤3的具体实施方式如下:
步骤a,对待判位加载点加载载荷,记录传感器的响应数据,根据传感器响应数据得到待判位加载点的应变极值坐标;
步骤b,寻找距离待判位加载点对应应变极值坐标最近的i个样本加载节点对应应变极值坐标,计算待判位加载点对应应变极值坐标与i个样本加载节点对应应变极值坐标的距离di,i为小于n的自然数;
步骤c,利用步骤a所述的i个样本加载节点对应应变极值坐标,反演步骤1得到的影射关系,得到i个样本加载节点加载坐标(xi,yi);利用反距离加权插值法确定待判位加载点坐标(x0,y0)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210551643.6A CN103049608B (zh) | 2012-12-18 | 2012-12-18 | 基于约束边应变极值坐标的载荷辨识系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210551643.6A CN103049608B (zh) | 2012-12-18 | 2012-12-18 | 基于约束边应变极值坐标的载荷辨识系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103049608A true CN103049608A (zh) | 2013-04-17 |
CN103049608B CN103049608B (zh) | 2015-10-28 |
Family
ID=48062245
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210551643.6A Active CN103049608B (zh) | 2012-12-18 | 2012-12-18 | 基于约束边应变极值坐标的载荷辨识系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103049608B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104019739A (zh) * | 2014-05-08 | 2014-09-03 | 南京航空航天大学 | 基于分布式光纤光栅传感网络的板结构冲击载荷定位方法 |
CN104776966A (zh) * | 2015-04-01 | 2015-07-15 | 南京航空航天大学 | 基于分形理论的板结构冲击监测方法 |
CN106592357A (zh) * | 2016-12-14 | 2017-04-26 | 中国人民解放军空军工程大学 | 一种机场道面弯道过渡面高程设计方法 |
CN107991060A (zh) * | 2017-11-20 | 2018-05-04 | 南京航空航天大学 | 基于自适应和迭代算法的载荷分布式光纤辨识方法 |
CN112948968A (zh) * | 2021-03-01 | 2021-06-11 | 南京航空航天大学 | 基于应变响应阈值计算的光纤传感器布局方法 |
CN113705034A (zh) * | 2021-07-19 | 2021-11-26 | 东风汽车集团股份有限公司 | 一种仿真结果处理方法和装置 |
CN116153842A (zh) * | 2023-04-20 | 2023-05-23 | 涌淳半导体(无锡)有限公司 | 用于测温的吸附装置及其测温方法、电子装置 |
CN116150551A (zh) * | 2023-04-19 | 2023-05-23 | 南京航空航天大学 | 基于等效附加载荷应变响应特征计算的结构损伤辨识方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020112548A1 (en) * | 2000-11-17 | 2002-08-22 | Pingsha Dong | Structural stress analysis |
CN102592017A (zh) * | 2011-12-31 | 2012-07-18 | 北京工业大学 | 一种双面锁紧刀柄/主轴联接性能仿真优化方法 |
-
2012
- 2012-12-18 CN CN201210551643.6A patent/CN103049608B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020112548A1 (en) * | 2000-11-17 | 2002-08-22 | Pingsha Dong | Structural stress analysis |
CN102592017A (zh) * | 2011-12-31 | 2012-07-18 | 北京工业大学 | 一种双面锁紧刀柄/主轴联接性能仿真优化方法 |
Non-Patent Citations (3)
Title |
---|
严刚 等: "加筋复合材料结构的冲击载荷识别", 《航空学报》 * |
周伟 等: "基于光纤智能夹层的车架有限元分析与试验研究", 《传感器与微系统》 * |
徐海伟 等: "基于光纤传感网络的可变体机翼应变场数值模拟及实验验证", 《航空学报》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104019739B (zh) * | 2014-05-08 | 2016-08-17 | 南京航空航天大学 | 基于分布式光纤光栅传感网络的板结构冲击载荷定位方法 |
CN104019739A (zh) * | 2014-05-08 | 2014-09-03 | 南京航空航天大学 | 基于分布式光纤光栅传感网络的板结构冲击载荷定位方法 |
CN104776966B (zh) * | 2015-04-01 | 2017-05-10 | 南京航空航天大学 | 基于分形理论的板结构冲击监测方法 |
CN104776966A (zh) * | 2015-04-01 | 2015-07-15 | 南京航空航天大学 | 基于分形理论的板结构冲击监测方法 |
CN106592357B (zh) * | 2016-12-14 | 2018-10-23 | 中国人民解放军空军工程大学 | 一种机场道面弯道过渡面高程设计方法 |
CN106592357A (zh) * | 2016-12-14 | 2017-04-26 | 中国人民解放军空军工程大学 | 一种机场道面弯道过渡面高程设计方法 |
CN107991060A (zh) * | 2017-11-20 | 2018-05-04 | 南京航空航天大学 | 基于自适应和迭代算法的载荷分布式光纤辨识方法 |
CN112948968A (zh) * | 2021-03-01 | 2021-06-11 | 南京航空航天大学 | 基于应变响应阈值计算的光纤传感器布局方法 |
CN112948968B (zh) * | 2021-03-01 | 2024-03-29 | 南京航空航天大学 | 基于应变响应阈值计算的光纤传感器布局方法 |
CN113705034A (zh) * | 2021-07-19 | 2021-11-26 | 东风汽车集团股份有限公司 | 一种仿真结果处理方法和装置 |
CN113705034B (zh) * | 2021-07-19 | 2023-12-19 | 东风汽车集团股份有限公司 | 一种仿真结果处理方法和装置 |
CN116150551A (zh) * | 2023-04-19 | 2023-05-23 | 南京航空航天大学 | 基于等效附加载荷应变响应特征计算的结构损伤辨识方法 |
CN116150551B (zh) * | 2023-04-19 | 2023-09-29 | 南京航空航天大学 | 基于等效附加载荷应变响应特征计算的结构损伤辨识方法 |
CN116153842A (zh) * | 2023-04-20 | 2023-05-23 | 涌淳半导体(无锡)有限公司 | 用于测温的吸附装置及其测温方法、电子装置 |
CN116153842B (zh) * | 2023-04-20 | 2023-06-23 | 涌淳半导体(无锡)有限公司 | 用于测温的吸附装置及其测温方法、电子装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103049608B (zh) | 2015-10-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103049608A (zh) | 基于约束边应变极值坐标的载荷辨识系统及方法 | |
Kim et al. | Pseudo-static approach for damage identification of bridges based on coupling vibration with a moving vehicle | |
Yang et al. | A damage identification approach for plate structures based on frequency measurements | |
Law et al. | Time-varying wind load identification from structural responses | |
CN109030202A (zh) | 一种快速确定岩石类脆性材料离散元模型参数的方法 | |
Xiang et al. | Identification of a crack in a beam based on the finite element method of a B-spline wavelet on the interval | |
Holl et al. | 3D multiscale crack propagation using the XFEM applied to a gas turbine blade | |
CN108875178B (zh) | 用于减小结构模态识别不确定性的传感器布置方法 | |
CN102254068A (zh) | 大跨度桥梁抖振响应的多尺度分析方法 | |
Maina et al. | Linear elastic analysis of pavement structure under non-circular loading | |
CN106354955A (zh) | 一种基于磨机振动模态参数的滑动轴承刚度识别方法 | |
CN103884593A (zh) | 基于少量静态测量数据的结构损伤定位与程度识别方法 | |
CN110807285B (zh) | 一种工程运输车车架疲劳分析方法、系统及装置 | |
CN102297784A (zh) | 一种建筑结构物抗风安全性能的评价方法 | |
CN104298823A (zh) | 一种高、低温天平分析方法及分析系统 | |
CN115659737A (zh) | 一种基于曲率模态和离散小波变换的梁构件损伤识别方法 | |
Pan et al. | Sparse regularization-based damage detection in a bridge subjected to unknown moving forces | |
CN110728094A (zh) | 一种新型整体结构拟静力试验方法 | |
Lan et al. | On the uniqueness and sensitivity issues in determining the elastic and plastic properties of power-law hardening materials through sharp and spherical indentation | |
CN103940453B (zh) | 一种提高传感器测量精度方法 | |
Li et al. | Structure design of six-component strain-gauge-based transducer for minimum cross-interference via hybrid optimization methods | |
CN203224874U (zh) | 基于约束边应变极值坐标的载荷辨识系统 | |
CN102270249B (zh) | 识别零部件特征频率的方法 | |
CN105510891A (zh) | 针对卫星合成孔径雷达的在轨模态辨识方法 | |
Sim | Estimation of flexibility matrix of beam structures using multisensor fusion |
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 |