用SPH和有限元方法研究鸟撞飞机风挡问题

合集下载

鸟撞飞机圆弧风挡实验的研究及数值模拟

鸟撞飞机圆弧风挡实验的研究及数值模拟

查些:坚三奎童丝圭垒塞耋堡丝塞一挡是甭失效。

1999年在太原理:『1人学应FH力学研究所的帮助h川I高速搬影机列歼八一II风捎的鸟撞过程进行了拍摄(如图I一9),对’;0撒n勺全]三[程仃了进步的认谚i。

图1.9高速摄影照片(5000幅/秒)随着有限元方法的完善及计算机技术的发展,以美国为首的围家又相继开发了其它可用于鸟撞分析的有限元计算软件,如LS—D'Y'NA¨…、DYTRAN等。

这些软件都是国际著名的结构瞬态动力非线性;1阳应分析软件,而且能进行流固结构的辎合汁算。

通过汁算机模拟非线性乌撞过槲,从lm大大节省研制新风挡和座舵盖的时削和经费。

Boroughs、RobertR.¨…(1998年)利用DYNA3D计算了乌撞[,earjet45飞机风挡玻璃的过程,他们用一个较详细的有限元模型描述了胍挡玻璃和其附属结构,并且与以自口的专门用于分析鸟撞风挡的控序进行了比较,阿到用这个模型来求解鸟撞Leauet45飞机风挡玻璃更合理的结论。

王爱俊等【^l胁1(1998年、1999年)利用LS.DYNA3D程J手作为i三安分析丁.具,采用碰撞接触有限元算泫,刘层合胍挡进行了鸟撺模拟。

采用AI,E天踩理5-k母碗士辱^j一论叉§2.1实验方法本实验采用幽际通用测试飞机风挡玻璃抗岛撞的实验方法(如图2-1所示)。

将规定质量的鸟弹装入鸟弹利壳,通过空气炮发射“呜弹”(吗墩活鸡代替,质量为I8Kg),撞击安装于台架上的-|毛机全尺、J‘风挡,水模拟空中的鸟撞。

划国产、进【J两种型号的圳弧jxb}"‘i驶璃进行,全』0、¨々撼试验。

试件参照飞机上的安簧角度安装于试验台架上,呈剁撞击姿势。

圈2—1鸟撞圆9i风挡的实验圈在实验中,位移传感器、加速发传感器安装于风挡内衣西I下方测量位移、加速度。

采用超动态应变仪测量风挡玻璃典型位戳的应变。

,本实验采用断丝法测量鸟速,用高速摄影机和高速摄像机从不同角度拍摄鸟撞的全部过程。

飞机尾翼前缘结构鸟撞模型与试验验证

飞机尾翼前缘结构鸟撞模型与试验验证

飞机尾翼前缘结构鸟撞模型与试验验证发表时间:2019-12-30T13:08:15.447Z 来源:《科学与技术》2019年 15期作者:倪磊[导读] 鸟撞是飞机在飞行中遇到的重要危险之一摘要:鸟撞是飞机在飞行中遇到的重要危险之一,同时也是一种突发性和多发性的飞行事故,因此,结构抗鸟撞设计成为飞机设计必须考虑的要素。

本文针对某飞机尾翼前缘结构基于光滑粒子流模拟了结构遭受鸟撞时结构损伤的全过程;以前缘结构鸟撞试验为基础,研究了前缘结构受鸟撞击的破坏模式;分析表明:模型的数值结果与试验结果在前缘结构的变形、破坏模式和应变曲线趋势吻合较好。

关键词:鸟撞;SPH方法;前缘结构中图分类号:V215.9, V216.5 文献标志码:A0 引言飞机飞行期间,飞机的迎风面对于鸟撞是敏感部位,这些结构中常设有各种管路系统等重要设施,一旦遭到破坏,飞机的安全性能就得不到保障,因此结构的抗鸟撞能力是设计必须考虑的。

过去通过昂贵的全尺寸鸟撞试验来验证,如果鸟撞验证试验未通过,这就会大大增加飞机的研制风险。

因此,为了降低设计阶段成本和风险,数值模拟方法被用来支持取证过程[1]。

一个精确的数值模拟模型能够让设计者对鸟撞撞击事件中所表现出的力学行为有更加明确的了解,本文采用大型非线性有限元方法,基于光滑粒子流建立了尾翼前缘结构鸟撞的数值分析模型,详细模拟了结构遭受鸟撞时结构损伤的全过程;同时,开展了飞机尾翼前缘结构鸟撞试验,研究前缘结构受鸟撞击的破坏模式。

最后,对比分析了模型和试验的前缘结构变形、破坏模式和应变曲线趋势,吻合比较好。

1 尾翼前缘结构根据某飞机平尾前缘危险部位筛选情况来看,前缘蒙皮、隔板和前梁腹板的对接位置为最危险部位,因该处的对接均在展向同一位置。

因此,从对接位置向左向右各延伸4个肋站位截取飞机平尾前缘结构中一段来作为分析和试验的目标,试件总长度1523mm,8个肋站位保证了试件遭受鸟撞时,能量能够充分的扩散。

飞机翼面结构抗鸟撞设计研究

飞机翼面结构抗鸟撞设计研究

飞机翼面结构抗鸟撞设计研究首先,文章针对抗鸟撞的研究活动,提出了抗鸟撞设计的计算模型,并分别对鸟体及结构的几何模型及材料参数设置进行了说明,提供了相应的计算公式;接下来,又针对有限元模型的建立,提供了适合研究鸟撞的SPH方法,同时阐述了网格划分及边界条件处理;最后,则在模拟结果的基础上展开了讨论,分析了撞击过程中各个结构所起到的作用。

标签:翼面结构;抗鸟撞设计;SPH方法前言随着经济的快速发展,民用飞机正在被大量使用,但随之也引发了大量的问题。

飞机在迎风飞行状态中,极易遭受鸟体撞击,尤其是飞机的翼面结构、机翼前缘等部位。

如果这些结构遭到破坏,就无法保障飞机的安全飞行。

因此,在飞机一翼面结构的设计上,不仅应符合空气动力学的相关标准,同时还应充分考虑鸟撞的发生,从而提升结构的强度。

但如果一味对强度进行提升,就会使结构重量增加,不利于飞机的性能维护,同时增加成本。

因此,为了使这一矛盾得以圆满地解决,就必须更多地运用新材料、新设计,最大程度地保证飞机的安全性能。

1 抗鸟撞设计的计算模型1.1 鳥体及结构的几何模型通过建立鸟体的几何模型,可知在整个缝翼结构中,分别由前后蒙皮、肋、梁及其他角材连接构成。

在计算过程中,前蒙皮将被视为均匀体,设定值为1.6mm;肋厚度为1.02mm;梁的厚度为1.8mm。

在鸟体形状上,采用实心的圆柱体进行模拟,两端均设为半球状,长径比为2:1,重量为 1.8kg,密度为900kg/m3,由此可以确定,圆柱直径,即半球体的直径为115mm[1]。

1.2材料参数设置在本次研究中,假设鸟体冲击速度恒定,设为150m/s,同时利用SPH来建构鸟体模型,可得到如下的本构关系:P=P0+B{{■}γ-1} (1)对正撞击,B=1.12×108Pa,γ=6.77;发生撞击时,如果角度恰在90°-45°之间,则B=1.28×108Pa,γ=7.98。

本次研究中将遵循这一情况,即取B=1.28×108Pa,γ=7.98。

基于abaqus的SPH算法鸟撞分析研究

基于abaqus的SPH算法鸟撞分析研究

基于abaqus 的SPH 算法鸟撞分析研究龙思海,滕春明,张海东(航空工业洪都,江西南昌,330024)要:本文采用abaqus 的SPH 算法进行了风挡鸟撞性能分析研究。

分析结果表明,基于的分析结果与试验结果趋势一致,风挡鸟撞临界速度与试验结果一致。

:风挡;抗鸟撞;abaqusAnalysis and Study of Bird Strike Based on SPH Algorithm of AbaquesLong Sihai,Teng Chunming,Zhang Haidong (AVIC Hongdu ,Nanchang ,Jiangxi ,330024)Abstract:This paper presents SPH algorithm of abaques for carrying out the analysis and study of wind strike.The analysis result shows that the analysis result obtained based on SPH identical with the trend of the test result,the critical velocity of wind shield bird strike is test result.Wind shield;Anti bird strike;Abaque0引言随着飞机低空高速飞行任务的增加,鸟撞飞机的概率也随之增加。

鸟撞飞机带来的危害已越来越为人们所认识。

由于工程试验法试验周期长、次数多,导致费用很高,例如在某型飞机风挡的研制过程中,全尺寸鸟撞试验就进行了多次,制造了大量的试验件,试验费用极高。

因此,鸟撞动响应仿真分析成为重点的研究方向,型号设计前期准确的鸟撞动响应仿真分析能够节省大量的试验费用。

张志林[1]等着重考虑了应变率对透明件材料性能的影响和几何非线性对刚度矩阵的影响,结果显示:应变率对位移、应变影响较大,考虑应变率相关性分析所得结果比不考虑应变率相关性分析结果更接近试验结果;几何非线性分析所获得的风挡最大法向位移比线性分析得到的值大,几何非线性对飞机风挡鸟撞动响应分析结果的影响不可忽略,并以此提出了鸟撞击载荷柔性靶理论。

基于SPH方法的LY12-CZ铝合金平板鸟撞模型

基于SPH方法的LY12-CZ铝合金平板鸟撞模型

基于SPH方法的LY12-CZ铝合金平板鸟撞模型张鼎逆;上官倩芡;刘富【摘要】To obtain the accurate analysis model of bird impact for aircraft structures,the numerical model of bird impact on LY1 2-CZ aluminum alloy flat plates was established based on the explicit finite element code of PAM-CRASH.Johnson-Cook equation was adopted to describe the constitutive model of LY1 2-CZ aluminum alloy.The stress-strain curves at 4 different strain rates of LY1 2-CZ were measured by dynamic tensile tests using electronic universal testing machine and split Hopkinson tensile bar (SHTB),and four constants in Johnson-Cook equation were fitted.Bird was modeled by smooth particle hydrodynamics (SPH ) method,and Monaghan EOS equation was introduced to simulate the bird material.Two corresponding bird impact experiments were implemented,and the strains of testing points were obtained.The calculation results were compared with the bird impact experimental results.The results show that the good agreement between numerical calculation and experimental results is obtained, which illuminates the bird constitutive model,the aluminum alloy constitutive model and the calculation method for high speed are reasonable and reliable.%为了得到更准确的飞机鸟撞分析模型,基于显式有限元分析程序PAM-CRASH建立了LY12-CZ铝合金平板的鸟撞数值模型.采用Johnson-Cook方程表述LY12-CZ材料的本构.LY12-CZ铝合金在4种不同应变率下的应力应变曲线通过电子万能试验机和分离式霍普金森拉杆(SHTB)拉伸试验获得,对曲线进行拟合得到Johnson-Cook 方程中的4个常数.基于光滑粒子流体动力学(SPH)方法建立了鸟体模型,引入Monaghan EOS 方程来描述鸟体材料.针对所建立的鸟撞数值计算模型,开展相对应的鸟撞试验,获得测试点的应变,并将数值计算结果和试验结果进行对比.结果表明:计算应变与试验测得的应变吻合较好,验证了鸟体、铝合金本构模型以及鸟体高速冲击计算分析模型的合理性、可靠性.【期刊名称】《江苏大学学报(自然科学版)》【年(卷),期】2016(037)004【总页数】5页(P418-422)【关键词】鸟撞;本构模型;SHTB;光滑粒子流体动力学;试验;数值计算【作者】张鼎逆;上官倩芡;刘富【作者单位】上海师范大学信息与机电工程学院,上海200234;上海师范大学信息与机电工程学院,上海200234;中国商用飞机有限责任公司上海飞机设计研究院,上海201210【正文语种】中文【中图分类】O347.3在飞机起飞着陆阶段,鸟撞问题严重威胁着乘客的生命安全[1].因此,民用航空规章已经针对飞机的鸟撞事件制定了相应的标准.近年来,诸多学者对飞机结构的抗鸟撞问题进行了大量研究.研究表明在鸟体高速撞击的过程中,存在着明显的大变形特性和流固耦合现象.鸟体被瞬间撞击,鸟体粉碎、飞溅过程类似于液体的流动.所以,建立结构尤其是鸟体本构模型是鸟撞非线性分析的难点和重点.对于合金和玻璃等材料,受到高速冲击时,材料的应变率效应往往是不可忽略的.结构的鸟撞数值分析方法主要有Lagrangian有限元法(finite element, FE)、任意拉格朗日-欧拉法(arbitrary Lagrange-Euler, ALE)和光滑粒子流体动力学方法(smooth particle hydrodynamics, SPH).Lagrangian有限元法不能解决大变形引起的网格畸变问题.任意拉格朗日-欧拉法很难模拟靶板穿透问题.光滑粒子流体动力学方法采用的是无网格拉格朗日技术,由于粒子在空间相互独立,因此比起以上2种方法更适用于鸟撞等高速碰撞问题.J. A. REGLERO等[2]提出了一种泡沫铝填充的抗鸟撞前缘,试验结果表明泡沫铝材料填充的结构,其支反力相对于其他结构较小,吸能效果较好.A.G.HANSSEN等[3]研究了AlSi7Mg0.5-泡沫夹芯板和2024-T3铝合金平板的鸟撞动态响应,2024-T3的本构模型通过力学性能试验获得.I.SMOJVER等[4]基于ABAQUS软件,分析了不同尺寸鸟体、撞击位置和速度对飞机复杂襟翼结构的冲击损伤影响.M. GUIDA等[5]采用Lsdyna和MSC-Dytran模拟了由铝-蜂窝结构和GLARE材料组成的复合材料前缘盒段鸟撞过程,计算了整个撞击过程结构的动态响应和破坏.刘军等[6-7]则对鸟体本构参数进行了一系列的研究,他们利用鸟撞平板试验,并使用神经网络方法对试验中的鸟体参数进行了反演,将反演得出的参数代入模型加以计算;模拟结果与试验结果比较吻合,所得参数为以后的鸟撞计算提供了参考.文中拟基于显式有限元分析程序PAM-CRASH,通过耦合SPH法和FE法建立LY12-CZ铝合金平板的鸟撞数值模型,并开展相应的鸟撞试验.进而对比计算和试验数据,以验证鸟撞数值分析模型的合理性和可靠性.鸟撞是发生在毫秒量级的瞬态动力学问题,此时一般都要考虑材料的应变率影响.Johnson-Cook方程是能够表征材料硬化效应,应变率强化效应和温度软化效应的弹塑性模型[8].LY12-CZ铝合金的本构模型由如下的Johnson-Cook方程表示:式中: σ为等效应力;ε为等效塑性应变0为量纲一化的塑性应变率0为参考应变率;T*=(T-Tr)/(Tm-Tr),其中Tr为室温,Tm为材料熔点.材料硬化效应,应变率强化效应和温度软化效应分别由和[1-(T*)m]描述.图1为霍普金森拉伸拉杆加载的示意图.工程应力、工程应变和应变率表示为式中: E为加载杆的弹性模量;A为加载杆的横截面积;As为试件的横截面积;L为试件的厚度;C0为波速;εI(t),εR(t)和εT(t)分别为入射应变、反射应变和透射应变,且有εI(t)=εR(t)+εT(t).开展4个应变率下的材料力学性能试验,准静态(应变率为0.001 s-1)试验采用电子万能试验机,高应变率(1 500,1 900,2 300 s-1)试验采用霍普金森拉杆.通过动态拉伸力学试验,获得材料的4条应力-应变曲线.不考虑材料的温度软化效应.图2为动态拉伸试验前后的试验件.图3为4种应变率下的σ-ε曲线.由图3可见,应力-应变曲线表明从准静态(0.001 s-1)到高应变率(2 300 s-1),屈服应力显著增加,但同一数量级应变率下的应力-应变关系基本接近.文中参考应变率取0.001 s-1,则通过应变率0.001 s-1的应力-应变曲线,计算出A,B和n.最后,通过4条应力-应变曲线,拟合出Johnson-Cook方程里的4个常数: A=380 MPa,B=520 MPa,C=0.03,n=0.465.2.1 鸟体模型鸟体形状为两端半球、中间圆柱,长径比为2,质量为1.8 kg,密度为950 kg·m-3.传统的基于网格的计算方法很难处理鸟体在高速撞击下的液体飞溅问题,而基于无网格的光滑粒子流体动力学方法则可以有效地解决此问题.因此,采用SPH单元来模拟鸟体.鸟体的本构通过Monaghan EOS方程[9-10]来描述,定义为式中: p0为参考压力;ρ0为初始密度;常数B和γ借助鸟撞试验后的参数反演来确定,文中取B=128 MPa,γ=7.98.2.2 计算模型图4为LY12-CZ铝合金平板固定在试验台上,进行了2种不同撞击速度的鸟撞试验,速度分别为70,120 m·s-1;平板长和宽均为600 mm,厚度为10 mm,撞击点位于平板的中心.应变测试点如图5所示.由图5可见,考虑到应变片在高速撞击中容易破坏,为了防止应变片损坏造成的数据丢失,在对称位置布置了另外4个应变片(S5,S6,S7,S8).夹具采用实体单元模拟,平板采用壳单元建模.定义鸟体(SPH)和平板(FE)之间的耦合接触,平板和夹具之间定义一般接触.夹具上的6个位置固支约束模拟夹具与台架的连接,板和夹具之间采用刚性单元连接来模拟螺栓.整个鸟撞计算模型如图6所示.2.3 结果分析图7为鸟撞速度v=70 m·s-1时铝合金平板等效应力云图.从撞击时刻起,鸟体被逐渐压缩,随后鸟体解体并像液体一样飞溅.应力集中区域位于平板撞击中心及螺栓连接处.图8和9给出了2种鸟撞速度下,测试点处计算值与试验值的应变-时间曲线.由图8-9可见,无论是应变峰值还是应变的变化趋势均与试验值吻合较好,表明文中的模型是比较准确的.随着撞击速度的提高,测试点的应变显著增加;对于特定撞击速度下,越靠近中心,应变值越大.S1测试点靠近靶板边缘,在撞击的初始时刻,边界条件对测试点的影响比鸟体更大,计算模型中的边界刚度要大于真实试验.因此,应变-时间曲线上第1个峰值的计算误差相对于其他测试点更大,且计算值高于试验值.随着时间的推移,鸟撞引起的应力波传播至平板边缘,并起主导作用,计算误差减小.图9中,S2,S3,S4测试点处的试验数据不完整,因为撞击速度120 m·s-1过高,从而S2,S3,S4,S5,S6,S7测试点处的应变片发生损坏,故而未采集到数据.1) 计算所得的测试点的应变响应与试验值吻合较好,表明拟合的LY12-CZ 铝合金Johnson-Cook本构方程、鸟体EOS方程以及瞬态冲击计算方法是合理、可靠的.2) 当鸟撞速度相对较低时,在撞击初始阶段,靠近平板边缘的测试点的动态响应受平板边界条件影响较大,随着应力波从撞击中心向平板四周扩散,边界条件的影响削弱.通过体单元模拟平板和夹具的螺栓连接取代文中的刚性连接能够有效减少边界刚度的影响.在后期的飞机结构鸟撞分析中将对边界条件进行修正和细化.3) 文中基于积木式试验验证的思路,从试片级的材料力学性能试验到单元级的平板鸟撞分析及试验验证为飞机结构的抗鸟撞研究提供了有力的技术支持.【相关文献】[1] WASHBURN B E, CISAR P J, DEVAULT T L. Wildlife strikes to civil helicopters in the US, 1990—2011[J]. Transportation Research Part D: Transport and Environment, 2013, 24: 83-88.[2] REGLERO J A, RODRGUEZ-PÉREZ M A, SOLRZANO E, et al. Aluminium foams as fillerfor leading edges: improvements in the mechanical behavior under bird strike impacttests[J]. Materials and Design, 2011, 32(2):907-910.[3] HANSSEN A G, GIRARD Y, OLOVSSON L, et al. A numerical model for bird strike of aluminium foam-based sandwich panels[J]. International Journal of Impact Engineering, 2006, 32(7): 1127-1144.[4] SMOJVER I, D. Numerical simulation of bird strike damage prediction in airplane flap structure[J]. Composite Structures, 2010, 92(9): 2016-2026.[5] GUIDA M, MARULO F, MEO M, et al. SPH-Lagran-gian study of bird impact on leading edge wing[J]. Composite Structures, 2011, 93(3): 1060-1071.[6] 刘军,李玉龙,郭伟国,等.鸟体本构模型参数反演Ⅰ:鸟撞平板试验研究[J].航空学报,2011,32(5):802-811.LIU J, LI Y L, GUO W G, et al. Parameters inversion on bird constitutive model part I: study on experiment of bird striking on plate[J]. Acta Aeronautica et Astronautica Sinica, 2011,32(5):802-811. (in Chinese)[7] 刘军,李玉龙,石霄鹏,等. 鸟体本构模型参数反演II:模型参数反演研究[J].航空学报,2011,32(5):812-821.LIU J, LI Y L, SHI X P, et al. Parameters inversion on bird constitutive model part II:study on model parameters inversion[J]. Acta Aeronautica et Astronautica Si-nica, 2011, 32(5):812-821. (in Chinese)[8] ZHANG D N, SHANGGUAN Q Q, XIE C J, et al. A modified Johnson-Cook model of dynamic tensile beha-viors for 7075-T6 aluminum alloy[J]. Journal of Alloys and Compounds, 2015, 619: 186-194.[9] MONAGHAN J J, RAFIEE A. A simple SPH algorithm for multi-fluid flow with high density ratios[J]. International Journal for Numerical Methods in Fluids, 2013, 71(5): 537-561.[10] VALIZADEH A, MONAGHAN J J. A study of solid wall models for weakly compressible SPH[J]. Journal of Computational Physics, 2015, 300: 5-19.。

