动态断裂力学的无限相似边界元法

对弹性动力学的相似边界元法进行了进一步研究,推导了相应的计算公式,并在此基础上提出了动态断裂力学的无限相似边界元法.与传统的边界元法相比,相似边界元法由于只需在少数单元上进行数值积分,大大减少了计算量.对动态断裂力学问题,无限相似边界元法由于在裂纹尖端的边界上设置了逼近于裂纹尖端的无限个相似边界单元,可直接得到裂纹尖端具有奇异性的应力,而不需要设置奇异单元,从而突破了奇异单元对应力奇异性阶次的局限.另外,还讨论了无限相似边界元法得到的无限阶的线性代数方程组的求解方法.     

Boundary element analysis of contact problems using artificial boundary node approach

An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary element and artificial boundary node approaches. The exisiting and derived approaches are tested using some case studies. The results of the artificial boundary node approach are compared with those of the existing boundary element program, the regular element approach, ANSYS and analytical solution whenever possible. The results show the effectiveness of the artificial boundary node approach for a wider range of boundary offsets.     

High accuracy non-equidistant method for singular perturbation reaction-diffusion problem

Singular perturbation reaction-diffusion problem with Dirichlet boundary condition is considered. It is a multi-scale problem. Presence of small parameter leads to boundary layer phenomena in both sides of the region. A non-equidistant finite difference method is presented according to the property of boundary layer. The region is divided into an inner boundary layer region and an outer boundary layer region according to transition point of Shishkin. The steps sizes are equidistant in the outer boundary layer region. The step sizes are gradually increased in the inner boundary layer region such that half of the step sizes are different from each other. Truncation error is estimated. The proposed method is stable and uniformly convergent with the order higher than 2. Numerical results are given, which are in agreement with the theoretical result.     

近场波动模拟的一种应力人工边界

采用平面波和远场散射波混合透射，引入无限介质线弹性本构关系建立了一种应力人工 边界条件. 其优点在于边界结点反应与内部有限元结点反应采用相同的积分格式计算，有限 元积分方法稳定时不存在人工边界失稳问题. 数值算例表明：边界精度高于现有黏性边 界、黏弹性人工边界，以及一、二阶透射人工边界.     

Effect of boundary condition nonlinearities on free large-amplitude vibrations of rectangular plates

The effect of boundary condition nonlinearities on free nonlinear vibrations of thin rectangular plates is analyzed. The method for analysis of the plate vibrations with geometrical nonlinearity and the boundary condition nonlinearity is suggested. The nonlinear boundary conditions for membrane forces are transformed into linear ones using the in-plane stress function. Additional boundary conditions for the in-plane displacements vanishing on the clamped edge of the plate are imposed on the stress function. Simply supported and cantilever plates are analyzed. The backbone curves obtained by satisfying linear and nonlinear boundary conditions are compared. It is shown that the results of the calculations with nonlinear boundary conditions differ essentially from the data obtained without these boundary conditions.     

浸入边界法及其应用

浸入边界法主要用于模拟存在复杂外形结构的流场的运动情况和处理各种动边界问题,目前已广泛应用于计算流体力学领域.浸入边界法既是数学建模方法又是数值离散方法,它将物体边界与流体的相互作用通过在流体运动方程中加体积力项来体现,并在数值计算中采用简单的笛卡尔网格,避免了按照物体边界形状生成贴体网格时所遇到的各种问题.浸入边界法分为连续力法和离散力法:连续力法主要用于处理弹性边界问题,它的力源项满足特定的力学关系式;离散力法主要用于处理固体界面问题,它的力源项由边界条件推导得到.着重阐述了浸入边界法的基本原理和数学构造,对目前已有的几种不同的浸入边界法做了简单地介绍,并给出了一些应用实例,最后提出了浸入边界法未来的发展方向.     

各向异性弹性力学一般边值问题的广义Stroh公式

当边值问题是简单的,即是应力边值问题时,Stroh公式是很有效的。对于混合边值问题,倒如滑动边界条件,Stroh公式中的简洁的矩阵表达式就失效了。我们提出了一个广义的Stroh公式,它可应用于一大类一般的边界条件。简单的边界条件和滑动边界条件是这一类一般边界条件的特殊情形。值得指出的是,这个关于Stroh公式所作的修正并不大。广义的公式和最后的解答看起来很类似于未修正的原公式和原来的解。然而这个修正却可应用于相当广的边界条件。     

