低重力环境下三维液体非线性晃动的数值模拟

本文讨论低重力环境下三维液体非线性晃动问题。采用ＡＬＥ（任意的拉格朗日－欧拉）运动学描述跟踪自由液面；对于ＡＬＥ描述的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，在时间域上采用分步离散方法中的速度修正格式进行数值离散，在空间域上利用有限元方法进行了数值离散；推导了考虑表面引力效应时有限元边界条件的弱积分形式；给出了表面张力的数值计算公式。模拟了考虑表面张力情况下圆筒型贮腔中液体的非线性晃动，并得到了贮腔壁面处     

低重环境下圆柱贮箱内液体晃动特性研究的一种解析法

考虑低重环境下圆柱贮箱内由于表面张力影响而呈现弯曲自由液面的情况,利用Fourier Bessel级数对贮箱受横向激励时的自由液面处的边界条件进行展开,得到液体晃动系统的广义状态方程,给出了固有频率、晃动波高、晃动力和晃动力矩等晃动特性的计算公式。通过数值算例与文献结果对比,在验证文中方法正确性的同时,具体研究了各晃动特性随充液深度、外激励频率和Bond数等参数变化的规律。     

悬挂工装对运载火箭模态试验影响分析

为分析悬挂工装对运载火箭模态影响,通过系统能量泛函建立了悬挂梁系统的振动方程和边界条件,得到了系统固有频率特征方程,并证明自由-自由边界模型是其蜕化情况。通过数值求解特征方程,得到了系统频率和振型相对于悬挂工装的长度、密度变化的规律。采用此方法解释了某型火箭模态试验出现的多组一阶横向模态现象。     

基于HTML语言的多媒体E—MAIL系统

介绍了一个在ＩＮＴＥＲＮＥＴ上实现的基于ＨＴＭＬ语言的多媒体Ｅ－ＭＡＩＬ系统，提出了一个适合实现多媒体Ｅ－ＭＡＩＬ的体系结构，通过将ＨＴＭＬ语言和Ｅ－ＭＡＩＬ有机结合起来，能达到系统的标准化和开放性等要求。本文在ＳＧＩ ＩＮＤＹ上实现。     

板式贮箱内部流动的定常和非定常数值模拟

为了研究微重力环境下板式表面张力贮箱的内部流动规律,对板式贮箱导流板的驱动机理进行了分析,得到了微重力环境下板式贮箱内部定常和非定常流动的基本方程。应用四阶龙格-库塔法对定常流动方程进行求解,得到定常流动时导流板上液带的分布规律;应用MacCormack数值离散方法,对非定常流动方程进行求解,得到了不同时刻导流板上液带的分布规律及液带恢复稳定所需的时间。计算结果为板式贮箱的设计提供了数值依据,奠定了板式贮箱的设计基础。     

卫星液体晃动气浮台仿真试验系统的运动分析

本文通过边界元数值法求解了部分充液自旋球腔内的液体晃动问题。以流体运动的基本方程和系统运动的Ｅｕｌｅｒ动力学方程为基础，考虑了贮箱偏置、涡旋、重力及Ｃｏｒｉｏｌｉｓ力等因素对流体晃动和系统运动状态的影响，求解出液体的速度场，并在此基础上估算液体的能量耗散率和系统的章动时间常数     

边界元法计算厚板表面裂纹问题的应力强度因子

建立了一种能反映裂纹尖端应力、应变奇异性的三维等参奇性边界元，应用边界元分析方法求得了具有半椭圆形和半圆形表面裂纹厚板的应力强度因子。结果表明，采用这种计算方法得到了满意的计算效果。     

全失重液体晃动的固有频率

本文研究失重条件下液滴，带核液体和球腔内液体晃动的固有频率，对于几何球对称系统，本文给出了其解析解；而对于非球对系系统，则采用边界元方法计算得到晃动频率的数值解。     

超音速射流与弹体绕流干扰流场数值分析

绕弹体的超音速气流与火箭发动机燃气喷流相互作用，在弹体底部形成了复杂的干扰流场。本文应用了四阶的ＭＵＳＣＬ ＴＶＤ格式，求解雷诺平均的Ｎ－Ｓ方程，对这一复杂流动进行了数值模拟，得到了流场结构随不同喷口压力的变化规律。     

典型结构的高频响应研究

利用ＮＡＳＴＲＡＮ软件得到一简支梁从低频到高频的频响特性，研究了不同有限元网格划分对结构高频特性的影响，并通过ＳＥＡ法、有法元法以及解析解的对比揭示了ＳＥＡ法在高频响应预示中的有效性。然后以耦合梁为例，通过有限元法得到两个梁之间的能量传递，并进而得到ＳＥＡ所需的耦合损耗因子；将有限元计算结果和ＳＥＡ结果对比验证了该方法的有效性，以其建立一种用有限元法获得ＳＥＡ参数的新途径。此外，研究了集中质量和不     

微重力下球形贮箱内液体晃动有限元模拟

