尖锐夹杂角端部局部应力场杂交有限元分析

本文首先利用作者曾提出的一维有限元特征分析方法计算所得到的尖锐夹杂角端部应力奇异指数和奇异应力场、位移场角分布函数,并依据Hellinger-Reissner原理,开发出了一个特殊的、能够反映夹杂角端部局部弹性现象的n结点多边形超级角端部单元,然后将该超级单元与标准的4结点杂交应力单元耦合在一起构建了一种分析异形夹杂角端部奇异弹性场的新型特殊杂交应力有限元方法.文中给出了两个应用算例,算例结果表明:本文方法不仅使用单元少、计算结果精度高,而且适用范围广,可拓展应用于分析复合材料微结构组织与力学行为关系.     

复合材料中矩形夹杂角端部力学行为分析

提出了一种分析矩形夹杂角端部奇异应力场的新型杂交有限元方法,该方法在分析矩形夹杂角端部奇异应力场时,需要在夹杂端部构造一个超级单元。超级单元的刚度矩阵可以通过夹杂端部特征问题数值解建立。我们用这种方法计算了单向载荷作用下无限大均质板中单个矩形夹杂角端部奇异应力场,并与现有的数值解进行了比较。比较结果表明:本文提出的方法是可行的、有效的,而且数值结果精度高。为说明本文方法适用范围更广,文章最后讨论了各向异性弹性材料和横观各向同性压电材料中矩形夹杂角端部电弹性场行为。     

夹杂角端部奇异应力场分析

提出一种分析夹杂角端部奇异应力场的新型杂交有限元方法.构造了一个角端部奇异单元,该单元刚度建立不依赖数学解析解.用这种方法计算了单向载荷作用下无限大板含单个方形夹杂和菱形夹杂角端部奇异应力场,并与现有的数值解进行了比较,结果表明:目前的数值方法是可行的、有效的、数值结果精度高,适用范围广.作为应用讨论了双方形夹杂刚度和位置对夹杂角端部奇异应力场的影响.     

周期分布多边形夹杂奇异性应力干涉问题的研究

采用一种新型的杂交元模型和一种单胞模型来解决周期分布多边形夹杂角部的奇异性应力相互干涉的问题。新型杂交元模型是基于广义Hellinger-Reissner变分原理建立的,其中奇异性应力场分量和位移场分量是采用有限元特征分析法的数值特征解得到的。使用当前的新型杂交元模型,只需要在夹杂角部邻域的周界上划分一维单元,避免了像传统有限元模型那样需要划分高密度二维单元。文中给出了代表奇异性应力场强度的夹杂角部广义应力强度因子数值解,并考虑材料属性、夹杂尺寸和夹杂位置关系的影响。算例中,考虑了夹杂和基体完全接合的情况,并给出了考核例。结果表明:当前模型能得到高精度数值解,且收敛性好;与传统有限元法和积分方程方法相比,该模型更具有通用性,为非均质材料的细观力学分析打下了基础。     

热-机载荷下不规则夹杂的奇异性应力场分析

基于有限元特征分析法得到的夹杂角部场数值特征解开发了一种超级奇异单元模型,并将其与普通四节点单元紧密结合,用于热-机载荷下夹杂角端部的应力场分析。在数值计算中,考察了热-机载荷下不同弹性比和不同夹杂尺寸的应力强度因子,并将所得结果与文献解和传统有限元方法解比对。结果表明,本文方法对热-机耦合条件下的不规则夹杂角端部的热弹性应力分析极为有效,可避免局部网格的高度加密,并提高计算效率。模型在复合材料夹杂的局部强度问题分析方面具有很好的实用性。     

椭圆类夹杂分析的杂交应力有限元方法

