有限厚度板穿透裂纹前缘附近三维弹性应力场分析

通过三维有限元计算来研究有限宽度、有限厚度含有穿透裂纹板的裂纹前缘应力场，从中找出应力强度因子与板的厚度、裂纹长度之间的关系，同时还分析了裂尖的三维约束程度和三维约束区的大小。分析结果表明：应力强度因子沿厚度的分布是不均匀的，应力强度因子的最大值及其位置与厚度有关；有限厚度板中面应力强度因子(KI)m-p及最大应力强度因子(KI)max均大于平面应力或平面应变的应力强度因子。对有限厚度裂纹问题，按平面应力或平面应变来考虑是不安全的；板中面的应力强度因子(KI)m-p及最大应力强度因子(KI)max是厚度B／a的函数；板的中面离面约束系数Tx最大，自由面(z=B)Tx=0。沿厚度方向裂尖附近的离面约束系数Tx也是z／B和B／a的函数，随着厚度的增加离面约束系数Tx增大，离中面越近离面约束系数Tx越大。Tx随着x的增大急剧减小，三维约束影响区域大小大约为板厚的一半，且裂纹长度a／W对应力强度因子沿厚度变化规律及Tx影响区域大小影响较小。     

有限厚双边U型切口板的三维弹塑性有限元分析

针对单向拉伸载荷作用下的有限厚双边 U型切口板 ,本文对其切口根部的三维效应进行了详细的三维弹塑性大变形有限元分析 ,得到了不同于平面切口问题的新结果 :( 1 )在线弹性阶段 ,三维切口根部的应力集中因子 Kt高于二维情形 ,且 Kt在板厚 ( B)约为 8倍切口根部半径 (ρ)时达到峰值 ;进入塑性后 ,切口根部的 Kt随载荷水平的增加逐渐下降 ,并介于平面应力和平面应变情形之间 ;( 2 )切口前缘等效应变εe的最大值出现在约束最高的中面 ,且可比约束最低的自由表面处的应变值高一倍 ;( 3)薄板中 ,塑性区尺寸 xp 的最大值不是出现在约束较低的自由表面 ,而是出现在约束较高的中面 ,但随 B/ρ的增大逐渐向自由表面外推 ;( 4 )面内应力比 Tx在切口根部塑性区内的分布对板厚和载荷水平不敏感 ,可以用 Hill的滑移线解很好地近似。另外 ,对三维影响区内一些重要的三维约束参量进行了详细分析 ,总结了它们的三维分布规律。     

紧凑拉剪试样中Ⅰ/Ⅱ复合型穿透裂纹的三维特性分析

进行了Ⅰ/Ⅱ复合型载荷作用下航空结构铝合金CTS试样线弹性的全场三维有限元计算,分析了复合型离面应力约束因子Tz和面内约束T应力的分布特性,研究了厚度和载荷条件对应力各分量及应力三轴性水平Rσ、应变能密度U0的影响以及这些量在实验中观察到的裂纹起裂方向上的特点。结果表明：（1）平面内约束T应力在此种试样形式下为零;（2）复合型裂纹的三维效应区与厚度成正比,为0.4-0.5B,与载荷条件基本无关,但是厚度效应随着载荷角的减小逐渐变小,到II型载荷时基本消失;（3）离面约束因子Tz随着径向和厚度尺寸的增加逐渐减小,但在周向基本不发生变化;（4）最小应变能密度U0min的方位能够表征此种材料三维复合型断裂的裂纹起裂方向。研究结果为建立三维复合型断裂准则提供了理论基础。     

紧凑拉剪试样中Ⅰ/Ⅱ复合型穿透裂纹的三维特性分析

进行了Ⅰ/Ⅱ复合型载荷作用下航空结构铝合金CTS试样线弹性的全场三维有限元计算,分析了复合型离面应力约束因子Tz和面内约束T应力的分布特性,研究了厚度和载荷条件对应力各分量及应力三轴性水平Rσ、应变能密度U0的影响以及这些量在实验中观察到的裂纹起裂方向上的特点。结果表明：（1）平面内约束T应力在此种试样形式下为零;（2...     

三维约束下幂硬化双材料界面裂尖场分析

本文用Ｊ２形变理论分析了三轴应力约束下，由具有不同硬化指数和硬化系数的幂硬化材料构成的界面裂纹问题．通过在本构方程中引入三轴应力约束因子Ｔｚ，我们考虑了三维约束效应对裂尖场的影响．计算结果表明，三轴应力约束对裂尖场的分布，尤其是对径向应力σｒ和有效应力σｅ有重要影响．本文还探讨了三轴约束和材料性能对裂尖应力奇异性的影响     

V型缺口裂端的三维应力状态及约束分析

利用三维有限元方法,研究了有限厚度板中V型缺口根部穿透裂纹前沿的三维弹性应力场。对不同厚度、不同缺口张开角和裂纹长度对应力强度因子及裂尖附近的三维约束程度的影响进行了分析,同时还讨论了三维约束区的大小。研究结果显示：当缺口张开角小于60度时,不同缺口的应力强度因子和离面约束因子的分布基本一致,角度的影响不明显。应力强度因子是厚度的函数,板中面的应力强度因子随厚度的增加逐渐减小趋近干平面值,最大为1．08倍的平面值。当板厚超过15倍的缺口深度时,应力强度因子最大值将从中面转移至接近自由表面位置,距中面约0．4倍板厚。三维约束非常明显的区域在裂尖前沿约0.45倍厚度的范围内．二维到三维的过渡区在裂尖前沿1.5倍厚度的区域内;在中面上三维效应影响区最大,随着离中面距离的增加逐渐减小,在自由表面上降为0。     

过载作用下聚乙烯疲劳裂纹扩展行为研究

论文针对中密度聚乙烯材料(MDPE)，采用平板试样进行了I型疲劳裂纹扩展和单次过载下裂纹扩展试验．发现与金属材料类似，单次拉伸过载对聚乙烯(PE)的疲劳裂纹扩展有明显的迟滞作用，降低了裂纹扩展速率．试验还通过变载荷刻线法获取疲劳裂纹扩展前缘的实际形貌和变化规律，对常规变载荷刻线方法进行了调整和验证，其修正方法对高分子材料的疲劳裂纹扩展前缘刻线具有较好的效果．通过观察发现含楔形塑性区的裂尖钝化是裂纹迟滞的主要原因．过载引入的塑性区内残余应力对裂纹迟滞也起了重要作用．论文利用Dugdale模型计算了塑性区尺寸，使用基于残余应力的Wheeler模型对过载迟滞进行了很好的拟合．     

I型裂纹尖端约束应力区模型及其解析解

考虑了I型裂纹尖端损伤区域内三种不同的约束应力分布形式,即右三角分布形式(情况A)、均匀分布形式(情况B)、左三角分布形式(情况C),并采用复变函数方法求得了应力强度因子与裂纹张开位移的解析解;在此基础上,通过数值计算得到了应力强度因子和裂纹张开位移随约束应力区长度、约束应力大小以及分布形式的变化规律。研究结果表明:随裂尖材料损伤程度的增加,裂尖损伤区内约束应力减小,应力强度因子和裂纹张开位移增大;约束应力的分布形式对应力强度因子和裂纹张开位移有显著影响;相对于其他区域,约束应力对裂纹尖端区域裂纹张开位移的影响较大。然而,对于裂尖损伤区域的形成与作用荷载、材料性质、构件几何尺寸之间的关系,还需要进行更为深入的研究。     

裂尖超钝化:实验、理论与数值模拟

本文报道了在改性聚丙烯裂尖平面应变段发生超钝化的实验发现。利用环向冷拉模型,可说明裂尖超钝化的形成取决于冷拉传播速率与加载率之比。裂尖超钝化可使材料韧性,尤其是冲击韧性的改善,达到一个新的水平。本文在超弹性-粘塑性各向异性损伤的有限变形本构模型下进行了详尽的数值计算。模拟结果成功地再现了裂尖超钝化的基本特征。     

金属韧性断裂的细观研究

对三点弯曲试件裂尖及断裂过程分别用 Prandtl-Reuss本构、Gurson本构方程作了大变形弹塑性有限元分析 ,并比较了二者的结果。对裂尖应力分布 ,应变分布和裂尖处空穴演化作了初步研究。计算模拟结果表明裂尖空穴化区域是一个很小量级的区域 ,采用 Gurson本构方程来研究裂尖“过程区”比采用 Prandtl-Reuss本构方程分析“过程区”来的合理 ,更接近材料实际断裂过程。     

CONSTRAINTS ON THE CRACK TIP PLASTIC FIELDS UNDER LARGE SCALE YIELD

