表面裂纹疲劳扩展的数值模拟(Ⅱ)

根据Paris疲劳裂纹扩展规律,对拉伸和纯弯曲疲劳载荷下表面裂纹扩展进行了数值模拟。数值模型中,用三次样条函数曲线拟合裂纹尖端,在裂纹扩展增量计算中考虑了裂纹闭合影响。裂纹形状演化的模拟结果与Newman和Raju经验公式预测结果进行了比较,表明了所采用的数值模拟方法的实用性。研究发现,裂纹闭合对疲劳裂纹扩展过程中的裂纹形状演化以及裂纹尖端的应力强度因子(SIF)分布都有明显影响。同裂纹形状演化一样,疲劳裂纹扩展过程中裂纹尖端的SIF分布表现出明显的特征。最后,建议了一个简单函数来统一描述表面裂纹尖端的SIF分布。     

巴西圆盘试样中心斜裂纹疲劳扩展轨迹的边界元模拟

提出了一个平面弹性体多裂纹疲劳扩展模型. 它主要涉及到复合型加载情况下多裂纹尖端疲劳扩展的数学模型及杂交位移不连续法(一种边界元法). 在数值模拟中, 对每一裂纹扩展增量分析时,在其先前的边界上增添裂纹扩展增量, 且只对新增添的裂纹扩展增量划分单元, 同时, 按照这种边界元法的实施方法对一些单元特征进行调整, 就可以方便地模拟裂纹扩展. 用这种数值方法模拟了巴西圆盘试样中心斜裂纹疲劳扩展轨迹,数值结果说明了预报模型的有效性, 揭示了裂纹体几何对疲劳扩展的影响.      

基于FE/EFG耦合方法的沥青路面裂纹扩展模拟

沥青路面裂纹在行车荷载下的裂纹扩展规律对于沥青路面的设计、维护具有重要的指导意义.通过预切口沥青混合料小梁疲劳试验,利用无网格伽辽金/有限元耦合方法,对沥青路面表面裂纹扩展进行了数值模拟,研究了面层厚度、面层模量、基层厚度、基层模量对裂纹扩展的影响规律.结果显示:随着预切口距离加载中心距离的增大,沥青混合料小梁的疲劳寿命增大,裂纹尖端的Ⅰ型应力强度因子减小,裂纹所表现出的Ⅱ型特性增 强;在沥青路面表面裂纹扩展过程中,应力强度因子经历一个急剧上升,达到一个极值后缓慢下降,然后又缓慢上升,最后急剧增加的过程;随面层、基层厚度的增加,表面裂纹尖端的应力强度因子降低;面层、基层模量对裂纹扩展路径的影响较大.     

考虑T-stress的非正态随机参数疲劳裂纹扩展寿命的可靠性分析

基于断裂力学的疲劳裂纹扩展寿命问题的研究常常将裂纹尖端应力展开项的高次项忽略,引起了裂纹扩展模拟的误差,本文考虑高次项T-stress对裂纹扩展角的影响,对裂纹扩展过程做了数值模拟,结果显示相同裂纹扩展长度下,考虑T-stress会延长裂纹扩展寿命。文章首先采用修正的Paris-Erdogan 公式计算了两端承受均布拉伸载荷的边缘斜裂纹板的疲劳裂纹扩展寿命,裂纹扩展方向采用两参数修正的最大拉应力准则。由于结构尺寸,材料特性和载荷等因素具有不确定性,导致疲劳裂纹扩展过程带有一定的随机性,本文以材料属性和载荷为随机变量,在随机有限元法的基础上,结合计算可靠度的四阶矩法,Edgeworth级数展开方法,提出随机参数服从任意分布时的结构疲劳裂纹扩展寿命可靠度的计算方法。分析了参数为非正态分布时的平板裂纹扩展寿命可靠度随裂纹扩展的变化过程。本文方法可预测工程中板裂纹的扩展行为,以及预测裂纹板的可靠度。     

