Nonlinear phase field theory for fracture and twinning with analysis of simple shear

A phase field theory for coupled twinning and fracture in single crystal domains is developed. Distinct order parameters denote twinned and fractured domains, finite strains are addressed and elastic nonlinearity is included via a neo-Hookean strain energy potential. The governing equations and boundary conditions are derived; an incremental energy minimization approach is advocated for prediction of equilibrium microstructural morphologies under quasi-static loading protocols. Aspects of the theory are analysed in detail for a material element undergoing simple shear deformation. Exact analytical and/or one-dimensional numerical solutions are obtained in dimensionless form for stress states, stability criteria and order parameter profiles at localized fractures or twinning zones. For sufficient applied strain, the relative likelihood of localized twinning vs. localized fracture is found to depend only on the ratio of twin boundary surface energy to fracture surface energy. Predicted criteria for shear stress-driven fracture or twinning are often found to be in closer agreement with test data for several types of real crystals than those based on the concept of theoretical strength.     

XPS study of Cr segregation in a martensitic stainless steel

Cr martensitic steels are promising materials for structural applications in future nuclear fusion reactors. Because the embrittlement after tempering treatments can be a serious problem, the fracture mode of a steel with 10.5 wt% of Cr treated at 700°C for 18 h has been investigated through Charpy tests in the temperature range from −100°C to +150°C. X-ray photoelectron spectroscopy (XPS) analyses carried out on the fracture surfaces evidenced the segregation of Cr in both ductile and brittle (quasicleavage) fields. The unexpected result indicates that Cr segregation weakens the atomic bonds; thus, the fracture path in both the cases corresponds to the zones with higher Cr content.     

Dynamic fracture testing of polymethyl-methacrylate (PMMA) single-edge notched beam

Dynamic fracture in single-edge notched polymethyl-methacrylate (PMMA) beams have been investigated by three-point-bending impact testing with a drop-weight machine. A high-speed camera combined with the digital image correlation (DIC) method is used to capture the impact-induced crack initiation and propagation, as well as the beam deformation fields and the open mode strain at the original notch tip. The crack propagation length is recorded and the instantaneous crack velocity is calculated. Furthermore, the dynamic fracture toughness K_Id is quantified from the loading-displacement relations at different impact velocities. The effects of the impact velocity and impact energy on dynamic fracture toughness, fracture initiation strain, as well as the corresponding influences on the fracture propagation velocity, are discussed.     

含中心裂纹梯度纳晶金属的数值模拟

梯度纳晶金属由于其微观组织的梯度分布,力学属性也呈现梯度变化,这使得其表现出不同于传统均匀材料的断裂行为．利用材料力学参数的梯度分布来表征梯度纳晶金属中晶粒尺寸的梯度变化,并编写ABAQUS和MATLAB脚本程序建立分层有限元模型．通过数值模拟计算了含有初始中心裂纹的梯度纳晶金属在受远端均匀拉应力作用下的裂尖J积分,分别研究了屈服应力梯度、裂纹角度和裂纹长度对金属材料断裂韧性的影响,并与传统粗晶进行了对比．结果表明梯度纳米结构的存在导致梯度纳晶金属内部的中心裂纹两端表现出不同的断裂韧性,小晶粒一侧裂尖的抗裂韧性优于大晶粒一侧裂尖,且屈服应力梯度绝对值越大,两者差距越大．梯度纳晶金属的断裂韧性受中心裂纹角度和长度变化的影响与传统粗晶金属基本一致,同时在晶粒尺寸梯度的作用下梯度纳晶的裂尖J积分略低于粗晶,即整体上拥有更好的抗裂韧性．     

高韧性合金钢断裂韧性及疲劳裂纹扩展门槛值的联合试验方法

在柴油机曲轴、连杆等关键零部件的可靠性设计和失效评估中,断裂韧性及疲劳裂纹扩展门槛值分别是衡量材料抵抗裂纹失稳扩展和裂纹开始扩展的重要指标．但是,对于高韧性合金材料,难以通过常规试验所推荐的厚度确定平面应变断裂韧性,而门槛值的测定通常不但非常耗时,且难以直接应用于不同循环特性的实际结构．本文针对高韧性合金钢34CrNi3MoA,提出一种将断裂韧性和疲劳裂纹扩展门槛值试验合二为一的试验方法,即用同一个试件可以同时测定门槛值和断裂韧性．利用断裂韧性关于试件厚度的渐近特性,以几种较薄试件的试验,确定平面应变状态下的断裂韧性．试验结果还表明,裂纹扩展门槛值的试件厚度依存性可以忽略,并给出了任意循环特性（应力比）下的门槛值计算公式．     