本文给出了一种分析椭圆类夹杂周边应力场的新型杂交应力有限元方法。基于弹性力学中平面问题的Muskhelishvili复势方法,应用保角变换映射技术,以Laurent级数和Faber级数为工具,借助Hellinger-Reissner原理构建一个能够反映椭圆类夹杂周边弹性现象同时包含椭圆夹杂的多边形超级单元。将该超级单元与标准的4节点杂交应力单元耦合在一起即可建立一种分析椭圆类夹杂周边弹性场的新型特殊杂交应力有限元方法。文中考核算例表明：本文方法不但使用简单、有效,而且精度高、单元少。作为本文方法的一个拓展应用,文章最后给出了一个分析含二个椭圆夹杂无限大各向同性板在远场均布载荷作用下椭圆夹杂周边弹性场的算例,并讨论了椭圆夹杂间距和弹性刚度比对应力集中系数的影响。     

Stress intensity factor analysis of a three-dimensional interfacial corner between anisotropic bimaterials under thermal stress

A numerical method using a path-independent H-integral based on the conservation integral was developed to analyze the singular stress field of a three-dimensional interfacial corner between anisotropic bimaterials under thermal stress. In the present method, the shape of the corner front is smooth. According to the theory of linear elasticity, asymptotic stress near the tip of a sharp interfacial corner is generally singular as a result of a mismatch of the materials’ elastic constants. The eigenvalues and the eigenfunctions are obtained using the Williams eigenfunction method, which depends on the anisotropic materials’ properties and the geometry of an interfacial corner. The order of the singularity related to the eigenvalue is real, complex or power-logarithmic. The amplitudes of the singular stress terms can be calculated using the H-integral. The stress and displacement around an interfacial corner for the H-integral are obtained using finite element analysis. In this study, a proposed definition of the stress intensity factors of an interfacial corner, which includes those of an interfacial crack and a homogeneous crack, is used to evaluate the singular stress fields. Asymptotic solutions of stress and displacement around an interfacial corner front are uniquely obtained using these stress intensity factors. To prove the accuracy of the present method, several different kinds of examples are shown such as interfacial corners or cracks in three-dimensional structures.     

受热载荷两相材料界面端应力场的新型有限元分析

基于一种特殊有限元特征分析方法获得两相材料界面端奇异性应力和位移场数值特征解, 据此开发了一种新型超级单元模型, 用于分析热载荷作用下两相材料界面端的应力场. 与机械载荷作用下超级单元模型的区别在于, 该模型在能量泛函中考虑了热-机耦合的影响, 将应力场分为奇异项和非奇异项, 而奇异性项又可分解为热致部分和力致部分. 模型的有效性通过了经验解和传统有限元方法的验证;模型可以避免在界面端邻域网格高度加密, 提高了计算速度, 对于分析多奇异性点应力干涉问题有重要意义.     

Order of singularity and singular stress field about an axisymmetric interface corner in three-dimensional isotropic elasticity

The spatial axisymmetric problem of an isotropic, elastic, homogeneous body containing an inclusion of a different material with an interface corner point (along a circular contour) and arbitrary joining angles is considered in this paper. It is found that the order of the stress singularity at the interface corner coincides with that of the corresponding plane strain problem, but the dependence of the singular stress field on the material properties depends on both the Poisson ratios (of the two isotropic materials) as well as on the ratio of their shear moduli. The theoretical results have been confirmed by numerical, finite-element-based results in a special bimaterial case.     

A unified definition for stress intensity factors of interface corners and cracks

Based upon linear fracture mechanics, it is well known that the singular order of stresses near the crack tip in homogeneous materials is a constant value −1/2, which is nothing to do with the material properties. For the interface cracks between two dissimilar materials, the near tip stresses are oscillatory due to the order of singularity being −1/2 ± iε and −1/2. The oscillation index ε is a constant related to the elastic properties of both materials. While for the general interface corners, their singular orders depend on the corner angle as well as the elastic properties of the materials. Owing to the difference of the singular orders of homogeneous cracks, interface cracks and interface corners, their associated stress intensity factors are usually defined separately and even not compatibly. Since homogenous cracks and interface cracks are just special cases of interface corners, in order to build a direct connection among them a unified definition for their stress intensity factors is proposed in this paper. Based upon the analytical solutions obtained previously for the multibonded anisotropic wedges, the near tip solutions for the general interface corners have been divided into five different categories depending on whether the singular order is distinct or repeated, real or complex. To provide a stable and efficient computing approach for the general mixed-mode stress intensity factors, the path-independent H-integral based on reciprocal theorem of Betti and Rayleigh is established in this paper. The complementary solutions needed for calculation of H-integral are also provided in this paper. To illustrate our results, several different kinds of examples are shown such as cracks in homogenous isotropic or anisotropic materials, central or edge notches in isotropic materials, interface cracks and interface corners between two dissimilar materials.     

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.     