粘聚裂纹扩展的强化有限元h型网格自适应模拟

非连续变形分析和非规则节点处理是基于单元细划的粘聚裂纹扩展网格自适应模拟的关键。首先,利用强化有限单元法中数学单元和物理单元分离的特点,通过引入过渡单元,将适用于非连续变形描述的数学模式覆盖法和方便处理非规则节点的物理模式重构法结合,提出了强化有限单元法的统一关联法则,并导出了相应的单元列式。其次,基于数学裂纹尖端影响域和裂尖单元尺寸,提出了基于强化有限单元法的粘聚裂纹扩展过程模拟的h型网格自适应策略。最后,通过两个算例验证了本文方法的合理性和有效性。     

数学网格和物理网格分离的有限单元法(II):粘聚裂纹扩展问题中的应用

强化有限单元法将物理网格与数学网格分离开来,可以方便地描述非连续变形;粘聚区域模型是模拟断裂过程区作用最简单有效的方法,且可以避免裂纹尖端的应力奇异性.本文以平面问题为例,将强化有限单元法与粘聚区域模型相结合,利用富集数学节点描述任意粘聚裂纹扩展过程中的非连续变形问题,提出了裂纹扩展过程中数学节点富集和数学单元定义的方法.本文还导出了与平面4～8节点平面等参单元对应的8～16节点粘聚裂纹单元列武.最后,通过三点弯梁的裂纹扩展过程模拟验证了本文提出的粘聚裂纹扩展模拟方法的有效性.     

基于光滑边域有限元法的二维裂纹扩展分析

为了提高预测裂纹扩展路径的准确性和效率,本文将光滑边域有限元法和虚拟裂纹闭合法相结合,利用常应变三角形单元,获得裂纹尖端处的断裂控制参量应力强度因子,并运用最大拉应力准则求得裂纹在荷载作用下的启裂方向,对裂纹扩展轨迹给出自动跟踪方法;对三个典型二维裂纹扩展模型,预测了裂纹扩展路径,并将结果与参考文献中的结果进行对比,验证方法的有效性与准确性。数值结果表明:该方法具有单元简单、网格尺寸要求低、裂纹尖端处单元不需特殊处理等优点,是分析裂纹扩展问题简洁高效的数值计算方法。     

三维裂纹扩展轨迹的边界元数值模拟

提出了一种对三维裂纹扩展轨迹进行数值模拟的新方法。采用一种新的具有C^1连续性、高精度的单节点二次边界单元，使边界元(BEM)的分析效率和裂纹张开位移(COD)、应力强度因子(SIF)的精度大大提高。采用裂纹张开位移全场拟合法(GCDFP)求出裂纹面前缘的SIF，所得到的SIF达到与所用的COD资料同样的精度。使用Paris公式求出裂纹前缘各点的裂纹扩展增量，并用三次B样条函数对这些增量进行拟合，得到新的光滑裂纹前缘。根据以上思想方法，开发了具有较高的计算效率和精度的数值模拟软件。此软件可以自动跟踪裂纹扩展，得到裂纹扩展的轨迹。运用该软件对椭圆和矩形裂纹的扩展轨迹进行了数值模拟。其结果与理论上的预言完全一致，裂纹最后都趋于一个圆裂纹，具有实际指导意义。     

扭力轴三维裂纹扩展寿命仿真研究

对疲劳载荷谱作用下三维表面裂纹,采用双重边界元理论求解裂纹前沿的应力应变场,运用Forman理论、最小应变能密度法和Elber模型,计算裂纹前沿各点的扩展长度、扩展方向和应力强度因子等特征量.根据增量步下裂纹几何形状的改变,对裂纹面进行网格重划分和迭代计算,模拟三维裂纹的扩展,预测裂纹扩展寿命.扭力轴表面裂纹扩展的仿真结果表明该方法合理可行.     

单向轴压条件下内置椭圆三维裂纹扩展无网格方法的研究

