非轴对称贮箱液体晃动的等效力学模型

针对非轴对称贮箱液体晃动问题,从一般意义上推导了作三维运动的任意三维贮箱内液体晃动的等效力学模型,采用有限元方法建立计算等效力学模型参数的数值算法,利用矩阵相似变换的性质分离液体零频以消除液体刚度矩阵的奇异性,从而将等效力学模型参数的计算归结为广义特征值问题和线性代数方程组的求解,具有较高的计算效率。通过非轴对称贮箱液体晃动的算例说明,不同的等效晃动方向可能不正交,晃动平面可能不过液体质心,因此非轴对称贮箱液体晃动的等效力学模型比轴对称贮箱液体晃动的更加复杂。最后以三棱柱形充液贮箱为例,采用商用CFD（计算流体动力学）分析软件FLOW-3D验证了在贮箱作典型三维运动时等效力学模型的有效性。
     

卫星载荷磁场优化计算方法研究

针对传统计算方法存在的计算复杂、计算成本高等不足,利用某处地磁场下卫星载荷的磁场测量数据,逆计算出卫星载荷的等效磁化率,再根据等效磁化率就可计算在不同磁化方向、不同地磁场大小下卫星载荷产生的感应磁场。根据均匀磁化卫星载荷的特点,将积分方程法中系数矩阵元素的体积分形式化为面积分形式,并消除积分计算的奇异性,建立基于单元表面积分的感应磁场计算模型,从而降低计算成本。最后,实验结果表明:磁场测量值和计算值吻合得较好,该算法的有效性和可行性得以验证。     

介质体散射分析中奇异积分处理

为消除高阶矩量法（MOM）求解表面积分方程时的奇异积分,提出了一种采用混合域基函数精确计算奇异积分的方法。通过分离奇异积分项,用变量替换求解析解,与高斯积分法获得的非奇异积分值相加,求得阻抗元素。仿真计算表明：该法的结果和文献中低价MOM吻合。     

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

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

电动伺服机构振荡问题研究

针对电动伺服机构中传动间隙引起的极限环振荡和较大惯量的负载加剧这一振荡问题,建立了系统间隙极限环模型和非线性振动模型,并采用摄动法求解出非线性振动规律,得出了间隙、转动惯量、伺服刚度及干扰力对系统的影响规律。仿真和试验结果证明这一理论是正确的。     

复合材料的二维轴对称瞬态导热分析

复合材料的瞬态导热问题在卫星、航天器和微电子等工程领域普遍存在,在热控制、热设计中起着非常重要的作用。本文从一般理论问题出发,采用汉克尔积分变换、分离变量法及矩阵运算方法,对复合材料的二维轴对称瞬态导热问题进行理论求解,并导出它们的通解。     

一种应用小波变换提高SAR距离向分辨率的设想

在分析信号和噪声的Lipschitz奇异性的基础上,提出了采用小波变换检测信号奇异性来估计SAR中线性调频脉冲回波的到达时间的新方法。由于小波变换能很好地检测奇异性,因此,为有效提高SAR距离向分辨率提供了一个新的途径。最后,给出了计算机模拟的结果,表明该方法是十分有效的。     

三维横向流喷与超声速主流干扰流场的数值模拟

本文对三维隅角机翼声速喷流与超声速主流的干扰流扬进行了数值模拟。三维欧拉方程的求解采用非结构网络有限体积伽辽金法（Finite Volume Galerkin Method）。引入了总体结点积分域的概念,简化了从单元矩阵到总体矩阵的汇总过程。通量的分裂采用Osher格式,通过外差使其由一阶精度上升为二阶精度。发展了一种基于线化流量的逆风非结构网格隐式有限元格式以提高求解精度及效率。最后给出了三维隅角机翼流场的算例。     

扭簧驱动构架式空间展开天线结构分析

首先仔细分析了扭簧驱动构架式空间展开天线结构的工作机理与构造特点，据此提出分析建模假设。采用静力缩减方法，由平面梁单元弯曲平衡方程推导出两端扭簧的杆件宏单元缩减刚度矩阵。根据梁的弯曲理论和位移形函数，直接导出宏单元的相容一致质量矩阵，这两个矩阵可用于这类结构在各工况下的结构特性分析。最后给出了一个数值分析算例，证明方法的有效性。     

铰接桁架结构动力学问题研究

以铰接桁架结构为对象,采用谱元法研究铰对桁架结构动力学行为的影响,拓展了谱元法的应用领域。在频域下将铰结构考虑为一个谱单元,分别建立铰、杆和梁单元的动力学刚度矩阵,并加以整合得到整体结构的动力学刚度阵,进而建立整体结构的动力学方程。通过求解整体结构的动力学方程,获得结构的固有频率、频响曲线和时间响应历程曲线,并将谱元法的计算结果同有限元法的相应结果进行了比较。分析结果表明,谱元法在求解铰接桁架结构的动力学问题上具有较高的精确性,并且铰的存在对桁架结构动力学行为具有显著影响。     

航天器着陆碰撞动力学CUDA并行计算技术

航天器着陆接触碰撞过程中的瞬态响应计算,涉及到几何非线性、接触非线性和材料非线性,计算过程复杂,计算量大。为此,文章将统一计算架构(CUDA)技术应用到接触碰撞问题的并行显式积分计算,研究了质量矩阵和刚度矩阵的并行生成方法,并将基于中心差分法的求解过程并行化实施。仿真算例表明:并行算法计算结果准确,较大地提高了计算仿真分析的效率,可为航天器着陆过程中冲击动力学计算的快速和高精度实现奠定基础。     

