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

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

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

本文提出了一种分析热载荷作用下夹杂角端部奇异应力场的简便方法。该方法是通过三个物理过程来完成的,研究结果发现夹杂角端部奇异热应力场问题可以等效地表示成一个双向等拉伸载荷作用下的奇异应力场弹性问题。换句话说,我们可以直接利用双向拉伸载荷作用下夹杂角端部奇异应力场现有数值结果计算奇异热应力场数值解。与有限元热弹性分析结果表明,本文提出的等效方法是正确的、有效的     

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

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

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

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

热机载荷下多边形夹杂角端部应力场的杂交元分析

开发了一种多边形夹杂角端部超级单元模型,并将其与传统四节点单元组装,用于分析热-机载荷下结构中多边形夹杂角端部的应力场.与机械载荷作用下超级单元模型的区别在于,该模型将夹杂角部邻域应力场分为奇异项和非奇异项,而奇异性项又可分解为热致部分和力致部分.在数值计算中,首先分析了单正方形夹杂问题,验证了模型的有效性,考察了材料匹配的影响;然后分析了双正方形夹杂的干涉问题,考察了夹杂间距的影响.结果表明,当前模型可避免局部网格的高度加密,从而提高有限元计算的效率.     

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

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

黏弹性接触界面端附近的奇异应力场

研究蠕变加载条件下线黏弹性材料接触界面端附近的奇异应力场问题.考虑接触界面的摩擦,假设界面端的滑移方向不改变,相对滑移量微小,且其与位移同量级,由此线性化局部边界条件,根据对应原理得到Laplace变换域中的界面端应力场,导出时域中奇异应力场的卷积积分表达式.对卷积积分核函数进行数值反演,考虑接触材料的两类组合,一是持久模量具有量级上的差异,另一是持久模量接近相同.算例结果证实核函数可以用准弹性法求得的解析式较准确地近似.在此基础上,利用积分中值定理,并引入各应力分量的修正系数,得到黏弹性奇异应力场的简化式.结合核函数的数值反演结果分析修正系数表达式的取值范围,得到如下结论,若两相接触材料的持久模量相差很大,可以采用准弹性解的解析式较准确地描述界面端的奇异应力场;一般情况下,应力场不存在统一的奇异值和应力强度系数,当采用类似于准弹性解的表达式近似给出黏弹性应力场时,可以估计此近似描述的误差限.文中最后采用有限元分析黏弹性板端部嵌入部位的应力场,算例包括了黏弹性板与弹性金属支承、黏弹性板与黏弹性垫层所形成的滑移接触界面端,利用黏弹性有限元的数值结果验证理论分析所得结论的有效性.      

黏弹性接触界面端附近的奇异应力场

研究蠕变加载条件下线黏弹性材料接触界面端附近的奇异应力场问题.考虑接触界面的摩擦,假设界面端的滑移方向不改变,相对滑移量微小,且其与位移同量级,由此线性化局部边界条件,根据对应原理得到Laplace变换域中的界面端应力场,导出时域中奇异应力场的卷积积分表达式.对卷积积分核函数进行数值反演,考虑接触材料的两类组合,一是持久模量具有量级上的差异,另一是持久模量接近相同.算例结果证实核函数可以用准弹性法求得的解析式较准确地近似.在此基础上,利用积分中值定理,并引入各应力分量的修正系数,得到黏弹性奇异应力场的简化式.结合核函数的数值反演结果分析修正系数表达式的取值范围,得到如下结论,若两相接触材料的持久模量相差很大,可以采用准弹性解的解析式较准确地描述界面端的奇异应力场;一般情况下,应力场不存在统一的奇异值和应力强度系数,当采用类似于准弹性解的表达式近似给出黏弹性应力场时,可以估计此近似描述的误差限.文中最后采用有限元分析黏弹性板端部嵌入部位的应力场,算例包括了黏弹性板与弹性金属支承、黏弹性板与黏弹性垫层所形成的滑移接触界面端,利用黏弹性有限元的数值结果验证理论分析所得结论的有效性.     

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

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

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

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

A novel hybrid finite element analysis of inplane singular elastic field around inclusion corners in elastic media

This paper deals with the inplane singular elastic field problems of inclusion corners in elastic media by an ad hoc hybrid-stress finite element method. A one-dimensional finite element method-based eigenanalysis is first applied to determine the order of singularity and the angular dependence of the stress and displacement field, which reflects elastic behavior around an inclusion corner. These numerical eigensolutions are subsequently used to develop a super element that simulates the elastic behavior around the inclusion corner. The super element is finally incorporated with standard four-node hybrid-stress elements to constitute an ad hoc hybrid-stress finite element method for the analysis of local singular stress fields arising from inclusion corners. The singular stress field is expressed by generalized stress intensity factors defined at the inclusion corner. The ad hoc finite element method is used to investigate the problem of a single rectangular or diamond inclusion in isotropic materials under longitudinal tension. Comparison with available numerical results shows the present method is an efficient mesh reducer and yields accurate stress distribution in the near-field region. As applications, the present ad hoc finite element method is extended to discuss the inplane singular elastic field problems of a single rectangular or diamond inclusion in anisotropic materials and of two interacting rectangular inclusions in isotropic materials. In the numerical analysis, the generalized stress intensity factors at the inclusion corner are systematically calculated for various material type, stiffness ratio, shape and spacing position of one or two inclusions in a plate subjected to tension and shear loadings.     

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

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

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.     

Torsion problems for a cylinder with a rectangular hole and a rectangular cylinder with a crack

Using the method of singular integral equation and the crack-cutting technique, the rigorous solutions are obtained for a cylinder with a rectangular hole and a rectangular cylinder with a crack, which exactly satisfy the boundary conditions and the conditions at the corner points. After that the torsional rigidities and the stress intensity factors at the crack tip are determined. Next, for the doubly connected circular cylinder with a rectangular hole the expressions for the singular stresses around the concave corner points are derived and the generalized stress intensity factors are then defined. Since the crack-cutting technique is used in this paper, the solution of the matching rectangular cylinder is also obtained and its numerical results coincide with those in references. Thus the method proposed here is verified. The project supported by National Natural Science Foundation of China     

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.     

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

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