高超声速飞行器的气动力工程计算

为满足高超声速飞行器在概念设计和优化设计中对气动力快速计算的需求,基于C/C++研究和开发了高超声速飞行器的气动力快速计算程序。使用该程序对HL-20升力体和双椭球模型的气动力特性进行计算,将计算结果与风洞实验数据进行对比分析。结果表明,计算结果与风洞实验数据吻合良好,该程序能够准确计算高超声速飞行器的升阻特性。通过修改工程计算方法,对比计算数据,验证本文使用的计算方法的精度,同时也体现了该程序可以灵活选择适合的计算方法的特性。此外,该程序不仅可以自由的更换计算方法,还能够作为C++的库文件链接成为优化程序的一部分,实现高超声速飞行器的概念设计和优化设计。     

橡胶密封圈在回弹过程中的密封性能分析

本文通过非线性有限元软件MSC.Marc,计算了O形密封圈在密封槽产生短时间隙时的回弹过程。计算分别采用动力学计算和静力计算,并将两种计算结果进行了比较。本文还通过用户子程序考虑了外界压力对回弹过程的影响。     

空天飞机的计算流体力学分析

本文分析了空天飞机空气动力学的特点和计算流体力学在发展空天飞机中的作用,讨论了发展空天飞机对计算流体力学的要求。强调了计算流体力学、地面试验和飞行试验应该相互结合成一个整体。在这基础上,探讨了利用地面试验对计算程序进行“确认”、“校准”和“鉴定”以及鉴定计算程序对地面试验的要求。本文还从计算方法、外流、进气道流动、喷管流动、燃烧室流动等方面,综述了空天飞机计算流体力学分析的最近进展和尚待解决的问题。最后对发展我国空天飞机计算流体力学提出了建议。     

一种飞机升降舵需用功率的计算方法

为提供舵面需用功率的计算依据,对工程作动器的设计提供理论支撑,提出了一种飞机升降舵需用功率的计算方法,通过不同类型飞机实例进行计算。首先,针对不同种类的工况进行阐释,并详述速率计算所用的评价准则;然后,提出带载荷平均速率计算方法步骤,分析铰链力矩特性;最后,通过计算空载速率,分析舵面需用功率,并进行计算分析。计算结果表明,该方法计算出的需用功率不仅保证稳定裕度,且与舵面载荷具有理论相关性,设计结果安全合理,对工程设计具有实用指导意义。     

径向基点插值法计算效率的改进方法

伽辽金弱形式和径向基点插值法(Radial basis point interpolation method,RPIM)的无网格法在解决偏微分方程问题中表现出良好的性能,但是在同时提高计算效率和精度方面存在困难。为了提高此类无网格法的计算效率,本文定义了一种基于背景网格的定义域,在计算定义域内的积分点插值时采用同一批节点,在插值计算过程中减少了部分矩阵计算次数,降低了RPIM无网格法的计算时间。在提高计算精度方面,本文提出一种杂交应力的无网格方法,用Hellinger-Reissner(H-R)变分原理推导求解方程,采用无网格方法求解。数值算例表明,本文方法计算二维固体力学时,在具备良好的计算精度的同时提高了计算速度,具有较高的实际应用价值。     

亚、超音速细长翼身组合体大迎角气动特性的计算研究

本文提出了一种适于初步设计使用、具有良好精度的亚、超音速细长翼身组合体大迎角气动特性的综合性计算方法。对大迎角情况下的涡升力,采用吸力比拟原理计算;位流升力的计算,采用基本解的数值计算方法。关于机翼翼剖面头部圆度和涡破碎对涡升力的影响,进行经验性修正。翼身干扰的贡献,通过翼身干扰系数进行计算。并按文[4]原理,将亚音速计算方法推广到亚音速前缘的超音速情况。对几种机翼与翼身组合体的计算结果表明,本文方法具有方法简便、计算快速和符合设计精度要求的优点。     

安全起飞和中断起飞的计算及实例分析

起飞性能的分析是确保航空公司安全和高效运行的关键。安全起飞和中断起飞的计算是起飞性能分析的主要内容。本文依据民航相关法规,提出了计算不同可用起飞距离和所需起飞距离的分段计算方法,给出了满足安全起飞和中断起飞的要求以及计算场长限重的工程方法。并以某机型在某机场跑道起飞为例,计算了各个可用起飞距离和所需起飞距离,给出了某条件下场长限重计算结果。     

