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

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

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

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

径向非均匀介质中圆形夹杂的动应力分析

基于复变函数理论,研究了径向非均匀弹性介质中均匀圆夹杂对弹性波的散射问题. 介质的非均匀性体现在介质密度沿着径向按幂函数形式变化且剪切模量是常数. 利用坐标变换法将变系数的非均匀波动方程转为标准亥姆霍兹(Helmholtz) 方程. 在复坐标系下求得非均匀基体和均匀夹杂同时存在的位移和应力表达式. 通过具体算例分析了圆夹杂周边的动应力集中系数(DSCF). 结果表明:基体与夹杂的波数比和剪切模量比,基体的参考波数和非均匀参数对动应力集中有较大的影响.      

层状压电半空间中的螺位错

研究了层状压电半空间中的螺位错与界面的相互作用,应用傅立叶变换的方法,推导了压电层和基体的弹性场和电场,并在此基础上给出了像力的表达式,分析了几何尺寸和材料的相对刚度对像力和界面应力的影响.结果表明:压电效应的存在可以改变位错与界面的相互作用机理;弹性常数的差异对界面应力有较大的影响.     

宏微观渐进展开损伤本构模型及其在碳纤维增强复合材料中的应用

基体开裂、纤维拔出、界面剥离等是碳纤维增强复合材料常出现的局部各向异性损伤现象,这些损伤逐渐扩展,削弱了材料的强度和刚度,影响材料的承载能力.对此利用宏微观摄动理论对位移进行双范围渐进展开,在微观位移中引入损伤应变,通过计算损伤应变集中因子,得到了含损伤的均质化损伤弹性常数(宏观有效刚度矩阵),用平均法和混合法检验了无...     

平行微裂纹削弱的弹性体的各向异性有效性能与双周期裂纹模型

平行微裂纹损伤模型被用于构建各向异性损伤理论.当施加在代表性体积单元上的边界条件满足Hill条件时,基于平均场理论论证了由平行穿透裂纹损伤的弹性体仅有6个独立有效弹性常数.除了原各向同性基体的2个弹性常数外,与损伤相关的另外4个常数中,3个描述有效弹性常数的折减,1个描述损伤导致的拉剪耦合效应.结合单胞模型和有限元方法分析了双周期阵列平行裂纹问题,数值结果显示:裂纹呈一般双周期阵列时,拉剪耦合参数相比其它模量小很多;当裂纹密度一定时,改变裂纹的排列形式,面内剪切模量和面外剪切模量的折减呈现出不同的规律.     

混凝土微观界面纳米划痕试验表征研究

鉴于混凝土微观结构对宏观性能的重要影响,利用纳米划痕试验方法,对常规水灰比配比的混凝土样品骨料-水泥浆体微观界面及界面过渡区进行了表征研究．基于纳米划痕试验方法和扫描电镜显微技术,对于样品微观情况下的骨料-水泥浆体界面进行了表征分析,并对表征结果经过多个高斯分布拟合的统计学分析方法得到基于界面位置的结构划痕硬度分布情况．通过划痕硬度与骨料-浆体界面的距离变化,定义了骨料和水泥浆体之间的界面过渡区,并得到界面过渡区厚度约为40μm,界面过渡区划痕硬度约为水泥浆料硬度的75％左右．通过比较划痕断裂韧性与宏观试验结果,证明利用纳米划痕试验可以初步分析材料断裂韧性．本研究为更好地研究混凝土微观性能对宏观力学性能的影响提供了可靠的数据支持,也为评价混凝土微观界面力学性能提供了一种新的可行方法．     

反平面剪切下电磁弹性纤维增强复合材料界面相效应研究

研究具有界面相电磁弹性纤维增强复合材料的反平面剪切问题,利用复变函数方法,获得了无穷域中带界面相纤维问题在远场力、电、磁多场作用下的闭合解,得到了复合材料内部各区域电磁弹性物理量的精确表达式.利用所得结果,考虑纤维和基体间的界面相效应,研究了界面相厚度及弹性模量对复合材料内部应力场、电场强度和磁场强度的影响,数值结果给出了复合材料电磁弹性物理量随界面相参数变化的规律,为该类复合材料的设计与计算提供了有价值的参考.     