以卫星球形贮箱为研究对象,用任意拉格朗日-欧拉(ALE)法建立了数学模型,给出了计算的基本方程,并基于MSC-Dytran有限元程序,对常重力和微重力环境中贮箱内推进剂的晃动进行了数值模拟。结果表明,微重力环境中液体更易发生大幅晃动,甚至产生飞溅,工程中常用的等效力学模型已不再适用,而有限元数值模拟方法可较好地解决该问题,为卫星及大充液复杂结构航天器控制系统的设计和分析提供理论参考。     

半直接配点法在航天器追逃问题求解中的应用

采用半直接配点法求解时间固定两航天器追逃问题,提出一种新的数值求解追逃双方最优控制策略的方式,避免了求解非线性两点边值问题。在两航天器均为连续小推力假设条件下,以终端距离为支付函数,给出了半直接配点法求解此追逃问题的过程。在此数值方法中,根据半直接转换将微分对策问题转化为一个最优控制问题,由Gauss-Lobbato配点法最终将此最优问题转化为非线性规划问题,继而通过序列二次规划方法求解。这种半直接配点法避免微分对策问题最优策略的必要条件(两点边值问题)求解,并且数值稳定性好。数值仿真给出了追逃双发的最优控制策略和相应的追逃轨迹。     

正交各向异性复合材料中裂纹问题的边界元法分析

在基本解的基础上，研究了正交各异性复合材料平面问题边界元法的有关基本问题，包括Ｃ矩阵，Ｇｉｊ矩阵，Ｈｉｊ矩阵和域内应力表达式等，并建立了常值边界单元和线性边界单元的计算公式（适用于各类边值问题）。用建立的公式对正交各向异性板含中心裂纹问题和含圆孔的无限大正交各向异性板问题进行了分析与计算，取得了较好的计算效果。     

包带连接力学分析与建模

研究了卫星包带连接的有限元力学建模。基于对包带预紧力、卡块与连接环间的接触力,以及两连接环间的接触力的分析,将接触力等效为有限单元节点力。对所建模型中的非线性项,用Newton-Raphson法求得给定包带顸紧力的平衡位置,再以该平衡位置处的切线刚度对模型作线性化,分析了模态和固有频率等动力学特性。某卫星的算例验证了模...     

基于有限元法的大气层外弹道式目标表面温度场分析

获取目标表面温度场是进行红外特征分析的重要前提。为确定大气层外弹道式目标的表面温度场分布,建立了有限元模型。根据目标温度场的轴对称分布特点,在柱坐标系内建立了二维瞬态热传递模型,从而降低了计算量。与当前大部分的有限元软件和文献不同,严格分析了温度非均匀性有限单元的辐射热损失,并推导出轴对称三结点三角形单元的辐射热损失公式。为了能够使用Galerkin法求解时间微分方程组,将辐射热损失视为与时间相关的热载荷项,并用Newton\|Raphson法迭代求解处理后得到的非线性方程组。一个简单的数值实验表明了所提出方法的有效性,最后应用该法求解了目标飞行全程的表面温度场分布。
     

航天器再入全过程轴对称烧蚀热防护数值仿真研究

对航天器再入全过程轴对称烧蚀热防护进行了全过程数值仿真研究.采用修正Lees驻点热流密度方法和参考焓方法计算再入热流密度.采用JANAF模型计算烧蚀率.利用有限元法计算钝锥体再入航天器烧蚀层在移动边界条件下的轴对称温度场.采用碳化层-热解面-原始材料的轴对称碳化烧蚀模型;推导了热解气体流量计算方法.针对再入飞行大热流密度条件下,用有限元方法求解瞬态温度场时会产生的时间和空间上解的振荡问题.通过分析温度振荡现象产生的原因,采用集中热容矩阵向后差分方法解决振荡问题.计算结果表明,在时间步长选择合适的情况下,求解集中热容矩阵能够很好地解决数值振荡问题,同时烧蚀率和温度场计算比较准确.     

纤维缠绕复合材料壳体的非线性大变形有限元分析

叙述了纤维缠绕复合材料固体火箭发动机旋转壳体与回转体构成的组合结构在轴对称大变形状态下的非线性有限元分析方法。给出了在轴对称大变形条件下的八结点等参单元和拟协调双曲壳单元的有限元列式。算例结果与已知的的理论符合较好。     

基于虚温定义的粘弹性有限元方法

基于NARAYANASWAMY提出的虚温定义,用虚温表示的温度变化产生的折算应变率替代Swanson模型中的应变率修正函数,获得了基于虚温定义的粘弹性材料的应力应变本构方程,并导出基于虚温定义的粘弹性非线性有限单元法计算公式。计算结果验证了此基于虚温定义的粘弹性分析方法的正确性。     

Adams--Cowel方法的扩充及其在人卫中的应用

本文将现有的Ａｄａｍｓ－Ｃｏｗｅｌ方法由求初值问题扩展到求一类边值问题，并将之应用于对同步卫星的轨道测定和预报，仿真与实算结果表明：方法具有独特的使用价值