基于CBS有限元的流动-传热-变形耦合计算方法

研究了基于CBS有限元法和常规有限元法相结合的流动-传热-变形耦合计算方法。在该方法中,流体流动和传热采用CBS有限元方法计算,固体变形采用常规的有限元方法计算,实现了流体域和固体域统一的有限元网格划分,简化了变形过程中的网格生成和不同网格间的数据交换。然后,依据此方法发展了计算程序,并通过算例分析,校验了计算方法的可行性和程序的计算能力。     

轴测投影基本定理的计算理论和方法

波克－许华尔兹定理是轴测投影的基本定理。本文研究了轴测投影基本定理的计算问题,给出了从投影图形计算空间四面体位置及运动参数的计算方法和公式;同时给出了从空间四面体计算投影图形位置及运动参数的计算方法和公式,据此建立了新的、含有定量关系的轴测投影基本定理。     

弹丸气动力特性计算方法

本文介绍了制导与非制导单独强身亚、跨、超声速有攻角时的升力、阻力特性及压力中心的计算方法。其中超声速波阻采用 Van Dyke 二级扰动理论计算,跨声速波阻用半经验方法处理。其余如摩阻、底阻、粘性分离和旋转带的阻力均用经验方法计算。无粘升力特性,在超声速时用一级横流理论计算,亚、跨声速时采用半经验方法处理,粘性升力特性则全部用经验方法计算。这套方法已软件化,用一个源程序表出。计算结果表明,该方法对尖头、截平头和半球形头部的实际弹形均有相当好的升力、阻力特性计算精度,并能给出合理的压力中心计算结果。     

再入端头帽热应力计算方法的改进

针对位移法求解热应力问题时应力计算精度较低的缺点,本文从广义变分原理出发,导出了应变加权积分法。它在位移法的基础上可以显著地改善应力计算精度,可以实现用较粗的网格划分获得较高的热应力计算精度。本文提出的方法可以看成是位移法的一种新的应力计算方法,因此可以很方便地引入现有位移法端头帽热应力计算程序。     

基于Euler方程解摄动的超声速旋转导数计算方法

基于笛卡尔网格和采用Euler方程基本流场解,通过Riemann不变量摄动,对旋转扰动引起的物面法向速度构建摄动关系,发展了一种超声速旋转导数快速计算方法。文中给出了具体的推导过程,可以看出由此得到的超声速旋转导数,仅与定常欧拉解基本流相关。因此,一旦CFD计算得到定常欧拉解,就可一次性地快速计算出感兴趣的超声速旋转导数。为了验证该计算方法,本文选用国际动导数标准模型Basic Finner算例进行了计算,计算结果与试验数据、文献数据具有很好的一致性,展示出本文方法计算耗时少、精度高、工程实用性强等优点,且该算法基于笛卡尔网格,特别适合复杂气动外形飞行器超声速旋转导数的快速计算。     

新复合材料夹芯壁板式组合结构等参元计算

飞机副翼采用碳纤维/环氧,Nomex纸蜂窝夹芯壁板与碳纤维/环氧实心梁、肋空间组合结构在国内还是首次。作者结合某改型机副翼的设计进行了应力应变有限元的分析计算,编制了一套通用于翼面的夹心壁板组合结构八节点等参元计算程序,按用于微机IBM-PC及其兼容机编制,从而使计算费用大为降低,将计算结果与原结构数据分析比较证明,计算所得的应力应变分布规律正确,适用于高性能的新型复合材料夹芯结构计算。     

防范秘密攻击的安全计算的博弈论实现

在博弈论中,惩罚博弈模拟了参与者试图欺骗但又不想被抓住,即安全计算中秘密攻击者的情形。针对密码学的计算博弈模型,本文对Halpern与Rafael提出的能否在计算具有成本的惩罚博弈与具有一定威慑度的防范秘密攻击的安全计算之间建立联系的问题给出肯定的回答,提出威慑度为1/2的防范秘密攻击的安全是计算博弈中错误可忽略的调解人的通用实现。     

飞机结冰中水滴撞击特性欧拉法的网格影响分析（英文）