含界面相桥联矩阵

桥联模型能够有效分析复合材料的弹塑性性能,其在World-Wide Failure Exercise(WWFE)评比中被认定为精度最高的细观力学理论.桥联模型的核心是桥联矩阵,现有桥联模型是在两相复合材料理想界面条件下建立的,而对于实际的复合材料而言,纤维和基体之间总或多或少的出现界面非理想现象.学者们常用一个存在于基体和纤维之间的界面相来描述这种非理想界面情形,这时就需要建立含界面相的三相桥联矩阵.论文采用三相CCA模型,求得六种不同边界载荷下的应力场,再将各相体积平均后的应力场代入三相桥联模型方程,求解线性方程组即可得到三相桥联矩阵元素的精确表达式.进一步,将纤维相和界面相看作一个局部复合材料或等效横观各向同性纤维,其力学性能利用两相桥联模型确定,与外部的基体相构成的复合材料依然可以通过两相桥联模型分析.利用这种等效纤维的概念,论文将三相隐式桥联矩阵进行了有效简化,通过对复合材料等效弹性常数的数值计算,证明简化后的三相显式桥联矩阵与精确的三相隐式桥联矩阵有极高的近似度,从而给实际应用带来极大方便.     

压电多层材料中的电极-陶瓷界面裂纹和椭圆夹杂

首先解析研究了压电多层材料中的电极-陶瓷界面裂纹,它是对Ru最近研究工作（见ASMB J·APPI.Mech.,2000年,67卷）的补充和完善.工作表明与 Ru的结论不同,对于一般的两相压电介质,仍可获得对于这种电极-陶瓷界面裂纹的精确解答.分析显示在界面处的电弹性场仍可显现出两类奇异性:振荡型奇异性一1/2土iε。和实指数型奇异性一1/2士 k,其中ε和k 由上下两相压电材料的本构常数加以确定.获得了界面电弹性场以及裂尖处能量释放率的显式和实形式解答.也讨论了在坐标变换下界面上物理量的一些不变性质.其次探讨了压电多层材料中的椭圆夹杂问题,并获得了当压电复合系统受到远场均布机电载荷时的通解.分析表明当压电基体受到远场均布机电载荷作用时,应力场和电场在压电夹杂体上仍然均匀分布井且整个夹杂体都为等电势体     

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.     

基于纳米压痕试验的混凝土微观耐久性表征

利用纳米压痕和扫描电镜研究了不同冻融循环下混凝土的微观性能,揭示了冻融循环所导致混凝土退化的机理．纳米压痕结果显示,混凝土浆体的压痕模量与硬度的相关系数很大,并且都不满足正态分布．本文利用核密度估计来计算浆体微观性能的概率密度函数．根据浆体和骨料的性能差异,用阶跃函数来拟合界面过渡区的微观力学性能,同时利用扫描电镜观察界面过渡区的形貌．采用以上两种方法研究不同冻融循环下的混凝土,发现在不同的尺度下,冻融循环对混凝土的影响不同：在纳米尺度下,由于水的参与,高密度水凝胶会进一步形成;而在微米尺度下,由于冻融循环中水结冰膨胀和材料多相性,微裂纹进一步扩展,导致断裂能量降低,混凝土表现出宏观退化．     

Multiscale lattice Boltzmann-finite element modelling of chloride diffusivity in cementitious materials. Part II: Simulation results and validation