无人机飞鸟撞击机翼损伤程度预测仿真

无人机飞鸟撞击机翼损伤程度预测仿真
第35卷 第9期
文章编号 :1006—9348(2018)09—0042—04
计 算 机 仿 真
2018年9月
无 人 机 飞 鸟撞 击 机 翼 损伤 程度 预 测 仿 真
王永 虎 ,昊 志坚 ,杨 敏
(中国民航飞行学院 飞行技术学院 ,四川 广汉 618307)
摘要 :日益普及的无人机“黑 飞”现象频频发生,严重危及航空安全和正常航 班运行 ,无人机撞击 民机关键 部位损伤程度的研
(Flight Technology Institute,Civil Aviation Hight University of China,Guanghan Sichuan 618307,China)
A BSTRA CT :The ”unregistered flight” phenomenon OCCUr s frequently,which ser iously endanger s aviation safety and normal flight operation.The research on damage degree of important part of civil aircra ft a fter the impact of u机翼 ;损伤 ;仿真
中 图 分 类号 :TP319.9 文献 标 识 码 :B
The Dam age Prediction and Sim ulation for the UAV and Birdstrike Im pact On W ing
WANG Yong-hu,WU Zhi-jian,YANG Min
energ y cha n ged with time a nd the deformation law of leading edge were obtained .Through the comparison,the dam-

