二维骨料随机分布混凝土的动态力学性能数值模拟北大核心CSCD

将混凝土看成粗骨料和水泥砂浆组成的二相非均质复合材料。根据富勒级配曲线和瓦拉文平面转化公式,编写了二维圆形骨料随机分布程序,对混凝土在冲击荷载作用下的动力响应进行了数值模拟。分析冲击速度、试件尺寸、粗骨料大小及分布和粗骨料体积分数对混凝土动态力学性能的影响,讨论了混凝土的冲击破坏模式。数值模拟表明:混凝土的峰值应力随着冲击速度的增大而增大,具有明显的率效应;随着模型尺寸的增大而减小,表现出明显的尺寸效应;随着粗骨料体积分数增大,冲击荷载作用下混凝土的峰值应力呈现先增大后减小的趋势,粗骨料体积分数为40%时混凝土峰值应力最大;保持粗骨料最大粒径不变,随着粗骨料最小粒径的增大,混凝土的峰值应力逐渐减小;保持粗骨料最小粒径不变,随着粗骨料最大粒径的增大,混凝土的峰值应力呈现先增大后减小的趋势。数值模拟结果为混凝土的工程应用提供了理论依据和技术支撑。     

混凝土材料动态力学性能实验与数值模拟研究

 利用轻气炮对素混凝土材料进行了3种不同冲击速度下的动态力学性能实验,得到了时间-应力曲线,随后利用软件进行数值模拟,将模拟结果与实验结果进行比较,解释了飞片和靶板破坏过程中波的传播过程。通过分析发现：混凝土材料表现出明显的率相关性、滞后效应和应力波的衰减特征,这些与材料内部损伤的演化密切相关;混凝土材料在压缩波、边侧稀疏效应和左、右自由表面反射拉伸波的共同作用下,最终形成宏观破碎;实验中靶板尺寸设计合理,有效避免了边侧稀疏效应对锰铜压阻计测试信号的影响。     

随机介质背散射二维Mueller矩阵的数值模拟与分析

推导了随机介质背散射Mueller矩阵的直接计算公式,并运用矢量Monte Carlo方法进行了数值模拟.结果表明随机介质背散射二维Mueller矩阵方位关系随散射系数的减小而增强,而与微粒大小关系不大;Mueller矩阵元素绝对值的空间分布随径向呈近似指数规律衰减,矩阵元素的方位变化具有周期性.对称系统的二维Mueller矩阵的花样图中仅有7幅独立,其余9幅可通过对称、旋转变换得到.     

二维磁通随机钉扎系统动力学行为的数值模拟

本文利用数值模拟的方法,研究在Ｌｏｒｅｎｔｚ 驱动力作用下,随机钉扎环境中运动的二维磁通点阵的动力学行为．通过对不同随机分布点钉扎中心密度和不同磁通间相互作用强度的研究,发现磁通点阵存在从弹性区向塑性区的渡越（ｃｒｏｓｓｏｖｅｒ） , 并且这种渡越强烈依赖于磁通之间及磁通与钉扎之间的相互作用．随着钉扎中心密度的增加及磁通间相互作用强度的减弱,临界钉扎力增强．模拟结果与最近实验结果定性一致．部分结果偏离集体钉扎理论     

模糊随机参数二维声场数值计算

解决声场参数同时具有模糊性和随机性的问题,实现模糊随机声场声压响应的预测,引入了信息熵理论,利用信息熵的等效转换,将模糊随机声场转化为纯随机声场或者纯模糊声场进行求解,推导了基于摄动法的二维随机声场和模糊声场的有限元计算公式。以模糊随机参数下的二维管道声场模型和某轿车二维声腔模型为例进行了数值计算,所得结果与蒙特卡洛法(Monte Carlo Method)所预测声压变化范围基本一致,同时,转化为纯随机声场和纯模糊声场所求得声压响应变化范围也基本一致,说明了本文方法计算结果的准确性。因此本文方法能很好地应用于模糊随机参数下二维声场的预测,具有重要的工程应用价值。     

模糊随机参数二维声场数值计算

