改进型扩展比例边界有限元法

结合了扩展有限元法(extended finite elementmethods,XFEM)和比例边界有限元法(scaled boundary finite elementmethods,SBFEM)的主要优点,提出了一种改进型扩展比例边界有限元法(improvedextended scaled boundary finite elementmethods,$i$XSBFEM),为断裂问题模拟提供了一条新的途径.类似XFEM,采用两个正交的水平集函数表征材料内部裂纹面,并基于水平集函数判断单元切割类型;将被裂纹切割的单元作为SBFE的子域处理,采用SBFEM求解单元刚度矩阵,从而避免了XFEM中求解不连续单元刚度矩阵需要进一步进行单元子划分的缺陷;同时,借助XFEM的主要思想,将裂纹与单元边界交点的真实位移作为单元结点的附加自由度考虑,赋予了单元结点附加自由度明确的物理意义,可以直接根据位移求解结果得出裂纹与单元边界交点的位移;对于含有裂尖的单元,选取围绕裂尖单元一圈的若干层单元作为超级单元,并将此超级单元作为SBFE的一个子域求解刚度矩阵,超级单元内部的结点位移可通过SBFE的位移模式求解得到,应力强度因子可基于裂尖处的奇异位移(应力)直接获得,无需借助其他的数值方法.最后,通过若干数值算例验证了建议的$i$XSBFEM的有效性,相比于常规XFEM,$i$XSBFEM的基于位移范数的相对误差收敛性较好;采用$i$XSBFEM通过应力法和位移法直接计算得到的裂尖应力强度因子均与解析解吻合\较好.      

裂纹面受荷载作用的应力强度因子的计算

基于比例边界有限元法计算了裂纹面有荷载作用情况下裂纹尖端的应力强度因子,给出了有限介质裂纹面作用荷载的比例边界有限元方程的基本求解过程.对于随径向坐标任意变化的一类面荷载的积分能够显式计算,不需要引入额外的近似;并将计算结果与解析解和数值结果进行对比,结果表明比例边界有限元法在计算裂纹面作用荷载时的应力强度因子是有效且精确的.此外,该方法可方便地处理各向异性材料裂纹问题,本文给出了正交各向异性矩形盘裂纹面受均布荷载情况的应力强度因子.     

Analysis method of planar cracks of arbitrary shape in the isotropic plane of a three-dimensional transversely isotropic magnetoelectroelastic medium

The hyper-singular boundary integral equation method of crack analysis in three-dimensional transversely isotropic magnetoelectroelastic media is proposed. Based on the fundamental solutions or Green’s functions of three-dimensional transversely isotropic magnetoelectroelastic media and the corresponding Somigliana identity, the boundary integral equations for a planar crack of arbitrary shape in the plane of isotropy are obtained in terms of the extended displacement discontinuities across crack faces. The extended displacement discontinuities include the displacement discontinuities, the electric potential discontinuity and the magnetic potential discontinuity, and correspondingly the extended tractions on crack face represent the conventional tractions, the electric displacement and the magnetic induction boundary values. The near crack tip fields and the intensity factors in terms of the extended displacement discontinuities are derived by boundary integral equation approach. A solution method is proposed by use of the analogy between the boundary integral equations of the magnetoelectroelastic media and the purely elastic materials. The influence of different electric and magnetic boundary conditions, i.e., electrically and magnetically impermeable and permeable conditions, electrically impermeable and magnetically permeable condition, and electrically permeable and magnetically impermeable condition, on the solutions is studied. The crack opening model is proposed to consider the real crack opening and the electric and magnetic fields in the crack cavity under combined mechanical-electric-magnetic loadings. An iteration approach is presented for the solution of the non-linear model. The exact solution is obtained for the case of uniformly applied loadings on the crack faces. Numerical results for a square crack under different electric and magnetic boundary conditions are displayed to demonstrate the proposed method.     

MODE III 2-D FRACTURE ANALYSIS BY THE SCALED BOUNDARY FINITE ELEMENT METHOD

The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique that combines the advantages of the finite element method and the boundary element method with unique properties of its own. This method has proven very efficient and accurate for determining the stress intensity factors （SIFs） for mode I and mode II two-dimensional crack problems. One main reason is that the SBFEM has a unique capacity of analytically representing the stress singularities at the crack tip. In this paper the SBFEM is developed for mode III （out of plane deformation） two-dimensional fracture anMysis. In addition, cubic B-spline functions are employed in this paper for constructing the shape functions in the circumferential direction so that higher continuity between elements is obtained. Numerical examples are presented at the end to demonstrate the simplicity and accuracy of the present approach for mode Ⅲ two-dimensional fracture analysis.     

