A modified Stillinger–Weber potential-based hyperelastic constitutive model for nonlinear elasticity

The Stillinger–Weber (SW) potential, which is a combination of the two- and three-body interaction, states that the bond energy is not only related to the distance between atoms, but also related to the bond angles subtended by this given bond and other bonds. The bond energy mechanism presented by the SW potential is different from that by the classical potentials, such as the Lennard–Jones, Tersoff and Embedded Atom potentials. Different micro energy mechanism reveals different micro fracture mechanism. The original SW potential takes the ‘ideal’ tetrahedral angle as the reference value of each bond angle in the current configuration, which makes it only applicable to the silicon materials. However, the micro fracture mechanism revealed by the SW potential should not be confined to the silicon materials. To extend the SW potential to a wider range of materials, the value of the bond angle in the reference configuration, not the ‘ideal’ tetrahedral angle, is taken as the reference value of this bond angle in the current configuration. Based on this modified SW potential, a constitutive model is developed. By this way, the micro fracture mechanism invoked by the SW potential is incorporated into the constitutive relation. Through this proposed constitutive model, it is found that the Hookean matrix derived from the SW potential matches that of a linear elastic continuum, which suggests that there exists a corresponding relationship between the micro physical parameters of SW potential and the macro material constants. The corresponding micro–macro parameter relationship is derived in this paper. To examine the application of this method to other materials, it is used to simulate the mixed fracture growth in concrete under static and dynamic load. The simulation results suggest that the present method can capture the characters of fracture growth in the quasi-brittle materials. It suggests that the constitutive model based on the modified SW potential can be applied to materials other than silicon. Because the interatomic potential-incorporated constitutive model makes the fracture criterion as the intrinsic property of a constitutive relation, it presents many advantages in fracture simulation. This paper enriches the constitutive relation with the micro fracture mechanism presented in SW potential, providing a new micro constitutive model for materials. Besides this, it also provides a feasible approach to calibrating the parameters of the SW potential for a certain material.     

Discretized virtual internal bond model for nonlinear elasticity

The virtual internal bond (VIB) is a micro–macro constitutive model. Although this model is based on a postulated discrete microstructure, it ultimately returns to a continuum constitutive relation through a homogenization process. The homogenization process can reduce the internal degrees of freedom, but it omits the effect of the individual micro bond that may play an important role in the fracture process. The present research develops a discrete system to represent the nonlinear elasticity by discretizing the continuous VIB. This discrete system is composed of unit cells, which can adopt any geometry with any number of bonds. The system is characterized by the force–displacement, not the stress–strain constitutive relationship. The nonlinear properties of this discrete system are governed by the micro-bond potential. The micro bond properties are related to Young’s modulus of the material, the volume and the bond number of the unit cell. For a given material, the unit cell has a certain topological structure and configuration. A discussion of two specific cases (the 2D triangular and 3D tetrahedral unit cells) suggests that the discrete system converges with decreasing unit cell size. In the unstructured unit cell scheme, the discrete system can almost precisely represent the initial Young’s modulus and the Poisson ratio of a nonlinear continuum. A mixed fracture example demonstrates that the present method can efficiently simulate the fracture propagation. The present paper provides a theory for developing a lattice-type mechanical model for nonlinear elasticity and provides new method for the fracture simulation of a nonlinear elastic material.     

A new quasi-continuum constitutive model for crack growth in an isotropic solid

A new quasi-continuum constitutive model is established based on the randomized cohesive bonds model proposed by Gao and Klein (1998). This model bridges the microscopic discrete constitution characters and the macroscopic mechanical properties of material. In the presented constitutive model, both the bond stretch energy potential and the rotation energy potential are considered, which makes the presented constitutive model applicable to different Poisson-ratio and Young's modulus materials. By establishing a phenomenological bond stiffness function according to the complete stress–strain relationship of uniaxial tension test, the fracture criterion is directly incorporated into the constitutive model. The method requires no external fracture criterion when simulating fracture initiation and propagation, which brings convenience in the numerical simulation. At last, the presented constitutive model is applied to an example of crack growth in an isotropic solid.     

Modified virtual internal bond model based on deformable Voronoi particles