基于SPH方法的鸟撞复合材料层合板数值分析

基于SPH方法的鸟撞复合材料层合板数值分析

基于SPH方法的鸟撞复合材料层合板数值分析
于永强;李成;铁瑛
【期刊名称】《玻璃钢/复合材料》
【年(卷),期】2017(000)005
【摘要】采用光滑粒子流体动力学法(SPH)耦合有限元法对复合材料层合板受鸟撞击的过程进行了数值模拟.复合材料层合板采用渐进损伤模型,鸟体采用SPH粒子建立模型,利用ANSYS/LS-DYNA显示动力分析模块分析了复合材料层合板结构非线性接触.分析了鸟撞层合板过程中鸟体损伤及层合板单层纤维失效和基体失效情况,分析了鸟体的入射角方向及层合板采用不同铺层时对层合板吸能效果的影响.计算结果表明,合理设计层合板铺层可以提高层合板的吸能效果.
【总页数】5页(P48-52)
【作者】于永强;李成;铁瑛
【作者单位】郑州大学,郑州 450000;郑州大学,郑州 450000;郑州大学,郑州450000
【正文语种】中文
【中图分类】TB332
【相关文献】
1.复合材料加筋层合板准静态压痕实验研究及数值分析 [J], 黄飞生;周仕刚;宫占锋;罗腾腾;薛元德
2.采用流固耦合方法的复合材料层合板鸟撞分析 [J], 李成;郭文辉;铁瑛
3.基于SPH方法鸟撞航空发动机进气风扇的数值分析 [J], 杨瑞进;姜楠;池剑虹;孟宪国;邹宜勇
4.基于非线性接触的修理后层合板鸟撞损伤特性研究 [J], 徐建新;马行驰;赵鲁春
5.复合材料层合板鸟撞损伤及吸能影响因素数值分析 [J], 于连超;陈伟;关玉璞;温海涛
因版权原因,仅展示原文概要,查看原文内容请购买。