一种XFEM断裂分析的裂尖单元新型改进函数

提出了一种适用于裂尖改进单元的新型改进函数, 基于三角变换的方法, 保留裂纹尖端场的应力奇异性和裂纹上、下表面的位移不连续性, 将常规扩展有限元法裂尖改进单元的4 项改进函数缩减为2 项, 裂尖改进单元的结点由常规的8 个改进自由度减少为4 个. 采用2 个正交的水平集函数表征材料内部裂纹面, 详细阐述了改进单元类型的判别方法, 给出一种改进单元的分区域积分方案. 最后, 若干断裂力学问题经典算例的数值计算结果表明:建议的裂尖改进函数具有较高的数值精度, 该方法是十分有效的.     

Crack at the apex of a loaded notch

This paper deals with the plane elastostatic problem of an infinite wedge, subjected to arbitrary surface tractions, and cracked along the wedge angle bisector. The problem is reduced to a single Fredholm integral equation, which is solved numerically for normal loads on the crack faces and various loads on the wedge faces. It is shown that the crack tip intensity factor depends strongly on the wedge angle. An approximation to a half plane with a notch of finite angle, cracked at its apex, is also obtained.     

双轴载荷作用下源于椭圆孔的分支裂纹的一种边界元分析

利用一种边界元方法来研究双轴载荷作用下无限大板中源于椭圆孔的分支裂纹．该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者提出的裂尖位移不连续单元构成．在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界,文中算例说明本数值方法对计算平面弹性裂纹的应力强度因子是非常有效的。该文对双轴载荷作用下无限大板中源于椭圆孔的分支裂纹的数值结果进一步证实本数值方法对计算复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示双轴载荷及裂纹体几何对应力强度因子的影响。     

On the crack face boundary conditions in electromechanical fracture and an experimental protocol for determining energy release rates

The fracture mechanics of electromechanical materials has been investigated for well over a decade, yet there still exists controversy over the appropriate crack face boundary conditions for non-conducting cracks. In this paper an experimental protocol for measuring the energy release rate in a non-linear reversible electromechanical body is proposed and summarized. The potential results from the proposed experimental approach are capable of shedding light on the true physical nature of the conditions prevailing at the crack surface and in the space within the crack. The experimental procedure is simulated numerically for a linear piezoelectric specimen in a four point bending configuration subjected to electrical loading perpendicular to the crack. The focus of these investigations is on a comparison between the commonly used exact crack face boundary condition and the recently proposed energetically consistent boundary conditions. To perform the numerical calculation with a wide range of electrical and mechanical loadings, two efficient finite element formulations are presented for the general analysis of crack problems with non-linear crack face boundary conditions. Methods for the numerical determination of the crack tip energy release rate and the simulation of the experimental method for obtaining the total energy release rate are developed. Numerical results for the crack tip and total energy release rate are given for both the exact and energetically consistent boundary conditions. It is shown that the crack tip energy release rate calculated under energetically consistent boundary conditions is equal to the total energy release rate generated from the simulated experimental method. When the exact boundary conditions are used, there is no such agreement.     

An effective boundary element method for analysis of crack problems in a plane elastic plate

A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.     

Reduced augmented Lagrangian bi-conjugate gradient method for impact-contact problems

To avoid the numerical oscillation of the penalty method and non-compatibility with explicit operators of conventional Lagrange multiplier methods used in transient contact problems to enforce surface contact conditions, a new approach to enforcing surface contact constraints for the transient nonlinear finite element problems, referred to as “the reduced augmented Lagrangian bi-conjugate gradient method (ALCG)”, is developed in this paper. Based on the nonlinear constrained optimization theory and is compatible with the explicit time integration scheme, this approach can also be used in implicit scheme naturally. The new surface contact constraint method presented has significant advantages over the widely adopted penalty function methods and the conventional Lagrangian multiplier methods. The surface contact constraints are satisfied more accurately for each step by the algorithm, so the oscillation of numerical solution for the explicit scheme is depressed. Through the development of new iteration strategy for solving nonlinear equations, ALCG method improves the computational efficiency greatly. Project supported by State Education Commission Doctoral Foundation and Natural Science Foundation of Liaoning Province.     