Chloride diffusivity in cementitious materials depends on the underlying microstructure and environmental conditions. The algorithms and implementation of the multiscale lattice Boltzmann-finite element scheme for prediction of chloride diffusivity in cementitious materials was described in detail in Part I (Zhang et al., 2013). Based on the obtained microstructures and the developed multiscale modelling scheme, chloride diffusivity in cementitious materials at the micro- and meso-scales, i.e. cement paste, mortar and concrete, are estimated and presented in Part II. The influences of w/c ratio, age, chloride binding, degree of water saturation, interfacial transition zone (ITZ) and aggregate content on chloride diffusivity are investigated in a quantitative manner. The simulations are validated with experimental data obtained from literature. The results indicate that the simulated chloride diffusivity in cementitious materials at each scale shows a good agreement with experimental data. In addition, the chloride binding, degree of water saturation, ITZ and aggregate content play significant roles in the chloride diffusivity in cementitious materials. The estimated chloride diffusivity in cementitious materials in this study accounting for the evolution of microstructure and environmental conditions can be directly used as input for the service life prediction of reinforced concrete structures.     

钢筋锈胀引发混凝土保护层开裂破坏的细观数值研究

钢筋锈蚀膨胀引起保护层混凝土开裂是影响钢筋混凝土结构耐久性和服役寿命的重要因素。考虑到混凝土细观结构组成对保护层破坏模式的影响,从细观角度出发,将混凝土看作由骨料、砂浆基质及两者间界面过渡区组成的三相复合材料,建立了描述钢筋锈胀力学行为的混凝土随机骨料模型。采用塑性损伤本构关系模型来表征砂浆基质和过渡区界面的力学行为,假定钢筋均匀锈蚀,对钢筋锈胀引起的混凝土保护层开裂破坏过程进行了细观数值研究。对比了宏观均匀模型与细观非均质模型下获得的保护层破坏模式,探讨了径厚比(c/d)、钢筋位置(中部和角区)及混凝土拉伸强度对保护层破坏模式及保护层胀裂时钢筋锈蚀水平的影响,得到了一些有益结论。     

A micromechanical model for effective elastic properties of particulate composites with imperfect interfacial bonds

A micromechanical model for effective elastic properties of particle filled acrylic composites with imperfect interfacial bonds is proposed. The constituents are treated as three distinct phases, consisting of agglomerate of particles, bulk matrix and interfacial transition zone around the agglomerate. The influence of the interfacial transition zone on the overall mechanical behavior of composites is studies analytically and experimentally. Test data on particle filled acrylic composites with three different interfacial properties are also presented. The comparison of analytical simulation with experimental data demonstrated the validity of the proposed micromechanical model with imperfect interface. Both the experimental results and analytical prediction show that interfacial conditions have great influence on the elastic properties of particle filled acrylic composites.     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Micromechanical models to guide the development of synthetic ‘brick and mortar’ composites

This paper describes a micromechanical analysis of the uniaxial response of composites comprising elastic platelets (bricks) bonded together with thin elastic perfectly plastic layers (mortar). The model yields closed-form results for the spatial variation of displacements in the bricks as a function of constituent properties, which can be used to calculate the effective properties of the composite, including elastic modulus, strength and work-to-failure. Regime maps are presented which indicate critical stresses for failure of the bricks and mortar as a function of constituent properties and brick architecture. The solution illustrates trade-offs between elastic modulus, strength and dissipated work that are a result of transitions between various failure mechanisms associated with brick rupture and rupture of the interfaces. Detailed scaling relationships are presented with the goal of providing material developers with a straightforward means to identify synthesis targets that balance competing mechanical behaviors and optimize material response. Ashby maps are presented to compare potential brick and mortar composites with existing materials, and identify future directions for material development.     

ANALYSIS OF MECHANICAL PROPERTIES OF CEMENT CLINKER UNDER DIFFERENT TEMPERATURES BASED ON HERTZIAN INDENTATION METHOD1)

