炸药冲击响应的二维细观离散元模拟

用离散元法在颗粒尺度模拟了塑料粘接炸药的冲击响应。炸药的细观结构基于Voronoi拼图建成。以含孔洞和不含孔洞的炸药模型为例,说明细观结构对热点的形成和分布有重要影响。在这些算例中,粘塑性变形是温度局域化的主要原因,而摩擦和体积粘性等机制可以忽略。     

多孔极化PZT95/5铁电陶瓷单轴压缩力学响应与放电特性

采用添加造孔剂的方法制备了四种不同孔隙率PZT95/5铁电陶瓷,对其进行电场极化,随后开展了准静态单轴压缩实验,讨论了畴变、相变以及孔隙率对极化PZT95/5铁电陶瓷的力学响应与放电特性的影响.研究结果表明:(1)多孔极化PZT95/5铁电陶瓷非线性力学响应行为主要归因于畴变和相变的共同作用,与内部孔洞变形和坍塌基本无关;(2)在准静态单轴压缩下极化PZT95/5铁电陶瓷的去极化机制是畴变和相变的共同作用;(3)孔隙率对极化PZT95/5铁电陶瓷的弹性模量、压缩强度有明显的影响,而对断裂应变的影响较小;(4)极化PZT95/5铁电陶瓷畴变和相变开始的临界应力都随着孔隙率的增大而线性衰减,但相变开始的临界体积应变却不依赖孔隙率;(5)极化PZT95/5铁电陶瓷电荷饱和释放量随着孔隙率呈线性减小,但孔隙率对电荷释放速率基本没有影响.     

损伤统计演化方程的性质和数值模拟

通过对一种含成核尺寸效应的损伤统计演化方程性质的分析和数值模拟,揭示了损伤率主要是由微损伤在二维相空间中的前沿的运动所决定的这也就是Kachanov提出的损伤率演化方程的物理基础数值结果进一步显示了含成核尺寸效应模型在损伤发展上与-维模型的区别而且,由几种形式的细观动力学算出的损伤率与损伤的关系简单,可近似拟合为宏观上封闭的形式     

沥青混合料车辙变形的离散元数值模拟

采用Burgers模型描述沥青基质的粘弹性,根据实验数据确定Burgers模型参数.借助数字图像处理方法建立试件几何模型,利用离散元法开展了车辙变形的离散元数值模拟,研究了沥青混合料在不同骨料分布、不同载荷和不同温度条件下的车辙深度变化情况.结果表明:温度对车辙变形的影响在沥青软化点附近较大,在离开软化点的范围较小;当温度较低时,混合料的整体性比较好,骨料分布对局部车辙变形的影响比较小,超载引起的车辙增量不太明显,而当温度达到或超过沥青软化点时,混合料的整体性明显下降,骨料分布对局部车辙变形的影响非常突出,超载引起的车辙增量也非常显著.因此,当超载和高温组合到一起时,车辙变形会大幅度增加,将给路面安全带来严重隐患.     

离散元数值模拟中查找邻居元关系的改进算法

离散元方法是近年发展起来的有前途的新概念数值模拟方法,然而查找邻居元关系耗费大量机时和内存。本文引入一种改进算法-“区域法”,并对其原理、算法、运行时间和空间进行了分析。和原算法中的“窗口法”的比较结果表明,新算法在运行时间和空间方面均明显优于原算法。     

基于离散元的杆件断裂全过程数值模拟

杆件的断裂会涉及到大变形、非线性以及不连续等问题,通常的数值计算方法模拟这种复杂力学行为具有局限性。本文基于颗粒离散元法DEM,将接触粘结处的分布式弹簧用梁纤维进行等效,提出了一种适于结构弹塑性分析的DEM纤维梁模型,然后在此基础上构建了构件断裂模拟算法以及纤维破环准则。将该模型应用于悬臂梁结构,模拟了悬臂梁从弹性到弹塑性阶段,再到断裂破坏的全过程,数值模拟得到的结构响应和截面开裂破坏形态均较合理。最后将该方法应用于单层网壳倒塌破坏模拟,并与网壳振动台倒塌试验进行对比,结果表明,数值模拟得到的杆件断裂过程及结构倒塌模式与试验现象一致,验证了该模型的正确性和适用性。     

