鸟击解析
鸟击的理论基础
鸟击是什么
老师，鸟击解析是指航空器与鸟碰撞的问题吗？
没错。FAR 25.631 / CS 25.631中规定了航空器的耐鸟击性能要求。发动机吸入、风挡冲击、机翼前缘冲击是主要评估对象。
鸟的建模
鸟在碰撞时会流体式变形（表现得更像流体而不是固体）。在有限元法中：
| 模型 | 特点 | 用途 |
|---|---|---|
| SPH粒子 | 无网格。对大变形有强耐受能力 | 最常见 |
| ALE | 将鸟作为欧拉流体处理 | 高精度 |
| 拉格朗日实体 | 常规实体单元 | 低速冲击 |
用SPH粒子来建模鸟！
SPH粒子之间没有固定连接，所以鸟在碰撞时"散开"的行为可以自然地表现出来。LS-DYNA的SPH＋拉格朗日（机体结构）耦合是业界标准。
碰撞条件
FAR 25.631的典型条件：
- 鸟的质量：1.8 kg（4磅）——中等大小的鸟
- 碰撞速度：V_c（巡航速度）等级。180～250 m/s
- 冲击能量：$E = mv^2/2$ ≈ 30～60 kJ
以180 m/s的速度碰撞…能量太惊人了。
60 kJ虽然比汽车碰撞（180 kJ at 56 km/h）要小，但碰撞面积极其有限，因此会导致局部穿透。风挡和发动机进气口的结构完整性受到质疑。
总结
关键要点：
- FAR/CS 25.631规定义务——耐鸟击性能
- 用SPH粒子建模鸟——表现流体式变形
- LS-DYNA的SPH＋拉格朗日耦合——业界标准
- 180～250 m/s的高速冲击——注意局部穿透
鸟击的冲击力超乎想象
1.8kg的鸟以时速800km撞上航空器时，其冲击荷重峰值可达约150kN。FAR 25.571条规定了1.8kg鸟在飞行速度下的试验义务，2000年代以后，欧美主要厂商在CFRP风挡认证试验中以Abaqus为基础的SPH法模拟已成标准。
鸟击的数值计算方法
SPH法建模鸟
LS-DYNA中的SPH鸟模型：
```
*SECTION_SPH
1, 1.0, 0.0, 0.0, 0, 0
*MAT_NULL
1, 950. $ 密度950 kg/m3（鸟≈水）
*EOS_GRUNEISEN
1, 1480., ... $ 水的状态方程
*INITIAL_VELOCITY_SET
bird, 200000., 0., 0. $ 200 m/s (mm/ms)
```
鸟的密度是950 kg/m³，和水一样？
鸟的身体大部分是水分（70～80%）。在高速变形碰撞时，鸟表现得像流体。用水的状态方程（EOS_GRUNEISEN）来描述压力-密度关系。
结构侧建模
总结
SPH法改变了鸟击解析
传统的拉格朗日有限元法因鸟体网格过度变形而计算中断，这个问题在1990年代后期LS-DYNA实装SPH（粒子法）后得到解决。将鸟体用水等同的状态方程（Mie-Grüneisen模型）建模，以约5000个粒子径0.003～0.005m的粒子，可以再现现实的压力波形。
鸟击的实务应用
鸟击的实务
航空器型式证明（TC：Type Certificate）中鸟击试验是必需的。用有限元法进行事先验证。
实务检查清单
不穿透就合格了吗？
风挡的标准是"不穿透"（乘员安全）。机翼前缘的标准是"能够安全继续飞行"。即使不穿透，大变形也会损伤液压管路等，就会不合格。
发动机风扇试验全球统一
基于EASA CS-E 800规范的鸟击认证中，涡轮风扇发动机一级风扇要射入1.8kg鸟，确认发动机能安全关闭。Boeing 787的GEnx-1B发动机认证前，用ANSYS LS-DYNA进行了200多种入射角、速度条件的参数化解析。
鸟击软件比较
鸟击的工具
选型指南
主要求解器鸟击对应比较
LS-DYNA以SPH、ALE、结冰功能的统一实现在鸟击解析中成为业界标准，在Airbus、Boeing认证解析中事实标准化。PAM-CRASH在薄壁结构接触处理精度高，曾用于欧洲战斗机EF-2000风挡认证。MSC Nastran的非线性SOL 700也采用LS-DYNA技术许可，后发制人提供集成环境。
鸟击的先端研究
鸟击的先端研究
无人机碰撞更加危险
DJI Phantom 4（约1.4kg）以时速900km撞上旅客机，其能量相当于同质量鸟的约1.6倍。机体接近刚体，能量吸收少，2020年以后FAA、EASA为应对无人机普及而加快了无人机碰撞HVI（高速冲击）分析手法标准化的进度。
鸟击故障排除
鸟击的故障
不要漏掉沙漏模式发散
SPH-有限元罚函数接触在鸟击解析中，当荷载振荡增长时往往是沙漏控制不足所致。LS-DYNA设定HOURGLASS TYPE=4（Belytschko法向补强）且系数≥0.1，加上SPH粒子密度保持在结构网格的1/3以下，稳定性可大幅改善。
更详细
错误