Modelling matrix damage and fibre–matrix interfacial decohesion in composite laminates via a multi-fibre multi-layer representative volume element (MRVE)

A three-dimensional multi-fibre multi-layer micromechanical finite element model was developed for the prediction of mechanical behaviour and damage response of composite laminates. Material response and micro-scale damage mechanism of cross-ply, [0/90]_ns, and angle-ply, [±45]_ns, glass-fibre/epoxy laminates were captured using multi-scale modelling via computational micromechanics. The framework of the homogenization theory for periodic media was used for the analysis of the proposed ‘multi-fibre multi-layer representative volume element’ (M²RVE). Each layer in M²RVE was represented by a unit cube with multiple randomly distributed, but longitudinally aligned, fibres of equal diameter and with a volume fraction corresponding to that of each lamina (equal in the present case). Periodic boundary conditions were applied to all the faces of the M²RVE. The non-homogeneous stress–strain fields within the M²RVE were related to the average stresses and strains by using Gauss’ theorem in conjunction with the Hill–Mandal strain energy equivalence principle. The global material response predicted by the M²RVE was found to be in good agreement with experimental results for both laminates. The model was used to study effect of matrix friction angle and cohesive strength of the fibre–matrix interface on the global material response. In addition, the M²RVE was also used to predict initiation and propagation of fibre–matrix interfacial decohesion and propagation at every point in the laminae.     

Particle fracture in high-volume-fraction ceramic-reinforced metals: Governing parameters and implications for composite failure

Weibull parameters of angular alumina particles are determined from experimental tensile test data on high-ceramic-content metal matrix composites using a micromechanical model that accounts for internal damage in the form of particle cracking, the dominant damage mode in these composites. The fraction of broken particles is assessed from the drop of Young's modulus and particle fracture is assumed to be stress controlled. Two extreme load-sharing modes, namely a purely local and a global load-sharing mode, are considered to account for the load redistribution due to particle fracture. Consistent powder strength parameters can be thus “back-calculated” for particles that are embedded in different Al-Cu matrices. On the other hand, this calculation fails for pure Al matrix composites, which exhibit a much larger strain to failure than Al-Cu matrix composites. It is shown that for Al matrix composites, the role of plastic (composite) strain on particle fracture constitutes a second parameter governing particle damage. This finding is rationalized by particle-particle interactions in these tightly packed ceramic particle-reinforced composites, and by the increase of matrix stress heterogeneity that is brought with increasing plastic strain. Failure of the alloyed matrix composites is well described by the (lower bound) local load-sharing micromechanical model, which predicts a catastrophic failure due to an avalanche of damage. The same model predicts failure of pure aluminium matrix composites to occur at the onset of tensile instability, also in agreement with experimental results once the role of plastic strain on damage accumulation is accounted for.     

Progressive failure analysis of tow-placed, variable-stiffness composite panels

The past developments on tow-placement technology led to the production of machines capable of controlling fibre tows individually and placing them onto the surface of a laminate with curvilinear topology. Due to the variation of properties along their surface, such structures are termed variable-stiffness composite panels.In previous experimental research tow-steered panels have shown increased buckling load capacity as compared with traditional straight-fibre laminates. Also, numerical analyses by the authors showed that first-ply failure occurs at a significant higher load level. The focus of this paper is to extend those analyses into the postbuckling progressive damage behaviour and final structural failure due to accumulation of fibre and matrix damage. A user-developed continuum damage model implemented in the finite element code ABAQUS^® is employed in the simulation of damage initiation and material stiffness degradation.In order to correctly predict the buckling loads of tow-steered panels under compression, it is of crucial importance to take into account the residual thermal stresses resulting from the curing process. Final failure of tow-steered panels in postbuckling is predicted to within 10% difference of the experimental results. Curvilinear-fibre panels have up to 56% higher strength than straight-fibre laminates and damage initiation is also remarkably postponed. Tow-steered designs also show more tolerance to central holes than traditional laminates.     

菱形颗粒冲击材料表面冲蚀磨损特性分析