含热冲击预损伤的陶瓷基复合材料损伤本构模型

陶瓷基复合材料结构在服役过程中不可避免地经受热冲击(较高的热应力梯度)而产生热机械损伤, 因此, 建立含循环热冲击预损伤材料的损伤本构模型, 以描述材料在热机械载荷作用下的力学行为, 对材料结构损伤容限设计与结构完整性评估非常重要. 本文首先对经历了循环热冲击的材料进行单调拉伸损伤实验, 发现对于含循环热冲击预损伤的材料, 其弹性模量的下降与所施加的应变直接相关. 然后在连续介质损伤力学的框架下, 基于平面应力假设, 建立了含循环热冲击预损伤材料的损伤演化模型, 该模型所涉及的参数可通过一个偏轴(45$^\circ$)以及两个正轴(平行于两个主方向)的单调拉伸试验获得. 最后, 采用经典塑性理论对由基体损伤引起的非弹性应变进行了描述. 本文所提出的应变损伤宏观模型可以描述陶瓷基复合材料在热机械载荷作用下的损伤演化, 同时弥补了含预损伤的陶瓷基复合材料在机械载荷下损伤本构模型在理论及实验研究方面的不足.      

半结晶聚合物损伤演化的实验表征与数值模拟

损伤本构模型对研究材料的断裂失效行为有重要意义, 但聚合物材料损伤演化的定量表征实验研究相对匮乏. 通过4种高密度聚乙烯(high density polythylene, HDPE)缺口圆棒试样的单轴拉伸实验获得了各类试样的载荷-位移曲线和真应力-应变曲线, 采用实验和有限元模拟相结合的方法确定了HDPE材料不同应力状态下的本构关系, 并建立了缺口半径与应力三轴度之间的关系;采用两阶段实验法定量描述了4种HDPE试样单轴拉伸过程中的弹性模量变化, 并建立了基于弹性模量衰减的损伤演化方程, 结合中断实验和扫描电子显微镜分析了应力状态对HDPE材料微观结构演化的影响. 结果表明缺口半径越小, 应力三轴度越大, 损伤起始越早、演化越快; 微观表现为: 高应力三轴度促进孔洞的萌生和发展, 但抑制纤维状结构的产生;基于实验和有限元模拟获得的断裂应变、应力三轴度、损伤演化方程等信息提出了一种适用于聚合物的损伤模型参数确定方法, 最后将本文获得的本构关系和损伤模型用于HDPE平板的冲压成形模拟, 模拟结果与实验结果吻合良好.      

高速铁路碎石道砟振动的离散元模拟

道砟振动对其磨损、破碎和道床累积变形有显著影响,为揭示高速车辆移动荷载作用下道砟动态响应特性,建立有砟道床离散元模型,开展车辆-轨道耦合动力学计算得到离散元模型输入荷载,模拟分析高速车辆以不同速度通过时有砟道床的振动响应,并与车辆-轨道耦合动力学计算结果进行对比分析。结果表明,轨枕、道砟和道床块振动位移波形相似,位移幅值沿道床深度方向减小,道床块振动位移与轨枕底面以下0.3m处道砟的振动位移相当;轨枕、道床块振动速度与加速度随行车速度提高而增大;受道砟颗粒间复杂相互作用的影响,道砟振动加速度会出现突变。道床离散元模型能合理反映道砟颗粒的振动响应特性,道床块模型体现了道床层在有砟轨道结构中的动力传递与减振特性,两种道床模型的计算结果具有一定的相似性。     

不排水条件下土石混合料双轴压缩离散元模拟

