CN112711808A - 一种内嵌式集中承载结构螺栓群载荷确定方法 - Google Patents

一种内嵌式集中承载结构螺栓群载荷确定方法 Download PDF

Info

Publication number
CN112711808A
CN112711808A CN202011612800.0A CN202011612800A CN112711808A CN 112711808 A CN112711808 A CN 112711808A CN 202011612800 A CN202011612800 A CN 202011612800A CN 112711808 A CN112711808 A CN 112711808A
Authority
CN
China
Prior art keywords
bolt group
composite material
bolt
metal pulley
composite
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
Application number
CN202011612800.0A
Other languages
English (en)
Other versions
CN112711808B (zh
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.)
Xian Aircraft Design and Research Institute of AVIC
Original Assignee
Xian Aircraft Design and Research Institute of AVIC
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 Xian Aircraft Design and Research Institute of AVIC filed Critical Xian Aircraft Design and Research Institute of AVIC
Priority to CN202011612800.0A priority Critical patent/CN112711808B/zh
Publication of CN112711808A publication Critical patent/CN112711808A/zh
Application granted granted Critical
Publication of CN112711808B publication Critical patent/CN112711808B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • 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
    • 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/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • 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)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明属于结构强度分析技术领域。本发明提供一种内嵌式集中承载结构螺栓群载荷确定方法，考虑结构非线性，依据工程受力分析，真实模拟连接基体刚度及连接刚度，考虑结构大集中力载荷，对连接钉载进行准确非线性模拟，准确表征金属接头、复材基体的应力应变水平，形成内嵌式结构集中承载结构钉载分配分析方法。

Description

