求解最优月球软着陆轨道的隐式打靶法

基于间接法思想推导出一种隐式打靶法对月球软着陆轨道优化问题进行了研究。建立月球软着陆轨道归一化系统模型,利用庞特亚金极大值原理将月球软着陆轨道优化问题转化为满足最优必要条件的两点边值问题（TPBVP）,采用一种新的时间变量使两点边值问题的终端时刻固定,同时将终端时刻看作状态变量并引入终端时刻的哈密尔顿函数值作为隐式终端条件,提出一种隐式打靶法对含有隐式终端条件的两点边值问题进行迭代求解,从而得出燃料消耗最优的月球软着陆轨道。仿真结果表明,与直接法和混合法相比,隐式打靶法优化精度高,收敛速度快,实现了月球软着陆过程燃料消耗最优。同时应用隐式打靶法求解不同发动机推力值的最优月球软着陆问题,得到燃料消耗最小的最优推重比,可为月面着陆器下降级发动机选型提供参考。     

连续空间离散时间—Riccati传递矩阵积分法

本文提出了分析非线性转子 -支承系统瞬态响应的一种传递矩阵积分法,它借助Newmark差分公式,建立传递关系,使用的传递向量是 {fe} T,而不是以往的 {fe} T;它用Riccati法计算的是各站位的加速度,而不是以往的位移,提高了方法的数值稳定法。本法具有Riccati法的优点,尤其适合于计算含有非线性元件的复杂转子系统的瞬态响应。      

水下爆源位置诱导舰艇总强度鞭状响应规律

基于速度势函数和 Morison相对速度公式,提出了气泡载荷水动力的计算方法,给出了用 Newmark数值求解法求解水下爆炸气泡载荷作用下船体梁动态响应的方法。经模型试验验证后,针对药包在船体梁模型中部位置的情况,研究了炸药深度、爆距和爆点轴向位置对鞭状运动响应的影响规律。研究发现:在其他条件不变的情况下,结构鞭状运动应变幅值随着周期比先增加后减小;鞭状运动应变幅值和长度比呈反比例函数关系;船体梁鞭状运动响 应应变幅值随着长度比先增大后减小。     

转子系统瞬态响应计算的精细迭代积分方法

提出了计算转子系统瞬态响应的精细迭代积分方法。采用精细积分法计算状态转移矩阵,可以得到齐次方程的精确解;结合轴承所提供的油膜力的特点,提出一种隐式积分法,允许使用较大的步长,并对各矩阵及向量适当分块,使非线形方程组的维数大为减小。此方法无条件稳定,可以保证较高的精度与效率。     

电控旋翼气弹动力学建模研究

研究了电控旋翼的气弹动力学建模方法.以有限元旋翼气弹分析程序LORA01为基础,首先给出了一种能计及襟翼移轴补偿影响的时域综合非定常气动力模型,并将伺服襟翼作为额外引入的质量体计入结构动力学模型;之后结合有限元法,根据Hamilton原理导出由广义力表示的带襟翼桨叶非线性运动方程,采用基于Newmark方法的隐式数值积分法对桨叶运动方程进行求解.最后,利用主动控制襟翼旋翼的试验数据对模型进行了比较验证,证明了所建模型的正确性.     

三浮陀螺仪有源磁悬浮加力系统研究

三浮陀螺仪输出轴采用有源磁悬浮,系统工作在分时控制模式下,控制变量为总周期和周期内加力占空比.加力电源类型的选取以及加力电压幅值的确定对电磁力的大小、磁悬浮的功耗有很大影响,进而影响陀螺温度场分布和陀螺精度.在研究了磁悬浮元件电磁特性之后,给出了确定加力电源的基于磁路分析的Flux电磁仿真法.分析给出直流、正弦和方波这三种电源的加力幅值范围为1.7V～ 2V,并结合试验确定1.8V的方波电压加力为最佳方案.实验表明采用该方法分析加力电压是有效的,同时节省了大量实验探索的时间.     

基于非结构网格的涡轮气热耦合数值方法研究

利用课题组自主开发的三维非结构隐式N-S计算软件CU_Turbo,采用气热耦合计算方法,对MarkⅡ内冷径向涡轮导向叶片、带气膜冷却涡轮导叶MT1的流场和温度场进行了数值模拟。计算过程中,隐式时间推进中Jacobi ans矩阵采用对Roe通量的一种近似方法求解。结果表明,计算值与试验值吻合良好,验证了气热耦合计算方法的实用性和有效性,为涡轮工程设计提供了一种新的计算分析方法;涡轮叶片通道内附面层的不同流动状态及气膜冷却,对当地换热都有很大的影响。     