根据三维裂纹前沿局部坐标系中主法平面内的周向应力计算公式,可以确定主法平面内的最大周向应力σθ,并由最大周向拉应力裂纹破裂准则确定裂隙岩体的破裂载荷和裂纹的破裂方向以及破裂步长。在三维无网格方法的节点的影响域中,选择合适的节点密度,以确保计算精度,并对单向轴压条件下内置三维裂纹扩展进行了数值模拟研究。根据数值计算结果,分析了破裂载荷、裂纹扩展角和扩展步长,揭示了内置三维裂纹的破裂过程,结合插值法清楚地展示了三维裂纹的破裂曲面形态。研究结果表明,针对三维裂纹破裂过程和破裂曲面形态无网格方法数值计算结果与FLAC-3D计算结果和试验结果吻合较好。     

海洋平台K型管节点的疲劳裂纹扩展分析 II:数值分析

对于已含初始裂纹平台管节点的寿命预测很大程度上依靠应力强度因子的精确值,而复杂载荷条件下的节点应力强度因子的计算尚无参数方程直接确定。本文提出了一种含表面裂纹的K节点的有限元网格产生方法,即把整个K节点划分为几个子区域,每个子区域的网格具有不同类型的单元和不同的密度。这种方法在控制网格密度,尤其是控制沿着裂纹边缘单元的边长比方面有其独特的优越性,当所有子区域的网格自动产生后,容易得到整个结构的有限元模型。同时用J积分和位移外推插值法分别计算了一个K型节点沿着裂纹前缘的应力强度因子值,发现:试验得到的应力强度因子值和提出的模型计算结果非常吻合,证明了所提有限元模型的准确性。     

Experimental and numerical analysis of fretting crack formation based on 3D X-FEM frictional contact fatigue crack model

Nowadays, numerical simulation of 3D fatigue crack growth is easily handled using the eXtended Finite Element Method coupled with level set techniques. The finite element mesh does not need to conform to the crack geometry. Most difficulties associated to complex mesh generation around the crack and the re-meshing steps during the possible propagation are hence avoided. A 3D two-scale frictional contact fatigue crack model developed within the X-FEM framework is presented in this article. It allows the use of a refined discretization of the crack interface independent from the underlying finite element mesh and adapted to the frictional contact crack scale. A stabilized three-field weak formulation is also proposed to avoid possible oscillations in the local solution linked to the LBB condition when tangential slip is occurring. Two basic three-dimensional numerical examples are presented. They aim at illustrating the capacities and the high level of accuracy of the proposed X-FEM model. Stress intensity factors are computed along the crack front. Finally an experimental 3D ball/plate fretting fatigue test with running conditions inducing crack nucleation and propagation is modeled. 3D crack shapes defined from actual experimental ones and fretting loading cycle are considered. This latter numerical simulation demonstrates the model ability to deal with challenging actual complex problems and the possibility to achieve tribological fatigue prediction at a design stage based on the fatigue crack modeling.     

A domain-independent integral for computation of stress intensity factors along three-dimensional crack fronts and edges by BEM

The present work deals with an evaluation of stress intensity factors (SIFs) along straight crack fronts and edges in three-dimensional isotropic elastic solids. A new numerical approach is developed for extraction, from a solution obtained by the boundary element method (BEM), of those SIFs, which are relevant for a failure assessment of mechanical components. In particular, the generalized SIFs associated to eigensolutions characterized by unbounded stresses at a neighbourhood of the crack front or a reentrant edge and also that associated to T-stress at the crack front can be extracted. The method introduced is based on a conservation integral, called H-integral, which leads to a new domain-independent integral represented by a scalar product of the SIF times some element shape function defined along the crack front or edge. For sufficiently small element lengths these weighted averages of SIFs give reasonable pointwise estimation of the SIFs. A proof of the domain integral independency, based on the bi-orthogonality of the classical two-dimensional eigensolutions associated to a corner problem, is presented. Numerical solutions of two three-dimensional problems, a crack problem and a reentrant edge problem, are presented, the accuracy and convergence of the new approach for SIF extraction being analysed.     

