基于等效黏性弹簧的黏弹性Timoshenko裂纹梁弯曲解析解

考虑裂纹的缝隙和黏性效应,将梁中横向裂纹等效为黏弹性扭转弹簧,利用广义Delta函数,给出了Laplace变换域内裂纹梁的等效抗弯刚度,得到了具有任意开闭裂纹数目且满足标准线性固体黏弹性本构的Timoshenko梁在时间域内的弯曲变形显式解析通解．在此基础上,通过两个数值算例,分析了时间、梁跨高比和裂纹深度等参数对黏弹性Timoshenko开裂纹梁弯曲变形的影响．结果表明：裂纹黏性对Timoshenko裂纹梁的弯曲具有显著的影响．相比于裂纹的弹性扭转弹簧模型,考虑裂纹黏性效应的黏弹性Timoshenko裂纹梁在裂纹处挠度尖点和转角跳跃现象十分明显．另外,由于横向剪切引起的附加变形,Timoshenko裂纹梁的稳态挠度与Euler-Bernoulli梁挠度的差值为常数,其大小与裂纹模型、梁跨高比或裂纹深度无关,这些结果对梁裂纹无损检测具有指导意义．     

黏弹性界面断裂分析的增量“加料”有限元法

提出了一种适用于黏弹性界面裂纹问题的增量“加料” 有限元方法. 利用弹性界面裂纹尖端位移场的解答,通过对应原理和拉普拉斯逆变换近似方法,得到了黏弹性界面裂纹的尖端位移场. 用该位移场构造了黏弹性界面裂纹“加料” 单元和过渡单元位移模式,推导了增量“加料” 有限元方程,求解有限元方程可获得应力强度因子和应变能释放率等断裂参量. 建立了典型黏弹性界面裂纹平面问题“加料” 有限元模型,计算结果表明,对于弹性/黏弹性界面裂纹和黏弹性/黏弹性界面裂纹,该方法都能得到相当精确地断裂参量,并能很好地反映蠕变和松弛特性,可推广应用于黏弹性界面断裂问题的计算分析.      

Viscoelastic Deformation of a Reinforced Plate with a Crack

The deformation of a viscoelastic reinforced composite is studied. The composite has an axis of elastic symmetry and consists of transversally isotropic fibers and a viscoelastic matrix, which differ by the volume concentration and mechanical characteristics. The material is modeled by a transversely isotropic homogeneous linearly viscoelastic medium with some averaged characteristics. A plate fabricated from the composite in question is weakened by a through mode I crack and is subjected to constant tensile forces. The viscoelastic properties of the matrix material are described by a convolution operator. The Volterra principle is used to derive expressions for the viscoelastic characteristics and crack opening. The irrational function of the integral operator that describes the crack opening is expanded into an operator continued fraction and is represented as the sum of base operators     

功能梯度材料的黏弹性断裂问题

功能梯度材料(FGM)是一种不同于传统复合材料的新型工程复合材料 [1], 国内外关于FGM的断裂力学方面的研究发展非常迅速. 关于FGM静态裂纹问题,学者们研究了不同类型裂纹尖端场的应力强度因子 [2-5], 探讨了有限长裂纹在不用载荷作用下的传播等问题. 而关于动态裂纹问题,也已经取得很大成就 [6-9]. FGM一个很重要的应用是高温结构材料,在强大的热环境中,很多材料都呈现出黏弹性. 因此,研究FGM的黏弹性断裂力学非常具有实际价值.对此,众多研究 [10-14]提出不同的分析模型,并在不同受载条件,通过理论计算,分析了黏弹性裂纹尖端场的力学 行为.本文考查了功能梯度材料板条中界面裂纹垂直于梯度方向时的黏弹性断裂问题,首先利用有限元法求解线弹性功能梯度材料板条的裂纹尖端场,然后根据黏弹性的对应性原理,求解出黏弹性功能梯度材料板条裂纹问题的应力场强度因子.      

Crack propagation in the power-law nonlinear viscoelastic material

I.IntroductionCrazingdamageisacommonphenon1enonoffractureofpolymericmaterials.Theformationofcrazezoneisamid-stateinthefractureprocessofthematerialsfromperfectstatetofaiIurc.Microscopically,inthisregionthereexistssomefibrilslinkingthetwocracksurfacesandres…     

