CN115371627A - 一种轻量化舱体结构应变场重构方法及*** - Google Patents

一种轻量化舱体结构应变场重构方法及*** Download PDF

Info

Publication number
CN115371627A
CN115371627A CN202210957162.9A CN202210957162A CN115371627A CN 115371627 A CN115371627 A CN 115371627A CN 202210957162 A CN202210957162 A CN 202210957162A CN 115371627 A CN115371627 A CN 115371627A
Authority
CN
China
Prior art keywords
strain
cabin
lightweight
strain response
response
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
Application number
CN202210957162.9A
Other languages
English (en)
Inventor
魏钧涛
李昭华
孙玲玉
张雷
张法业
姜明顺
隋青美
贾磊
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University
Original Assignee
Shandong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shandong University filed Critical Shandong University
Priority to CN202210957162.9A priority Critical patent/CN115371627A/zh
Publication of CN115371627A publication Critical patent/CN115371627A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
    • G01B21/32—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring the deformation in a solid
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M7/00—Vibration-testing of structures; Shock-testing of structures
    • G01M7/02—Vibration-testing by means of a shake table
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • G01N3/10—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces generated by pneumatic or hydraulic pressure
    • G01N3/12—Pressure testing
    • 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
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/0202—Control of the test
    • G01N2203/0212—Theories, calculations
    • G01N2203/0216—Finite elements
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/06—Indicating or recording means; Sensing means
    • G01N2203/067—Parameter measured for estimating the property
    • G01N2203/0676—Force, weight, load, energy, speed or acceleration
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/06—Indicating or recording means; Sensing means
    • G01N2203/067—Parameter measured for estimating the property
    • G01N2203/0682—Spatial dimension, e.g. length, area, angle
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Biochemistry (AREA)
  • Chemical & Material Sciences (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本公开属于应变重构技术领域，具体涉及一种轻量化舱体结构应变场重构方法及***，包括：构建轻量化舱体结构模型；对所构建的轻量化舱体结构模型进行不同载荷下的分析，得到不同载荷下的应变响应；根据所得到的不同载荷下的应变响应，计算轻量化舱体的载荷应变响应；对所构建的轻量化舱体结构模型进行不同模态下的振动分析，得到模态应变响应；根据所得到的模态应变响应，计算轻量化舱体的振动应变响应；基于所得到的轻量化舱体的载荷应变响应和振动应变响应，计算轻量化舱体的应变响应，完成轻量化舱体的结构应变场重构。

Description

一种轻量化舱体结构应变场重构方法及***
技术领域
本公开属于应变重构技术领域，具体涉及一种轻量化舱体结构应变场重构方法及***。
背景技术
本部分的陈述仅仅是提供了与本公开相关的背景技术信息，不必然构成在先技术。
航天器密封舱是航空航天技术领域中的关键结构之一，主要的功能是在轻量化设计的前提下保持舱体的密封状态和承受载荷。它在空间运行时，一方面要承受一定的内压环境，另一方面也要承受和传递运载发射加速度载荷等载荷的作用。在运行过程中，舱内需要维持一定的压力，以满足舱内航天员和仪器设备的需要。由于服役环境恶劣，密封舱在运行时各种复杂的载荷容易引起结构的损伤，损伤累积容易造成密封舱的功能失效，直接影响密封舱的运行安全，而应变数据是结构疲劳强度和损伤预测的重要依据，因此开展针对密封舱结构的应变监测意义重大。
应变测量可以分为非接触式测量(如摄影测量、相机测量和激光扫描)和接触式测量，其中非接触式测量是通过获得变形前后结构表面的图像，以获得坐标变形前后的匹配点，后续通过计算得到结构的应变和位移，但该方法容易受环境因素影响，例如摇晃和振动，对设备的安装也有空间要求。接触式测量方法包括电阻应变传感器、加速度传感器和光纤布拉格光栅(FBG)传感器。其中，FBG传感器不仅具有灵敏度高、体积小、抗电磁干扰等优点，而且可以进行空分复用和波分复用，易于形成分布式传感网络。因此，本发明将采用FBG传感器应变测量技术。
但是由于密封舱结构复杂，难以在部分关键位置甚至所有的位置上均安装传感器，应变场重构算法为应变监测提供了一种解决思路。应变场重构是指通过获取结构上少量测点的应变值，通过算法，计算出结构上其余位置的应变。目前主要有两种算法，一种是插值法，一种是力学分析法。插值法就是对测量的少量点的应变值进行插值，通过拟合函数计算获得结构上其余位置的应变值。力学分析法是对结构的力学模型进行分析，根据结构的力学特性及受载情况，获得结构的应变场。插值法虽然计算简单，但是该方法需要确定准确的插值函数，而复杂函数的插值函数确定起来较为困难。所以对于复杂结构来说，力学分析法更为适用。
依据力学分析的应变场重构算法主要包括模态法(MM)、逆有限元法(iFEM)和Ko位移理论。MM是将结构变形视为各模态下变形的叠加，重建方法对传感器布局、构件形状和模态顺序选择更敏感，在风力涡轮机叶片和简支梁中有更多的应用。最近，MM在中通过实验验证用于重建机翼的静态和动态位移。iFEM遵循有限元离散化原理，通过将结构离散为单元并使用位移插值函数和形式函数构造逆单元，通过测量应变和理论应变的最小二乘插值来求解节点位移。由于求解速度快，算法不依赖于材料性质和外载荷性质，近五年来在海洋工程、航空航天、复合材料等领域得到更广泛的应用，如逆有限元单元迄今为止，iQS4已用于船体加筋和复合材料机翼，均具有令人满意的重建精度。iFEM领域的最新发展体现在引入可变截面梁、弯曲壳结构的等距公式，以及用于阶梯梁结构形状传感的新型等几何分析方法。Ko位移理论是基于欧拉-伯努利梁的力模型，在被测构件受力不明确的情况下，可以实现从构件的实测应变信息到未知位移信息的转换，也不需要具体的被测构件的材料性能参数，对梁结构的变形有较好的影响，但传感器的布置间隔对精度的影响较大。
据发明人了解，虽然许多学者在应变测量和形状重构方面已经做了很多工作，但关于使用直接测量获得的应变数据对整个结构应变场的反演的研究较少。其中，应变场重构的研究对象大部分是简单结构，针对复杂结构的研究较少，大部分只针对应变的单一影响因素进行了研究，比如说载荷造成的应变，但空间舱舱体应变产生因素较多，仅考虑单一因素会造成较大的误差，因此选择恰当的方法对密封舱舱体进行应变场重构是一个有待于解决的难点。
发明内容
为了解决上述问题，本公开提出了一种轻量化舱体结构应变场重构方法及***，以密封舱为研究对象，基于密封舱的服役过程中的受载情况，考虑了密封舱的静力学特性与动力学特性，将载荷/模态—应变矩阵代替了模态叠加法中的应变振型矩阵，避免了模态法中模态测试困难、模态阶数难以确定等问题，能够实现复杂结构的应变场重构；所述方案在线弹性范围内应变场重构的精度更高，重构效果更好。。
根据一些实施例，本公开的第一方案提供了一种轻量化舱体结构应变场重构方法，采用如下技术方案：
一种轻量化舱体结构应变场重构方法，包括：
构建轻量化舱体结构模型；
对所构建的轻量化舱体结构模型进行不同载荷下的分析，得到不同载荷下的应变响应；
根据所得到的不同载荷下的应变响应，计算轻量化舱体的载荷应变响应；
对所构建的轻量化舱体结构模型进行不同模态下的振动分析，得到模态应变响应；
根据所得到的模态应变响应，计算轻量化舱体的振动应变响应；
基于所得到的轻量化舱体的载荷应变响应和振动应变响应，计算轻量化舱体的应变响应，完成轻量化舱体的结构应变场重构。
作为进一步的技术限定，所述不同载荷包括充气载荷、加速度载荷和重力及支撑载荷；所述不同载荷下的应变响应包括充气载荷所产生的应变响应、射过程中加速度载荷所产生的应变响应和垂直转运过程中重力和支撑载荷所产生的应变响应。
作为进一步的技术限定，对所构建的轻量化舱体结构模型进行有限元模态分析，根据模态有效质量原则，确定所截取的轻量化舱体结构模型的模态数量，到不同模态下的应变响应。
作为进一步的技术限定，在计算轻量化舱体的载荷应变响应的过程中，对所得到的不同载荷下的应变响应进行线性叠加，得到轻量化舱体的载荷应变响应；
所述轻量化舱体的载荷应变响应ε'M为
Figure BDA0003791819360000051
其中，M表示监测点，ε'1,ε'2…,ε'M分别为矩阵ε'M中的矩阵元素，
Figure BDA0003791819360000052
Figure BDA0003791819360000053
和
Figure BDA0003791819360000054
分别表示不同载荷的应变响应，ω'1，ω'2和ω'3分别表示不同载荷的应变响应的权重；
Figure BDA0003791819360000055
表示在不同载荷分别作用下M个监测点的应变响应矩阵，且
Figure BDA0003791819360000056
ω'3×1表示不同载荷分别作用下M个监测点的应变响应权重，且
Figure BDA0003791819360000057
进一步的，在计算轻量化舱体的振动应变响应的过程中，对所得到的模态应变响应进行线性叠加，得到轻量化舱体的振动应变响应；
所述轻量化舱体的振动应变响应ε”M为多个载荷同时作用在轻量化舱体上时的振动应变响应，即
Figure BDA0003791819360000061
其中，M表示监测点，ε”1,ε”2…,ε”M分别为矩阵ε”M中的矩阵元素，
Figure BDA0003791819360000062
Figure BDA0003791819360000063
和
Figure BDA0003791819360000064
分别表示不同模态下M个监测点的模态应变响应，ω”1，ω”2和ω”b分别表示不同模态下M个监测点的模态应变响应的权重；
Figure BDA0003791819360000065
表示在不同模态下M个监测点的模态应变响应矩阵，且
Figure BDA0003791819360000066
ω'b×1表示不同模态下M个监测点的模态应变响应权重，且
Figure BDA0003791819360000067
进一步的，在计算轻量化舱体的应变响应的过程中，对所得到的轻量化舱体的载荷应变响应和振动应变响应进行线性叠加，得到轻量化舱体的应变响应；
所述轻量化舱体的应变响应εM为
Figure BDA0003791819360000068
其中，ε1,ε2…,εM分别为矩阵εM中的矩阵元素；
Figure BDA0003791819360000069
表示轻量化舱体上M个监测点的应变响应，且
Figure BDA0003791819360000071
ω(3+b)×1表示轻量化舱体上M个监测点的应变响应的权重，且
Figure BDA0003791819360000072
进一步的，根据最小二乘法可得
Figure BDA0003791819360000073
结合
Figure BDA0003791819360000074
即可得到监测点的轻量化舱体的应变响应；基于监测点进行轻量化舱体结构上未监测点的应变场重构，得到轻量化舱体结构所有点的应变响应，完成轻量化舱体的结构应变场重构。
根据一些实施例，本公开的第二方案提供了一种轻量化舱体结构应变场重构***，采用如下技术方案：
一种轻量化舱体结构应变场重构***，包括：
建模模块，其被配置为构建轻量化舱体结构模型；
载荷应变响应模块，其被配置为对所构建的轻量化舱体结构模型进行不同载荷下的分析，得到不同载荷下的应变响应；根据所得到的不同载荷下的应变响应，计算轻量化舱体的载荷应变响应；
模态振动应变响应模块，其被配置为对所构建的轻量化舱体结构模型进行不同模态下的振动分析，得到模态应变响应；根据所得到的模态应变响应，计算轻量化舱体的振动应变响应；
结构应变场重构模块，其被配置为基于所得到的轻量化舱体的载荷应变响应和振动应变响应，计算轻量化舱体的应变响应，完成轻量化舱体的结构应变场重构。
根据一些实施例，本公开的第三方案提供了一种计算机可读存储介质，采用如下技术方案：
一种计算机可读存储介质，其上存储有程序，该程序被处理器执行时实现如本公开第一方面所述的轻量化舱体结构应变场重构方法中的步骤。
根据一些实施例，本公开的第四方案提供了一种电子设备，采用如下技术方案：
一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的程序，所述处理器执行所述程序时实现如本公开第一方面所述的轻量化舱体结构应变场重构方法中的步骤。
与现有技术相比，本公开的有益效果为：
本公开通过有限的应变测点，实现密封舱整体的应变场重构，避免了在密封舱大量布置传感器进行应变监测的弊端；
综合考虑转运载荷、发射载荷、充气载荷和结构振动的影响，对多种因素引起的应变场进行同时求解，避免了以往应变场重构算法中仅考虑单一因素进行应变场重构的缺点。
附图说明
构成本公开的一部分的说明书附图用来提供对本公开的进一步理解，本公开的示意性实施例及其说明用于解释本公开，并不构成对本公开的不当限定。
图1是本公开实施例一中的轻量化舱体结构应变场重构方法的流程图；
图2是本公开实施例一中的轻量化舱体结构应变场重构方法的整体框图；
图3(a)是本公开实施例一中的轻量化舱体结构的1阶模态应变图；
图3(b)是本公开实施例一中的轻量化舱体结构的2阶模态应变图；
图3(c)是本公开实施例一中的轻量化舱体结构的3阶模态应变图；
图3(d)是本公开实施例一中的轻量化舱体结构的4阶模态应变图；
图3(e)是本公开实施例一中的轻量化舱体结构的5阶模态应变图；
图3(f)是本公开实施例一中的轻量化舱体结构的6阶模态应变图；
图4是本公开实施例一中的重力及支撑载荷作用下的应变图；
图5是本公开实施例一中的充气荷作用下的应变图；
图6是本公开实施例一中的加速度载荷作用下的应变图；
图7是本公开实施例二中的轻量化舱体结构应变场重构***的结构框图。
具体实施方式
下面结合附图与实施例对本公开作进一步说明。
应该指出，以下详细说明都是例示性的，旨在对本公开提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本公开所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本公开的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
在本公开中，术语如“上”、“下”、“左”、“右”、“前”、“后”、“竖直”、“水平”、“侧”、“底”等指示的方位或位置关系为基于附图所示的方位或位置关系，只是为了便于叙述本公开各部件或元件结构关系而确定的关系词，并非特指本公开中任一部件或元件，不能理解为对本公开的限制。
本公开中，术语如“固接”、“相连”、“连接”等应做广义理解，表示可以是固定连接，也可以是一体地连接或可拆卸连接；可以是直接相连，也可以通过中间媒介间接相连。对于本领域的相关科研或技术人员，可以根据具体情况确定上述术语在本实公开中的具体含义，不能理解为对本公开的限制。
在不冲突的情况下，本公开中的实施例及实施例中的特征可以相互组合。
实施例一
本公开实施例一介绍了一种轻量化舱体结构应变场重构方法。
如图1和图2所示的一种轻量化舱体结构应变场重构方法，包括：
构建轻量化舱体结构模型；
对所构建的轻量化舱体结构模型进行不同载荷下的分析，得到不同载荷下的应变响应；
根据所得到的不同载荷下的应变响应，计算轻量化舱体的载荷应变响应；
对所构建的轻量化舱体结构模型进行不同模态下的振动分析，得到模态应变响应；
根据所得到的模态应变响应，计算轻量化舱体的振动应变响应；
基于所得到的轻量化舱体的载荷应变响应和振动应变响应，计算轻量化舱体的应变响应，完成轻量化舱体的结构应变场重构。
航天器铝合金轻量化加筋封闭舱舱体在服役过程中，首先，在垂直转运阶段，由于舱体重力以及支撑载荷产生应变，在发射过程中，由于燃料助推等承受的加速度载荷产生应变，在服役过程中，舱体内部气压和外部气压不同，会由于承受的充气载荷产生应变，以上应变均为由于舱体所受载荷所产生的应变；在舱体发射及服役过程中，也会由于舱体的振动产生应变，振动的来源一般认为是发动机工作产生的振动，以及各种噪声和冲击。因此，当考虑到舱体的整体应变时，可以将其应变分为两种，一种是由于载荷产生的应变，一种是由于舱体的振动产生的应变。根据***动力学原理，可将舱体的应变表示为两种应变的叠加。
如图4、图5和图6所示，舱体由于载荷所受的应变主要有以下三种，分别是垂直转运过程中重力和支撑载荷所产生的应变，充气载荷所产生的应变，发射过程中加速度载荷所产生的应变。搭建铝合金轻量化加筋密封舱的高保真模型，对结构进行有限元仿真，分析密封舱的服役工况，确定密封舱的受载类型，设计载荷的加载方法，获得载荷应变响应。
根据***动力学原理，在弹性范围内，***产生的总应变可以认为是***在不同载荷作用下的应变的线性叠加。当结构承受外力时，针对监测点，轻量化舱体结构的监测点的应变值可以分别表示为：
轻量化舱体的载荷应变响应ε'M为
Figure BDA0003791819360000121
其中，M表示监测点，ε'1,ε'2…,ε'M分别为矩阵ε'M中的矩阵元素，
Figure BDA0003791819360000122
和
Figure BDA0003791819360000123
分别表示不同载荷的应变响应，ω'1，ω'2和ω'3分别表示不同载荷的应变响应的权重；
Figure BDA0003791819360000124
表示在不同载荷分别作用下M个监测点的应变响应矩阵，且
Figure BDA0003791819360000131
ω'3×1表示不同载荷分别作用下M个监测点的应变响应权重，且
Figure BDA0003791819360000132
同理，能够得到轻量化舱体结构上的未监测点的轻量化舱体的载荷应变响应ε'N为
Figure BDA0003791819360000133
其中，N表示未监测点，ε'1,ε'2…,ε'N分别为矩阵ε'N中的矩阵元素，
Figure BDA0003791819360000134
和
Figure BDA0003791819360000135
分别表示不同载荷的应变响应，ω'1，ω'2和ω'3分别表示不同载荷的应变响应的权重；
Figure BDA0003791819360000136
表示在不同载荷分别作用下N个未监测点的应变响应矩阵，且
Figure BDA0003791819360000137
ω'3×1表示不同载荷分别作用下N个未监测点的应变响应权重，且
Figure BDA0003791819360000138
根据模态叠加原理，当轻量化舱体结构产生振动响应时，结构在任意初始条件下的振动响应都是所有固有振型的线性叠加。当结构产生振动响应时，针对监测点，轻量化舱体结构的振动应变响应为多个载荷同时作用在轻量化舱体上时的振动应变响应，即：
Figure BDA0003791819360000141
其中，M表示监测点，ε”1,ε”2…,ε”M分别为矩阵ε”M中的矩阵元素，
Figure BDA0003791819360000142
和
Figure BDA0003791819360000143
分别表示不同模态下M个监测点的模态应变响应，ω”1，ω”2和ω”b分别表示不同模态下M个监测点的模态应变响应的权重；
Figure BDA0003791819360000144
表示在不同模态下M个监测点的模态应变响应矩阵，且
Figure BDA0003791819360000145
ω'b×1表示不同模态下M个监测点的模态应变响应权重，且
Figure BDA0003791819360000146
在本实施例中采用6阶模态分析，其模态应变图分别如图3(a)、图3(b)、图3(c)、图3(d)、图3(e)和图3(f)所示。
同理，能够得到轻量化舱体结构上的未监测点的轻量化舱体的振动应变响应ε”N为
Figure BDA0003791819360000147
其中，N表示未监测点，ε”1,ε”2…,ε”N分别为矩阵ε”N中的矩阵元素，
Figure BDA0003791819360000148
和
Figure BDA0003791819360000149
分别表示不同模态下N个未监测点的模态应变响应，ω”1，ω”2和ω”b分别表示不同模态下N个未监测点的模态应变响应的权重；
Figure BDA0003791819360000151
表示在不同模态下N个监测点的模态应变响应矩阵，且
Figure BDA0003791819360000152
ω'b×1表示不同模态下N个未监测点的模态应变响应权重，且
Figure BDA0003791819360000153
当***由于载荷产生应变或者因为自身振动应变时，不同的载荷会对应其载荷-应变向量，不同的模态对应其模态-应变向量，那么所有的载荷引起的应变可表示为载荷-应变矩阵，所有的模态-应变向量可以表示为模态-应变矩阵，这两个矩阵的列与列之间不存在线性关系，那么，其总应变可以以结构所承担载荷造成的应变与模态应变的叠加表示。因此本实施例将结构的载荷引起的应变与振动引起的应变进行叠加，针对监测点，轻量化舱体的应变响应εM为
Figure BDA0003791819360000154
其中，ε1,ε2…,εM分别为矩阵εM中的矩阵元素；
Figure BDA0003791819360000155
表示轻量化舱体上M个监测点的应变响应，且
Figure BDA0003791819360000161
ω(3+b)×1表示轻量化舱体上M个监测点的应变响应的权重，且
Figure BDA0003791819360000162
针对未监测点，轻量化舱体的应变响应εN为
Figure BDA0003791819360000163
其中，ε1,ε2…,εN分别为矩阵εN中的矩阵元素；
Figure BDA0003791819360000164
表示轻量化舱体上N个未监测点的应变响应，且
Figure BDA0003791819360000165
ω(3+b)×1表示轻量化舱体上N个未监测点的应变响应的权重，且
Figure BDA0003791819360000166
由于在实际工程中，载荷及模态振型的数量一般大于测点的数量，所以依据最小二乘原理求解可得，
Figure BDA0003791819360000171
再将ω(3+b)×1及通过有限元分析获得的
Figure BDA0003791819360000172
带入公式(6)中，即可求得监测点的应变值，实现应变场的重构。
在具体实施阶段，首先进行结构的静力学仿真，首先依据服役的密封舱结构，创建结构的高保真模型，明确结构的材料属性及边界条件等。根据密封舱的服役工况，在运行过程中，主要重力及支撑力载荷导致的应变，加速度载荷导致的应变，充气载荷导致的应变，以及舱体振动造成得应变，根据上述四种应变情况进行有限元分析，获得对应工况下的应变响应。根据工程经验，确定测点位置，测点的位置以能够安装传感器以及应力集中处为主，也就是易产生损伤的地方。确定测点位置后，从有限元仿真获得的应变响应中提取对应工况的测点的载荷-应变响应数据，即
Figure BDA0003791819360000173
与全场的载荷-应变数据
Figure BDA0003791819360000174
对结构进行模态有限元分析，通过截取的模态有效质量应占全部有效质量的85％为准则，截取模态，确定模态阶数为b，并提取对应模态阶数下的测点以及全部点的模态-应变矩阵即
Figure BDA0003791819360000175
与
Figure BDA0003791819360000176
将受外力作用下的载荷-应变矩阵与振动响应下的模态-应变矩阵结合为载荷/模态-应变矩阵，即
Figure BDA0003791819360000177
随后设计应变测量***，获得结构在服役时下的所述测点的应变，进而重构密封舱的应变场。***主要由传感器、解调仪、光缆组成。根据有限元仿真的结果在结构上安装FBG光纤光栅应变传感器，根据有限元分析结果，以应变最大的方向作为传感器的粘贴方向，当密封舱在服役过程中，结构的应变发生变化时，光纤光栅传感器的中心波长发生改变，传感器将波长变化传递给解调仪进行解调，通过对解调仪的输出数据进行处理，获得结构上传感器位置处即所述测点的应变数据εM。
结合有限元分析得到的载荷/模态-应变矩阵
Figure BDA0003791819360000181
与实验获得的测点的在多种因素共同作用下的实际的应变εM数据，根据式(6)和式(7)即可重构密封舱的应变场。
本实施例基于有限元分析技术，获得结构的标定矩阵即模态-应变矩阵与在基本载荷系下的载荷-应变矩阵，通过载荷/模态-应变线性叠加原理的逆过程求解，实现铝合金轻量化加筋密封舱舱体结构的应变场重构；
首先，搭建铝合金轻量化加筋密封舱的高保真模型，对结构进行有限元仿真，分析密封舱的服役工况，确定密封舱的受载类型，设计载荷的加载方法，获取结构的载荷-应变矩阵；其次，对铝合金轻量化加筋密封舱进行有限元模态分析，根据模态有效质量原则，确定截取的模态数量，获取密封舱的模态-应变矩阵；再次，结合载荷-应变矩阵与模态振型矩阵，构建结构的载荷/模态-应变矩阵；搭建光纤光栅传感器应变测量***，依据有限元数值仿真结果，在舱体上安装传感器并测量应变数据；最后，基于有限元仿真基准矩阵与测量***测得的传感器的应变值，根据叠加原理的逆过程求解，实现密封舱的应变场重构。
实施例二
本公开实施例二介绍了一种轻量化舱体结构应变场重构***。
如图7所示的一种轻量化舱体结构应变场重构***，包括：
建模模块，其被配置为构建轻量化舱体结构模型；
载荷应变响应模块，其被配置为对所构建的轻量化舱体结构模型进行不同载荷下的分析，得到不同载荷下的应变响应；根据所得到的不同载荷下的应变响应，计算轻量化舱体的载荷应变响应；
模态振动应变响应模块，其被配置为对所构建的轻量化舱体结构模型进行不同模态下的振动分析，得到模态应变响应；根据所得到的模态应变响应，计算轻量化舱体的振动应变响应；
结构应变场重构模块，其被配置为基于所得到的轻量化舱体的载荷应变响应和振动应变响应，计算轻量化舱体的应变响应，完成轻量化舱体的结构应变场重构。
详细步骤与实施例一提供的轻量化舱体结构应变场重构方法相同，在此不再赘述。
实施例三
本公开实施例三提供了一种计算机可读存储介质。
一种计算机可读存储介质，其上存储有程序，该程序被处理器执行时实现如本公开实施例一所述的轻量化舱体结构应变场重构方法中的步骤。
详细步骤与实施例一提供的轻量化舱体结构应变场重构方法相同，在此不再赘述。
实施例四
本公开实施例四提供了一种电子设备。
一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的程序，所述处理器执行所述程序时实现如本公开实施例一所述的轻量化舱体结构应变场重构方法中的步骤。
详细步骤与实施例一提供的轻量化舱体结构应变场重构方法相同，在此不再赘述。
以上所述仅为本公开的优选实施例而已，并不用于限制本公开，对于本领域的技术人员来说，本公开可以有各种更改和变化。凡在本公开的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本公开的保护范围之内。
上述虽然结合附图对本公开的具体实施方式进行了描述，但并非对本公开保护范围的限制，所属领域技术人员应该明白，在本公开的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本公开的保护范围以内。

Claims (10)

1.一种轻量化舱体结构应变场重构方法，其特征在于，包括：
构建轻量化舱体结构模型；
对所构建的轻量化舱体结构模型进行不同载荷下的分析，得到不同载荷下的应变响应；
根据所得到的不同载荷下的应变响应，计算轻量化舱体的载荷应变响应；
对所构建的轻量化舱体结构模型进行不同模态下的振动分析，得到模态应变响应；
根据所得到的模态应变响应，计算轻量化舱体的振动应变响应；
基于所得到的轻量化舱体的载荷应变响应和振动应变响应，计算轻量化舱体的应变响应，完成轻量化舱体的结构应变场重构。
2.如权利要求1中所述的一种轻量化舱体结构应变场重构方法，其特征在于，所述不同载荷包括充气载荷、加速度载荷和重力及支撑载荷；所述不同载荷下的应变响应包括充气载荷所产生的应变响应、射过程中加速度载荷所产生的应变响应和垂直转运过程中重力和支撑载荷所产生的应变响应。
3.如权利要求1中所述的一种轻量化舱体结构应变场重构方法，其特征在于，对所构建的轻量化舱体结构模型进行有限元模态分析，根据模态有效质量原则，确定所截取的轻量化舱体结构模型的模态数量，到不同模态下的应变响应。
4.如权利要求1中所述的一种轻量化舱体结构应变场重构方法，其特征在于，在计算轻量化舱体的载荷应变响应的过程中，对所得到的不同载荷下的应变响应进行线性叠加，得到轻量化舱体的载荷应变响应；
所述轻量化舱体的载荷应变响应ε'M为
Figure FDA0003791819350000021
其中，M表示监测点，ε'1,ε'2…,ε'M分别为矩阵ε'M中的矩阵元素，
Figure FDA0003791819350000022
Figure FDA0003791819350000023
和
Figure FDA0003791819350000024
分别表示不同载荷的应变响应，ω'1，ω'2和ω'3分别表示不同载荷的应变响应的权重；
Figure FDA0003791819350000025
表示在不同载荷分别作用下M个监测点的应变响应矩阵，且
Figure FDA0003791819350000026
ω'3×1表示不同载荷分别作用下M个监测点的应变响应权重，且
Figure FDA0003791819350000027
5.如权利要求4中所述的一种轻量化舱体结构应变场重构方法，其特征在于，在计算轻量化舱体的振动应变响应的过程中，对所得到的模态应变响应进行线性叠加，得到轻量化舱体的振动应变响应；
所述轻量化舱体的振动应变响应ε”M为多个载荷同时作用在轻量化舱体上时的振动应变响应，即
Figure FDA0003791819350000028
其中，M表示监测点，ε”1,ε”2…,ε”M分别为矩阵ε”M中的矩阵元素，
Figure FDA0003791819350000031
Figure FDA0003791819350000032
和
Figure FDA0003791819350000033
分别表示不同模态下M个监测点的模态应变响应，ω”1，ω”2和ω”b分别表示不同模态下M个监测点的模态应变响应的权重；
Figure FDA0003791819350000034
表示在不同模态下M个监测点的模态应变响应矩阵，且
Figure FDA0003791819350000035
ω'b×1表示不同模态下M个监测点的模态应变响应权重，且
Figure FDA0003791819350000036
6.如权利要求5中所述的一种轻量化舱体结构应变场重构方法，其特征在于，在计算轻量化舱体的应变响应的过程中，对所得到的轻量化舱体的载荷应变响应和振动应变响应进行线性叠加，得到轻量化舱体的应变响应；
所述轻量化舱体的应变响应εM为
Figure FDA0003791819350000037
其中，ε1,ε2…,εM分别为矩阵εM中的矩阵元素；
Figure FDA0003791819350000038
表示轻量化舱体上M个监测点的应变响应，且
Figure FDA0003791819350000041
ω(3+b)×1表示轻量化舱体上M个监测点的应变响应的权重，且
Figure FDA0003791819350000042
7.如权利要求6中所述的一种轻量化舱体结构应变场重构方法，其特征在于，根据最小二乘法可得
Figure FDA0003791819350000043
结合
Figure FDA0003791819350000044
即可得到监测点的轻量化舱体的应变响应；基于监测点进行轻量化舱体结构上未监测点的应变场重构，得到轻量化舱体结构所有点的应变响应，完成轻量化舱体的结构应变场重构。
8.一种轻量化舱体结构应变场重构***，其特征在于，包括：
建模模块，其被配置为构建轻量化舱体结构模型；
载荷应变响应模块，其被配置为对所构建的轻量化舱体结构模型进行不同载荷下的分析，得到不同载荷下的应变响应；根据所得到的不同载荷下的应变响应，计算轻量化舱体的载荷应变响应；
模态振动应变响应模块，其被配置为对所构建的轻量化舱体结构模型进行不同模态下的振动分析，得到模态应变响应；根据所得到的模态应变响应，计算轻量化舱体的振动应变响应；
结构应变场重构模块，其被配置为基于所得到的轻量化舱体的载荷应变响应和振动应变响应，计算轻量化舱体的应变响应，完成轻量化舱体的结构应变场重构。
9.一种计算机可读存储介质，其上存储有程序，其特征在于，该程序被处理器执行时实现如权利要求1-7中任一项所述的轻量化舱体结构应变场重构方法中的步骤。
10.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的程序，其特征在于，所述处理器执行所述程序时实现如权利要求1-7中任一项所述的轻量化舱体结构应变场重构方法中的步骤。
CN202210957162.9A 2022-08-10 2022-08-10 一种轻量化舱体结构应变场重构方法及*** Pending CN115371627A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210957162.9A CN115371627A (zh) 2022-08-10 2022-08-10 一种轻量化舱体结构应变场重构方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210957162.9A CN115371627A (zh) 2022-08-10 2022-08-10 一种轻量化舱体结构应变场重构方法及***

Publications (1)

Publication Number Publication Date
CN115371627A true CN115371627A (zh) 2022-11-22

Family

ID=84066452

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210957162.9A Pending CN115371627A (zh) 2022-08-10 2022-08-10 一种轻量化舱体结构应变场重构方法及***

Country Status (1)

Country Link
CN (1) CN115371627A (zh)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107103111A (zh) * 2017-03-02 2017-08-29 西安电子科技大学 基于应变传感器的电子装备功能形面特征点位移场重构方法
CN110059373A (zh) * 2019-04-01 2019-07-26 南京航空航天大学 基于模态叠加原理的机翼应变场重构分布式光纤计算方法
CN112989491A (zh) * 2021-03-30 2021-06-18 山东大学 基于载荷应变线性叠加的应变场重构可视化方法及***
US20210389207A1 (en) * 2019-03-22 2021-12-16 Xi'an Jiaotong University Non-Contact Dynamic Strain Field Measuring Method and System for Rotating Blade
CN114580080A (zh) * 2022-03-03 2022-06-03 山东大学 一种基于贝叶斯有限元模型修正的应变场重构方法及***

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107103111A (zh) * 2017-03-02 2017-08-29 西安电子科技大学 基于应变传感器的电子装备功能形面特征点位移场重构方法
US20210389207A1 (en) * 2019-03-22 2021-12-16 Xi'an Jiaotong University Non-Contact Dynamic Strain Field Measuring Method and System for Rotating Blade
CN110059373A (zh) * 2019-04-01 2019-07-26 南京航空航天大学 基于模态叠加原理的机翼应变场重构分布式光纤计算方法
CN112989491A (zh) * 2021-03-30 2021-06-18 山东大学 基于载荷应变线性叠加的应变场重构可视化方法及***
CN114580080A (zh) * 2022-03-03 2022-06-03 山东大学 一种基于贝叶斯有限元模型修正的应变场重构方法及***

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ZHAOHUA LI 等: "Strain Field Reconstruction of High-Speed Train Crossbeam Based on FBG Sensing Network and Load-Strain Linear Superposition Algorithm", 《IEEE SENSOR JOURNAL》, vol. 22, no. 4, 30 December 2021 (2021-12-30), XP011900459, DOI: 10.1109/JSEN.2021.3139648 *
司亚文等: "机翼模型应变场分布式光纤监测与重构方法", 《振动、测试与诊断》, vol. 40, no. 4, 31 August 2020 (2020-08-31) *
孙玲玉等: "基于光纤布拉格光栅阵列的铝合金疲劳裂纹预测方法", 《中国激光》, vol. 48, no. 13, 24 February 2021 (2021-02-24) *
闫洁等: "基于光纤光栅传感器的板状结构形态感知与三维重构技术", 《中国激光》, vol. 47, no. 11, 30 November 2020 (2020-11-30) *

Similar Documents

Publication Publication Date Title
Gherlone et al. Shape sensing methods: Review and experimental comparison on a wing-shaped plate
Pak Wing shape sensing from measured strain
Papa et al. Health structure monitoring for the design of an innovative UAS fixed wing through inverse finite element method (iFEM)
Cerracchio et al. Real-time displacement monitoring of a composite stiffened panel subjected to mechanical and thermal loads
Tessler et al. Inverse FEM for full-field reconstruction of elastic deformations in shear deformable plates and shells
CN108846149B (zh) 一种基于多源不确定性的结构分布式动态载荷识别的方法
Tessler Structural analysis methods for structural health management of future aerospace vehicles
Robinson et al. Design and implementation of an internal stress wave force balance in a shock tunnel
CN103049608B (zh) 基于约束边应变极值坐标的载荷辨识***及方法
CN112613129A (zh) 一种基于逆有限元与微元动态响应法的复合材料板壳结构健康监测方法
CN115879346A (zh) 基于改进型四节点逆有限元理论的结构应变场反演方法
Yuan et al. An equivalent modeling method for honeycomb sandwich structure based on orthogonal anisotropic solid element
CN108038277A (zh) 一种航天器有限元模型的二次缩聚方法
Alioli et al. Membrane shape and load reconstruction from measurements using inverse finite element analysis
Tsushima et al. Monitoring multi-axial vibrations of flexible rockets using sensor-instrumented reference strain structures
CN116205115B (zh) 基于逆元法与虚实结合技术的结构形态反演精度提升方法
Bottasso et al. Estimation of blade structural properties from experimental data
Colombo et al. Numerical and experimental verification of an inverse‐direct approach for load and strain monitoring in aeronautical structures
CN115371627A (zh) 一种轻量化舱体结构应变场重构方法及***
Liu et al. Deformation reconstruction for a heavy-duty machine column through the inverse finite element method
CN105404740B (zh) 一种基于联合指标与有限信息的修正灵敏度结构损伤识别方法
Warwick et al. Measurement of aeroelastic wing deflections using modal shapes and strain pattern analysis
Heaney et al. Distributed sensing of a cantilever beam and plate using a fiber optic sensing system
Miglani et al. Theoretical and Experimental Analyses of Inflatable Structures
Kobayashi et al. Three-dimensional shape sensing by inverse finite element method based on distributed fiber-optic sensors

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