CN113358308B - 基于有限测点和全局模态的组合结构横向位移确定方法 - Google Patents
基于有限测点和全局模态的组合结构横向位移确定方法 Download PDFInfo
- Publication number
- CN113358308B CN113358308B CN202110619854.8A CN202110619854A CN113358308B CN 113358308 B CN113358308 B CN 113358308B CN 202110619854 A CN202110619854 A CN 202110619854A CN 113358308 B CN113358308 B CN 113358308B
- Authority
- CN
- China
- Prior art keywords
- global
- combined structure
- displacement
- measuring points
- sensor
- 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
Images
Classifications
-
- 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
- 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
- G01M7/025—Measuring arrangements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
基于有限测点和全局模态的组合结构横向位移确定方法,属于组合结构的横向振动位移确定方法领域。本发明解决了实验或实际工程中对组合结构横向位移、速度等物理场的有效测量与实验设备复杂,成本高,数据结果精度低等问题,在组合结构上合理布置传感器,测量并记录各测点的位移和速度,采用理论计算或有限元仿真获取组合结构的全局模态,获取可逆M方阵,计算解得组合结构中任意结构上任意点的位移和速度。本发明所使用的装置简单,成本低,计算效率高,相对精度高,通过有限个测点的横向位移数据和组合结构全局模态,便可实现组合结构全位移场的构建,有着丰富的应用前景和现实意义。
Description
技术领域
本发明属于组合结构的横向振动位移确定方法领域,具体涉及一种基于有限测点和全局模态的组合结构横向位移确定方法。
背景技术
组合结构是由多个物体(包括柔体、刚体等)以某种方式相联接而构成的多体***,被广泛应用在航空航天、机械、土木等工程领域中,如带太阳能帆板的空间组合体、空间柔性机械臂、带细长叶片的大型旋转机械等。在组合结构中,由于各部件之间的振动耦合效应,***中单一柔性部件在悬臂、简支和自由等静定边界下的动力学行为与整体结构的动力学行为有很大差异。研究组合结构的振动特性对于满足工程愈加增长的结构需求和提升***安全性有重大意义。
在组合结构的预研实验或实际工程健康监测中,组合结构振动的位移、速度等物理量是必要的数据。而对于一些特别的组合结构,受工作环境苛刻、载荷复杂等因素的影响,工件变形及振动具有时空不均匀性和强时变性,仅有单点或有限点测量信息不足以描述实际工况,更需要提供横向位移、速度等物理量的场信息。如在组合结构形式飞行器气动弹性分析地面实验中,结构各点所受到的气动力是对应结构各点的横向位移和速度的函数,故只有确定结构的横向位移场和速度场才能提供精确的气动力场。
当前确定结构横向位移场的方式有三种:一是先测量结构上部分点横向位移,再通过插值法确定整个结构的位移场。尽管该方法易操作,计算量不大,但是难以对误差进行估计,且插值效果依赖于采样点数量,即若想要得到精确的结果需要大量的采样点。此外,由于算法特征,不同的插值法存在其特有的缺陷,如基于曲率的位移插值方法(CBDI)在边点区域内位移值会有震荡效应而使得到的结果精度较差;二是利用激光测振仪实现对结构全覆盖的横向位移测量。激光测振具有精度高,响应快,动态范围大等优点。但该设备对使用环境条件极其苛刻,测量过程受其他杂散光影响较大,且设备复杂,造价高。此外,在结构振幅较大时,激光测振仪对结构边角处的测量会出现漏采、失真问题;三是基于有限测点和振动模态的板/梁结构横向位移确定方法,这是一种简单且精度效果好的测量板/梁结构位移、速度场的方法,但应用在组合结构上可能会出现严重的误差。在组合结构中,由于各部件之间的振动耦合效应,***中单一柔性部件在悬臂、简支和自由等静定边界下的模态与整体结构的模态有很大差异。故采用部件级模态可能会导致所得到位移场和速度场准确度下降。
综上所述,在组合结构横向振动实验或实际工程应用中,存在着对组合结构横向位移、速度等物理场的有效测量与实验设备复杂,成本高,数据结果精度低等问题。
发明内容
本发明解决了实验或实际工程中对组合结构横向位移、速度等物理场的有效测量与实验设备复杂,成本高,数据结果精度低等问题,提供一种基于有限测点和全局模态的组合结构横向位移确定方法。
为实现上述目的,本发明采取的技术方案如下:
基于有限测点和全局模态的组合结构横向位移确定方法,所述方法具体过程为:
在组合结构振动实验中,针对组合结构横向位移w(x,y,t)的测量,在结构上表面或下表面布置m个传感器,其m个传感器铺设位置为sn(xn,yn),测量位移wn和速度其中n=1…m,xn为传感器的横坐标,yn为传感器的纵坐标;
对于确定的组合结构,第k个结构上任一点的横向位移可表示为
其中,Φik(x,y)为满足组合结构匹配条件和边界条件的全局模态振型函数,qi(t)是对应于Φi(x,y)的广义坐标;
其中,M为由组合结构前m阶全局模态向量组成的矩阵,每一阶全局模态向量包含m个测点的模态;在这里需要注意的是,需要截取与测点数目相等的,且互不相关的全局模态,以保证M是可逆方阵;M表达式如下:
进而计算组合结构上任意一点(xr,yr)的横向位移和速度
本发明相对于现有技术的有益效果为:本发明所使用的装置简单,成本低,计算效率高,相对精度高,通过有限个测点的横向位移数据和组合结构全局模态,便可实现组合结构全位移场的构建,有着丰富的应用前景和现实意义。
附图说明
图1是本发明原理流程图;
图2是实施例1中铰链连接多梁结构及电涡流传感器布置示意图;
图3是实施例2中铰链连接多板结构及电涡流传感器布置示意图。
具体实施方式
下面将参照附图1至附图3更详细地描述本发明的具体实施例。然而应当注意,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本发明,并且能够将本发明的范围完整地传达给本领域的技术人员。
全局模态指组合结构的包含***刚体运动和全部柔性部件弹性振动的模态。
具体实施方式一:本实施方式记载的是一种基于有限测点和全局模态的组合结构横向位移确定方法,所述方法具体过程为:
在组合结构振动实验中,针对组合结构横向位移w(x,y,t)的测量,其中w(x,y,t)表示在t时刻组合结构上(x,y,z)处的横向位移,w表示横向位移,x,y是组合结构上的坐标,t表示时间。在结构上表面或下表面布置m个传感器,其m个传感器铺设位置为sn(xn,yn),测量位移wn和速度其中n=1…m,xn为传感器所在位置的横坐标,yn为传感器所在位置的纵坐标;
对于确定的组合结构,第k个结构上任一点的横向位移可表示为
其中,Φik(x,y)为满足组合结构匹配条件和边界条件的全局模态振型函数,qi(t)是对应于Φi(x,y)的广义坐标;举例说明:表示的是Φ1(x,y)q1(t)+Φ2(x,y)q2(t)+Φ3(x,y)q3(t))。
其中,变量上带点表示该变量对时间的导数,M为由组合结构前m阶全局模态向量组成的矩阵,每一阶全局模态向量包含m个测点的模态;组合结构的全局模态通过理论计算、有限元仿真或模态试验获得。在这里需要注意的是,需要截取与测点数目相等的,且互不相关的全局模态,以保证M是可逆方阵;M表达式如下:
进而计算组合结构上任意一点(xr,yr)的横向位移和速度
具体实施方式二:具体实施方式一所述的基于有限测点和全局模态的组合结构横向位移确定方法,传感器等距均匀地布置在结构上。
具体实施方式三:具体实施方式二所述的基于有限测点和全局模态的组合结构横向位移确定方法,根据结构特性改善传感器的布置方式,避开结构振型的节点位置。
实施例1:
应用在铰链连接的多梁结构横向振动实验中,结合图1说明本实施方式,本发明具体按以下操作实现:
第k根梁上任一点的横向位移可表示为
其中Φik(x)=Ak cos(βkx)+Bk sin(βkx)+Ck cosh(βkx)+Dk sinh(βkx),k=1,2,3为多梁结构满足匹配条件和边界条件的全局模态振型函数。A,B,C,D是由梁的边界条件所决定的常数,β是由梁的参数所确定的常数;
其中M为由梁结构前3阶模态向量组成的矩阵,每一阶模态向量包含3个测点的模态。梁的模态可以通过理论计算、有限元仿真或模态试验获得。在这里需要注意的是,需要截取与测点数目相等的,且互不相关的模态,以保证M是可逆方阵。M表达式如下:
进而可以计算解得第k根梁上任意一点sr(xr)的横向位移和速度
实施例2:
应用在铰链连接多板结构振动实验中,结合图1说明本实施方式,本发明具体按以下操作实现:
第k块板上任一点的横向位移可表示为
其中M为由结构前6阶模态向量组成的矩阵,每一阶模态向量包含6个测点的模态。多板的全局模态可以通过理论计算、有限元仿真或模态试验获得。在这里需要注意的是,需要截取与测点数目相等的,且互不相关的全局模态,以保证M是可逆方阵。M表达式如下:
进而可以计算解得第k块板上任意一点sr(xr,yr)的横向位移和速度
尽管以上结合附图对本发明的具体实施方案进行了描述,但本发明不局限于以上阐述的两个具体实施方案和应用领域,上述的具体实施方案仅仅是示意性的和指导性的,而不是限制性的。本领域的普通技术人员在本说明书的启示下和在不脱离本发明权利要求所保护的范围的情况下,还可以做出很多种的形式和应用,这些均属于本发明保护之列。
Claims (4)
1.基于有限测点和全局模态的组合结构横向位移确定方法,其特征在于:所述方法具体过程为:
在组合结构振动实验中,针对组合结构横向位移w(x,y,t)的测量,在结构上表面或下表面布置m个传感器,其m个传感器铺设位置为sn(xn,yn),测量位移wn和速度其中n=1…m,xn为传感器的横坐标,yn为传感器的纵坐标;
对于确定的组合结构,第k个结构上任一点的横向位移表示为
对于具有n个梁的多梁结构,其满足多梁匹配条件和边界条件全局模态振型函数为
Φik(x)=Ak cos(βkx)+Bk sin(βkx)+Ck cosh(βkx)+Dk sinh(βkx),k=1,2,…,n (2)
A,B,C,D是由梁的边界条件所决定的常数,β是由梁的参数所确定的常数;
对于具有n个板的多板结构,其满足多板匹配条件和边界条件全局模态振型函数为
其中,Φik(x,y)为满足组合结构匹配条件和边界条件的全局模态振型函数,qi(t)是对应于Φik(x,y)的广义坐标;
其中,M为由组合结构前m阶全局模态向量组成的矩阵,每一阶全局模态向量包含m个测点的模态;在这里需要注意的是,需要截取与测点数目相等的,且互不相关的全局模态,以保证M是可逆方阵;M表达式如下:
进而计算组合结构上任意一点(xr,yr)的横向位移和速度
2.根据权利要求1所述的基于有限测点和全局模态的组合结构横向位移确定方法,其特征在于:所述传感器不局限于常规传感器,还包括电涡流传感器或激光传感器等非接触式传感器。
3.根据权利要求1所述的基于有限测点和全局模态的组合结构横向位移确定方法,其特征在于:传感器等距均匀布置在结构上。
4.根据权利要求1所述的基于有限测点和全局模态的组合结构横向位移确定方法,其特征在于:根据结构特性改善传感器的布置方式,避开结构振型的节点位置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110619854.8A CN113358308B (zh) | 2021-06-03 | 2021-06-03 | 基于有限测点和全局模态的组合结构横向位移确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110619854.8A CN113358308B (zh) | 2021-06-03 | 2021-06-03 | 基于有限测点和全局模态的组合结构横向位移确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113358308A CN113358308A (zh) | 2021-09-07 |
CN113358308B true CN113358308B (zh) | 2022-10-25 |
Family
ID=77531663
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110619854.8A Active CN113358308B (zh) | 2021-06-03 | 2021-06-03 | 基于有限测点和全局模态的组合结构横向位移确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113358308B (zh) |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1991019173A1 (en) * | 1990-06-01 | 1991-12-12 | Technology Integration And Development Group Incorporated | Method for assessing structural integrity of composite structures |
CN104589359A (zh) * | 2014-12-25 | 2015-05-06 | 中国矿业大学 | 一种基于振动观测器的柔性机械臂振动控制方法 |
CN105631090A (zh) * | 2015-12-02 | 2016-06-01 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种有限元模型优化装置及方法 |
CN107292054A (zh) * | 2017-07-12 | 2017-10-24 | 哈尔滨工程大学 | 一种复杂船舶推进轴系扭转振动特性分析方法 |
CN107345786A (zh) * | 2017-08-29 | 2017-11-14 | 华南理工大学 | 一种测量柔性梁横向振动位移和应变关系的装置及方法 |
CN108228958A (zh) * | 2017-11-08 | 2018-06-29 | 中国舰船研究设计中心 | 基于模态强度的振动场监测优化方法 |
CN110375690A (zh) * | 2019-03-22 | 2019-10-25 | 西安交通大学 | 一种旋转叶片非接触式位移场测量方法及其*** |
KR20210006832A (ko) * | 2019-07-09 | 2021-01-19 | 한국전자통신연구원 | 기계고장 진단 방법 및 장치 |
CN112383242A (zh) * | 2020-12-09 | 2021-02-19 | 哈尔滨工业大学 | 一种带尾翼的薄板框架结构直线超声电机定子及其激励方法 |
CN112629786A (zh) * | 2020-12-03 | 2021-04-09 | 华侨大学 | 一种工作模态参数识别方法及设备故障诊断方法 |
CN112720456A (zh) * | 2020-12-02 | 2021-04-30 | 重庆大学 | 基于双时间尺度的柔性关节柔性杆机械臂控制方法及*** |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0507618D0 (en) * | 2005-04-15 | 2005-05-25 | Lms Internat Nv | Method and system for dynamic analysis of complex systems |
TW201035416A (en) * | 2009-03-26 | 2010-10-01 | Accuracy Structural Technology Corp | Method for evaluating bridge safety with vibration measurement |
US10473553B2 (en) * | 2015-02-26 | 2019-11-12 | City University Of Hong Kong | System and a method for performing modal analysis on a structure |
US10318798B2 (en) * | 2015-07-10 | 2019-06-11 | Booz Allen Hamilton Inc. | Device and method for detecting non-visible content in a non-contact manner |
CN106777825B (zh) * | 2017-01-24 | 2020-03-27 | 湖南科技大学 | 一种基于谱有限元的矩形板振动模态计算方法 |
IT201800006499A1 (it) * | 2018-06-20 | 2019-12-20 | Procedimento per la diagnostica di una struttura sottoposta a carichi basato sulla misura di spostamenti, e sistema per l'attuazione di detto procedimento. | |
CN109918614B (zh) * | 2019-03-14 | 2022-10-14 | 合肥工业大学 | 一种基于模态学习的全局动应变测量方法 |
CN109933916B (zh) * | 2019-03-18 | 2020-07-31 | 山东大学 | 基于轴系纵向振动响应测量反演螺旋桨纵向激励的方法及*** |
-
2021
- 2021-06-03 CN CN202110619854.8A patent/CN113358308B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1991019173A1 (en) * | 1990-06-01 | 1991-12-12 | Technology Integration And Development Group Incorporated | Method for assessing structural integrity of composite structures |
CN104589359A (zh) * | 2014-12-25 | 2015-05-06 | 中国矿业大学 | 一种基于振动观测器的柔性机械臂振动控制方法 |
CN105631090A (zh) * | 2015-12-02 | 2016-06-01 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种有限元模型优化装置及方法 |
CN107292054A (zh) * | 2017-07-12 | 2017-10-24 | 哈尔滨工程大学 | 一种复杂船舶推进轴系扭转振动特性分析方法 |
CN107345786A (zh) * | 2017-08-29 | 2017-11-14 | 华南理工大学 | 一种测量柔性梁横向振动位移和应变关系的装置及方法 |
CN108228958A (zh) * | 2017-11-08 | 2018-06-29 | 中国舰船研究设计中心 | 基于模态强度的振动场监测优化方法 |
CN110375690A (zh) * | 2019-03-22 | 2019-10-25 | 西安交通大学 | 一种旋转叶片非接触式位移场测量方法及其*** |
KR20210006832A (ko) * | 2019-07-09 | 2021-01-19 | 한국전자통신연구원 | 기계고장 진단 방법 및 장치 |
CN112720456A (zh) * | 2020-12-02 | 2021-04-30 | 重庆大学 | 基于双时间尺度的柔性关节柔性杆机械臂控制方法及*** |
CN112629786A (zh) * | 2020-12-03 | 2021-04-09 | 华侨大学 | 一种工作模态参数识别方法及设备故障诊断方法 |
CN112383242A (zh) * | 2020-12-09 | 2021-02-19 | 哈尔滨工业大学 | 一种带尾翼的薄板框架结构直线超声电机定子及其激励方法 |
Non-Patent Citations (2)
Title |
---|
复杂柔性结构全局模态方法及其在空间组合结构中的应用;魏进;《中国博士学位论文全文数据库-基础科学辑》;20200215(第02期);全文 * |
多环盘组合结构的横向弯曲耦合振动建模与求解;付俊帆等;《工程设计学报》;20170628(第03期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113358308A (zh) | 2021-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9989973B2 (en) | Measurement control system for multi-shaft supported air floatation platform | |
CN101746510B (zh) | 一种基于激光测量技术的前缘襟翼的装配方法 | |
CN106289246B (zh) | 一种基于位置和姿态测量***的柔性杆臂测量方法 | |
Chwalowski et al. | Preliminary computational analysis of the (hirenasd) configuration in preparation for the aeroelastic prediction workshop | |
CN101738203B (zh) | 挠性陀螺仪静态漂移零次和一次加速度相关项误差模型最优位置标定方法 | |
CN104807414A (zh) | 基于分布式光纤传感技术的地铁隧道沉降变形监测方法 | |
CN105424024B (zh) | 一种基于全站仪的空间目标的位置和朝向标定方法 | |
CN109282785B (zh) | 一种基于应变监测的弹性支承板的变形监测方法 | |
KR20130094071A (ko) | 부유식 해상풍력발전 하부구조의 동적피로하중산출 시뮬레이션 방법 | |
Attar et al. | Modeling delta wing limit-cycle oscillations using a high-fidelity structural model | |
CN103808286A (zh) | 一种基于全站仪的钢结构三维精度检测分析方法及其应用 | |
CN105547235A (zh) | 变截面梁结构拉压与弯曲复合变形场的测量方法 | |
CN113358308B (zh) | 基于有限测点和全局模态的组合结构横向位移确定方法 | |
CN109543227A (zh) | 大型反射面天线的风致指向误差实时测量方法、雷达天线 | |
CN113587819A (zh) | 一种基于三路激光跟踪的大尺度空间位姿动态测量方法与测量精度验证方法 | |
Chwalowski et al. | FUN3D analyses in support of the First Aeroelastic Prediction Workshop | |
Carta | Unsteady aerodynamics and gapwise periodicity of oscillating cascaded airfoils | |
CN113358311B (zh) | 基于有限测点和振动模态的板/梁结构横向位移确定方法 | |
Gang et al. | Optimal load balancing leveling method for multi-leg flexible platforms | |
CN103177192A (zh) | 五孔探针角度特性曲线数据处理方法 | |
CN104330078B (zh) | 一种基于三点后方交会模型的联合测量方法 | |
CN106323587A (zh) | 机翼风洞试验模型弹性变形的单目视频高精度测量方法 | |
Chwalowski et al. | Collaborative HIRENASD analyses to eliminate variations in computational results | |
CN114659709A (zh) | 一种用于大型带翼航天飞行器的三维质心测量方法 | |
Pham et al. | Real-time aerodynamic load estimation for hypersonics via strain-based inverse maps |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |