基于扩展有限元法的裂尖场精度研究

扩展有限元方法基于单元分解的基本思想,通过引入位移加强函数来表征裂纹的不连续性和裂尖的奇异性。在裂尖加强单元与常规单元之间有一层混合单元,当对裂尖特定区域进行加强时,混合单元个数相应增加,混合单元个数与计算精度存在一定联系。本文提出一种正方形裂尖加强区域的选择方式,可得到较单个加强和圆形加强精度更高、更稳定的计算结果。对于不同长度的裂纹,表征裂尖场奇异性所需的裂尖加强范围存在较大差异,以正方形裂尖加强方式进行计算,得到了不同裂纹长度下最优的加强尺寸。     

A NEW CRITERION OF COMBINED TYPE CRACK—THE CRITERION OF MINIMUM DISTANCE(r_(min))IN PLASTIC REGION

The existence of the plastic area around the crack tip is an important factor against thecracking.In this plastic region.the cracking most likely develops in such a direction alongwhich the distance from the existant crack tip to the plastic region edge is the shortest one.Substituting the stress components at the crack tip for a combined cracking into R.Von.Mises equation of yielding condition,the boundary equation of the plastic regionaround the crack tip can be expressed as:where r is the distance from the crack tip to the boundary of plastic region.The direction inwhich the cracking develops can be determined by following conditions:One form of the material failures is the slipping under the action of shear stress.Underthe conditions of complex stress,the octahedral shear stress on a point around the crack tipis here.B devotes the.factor of the octahedral shear stress,which relates to the size. pattern ofcrack and the exerted stress.Its magnitude reflects the strength of the stress field aroundthe c     

Failure analysis of stainless steel at elevated temperatures

In this paper, particular emphasis has been put on gathering information on the phenomena that take place at the crack tip of a crack propagating at 1100°F. Since the experimental program was directed toward studying crack propagation in tubing, the tests were conducted on rings.From the experimentally obtained data and from the correlation with the theoretically predicted values, the following picture emerges for the fracture behavior with full plasticity present. There is a region surrounding the crack tip where very large plastic deformations take place. This region is surrounded by a much larger region where the loading is nearly proportional and the behavior can be predicted well by the results of the deformation theory of plasticity and the theory of singularity fields. As the crack propagation initiates, there is a drastic change in the crack-tip configuration. The crack tip does not blunt and a fairly sharp crack-tip region is observed. The crack tip carries a large deformation field of a far more localized nature than that observed at the initiation of the crack growth.Paper was presented at 1978 SESA Spring Meeting held in Wichita, KS on May 14–19.     

三维切口/裂纹结构的扩展边界元法分析

在线弹性理论中,三维 V 形切口/裂纹结构尖端区域存在多重应力奇异性,常规数值方法不易求解. 本文提出和建立了三维扩展边界元法 (XBEM),用于分析三维线弹性 V 形切口/裂纹结构完整的位移和应力场. 先将三维线弹性 V 形切口/裂纹结构分为尖端小扇形柱和挖去小扇形柱后的外围结构. 尖端小扇形柱内的位移函数采用自尖端径向距离 $r$ 的渐近级数展开式表达,其中尖端区域的应力奇异指数、位移和应力特征角函数通过插值矩阵法获得. 而级数展开式各项的幅值系数作为基本未知量. 挖去扇形域后的外围结构采用常规边界元法分析. 两者方程联立求解可获得三维 V 形切口/裂纹结构完整的位移和应力场,包括切口/裂纹尖端区域精细的应力场. 扩展边界元法具有半解析法特征,适用于一般三维 V 形切口/裂纹结构完整位移场和应力场的分析,其解可精细描述从尖端区域到整体结构区域的完整应力场. 作者研制了三维扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明了扩展边界元法求解三维 V 形切口/裂纹结构完整应力场的准确性和有效性.      