Analysis of stress concentrations in plates with rectangular openings by a combined conformal mapping – Finite element approach

Plates with rectangular openings develop stress concentrations under bending. While these stresses can be determined using finite elements, in many problems this would be difficult because a high density mesh would be needed in the neighborhood of every opening corner. In this paper, it is shown how a complex-variable conformal mapping approach can be numerically coupled with the finite element method to analyze these corner stresses. This approach can be used even with relatively coarse meshes where the finite element results, by themselves, do not resolve the stress concentrations. In essence, the method relies on finite element analysis to obtain information on the stress field in a region surrounding the plate opening; this information is subsequently used to set up the parameters of the conformal mapping approach to obtain the near field stresses at the opening corners.     

THE INTERACTION PROBLEM BETWEEN THE ELASTIC LINE INCLUSIONS

ＩｎｔｒｏｄｕｃｔｉｏｎＯｆａｌｌｔｈｅｆｉｂｅｒ＿ｒｅｉｎｆｏｒｃｅｄｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌｓ,ｔｈｅｓｈｏｒｔ＿ｆｉｂｅｒｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌｎｏｔｏｎｌｙｓｔｒｅｎｇｔｈｅｎｓｔｈｅｍａｔｒｉｘｂｕｔａｖｏｉｄｓｄｅｆｅｃｔｉｏｎｓｏｆｔｈｅｌｏｎｇ＿ｆｉｂｅｒｃｏｍｐｏｓｉｔｅｍａｔｅｒｉａｌａｓｗｅｌｌ.Ｔｈｅｍｉｃｒｏ＿ｍｅｃｈａｎｉｃｓａｂｏｕｔｉｔｓｕｃｈａｓｆｒａｃｔｕｒｅ ,ｆａｔｉｇｕｅａｎｄｄａｍａｇｅｉｓｖｅｒｙｃｏｍｐｌｅｘ .Ｉｎｔｈｅｐｒｅｖｉ…     

A semi-analytical elastic solution for stress field of lined non-circular tunnels at great depth using complex variable method

In this paper, a semi-analytical elastic plane strain solution was provided for stress field around a lined non-circular tunnel subjected to uniform ground load. Concrete lining and the surrounding rock mass were assumed as linearly elastic materials. Due to complexity of the problem for non-circular geometric configurations, complex variable method introduced by Muskhelishvili and conformal mapping functions were used to determine stress components within concrete lining and the surrounding rock mass. Finally, the solution was validated by ABAQUS finite element software through an example. Very good agreement was demonstrated between semi-analytical and numerical solution although some discrepancies were found at tunnel corners where large curvature existed. It was demonstrated that the solution predicted stress components more accurately around the tunnels, especially the corners with large stress concentration. Practical significance of the solution was placed in the fact that it could be used as a quick-solver with high accuracy.     

压电材料中切口/接头端部平面电弹性场奇异性有限元分析

为了对平面载荷作用下压电材料中切口或接头端部附近电弹性场奇异性问题进行分析，首先以应力平衡方程、Maxwell方程和和边界条件为基础，得到一种求解压电材料特征问题的弱式方程；其次，假定楔形切口或接头端部附近单元内位移和电势沿径向分布为指数形式，而周向方向分布则采用泡函数插值，将其代入弱式方程，建立一种只需对楔形切口或接头端部附近周边进行离散的一维简单有限元方法．压电材料的极化轴可以是任意方向．利用该有限元模型讨论了楔形切口角度、极化轴方向和边界条件对奇性场的影响．通过和其它特定情况下的现有解相比，证实了该文有限元数值方法的有效性，而且精度很高．     

