Numerical analysis of macrocrack propagation along a bimaterial interface under dynamic loading processes

Advances in computing as well as measurement instrumentation have recently allowed for the investigation of a wider spectrum of physical phenomena in dynamic failure than previously possible. With increasing demand for specialized lightweight, high strength structures, failure of inhomogeneous solids has been receiving increased attention. Such inhomogeneous solids include structural composites such as bonded and sandwich structures, layered and composite materials as well as functionally graded solids. Many of such solids are composed of brittle constituents possessing substantial mismatch in wave speeds, and are bonded together with weak interfaces, which may serve as sites for catastrophic failure (Rosakis and Ravichandran (2000)).In the present study numerical analysis of macrocrack propagation along a bimaterial interface under dynamic loading processes is presented. A general constitutive model of elasto-viscoplastic damaged polycrystalline solids is developed within the thermodynamic framework of the rate type covariance structure with finite set of the internal state variables. A set of the internal state variables is assumed and interpreted such that the theory developed takes account of the effects as follows: (i) plastic non-normality; (ii) softening generated by microdamage mechanisms; (iii) thermomechanical coupling (thermal plastic softening and thermal expansion); (iv) rate sensitivity.To describe suitably the time and temperature dependent effects observed experimentally during dynamic loading processes the kinetics of microdamage has been modified. The relaxation time is used as a regularization parameter. By assuming that the relaxation time tends to zero, the rate independent elastic–plastic response can be obtained. The identification procedure is developed basing on the experimental observations. The finite difference method for regularized elasto-viscoplastic model is used. The edge-cracked bimaterial specimen is considered. In the initial configuration, the height of the specimen is equal to 30 cm, width is 12.5 cm and the length of the initial crack is equal to 2.5 cm. The length of the boundary over which impact is applied is equal to 5 cm, the rise time is fixed at 0.1 μs and the impact velocity is varied. The impact area is localized symmetrically or asymmetrically to the shorter axis of the specimen (symmetry axis of the cohesive band). Basing on the available data of recent experimental observation Rosakis et al. (1999) that have been carried out for relatively thin specimens both the plane stress and plane strain conditions are considered. The material of the specimen is AISI 4340 steel, while PMMA is the cohesive band, both modelled by thermo-elasto-viscoplastic constitutive equations with effects of isotropic hardening and softening generated by microdamage mechanisms and thermomechanical coupling. Fracture criterion based on the evolution of microdamage is assumed. Both, isothermal and adiabatic processes are considered.Particular attention is focused on the investigation of the interactions and reflections of stress waves and the influence of these waves on the propagation of macrocrack within the interface band. The propagation of the macroscopic crack within the material of the interface band for both symmetrical and asymmetrical impact cases has been investigated. It has been found that macrocrack-tip speeds vary from the shear wave speed to the dilatational wave speed of the material and is higher than the Rayleigh surface wave speed. This result is in accord with the experimental observations performed by Rosakis et al. (1999).     

应力波载荷作用下线弹性断裂过程的动态分析方法研究

利用Hopkinson单压杆实验装置 ,对材料的线弹性动态断裂特性进行了研究 ,建立了应力波载荷作用下动态裂纹起裂及扩展过程的动态分析方法 ,采用该方法可同时测得材料的动态裂纹起裂时间、断裂韧性和裂纹扩展速度。40Cr钢三点弯曲试样的实验结果表明 :该钢的动态裂纹扩展过程主要是减速过程 ,在2 2 5TPam /s的加载速率下 ,起裂时间为 2 8 0 0 s,最大裂纹扩展速度为 478 91m/s ,动态断裂韧性为6 3 12MPam。     

基于耦合算法的三维复杂结构冲击动力学特性

光滑粒子流体动力学-有限元耦合算法（SPH-FEM）较好地结合了SPH和FEM的优势,近年来逐渐被引入冲击动力学相关问题研究中。然而早期的研究对象多为单一材料的简单结构,所取得的研究成果距离实际工程应用仍有一定差距。为此,在总结前人工作的基础上,对SPH-FEM耦合算法进行适当改进,通过引入复合材料损伤模型,对复合材料蒙皮结构飞行器舱段结构进行建模计算,分析其在爆炸冲击激励下的冲击动力学特性。将数值计算结果与试验结果进行对比分析,验证该算法和模型的有效性和准确性,初步实现SPH-FEM的工程实际应用。最后总结了复合材料蒙皮结构飞行器在爆炸冲击激励下的一系列结构动态响应规律,以期为航天飞行器结构设计与防护提供参考。     

Damage localization using transient non-linear elastic wave spectroscopy on composite structures

