边界伸缩原理

文[5]提出边界伸缩原理及边界伸缩法.本文补充叙述了边界伸缩原理,并根据已有研究成果给予了较严格地证明,进一步完善了边界伸缩法理论基础.     

无界区域上Stokes问题的自然边界元与有限元耦合法

§1.引言 对于用有限元方法求解平面有界区域上的Stokes问题,国内外已有大量工作,例如可见[2]、[9]及其所引文献.但对无界区域上的这一问题,由于区域的无界性给有限元方法带来了困难,边界元方法及边界元与有限元的耦合法便显示其优越性.本文提出用自然边界元与有限元的耦合法求解无界区域上的Stokes问题.这一耦合法早在作者以前的工作中被应用于求解调和问题、重调和问题和平面弹性问题,但将它用于求解     

一种改进边界元法在弹性扭转问题中的应用

本文用一种改进边界元法分析与计算了椭圆截面等直杆的扭转问题.并与边界元法的解进行比较,其结果极为符合.然而,改进边界元法较边界元法所需要的数据量少得多,计算时间也将大大减少了.因此,本文方法对求解Poisson方程问题是一种经济而行之有效的数值计算方法.     

三维瞬态弹性动力场的基本解及边界积分方程

三维瞬态动力场的求解,一直为世人所瞩目,本文拟采用边界积分方程——边界单元法进行求解。与其它数值方法相比较,本方法由于采用了动力场区域的基本解,因此在处理无限域或半无限域的动力场问题时可以毋须人为地去划定边界,附加相应的边界约束条件。同时在求解过程中,采用了权函数与奇异格林函数,使得域积分可以转化为边界积     

利用边界元法求解一类重调和方程

边界元法(BEM)和多重互易法(MRM)相结合求解一类重调和方程.通过重调和基本解序列给出的MRM-方法和BEM, 推导出该类问题的MRM-边界变分方程, 用边界元法求解该变分方程, 从而得到重调和方程的近似解, 并给出了解的存在唯一性证明.通过数值算例说明了MRM-方法具有收敛速度快、计算精度高, 易编程等优点, 为使用边界元法数值求解重调和方程提供了方法和理论依据.适合于工程中的实际运算.     

发展方程初边值问题的高阶差分-边界积分方程法及误差分析

本文结合差分方法与边界积分方程方法,提出并研究了一类新的求解发展型方程初边值问题的高阶差分与边界积分方程耦合数值方法.对于有界区域问题与无界区域问题给出了数值计算格式及其误差的先验估计.     

关于拟共形映照边界伸缩商的一个注记

本文给出了拟共形映照边界伸缩商与无限小边界伸缩商的一个等式h([μ])=inf_(μ1∈[μ])b([μ1]B)；并给出了一个关于T_0空间的推论．     

二维抛物型方程初边值问题的高阶差分边界元法及误差估计

本文讨论二维抛物型方程初边值问题的求解问题，提出了一种差分与边界元耦合的高阶求解方法，并给出了先验误差估计。     

含开边界二维Stokes问题的Galerkin边界元解法

本文推导了含有开边界的二维有限域上Stokes问题的边界积分方程, 得出基于单层位势的第一类间接边界积分方程.对与之等价的边界变分方程用Galerkin边界元求解以得出单层位势的向量密度. 对于含有开边界端点的边界单元,采用特别的插值函数, 以模拟其固有的奇异性.论文用若干数值算例模拟了含有开边界的有限区域上不可压缩粘性流体的绕流.       

弹性动力问题的有限元与特解边界元耦合的数值方法

导出了特解边界元法与有限元法的耦合方程。并应用自由度缩减技术,使耦合方程的自由度缩减到有限元域及其和边界元域的耦合边界上。这样得到的耦合方程不增加原有限元方程的带宽和阶数。耦合方程的求解可以引用求解有限元方程的所有方法,易于程序实现。数值算例结果表明,本文所提出的方法是正确的,是一种较为理想的耦合方法。     