Fatigue crack shapes development of plate-to-plate welded joints

This paper presents crack shapes development of 3D surface semielliptical cracks which are always found at the weld toe of non-load-carrying cruciform welded joints. The alternating current potential drop technique has been used to measure the eight crack depths along the weld toe where probes are placed at 10-mm intervals. 3D crack shapes can be obtained at any particular cycle during the fatigue test. Subregion boundary element methods incorporating the transition and quarter-point elements along the crack front are used to validate the experimental stress intensity factors along the crack front. A method for evaluating the effective stress intensity factors is proposed. It is found that the numerical results generally agree with the experimental results.     

Evaluation of finite-element calculations in a part-circular crack by coherent optics techniques

Two techniques, speckle photography and holographic interferometry, were used to test three-dimensional finite-element calculations in an internally pressurized cylinder with an external part-circular crack. Opening displacements along the crack line were measured by speckle photography. Radial displacements were obtained from holographic fringe patterns. Good agreement between experimental and numerical data is obtained. Stress-intensity factor variations along the crack front are calculated from numerical results.     

裂纹扩展独立于网格的一种有限元单元子划分方法

提出了一种有限元模拟裂纹扩展的单元子划分结合子结构的方法。本方法中,裂纹可以进入或穿过一个单元,或沿单元的边界扩展,因此裂纹可以沿任意路径扩展而不受初始网格的限制。对上述几类包含裂纹的单元按照裂纹的路径进行子划分,覆盖一条裂纹的所有子划分单元就组成了一个子结构,子结构规模随裂纹的扩展而增大。子结构中因单元子划分而新增的结点自由度,通过自由度的凝聚用初始网格结点的自由度表示,因此结构整体分析的总自由度不变。以上述方法为基础建立了裂纹萌生和扩展的准则。用本文的方法分析了单（双）材料无限大平面中心（界面）裂纹的裂尖场,验证了本文方法的精度,并模拟了颗粒复合材料中微裂纹在颗粒、基体和界面中逐步扩展的过程,考核了本文方法对复杂裂纹扩展问题模拟的适用性。     

Simulation of crack extension in shell structures and prediction of residual strength

Two models for the numerical simulation of ductile crack extension in shell structures are presented and compared. They are based on the crack tip opening angle and a cohesive zone approach. After identification of the model parameters and investigations of the mesh dependence, the models are applied to various specimen configurations and structural components. Their excellent numerical performance favors their application for predicting the residual strength of lightweight components like aircraft fuselages.     

基于网格生成的随机凹凸型混凝土骨料细观建模方法

基于目前混凝土细观数值分析多采用二维模型,提出了一种基于网格预先生成的二维随机凹凸型骨料的混凝土细观建模方法。该方法利用已生成的网格状态进行骨料投放,有效解决了在细观三相介质中的网格质量问题。采用随机凹凸型骨料使骨料模型更加接近真实形态,避免了该方法对骨料形貌描述不准确的问题。该方法可保证投放的随机性以及级配分布的合理性,二维投放骨料含量可达75%。最后,采用混凝土损伤塑性模型模拟了混凝土保护层锈胀开裂的细观断裂破坏过程,分析了混凝土裂纹萌生、扩展过程及破坏形态,数值结果与实验结果吻合良好。     

Energy domain integral applied to solve center and double-edge crack problems in three-dimensions

A three-dimensional Boundary Element Method (BEM) implementation of the energy domain integral for the numerical computation of the crack energy release rate is presented in this paper. The domain expression of the energy release rate is naturally compatible with the BEM, since stresses, strains and derivatives of displacements at internal points can be evaluated using the appropriate boundary integral equations. The pointwise crack energy release rate is evaluated along the three-dimensional crack front over a cylindrical domains that surround a segment of the crack front. The accuracy of the implementation is demonstrated by solving several problems, which include geometries containing straight as well as curved crack fronts.