基于弹射试验装置,借助高速摄像机捕捉不同入射条件下单个菱形颗粒冲击金属表面的动态过程,同时结合试验过程建立菱形颗粒冲击金属表面的FEM-SPH耦合数值模型,通过对比试验现象与仿真结果优化数值模型参数,最后借助数值模型进一步分析菱形颗粒在临界冲击、自身初始旋转以及重复冲击等工况下的运动行为及预测的凹坑轮廓形态. 结果表明:优化后的模型能够很好地捕捉颗粒冲击过程中金属表面凹坑的产生及演化规律,并能详细记录颗粒的入射行为及反弹规律,测得颗粒反弹速度和反弹角度误差均在14%以内. 临界冲击工况下颗粒动能损失最大,且冲击角越高,残余动能越少;颗粒初始旋转能够改变其反弹后的运动行为及金属表面材料的失效方式;颗粒重复冲蚀对材料表面的作用机制与后续颗粒的入射条件有密切关系,模型成功捕捉到重复冲蚀导致的材料破坏加深和破坏减缓两种特殊现象.      

磨料水射流冲蚀金属损伤形貌特征及机理研究

针对前混合磨料水射流冲蚀过程,进行了冲蚀试验及仿真研究.从试验中得到了材料冲蚀损伤形貌特征.采用SPH耦合FEM方法建立相应的冲蚀模型,对材料冲蚀过程进行模拟分析,揭示了材料损伤形貌特征产生的机理.结果表明:随着冲蚀深度的增加,磨料颗粒的冲蚀动能逐渐减小,冲蚀角度逐渐增大,材料冲蚀损伤由微切削和微犁削逐渐变为冲击变形,材料冲蚀损伤断面的形貌特征逐渐恶化,具体表现为拖尾角和表面粗糙度逐渐增大.不同金属材料损伤断面形貌特征具有相似性,但是金属材料属性的差异会对磨料冲蚀过程产生影响,导致条纹在材料损伤断面上的分布和条纹角度出现差异.由于材料损伤形貌特征受控于磨料颗粒运动特性,因此材料冲蚀断面质量改善应该从改变磨料颗粒运动特性角度出发,这为进一步研究材料冲蚀断面质量改进奠定了理论基础.     

双滑移转动率模型对椭圆散粒体的适用性

首先基于椭圆颗粒接触点的运动关系推导出新的平均纯转动率（APR）表达式,并将其引入到已有的描述圆形散粒体流动特性的双滑移转动率模型（DSR²模型）中;而后采用改进的NS2D离散元程序对长短轴比例分别为1.4和1.7的椭圆颗粒堆积体进行一系列不排水单剪试验,验证以椭圆颗粒为基础的离散元方法模拟砂土流动特性的可行性及DSR²模型的合理性。数值试验及已有成果表明,以椭圆颗粒为基础的NS2D程序能够模拟砂土的流动特性,对圆形和椭圆形颗粒体系,DSR²模型均能很好地预测运动模型中转动率参量的变化情况;APR是联系连续介质力学和离散介质力学的重要参数,它将二者有机结合成一个整体。     

基于Rayleigh-Ritz算法的梁裂缝损伤识别及试验研究

结构动态特性的变化则预示结构出现损伤,基于Ritz线性近似,文中提出一种裂缝损伤识别的方法,用以识别结构裂缝损伤位置及损伤程度。该方法分两步：首先用模态子空间近似关系．消去模态应变能表达式中的整体刚度矩阵和整体质量矩阵,避免模型误差对识别结果的影响。再采用计算向量空间夹角的方法分离损伤位置与损伤程度的耦合影响,进而识别出单元裂缝损伤位置;其次,识别损伤程度采用二次线性规划方法,不再计算特征值灵敏度。线性约束条件保证了二次规划问题的解是唯一的。模拟筒支梁几种裂缝损伤情况进行数值计算与模态试验,利用所得模态参数对该算法程序进行了验证,识别出了裂缝损伤的确切位置及损伤程度,并进行了误差对比。结果表明,该算法由于不用整体结构的数值模型,从而避免了边界条件、连接条件及材料特性参数等因素对识别结果的干扰,识别精度得到提高,将其用于结构损伤识别是可行的。     

The influence of particle shape,volume fraction and distribution on post-necking deformation and fracture in uniaxial tension of AA5754 sheet materials