基于SPH算法的襟翼鸟撞分析研究

基于SPH算法的襟翼鸟撞分析研究

过内部转换器由六面体单元转换生成。

图1内襟翼结构图2内襟翼结构(隐藏蒙皮)图3内襟翼边界条件Science&Technology Vision科技视界图4鸟体模型考虑速度、迎角、卡位、撞击位置等伤的影响,这里选取四个典型的撞中间、肋与后梁连接处、盒段腹板如图5所示。

鸟体撞击速度选取am机迎角选取三种情形d°、e°、f°,襟g°。

这里假设以上所设置的工况情到所选位置。

:a<b<c;d<e<f。

图5鸟撞击位置SPH算法SPH方法(Smoothed Particle Hydrodynamics)的核心是一种插值技术。

每一个粒子与其相距设定距离范围内的所有其他粒子发生相互作用。

它们间的相互作用是由未知函数来衡量的,设定距离为光滑长度的两倍。

鸟体本构模型采用Murnaghan状态方程[5]。

此模型中状态方程为()[]1撞击过程图6所选工况下鸟撞结构破坏过程体的撞击过程为例,说明如下:图7模型能量变化曲线从图7中可以看出,计算的最后时刻系统的动能并未完全耗散,动能的耗散主要发生在撞击之后的5ms内,占总耗散动能的35%左右,当然这与初始速度密切相关。

该工况下,结构损伤情况见图8-图9。

结构损伤最大值发生在后梁腹板处。

图8撞击后整体结构损伤云图图9撞击后前后梁局部损伤云图3.3不同工况下鸟体的撞击结果对比分析这里着重关注下内襟翼前后梁、肋和蒙皮上的结构最大损伤,图10、图11、图12为飞机迎角分别为d°、e°、f°情形下对应不同速度和不同撞击位置下的结构损伤最大值对比。

图13为撞击位置3对应的不同飞. All Rights Reserved.Science&Technology Vision科技视界结构损伤更为明显。

图10迎角d°图11迎角e°图12迎角f°3.4沙漏现象使用单点积分的有限元分析中,显示积分算法因大变形、单元畸变等原因可能会出现沙漏现象,观察计算结果中各部件的沙漏现象。

基于响应面法的鸟撞风扇叶片损伤预测

基于响应面法的鸟撞风扇叶片损伤预测

第 50 卷第 2 期2024 年 4 月Vol. 50 No. 2Apr. 2024航空发动机Aeroengine基于响应面法的鸟撞风扇叶片损伤预测刘洋1,罗忠2,张海洋1,沈丽娟1,赵凤飞1,张宗锴1(1.中国航发沈阳发动机研究所,沈阳 110015; 2.东北大学机械工程与自动化学院,沈阳 110819)摘要:有限元模拟鸟撞风扇叶片损伤成本高,为解决工程问题,采用经典叶栅鸟撞切割模型建立了鸟撞风扇叶片动载荷数学模型,结合鸟撞部件试验结果,以拟合技术明确风扇叶片损伤程度与最大关键动载荷计算值间的函数关系,形成叶片损伤预测响应面,实现对鸟撞风扇叶片损伤的快速预测,并建立基于响应面法的鸟撞风扇叶片损伤预测工作流程。

结合涡扇发动机吞鸟试验技术要求、风扇结构设计特征及已开展的鸟撞部件试验结果,建立叶片损伤预测响应面,初步识别2种鸟撞方案的径向弯曲、弦向弯曲,并计算撕裂范围分别不超过0.3867和0.3941,撕裂与弦向弯曲相关性显著,呈抛物线变化趋势。

结果表明:预测的损伤在可接受的安全性水平范围内,预测方法能够识别损伤范围及趋势,可为后续鸟撞有限元模拟、试验策划、安全性分析、风扇叶片抗鸟撞设计等工作提供量化的技术支持。

关键词:鸟撞;风扇叶片;损伤;响应面法;航空发动机中图分类号:V232.4文献标识码:A doi:10.13477/ki.aeroengine.2024.02.012 Prediction of Damage to Fan Blade from Bird Strike Based on Response Surface MethodLIU Yang1, LUO Zhong2, ZHANG Hai-yang1, SHEN Li-juan1, ZHAO Feng-fei1, ZHANG Zong-kai1(1. AECC Shenyang Engine Research Institute,Shenyang 110015;2. School of Mechanical Engineering & Automation,Northeastern University,Shenyang 110819)Abstract:The cost of finite element simulation of bird strike fan blade damage is high. In order to solve engineering problems, the classic bird strike cascade slicing model is used to establish the mathematical model for the dynamic load of bird impact. Combined with the results of bird strike component tests, the fitting technology was used to clarify the functional relationship between the degree of fan blade damage and the calculated value of the maximum critical dynamic load, and the blade damage prediction response surface was formed. Rapid prediction of fan blade damage due to bird strike was realized, and the prediction workflow based on the response surface method was established. Combined with the requirements of the turbofan engine bird ingestion test, the fan structural features, and the results of the conducted bird strike component tests, the blade damage prediction response surface was established, the radial bending, and chordwise bending of two bird-strike schemes were preliminarily identified, and the tearing ranges do not exceed 0.3867 and 0.3941, respectively. The tearing and the chordwise bending are significantly correlated, showing a parabolic trend. The results show that the predicted damage is within the acceptable safety level. The prediction method can identify the damage range and trend. It can provide quantitative technical support for subsequent work such as bird strike finite element simulation, test planning, safety analysis, fan blade bird strike resistant design, etc.Key words:bird strike; fan blade; damage; response surface method; aeroengine0 引言涡扇发动机运行中可能遇到不同类型的外物侵入进而造成损伤,称为可能损伤飞机/发动机的外来物质、碎屑或物体(Foreign Object Debris,FOD),分为软体和硬体2类。

2024-T3铝合金动力学实验及其平板鸟撞动态响应分析

2024-T3铝合金动力学实验及其平板鸟撞动态响应分析

2024-T3铝合金动力学实验及其平板鸟撞动态响应分析刘富;张嘉振;童明波;胡忠民;郭亚洲;臧曙光【摘要】通过电子万能试验机和分离式霍普金森拉杆(SHTB)拉伸试验分别获得2024-T3铝合金材料准静态和高应变率两种应变率下的应力-应变曲线。

铝合金材料的本构关系由能够反映材料硬化效应和应变率强化效应的John-son-Cook材料模型描述,方程中的4个参数通过不同应变率下的应力-应变曲线拟合得到。

