主裂纹与微裂纹相互作用问题的有效数值计算方法

提出了平面弹性介质中主裂纹与微裂纹相互作用问题的有效数值计算 方法. 通过把适于单一裂纹的Bueckner原理扩充到含有多裂纹的一般体系，将原问题分解 为承受远处载荷不含裂纹的均匀问题，和在远处不承受载荷但在裂纹面上承受面力的多裂纹 问题. 于是，以应力强度因子作为参量的问题可以通过考虑后者(多裂纹问题)来解决，而 利用提出的杂交位移不连续法，这种多裂纹问题是容易数值求解的. 列举 Cai和 Faber为评价主裂纹与微裂纹相互作用问题的近似方法而列举的算例，说明 该数值方法对分析平面弹性介质中主裂纹与微裂纹相互作用问题既简单又非常有效.     

带微裂纹物体的有效断裂韧性

按照等效介质的思想，引进有效表面能密度的概念，建立了带微裂纹物体有效断裂韧性的公式．具体计算了微裂纹群分别平行和垂直于宏观裂纹两种情况的减韧比．表明微裂纹群在产生应力屏蔽（或反屏蔽）效应的同时，也降低了材料的有效断裂韧性，减小了对宏观裂纹的扩展阻力．     

二维微裂纹体有效性质之间的关联研究

通过解析和数值的方法研究了二维微裂纹体有效模量和有效导电系数之间的关联，针对二维含有任意取向和平行的微裂纹介质，通过基于象素的有限元方法分别计算了它们的有效杨氏模量和有效导电系数，并建立了两种不同物理量之间的关联．研究结果表明利用简单的解析Mori-Tanaka方法也可以建立上述关联，并且与数值结果吻合很好。     

损伤体有效弹性性质的细观分析和不变性描述—— 一个考虑微缺陷相互作用的一般理论模式

通过引进微缺陷相互作用张量，建立了一个二维情况下考虑微缺陷（微裂纹或微孔洞）间相互作用的损伤固体有效弹性性质的一般理论模式模型中考虑了微缺陷的几何形状、取向分布和空间分布所造成的有效柔度张量的各向异性和材料中微缺陷之间的相互作用所引起的损伤柔度张量的高阶效应针对微椭圆孔、微圆孔和微裂纹问题，求得了相互作用张量的解析形式     

损伤体有效弹性性质的细观分析和不变性描述—— 一个考虑微缺陷相互作用的一般理论模式

通过引进微缺陷相互作用张量，建立了一个二维情况下考虑微缺陷（微裂纹或微孔洞）间相互作用的损伤固体有效弹性性质的一般理论模式模型中考虑了微缺陷的几何形状、取向分布和空间分布所造成的有效柔度张量的各向异性和材料中微缺陷之间的相互作用所引起的损伤柔度张量的高阶效应针对微椭圆孔、微圆孔和微裂纹问题，求得了相互作用张量的解析形式     

PZT-4紧凑拉伸试样的断裂分析

基于线性压电材料的复势理论，通过解析分析，导出了一种分析有限压电板裂纹问题的 解析数值方法. 首先，计算了含中心裂纹有限板的断裂参数，与Woo和Wang的解析数值 法(Int J Fract, 1993, 62: 203$\sim$218)相比较，表明该方法具有很高的精度和 很好的计算效率. 随后，采用该方法和有限元 法计算了PZT-4紧凑拉伸试样在绝缘裂纹面边界条件下断裂时的断裂参数，发现各断裂参数 的临界值分散性很大，不能作为压电材料的单参数断裂准则. 进而，针对试样真实的裂隙形 状，采用有限元法计算了裂隙尖端的应力、电位移场，比较了裂隙内介质的介电性能对裂隙 尖端场的影响，计算了带微裂纹的真实裂隙模型的断裂参数并进行了理论分析.     

多相材料有效性能预测的高精度方法

有效介质方法是常用的细观力学方法之一.其可用于计算多相材料的有效性能,并建立材料微细观结构和宏观性能的定量关系;有助于指导新材料设计,减少试验工作量等.然而,当夹杂含量升高时,传统有效介质方法的计算精度下降.本文以两相材料为研究对象,提出一种新的参考介质,即:为更合理考虑不同夹杂颗粒间的相互作用,假定参考介质的应变是基体相平均应变和某一修正张量的双点积.在此基础上,推导了新参考介质下两相材料的有效模量表达式,并给出该修正张量的近似计算方法;通过反复更新参考介质,采用多层次均匀化思路,将本文方法进一步用于多相材料性能的预测.为验证方法的有效性,将预测结果与已有模型结果和试验数据进行对比.结果表明本文方法较已有方法更为合理、有效.当夹杂含量升高时,本文方法较传统有效介质方法的计算精度有所提升.     

三维弹性问题Taylor展开多极边界元法的误差分析

