CN103049608B - 基于约束边应变极值坐标的载荷辨识系统及方法 - Google Patents
基于约束边应变极值坐标的载荷辨识系统及方法 Download PDFInfo
- Publication number
- CN103049608B CN103049608B CN201210551643.6A CN201210551643A CN103049608B CN 103049608 B CN103049608 B CN 103049608B CN 201210551643 A CN201210551643 A CN 201210551643A CN 103049608 B CN103049608 B CN 103049608B
- Authority
- CN
- China
- Prior art keywords
- load
- extreme coordinates
- strain
- ing
- point
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000013507 mapping Methods 0.000 claims abstract description 12
- 230000004044 response Effects 0.000 claims description 12
- 238000012850 discrimination method Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 3
- 238000004088 simulation Methods 0.000 claims description 2
- 239000000463 material Substances 0.000 abstract description 4
- 230000008901 benefit Effects 0.000 abstract description 3
- 238000012544 monitoring process Methods 0.000 abstract description 3
- 229910000838 Al alloy Inorganic materials 0.000 description 7
- 229920000049 Carbon (fiber) Polymers 0.000 description 7
- 239000004917 carbon fiber Substances 0.000 description 7
- 239000002131 composite material Substances 0.000 description 7
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 7
- 238000010276 construction Methods 0.000 description 5
- 230000003068 static effect Effects 0.000 description 5
- 238000007906 compression Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000009977 dual effect Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001235 sensitizing effect Effects 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
本发明公开了基于约束边应变极值坐标的载荷辨识系统及方法,属于板壳结构监测技术领域。所述系统包括:用于测量板壳结构应变值的传感器、数据处理器、计算机。所述方法首先建立加载点坐标与对应应变极值坐标的映射关系;接着沿板壳结构边缘均匀布局传感器,测量待判位加载点对应应变极值坐标;然后需找距离待判位加载点对应应变极值坐标最小的样本加载节点;最后根据反距离加权差值法计算得到待判位加载点坐标。本发明所述系统具有结构简单,传感器布局简单的优点。本发明所述方法通过采集应变极值坐标,简化了加载点判位方法;对不同材料参数的结构具有通用性。
Description
技术领域
本发明公开了基于约束边应变极值坐标的载荷辨识系统及方法,属于板壳结构监测技术领域。
背景技术
板壳类结构是飞机中一种重要的典型结构,南京航空航天大学陆观等人的研究成果显示,传感网络信号的有效性和传感网络覆盖位置的正确性依赖于传感网络敏感区的分布,首先需要对整体的结构状态特征进行分析,采用与结构相匹配的算法,优化得到传感器的分布。此外,针对传感器分布常用的构建方法包括有效独立、Guvan模型缩减法。有效独立法从所有可能测点出发,利用复模态矩阵的幂等型,计算有效独立向量,按照目标模态矩阵独立性排序,删除对其秩贡献最小的自由度,从而优化Fisher信息阵而使感兴趣的模态向量尽可能保持线性无关。Guyan模型缩减法也是一种常用的测点选择方法,通过刚度(静力缩减)或质量(动力缩减)子矩阵构成的转换矩阵,可以把那些对模态反应起主要作用的自由度保留下来作为测点的位置。以上传感器分布方法都要对结构整体先进行分析研究,在此基础上优化得到传感器的分布,确定实验系统。
飞机在工作环境中由于承受多各种形式载荷共同作用,板壳类结构会产生不可以预见的损伤。由于飞行任务模式中强过载和加载形式日趋复杂化,容易造成结构强度失效,如复合材料板会造成板结构内部纤维断裂,金属结构会造成局部屈服变形,使得飞机典型承载部件如机体、机翼力学性能大幅下降,直接危及飞行安全。因此必须及时获取载荷加载信息,从而为结构健康状态准确评估和视情维修提供可靠依据。
结构损伤识别的方法主要有:模式匹配法、敏感因子法、模型修正法和参数识别法。由于基于结构静力响应的损伤识别具有稳定性高、数据处理简单的优点,使得目前针对结构静力响应特性的研究逐渐增加。
传统板壳结构载荷加载点判位方法对传感器布置缺少统一的布置规则,因此采用不同布置方法时,传感器测量到的应变/应力数值也不相同。,因此需要以从实际目标结构获取基于样本加载节点的先验结构响应数据作为静载判位的参考基准,根据先验数据布局传感器,这在一定程度上使得判位方法只适用于特定结构,对不同形状的结构需要重新设计传感器测点位置,并再次获取先验数据。随着实际监测场合和承载形式日趋复杂多变,不仅难以直接获取这些实测数据,而且还会使得相关传感器网络配置形式难以确定。
发明内容
本发明所要解决的技术问题是针对上述背景技术的不足,提供了基于约束边应变极值坐标的载荷辨识系统及方法。
本发明为实现上述发明目的采用如下技术方案:
基于约束边应变极值坐标的载荷辨识方法,建立包括:用于测量板壳结构应变值的传感器、数据处理器、计算机的系统,所述数据处理器输入端与用于测量板壳结构应变值的传感器连接,输出端与计算机连接,
载荷辨识方法具体包括如下步骤:
步骤1,以板壳结构上的一点为原点建立二维坐标系,建立加载点坐标与对应应变极值坐标的映射关系:
步骤A,在板壳结构上取n个样本加载节点,在第k个样本加载节点上加载载荷,利用有限元仿真方法得到第k个样本加载节点对应的应变极值坐标;其中,n为自然数,k为小于n的自然数,
步骤B,构建由每个样本加载节点坐标及其相应应变极值坐标组成的矩阵,并建立样本加载节点坐标及其相应应变极值坐标间映射关系;
步骤2,在所述板壳结构上布局用于测量板壳结构应变值的传感器:沿板壳结构约束边均匀布置传感器,所述传感器粘贴方向与约束边垂直;
步骤3,对加载点加载载荷,根据按照步骤2所述方式布局的传感器的响应数据,对步骤1确定的映射关系反演得到加载点目标。
所述基于约束边应变极值坐标的载荷辨识方法,步骤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的自然数:
反距离加权插值法公式是:
式(2)中,(xi,yi)为四个样本加载节点坐标,为4个样本加载节点中一个样本节点分。
距离越小表示待判位加载点与已知样本加载节点距离越近,选取距离最近的4个节点,按距离的反比分配权重,根据公式(2)计算实际加载坐标(x0,y0)。
图7是各向同性铝合金四边固支板不同泊松比μ对应的x方向应变坐标。图8是各向同性铝合金四边固支板不同弹性模量E对应的x方向应变坐标。图9是各向异性复合材料碳纤维四边固支板不同泊松比PRYZ对应的x方向应变坐标。图10是各向异性复合材料碳纤维四边固支板不同弹性模量MOEX对应的x方向应变坐标。从图中可以看出在结构参数(泊松比、弹性模量)改变较大的时候,测试的板壳结构各点应变值大小会发生改变,但是应变极值对应的坐标基本没有变化。可见,利用本发明所述方法判位同形状不同材料的板壳结构载荷加载点时,应变坐标与加载坐标之间的映射f保持不变。
综上所述,本发明通过采集应变极值坐标,简化了加载点判位方法;对不同材料参数的结构具有通用性。
Claims (2)
1.基于约束边应变极值坐标的载荷辨识方法,建立包括:用于测量板壳结构应变值的传感器、数据处理器、计算机的系统,所述数据处理器输入端与用于测量板壳结构应变值的传感器连接,输出端与计算机连接,
其特征在于,载荷辨识方法具体包括如下步骤:
步骤1,以板壳结构上的一点为原点建立二维坐标系,建立加载点坐标与对应应变极值坐标的映射关系:
步骤A,在板壳结构上取n个样本加载节点,在第k个样本加载节点上加载载荷,利用有限元仿真方法得到第k个样本加载节点对应的应变极值坐标;其中,n为自然数,k为小于n的自然数,
步骤B,构建由每个样本加载节点坐标及其相应应变极值坐标组成的矩阵,并建立样本加载节点坐标及其相应应变极值坐标间映射关系;
步骤2,在所述板壳结构上布局用于测量板壳结构应变值的传感器:沿板壳结构约束边均匀布置传感器,所述传感器粘贴方向与约束边垂直;
步骤3,对加载点加载载荷,根据按照步骤2所述方式布局的传感器的响应数据,对步骤1确定的映射关系反演得到加载点目标。
2.根据权利要求1所述的基于约束边应变极值坐标的载荷辨识方法,其特征在于,所述步骤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 CN103049608A (zh) | 2013-04-17 |
CN103049608B true 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) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104019739B (zh) * | 2014-05-08 | 2016-08-17 | 南京航空航天大学 | 基于分布式光纤光栅传感网络的板结构冲击载荷定位方法 |
CN104776966B (zh) * | 2015-04-01 | 2017-05-10 | 南京航空航天大学 | 基于分形理论的板结构冲击监测方法 |
CN106592357B (zh) * | 2016-12-14 | 2018-10-23 | 中国人民解放军空军工程大学 | 一种机场道面弯道过渡面高程设计方法 |
CN107991060B (zh) * | 2017-11-20 | 2019-08-20 | 南京航空航天大学 | 基于自适应和迭代算法的载荷分布式光纤辨识方法 |
CN112948968B (zh) * | 2021-03-01 | 2024-03-29 | 南京航空航天大学 | 基于应变响应阈值计算的光纤传感器布局方法 |
CN113705034B (zh) * | 2021-07-19 | 2023-12-19 | 东风汽车集团股份有限公司 | 一种仿真结果处理方法和装置 |
CN116150551B (zh) * | 2023-04-19 | 2023-09-29 | 南京航空航天大学 | 基于等效附加载荷应变响应特征计算的结构损伤辨识方法 |
CN116153842B (zh) * | 2023-04-20 | 2023-06-23 | 涌淳半导体(无锡)有限公司 | 用于测温的吸附装置及其测温方法、电子装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102592017A (zh) * | 2011-12-31 | 2012-07-18 | 北京工业大学 | 一种双面锁紧刀柄/主轴联接性能仿真优化方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1337942B1 (en) * | 2000-11-17 | 2016-10-12 | Battelle Memorial Institute | Method and system for structural stress analysis |
-
2012
- 2012-12-18 CN CN201210551643.6A patent/CN103049608B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102592017A (zh) * | 2011-12-31 | 2012-07-18 | 北京工业大学 | 一种双面锁紧刀柄/主轴联接性能仿真优化方法 |
Non-Patent Citations (3)
Title |
---|
严刚 等.加筋复合材料结构的冲击载荷识别.《航空学报》.2008,第29卷(第5期), * |
周伟 等.基于光纤智能夹层的车架有限元分析与试验研究.《传感器与微系统》.2008,第27卷(第12期), * |
基于光纤传感网络的可变体机翼应变场数值模拟及实验验证;徐海伟 等;《航空学报》;20111025;第32卷(第10期);第3页左栏第1段至第4页右栏最后1段,图1至图4 * |
Also Published As
Publication number | Publication date |
---|---|
CN103049608A (zh) | 2013-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103049608B (zh) | 基于约束边应变极值坐标的载荷辨识系统及方法 | |
Dahak et al. | Damage detection in beam through change in measured frequency and undamaged curvature mode shape | |
Gravenkamp et al. | A numerical approach for the computation of dispersion relations for plate structures using the scaled boundary finite element method | |
Xiang et al. | Identification of a crack in a beam based on the finite element method of a B-spline wavelet on the interval | |
CN108875178B (zh) | 用于减小结构模态识别不确定性的传感器布置方法 | |
CN106289947B (zh) | 一种轻质高强梁结构损伤识别方法 | |
Khiem et al. | A method for crack identification in functionally graded Timoshenko beam | |
CN106529049A (zh) | 一种用于模型修正的应变模态相关性评价方法 | |
CN106528945A (zh) | 基于加载夹角归一化的板结构载荷位置与大小辨识方法 | |
CN110375690A (zh) | 一种旋转叶片非接触式位移场测量方法及其系统 | |
Lecieux et al. | Monitoring of a reinforced concrete wharf using structural health monitoring system and material testing | |
Zhou et al. | The magneto-electro-elastic multi-physics coupling element free Galerkin method for smart structures in statics and dynamics problems | |
Kumar et al. | Development and characterization of a modified ring shaped force transducer | |
Ştefănescu | Strain sensitivity as selection criterion for elastic elements of force transducers: A brief review | |
CN115659737A (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 | |
CN114021403A (zh) | 基于应变模态的承力结构件损伤识别方法及系统 | |
CN103940453B (zh) | 一种提高传感器测量精度方法 | |
CN113639941A (zh) | 一种测试桥梁应变柔度矩阵的方法 | |
Li et al. | Health monitoring of flexural steel structures based on distributed fibre optic sensors | |
Sim | Estimation of flexibility matrix of beam structures using multisensor fusion | |
Kaul | Crack diagnostics in beams using wavelets, kurtosis and skewness | |
CN203224874U (zh) | 基于约束边应变极值坐标的载荷辨识系统 | |
CN102270249B (zh) | 识别零部件特征频率的方法 | |
Yan | An effective numerical approach for multiple void-crack interaction |
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 |