一种内嵌式集中承载结构螺栓群载荷确定方法
技术领域
本发明属于结构强度分析技术领域，特别涉及集中承载结构钉群载荷分配及精确应力分析。
背景技术
集中传载结构广泛应用于飞机各个部位，这些部位集中载荷的传递多采用金属接头_复材层板混合多钉连接形式。鉴于飞机设计空间限制，金属接头设计成内嵌于翼面的结构形式。
此结构约束条件刚硬，导致集中载荷较大，结构受力复杂，连接金属接头_复材层板钉载分配受结构刚度影响，往往出现载荷分配不均匀的现象。目前，无论是解析法还是有限元素法，针对多钉连接钉载分配的计算基本局限于工程线弹性理论，此方法分析精确低，分析可靠性差。
发明内容
发明目的：提供一种内嵌式集中承载结构钉群载荷确定方法，考虑结构非线性，依据工程受力分析，真实模拟连接基体刚度及连接刚度，考虑结构大集中力载荷，对连接钉载进行准确非线性模拟，准确表征金属接头、复材基体的应力应变水平，形成内嵌式结构集中承载结构钉载分配分析方法。
本发明的技术方案：提供一种内嵌式集中承载结构螺栓群载荷确定方法，所述集中承载结构包括金属滑轮架1、复合材料梁2、前缘层压板复材蒙皮5、盒段复材夹芯壁板3；金属滑轮架1与两侧的复合材料梁2、前缘层压板复材蒙皮5通过螺栓群4机械连接；复合材料梁2与前缘层压板复材蒙皮5通过螺栓群4机械连接；盒段复材夹芯壁板3位于复合材料梁2、滑轮架1围成的框架内；所述螺栓群载荷确定方法包括：
步骤1：对集中承载结构进行工程受力分析，确定集中承载结构的螺栓群4受力状况分类，将集中承载结构的螺栓群4分为关键螺栓群和一般螺栓群；
步骤2：对关键螺栓群，利用关键螺栓群与金属滑轮架1、复合材料梁2的接触非线性，将不连续的接触边界条件运用带有罚函数的摩擦接触对进行表征；
步骤3：对一般螺栓群，不考虑螺栓连接件的连接刚度，只考虑金属滑轮架1、前缘层压板复材蒙皮5的刚度，建立一般螺栓群的刚性连接关系；
步骤4：螺栓连接件、金属滑轮架1受力过程中产生局部塑性流动；利用螺栓连接件及对应金属滑轮架1的材料非线性，采用非线性塑性段线性化等效处理方法表征螺栓连接件、金属滑轮架1材料的真实塑性流动过程；
步骤5：对集中承载结构进行有限元网格剖分，采用渐进式收敛方法确定集中承载结构的有限元模拟网格尺寸及单元类型；
步骤6：根据表征的摩擦接触对、一般螺栓群的刚性连接关系、真实塑性流动过程、有限元网格，采用劲度切线法迭代求解器进行求解，得到关键螺栓群和普通螺栓群的载荷、金属滑轮架的应力分布、前缘层压板复材蒙皮和复合材料梁2的应变分布。
进一步地，步骤1中，金属滑轮架1承受的集中载荷P1、P2，对复合材料梁2平面的力矩M，由复合材料梁2的上下缘条处的螺栓群实现平衡，受力较大；将金属滑轮架1与两侧复合材料梁2上下缘条处的螺栓群定义为关键螺栓群；定义除关键螺栓群外的其他螺栓群为一般螺栓群。
进一步地，步骤2中，基于螺栓连接件与金属滑轮架1、复合材料梁2的切向接触行为，采用的罚函数是指利用摩擦系数模拟所述切向接触行为。
进一步地，步骤2中，摩擦接触对包括：螺栓头与复合材料梁、螺母与金属滑轮架、螺柱与复合材料梁、螺柱与金属滑轮架、金属滑轮架与复合材料梁的接触对。
进一步地，步骤3中，一般螺栓群的刚性连接关系的建立包括：前缘层压板复材蒙皮5与金属滑轮架1、金属滑轮架1与复合材料梁2的腹板的刚性连接关系。
进一步地，步骤5中，有限元网格剖分包括：关键螺栓群及对应的金属滑轮架部位采用2次四面体单元进行剖分，单元网格尺寸为3-5mm；关键螺栓群及对应的复合材料梁2腹板部位采用4-6mm连续体壳元；一般螺栓群对应的复合材料梁、前缘层压板复材蒙皮5部位采用6-7mm的弯曲板元；一般螺栓群对应的金属滑轮架1部位采用12-14mm一次四面体单元；盒段复材夹芯壁板3的芯子采用8-9mm的一次六面体单元。
进一步地，步骤6中，劲度切线法是指：将非线性迭代求解转化为若干个分段线性化的求解。
进一步地，步骤2中，采用非线性塑性段线性化等效处理方法包括：对螺栓连接件及对应的金属滑轮架部位的非线性塑性段采用若干个渐进逼近点进行无限逼近。
本发明的技术效果：本文旨在考虑结构非线性，依据工程受力分析，真实模拟连接基体刚度及连接刚度，考虑结构大集中力载荷，对连接钉载进行准确非线性模拟，准确表征金属接头_复材基体的应力应变水平，形成一整内嵌式结构集中承载结构钉载分配分析方法。
附图说明：
图1为内嵌式集中传载结构图；
图2为工程计算钉载分析图；
图3为非线性硬化应力应变曲线；
图4为螺栓与滑轮架、梁接触示意图；
图5为切向接触行为示意图；
图6为普通区域连接单元示意图；
图7为有限元网格剖分图(左图为整体网格图，右图为A区网格图)；
图8为劲度切线法迭代求解；
图9金属滑轮架应力云图；
图10为复材层板应变云图。
具体实施方式
已知某型飞机内嵌式集中传载结构，由金属滑轮架1、复合材料梁2、前缘层压板复材蒙皮5、盒段复材夹芯壁板3组成，金属滑轮架与周围复材基体结构由螺栓群4连接，见图1。此结构传载为钉群传载，基体刚度及连接刚度决定了钉载分配。本实施例，提供一种内嵌式集中承载结构螺栓群载荷确定方法，具体包括以下步骤：
步骤1：对图1结构进行工程受力分析，确定集中传载结构的钉群分类：滑轮架主要承受集中载荷P1、集中载荷P2，其受力工程图见图2，滑轮架与梁缘条连接钉、滑轮架与梁腹板连接钉为关键螺栓群；
梁缘条连接的螺栓承受的剪力：
Figure BDA0002869835510000041
梁腹板连接的螺栓承受的剪力：Q＝P1-P2。
步骤2：基于步骤(1)依据受载情况对连接钉群分类：A区为关键受力钉，滑轮架其他连接区为基本传力钉，分区见图1。
步骤3：初步工程分析，静力分析中考虑材料非线性硬化：A区金属滑轮架进入塑性流动，图3非线性材料硬化段进行线性化分解近似处理，输入相应的塑性数据，以便对应力水平、载荷传递路径进行精确模拟。
步骤4：对A区关键钉建立非线性接触弹塑性接触模型：每个螺栓建立螺栓头与梁、螺母与复材板、螺柱与梁、螺柱与复材板、滑轮架与复材板共计5对弹塑性面接触对，见图4，将不连续的接触边界条件运用带有罚函数的摩擦切向接触对进行表征，见图5。
步骤5：一般螺栓群建立钉载连接单元：不考虑螺栓连接件的连接刚度，只考虑金属滑轮架、前缘层压板复材蒙皮的刚度建立滑轮架与复材层板结构的连接单元，连接单元属性为Hinge(释放绕转轴方向自由度)，并于连接点创建Faster约束，影响半径为钉直径D，见图6。
步骤6：进行合理的网格剖分采用渐进式收敛方法确定集中承载结构的有限元模拟网格尺寸及单元类型：螺栓、金属滑轮架采用2次四面体单元C3D10M进行剖分，A区网格尺寸大小为3mm，其他部位网格尺寸5mm；A区复材层板采用3mm连续体壳元SC8R以便建立接触模型；其他部位复材层板结构采用5mm普通弯曲板元S4R单元，芯子采用5mm一次六面体单元C3D8R，网格剖分见图7；
步骤7：选用Newton-ruphson劲度切线法迭代求解器进行求解：采用Newton-ruphson将非线性硬化问题转化为一种分段线性化的求解方法，见图8；
步骤8：通过非线性获得钉群载荷、滑轮架应力分布、复材应变分布，见图9、图10。