To reduce the costs related to maintenance of aircraft structures, there is the need to develop new robust, accurate and reliable damage detection methods. A possible answer to this problem is offered by newly developed non-linear acoustic/ultrasonic techniques, which monitor the non-linear elastic wave propagation behaviour introduced by damage, to detect its presence and location.In this paper, a new transient non-linear elastic wave spectroscopy (TNEWS) is presented for the detection and localization of a scattered zone (damage) in a composite plate. The TNEWS analyses the uncorrelations between two structural dynamic responses generated by two different pulse excitation amplitudes by using a time-frequency coherence function. A numerical validation of the proposed method is presented. Damage was introduced and modelled using a multi-scale material constitutive model (Preisach-Mayergoyz space).The developed technique identified in a clear manner the faulted zone, showing its robustness to locate and characterize non-linear sources in composite materials     

Asymmetrical dynamic fracture model of bridging fiber pull-out of unidirectional composite materials

An elastic analysis of an internal crack with bridging fibers parallel to the free surface in an infinite orthotropic anisotropic elastic plane is studied, and asymmetrical dynamic fracture model of bridging fiber pull-out of unidirectional composite materials is presented for analyzing the distributions of stress and displacement with the internal asymmetrical crack under the loading conditions of an applied non-homogenous stress and the traction forces on crack faces yielded by the bridging fiber pull-out model. Thus the fiber failure is ascertained by maximum tensile stress, the fiber ruptures and hence the crack propagation should also appear in the modality of self-similarity. The formulation involves the development of a Riemann-Hilbert problem. Analytical solution of an asymmetrical propagation crack of unidirectional composite materials under the conditions of two increasing loads given is obtained, respectively. In terms of correlative material properties, the variable rule of dynamic stress intensity factor was depicted very well. After those analytical solutions were utilized by superposition theorem, the solutions of arbitrary complex problems could be gained.     

An asymmetrical dynamic model for bridging fiber pull-out of unidirectional composite materials

An elastic analysis of an internal crack with bridging fibers parallel to the free surface in an infinite orthotropic elastic plane is studied. An asymmetrical dynamic model for bridging fiber pull-out of unidirectional composite materials is presented for analyzing the distributions of stress and displacement with the internal asymmetrical crack under the loading conditions of an applied non-homogenous stress and the traction forces on crack faces yielded by the bridging fiber pull-out model. Thus the fiber failure is determined by maximum tensile stress, resulting in fiber rupture and hence the crack propagation would occur in a self-similarity manner. The formulation involves the development of a Riemann-Hilbert problem. Analytical solution of an asymmetrical propagation crack of unidirectional composite materials under the conditions of two moving loads given is obtained, respectively. After those analytical solutions were utilized by superposition theorem, the solutions of arbitrary complex problems could be obtained.     

冻融循环对含纯Ⅰ型裂隙围岩的动态起裂特性影响规律

以寒区隧道为工程背景研究在冻融循环作用下围岩内Ⅰ型裂纹的动态起裂特性演化规律,采用隧道模型试件作为研究对象,开展冻融循环试验与大尺度落锤冲击试验,得到岩石试件经历不同冻融循环次数后的相关力学参数,并在裂纹尖端粘贴裂纹扩展计(crack propagation gauge, CPG)测量预制裂纹的动态起裂时间。采用有限元软件建立相应的数值模型计算裂纹尖端的动态应力强度因子,采用试验-数值法计算动态起裂韧度,随后采用电镜对冻融循环后的试样进行扫描,研究冻融循环对岩石材料的细观损伤机制。研究结果表明:随着冻融循环次数的增加,岩石材料的纵波、横波波速与弹性模量逐渐减小,而泊松比逐渐增大;砂岩材料的动态起裂韧度随着冻融循环次数的增加逐渐减小,表征线性反比例的特性;材料内部的胶结物质会由于冻融循环的影响而流失,材料的孔隙和裂纹也随着冻融循环次数的增加而变多变大。     

含偏置裂纹三点弯曲梁的动态断裂行为研究

采用动态焦散线方法，对含偏置裂纹三点弯曲梁承受横向冲击的弯曲断裂行为进行了一系列动态断裂力学实验研究，分析了无量纲量ａ／ｌ的改变（ａ——初始裂纹偏离梁中心线的距离；ｌ——梁长度的一半）对于裂纹动态扩展行为（裂纹起始状态、裂纹尖端的复合应力强度因子、裂纹扩展速度、裂纹扩展轨迹）的影响，并借助动态光弹性应力分析，对应力波与扩展裂纹的相互作用以及应力波传播规律进行探讨．给出了裂纹尖端复合应力强度因子、裂纹扩展速度的变化、裂纹曲裂轨迹以及方向与梁中应力波传播的相互关系     

冲击加载下巷道内裂纹的扩展特性及破坏行为