Plane strain elastic-plastic finite element analyses are used to study the stress, strain fields around a straight crack in power hardening plastic material. Center crack panel (CCP), single edge crack panel (SECP) and double edge crack panel (DECP) tension specimens are analyzed with various crack lengths. Two local constraint parameters, i.e. in-plane stress ratioT _x and out-of-plane constraintT are analyzed, which are defined as tangential stress dividing normal (open) stress and out-of-plane stress dividing the sum of tangential stress and normal stress respectively. Numerical analyses indicate that the two local constraint parameters are nearly independent of the specimen geometry and the loading level in the plastic zone. Methods of estimating two constraint parameters and ways to define the boundary of the blunt zone are given. By using these parameters, the stress distributions in the plastic zone, especially in the blunt zone can be accurately predicted. Project supported by the National Foundation of Distinguished Young Scientists of China (No. 59625510).     

In-plane and out-of-plane constraint effects under pressurized thermal shocks

Transferability of fracture toughness data obtained on small scale specimens to a full-scale cracked structure is one of the key issues in integrity assessment of engineering structures. In order to transfer fracture toughness under different constraints, both in-plane and out-of-plane constraint effect should be considered for the specimens and structures. In this paper both in-plane and out-of-plane constraint effects of a crack in a reference reactor pressure vessel (RPV) subjected to pressurized thermal shocks (PTSs) are analyzed by two-parameter and three-parameter methods. The comparison between elastic and elastic–plastic analysis shows that the constraint effect varies with the material property. T₁₁ (the second term of William’s extension acting parallel to the crack plane) generally displays a reversed relation to the stress intensity factor (SIF) with the transient time, which indicates that the loading (SIF) plays an important role on the in-plane constraint effect. The thickness at the crack tip contributes more than the loading to the out-of-plane constraint, such that T₃₃ (the second term of William’s extension acting along the thickness) displays a similar relation to ε₃₃ (strain along the thickness direction) and a different relation to T₁₁ during the transient. The results demonstrate that both in-plane and out-of-plane constraint effect should be analyzed separately in order to describe precisely the stress distribution ahead of the crack tip.     

Blind-hole residual stress determination using optical interferometry

The in-plane method and the out-of-plane method are used to analyze blind-hole residual stress as measured by optical interferometry. The in-plane method, which constructs a relation between the in-plane displacement field and the residual stress released from blind-hole drilling, is applicable when the sensitivity vector of the interferometer used in the measuring system is parallel to the object surface. Three in-plane displacements obtained from one interference pattern are sufficient to determine the residual stress. The out-of-plane method, which establishes a new relation between the out-of-plane displacement field and the released residual stress, is suggested when the sensitivity vector is perpendicular to the object surface. Two relative out-of-plane displacements extracted from one interference pattern are sufficient to determine the residual stress. With the adoption of these two methods, interpolating calculation is not needed to determine the fringe order of each data point, since the selections of the required data points are flexible using these two methods. Two experiments, one for the in-plane method and the other for the out-of-plane method, were carried out to illustrate the applicability and usefulness of these two methods.     

3D Strain Field Evolution and Failure Mechanisms in Anisotropic Paperboard

Background

Experimental analyses of the 3D strain field evolution during loading allows for better understanding of deformation and failure mechanisms at the meso- and microscale in different materials. In order to understand the auxetic behaviour and delamination process in paperboard materials during tensile deformation, it is essential to study the out-of-plane component of the strain tensor that is, in contrast to previous 2D studies, only achievable in 3D.

Objective

The main objective of this study is to obtain a better understanding of the influence of different out-of-plane structures and in-plane material directions on the deformation and failure mechanisms at the meso- and microscale in paperboard samples.

Methods

X-ray tomography imaging during in-situ uniaxial tensile testing and Digital Volume Correlation analysis was performed to investigate the 3D strain field evolution and microscale mechanical behaviour in two different types of commercial paperboards and in two material directions. The evolution of sample properties such as the spatial variation in sample thickness, solid fraction and fibre orientation distribution were also obtained from the images. A comprehensive analysis of the full strain tensor in paperboards is lacking in previous research, and the influence of material directions and out-of-plane structures on 3D strain field patterns as well as the spatial and temporal quantification of the auxetic behaviour in paperboard are novel contributions.

Results

The results show that volumetric and deviatoric strain, dominated by the out-of-plane normal strain component of the strain tensor, localize in the out-of-plane centre already in the initial linear stress-strain regime. In-plane strain field patterns differ between samples loaded in the Machine Direction (MD) and Cross Direction (CD); in MD, strain localizes in a more well-defined zone close to the notches and the failure occurs abruptly at peak load, resulting in angular fracture paths extending through the stiffer surface planes of the samples. In CD, strain localizes in more horizontal and continuous bands between the notches and at peak load, fractures are not clearly visible at the surfaces of CD-tested samples that appear to fail internally through more well-distributed delamination.

Conclusions

