动载下复合材料基体裂纹与夹杂相互作用的数值模拟

采用ABAQUS软件及粘聚裂纹模型对动态拉伸载荷下复合材料中垂直于基体-夹杂界面的基体裂纹与夹杂的相互作用进行了数值模拟.结果表明:在给定的界面强度下,当加载率(或应变率)低于某一临界值(临界应变率)时裂纹将沿基体-夹杂界面扩展,当高于该临界值时裂纹可以穿过界面在夹杂中扩展.此外,随着界面强度的提高,临界应变率降低;当界面强度超过一定值后,裂纹扩展方向将不受外部载荷影响,裂纹将沿自相似方向扩展;当界面强度低于该临界值时,裂纹将沿界面扩展;并且临界应变率随夹杂尺寸的增大而降低,即小夹杂更难于破坏.上述结果可为前人的混凝土动态实验数据提供合理的细观理论解释.     

夹杂和裂纹相互影响时的界面应力分析

本文结合无限域上单根夹杂和单根裂纹的基本解,将裂纹与夹杂相互作用的问题归结为解一组柯西型奇异积分的积分方程组,使问题得到解决。本文还使用夹杂两侧的未知界面应力差,进一步推导了夹杂两侧的界面应力,并做了数值计算。有关这方面的计算可以作为研究与设计纤维与基体的联结强度的工程参考。     

线夹杂和线夹杂的相互作用

以短纤维复合材料为工程背景，本文利用线夹杂的工程计算模型以及无限平面中单夹杂的基本解，导出了线夹杂和线夹杂相互作用的平面问题的奇异积分方程。给出了夹杂端点的应力强度因子和夹杂界面应力的表达式，并作了具体的数值计算。     

纳米夹杂复合材料的有效反平面剪切模量研究

基于Gurtin-Murdoch表面/界面理论模型,利用复变函数方法,获得了考虑夹杂界面应力时夹杂/基体/等效介质模型的全场精确解,发展了能够预测纳米夹杂复合材料有效反平面剪切模量的广义自洽方法,给出了复合材料有效反平面剪切模量的封闭形式解。数值结果显示：当夹杂尺寸在纳米量级时,复合材料的有效反平面剪切模量具有尺度相关性,随着夹杂尺寸的增大,本文结果趋近于经典弹性理论的预测值;夹杂尺寸对于有效反平面剪切模量(本文结果)的影响范围要小于其对有效体积模量与剪切模量(各向同性材料)的影响范围;有效反平面剪切模量受夹杂的界面性能和夹杂刚度影响显著。     

含界面相Voronoi单元的等效弹性常数的计算

采用含界面相Voronoi单元有限元法,根据广义胡克定律,计算了在给定边界条件下,颗粒增强复合材料的等效弹性常数。建立了含多个随机分布的椭圆形夹杂及界面相的VCFEM模型,分析了夹杂体分比,界面相厚度和界面相弹性模量等因素对颗粒增强复合材料等效弹性常数的影响,并利用普通有限元方法对比验证。结果表明,当界面相弹性模量小于基体与夹杂时,材料的等效弹性模量会随着界面相厚度的增大而减小,随着夹杂体分比的增大而减小,并且界面过薄时,材料的等效弹性模量会随着夹杂体分比的增大而增大;当界面相弹性模量大于基体或夹杂时,材料的等效弹性模量会随着夹杂体分比和界面相厚度的增大而增大。而界面相的厚度和弹性模量对材料的等效泊松比的影响较小,材料的等效泊松比主要受夹杂体分比的影响,与其呈反比关系。     

含弧形裂纹复合陶瓷强度的尺度效应

认为含弧形裂纹复合陶瓷由随机方向的三相胞元与有效介质构成,用细观力学的方法研究了复合陶瓷的损伤失效和强度。首先确定三相胞元的外载应变,再依据复合陶瓷在损伤过程中的细观应力场和广义热力学力,计算出三相胞元内基体和颗粒的损伤等效应力,当基体和颗粒的损伤等效应力分别等于两者的极限应力时,得到基体和颗粒的破坏应力。然后,根据混合型应力强度因子计算弧形裂纹扩展时的能量释放率,进而得到界面的破坏应力。最后综合考虑基体、颗粒和和界面损伤影响,获得含弧形裂纹复合陶瓷的宏观强度及其尺度效应。     