本文通过高温球压法得到各种脆性和准脆性材料的表面压痕应力与应变之间的关系曲线。通过压痕应力-应变曲线的分析既可比较方便地确定出材料的压痕弹性模量、剪切模量和布氏硬度,又可比较不同温度下水泥熟料的变形性能。在不同温度(25 ${^\circ}\!$C$\sim $1400 ${^\circ}\!$C)处理下,球压应力松弛试验载荷松弛,在载荷峰值为100 N时,随着温度的升高,水泥熟料载荷松弛更明显。随着温度从500 ${^\circ}\!$C升高到1400 ${^\circ}\!$C,载荷松弛非常明显,尤其温度高于1200 ${^\circ}\!$C,水泥熟料样品内部的硅酸三钙(Ca$_3$SiO$_5$, 简称C3S)分解以及有部分液相的出现引起的应力松弛现象最为明显,在1275 ${^\circ}\!$C时熟料基本上已经软化,载荷急速松弛,所以认为1275 ${^\circ}\!$C为熟料的脆延转化温度。通过水泥熟料高温球压松弛试验可以确定水泥熟料在二次加热过程中的脆-延转化温度,测定熟料弹性模量和抗压强度急剧变化的温度范围。研究水泥熟料在不同温度下的力学行为和力学特性,探索提高粉磨效率的新途径,实现高温下的低能耗粉碎。     

Modelling of compressive behaviour of concrete-like materials at high strain rate

Uniaxial compression tests are the most common tests for characterizing the strength of concrete-like materials. The dynamic compression strength of concrete-like material is typically obtained by Split Hopkinson Pressure Bar (SHPB) tests. The increase in material strength under dynamic loading is usually attributed to the strain rate effect and modelled with a dynamic increase factor (DIF). However, it was observed by some researchers that the radial inertial confinement caused apparent increase of dynamic strength of concrete-like specimen in SHPB tests. They attributed the material strength increase to this inertial effect, instead of the strain rate effect. In the present study, numerical analyses are performed to investigate the compressive behaviour of concrete-like material at high strain rates. A homogeneous macroscale model and a heterogeneous mesoscale model are developed in the study. In the macroscale model, the material is assumed to be homogeneous and isotropic. In the mesoscale model, the test sample is modelled as a three-phase composite consisting of aggregate, mortar matrix and interfacial transaction zone (ITZ) between the aggregate and the mortar matrix. The aggregate is assumed to be circular and the ITZ is modelled as a thin boundary around the aggregate. In the both models, the materials are assumed to be insensitive to the strain rate first. Therefore, the obtained strength enhancement is only due to the inertial confinement. Strain rate sensitive material properties are then used in the two models in the calculations. Numerical simulations of the concrete samples under compression at different strain rates are carried out. The relative contribution of the inertial effect and the strain rate effect on the compressive strength DIF is examined based on the numerical results. The failure process of concrete specimen is also studied.     

Nonlinear analysis of compressed elastic thin films on elastic substrates: From wrinkling to buckle-delamination

Nonlinear buckling of elastic thin films on compliant substrates is studied by modeling and simulations to reveal the roles of pre-strain, elastic modulus ratio, and interfacial properties in morphological transition from wrinkles to buckle-delamination blisters. The model integrates an interfacial cohesive zone model with the Föppl–von Kármán plate theory and Green function method within the general framework of energy minimization. A kinetics approach is developed for numerical simulations. Subject to a uniaxial pre-strain, the numerical simulations confirm the analytically predicted critical conditions for onset of wrinkling and wrinkle-induced delamination, with which a phase diagram is constructed. It is found that, with increasing pre-strain, the equilibrium configuration evolves from flat to wrinkles, to concomitant wrinkles and buckle-delamination, and to an array of parallel straight blisters. The height and width of the buckle-delamination blisters can be approximately described by a set of scaling laws with respect to the pre-strain and interfacial toughness. Subject to an equi-biaxial pre-strain, the critical conditions are determined numerically to construct a similar phase diagram for the buckling modes. Moreover, by varying the pre-strain, modulus ratio, and interfacial toughness, a rich variety of equilibrium configurations are simulated, including straight blisters, and network blisters with or without wrinkles. These results provide considerable insight into diverse surface patterns in layered material systems as a result of the mechanical interactions between the film and the substrate through their interface, which suggests potential control parameters for designing specific surface patterns.