CN109163835A - 基于应变识别多接头的接头载荷测量方法 - Google Patents
基于应变识别多接头的接头载荷测量方法 Download PDFInfo
- Publication number
- CN109163835A CN109163835A CN201811205768.7A CN201811205768A CN109163835A CN 109163835 A CN109163835 A CN 109163835A CN 201811205768 A CN201811205768 A CN 201811205768A CN 109163835 A CN109163835 A CN 109163835A
- Authority
- CN
- China
- Prior art keywords
- load
- strain
- joint load
- transfer coefficient
- polylinker
- 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
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;
假定每一个接头载荷分别为:F′1,F′2,…,F′3m,因结构在线性范围内,根据线性叠加原理具有如下关系:
于是可得到:
简记为:
{ε′j}=[kij]{F′i} (5)
因此,得到关于载荷和应变之间的一个过定义线性方程组式(5),再采用奇异值分解求解该过定义线性方程组式(5)的解:
其中:若传递系数矩阵[kij]是病态,需增加贴片点,再一次标定传递系数矩阵[kij],并检查传递系数矩阵的条件数;
若始终无法获得一个良态的传递系数矩阵[kij],或传递系数矩阵[kij]本质就是一个病态传递系数矩阵,对于基于测量应变反演识别接头载荷问题,常见的情况是问题的本质是一个病态问题,因此传递系数矩阵是病态,此时不能直接求解,故将病态矩阵的正则化方法引入至基于测量应变反演识别接头载荷问题中,于是通过下式反演计算:
λ为截断的奇异值的阶数,检查和σi之间的衰减关系,若衰减比σi快,即满足Picard条件,那么使用L-curve法选取截断点λ;若不满足Picard条件,通过对不同类型病态问题的大量数值试验,通过:
确定截断的奇异值阶数,将取得更好的效果,此时Δ一般取值0.02,具体算法如下:
第一步,对传递系数矩阵[kij]作奇异值分解,并检查传递系数矩阵[kij]的条件数;
第二步,计算
第三步,取绝对值并找出
第四步,找到满足最大的下标,记为λ;
第五步,计算
有益效果:本发明通过测量应变以识别不同结构件之间的多个连接铰点的载荷,尤其是飞机机翼接头的连接载荷,通过测量标定应变与接头载荷之间的传递关系,而后反演具有过约束系统的接头载荷,从而实现对结构内部接头载荷的测量,有效克服传统通过力传感器来测量载荷的不足,测量出接头载荷对结构的设计优化有着重要意义。
附图说明
图1是本发明的较佳实施例中的具有m个接头的任意弹性体标定载荷示意图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体图示,进一步阐述本发明。
基于应变识别多接头的接头载荷测量方法,具体步骤如下:
1)标定载荷,首先拆开对接的结构件,而后选择其中一个结构件标定,如图1所示,选择第一个接头A1沿x方向逐级加载,一共加载8级,记录载荷F1与每一个应变εi1的数值,线性拟合出载荷F1与每一个应变εi1之间的关系直线,得到3m个斜率ki1,于是:
对于每一个接头任一个方向都按照相同的方式标定,形成应变与载荷之间的传递系数矩阵
标定传递系数矩阵需要注意,由于是反演识别载荷,故应变片的数目要远远大于需要测量的接头载荷数目,只有这样才能够获得比较好的反演解;
2)将两个部件,比如机翼机身,按照实际情况连接,然后在机翼上加载,测量应变并记为ε′1,ε′2,…,ε′3n;
假定每一个接头载荷分别为:F′1,F′2,…,F′3m,因结构在线性范围内,根据线性叠加原理具有如下关系:
于是可得到:
简记为:
{ε′j}=[kij]{F′i} (5)
因此,得到关于载荷和应变之间的一个过定义线性方程组式(5),
再采用奇异值分解求解该过定义线性方程组式(5)的解:
其中:若传递系数矩阵[kij]是病态,需增加贴片点,再一次标定传递系数矩阵[kij],并检查传递系数矩阵的条件数;
若始终无法获得一个良态的传递系数矩阵[kij],或问题的本质就是一个病态问题,对于基于测量应变反演识别接头载荷问题,常见的情况是问题的本质是一个病态问题,因此传递系数矩阵是病态,此时不能直接求解,故将病态矩阵的正则化方法引入至基于测量应变反演识别接头载荷问题中,于是通过下式反演计算:
λ为截断的奇异值的阶数,检查和σi之间的衰减关系,若衰减比σi快,即满足Picard条件,那么使用L-curve法选取截断点λ;若不满足Picard条件,通过对不同类型病态问题的大量数值试验,通过:
确定截断的奇异值阶数,会取得更好的效果,此时Δ一般取值0.02,具体算法如下:
第一步,对传递系数矩阵[kij]作奇异值分解,并检查传递系数矩阵[kij]的条件数;
第二步,计算
第三步,取绝对值并找出
第四步,找到满足最大的下标,记为λ;
第五步,计算
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (8)
1.基于应变识别多接头的接头载荷测量方法,其特征在于,具体步骤如下:
1)标定载荷,首先拆开对接的结构件,而后选择其中一个结构件标定,再将第一个接头A1沿x方向逐级加载,记录载荷F1与每一个应变εi1的数值,而后得到3m个斜率ki1,于是:
对于每一个接头任一个方向都按照相同的方式标定,形成应变与载荷之间的传递系数矩阵
2)将两个部件,按照实际情况连接,在第二个部件上加载,测量应变并记为ε′1,ε′2,…,ε′3n;
假定每一个接头载荷分别为:F1',F′2,…,F′3m,因结构在线性范围内,根据线性叠加原理具有如下关系:
于是得到:
简记为:
{ε′j}=[kij]{F′i} (5)
因此,得到关于载荷和应变之间的一个过定义线性方程组式(5),
再采用奇异值分解求解该过定义线性方程组式(5)的解:
其中:
若始终无法获得一个良态的传递系数矩阵[kij],或传递系数矩阵[kij]本质就是一个病态矩阵,此时将病态矩阵的正则化方法引入至基于测量应变反演识别接头载荷问题中,于是通过下式反演计算:
λ为截断的奇异值的阶数,检查和σi之间的衰减关系,若衰减比σi快,即满足Picard条件,选取截断点λ,从而实现对结构内部接头载荷的测量。
2.根据权利要求1所述的基于应变识别多接头的接头载荷测量方法,其特征在于,步骤1)中,标定传递系数矩阵中,应变片的数目远大于需要测量的接头载荷数目。
3.根据权利要求1所述的基于应变识别多接头的接头载荷测量方法,其特征在于,步骤1)中,3m个斜率ki1由线性拟合出载荷F1与每一个应变εi1之间的关系直线得到。
4.根据权利要求1所述的基于应变识别多接头的接头载荷测量方法,其特征在于,步骤2)中,若传递系数矩阵[kij]是病态,需增加贴片点。
5.根据权利要求4所述的基于应变识别多接头的接头载荷测量方法,其特征在于,再一次标定增加的贴片点传递系数矩阵[kij],并检查传递系数矩阵的条件数。
6.根据权利要求1所述的基于应变识别多接头的接头载荷测量方法,其特征在于,步骤2)中,采用L-curve法选取截断点λ。
7.根据权利要求1所述的基于应变识别多接头的接头载荷测量方法,其特征在于,步骤2)中,若衰减比σi慢,即不满足Picard条件,通过以下:
确定截断的奇异值阶数。
8.根据权利要求7所述的基于应变识别多接头的接头载荷测量方法,其特征在于,Δ取值为0.02,具体算法如下:
第一步,对传递系数矩阵[kij]作奇异值分解,并检查传递系数矩阵[kij]的条件数;
第二步,计算
第三步,取绝对值并找出
第四步,找到满足最大的下标,记为λ;
第五步,计算
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 true CN109163835A (zh) | 2019-01-08 |
CN109163835B 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) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110243564A (zh) * | 2019-06-12 | 2019-09-17 | 南京航空航天大学 | 拦阻钩挂索冲击载荷测试试验装置及试验方法 |
CN112052518A (zh) * | 2020-08-17 | 2020-12-08 | 大连理工大学 | 一种极地船舶结构冰载荷的远场测量识别方法 |
CN112109919A (zh) * | 2020-04-30 | 2020-12-22 | 中国飞机强度研究所 | 一种强度试验加载点布局方法 |
CN112699480A (zh) * | 2020-12-29 | 2021-04-23 | 中国航空工业集团公司西安飞机设计研究所 | 一种动翼面接头布置分析方法 |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103278272A (zh) * | 2013-05-31 | 2013-09-04 | 哈尔滨工业大学 | 基于激光散斑技术的焊接接头应力集中测试的装置及应力集中测试方法 |
CN103455671A (zh) * | 2013-08-27 | 2013-12-18 | 西北工业大学 | 电磁铆接接头疲劳寿命预测方法 |
JP2014163877A (ja) * | 2013-02-27 | 2014-09-08 | Yamato Scale Co Ltd | ロードセル |
CN104536941A (zh) * | 2015-01-15 | 2015-04-22 | 西安交通大学 | 一种基于Tikhonov正则化的频域载荷识别方法 |
US20150185067A1 (en) * | 2013-12-26 | 2015-07-02 | Aisin Seiki Kabushiki Kaisha | Cover for load detection sensor and load detection device |
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 | 中国直升机设计研究所 | 一种接头约束载荷的分析方法 |
-
2018
- 2018-10-17 CN CN201811205768.7A patent/CN109163835B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2014163877A (ja) * | 2013-02-27 | 2014-09-08 | Yamato Scale Co Ltd | ロードセル |
CN103196591B (zh) * | 2013-03-07 | 2015-10-21 | 同济大学 | 一种基于正则化和奇异值分解的结构载荷识别方法 |
CN103278272A (zh) * | 2013-05-31 | 2013-09-04 | 哈尔滨工业大学 | 基于激光散斑技术的焊接接头应力集中测试的装置及应力集中测试方法 |
CN103455671A (zh) * | 2013-08-27 | 2013-12-18 | 西北工业大学 | 电磁铆接接头疲劳寿命预测方法 |
US20150185067A1 (en) * | 2013-12-26 | 2015-07-02 | Aisin Seiki Kabushiki Kaisha | Cover for load detection sensor and load detection device |
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 (2)
Title |
---|
孙滢: "《某型飞机翼身主接头载荷分配研究》", 《贵州工业大学学报(自然科学版)》 * |
阎国良: "《复杂接头机械连接载荷分布工程化算法研究》", 《中国优秀硕士学位论文全文数据库(电子期刊)》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110243564A (zh) * | 2019-06-12 | 2019-09-17 | 南京航空航天大学 | 拦阻钩挂索冲击载荷测试试验装置及试验方法 |
CN112109919A (zh) * | 2020-04-30 | 2020-12-22 | 中国飞机强度研究所 | 一种强度试验加载点布局方法 |
CN112109919B (zh) * | 2020-04-30 | 2024-04-19 | 中国飞机强度研究所 | 一种强度试验加载点布局方法 |
CN112052518A (zh) * | 2020-08-17 | 2020-12-08 | 大连理工大学 | 一种极地船舶结构冰载荷的远场测量识别方法 |
CN112699480A (zh) * | 2020-12-29 | 2021-04-23 | 中国航空工业集团公司西安飞机设计研究所 | 一种动翼面接头布置分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109163835B (zh) | 2020-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109163835A (zh) | 基于应变识别多接头的接头载荷测量方法 | |
Nakamura et al. | Inverse identification of continuously distributed loads using strain data | |
CN103033314B (zh) | 船模波浪载荷试验用测量梁的扭矩标定方法 | |
CN109084687A (zh) | 基于无人机激光雷达测量技术的弧垂计算方法及其装置 | |
CN108195554A (zh) | 六分量光纤气动力测量天平及输出信号组合方法 | |
Moghaddasie et al. | Equilibria and stability boundaries of shallow arches under static loading in a thermal environment | |
CN103018006A (zh) | 船模波浪载荷试验用测量梁的弯矩标定装置及方法 | |
Hu et al. | Performance of wind-excited linked building systems considering the link-induced structural coupling | |
CN109100073A (zh) | 一种基于应变反演的六维力传感器及其测量方法 | |
CN107766612A (zh) | 一种连接翼结构形式机翼载荷测量方法 | |
Afshin et al. | Interlaminar stresses of laminated composite beams resting on elastic foundation subjected to transverse loading | |
US8474326B2 (en) | Load pin with increased performance | |
CN102570452A (zh) | 一种基于模态级数法的facts交互影响程度评估方法 | |
Klump et al. | A new method for finding low-voltage power flow solutions | |
CN103474993B (zh) | 基于映射弹性势能的电网有功承载能力定量分析指标 | |
Tahani et al. | Three-dimensional interlaminar stress analysis at free edges of general cross-ply composite laminates | |
CN105277130A (zh) | 折叠展开结构褶皱区域轴向变形光纤快速计算方法 | |
Samy et al. | Unmanned air vehicle air data estimation using a matrix of pressure sensors: a comparison of neural networks and look-up tables | |
US20100250212A1 (en) | Computer-aided method for a cost-optimized calculation of variable distributions over an aerodynamic surface | |
CN107701379A (zh) | 用于风力发电机组叶片的传感器标定方法和装置 | |
CN205748253U (zh) | 一种盾构机刀具磨损检测装置及盾构机刀具系统 | |
Gangone et al. | Experimental characterization and diagnostics of the early‐age behavior of a semi‐integral abutment FRP deck bridge | |
CN106802205B (zh) | 一种消除剖面内载荷影响的飞行载荷测量方法 | |
CN207675407U (zh) | 六分量光纤气动力测量天平 | |
CN109241688A (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 | ||
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. |