基于Voronoi建模和Cohesive单元的PBX细观力学模型

针对PBX的细观力学模型存在模型结构与真实PBX细观形貌不符、破坏行为描述不全面(仅考虑界面脱粘)、研究对象单一(大多针对PBX-9501)等问题,建立了基于Voronoi建模和Cohesive单元的PBX细观力学模型。首先,基于光学显微镜下的细观结构图,采用Voronoi方法构建了更贴近实际形貌的PBX炸药细观模型;然后,采用内聚力有限元方法在炸药晶体内部、黏结剂内部及炸药颗粒/黏结剂界面处都引入了Cohesive单元,以在二维尺度上实现任意路径破坏行为的模拟;最后,针对某HMX基-F2311黏结剂炸药和某HMX基-F2313黏结剂炸药标定了模型参数,模拟了这两种炸药的准静态拉伸行为。结果表明,两种炸药数值模拟的拉伸变形曲线、破坏应力、破坏应变都与试验数据吻合良好,表明该模型不仅适用于不同PBX炸药材料(线弹性PBX、非线性PBX)的变形行为描述,还可以分析PBX炸药在拉伸破坏过程中的细观机理,刻画不同PBX炸药的裂纹萌生、扩展和贯穿的过程。

Finite Element Simulations of the Localized Failure and Fracture Propagation in Cohesive Materials with Friction

Strain localization frequently occurs in cohesive materials with friction(e.g.,composites,soils,rocks)and is widely recognized as a fundamental cause of progressive structural failure.Nonetheless,achieving high-fidelity simulation for this issue,particularly concerning strong discontinuities and tension-compression-shear behaviors within localized zones,remains significantly constrained.In response,this study introduces an integrated algorithmwithin the finite element framework,merging a coupled cohesive zone model(CZM)with the nonlinear augmented finite elementmethod(N-AFEM).The coupledCZMcomprehensively describes tension-compression and compressionshear failure behaviors in cohesive,frictional materials,while the N-AFEM allows nonlinear coupled intraelement discontinuities without necessitating extra nodes or nodal DoFs.Following CZM validation using existing experimental data,this integrated algorithm was utilized to analyze soil slope failure mechanisms involving a specific tensile strength and to assess the impact of mechanical parameters(e.g.,tensile strength,weighting factor,modulus)in soils.

Extended finite element-based cohesive zone method for modeling simultaneous hydraulic fracture height growth in layered reservoirs

In this study,a fully coupled hydromechanical model within the extended finite element method(XFEM)-based cohesive zone method(CZM)is employed to investigate the simultaneous height growth behavior of multi-cluster hydraulic fractures in layered porous reservoirs with modulus contrast.The coupled hydromechanical model is first verified against an analytical solution and a laboratory experiment.Then,the fracture geometry(e.g.height,aperture,and area)and fluid pressure evolutions of multiple hydraulic fractures placed in a porous reservoir interbedded with alternating stiff and soft layers are investigated using the model.The stress and pore pressure distributions within the layered reservoir during fluid injection are also presented.The simulation results reveal that stress umbrellas are easily to form among multiple hydraulic fractures’tips when propagating in soft layers,which impedes the simultaneous height growth.It is also observed that the impediment effect of soft layer is much more significant in the fractures suppressed by the preferential growth of adjoining fractures.After that,the combined effect of in situ stress ratio and fracturing spacing on the multi-fracture height growth is presented,and the results elucidate the influence of in situ stress ratio on the height growth behavior depending on the fracture spacing.Finally,it is found that the inclusion of soft layers changes the aperture distribution of outmost and interior hydraulic fractures.The results obtained from this study may provide some insights on the understanding of hydraulic fracture height containment observed in filed.

Investigation of FRP and SFRC Technologies for Efficient Tunnel Reinforcement Using the Cohesive Zone Model