Finite element analysis was performed over a small particle field, edge constraint plane strain post-necking model. The aim is to understand the roles of particle shape, volume fraction and distribution over the post-necking deformation and fracture of AA5754-O sheet materials. For models containing one single particle, the post-necking deformation decreases when the particle varies from circular to elliptical. The inter-particle spacing, the major parameter of distribution to determine whether a pair of particles belongs to a stringer or not, was varied for models with two particles of circular or elliptical shape. The general trend is that the post-necking deformation and fracture strains decrease with decreasing spacing between particles. There is considerable difference in terms of both fracture topographies and strains for models containing 16 particles when distributions varied from random/uniform to stringer distributions. The post-necking deformation and fracture strains monotonically decrease with particle volume fractions for models with 4–64 particles of random or stringer distribution. This indicates that the post-necking behavior for AA5754-O alloys where the matrix material is rather ductile is not solely controlled by a single or pair of particles although they may become initiation places of damage. Multiple damaging sources such as stringers or large particles can act cooperatively and speed up the damaging propagation of the material, and therefore produce small post-necking deformation and early fracture. The center clustering of particles can be beneficial for post-necking behavior and bendability of sheet materials.     

On the extension of the Gurson-type porous plasticity models for prediction of ductile fracture under shear-dominated conditions

One of the major drawbacks of the Gurson-type of porous plasticity models is the inability of these models to predict material failure under low stress triaxiality, shear dominated conditions. This study addresses this issue by combining the damage mechanics concept with the porous plasticity model that accounts for void nucleation, growth and coalescence. In particular, the widely adopted Gurson–Tvergaard–Needleman (GTN) model is extended by coupling two damage parameters, representing the volumetric damage (void volume fraction) and the shear damage, respectively, into the yield function and flow potential. The effectiveness of the new model is illustrated through a series of numerical tests comparing its performance with existing models. The current model not only is capable of predicting damage and fracture under low (even negative) triaxiality conditions but also suppresses spurious damage that has been shown to develop in earlier modifications of the GTN model for moderate to high triaxiality regimes. Finally the modified GTN model is applied to predict the ductile fracture behavior of a beta-treated Zircaloy-4 by coupling the proposed damage modeling framework with a recently developed J₂–J₃ plasticity model for the matrix material. Model parameters are calibrated using experimental data, and the calibrated model predicts failure initiation and propagation in various specimens experiencing a wide range of triaxiality and Lode parameter combinations.     

Damage in a viscoplastic material part I: Cavity growth

The exact velocity, stress and strain rate fields around a spheroidal cavity in an infinite linear viscoplastic compressible matrix are derived analytically by the ‘three function approach’. The perturbation of the velocity field due to the cavity is the superposition of three independent modes, inducing homothetic growth, pure distortion and both volume and shape changes, respectively. This solution is then used to investigate the velocity field around a spheroidal cavity in a nonlinear viscous compressible material by means of a variational principle. The behaviour of such damaged linear and nonlinear materials will be discussed in a forthcoming companion paper.The importance of the reference strain, while studying void growth in a compressible material, is emphasized. If the axial strain is chosen as a reference, void growth is found to be enhanced at low triaxiality ratios, but lowered at high triaxiality ratios in a compressible matrix relative to an incompressible one. Finally, the transition from a power law to a linear law with intercept, at increasing strain rates, is shown to reduce damage growth rate.     

复合材料层合板低速冲击损伤分析的连续介质损伤力学模型

针对复合材料层合板低速冲击损伤问题,提出了一种各向异性材料连续介质损伤力学模型,模型涵盖损伤表征、损伤起始判定和损伤演化法则3 个方面. 通过材料断裂面坐标下的损伤状态变量矩阵完成损伤表征,并考虑断裂面角度的影响,建立了主轴坐标系下的材料损伤本构关系. 损伤起始由卜克(Puck) 失效准则预测,损伤演化由断裂面上的等效应变控制,服从基于材料应变能释放的线性软化行为. 模型区分了纤维损伤和基体损伤,并根据冲击载荷下层内产生多条基体裂纹继而扩展至界面形成层间裂纹(分层) 的试验观察,引入基体裂纹饱和密度参数表征层间分层. 以[0₃/45/-45]_S 和[45/0/-45/90]_4S 两种铺层的复合材料层合板为例,预测了不同冲击能量下复合材料层合板的低速冲击损伤响应参数,试验结果证明了连续介质损伤力学模型的有效性.模型在不同网格密度下的计算结果表明单元特征长度的引入可以在一定程度上降低损伤演化阶段对网格密度的依赖性.      