基于离散单元法构建了考虑碎石形状的柔性边界不排水双轴压缩模型,研究了11组不同含石量土石混合料在不排水条件下的宏细观力学特性。研究发现,(1)土石混合料强度随含石量增加先增后减,含石量超过约60%时形成碎石骨架,强度快速增长;含石量约80%时强度和内摩擦角最大;含石量增加,混合料黏聚力逐渐减小。(2)不排水受载下,土石混合料孔隙水压力先正后负、先增后减。含石量较小时,含石量增加,混合料孔隙水压力变化较小,混合料在碎石骨架形成初期孔隙水压力增加;碎石骨架完全形成后最大孔隙水压力减小。(3)含石量小于约30%时,混合料应变局部化呈规则剪切带,随后含石量增加破坏规则剪切带,当含石量约为80%时,混合料再次形成规则剪切带。含石量增加,配位数逐渐减小,颗粒间接触力在含石量小于60%时随含石量增加缓慢增长,超过70%时则快速增长。     

Mechanical behaviors and damage constitutive model of ceramics under shock compression

One-stage light gas gun was utilized to study the dynamic mechanical properties of AD90 alumina subjected to the shock loading. Manganin gauges were adopted to obtain the stress-time histories. The velocity interferometer system for any reflector （VISAR） was used to obtain the free surface velocity profile and determine the Hugoniot elastic limit. The Hugoniot curves were fitted with the experimental data. From Hugoniot curves the compressive behaviors of AD90 alumina were found to change typically from elastic to ＂plastic＂. The dynamic mechanical behaviors for alumina under impact loadings were analyzed by using the path line principle of Lagrange analysis, including the nonlinear characteristics, the strain rate dependence, the dispersion and declination of shock wave in the material. A damage model applicable to ceramics subjected to dynamic compressive loading has been developed. The model was based on the damage micromechanics and wing crack nucleation and growth. The effects of parameters of both the micro-cracks nucleation and the initial crack size on the dynamic fracture strength were discussed. The results of the dynamic damage evolution model were compared with the experimental results and a good agreement was found.     

单轴压缩下页岩类脆性材料细观损伤演化的离散元分析

了解广泛存在的类似页岩的脆性材料各向异性对工程安全具有重要意义。本研究将页岩视为粘结颗粒材料,基于离散单元方法研究了横观各向同性脆性页岩的损伤演化。再现了不同层理角的页岩试样的破坏模式,并对比了实验和数值模拟的抗压强度和弹性模量。引入微裂纹的概念,通过定义裂纹密度函数,系统地研究了单轴压缩条件下,页岩层理角对细观结构的影响。此外基于平均配位数建立了配位数变化与细观损伤的联系,并根据配位数的变化与裂纹数量将加载过程分为三个阶段,分析了不同阶段配位数与裂纹数量的对应变化关系。研究表明,页岩的裂纹密度随着层理角的增加而增加,而试样的平均配位数在加载过程中先上升后剧烈下降,颗粒集合体在单轴压缩条件下的应力应变及裂纹数量曲线与平均配位数曲线有良好的一致性。该研究揭示了横观各向同性脆性岩石的破坏过程和内在机理,将为页岩类脆性材料的工程应用提供理论指导。     

Discrete element simulation of shock wave propagation in polycrystalline copper

The implementation of the characteristic of compressive plasticity into the Discrete Element Code, DM2, while maintaining its quasi-molecular scheme, is described. The code is used to simulate the shock compression of polycrystalline copper at 3.35 and 11.0 GPa. The model polycrystal has a normal distribution of grain sizes, with mean diameter 14 μm, and three distinct grain orientations are permitted with respect to the shock direction; 〈1 0 0〉, 〈1 1 0〉, and 〈1 1 1〉. Particle velocity dispersion (PVD) is present in the shock-induced flow, attaining its maximum magnitude at the plastic wave rise. PVD normalised to the average particle velocity of and are yielded for the 3.35 and 11.0 GPa shocks, respectively, and are of the same order as those seen in the experiment. Non-planar elastic and plastic wave fronts are present, the distribution in shock front position increasing with propagation distance. The rate of increase of the spread in shock front positions is found to be significantly smaller than that seen in probabilistic calculations on nickel polycrystals, and this difference is attributed, in the main, to grain interaction. Reflections at free surfaces yield a region of tension near to the target free surface. Due to the dispersive nature of the shock particle velocity and the non-planarity of the shock front, the tensile pressure is distributed. This may have implications for the spall strength, which are discussed. Simulations reveal a transient shear stress distribution behind the shock front. Such a distribution agrees with that put forward by Lipkin and Asay to explain the quasi-elastic reloading phenomenon. Simulation of reloading shocks show that the shear stress distribution can give rise to quasi-elastic reloading on the grain scale.     

