A directional crack path criterion for crack growth in ductile materials subjected to shear and compressive loading under plane strain conditions

A directional crack growth criterion in a compressed elastic perfectly plastic material is considered. The conditions at the crack-tip are evaluated for a straight stationary crack with a small incipient kink. Remote load is a combined hydrostatic pressure and pure shear applied via a boundary layer. Crack surfaces in contact are assumed to develop homogenous Coulomb friction.The crack opening displacement of an extended kink is examined in a finite element analysis to judge the risk of opening mode failure. It has been found that the direction that maximizes the crack opening displacement of an extended kink tip coincides very well with a prediction of the crack growth direction obtained by using a criterion for continued crack growth direction discussed by the authors elsewhere [Int. J. Fract. 108 (2001) 351].Moreover, the by the model predicted incipient crack growth directions are qualitatively comparable with reported crack paths obtained in ductile materials in a limited number of experiments performed under a combined load of in-plane shear and compression.     

A MIXED MODE FRACTURE CRITERION BASED ON THE MAXIMUM TANGENTIAL STRESS IN BRITTLE INCLUSION

A closed-form solution for predicting the tangential stress of an inclusion located in mixed mode Ⅰ and Ⅱ crack tip field was developed based on the Eshelby equivalent inclusion theory. Then a mixed mode fracture criterion, including the fracture direction and the critical load, was established based on the maximum tangential stress in the inclusion for brittle inclusioninduced fracture materials. The proposed fracture criterion is a function of the inclusion fracture stress, its size and volume fraction, as well as the elastic constants of the inclusion and the matrix material. The present criterion will reduce to the conventional one as the inclusion having the same elastic behavior as the matrix material. The proposed solutions are in good agreement with detailed finite element analysis and measurement.     

Adhesion-induced instabilities in elastic and elastic-plastic contacts during single and repetitive normal loading

Adhesive interaction in spherical contacts was modeled with the Lennard-Jones (L-J) potential. Elastic adhesive contact was analyzed by the equivalent system of a rigid sphere with reduced radius of curvature and a half-space of effective elastic modulus. The critical gap at the instant of abrupt surface contact (jump-in) and separation (jump-out) was determined from the deformed surface profile of the elastic half-space and geometrical relationships. A finite element model of a rigid sphere and an elastic-plastic half-space was used to examine elastic-plastic adhesive contact. Surface adhesion was modeled by nonlinear springs with a force-displacement relationship governed by the L-J potential. The evolution of the interfacial force and the central gap distance as well as the occurrence of jump-in and jump-out instabilities were investigated in terms of the Tabor parameter, plasticity parameter, and dimensionless maximum normal displacement. The force-displacement response due to several approach-retraction cycles was interpreted in the context of elastic and plastic shakedown behaviors using dimensionless parameters.     

The Influence of Indenter Tip Imperfection and Deformability on Analysing Instrumented Indentation Tests at Shallow Depths of Penetration on Stiff and Hard Materials

We report on the difficulties of extracting plastic parameters from constitutive equations derived by instrumented indentation tests on hard and stiff materials at shallow depths of penetration. As a general rule, we refer here to materials with an elastic stiffness more than 10 % of that of the indenter and a yield strain higher than 1 %, as well as to penetration depths less than ～ 5 times the characteristic tip defect length of the indenter. We experimentally tested such a material (an amorphous alloy) by nanoindentation. To describe the mechanical response of the test, namely the force-displacement curve, it is necessary to consider the combined effects of indenter tip imperfections and indenter deformability. For this purpose, an identification procedure has been carried out by performing numerical simulations (using Finite Element Analysis) with constitutive equations that are known to satisfactorily describe the behaviour of the tested material. We propose a straightforward procedure to address indenter tip imperfection and deformability, which consists of firstly taking account of a deformable indenter in the numerical simulations. This procedure also involves modifying the experimental curve by considering a truncated length to create artificially the material’s response to a perfectly sharp indentation. The truncated length is determined directly from the loading part of the force-displacement curve. We also show that ignoring one or both of these issues results in large errors in the plastic parameters extracted from the data.     

纳米压痕过程的三维有限元数值试验研究

采用有限元方法模拟了纳米压痕仪的加、卸载过程，三维有限元模型考虑了纳米压痕仪的标准Berkovich压头．介绍了有限元模型的几何参数、边界条件、材料特性与加载方式，讨论了摩擦、滑动机制、试件模型的大小对计算结果的影响，进行了计算结果与标准试样实验结果的比较，证实了模拟的可靠性．在此基础上，重点研究了压头尖端曲率半径对纳米压痕实验数据的影响．对比分析了尖端曲率半径r＝0与r＝100nm两种压头的材料压痕载荷—位移曲线．结果表明，当压头尖端曲率半径r≠0时，基于经典的均匀连续介质力学本构理论、传统的实验手段与数据处理方法，压痕硬度值会随着压痕深度的减小而升高．     

