砌体结构长脉宽爆炸荷载损伤等级评估方法

随着百千吨级当量爆炸工业事故的频繁发生,建筑结构的损伤评估和抗爆安全性更受到关注。目前,构件级的评估方法相对成熟,而大当量冲击波作用下的建筑结构整体毁伤评估依旧是个开放性问题。本文中,面向结构级的毁伤评估,提出了一种新的评估方法?构件损伤加权。该方法以承重构件损伤程度为基础,通过基于应变能的构件权系数加权,进而评估结构级损伤破坏程度。为了验证评估方法的有效性,以典型砌体结构为例,利用自主研发的冲击波结构毁伤模拟有限元程序,开展了百毫秒脉宽爆炸冲击波荷载下结构动力学响应数值模拟。根据数值模拟结果,结合构件损伤加权的评估方法,获取砌体结构损伤等级与冲击波超压的关系。预测的超压值的相对误差为?16.9%～26.2%,验证了评估方法的有效性。该评估方法为获取砌体结构的超压-冲量曲线提供了可行的途径,可为结构的抗爆安全设计提供参考。     

改进型XFEM进展

扩展有限元法(XFEM)在诞生后的十几年时间里,引起学术界和工业界的广泛关注,并已经成为目前裂纹分析的主流数值方法。然而,在实际应用中该法一直受到两方面的困扰,(1)总体方程高度病态;刚度阵条件数随网格尺寸呈h-6变化(普通有限元为h-2)。(2)裂尖强化插值由于能量一致性问题无法直接推广应用于动力学计算。前者表现在XFEM稳态问题的迭代求解收敛慢或难以收敛,后者长期以来导致XFEM裂纹扩展动力学计算实施困难。本文认为XFEM目前遇到的种种困难,均与单位分解引入的额外自由度相关。为此,提出了无额外自由度的单位分解插值格式,基于此格式,进一步构造出改进型扩展有限元方法。改进型XFEM具有如下特点,(1)可以消除原有XFEM的线性依赖性和总体方程病态的问题。(2)避免动力学问题中额外自由度引起的质量集中、零临界时间步长问题以及裂纹扩展过程中的能量一致性问题。本文结合静动力学测试问题综述上述改进。     

基于VENUS-III国际基准模型的JMCT程序验证

基于自主研制的三维中子-光子耦合输运蒙特卡罗通用程序JMCT（J Monte Carlo Transport Code）,采用连续点截面,对国际基准屏蔽VENUS-III模型开展精细建模和中子输运临界及屏蔽计算.临界计算得到系统k_eff、重要区域的通量及能谱.结果表明,JMCT和MCNP程序的重要区域体通量计数吻合较好,偏差均在1%以内.深穿透屏蔽计算采用外源模式,点探测器计数,JMCT计算值与实验测量值偏差在15%以内,满足屏蔽设计对误差的要求.初步验证了JMCT程序临界及屏蔽计算的可用性.     

嵌入式薄片模型在时域有限差分算法中的应用

将一等效薄片模型嵌入到时域有限差分算法（FDTD）中,以快速而有效地解决复合材料薄片在电磁计算中的多尺度问题。在该嵌入式薄片模型中,薄片不需要被剖分网格,而是被嵌入到相邻的网格间,从而可以使用相对较大的网格剖分周围物体,进而可节省大量的计算资源。在这一模型中,薄片被等效为一段传输线,并用其频域的导纳矩阵代替。使用数字滤波器理论以及逆Z变换可将频域的导纳矩阵转换到时域,并将其嵌入到时域有限差分算法中。该嵌入式薄片模型被用来计算一单层碳纤维复合材料薄片的反射以及透射性能,并与其解析解进行对比,从而验证该模型的准确性、收敛性以及高效性。该模型被用来计算三种具有不同电参数的单层碳纤维复合材料薄片的屏蔽性能,以研究各电参数对其屏蔽性能的影响。     

A NEW INTERFACE CAPTURING METHOD BASED ON DOUBLE INTERFACE FUNCTIONS