Photoelastic modeling of the fracture of viscoelastic orthotropic plates with a crack

The well-known equations of photoelasticity of linear viscoelastic bodies are used to describe the photoelastic behavior of a viscoelastic orthotropic plate with a crack. Expressions for the stress intensity factors (SIFs) at the crack tip are obtained using photoelastic measurements. The time dependence of the SIFs is analyzed and shown to be determined by the angles between directions of the crack and tension     

Some basic problems in anisotropic bimaterials with a viscoelastic interface

This research is devoted to the study of anisotropic bimaterials with Kelvin-type viscoelastic interface under antiplane deformations. First we derive the Green’s function for a bimaterial with a Kelvin-type viscoelastic interface subjected to an antiplane force and a screw dislocation by means of the complex variable method. Explicit expressions are derived for the time-dependent stress field induced by the antiplane force and screw dislocation. Also presented is the time-dependent image force acting on the screw dislocation due to its interaction with the Kelvin-type viscoelastic interface. Second we investigate a rectangular inclusion with uniform antiplane eigenstrains embedded in one of the two bonded anisotropic half-planes by virtue of the derived Green’s function for a line force. The explicit expressions for the time-dependent stress field induced by the rectangular inclusion are obtained in terms of the simple logarithmic and exponential integral functions. It is observed that in general the stresses exhibit the logarithmic singularity at the four corners of the rectangular inclusion. Our results also show that when one side of the rectangular inclusion lies on the viscoelastic interface, the interfacial tractions are still regular at the two corners of the inclusion which are located on the interface. Last we address a finite Griffith crack normal to the viscoelastic interface by means of the obtained Green’s function for a screw dislocation. The crack problem is formulated in terms of a resulting singular integral equation which is solved numerically. The time-dependent stress intensity factors at the two crack tips are obtained and some interesting features are discussed.     

Specification of crack opening and surface sliding in orthotropic viscoelastic material

Modelling of crack opening and surface sliding in an orthotropic viscoelastic material is made by introducing two coefficients: one for the surface displacement and surface friction. The material possesses orthotropy in two dimensions and viscoelastic property consisting of a Kelvin element in series with a spring. The method of Laplace transform is applied to obtain a closed form solution to the problem. Explicit expressions of Mode I and II stress intensity factors are computed together with crack surface opening. Trade-off between the Mode I and II stress intensity factors depends on the nature of material orthotropy.     

Subsonic and intersonic mode II crack propagation with a rate-dependent cohesive zone

A recent experimental study has demonstrated the attainability of intersonic shear crack growth along weak planes in otherwise homogeneous, isotropic, linear elastic solids subjected to remote loading conditions (Rosakis et al., Science 284 (5418) (1999) 1337). The relevant experimental observations are summarized briefly here and the conditions governing the attainment of intersonic crack speeds are examined. Motivated by experimental observations, subsonic and intersonic mode II crack propagation with a rate-dependent cohesive zone is subsequently analyzed. A cohesive law is assumed, wherein the cohesive shear traction is either a constant or varies linearly with the local sliding rate. Complete decohesion is assumed to occur when the crack tip sliding displacement reaches a material-specific critical value. Closed form expressions are obtained for the near-tip fields. With a cohesive zone of finite size, it is found that the dynamic energy release rate is finite through out the intersonic regime. Crack tip stability issues are addressed and favorable speed regimes are identified. The influence of shear strength of the crack plane and of a rate parameter on crack propagation behavior is also investigated. The isochromatic fringe patterns predicted by the analytical solution are compared with the experimental observations of Rosakis et al. (1999) and comments are made on the validity of the proposed model.     

The scattering of general SH plane wave by interface crack between two dissimilar viscoelastic bodies