Claims (8)

1.一种内嵌式集中承载结构螺栓群载荷确定方法，所述集中承载结构包括金属滑轮架(1)、复合材料梁(2)、前缘层压板复材蒙皮(5)、盒段复材夹芯壁板(3)；金属滑轮架(1)与两侧的复合材料梁(2)、前缘层压板复材蒙皮(5)通过螺栓群(4)机械连接；复合材料梁(2)与前缘层压板复材蒙皮(5)通过螺栓群(4)机械连接；盒段复材夹芯壁板(3)位于复合材料梁(2)、滑轮架(1)围成的框架内；其特征在于，所述螺栓群载荷确定方法包括：
步骤1：对集中承载结构进行工程受力分析，确定集中承载结构的螺栓群(4)受力状况分类，将集中承载结构的螺栓群(4)分为关键螺栓群和一般螺栓群；
步骤2：对关键螺栓群，利用关键螺栓群与金属滑轮架(1)、复合材料梁(2)的接触非线性，将不连续的接触边界条件运用带有罚函数的摩擦接触对进行表征；
步骤3：对一般螺栓群，不考虑螺栓连接件的连接刚度，只考虑金属滑轮架(1)、前缘层压板复材蒙皮(5)的刚度，建立一般螺栓群的刚性连接关系；
步骤4：螺栓连接件、金属滑轮架(1)受力过程中产生局部塑性流动；利用螺栓连接件及对应金属滑轮架(1)的材料非线性，采用非线性塑性段线性化等效处理方法表征螺栓连接件、金属滑轮架(1)材料的真实塑性流动过程；
步骤5：对集中承载结构进行有限元网格剖分，采用渐进式收敛方法确定集中承载结构的有限元模拟网格尺寸及单元类型；
步骤6：根据表征的摩擦接触对、一般螺栓群的刚性连接关系、真实塑性流动过程、有限元网格，采用劲度切线法迭代求解器进行求解，得到关键螺栓群和普通螺栓群的载荷、金属滑轮架的应力分布、前缘层压板复材蒙皮和复合材料梁(2)的应变分布。
2.根据权利要求1所述的螺栓群载荷确定方法，其特征在于，步骤1中，
金属滑轮架(1)承受的集中载荷P1、P2，对复合材料梁(2)平面的力矩M，由复合材料梁(2)的上下缘条处的螺栓群实现平衡，受力较大；将金属滑轮架(1)与两侧复合材料梁(2)上下缘条处的螺栓群定义为关键螺栓群；
定义除关键螺栓群外的其他螺栓群为一般螺栓群。
3.根据权利要求2所述的螺栓群载荷确定方法，其特征在于，步骤2中，基于螺栓连接件与金属滑轮架(1)、复合材料梁(2)的切向接触行为，采用的罚函数是指利用摩擦系数模拟所述切向接触行为。
4.根据权利要求3所述的螺栓群载荷确定方法，其特征在于，步骤2中，摩擦接触对包括：螺栓头与复合材料梁、螺母与金属滑轮架、螺柱与复合材料梁、螺柱与金属滑轮架、金属滑轮架与复合材料梁的接触对。
5.根据权利要求4所述的螺栓群载荷确定方法，其特征在于，步骤3中，一般螺栓群的刚性连接关系的建立包括：前缘层压板复材蒙皮(5)与金属滑轮架(1)、金属滑轮架(1)与复合材料梁(2)的腹板的刚性连接关系。
6.根据权利要求5所述的螺栓群载荷确定方法，其特征在于，步骤5中，有限元网格剖分包括：关键螺栓群及对应的金属滑轮架部位采用2次四面体单元进行剖分，单元网格尺寸为3-5mm；关键螺栓群及对应的复合材料梁(2)腹板部位采用4-6mm连续体壳元；一般螺栓群对应的复合材料梁、前缘层压板复材蒙皮(5)部位采用6-7mm的弯曲板元；一般螺栓群对应的金属滑轮架(1)部位采用12-14mm一次四面体单元；盒段复材夹芯壁板(3)的芯子采用8-9mm的一次六面体单元。
7.根据权利要求6所述的螺栓群载荷确定方法，其特征在于，步骤6中，劲度切线法是指：将非线性迭代求解转化为若干个分段线性化的求解。
8.根据权利要求7所述的螺栓群载荷确定方法，其特征在于，步骤2中，采用非线性塑性段线性化等效处理方法包括：对螺栓连接件及对应的金属滑轮架部位的非线性塑性段采用若干个渐进逼近点进行无限逼近。
CN202011612800.0A 2020-12-29 2020-12-29 一种内嵌式集中承载结构螺栓群载荷确定方法 Active CN112711808B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011612800.0A CN112711808B (zh) 2020-12-29 2020-12-29 一种内嵌式集中承载结构螺栓群载荷确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011612800.0A CN112711808B (zh) 2020-12-29 2020-12-29 一种内嵌式集中承载结构螺栓群载荷确定方法

