Dynamic Response Analysis of Bird Strike on Aircraft Windshield Based on Damage-modified Nonlinear Viscoelastic Constitutive Relation

Damage-modified nonlinear viscoelastic constitutive equation and failure criterion are introduced and the three-dimensional incremental forms are deduced based on the updated Lagrangian approach. A simple tensile test model and a split Hopkinson pressure bar model are built to verify the accuracy of the subroutine implemented within the non-linear finite element program LS-DYNA. A numerical model of bird strike on windshield is established to study the responses of windshield under three different bird velocities at three sites. The bird is represented by a cylinder with a hemisphere at each end and the contact-impact coupling algorithm is used in this study. It is found that the implemented subroutine can properly describe the mechanical behavior of polymethyl methacrylate under low and high strain rates and large deformation, and can be used validly.     

非结构网格上Euler方程的区域分裂算法及并行计算

改进了“波阵面”区域分裂算法,并应用于流场区域的划分;对于子区域边界的不光滑现象,为尽量减少通讯消耗,提出了一种边界并行优化策略。利用PVM并行环境,探讨了非结构网格上求解Euler方程的分区并行算法。根据改进的区域分裂算法及优化策略,运用Jameson有限体积法,对二维翼型流场进行了分区并行求解,多区计算的结果与单区计算的结果作了比较,表明了本文研究方法的有效性。     

隧道内危险性介质泄漏扩散的数值模拟

化学危险品的意外泄漏已经造成很多事故,尤其是在隧道这样相对封闭的空间内,事故发生的概率更大。为了有效地控制和预防事故的发生,本文以RANS方程为控制方程,采用基于simple算法的有限体积法,对隧道内危险性介质泄漏扩散过程进行了数值模拟。讨论隧道内危险品(乙炔)运载车在行驶过程中发生小孔泄漏扩散的情况,计算得到了危险性介质浓度的时间和空间分布规律。同时,分析了运载车不同行驶速度和危险品泄漏释放方式(水平释放和垂直释放)对其扩散过程的影响,为进一步了解危险性介质的扩散行为和控制方式提供了理论基础。     

耦合Navier-Stokes方程的可压缩边界层数值求解

为了在Reynolds-averaged Navier-Stokes(RANS)方程计算中耦合eN方法进行转捩判断,在RANS方程求解过程中耦合求解了可压缩层流边界层方程为判断转捩提供了精确、可靠的边界层信息.利用等熵关系由RANS方程求出的物面压力分布确定边界层外边界的速度分布,进一步确定出边界层外边界.边界层方程的求解使用Keller提出的二阶BOX方法.为了验证方法求解边界层方程的正确性,在低速流动状态下将计算结果和XFOIL的计算得到的边界层解进行了比较,二者吻合较好.     

基于RBF插值边界条件的笛卡尔网格方法

介绍了一种基于RBF(Radial Basis Function,径向基函数)插值处理物面边界条件的笛卡尔网格方法.整个流场采用有限体积方法计算,但是对于物面附近切割的网格计算域不完整,因此通过RBF插值来得到物体内部实单元上的值,从而使得流场计算顺利的进行推进.相比于其他插值方法来说,RBF插值特别适合于多维离散点的插值,不会因为插值点的位置关系而产生病态的系数矩阵,而且使用比较少的网格点就可以得到同最小二乘方法相当的精度.将这种方法应用于圆柱和翼型绕流中,得到了满意的结果,证明了方法的有效性,还对边界条件的各种不同形式对结果的影响作了对比.     

NEW APPROACH TO OBTAIN SPACE-TIME CONSERVATION SCHEMES

Inordertoovercometheseverallimitationsofthetraditionalmetheds,anewnumericalmethodofspace-timeconservationelementandsolutionelement(CE/SE)forsolvingconser-vationlawswasdevelopedbyChang['j.Thismethoddifferssubstantiallyinbothconceptandmethedologyfromthewell-establishedmethodslsuchasfinitedifference,finiteelement,etc.,andithasmanynontraditionalfeatures.First,spaceandtimeareunifiedandtreatedonthesamefooting,andbytheintroductionofconservationelementandsolutionelement,bothlocalandglobalfluxconserv…     

一类电磁脉冲波方程的初值问题

对是民磁脉冲波方程的数学征物理模型之初值问题进行了研究，应用细致的先验估计，建立了该问题解的整体存在性与唯一性，并对解趋于│ｘ│→∞处的渐近性态做了讨论。     

自适应网格局部加密方法在可压缩流场模拟中的应用

本文基于M.J.Berger和J.Oliger发展的自适应网格局部加密算法,它采用多个组合网格的思想,用于求解双曲型方程组。它结合显式有限差分格式,采用Richardson外推技术自动进行局部截断误差估计,对精度低的区域产生新细网格局部加密,或去除无需再加密的旧细网格,以最小的运算量达到给定的精度要求。网格可一层层加密下去,按层次覆盖,每个网格是有任意方向的矩形均匀网格。这套算法独立于求解所用的差分格式,很容易和各种格式结合。 我们从M.J.Berger和J.Oliger程序出发,参考M.J.Berger和P.Colella,对守恒型方程组,实现了粗细网格交接面上数值通量守恒,从而能计算间断解。对网格边界处理、误差估计和网格重新生成等做了一些改进,并且可以灵活地在任何感兴趣区域指定加密。 我们用不定常Euler方程组计算了带前台阶的二维管道中Mach-3流动,采用MacCormack显式附加人工粘性差分格式。共用了四层网格结构:Δ=1/10、1/20、1/80和1/320。其中基本网格是Δ=1/10;到台阶上第一次激波反射区域从Δ=1/20加密至Δ=1/80,计算到定常所需运算量仅是同等均匀网格△=1/80的28％;我们还用Δ=1/320网格对台阶拐角和切向间断处指定加密,结果表明可明显改善台阶上的激波Mach反射,切向间断清晰可辨,所需的计算量是上述均匀网格Δ=1/     

卫星太阳阵展开的动力学方程

研究了卫星太阳阵展开的动力学问题，将卫星太阳阵的结构视为树形多体系统，利用给出的自由形式树形多刚体系统运动微分方程建立太阳阵系统的动力学方程，应用为该系统动力学方程编制的计算程序，可对太阳阵展开过程中卫星本体姿态的改变量以及太阳帆板在展开过程中的运动参数进行数值计算。     

强粘性层流理论及在粘性流计算中的应用

对高雷诺数流动计算，为了解决部分粘性项甚至全部惯性粘性项落入误差与流动物理尺度的关系问题，本文提出强粘性流动理论。强粘性流中至少有一个粘性项与惯性项同量阶，理论包含了物理尺度各向极限、经典边界层和多层边界层理论为其特例，给出了从经典边界层向物理尺度各向相同极限演化的尺度规律和粘性惯性诸项变化的量阶关系，阐明了粘性与惯性力强相互作用将在剪切层的法向以流向同时“激发”涉尺度结构。对粘性流计算，利用强粘性剪切流尺度律重新标度NS格式的修正微分方程，给出临界网格尺度与流动物理尺度和差分格式精度的关系，得到部分粘性项落入误差和计算结果为非物理数值粘性解的二个判据。并以流场中的边界层、驻点和分离点领域计算为例说明理论的应用，对强粘性剪切流计算、证实部分粘性项甚至全部粘性惯性项落入误差的问题值得重视。