In last time, the series of virtual internal bond model was proposed for solving rock mechanics problems. In these models, the rock continuum is considered as a structure of discrete particles connected by normal and shear springs(bonds). It is well announced that the normal springs structure corresponds to a linear elastic solid with a fixed Poisson ratio, namely, 0.25 for threedimensional cases. So the shear springs used to represent the diversity of the Poisson ratio.However, the shearing force calculation is not rotationally invariant and it produce difficulties in application of these models for rock mechanics problems with sufficient displacements. In this letter, we proposed the approach to support the diversity of the Poisson ratio that based on usage of deformable Voronoi cells as set of particles. The edges of dual Delaunay tetrahedralization are considered as structure of normal springs(bonds). The movements of particle's centers lead to deformation of tetrahedrals and as result to deformation of Voronoi cells. For each bond, there are the corresponded dual face of some Voronoi cell. We can consider the normal bond as some beam and in this case, the appropriate face of Voronoi cell will be a cross section of this beam. If during deformation the Voronoi face was expand, then, according Poisson effect, the length of bond should be decrees. The above mechanism was numerically investigated and we shown that it is acceptable for simulation of elastic behavior in 0.1–0.3 interval of Poisson ratio. Unexpected surprise is that proposed approach give possibility to simulate auxetic materials with negative Poisson's ratio in interval from –0.5 to –0.1.     

含初始裂纹巴西圆盘劈裂问题的非局部近 场动力学建模

巴西圆盘劈裂是弹性力学及岩石力学与工程中的经典问题。在非局部键型近场动力 学理论的基础上,引入物质点对的转动自由度构建双参数微观弹脆性近场动力学本构力模型 以突破常规模型的应用范围限制,并考虑岩石混凝土类材料的宏观拉压异性和断裂特征。引 入动态松弛、粒子系统力边界条件和系统平衡弛豫等算法,实现了含不同倾角中心裂纹巴西 圆盘受压劈裂破坏全过程的近场动力学数值模拟,裂纹扩展路径及破坏形式均与试验结果高 度吻合,为裂纹扩展和断裂破坏问题的数值模拟提供了新的选择。     

岩石边坡裂隙渗流的流形元模拟

为了能够更好地考虑裂隙渗流在岩石边坡破坏中的作用,在考虑裂纹扩展的数值流形方法程序基础上,从最小势能原理出发,阐述了渗流与变形的耦合作用机理并推导了相应的耦合方程,在程序中实现了对渗流与断裂的模拟。最后利用该程序对含初始裂隙的岩石边坡在渗流作用下的破坏过程进行了模拟。模拟结果很好地再现了边坡在渗流作用下的实际破坏过程。同时通过对裂纹扩展过程中应力强度因子的计算发现,第一应力强度因子始终大于第二应力强度因子,且随着计算时间增加,二者的差值逐渐增大,裂纹扩展以Ⅰ型为主。     

综放顶煤爆破能量的分形研究

基于综放开采预爆破弱化坚硬顶工艺,利用分形理论建立岩石断耗能量与块度分维值之间的理论关系,推导出岩石爆破块度分数维与炸药单耗之间关系,进行了室内相似材料模型爆破实验,验证了上述理论关系,利用爆破岩石表面裂隙分布的分形维数值作了衡量指标,优化爆破孔网参数。     

考虑粘结面厚度的局部粘结单元法

传统无厚度粘结单元法CFEM (Cohesive finite element method)在模拟脆性材料断裂方面具有很强的优势,但也存在很大问题.一是单元尺寸增大,收敛性变差;二是单元尺寸变小,模型刚度发生折减.为了克服这两个问题,发展了考虑厚度的局部粘结单元法,即在裂纹可能扩展区插入具有一定厚度的粘结面单元.粘结面单元采用拓展虚内键本构(Augmented virtual internal bond)描述.由于考虑了厚度,粘结面交叉处会形成多边形空缺.为了弥补这一空缺,将其看作多边形键元胞,采用离散虚内键模型(Discretized virtual internal bond)对其建模,保证了模型的几何完整性.模拟结果表明,本文方法有效,克服了传统CFEM方法的刚度折减问题,提高了计算稳定性和收敛性.     

块石形状对土石混合体力学行为影响的颗粒流模拟

考虑块石形状为球体、正方体和长方体三种情况,通过正方体与球体相比较来探究块石棱角度不同对土石混合体力学特性的影响,通过长方体与正方体相比较来探究块石球度不同对土石混合体力学特性的影响。首先,提出特定形状块石三维离散元精细建模的方法;接着建立含石量为30%和80%的块石形状分别为球体、正方体和长方体的土石混合体三维颗粒流数值模型;然后,对土石混合体大三轴试验进行颗粒流模拟,获得了不同含石量及不同块石形状的土石混合体试样的宏观力学特征;最后,详细分析了块石形状对土石混合体力学行为影响的细观机理。结果表明,块石形状会影响土石混合体的力学行为,其影响的程度与含石量密切相关;配位数、块石颗粒平均旋转量和摩擦功的演化规律能够很好地从细观水平上反映块石形状的作用机理。     