Publications (2)

Publication Number Publication Date
CN112711808A true CN112711808A (zh) 2021-04-27
CN112711808B CN112711808B (zh) 2024-05-24

Family

ID=75547346

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011612800.0A Active CN112711808B (zh) 2020-12-29 2020-12-29 一种内嵌式集中承载结构螺栓群载荷确定方法

Country Status (1)

Country Link
CN (1) CN112711808B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004212210A (ja) * 2002-12-27 2004-07-29 National Aerospace Laboratory Of Japan 締結具の剪断荷重測定方法
CN107273646A (zh) * 2017-08-04 2017-10-20 中国航空工业集团公司西安飞机设计研究所 变温度下金属与复材混合结构中螺栓连接内力计算方法
CN107451365A (zh) * 2017-08-04 2017-12-08 中国航空工业集团公司西安飞机设计研究所 一种集中传载的长螺栓强度及刚度分析方法
CN109726438A (zh) * 2018-12-04 2019-05-07 中国航空工业集团公司西安飞机设计研究所 一种混合结构温度应力计算方法
CN109753686A (zh) * 2018-12-04 2019-05-14 中国航空工业集团公司西安飞机设计研究所 一种集中传载结构钉群载荷分配计算方法
CN111159879A (zh) * 2019-12-25 2020-05-15 中国航空工业集团公司西安飞机设计研究所 一种集中传载结构的钉载处理方法
CN111159943A (zh) * 2019-12-25 2020-05-15 中国航空工业集团公司西安飞机设计研究所 一种动翼面封严结构的屈曲处理方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004212210A (ja) * 2002-12-27 2004-07-29 National Aerospace Laboratory Of Japan 締結具の剪断荷重測定方法
CN107273646A (zh) * 2017-08-04 2017-10-20 中国航空工业集团公司西安飞机设计研究所 变温度下金属与复材混合结构中螺栓连接内力计算方法
CN107451365A (zh) * 2017-08-04 2017-12-08 中国航空工业集团公司西安飞机设计研究所 一种集中传载的长螺栓强度及刚度分析方法
CN109726438A (zh) * 2018-12-04 2019-05-07 中国航空工业集团公司西安飞机设计研究所 一种混合结构温度应力计算方法
CN109753686A (zh) * 2018-12-04 2019-05-14 中国航空工业集团公司西安飞机设计研究所 一种集中传载结构钉群载荷分配计算方法
CN111159879A (zh) * 2019-12-25 2020-05-15 中国航空工业集团公司西安飞机设计研究所 一种集中传载结构的钉载处理方法
CN111159943A (zh) * 2019-12-25 2020-05-15 中国航空工业集团公司西安飞机设计研究所 一种动翼面封严结构的屈曲处理方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
郑茂亮;范瑞娟;李健;魏洪;: "大间隙集中传载结构强度分析方法研究", 航空计算技术, no. 02 *
郑茂亮;范瑞娟;李健;魏洪;: "大间隙集中传载结构强度分析方法研究", 航空计算技术, no. 02, 25 March 2016 (2016-03-25) *

