求解双材料裂纹结构全域应力场的扩展边界元法

在线弹性理论中,复合材料裂纹尖端具有多重应力奇异性,常规数值方法不易求解.该文建立的扩展边界元法(XBEM)对围绕尖端区域位移函数采用自尖端径向距离r的渐近级数展开式表达,其幅值系数作为基本未知量,而尖端外部区域采用常规边界元法离散方程.两方程联立求解可获得裂纹结构完整的位移和应力场.对两相材料裂纹结构尖端的两个材料域分别采用合理的应力特征对,然后对其进行计算,通过计算结果的对比分析,表明了扩展边界元法求解两相材料裂纹结构全域应力场的准确性和有效性.     

Evaluation of the stress singularities of plane V-notches in bonded dissimilar materials

According to the linear theory of elasticity, there exists a combination of different orders of stress singularity at a V-notch tip of bonded dissimilar materials. The singularity reflects a strong stress concentration near the sharp V-notches. In this paper, a new way is proposed in order to determine the orders of singularity for two-dimensional V-notch problems. Firstly, on the basis of an asymptotic stress field in terms of radial coordinates at the V-notch tip, the governing equations of the elastic theory are transformed into an eigenvalue problem of ordinary differential equations (ODEs) with respect to the circumferential coordinate θ around the notch tip. Then the interpolating matrix method established by the first author is further developed to solve the general eigenvalue problem. Hence, the singularity orders of the V-notch problem are determined through solving the corresponding ODEs by means of the interpolating matrix method. Meanwhile, the associated eigenvectors of the displacement and stress fields near the V-notches are also obtained. These functions are essential in calculating the amplitude of the stress field described as generalized stress intensity factors of the V-notches. The present method is also available to deal with the plane V-notch problems in bonded orthotropic multi-material. Finally, numerical examples are presented to illustrate the accuracy and the effectiveness of the method.     

Elastic stress analysis of blunt V-notches under mixed mode loading by considering higher order terms

This paper presents the in-plane asymptotic displacement and stress fields for blunt V-notched components based on Kolosov–Muskhelishvili's approach. In the first part, the displacement and stress components in the polar coordinate system are determined by choosing appropriate complex potential functions. In order to construct the notch geometry, the Neuber's mapping relation is utilized. Then, the notch boundary conditions are imposed to calculate the free parameters of the stress distribution. Eventually, the stress and displacement components are calculated in the Cartesian and polar coordinates in the forms of series expansion. In the second part, the coefficients of series expansions are computed by using the least square method (LSM). The blunt V-notched Brazilian disk (BV-BD) specimen under mixed mode loading is used as an example to verify the proposed procedure. The stress components in arbitrary distances and directions are determined for different blunt V-notches in order to evaluate the accuracy of the calculated stress series solutions and their associated coefficients. The numerical results indicate that a single-term solution can lead to considerable errors, and to achieve good accuracy in the stress field calculation, one should take account of at least three terms in the stress series solution.     

三维切口尖端应力应变场

本文利用双重幂级数展开法分析三维切口尖端应力应变奇异性，通过切口边界条件导出切口特征方程，进而求得不同切口内角下特征值序列解答，最后推得切口尖端应力应变场。     

Complete elastic–plastic stress asymptotic solutions near general plane V-notch tips

An efficient method is developed to determine the multiple term eigen-solutions of the elastic–plastic stress fields at the plane V-notch tip in power-law hardening materials. By introducing the asymptotic expansions of stress and displacement fields around the V-notch tip into the fundamental equations of elastic–plastic theory, the governing ordinary differential equations (ODEs) with the stress and displacement eigen-functions are established. Then the interpolating matrix method is employed to solve the resulting nonlinear and linear ODEs. Consequently, the first four and even more terms of the stress exponents and the associated eigen-solutions are obtained. The present method has the advantages of greater versatility and high accuracy, which is capable of dealing with the V-notches with arbitrary opening angle under plane strain and plane stress. In the present analysis, both the elastic and the plastic deformations are considered, thus the complete elastic and plastic stress asymptotic solutions are evaluated. Numerical examples are shown to demonstrate the accuracy and effectiveness of the present method.     

反平面塑性V形切口尖端应力和位移渐近解