脆塑性岩石破坏后区应力跌落效应数值模拟

采用三轴试验和数值模拟研究了岩石类脆塑性材料的应力跌落效应,并用塑性流动因子λ来描述应力跌落效应。为简化计算,给出了一种应力陡降过程中伴随的非零应变增量的工程近似处理方法,基于此针对性地编制了处理脆塑性材料应力跌落的有限元分析程序代码。数值模拟结果验证了塑性流动因子λ以及该近似处理方法的有效性。     

Numerical simulation of rock fragmentation during cutting by conical picks under confining pressure

In this article, the effect of confining pressure on rock fragmentation process during cutting was investigated by numerical simulation with a discrete element method (DEM). Four kinds of sandstones with different physical properties were simulated in the rock cutting models under different confining pressures. The rock fragmentation process, the cutting force, and the specific energy under different confining pressures were analyzed. With the increase in confining pressure and rock strength, the vertical propagation of cracks was restrained. Rock samples were compacted and strengthened by confining pressure resulting in the increase of the cutting force. The specific energy of rock cutting linearly increased with the increase of the confining pressure ratio.     

Determination of Bond Strengths in Non-woven Fabrics: a Combined Experimental and Computational Approach

Interfiber bonds are important structural components in non-woven fabrics. Bond fracture greatly affects the strength and damage progression in a fiber network structure. Here, we present a novel combined experimental and computational approach to extract bond strengths in non-wovens. In this method, a small specimen is imaged and the obtained 3D geometry of the network is directly modeled in a finite element framework. Bond properties are determined by matching finite element simulation predicted mechanical response to the experimental data. This method is demonstrated by applying it to six specimens of a commercial polypropylene non-woven. A four parameter bi-linear interface law is used with normal stiffness k, shear stiffness βk, separation at the start of damage d ₁, and separation at total loss of bond stiffness d ₂. The determined normal strength (kd ₁)and shear strength (βkd ₁) are (1.3 ± 0.3) × 10² MPa and (1.0 ± 0.2) × 10² MPa, respectively. To show that the obtained bond parameters can be applied to a new specimen, a cross validation is conducted whereby parameters are fit from five specimens and then evaluated on the sixth. Additional validation of the obtained bond strength parameters was conducted with larger size artificial network simulations and peel tests. The proposed method in this work carries the dual advantages of characterizing actual bonds in a non-woven and characterizing hundreds of bonds simultaneously. The method can be applied to a variety of non-woven fabrics that are bonded at fiber-fiber intersections.     

基于XFEM-MBEM的嵌入式离散裂缝模型流固耦合数值模拟方法

离散缝网的表征与模拟是目前国内外研究的热点. 在非常规油气开发过程中, 由于地应力场的存在会对裂缝的流动属性产生显著影响, 若将裂缝视为静态对象, 与矿场数据会出现极大偏差, 因此要基于动态裂缝做更深入的研究. 本文针对致密油藏应力场?渗流场耦合力学问题, 提出了一种高效的混合数值离散化方法, 其中采用扩展有限元法 (XFEM) 求解岩石的弹性形变, 采用了混合边界元法 (MBEM) 精确计算基岩与裂缝间的非稳态窜流, 这两种数值格式是完全耦合的, 并对整体计算格式的时间项进行了全隐式求解, 可准确表征致密油藏开采过程中的裂缝变形及流体流动机理. 此外, 本文采用了嵌入式离散裂缝前处理算法显式表征大尺度水力压裂缝, 并考虑了支撑剂的作用; 采用了双孔有效应力原理和双重介质隐式裂缝表征方法, 可捕捉基质与小尺度天然裂缝的动态信息; 由此, 本文所提出的混合模型综合表征了基质?天然裂缝?水力压裂缝共同组成的致密油藏复杂渗流环境, 并通过几个实例论证了模型的准确性, 研究表明: 对致密油藏压裂水平井进行产能评价时, 应力场所引起渗流参数的改变及裂缝开度降低的影响不可忽略. 本文研究可为非常规油气资源的开发提供理论指导.      

基于损伤-虚拟张拉裂纹模型的地下爆炸围岩破坏规律研究