增强相/夹杂内螺位错偶极子与含共焦钝裂纹椭圆夹杂的干涉效应

运用弹性力学的复势方法,研究了纵向剪切下增强相/夹杂内螺型位错偶极子与含共焦钝裂纹椭圆夹杂的干涉效应,得到了该问题复势函数的封闭形式解答,由此推导出了夹杂区域的应力场、作用在螺型位错偶极子中心的像力和像力偶矩以及裂纹尖端应力强度因子级数形式解。并分析了位错偶极子倾角 、钝裂纹尺寸和材料常数对位错像力、像力偶矩以及应力强度因子的影响。数值计算结果表明：位错像力、像力偶矩以及应力强度因子均随位错偶极子倾角做周期变化;夹杂内部的椭圆钝裂纹明显增强了硬基体对位错的排斥,减弱了软基体对位错的吸引,且对于硬夹杂,位错出现了一个不稳定平衡位置,该平衡位置随钝裂纹曲率的增大不断向界面靠近;变化 值将出现改变位错偶极子对应力强度因子作用方向的临界值。     

含微裂纹和椭球颗粒介质的强度及本构关系

针对含随机分布微裂纹及椭球颗粒的复合材料，通过考虑椭球颗粒内的本征应变及其与微裂纹的相互作用，利用等效夹杂方法研究了微裂纹损伤对材料有效模量和强度的影响，推导了复合材料的细观应力场及本构关系，并导出了材料破坏的临界条件．     

陶瓷颗粒增强金属基复合材料的细观强度分析

陶瓷颗粒增强金属基复合材料的失效主要有界面脱粘、增强粒子开裂等新的细观结构损伤机制。为了减小这些不足并对细观失效过程有一个清晰的了解,近来人们对金属基复合材料进行了大量研究,在此基础上,本文用细观力学的方法和损伤模型研究了陶瓷颗粒增强金属基复合材料的强度和损伤失效。为了计算方便,陶瓷颗粒简化为在复合材料中随机分布的椭球形粒子,然后以二相胞元模型计算分析了金属基体、颗粒中的应力应变分布情况,结果表明,基体中应力极不均匀,界面区存在应力集中,并计算了界面弧形裂纹扩展时的能量。最后分别提出了基体,颗粒和界面的失效强度准则,本文结果对于颗粒增强金属基复合材料具有普遍的实用性。     

开裂孔隙材料渗透率的细观力学模型研究

采用细观力学方法对含随机裂纹网络的孔隙材料渗透性进行研究.开裂孔隙材料渗透性的影响因素包括裂纹网络的密度、连通度、裂纹的开度以及孔隙材料基体渗透性.对于不连通的裂纹网络,该文采用已有的相互作用直推法(interaction direct derivative,IDD)的理论框架,引入裂纹的密度$\rho$和裂纹开度比$b$,提出了裂纹夹杂$\!$-$\!$-$\!$基体两相复合材料渗透率的IDD理论解.对于部分连通裂纹网络,考虑局部裂纹团内部各个裂纹对有效渗透率的相互放大作用,引入裂纹网络的连通度$f$,定义与连通度相关的水平裂纹密度$\rho^{h}$,按照增量法将表征连通特征的水平裂纹嵌入有效基体中,以此方式来考虑裂纹夹杂间的相互搭接,提出了考虑裂纹连通特征的扩展IDD理论解,分别考虑了基体材料渗透率$K_{m}$、裂纹密度$\rho $、裂纹开度比$b$以及与连通度$f$相关的$\rho ^{\rm h}$.最后通过对有限区域内含随机裂纹网络孔隙材料渗透过程的有限元模拟分别验证了不连通和部分连通裂纹网络扩展IDD模型的适用性:(1)当裂纹不连通时,由于基体对流体渗透的阻隔作用,裂纹的开度对有效渗透率影响不大;(2)当裂纹部分连通时,裂纹密度分别小于1.1(无关联裂纹网络,分形维数为2.0)、1.2(关联裂纹网络,分形维数为1.75)时,扩展IDD模型能够很好地估计开裂孔隙材料的有效渗透率,但是随着裂纹进一步扩展,最大裂纹团主导作用凸显,扩展IDD模型不再适用.      