裂缝分析的比例边界有限元与有限元耦合的虚拟结构面模型

比例边界有限元法SBFEM(Scaled Boundary Finite Element Method)是一种半解析数值方法,在裂缝分析特别是强度因子计算上具有相当高的精度。本文提出了一种用于裂缝分析的基于虚拟结构面的SBFEM与常规FEM的耦合分析方法。首先选取裂缝周边一定范围的计算域,并将结构分成不含裂缝区域和含裂缝区域两部分。然后,对不含裂缝区域,采用FEM进行网格离散;对含裂缝区域,采用SBFEM进行网格离散;两者相互独立,在这两个域内,分别采用各自相应的位移模式。最后通过在SBFEM网格的外边界设置虚拟耦合结构面的模式,实现有限元网格和比例边界有限元网格的耦合。通过两个经典的含裂缝平板的算例研究,探讨了本文方法在I型开裂和混合型开裂分析中,影响应力强度因子精度的因素。算例表明,SBFEM具有的降维和半解析性质,使本文方法在裂缝分析中的前处理简单易行,且计算结果具有相当高的计算精度。     

Analysis of electric boundary condition effects on crack propagation in piezoelectric ceramics

There are three types of cracks: impermeable crack, permeable crack and conducting crack, with different electric boundary conditions on faces of cracks in piezoelectric ceramics, which poses difficulties in the analysis of piezoelectric fracture problems. In this paper, in contrast to our previous FEM formulation, the numerical analysis is based on the used of exact electric boundary conditions at the crack faces, thus the common assumption of electric impermeability in the FEM analysis is avoided. The crack behavior and elasto-electric fields near a crack tip in a PZT-5 piezoelectric ceramic under mechanical, electrical and coupled mechanical-electrical loads with different electric boundary conditions on crack faces are investigated. It is found that the dielectric medium between the crack faces will reduce the singularity of stress and electric displacement. Furthermore, when the permittivity of the dielectric medium in the crack gap is of the same order as that of the piezoelectric ceramic, the crack becomes a conducting crack, the applied electric field has no effect on the crack propagation. The project supported by the National Natural Science Foundation of China (19672026, 19891180)     

Linear and Nonlinear Algorithms for Stress Separation in Photoelasticity

An experimental-numerical hybrid method for the stress separation in photoelasticity is proposed in this study. In the proposed method, boundary conditions for a local finite element model, that is, tractions along boundaries are inversely determined from photoelastic fringes. Two algorithms are proposed for determining the boundary condition. One is a linear algorithm in which the tractions are obtained by the method of linear least-squares from both principal stress difference and principal direction. Another is the nonlinear algorithm in which the tractions are determined only from the principal stress difference. After determining the boundary conditions for the local finite element model, the stresses can be obtained by finite element direct analysis. The effectiveness is demonstrated by applying the proposed method to a perforated plate under tension and contact problems. Results show that the boundary conditions of the local finite element model can be determined from the photoelastic fringes and then the individual stresses can be obtained by the proposed method. Furthermore, the stresses can be evaluated even if the boundary condition is complicated such as at the contact surface. It is expected that the proposed method can be powerful tool for stress analysis.     

Influence of non-singular stress terms and specimen geometry on small-scale yielding at crack tips in elastic-plastic materials

The plane strain elastic-plastic state at a crack tip is determined for compact tension, bend, double edge-cracked and centre-cracked specimens using a finite element method with triangular constant-strain elements. The solutions are found to differ by 10 to 30 per cent at the ASTM-limit as regards fracture surface displacement, normal stress and plastic zone size. In order to bring the boundary layer solution for the crack problem into agreement with the solution for a specific specimen one has to modify this solution. The modification consists of an addition to the boundary tractions for the boundary layer problem of tractions corresponding to the non-singular, constant second term in a series expansion of the normal stress parallel to the crack plane.     

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method (SBFEM) is a recently developed numerical method combining advantages of both finite element methods (FEM) and boundary element methods (BEM) and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors (SIFs) directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.The project supported by the National Natural Science Foundation of China (50579081) and the Australian Research Council (DP0452681)The English text was polished by Keren Wang.     

2D dynamic analysis of cracks and interface cracks in piezoelectric composites using the SBFEM