基于均匀化理论的页岩微观多孔黏土强度特性

页岩强度是页岩油气开发所必需的基础技术参数之一,对页岩强度的研究贯穿于钻完井、压裂工艺施工的全过程.常规宏观室内实验存在试样获取困难、耗时较长,受井下工矿制约,地球物理方法获取资料品质欠佳且增加了井下设备卡、埋风险.因此,提出基于均匀化理论评价页岩微观多孔黏土强度的方法,进行多孔黏土组成与力学分析.基于耗散能原理和Drucker-Prager准则,开展了微观多孔黏土的强度与$\pi$函数的应变求解分析;讨论黏土颗粒与粒间孔隙的力学特性,建立多孔黏土的均匀化应变能;采用强度均匀化理论构建微观非线性函数模型,建立与多孔黏土组成、摩擦系数、内聚系数等参数相关的均匀化函数模型;基于纳米力学实验、量纲分析和有限元模拟,分析多孔黏土硬度、强度与组成的内在关系.研究结果表明,页岩微观多孔黏土的弹性模量和硬度与黏土堆积密度正相关,当黏土堆积密度一定时,硬度与内聚系数的比值受摩擦系数影响较大,为非线性递增;通过量纲分析和有限元模拟,求解页岩微观多孔黏土关于硬度--强度--堆积密度的$\pi$函数,揭示页岩微观黏土矿物的组成与力学性质的关系,为进一步深入研究页岩细观强度参数和宏观强度预测奠定基础.      

准静态条件下金属材料的临界断裂准则研究

本文针对9种金属材料完成了具有不同约束程度的10类试样的延性断裂试验, 获得了发生拉、压、扭和裂尖断裂等破坏形式构型试样的载荷-位移试验关系; 基于圆棒漏斗试样拉伸试验所得直至破坏的载荷-位移曲线, 采用有限元辅助试验(finite-element-analysis aided testing, FAT)方法得到了9种材料直至破坏的全程等效应力-应变曲线, 以此作为材料本构关系通过有限元分析获得了各类试样直至临界破坏的载荷-位移关系模拟. 载荷-位移关系模拟结果与试验结果有较好的一致性, 表明用于解决试样颈缩问题的FAT方法所获得的全程材料本构关系针对各向同性材料具有真实性和普适性. 对应9种材料、10类试样的36 个载荷-位移临界断裂点, 通过有限元分析获得了对应的材料临界断裂应力、应变与临界应力三轴度, 研究表明, 第一主应力在延性变形过程中为主控断裂的主导参量; 通过研究光滑、缺口、裂纹等构型试样的断裂状态, 提出了$-1$至3范围的应力三轴度下由第一主应力主控的统一塑性临界断裂准则.      

Crushing of particles in idealised granular assemblies

Four idealised assemblies of equally sized spherical particles are subjected to a range of macroscopic compressive principal stresses and the contact forces on individual particles are determined. For each set of contact forces the stress fields within individual particles are studied. A failure criterion for brittle materials is imposed and indicates that crushing (or rupture) occurs when the maximum contact force reaches a threshold particle strength value, irrespective of the presence and magnitude of other lesser contact forces acting on the particle and the material properties of the particle. Combining the crushing mechanism with an assembly instability mechanism enables failure surfaces to be drawn in the three-dimensional stress space. A simple spatial averaging technique has been applied to the failure surfaces to remove the effects of assembly anisotropies. Sections of the failure surfaces on π planes have similarities to those commonly used in sand modelling.     

Micromechanical study of elastic moduli of three-dimensional granular assemblies

In micromechanics of granular materials, relationships are investigated between micro-scale characteristics of particles and contacts and macro-scale continuum characteristics. For three-dimensional isotropic assemblies the macro-scale elastic characteristics are described by the bulk and the shear modulus, which depend on the micro-scale characteristics of the coordination number (i.e. the average number of contacts per particle) and the interparticle contact stiffnesses in directions normal and tangential to the contact.     

Two new expanding cavity models for indentation deformations of elastic strain-hardening materials

