冲击载荷下含表面裂纹圆柱壳体的动态断裂

动态载荷下含表面裂纹的有限尺寸构件的断裂问题在工程实践中有着重要意义,但由于此类问题非常复杂,目前还不能求得解析解。本文针对含轴向半椭圆盘状表面裂纹的圆柱壳体,应用有限元法研究了动态载荷下其断裂问题,计算了动态应力强度因子与静态应力强度因子的比值KdIyn(t)/KsIta。从计算结果可以得出,比值KdIyn(t)/KsIta与结构和裂纹的尺寸有关,而与冲击载荷的大小无关。本文所得结果在一定程度上揭示了圆柱壳体表面、裂纹面、物质惯性和弹性波的相互作用及其对动态断裂的影响。     

三点弯曲试样动态冲击特性的有限元分析

本文使用动态有限元技术，对于两种不同几何尺寸，两种不同材料的三点弯曲试样在三类七种不同冲击载荷作用下的动态响应进行了分析，求得了动态应力强度因子随时间的变化规律。并与准静态应力强度因子进行了比较。计算结果表明：将冲击载荷历史代入静态公式确定动态应力强度因子的做法是不正确的，要求得动态应力强度因子，必须对试样进行完全的动态分析。当材料的Ｅ／ρ值相同时，动态应力强度因子的响应曲线完全相同。而动态应力强度因子分别与加载点的位移及裂纹的张开位移之间存在着与准静态情况下各自相同的线性关系。这与资料［５］［６］中的结论完全相同。     

界面裂纹问题中的通用权函数法

热载荷和机械载荷共同作用下复合材料中的裂纹扩展往往发生在界面处.传统求解热冲击及机械载荷共同作用下界面裂纹尖端的应力强度因子的数值方法(如有限元、边界元法等),计算工作量大、效率低.通用权函数与时间无关,运用通用权函数法可以免除对每个时刻的应力分析,计算效率可得到很大提高.本文将通用权函数法推广到求解热载荷和机械载荷共同作用下界面裂纹尖端的应力强度因子过渡过程的问题中,推导出求解平面双材料界面裂纹问题应力强度因子的通用权函数法计算格式.基于此格式,计算热载荷和机械载荷共同作用下界面裂纹尖端的应力强度因子.通过实例计算比较,表明此方法得到的结果可以达到与相互作用积分法相当的工程应用精度.最后,应用此方法研究了热障涂层受热冲击及表面力共同作用时裂纹长度以及涂层厚度对应力强度因子的影响.结果表明:在一定边界条件下,当热障涂层中存在边缘裂纹时,随着涂层厚度的增加,更容易导致裂纹的扩展和涂层的剥落.     

裂纹表面承受均布载荷时的应力强度因子及裂纹张开位移的边界配…

本文针对裂纹表面承受载荷时的应力条件，提出了新的应力函数，对于各种裂纹模型，各种边界条件，各种边界形状，裂纹表面自由或承受均布载荷等均适用。并利用边界配位法，计算了裂纹表面承受均布载荷的方型板内中心裂纹的应力强度因子及裂纹的张开位移。     

裂纹表面承受均布载荷时的应力强度因子及裂纹张开位移的边界配位法

本文针对裂纹表面承受载荷时的应力条件，提出了新的应力函数，对于各种裂纹模型、各种边界条件、各种边界形状、裂纹表面自由或承受均布载荷等均适用。并利用边界配位法，计算了裂纹表面承受均布载荷时的方型板内中心裂纹的应力强度因子（ＳＩＦ）及裂纹的张开位移（ＣＯＤ）。     

冲击载荷作用下裂纹动态响应的数值模拟

对垂直、剪切以及斜向等各种冲击载荷作用下裂纹的动态响应进行了数值模拟,得到了一系列随时间变化的动态应力场以及应变场图;根据其定义,计算出了相应的动态应力强度因子,进而分析了斜向载荷作用下裂纹起裂情况,并对最优断裂问题进行了阐述。     

复合材料的Ⅲ型动态断裂力学分析

本文综合利用 Laplace 变换和 Fourier 变换,从理论上分析了含有贯穿裂纹的层合复合材料在动载荷作用下裂纹尖端的Ⅲ型动态应力强度因子.并以冲击载荷为算例,考查了外层的厚度、拉、剪模量的变化对裂纹尖端应力强度因子的影响,为改善复合材料的断裂动力学性能提供了参考依据.     

