A boundary element analysis for stress intensity factors of multiple circular arc cracks in a plane elasticity plate

In this paper, a numerical approach for analyzing interacting multiple cracks in infinite linear elastic media is presented. By extending Bueckner’s principle suited for a crack to a general system containing multiple interacting cracks, the original problem is divided into a homogeneous problem (the one without cracks) subjected to remote loads and a multiple crack problem in an unloaded body with applied tractions on the crack surfaces. Thus, the results in terms of the stress intensity factors (SIFs) can be obtained by considering the latter problem, which is analyzed easily by means of the displacement discontinuity method with crack-tip elements proposed recently by the author. Test examples are given to illustrate that the numerical approach is very accurate for analyzing interacting multiple cracks in an infinite linear elastic media under remote uniform stresses. In addition, the displacement discontinuity method with crack-tip elements is used to analyze a multiple crack problem in a finite plate. It is found that the boundary element method is also very accurate for investigating interacting multiple cracks in a finite plate. Specially, a generalization of Bueckner’s principle and the displacement discontinuity method with crack-tip elements are used to analyze multiple circular arc crack problems in infinite plate in tension (including: Two Collinear Circular Arc Cracks, Three Collinear Circular Arc Cracks, Two Parallel Circular Arc Cracks, Three Parallel Circular Arc Cracks and Two Circular Arc Cracks) in a plane elasticity plate. Many results are given.     

平面弹性裂纹分析的一种有效边界元方法

提出了一种简单而有效的平面弹性裂纹应力强度因子的边界元计算方法．该方法由Crouch与Starfield建立的常位移不连续单元和闫相桥最近提出的裂尖位移不连续单元构成A·D2在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界．算例(如单向拉伸无限大板中心裂纹、单向拉伸无限大板中圆孔与裂纹的作用)说明平面弹性裂纹应力强度因子的边界元计算方法是非常有效的．此外,还对双轴载荷作用下有限大板中方孔分支裂纹进行了分析．这一数值结果说明平面弹性裂纹应力强度因子的边界元计算方法对有限体中复杂裂纹的有效性,可以揭示双轴载荷及裂纹体几何对应力强度因子的影响．     

Simulation of short crack propagation using a hybrid boundary element technique

Service life of cyclically loaded components is often determined by the propagation of short fatigue cracks, which is highly influenced by microstructural features such as grain boundaries. A two-dimensional model to simulate the growth of such stage I-cracks is presented. The crack is discretised by dislocation discontinuity boundary elements and the direct boundary element method is used to mesh the grain boundaries. A superposition procedure couples these different boundary element methods to employ them in one model. Varying elastic properties of the grains are considered and their influence on short crack propagation is studied. A change in crack tip slide displacement determining short crack propagation is observed. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

3D modeling of crack growth in electro-magneto-thermo-elastic coupled viscoplastic multiphase composites

This work presents a time-domain hypersingular integral equation (TD-HIE) method for modeling 3D crack growth in electro-magneto-thermo-elastic coupled viscoplastic multiphase composites (EMTE-CVP-MCs) under extended incremental loads rate through intricate theoretical analysis and numerical simulations. Using Green’s functions, the extended general incremental displacement rate solutions are obtained by time-domain boundary element method. Three-dimensional arbitrary crack growth problem in EMTE-CVP-MCs is reduced to solving a set of TD-HIEs coupled with boundary integral equations, in which the unknown functions are the extended incremental displacement discontinuities gradient. Then, the behavior of the extended incremental displacement discontinuities gradient around the crack front terminating at the interface is analyzed by the time-domain main-part analysis method of TD-HIE. Also, analytical solutions of the extended singular incremental stresses gradient and extended incremental integral near the crack fronts in EMTE-CVP-MCs are provided. In addition, a numerical method of the TD-HIE for a 3D crack subjected to extended incremental loads rate is put forward with the extended incremental displacement discontinuities gradient approximated by the product of time-domain basic density functions and polynomials. Finally, examples are presented to demonstrate the application of the proposed method.     

Three-dimensional finite/infinite elements analysis of fluid flow in porous media

Fluid flow in naturally fractured porous media can always be regarded as an unbounded domain problem and be better solved by finite/infinite elements. In this paper, a three-dimensional two-direction mapped infinite element is generated and combined with conventional finite elements and one direction infinite element to simulate poroelasticity. Therefore, the entire semi-infinite domain can be included in the numerical analysis. Both single- and dual-porosity porous media are considered. For the purpose of validation, we compare the results of finite/infinite elements with those of finite elements under two extreme boundary conditions. The comparison indicated that mapped infinite element is an appropriate approach to model fluid flow in porous media and provides an intermediate solution.     