Amid urbanization and the continuous expansion of transportation networks,the necessity for tunnel construction and maintenance has become paramount.Addressing this need requires the investigation of efficient,economical,and robust tunnel reinforcement techniques.This paper explores fiber reinforced polymer(FRP)and steel fiber reinforced concrete(SFRC)technologies,which have emerged as viable solutions for enhancing tunnel structures.FRP is celebrated for its lightweight and high-strength attributes,effectively augmenting load-bearing capacity and seismic resistance,while SFRC’s notable crack resistance and longevity potentially enhance the performance of tunnel segments.Nonetheless,current research predominantly focuses on experimental analysis,lacking comprehensive theoretical models.To bridge this gap,the cohesive zone model(CZM),which utilizes cohesive elements to characterize the potential fracture surfaces of concrete/SFRC,the rebar-concrete interface,and the FRP-concrete interface,was employed.A modeling approach was subsequently proposed to construct a tunnel segment model reinforced with either SFRC or FRP.Moreover,the corresponding mixed-mode constitutive models,considering interfacial friction,were integrated into the proposed model.Experimental validation and numerical simulations corroborated the accuracy of the proposed model.Additionally,this study examined the reinforcement design of tunnel segments.Through a numerical evaluation,the effectiveness of innovative reinforcement schemes,such as substituting concrete with SFRC and externally bonding FRP sheets,was assessed utilizing a case study from the Fuzhou Metro Shield Tunnel Construction Project.

Virtual strain loading method for low temperature cohesive failure of asphalt binder

Cohesive failure is one of the primary reasons for low-temperature cracking in asphalt pavements.Understanding the micro-level mechanism is crucial for comprehending cohesive failure behavior.However,previous literature has not fully reported on this aspect.Moreover,there has been insufficient attention given to the correlation between macroscopic and microscopic failures.To address these issues,this study employed molecular dynamics simulation to investigate the low-temperature tensile behavior of asphalt binder.By applying virtual strain,the separation work during asphalt binder tensile failure was calculated.Additionally,a correlation between macroscopic and microscopic tensile behaviors was established.Specifically,a quadrilateral asphalt binder model was generated based on SARA fractions.By applying various combinations of virtual strain loading,the separation work at tensile failure was determined.Furthermore,the impact of strain loading combinations on separation work was analyzed.Normalization was employed to establish the correlation between macroscopic and microscopic tensile behaviors.The results indicated that thermodynamic and classical mechanical indicators validated the reliability of the tetragonal asphalt binder model.The strain loading combination consists of strain rate and loading number.All strain loading combinations exhibited the similar tensile failure characteristic.The critical separation strain was hardly influenced by strain loading combination.However,increasing strain rate significantly enhanced both the maximum traction stress and separation work of the asphalt binder.An increment in the loading number led to a decrease in separation work.The virtual strain combination of 0.5%-80 provided a more accurate representation of the actual asphalt's tensile behavior trend.

基于cohesive单元海域天然气水合物储层水力压裂模拟

天然气水合物作为一种高效清洁能源,广泛分布于我国南海海域的沉积地层中。我国先后于2017年和2020年成功开展了2次试开采,但由于海域天然气水合物特殊的赋存条件,单井水合物试采仍然面临着开采范围小、高产稳产时间短等问题。为了提高水合物的开采范围,基于cohesive单元进行了水合物储层二维水力压裂数值模型研究,比较了100 m×100 m和20 m×20 m两种模型的裂缝半长和宽度,得出了当注入压力为25 MPa时,压裂裂缝半长均为6 m,最大宽度分别为5.8、5.5 mm,构建尺寸较大的模型得出的实验结果更加准确。并且研究了裂缝宽度随注入时间的变化规律,随着注入压力和注入量的不断增加,初期裂缝宽度急速变大,后续在地应力和注入流体压力的共同作用下裂缝出现“阶梯式”的扩展规律。该研究在页岩气和煤层气等非常规能源储层水力压裂模型分析中得到了成功运用,为海域天然气水合物储层水力压裂提供一定的理论指导。

Separation work analysis of cohesive law and consistently coupled cohesive law

An appropriate coupled cohesive law for predicting the mixed mode failure is established by combining normal separation and tangential separation of surfaces in the cohesive zone model （CZM） and the cohesive element method. The Xu-Needleman exponential cohesive law with the fully shear failure mechanism is one of the most popular models. Based on the proposed consistently coupled rule/principle, the Xu-Needleman law with the fully shear failure mechanism is proved to be a non-consistently coupled cohesive law by analyzing the surface separation work. It is shown that the Xu-Needleman law is only valid in the mixed mode fracture when the normal separation work equals the tangential separation work. Based on the consistently coupled principle and the modification of the Xu-Needleman law, a consistently coupled cohesive （CCC） law is given. It is shown that the proposed CCC law has already overcome the non-consistency defect of the Xu-Needleman law with great promise in mixed mode analyses.