一种新的计算各向异性材料裂纹尖端应力强度因子杂交元法

利用有限元特征分析法研究了平面各向异性材料裂纹端部的奇性应力指数以及应力场和位移场的角分布函数,以此构造了一个新的裂纹尖端单元。文中利用该单元建立了研究裂纹尖端奇性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型。与四节点单元相结合,由此提出了一种新的求解应力强度因子的杂交元法。最后给出了在平面应力和平面应变下求解裂纹尖端奇性场的算例。算例表明,本文所述方法不仅精度高,而且适应性强。     

刚性夹杂角点处奇异场的数值分析及双夹杂间的相互影响

本文采用编制的有限元程序求解了无限大板中含有单个具有角点的刚性夹杂在无限远处均布外载作用下夹杂角点附近的位移场，并结合Ｉｎｓｈｉｋａｗａ的工作，求得了相应的应力强度因子。在此基础上，着重分析了两个刚性夹杂间的相互作用，得到了双夹杂左右角点附近的ＫⅠ，Ⅱ值随夹寻间距离改变而变化的规律。     

Interaction of plane nonstationary waves with a thin elastic inclusion under smooth contact conditions

We solve the problem of determining the stress state near a thin elastic inclusion in the form of a strip of finite width in an unbounded elastic body (matrix) with plane nonstationary waves propagating through it and with the forces exerted by the ambient medium taken into account. We assume that the matrix is in the plane strain state, and the smooth contact conditions are realized on both sides of the inclusion. The method for solving this problem consists in using the integral Laplace transform with respect to time and in representing the stress and displacement images in terms of the discontinuous solution of Lamé equations in the case of plane strain. As a result, the initial problem is reduced to a system of singular integral equations for the transforms of the unknown stress and displacement jumps. To invert the Laplace transform, we use a numerical method based on replacing the Mellin integral by the Fourier series. As a result, we obtain approximate formulas for calculating the stress intensity factors (SIF) for the inclusion, which are used to study the SIF time-dependence and its influence on the values of the inclusion rigidity. We also studied the possibility of considering the inclusions of higher rigidity as absolutely rigid inclusions.     

基于新型裂尖杂交元的压电材料断裂力学研究

提出了一种裂尖邻域杂交元模型，将其与标准杂交应力元结合来求解压电材料裂纹尖 端的奇性电弹场和断裂参数的数值解．裂纹尖端杂交元的建立步骤为：1) 利用高次内插有限元特征法求解特征问题，得到反映裂尖奇异性电弹场状况的特 征值和特征角分布函数；2) 利用广义Hellinger-Reissner变分泛函以及特征问题的解来建立裂尖邻域杂交元模型．该 方法求解电弹场时，摒弃了传统有限元方法中裂尖奇异性场需要借助解析解的做法，也避免 了单纯有限元方法中需要在裂尖端部进行高密度单元划分．采用PZT5板中心裂纹问题 作为考核例，数值结果显示了良好的精确性．作为进一步应用，求解了含中心界面裂纹 的PZT4-PZT5两相压电材料的应力强度因子和电位移强度因子．所有的算例都考虑 了3种裂纹面电边界条件．     

Variational bounds for anisotropic elastic multiphase composites with different shapes of inclusions

In the present work, unified formulae for the overall elastic bounds for multiphase transversely isotropic composites with different geometrical types of inclusions embedded in a matrix are calculated, including the spherical and long or short continuous cylindrical fiber cases. The influence of the different geometrical configurations of the inclusions on the composites is studied. The transversely isotropic effective bounds are obtained by applying the variational formulation for anisotropic composites developed by Willis, which relies on expressions for the static transversely isotropic Green’s function. Some numerical calculations and comparisons with the effective coefficients derived from the self-consistent approach, asymptotic homogenization method, and finite element method (FEM) are shown for different aspect ratio values, exhibiting good agreement.