Multiscale modeling of the damaged plastic behavior and failure of Al/SiCp composites

We studied the tensile behavior and damage of an aluminium X2080 reinforced with different volume fractions of silicon carbide particles. The main damage mechanism is particle failure. Regions of the matrix adjacent to broken particles are sites with high hydrostatic tension and hence the nucleation of cavities is expected. Using J integral and HRR modified solution it is possible to calculate the growth of these voids. Macroscopic failure is governed by a critical volume fraction of voids. The originality of this work lies in the modeling of the composite using a micromechanical approach.     

Experimental simulation of matrix cracking and debonding in a model brittle matrix composite

A model composite material system was designed to simulate typical damage mechanisms in unidirectional fiber reinforced brittle matrix composites. Experiments were performed at low to high quasistatic, macroscopic loading rates . At all loading rates reversal of the transverse strain was observed and was correlated to matrix cracking and debonding. The optical method of coherent gradient sensing (CGS) was used to obtain qualitative information regarding the stress fields and to observe the progression of damage. It was found that the sequence of damage formation (damage path) depended on the macroscopic loading rate. At lower loading rates periodic matrix cracks developed; minimal debonding of the reinforcement-matrix interface occurred only much later in the experiment. At higher loading rates extensive debonding followed propagation of the initial matrix crack, and periodic cracking was not observed. Several features of the material response of the model material system were also observed in a previously studied unidirectional ceramic matrix composite.     

Deep Bed Filtration Modelling of Formation Damage Due to Particulate Invasion from Drilling Fluids

A Deep Bed Filtration model has been developed to quantify the effect of solids invasion from drilling fluids on the permeability of rock formations. The calculated particle-trapping profiles are compared directly with experimental profiles from scanning electron microscopy and synchrotron X-ray diffraction tomography mapping. The computed permeability reduction as a consequence of particle invasion is in broad agreement with experiment. Backflow was modelled by reversing the flow rate, starting off with a situation where all particles either remain trapped or are all released. It appears that the experimentally observed 30% release of particles upon backflow is reproducible within the limits of the two extreme cases. When erosion is included in the model, a peak in the backflow pressure time series can be observed. This peak may be correlated with the experimentally observed flow initiation pressure, which is the backflow pressure needed to initiate flow after initial inflow filtration. Finally, we conclude that internal reservoir damage, within the limits of our 1-D single phase DBF model, may contribute to the experimentally observed flow initiation pressure.     

A thermodynamics based damage mechanics model for particulate composites

A micro-mechanical damage model is proposed to predict the overall viscoplastic behavior and damage evolution in a particle filled polymer matrix composite. Particulate composite consists of polymer matrix, particle fillers, and an interfacial transition interphase around the filler particles. Yet the composite is treated as a two distinct phase material, namely the matrix and the equivalent particle-interface assembly. The CTE mismatch between the matrix and the filler particles is introduced into the model. A damage evolution function based on irreversible thermodynamics is also introduced into the constitutive model to describe the degradation of the composite. The efficient general return-mapping algorithm is exploited to implement the proposed unified damage coupled viscoplastic model into finite element formulation. Furthermore, the model predictions for uniaxial loading conditions are compared with the experimental data.     

Effect of cavity aspect ratio on flow and heat transfer characteristics in pipes: a numerical study

Using the standard k–ε turbulence model, a two-dimensional turbulent pipe flow was simulated with and without square cavities. Effect of cavity aspect ratio on flow and heat transfer characteristics was investigated. Uncertainty was approximated through experimental validation and grid independence. The simulation revealed circulation inside the cavities. Cavity boundaries were shown to contribute significantly toward turbulence production. Cavity presence was shown to enhance overall heat transfer through the wall, while increasing pressure drop significantly across the pipe. It was predicted that cavities with higher aspect ratio enhance heat transfer more while increasing pressure drop.     

Plastic potentials for anisotropic porous solids