The dynamic stress and electric displacement intensity factors of impermeable cracks in homogeneous piezoelectric materials and interface cracks in piezoelectric bimaterials are evaluated by extending the scaled boundary finite element method (SBFEM). In this method, a piezoelectric plate is divided into polygons. Each polygon is treated as a scaled boundary finite element subdomain. Only the boundaries of the subdomains need to be discretized with line elements. The dynamic properties of a subdomain are represented by the high order stiffness and mass matrices obtained from a continued fraction solution, which is able to represent the high frequency response with only 3–4 terms per wavelength. The semi-analytical solutions model singular stress and electric displacement fields in the vicinity of crack tips accurately and efficiently. The dynamic stress and electric displacement intensity factors are evaluated directly from the scaled boundary finite element solutions. No asymptotic solution, local mesh refinement or other special treatments around a crack tip are required. Numerical examples are presented to verify the proposed technique with the analytical solutions and the results from the literature. The present results highlight the accuracy, simplicity and efficiency of the proposed technique.     

A contact force model in the dynamic analysis of contactable cracks

In this paper a finite element method in the time domain is proposed to analyze the dynamic response of contactable cracks. By this method, the contact normal stress caused by impact or contact of the crack faces can be directly obtained and then the motion of the model containing the contactable cracks can be calculated, and the processes of contact, slip and separation of the crack faces can also be correctly determined and simulated. The method proposed can be used to study the scattering problems of an elastic wave caused by the crack and the dynamic extension problems of the crack arising from impact or cyclic loading. Numerical examples indicate that the method is efficient and accurate enough. It is also shown that the effect of dynamic contact of crack faces on the motion of itself and the area near the crack is significant by comparing the results with the contact considered and neglected.     

A numerical analysis of cracks emanating from an elliptical hole in a 2-D elasticity plate

This paper is concerned with the stress intensity factors (SIFs) of cracks emanating from an elliptical hole in an infinite or a finite plate under biaxial loads by using a boundary element method, which consists of the non-singular displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements due to the author. In the boundary element implementation the left or the right crack-tip element is placed locally at the corresponding left or right crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. A few numerical examples are included to show that the present approach is very efficient and accurate for the calculating the SIFs of crack problems in an infinite or a finite plate. The present numerical results of cracks emanating from an elliptical hole under biaxial loads can reveal the effect of the elliptical aspect ratio and the transverse load on the SIFs.     

基于SBFEM和深度学习的裂纹状缺陷反演模型

结构内部缺陷的识别是结构健康监测的重要研究内容, 而当前以无损检测为主的结构安全检测多以定性分析为主, 定量识别缺陷的尺度较困难. 本文将比例边界有限元法(scaled boundary finite element methods, SBFEM)和深度学习相结合, 提出了基于Lamb波在结构中传播时的反馈信号定量识别结构内部裂纹状缺陷的反演模型. 通过随机生成缺陷信息(位置、大小), 采用SBFEM模拟Lamb波在含不同缺陷信息的结构中的信号传播过程, SBFEM仅需对结构边界离散可最小化网格重划分过程, 大大提高了计算效率. Lamb波在含裂纹状缺陷结构中传播时观测点的反馈信号包含大量的裂纹信息, 基于这一特性可为深度学习模型提供足够多的反映问题特性的训练数据. 建议的缺陷反演模型规避了传统反分析问题的目标函数极小化迭代过程, 在保证计算精度的前提下大大减少了计算成本. 对含单裂纹和多裂纹板的数值算例进行分析, 结果表明: 建立的缺陷识别模型能够准确地量化结构内部的缺陷, 对浅表裂纹亦有很好的识别效果, 且对于含噪信号模型仍具有较好的鲁棒性.      

THE STRESS FIELD AT MODE I CRACK TIP ON STEEL PIPE UNDER AXIAL TENSION

采用弹性断裂力学Westergaard 的方法,分析在轴向拉伸载荷作用下钢管的Ⅰ型裂纹尖端处的应力场. 基于Ⅰ型裂纹钢管应力场的特征,设定尖端处的应力艾雷函数,给出其应力场模型边界条件和应力场的解析函数,并利用裂纹尖端处的切平面研究裂纹尖端局部应力场,建立了钢管裂纹尖端应力场模型. 通过钢管与平板Ⅰ型裂纹应力场模型的对比,结果表明二者明显不同,钢管裂纹尖端处应力峰值影响范围仅与裂纹长度、拉应力相关.