为了开展含预制裂纹的巷道模型试样在冲击载荷下的动态断裂响应实验,选用青砂岩作为模型材料制作巷道模型试样,以可调速落锤冲击实验机作为冲击加载装置进行试样的动态断裂实验,分析冲击载荷作用下的巷道内裂纹的扩展规律。采用裂纹扩展计及应变片测试系统监测裂纹的起裂时间、扩展速度及止裂时间,并借助于AUTODYN、ABAQUS有限元数值分析软件对实验模型进行数值模拟,计算裂纹的动态起裂韧度、动态扩展韧度、动态止裂韧度等断裂参数。结果表明:巷道内裂纹在扩展路径过程中存在着明显的止裂现象;采用实验-数值方法能够较好地得出裂纹的起裂韧度、扩展韧度和止裂韧度等参数。另外,对止裂现象进行了初步的分析,讨论了试样内应力反射波与透射波对止裂问题的影响。     

DYNAMIC CRACK MODELS ON PROBLEM OF BRIDGING FIBER PULL-OUT OF COMPOSITE MATERIALS

An elastic analysis of an internal central crack with bridging fibers parallel to the free surface in an infinite orthotropic anisotropic elastic plane was performed. A dynamic model of bridging fiber pull-out of composite materials was presented. Resultingly the fiber failure is governed by maximum tensile stress, the fiber breaks and hence the crack extension should occur in self-similar fashion. By the methods of complex functions, the problem studied can be transformed into the dynamic model to the Reimann-Hilbert mixed boundary value problem, and a straightforward and easy analytical solution is presented. Analytical study on the crack propagation subjected to a ladder load and an instantaneous pulse loading is obtained respectively for orthotropic anisotropic body. By utilizing the solution, the concrete solutions of this model are attained by ways of superposition.     

The propagation of elastic stress waves in a conical shell under axial impulsive loading

The propagation of elastic stress waves in a conical shell subjected to axial impulsive loading is studied in this paper by means of the finite element calculation and model experiments. It is shown that there are two axisymmetrical elastic stress waves propagating with different velocities, i.e., the longitudinal wave and the bending wave. The attenuation of these waves while propagating along the shell surface is discussed. It is found in experiments that the bending wave is also generated when a longitudinal wave reflects from the fixed end of the shell, and both reflected waves will separate during the propagation due to their different velocities. Southwest Institute of Structural Mechanics     

半空间饱和土中圆柱形衬砌内荷载引起的动力响应解答

地下衬砌结构经常会受到内部动荷载的作用,内荷载引起的衬砌结构的动应力集中备受关注。采用Laplace变换和波函数展开法,对半空间饱和土中均布突加荷载作用下圆柱形衬砌结构的动力响应进行了研究。对以大半径凸圆弧来近似半空间表面,将半空间饱和土中的波动方程展开成无穷级数的形式,应用Graff加法公式进行坐标变换,根据连续条件和边界条件,确定波函数展开式中的未知系数,求得了半空间饱和土中均布突加荷载作用下圆柱形衬砌的动力响应解答,讨论了衬砌埋深对圆柱形衬砌动力响应的影响。该问题的求解,为地下结构的动力分析提供了一种有效方法。     

A dynamic model of bridging fiber pull-out of composite materials

An elastic analysis of an internal central crack with bridging fibers parallel to the free surface in an infinite orthotropic anisotropic elastic plane is carried out. In this paper a dynamic model of bridging fiber pull-out is presented for analyzing the distributions stress and displacement of composite materials with the internal central crack under the loading conditions of an applied non-uniform stress and the traction forces on crack faces yielded by the fiber pull-out model. Thus the fiber failure is determined by maximum tensile stress, the fiber breaks and hence the crack propagation should occur in self-similar fashion. By reducing the dynamic model to the Keldysh–Sedov mixed boundary value problem, a straightforward and easy analytical solution can be attained. When the crack extends, its fibers continue to break. Analytical study on the crack extension under the action of an inhomogeneous point force Px/t, Pt is obtained for orthotropic anisotropic body, respectively; and it can be utilized to attain the concrete solutions of the model by the ways of superposition.     

Diffraction of a shear plane wave in a compound elastic space with a half-infinite crack parallel to the inhomogeneity line

Problems of stress wave propagation and diffraction in elastic inhomogeneous media are undoubtedly of interest to scientists from the viewpoint of investigation of fundamental laws of dynamic processes and of the use of the results in technical and technological applications. The paper deals with the dynamic contact problem of shear plane wave diffraction at the edge of a semi-infinite crack in a compound space consisting of two elastic half-spaces. The questions related to the onset of surface waves and the wave field behavior in far-field regions are also considered.     

