CN109163835B - 基于应变识别多接头的接头载荷测量方法 - Google Patents
基于应变识别多接头的接头载荷测量方法 Download PDFInfo
- Publication number
- CN109163835B CN109163835B CN201811205768.7A CN201811205768A CN109163835B CN 109163835 B CN109163835 B CN 109163835B CN 201811205768 A CN201811205768 A CN 201811205768A CN 109163835 B CN109163835 B CN 109163835B
- Authority
- CN
- China
- Prior art keywords
- joint
- load
- strain
- coefficient matrix
- method based
- 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
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L5/00—Apparatus for, or methods of, measuring force, work, mechanical power, or torque, specially adapted for specific purposes
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
- B64F5/60—Testing or inspecting aircraft components or systems
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Force Measurement Appropriate To Specific Purposes (AREA)
Abstract
基于应变识别多接头的接头载荷测量方法,首先拆开对接的结构件,而后选择其中一个结构件标定,选择第一个接头A1沿x方向逐级加载,一共加载8级,记录载荷F1与每一个应变εi1的数值,线性拟合出载荷F1与每一个应变εi1之间的关系直线,得到3m个斜率ki1,对每一个接头任一个方向都按照相同的方式标定,形成应变与载荷之间的传递系数矩阵,并对传递系数矩阵进行线性叠加处理,通过测量应变以识别不同结构件之间的多个连接铰点的载荷,尤其是飞机机翼接头的连接载荷,通过测量标定应变与接头载荷之间的传递关系,而后反演具有过约束系统的接头载荷,从而实现对结构内部接头载荷的测量,有效克服通过力传感器来测量载荷的不足,测量出接头载荷对结构设计优化有着重要意义。
Description
技术领域
本发明涉及接头载荷测量技术领域,尤其涉及一种基于应变识别多接头的接头载荷测量方法。
背景技术
飞机机翼与机身或其他一些需要通过多个接头连接在一起的结构件,由于接头较多而形成过约束,其接头载荷的大小无法通过力传感器直接测量,也难以根据实际的接头连接情况建立准确的有限元模型,而具体每一个接头传递多少载荷对结构设计和优化有着重要的意义。国内外已经有学者发表论文研究边界条件的识别问题,是将整个边界Γ划分为未知部分Γ1和已知部分Γ2,通过有限元或边界元利用Γ2的边界条件(力和位移)反演Γ1的边界条件(力和位移),从理论上取得了一定的进展,但也存在一定不足:首先,Γ2的边界条件在实际条件下往往是没有办法获取;其次,研究的载荷形式往往都是均匀分布的载荷,但是实际情况中的接头载荷是集中载荷,且每一个接头的载荷一般都不同。中国专利ZL201410154031.2公开一种载荷测量方法,该方法的缺点是难以测量对多个接头的传载;中国专利ZL201710804690.X公开一种连接翼结构形式机翼载荷测量方法,该方法需要事先根据结构的受力特点来布置应变片,一旦结构连接形式过于复杂就难以实施;中国专利ZL201510021235.3公开一种基于Tikhonov正则化的频域载荷识别方法,中国专利ZL201310073284.2公开一种基于正则化和奇异值分解的结构载荷识别方法,这两种方法都是基于Tikhonov正则化方法的载荷识别方法,本质上都要求方程满足Picard条件。对于载荷识别问题或者一般的病态问题可能会出现方程不满足Picard条件。
发明内容
本发明所解决的技术问题在于提供基于应变识别多接头的接头载荷测量方法,以解决上述背景技术中的缺点。
本发明所解决的技术问题采用以下技术方案来实现:
基于应变识别多接头的接头载荷测量方法,具体步骤如下:
1)标定载荷,首先拆开对接的结构件,而后选择其中一个结构件标定,再将第一个接头A1沿x方向逐级加载,一共加载8级,记录载荷F1与每一个应变εi1的数值,线性拟合出载荷F1与每一个应变εi1之间的关系直线,得到3m个斜率ki1,于是:
对于每一个接头任一个方向都按照相同的方式标定,形成应变与载荷之间的传递系数矩阵
标定传递系数矩阵中,应变片的数目远大于需要测量的接头载荷数目;
2)将两个结构件,按照实际情况连接,在第二个结构件上加载,测量应变并记为ε′1,ε′2,…,ε′3n;
假定第二个结构件上的每一个接头载荷分别为:F1',F2',…,F′3m,因结构在线性范围内,根据线性叠加原理具有如下关系:
于是可得到:
简记为:
{ε′j}=[kij]{F′i} (5)
因此,得到关于载荷和应变之间的一个过定义线性方程组式(5),再采用奇异值分解求解该过定义线性方程组式(5)的解:
若始终无法获得一个良态的传递系数矩阵[kij],或传递系数矩阵[kij]本质就是一个病态传递系数矩阵,对于基于测量应变反演识别接头载荷问题,常见的情况是问题的本质是一个病态问题,因此传递系数矩阵是病态,此时不能直接求解,故将病态矩阵的正则化方法引入至基于测量应变反演识别接头载荷问题中,于是通过下式反演计算:
确定截断的奇异值阶数,将取得更好的效果,此时Δ一般取值0.02,具体算法如下:
第一步,对传递系数矩阵[kij]作奇异值分解,并检查传递系数矩阵[kij]的接头载荷数目;
有益效果:本发明通过测量应变以识别不同结构件之间的多个连接铰点的载荷,尤其是飞机机翼接头的连接载荷,通过测量标定应变与接头载荷之间的传递关系,而后反演具有过约束系统的接头载荷,从而实现对结构内部接头载荷的测量,有效克服传统通过力传感器来测量载荷的不足,测量出接头载荷对结构的设计优化有着重要意义。
附图说明
图1是本发明的较佳实施例中的具有m个接头的任意弹性体标定载荷示意图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体图示,进一步阐述本发明。
基于应变识别多接头的接头载荷测量方法,具体步骤如下:
1)标定载荷,首先拆开对接的结构件,而后选择其中一个结构件标定,如图1所示,选择第一个接头A1沿x方向逐级加载,一共加载8级,记录载荷F1与每一个应变εi1的数值,线性拟合出载荷F1与每一个应变εi1之间的关系直线,得到3m个斜率ki1,于是:
对于每一个接头任一个方向都按照相同的方式标定,形成应变与载荷之间的传递系数矩阵
标定传递系数矩阵需要注意,由于是反演识别载荷,故应变片的数目要远远大于需要测量的接头载荷数目,只有这样才能够获得比较好的反演解;
2)将两个结构件,比如机翼机身,按照实际情况连接,然后在机翼上加载,测量应变并记为ε′1,ε′2,...,ε′3n;
假定机翼上的每一个接头载荷分别为:F1',F2',...,F′3m,因结构在线性范围内,根据线性叠加原理具有如下关系:
于是可得到:
简记为:
{ε′j}=[kij]{F′i} (5)
因此,得到关于载荷和应变之间的一个过定义线性方程组式(5),再采用奇异值分解求解该过定义线性方程组式(5)的解:
若始终无法获得一个良态的传递系数矩阵[kij],或问题的本质就是一个病态问题,对于基于测量应变反演识别接头载荷问题,常见的情况是问题的本质是一个病态问题,因此传递系数矩阵是病态,此时不能直接求解,故将病态矩阵的正则化方法引入至基于测量应变反演识别接头载荷问题中,于是通过下式反演计算:
确定截断的奇异值阶数,会取得更好的效果,此时Δ一般取值0.02,具体算法如下:
第一步,对传递系数矩阵[kij]作奇异值分解,并检查传递系数矩阵[kij]的接头载荷数目;
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (8)
1.基于应变识别多接头的接头载荷测量方法,其特征在于,具体步骤如下:
1)标定载荷,首先拆开对接的结构件,而后选择其中一个结构件标定,再将第一个接头A1沿x方向逐级加载,记录载荷F1与每一个应变εi1的数值,而后得到3m个斜率ki1,于是:
对于每一个接头任一个方向都按照相同的方式标定,形成应变与载荷之间的传递系数矩阵
2)将两个结构件,按照实际情况连接,在第二个结构件上加载,测量应变并记为ε′1,ε′2,…,ε′3n;
假定在第二个结构件上的每一个接头载荷分别为:F1',F2',…,F′3m,因结构在线性范围内,根据线性叠加原理具有如下关系:
于是得到:
简记为:
{ε′j}=[kij]{F′i} (5)
因此,得到关于载荷和应变之间的一个过定义线性方程组式(5),再采用奇异值分解求解该过定义线性方程组式(5)的解:
若始终无法获得一个良态的传递系数矩阵[kij],或传递系数矩阵[kij]本质就是一个病态矩阵,此时将病态矩阵的正则化方法引入至基于测量应变反演识别接头载荷问题中,于是通过下式反演计算:
2.根据权利要求1所述的基于应变识别多接头的接头载荷测量方法,其特征在于,步骤1)中,标定传递系数矩阵中,应变片的数目远大于需要测量的接头载荷数目。
3.根据权利要求1所述的基于应变识别多接头的接头载荷测量方法,其特征在于,步骤1)中,3m个斜率ki1由线性拟合出载荷F1与每一个应变εi1之间的关系直线得到。
4.根据权利要求1所述的基于应变识别多接头的接头载荷测量方法,其特征在于,步骤2)中,若传递系数矩阵[kij]是病态,需增加贴片点。
5.根据权利要求4所述的基于应变识别多接头的接头载荷测量方法,其特征在于,再一次标定增加的贴片点传递系数矩阵[kij],并检查传递系数矩阵的接头载荷数目。
6.根据权利要求1所述的基于应变识别多接头的接头载荷测量方法,其特征在于,步骤2)中,采用L-curve法选取截断点λ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811205768.7A CN109163835B (zh) | 2018-10-17 | 2018-10-17 | 基于应变识别多接头的接头载荷测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811205768.7A CN109163835B (zh) | 2018-10-17 | 2018-10-17 | 基于应变识别多接头的接头载荷测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109163835A CN109163835A (zh) | 2019-01-08 |
CN109163835B true CN109163835B (zh) | 2020-11-06 |
Family
ID=64878380
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811205768.7A Active CN109163835B (zh) | 2018-10-17 | 2018-10-17 | 基于应变识别多接头的接头载荷测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109163835B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110243564A (zh) * | 2019-06-12 | 2019-09-17 | 南京航空航天大学 | 拦阻钩挂索冲击载荷测试试验装置及试验方法 |
CN112109919B (zh) * | 2020-04-30 | 2024-04-19 | 中国飞机强度研究所 | 一种强度试验加载点布局方法 |
CN112052518B (zh) * | 2020-08-17 | 2022-11-18 | 大连理工大学 | 一种极地船舶结构冰载荷的远场测量识别方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104536941A (zh) * | 2015-01-15 | 2015-04-22 | 西安交通大学 | 一种基于Tikhonov正则化的频域载荷识别方法 |
CN104794299A (zh) * | 2015-04-29 | 2015-07-22 | 西北工业大学 | 一种复合材料干涉配合接头应力分布计算方法 |
CN103196591B (zh) * | 2013-03-07 | 2015-10-21 | 同济大学 | 一种基于正则化和奇异值分解的结构载荷识别方法 |
CN105043728A (zh) * | 2015-06-23 | 2015-11-11 | 中国航空工业集团公司西安飞机设计研究所 | 一种计算待测接头的拉伸承载能力的方法 |
CN105975704A (zh) * | 2016-05-13 | 2016-09-28 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机悬挂接头结构疲劳试验的设计方法 |
CN106709183A (zh) * | 2016-12-26 | 2017-05-24 | 中国航空工业集团公司西安飞机设计研究所 | 采用金属板替换的复合材料混杂接头的细节应力仿真方法 |
CN107766612A (zh) * | 2017-09-08 | 2018-03-06 | 中国飞行试验研究院 | 一种连接翼结构形式机翼载荷测量方法 |
CN108090260A (zh) * | 2017-12-01 | 2018-05-29 | 中国直升机设计研究所 | 一种接头约束载荷的分析方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6052880B2 (ja) * | 2013-02-27 | 2016-12-27 | 大和製衡株式会社 | ロードセル |
CN103278272A (zh) * | 2013-05-31 | 2013-09-04 | 哈尔滨工业大学 | 基于激光散斑技术的焊接接头应力集中测试的装置及应力集中测试方法 |
CN103455671B (zh) * | 2013-08-27 | 2016-02-17 | 西北工业大学 | 电磁铆接接头疲劳寿命预测方法 |
JP6237222B2 (ja) * | 2013-12-26 | 2017-11-29 | アイシン精機株式会社 | 荷重検出センサ用カバー及び荷重検出装置 |
-
2018
- 2018-10-17 CN CN201811205768.7A patent/CN109163835B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103196591B (zh) * | 2013-03-07 | 2015-10-21 | 同济大学 | 一种基于正则化和奇异值分解的结构载荷识别方法 |
CN104536941A (zh) * | 2015-01-15 | 2015-04-22 | 西安交通大学 | 一种基于Tikhonov正则化的频域载荷识别方法 |
CN104794299A (zh) * | 2015-04-29 | 2015-07-22 | 西北工业大学 | 一种复合材料干涉配合接头应力分布计算方法 |
CN105043728A (zh) * | 2015-06-23 | 2015-11-11 | 中国航空工业集团公司西安飞机设计研究所 | 一种计算待测接头的拉伸承载能力的方法 |
CN105975704A (zh) * | 2016-05-13 | 2016-09-28 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机悬挂接头结构疲劳试验的设计方法 |
CN106709183A (zh) * | 2016-12-26 | 2017-05-24 | 中国航空工业集团公司西安飞机设计研究所 | 采用金属板替换的复合材料混杂接头的细节应力仿真方法 |
CN107766612A (zh) * | 2017-09-08 | 2018-03-06 | 中国飞行试验研究院 | 一种连接翼结构形式机翼载荷测量方法 |
CN108090260A (zh) * | 2017-12-01 | 2018-05-29 | 中国直升机设计研究所 | 一种接头约束载荷的分析方法 |
Non-Patent Citations (1)
Title |
---|
《某型飞机翼身主接头载荷分配研究》;孙滢;《贵州工业大学学报(自然科学版)》;20051031;第100-104页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109163835A (zh) | 2019-01-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109163835B (zh) | 基于应变识别多接头的接头载荷测量方法 | |
WO2020192623A2 (zh) | 一种旋转叶片非接触式动应变场测量方法及其系统 | |
Glaser et al. | Shape monitoring of a beam structure from measured strain or curvature | |
Nakamura et al. | Inverse identification of continuously distributed loads using strain data | |
CN108254153B (zh) | 光纤气动力测量天平温度补偿方法 | |
Shen et al. | Modeling and analysis of the electrical resistance measurement of carbon fiber polymer–matrix composites | |
DE112018000571T5 (de) | Temperaturkompensation eines kraft-/drehmoment-sensors | |
WO2019032470A1 (en) | AXIAL FORCE PRESSURE TRANSDUCER | |
TWI438450B (zh) | Correction method of measurement error and electronic component characteristic measuring device | |
US7620507B2 (en) | Impedance-estimation methods, modeling methods, articles of manufacture, impedance-modeling devices, and estimated-impedance monitoring systems | |
US20210278330A1 (en) | Quantification of Shear Stress Field Using Tuft Visualization and Deep Learning | |
Iriarte et al. | Optimal strain-gauge placement for mechanical load estimation in circular cross-section shafts | |
CN111232239B (zh) | 曲面挠变位移场重构方法、装置及设备 | |
KR20160097524A (ko) | 동특성 분석에 의한 사장교 케이블의 손상추정방법 | |
CN111220346A (zh) | 一种压电天平校准与使用不确定度的评估方法 | |
CN107766612B (zh) | 一种连接翼结构形式机翼载荷测量方法 | |
Zhao et al. | The application research of inverse finite element method for frame deformation estimation | |
US9274015B2 (en) | Device and method for measuring sectional forces | |
BR102017023377B1 (pt) | Estrutura de compósito, e, sistema e método para medir tensão em uma estrutura de compósito | |
RU2307317C1 (ru) | Косвенный способ настройки тензорезисторных датчиков с мостовой измерительной цепью по аддитивной температурной погрешности | |
CN111380476B (zh) | 一种基于应变测量数据的梁式结构变形测量方法和装置 | |
CN114486302B (zh) | 考虑纵向加载附加弯矩的车钩力测量方法及系统 | |
Rus et al. | Boundary integral equation for inclusion and cavity shape sensitivity in harmonic elastodynamics | |
Liu et al. | Three components strain-gauge type aircraft surface friction resistance sensor: Design, manufacturing, and calibration | |
CN106802205B (zh) | 一种消除剖面内载荷影响的飞行载荷测量方法 |
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 | ||
CP02 | Change in the address of a patent holder |
Address after: Nanchang high tech Industrial Development Zone, Jiangxi Province Patentee after: JIANGXI HONGDU AVIATION INDUSTRY GROUP Co.,Ltd. Address before: 330000 Jiangxi city in Nanchang Province, the new bridge box 460 box 5001 Patentee before: JIANGXI HONGDU AVIATION INDUSTRY GROUP Co.,Ltd. |
|
CP02 | Change in the address of a patent holder |