Discrete element modeling on the crack evolution behavior of brittle sandstone containing three fissures under uniaxial compression

Based on experimental results of brittle, intact sandstone under uniaxial compression, the micro-parameters were firstly confirmed by adopting particle flow code(PFC2D). Then, the validation of the simulated models were cross checked with the experimental results of brittle sandstone containing three parallel fissures under uniaxial compression. The simulated results agreed very well with the experimental results, including the peak strength, peak axial strain, and ultimate failure mode. Using the same microparameters, the numerical models containing a new geometry of three fissures are constructed to investigate the fissure angle on the fracture mechanical behavior of brittle sandstone under uniaxial compression. The strength and deformation parameters of brittle sandstone containing new three fissures are dependent to the fissure angle. With the increase of the fissure angle, the elastic modulus, the crack damage threshold,and the peak strength of brittle sandstone containing three fissures firstly increase and secondly decrease. But the peak axial strain is nonlinearly related to the fissure angle. In the entire process of deformation, the crack initiation and propagation behavior of brittle sandstone containing three fissures under uniaxial compression are investigated with respect tothe fissure angle. Six different crack coalescence modes are identified for brittle sandstone containing three fissures under uniaxial compression. The influence of the fissure angle on the length of crack propagation and crack coalescence stress is evaluated. These investigated conclusions are very important for ensuring the stability and safety of rock engineering with intermittent structures.     

A95陶瓷材料的动态压缩测试研究

为了得到陶瓷材料有效的实验结果,必须采用多种改进措施(PulseShaper、高强度和高刚度垫块、直接测试试样应变以及使用较大长径比的试样)来满足SHPB装置所要求的基本假设条件。研究以A95陶瓷的测试为例,通过数值模拟分析证明:针对陶瓷材料改进的SHPB装置,可以得到陶瓷材料的有效压缩应力 应变曲线和压缩破坏强度。但是在分析材料的压缩破坏强度时,必须排除劈开破坏的干扰。数值模拟研究发现,压缩测试中出现的劈开破坏是由卸载过程中出现的局部横向(沿直径)拉伸造成的拉伸破坏。在修正了弥散对压缩破坏强度的影响后,研究得到了A95陶瓷的压缩破坏强度为3.5GPa。A95陶瓷的高应变率压缩应力 应变曲线显示了良好的线弹性,并且曲线的斜率与杨氏模量相同。     

A theory of the densification-induced fragmentation in glasses and ceramics under dynamic compression

A boundary of failure which follows a shock front is observed in glasses and ceramics above a critical compressive shock load. This boundary, called the failure wave leaves behind a damaged material with newly evolved properties. These include the Tresca yield behavior and the reduction in sound speed. The evolution of the Mescall zone during long-rod penetration of these materials is associated with this wave. But, the failure wave and this accompanying process of fragmentation under dynamic compression are still not understood. It is known that some brittle solids undergo an irreversible density increase when subjected to high compression. This phenomenon, called densification is linked to the formations of slip lines and cracks in intensely compressed regions of silica glass. It also corresponds directly to the losses in shock wave speeds. Once densified, a region tends to shrink, straining the interface between it and the original solid. Stressed interfaces are unstable and may roughen, causing local cracking. On this basis, the failure wave is idealized as a propagating fracture boundary layer where the solid is comminuted by a process of densification interface roughness. The kinetics for this process are established using the fluctuation dissipation theorem. Shear and tensile modes of fragmentation are studied in plane stress. The theory predicts the powder size in the Mescall zone of silica glass. Nevertheless, this theory still needs experimental verification.