Taylor展开多极边界元法有效的提高了边界元法的求解效率,使之可用于大规模问题的计算。然而,由于计算中对基本解进行了Taylor级数展开,与传统边界元方法相比计算精度有所下降。本文主要针对三维弹性问题Taylor展开多极边界元法的计算精度和误差进行研究。文中对两种方法的计算精度进行了比较;研究了核函数的Taylor展开性质;推导了三维弹性问题基本解的误差估计公式;给出了Taylor展开多极边界元法中远近场的划分原则。通过具体的算例,证明了该方法的正确性和误差估计公式的有效性,说明了影响Taylor展开多极边界元法求解精度的因素。     

粉末冶金材料各向异性连续型损伤模型研究

在采用微裂纹扩展区描述脆性材料损伤状态的基础上,通过引入颗粒度大小和微裂纹概率分布密度函数来研究粉末冶金材料特征体积单元的损伤演化.同时结合Hill自洽模型采用空间混合平均的方法,可考虑粉末冶金材料的微结构特征及微裂纹间相互作用的影响,可描述粉末冶金类材料的各向异性连续型损伤.     

多孔介质中流固作用力的动量交换计算

流经多孔介质的流动广泛存在于化工生产、多孔颗粒悬浮流等领域,如何准确计算流体与多孔介质之间的作用力是研究此类流动的一个关键因素.作为一种有效的流体计算方法,格子波尔兹曼方法(latticeBoltzmann method,LBM)常采用动量交换法计算流体与固体之间相互作用力.分析了流体流经多孔介质时两者的动量交换过程,提出了一种高效的动量交换法来计算流固作用力,并在孔隙尺度下对其进行了验证,结果表明该方法是可行的.进而将该方法推广到计算表征体元(表征体元)尺度下的流--固相互作用,并对不同雷诺数(Re)下的多孔方柱绕流问题进行了模拟和验证.     

A simple method for calculating interaction of numerous microcracks and its applications

The effects of microcrack interaction on the failure behavior of materials present one problem of considerable interest in micromechanics, which has been extensively argued but has not been resolved as yet. In the present paper, a simple and effective method is presented based on the concept of the effective field to analyze the interaction of microcracks of a large number or of a high density. To determine the stress intensity factors of a microcrack embedded in a solid containing numerous or even countless microcracks, the solid is divided into two regions. The interaction of microcracks in a circular or elliptical region around the considered microcrack is calculated directly by using Kachanov’s micromechanics method, while the influence of all other microcracks is reflected by modifying the stress applied in the far field. Both the cases of tensile and compressive loading are considered. This simplified scheme may yield an estimate for stress intensity factors of satisfactory accuracy, and therefore provide a potential tool for elucidating some phenomena of material failure associated with microcracking. As two of its various promising applications, the above scheme is employed to investigate the size effects of material strength due to stochastic distribution of interacting microcracks and to calculate the effective elastic moduli of elastic solids containing distributed microcracks. Some conventional micromechanics methods for estimating the effective moduli of microcracked materials are evaluated by comparing with the numerical results. Only two-dimensional problems have been considered here, though the three-dimensional extension of the present method is of greater interest.     

Estimate of effective elastic moduli with microcrack interaction effects

A method is developed for estimating the effects of microdefect interaction on the effective elastic properties of heterogeneous solids. An effective medium is defined to calculate the global effective elastic moduli of brittle materials weakened by distributed microcracks. Each microcrack is assumed to be embedded in an effective medium, the compliance of which is obtained from the dilute concentration method without accounting for interaction. The present scheme requires no iteration; it can account for microcrack interaction with sufficient accuracy. Analytical solutions are given for several two- and three-dimension problems with and without anisotropy.     

Microcrack interaction brittle rock subjected to uniaxial tensile loads

A micromechanics-based model is proposed to describe unstable damage evolution in microcrack-weakened brittle rock material. The influence of all microcracks with different sizes and orientations are introduced into the constitutive relation by using the statistical average method. Effects of microcrack interaction on the complete stress–strain relation as well as the localization of damage for microcrack-weakened brittle rock material are analyzed by using effective medium method. Each microcrack is assumed to be embedded in an approximate effective medium that is weakened by uniformly distributed microcracks of the statistically-averaged length depending on the actual damage state. The elastic moduli of the approximate effective medium can be determined by using the dilute distribution method. Micromechanical kinetic equations for stable and unstable growth characterizing the ‘process domains’ of active microcracks are taken into account. These ‘process domains’ together with ‘open microcrack domains’ completely determine the integration domains of ensemble averaged constitutive equations relating macro-strain and macro-stress. Theoretical predictions have shown to be consistent with the experimental results.     

Constitutive relationship of brittle rock subjected to dynamic uniaxial tensile loads with microcrack interaction effects

