基于界面层模型的双周期纳米夹杂复合材料反平面问题研究

利用复变函数方法并结合双准周期Riemann边值问题理论,获得了含双周期分布非均匀相（夹杂/界面层）的复合材料在远场均匀反平面应力下弹性场的全场解答．该解答可用于对纳米夹杂复合材料的应力进行分析,结合平均场理论也用于预测纳米夹杂复合材料的有效性能．计算结果表明：当夹杂尺度在纳米量级时,应力和有效反平面剪切模量具有明显的尺度依赖性,并且随着夹杂尺寸的增加,趋近于不考虑界面效应时的结果;界面层厚度和性能对应力和有效反平面剪切模量明显变化时所对应的夹杂尺度范围和趋近于无界面效应结果的快慢有显著影响;当界面厚度足够薄时,界面层模型可用于模拟零厚度界面情况．     

圆形刚性夹杂双周期分布的反平面问题

研究含双周期分布的圆形刚性夹杂在无穷远受纵向剪切的弹性平面问题，遵循复合材料中各夹杂相互影响的重要条件。采用复变函数方法。构造相应模型的复应力函数。通过坐标变换，同时满足夹杂边界位移条件，再利用围线积分将求争方程组化为线性代数方程组。导出了圆形刚性夹杂双周期分布的界面应力解析表达式。算例给出了界面应力最大值与夹杂间距的变化规律。求出了刚性夹杂的合理间距问题，本文发展的分析方法为研究夹杂材料的细观机理探索了一条有效的分析途径。     

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

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

横观各向同性介质中椭圆夹杂受非弹性剪切变形引起的弹性场的确定

本文求解了横观各向同性介质中椭圆夹杂内受非弹性剪切变形引起的弹性场。采用各向异性弹性力学平面问题的复变函数解法,结合保角变换,获得夹杂内应变能和基体内边界的应力分布和相应的应变能的表达式。进一步,根据最小应变能原理,获得表征夹杂平衡边界的两个特征剪切应变,从而得到了弹性场的解析解。通过应力转换关系,验证了应力解满足夹杂边界上法向正应力和剪应力的连续条件,表明了该解的正确性。本文解可用于复合材料断裂强度的分析中。     

考虑夹杂相互作用的复合陶瓷夹杂界面的断裂分析

复合材料中夹杂含量较高时,夹杂间的相互作用能显著改变材料细观应力应变场分布,基体和夹杂中的平均应力应变水平也会发生较大变化,导致复合材料强度等力学性能发生显著变化. 为修正单一夹杂模型运用在实际材料中的误差,基于相互作用直推估计法,建立一种考虑含夹杂相互作用的夹杂界面裂纹开裂模型. 首先根据相互作用直推估计法,得到残余应力和外载应力共同作用下夹杂中的平均应力,再计算无限大基体中相同的夹杂达到相同应力场时的等效加载应力,将此加载应力作为含界面裂纹夹杂的等效应力边界条件,在此边界条件下求得界面裂纹尖端的应力强度因子,进而得到界面裂纹开裂的极限加载条件,并分析了夹杂弹性性能、含量、热残余应力、夹杂尺寸等因素对界面裂纹开裂条件的影响. 结果表明,方法能够有效修正单夹杂模型运用在实际材料中的误差,较大的残余应力对界面裂纹开裂有重要的影响,夹杂刚度的影响并非单调且比较复杂;在残余应力较小时,降低柔性夹杂刚度或者增大刚性夹杂刚度都有利于提高材料强度;扩大夹杂尺寸将导致裂纹开裂极限应力显著降低,从而降低材料强度.      

含刚性椭圆夹杂的弹性平面奇异解

运用复变函数方法,求解了含刚性椭圆夹杂的无限弹性平面在任意位置作用集中力和集中力偶的问题,导出了界面应力公式,绘出了应力分布曲线。     

圆内旋轮线形状异质夹杂热弹性问题的解析解

论文研究具有圆内旋轮线型形状夹杂域的平面热弹性体在夹杂域内非均匀温度场作用下对弹性场所产生的影响,其中考虑的夹杂与基体的材料不同但是具有相同的剪切模量。借助黎曼映射理论,将平面光滑闭合曲线外部区域映射到单位圆外部区域,进而利用解析函数性质,结合柯西型积分与Faber多项式,求解得到夹杂域内外场势函数的显式解。通过得到的势函数求出内外场应力,并对应力分布进行分析,结果表明：一般形状异质夹杂时,内场应力值与有限元计算值相吻合;退化到椭圆同质夹杂时,与相关文献中的结果相同,但是更具一般性与实际可操作性。     

多圆形刚性夹杂的反平面问题

构造任意分布且相互影响的多个圆形刚性夹杂模型的复应力函数，采用复变函数方法，达到满足各个夹杂的边界条件，利用坐标变换和围线积分将求解方程组化为线性代数方程组，推导出了圆形刚性夹杂任意分布的界面应力解析表达式，算例对多夹杂与单夹杂两种模型的界面应力最大值进行了对比，同时还给出了界面应力最大值随夹杂间距的变化规律，求出了刚性夹杂的合理间距。本文发展的分析方法为研究夹杂材料的细观机理探索了一条有效的分析途径。     

