Some details of the Galerkin finite element method

The incorporation of the Galerkin technique in the finite element method has removed the constraint of finding a variational formulation for many problems of mathematical physics. The method has been successfully applied to many areas and has received wide acceptance. However, in the process of transplanting the concept from the Galerkin method for the entire domain to the Galerkin finite element method, some formal details have been overlooked or glossed over in the literature. This paper considers some of these details, including a possible reason for integration by parts and the contribution of interelement discontinuity terms.     

Numerical analysis of quasi-static crack branching in brittle solids by a modified displacement discontinuity method

Mechanism of quasi-static crack branching in brittle solids has been analyzed by a modified displacement discontinuity method. It has been assumed that the pre-existing cracks in brittle solids may propagate at the crack tips due to the initiation and propagation of the kink (or wing) cracks. The originated wing cracks will act as new cracks and can be further propagated from their tips according to the linear elastic fracture mechanics (LEFM) theory. The kink displacement discontinuity formulations (considering the linear and quadratic interpolation functions) are specially developed to calculate the displacement discontinuities for the left and right sides of a kink point so that the first and second mode kink stress intensity factors can be estimated. The crack tips are also treated by boundary displacement collocation technique considering the singularity variation of the displacements and stresses near the crack tip. The propagating direction of the secondary cracks can be predicted by using the maximum tangential stress criterion. An iterative algorithm is used to predict the crack propagating path assuming an incremental increase of the crack length in the predicted direction (straight and curved cracks have been treated). The same approach has been used for estimating the crack propagating direction and path of the original and wing cracks considering the special crack tip elements. Some example problems are numerically solved assuming quasi-static conditions. These results are compared with the corresponding experimental and numerical results given in the literature. This comparison validates the accuracy and applicability of the proposed method.     

弱不连续问题扩展有限元法的数值精度研究

主要研究了扩展有限元法(extended finite element method, XFEM)在处理弱不连续问题时不同改进函数形式对XFEM数值求解精度的影响,阐述了各种改进函数影响XFEM求解精度的关键因素,指出校正的扩展有限元法(corrected-XFEM)能够提高数值求解精度的实质在于它拓展了改进结点域,即将常规扩展有限元法(standard-XFEM)的改进结点域增加一层作为corrected-XFEM的改进结点域,文中建议延拓corrected-XFEM的改进结点域,即在corrected-XFEM的改进结点域基础上再增加一层改进结点. 利用水平集函数表征材料内部的不连续界面,推导了XFEM求解的支配方程,给出了一种改进单元的数值积分方案以及改进单元处高精度应力的求解方法. 含夹杂问题的数值计算结果表明:建议的延拓corrected-XFEM改进结点域的方法能够明显提高XFEM的数值求解精度.      

A duality recipe for non-convex variational problems

The aim of this paper is to present a general convexification recipe that can be useful for studying non-convex variational problems. In particular, this allows us to treat such problems by using a powerful primal–dual scheme. Possible further developments and open issues are given.     

Regularizing algorithms for detecting discontinuities in ill-posed problems

The problem of detecting singularities (discontinuities of the first kind) of a noisy function in L ₂ is considered. A wide class of regularizing algorithms that can detect discontinuities is constructed. New estimates of accuracy of determining the location of discontinuities are obtained and their optimality in terms of order with respect to the error level δ is proved for some classes of functions with isolated singularities. New upper bounds for the singularity separation threshold are obtained.     

A finite strain, finite band method for modeling ductile fracture

We present a finite deformation generalization of the finite thickness embedded discontinuity formulation presented in our previous paper [A.E. Huespe, A. Needleman, J. Oliver, P.J. Sánchez, A finite thickness band method for ductile fracture analysis, Int. J. Plasticity 25 (2009) 2349-2365]. In this framework the transition from a weak discontinuity to a strong discontinuity can occur using a single constitutive relation which is of importance in a range of applications, in particular ductile fracture, where localization typically precedes the creation of new free surface. An embedded weak discontinuity is introduced when the loss of ellipticity condition is met. The resulting localized deformation band is given a specified thickness which introduces a length scale thus providing a regularization of the post-localization response. The methodology is illustrated through several example problems emphasizing finite deformation effects including the development of a cup-cone failure in round bar tension.     

非线性兰姆波在阶梯板中传播的有限元仿真

在非线性兰姆波检测中,累积二次谐波是重要的检测信号.该文利用有限元法仿真了非线性兰姆波在阶梯板中的产生和传播过程,通过二维离散傅里叶变换分析了兰姆波的传播模态和各个模态通过波导横截面的能流大小,给出了各个模态在阶梯板中的反射、透射系数,总结了其随阶梯板厚度的变化规律.研究发现,随着台阶增高,散射效应使得兰姆波更难入射波...