Numerical simulation of the mechanical behavior of superconducting tape in conductor on round core cable using the cohesive zone model

Cables composed of rare-earth barium copper oxide(REBCO)tapes have been extensively used in various superconducting devices.In recent years,conductor on round core(CORC)cable has drawn the attention of researchers with its outstanding current-carrying capacity and mechanical properties.The REBCO tapes are wound spirally on the surface of CORC cable.Under extreme loadings,the REBCO tapes with layered composite structures are vulnerable,which can lead to degradation of critical current and even quenching of superconducting devices.In this paper,we simulate the deformation of CORC cable under external loads,and analyze the damage inside the tape with the cohesive zone model(CZM).Firstly,the fabrication and cabling of CORC are simulated,and the stresses and strains generated in the tape are extracted as the initial condition of the next step.Then,the tension and bending loads are applied to CORC cable,and the damage distribution inside the tape is presented.In addition,the effects of some parameters on the damage are discussed during the bending simulations.

基于XFEM与Cohesive模型分析PBX裂纹产生与扩展

利用扩展有限元法(XFEM)分析PBX-9502带孔板状试件在整体压缩下由局部裂纹萌生到裂纹扩展全过程的开裂破坏机理。采用应力状态相关的强度面、非关联流动法则及Cohesive模型,描述了材料在复杂应力状态下的非线性本构行为以及材料的破坏行为。进行了数值模拟结果与美国洛斯阿拉莫斯国家实验室(LANL)试验结果的对比。结果表明,含孔洞的平板在整体压应力环境下孔洞周围产生局部拉伸应力,这种拉伸条件导致局部裂纹萌生。数值模拟的裂纹发展趋势与试验结果相吻合,包括裂纹时程的整体走势和拐点、启裂时刻、裂纹初期扩展速度等。基于扩展有限元方法和内聚模型法,可模拟高聚物粘结炸药(PBX)含能材料的裂纹萌生、扩展。

基于Cohesive单元法的弃置井水泥塞-套管界面胶结失效数值模拟

水泥塞-套管界面胶结失效对弃置井井筒完整性提出了严峻挑战。考虑地下流体与水泥塞的流固耦合作用,基于cohesive单元方法,建立水泥塞-套管-水泥环-地层系统三维有限元模型,模拟垂直井水泥塞-套管界面裂缝剥离过程,研究地应力对界面裂缝损伤演化的影响,并分析水平地应力、水泥塞力学参数及界面性质对裂缝剥离高度的影响。结果显示:水平地应力均匀时剥离裂缝沿着整个界面圆周延伸且高度相等,水平地应力非均匀时剥离裂缝倾向于沿着界面某一圆周角扩展且在最大水平地应力方向具有较大的高度;水泥塞弹性模量从30 GPa减小到1 GPa,裂缝剥离高度降低9.3 m,临界法向强度从0.25 MPa增大到2.0 MPa,裂缝剥离高度降低6.5 m,表明较低的弹性模量及较大的临界法向强度有利于减小水泥塞-套管界面胶结失效的风险;水泥塞泊松比从0.35减小到0.10,裂缝剥离高度仅降低2.0 m,临界剪切强度从0.5 MPa增加到4.0 MPa,裂缝剥离高度仅降低3.3 m,表明泊松比和临界剪切强度对界面胶结失效影响较小。建立的模型能够为水泥浆配方优选和井筒弃置工艺优化提供指导。

Numerical study of fatigue damage of asphalt concrete using cohesive zone model

In order to investigate the fatigue behavior of asphalt concrete, a new numerical approach based on a bi-linear cohesive zone model （CZM） is developed. Integrated with the CZM, a fatigue damage evolution model is established to indicate the gradual degradation of cohesive properties of asphalt concrete under cyclic loading. Then the model is implemented in the finite element software ABAQUS through a user-defined subroutine. Based on the proposed model, an indirect tensile fatigue test is finally simulated. The fatigue lives obtained through numerical analysis show good agreement with laboratory results. Fatigue damage accumulates in a nonlinear manner during the cyclic loading process and damage initiation phase is the major part of fatigue failure. As the stress ratio increases, the time of the steady damage growth stage decreases significantly. It is found that the proposed fatigue damage evolution model can serve as an accurate and efficient tool for the prediction of fatigue damage of asphalt concrete.

基于Cohesive单元的石拱桥主拱圈二相数值模拟方法