Two expanding cavity models (ECMs) are developed for describing indentation deformations of elastic power-law hardening and elastic linear-hardening materials. The derivations are based on two elastic–plastic solutions for internally pressurized thick-walled spherical shells of strain-hardening materials. Closed-form formulas are provided for both conical and spherical indentations, which explicitly show that for a given indenter geometry indentation hardness depends on Young’s modulus, yield stress and strain-hardening index of the indented material. The two new models reduce to Johnson’s ECM for elastic-perfectly plastic materials when the strain-hardening effect is not considered. The sample numerical results obtained using the two newly developed models reveal that the indentation hardness increases with the Young’s modulus and strain-hardening level of the indented material. For conical indentations the values of the indentation hardness are found to depend on the sharpness of the indenter: the sharper the indenter, the larger the hardness. For spherical indentations it is shown that the hardness is significantly affected by the strain-hardening level when the indented material is stiff (i.e., with a large ratio of Young’s modulus to yield stress) and/or the indentation depth is large. When the indentation depth is small such that little or no plastic deformation is induced by the spherical indenter, the hardness appears to be independent of the strain-hardening level. These predicted trends for spherical indentations are in fairly good agreement with the recent finite element results of Park and Pharr.     

考虑等效曲率的超二次曲面单元非线性接触模型

基于连续函数包络的超二次曲面单元可有效地描述自然界和工业生产中的非球体颗粒形态, 并通过非线性迭代方法精确计算单元间的接触力. 对于具有复杂几何形态的超二次曲面单元, 线性接触模型不能准确地计算不同接触模式下的作用力. 考虑超二次曲面单元相互作用时不同颗粒形状及表面曲率的影响, 本文发展了相应的非线性黏弹性接触模型. 该模型将不同接触模式下的法向刚度和黏滞力统一表述为单元间局部接触点处等效曲率半径的函数; 切向接触作用则借鉴基于Mohr-Coulomb摩擦定律的球体单元非线性接触模型的计算方法. 为检验超二次曲面单元接触模型的可靠性, 对球形颗粒间的法向碰撞、椭球体颗粒间的斜冲击过程、圆柱体的静态堆积和椭球体的动态卸料过程进行离散元模拟, 并与有限元数值结果及试验结果进行对比验证. 计算表明, 考虑接触点处等效曲率半径的超二次曲面非线性接触模型可准确地计算单元间的接触碰撞作用, 并合理地反映非球形颗粒体系的运动规律. 在此基础上进一步分析了不同长宽比和表面尖锐度对卸料过程中颗粒流动特性的影响, 为非球形颗粒材料的流动特性分析提供了一种有效的离散元方法.      

Influence of penetration depth and mechanical properties on contact radius determination for spherical indentation

Knowledge of the relationship between the penetration depth and the contact radius is required in order to determine the mechanical properties of a material starting from an instrumented indentation test. The aim of this work is to propose a new penetration depth–contact radius relationship valid for most metals which are deformed plastically by parabolic and spherical indenters. Numerical simulation results of the indentation of an elastic–plastic half-space by a frictionless rigid paraboloïd of revolution show that the contact radius–indentation depth relationship can be represented by a power law, which depends on the reduced Young’s modulus of the contact, on the strain hardening exponent and on the yield stress of the indented material. In order to use the proposed formulation for experimental spherical indentations, adaptation of the model is performed in the case of a rigid spherical indenter. Compared to the previous formulations, the model proposed in the present study for spherical indentation has the advantage of being accurate in the plastic regime for a large range of contact radii and for materials of well-developed yield stress. Lastly, a simple criterion, depending on the material mechanical properties, is proposed in order to know when piling-up appears for the spherical indentation.     

Effects of large deformation and material nonlinearity on spherical indentation of hyperelastic soft materials

Up to now, the indentation of hyperelastic soft materials has not been completely understood. In this paper, the spherical indentation on hyperelastic soft solids was systematically investigated through theoretical analysis and finite element method (FEM). The validation and application of the Hertzian load-displacement relation for indentation of hyperelastic soft materials were clarified, the effects of large deformation and material nonlinearity on spherical indentation of hyperelastic soft materials were analyzed and discussed. It was found that the complicated indentation behaviors of hyperelastic soft solids mainly depended on the coupling interactions of large deformation and material nonlinearity. Besides, we proposed two new nonlinear elastic contact models to separate the effects of large deformation and material nonlinearity on spherical indentation responses of hyperelastic soft solids. Our efforts might help to enhance the understanding of hyperelastic indentation problems and provided necessary instructions for the mechanical characterization of hyperelastic soft materials.     

A MOVING CRACK IN A NONHOMOGENEOUS MATERIAL STRIP

This paper considers an anti-plane moving crack in a nonhomogeneous material strip of finite thickness. The shear modulus and the mass density of the strip are considered for a class of functional forms for which the equilibrium equation has analytical solutions. The problem is solved by means of the singular integral equation technique. The stress field near the crack tip is obtained. The results are plotted to show the effect of the material non-homogeneity and crack moving velocity on the crack tip field. Crack bifurcation behaviour is also discussed. The paper points out that use of an appropriate fracture criterion is essential for studying the stability of a moving crack in nonhomogeneous materials. The prediction whether the unstable crack growth will be enhanced or retarded is strongly dependent on the type of the fracture criterion used. Based on the analysis, it seems that the maximum 'anti-plane shear' stress around the crack tip is a suitable failure criterion for moving cracks in nonhomogeneous materials.     