理想弹塑性Ⅲ型扩展裂纹的全新和精确分析

本文采用线场分析方法对理想弹塑性Ⅲ型准静态扩展裂纹进行了分析.本文的意义在于突破了小范围屈服理论的限制.通过求得裂纹线附近塑性区应力和位移率的通解,并将此通解(而不是过去一直采用的特解)与弹性场的精确解(而不是线弹性奇异K场)在裂纹线附近的弹塑性边界上匹配,本文得出了裂纹线附近塑性区的应力变形场、塑性区的长度及弹塑性边界的单位法向量的全新和精确解答.本文的分析放弃了小范屈服理论的所有近似假定并且不再附加任何其它的近似假定,本文的结果在裂纹线附近是足够精确的.本文的结果表明:对理想弹塑性Ⅲ型准静态扩展裂纹,不存在“定常扩展状态”,且裂纹线附近塑性应变不存在奇异性.本文还对裂纹稳定扩展过程讨论了两种重要情形.     

Ⅱ型平面应力裂纹线场的弹塑性精确解

本文采用线场分析方法对理想弹塑性Ⅱ型平面应力裂纹裂纹线附近的应力场及弹塑性边界进行了精确分析。本文完全放弃了小范围屈服条件,探讨了弹塑性边界上弹塑性应力场匹配条件的正确提法,通过将裂纹线附近塑性区应力场的通解(而不是过去采用的特解)与弹性应力场的精确解(而不是通常的裂尖应力强度因子K场)在裂纹线附近的弹塑性边界上匹配,本文得出了塑性区应力场,塑性区长度及弹塑性边界的单位法向量在裂纹线附近的足够精确的表达式。     

对称迭层矩形板的平面应力分析

常用的对称迭层板为各向异性板．根据平面应力问题的基本方程精确地用应力函数解法求得了各向异性板的一般解析解．推导出平面内应力和位移的一般公式,其中积分常数由边界条件来决定．一般解包括三角函数和双曲函数组成的解,它能满足4个边为任意边界条件的问题．还有代数多项式解,它能满足4个角的边界条件．因此一般解可用以求解任意边界条件下的平面应力问题．以4边承受均匀法向和切向载荷以及非均匀法向载荷的对称迭层方板为例,进行了计算和分析．     

理想弹塑性I型平面应力裂纹线场的精确解

本文纠正了过去在裂纹弹塑性场匹配上存在的问题,采用线场分析方法,通过求得塑性区应力场的合理解答,使之与弹性精确场在裂纹线附近的弹塑性边界上匹配。本文就远场受单向拉伸及双向拉伸的理想弹塑性平面应力裂纹无限板,在完全放弃了小范围屈服条件的情况下求得了塑性区应力场、塑性区长度以及弹塑性边界的单位法向量在裂纹线附近足够精确的表达式。结果表明,无论单向拉伸和双向拉伸,塑性区应力分量σ_y,σ_xy,塑性区长度以及弹塑性边界的单位法向量在裂纹线附近的表达式完全相同,但塑性区沿X方向的正应力σ_x存在差别。     

任意形状孔口双边裂纹平板的应力强度因子计算

本文采用Muskhclishvili弹性力学的复变函数和边界配位方法对不同形状孔口双边裂纹问题进行了研究,计算了圆孔、椭圆孔、矩形孔、菱形孔等不同形状孔口双边裂纹,以及Ⅰ型和复合型等不同类型断裂试件的应力强度因子,本文方法简单方便,精度较高,与某些已有计算结果的问题比较,本文方法所得的结果是令人满意的.同时,本方法可以应用于不同几何形状和加载条件下的孔口双边裂纹有限大板的计算,是解这一类问题的一致有效方法.     

Boundary element analysis of uncoupled transient thermo-elastic problems with time- and space-dependent heat sources