ANALYSIS OF 3-D NOTCHED/CRACKED STRUCTURES BY USING EXTENDED BOUNDARY ELEMENT METHOD 1)

在线弹性理论中,三维 V 形切口/裂纹结构尖端区域存在多重应力奇异性,常规数值方法不易求解. 本文提出和建立了三维扩展边界元法 (XBEM),用于分析三维线弹性 V 形切口/裂纹结构完整的位移和应力场. 先将三维线弹性 V 形切口/裂纹结构分为尖端小扇形柱和挖去小扇形柱后的外围结构. 尖端小扇形柱内的位移函数采用自尖端径向距离 $r$ 的渐近级数展开式表达,其中尖端区域的应力奇异指数、位移和应力特征角函数通过插值矩阵法获得. 而级数展开式各项的幅值系数作为基本未知量. 挖去扇形域后的外围结构采用常规边界元法分析. 两者方程联立求解可获得三维 V 形切口/裂纹结构完整的位移和应力场,包括切口/裂纹尖端区域精细的应力场. 扩展边界元法具有半解析法特征,适用于一般三维 V 形切口/裂纹结构完整位移场和应力场的分析,其解可精细描述从尖端区域到整体结构区域的完整应力场. 作者研制了三维扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明了扩展边界元法求解三维 V 形切口/裂纹结构完整应力场的准确性和有效性.     

非晶态高聚物裂尖银纹损伤场

在非晶态高聚物裂端塑性区引入银纹损伤，提出了一个修正的Dugdale模型，并由它确定塑性区的大小，算例表明，银纹损伤导致了塑性区深度的增加。此外，在Dugdale模型下研究了稳态扩展的端塑性区，在裂纹缓慢扩展、小范围银纹化以及材料为线粘弹性等假定下，由能量原理得到了控制裂纹生长的一个非线性微分一积分方程，从中看到银纹损伤的影响十分显著。     

Photoelastic study of the stresses near a cleavage crack

Photoelastic analysis has been used to determine the stress distribution near the tip of a cleavage crack. The crack bisects a slender rectangular bar along its length, and forces are applied at one end of it tending to open it out. The region of high stresses is localized in a region of the order of the half-width of the bar in size. Very near the crack tip, the stresses decrease in proportion to the inverse square root of the distance from the crack tip. The maximum gradient of the principal tensile stress lies at an angle of 70 to 80 deg from the plane of the crack. Contour maps of the principal stresses, the 45-deg shear stress and the maximum tensile stress are presented.     

Finite deformation analysis of COD,J-integral and crack tip intense strain region in plane strain large-scale yielding

A finite element analysis was performed to simulate crack tip blunting and the development of the intense strain region in a small compact tension specimen (0.4 T CT) of SA533B-1 under plane strain large-scale yielding, with the condition of large-geometry change around the crack tip taken into consideration. The region where the equivalent plastic strain $\text{}\text{?}$g3_p is greater than 0.15 was defined as the intense strain region, which corresponded to the recrystallized-etched zone delineated experimentally around the blunting crack tip. The development of the intense strain region was discussed as a function of the J-integral and the crack opening displacement. A linear relationship was obtained between the plastic work W_p dissipated within the intense strain region and (J/σ_y)² or b², where b is the crack opening displacement, defined as the separation of the two points at which the boundary of the intense strain region surrounding the crack tip intersects with the free surfaces of the crack.     

Plastic flow at a crack tip. Energy failure criterion and its relationship to the J integral

Plastic flow at the tip of a crack moving in an elastic body is considered using the theory of an ideally rigid-plastic body. The material at the crack tip is treated as a body consisting of an elastic outer region and a rigid-plastic inner region. It is shown that this representation is energetically justified for small plastic regions. The distribution of the specific dissipation of the work of internal forces and deformations along the particle trajectory at the crack tip is obtained. The relationship between the specific dissipation of the work of internal forces and the J integral under plane-strain conditions is established.     

Multiscale modeling of crack initiation and propagation at the nanoscale