倾斜界面耦合弹性层中的渡越辐射能

物理学中,摄动源在非均匀介质中或非均匀介质附近匀速直线运动所产生的能量辐射现象称为渡越辐射.列车沿轨道运行,由轮轨接触产生的弹性波在非均匀轨道和基础中传播将发生渡越辐射,而轨道和基础的非均匀性集中体现在不同轨道基础之间的过渡段(如路桥过渡段、桥隧过渡段或有砟-无砟轨道过渡段).为研究车致弹性波在过渡段中引发的渡越辐射现象,本文以典型高速铁路路桥过渡段结构形式为依据,建立了二维平面应力渡越辐射能计算模型.其中,两个材料参数不同的半无限弹性层由一倾斜界面耦合,底端固定,上表面自由,一个集中载荷在自由表面上匀速运动.界面两侧弹性体中的波动方程均分解为本征场、自由场两个部分分别求解,其中自由场波动方程采用分离变量法数值求解.通过模型求解得到了不同载荷移动速度和界面倾斜角度条件下的渡越辐射能及界面附近应变能密度.结果表明,渡越辐射能的大小随载荷移动速度增大单调非线性增大,移动载荷速度达到刚度较大一侧介质表面波速的74%时产生的渡越辐射能就将超过载荷本身激发的本征场应变能;界面倾斜角度越大,即两侧介质刚度过渡距离越短,渡越辐射能与本征场应变能比值越大.     

正交弹性材料中双周期裂纹反平面问题的封闭解

研究了无限大正交弹性材料中含双周期裂纹的反平面问题,其基本胞元含有三条裂纹,且三条裂纹的中心恰好位于一等腰三角形顶点。运用椭圆函数、保角变换理论、施瓦兹公式获得了该问题应力场的封闭解,并得到了裂纹尖端处的应力强度因子。该问题结果取特殊情形退化对应于单排共线周期裂纹的解答。通过数值算例,分析了双周期裂纹归一化的应力强度因子随双周期裂纹的横向间距和纵向间距之比/b a 分别取10、5、2、1时的变化曲线。结果表明：对于一定的横向间距,应力奇异因子随纵向间距的增大而减小,但随着纵向间距的增大,纵向间距对应力奇异因子的影响变得不明显;对于一定的纵向间距,应力奇异因子随横向间距的减小而减小,但随着横向间距的减小,横向间距对应力奇异因子的影响变得不明显。     

Interfacial Energies for Incoherent Inclusions

We study a variational problem describing an incoherent interface between a rigid inclusion and a linearly elastic matrix. The elastic material is allowed to slip relative to the inclusion along the interface, and the resulting mismatch is penalized by an interfacial energy term that depends on the surface gradient of the relative displacement. The competition between the elastic and interfacial energies induces a threshold effect when the interfacial energy density is non-smooth: small inclusions are coherent (no mismatch); sufficiently large inclusions are incoherent. We also show that the relaxation of the energy functional can be written as the sum of the bulk elastic energy functional and the tangential quasiconvex envelope of the interfacial energy functional.     

On multiple circular inclusions in plane magnetoelasticity

A general series solution to the magnetoelastic problem of interacting circular inclusions in plane magnetoelasticity is provided in this paper. By the use of complex variable theory and Laurent series expansion method, the general expression of the magnetic and the magnetoelastic complex potentials for the circular inclusion problem is derived. Expanding the definition of the Airy’s stress function of pure elastic field into the magnetoelastic field and applying the superposition method, the general expression then can be reduced to a set of linear algebraic equations and solved in a series form. An approximate closed form solution for the case of two arbitrarily located inclusions is also provided. For illustrating the effect of the pertinent parameters, the numerical results of the interfacial magnetoelastic stresses are displayed in graphic form.     

Intensity of singular stress fields causing interfacial debonding at the end of a fiber under pullout force and transverse tension

In this study, singular stress fields at the ends of fibers are discussed by the use of models of rectangular and cylindrical inclusions in a semi-infinite body under pullout force. Those singular stresses have not been discussed yet in the previous studies for pullout problems although they are important for causing interfacial initial debonding. The body force method is used to formulate those problems as a system of singular integral equations where unknowns are densities of the body forces distributed in a semi-infinite body having the same elastic constants as those of the matrix and inclusions. In order to compare the results with the previous solutions, tension problems of a fiber in a semi-infinite body are also considered. Then, generalized stress intensity factors at the corner of rectangular and cylindrical inclusions are systematically calculated for various geometrical conditions with varying the elastic ratio, length, and spacing of the location from edge to inner of the body. The effects of elastic modulus ratio and aspect ratio of inclusion upon the stress intensity factors are discussed for pullout problems.     

Stress concentrations in the particulate composite with transversely isotropic phases