提出了一种确定幂硬化材料反平面V形切口尖端应力和位移渐近解的主导项和高阶项的有效方法。首先通过在弹塑性理论基本方程中引入V形切口尖端应力场和位移场的渐近级数展开,建立以应力和位移为特征函数的非线性和线性常微分方程组。然后采用插值矩阵法求解常微分方程组,可得到多阶应力特征指数和其相对应的特征函数。该方法具有通用性强、精度高等优点,可处理任意开口角度和应变硬化指数的V形切口。典型算例验证了该方法的准确性和有效性。     

Analysis of the singular function boundary integral method for a biharmonic problem with one boundary singularity

In this article, we analyze the singular function boundary integral method (SFBIM) for a two‐dimensional biharmonic problem with one boundary singularity, as a model for the Newtonian stick‐slip flow problem. In the SFBIM, the leading terms of the local asymptotic solution expansion near the singular point are used to approximate the solution, and the Dirichlet boundary conditions are weakly enforced by means of Lagrange multiplier functions. By means of Green's theorem, the resulting discretized equations are posed and solved on the boundary of the domain, away from the point where the singularity arises. We analyze the convergence of the method and prove that the coefficients in the local asymptotic expansion, also referred to as stress intensity factors, are approximated at an exponential rate as the number of the employed expansion terms is increased. Our theoretical results are illustrated through a numerical experiment. © 2011 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2011     

Asymptotic expansions using blow-up

The method of matched asymptotic expansions and geometric singular perturbation theory are the most common and successful approaches to singular perturbation problems. In this work we establish a connection between the two approaches in the context of the simple fold problem. Using the blow-up technique [5], [12] and the tools of geometric singular perturbation theory we derive asymptotic expansions of slow manifolds continued beyond the fold point. Our analysis explains the structure of the expansion and gives an algorithm for computing its coefficients.     

Superelements for the finite element solution of two-dimensional elliptic problems with boundary singularities

A novel singular superelement (SSE) formulation has been developed to overcome the loss of accuracy encountered when applying the standard finite element schemes to two-dimensional elliptic problems possessing a singularity on the boundary arising from an abrupt change of boundary conditions or a reentrant corner. The SSE consists of an inner region over which the known analytic form of the solution in the vicinity of the singular point is utilized, and a transition region in which blending functions are used to provide a smooth transition to the usual linear or quadratic isoparametric elements used over the remainder of the domain. Solution of the finite element equations yield directly the coefficients of the asymptotic series, known as the flux/stress intensity factors in linear heat transfer or elasticity theories, respectively. Numerical examples using the SSE for the Laplace equation and for computing the stress intensity factors in the linear theory of elasticity are given, demonstrating that accurate results can be attained for a moderate computational effort.     

Asymptotic expansions using blow-up

The method of matched asymptotic expansions and geometric singular perturbation theory are the most common and successful approaches to singular perturbation problems. In this work we establish a connection between the two approaches in the context of the simple fold problem. Using the blow-up technique [5], [12] and the tools of geometric singular perturbation theory we derive asymptotic expansions of slow manifolds continued beyond the fold point. Our analysis explains the structure of the expansion and gives an algorithm for computing its coefficients.*Research supported by the Austrian Science Foundation under grant Y 42-MAT.Received: February 1, 2001; revised: November 22, 2002     

关于边界层方法

本文指出传统的边界层方法(包括匹配法和Vi?ik—Lyusternik方法)的不足:不能作出边界层项的渐近展开式.提出多重尺度构造边界层项的方法,得到符合实情的结果.又与Levinson所用的方法比较,本方法能更简单地导出后一方法给出的边界层项的渐近展开式.又应用此方法研究现有的关于奇异摄动的某些成果,指出这些成果的局限性,并在一般情况下作出解的渐近展开式.     

On the extraction technique in boundary integral equations

In this paper we develop and analyze a bootstrapping algorithm for the extraction of potentials and arbitrary derivatives of the Cauchy data of regular three-dimensional second order elliptic boundary value problems in connection with corresponding boundary integral equations. The method rests on the derivatives of the generalized Green's representation formula, which are expressed in terms of singular boundary integrals as Hadamard's finite parts. Their regularization, together with asymptotic pseudohomogeneous kernel expansions, yields a constructive method for obtaining generalized jump relations. These expansions are obtained via composition of Taylor expansions of the local surface representation, the density functions, differential operators and the fundamental solution of the original problem, together with the use of local polar coordinates in the parameter domain. For boundary integral equations obtained by the direct method, this method allows the recursive numerical extraction of potentials and their derivatives near and up to the boundary surface.

     