Fracture occurs on multiple interacting length scales; atoms separate on the atomic scale while plasticity develops on the microscale. A dynamic multiscale approach (CADD: coupled atomistics and discrete dislocations) is employed to investigate an edge-cracked specimen of single-crystal nickel, Ni, (brittle failure) and aluminum, Al, (ductile failure) subjected to mode-I loading. The dynamic model couples continuum finite elements to a fully atomistic region, with key advantages such as the ability to accommodate discrete dislocations in the continuum region and an algorithm for automatically detecting dislocations as they move from the atomistic region to the continuum region and then correctly “converting” the atomistic dislocations into discrete dislocations, or vice-versa. An ad hoc computational technique is also applied to dissipate localized waves formed during crack advance in the atomistic zone, whereby an embedded damping zone at the atomistic/continuum interface effectively eliminates the spurious reflection of high-frequency phonons, while allowing low-frequency phonons to pass into the continuum region.The simulations accurately capture the essential physics of the crack propagation in a Ni specimen at different temperatures, including the formation of nano-voids and the sudden acceleration of the crack tip to a velocity close to the material Rayleigh wave speed. The nanoscale brittle fracture happens through the crack growth in the form of nano-void nucleation, growth and coalescence ahead of the crack tip, and as such resembles fracture at the microscale. When the crack tip behaves in a ductile manner, the crack does not advance rapidly after the pre-opening process but is blunted by dislocation generation from its tip. The effect of temperature on crack speed is found to be perceptible in both ductile and brittle specimens.     

Stress and energy density fields near a blunted crack front in power hardening material

Using Jaumann and Dienes rates of Euler stress in elastic-plastic constitutive equations of finite deformation, plane strain finite element analysis for a compact tension specimen with a blunted crack front is made. The Euler stress, Kirchhoff stress and volume strain energy density near a blunted crack tip are computed. Constitutive relations with different deformation rates affect the the near crack tip solution in a region within an order of magnitude of the crack opening displacement. The results differed from the corresponding solution of deformation plasticity (or nonlinear elasticity) with increasing deformation. They are smaller in a local region of about 2 to 10 times of the crack opening distance.The volume energy density near the crack tip is computed, the stationary values of which determine the locations of extensive yielding and possible sites of crack initiation. It remained nearly constant with increasing deformation. Such a character tends to support the volume energy density criterion as a means for quantifying the ductile fracture behavior of metals.     

NiTi形状记忆合金裂纹尖端热力耦合行为研究

本文对NiTi形状记忆合金I型裂纹尖端热力耦合行为进行了数值仿真分析和实验验证。建立了包含相变和热力耦合的本构模型,通过有限元计算得到了裂纹尖端附近的纵向应变、马氏体体积分数和温度场分布,依据马氏体相变情况对裂纹尖端有效应力强度因子进行了修正,揭示了加载速率对形状记忆合金裂纹尖端有效应力强度影子的影响规律。参数研究表明,随着加载频率的增加,裂纹尖端附近温度逐渐升高,马氏体相变区域逐渐缩小,有效应力强度因子呈下降趋势,形状记忆合金表现出增韧效应,有助于减缓裂纹扩展。本研究结果对于揭示热力耦合作用下超弹性形状记忆合金疲劳裂纹扩展规律具有重要参考意义。     

Nonlinear response of anti-plane shear crack

Investigated in this communication is the nonlinear response of an anti-plane shear crack. A nonlinear hardening/softening relation is assumed between the resultant shear stress and shear strain for determining the state of affairs near the crack tip. Amplification of the local stress fields is obtained in terms of the material parameters by assuming that nonlinearity is confined to the immediate vicinity of the crack tip region.     

扩展边界元法分析切口和裂纹结构应力场的准确性