解决声场参数同时具有模糊性和随机性的问题,实现模糊随机声场声压响应的预测,引入了信息熵理论,利用信息熵的等效转换,将模糊随机声场转化为纯随机声场或者纯模糊声场进行求解,推导了基于摄动法的二维随机声场和模糊声场的有限元计算公式。以模糊随机参数下的二维管道声场模型和某轿车二维声腔模型为例进行了数值计算,所得结果与蒙特卡洛法(Monte Carlo Method)所预测声压变化范围基本一致,同时,转化为纯随机声场和纯模糊声场所求得声压响应变化范围也基本一致,说明了本文方法计算结果的准确性。因此本文方法能很好地应用于模糊随机参数下二维声场的预测,具有重要的工程应用价值。     

Z箍缩丝阵动态烧蚀的二维数值模拟

采用FOI-PERFECT程序对局域欧姆加热主导的动态烧蚀进行了数值模拟研究,定量诊断并比较了烧蚀率、先驱电流和烧蚀时间,考查烧蚀阶段之初的丝等离子体的直径和丝阵的丝数对烧蚀的影响,结果表明,在等离子体流建立后,随丝数或者丝径的增加,烧蚀率增大,先驱电流减少,烧蚀时间减小。     

沙漏计时原理二维数值模拟

提出沙漏计时二维数值模型,并用元胞自动机的方法对该模型进行计算机模拟.结果表明由于开口处沙粒的流动而引起的沙堆崩塌行为是一个准周期性振荡过程,开口大小与崩塌总时间近似成线性关系.特殊的开口位置会引起崩塌"拍",且随沙堆高度呈周期性变化.     

二维弹塑性磁流体力学数值模拟

为了研究物质弹塑性对磁驱动实验运动过程、不稳定性发展等的影响,在MDSC2程序的基础上,增加了弹塑性模块,研制了包括弹塑性的磁流体力学程序,并进行了弹塑性项影响的数值模拟和分析。数值模拟表明：没有初始扰动时,弹塑性项几乎不影响套筒内外半径的运动轨迹;有初始扰动时,弹塑性项对磁驱动固体套筒的Rayleigh-Tayor不稳定性有明显的抑制作用。     

Z箍缩丝阵动态烧蚀的二维数值模拟

采用FOI-PERFECT程序对局域欧姆加热主导的动态烧蚀进行了数值模拟研究, 定量诊断并比较了烧蚀率、先驱电流和烧蚀时间, 考查烧蚀阶段之初的丝等离子体的直径和丝阵的丝数对烧蚀的影响, 结果表明, 在等离子体流建立后, 随丝数或者丝径的增加, 烧蚀率增大, 先驱电流减少, 烧蚀时间减小。     

自由电子激光振荡器的二维数值模拟

在柱二维坐标系中,用数值计算方法求解了摇摆器内电子模拟方程组、光场方程组和无源区(不含摇摆器区)经非自适应方法变换后的无源旁轴波动方程并设计了相应的程序R_2D。用该程序计算的数值结果表明,在有源条件下二维数值模拟结果与高斯光束近似下的数值模拟结果符合较好;取完全相同的计算条件R-2D程序与洛斯阿拉莫斯实验室(LANL)的FELEX程序的计算结果符合较好;初步认为LANL实验数据的物理图象合理,数值结果基本可信;最后给出北京自由电子激光器(BFEL)的一些物理数据。     

二维二相不规则网格的油藏数值模拟

本文采用不规则网格的差分方法,对平面油水模型进行了数值模拟,并用IMPES方法计算了两个实例。例1的计算结果与精确解吻合良好,比Pedrosa的局部加密网格法使用的网格数少而且计算更简便。例2为矩形网格难以计算的问题,应用本文方法也得出了较好的结果。这表明本文的不规则网格差分方法是可靠而简便的,它可以有效地应用于复杂外边界和复杂地层结构的油藏模拟问题。     

锥形二维光子晶体太阳电池数值模拟

采用时域有限差分方法,模拟研究在本征吸收层引入锥形二维光子晶体(2D PC)后,其结构参数变化对单结微晶硅电池各膜层吸收的影响规律.研究表明,2D PC的纵横比(高度与周期之比)对电池本征吸收具有决定性影响.周期小于1μm时,本征吸收随着纵横比的增大先上升后下降,纵横比为1时达到最大值;周期大于1μm时,本征吸收达到最大值的纵横比小于1,且周期越大,实现本征吸收最大化的纵横比越小.当周期为0.5μm,纵横比为1时,锥形2D PC电池的本征吸收达到峰值,短路电流密度为27.8 mA/cm2;与平面结构相比,短路电流密度提升5.8 mA/cm~2,相对增加27%.该研究突破了以往认为绒面陷光效果主要取决于绒面形貌横向特征尺寸的观点,对实验获取最佳的周期或随机绒面陷光结构具有指导意义.     