In-plane strain localization preceded a local increase of sample thickness, i.e. the initiation of the delamination process, and at peak load, a dramatic increase in average sample thickening occurred. Different in-plane material directions affected the angles and continuity of the in-plane strain patterns as well as the sample and fracture properties at failure, while the out-of-plane structure affected how the strain fields distributed within the samples.

     

冲击载荷作用下柱壳链中的弹性波传播简化模型及其解析解

柱壳链能引起波形的弥散,具备操控波形的潜力。建立了柱壳链结构的等效连续介质模型和细观有限元模型,研究了质量块冲击作用下柱壳链中的弹性应力波传播过程及其几何弥散特性。基于考虑横向惯性修正的Rayleigh-Love波动方程,建立了柱壳链在质量块冲击下的控制方程,采用Laplace变换及其逆变换获得了位移场、速度场和应变场的解析解,所得结果与细观有限元模拟结果较好吻合。结果表明,在冲击过程中应变和速度峰值均逐渐减小,应变峰值、振荡幅度和波形前沿宽度与泊松比和惯性半径相关,泊松比和惯性半径越大,应变峰值越小,应变分布振荡越剧烈,波形前沿宽度越宽。     

基体开裂的正交复合材料对称层合板的粘弹性力学行为

本文讨论了含有均匀基体裂纹的正交复合材料对称层合板的线性粘弹性力学行为.采用二维剪切滞后模型并对其层间剪应力在厚度方向进行线性假设分布,求得层合板的平均应力应变的线弹性解,利用等效约束模型和经典层合理论可得到层合板因为含有基体裂纹而所引起的刚度退化现象.在弹性-粘弹性对应原理的基础上对其层合板的线粘弹性进行了讨论和研究.结果表明:层合板的松弛模量和蠕变泊松比随着时间的增加而减少,到达稳态后其值基本上是恒值.并跟Zocher的解析解和有限元数值解作了比较,发现结果非常吻合.     

Planar and non-planar non-linear dynamics of suspended cables under random in-plane loading-I. Single internal resonance

The non-linear interaction of the in-plane and out-of-plane motions of a suspended cable in the neighbourhood of 2:1 internal resonance under random loading is studied. The random loading acts externally on the in-plane mode, while the out-of-plane mode is non-linearly coupled with the in-plane mode. Any non-trivial motion of the out-of-plane mode is mainly due to this non-linear coupling, which becomes significant in the neighbourhood of internal resonance. The response statistics are estimated by employing the Fokker-Planck equation together with Gaussian and non-Gaussian closures. Monte-Carlo simulation is also used for numerical verification. Away from the internal resonance condition, the response is governed by the inplane motion, and the non-Gaussian closure solution is found to be in good agreement with numerical simulation results. The stochastic bifurcation of the out-of-plane mode is predicted by Gaussian and non-Gaussian closures, and by Monte-Carlo simulation. The non-Gaussian closure can only predict bounded solutions within a limited region. The on-off intermittency of the second mode is observed in the Monte-Carlo simulation over a small range of excitation level. The influence on response statistics of excitation level and cable parameters, such as internal detuning, damping ratios, and sag-to-span ratio, is reported.     

高温结构蠕变裂尖拘束效应

如何对蠕变裂纹扩展寿命进行准确预测和评价是高温结构完整性评定、寿命设计和运行维护中需要解决的核心问题.基于宏观单参数C?的蠕变断裂行为的评价方法,未有效纳入裂尖拘束效应的影响,因而其评价结果过于保守或非保守.目前国内外还未建立起有效纳入裂尖拘束效应的高温结构蠕变寿命评价的理论体系和技术方法,还没有纳入蠕变拘束效应的高温结构完整性评定规范.本文综述了作者近年来在高温蠕变断裂拘束效应方面的研究工作.主要包括:裂尖拘束对材料蠕变裂纹扩展行为的影响及机理;蠕变裂尖场和拘束参数R的定义和影响因素;载荷无关的蠕变拘束参数R? 的提出及其应用基础;承压管道表面裂纹的拘束参数R? 解及纳入裂尖拘束的蠕变寿命评价方法;试样与管道轴向裂纹蠕变裂尖拘束的关联;基于裂尖等效蠕变应变的面内与面外蠕变裂尖拘束的统一表征参数Ac的研究;材料拘束相关的蠕变裂纹扩展速率的建立;宽范围C? 区蠕变裂纹扩展速率及其拘束效应的数值预测;材料拘束对焊接接头蠕变裂纹扩展行为的影响及机理等.这些研究为建立纳入裂尖拘束效应的高温部件的蠕变裂纹扩展寿命评价方法奠定了理论和技术基础.论文对后续拟开展的工作也进行了展望.