在线弹性理论中,切口/裂纹结构尖端区域存在奇异应力场,数值方法不易求解。本文建立的扩展边界元法(XBEM)对围绕尖端区域位移函数采用自尖端径向距离 的渐近级数展开式表达,其级数项的幅值系数作为基本未知量,而外部区域采用常规边界元法离散方程。两者方程联立求解可获得切口和裂纹结构完整的位移和应力场。扩展边界元法具有半解析法特征,适用于一般的切口和裂纹结构应力场分析,其解可精细描述从尖端区域到整体结构区域的应力场。作者研制了扩展边界元法程序,文中给出了两个算例,通过计算结果分析,表明扩展边界元法求解切口和裂纹结构应力场的准确性和有效性。     

Logarithm strain singularity at tip of mode I growing crack in elasticand perfectly-plastic material

Reanalyzed in detail is the stress and strain distribution near the tip of a Mode I steadily growing crack in an elastic and perfectly-plastic material. The crack tip region is divided into five angular sectors, one of which is singular in character and represents a rapid transition zone that becomes a line of strain discontinuity in the limit as crack tip is approached. It is shown for an incompressible material (ν=0.5) under plane strain that the local strain in all the angular sectors possesses the same logarithm singularity, i.e., In r where r is the radial distance measured from the crack tip. This result also prevails for the compressible material ( v < 0.5) and resolves a long standing controversy concerning the strain singularity in the sector just ahead of the crack tip.     

热辐射应力图像分析确定裂纹尖端塑性区

本文借助于热辐应力图像分析技术记录了铝质紧养伤拉伸试件的热辐射应力图像。根据裂尖周围温度变化曲线的分布规律确定其特征区的,从而得以裂尖塑性区的开头曲线,并与理论值比较,符合良好。对裂纹扩展面进行ＳＥＭ分析,得出该试件的断裂特性,对于进行细观损伤研究有重要的意义。     

Crack tip asymptotic character of anti-plane stress and strain rate for linear fractional constitutive relations

Linear-fractional strain rate and stress relations are used to simulate materials undergoing steady state creep. The crack tip asymptotic character of the stress and strain rate field is obtained in exact and approximate form. In the limit as the radial distance emanating from the crack tip approaches zero, the stress field corresponds to that for an ideal plastic material while the exact and approximate solutions tend to coincide. Discussed is the nonhomogeneous singular character of the strain rate field that possess different orders of singularities in a circular region around the crack tip.     

复合型裂纹断裂的新准则

以Ⅰ-Ⅱ复合型裂纹为研究对象,对裂纹尖端的塑性区分布规律进行了理论分析。引入两组评价裂纹尖端应力场对裂纹扩展影响的参数,考虑裂纹尖端存在的局部塑性变形,并采纳如下两个假设,(1)裂纹沿最短路径穿过塑性区向弹性区扩展,(2)当在扩展方向上的弹塑性边界极半径r大于其临界极半径rC时,裂纹开始扩展。在此基础上,导出了新的复合型裂纹断裂准则,并与现有部分断裂准则及实验结果进行了对比。结果表明：新建立的复合裂纹断裂准则与实验结果吻合程度非常高。最后,阐明了以单一KIC或KIIC建立的复合裂纹扩展准则的局限性以及考虑裂纹尖端当前应力场的必要性。     

A new criterion for mixed mode fracture initiation based on the crack tip plastic core region

In this work, we propose a new criterion for mixed mode I-II crack initiation angles based on the characteristics of the plastic core region surrounding the crack tip. The shape and size of the plastic core region are thoroughly analyzed under different loading conditions and a new formulation for the non-dimensional variable radius of the core region is presented for mixed mode (K_I−K_II) fracture. The proposed criterion states that the crack extends in the direction of the local or global minimum of the plastic core region boundary depending on the resultant stress state at the crack tip. The results show a well-defined correlation between the plastic core region characteristics and crack extension angles predicted by other criteria. The proposed criterion is formulated for various loading conditions and is compared with other available criteria against the limited available experimental data. It is shown that the proposed criterion provides a better agreement with the experimental data.