A micromechanical model for predicting the fracture toughness of functionally graded foams

In this paper, finite element method based micromechanical analysis is used to understand the fracture behavior of functionally graded foams. The finite element analysis uses a micro-mechanical model in conjunction with a macro-mechanical model in order to relate the stress intensity factor to the stresses in the struts of the foam. The stress intensity factor at the crack tip of the macro-mechanical model can be evaluated using either the J-contour integral or the stresses in the singularity-dominated zone. The fracture toughness is evaluated for various crack positions and length within the functionally graded foam. Then the relationship between the fracture toughness of the graded foam and the local density at the crack tip is studied. Convergence tests for both macro-mechanical and micro-mechanical model analysis were conducted in order to maintain adequate accuracy with reasonable computational time. Fracture toughness of homogenous foams and functionally graded foams for various cases are presented as a function of relative density. This study indicates that the fracture toughness of functionally graded foams mainly depends on the relative density at the crack-tip.     

Analytical representation of the non-square-root singular stress field at a finite angle sharp notch

The stress field near the tip of a finite angle sharp notch is singular. However, unlike a crack, the order of the singularity at the notch tip is less than one-half. Under tensile loading, such a singularity is characterized by a generalized stress intensity factor which is analogous to the mode I stress intensity factor used in fracture mechanics, but which has order less than one-half. By using a cohesive zone model for a notional crack emanating from the notch tip, we relate the critical value of the generalized stress intensity factor to the fracture toughness. The results show that this relation depends not only on the notch angle, but also on the maximum stress of the cohesive zone model. As expected the dependence on that maximum stress vanishes as the notch angle approaches zero. The results of this analysis compare very well with a numerical (finite element) analysis in the literature. For mixed-mode loading the limits of applicability of using a mode I failure criterion are explored.     

求解I型裂纹构元J积分的半解析方法

表征裂纹尖端应力应变场程度的J积分是一个定义明确、理论严密的弹塑性断裂力学基础参量. 目前J积分的计算主要是依靠塑性因子法和有限元法,但对各类裂纹构元获得J积分以及载荷-位移关系的解析公式以实现材料断裂韧性理论预测和材料测试是断裂力学的重要和困难的任务. 以J积分为参量的材料断裂测试中应用最广的是I型裂纹试样的断裂韧性测试. 本文在平面应变条件下,针对断裂韧性测试中使用的6种I型裂纹构元,基于能量等效假设,提出了J积分-载荷和载荷-位移的工程半解析统一表征方法,进而结合有限元分析的少量计算获得J积分-载荷和载荷-位移关系的半解析公式待定参数. 分析表明,6种I型裂纹构元的J积分-载荷和载荷-位移统一公式的预测结果与有限元结果吻合良好. 新提出的J积分-载荷工程半解析公式包含了材料的弹性模量、应力强度系数和应变硬化指数,能够广泛适应不同的材料,且运用该公式能够方便获取任意载荷点对应的J积分值. 应用新方法可便于获得各类I型裂纹构元的J积分-载荷和载荷-位移工程半解析公式.      

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

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

随机有限元方法在断裂分析中的应用

在幂律非线性随机有限元基础上,以单边裂纹板为例给出计算含量钢继裂参数,J（J积分）,δ（裂纹张开位移）,Δ（由裂纹引起的裂纹板上下底面相对位移）,θ（由裂纹引起的裂纹板上下底在相对转角）及其对基本随机变量变化率的方法和分析算例。     

基于临界距离点法的混凝土Ⅰ型断裂韧度的预测

混凝土的断裂韧度是极为重要的断裂力学材料参数，本文利用临界距离理论(TCD)中的点法和混凝土标准断裂梁试样裂缝尖端应力场的近似解，提出了一种对混凝土断裂问题按临界距离的点法进行分析的方法．基于已有文献的试验材料及其实验结果，分析了这种断裂研究分析方法的可行性．将临界距离点法与双K断裂准则应用于几何相似的混凝土梁的断裂韧性和断裂过程区分析，验证了TCD点法的可靠性．实验表明：临界距离理论的点法能够获得比双K断裂准则相对安全的失稳断裂韧度，利用临界距离法还可得到相应断裂过程区长度的极限估算值．     

The analysis of crack problems with non-local elasticity