A micromechanical model is proposed to describe both stable and unstable damage evolution in microcrack-weakened brittle rock material subjected to dynamic uniaxial tensile loads. The basic idea of the present model is to classify the constitution relationship of rock material subjected to dynamic uniaxial tensile loads into four stages including some of the stages of linear elasticity, pre-peak nonlinear hardening, rapid stress drop, and strain softening, and to investigate their corresponding micromechanical damage mechanisms individually. Special attention is paid to the transition from structure rearrangements on microscale to the macroscopic inelastic strain, to the transition from distribution damage to localization of damage and the transition from homogeneous deformation to localization of deformation. The influence of all microcracks with different sizes and orientations are introduced into the constitutive relation by using the statistical average method. Effects of microcrack interaction on the complete stress-strain relation as well as the localization of damage for microcrack-weakened brittle rock material are analyzed by using effective medium method. Each microcrack is assumed to be embedded in an approximate effective medium that is weakened by uniformly distributed microcracks of the statistically-averaged length depending on the actual damage state. The elastic moduli of the approximate effective medium can be determined by using the dilute distribution method. Micromechanical kinetic equations for stable and unstable growth characterizing the ‘process domains’ of active microcracks are taken into account. These ‘process domains’ together with ‘open microcrack domains’ completely determine the integration domains of ensemble averaged constitutive equations relating macro-strain and macro-stress. Theoretical predictions have shown to consistent with the experimental results.     

On collinear microcrack–inclusion arrays interaction and related problems

Investigated is a crack problem for an array of collinear microcracks in composite matrix. Inclusions are situated in between the neighbouring microcracks tips and exhibit different elastic properties than matrix. The problem is solved using the technique of distributed dislocations. A developed approximate fundamental solution for a single dislocation lying in a general point between inclusions is employed in the distribution of continuously distributed dislocation to cracks modelling. Stress intensity factor is calculated for various cracks/inclusions geometries and elastic moduli mismatches. Stability and/or instability of the straight microcrack paths is investigated for slowly growing microcracks with inclusions located in between the neighbouring microcracks tips. Applications to periodic microcrack tunnelling and microcracks weakening ahead of the main crack are discussed.     

A unified energy approach to a class of micromechanics models for composite materials

Several micromechanics models for the determination of composite moduli are investigated in this paper, including the dilute solution, self-consistent method, generalized self-consistent method, and Mori-Tanaka's method. These micromechanical models have been developed by following quite different approaches and physical interpretations. It is shown that all the micromechanics models share a common ground, the generalized Budiansky's energy-equivalence framework. The difference among the various models is shown to be the way in which the average strain of the inclusion phase is evaluated. As a bonus of this theoretical development, the asymmetry suffered in Mori-Tanaka's method can be circumvented and the applicability of the generalized self-consistent method can be extended to materials containing microcracks, multiphase inclusions, non-spherical inclusions, or non-cylindrical inclusions. The relevance to the differential method, double-inclusion model, and Hashin-Shtrikman bounds is also discussed. The application of these micromechanics models to particulate-reinforced composites and microcracked solids is reviewed and some new results are presented.     

Electroelastic behavior modeling of piezoelectric composite materials containing spatially oriented reinforcements

In this work, a modeling of electroelastic composite materials is proposed. The extension of the heterogeneous inclusion problem of Eshelby for elastic to electroelastic behavior is formulated in terms of four interaction tensors related to Eshelby’s electroelastic tensors. Analytical formulations of interaction tensors are presented for ellipsoidal inclusions. These tensors are basically used to derive the self-consistent model, Mori–Tanaka and dilute approaches. Numerical solutions are based on numerical computations of these tensors for various types of inclusions. Using the obtained results, effective electroelastic moduli of piezoelectric multiphase composites are investigated by an iterative procedure in the context of self-consistent scheme. Generalised Mori–Tanaka’s model and dilute approach are re-formulated and the three models are deeply analysed. Concentration tensors corresponding to each model are presented and relationships of effective coefficients are given. Numerical results of effective electroelastic moduli are presented for various types of piezoelectric inclusions and for various orientations and compared to existing experimental and theoretical ones.     

Prediction of the overall moduli of a cylindrical short-fiber reinforced composite

With respect to obtaining the effective elastic moduli of the composite, the present theory differs from both Eshelby's equivalent inclusion method and Hill's self-consistent one, both of which only consider the mechanical properties of the matrix and inclusions (fibers). In fact, the inclusion-inclusion interaction is more pronounced when the volume fraction of inclusions of the composite increases. Hence, in this paper the effective elastic moduli of the composite are derived by taking into account the shapes, sizes and distribution of inclusions, and the interactions between inclusions. In addition, it is more convincing to assume short-fibers as cylindrical inclusions as in the present paper than as ellipsoidal ones as in others^[7,8]. Finally, numerical results are given.     

A UNIFIED ENERGY APPROACH TO A CLASS OF MICROMECHANICS MODELS FOR MICROCRACKED SOLIDS

Micromechanics models have been developed for the determination ofthe elastic moduli of microcracked solids based on different approaches andinterpretations,including the dilute or non-interacting solution,the Mori-Tanakamethod,the self-consistent method,and the generalized self-consistent method.It isshown in the present study that all these micromechanics models can be unified withinan energy-equivalence framework,and that they differ only in the way in which themicrocrack opening and sliding displacements are evaluated.Relevance to thedifferential methods and the verification of these models are discussed.