基于瞬态动力学软件PAM-CRASH,结合材料动态力学性能试验所获得的2024-T3铝合金Johnson-Cook模型方程,耦合光滑粒子流体动力学(SPH)方法和有限元(FE)方法建立2024-T3铝合金平板的鸟撞数值模型,数值计算所得动态响应与鸟撞试验结果吻合较好,表明建立的鸟撞数值计算模型是合理、可靠的,整个分析流程从材料动态力学性能试验、鸟撞数值计算到最终的鸟撞试验验证为飞机结构的抗鸟撞设计与分析提供了有力的参考。

%Stress-strain curves at quasi-static and high strain rates of 2024 -T3 aluminum alloy were measured using dynamic tensile tests with an electronic universal testing machine and a split Hopkinson tensile bar (SHTB).The constitutive model of this aluminum alloy was described with Johnson-Cook model reflecting the material hardening effect the strain rate strengthening effect.4 parameters in Johnson-Cook equation were obtained by using stress-strain curves fitting at different strain rates.By combining the smooth particle hydrodynamics (SPH)method and the finite element (FE)method,and using Johnson-Cook equation of 2024-T3 aluminum alloy,the numerical model of bird impact on a 2024-T3 aluminum alloy square plate was established with the transient dynamic software PAM-CRASH.Dynamic responses of the platewith numerical calculation agreed well with the test results of bird impact.It was shown that the numerical calculation model is reasonable and reliable.The whole analysis process including material dynamic tests, numerical calculation and bird impact tests provided a reference for anti-bird impact design and analysis of aircraft structures.【期刊名称】《振动与冲击》【年(卷),期】2014(000)004【总页数】6页(P113-118)【关键词】动态响应;数值计算;本构模型;鸟撞;试验【作者】刘富;张嘉振;童明波;胡忠民;郭亚洲;臧曙光【作者单位】中国商飞北京民用飞机技术研究中心,北京 102211; 南京航空航天大学飞行器先进设计技术国防重点学科实验室,南京 210016;中国商飞北京民用飞机技术研究中心,北京 102211;南京航空航天大学飞行器先进设计技术国防重点学科实验室,南京 210016;中国商飞北京民用飞机技术研究中心,北京 102211;西北工业大学航空学院,西安 710012;中国建筑材料检验认证中心,北京100024【正文语种】中文【中图分类】O347.3鸟撞是发生在毫秒量级内的瞬态冲击现象。

鸟撞风挡问题分析方法研究

鸟撞风挡问题分析方法研究
SongBotao,MiaoZhitao,ChiFujian.StudyonBirdStrikeAnalysisMethodforAircraftWindshield[J].AdvancesinAeronauti灢 calScienceandEngineering,2019,10(S1):50灢54.(inChinese)
文 章 编 号 :1674灢8190(2019)S1灢050灢05
Vol灡10 Suppl灡1 May 2019
鸟撞风挡问题分析方法研究
宋 波 涛 ,苗 志 桃 ,池 福 俭
(航 空 工 业 西 安 飞 机 工 业 (集 团 )有 限 责 任 公 司 ,西 安 暋710089)
摘暋要:因鸟撞问题引起的飞机结构破损等损失问题需要解决,采用对比分 析 的 方 法,对 目 前 鸟 撞 风 挡 问 题 的 地面试验法、工程计算法及有限元仿真法进行详细介 绍,比 较 不 同 方 法 的 适 用 条 件、范 围 及 其 优 劣。 选 取 某 飞 机风挡层合玻璃作为研究对象,采用有限元仿真法,建 立 风 挡 鸟 撞 模 型,利 用 任 意 拉 格 朗 日 欧 拉 耦 合 法 (ALE) 完成相关分析,结果得到鸟撞风挡的变形及其速度、加速 度 等 参 数,并 总 结 鸟 撞 风 挡 问 题 的 三 种 分 析 方 法 的 优 劣 性 ,为 飞 机 风 挡 设 计 提 供 参 考 。 关 键 词 :风 挡 ;鸟 撞 ;地 面 试 验 法 ;任 意 拉 格 朗 日 欧 拉 耦 合 法 中 图 分 类 号 :v215.2暋 暋 暋 暋 文 献 标 识 码 :A暋 暋 暋 暋 暋 暋 暋DOI:10.16615/ki.1674灢8190.2019.S1.010
增 刊 1暋 暋 暋 暋 暋 暋 暋 暋 暋 暋 暋 暋 暋 暋 宋 波 涛 等 :鸟 撞 风 挡 问 题 分 析 方 法 研 究

某型民机风挡结构抗鸟撞性能分析

某型民机风挡结构抗鸟撞性能分析

某型民机风挡结构抗鸟撞性能分析作者:金日浩来源:《科技视界》2017年第09期【摘要】运用LS-DYNA软件的显式中心差分法求解器进行了大规模数值计算,考察了风挡结构抗鸟撞特性,数值分析结果对于易受鸟撞的部位有一定的工程价值,为结构的鸟撞地面模型试验验证提供有力技术支持。

【关键词】鸟撞;风挡;有限元分析0 引言飞机在飞行运营中会受到外部环境的种种威胁,其中在运营中由鸟撞引起的问题长期以来受到研究人员的关注。

在CCAR-25中,对于风挡的抗鸟撞性能提出了明确的要求:风挡结构在受到1.8公斤重的鸟的撞击后能够完成飞行。

飞机与鸟沿着飞机飞行航迹的相对速度取海平面VC或2450米(8000英尺)0.85VC,两者中的较严重者[1]。

本文针对某型民机的风挡结构,运用有限元分析对风挡结构的抗鸟撞性能进行评估。

1 数值计算模型1.1 鸟体模型在鸟撞模拟中,通常采用中间为圆柱体,两端为半球体的胶囊型结构来模拟鸟体形状,鸟体密度为为0.9~0.95g/cm3,鸟体质量为1.81kg,速度为125m/s,具体尺寸见图1。

由于鸟体在撞击过程中变形极大,会出现鸟体四溅的状况。

为了避免数值计算困难,采用SPH粒子技术来建立鸟体模型,粒子数为35808个。

采用带失效应变的弹塑性流体动力学本构模型模拟鸟体材料,具体参数见文献[2]。

1.2 风挡模型在大型问题的数值模拟计算中,一般采用壳单元来模拟结构。

选用壳单元可以有效的降低求解的规模,在保证计算精度的前提下节省计算时间。

在文献[3]中考察了不同的壳单元尺度对于有限元计算精度的影响。

本文采用的壳单元尺寸为5mm。

本文计算的风挡结构包括:机头蒙皮、风挡玻璃、风挡骨架及顶部纵向加强结构等。

风挡玻璃采用3+2的结构形式,总厚度为28.5mm。

风挡有限元模型中除风挡玻璃用体单元建模外,其它构件(如蒙皮等)均采用壳单元建模。

壳单元网格大小为5mm,体单元在长度和宽度方向的尺寸大小为6mm,厚度方向的尺寸大小与结构的实际厚度有关,共计壳单元844139个,体单元269208个。

基于SPH的航空发动机风扇鸟撞分析

基于SPH的航空发动机风扇鸟撞分析

科学技术创新2019.191概述随着航空技术的迅猛发展,飞机数量和飞行航次急速增多,飞机鸟撞事故的数量也呈现上升趋势[1]。

鸟撞事故一般发生在飞机起飞降落阶段,以及军用飞机低空高速飞行时。

飞机鸟撞事故的严重程度取决于所撞飞机部位、鸟体质量以及鸟与飞机相对撞击速度[2]。

根据统计,发动机风扇叶片和风挡是受鸟撞击概率最大的两个部位[3]。

由于鸟体的冲击力可能会打碎发动机叶片,而鸟在被搅碎之后,遗骸也可能堵塞发动机的管道,在撞鸟后,发动机往往会出现喘振起火,甚至自行停车,因此鸟撞发动机叶片的危害极大[4]。

鸟撞发动机的研究主要有实验和数值仿真方法两种。

早期主要通过实验进行,但这类试验成本很高。

