【摘 要】
:
轨道动力学快速计算是航天工程中的基础问题,广泛存在于轨道设计、空间抓捕以及深空探测等任务中.基于有限差分原理的经典数值积分算法,由于精度严重依赖小积分步长,难以满足航天器在轨快速计算需求.针对该问题,提出一种局部配点反馈迭代算法,该算法能高效解算受到初值约束和两点边值约束的轨道动力学方程.基于Picard迭代公式建立数值算法以避免计算雅可比矩阵,引入误差反馈以加快迭代收敛速度.通过时域配点法消除迭代公式推导中的复杂符号运算,使迭代公式更为简洁.结合拟线性化法及叠加法,算法能对两点边值约束下的Lambert
【机 构】
:
西北工业大学航天学院, 西安 710072;航天飞行动力学技术国家级重点实验室, 西安 710072
论文部分内容阅读
轨道动力学快速计算是航天工程中的基础问题,广泛存在于轨道设计、空间抓捕以及深空探测等任务中.基于有限差分原理的经典数值积分算法,由于精度严重依赖小积分步长,难以满足航天器在轨快速计算需求.针对该问题,提出一种局部配点反馈迭代算法,该算法能高效解算受到初值约束和两点边值约束的轨道动力学方程.基于Picard迭代公式建立数值算法以避免计算雅可比矩阵,引入误差反馈以加快迭代收敛速度.通过时域配点法消除迭代公式推导中的复杂符号运算,使迭代公式更为简洁.结合拟线性化法及叠加法,算法能对两点边值约束下的Lambert问题实现高效解算.基于ph网格细化法建立计算参数自适应调节算法,能进一步增强局部配点反馈迭代法的大步长计算优势.通过求解二体动力学模型下的递推轨道,摄动Lambert问题以及限制性三体动力学模型下的转移轨道验证了算法有效性.仿真结果表明,在相同计算精度下,局部配点反馈迭代算法计算速度比拟线性化-局部变分迭代法提高1.5倍以上,引入变参数方案能够进一步将算法计算速度提高6倍以上.
其他文献
增材制造TC4钛合金是具有优良的力学性能和工艺性能的金属材料,在航空航天领域已得到重要应用.近年来,在塑性力学的研究中,探究应力状态对金属材料变形和失效行为的影响得到广泛关注,然而大部分的研究都是在准静态下完成的,对于中高应变率下金属材料变形失效的研究较少.本文从增材制造TC4钛合金的基本力学性能出发,考虑应力状态和应变率对其变形和失效行为的影响,采用应力三轴度η和罗德角参数(θ)表征应力状态,设计了相应的试样形式和实验方法.利用电子万能实验机、高速液压伺服实验机以及分离式Hopkinson杆,结合数字图
采用基于谱元法线性稳定性分析方法,研究了高径比对GaAs熔体(Pr=0.068)液桥热毛细对流失稳的影响,同时结合能量分析揭示了热毛细对流的失稳机制.研究结果表明:与典型低普朗特数(例如Pr=0.011)熔体静态失稳模式和典型高普朗特数(例如Pr>1)熔体振荡失稳模式不同,GaAs熔体热毛细对流失稳模式依赖于液桥高径比(As).随高径比的变化,GaAs熔体热毛细对流存在两种失稳模式.高径比As在0.4≤As≤1.18范围内,热毛细对流失稳是从二维轴对称定常对流转变为三维周期性振荡对流(振荡失稳);高径比在
由复合材料构成的板结构一直以来受到很大关注,其中功能梯度碳纳米管增强复合材料(functionally graded carbon nanotube-reinforced composite,FG-CNTRC)具有异常优越的力学性能,使得诸多学者展开了对功能梯度碳纳米管增强复合材料板结构力学行为的研究.本文以FG-CNTRC板为研究对象,将一种新型的区域型无网格方法——广义有限差分法应用于求解基于一阶剪切变形的FG-CNTRC板结构的静态线性弯曲和自振模态问题.广义有限差分法(generalized fi
基于Pope修正的有效黏度假设,张量基神经网络(tensor based neural network,TBNN)构建了从雷诺平均方程湍流模型(RANS)的平均应变率张量和平均旋转率张量到高精度数值解的雷诺应力各向异性张量的映射.将高精度数值解用于TBNN的训练,从而使TBNN根据RANS求解的湍动能、湍流耗散率和速度梯度预测其雷诺应力各向异性张量,并与对应的高精度数值模拟结果以及风洞实验结果对比以评估TBNN的预测能力.本工作将TBNN的预测能力从低速域拓展至高超声速工况,分别对低速槽道流、低速NACA
由于外部周期性的振动而在液滴表面产生的Faraday不稳定效应广泛存在于超声雾化、喷涂加工等应用中,对Faraday不稳定性进行分析对研究振动液滴的表面动力学有着重要意义.本文将Faraday不稳定性问题从径向振动拓展到竖直振动,研究了竖直振动无黏液滴表面波的不稳定性.竖直方向的振动使得液滴动量方程为含有空间相关项和时间周期系数的Mathieu方程.采用Floquet理论进行求解,得到了竖直振动液滴表面波线性增长率与模态数以及流动参数之间的色散关系.通过求解一个关于表面变形模态的特征值问题,得到了竖直振动
针对同轴气流式液体射流分裂液滴粒径预测模型缺乏的现状,结合射流线性稳定性理论,建立了基于临界模数的同轴气流式黏性液体射流分裂液滴粒径表达式,在此基础上,分别研究了气流旋拧(气流同时存在轴向和周向运动)及流体物性(气体可压缩性、液体黏性、气液密度比和表面张力)对液滴粒径的影响规律.研究发现:周围气流轴向引射作用和同轴旋转作用均会导致分裂液滴粒径整体呈先增大后减小的趋势;且在气流仅作同轴旋转运动时,相同临界模数下气流旋转对分裂液滴粒径的影响较小.在本文的研究参数范围内,分裂液滴粒径随气体可压缩性和气液密度比的
对称性是振动理论中5大美学特征之一,然而对称性破缺又难以避免.本文以工程中常见的易损结构—悬索为例,探究当该系统遭遇非对称性损伤时,对称性破缺对其面内耦合振动特性影响.首先建立受损悬索面内非线性动力学模型,并采用Galerkin法得到离散的无穷维微分方程.利用多尺度法计算该非线性系统发生面内耦合共振响应的调谐方程.截取前9阶模态,利用数值计算方法得到无损和受损悬索的各类共振曲线及其稳定性,通过计算最大李雅普诺夫指数来确定系统的混沌运动.研究结果表明:已有研究常采用抛物线模拟悬索静态构形,然而一旦发生不对称
由于单元应力属于局部性能约束,导致相应的结构拓扑优化存在难以承受的大量约束条件;尽管化整方法极大地减少了约束数量,但是优化结果中有少数应力超限现象.为此,本文在应力约束的结构拓扑优化中,瞄准克服应力超限和提高求解效率两个目标,进行了探索.提出了乘子法及序列二次规划(SQP)法两种解法,首先在化整交融(即化整-集成)解法中的m方集成模型应用,与一阶近似的移动渐近线(MMA)解法进行了求解效率对比.然后,在此基础上采用了效果最好的m方集成模型的SQP解法,建立了求解应力约束下结构体积极小化模型(即s方模型),
深埋软岩隧道围岩表现出显著的塑性软化与剪胀特性,而当下的理论分析很少同时考虑这两点,导致预测结果与隧道实际变形行为存在一定误差.为解决该问题,本文基于Kelvin-Voigt流变模型和Mohr-Coulomb强度准则,考虑了塑性阶段时围岩软化与剪胀特征,并引入了掌子面空间约束效应,建立了深埋软岩隧道黏弹-塑性计算分析模型.进一步,为考虑锚杆对隧道围岩的支护作用,在理论模型中,利用等效刚度法建立了加固围岩的力学模型.结合围岩塑性半径与锚杆长度相对关系,给出了6种工况下考虑锚杆加固作用的隧道黏弹塑性力学响应的
黏弹性材料在航空、机械、土木等领域具有广阔的应用前景,而具有1.5自由度的非线性Zener模型能更好地描述其特性.因此,研究多尺度法的推广和应用具有重要意义.在传统多尺度法的基础上,推广并利用多尺度法对非线性奇数阶微分方程进行研究,解决非线性奇数阶系统的动力学求解问题.以非线性Zener模型为例,首先通过推广的多尺度法对该模型进行近似解析求解,并通过数值方法对解析解进行数值验证,结果吻合良好,证明该推广方法求解过程的正确性.然后,从解析解中推导出稳态响应的幅频方程和相频方程,从幅频曲线中发现对于弱阻尼系统