The scattering of general SH plane wave by an interface crack between two dissimilar viscoelastic bodies is studied and the dynamic stress intensity factor at the crack-tip is computed. The scattering problem can be decomposed into two problems: one is the reflection and refraction problem of general SH plane waves at perfect interface (with no crack); another is the scattering problem due to the existence of crack. For the first problem, the viscoelastic wave equation, displacement and stress continuity conditions across the interface are used to obtain the shear stress distribution at the interface. For the second problem, the integral transformation method is used to reduce the scattering problem into dual integral equations. Then, the dual integral equations are transformed into the Cauchy singular integral equation of first kind by introduction of the crack dislocation density function. Finally, the singular integral equation is solved by Kurtz's piecewise continuous function method. As a consequence, the crack opening displacement and dynamic stress intensity factor are obtained. At the end of the paper, a numerical example is given. The effects of incident angle, incident frequency and viscoelastic material parameters are analyzed. It is found that there is a frequency region for viscoelastic material within which the viscoelastic effects cannot be ignored. This work was supported by the National Natural Science Foundation of China (No.19772064) and by the project of CAS KJ 951-1-20     

Superposition principles for small amplitude oscillatory shearing of nematic mesophases

The linear viscoelasticty of Leslie-Ericksen monodomain liquid crystals subjected to a bend distortion through a small amplitude oscillatory shear flow driven by harmonic wall stress is analyzed, using numerical and asymptotic methods. The viscoelastic material functions were derived using a new scaling approach that extracts the material parameters that control superposition. Small and high frequency superposition schemes for linear viscoleasticity were derived. The schemes were successfully applied to collapse the predicted loss and storage linear viscoelastic moduli of seven experimental data sets. Comparisons between different shear flows (simple shear and capillary Poiseuille) and different director distortion modes (splay and bend) shows that the superposition schemes are applicable to shear flows in any single director distortion mode.     

Subcritical stable growth of a penny-shaped crack in an aging viscoelastic body with cylindrical anisotropy

The paper studies the subcritical stable growth of a penny-shaped mode I crack in an aging linear viscoelastic body with cylindrical anisotropy under tensile loads applied at infinity. The equations of quasistatic crack growth are derived using the Volterra principle for a modified Dugdale model and the critical crack opening displacement as a failure criterion. The kernel of the operator that describes the viscoelastic crack opening is found using operator-valued continued fractions. As an example, the subcritical stable growth of a crack in a specific material is studied numerically __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 1, pp. 82–92, January 2007.     

粘弹性Timoshenko梁的变分原理和静动力学行为分析

从线性，各向同性，均匀粘弹性材料的Boltzmann本构定律出发，通过Laplace变换及其反变换，由三维积分型本构关系给出了Timoshenko梁的本构关系，并由此建立了小挠度情况下粘弹性Timoshenko梁的静动力学行为分析的数学模型，一个积分－偏微分方程组的初边值问题。同时，采用卷积，建立了相应的简化Gurtin型变分原理。给出了两个算例，考查了梁的厚度h与梁的长度l之比β对梁的力学行为的影响。     

粘弹性界面裂纹奇异场

对于许多粘弹性问题,通常可以利用对应性原理,即由弹性问题的结果得到对应的粘弹性问题在拉普拉斯变换域内的解,再通过反演变换求得最终时域中的解.但是,由于界面裂纹场存在着振荡奇异性,弹性问题解的形式就已经非常复杂,对应的粘弹性问题要通过反演变换直接求得准确的解析解几乎是不可能的.本文在利用对应性原理时做了更简单的准静态处理,即将弹性结果中的材料参数用粘弹性材料参数做对应替代,得到了粘弹性界面裂纹场近似的经典解,并与有限元分析结果作了比较.同时,利用Comninou接触模型,对粘弹性界面裂纹在远场拉剪混合加载情况下的裂尖应力场和接触区做了考察,并与经典解作了比较.     

Fracture mechanics of plates and shells applied to fail-safe analysis of fuselage Part I: Theory