Generalized stroh formalism for anisotropic elasticity for general boundary conditions

The Stroh formalism is most elegant when the boundary conditions are simple, namely, they are prescribed in terms of traction or displacement. For mixed boundary conditions such as there for a slippery boundary, the concise matrix expressions of the Stroh formalism are destroyed. We present a generalized Stroh formalism which is applicable to a class of general boundary conditions. The general boundary conditions include the simple and slippery boundary conditions as special cases. For Green's functions for the half space, the general solution is applicable to the case when the surface of the half-space is a fixed, a free, a slippery, or other more general boundary. For the Griffith crack in the infinite space, the crack can be a slit-like crack with free surfaces, a rigid line inclusion (which is sometimes called an anticrack), or a rigid line with slippery surface or with other general surface conditions. It is worth mention that the modifications required on the Stroh formalism are minor. The generalized formalism and the final solutions look very similar to those of unmodified version. Yet the results are applicable to a rather wide range of boundary conditions.     

On the existence of solutions of boundary value problems of duffing type systems

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎ [1 ] ,ｕｎｄｅｒａｇｒｏｕｐｏｆｖｅｒｙｅｘｔｅｎｓｉｖｅｃｏｎｄｉｔｉｏｎｓ,ＳＨＥＮＺｕ＿ｈｅｓｔｕｄｉｅｄｔｈｅｅｘｉｓｔｅｎｃｅｏｆａｕｎｉｑｕｅ 2π＿ｐｅｒｉｏｄｉｃｓｏｌｕｔｉｏｎｏｆｔｈｅｓｙｓｔｅｍｏｆｏｒｄｉｎａｒｙｄｉｆｆｅｒｅｎｔｉａｌｅｑｕａｔｉｏｎｓｕ″(ｔ) Ｇ(ｕ(ｔ) ) =ｐ(ｔ) ,(1 )ｗｈｅｒｅＧ :Ｒｎ →Ｒｈａｓａｃｏｎｔｉｎｕｏｕｓｓｅｃｏｎｄｐａｒｔｉａｌｄｅｒｉｖａｔｉｖｅｓ,ａｎｄｐ:Ｒ→Ｒｎｉｓｃｏｎｔｉｎｕｏｕｓａｎｄ2π＿ｐｅｒ…     

Boundary eigensolutions in elasticity II. Application to computational mechanics

The theory of fundamental boundary eigensolutions for elastostatic problems, developed in Part I, is applied to formulate methods for computational mechanics. This theory shows that every elastic solution can be written as a linear combination of some fundamental boundary orthogonal deformations, thus providing a generalized Fourier expansion. One finds that traditional boundary element and finite element methods are largely consistent with this theory, but do not harness its full power. This theory shows that these computational methods are indirectly a generalized discrete Fourier analysis. Furthermore, by utilizing suitable boundary weight functions, boundary element and finite element formulations may be written exclusively in terms of bounded quantities, even for non-smooth problems involving notches, cracks, mixed boundary conditions and bi-material interfaces. The close relationship between the resulting boundary element and finite element methods also becomes evident. Both use displacement and surface traction as primary variables. A new degree-of-freedom concept is introduced, along with a stiffness tensor that enables one to visualize a finite element method via a boundary discretization process, just as in a boundary element approach. Global convergence characteristics of the traction-oriented finite element method are also developed. Comparisons with closed-form fundamental boundary eigensolutions for a circular elastic disc are presented in order to provide a means for assessing the numerical methods. Several other numerical examples are solved efficiently by using the concept of boundary eigensolutions in an indirect fashion. The results indicate that the algorithms follow the underlying theory and that solutions to non-smooth problems can be obtained in a systematic manner. Beyond this, the concept of boundary eigensolutions provides an alternative view of computational continuum mechanics that may lead to the development of other non-traditional approaches.     

Approximate nonreflecting inflow-outflow boundary conditions for the calculation of navier-stokes equations in internal flows

In this paper, the nonreflecting boundary conditions based upon fundamental ideas of the linear analysis are developed for gas dynamic equations, and the modified boundary conditions for Navier-Stokes equations are proposed as a substitute of the nonreflecting boundary conditions inside boundary layers near rigid walls. These derived boundary conditions are then applied to calculations both for the Euler equations and the Navier-Stokes equations to determine if they can produce acceptable results for the subsonic flows in channels. The numerical results obtained by an implicit second-order upwind difference scheme show the effectiveness and generality of the boundary conditions. Furthermore, the formulae and the analysis performed here may be extended to three dimensional problems. recommended by Prof. Cui Erjie     