Simulation of proppant transport at intersection of hydraulic fracture and natural fracture of wellbores using CFD-DEM

Proppants transport is an advanced technique to improve the hydraulic fracture phenomenon, in order to promote the versatility of gas/oil reservoirs. A numerical simulation of proppants transport at both hydraulic fracture (HF) and natural fracture (NF) intersection is performed to provide a better understanding of key factors which cause, or contribute to proppants transport in HF–NF intersection. Computational fluid dynamics (CFD) in association with discrete element method (DEM) is used to model the complex interactions between proppant particles, host fluid medium and fractured walls. The effect of non-spherical geometry of particles is considered in this model, using the multi-sphere method. All interaction forces between fluid flow and particles are considered in the computational model. Moreover, the interactions of particle–particle and particle–wall are taken into account via Hertz–Mindlin model. The results of the CFD-DEM simulations are compared to the experimental data. It is found that the CFD-DEM simulation is capable of predicting proppant transport and deposition quality at intersections which are in agreement with experimental data. The results indicate that the HF–NF intersection type, fluid velocity and NF aperture affect the quality of blockage occurrence, presenting a new index, called the blockage coefficient which indicates the severity of the blockage.     

基于FDM-DEM耦合的冲击损伤大理岩静态断裂力学特征研究

为探究循环冲击损伤后大理岩的静态断裂力学特征,基于有限差分（finite difference method,FDM）-离散元（discrete element method,DEM）耦合的建模技术构建了三维分离式霍普金森压杆（split Hopkinson pressure bar,SHPB）数值模型,其中杆件系统和岩石试件分别采用FLAC3D和PFC3D程序建模。利用该模型对中心直切槽半圆盘（NSCB）试样进行了恒定子弹速度下的循环冲击,随后对受损试样进行静态三点弯曲断裂实验。通过编写Fish程序,提取试样断裂面数据,对断裂面进行重构并定量计算表面粗糙度。通过与相关室内实验结果的对比分析,验证了本文数值分析的合理性与可靠性。模拟结果表明,随着循环冲击次数的增加,试样内部微裂纹、破碎颗粒均增加。连接力场分布混乱,部分力链发生断裂。力链的变化是试样力学性能劣化的根本原因。在静态三点弯曲断裂实验中,冲击5次后试样的静态断裂韧度较天然试样产生一定程度的降低。试样在静载过程中产生的微裂纹和碎块的数量随循环冲击次数的增加而增加,断裂面粗糙度随循环冲击次数的增加而增加。     

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

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

直通型CT试样COD弹塑性换算

基于刚性转动假设和功等效原理,推导了直通型CT试样的裂纹嘴张开位移和加载线张开位移弹塑性换算公式,并采用Cr2Ni2MoV钢对公式进行试验验证.根据弹塑性有限元分析,刚性转动假设得到了很好的验证,刚性转动中心靠近裂尖位置且受材料本构关系的影响甚微,转动半径R仅与裂纹长度和试样宽度之比a/W有关.有限元分析和试验结果表明,直通型CT试样的COD弹塑性换算公式用于塑性功计算具有更好的精度.     

X100直缝埋弧焊管断裂韧性的裂纹尖端张开位移试验研究

为了检验X100直缝埋弧焊管的断裂韧性,采用三点弯曲多试样法,在室温、-10℃和-30℃下对该焊管的母材、热影响区和焊缝部位进行了裂纹尖端张开位移(CTOD)试验.测试结果表明:在同一温度下,CTOD试验的表观启裂值δ0.05、条件启裂值δi和最大载荷值δm的分布规律为母材最好,焊缝次之,热影响区最差;采用裂纹扩展阻力和撕裂模量表征的分布规律为母材最好,热影响区次之,焊缝最差;焊接热影响区组织中含有大量呈块状分布的铁素体,而粒状贝氏体等增强相的含量较少,这会造成其韧性下降,使焊管出现局部失稳现象.