脆性岩石断裂破坏机理的边界配位法分析

针对裂纹表面承受载荷时的应力条件，提出了新的应力函数，该应力函数对于各种裂纹模型、各种边界条件、各种边界形状、裂纹表面自由或承受均布载荷等均适用．并利用边界配位法，计算了在压缩载荷下，岩石内部裂纹的应力强度因子（ＳＩＦ），给出了关于岩石断裂破坏的一些新结论     

三点弯曲试样动态冲击特性的有限元分析

本文使用动态有限元技术，对两种不同几何尺寸，两种不同材料的三点弯曲试样在三类七种不同冲击载荷作用下的动态响应进行了分析，求得了动态应力强度因子随时间的变化规律，并与准静态应力强度因子进行了比较，计算结果表明：半冲击载荷历史代入静态公式确定动态应力强度因子的做法是不正确的，要求得动态应力强度因子，必须对试样进行完全的动态分析，当材料的Ｅ／ρ值相同时，动态应力强度因子的响应曲线完全相同，而动态应力强度     

含双共线裂纹板的弹性分析

本文研究含双共线裂纹无限大平板的平面弹性问题,目的是要确定当裂纹上下表面受到任意的张开型载荷(正应力)作用时,板内的应力状态,特别是裂纹尖端上的应力强度因子.已有的若干文献考虑的是均布正应力情况下的问题,只是本文所处理的情况的一个特例.1.问题和基本解答考虑图1所示的含裂纹无限大平板.用t 表示裂纹位置线L 上的坐标x,我们要确定当裂纹上下表面受到任意的正应力σ(t)作用时,板内的应力和应力强度因子的计算公式.为此,就要解出由下列边界条件确定的一个平面弹性问题:     

Full field analysis of an anti-plane interface crack subjected to dynamic body forces

In this study, the transient full field response of an interface crack between two different media subjected to dynamic body force at one material is investigated. For time t < 0, the bimaterial medium is stress free and at rest. At t = 0, a concentrated anti-plane dynamic point loading is applied at the medium as shown in Fig. 1. The total wave field is due to the effect of this point loading and the scattering of the incident waves by the interface crack. An alternative methodology that is different from the conventional superposition method is used to construct the reflected, refracted and diffracted wave fields. A useful fundamental solution is proposed in this study and the full field solution is determined by superposition of the fundamental solution in the Laplace transform domain. The proposed fundamental problem is the problem of applying an exponentially distributed traction (in the Laplace transform domain) on the interfacial crack faces. The Cagniard–de Hoop method of Laplace inversion is used to obtain the transient solution in time domain. Exact transient closed form solutions for stresses and stress intensity factors are obtained. Numerical results for the time history of stresses and stress intensity factors during the transient process are discussed in detail.     

The Stress Intensity Factor for an Elliptic Crack Under a Piecewise-Constant Load

A closed approximate solution is given for the stress intensity factor (SIF) at the front of an elliptic crack under uniform loading on the elliptic ring. The solution is constructed by modifying Rice's variational formula that integrally relates the SIFs for two different states on the crack faces. The first state is a combination of the second and principal states such that the total SIF is zero. In this case, Rice's formula reduces to the virtual-displacement principle for a combined loading on the crack faces.     

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.     

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