The asymptotic solution of the contact problem for a three-dimensional elastic body of finite dimensions

The asymptotic form of Green's vector function with a pole on the boundary is calculated by the method of matched asymptotic expansions. The expansion obtained is used to construct the asymptotic form of the contact pressure. The equations of the contact problem are derived with integral corrections, which take into account the nature of the attachment and the geometry of the elastic body. Examples of calculations for an elliptic punch are given.     

一类具非线性边界条件的高阶方程的奇摄动问题

通过引入伸展变量和非常规的渐近序列{∈}）,运用合成展开法,对一类具非线性边界条件的非线性高阶微分方程的奇摄动问题构造了形式渐近解,再运用微分不等式理论证明了原问题解的存在性及所得渐近近似式的一致有效性．     

Corner asymptotics of the magnetic potential in the eddy‐current model

In this paper, we describe the magnetic potential in the vicinity of a corner of a conducting body embedded in a dielectric medium in a bidimensional setting. We make explicit the corner asymptotic expansion for this potential as the distance to the corner goes to zero. This expansion involves singular functions and singular coefficients. We introduce a method for the calculation of the singular functions near the corner, and we provide two methods to compute the singular coefficients: the method of moments and the method of quasi‐dual singular functions. Estimates for the convergence of both approximate methods are proven. We eventually illustrate the theoretical results with finite element computations. The specific nonstandard feature of this problem lies in the structure of its singular functions: They have the form of series whose first terms are harmonic polynomials, and further terms are genuine nonsmooth functions generated by the piecewise constant zeroth order term of the operator. Copyright © 2013 John Wiley & Sons, Ltd.     

非Fourier温度场分布的奇摄动解

应用非Fourier热传导定律构建了单层材料中温度场模型，即一类在无界域上带小参数的奇摄动双曲方程，通过奇摄动展开方法，得到了该问题的渐近解.首先应用奇摄动方法得到了该问题的外解和边界层矫正项，通过对内解和外解的最大模估计和关于时间导数的最大模估计以及线性抛物方程理论，得到了内外解的存在唯一性，从而得到了解的形式渐近展开式.通过余项估计，给出了渐近解的L2估计，得到了渐近解的一致有效性，从而得到了无界域上温度场的分布.通过奇摄动分析，给出了非Fourier 温度场与Fourier 温度场的关系，描述了非Fourier温度场的具体形态.     

一类二阶非线性奇摄动边值问题解的渐近展开

本文研究一类非线性微分方程的非线性边值问题的奇摄动,应用边界层校正法构造出解的形式渐近展开式,并借助于上,下解及微分不等式理论研究解及其一阶导数的有关余项估计。     

An asymptotic analysis of the stress—strain state of a strip reinforced with ribs

An approximate solution of the problem of the stress—strain state of an anisotropic strip reinforced with two-dimensional ribs is constructed using the method of asymptotic expansion of generalized functions, the averaging method and the method of singular expansions.     

蠕变损伤材料中缺口尖端稳定扩展的应力场

研究了在拉伸载荷和反平面载荷作用下蠕变损伤材料缺口尖端稳定扩展的应力场．假设材料的应力及位移势函数,得到了缺口尖端场的各参数分量,进而在小范围蠕变条件下,建立了缺口尖端稳定扩展的蠕变损伤控制方程,并考虑缺口尖端蠕变钝化作用和问题的边界条件,对控制方程进行了数值分析,得到了缺口尖端的应力场,并讨论了缺口尖端应力场随各影响参数的变化规律．结果表明,缺口尖端的应力具有r1/(1-n)的奇异性,应力率具有rn/(1-n)的奇异性,n是蠕变指数．     

微分求积法分析平面接头应力奇异性北大核心CSCD

对于双材料平面接头问题提出了一个分析应力奇性指数的新方法:微分求积法(DQM).首先,将平面接头连接点处位移场的径向渐近展开格式代入平面弹性力学控制方程,获得了关于应力奇性指数的常微分方程组(ODEs)特征值问题.然后,基于DQM理论,将ODEs的特征值问题转化为标准型广义代数方程组特征值问题,求解之可一次性地计算出双材料平面接头连接点处应力奇性指数,同时,一并求出了接头连接点处相应的位移和应力特征函数.数值计算结果说明该文DQM计算平面接头连接点处应力奇性指数的结果是正确的.