为解决石砌体材料非均质的描述问题,提出一种基于Cohesive单元的石拱桥主拱圈两相数值模拟方法。视石砌体为两相材料(砌块和砌缝),采用实体单元模拟砌块并引入非线性本构描述其破坏行为,在相邻砌块间插入Cohesive单元考虑砌缝砂浆的剪切和拉伸破坏。通过室内试验与数值模拟对比验证方法的有效性及适用性,分析了砌缝抗剪摩擦系数μ、加载位置等敏感参数对拱桥承载力的影响。结果表明:基于Cohesive单元的石砌体两相数值模型,可以有效描述石砌体材料的非均匀性及石拱桥的破坏过程(尤其是砌缝剪切滑移破坏行为),可为石拱桥极限承载力评估提供重要信息,如荷载-位移曲线、破坏模式等。此外,研究结果还发现主拱圈破坏机制由拱的受弯、受剪特性决定,并与砌缝抗剪摩擦系数μ强相关。

Delamination analysis of woven fabrication laminates using cohesive zone model

A new test method was proposed to evaluate the cohesive strength of composite laminates. Cohesive strength and the critical strain energy for Mode-II interlamiar fracture of E-glass/epoxy woven fabrication were determined from the single lap joint(SLJ) and end notch flexure(ENF) test, respectively. In order to verify their adequacy, a cohesive zone model simulation based on interface finite elements was performed. A closed form solution for determination of the penalty stiffness parameter was proposed. Modified form of Park-Paulino-Roesler traction-separation law was provided and conducted altogether with trapezoidal and bilinear mixed-mode damage models to simulate damage using Abaqus cohesive elements. It was observed that accurate damage prediction and numerical convergence were obtained using the proposed penalty stiffness. Comparison between three damage models reveals that good simulation of fracture process zone and delamination prediction were obtained using the modified PPR model as damage model. Cohesive zone length as a material property was determined. To ensure the sufficient dissipation of energy, it was recommended that at least 4 elements should span cohesive zone length.

Survey Research into Cohesive Device Misuses in Business English Writing

Cohesive devices in students’business English writing are regarded as the object of the research.Based on Haliday and Hasan’s cohesion theory,this paper introduces commonly-used cohesive devices in English writing.With the method of quantitative data,use and misuse frequency of cohesive devices in students’writings can be known.The paper will also analyze why misuses happen through qualitative data analysis and explore right ways of using cohesive devices.

Analysis on the Cohesive and Coherent Features of John Kennedy's First Inaugural Address

John Kennedy's first inaugural address is one of the widely appreciated speeches worldwide. It is famous not only for calling up the American people to well serve the country, but also for its extraordinary linguistic power to arouse the listeners' emotions, which lies to a great extent in the marvelous employment of the cohesive and coherent devices in the process of its delivery. Cohesion and coherence are two elementary and significant concepts in the theoretical system of discourse analysis. There-fore, they play an important role in the structuring, arrangement, interpretation and analysis of a discourse. In this sense, it is significant to analyze the cohesive and coherent features of John Kennedy's first inaugural address in order to obtain a penetrating comprehension of the speech in many aspects. A detailed analysis on the cohesive and coherent features of the speech has been conducted in this paper. In the aspect of cohesion in the address, the devices employed fall into two categories: structural cohesion and non-structural cohesion. Structural cohesive devices used in the discourse are mainly grammatical cohesion and lexical cohesion like repetition, ellipsis, conjunction, etc. Non-structural methods adopted in the speech are transitivity, mood and modality, thematic progression, parallel structure and so on. In the aspect of coherence, five levels of coherent methods have been employed, namely, lexical level, syntax level, semantic level, phonological level and social semiotic level. The neat intermingling of the cohesive and coherent methods function cooperatively and lead to the smooth going of the text.

界面应力传递重新分析及Cohesive模型参数的确定

在经典剪滞理论中引入双线性cohesive模型表征纤维/基体之间的非理想界面,重新分析了纤维增强复合材料中的应力传递机理,得到了考虑界面因素的应力分布。用上述结果解释了单丝段裂实验过程中的现象,讨论了界面参数和材料性能对应力分布的影响。基于上述理论,建立了用cohesive单元表征界面的模拟单丝段裂实验的三维有限元模型,结合单丝段裂实验结果,提出了一种估测cohesive界面刚度参数的新方法。数值和理论分析结果与实验结果对比,吻合良好,可以为材料的界面性能分析和材料设计提供参考依据。