Numerical electroseismic modeling: A finite element approach

Electroseismics is a procedure that uses the conversion of electromagnetic to seismic waves in a fluid-saturated porous rock due to the electrokinetic phenomenon. This work presents a collection of continuous and discrete time finite element procedures for electroseismic modeling in poroelastic fluid-saturated media. The model involves the simultaneous solution of Biot’s equations of motion and Maxwell’s equations in a bounded domain, coupled via an electrokinetic coefficient, with appropriate initial conditions and employing absorbing boundary conditions at the artificial boundaries. The 3D case is formulated and analyzed in detail including results on the existence and uniqueness of the solution of the initial boundary value problem. Apriori error estimates for a continuous-time finite element procedure based on parallelepiped elements are derived, with Maxwell’s equations discretized in space using the lowest order mixed finite element spaces of Nédélec, while for Biot’s equations a nonconforming element for each component of the solid displacement vector and the vector part of the Raviart-Thomas-Nédélec of zero order for the fluid displacement vector are employed. A fully implicit discrete-time finite element method is also defined and its stability is demonstrated. The results are also extended to the case of tetrahedral elements. The 2D cases of compressional and vertically polarized shear waves coupled with the transverse magnetic polarization (PSVTM-mode) and horizontally polarized shear waves coupled with the transverse electric polarization (SHTE-mode) are also formulated and the corresponding finite element spaces are defined. The 1D SHTE initial boundary value problem is also formulated and approximately solved using a discrete-time finite element procedure, which was implemented to obtain the numerical examples presented.     

分析表面裂纹的一种新方法

本文作者综合了线弹簧模型及边界元法的优点,开发了一种新的线弹簧边界元法.该方法把表面裂纹这一三维问题简化为拟一维问题,可用于分析受到多种载荷作用的含表面裂纹的板.本文对该方法进行了理论分析和数值验证,报告了计算结果.结果表明,该方法经济有效.利用该方法仅使用个人计算机就可以分析表面裂纹问题.     

A strong discontinuity based adaptive refinement approach for the modeling of crack branching

In the current work, the physical phenomena of dynamic fracture of brittle materials involving crack growth, acceleration and consequent branching is simulated. The numerical modeling is based on the approach where the failure in the form of cracks or shear bands is modeled by a jump in the displacement field, the so called ‘strong discontinuity’. The finite element method is employed with this strong discontinuity approach where each finite element is capable of developing a strong discontinuity locally embedded into it. The focus in this work is on branching phenomena which is modeled by an adaptive refinement method by solving a new sub-boundary value problem represented by a finite element at the growing crack tip. The sub-boundary value problem is subjected to a certain kinematic constraint on the boundary in the form of a linear deformation constraint. An accurate resolution of the state of material at the branching crack tip is achieved which results in realistic dynamic fracture simulations. A comparison of resulting numerical simulations is provided with the experiment of dynamic fracture from the literature. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

无穷凹角区域各向异性问题的自然边界元与有限元耦合法

本文研究无穷凹角区域上一类各向异性问题的自然边界元与有限元耦合法.利用自然边界归化原理,获得圆弧或椭圆弧人工边界上的自然积分方程,给出了耦合的变分形式及其数值方法,以及逼近解的收敛性和误差估计,最后给出了数值例子,以示方法的可行性和有效性.     

一种新的声椭球无限单元

提出了一种新的声椭球无限单元．这种声无限单元基于一种新的声压表达式,这种声压表达式能够更准确地代表着椭球声场的声传播模式．这种新方法的形函数类似于Burnett方法,而权函数定义为形函数和一个附加因子的乘积．因为仅需要一维的数值积分,这种新方法的代码生成十分容易,就像处理一维单元一样．耦合标准的有限元程序,这种声无限单元理论上能够高效地求解任何形状的声源的声辐射和声散射现象．简要地推导了这种新方法,并给出了这种方法详尽的推导结果．为更有效地检验该无限元方法的可行性,文中例子仅考虑无限元求解的精度,而不包括相应的有限元．使用这种新方法,精确地推导出了摆动球的理论计算公式．而长旋转椭球的例子则表明了这种方法优于边界元方法和其他声椭球无限元方法．这些例子表明了这种新方法是切实可行的．     

描述低周疲劳裂纹扩展速率的循环J积分新参量

探讨了低周疲劳加载条件下的应力增量．应变增量关系，提出了模拟裂纹疲劳扩展的二维模型以建立新的循环．积分参量，详细阐述了该积分参量的定义、主要特点、物理意义以及数值计算方法，并通过紧凑拉伸试样的疲劳试验检验该积分参量的有效性．结果表明：该积分参量能够较好描述恒幅低周疲劳裂纹的扩展速率．此外，基于积分参量体系，从能量的角度解释了疲劳迟滞现象．     