过冷水滴撞击特性计算研究是飞机结冰预测与防除冰系统性能分析的基础,欧拉法是过冷水滴运动及撞击特性计算的常用方法之一,分析欧拉法在计算复杂表面水滴运动撞击时的准确性具有重要意义。以NACA0012翼型、冰风洞风道、S形进气道及三段翼为研究对象,利用欧拉法采用不同网格计算获得了水滴运动及局部水收集系数,结果表明：当水滴运动没有受到上游部件的导流或者遮挡时,欧拉法计算结果具有网格无关性;当水滴运动受到上游部件影响发生轨迹偏转时,欧拉法计算结果具有网格相关性,不同网格计算得到的液态水含量及局部水收集系数不一致,需要进一步考虑水滴运动受到上游效应影响时网格对欧拉法计算结果的影响。     

大型太阳帆板展开状态动态特性研究

利用有限元技术及自编的 FORTRAN程序计算了单翼太阳帆板约束模态的振型和频率 ,计算结果和实验结果比较表明 ,计算精度在工程上是可以接受的。在此基础上 ,以飞船带有双翼对称太阳帆板为模型 ,计算了太阳帆板在转动和不转动情况下非约束模态的振型和频率。上述计算结果为研究带有挠性太阳帆板的航天器的动力学与控制问题打下了基础     

超长大直径钻孔灌注桩桩端承载力研究

阐述了超长大直径钻孔灌注桩桩端承载力的影响因素,包括施工工艺、深度效应和尺寸效应等。分析和对比了目前国内常用的各规范所采用的计算方法。文中算例的计算结果表明,对于同等条件的桩,各规范的计算结果相差较大。采用自平衡试桩法获得了几个重大工程的超长大直径钻孔灌注桩的桩端承载力的实测值。将工程实例的测试结果与各规范的计算结果进行对比,证实了目前规范的计算方法不尽合理。目前各规范中桩端承载力的计算方法是在针对小直径的中短桩的研究基础上建立的。如今,大型桥梁工程中普遍采用超长大直径钻孔灌注桩的情况下,对其桩端承载力的计算应当进行相应的调整,使设计更为合理。     

布撒器雷达截面计算与测量

用物理光学法计算了布撒器的雷达截面,并与外场测量结果进行了比较。计算的RCS～θ曲线与测量的RCS～θ曲线变化趋势相同,在典型方向上的计算结果与测量结果比较吻合。表明镜面回波在强度上占绝对优势的目标RCS计算中,物理光学法可以给出工程上可用的结果。     

大机动无人机飞行轨道计算中两个问题的分析

本文讨论无人机作大机动飞行时,轨道计算中必须考虑的两个问题。 一个问题是按照大机动无人机的油箱结构,考虑燃料消耗和前后左右的波动所引起的无人机重心位置的变化。这个变化引起的作用在无人机上的力矩的变化对无人机进入和退出大机动时的飞行参数有明显的影响。另一个问题是da/dt的计算,文中推导了准确计算的公式,按此公式考虑无人机转动受下洗气流的影响,可以准确地计算出阻尼力矩的影响。 在轨道计算中计及这两个问题后,不仅可以计算出无人机在大机动飞行中高度和速度的波动,而且可以计算出滚转角的波动,这和实际飞行情况是比较接近的。     

航天飞机在再入段的亚音速空气动力特性的一种有效算法

本文采用吸力比拟原理,结合基本解的数值计算方法,用来计算航天飞机机翼从小迎角到大迎角(a=0°～30°)的亚音速纵向气动特性;而对零升阻力和机身气动特性,则用工程估算方法计算。由于目前的航天飞机,一般为下单翼的复杂外形翼-身组合体,根据文[9]的原理,可忽略翼-身干扰对纵向气动特性的影响。 本文导得可以计及涡效应的任意平面形状边条机翼的亚音速气动特性的计算公式,亦可计算尖梢机翼的展向升力分布。公式中所需的位流系数可采用涡格面元法进行数值计算来获得,压缩性效应则通过位流系数来计及。 本文计算了多种机翼和航天飞机的气动特性。与实验数据比较表明,本方法具有方法简便、计算快速和计算结果具有设计精度的优点,是计算航天飞机亚音速气动特性的一种有效方法。可供航天飞机初步设计使用,亦可作为航天飞机气动优化设计系统中的子系统。经过适当推导,本方法可推广应用于亚音速前缘的超音速情况。