Mesh-dependence of Material with Softening Behavior

The issue of mesh-dependence emerges when the conventional continuum damage model is applied to handling the softening behavior. In order to circumvent the mesh-dependence, the non-local theory is introduced into the conventional damage model and the finite element formulas are derived for two-dimensional gradient-enhanced damage model. A new element is proposed in which the basic unknown quantities are displacement, non-local equivalent strain and the gradient of non-local equivalent strain. The element and constitutive equation proposed in this article are added to the finite element software ABAQUS through user subroutine UEL. Numerical results show that the gradient-enhanced damage model can eliminate the mesh-dependence and is effective for dealing with the issue of softening behavior.     

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.     

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

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

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

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

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

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

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

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

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

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

遥感卫星扰振对TDICCD相机图像影响的仿真与分析

随着传输型遥感卫星结构的扩大化及灵活性的增强,对图像品质产生影响的扰振越来越受到重视。当遥感卫星在轨进行成像时,相机的图像品质会由于星上活动部件的运动的影响而下降。文章从调制传递函数（Modulation Transfen Function,MTF）角度出发,对不同工况下太阳翼转动驱动机构（BAPTA）和动量轮对图像品质的影响进行了仿真和评估。同时,针对高分辨率相机成像品质影响,从动力学仿真方面对星体扰振的幅值和频率进行了有限元分析,确定了扰振源的频率和幅值。最后,根据仿真分析结果,提出了星上扰振控制方案建议。     

定向破坏复合材料层合壳体结构铺层的分析与设计

以复合材料层合板为基本材料设计了一种结构新颖的定向破坏复合材料层合壳体结构,在导弹与

该壳体结构的初始接触区域预制缺陷区,导弹发射时以四瓣形式开裂,大幅度提高了导弹的发射效率。利用ANSYS

有限元模拟分析软件,通过强度分析确定了最佳铺层设计方案,单层板厚度为0.2 mm,采用基体本身加强预

置缺陷区的[0°/90°/ ±45°]2S 铺层方式。该种铺层方式的爆破载荷为0.21 MPa,最小顶破载荷为1 kN,完全能够

满足其在工作环境下既能承受一定的压力不漏气又能在较小的顶破载荷下自行破裂的设计要求。

     

基于LBM-LES方法的多段翼型流场数值模拟

运用Lattice Boltzmann(LB)方法与Large Eddy Simulation(LES)结合的方法数值模拟研究了二维翼型绕流流场,其中,粒子速度模型采用D2Q9模型,松弛时间中对应的湍流粘性系数由Smagorinsky模型计算,网格生成方法采用无网格技术。对NACA0012翼型进行数值模拟,验证其有效性;为进一步研究其模拟复杂物体绕流的能力,以三段翼型30P30N为研究对象进行数值研究,计算出了给定攻角条件下的气动特性,并与实验值进行比较。计算结果表明,LBM-LES方法的计算结果能够与实验值一致,是一种可行、有效且较新的数值模拟方法。