双方程 k-ω SST 湍流模型的显式耦合求解及其在叶轮机械中的应用简

 双方程k-ω剪切应力输运（SST）湍流模型通常以隐式耦合方式或者显式半耦合/解耦的方式来求解。本文提出了该模型的一种显式耦合应用方法,即通过点隐的方式来处理湍流源项的刚性,并与混合Runge-Kutta时间推进以及当地时间步长、隐式残差光顺等加速收敛技术相结合,从而使得湍流方程可以与流动方程同时求解。为了增强计算的鲁棒性,进一步对湍流变量进行了限制。将所发展的方法用于DLR平面叶栅算例,确认了求解结果的正确性以及刚性的来源。通过对三维NASA Rotor 67的模拟,验证了SST模型的精度;进一步将其与Badwin-Lomax（BL）模型、Spalart-Allmaras（SA）模型对比,发现三者都能正确地捕捉出口参数分布,且SST与SA模型的模拟结果比较一致;对于该算例,SST模型在总温模拟上更具优势,而BL模型在总压分布上与试验值更加接近。     

动能BGK格式模拟高雷诺数定常流的收敛性研究

为克服动能BGK格式计算耗时多、收敛慢的缺点,在格式中引入当地时间步长、隐式LU-SGS方法和多重网格技术,基于RAE2822翼型粘性跨声速绕流的数值模拟,对其收敛性的改善进行了研究。结果表明:上述3种加速手段对求解定常问题时的收敛性有明显的改善,达到定常解所需的迭代步数随着CFL数的提高而减少;因而求解某类问题时,隐式格式比显式格式更具优势;多重网格既适用于显式格式,也适用于隐式格式,加速性更加突出。     

求解N-S方程的一个简单隐式算法

1.引言 近十几年来,求解分离流动的N-S方程的数值解法,取得了很大进展。显式有限差分解已发展得较为成熟。但由于稳定条件的限制,获得稳定解的计算时间较长。为了克服这一困难,发展了隐式方法,如隐式因式分解法和交替方向的隐式方法。但是这些方法都要求解三对角矩阵方程,同显式方法比较,程序复杂,每步计算时间较长。为了避免三对角矩阵迭代的困难,MacCormack提出了新的隐式算法,其计算迭代过程为两对角矩阵追杆。本文在[1]的基础上,提出了一种更为简单的隐式算法,其迭代过程     

梁/液体耦合振动的边界积分对称列式法

 利用格林积分导出的液体内部扰动压力函数作为测试函数,通过变分得到一组具有对称系数的代数方程,将该方程同梁的有限元振动方程联立后,即可求解梁在液体中的耦合振动问题。本方法具有未知数目少,数据存储和计算效率高的特点。     

疲劳裂纹扩展速率统计分析

给出了4种常幅载荷下疲劳裂纹扩展速率的统计处理方法,获得了在规定存活率P下的“P-da/dN-△K”关系,并对LC9中心裂纹试样试验数据进行了统计分析,给出了存活率为0.50、0.90、0.95和0.99时Paris方程中的c和n值。最后对4种方法进行了比较。     

Motycka估算进气最大瞬时畸变方法的改进

一、引言 用Motycka估算方法去估算最大瞬时畸变值很费机时,特别是当采样频率很高,估算时间较长时。为了工程应用,有必要改进这种估算最大瞬时畸变值的方法,以减少估算机时。     

结构动力响应精细积分级数解的并行计算

讨论了精细积分法的存储需求问题,指出级数解更适合计算大型结构的动力响应分析.研究了精细积分法级数解的并行计算,精细积分法的时程积分公式包含两项,采用第一项串行、第二项并行的并行算法对精细积分级数解实施并行计算,给出了相应的实现流程.讨论了并行算法的负载分配策略,为减少处理器等待初值的时间,对时间步数实施非平均分配.算例表明,并行算法具有好的加速比.     

圆管孔边角裂纹剪切型应力强度因子的封闭解

 在飞机全动平尾大轴中,经常出现孔边角裂纹。本文提出了一种新的工程-解析方法,导出了飞机平尾大轴孔边角裂纹剪切型应力强度因子的封闭解。这种能量释放方法非常节省机时与人力,可以获得大量系统的新结果。     

铆接加筋板应力强度因子计算中几个问题的研究

给出了采用J积分法计算加筋板应力强度因子时选择积分回路的方法;同时给出了采用刚度导数法计算加筋板应力强度因子时刚度导数的计算方法;另外,还对模拟铆钉连接作用的3种方法进行了初步研究。     

口径天线绕射场的研究

介绍口径天线近场绕射计算的卷积积分法 ,该方法可用 FFT实现快速计算 ,故它是比传统的直接积分或球面波展开法更有效的一种计算模式。在此基础上 ,对圆形和矩形口径天线的近场绕射特性进行了研究。由计算结果可见 ,口径天线在观察面上的绕场分布与口径形状、口径尺寸 ( D/λ)和观察面至口径有关。还对减小口径边缘的近场绕射问题进行了分析     

基于一种新的弯曲三角形元素的薄板极限分析方法

 本文提出了平面薄板极限分析中的一种新的弯曲三角形元素,简化了元素刚阵的形成,大大降低了结构刚阵阶数。用塑性应变模态思想对平面薄板结构进行了极限分析。本文方法具有物理意义明显、数学形式简洁、元素简单、计算量小、实用性强等优点,在若干算例中均取得了满意的结果。