二维翼型跨音速绕流的数值模拟

本文结合高精度TVD格式的数值通量和时间进展多步法给出了一种求解定常流问题的数值方法。同时给出一些特殊处理来加快数值解的收敛速度。本文用以上方法计算了翼型跨音速绕流问题,结果表明此方法具有分辨率高,收敛速度较快之性质。     

二维颗粒运动以及传热的数值模拟

本文提出了直接热源算法(Direct heat source)使得颗粒壁面温度满足Dirichlet边界条件,同时采用Wang等提出的内嵌边界方法(Immersed boundary method)和多重直接力算法(Multi-direct forcing)全尺度模拟二维颗粒的运动.本文计算了运动的冷颗粒在不同Grashof数时的运动状态.本文的计算结果与其他学者的计算结果符合得很好,从而验证了本文采用的方法的精确性与有效性.     

超音速连续波COIL的二维半数值模拟

给出了超音速连续波氧碘化学激光器的理论模型，利用此模型对典型实验装置进行了二维半几何近似下的初步数值模拟研究，对所得结果进行了讨论，并与实验或文献结果作了比较。     

自锻破片形成的二维数值模拟

自锻破片装药(简称SFF)是从聚能破甲装药派生出来的一种新的反坦克装药结构。描述SFF的形成过程,我们采用二维轴对称弹性-理想塑性-流体力学方程组。基本数值方法采用M.L.Wilkins的迴路积分法和Legrange四边形网格。数值模拟得到的SFF完整的形成过程和试验中x射线高速摄影照片基本相符。数值模拟中的若干方法简化了计算,节省了机时。     

非牛顿流体二维溃坝问题的数值模拟

移动粒子半隐式法(Moving Particle Semi-implicit,MPS)是一种广泛应用于不可压缩自由表面流动的粒子方法。本文通过把非牛顿流体转化为变黏度牛顿流体的处理方法将MPS方法拓展至非牛顿自由表面流动,以Casson流体和Cross流体为例计算了非牛顿流体的二维溃坝问题。将计算结果与前人数据进行了对比,结果吻合较好。同时对比了非牛顿流体与牛顿流体的计算结果,发现非牛顿流体的溃坝前端发展速度较慢。     

水下爆炸特性的二维数值模拟研究

 对带有3 mm左右厚铝壳或钢壳的柱形TNT、RS-211、T/Γ（35/65）以及ROT901装药在水下采用一端点起爆后,水中冲击波超压分布进行了二维数值模拟研究。对该四种炸药的每一种,分别考虑了药量为1 kg、3 kg、5 kg和8 kg四种TNT当量的柱形装药情况。柱形装药的长径比为L/D=1.5左右。计算使用的程序为二维欧拉多流体网格法流体力学计算程序MFIC。给出了水中冲击波超压的空间和时间分布;对铝壳和钢壳两种情况的结果进行了分析、比较;同时给出了缩比模型与1∶1模型之间相应物理量之间的缩比关系。     

二维随机蜂巢网格熔断动力学过程和熔断面标度性质的数值模拟

石墨烯等材料具有典型的二维蜂巢结构,而随机电阻丝模型则是研究非均匀材料断裂十分有效的统计物理学模型.本文尝试对二维蜂巢结构随机电阻丝网络熔断的动力学过程及熔断面性质进行数值模拟分析,以此来研究二维非均质蜂窝材料熔断的动力学性质和熔断面的动力学标度性质.模拟研究表明,二维随机蜂窝网格的熔断动力学过程和熔断面具有明显的标度性质,得到的熔断面整体和局域粗糙度指数分别为α=0.911±0.005和α_(loc)=0.808 ± 0.003,这两者之间的明显差异表明熔断面具有奇异标度性.通过对熔断面极值高度的分析发现,熔断面高度的极值统计分布能很好地满足Asym2sig型分布,而不是最常见的三种极值统计分布.本文的研究表明,随机电阻丝模型在模拟非均匀材料的电流熔断过程和熔断表面标度性的分析中同样适用和有效.