IntroductionTheclassicalconhnuummechanicshasbeenusedtosolvemanyproblemsinmacrofracturemechanics,butencountersdifficulheswhentheeffectofITilcrocharacteristicdimensionshouldbetakenintoaccount.Thestressfieldverynearthecracktipisstillnotclear.Somephenomenaofshortcrackscannotbeexplained["']andsomemechanismoffracturehasnotbeensolvedyet.Thenon-localelashcitytheoryseemsattractivetotheseproblems.Thetheoryofnon-localelasticity,establishedanddevelopedbyEringenetal[3),connectstheclassicalcontinuummechan…     

40Cr材料动态起裂韧性KId（）的实验测试

描述了利用Ｈｏｐｋｉｎｓｏｎ压杆技术加载三点弯曲试样测试４０Ｃｒ，材料动态起裂韧性ＫＩｄ（）的试验方法。试样上的动态载荷历程由Ｈｏｐｋｉｎｓｏｎ杆直接测得，并分别代入动态有限元程序及近似公式求得动态应力强度因子历史；由贴在试样裂尖附近的应变片确定起裂时间，最终确定起裂时的动态应力强度因子值，即动态起裂韧性ＫＩｄ（）。试验结果表明:利用Ｈｏｐｋｉｎｓｏｎ压杆技术加载三点弯曲试样测试材料动态起裂韧性的方法是可行的，起裂时，动态有限元的位移法、应力法及近似公式法求得的动态应力强度因子值比较吻合；在本文的载荷速率下，４０Ｃｒ材料动态起裂韧性ＫＩｄ（）与准静态裂韧性ＫＩｄ（）相比，降低了约２８％。     

The damage tolerance of elastic-brittle, two-dimensional isotropic lattices

The fracture toughness of elastic-brittle 2D lattices is determined by the finite element method for three isotropic periodic topologies: the regular hexagonal honeycomb, the Kagome lattice and the regular triangular honeycomb. The dependence of mode I and mode II fracture toughness upon relative density is determined for each lattice, and the fracture envelope is obtained in combined mode I-mode II stress intensity factor space. Analytical estimates are also made for the dependence of mode I and mode II toughness upon relative density. The high nodal connectivity of the triangular grid ensures that it deforms predominantly by stretching of the constituent bars, while the hexagonal honeycomb deforms by bar bending. The Kagome microstructure deforms by bar stretching remote from the crack tip, and by a combination of bar bending and bar stretching within a characteristic elastic deformation zone near the crack tip. This elastic zone reduces the stress concentration at the crack tip in the Kagome lattice and leads to an elevated macroscopic toughness.Predictions are given for the tensile and shear strengths of a centre-cracked panel with microstructure given explicitly by each of the three topologies. The hexagonal and triangular honeycombs are flaw-sensitive, with a strength adequately predicted by linear elastic fracture mechanics (LEFM) for cracks spanning more than a few cells. In contrast, the Kagome microstructure is damage tolerant, and for cracks shorter than a transition length its tensile strength and shear strength are independent of crack length but are somewhat below the unnotched strength. At crack lengths exceeding the transition value, the strength decreases with increasing crack length in accordance with the LEFM estimate. This transition crack length scales with the parameter of bar length divided by relative density of the Kagome grid, and can be an order of magnitude greater than the cell size at low relative densities. Finally, the presence of a boundary layer is noted at the free edge of a crack-free Kagome grid loaded in tension and in shear. Deformation within this boundary layer is by a combination of bar bending and stretching whereas remote from the free edge the Kagome grid deforms by bar stretching (with a negligible contribution from bar bending). The edge boundary layer degrades both the macroscopic stiffness and strength of the Kagome plate. No such boundary layer is evident for the hexagonal and triangular honeycombs.     

功能梯度材料的黏弹性断裂问题

功能梯度材料(FGM)是一种不同于传统复合材料的新型工程复合材料 [1], 国内外关于FGM 的断裂力学方面的研究发展非常迅速. 关于FGM静态裂纹问题，学者们研究了不同类型裂纹 尖端场的应力强度因子 [2-5], 探讨了有限长裂纹在不用载荷作用下的传播等问题. 而关于动 态裂纹问题，也已经取得很大成就 [6-9]. FGM一个很重要的应用是高温结构材料，在强大的 热环境中，很多材料都呈现出黏弹性. 因此，研究FGM的黏弹性断裂力学非常具有实际价值. 对此，众多研究 [10-14]提出不同的分析模型，并在不同受载条件，通过 理论计算，分析了黏弹性裂纹尖端场的力学 行为. 本文考查了功能梯度材料板条中界面裂纹垂直于梯度方向时的黏弹性断裂问题，首先利用有 限元法求解线弹性功能梯度材料板条的裂纹尖端场，然后根据黏弹性的对应性原理，求解出 黏弹性功能梯度材料板条裂纹问题的应力场强度因子.     

BRITTLE-DAMAGE ANALYSIS OF CRACKS UNDER CONDITIONS OF PLANE-STRAIN TENSILE LOADING

An abrupt damage model, taking full account of finite geometrychanges, is used to study both the shapes of damage zones and the stress strain fieldsfor a plane-strain tensile crack under the small-scale yielding condition. Two typicalcrack-tip damage profiles are simulated by the element vanish technique. The fracturetoughness increment due to damage dissipation is evaluated in terms of the energyrelease rate.     

Finite deformation analysis of crack-tip opening in elastic-plastic materials and implications for fracture

Analyses of the stress and strain fields around smoothly-blunting crack tips in both non-hardening and hardening elastic-plastic materials, under contained plane-strain yielding and subject to mode I opening loads, have been carried out by use of a finite element method suitably formulated to admit large geometry changes. The results include the crack-tip shape and near-tip deformation field, and the crack-tip opening displacement has been related to a parameter of the applied load, the J-integral. The hydrostatic stresses near the crack tip are limited due to the lack of constraint on the blunted tip, limiting achievable stress levels except in a very small region around the crack tip in power-law hardening materials. The J-integral is found to be path-independent except very close to the crack tip in the region affected by the blunted tip. Models for fracture are discussed in the light of these results including one based on the growth of voids. The rate of void-growth near the tip in hardening materials seems to be little different from the rate in non-hardening ones when measured in terms of crack-tip opening displacement, which leads to a prediction of higher toughness in hardening materials. It is suggested that improvement of this model would follow from better understanding of void-void and void-crack coalescence and void nucleation, and some criteria and models for these effects are discussed. The implications of the finite element results for fracture criteria based on critical stress or strain, or both, is discussed with respect to transition of fracture mode and the angle of initial crack-growth. Localization of flow is discussed as a possible fracture model and as a model for void-crack coalescence.     

Stress functions for central straight cracked anisotropic plates

In this paper, the plane problem for an anisotropic plate with a central straight crack in any direction is solved. The stress functions are given to represent the finite stress concentrations near the crack tips by the weight integral method. It shows that there is no stress singularity at the crack tip. The model can be used to appropriate to fracture mechanics for non-metallic materials.     

混凝土双K断裂参数计算的半解析有限元法

混凝土裂缝扩展的双$K$断裂准则，用于描述混凝土结构裂缝的起裂、稳定扩 展和失稳断裂. 其相应的双$K$断裂参数(起裂断裂韧度$K_{\rm IC}^{\rm ini} $和失 稳断裂韧度$K_{\rm IC}^{\rm un}$)一般通过简便的试验和基于虚拟裂缝扩展粘 聚力的解析方法确定. 利用平面扇形域哈 密顿体系的方程，通过分离变量法及共轭辛本征函数向量展开法，以解析的方法推导出基于混 凝土虚拟裂缝扩展线性粘聚力模型的平面裂缝解析元列式. 将该解析元与有限元相结合，构成 半解析的有限元法，可求解任意结构几何形状的混凝土平面裂缝双$K$断裂参数的计算问题. 数值计算结果表明半解析有限元法对该类问题的求解是十分有效的.     

On the uses of special crack tip elements in numerical rock fracture mechanics

Numerical methods such as boundary element methods are widely used for the stress analysis in solid mechanics. These methods are also used for crack analysis in rock fracture mechanics. There are singularities for the stresses and displacements at the crack tips in fracture mechanics problem, which decrease the accuracy of the numerical results in areas very close to the crack ends. To overcome this, higher order elements and isoperimetric higher order elements have been used. Recently, special crack tip elements have been proposed and used in most of the numerical fracture mechanics models. These elements can drastically increase the accuracy of the results near the crack tips, but in most of the models only one special crack tip element has been used for each crack end. In this study the uses of higher order crack tip elements are discussed and a higher order displacement discontinuity method is used to investigate the effect of these elements on the accuracy of the results in some crack problems. The useful shape functions for two special crack tip elements, are derived and given in the text and appendix for both infinite and semi-infinite plane problems. In this analysis both Mode I and Mode II stress intensity factors are computed . Some example problems are solved and the computed results are compared with the results given in the literature. The numerical results obtained here are in good agreement with those cited in the literature. For the curved crack problem, the strain energy release rate, G can be calculated accurately in the vicinity of the crack tips by using the higher order displacement discontinuity method with a quadratic variation of displacement discontinuity elements and with two special crack tip elements at each crack end.     

断裂力学判据的评述

从Inglis和Griffith的著名论文到Irwin和Rice等的奠基性贡献,对断裂力学中的线弹性断裂力学的K判据,界面断裂力学的G判据,和弹塑性断裂力学的J判据作了扼要的综述.介绍了在界面断裂力学G判据的基础上提出的界面断裂力学的K判据,以说明断裂力学的判据存在改进的可能性.在综述中归纳出断裂力学判据中目前还没有较好解决的几个问题.在总结以往断裂力学研究经验的基础上,指出裂纹端应力奇异性的源是对断裂力学判据存在的问题作进一步研究的切入点.探讨了裂纹端应变间断的奇点是裂纹端应力奇异性的源的问题,从而对裂纹端应力强度因子的物理意义进行了讨论.最后,阐述了进行可靠的裂纹端应力场的弹塑性分析是改进弹塑性断裂力学判据的关键,而进行可靠的裂纹端应力场的弹塑性分析的前提是要通过裂纹端应力奇异性的源的研究来获得作用在裂纹端的造成裂纹端应变间断的有限值应力.     

Simulation of glass cutting with an impinging hot air jet

Thermal cutting of glass sheet due to an impinging hot air jet is simulated and analyzed. Induced thermal stresses due to the moving heat source can be used to stably initiate and attract a crack toward the jet axis. Relative motion between the jet and glass sheet then can be used to cut the glass sheet. This paper presents a theoretical study of this process for straight cuts. Process simulation is accomplished by analyzing the coupled temperature and stress fields together with the fracture mechanics criteria for the crack growth.A finite element remeshing technique is employed for the analysis and singular elements are used around the crack tip for a more precise computing of the stress intensity factor. It is shown that a certain minimum air jet temperature for a given nozzle velocity and a certain maximum air jet velocity for a given temperature are required for continuous cutting. The results of the simulation show a good agreement with the published results in the literature. However a variating nature is detected for the distance between the crack tip and the air jet nozzle from a starting value to the steady-state one.     

爆炸荷载下青砂岩动态起裂韧度的测试方法

为了研究爆炸荷载下青砂岩I型裂纹动态断裂韧度的测试方法，利用内部中心单裂纹圆盘（internal center single crack disc，ICSCD）试样进行了爆炸试验研究。试样由外径为400 mm、内部加载孔径为40 mm、预制裂纹长为60 mm的青砂岩制成。利用同步触发器实现圆盘中心起爆，并同步触发超动态应变仪，通过径向应变片获取爆炸应变曲线、裂纹尖端的环向应变片获取裂纹起裂时刻。以实测爆炸应变曲线为参量，应用Laplace变换推导出试样加载孔壁应力时程曲线表达式，并用数值反演法得出其数值解。利用ANSYS有限元软件，建立数值计算模型，通过相互作用积分法得出了在爆炸荷载作用下砂岩的I型动态应力强度因子曲线。研究结果表明:（1）ICSCD试件能够很好地用来测试岩石的动态起裂韧度；（2）炮孔周边的应力可以通过拉普拉斯变换的数值反演方法得到；（3）通过试验-数值法能稳定计算出ICSCD砂岩构型的动态起裂韧度，其最大误差仅为7%。     

断裂过程的有限元模拟

讨论了材料断裂过程的有限元模拟技术。基于自适应有限元的一般原理，并针对多相材料的裂纹扩展的特点，提出了一种简化的高精度和高效率有限元网格的动态重新划分策略。裂纹被假设沿着单元之间的路径连续扩展，利用节点力释放技术生成新的裂纹自由表面，发展了一种可随裂尖连续移动的网格动态加密和释放方法。这种方法已在各种裂纹问题中得以实现与应用。