On stresses in a strip with a crack

The solution of the problem of a loaded crack in an infinite strip is given using the method of superposition of three problems (a loaded crack in the infinite plane; an infinite homogeneous strip with normal and tangent stresses that are given on nonhomogeneous boundaries; an infinite strip with longitudinal generators which are free from load and an arbitrary load at the end), which makes it possible to satisfy the boundary conditions exactly.Translated from Dinamicheskie Sistemy, No. 9, pp. 65–71, 1990.     

Some integral relations of Hankel transform type and applications to elasticity theory

The dominant part of an integral equation arising in connection with boundary value problems for the circular disc is evaluated in terms of orthogonal polynomials. This relation leads to an efficient method for numerical solution of the complete integral equation even in the presence of a complicated bounded kernel. The static problem of a circular crack in an infinite elastic body under general loads is used to illustrate vector boundary conditions leading to two coupled integral equations, while the problem of a vibrating flexible circular plate in frictionless contact with an elastic half space is solved by use of the associated numerical method.     

Boundary collocation method for a cracked rectangular plate with double external tension

In this article, the boundary collocation method is employed to investigate the problems of a central crack in a rectangular plate which applied double external tension on the outer boundary under the assumption that the dimensions of the plate are much larger than that of the crack. A set of stress functions has also been proposed based on the theoretical analysis which satisfies the condition that there is no external force on the crack surfaces. It is only necessary to consider the condition on the external boundary. Using boundary collocation method, the linear algebra equations at collocation points are obtained. The least squares method is used to obtain the solution of the equations, so that the unknown coefficients can be obtained. According to the expression of the stress intensity factor at crack tip, we can obtain the numerical results of stress intensity factor. Numerical experiments show that the results coincide with the exact solution of the infinite plate. In particular, this case of the double external tension applied on the outer boundary is seldom studied by boundary collocation method.     

有一条裂纹的圆形焊接问题

本文讨论了在带圆孔的无限平面中焊接一个不同材料的带裂纹的近似圆板的问题.该问题化为求解解析函数边值问题然后又转化为求解沿裂纹的奇异积分方程.后者的数字解法也已给出.文末并对Ⅰ型、Ⅱ型情况得出了应力强度因子的公式以及数字结果.     

Efficient simulation of 3D fatigue crack growth based on boundary elements

The three-dimensional simulation of fatigue crack growth under the consideration of 3D effects is presented to follow complex crack paths as realistic as possible from the macroscopic point of view. This aim is mainly based on two aspects. Firstly, an accurate stress analysis is performed by the boundary element method in terms of the 3D Dual BEM. Secondly, a suitable 3D crack growth criterion based on experimental observations is utilized. Due to the non-linearity of crack growth the whole procedure is embedded in an advanced incremental crack growth algorithm. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Probabilistic modeling and simulation of multiple surface crack propagation and coalescence

In the low cycle fatigue (LCF) regime, fatigue failure of metallic materials with high strength and less impurities generally dominates by multiple surface crack propagation and coalescence, in which its final failure shows a stochastic nature on crack initiation, propagation and coalescence under cyclic loadings. According to this, the competing failure modes of multiple surface cracks and interior cracks are studied through coupling numerical simulations with fracture mechanics methods. In particular, a probabilistic procedure for modeling multiple surface crack propagation and coalescence is established by incorporating Monte Carlo simulation with experimental evidences, including surface crack density and crack length distributions measured from LCF replica tests of 30NiCrMoV12 steel. In addition, it calculates the probability of coalescence of neighboring cracks with allowance for their interactions and local plastic deformation at the crack tips. Finally, it estimates the remaining usage lives of specimens from initial state to critical cracks by propagation and coalescence of dispersed cracks.     

A stable boundary elements method for magnetohydrodynamic channel flows at high Hartmann numbers

The article is devoted to extension of boundary element method (BEM) for solving coupled equations in velocity and induced magnetic field for time dependent magnetohydrodynamic (MHD) flows through a rectangular pipe. The BEM is equipped with finite difference approach to solve MHD problem at high Hartmann numbers up to 10⁶. In fact, the finite difference approach is used to approximate partial derivatives of unknown functions at boundary points respect to outward normal vector. It yields a numerical method with no singular boundary integrals. Besides, a new approach is suggested in this article where transforms 2D singular BEM's integrals to 1D nonsingular ones. The new approach reduces computational cost, significantly. Note that the stability of the numerical scheme is proved mathematically when computational domain is discretized uniformly and Hartmann number is 40 times bigger than length of boundary elements. Numerical examples show behavior of velocity and induced magnetic field across the sections.