Stress intensity factors for ring-shaped crack surrounded by spherical inclusions

Interaction of a ring-shaped crack with inhomogeneities such as inclusions is analyzed for the resulting three-dimensional stress field. Considered for the composite solid with a given volume fraction of inclusions are the two cases of (a) spherical voids and (b) spherical inclusions with elastic moduli different from the matrix. A ring-shaped crack is initiated at the equator of one of the voids or inclusions. A three-phase model is used to examine the interaction between the crack and surrounding inhomogeneities. Finite element method is then applied to calculate the stress intensity factor for different configurations. The effects of volume fraction of inhomogeneities, relative size of crack to inclusions, and material constants on crack behavior are discussed.     

弹塑性复合材料力学性能的细观研究

应用细观力学的Eshelby等效夹杂理论研究了复合材料的弹塑性问题。以铝基复合材料为例,建立了多轴载荷下复合材料弹塑性应力-应变关系,并且理论预报与实验结果符合较好,分析了夹杂形状、体积分数及加载路径对材料宏观性能的影响。同时,还研究了热塑性复合材料热膨胀系数与工艺温度之间的变化规律,分析了热残余应变对材料设计的影响。     

Prediction of overall tension behavior of short fiber-reinforced composites

The bridging stress of fibers along the crack surface plays an important role in analyzing the tension behavior of short or long fiber-reinforced composites. This paper uses the inclusion theory to obtain the expression of bridging stress for short fiber reinforced composite (SFRC) . A simplified model with periodically distributed fibers is proposed to estimate the average fiber spacings. The total fracture resistance is calculated as an energy summation including interface debonding energy dissipation, frictional sliding work between fibers and matrix, strain energy increment of fibers and matrix. The bend over point (BOP) stress is calculated by this fracture resistance. The necessary conditions of the fibers and matrix for the multiple cracking in SFRCs are discussed and the expression of ultimate external stress is derived. The critical fiber volume fraction for the strain hardening response is determined by an iteration method. In the meanwhile, the average spacing between two short fibers is proposed by a periodical distribution assumption. The theoretical prediction is compared with experimental data.     

Mechanical effects of circular liquid inclusions inside soft matrix: role of internal pressure change and surface tension

The mechanical effects of dilute liquid inclusions on the solid-liquid composite are explored,based on an analytical circular inclusion model incorporating the internal pressure change of the liquid and the surface tension of the interface.Several simple explicit dependences of the stress field and effective stiffness on the bulk modulus and the size of the liquid,the surface tension,and Poisson’s ratio of the matrix are derived.The results show that the stresses in the matrix are reduced,and the stiffness of the solid-liquid composite is enhanced with the consideration of either the surface tension or the internal pressure change.Particularly,the effective Young’s modulus predicted by the present model for either soft or stiff matrices agrees well with the known experimental data.In addition,according to the theoretical results,it is possible to stiffen a soft solid by pressured gas with the presence of the surface tension of the solid-gas interface.     

Effect of surface/interface stress on the plastic deformation of nanoporous materials and nanocomposites

Surface and interface play an important role on the overall mechanical behaviors of nanostructured materials. We investigate the effect of surface/interface stress on the macroscopic plastic behaviors of nanoporous materials and nanocomposites, where both the surface/interface residual stress and surface/interface elasticity are taken into account. A new second-order moment nonlinear micromechanics theory is developed and then reduced to macroscopically isotropic materials. It is found that the effect of surface/interface residual stress is much more prominent than that of the surface/interface elasticity, causing strong size effect as well as asymmetric plastic deformation for tension and compression. The variation of yield strength is more prominent with smaller pore/inclusion size or higher pore/inclusion volume fraction. For a representative nanoporous aluminum, the surface effect becomes significant when the pore radius is smaller than about 50 nm. When hard inclusions are embedded in a ductile metal matrix, the interface effect and resulting size effect are much smaller than that of nanoporous materials. The results may be useful for evaluating the mechanical integrity of nanostructured materials.     

