高超声速气动力及导数计算报告.doc

高超声速气动力及导数计算报告 撰写人： 学 号: 班 级： 2012年10月25日 一 引言 实验目的： ①熟练运用面元法中有关网格划分的方法； ②掌握高超音速气动力工程估算方法中面元法及牛顿法；同时，比较两者的计算结果，并分析差异产生原因。 实验条件： 编程计算该旋成体的升力、阻力、升阻比及俯仰力矩系数，还有导数作出曲线。 图（1） 计算条件： 几何参数： 飞行器运动状态： 二 计算方法 规定导弹的体坐标系为：轴沿着导弹纵轴向后，轴垂直于弹体向上，轴于其它两轴构成右手坐标系，即指向导弹左侧，原点位于导弹前缘点。取原点为参考点。 牛顿法： 牛顿法的假设如下： ①攻角α小于物面倾角θ； ②假设流体由大量均匀分布的，彼此独立无相互作用的质点所组成，它们排列整齐、平行地沿着直迹线流向物体。 ③流体质点流与物面碰撞时，流体质点将失去与物面垂直的法向动量，而保持原有的切向动量沿物面向下流下去。由于法向动量的变化从而引起流体作用在物体上的力。 ④流体对物面的压力只作用在物面能与流体质点相碰撞的表面上，而遮蔽区上压力为零。 牛顿公式： 其中为来流速度方向与物面切面的夹角。 由课本可知： 对于圆球部分，由于用于验证，我们假设αθ,: 对于圆锥部分，由于每一部分圆锥各自的一直为常值，因此直接带入积分进行计算，仅与有关，当时，.其余为零。由于用于验证，我们假设αθ,因此可得对于圆锥其相关系数为： 对于圆柱部分，由于，因此，仅计算其法向力系数，可得 于是 对于本实验： 面元法： 1）确定平面面元 小曲面元四个角点的坐标为。其中，两个对角点形成矢量，分别把它们标记为 这里 通过两个矢量的矢量积，可以得到一个新的矢量，把它记为： 的单位矢量定义为： 我们知道：由一个法向矢量和空间任一点可以确定一个平面。因此，把矢量作为法向矢量，并选取曲面的四个角点的平均值作为一个空间点,从而可以确定一个平面，这就是我们要选取的面元所在的平面。其中为面元的外法线方向单位矢量，由如下公式计算： 现在，将四个角点沿矢量的方向投影到面元平面上，可以求得四个投影点，其坐标分别为 式中， 2）建立面元坐标系 在面元形成后，需要建立一个面元坐标系。现在我们选取面元法向单位矢量作为一个单位矢量，再把的单位矢量作为一个单位矢量： 由和的矢量积可以决定第三的单位矢量： 至此，已经建立了面元坐标系，原点取，轴平行于，轴平行于，轴平行于。易知轴和轴在面元平面上。 3）面元坐标系与飞行器参考坐标系之间的变换 在具体计算中，需要把飞行器坐标系中的坐标值变换到面元坐标系中去，也需要把面元系中的值转换到飞行器坐标系中去。此时，转换公式为： 飞行器坐标转换到面元坐标 从面元坐标系转换到飞行器坐标系 4）面元的面积和质心坐标 现在把投影点在飞行器坐标系中的值变换到面元坐标系中去，在面元坐标系中的值用来表示，则 在面元平面上，由凸四边形面积公式可求得面元的面积为 这里由于平行于轴，因此，面积公式可以化简为 在面元坐标系中，面元质心的坐标为 利用面元坐标系转换到飞行器坐标系的转换关系，可求得质心在飞行器坐标中的质心坐标 5）面元冲击角 在气动力的计算中，我们必须要知道平面面元与自由流速度的夹角，即冲击角。而第个面元平面的冲击角： 式中为第面元平面的单位法向向量。 速度矢量 式中，总角速度 向径 ： 为来流速度，为物体滚转、偏航、俯仰角速度，为参考点的坐标，在编程过程中取（0，0，0）。 6）计算面元压力系数 根据冲击角是否大于零，可以把面元所处的流场分为迎风面和背风面。 若冲击角大于零，则可利用修正牛顿理论：，计算面元的压力系数。 其中为来流速度与面元内法线的夹角 若冲击角小于零，则可直接令。 7） 面元上气动力系数及动静导数计算 将飞行器表面分成若干面元后，用积分的方法可以得到它的气动力系数及其静、动导数。 参考坐标系中（轴沿着导弹纵轴向后，轴垂直于弹体向上，轴于其它两轴构成右手坐标系，即指向导弹左侧，原点位于导弹前缘点。），来流速度矢量，假设面元外法向矢量，则来流速度在面元内法向上的投影为： 分别对和求导，得 物体滚转、俯仰、偏航角速度为，由三个角速度引起的物体与气流相对运动的三个速度为： 其中表示面元质心到参考点的距离。 令在内法向的投影为，若记对的导数为 式中—参考面积。 总速度在面元法向的投影：，则修正的牛顿理论亦可表示为。 根据上述计算，压力系数的静、动导数计算公式为： 根据前面计算的面元上的几何参数及压力系数和静、动导数，可用下列公式计算面元上气动力、力矩及其静、动导数。 面元上所受的气动力： 面元上气动力对参考点的力矩： 面元上气动力的静导数： 面元上气动力产生的力矩的静导数： 面元上气动力的动导数: 面元上力矩的动导数: 求出每个面元上的气动力、力矩及其静动导