地下硐室作为爆炸危险物的隐蔽贮藏空间,有潜在的内爆炸风险。为研究内爆炸作用下硐室围岩的动态响应机制,提出了一种基于岩石HJC (Holmquist-Johnson-Cook)模型和节理内聚力单元的损伤-虚拟裂纹模型。分析了模拟方法的可靠性,并在此基础上,通过多物质ALE算法对球形硐室内爆炸过程进行数值模拟,分析了围岩损伤范围和分区破坏规律。研究表明:插入内聚力单元弥补了HJC模型无法模拟低静水压力下张拉破坏的不足,且尺寸效应易于处理。模拟方法同时考虑了岩体内张拉裂纹的扩展和岩石材料的塑性损伤,能够真实地反映岩石破坏的全过程。以红砂岩为例,根据数值模拟结果,填实（耦合装药）爆炸时围岩分区破坏规律明显,破碎区比例半径为0.26 m/kg1/3、裂隙区比例半径为0.47 m/kg1/3。随着硐室尺寸的增大,空气的间隔作用可以减小爆炸荷载对围岩的损伤作用,比例半径达到0.52 m/kg1/3时,可以实现爆炸荷载的完全解耦。     

基于连续-非连续单元方法的露天矿三维台阶爆破全过程数值模拟

爆破开采是露天矿采选总成本控制的首要环节，数值模拟是进行露天矿爆破开采优化设计及爆破效果分析的有效手段。利用连续-非连续单元方法（continuum-discontinuum element method，CDEM）对露天矿的三维台阶爆破过程进行了模拟，通过朗道爆炸模型实现了爆炸作用力的精确计算，通过弹性-损伤-断裂本构实现了岩体损伤破裂过程的描述，通过半弹簧-目标面及半棱-目标棱的联合接触算法实现了破碎岩块碰撞、飞散及堆积过程的高效模拟。开展了小尺度单自由面爆破过程的数值模拟，计算给出的块度分布曲线、爆破漏斗体积等参数与文献中模型实验的结果基本一致，证明了CDEM及本文所述各类模型在模拟爆炸破岩方面的精确性。以鞍千矿南采区的露天铁矿爆破开采为研究对象，建立了3排21炮孔的三维台阶爆破概化模型，模拟了从炸药起爆、岩体损伤破裂到最后爆堆形成的全过程；计算结果表明，除后缘拉裂槽外，数值计算给出的爆堆形态、顶部鼓起高度等与现场的测试结果基本一致，证明了利用CDEM开展三维露天台阶爆破全过程模拟的可行性。     

Bedding plane-emb e dde d augmente d virtual internal bonds for fracture propagation simulation in shale

To effectively simulate the fracture propagation in shale, the bedding plane(BP) effect is incorporated into the augmented virtual internal bond(AVIB) constitutive relation through BP tensor. Comparing the BP-embedded AVIB with the theory of transverse isotropy, it is found the approach can represent the anisotropic properties induced by parallel BPs. Through the simulation example, it is found that this method can simulate the stiffness anisotropy of shale and can represent the effect of BPs on hydraulic fracture propagation direction. Compared with the BP-embedded virtual internal bond(VIB), this method can account for the various Poisson's ratio. It provides a feasible approach to simulate the fracture propagation in shale.     

岩体裂隙扩展过程的数值模拟

本文根据断裂力学理论及节理岩体的等效连续模型,探讨了岩体裂隙扩展过程的数值模拟方法,文末算例说明了方法的有效性。     

下卧软弱岩带变形模量取值研究

软弱岩带的变形模量对于工程建设意义重大,当前变形试验方法难以表征原位条件下软弱岩带的变形模量,导致取值偏于保守。本文基于弹性力学推导了下卧软弱岩带变形模量的解析解,建立了一种新的软弱岩带变形模量的取值方法。在现场试验成果的基础上,采用有限元反演,数值模拟成果与试验计算成果较为吻合,证明该计算方法合理可行,较常规方法直接在软弱岩带上开展的试验成果大幅提高。     

A method for calculating the stress–strain state in the general boundary value problem of metal forming—part 2. Impact of a bar against a rigid obstacle

The method for calculating stress–strain state and fracture proposed by Kolmogorov, 1995and in Part 1 of this present paper is illustrated by the simple problem of a thin bar impacting a rigid obstacle. Known exact solutions are used to test the method. On the basis of the stability theory, the one-dimensional solution has been shown to be legitimate. Mathematical simulation of bar fragmentation resulting from impact has been carried out.     

岩石爆破分形损伤模型研究

在Ｋｕｓ和Ｔｈｏｒｎｅ等人研究基础上，利用分形理论及应用成果提出了一种新的岩石爆破损伤模型。新模型将岩体中各种结构弱面视为初始损伤，并将其分形维数作为岩石性质的主要参量，用损伤能量耗散率及与分形维数的关系表征爆破过程中的岩石损伤演化规律，从而弥补了现有模型未考虑初始损伤影响和将损伤演化归结为体积应变函数的不足。还通过模拟试验等手段验证了新模型的如上观点。