基于代数重构思想,发展了一种新的双界面函数重构方法,并采用双正弦函数构造了双正弦界面重构方法（double sine interface capturing,DSINC）.为验证不同界面函数对界面捕捉效果的影响,用数值方法求解了可压缩五方程模型,其中对流项的离散采用五阶WENO（weighted essentially non-oscillatory method）格式,时间积分采用三阶Runge——Kutta方法,通量计算分别考虑了HLL和HLLC方法,而状态方程采用Mie-Grüneisen状态方程.在数值计算中,在界面附近,采用DSINC来获得体积分数的重构,而在远离界面的区域采用WENO格式来获得高阶插值状态.相比采用单界面函数的方法,如双曲正切界面重构方法（tangent of hyperbola for interfacecapturing,THINC）,DSINC方法同样具有界面重构算法简单,在程序中添加方便等特点,两者区别在于,DSINC方法在重构过程中未知函数更易于求解,而无需求解复杂的非线性超越方程,这就使其具有易于向多维扩展的能力.一些典型的两相流动问题,如圆形水柱对流问题,两相三波点问题和激波——界面不稳定性问题等被用作不同界面函数对界面捕捉效果的影响对比.对比分析发现,DSINC与THINC在界面捕捉效果上大致保持一致,并在计算中表现出了较好的稳定性.双界面函数重构思想可以为多相流动界面的代数重构提供了一种新的思路.     

JEMS-FDTD软件在运输机电磁特性仿真中的应用

复杂电磁环境的研究包括多方面内容,其中飞行器的电磁环境效应研究受到了较多的关注。随着计算机技术的飞速发展,利用数值模拟技术获得飞行器电磁特性已经成为飞行器电磁环境效应研究的一个重要手段。我们将自主研发的大规模并行三维时域全波电磁模拟软件JEMS-FDTD应用于运输机整机电磁特性仿真,计算并分析了飞机的电磁耦合散射特性,获得了飞行器在电磁脉冲照射下的时域/频域、近场/远场等电磁信息。计算中为保证精度采用了非均匀以及高阶FDTD计算技术。     

JMCT与子通道程序耦合方法研究及验证

蒙特卡罗与热工水力的耦合计算是目前反应堆数值模拟的重要研究方向,在蒙特卡罗方法连续能量点截面的基础上结合热工程序的温度反馈,反应堆中子计算的准确性得到大幅提高。为了提高计算精度,堆芯模型分辨率也需进一步提高,相比于组件均匀化模型,pin-by-pin的建模方式能够获得更好的结果。利用蒙特卡罗程序JMCT与子通道程序COBRA-EN实现了蒙特卡罗-热工的内耦合,内耦合方式通过内存进行数据传递,其计算效率及安全性均优于外耦合方法。随后利用NURISP项目迷你堆的pin-by-pin模型对耦合程序进行验证。计算结果与同类耦合程序相似,验证了程序的准确性。同时,对耦合过程的收敛性问题进行了初步分析。     

基于神经网络的粒子输运问题高效计算方法

蒙特卡罗方法是求解粒子输运问题的有力工具之一,其局限性在于为达到精度要求需模拟大量粒子,计算耗时长,这阻碍了该方法的进一步应用,尤其在需快速响应的情形.本文结合神经网络和若干蒙特卡罗方法基本原理发展了一种计算方法,能够实现源分布可变,几何、材料和目标计数不变的中子输运问题的快速准确求解.首先,为高效生成用于神经网络训练的数据,利用重要性原理实现在同样模拟次数基础上有效扩充训练数据集容量,在一定程度上克服了使用蒙特卡罗计算获取训练数据耗时长的缺点.进而,基于目标计数是源分布与重要性函数乘积积分的事实,设计了利用神经网络实现快速输运计算的策略.该网络的输入是中子源项,输出是目标计数,在几何、材料和目标计数固定的情况下,该神经网络可重复使用,根据新的源项快速准确得到目标计数.本文所提出方法的原理和框架同样适用于其他种类粒子的同类型输运问题.基于若干基准模型的验证表明,训练得到的神经网络能在不到1 s的时间内得到目标计数,且与蒙特卡罗大样本模拟得到基准结果的平均相对偏差均低于5%.