COHESIVE ZONE FINITE ELEMENT-BASED MODELING OF HYDRAULIC FRACTURES

Hydraulic fracturing is a powerful technology used to stimulate fluid production from reservoirs. The fully 3-D numerical simulation of the hydraulic fracturing process is of great importance to the efficient application of this technology, but is also a great challenge because of the strong nonlinear coupling between the viscous flow of fluid and fracture propagation. By taking advantage of a cohesive zone method to simulate the fracture process, a finite element model based on the existing pore pressure cohesive finite elements has been established to investigate the propagation of a penny-shaped hydraulic fracture in an infinite elastic medium. The effect of cohesive material parameters and fluid viscosity on the hydraulic fracture behaviour has been investigated. Excellent agreement between the finite element results and analytical solutions for the limiting case where the fracture process is dominated by rock fracture toughness demonstrates the ability of the cohesive zone finite element model in simulating the hydraulic fracture growth for this case.

Rheological Properties and Incipient Motion of Cohesive Sediment in the Haine Estuary of China

The Haihe cohesive sediment, which is typical in China, is studied systematically for its basic physical and incipient motion properties. Following the requirements of dredging works in the Haihe Estuary, cohesive sediment samples were taken from three locations. Laboratory experiments were conducted to determine the rheological properties of these samples and to examine the incipient motion of the cohesive sediment. It is found that the cohesive sediment has an obvious yield stress tau(b), which increases with the mud density in a manner of an exponential function, and so does the viscosity parameter eta. The cohesive sediment behaves like a Bingham fluid when its density is below 1.38 similar to 1.40 g/cm(3), and when denser than these values, it may become a power-law fluid. The incipient motion experiment also revealed that the incipient velocity of the cohesive sediment increases with die density in an exponential manner. Therefore, the incipient motion is primarily related to the density, which is different from the case for non-cohesive sediment in which the incipient motion is con-elated with the diameter of sand particles instead. The incipient motion occurs in two different ways depending on the concentration of mud in the bottom. For sufficiently fine particles and a concentration lower than 1.20 g/cm(3), the cohesive sediment appears as fluidized mud, and the incipient motion is in the form of instability of an internal wave. For a higher concentration, the cohesive sediment appears as general quasi-solid-mud, and the incipient motion can be described by a series of extended Shields curves each with a different porosity for newly deposited alluvial mud.

Face stability analysis for a longitudinally inclined tunnel in anisotropic cohesive soils

A stability analysis approach of tunnel face considering a longitudinally inclined tunnel angle and anisotropic purely cohesive soils based on a continuous velocity field (CVF) is investigated in this study. Based on the kinematic approach of limit analysis and the discretization technique, an improved three-dimensional CVF model for longitudinally inclined tunnels driven by pressurized shields is proposed. With the proposed model, the critical support pressure acted on tunnel face is determined by the work-balance equation. A serial of finite element numerical models are conducted to validate the proposed model. Finally, the effects of tunnel inclination angles, several dimensionless parameters as well as soil anisotropy on the critical support pressure are investigated. The numerical results show that the effects of the soil anisotropy and the tunnel inclination angle on tunnel faces should be considered in the actual design of tunneling engineering.

Maximum Force of Inclined Pullout of A Torpedo Anchor in Cohesive Beds

Torpedo anchors have been used in mooring systems for deep-water oil and gas projects owing to their prominent advantages, such as low cost and easy installation. The maximum force of torpedo anchors is crucial not only to the safety and stability of vessels and other marine facilities, but also for an economical design. It is necessary to develop reliable formula for fast predicting their maximum inclined force of a torpedo anchor in cohesive beds. In this study, the maximum inclined force of a torpedo anchor vertically embedded in cohesive beds was extensively investigated. 316 sets of inclined pullout laboratory tests were carried out for 9 differently shaped torpedo anchors which were vertically embedded in different cohesive beds. The loading curves were automatically acquisitioned and their characteristics were analyzed. The load angle relative to the horizontal varied from 20° to 90°. A new formula for fast calculating the maximum inclined force of the torpedo anchor vertically embedded in cohesive beds was obtained based on force analysis and a nonlinear regression on the data from the present and other studies. Effect aspects on the tests are discussed and further studies are highlighted.