The accurate series solution have been obtained of the elasticity theory problem for a transversely isotropic solid containing a finite or infinite periodic array of anisotropic spherical inclusions. The method of solution has been developed based on the multipole expansion technique. The basic idea of method consists in expansion the displacement vector into a series over the set of vectorial functions satisfying the governing equations of elastic equilibrium. The re-expansion formulae derived for these functions provide exact satisfaction of the interfacial boundary conditions. As a result, the primary spatial boundary-value problem is reduced to an infinite set of linear algebraic equations. The method has been applied systematically to solve for three models of composite, namely a single inclusion, a finite array of inclusions and an infinite periodic array of inclusions, respectively, embedded in a transversely isotropic solid. The numerical results are presented demonstrating that elastic properties mismatch, anisotropy degree, orientation of the anisotropy axes and interactions between the inclusions can produce significant local stress concentration and, thus, affect greatly the overall elastic behavior of composite.     

Analysis of stress concentration around a spheroidal cavity under asymmetric dynamic loading

The fracture and fatigue properties of porous materials are strongly influenced by stress concentrations around the pores. In addition, failure of structural components initiates at locations of high stress concentration which is often caused by holes, inclusions or other discontinuities. In view of this, the stress concentration around a spheroidal cavity embedded in an elastic medium is studied under dynamic loading conditions. While solutions abound for static loads, only limited solutions exist for dynamic loads. The stress field around a spheroidal cavity is determined by using a hybrid methodology that combines the finite element technique with a spherical wave function expansion method. The stress concentrations within the matrix are found to be dependent on the frequency of excitation, aspect ratio of the cavity and the Poisson’s ratio of the matrix. The study reveals that dynamic stress concentrations can reach much higher values than those encountered under static loading.     

Spherical ceramic pebbles subjected to multiple non-concentrated surface loads

This paper presents an analytical solution for the stress distributions within spherical ceramic pebbles subjected to multiple surface loads along different directions. The method of solution employs a displacement approach together with the Fourier associated Legendre expansion for piecewise boundary loads. The solution corresponds to spherically isotropic elastic spheres. The classical solution for isotropic spheres subjected diametral point loads is recovered as a special case of our solution. For the isotropic pebbles under consideration, stresses within spheres are numerically evaluated. The results show that the number of loads does have significant influence on the maximum tensile stress inside the sphere. Moreover, the applicability of solutions using the series expansion method for stresses near surface load areas is also examined. The stresses evaluated with large enough number of terms agree quite well with those derived from FEM simulations, except around the edge of circle load area.     

Effect of geometrical and material parameters in nanoindentation of layered materials with an interphase

Indentation models for thin layer-substrate geometry with an interphase have been developed. The interphase can be modeled either as a nonhomogeneous layer or as a homogeneous layer. Between the two models of the interphase, contact depth and critical interfacial stresses are compared to find the effect of indentation area, film and substrate Young’s moduli, and the interphase and film thicknesses. Although contact depth is found not to be sensitive to the type of interphase model used, critical interfacial stresses are significantly different (up to 15%) for film to substrate elastic Young’s moduli ratios of more than 25. A formal sensitivity analysis based on design of experiments shows that on critical interfacial stresses, interphase to film thickness ratio and film to substrate Young’s moduli ratio has the most impact, while type of elastic moduli variation in the interphase and indentor width to film thickness ratio has the least impact.     

含空心球涂层粒子复合材料界面的脱粘应力分析

在球对称拉伸载荷作用下针对空心球涂层复合材料分析了空心球涂层粒子增强复合材料的局部应力场,得到了界面临界脱粘应力的解析表达式.讨论了各相几何参数对非均匀涂层空心球粒子临界脱粘应力的影响,比较了均匀涂层和非均匀涂层的脱粘应力.结果表明:在球对称拉伸下界面脱粘更容易发生在涂层相与基体相界面间,空心球的壁厚和涂层厚度是影响界面临界脱粘应力的重要因素,因而选择适当的空心球、涂层厚度和提高界面粘结能将有利于提高界面的临界脱粘应力.     

Hamilton体系下矩形薄板受抛物线压力载荷的屈曲分析

针对四边简支矩形薄板在两对边相向的非线性分布压力下的面内应力分布以及屈曲问题,应用弹性力学的Hamilton体系和Galerkin法进行了研究.基于弹性力学的平面矩形域Hamilton体系,根据辛本征向量展开解法,得到了对应于零本征值和非零本征值的含待定常数的实数型面内应力分布通解.依据必须满足的应力边界条件,导出了矩形薄板在抛物线分布载荷下的面内应力分布.考虑到应力分布表达式的复杂性,用完全的解析方法得到屈曲载荷是不可能的.因此,运用基于虚功原理的Galerkin法,根据四边简支矩形薄板弯曲的位移边界条件,给出了不同长宽比矩形薄板受抛物线分布载荷的屈曲临界载荷.通过与已有文献中DQ法给出的数值计算结果比较,表明了本文求解方法的有效性和正确性.基于所给出的结果,可望为解决矩形薄板在非线性分布载荷下的面内应力分布以及屈曲问题提供一种新的研究方法.