Dual-parameter elastic-plastic fracture criterion of metals and alloys in the plane strain case

In recent years, some investigators discussed the applicability of the HRR theory for engineering materials based on the results of numerical analyses and experimental studies. In the present paper, the finite element method is employed to analyze the crack tip fields of the engineering elastic-plastic material with a variety of geometry configurations of cracked specimens from elastic state to intensely general yielded state in the plane strain case. The results indicate that the HRR theory loses its validity of application for engineering elastic-plastic materials in the plane strain case. The reasons for this are analyzed. A dual-parameter fracture criterion is suggested for this case.     

Refined Series Methodology for the Fully Coupled Thin-Walled Laminated Beams Considering Foundation Effects

The refined power series solutions are presented for the coupled static analysis of thin-walled laminated beams resting on elastic foundation. For this purpose, the elastic strain energy considering the material and structural coupling effects and the energy including the foundation effects are constructed. The equilibrium equations and the force-displacement relationships are derived from the extended Hamilton's principle, and the explicit expressions for displacement parameters are presented based on power series expansions of displacement components. Finally, the member stiffness matrix is determined by using the force-displacement relationships. For comparison, the finite element model based on the Hermite cubic interpolation polynomial is presented. In order to verify the accuracy and the superiority of the laminated beam element developed by this study, the numerical solutions are presented and compared with results obtained from the regular finite beam elements and the ABAQUS's shell elements. The influences of the fiber angle change and the boundary conditions on the coupled behavior of laminated beams with mono-symmetric I-sections are investigated.     

Contact laws between solid particles

This paper presents a comprehensive study for the contact laws between solid particles taking into account the effects of plasticity, strain hardening and very large deformation. The study takes advantage of the development of a so-called material point method (MPM) which requires neither remeshing for large deformation problems, nor iterative schemes to satisfy the contact boundary conditions. The numerical results show that the contact law is sensitive to impact velocity and material properties. The contact laws currently used in the discrete element simulations often ignore these factors and are therefore over-simplistic. For spherical particles made of elastic perfectly plastic material, the study shows that the contact law can be fully determined by knowing the relative impact velocity and the ratio between the effective elastic modulus and yield stress. For particles with strain hardening, the study shows that it is difficult to develop an analytical contact law. The same difficulty exists when dealing with particles of irregular shapes or made of heterogeneous materials. The problem can be overcome by using numerical contact laws which can be easily obtained using the material point method.     

Non-slipping JKR model for transversely isotropic materials

A generalized plane strain JKR model is established for non-slipping adhesive contact between an elastic transversely isotropic cylinder and a dissimilar elastic transversely isotropic half plane, in which a pulling force acts on the cylinder with the pulling direction at an angle inclined to the contact interface. Full-coupled solutions are obtained through the Griffith energy balance between elastic and surface energies. The analysis shows that, for a special case, i.e., the direction of pulling normal to the contact interface, the full-coupled solution can be approximated by a non-oscillatory one, in which the critical pull-off force, pull-off contact half-width and adhesion strength can be expressed explicitly. For the other cases, i.e., the direction of pulling inclined to the contact interface, tangential tractions have significant effects on the pull-off process, it should be described by an exact full-coupled solution. The elastic anisotropy leads to an orientation-dependent pull-off force and adhesion strength. This study could not only supply an exact solution to the generalized JKR model of transversely isotropic materials, but also suggest a reversible adhesion sensor designed by transversely isotropic materials, such as PZT or fiber-reinforced materials with parallel fibers.     

Novel material tailoring method for internally pressurized FG spherical and cylindrical vessels

A new material tailoring method for spherical and cylindrical vessels made of functionally graded materials (FGMs) is presented. It is assumed that the FG material is composed of an Al–SiC metallic-matrix composite. A uniform ratio of in-plane shear stress to yield strength [\(\varphi \left( r \right) \)] is used as the design criterion to utilize the maximum capacity of the vessel. The aim is to find a distribution of SiC particles in the radial direction, i.e., \(f\left( r \right) \), that achieves a uniform index \(\varphi \left( r \right) =\hbox {const}.\) through the wall thickness of the internally pressurized spherical or cylindrical vessel. Both the Mori–Tanaka and rule-of-mixtures homogenization schemes are used to express the effective elastic module and Poisson’s ratio. Moreover, the strength of the composite is expressed based on the rule of mixtures. Besides, finite element simulation is carried out to verify the accuracy of the analytical solution. The effects of input parameters such as the internal pressure, strength of the SiC particles, ratio of in-plane shear stress to effective yield strength, and choice of homogenization scheme on the tailored distribution of the SiC volume fraction in the radial direction are also investigated.