CN104330253B - 一种加筋壁板损伤容限特性分析方法 - Google Patents
一种加筋壁板损伤容限特性分析方法 Download PDFInfo
- Publication number
- CN104330253B CN104330253B CN201410713040.0A CN201410713040A CN104330253B CN 104330253 B CN104330253 B CN 104330253B CN 201410713040 A CN201410713040 A CN 201410713040A CN 104330253 B CN104330253 B CN 104330253B
- Authority
- CN
- China
- Prior art keywords
- rib
- load
- stiffened panel
- analysis
- factor
- 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
Abstract
本发明涉及一种加筋壁板损伤容限特性分析方法,步骤为:首先选取典型截面位置对加筋壁板结构进行评估,确定影响因素;应用工程方法求取基本加筋影响系数β加筋;应用工程方法求取结构鼓胀影响系数β鼓胀;引入载荷影响系数,采用试验数据修正β载荷;分别计算得到每一个裂纹长度相应的应力强度因子:结合裂纹扩展分析、剩余强度分析,完成整个加筋壁板的损伤容限分析过程。本发明对于飞机加筋壁板的损伤容限特性的分析,提供了一种详细的、准确的、快速的计算方法,填补了目前国内对于该领域在实际应用中的空白。具有实际应用的前景。
Description
技术领域
本发明涉及一种加筋壁板损伤容限特性分析方法,属于结构损伤容限设计领域。
背景技术
在航空设计领域,现代大型运输机大都采用损伤容限设计理念进行设计,而在整个损伤容限分析中,机身、机翼蒙皮壁板损伤容限分析在整个损伤容限分析中占有极其重要地位,但是其分析难度高、工作量大,计算时间久,且由于分析结果的不可靠性,大多需要单独的研发或相关试验进行验证。同时国内在大型运输机设计方面还处于初期阶段,损伤容限分析方面还远未成熟,只能够进行简单的加筋分析。目前国内对于加筋壁板的损伤容限没有提供详细的分析方法或者运算工具,其分析仅仅是从一些设计手册中查找图表,而这些图表只是特定结构的结果,无法考虑实际结构在几何构型、加载情况或者加工工艺等发生的变化,其分析结果的完整性及可靠性远远不足,而从保守的角度设计角度出发,又导致了飞机结构的厚重从而带来了设计成本的增加。
发明内容
本发明的目的是针对现有方法的不足,提供一种快速、准确的加筋壁板损伤容限评估方法。本方法在设计初期及后期均有着广泛的应用空间,可以极大的提高分析的准确性,从损伤容限的角度保证飞机服役安全。
本发明的目的是通过下述技术方案实现的:
一种加筋壁板损伤容限特性分析方法,其特征在于,包括如下步骤:
步骤一:选取典型截面位置对加筋壁板结构进行评估,确定影响因素;
其中,所述的截面位置:包括典型的剖面位置及除开口和连接结构外的具有代表性的蒙皮曲率的结构位置,具体为:机身顶部、底部、左右两侧;影响因素包括,几何特性——即筋条对于蒙皮上裂纹扩展的作用、鼓胀效应及载荷因素;
步骤二:应用工程方法求取基本加筋影响系数β加筋;
具体过程如下:
将加筋板结构分为无筋条的壁板和筋条两部分,其中壁板受到单轴的拉应力即垂直于裂纹方向,以及一系列的钉载,筋条受到的是紧固件载荷的反作用力;该壁板可视为三种构型的叠加,包括:①受到均匀拉应力的有裂纹的板;②施加集中力,即钉载,但不包含裂纹的板;③考虑裂纹表面受到的分布载荷的板;
所述的第①个构型的应力强度因子计算公式为:
第②个构型的应力强度因子计算公式为:
K=0
第③个构型的应力强度因子计算公式为:
其中,
其中,B为板厚,σ为拉应力,a为半裂纹长度,s为筋件间距的一半,y0为筋条上紧固件间距,ν为裂纹区沿载荷方向位移,F为钉载,x为裂纹中心到裂尖距离;式中参数见图4;
当钉载F已知时,基于位移协调法,可以得到筋条的载荷集中系数Ls以及筋条影响系数β加筋,即加筋板的应力强度因子与未加筋板的应力强度因子之比:
其中Fmax为附加紧固件载荷后筋条上最大拉伸载荷,F∞为无裂纹情况下筋条上载荷,As为筋条面积;
筋条上的应力可由下式得到:
σ筋条=LSσ;
F1+F2…+Fj=(Ls-1)*σ*As
其中Fj代表加筋件上裂纹一侧各紧固件载荷,j取值一般小于等于7,筋条的刚度对于β加筋和Ls的数值有较明显的影响;
步骤三:应用工程方法求取结构鼓胀影响系数β鼓胀;
具体过程如下:
鼓胀系数仅仅是在增压情况时有曲率半径的蒙皮壁板上的纵向裂纹有影响,计算过程基于如下公式;
一跨裂纹时的计算公式:
两跨裂纹时筋条完整的计算公式:
两跨裂纹时中间筋条断裂的计算公式:
其中,c:半裂纹长度;R:曲率半径;L:筋条间距;Z:阻尼系数,即,位于筋条位置处的鼓胀与筋条之间全部鼓胀的比值;w:两筋条中点到裂纹尖端的距离;
步骤四:引入载荷影响系数,采用试验数据修正β载荷,
其修正公式为:
β载荷=1.6138*(a/L)6-5.1836*(a/L)5+6.4747*(a/L)4
-3.8982*(a/L)3+0.7567*(a/L)2+0.338*(a/L)+0.8256
其中,a为半裂纹长度,L为筋条间距;
步骤五:根据前述步骤,分别计算得到每一个裂纹长度相应的应力强度因子:结合裂纹扩展分析、剩余强度分析,完成整个加筋壁板的损伤容限分析过程。
本发明的有益的效果:
本发明对于飞机加筋壁板的损伤容限特性的分析,提供了一种详细的、准确的、快速的计算方法,填补了目前国内对于该领域在实际应用中的空白,目前国内对于飞机加筋壁板的损伤容限特性的分析仅仅依靠现有的设计手册中单一特定构型的图表进行查询,而对于实际结构的几何、载荷及加工工艺的变化无法体现,没有系统的对于该结构损伤容限分析的方法及计算工具,而这正是本发明方法具备的特点;本发明方法可以在设计初期进行结构方案评估,避免后期工作中的概念更改;也可在设计的中后期对结构进行损伤容限分析以及试验设计,提高结构安全性以及大量降低试验数量,极大的节约项目设计成本。通过经验公式及试验数据修正影响系数,例如载荷影响系数,可以使得原有工程分析方法更加精确,解决工程方法无法完全反应真实细节,只能将结构模型进一步简化而导致分析结果精度过低的问题;通过引入试验数据修正载荷系数,解决原有工程算法对特殊结构的影响系数无法真实求取的问题,如筋条结构对化铣蒙皮处应力水平的影响,工程方法近似的认为在两个筋条之间的蒙皮应力是不变的,其值为该区域的最大应力水平,而实际结构中,这种假设是保守的,越靠近筋条蒙皮上的应力越小。从而在实际分析中,免去了不必要的保守因素,使得计算结果与实际结果更加吻合;通过工程算法与经验公式及试验数据相结合,大大的提高了分析的精度及分析速度,极大的减少试验件数量,节约设计成本。通过本方法很容易得到典型位置的加筋壁板损伤容限特性的分析结果。而一般民用飞机通常包含大面积的加筋壁板结构,因此对其他位置的相似结构,只需更改几何构型等必要的输入参数,即可应用该分析结果对相似结构的损伤容限特性进行快速、准确的分析。经实践证明,该方法具有实际应用的前景。
附图说明
图1是建立细节模型选取的一组典型位置示意图,图中,包括顶部1L-1R长桁区域、上中部8L-9L长桁区域、与地板梁连接附近位置19L-20L及21L-22L长桁区域以及底部36L-36R长桁区域分别建立细节模型进行分析。
图2选取图1中部分区域进行分析的状态图。
图3是有裂纹的壁板的应力强度因子分解示意图。
图4是筋条对于蒙皮上裂纹扩展的影响的工程算法原理图。
图5是加筋系数随半裂纹尺寸变化曲线图。
图6是膨胀系数随半裂纹尺寸变化曲线图。
图7是应力分布系数随半裂纹尺寸变化曲线图。
图8是应力强度因子随半裂纹尺寸变化曲线图。
图9是修正后的工程算法和试验结果的裂纹扩展曲线对比图。
具体实施方式
以某飞机机身蒙皮壁板的损伤容限分析为例,说明对典型机身蒙皮加筋壁板的损伤容限分析方法的内容和步骤,附图均取自该部段分析计算结果,此分析方法全部采用本发明方法。
下面结合附图,对本发明具体实现过程作进一步详细描述。具体步骤描述过程以某飞机机身壁板纵向裂纹损伤容限分析为例。其基本实现过程如下:
步骤一:选取典型截面位置对加筋壁板结构进行评估,确定影响因素,如图1、2所示;
其中,所述的截面位置:包括典型的剖面位置及除开口和连接结构外的具有代表性的蒙皮曲率的结构位置,具体为:机身顶部、底部、左右两侧;影响因素包括,几何特性——即筋条对于蒙皮上裂纹扩展的作用、鼓胀效应及载荷因素;
步骤二:应用工程方法求取基本加筋影响系数β加筋;
具体过程如下:
一般情况下,受到不均匀载荷的壁板上有裂纹的应力强度因子等于没有裂纹的板上在裂纹处施加局部应力的应力强度因子,即如图3(a)和(b)所示,Ka=Kb;同时这种情况又可以将壁板分为图3(d)和(e)两者应力强度因子叠加的结果,即等把远端载荷和局部载荷分为两个模型,Ka=Kb=Kd+Ke。对于实际加筋壁板来说,可以将加筋板结构分为无筋条的壁板和筋条两部分,其中壁板结构又可以按照上述情况继续细分,考虑到钉载的作用,将其分为三部分,包括:①受到均匀拉应力的有裂纹的板;②施加集中力,即钉载,但不包含裂纹的板;③考虑裂纹表面受到的分布载荷的板。筋条受到的是紧固件载荷的反作用力,如图4所示。
所述的第①个构型的应力强度因子计算公式为:
第②个构型的应力强度因子计算公式为:
K=0
第③个构型的应力强度因子计算公式为:
其中,
其中,B为板厚,σ为拉应力,a为半裂纹长度,s为筋件间距的一半,y0为筋条上紧固件间距,ν为裂纹区沿载荷方向位移,F为钉载,x为裂纹中心到裂尖距离;式中参数见图4;
当钉载F已知时,基于位移协调法,可以得到筋条的载荷集中系数Ls以及筋条影响系数β加筋,即加筋板的应力强度因子与未加筋板的应力强度因子之比:
其中Fmax为附加紧固件载荷后筋条上最大拉伸载荷,F∞为无裂纹情况下筋条上载荷,As为筋条面积;
筋条上的应力可由下式得到:
σ筋条=LSσ;
F1+F2…+Fj=(Ls-1)*σ*As
其中Fj代表加筋件上裂纹一侧各紧固件载荷,j取值一般小于等于7,筋条的刚度对于β加筋和Ls的数值有较明显的影响;
上述数据为该工程算法计算结果,其中第一列数据为裂纹长度,第三列数据为对应的β加筋,其余分别为应力强度因子以及筋条受到的应力及钉载。计算结果曲线如图5所示。
步骤三:应用工程方法求取结构鼓胀影响系数β鼓胀;
具体过程如下:
鼓胀系数仅仅是在增压情况时有曲率半径的蒙皮壁板上的纵向裂纹有影响,计算过程基于如下公式;
一跨裂纹时的计算公式:
两跨裂纹时筋条完整的计算公式:
两跨裂纹时中间筋条断裂的计算公式:
其中,c:半裂纹长度;R:曲率半径;L:筋条间距;Z:阻尼系数,即,位于筋条位置处的鼓胀与筋条之间全部鼓胀的比值;w:两筋条中点到裂纹尖端的距离;
其计算结果如图6所示。
步骤四:引入载荷影响系数,引入试验数据修正β载荷,如图7所示。
其修正公式为:
β载荷=1.6138*(a/L)6-5.1836*(a/L)5+6.4747*(a/L)4
-3.8982*(a/L)3+0.7567*(a/L)2+0.338*(a/L)+0.8256
其中,a为半裂纹长度,L为筋条间距;
步骤五:根据前述步骤,得到修正后的应力强度因子,如图8所示。
通过上述发明方法得到每一个裂纹长度对应的应力强度因子后,即可继续进行裂纹扩展分析以及剩余强度分析,来完成整个结构的损伤容限分析。
本方法考虑了加筋的影响、鼓胀效应以及试验数据修正后的载荷分布,对于壁板结构的损伤容限分析来说,考虑全面且精度很高,可以快速并准确的得到分析结果,与试验结果吻合较好,误差控制在30%以内,如图9所示。
以上所述仅是本发明的实现方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明技术原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (1)
1.一种加筋壁板损伤容限特性分析方法,其特征在于,包括如下步骤:
步骤一:选取典型截面位置对加筋壁板结构进行评估,确定影响因素;
其中,所述的截面位置:包括典型的剖面位置及除开口和连接结构外的具有代表性的蒙皮曲率的结构位置,具体为:机身顶部、底部、左右两侧;影响因素包括,几何特性——即筋条对于蒙皮上裂纹扩展的作用、鼓胀效应及载荷因素;
步骤二:应用工程方法求取基本加筋影响系数β加筋;
具体过程如下:
将加筋板结构分为无筋条的壁板和筋条两部分,其中壁板受到单轴的拉应力即垂直于裂纹方向,以及一系列的钉载,筋条受到的是紧固件载荷的反作用力;该壁板可视为三种构型的叠加,包括:①受到均匀拉应力的有裂纹的板;②施加集中力,即钉载,但不包含裂纹的板;③考虑裂纹表面受到的分布载荷的板;
所述的第①个构型的应力强度因子计算公式为:
第②个构型的应力强度因子计算公式为:
K=0
第③个构型的应力强度因子计算公式为:
其中,
其中,B为板厚,σ为拉应力,a为半裂纹长度,s为筋件间距的一半,y0为筋条上紧固件间距,ν为裂纹区沿载荷方向位移,F为钉载,x为裂纹中心到裂尖距离;
当钉载F已知时,基于位移协调法,可以得到筋条的载荷集中系数Ls以及筋条影响系数β加筋,即加筋板的应力强度因子与未加筋板的应力强度因子之比:
其中Fmax为附加紧固件载荷后筋条上最大拉伸载荷,F∞为无裂纹情况下筋条上载荷,As为筋条面积;
筋条上的应力可由下式得到:
σ筋条=LSσ;
F1+F2…+Fj=(Ls-1)*σ*As
其中Fj代表加筋件上裂纹一侧各紧固件载荷,筋条的刚度对于β加筋和Ls的数值有较明显的影响;
步骤三:应用工程方法求取结构鼓胀影响系数β鼓胀;
具体过程如下:
鼓胀系数仅仅是在增压情况时有曲率半径的蒙皮壁板上的纵向裂纹有影响,计算过程基于如下公式;
一跨裂纹时的计算公式:
两跨裂纹时筋条完整的计算公式:
两跨裂纹时中间筋条断裂的计算公式:
其中,c:半裂纹长度;R:曲率半径;L:筋条间距;Z:阻尼系数,即,位于筋条位置处的鼓胀与筋条之间全部鼓胀的比值;w:两筋条中点到裂纹尖端的距离;
步骤四:引入载荷影响系数,引入试验数据修正β载荷,
其修正公式为:
β载荷=1.6138*(a/L)6-5.1836*(a/L)5+6.4747*(a/L)4
-3.8982*(a/L)3+0.7567*(a/L)2+0.338*(a/L)+0.8256
其中,a为半裂纹长度,L为筋条间距;
步骤五:根据前述步骤,分别计算得到每一个裂纹长度相应的应力强度因子:结合裂纹扩展分析、剩余强度分析,完成整个加筋壁板的损伤容限分析过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410713040.0A CN104330253B (zh) | 2014-11-28 | 2014-11-28 | 一种加筋壁板损伤容限特性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410713040.0A CN104330253B (zh) | 2014-11-28 | 2014-11-28 | 一种加筋壁板损伤容限特性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104330253A CN104330253A (zh) | 2015-02-04 |
CN104330253B true CN104330253B (zh) | 2017-03-08 |
Family
ID=52405023
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410713040.0A Active CN104330253B (zh) | 2014-11-28 | 2014-11-28 | 一种加筋壁板损伤容限特性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104330253B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105183926A (zh) * | 2015-05-07 | 2015-12-23 | 苏州纽威阀门股份有限公司 | 一种阀体优化设计方法 |
CN106354934B (zh) * | 2016-08-30 | 2019-07-02 | 中航沈飞民用飞机有限责任公司 | 一种基于有限元分析求解的加筋壁板损伤容限分析方法 |
CN107356420A (zh) * | 2017-07-28 | 2017-11-17 | 中航沈飞民用飞机有限责任公司 | 一种机身壁板搭接结构广布疲劳损伤分析方法 |
CN108036906B (zh) * | 2017-12-08 | 2019-08-30 | 合肥工业大学 | 一种裂纹转子刚度系数测量方法 |
CN108984909B (zh) * | 2018-07-19 | 2022-11-11 | 中航沈飞民用飞机有限责任公司 | 一种基于Mar-Lin模型的含大损伤飞机复合材料结构剩余强度分析方法 |
CN113109429B (zh) * | 2021-03-22 | 2023-09-19 | 国营芜湖机械厂 | 一种飞机缺陷评估用的对比试块及其使用方法 |
CN113443169B (zh) * | 2021-06-11 | 2022-08-02 | 航空工业第一飞机设计研究院 | 一种整体壁板结构损伤断裂主动控制的试验方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101561459A (zh) * | 2009-06-04 | 2009-10-21 | 中国航空工业集团公司西安飞机设计研究所 | 飞机外翼下壁板损伤容限试验方法 |
CN101561365A (zh) * | 2009-06-04 | 2009-10-21 | 中国航空工业集团公司西安飞机设计研究所 | 飞机整体翼梁损伤容限试验方法 |
CN202141630U (zh) * | 2011-05-04 | 2012-02-08 | 中国航空工业集团公司西安飞机设计研究所 | 一种整体加筋平板损伤容限试验件过渡段 |
GB201108170D0 (en) * | 2011-05-17 | 2011-06-29 | Rolls Royce Plc | Improving the damage tolerance of a rotor assembly |
KR20130125935A (ko) * | 2012-05-10 | 2013-11-20 | 한국항공우주산업 주식회사 | 항공기 구조물의 충격 후 압축 파손 강도를 예측하는 방법 |
CN103870613A (zh) * | 2012-12-10 | 2014-06-18 | 中国飞机强度研究所 | 一种加筋壁板承载能力计算方法 |
CN104034576B (zh) * | 2014-06-16 | 2016-11-02 | 中国飞机强度研究所 | 一种飞机结构裂纹扩展试验用载荷谱简化方法 |
-
2014
- 2014-11-28 CN CN201410713040.0A patent/CN104330253B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN104330253A (zh) | 2015-02-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104330253B (zh) | 一种加筋壁板损伤容限特性分析方法 | |
CN107506529B (zh) | 一种复合材料加筋壁板轴压稳定性计算方法 | |
CN105184390A (zh) | 一种壁板结构静强度、刚度、稳定性的综合优化方法 | |
CN109033526B (zh) | 一种翼肋与蒙皮铆钉连接载荷计算方法 | |
CN108153981A (zh) | 一种基于有限元分析的复合材料机身加筋壁板结构后屈曲分析方法 | |
CN103870613A (zh) | 一种加筋壁板承载能力计算方法 | |
CN103984803B (zh) | 一种有限元载荷施加方法 | |
CN105197253A (zh) | 一种机翼加筋壁板安全裕度计算方法 | |
CN106354934A (zh) | 一种基于有限元分析求解的加筋壁板损伤容限分析方法 | |
Trubina et al. | Geometric nonlinearity of the thin-walled profile under transverse bending | |
CN111046610B (zh) | 一种飞机整体翼梁无量纲应力强度因子的计算方法 | |
CN111914351B (zh) | 一种机身结构加筋壁板总体稳定性计算方法 | |
CN112784359B (zh) | 薄壁梁约束扭转极限承载力迭代计算方法 | |
CN107463746A (zh) | 一种气密载荷下机身隔框环向应力计算方法 | |
US9026406B2 (en) | Method, apparatus and computer program product for determining the strain induced at a selected point in a stiffened panel structure in response to a load, taking into account one or more out of plane (OOP) effects | |
CN111159943A (zh) | 一种动翼面封严结构的屈曲处理方法 | |
CN107526866B (zh) | 基于特征驱动的翼面结构拓扑优化方法 | |
CN107958126A (zh) | 一种双材料v型切口应力强度因子的计算方法 | |
CN104699873A (zh) | 一种飞机壁板开口结构数值分析方法 | |
CN109252441B (zh) | 一种变截面箱梁剪力滞效应的分析方法 | |
Dorbath et al. | Comparison of beam and shell theory for mass estimation in preliminary wing design | |
Dubovikov | Novel approach and algorithm for searching rational nonconventional airframe concepts of new generation aircrafts | |
CN105501462B (zh) | 一种机翼结构设计方法 | |
CN113051660B (zh) | 一种飞机机身门框区截面侧向弯曲刚度的计算方法 | |
Baker et al. | A practical approach to MDO and its application to an HSCT aircraft (multidisciplinary design optimization) |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |