CN103544402B - 一种疲劳开裂结构等效分析谱的构造方法 - Google Patents
一种疲劳开裂结构等效分析谱的构造方法 Download PDFInfo
- Publication number
- CN103544402B CN103544402B CN201310541197.5A CN201310541197A CN103544402B CN 103544402 B CN103544402 B CN 103544402B CN 201310541197 A CN201310541197 A CN 201310541197A CN 103544402 B CN103544402 B CN 103544402B
- Authority
- CN
- China
- Prior art keywords
- fatigue
- vacant lot
- stress
- load
- gag
- 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
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明属于疲劳试验技术领域,涉及一种疲劳开裂结构等效分析谱的构造方法。本发明基于外场数据,因此,得到的疲劳开裂结构等效分析谱能够反映外场实际,谱型为等幅谱,不需要再做载荷谱简化工作,即可应用于改进分析和对比试验,确定改进方案的有效性,大幅缩短试验周期。本发明解决了外场疲劳裂纹排故中结构疲劳载荷谱数据不全或过于复杂的问题。
Description
技术领域
本发明属于疲劳试验技术领域,涉及一种疲劳开裂结构等效分析谱的构造方法。
背景技术
飞机外场结构疲劳开裂排故,是飞机服役保障的一项很重要的工作内容,而在故障发生后,原因分析和结构改进都需要耗费大量的人力、物力和时间,其中从保障飞机飞行安全的角度,时间尤为关键,而发现原因和完成改进设计的前提,就是有相应结构的分析谱,现在面临的问题是,出现疲劳裂纹的结构往往疲劳载荷谱数据不完备,需要进行飞行实测,或者其疲劳载荷谱较为复杂,在做相关分析和改进验证性试验时,需要较长的时间周期。因此,使用较少的输入条件,构造简单的等效结构分析谱,用于飞机外场排故,节省经费和排故时间就显得特别的重要。
发明内容
本发明的目的是为了解决外场疲劳裂纹排故中结构疲劳载荷谱数据不全或过于复杂的问题,构造能用于故障分析的等效分析谱。
本发明的技术解决方案是,第一步:确定疲劳关键部位PSE的数据,疲劳关键部位包括疲劳开裂部位和相邻结构部位,数据包括飞机整机目标寿命N0和疲劳开裂部位寿命N1;
第二步:根据疲劳开裂结构部位,查阅相关疲劳分析手册、图册,获得同类飞机相似结构的地空地损伤比λj;
第三步:根据飞机疲劳开裂结构部位及周边部位的静强度气动载荷分布,得到相应部位段的静强度载荷P和作用点Z,用静强度载荷替代疲劳载荷,疲劳载荷P'=P/1.5;
第四步:计算疲劳关键部位PSE对应的细节疲劳额定值DFR,DFR=DFRbase·K·F·RC;其中,DFRbase是结构DFR基准值,K是材料系数,F是粗糙度系数,RC疲劳额定值系数;
第五步:按照飞机的目标寿命取N0、疲劳分析手册给定的疲劳可靠性系数FRF及计算得到的DFR,通过DFR法,得到不同应力比RGAG和地空地损伤比下的地空地循环最大应力σmaxGAG对应不同应力水平的疲劳寿命NR变化曲线;
第六步:建立疲劳关键部位段的有限元分析模型,施加疲劳载荷P',获得关键部位PSE的最大应力σ1;
第七步:选择一个应力比RGAG下的不同地空地损伤比的地空地循环最大应力σmaxGAG与NR的变化曲线,在这幅曲线中找出关键部位PSE的最大应力σ1和疲劳开裂部位寿命N1对应的点,若落在某个曲线上,则直接得到该应力比RGAG下的疲劳开裂部位的地空地损伤比λ',若落在两个曲线之间,则用插值的方式得到该应力比RGAG下的疲劳开裂部位的地空地损伤比λ';
第八步:依据疲劳载荷P'、地空地损伤比λ'得到构造的等效分析谱,谱型为等幅谱,该谱中地空地损伤比为λ',载荷循环中峰值等于疲劳载荷P',谷值为RGAG·P’。
所述的应力比RGAG设置为0。06,得到该应力比下的不同地空地损伤比下的地空地循环最大应力σmaxGAG对应不同应力水平的疲劳寿命NR变化曲线。
本发明具有的优点和有益效果,
本发明提出了一种疲劳开裂结构等效分析谱的构造方法。
本发明所需的数据包括飞机整机目标寿命N0,疲劳开裂部位寿命N1,静强度载荷P,同类飞机相似结构的地空地损伤比λj。构造过程包括获取疲劳载荷P′,计算疲劳开裂部位的细节疲劳额定值DFR,计算不同应力比RGAG和地空地损伤比下的地空地循环最大应力σmaxGAG对应不同应力水平的疲劳寿命NR变化曲线,计算疲劳开裂部位最大应力σ1,确定疲劳开裂部位的地空地损伤比λ',形成疲劳开裂部位等效分析谱。本发明不需要进行飞机载荷实测,亦不需要按照常规疲劳载荷谱编制过程进行复杂的当量折算和数据拟合,具有所需数据少,编制过程简单的优点,节省了人力、物力。
本发明基于外场数据,因此,得到的疲劳开裂结构等效分析谱能够反映外场实际,谱型为等幅谱,不需要再做载荷谱简化工作,即可应用于改进分析和对比试验,确定改进方案的有效性,大幅缩短试验周期。
本发明解决了外场疲劳裂纹排故中结构疲劳载荷谱数据不全或过于复杂的问题。
附图说明
图1是本发明在RGAG=0.06时,不同λ的σmaxGAG和NR变化曲线示意图。
具体实施方式:
下面结合附图对本发明作详细说明。
第一步:确定疲劳关键部位PSE的数据,疲劳关键部位包括疲劳开裂部位和相邻结构部位,数据包括飞机整机目标寿命N0和疲劳开裂部位寿命N1;
第二步:根据疲劳开裂结构部位,查阅相关疲劳分析手册、图册,获得同类飞机相似结构的地空地损伤比λj;
第三步:根据飞机疲劳开裂结构部位及周边部位的静强度气动载荷分布,得到相应部位段的静强度载荷P和作用点Z,用静强度载荷替代疲劳载荷,疲劳载荷P'=P/1.5;
第四步:计算疲劳关键部位PSE对应的细节疲劳额定值DFR,DFR=DFRbase·K·F·RC;其中,DFRbase是结构DFR基准值,K是材料系数,对于铝:K=1.0;对于钛:K=1.6;对于高强度钢(σb>1655MPa):K=1.9;对于中等强度钢(σb≤1655MPa):K=2.2,F是粗糙度系数,RC疲劳额定值系数;
第五步:按照飞机的目标寿命取N0、疲劳分析手册给定的疲劳可靠性系数FRF及计算得到的DFR,通过DFR法,其中单独作用地空地循环应力破坏时的飞行次数 S和σm0为疲劳特征参数,铝(全部):σm0=310MPa,S=2;钛(全部)σm0=620MPa,S=2;钢(中强):σm0=930MPa,S=1.8;钢(高强):σm0=1240MPa,S=1.8,得到不同应力比RGAG和地空地损伤比下的地空地循环最大应力σmaxGAG对应不同应力水平的疲劳寿命NR变化曲线;
第六步:建立疲劳关键部位段的有限元分析模型,施加疲劳载荷P',获得关键部位PSE的最大应力σ1;
第七步:选择一个应力比RGAG下的不同地空地损伤比的地空地循环最大应力σmaxGAG与NR的变化曲线,在这幅曲线中找出关键部位PSE的最大应力σ1和疲劳开裂部位寿命N1对应的点,若落在某个曲线上,则直接得到该应力比RGAG下的疲劳开裂部位的地空地损伤比λ',若落在两个曲线之间,则用插值的方式得到该应力比RGAG下的疲劳开裂部位的地空地损伤比λ';
第八步:依据疲劳载荷P'、地空地损伤比λ'得到构造的等效分析谱,谱型为等幅谱,该谱中地空地损伤比为λ',载荷循环中峰值等于疲劳载荷P',谷值为RGAG·P’。
该发明已应用到某型飞机腹鳍连接角材疲劳开裂排故,下面具体介绍该结构等效分析谱的构造过程。
实施例
第一步,确定疲劳关键部位包括疲劳开裂部位某型飞机腹鳍连接角材根部和相邻部位腹鳍连接角材螺栓连接孔,N0=5000起落数,N1=660起落数。
第二步:查飞机疲劳手册,λ1取波音737飞机垂尾与机身连接后梁端部接头的地-空-地损伤比λGAG=0.14,λ2取波音747飞机平尾结构的地-空-地损伤比λGAG=0.3,λ3亦取波音747飞机平尾结构的地-空-地损伤比λGAG=0.8。
第三步:根据静强度气动载荷分布,该结构主要承受侧向气动力引起的弯矩M。因此,分段计算侧面分布载荷,获得每段的平均弯矩,按照结构段,来分析对比找出,受载严重的关键部位段,确定加载点为713mm处时,疲劳载荷为8000N。
第四步:N0=5000飞行起落数,FRF取4,应力比RGAG取0.06,结构材料为铝合金,DFR=DFRbase·K·F·RC=133.8·1.0·0.95·1.44=183MPa。
第五步:结构为铝合金,疲劳特征参数σm0=310MPa,S=2,计算得到,RGAG=0.06时,三种地空地损伤比下,地空地循环最大应力σmaxGAG对应不同应力水平的疲劳寿命N变化曲线。
第六步:建立有限元模型,施加疲劳载荷P',应力分析得到疲劳开裂部位的最大应力σ1=270MPa。
第七步:查找关键部位的应力水平和疲劳寿命对应曲线中,对应σ1和N1的点,确认RGAG=0.06时,疲劳开裂部位的地空地损伤比λ'=λ1=0.14。
第八步:得到结构等效分析谱,谱型等幅谱,应力比RGAG=0.06,地空地损伤比λ=0.14,载荷循环为(8000N,480N)。
构造得到的分析谱,已成功用于结构疲劳分析和原结构与改进结构的对比试验,解决了某型飞机腹鳍连接角材疲劳开裂问题。
Claims (2)
1.一种疲劳开裂结构等效分析谱的构造方法,其特征是,第一步:确定疲劳关键部位PSE的数据,疲劳关键部位包括疲劳开裂部位和相邻结构部位,数据包括飞机整机目标寿命N0和疲劳开裂部位寿命N1;
第二步:根据疲劳开裂结构部位,查阅相关疲劳分析手册、图册,获得同类飞机相似结构的地空地损伤比λj;
第三步:根据飞机疲劳开裂结构部位及周边部位的静强度气动载荷分布,得到相应部位段的静强度载荷P和作用点Z,用静强度载荷替代疲劳载荷,疲劳载荷P′=P/1.5;
第四步:计算疲劳关键部位PSE对应的细节疲劳额定值DFR,DFR=DFRbase·K·F·RC;其中,DFRbase是结构DFR基准值,K是应力集中系数,F是粗糙度系数,RC疲劳额定值系数;
第五步:按照飞机的目标寿命取N0、疲劳分析手册给定的疲劳可靠性系数FRF及计算得到的DFR,通过DFR法,得到不同应力比RGAG和地空地损伤比下的地空地循环最大应力σmaxGAG对应不同应力水平的疲劳寿命NR变化曲线;
第六步:建立疲劳关键部位段的有限元分析模型,施加疲劳载荷P',获得关键部位PSE的最大应力σ1;
第七步:选择一个应力比RGAG下的不同地空地损伤比的地空地循环最大应力σmaxGAG与NR的变化曲线,在这幅曲线中找出关键部位PSE的最大应力σ1和疲劳开裂部位寿命N1对应的点,若落在某个曲线上,则直接得到该应力比RGAG下的疲劳开裂部位的地空地损伤比λ',若落在两个曲线之间,则用插值的方式得到该应力比RGAG下的疲劳开裂部位的地空地损伤比λ';
第八步:依据疲劳载荷P'、地空地损伤比λ'得到构造的等效分析谱,谱型为等幅谱,该谱中地空地损伤比为λ',载荷循环中峰值等于疲劳载荷P',谷值为RGAG·P。
2.根据权利要求1所述的一种疲劳开裂结构等效分析谱的构造方法,其特征是,所述的应力比RGAG设置为0.06,得到该应力比下的不同地空地损伤比下的地空地循环最大应力σmaxGAG对应不同应力水平的疲劳寿命NR变化曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310541197.5A CN103544402B (zh) | 2013-11-05 | 2013-11-05 | 一种疲劳开裂结构等效分析谱的构造方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310541197.5A CN103544402B (zh) | 2013-11-05 | 2013-11-05 | 一种疲劳开裂结构等效分析谱的构造方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103544402A CN103544402A (zh) | 2014-01-29 |
CN103544402B true CN103544402B (zh) | 2016-08-24 |
Family
ID=49967847
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310541197.5A Active CN103544402B (zh) | 2013-11-05 | 2013-11-05 | 一种疲劳开裂结构等效分析谱的构造方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103544402B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103902820B (zh) * | 2014-03-22 | 2016-09-14 | 淮阴工学院 | 一种组合路况下车辆结构疲劳损伤计算的当量系数法 |
CN104034576B (zh) * | 2014-06-16 | 2016-11-02 | 中国飞机强度研究所 | 一种飞机结构裂纹扩展试验用载荷谱简化方法 |
CN104679933A (zh) * | 2014-07-30 | 2015-06-03 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机疲劳载荷设计方法 |
CN105574307B (zh) * | 2014-10-11 | 2018-05-29 | 中国航空工业集团公司西安飞机设计研究所 | 一种展向连接结构dfr双向修正临界点的判断方法 |
CN112268799B (zh) * | 2020-10-16 | 2022-09-09 | 中国直升机设计研究所 | 一种复合材料结构静强度和疲劳强度一体化试验验证方法 |
CN114428026A (zh) * | 2021-12-30 | 2022-05-03 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机机身蒙皮环向对缝带板连接结构疲劳分析方法 |
-
2013
- 2013-11-05 CN CN201310541197.5A patent/CN103544402B/zh active Active
Non-Patent Citations (6)
Title |
---|
军用飞机结构耐久性设计的细节疲劳额定值方法;董彦民;《航空学报》;20101231;第31卷(第12期);全文 * |
小型无人机机体寿命分析研究;李光超;《中国优秀硕士学位论文全文数据库》;20070615;第2007年卷(第6期);全文 * |
某型飞机前墙6-8肋段连接区疲劳分析;汤庆辉;《中国优秀硕士学位论文全文数据库》;20030315;第2003年卷(第1期);全文 * |
某航空发动机使用载荷对其整机寿命和可靠性的影响研究;公杰;《中国优秀硕士学位论文全文数据库》;20130115;第2013年卷(第1期);全文 * |
飞机疲劳载荷谱分析的一种新方法;李天亮 等;《机械科学与技术》;20060630;第25卷(第6期);全文 * |
飞机结构WFD发生的概率分析;雷晓欣 等;《机械科学与技术》;20111130;第30卷(第11期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103544402A (zh) | 2014-01-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103544402B (zh) | 一种疲劳开裂结构等效分析谱的构造方法 | |
CN109033526B (zh) | 一种翼肋与蒙皮铆钉连接载荷计算方法 | |
RU2599108C1 (ru) | Способ мониторинга нагрузок и накопленной усталостной повреждаемости в условиях эксплуатации самолета | |
CN109902330B (zh) | 一种无长桁和肋腹板的翼面及其强度计算方法 | |
CN113837639A (zh) | 一种基于多维指标的航空器遭遇尾流风险评估系统 | |
CN108984909B (zh) | 一种基于Mar-Lin模型的含大损伤飞机复合材料结构剩余强度分析方法 | |
CN102944211A (zh) | 一种飞机整体翼梁止裂筋条面积的确定方法 | |
CN105334032B (zh) | 一种机型结构寿命的计算方法 | |
Sitaraman et al. | Coupled CFD/CSD prediction of rotor aerodynamic and structural dynamic loads for three critical flight conditions | |
Molent et al. | Damage tolerance demonstration testing for the Australian F/A-18 | |
Sarath et al. | Analysis of the wingbox with spliced skin and estimation of the fatigue life for the wingbox | |
Vargas-Jimenez et al. | Development of a Wave Drag Prediction Method for the Conceptual Design Phase | |
CN105823688A (zh) | 一种翼面前缘曲面结构的承载试验方法 | |
Atmeh et al. | Design and stress analysis of a general aviation aircraft wing | |
Woodson et al. | Understanding abrupt wing stall with CFD | |
Cunningham Jr et al. | A system analysis study comparing reverse engineered combinatorial testing to expert judgment | |
Breitsamter | Aerodynamic efficiency of high maneuverable aircraft applying adaptive wing trailing edge section | |
Snider et al. | Residual strength and crack propagation tests on C-130 airplane center wings with service-imposed fatigue damage | |
CN106769533A (zh) | 一种复合材料隔框结构四点弯曲试验方法 | |
Chipman et al. | Numerical computation of aeroelastically corrected transonic loads | |
Castells et al. | A numerical investigation of the influence of the blade-vortex interaction on the dynamic stall onset in simplified rotor environment | |
Lefebvre | Cylinder-Airfoil Interactions and the Effect on Airfoil Performance | |
Pickhaver et al. | A technique to predict the aerodynamic effects of battle damage on an aircraft’s wing | |
Himisch et al. | The influence of target function selection on the optimization of winglets for the glider SB 15 | |
Giglio et al. | Comparison of fatigue crack propagation behavior of Al 2024 and Al–Li 8090 helicopter fuselage panels |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |