基于EEP超收敛解的自适应有限元法特性分析

基于EEP （单元能量投影）超收敛计算的自适应有限元法,已对一系列问题取得成功,但其自适应特性尚缺乏相关研究。该文以二阶常微分方程为模型问题,同时考察基于EEP和SPR （超收敛分片恢复）超收敛解的自适应分析方法,与有限元最优网格进行了比较分析,进而提出反映自适应有限元收敛特性的估计式,并给出了自适应收敛率β的定义。该文给出的数值试验表明:采用m次单元,对于解答光滑的问题,SPR法与EEP法均可有效用于自适应求解,其位移可按最大模获得m+1的自适应收敛率;对于奇异因子为α（<1）的奇异问题,SPR法失效,而基于EEP法的自适应求解,其位移按最大模可获得m+α的自适应收敛率,远高于α的常规有限元收敛率。     

变截面变曲率梁振型的有限元超收敛拼片恢复解和网格自适应分析

该文提出变截面变曲率梁振型的有限元后处理超收敛拼片恢复方法,建立各阶振型的超收敛解,并基于振型超收敛解进行变截面曲梁面内和面外自由振动的自适应分析。在位移型有限元后处理阶段,引入超收敛拼片恢复方法和高阶形函数插值技术,得到振型(位移)的超收敛解。利用振型超收敛解估计当前网格下振型有限元解的能量模形式下的误差,并指导网格进行自适应细分加密分析,获得优化的网格和满足预设误差限的高精度解答。数值算例表明该算法适于求解不同曲线形态、多类边界条件、变截面、变曲率形式的曲梁面内和面外自由振动连续阶频率和振型,解答精确、分析过程高效可靠。     

基于应力超收敛恢复技术的广义特征值问题后验误差估计

根据有限元解的超收敛特性提出了一种基于应力超收敛恢复技术的广义特征值问题后验误差估计。通过对单元内的应力超收敛点以及相邻单元的应力超收敛点进行插值或外推处理,得到单元内其它点处更高精度的应力解。通过高精度的应力值可以得到结构处理后改进的势能。将改进的势能代入瑞利商,最终得到比原始有限元解更高精度的特征解。将后处理特征解作为“准精确解”代替误差估计因子中未知的精确解,实现后验误差估计过程。数值计算结果表明,所提出的后验误差估计是渐进精确的,因此可作为结构广义特征值问题自适应有限元方法的误差估计因子。     

基于p型超收敛计算的一维有限元自适应分析

基于提高单元阶次的p型超收敛算法,可以在有限元解答基础上求得超收敛解。用该超收敛解代替精确解可以对有限元解答进行可靠的误差估计。对Zienkiewicz网格划分策略进行一定的改进,得到一种更有效的网格划分策略。基于可靠的误差估计和高效的网格划分,可以进行有限元自适应求解。数值试验表明,该文的自适应求解方案能够得到较优的网格和满足误差限的解答。     

流体固体动力耦合分析的有限元法

应用有限元法探讨了流体、固体接触界面由无限接触点对组成,并以接触点对的瞬态接触内力作为待定变量的流体固体动力耦合模型的数值求解方法。分析了流体、固体域插值函数的特点,用二维八节点等参元及流固接触面上的接触点对单元,对流固耦合系统进行了离散化处理;并采用变分原理推导了反映流体固体动力相互作用机理的接触约束矩阵(或称动力耦合矩阵),建立了有限元控制方程,给出了完整的数值计算方法,研编了动力耦合系统的分析程序。数值计算结果与经典理论解误差很小,验证了动力耦合模型和有限元求解方法的正确性及其较高的计算精度。     

横观各向同性岩石力学-化学-损伤耦合有限元分析

该文基于层理性岩石的水化作用和连续损伤特性,建立力学-化学-损伤耦合的有限元（FEM）求解方法,开展含钻井孔岩石的井壁应力和围岩损伤演化分析。该文发展横观各向同性Biot本构关系,采用Weibull分布函数表征岩石的非均质性;考虑水化作用引起的损伤,结合当前应力状态的应力损伤得到损伤张量,对弹性模量和渗透率进行损伤分析,实现层理性岩石在水化和荷载作用下的连续损伤演化,形成一套渗流-应力-化学-损伤（HMCD）耦合分析方法。该文给出数值算例,将含损伤横观各向同性模型用于研究层理性岩石的水化特性,表明岩石的非均匀性和水化作用对井壁应力解答具有重要影响,该有限元求解方法可对岩石水化、损伤进行可靠、有效的数值分析。     

二维变分不等式问题的自适应有限元分析

有限元后处理超收敛计算的EEP(单元能量投影)法以及基于该法的自适应有限元分析已在一维变分不等式问题的求解中取得显著成功。以此为基础,该文对二维变分不等式问题成功地实现了自适应有限元分析。该文提出二维区域二分法和二维C 检验技术,有效地提升了松弛迭代的收敛速率,进而应用EEP 超收敛公式计算超收敛解答,用其检验误差并指导网格细分。该文给出的典型数值算例表明该文算法高效、稳定、精确,解答可逐点以最大模度量满足精度要求,堪称为数值精确解。     

一个高效的一维有限元自适应求解的新方案 第十三届全国结构工程学术大会特邀报告

基于新近提出的一维有限元后处理超收敛算法——单元能量投影(EEP)法,将有限元自适应求解问题转化为对超收敛解答的自适应分段多项式插值问题,一步便可获得最优的有限元网格划分,在该网格上再次进行有限元计算,即可获得满足用户给定的误差限的有限元解答。该法简单实用、快速高效,是一个颇具优势和潜力的自适应方法。文中以二阶常微分方程模型问题为例,对该法的形成思路和实施策略做一介绍,并给出有代表性的数值算例用以展示该法的优良性能和效果。     

基于X-SBFEM的裂纹体非网格重剖分耦合模型研究

扩展有限元法利用了非网格重剖分技术,但需要基于裂尖解析解构造复杂的插值基函数,计算精度受网格疏密和插值基函数等因素影响。比例边界有限元法则在求解无限域和裂尖奇异性问题优势明显,两者衔接于有限元法理论内,可建立一种结合二者优势的断裂耦合数值模型。该文从虚功原理出发,利用位移协调与力平衡机制,提出了一种断裂计算的新方法X-SBFEM,达到了扩展有限元模拟裂纹主体、比例边界有限元模拟裂尖的目的。在数值算例中,通过边裂纹和混合型裂纹的应力强度因子计算,并与理论解对比,验证了该方法的准确性和有效性。     

一维C~1有限元后处理超收敛计算的新型算法

该文针对一维C~1有限元提出一种新型后处理超收敛算法,由该法可求得全域超收敛的位移和内力。该法在单个单元上逐单元实施,通过将单元端部结点位移有限元解设为本质边界条件,在单元域上建立单元位移恢复的局部边值问题。对该局部边值问题,以单元内任一点为结点将单元划分为两个子单元进行有限元求解,子单元次数与原单元相同,由此获得该点位移的超收敛解。对单元内所有点均作这样的超收敛求解,可获得整个单元上位移的超收敛解。该位移超收敛解光滑、连续,通过对该位移超收敛解求导可获得转角和内力的超收敛解。数值结果表明,对于m次元,该法得到的挠度和转角具备与结点位移相同的h~(2m-2)阶的最佳收敛阶;弯矩和剪力则分别具备h~(2m-3)、h~(2m-4)阶的收敛阶,均比相应有限元解高出m-2阶。该法可靠、高效、易于实施,是一种颇具潜力的后处理超收敛算法。     

Rapid evaluation of notch stress intensity factors using the peak stress method: Comparison of commercial finite element codes for a range of mesh patterns

The peak stress method (PSM) is an engineering, finite element (FE)‐oriented method to rapidly estimate the notch stress intensity factors by using the singular linear elastic peak stresses calculated from coarse FE analyses. The average element size adopted to generate the mesh pattern can be chosen arbitrarily within a given range. Originally, the PSM has been calibrated under pure mode I and pure mode II loadings by means of Ansys FE software. In the present contribution, a round robin between 10 Italian universities has been carried out to calibrate the PSM with 7 different commercial FE codes. To this aim, several two‐dimensional mode I and mode II problems have been analysed independently by the participants. The obtained results have been used to calibrate the PSM for given stress analysis conditions in (i) FE software, (ii) element type and element formulation, (iii) mesh pattern, and (iv) criteria for stress extrapolation and principal stress analysis at FE nodes.     

Coupled meshfree and fractal finite element method for mixed mode two‐dimensional crack problems

This paper presents a coupling technique for integrating the element‐free Galerkin method (EFGM) with the fractal finite element method (FFEM) for analyzing homogeneous, isotropic, and two‐dimensional linear‐elastic cracked structures subjected to mixed‐mode (modes I and II) loading conditions. FFEM is adopted for discretization of the domain close to the crack tip and EFGM is adopted in the rest of the domain. In the transition region interface elements are employed. The shape functions within interface elements which comprise both the EFG and the finite element (FE) shape functions, satisfies the consistency condition thus ensuring convergence of the proposed coupled EFGM–FFEM. The proposed method combines the best features of EFGM and FFEM, in the sense that no special enriched basis functions or no structured mesh with special FEs are necessary and no post‐processing (employing any path independent integrals) is needed to determine fracture parameters, such as stress‐intensity factors (SIFs) and T‐stress. The numerical results show that SIFs and T‐stress obtained using the proposed method are in excellent agreement with the reference solutions for the structural and crack geometries considered in the present study. Also, a parametric study is carried out to examine the effects of the integration order, the similarity ratio, the number of transformation terms, and the crack length to width ratio on the quality of the numerical solutions. A numerical example on mixed‐mode condition is presented to simulate crack propagation. As in the proposed coupled EFGM–FFEM at each increment during the crack propagation, the FFEM mesh (around the crack tip) is shifted as it is to the new updated position of the crack tip (such that FFEM mesh center coincides with the crack tip) and few meshless nodes are sprinkled in the location where the FFEM mesh was lying previously, crack‐propagation analysis can be dramatically simplified. Copyright © 2010 John Wiley & Sons, Ltd.     

A finite element algorithm to study creep cracks based on the creep hardening surface

This work addresses finite element (FE) modelling of creep cracks under reversed and cyclic loads in steels. A constitutive model based on the creep hardening surface developed by Murakami and Ohno has been selected for this purpose. This model is particularly accurate for describing creep under reversed and cyclic loads and requires no additional material constants. An FE algorithm for this model has been derived and implemented into a research code FVP. The algorithm is verified by comparing the numerical predictions with closed form solutions for simple geometries and loading configurations. FE predictions are compared with experimental data for a stationary crack in a compact tension specimen. The stress and strain fields in the vicinity of a crack under a sustained load are compared with those for the intermediate unloading case. Several integral fracture parameters are investigated as to their appropriateness for describing creep cracks under reversed and cyclic loads.     

Stress intensity factors and crack opening displacements for slanted axial through-wall cracks in pressurized pipes

This paper proposes elastic stress intensity factors and crack opening displacements (CODs) for a slanted axial through-wall cracked cylinder under an internal pressure based on detailed three-dimensional (3D) elastic finite element (FE) analyses. The FE model and analysis procedure were validated against existing solutions for both elastic stress intensity factor and COD of an idealized axial through-wall cracked cylinder. To cover a practical range, four different values of the ratio of the mean radius of cylinder to the thickness ( R _m/ t ) were selected. Furthermore, four different values of the normalized crack length and five different values of the ratio of the crack length at the inner surface to the crack length at the outer surface representing the slant angle were selected. Based on the elastic FE results, the stress intensity factors along the crack front and CODs through the thickness at the centre of the crack were provided. These values were also tabulated for three selected points, that is, the inner and outer surfaces and at the mid-thickness. The present results can be used to evaluate the crack growth rate and leak rate of a slanted axial through-wall crack due to stress corrosion cracking and fatigue. Moreover, the present results can be used to perform a detailed leak-before-break analysis considering more realistic crack shape development.     

Local/global non-intrusive crack propagation simulation using a multigrid X-FEM solver

A Local/global non-intrusive coupling algorithm is proposed for the analysis of mixed-mode crack propagation. It is based on a three scale multigrid and extended finite element method, that was proposed recently for the direct estimation of stress intensity factors of static cracks. The algorithm couples a linear elastic global model (possibly performed by a industrial software) with an enhanced local model capable of modeling a crack and accurately estimating SIFs (performed by a separate research code). It is said non-intrusive since it does not modify the global mesh, its connectivity and solver. For the global model, the contribution of the local patch consists in additional nodal efforts near the crack, which makes it compatible with most softwares. Further the shape of the domain over which the local model is applied is automatically adapted during propagation.     

Fatigue crack growth in residual stress fields

A fatigue crack growth (FCG) model for specimens with well-characterized residual stress fields has been studied using experimental analysis and finite element (FE) modeling. The residual stress field was obtained using four point bending tests performed on 7050-T7451 aluminum alloy rectangular specimens and consecutively modeled using the FE method. The experimentally obtained residual stress fields were characterized using a digital image correlation technique and a slitting method, and a good agreement between the experimental residual stress fields and the stress field in the FE model was obtained. The FE FCG models were developed using a linear elastic model, a linear elastic model with crack closure and an elastic–plastic model with crack closure. The crack growth in the FE FCG model was predicted using Paris–Erdogan data obtained from the residual stress free samples, using the Harter T-method for interpolating between different baseline crack growth curves, and using the effective stress intensity factor range and stress ratio. The elastic–plastic model with crack closure effects provides results close to the experimental data for the FCG with positive applied stress ratios reproducing the FCG deceleration in the compressive zone of the residual stress field. However, in the case of a negative stress ratio all models with crack closure effects strongly underestimate the FCG rates, in which case a linear elastic model provides the best fit with the experimental data. The results demonstrate that the negative part of the stress cycle with a fully closed crack contributes to the driving force for the FCG and thus should be accounted for in the fatigue life estimates.     

二维泊松方程问题Lagrange型有限元p型超收敛算法

该文针对二维泊松方程问题的Lagrange型有限元法提出了一种p型超收敛算法。该法受有限元线法对二维问题降维思想的启发,基于网格结点位移的天然超收敛性,通过从网格中取出一行对边相邻的单元作一子域,将子域内各单元另一对边解答取为原有限元解答,在子域上建立真解近似满足的局部偏微分方程边值问题,对该局部边值问题,沿对边方向单向提高单元阶次进行有限元求解获得单元对边上的超收敛解。单元另一对边上的超收敛解可通过另一方向的单元行类似获得。在单元边超收敛解的基础上,依次取出各个单元,以单元边位移超收敛解为Dirichlet边界条件,双向提高单元阶次对原泊松方程问题进行有限元求解即可获得全域超收敛解。数值算例表明,通过简单的后处理计算本法可显著提高解答的精度和收敛阶。     

Elastic T stress estimates for circumferential surface‐cracked cylinders*

On the basis of detailed three‐dimensional (3D) elastic finite element (FE) analyses, this paper provides tractable approximations for elastic T stress solutions for circumferential inner‐surface cracks in cylinders. Internal pressure and global bending moment were considered. The FE model and analysis procedure employed in the analysis were verified using existing solutions for both elastic stress intensity factor and T stress. To cover a practical range, three different values of the ratio of the mean radius of cylinder to the thickness, R_m/t, were selected; furthermore, four different values of the ratio of the crack depth to the thickness, a/t, ranging from 0.1 to 0.75 and three different values of θ/π ranging from 0.1 to 0.4 were selected. On the basis of FE analyses results, polynomial approximations were proposed at three different locations: surface point, middle point and deepest point. On the basis of the detailed 3D elastic FE analysis, the solutions presented are believed to be the most accurate, and thus provide valuable information for structural integrity assessment considering a crack‐tip constraint.     

Non-local regularization for FE simulation of damage in ductile materials

The paper presents an application of the non-local regularization to the finite element modelling of ductile damage and tearing. In order to model the damage growth in ductile materials of structural components, integral limiters have been introduced. These integral limiters, which are spatial averaging operators, can prevent the problem of mesh-sensitivity of the finite element computations. Hence, it was important to establish the characteristic length l_c for spatial averaging operators of the non-local regularization. More formally, the in-plane distance l_c and the out-of-plane FE dimension have been specified, characterizing the volume over which averaging of stress and strain was carried out to ensure that the continuum theory can represent the physical process of damage.

In order to check the reliability and transferability of the method, FE simulations of various testing examples have been carried out, namely ductile fracture in notched tensile specimens, ductile crack growth in C(T) specimens and application to the pipe-bending test.