CN112699468A - 主翼面后缘结构及避免后缘卷曲及连接破坏强度分析方法 - Google Patents
主翼面后缘结构及避免后缘卷曲及连接破坏强度分析方法 Download PDFInfo
- Publication number
- CN112699468A CN112699468A CN202011602082.9A CN202011602082A CN112699468A CN 112699468 A CN112699468 A CN 112699468A CN 202011602082 A CN202011602082 A CN 202011602082A CN 112699468 A CN112699468 A CN 112699468A
- Authority
- CN
- China
- Prior art keywords
- stress
- panel
- trailing edge
- sealing structure
- outer panel
- 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
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 13
- 238000007789 sealing Methods 0.000 claims abstract description 36
- 238000000034 method Methods 0.000 claims abstract description 26
- 239000002131 composite material Substances 0.000 claims abstract description 16
- 239000011162 core material Substances 0.000 claims description 28
- 238000009826 distribution Methods 0.000 claims description 11
- 238000001125 extrusion Methods 0.000 claims description 11
- 238000005457 optimization Methods 0.000 claims description 9
- 239000000463 material Substances 0.000 claims description 8
- 238000010008 shearing Methods 0.000 claims description 5
- 230000003068 static effect Effects 0.000 claims description 4
- 230000008014 freezing Effects 0.000 claims description 2
- 238000007710 freezing Methods 0.000 claims description 2
- 238000004519 manufacturing process Methods 0.000 claims description 2
- 238000005452 bending Methods 0.000 description 8
- 230000006378 damage Effects 0.000 description 6
- 238000003825 pressing Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005381 potential energy Methods 0.000 description 2
- 229920000049 Carbon (fiber) Polymers 0.000 description 1
- 239000004917 carbon fiber Substances 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000002788 crimping Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 1
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/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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及结构强度分析技术领域,特别涉及主翼面后缘结构及避免后缘卷曲及连接破坏强度分析方法。该结构包括:主翼面后缘壁板1、密封结构2、动翼面3;其中,密封结构2包括连接区B和夹芯区A;连接区B与主翼面后缘壁板1机械连接,夹芯区A的顶端与夹芯区B接触或存在间隙;连接区B是复材承压板结构,夹芯区A包括内面板A2、外面板A1和橡胶芯子A3;内面板A2和外面板A1之间填充橡胶芯子A3。
Description
技术领域
本发明涉及结构强度分析技术领域,特别涉及主翼面后缘结构及避免后缘卷曲及连接破坏强度分析方法。
背景技术
后缘结构作为飞机的次承力结构,其结构构型设计主要参照已有型号设计经验。在动翼面重复大行程操纵下,载荷形式复杂,结构屈曲严重,局部连接产生破坏,现有分析主要局限在工程经验方面,已不能满足操纵面型号设计需求。
发明内容
发明目的:适用于大型飞机主翼面后缘密封结构的强度设计,以避免其卷曲强度破坏和连接强度破坏。此发明具有较大应用价值,已在某型号主翼面后缘密封结构设计中使用。
技术方案:
一种主翼面后缘结构,包括:主翼面后缘壁板1、密封结构2、动翼面3;
其中,密封结构2包括连接区B和夹芯区A;连接区B与主翼面后缘壁板1机械连接,夹芯区A的顶端与夹芯区B接触或存在间隙;连接区B是复材承压板结构,夹芯区A包括内面板A2、外面板A1和橡胶芯子A3;内面板A2和外面板A1之间填充橡胶芯子A3。
进一步的,橡胶芯子A3为蜂窝结构,其分布厚度与总体外形布置相关。
一种避免后缘卷曲及连接破坏强度分析方法,包括:
对密封结构进行静态应力分析,得到应力以及变形分布趋势;
依据应力以及变形分布趋势,对密封结构的内外面板在预设约束条件下进行结构参数优化,得到内外面板的最优结构参数;预设约束条件包括:体积比、结构最大变形、制造工艺。
进一步的,对密封结构进行静态应力分析,得到应力以及变形分布趋势,包括:
依据总体外形、内外面板的最大厚度,确定橡胶芯子的最小厚度;将橡胶芯子的最小厚度和内外面板的最大厚度作为密封结构的初始结构参数;
采用“组合元素法”对具有初始结构参数的密封结构进行有限元建模;
将内外面板材料的表征、橡胶芯子材料特性的表征输入模型;
将内外面板、内外面板三部分的刚度矩阵进行叠加,并进行有限元线性迭代求解,得到其应力云图及变形云图。
进一步的,依据应力以及变形分布趋势,对密封结构的内外面板在预设约束条件下进行结构参数优化,得到内外面板的最优结构参数,包括:
以密封结构重量最小化为目标函数,以预设的体积比、结构最大变形、面板最大/最小厚度尺寸限制作为约束条件,冻结模型的橡胶芯子区域的优化,对模型的内外面板进行结构参数优化。
进一步的,得到内外面板的最优结构参数之后,所述方法还包括:
采用“组合元素法”对具有初始优化后的结构参数的密封结构进行有限元建模,得到应力数值;
根据应力数值中最大值计算螺栓剪切载荷;
根据螺栓剪切载荷和最优结构参数,计算内外面板钉孔挤压应力;
根据内外面板钉孔挤压应力和许用挤压应力,确定内外面板钉孔挤压强度的安全裕度。
进一步的,采用“组合元素法”对具有初始结构参数的密封结构进行有限元建模,包括:
将内外面板分别等效为二维各向异性层压板单元,单元选取弯曲板元S4R,网格尺寸10mm;将橡胶芯子等效为各向异性三维体单元,单元选取弯曲板元C3D8R,网格尺寸10mm。
进一步的,夹芯区的厚度是决定卷曲破坏和连接破坏的关键参数,其中,厚度越大,其结构复合应力越低。
进一步的,连接区运选用大弹性模量的复合材料,以有效提高连接区临界连接应力,以防止连接区发生连接破坏。
有益效果:
创造性地对后缘密封结构的截面形式、密封橡胶弹性、压紧干涉量进行结构再次优化设计,有效解决了主翼面后缘与动翼面之间的卷曲干涉破坏和连接强度破坏,保证了飞机安全。
附图说明
图1为后缘密封结构受力图;
图2为本发明夹芯区结构平面内受弯示意图;
图3为本发明后缘密封结构有限元网格图;
图4为橡胶的典型应力-应变曲线;
图5a为后缘密封结构应力图;
图5b为后缘密封结构变形图;
图6为优化后内外面板结构参数的示意图;
图7为本发明提供的机翼后缘密封结构的示意图。
具体实施方式
已知某大型飞机主翼面后缘结构,主要有主翼面后缘壁板、密封结构组成,后缘壁板与复材密封结构通过连接区机械连接,此后缘结构主要保证动翼面运动过程中密封性能。由于动翼面运动过程中后缘结构产生复材密封结构的卷曲破坏及连接区的连接破坏,这使得后缘复材密封结构的优化设计成为必然。
a)对密封结构进行工程受力分析如图1所示:
A端为密封结构与动翼面密封接触区,B段为密封结构胡连接区,A段承受动翼面结构的压紧力F、压紧力F产生的摩擦力Q=μF及摩擦力对A段偏心力矩M1=Qδ/2。
式中μ为摩擦系数,δ为A段厚度。
一定的壁板宽度,F一定情况下增加A段厚度δ,则有效降低A段的复合应力,有效防止A段卷曲破坏。
由于力矩M2影响占主,提高材料刚度及高度有效降低B端复材应变。
b)对B段进行结构选型优化
基于步骤(1),增加B段厚度δ,有效降低B段复合应力,并运用大弹性模量的复合材料结构,有效提高B段临界连接应力σxcr,以防止B段发生连接破坏。
式中,D11、D12、D22、D66是指层压板的弯曲刚度系数,m板长度方向的半波数,a、b、t是指板的长度、宽度和厚度。
c)对密封结构A段进行有限元分析并对结构参数进行优化:
为有效承载A段弯矩,降低其弯曲应力,将A段结构设计为夹芯结构,作为夹芯区,如图2所示:由外面板A1、内面板A2和橡胶芯子A3组成。弯曲载荷作用下,内外面板以拉压载荷形式平衡整体弯矩。
(1)采用“组合元素法”的进行有限元建模:
该方法的总体建模思路是将蜂窝夹层结构外面板、中间蜂窝芯子、内面板分别进行建模,将内外面板分别等效为二维各向异性层压板单元,单元选取弯曲板元S4R,网格尺寸10mm;将橡胶芯子等效为各向异性三维体单元,单元选取弯曲板元C3D8R,网格尺寸10mm,有限元网格见图3;
(2)合理表征内外面板、橡胶芯子材料特性:
上下面板为碳纤维复合材料,采用正交各项层压板表述;橡胶芯子为一种高度超弹性材料(非线性弹性),见图4,其应力应变关系采用应变势能U来表达超弹性材料的应力-应变关系,多项式应变势能U表达为,Cij描述肋材料的剪切性能,Di代表肋材料的可压缩性;
(3)矩阵迭代求解:
将三部分的刚度矩阵进行叠加,得到组合元素的刚度矩阵K=K1+K2+K3,再进行有限元线性迭代求解,得到其应力云图及变形云图,见图5a和图5b;
式中,K1、K2、K3指外面板、中间蜂窝芯子、内面板刚度。
(4)优化内外面板铺层参数,得到满足约束条件的结构参数:
内外面板为主要承力件,对其进行结构参数优化;以结构重量最小化为目标函数,以体积比(体积比0.1)和结构最大变形(最大变形13mm)作为约束条件,冻结橡胶芯子区域优化,考虑加工工艺影响,设置最大/最小原件尺寸限制(单元最小尺寸1.296mm、最大尺寸3.5mm),数学表述如下;
u≤uf
δ1≤δ≤δ2
式中,wij(x1,…,xk)、V0、Vf、uf、δ1、δ、δ2是指重量函数、优化前体积、优化后体积、优化后变形数值、优化最小厚度、优化厚度、优化最大厚度。
其优化后内外面板结构参数见图6。
d)对图6结构参数进行连接强度计算,得到满足连接强度最优结构参数,具体结构见图7。
内外面板钉孔挤压应力为:
内外面板钉孔挤压强度的安全裕度为:
式中,α、Q、L、n、D、t、[σbr]是指钉载分配系数、面板剪流、钉间距、钉排数、钉孔直径、面板厚度、面板许用挤压应力。
Claims (9)
1.一种主翼面后缘结构,其特征在于,包括:主翼面后缘壁板(1)、密封结构(2)、动翼面(3);
其中,密封结构(2)包括连接区(B)和夹芯区(A);连接区(B)与主翼面后缘壁板(1)机械连接,夹芯区(A)的顶端与夹芯区(B)接触或存在间隙;连接区(B)是复材承压板结构,夹芯区(A)包括内面板(A2)、外面板(A1)和橡胶芯子(A3);内面板(A2)和外面板(A1)之间填充橡胶芯子(A3)。
2.根据权利要求1所述的结构,其特征在于,橡胶芯子(A3)为蜂窝结构,其分布厚度与总体外形布置相关。
3.一种避免后缘卷曲及连接破坏强度分析方法,其特征在于,包括:
对密封结构进行静态应力分析,得到应力以及变形分布趋势;
依据应力以及变形分布趋势,对密封结构的内外面板在预设约束条件下进行结构参数优化,得到内外面板的最优结构参数;预设约束条件包括:体积比、结构最大变形、制造工艺。
4.根据权利要求3所述的方法,其特征在于,对密封结构进行静态应力分析,得到应力以及变形分布趋势,包括:
依据总体外形、内外面板的最大厚度,确定橡胶芯子的最小厚度;将橡胶芯子的最小厚度和内外面板的最大厚度作为密封结构的初始结构参数;
采用“组合元素法”对具有初始结构参数的密封结构进行有限元建模;
将内外面板材料的表征、橡胶芯子材料特性的表征输入模型;
将内外面板、内外面板三部分的刚度矩阵进行叠加,并进行有限元线性迭代求解,得到其应力云图及变形云图。
5.根据权利要求4所述的方法,其特征在于,依据应力以及变形分布趋势,对密封结构的内外面板在预设约束条件下进行结构参数优化,得到内外面板的最优结构参数,包括:
以密封结构结构重量最小化为目标函数,以预设的体积比、结构最大变形、面板最大/最小厚度尺寸限制作为约束条件,冻结模型的橡胶芯子区域的优化,对模型的内外面板进行结构参数优化。
6.根据权利要求5所述的方法,其特征在于,得到内外面板的最优结构参数之后,所述方法还包括:
采用“组合元素法”对具有初始优化后的结构参数的密封结构进行有限元建模,得到应力数值;
根据应力数值中最大值计算螺栓剪切载荷;
根据螺栓剪切载荷和最优结构参数,计算内外面板钉孔挤压应力;
根据内外面板钉孔挤压应力和许用挤压应力,确定内外面板钉孔挤压强度的安全裕度。
7.根据权利要求4所述的方法,其特征在于,采用“组合元素法”对具有初始结构参数的密封结构进行有限元建模,包括:
将内外面板分别等效为二维各向异性层压板单元,单元选取弯曲板元S4R,网格尺寸10mm;将橡胶芯子等效为各向异性三维体单元,单元选取弯曲板元C3D8R,网格尺寸10mm。
8.根据权利要求3所述的方法,其特征在于,夹芯区的厚度是决定卷曲破坏和连接破坏的关键参数,其中,厚度越大,其结构复合应力越低。
9.根据权利要求3所述的方法,其特征在于,连接区运选用大弹性模量的复合材料,以有效提高连接区临界连接应力,以防止连接区发生连接破坏。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011602082.9A CN112699468A (zh) | 2020-12-29 | 2020-12-29 | 主翼面后缘结构及避免后缘卷曲及连接破坏强度分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011602082.9A CN112699468A (zh) | 2020-12-29 | 2020-12-29 | 主翼面后缘结构及避免后缘卷曲及连接破坏强度分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN112699468A true CN112699468A (zh) | 2021-04-23 |
Family
ID=75512144
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011602082.9A Pending CN112699468A (zh) | 2020-12-29 | 2020-12-29 | 主翼面后缘结构及避免后缘卷曲及连接破坏强度分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112699468A (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105678015A (zh) * | 2016-02-04 | 2016-06-15 | 北京航空航天大学 | 一种高超声速三维机翼的非概率可靠性气动结构耦合优化设计方法 |
US20170129596A1 (en) * | 2015-04-21 | 2017-05-11 | Agustawestland Limited | Helicopter aerofoil with trailing edge flap |
CN110955938A (zh) * | 2019-11-06 | 2020-04-03 | 华中科技大学 | 一种具有梯度多孔夹芯的夹层结构拓扑优化方法 |
CN111159943A (zh) * | 2019-12-25 | 2020-05-15 | 中国航空工业集团公司西安飞机设计研究所 | 一种动翼面封严结构的屈曲处理方法 |
-
2020
- 2020-12-29 CN CN202011602082.9A patent/CN112699468A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170129596A1 (en) * | 2015-04-21 | 2017-05-11 | Agustawestland Limited | Helicopter aerofoil with trailing edge flap |
CN105678015A (zh) * | 2016-02-04 | 2016-06-15 | 北京航空航天大学 | 一种高超声速三维机翼的非概率可靠性气动结构耦合优化设计方法 |
CN110955938A (zh) * | 2019-11-06 | 2020-04-03 | 华中科技大学 | 一种具有梯度多孔夹芯的夹层结构拓扑优化方法 |
CN111159943A (zh) * | 2019-12-25 | 2020-05-15 | 中国航空工业集团公司西安飞机设计研究所 | 一种动翼面封严结构的屈曲处理方法 |
Non-Patent Citations (4)
Title |
---|
刘混举: "机械可靠性设计", 31 January 2009, 国防工业出版社, pages: 137 * |
李新德: "工程力学", 31 July 2006, 中国商业出版社, pages: 189 - 195 * |
范瑞娟;冯翃翎;魏洪: "基于工程应用的蜂窝夹层结构有限元建模方法", 航空计算技术, vol. 46, no. 3, pages 46 - 53 * |
郭学伟: "复合材料无人机机翼结构的轻量化优化研究", 中国优秀硕士学位论文全文数据库工程科技Ⅱ辑, no. 11, pages 36 - 60 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108153981B (zh) | 一种基于有限元分析的机身加筋壁板后屈曲分析方法 | |
Dahmani et al. | Crack identification in reinforced concrete beams using ANSYS software | |
Ovesy et al. | Compressional stability behavior of composite plates with through-the-width and embedded delaminations by using first order shear deformation theory | |
CN108345742B (zh) | 一种适用于波纹夹芯结构的动力学等效建模方法 | |
CN110837690B (zh) | 蜂窝结构非线性本构关系的建立方法、介质和设备 | |
CN109684724B (zh) | 一种客车多材料车身粘铆复合钢骨架结构分析方法 | |
Zhang et al. | Investigation of bird-strike resistance of composite sandwich curved plates with lattice/foam cores | |
Lam et al. | Non‐linear analysis of rectangular laminates under end shortening, using shear deformation plate theory | |
CN112699468A (zh) | 主翼面后缘结构及避免后缘卷曲及连接破坏强度分析方法 | |
CN112926244A (zh) | 一种复合材料层合板开孔件极限载荷确定方法 | |
CN113158508A (zh) | 一种复合材料层合板贴补修复后极限载荷确定方法 | |
CN115544833A (zh) | 一种复合材料蜂窝夹芯板冲击损伤修理方案评估方法 | |
Ovesy et al. | Large deflection finite strip analysis of functionally graded plates under pressure loads | |
CN113408024B (zh) | 装配式地下结构注浆式榫槽接头抗弯承载能力计算方法 | |
CN113343519B (zh) | 一种复合材料夹芯板抗压极限载荷计算方法 | |
Tong et al. | Adhesive elements for stress analysis of bonded patch to curved thin-walled structures | |
CN113722861A (zh) | 一种复合材料螺栓连接结构的强度和失效模式的预测方法 | |
Mikhasev et al. | Flexural deformations and vibrations of a three-layer beam-strip with a stiff core and soft skins | |
Sun et al. | Progressive failure analysis of laminated plates with delamination | |
CN112711808B (zh) | 一种内嵌式集中承载结构螺栓群载荷确定方法 | |
Saraswathy et al. | Dynamic analysis of honeycomb sandwich beam with multiple debonds | |
Arslan et al. | Elasto-plastic behavior of thermoplastic matrix roller chain link plates reinforced with steel fibers | |
Kim et al. | Global stiffness determination of a cross-ply laminated composite plate with distributed delaminations and matrix cracks and its application | |
Ouinas et al. | Full-width disbonding effect on repaired cracks in structural panels with bonded composite patches | |
Pardo et al. | Finite element analysis of the cross-section of wind turbine blades; a comparison between shell and 2D-solid models |
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 |