加密，仅对导弹轨迹区域附近进行了网格局部加密。具体网格示意图如图8所示。
f龙源期刊网httpwwwqika
comc
本次计算来流马赫数为095，单位雷诺数为7874×106，迎角为0°，飞行高度是8000m。流场求解采用欧拉方程准定常计算，时间步长△t为0002s，求解总时长1s。
423模拟结果
通过网格局部重构的方法生成的网格如图9所示，为了保证导弹分离过程中对导弹流体计算的准确性，需保证计算过程中导弹周围网格的疏密度，因此，在设置时，将导弹周围加密网格区域在i
erbox的范围内，如图9所示导弹周围的网格为导弹初始网格的保留。
初始时刻定常流场计算的压力（压强）分布如图10所示，由于流动是跨声速的，因此，在外挂物中部和尾部区域出现激波。由于吊舱和挂架很近，二者之间有较强的气动干扰，使这一区域的流动比较复杂。图11为质心位移与试验值的对比图。从图11可以看出，采用准定常流动计算方法，实现了AEDC机翼导弹分离标准模型轨迹模拟，Is时间内导弹投放迹计算结果与风洞试验数据基本吻合。
5结论
本文通过采用网格变形与局部重构相结合的网格变形技术，以NACA0012和AFDC为例，耦合求解Euler方程与六自由度方程，采用准定常的方法实现了机翼外挂物分离轨迹数值模拟，证明本文网格变形局部网格重构方法有很好的实用性。AEDC仿真结果与试验值吻合较好，验证了本文多体分离计算方法的准确性。
参考文献
1唐志共，李彬，郑鸣，等飞行器外挂投放数值模拟J空气动力学学报，2009，27（5）：592596
2曾铮，王，叶正寅RBF整体网格变形技术与多体轨迹仿真J，空气动力学学报，2015，33（2）：170177
3Torste
B，LarsTTimeaccurateCFDapproachto
umericalsimulatio
ofstoreseparatio
trajectorypredictio
RAIAA20113958，2011
4MahmoodT，AizudMN，ZahirSAerody
amiceffectsofthestorereleaseo
therollattitudeofawi
gco
figuratio
CTra
so
icFlightProceedi
gsofI
ter
atio
alBhurba
Co
fere
ceo
AppliedScie
cesTech
ology，Islamabad，Pakista
，2011
5聂璐，向锦武，飞机外挂物投放安全性的参数影响分析阴飞行力学，2011，29（2）：2528
6黄鹏飞拉普拉斯加权聚类算法的研究D南京：南京航空航天大学，2009
f龙源期刊网httpwwwqika
comc
7胡勇全基于自动生成控制点的拉普拉斯变形D大连：大连理工大学，2012
8EliasEP，Spyrido
DKCIDtra
so
icstoreseparatio
trajectorypredictio
swithcompariso
towi
dtu
eli
vestigatio
sJI
ter
atio
alJour
alofE
gi
eeri
g（UE），2010，3（6）：538553
9范立钦飞机飞行力学研究工作的回顾与展望J飞行力学，1987（01）r