20世纪随着计算机和仿真技术的发展,数值仿真在鸟撞发动机的研究中得到了广泛应用[5]。

鸟撞发动机问题属于高度非线性冲击动力学问题,撞击过程中叶片会产生大变形,而鸟体会呈现碎裂、流变现象[6]。

因此对鸟体建立准确地数值模型是鸟撞数值分析中的难点[7]。

目前主要的建模方法有拉格朗日(Lagrange )法和任意拉格朗日-欧拉(Arbitrary Lagrange Euler,ALE )法等[8]。

光滑粒子流体动力学(Smoothed Particle Hydrodynamics,SPH )法是一种无网格化Lagrange 算法,它将材料处理成一组流体粒子的集合,每个粒子具有自己的速度、能量和质量特征,并描述为一个与流体性质相关的插值点[9]。

该方法避免了Lagrange 方法中的网格畸变、失效及ALE 方法中的传输和重复计算,因而逐渐在流体动力学、侵彻、碰撞等领域得到广泛应用[10]。

根据鸟撞发动机风扇叶片动态响应的特点,本文混合使用SPH 方法和有限元方法,鸟体采用SPH 方法建模,用流动的粒子描述鸟体的大变形、破碎及飞散,其它区域使用有限元Lagrange 方法。

2数值计算模型2.1发动机风扇模型涡扇发动机风扇由叶片和轮毂组成。

基于SPH算法的某型飞机座舱盖的鸟撞数值分析_王意

基于SPH算法的某型飞机座舱盖的鸟撞数值分析_王意

基于SPH算法的某型飞机座舱盖的鸟撞数值分析王意 王冰上海飞机设计研究院结构设计研究部 上海 200232摘 要:在HyperMesh中建立座舱盖和鸟体的数值简化模型,并用SPH(Smoothed Particle Hydrodynamics)算法对不同条件下的鸟撞进行对比分析,得出了以下结论:合理地采用材料和结构布置,使结构刚度协调分布,可使鸟体顺利地滑出结构,带走冲击能量;座舱盖结构的局部刚度过大,反而降低了其抗鸟撞性。

关键词: SPH,座舱盖,鸟撞0 前言鸟撞是指飞机等飞行器与在天空中飞行的鸟类相撞造成飞行事故的简称,它是随着飞机的诞生。

就出现的一种飞行事故。

随着飞机的广泛使用,鸟撞问题也日益突出。

美国交通部、农业部及联邦航空局1990-2008年的统计资料表明,从1990年到2008年之间,美国民用飞机统计到的鸟撞事件共发生89727起,造成经济损失约3.5亿美元[1][2]鸟撞威胁飞行人员的安全,在经济上造成重大损失。

飞机的座舱盖是飞机遭受鸟撞概率较高的部件,其设计要求必须达到能承受一定的鸟撞载荷。

传统的飞机的抗鸟撞能力的研究,主要是基于试验研究,但试验研究周期长、试验费高。

例如在某型军机风挡的研制过程中,全尺寸鸟撞试验就进行了200多次[3],因此 试验研究通常是用于设计验证。

随着算法的不断完善,求解精度的不断提高,在工程实际中,逐渐形成了以有限元数值模拟与鸟撞试验相结合的方法来进行抗鸟撞设计研究。

1 座舱盖结构的鸟撞分析1.1 鸟撞的分析方法鸟撞问题的实质是流体和结构耦合动力响应问题,载荷与变形是相互影响且不断变化的。

鸟撞是一个伴随瞬时冲击载荷、高应变率、大变形、耦合性很强的过程[4]。

在没有商用软件的条件下,可用解耦的方法求解,解耦法认为鸟撞分析中,鸟体不是关心的对象,鸟撞击的结构的响应才是重点研究的对象。

通过各种假设提取鸟体载荷情况,然后当做已知条件加载到鸟撞结构上,即先确定鸟撞的载荷,再求在此载荷作用下的结构响应。

大型民用飞机风挡鸟撞的适航分析与数值仿真

大型民用飞机风挡鸟撞的适航分析与数值仿真
添加标题
具体 较 细 较 粗
添加标题
具体 较 细 较 粗
添加标题
具体 较 细 较 粗
添加标题
滚动计划法
二、甘特图法
甘特图(Gantt chart)是在20世纪初由亨利·甘特创立的,也叫线条图或横道图,是一种用线条来安排生产或工程进度的方法。
1、目标管理(MBO)的涵义
目标管理:用系统方法,通过科学制定目标、实施目标,依据目标进行考核评价来实施管理任务的过程(管理方法)。
从形式上看,目标管理是一种程序和过程。一切管理活动以制定目标开始、以目标为导向,以目标情况作为管理依据。目标贯穿于组织管理活动的全过程。
滚动计划法的基本思想
滚动计划法是一种将短期计划、中期计划和长期计划有机的结合起来。根据近期计划的执行情况和环境变化情况,定期修订未来计划并逐期向前推移的方法。
具体做法:在计划制订时同时制订出未来若干期的计划但计划内容采用:“近细远粗”的办法,近期计划尽可能的详尽,远期计划的内容则较粗,在计划期的第一阶段结束时,根据该阶段计划执行情况的内外环境变化情况,对原计划进行修订,并将整个计划向前滚动一个阶段,以后根据同样的原则逐期滚动。
1
2
3
4
目标
实际进度
报告日期
月份
设计封面
印刷校样
打印长条校样
制图
设计版式
活动 编辑加工
图书出版甘特图
资料来源:斯蒂芬·P·罗宾斯《管理学》,中国人民大学出版社1997年出版
时间
1

3
4
6
5
7
8

1
2
2
1
4
1
2

无人机飞鸟撞击机翼损伤程度预测仿真

无人机飞鸟撞击机翼损伤程度预测仿真

无人机飞鸟撞击机翼损伤程度预测仿真王永虎;吴志坚;杨敏【期刊名称】《计算机仿真》【年(卷),期】2018(035)009【摘要】日益普及的无人机“黑飞”现象频频发生,严重危及航空安全和正常航班运行,无人机撞击民机关键部位损伤程度的研究是航空冲击动力学的一个新型热点问题.此前,各民航大国都致力于外来物(Foreign Object Debris,FOD)尤其是鸟撞和冰雹撞击等,而对无人机撞击的研究成果屈指可数.依据鸟撞适航研究思路,采用有限元仿真软件ABAQUS分别用SPH法和拉格朗日法建立鸟撞和无人机撞击机翼的有限元模型,对二者撞击机翼过程的某个工况下进行损伤程度预测仿真.在相同的冲击初始条件下进行计算,得出了鸟体和无人机冲击机翼损伤耗散能、应变能随时间变化规律和机翼前缘变形规律.通过对比分析二者撞击机翼的损伤机理和失效过程,相同条件下鸟体撞击机翼蒙皮产生塑性变形,而无人机冲击机翼蒙皮会产生撕裂,甚至会穿透蒙皮造成更严重后果,得出了无人机撞击飞机比鸟体撞击飞机具有破坏性的结论.【总页数】5页(P42-45,83)【作者】王永虎;吴志坚;杨敏【作者单位】中国民航飞行学院飞行技术学院,四川广汉618307;中国民航飞行学院飞行技术学院,四川广汉618307;中国民航飞行学院飞行技术学院,四川广汉618307【正文语种】中文【中图分类】TP319.9【相关文献】1.基于ABAQUS的某无人机机翼肋板裂纹扩展仿真 [J], 王海鹏;吴永东2.基于神经网络的机翼气动参数预测仿真研究 [J], 刘昕3.关于液滴撞击弹性固体变形预测仿真 [J], 马小晶; 买买提明·艾尼; 闫亚岭; 金阿芳4.扇翼无人机二维机翼气动特性仿真研究 [J], 王旭;刘明鑫5.基于ANSYS Workbench新型折叠翼无人机机翼仿真优化分析 [J], 左晓军;齐元胜;马克西姆因版权原因,仅展示原文概要,查看原文内容请购买。

平板模拟叶片不同角度鸟撞响应数值研究