三维正交异性材料旋转体热结构分析方法的改进

在轴对称温度场和轴对称分布载荷作用下,由三维正交异性材料制造的旋转体的应力、应变计算是一种非轴对称问题。采用二维轴对称力学模型求解是比较粗糙的,不能令人满意。采用三维有限元法求解则计算工作量很大,划分单元体等方面也比较复杂,因此也不够理想。本文利用物体形状的轴对称性,采用半解析有限元法来求解这个非轴对称问题。根据三维正交异性材料的特点,本文分析了方程组的偶合特征,提出了大幅度降低节点自由度的方法,从而简化了计算工作。为三维正交异性材料旋转体的热结构分析提供了一个较为满意的方法。     

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

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

再入飞行器多变量稳定裕度研究

回差矩阵奇异值法是一种计算多变量系统稳定裕度的方法,由于其计算简单,易于计算机实现,已得到了较广泛的应用。但是按照这种方法计算的系统稳定裕度存在一定保守性,本文提出了基于回差矩阵奇异值的摄动矩阵迭代法,改善了计算多变量系统稳定裕度的保守性,并通过定点仿真验证,取得了较好的效果。     

固体火箭发动机界面脱粘裂纹分析

使用有限元法,在裂纹尖端周围布置有限奇异裂纹单元以模拟裂纹尖端附近的奇异性。针对轴对称发动机头部的界面脱粘裂纹,计算了点火内压作用下,发动机衬层/药柱、壳体/绝热层界面不同深度脱粘裂纹尖端的应力强度因子,指出应力强度因子随裂纹深度的发展规律。结果表明,当裂纹深度较小时,衬层/药柱界面处于闭合状态,应力强度因子几乎不发生变化,随着裂纹深度的增加,裂纹呈张开状态,裂纹尖端的应力强度因子不断增大;壳体/绝热层界面裂纹总是处于张开状态,且应力强度因子随裂纹深度的增加而增大。     

相变及炭黑杂质对氧化铝粒子流辐射特性的影响

基于Bruggeman有效介质理论与反向蒙特卡罗法计算了轴对称高温氧化铝粒子流在轴线正侧向辐射强度分布。采用Bruggeman有效介质理论计算多相态共存氧化铝粒子以及含少量炭黑杂质的氧化铝粒子的等效光学常数,以此为基础,通过反向蒙特卡罗法考察氧化铝粒子相变过程以及粒子中少量炭黑杂质对粒子流场轴线上辐射特性的影响。发现粒子相变过程中多相共存状态对高温粒子流场紫外波段辐射信号的影响不可忽略,氧化铝粒子中含有少量炭黑杂质就能够明显降低粒子流高温区间轴线上光谱辐射强度,当波长增加到近红外波段时,离喷口较近处轴线正侧向相对光谱辐射强度随着含炭黑量的增加而增大。     

基于基元反应的总包机理建模及算法优化

以煤油/氧的化学动力学体系为例,介绍了基于基元反应的总包机理建模、反应速率计算等国际先进方法,归纳了针对化学动力学过程由建模到求解的一套完整流程。首先选用正十二烷作为煤油的替代燃料,通过直接关系图法(DRG)结合计算奇异值摄动法(CSP)的机理简化方法将正十二烷203组分、738步的详细反应机理化简为了32组分、36步的总包反应机理。为了计算总包反应速率,采用线性准稳态假设(LQSSA)计算准稳态组分浓度,解耦了准稳态组分间的非线性耦合,在保证计算精度的同时提高算法稳定性和效率。最终,通过VODE刚性积分求解器求解化学动力学方程组,并形成了一套可与CFD主程序耦合使用的化学动力学计算子程序。     

考虑前屈曲几何非线性影响的结构稳定性分析

本文采用有限元法分析考虑前屈曲几何非线性影响的结构稳定性。应用伴有线性搜索和各种牛顿型解法的组合弧长法可以顺利地描述包含有极值点的载荷一位移曲线,与此同时,若在加载点处计算结构刚度矩阵的行列式值、临界载荷因子和特性刚度参数等值,则可根据稳定准则,通过插值确定临界点的位置。从而为研究复杂结构稳定性提供简便而有效的方法。     

航天器发动机羽流对敏感器热效应仿真研究

采用差分求解N-S方程与DSMC方法相结合的方法,研究了航天器单台发动机连续工作情况下真空羽流对航天器敏感器的热效应。首先通过求解N-S方程,获得发动机喷管的内流场,然后应用DSMC方法对喷管出口外轴对称羽流场进行计算,最后将轴对称羽流场计算结果作为模拟粒子入口边界条件,在并行计算机平台上进行三维羽流场和热效应计算得到航天器单发动机连续工作情况下羽流场对敏感器的热效应。以两个敏感器为例,对仿真结果进行了分析和比较,并得出了相应结论。     

索网结构非线性有限元分析中初始张力状态的快速构造方法

预应力索网结构是以柔性拉索为主要承重构件或辅助承重构件的一种结构体系,用于可展开天线等空间结构。为在有限元数值分析时准确、高效地施加张拉预应力,综合考虑索网的几何非线性、刚度矩阵奇异性、索网与支撑结构存在耦合变形等问题,利用索网式可展开天线的特性,在非线性有限元计算时忽略初位移刚度矩阵,将问题线性化,并详细介绍了以降温法模拟预应力时降温值的计算过程。此方法无需迭代,且得到的降温值在用于索网非线性计算时具有足够的精度;此方法能充分利用通用有限元计算程序,可进一步开展模态和响应等力学计算,便于在工程中实际应用。