非齐次边界条件下轴向运动梁的非线性振动

轴向运动系统的横向非线性振动一直是国内外研究的热点课题之一.目前相关研究大都是针对齐次边界条件的.但是在工程实际中,非齐次边界条件更为常见,而针对非齐次边界条件的研究相对较少.为深入研究非齐次边界条件对轴向运动系统横向非线性振动的影响,本文以轴向变速运动黏弹性Euler梁为例,引入由黏弹性引起的非齐次边界条件,同时还引入由轴向加速度引起的径向变化张力,建立梁横向振动的积分-偏微分型运动方程,并导出了相应的非齐次边界条件.采用直接多尺度法分析了梁的次谐波参数共振.由可解性条件得到了梁的稳态响应,并根据Routh-Hurvitz判据确定了系统稳态响应的稳定性.通过数值例子讨论了黏弹性系数,轴向运动速度,轴向速度脉动幅值和非线性系数对幅频响应的影响,并详细对比分析了非齐次边界条件和齐次边界条件对幅频响应的影响.结果表明:随着黏弹性系数的增大,非齐次边界条件下的零解失稳区域和稳态响应幅值比齐次边界条件下的失稳区域和幅值大,非齐次边界条件对高阶次谐波参数共振的影响更加显著.最后,引入微分求积法来验证直接多尺度法的近似解结果.      

Stratified Atmospheric Boundary Layers and Breakdown of Models

The goal of this study is to assess complications in atmospheric stable boundary layers which are not included in numerical models of the stably stratified boundary layer and to provide a formulation of surface fluxes for use in numerical models. Based on an extensive interpretive literature survey and new eddy correlation data for the stable boundary layer, this study defines two prototype stable boundary layers: the weakly stable case and the very stable case. The weakly stable boundary layer is amenable to existing models. The very stable boundary layer eludes modeling attempts due to breakdown of existing formulations of turbulence and due to features found in the atmosphere which are not normally included in models. The latter includes clear-air radiative cooling, low-level jets, surface heterogeneity, gravity waves, meandering motions, and other mesoscale motions which propagate from outside the local domain. While these mechanisms are not essential to understanding idealized or laboratory versions of the stable boundary layer, they complicate comparisons of numerical models and theories with actual atmospheric boundary layers. Statistics which describe various features of the stable boundary layer are offered for future comparison with modeling results. Received 13 June 1997 and accepted 19 November 1997     

弹性力学的一种边界无单元法

首先对移动最小二乘副近法进行了研究，针对其容易形成病态方程的缺点，提出了以带权的正交函数作为基函数的方法－改进的移动最小二乘副近法，改进的移动最小二乘逼近法比原方法计算量小，精度高，且不会形成病态方程组，然后，将弹性力学的边界积分方程方法与改进的移动最小二乘逼近法结合，提出了弹性力学的一种边界无单元法，这种边界无单元法法是边界积分方程的无网格方法，与原有的边界积分方程的无网格方法相比，该方法直接采用节点变量的真实解为基本未知量，是边界积分方程无网格方法的直接解法，更容易引入界条件，且具有更高的精度，最后给出了弹性力学的边界无单元法的数值算例，并与原有的边界积分方程的无网格方法进行了较为详细的比较和讨论。     

Exterior statistics based boundary conditions for representative volume elements of elastic composites

Statistically equivalent representative volume elements or SERVEs are representations of the microstructure that are used for micromechanical simulations to generate homogenized material constitutive responses and properties (Swaminathan et al., 2006a, Ghosh, 2011). Typically, a SERVE is generated from the parent microstructure as a statistically equivalent region, whose size is determined from the requirements of convergence of macroscopic properties. Standard boundary conditions, such as affine transformation-based displacement boundary conditions (ATDBCs), uniform traction boundary conditions (UTBCs) or periodic boundary conditions (PBCs) are conventionally applied on the SERVE boundary for micromechanical simulations. However, when the microstructure is characterized by arbitrary, nonuniform distributions of heterogeneities, these simple boundary conditions do not represent the effect of regions exterior to the SERVE. Improper boundary conditions can result in significantly larger than optimal SERVE domains, needed for converged properties. In an attempt to overcome the limitations of the conventional boundary conditions on the SERVE, this paper explores the effect of boundary conditions that incorporate the statistics of the exterior region on the SERVE of elastic composites. Using Green's function based interaction kernels, coupled with statistical functions of the microstructural characteristics like one-point and two-point correlation functions, a novel exterior statistics-based boundary condition or ESBC is derived for the SERVE. The advantages of the ESBC are established by comparing with results of simulations using conventional boundary conditions. Results of the SERVE simulations subjected to ESBCs are also compared with those from other popular methods like statistical volume element (SVE) and weighted statistical volume element (WSVE). The proposed ESBCs offer significant advantages over other methods in the SERVE-based analysis of heterogeneous materials.     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

