华南理工大学学报(自然科学版) ›› 2009, Vol. 37 ›› Issue (12): 146-151.

• 力学 • 上一篇    

基于SPH方法的鸟撞飞机风挡的数值模拟

李志强1  韩强1  杨建林2  赵隆茂2  姚小虎1   

  1. 1.华南理工大学 土木与交通学院, 广东 广州 510640; 2.太原理工大学 应用力学与生物医学工程研究所, 山西 太原 030024
  • 收稿日期:2008-10-28 修回日期:2009-03-09 出版日期:2009-12-25 发布日期:2009-12-25
  • 通信作者: 李志强(1973-),男,博士,讲师,主要从事冲击动力学和数值模拟研究. E-mail:zqli_tyut@163.com
  • 作者简介:李志强(1973-),男,博士,讲师,主要从事冲击动力学和数值模拟研究.
  • 基金资助:

    中国博士后科学基金资助项目(20090450869);山西省青年科技基金资助项目(200703005)

SPH-Based Numerical Simulation of Aircraft Windshield Under Bird Impact

Li Zhi-qiang Han Qiang1  Yang Jian-lin2  Zhao Long-mao Yao Xiao-hu1   

  1. 1. School of Civil Engineering and Transportation, South China University of Technology, Guangzhou 510640, Guangdong, China ; 2. Institute of Applied Mechanics and Biomedical Engineering, Taiyuan University of Technology, Taiyuan 030024, Shanxi, China
  • Received:2008-10-28 Revised:2009-03-09 Online:2009-12-25 Published:2009-12-25
  • Contact: 李志强(1973-),男,博士,讲师,主要从事冲击动力学和数值模拟研究. E-mail:zqli_tyut@163.com
  • About author:李志强(1973-),男,博士,讲师,主要从事冲击动力学和数值模拟研究.
  • Supported by:

    中国博士后科学基金资助项目(20090450869);山西省青年科技基金资助项目(200703005)

摘要: 风挡抗鸟撞是飞机安全飞行的重要保证.文中基于飞机圆弧风挡受鸟体撞击的实验观察,建立了国产某型军用飞机圆弧风挡及鸟体的有限元分析模型,利用光滑粒子流体动力学法(SPH)耦合有限元法对圆弧风挡受鸟撞击的过程进行了数值模拟,计算得到风挡结构的变形、位移、应变、撞击力、应力、临界撞速、发生破坏的可能位置及其破坏方式等几方面的数据,并考察了SPH粒子疏密对计算结果的影响.研究表明,数值模拟结果与实验结果基本吻合;鸟撞整个过程约4ms,撞击中点、前1/3处和后1/3处,风挡发生破坏(包括安全破坏)的临界撞速分别约为(540±5)、(600±5)和(470±5)km/h;鸟撞过程中,风挡的位移与其厚度是同一量级,风挡的最大应变已达到1012量级;风挡首先发生破坏的位置在后弧框附近,然后向与风挡中线成45。角的方向发展;SPH粒子数越多,鸟体变形模态越好.

关键词: 鸟撞, 飞机风挡, 数值模拟, 光滑粒子流体动力学法, 有限元法

Abstract:

Bird impact resistance of windshield is most important for flight safety. In this paper, based on the experimental observation of aircraft arc windshield under bird impact, a finite element analysis model of arc windshield and bird body for a certain aircraft made in China is established, and the impact process is numerically simulated by means of the Smoothed Particle Hydrodynamics (SPH) coupling with the finite element method. Then, a series of data such as deformation, displacement, strain, impact force, stress, critical impact velocity, possible damage location, and mode of arc windshield are obtained. Moreover, the effect of SPH particle density on the calculation results is analyzed. It is found that ( 1 ) the simulated results basically accord with the experimental ones ;(2) the whole process of bird impact lasts about 4ms, and the critical impact velocities of windshield damage ( including safety damage) are (540 ±5), (600 ±5) and (470 ±5) km/h at the middle point, the front 1/3 and the back 1/3 of the windshield, respectively; (3) the displacement and thickness of windshield in the impact process are of the same order of magnitude, and the maximum strain of windshield reaches 10-2 order of magnitude; (4) the damage of windshield first locates near the back arc frame, then expands towards 45°direction with respect to the windshield midline ; and (5) the deformation mode of bird body becomes better with the increase of SPH particles.

Key words: bird impact, aircraft windshield, numerical simulation, Smoothed Particle Hydrodynamics, finite element method