EBE技术在结构分析中的应用（Ⅰ）─EBE矩阵向量乘积方法及其实现

矩阵向量乘积运算是数值分析领域中的一种基本运算。本文基于单元分组技术，提出了在有限元结构分析中，当不形成总刚度矩阵而进行矩阵向量乘积的ＥＢＥ计算方法。该方法不仅预处理过程简单，实现容易，而且其并行化、向量化处理都极为直接。同时，数值试验结果还表明了它的有效性。     

The use of Gaussian equations of motions of a satellite for local gravity anomaly recovery

The orbital elements of a low Earth orbiting satellite and their velocities can be used for local determination of gravity anomaly. The important issue is to find direct relations among the anomalies and these parameters. Here, a primary theoretical study is presented for this purpose. The Gaussian equations of motion of a satellite are used to develop integral formulas for recovering the gravity anomalies. The behaviour of kernels of the integrals are investigated for a two-month simulated orbit similar to that of the Gravity field and steady-state ocean circulation explorer (GOCE) mission over Fennoscandia. Numerical investigations show that the integral formulas have neither isotropic nor well-behaved kernels. In such a case, gravity anomaly recovery is not successful due to large spatial truncation error of the integral formulas. Reformulation of the problem by combining the orbital elements and their velocities leads to an integral with a well-behaved kernel which is suitable for our purpose. Also based on these combinations some general relations among the orbital elements and their velocities are obtained which can be used for validation of orbital parameters and their velocities.     

可压缩流动激波装配在格心型有限体积法中的应用

发展了一种基于格心型有限体积方法(FVM)的激波装配算法。通过定义网格节点属性可以灵活调用激波装配和激波捕捉计算方法。在使用激波装配方法时,激波节点运动速度和下游运动速度通过Rankine-Hugoniot(R-H)关系式获得,同时采用非结构动网格技术描述激波的运动以及调整其他网格节点的位置。流过激波面元的通量为上游单元的基本通量,物理概念更加清晰,通量计算也更为准确。在计算过程中,网格节点属性可以发生变化,以此实现对带有拓扑变化流场的描述。数值试验表明:本文提出的计算方法不但具有较高的计算精度,同时能有效地避免由于捕捉激波而出现的数值问题。     

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

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

Bush有限单元原理及其在航天器结构建模中的应用

利用大量的数值计算确定了MSC.NASTRAN软件中广义弹簧单元（Bush）刚度阵的解析表达式,通过对比Bush元与剪切梁元的刚度阵,推导了Bush元参数与剪切梁刚度参数的对应关系,并利用展开状态的太阳翼结构和大型桁架结构的静力和动力分析验证了推导结果。该研究揭示了Bush有限单元的原理和所有输入参数的物理含义,有助于Bush元在结构有限元分析中的准确使用。     

热-机耦合分析的有限元法及其应用

基于有限元分析的基本原理,建立了热-机耦合问题分析的有限元基本方程,并对工件挤压过程及摩擦制动问题进行了分析。算例表明,热-机耦合分析方法能很好地模拟受热结构的温度场和应力应变场,为结构设计和工艺设计提供依据。     

运输环境中某战术导弹结构振动特性分析

应用三维有限单元和导弹实体简化模型对一种战术导弹的结构振动特性进行分析，并动用ＡＮＳＹＳ有限元分析程序，计算得到了“自由－固支”条件下，导弹竖立状态的固有模态，取得了与试验和数据基本一致的结果，证明了所建模型物正确性，进而对运输环境中导弹水平支承条件下的固有频率和振型进行了分析计算，所得结果为导弹的使用可靠性分析提供了依据。     

三维角联锁结构复合材料等效弹性性能有限元分析

在复合材料承载时变形协调的前提下,采用CAD软件绘制三维织物单元体几何模型,然后导入结构分析有限元软件包,对工程中常用的带衬经三维角联锁结构复合材料的弹性性能进行分析,并将0°和90°方向弹性性能预测值及实验值与取向平均法预测值进行比较。结果表明,采用的有限元方法所得预测值及实验值与取向平均法预测值之间具有很好的一致性。     

针栓式喷注单元膜束撞击雾化混合过程数值模拟

为了全面认识针栓式喷注器喷雾场结构，基于自适应网格加密技术和分三相计算的PLIC VOF（Piecewise Linear Interface Calculation Volume of Fluid）方法对针栓式喷注单元膜束撞击雾化混合过程进行了仿真分析，通过对两路推进剂分别进行界面追踪，获得了膜束撞击雾化混合过程的详细结构特征，与高速摄影试验结果定性定量对比均吻合较好，验证了数值方法的准确性。以此为基础对膜束撞击的喷雾场结构、撞击变形过程、流场涡结构、雾化破碎典型特征及破碎后的雾化混合分布特征进行了识别分析，结果表明：膜束撞击形成了液束未穿透液膜和液束穿透液膜2种不同的喷雾扇结构。膜束撞击形成的喷雾扇呈"Ω"形，膜束同时发生弯曲变形和横截面变形。另外，膜束撞击同时受到正压和剪切应力作用，导致了一系列复杂涡流现象，使得相互作用增强，雾化混合均增强，这也是膜束撞击喷注构型优于膜膜撞击的本质原因。最后，还发现膜束撞击喷雾场液滴分布呈现分区结构特征，分别是液束控制主导的上雾化区、液膜控制主导的下雾化区及夹在中间的混合区，实际中应兼顾雾化特性和混合特性，选取中等动量比膜束撞击，这可为针栓式喷注器的理论研究和工程设计提供重要参考。     

尾涡数值模拟的曲涡元技术

连接涡线端点形成涡段并计算其对空间任意点的影响系数是尾迹分析的基本环节,称为涡元技术。本文推出圆弧涡元与其应用程序,并引进抛物拱弧涡元,以代替传统的直涡元。提供的曲涡元程序可方便地用于尾迹分析,算例表明曲涡元在提高精度减少计算量方面的先进性。