Viscous boundary layers in flows through a domain with permeable boundary

The effect of small viscosity on nearly inviscid flows of an incompressible fluid through a given domain with permeable boundary is studied. The Vishik–Lyusternik method is applied to construct a boundary layer asymptotic at the outlet in the limit of vanishing viscosity. Mathematical problems with both consistent and inconsistent initial and boundary conditions at the outlet are considered. It is shown that in the former case, the viscosity leads to a boundary layer only at the outlet. In the latter case, in the leading term of the expansion there is a boundary layer at the outlet and there is no boundary layer at the inlet, but in higher order terms another boundary layer appears at the inlet. To verify the validity of the expansion, a number of simple examples are presented. The examples demonstrate that asymptotic solutions are in quite good agreement with exact or numerical solutions.     

Reflection of longitudinal displacement wave at the viscoelastically supported boundary of micropolar half-space

The reflection of elastic waves at the boundary of micropolar half-space with the viscoelastic support is studied in this paper. The spring-dashpot model is used to model the viscoelastic support. The boundary condition includes the force stress, couple stress, the displacements and the micro-rotation and the contribution from the spring and the dashpot on them. The amplitude ratios and phase shifts of reflection waves with respect to the incident wave are obtained from the visco-elastic boundary conditions. Further, the energy flux ratios of the reflection waves to the incident wave are estimated and the energy flux conservation with consideration of the energy dissipation of visco-elastic boundary is used to validate the numerical results. Based on the numerical results, the influences of spring and dashpot are studied respectively. It is found that the elastically supported boundary and the viscously supported boundary have evident different influences on the amplitude ratio and the phase shift. The causes resulting in these deviations are related with the instantaneous elasticity of elastic boundary and the time-delay effects of viscous boundary.     

A NOVEL BOUNDARY INTEGRAL EQUATION METHOD FOR LINEAR ELASTICITY--NATURAL BOUNDARY INTEGRAL EQUATION

The boundary integral equation (BIE) of displacement derivatives is put at a disadvantage for the difficulty involved in the evaluation of the hypersingular integrals. In this paper, the operators δij and εij are used to act on the derivative BIE. The boundary displacements, tractions and displacement derivatives are transformed into a set of new boundary tensors as boundary variables. A new BIE formulation termed natural boundary integral equation (NBIE) is obtained. The NBIE is applied to solving two-dimensional elasticity problems. In the NBIE only the strongly singular integrals are contained. The Cauchy principal value integrals occurring in the NBIE are evaluated. A combination of the NBIE and displacement BIE can be used to directly calculate the boundary stresses. The numerical results of several examples demonstrate the accuracy of the NBIE.     

Boundary integral equations of unique solutions in elasticity

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｂｏｕｎｄａｒｙｅｌｅｍｅｎｔｍｅｔｈｏｄ（ＢＥＭ）ｐｒｏｖｉｄｅｓａｎａｔｔｒａｃｔｉｖｅａｌｔｅｒｎａｔｉｖｅｆｏｒｔｈｅａｎａｌｙｓｉｓｏｆｅｎｇｉｎｅｅｒｉｎｇｐｒｏｂｌｅｍｓ．Ｉｔｓｍａｉｎａｄｖａｎｔａｇｅｓａｒｅｅｃｏｎｏｍｉｃａｌａｎｄｐａｒｔｉｃｕｌａｒｌｙｃｏｎｖｅｎｉｅｎｔｆｏｒｕｎｂｏｕｎｄｅｄｄｏｍａｉｎａｎｄｓｔｒｅｓｓｃｏｎｃｅｎｔｒａｔｉｏｎｐｒｏｂｌｅｍｓ．Ｔｈｅｂｏｕｎｄａｒｙｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎ（ＢＩＥ）ｉｓｔｈｅ…