The aim of this paper is to incorporate plastic anisotropy into constitutive equations of porous ductile metals. It is shown that plastic anisotropy of the matrix surrounding the voids in a ductile material could have an influence on both effective stress–strain relation and damage evolution. Two theoretical frameworks are envisageable to study the influence of plastic flow anisotropy: continuum thermodynamics and micromechanics. By going through the Rousselier thermodynamical formulation, one can account for the overall plastic anisotropy, in a very simple manner. However, since this model is based on a weak coupling between plasticity and damage dissipative processes, it does not predict any influence of plastic anisotropy on cavity growth, unless a more suitable choice of the thermodynamical force associated with the damage parameter is made. Micromechanically-based models are then proposed. They consist of extending the famous Gurson model for spherical and cylindrical voids to the case of an orthotropic material. We derive an upper bound of the yield surface of a hollow sphere, or a hollow cylinder, made of a perfectly plastic matrix obeying the Hill criterion. The main findings are related to the so-called ‘scalar effect’ and ‘directional effect’. First, the effect of plastic flow anisotropy on the spherical term of the plastic potential is quantified. This allows a classification of sheet materials with regard to the anisotropy factor h; this is the scalar effect. A second feature of the model is the plasticity-induced damage anisotropy. This results in directionality of fracture properties (‘directional effect’). The latter is mainly due to the principal Hill coefficients whilst the scalar effect is enhanced by ‘shear’ Hill coefficients. Results are compared to some micromechanical calculations using the finite element method.     

A constitutive model for cavitation and cavity growth in rubber-like materials under arbitrary tri-axial loading

In this paper, we attempted to construct a constitutive model to deal with the phenomenon of cavitation and cavity growth in a rubber-like material subjected to an arbitrary tri-axial loading. To this end, we considered a spherical elementary representative volume in a general Rivlin’s incompressible material containing a central spherical cavity. The kinematics proposed by [Hou, H.S., Abeyaratne, R., 1992. Cavitation in elastic and elastic-plastic solids. J. Mech. Phys. Solids 40, 571–722] was adopted in order to construct an approximate but optimal field. In order to establish a suitable constitutive law for this class of materials, we utilized the homogenisation technique that permits us to calculate the average strain energy density of the volume. The cavity growth was considered through a physically realistic failure criterion. Combination of the constitutive law and the failure criterion enables us to describe correctly the global behaviour and the damage evolution of the material under tri-axial loading. It was shown that the present models can efficiently reproduce different stress states, varying from uniaxial to tri-axial tensions, observed in experimentations. Comparison between predicted results and experimental data proves that the proposed model is accurate and physically reasonable. Another advantage is that the proposed model does not need special identification work, the initial Rivlin’s law for the corresponding incompressible material is sufficient to form the new law for the compressible material resulted from cavitation procedure.     

Calibration,validation, and verification including uncertainty of a physically motivated internal state variable plasticity and damage model

In this paper, we illustrate a formal calibration, validation, and verification process that includes uncertainty in an internal state variable plasticity-damage model that is implemented in a finite element code. The physically motivated continuum model characterizes damage evolution by incorporating material uncertainty due to microstructural spatial clustering. The uncertainty analysis is performed by introducing material variation through model validation and verification. The effect of variability in microstructural clustering and boundary conditions on the sensitivities and uncertainty of the plasticity-damage evolution for the 7075 aluminum alloy is characterized. The results show the potential of this methodology in the evaluation of material response uncertainty due to microstructure spatial clustering and its effect on damage evolution. For damage evolution, we have shown that the initial isotropic damage evolved into an anisotropic form as the deformation increased which is consistent with experimentally observed behavior for 7075 aluminum alloy in literature. Also, the sensitivities were found to be consistent with the physics of damage progression for this particular type of material. Through the sensitivity analysis, the initial defect size and number density of cracked particles are important at the beginning of deformation. As damage evolves, more voids are nucleated and grow and the sensitivity analysis illustrates this as well. Then, voids combine with each other and coalescence becomes the main driver, which is also confirmed by the sensitivity analysis. This work also shows that the microstructurally based damage evolution equations provide an accurate representation of the damage progression due to large intermetallic particles. Finally, we illustrate that the initial variation in the microstructure clustering can lead to about ±7.0%, ±8.1%, and ±9.75% variation in the elongation to failure strain for torsion, tensile, and compressive loading, respectively.