Modeling segmentation cracking of a brittle coating due to underneath periodic eigenstrains

The formation mechanism for a segmented Cr coating by the hybrid technique of laser pre-quenching steel substrate plus post-electroplating was investigated. The discrete laser quenched zones (LQZs) were modeled as multiple inclusions with prescribed eigenstrains. The stress field was determined to account for the onset of segmentation cracking. Also addressed was segmentation crack growth through the evaluation of the stress intensity factor (SIF). The computations were implemented by using the finite element (FE) method. The dependence of a wide range of dimensionless variables of interest on both the stress field and SIF was assessed through a detailed parametric study.     

Interface edge crack in a multiferroic semicylinder

The purpose of the present work is to study the interface edge cracking problem in a semicircular cylindrical multiferroic composite theoretically by the methods of infinite series and singular integral equation. Numerical results of the stress intensity factor (SIF) are obtained and the computational accuracy is demonstrated. Discussion on the numerical results indicates that: (a) in order to reduce the SIF the inner piezoelectric layer should be softer and the outer piezomagnetic layer be stiffer; (b) the piezoelectric stiffening affect the SIF notably, but the piezomagnetic stiffening can only influence it very little.     

The effect of interfacial properties on the cohesion of highly filled composite materials

It has long been recognized that the cohesion of composite materials, in low confinement, is strongly affected by the properties of the interfacial transition zone (ITZ) between inclusions and matrix. While the effect of the ITZ on the elasticity properties of composites has been studied by many authors in the context of linear homogenization methods, the upscaling of the cohesion strength of highly filled composite materials has not been addressed. This is the focus of the non-linear homogenization procedure developed in this paper, which is based on the separation of the heterogeneous material system in phases of constant strength properties, a non-linear elastic representation of the limit stress state in each phase, and the definition of appropriate effective strain quantities that capture the morphological features of the microstructure. Applied to a three phase composite model composed of rigid inclusion, interface zone and matrix, the model provides a quantitative means of studying the effect of the interface cohesion and the interface volume fraction on the composite cohesion. In particular, we identify a critical interface-to-matrix cohesion ratio, below which the composite cohesion is smaller than the one of the matrix. Furthermore, the model lends itself readily to the study of the degradation of the interfacial properties in composite materials. This is shown for non-degraded and chemically softened cement-based materials, for which we provide conclusive evidence (1) that the interface strength properties of mortar are far more affected by chemical degradation than the one of the cement paste matrix; and (2) that chemical degradation does affect the mechanical strength performance of the cement paste not only through a change of volume proportions (i.e. increase of porosity), but as well through a pure chemical softening of the solid’s cohesion.     

A Zener–Stroh crack interacting with a coated inclusion with generalized Irwin plastic zone correction

The elastic–plastic fracture behavior of a Zener–Stroh crack interacting with a coated inclusion in composite materials has been investigated with crack tip plastic zone corrections. With the distributed dislocation method, the crack problem is formulated into a set of singular integral equations which are solved numerically. The plastic zone sizes at the both crack tips are determined by a generalized Irwin model where Von Mises stress yielding criterion is used. The stress intensity factor (SIF), the plastic zone size (PZS), the crack tip opening displacement (CTOD) and the effective stress intensity factor have been evaluated. In the numerical examples, the influence of the inclusion shear modulus, the coating-layer thickness and shear modulus, as well as the distance between the crack and inclusion, on the SIF, the PZS and the CTOD are discussed in detail. Numerical examples show that increasing the shear modulus or the thickness of the coating phase, the influence of the inclusion on the normalized SIF and the normalized PZS will be shielded.