平板模拟叶片不同角度鸟撞响应数值研究

摘要:针对航空发动机风扇叶片的抗飞鸟撞击设计需求,建立了钛合金平板模拟叶片受模拟鸟撞击的有限元模型,分析了不同撞击角度下叶片的动态响应特征变化规律。

关键词:风扇叶片;不同角度;鸟撞响应;数值模拟0 引言航空发动机风扇叶片在工作中非常容易受到飞鸟撞击而导致叶片受损,直接危害飞行器的安全运行,而真实的鸟撞发动机试验费用十分昂贵],因此,在试验前采用数值模拟方式,对叶片鸟撞响应进行多种工况的分析和预估具有重要意义。

鸟撞击叶片时,鸟体和弯扭叶片间的相对作用位置随发动机转速和撞击速度变化而变化,可以将其简化为叶片前缘以不同角度受鸟体撞击和切割的工况。

本文建立了平板叶片受不同角度鸟撞击的有限元模型,分析了不同的鸟撞击工况下叶片的动态响应变化规律,为叶片的抗鸟撞击分析提供依据。

1 平板模拟叶片鸟撞建模在瞬态动力学分析软件LS-DYNA中,建立鸟撞TC4钛合金平板有限元模型。

模拟叶片在夹持端加厚,以防止在叶片根部产生过大的集中应力,撞击段为矩形平板,长300 mm、宽100 mm,夹持段采用固定支架约束,平板采用六面体网格划分,撞击区域适当加密,平板力学性能采用文献中的双线性塑性模型;模拟鸟采用长径比为2、直径40 mm的圆柱几何形状,采用SPH方法建立,鸟体力学性能参数采用文献中弹塑性水动力本构模型与多项式状态方程描述。

变化平板倾斜角度以模拟不同的撞击角度,共设计7组撞击角度,角度为20°~90°,鸟撞平板有限元模型如图1所示。

2 模拟鸟撞平板叶片过程以45°撞击为例,鸟撞叶片的过程如图2所示。

0.06 ms时,鸟体接触平板前缘,撞击导致叶片前缘开始产生变形,鸟体被叶片前缘切割成两部分;0.6 ms时,鸟体因撞击损耗过半,传递给平板的能量导致叶片前缘变形达到最大,同时应力波使叶片前缘发生波浪形振动,进而在振动节线位置产生较大的应力;1 ms以后,鸟撞击过程结束,叶片进入持续振动状态,在本身阻尼作用下逐渐停止,叶片前缘保持一定残余变形;50 ms时,叶片振动基本停止,叶片前缘保持一定残余变形。

鸟撞飞机风挡夹层结构数值分析

鸟撞飞机风挡夹层结构数值分析

鸟撞飞机风挡夹层结构数值分析
张博平;孙秦;李玉龙
【期刊名称】《航空计算技术》
【年(卷),期】2008(38)1
【摘要】用非线性显式动力分析软件LS-DYNA建立了风挡夹层结构有限元模型,鸟体用ALE格式和与应变率相关的随动硬化材料模型,风挡夹层用Lagrange格式和双线性随动硬化材料模型,对鸟撞击风挡进行了数值计算,分析了夹层结构应力分布,为飞机风挡设计提供理论依据.
【总页数】4页(P6-8,22)
【作者】张博平;孙秦;李玉龙
【作者单位】西北工业大学,航空学院,陕西,西安,710072;西北工业大学,航空学院,陕西,西安,710072;西北工业大学,航空学院,陕西,西安,710072
【正文语种】中文
【中图分类】O242
【相关文献】
1.鸟撞飞机风挡非线性数值分析 [J], 朱书华;王跃全;郭亮;张立圣;童明波
2.飞机风挡夹层结构耐撞性数值分析 [J], 张博平;孙秦;李玉龙
3.飞机风挡结构抗鸟撞动响应数值模拟 [J], 王新军;岳珠峰;王富生;冯震宙
4.飞机前风挡非对称结构的鸟撞数值模拟 [J], 王猛;黄德武;罗荣梅
5.飞机风挡结构抗鸟撞数值模拟 [J], 高磊;万小朋;赵美英;王文智
因版权原因,仅展示原文概要,查看原文内容请购买。

一种飞机鸟撞试验的有限元分析方法

一种飞机鸟撞试验的有限元分析方法

一种飞机鸟撞试验的有限元分析方法王卫伟【摘要】本文主要介绍飞机鸟撞试验有限元三维仿真模型的建立及有限元分析的方法,模型仿真计算结果与试验结果相比较,具有良好的一致性,从而验证了此分析方法的分析结果是具有一定的参考价值.可降低直接通过试验确定鸟体参数的困难度及高成本.【期刊名称】《科技风》【年(卷),期】2015(000)017【总页数】2页(P57,69)【关键词】鸟撞;SPH;有限元;仿真【作者】王卫伟【作者单位】中航飞机股份有限公司研发中心,陕西西安710089【正文语种】中文1 引言鸟撞是指飞机等飞行器与天空中飞行的鸟类相撞造成的飞行事故。

随着世界航空的快速发展,鸟撞事故的发生越来越频繁,已经成为世界航空运输业的三大灾害之一。

国际航空联合会已经把鸟撞事故升级为“A”级航空灾难。

对于鸟撞问题的研究受到了各国政府和学者的重视,一方面学者的研究集中在如何避免鸟撞事故的发生上,如在机场附近建立高效的驱鸟设施,研究鸟类的迁徙时间和路径以避免与航线冲突等。

另一方面的研究则侧重于如何提高飞行器结构抗击鸟撞的能力上。

然而直接通过试验确定鸟体参数非常困难而且代价较高。

随着计算机技术和有限元数值计算理论的发展,形成了以鸟撞数值仿真分析与试验相结合的方法来进行飞行器抗鸟撞的研究。

本文主要介绍飞机鸟撞试验有限元三维仿真模型的建立及分析,为飞机鸟撞试验的分析研究提供一定的借鉴与参考。

2 模型建立鸟体本构模型与参数作为鸟撞数值模拟计算的基本性能数据是开展抗鸟撞结构设计分析的基础,由于鸟体在碰撞的过程中发生了较大的变形,在这里,我们对鸟体采用SPH方法进行仿真计算。

SPH是一种无网格拉格朗日计算方法,与传统的基于网格有限元法不同,它是用一系列具有物理量的粒子质点来代表整个连续的物体,然后通过一个称为“核函数”的积分进行“核函数估值”,从而求得流场中不同位置在不同时刻的各种动力学量。

由于该方法没有网格,可以避免传统有限元方法中的网格缠绕及扭曲问题,因而被广泛的应用于大型结构断裂及大变形问题的数值计算分析中,本文中模拟的鸟体在撞击飞机后有较大的变形,适用于SPH方法。

基于PAM-CRASH的鸟撞飞机风挡动响应分析

基于PAM-CRASH的鸟撞飞机风挡动响应分析

基于PAM-CRASH的鸟撞飞机风挡动响应分析
刘军;李玉龙;徐绯
【期刊名称】《爆炸与冲击》
【年(卷),期】2009(029)001
【摘要】结合显式动力分析有限元软件PAM-CRASH及其提供的SPH算法,建立了鸟撞飞机风挡数值分析模型.对某飞机风挡进行了动响应分析,计算了风挡中轴线上四点位移时间曲线;鸟撞击过程的仿真结果表明,SPH鸟体模型能有效模拟撞击时鸟体溅射成碎片的情形;建立了鸟撞击作用下风挡破坏判据,对风挡在试验条件下是否破坏进行了模拟计算.计算结果和试验结果吻合较好,表明本文建立的鸟撞飞机风挡数值模型是有效的.
【总页数】5页(P80-84)
【作者】刘军;李玉龙;徐绯
【作者单位】西北工业大学航空学院航空结构工程系,陕西,西安,710072;西北工业大学航空学院航空结构工程系,陕西,西安,710072;西北工业大学航空学院航空结构工程系,陕西,西安,710072
【正文语种】中文
【中图分类】O347
【相关文献】
1.飞机风挡鸟撞动响应分析方法研究 [J], 彭迎风;滕春明
2.飞机前风挡鸟撞动力响应分析 [J], 臧曙光;武存浩;汪如洋;马眷荣
3.飞机风挡鸟撞动响应分析方法研究 [J], 张志林;姚卫星
4.飞机圆弧风挡鸟撞动响应分析 [J], 张志林;张启桥;李铭兴
5.基于PAM-CRASH的鸟撞风挡有限元分析 [J], 刁斌;左洪福;蔡景
因版权原因,仅展示原文概要,查看原文内容请购买。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