通过三维有限元计算来研究有限宽度、有限厚度含有穿透裂纹板的裂纹前缘应力场，从中找出应力强度因子与板的厚度、裂纹长度之间的关系，同时还分析了裂尖的三维约束程度和三维约束区的大小。分析结果表明：应力强度因子沿厚度的分布是不均匀的，应力强度因子的最大值及其位置与厚度有关；有限厚度板中面应力强度因子(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影响区域大小影响较小。     

Thermostressed state of a piezoelectric bodywith a plane crack under symmetric thermal load

The paper addresses a thermoelectroelastic problem for a piezoelectric body with an arbitrarily shaped plane crack in a plane perpendicular to the polarization axis under a symmetric thermal load. A relationship between the intensity factors for stress (SIF) and electric displacement (EDIF) in an infinite piezoceramic body with a crack under a thermal load and the SIF for a purely elastic body with a crack of the same shape under a mechanical load is established. This makes it possible to find the SIF and EDIF for an electroelastic material from the elastic solution without the need to solve specific problems of thermoelasticity. The SIF and EDIF for a piezoceramic body with an elliptic crack and linear distribution of temperature over the crack surface are found as an example __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 3, pp. 96–108, March 2008.     

Parameter identification for cracks in elastic plate structures based on remote strain fields

A method for the detection of cracks in plate structures is presented. In contrast to most of the common monitoring concepts taking advantage of the reflection of elastic waves at crack faces, the presented approach is based on the strain measured at different locations on the surface of the structure. This allows both the identification of crack position parameters, such as length, location and angles with respect to a reference coordinate system and the calculation of stress intensity factors (SIF). The solution of the direct problem is performed on the basis of the BFM (body force method). The inverse problem is solved applying the particle swarm optimization (PSO) algorithm. The BFM is based on the principle of linear superposition which allows the calculation of the strain field in a cracked body. The strain at an arbitrary point in the structure is replaced by the strain provided by body force doublets in the uncracked structure. The doublets as well as external loads are parameters which have to be determined solving the inverse problem by minimizing a fitness function, which is defined by a square sum of residuals between measured strain distributions and computed ones for an assumed crack. The PSO algorithm applied to the fitness function operates on the basis of a swarm of candidate solutions. Once knowing loading and crack parameters, the SIF can be determined.     

Magnetostrictive/electrostrictive fracture of the piezomagnetic and piezoelectric layers in a multiferroic composite: Anti-plane case

The main purpose of the present work is to study the influences of magnetostriction, electrostriction and piezomagnetic/piezoelectric stiffening on the fracture behavior of a layered multiferroic composite. For comparison, it is assumed that there is a crack, parallel to the interface, in each layer. Methods of cosine transform and Cauchy singular integral equations are used to solve the crack problem. Numerical results of the stress intensity factor (SIF) are provided and the computational accuracy is demonstrated. Discussion on the numerical results indicates that the multiferroic composite consisting of cobalt ferrite and barium titanate layers are more prone to fracture under electric loading than under magnetic loading. In the case of magnetostriction, to increase the shear modulus of the piezomagnetic layer would raise the SIF; but to increase that of the piezoelectric layer would reduce the SIF; in the case of electrostriction, inverse results are obtained. Piezomagnetic stiffening can affect the SIF when the composite is under electrostriction; piezoelectric stiffening can influence the SIF if the composite is under magnetostriction. In addition, it is also revealed that two parallel equal cracks may shield each other even if an interface exists between them.     

THE MIXED STATE HAMILTONIAN DYNAMIC ELEMENT AND A SEMI-ANALYTICAL SOLUTION FOR THE ANALYSIS OF THICK LAMINATED COMPOSITE PLATES

Through introducing the Laplace transformation in the time direction,the mixed state Hamilton canonical equation and a semi-analytical solution arepresented for analyzing the dynamic response of laminated composite plates.Thismethod accounts for the separation of variables,the finite element discretization can beemployed in the plane of laminar,and the exact solution in the thickness direction isderived by the state space control method.To apply the transfer matrix method,therelational expression at the top and bottom surface is established.So the generalsolution in transformation space is deduced by the spot method.By the application ofinversion of Laplace transformation,the transient displacements and stresses can bederived.     

Transient response of a piezoelectric ceramic with two coplanar cracks under electromechanical impact

The transient response of two coplanar cracks in a piezoelectric ceramic under antiplane mechanical and inplane electric impacting loads is investigated in the present paper. Laplace and Fourier transforms are used to reduce the mixed boundary value problems to Cauchy-type singular integral equations in Laplace transform domain, which are solved numerically. The dynamic stress and electric displacement factors are obtained as the functions of time and geometry parameters. The present study shows that the presence of the dynamic electric field will impede or enhance the propagation of the crack in piezoelectric ceramics at different stages of the dynamic electromechanical load. Moreover, the electromechanical response is greatly affected by the ratio of the space of the cracks and the crack length.     

Multiple axial cracks in a coated hollow cylinder due to thermal shock

The transient thermal stress problem of an inner-surface-coated hollow cylinder with multiple pre-existing surface cracks contained in the coating is considered. The transient temperature, induced thermal stress, and the crack tip stress intensity factor (SIF) are calculated for the cylinder via finite element method (FEM), which is exposed to convective cooling from the inner surface. As an example, the material pair of a chromium coating and an underlying steel substrate 30CrNi2MoVA is particularly evaluated. Numerical results are obtained for the stress intensity factors as a function of normalized quantities such as time, crack length, convection severity, material constants and crack spacing.