Also Published As

Publication number Publication date
CN112711808B (zh) 2024-05-24

Similar Documents

Publication Publication Date Title
CN108153981B (zh) 一种基于有限元分析的机身加筋壁板后屈曲分析方法
Grover et al. Analytical and finite element modeling of laminated composite and sandwich plates: An assessment of a new shear deformation theory for free vibration response
Overgaard et al. Structural collapse of a wind turbine blade. Part B: Progressive interlaminar failure models
CN107729603A (zh) 复合材料壁板螺栓连接的非线性分析方法
US8255174B2 (en) Method and device for determining critical buffeting loads on a structure of an aircraft
US10824778B2 (en) Computer-implemented method for simulating behavior of a screw connection between metal components of engineering structures
Bhardwaj et al. Numerical simulations of cracked plate using XIGA under different loads and boundary conditions
CN105184390A (zh) 一种壁板结构静强度、刚度、稳定性的综合优化方法
CN111159943B (zh) 一种动翼面封严结构的屈曲处理方法
Hissaria et al. Numerical investigation of eigenvalue characteristics (vibration and buckling) of damaged porous bidirectional FG panels
Demirtaş et al. Free vibration analysis of an aircraft wing by considering as a cantilever beam
Solanki et al. Meshfree approach for linear and nonlinear analysis of sandwich plates: A critical review of twenty plate theories
Sit et al. A third order nonlinear model to study the dynamic behaviour of composite laminated structures under thermal effect with experimental verification
CN111159879B (zh) 一种集中传载结构的钉载处理方法
Lokatt Aeroelastic flutter analysis considering modeling uncertainties
Wu et al. Free vibration analysis and multi-objective optimization of lattice sandwich beams
CN112711808A (zh) 一种内嵌式集中承载结构螺栓群载荷确定方法
AsadiGorgi et al. Effects of all-over part-through cracks on the aeroelastic characteristics of rectangular panels
Jam et al. A new improved high‐order theory for analysis of free vibration of sandwich panels
Beheshti Novel quadrilateral elements based on explicit Hermite polynomials for bending of Kirchhoff–Love plates
CN115391947A (zh) 轨道交通车辆复合材料结构虚拟分析方法及系统
Bayraktar et al. Free vibration analysis of an aircraft wing by considering as a cantilever beam
Murphy et al. Tailoring static strength performance of metallic stiffened panels by selective local sub-stiffening
JP4179095B2 (ja) 積層材料内の層間応力演算方法
Szyndler et al. Numerical and Experimental Investigation of the Innovatory Incremental-Forming Process Dedicated to the Aerospace Industry

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