( , , , ) 1 . S c h o o l o fA e r o n a u t i c sa n dA s t r o n a u t i c s N a n i n n i v e r s i t fA e r o n a u t i c sa n dA s t r o n a u t i c s N a n i n 1 0 0 1 6 C h i n a 2 j gU yo j g ( , , ) 2. I n s t i t u t eo fA l i e dM e c h a n i c s T a i u a nU n i v e r s i t fT e c h n o l o T a i u a n 3 0 0 2 4, C h i n a 0 p p y yo g y y ( , ) , , ) 3. R e a r c ha n dD e v e l o m e n tC e n t e r C h e n d uA i r c r a f t I n d u s t r i a l( G r o u C o . L t d . C h e n d u 1 0 0 9 2, C h i n a 6 p g p g , , ) ( 4. S c h o o l o fC i v i lE n i n e e r i n n dT r a n s o r t a t i o n S o u t hC h i n aU n i v e r s i t fT e c h n o l o G u a n z h o u 1 0 6 4 0, C h i n a 5 g ga p yo g y g 摘 要 :鸟与飞行中的飞机相撞是飞机结构损坏的重要因素 , 严重时会引发机毁人亡的灾 难 性 事 故 。 对 高 速 低空飞行的军用飞机而言 , 风挡部分抗鸟撞的研究对保证飞行安全尤其重要 。 基于飞机圆弧风挡受鸟体撞击 建立了国产某型军用飞机圆弧风挡及鸟体的 计 算 模 型 , 采用 L 的实验观察 , S D YNA 3 D 中有限元和光滑粒子 流体动力学 ( 耦合的数值分析方法 , 对某飞机圆弧风挡受鸟体撞击的过程进行了数值模拟 。 计算结果得 S P H) / 到了风挡结 构 的 变 形 、 位移和应变等几方面的数据, 与实验结果基本吻合。同时, 给出了5 0 0~6 5 0k m h速 风挡发生破坏的临界 撞 速 、 圆弧风挡经受鸟体撞击时发生破坏的可能位 度范围内的撞击力和应力时程曲线 、 置及其破坏方式 。 最后 , 与鸟 体 采 用 任 意 拉 格 朗 日( 和无网格伽辽金方法( 进 行 了 对 比, 验证了 A L E) E F G) S P H 方法在分析鸟撞问题中的优越性 。 研究结果为风挡的安全设计和研制新机型提供了有价值的数据 。 关键词 :鸟撞 ;飞机风挡 ;光滑粒子流体动力学 ;有限元方法 中图分类号 :V 2 1 5. 2; O 3 4 7. 3 文献标识码 :A : A b s t r a c t T h e c l a s hb e t w e e n f l i n i r d s a n da i r l a n e s f r e u e n t l e s u l t s i nd a m a eo f a i r c r a f t s t r u c t u r e sw h i c h y gb p q yr g c o u l de v e nl e a dt oa i r c r a f t c r a s ha n dd e a t ho f c r e wm e m b e r s i nt h em o s td i s a s t r o u sc a s e s . I n v e s t i a t i o no nt h e g w i n d s h i e l dr e s i s t a n c e t ob i r d i m a c t i s t h e r e f o r e c r u c i a l t oe n s u r e t h e s a f e f l i h t o f h i h s e e d l o w a l t i t u d em i l i p g g p , t a r i r c r a f t . B a s e do nt h ee x e r i m e n t a l o b s e r v a t i o no fb i r d i m a c t o na i r c r a f t a r cw i n d s h i e l d s ac o m u t a t i o n a l ya p p p m o d e lo f a na r cw i n d s h i e l da n db i r d f o r ac e r t a i nt eo f a i r c r a f tm a d e i nC h i n a i se s t a b l i s h e d . I s s u e so nt h ea r c y p w i n d s h i e l d i m a c t e db i r d sa r es t u d i e du s i n m o o t h e dp a r t i c l eh d r o d n a m i c s( S P H) a n dt h e f i n i t ee l e m e n t p yb gs y y , F EM) d i s l a c e m e n ta n ds t r a i no ft h e m e t h o d( r o v i d e db S D YNA 3 D. R e s u l t si n c l u d i n h ed e f o r m a t i o n p p yL gt ,w a r cw i n d s h i e l da r eo b t a i n e dt h r o u hn u m e r i c a l s i m u l a t i o n h i c ha r e ew e l lw i t ht h ee x e r i m e n t a ld a t a . M o r e g g p , u s e f u l c o n c l u s i o n s s u c ha st h ei m a c t f o r c eh i s t o r u r v ea n ds t r e s sh i s t o r u r v eu n d e r i m a c tv e l o c i t i e si n p yc yc p / , t h er a n eo f 5 0 0 6 5 0k m h, t h ec r i t i c a l i m a c tv e l o c i t h i c hc a nl e a dt od a m a eo f t h ew i n d s h i e l d o s s i b l e g p yw g p , d a m a e l o c a t i o n sa n dm o d e so f t h ew i n d s h i e l da f t e r i m a c t e db i r d s a r ed r a w nf r o mt h en u m e r i c a l s i m u l a g p yb , t i o n . F i n a l l ac o m a r i s o no f t h ea r b i t a r a r a n eE u l e r( A L E) a n de l e m e n t f r e eG a l e r k i n( E F G)m e t h o d s y p yL g g w i t ht h eS P Hf o rb i r d i m a c t sv e r i f i e s t h a t t h eS P Hm e t h o di sf e a s i b l ea n de f f i c i e n t f o ra n a l z i n i r di m a c t p y gb p r o b l e m s . T h er e s e a r c hr e s u l t sc a np r o v i d ev a l u a b l ed a t af o rt h ed e s i na n dm a n u f a c t u r eo fr e l i a b l ew i n d p g s h i e l d s . : ; ; K e o r d s b i r d i m a c t a i r c r a f tw i n d s h i e l d S P H; f i n i t ee l e m e n tm e t h o d p yw
第1期ቤተ መጻሕፍቲ ባይዱ
贾建东等 : 用S P H 和有限元方法研究鸟撞飞机风挡问题
1 3 7
1] 2] 研究工作 。 龚尧南 [ 和王礼立 [ 等提出了合理分
( 抗鸟撞的实验方法 , 通过空气炮发射 “ 鸟弹 ” 取活 鸡代 )撞击安装于台架上的全尺寸风挡来模拟鸟 撞过程 。 试件参照飞机上的安置角度安装于实验 , 台架上 , 呈斜 撞 击 姿 势 , 角度为2 安装状态保 7 ° 。鸟 体 以 不 持原 设 计 的 空 间 状 态 ( 如 图 2 所 示) / 同 速 度v 从 3 个 不 同 的 位 置 ( 风挡对称线前1 3 处、 中点以及后弧 框) 分 别 对 风 挡 进 行 撞 击。风 挡厚度 均 为 1 材 料 为 航 空 3 号 有 机 玻 璃, 8 mm, , 鸟重为 1 鸟的撞击速度范围为4 8k 3 9~6 2 7 g / k m h。 采用高速 摄 像 和 高 速 摄 影 装 置 从 不 同 角 度对鸟撞过程进行了拍摄 。 使用动态应变仪测量 激光位移传感器测量位移 。 应变 ,
它 风挡是现代 飞 机 上 极 为 关 注 的 关 键 部 件 , 关系到飞行员和整个飞机的安全 。 对高速低空飞
相关文档
最新文档