In this paper, a general theory on the asymptotic field near the crack tip for plates and shells with and without shear deformation effect is established. It is found that four stress intensity factors, two for symmetrical and antisymmetrical stretching and two for symmetrical and antisymmetrical bending, are required to describe arbitrary asymptotic fields near the crack tip for plates without shear deformation. An additional stress intensity factor is required for the transverse shearing force induced by antisymmetrical bending when the shear deformation is included in the analysis. It is also proven by means of the complex variable technique that for problems of plates with shear deformation, there exist similarities in the asymptotic expressions of moments and membrane forces and also in the asymptotic expressions of in-plane displacements and rotations of the mid-surface. The energy release rate associated with crack growth in the direction of the crack line can be expressed in terms of stress intensity factors by means of Irwin's method of work and energy associated with a virtual crack extension. A combined stress intensity factor can be defined through the total energy release rate. The theory of the fracture of plates is generalized and applied to the study of problems in the fracture of shells. An example of an infinitely long cylindrical shell with a circumferential crack subjected to remote axial tension is given to demonstrate the application of the theory and to test the accuracy of the numerical analysis used for the problem.     

Dynamic crack propagation in a viscoelastic strip

The dynamic propagation of a semi-infinite crack in a finite linear viscoelastic strip subjected to Mode I loading is investigated. Through the use of integral transforms the problem is reduced to solving a Wiener-Hopf equation. The asymptotic properties of the transforms are exploited to establish the stress intensity factor. Plane-stress and plane-strain stress intensity factors as a function of crack speed for both fully-clamped and shear-free lateral boundaries are presented for the standard linear viscoelastic solid. Comparisons are made with previously obtained asymptotic stress intensity factors and with stress intensity factors for the equivalent elastic strips.     

Quasi-static growth and dynamic extension of anti-plane shear crack in periodically nonhomogenous elastic material

The quasi-static growth of an anti-plane shear crack in a periodically nonhomogeneous elastic material is considered. An effective velocity dependent energy release rate is found for a semi-infinite crack. Approximate expressions of stress intensity factors for in-plane cracks in periodically nonhomogeneous material are also given.     

An outline of the non-linear viscoelastic behaviour of wheat flour dough in shear

Creep and creep recovery, stress relaxation and small- and large-amplitude oscillatory shear experiments have been used to study the steady-state flow behaviour and the transient viscoelastic response of wheat flour dough in shear over large ranges of time, stress and strain. The results are discussed with reference to the limited body of reliable literature data. Dough does display a linear viscoelastic domain. The complex character of its non-linear viscoelastic properties is essentially due to the extremely low shear rate limit of the initial Newtonian plateau and to the onset of time-dependent flow behaviour above a certain strain threshold, which explain qualitatively the discrepancies observed in certain cases on a part of the range of the rheological variables explored, despite global self-consistency of the results. Comparison of gluten and dough linear viscoelastic properties shows that dough cannot be viewed simply as a concentrated suspension of starch granules in the hydrated viscoelastic gluten matrix.Paper presented at the second Annual European Rheology Conference (AERC 2005) held in Grenoble, France on April 21–23, 2005.     

Mapping of the relaxation patterns of polymer melts with linear flexible molecules of uniform length

The linear viscoelastic material functions of linear flexible polymers of uniform length are calculated from the BSW spectrum (Baumgaertel et al., 1990, 1992), and explicit analytic expressions are presented for several of the most common material functions for transient and dynamic experiments. However, numerical calculations are presented whenever needed. The BSW spectrum was determined from experimental G, G data of two sets of molten polymers of narrow molecular weight distribution, polystyrene and polybutadiene. The purpose of the mapping is to show a wide range of viscoelastic behavior which otherwise is not available in such comprehensive form. Experimental check of these predictions is still needed in most cases. Also, some insight into the predictions for the non-linear (including the non-equilibrium) viscoelastic behavior is achieved by studying two particular experiments: the start-up of uniaxial extension at constant rate and the start-up of shear flow at constant rate.Dedicated to Prof. John D. Ferry on the occasion of his 85th birthday.     

黏弹性材料界面裂纹应力场奇性分析

研究两半无限大黏弹性体间Griffith界面裂纹在简谐载荷作用下裂纹尖端动应力场的奇异特性.通过引入裂纹张开位移和裂纹位错密度函数,相应的混合边值问题归结为一组耦合的奇异积分方程.渐近分析表明裂尖动应力场的奇异特征完全包含在奇异积分方程的基本解中.通过对基本解的深入分析发现黏弹性材料界面裂纹裂尖动应力场具有与材料参数和外载荷频率相关的振荡奇异特性.以标准线性固体黏弹材料为例讨论了材料参数和载荷频率对奇性指数和振荡指数的影响.