A boundary element method (BEM) for the analysis of two- and three-dimensional uncoupled transient thermo-elastic problems involving time- and space-dependent heat sources is presented. The domain integrals are efficiently treated using the Cartesian transformation and the radial integration methods without considering any internal cells. Similar to the dual reciprocity method (DRM), some internal points without any connectivity are considered; however, in contrast to the DRM, any arbitrary mesh-free interpolation method can be used in the present formulation. There is no need to find any particular solutions and the shape functions in the mesh-free interpolation method can be arbitrary and sufficiently complicated. Unlike the DRM, the generated system of equations contains the unknowns only on the boundary. After finding the primary unknowns on the boundary, the temperature, displacement, and stress components at all internal points can directly be found without solving any system of equations. Three examples with different forms of heat sources are presented to demonstrate the efficiency and accuracy of the proposed method. Although the proposed BEM is mathematically more complicated than domain methods, such as the finite element method (FEM), it is more efficient from a modelling viewpoint since only the surface mesh has to be generated in the presented method.     

两连域保角映射成环域的方法

本文论述一种把型区域保角映射成环域的方法.其要点是将问题转化为Dirichlet问题,并证明该映像函数之实部应满足本文所示的边界条件,进而依据两连域上定义的调和函数的单值特性确定环域的内半径.映像函数的虚部可由Cauchy-Riemann条件得到,由此产生的积分常数仅影响映像点的幅角,并可由一一对应的映像来确定.不失其一般性,本方法可将由矩形拼成的复杂两连域保角映射成环域.笔者还对本方法作了电算,证明本方法可靠、经济、结果附有表格.     

A new coupled approach high accuracy numerical method for the solution of 3D non-linear biharmonic equations

In this paper, we derive a new fourth order finite difference approximation based on arithmetic average discretization for the solution of three-dimensional non-linear biharmonic partial differential equations on a 19-point compact stencil using coupled approach. The numerical solutions of unknown variable u(x,y,z) and its Laplacian ∇²u are obtained at each internal grid point. The resulting stencil algorithm is presented which can be used to solve many physical problems. The proposed method allows us to use the Dirichlet boundary conditions directly and there is no need to discretize the derivative boundary conditions near the boundary. We also show that special treatment is required to handle the boundary conditions. The new method is tested on three problems and the results are compared with the corresponding second order approximation, which we also discuss using coupled approach.     

A singular perturbation approach to integral equations occurring in poroelasticity

Singular perturbation theory is used to solve the integral equationswhich occur when treating finite-length crack problems in porouselastic materials. The method provides the stress intensityfactors which characterize the near crack tip stress and displacementfields for small times. The method also gives the stress andpore pressure fields on the fracture plane for small times relativeto the diffusive time scale. In this paper, the authors treatcrack problems which are unmixed in the pore pressure boundarycondition on the fracture plane. The Abelian result that smalltimes correspond, in Laplace transform space, to large valuesof the transform variable is used to formulate the problemsin terms of a small parameter. Rescaling on this small parameterleads to inner problems which are eigensolutions of the semi-infiniteproblems treated earlier by the authors. The outer solutionsare given by elastic eigensolutions together with appropriatefluid dipole responses. These outer solutions give the completestress and pore pressure fields except in the neighbourhoodof the crack tips; in this region the outer solutions are asymptoticallymatched with inner solutions. The full outer solutions are givenhere as an asymptotic expansion for small times and enable thedevelopment of the outer fields to be followed in real time.A reciprocal theorem in Laplace transform space is used to checkthe small-time solutions. The inner problem is rescaled to asemi-infinite crack problem, so eigensolutions of this semi-infiniteproblem are used together with the known asymptotic behaviourof the real solution to identify the stress intensity factor.The stress intensity factor is then related to an integral involvingthe inner limit of the outer solution together with the eigensolutionof the semi-infinite problem. Using this integral, we recoverthe result for the stress intensity factor found using singularperturbation theory. A ‘nearly’ invariant integralanalogous to the invariant M integral used in elastostaticsis derived. Unfortunately, the poroelastic analogue is not invariant,although it is used to verify the small-time results.