金属铅在斜冲击波对碰加载下动态行为的数值模拟

高强度冲击加载作用下金属材料的动态物理行为是当前冲击波领域基础研究和工程应用最为关注的焦点。采用光滑粒子法(SPH)开展不同位置起爆诱发的斜冲击波对碰加载金属铅的二维数值模拟研究,得到了金属铅内入射斜冲击波的角度和强度,并利用极曲线方法理论上导出发生马赫反射时的临界入射角和入射马赫数关系。根据计算结果可知,金属铅内入射斜波对碰后将发生马赫反射。随着起爆位置与金属铅表面距离的增加,不仅金属铅内入射冲击波强度和入射角增加,而且形成的马赫杆宽度也在增加。由自由面速度剖面给出了马赫杆宽度及张角,结果与理论预测的结果吻合较好。     

球壳中球形夹杂对SV波的三维散射与动应力集中研究

夹杂将导致结构应力集中,是降低结构承载能力重要影响因素,尤其是动载作用情况下,弹性波衍射和叠加将加剧应力集中程度.弹性波衍射方程建立和求解非常复杂,目前主要研究对象集中在二维模型情况,三维有限域内夹杂引起的动应力集中现象在大型结构中比较常见,有界域边界不仅作为边界条件,同时也是散射波波源,提高了求解难度.一般通过近似方...     

The dynamic behaviour of a piezoelectric actuator bonded to an anisotropic elastic medium

Surface-bonded piezoelectric actuators can be used to generate elastic waves for monitoring damages of composite materials. This paper provides an analytical and numerical study to simulate wave propagation in an anisotropic medium induced by surface-bonded piezocermic actuators under high-frequency electric loads. Based on a one-dimensional actuator model, the dynamic load transfer between a piezoceramic actuator and an anisotropic elastic medium under in-plane mechanical and electrical loading is obtained. The wave propagation induced by the surface-bonded actuator is also studied in detail by using Fourier transform technique and solving the resulting integral equations in terms of the interfacial shear stress. Typical examples are provided to show effects of the geometry, the material combination, the loading frequency and the material anisotropy of the composite upon the load transfer and the resulting wave propagation.     

压剪联合冲击下K9玻璃中的失效波

对K9玻璃进行了冲击速度为150~400 m/s,倾斜角为10、15的纵剖试样斜撞击实验。结果表明，在加载剪切横波(S+)和卸载纵波(P-)之间有波速超过了纵波波速的波阵面存在，此波的存在对卸载纵波和卸载横波的幅值有一定的影响。由于失效波的产生与材料的表面性质密切相关，为消除纵剖试样中间界面的影响，模拟VISAR的实验条件，进行了冲击速度为70~300 m/s、倾斜角为10的横剖试样斜撞击实验。可以确定在此冲击范围失效波的波速在0.98~1.4 km/s，产生失效波的临界状态为:压应力0.86~1.01 GPa，对应剪应力0.053~0.071 GPa。表明剪应力的存在极大降低了失效波产生的阈值。在此基础上初步分析了撞击面的动摩擦因数和相对滑移速度。     

爆炸应力波与裂纹作用实验研究

为研究爆炸应力波与裂纹相互作用机理,利用透射式爆炸动态焦散线光学实验系统研究了预制水平静态裂纹和切缝药包炮孔爆破产生的水平运动裂纹受正入射爆炸动载作用后动态特性的变化规律。结果表明:正入射爆炸应力波与静止裂纹作用时,爆炸应力波P波使得裂纹先闭合后张开,S波在裂纹壁面形成波浪状散斑上下交替向外扩展;运动裂纹尖端应力场对静止裂纹的起裂和扩展有重要影响。后爆孔爆炸应力波对先爆孔产生的水平定向运动裂纹尖端动力学特性影响显著。当爆炸应力波与运动裂纹同向时,P波使得裂纹扩展速度和应力强度因子KI^d先减小后增大,S波促进了裂纹的扩展,波与裂纹作用之后,裂纹扩展速度增大;当爆炸应力波与运动裂纹反向时,P波抑制了运动裂纹的扩展,波与裂纹作用之后,裂纹扩展速度和应力强度因子KI^d均逐渐降低。     

SEMI-ELLIPTIC SURFACE CRACK IN AN ELASTIC SOLID WITH FINITE SIZE UNDER IMPACT LOADING

In this paper a semi-elliptic surface crack problem in an elastic solid of finite size under impact loading is investigated. An analysis is performed by means of fracture dynamics and the finite element method, and a three-dimensional finite element program is developed to compute the dynamic stress intensity factor. The results reveal that the effects of the solid's boundary surface, crack surface, material inertia and stress wave interactions play significant roles in dynamic fracture.