密度梯度柱壳链的弹性波传播特性研究

均匀圆柱壳链可以调控弹性波传播,引入密度梯度有望进一步提高波形调控能力.通过建立密度梯度柱壳链的细观有限元模型和连续介质模型,研究了密度梯度柱壳链的弹性波传播特性.通过将密度梯度柱壳链等效为变密度连续介质弹性杆,建立了其在应力脉冲作用下的控制方程.运用拉普拉斯积分变换方法,考虑杆中密度遵循线性分布,获得了方程的解析解.以三角形应力脉冲作用为例,通过与细观有限元模拟结果比较,发现解析解可以较好地预测梯度柱壳链中载荷的演化趋势.正梯度链中载荷峰值随着波传播逐渐增大,负梯度链中载荷峰值随着波传播逐渐减小.正梯度链支撑端峰值载荷高于均匀链,负梯度链支撑端峰值载荷低于均匀链,这表明相较于均匀柱壳链,密度梯度柱壳链可以在更大范围内对波形进行调控.线性密度梯度参数对梯度柱壳链的波形调控能力影响较大,梯度参数越小,传递到支撑端的峰值载荷越小;相反,梯度参数越大,支撑端的峰值载荷越大.建立的理论模型及其解析解为研究梯度柱壳链中应力波传播规律及揭示载荷调控机理提供了理论基础.     

近场动力学系统中波传播特性的探究

近场动力学是一类基于非局部思想的新固体力学方法,其采用积分形式的控制方程,自然地适用于极端载荷下材料破碎和裂纹发展的模拟,被广泛用于国防安全等领域的研究.但是,非局部性会引入色散效应,对波的传播产生不利影响,制约其对断裂等固体行为的捕捉能力.为此,采用谱分析方法,对近场动力学系统的色散行为进行了全面的研究.发现相比于低频成分,高频成分的色散关系呈现出振荡趋势和零能模式,色散问题更为严重.高频域的色散行为还随波的传播方向发生改变,呈现出沿45°的对称性.而近场动力学系统本身缺乏数值耗散,无法抑制色散问题带来的不利影响.因此,从引入数值耗散的角度出发,在合理保留传统近场动力学理论框架的基础上,建立了黏性引入的控制方程.并考虑固体中常见的体积变形和对高频成分的选择性抑制,构造了相应的黏性力态.最后,在数值研究中模拟了极端载荷下激波的产生,以探究波的间断性对色散行为的影响.发现间断性强的波表现出更为显著的色散行为,呈现出Gibbs不稳定性.这些均能有效地被黏性力态所抑制,验证了所提方法的正确性.这为在近场动力学系统中实现对波传播过程的正确捕捉,获得正确的固体行为提供了重要参考,从而为国防安全领...     

考虑阻尼的无限长小垂度缆索弹性波的传播

缆索结构被广泛应用于电气、土木、海洋和航空工程等领域,随着缆索在工程中的应用长度越来越长,高阶振动越来越明显,研究时应该考虑扰动沿着缆索的传播.现有对缆索弹性波传播的研究中,通常不考虑阻尼项,然而阻尼对于波的传播有着重要影响.文章考虑阻尼的影响,发展了包含阻尼项的三维弹性缆索运动方程.通过求解上述含阻尼项的运动方程,分别考察了面内面外弹性波的频率关系、相速度和群速度等自由传播特性,进而通过计算无限长弹性缆索在初始余弦型脉冲作用下的位移响应,分析扰动沿着该缆索的传播规律,考察波的色散现象以及阻尼对于缆索弹性波传播的影响.结果表明,考虑阻尼后,面内波和面外波均为色散波,面内波在曲率的作用下,为高度色散波.此外,在阻尼的影响下,波的峰值在传播过程不断减小,且波的后缘端点响应总是高于前缘端点响应.     

非均匀损伤介质中波传播的数值解

对弹性波在非均匀损伤介质中的传播理论进行了研究。通过将非均匀损伤区域离散成分层均匀的区域,结合相邻区域交界面处的连续条件,推导出了以右行波、左行波为状态向量的波动方程和传递矩阵。对几种非均匀损伤介质中波的传播进行了实例数值计算,并和其解析解的结果进行了比较,讨论了弹性波在非均匀损伤介质中传播的一般性质。     

基于Lord-Shulman理论的热弹性波在偶应力三明治固体中的反射和透射

在温度急剧变化、短时间极速加热等极端情况下,基于Fourier定律的热流矢量与温度梯度成正比关系的经典热传导理论不能准确描述其物理过程．经典热弹性理论的热传导方程是抛物型的,而广义热弹性理论包含双曲型方程,热将以具有有限传播速度的波动形式传播．本文基于Lord-Shulman广义热弹性理论和修正偶应力弹性理论,得到在偶应力热弹性固体中四种色散波,研究热弹性波的传播和在偶应力固体三明治结构中的反射透射问题,重点研究横波入射时偶应力参数和热弛豫时间对各种热弹性耦合波反射透射系数的影响．     

粘弹性损伤混凝土介质中波的传播

混凝土的损伤会引起弹性和粘性性质的改变,利用弹性和粘性的双参数来描述损伤,建立双参数损伤理论。根据不同的损伤程度对混凝土介质进行分区处理,建立基本方程式,考虑连续性条件和边界条件,求解波动方程。比较是否考虑粘性时的波幅和波传播时间与损伤的关系,表明粘性是损伤混凝土介质中波传播所必须要考虑的问题。分析了粘弹性混凝土介质中损伤区域长度、损伤度等对波传播的影响,给出了它们的关系曲线,可为波的反分析提供依据。     

一维声子晶体带隙的非局部辛分析

周期性弹性复合结构(声子晶体)中传播的弹性波存在特殊的色散关系:弹性波只能在某段频率范围内无损耗的传播,该频率范围称为通带.一维声子晶体的色散问题可以看作分层介质中弹性波的传播问题,利用二维弹性理论予以分析.为了研究非局部效应对声子晶体带隙特性的影响,将Eringen的二维非局部弹性理论引入到Hamilton体系下,利用精细积分与扩展的Wittrick Williams算法可获取任意频率范围内的本征解.通过对不同算例的数值计算,分析和对比了非局部理论方法与传统局部理论方法的差别.并进一步指出了该套算法的适用性和优势所在.     

关于热波理论的研究

当热传导服从经典的Fourier定律时,温度场由抛物型方程所控制,热扰动以无限大速度传播。在通常情况下,热波的迅速衰减掩盖了这种佯谬。但对于热爆炸,热核聚变,快速化学反应,强激光与物质相互作用这样一些时间尺度极短(与从非平衡态达到局部平衡态的时间相比)的情形,需对Fourier定律进行修正。本文从Cattaneo方程及其唯象修正,Boltzmann演化方程,分子动力学数值模拟以及连续介质热力学理论四个方面对热波理论近20年的进展进行了评述,并对热波的实验验证及数值分析的某些重要结果进行了介绍,且提出了关于热波非弹性理论的新认识,讨论了可能的应用前景。     

管内流体对单壁碳纳米管中弯曲波频散的影响

由于其管状结构,碳纳米管在纳机械系统中可望被用作输流管道.采用连续介质力学方法,研究管内有流体存在时碳纳米管中弯曲波的传播和频散.建立流体存在时考虑二阶应变梯度的非局部弹性Timoshenko梁方程.流体的存在,使相速度最低的解支相速度降低.当流速较低时,流速对碳纳米管中弯曲波传播的影响不大.当流速较高时,相速度最低的一支随流速增加相速度降低.当流速非常高时,该解支会消失.但流体的存在对其他解支影响不大.随着波数的升高,非局部弹性所描述的微结构对碳纳米管中弯曲波传播的影响越来越明显.     

功能梯度材料板中Lamb波传播特性研究

对材料性能参数沿厚度连续变化的横观各向同性热应力缓和型功能梯度材料板中Lamb波的传播问题,采用幂级数法,求得其相速度方程.借助数值算例,分析了参数梯度变化对Lamb波频散曲线的影响,并与相应陶瓷板和金属板中的频散曲线进行了对比.进一步研究了参数梯度变化对波结构的影响,揭示了Lamb波在这种非均质板中的传播行为,所得结果可以为功能梯度材料及结构的超声表征与检测提供理论依据.     

Strain localisation patterns under equal-channel angular pressing

Instabilities of plastic flow in the form of localised shear bands were experimentally observed to result from equal-channel angular pressing (ECAP) of magnesium alloy AZ31. The appearance of shear bands and their spacing were dependent on velocity of the pressing and applied back-pressure. A generic gradient plasticity theory involving second-order strain gradient terms in a constitutive model was applied to the case of AZ31 deformed by ECAP. Linear stability analysis was applied to the set of equations describing the deformation behaviour in the process zone idealised as a planar shear zone. A full analytical solution providing a dispersion relation between the rate of growth of a perturbation and the wave number was obtained. It was shown that the pattern of incipient localised shear bands exhibits a spectrum of characteristic lengths corresponding to admissible wave numbers. The interval of the spectrum of wave numbers of viable, i.e. growing, perturbations predicted by linear stability analysis was shown to be in good agreement with the experimentally observed spectrum. The effect of back-pressure applied during ECAP was also considered. The predicted displacement of the shear band spectrum towards lower wave numbers, shown to be a result of the decreased shear strain rate in the shear zone, was consistent with the experimentally observed increase of the band spacing with increased back-pressure. A good predictive capability of the general modelling frame used in conjunction with linear stability analysis was thus demonstrated in the instance of the particular alloy system and the specific processing conditions considered.     

Wave dispersion in gradient elastic solids and structures: A unified treatment

Analytical wave propagation studies in gradient elastic solids and structures are presented. These solids and structures involve an infinite space, a simple axial bar, a Bernoulli–Euler flexural beam and a Kirchhoff flexural plate. In all cases wave dispersion is observed as a result of introducing microstructural effects into the classical elastic material behavior through a simple gradient elasticity theory involving both micro-elastic and micro-inertia characteristics. It is observed that the micro-elastic characteristics are not enough for resulting in realistic dispersion curves and that the micro-inertia characteristics are needed in addition for that purpose for all the cases of solids and structures considered here. It is further observed that there exist similarities between the shear and rotary inertia corrections in the governing equations of motion for bars, beams and plates and the additions of micro-elastic (gradient elastic) and micro-inertia terms in the classical elastic material behavior in order to have wave dispersion in the above structures.     

Characteristics of elastic waves in FGM spherical shells,an analytical solution

In this paper, we investigate the propagation characteristics of elastic guided waves in FGM spherical shells with exponentially graded material in the radial direction. A new separation of variables technique to displacements is proposed to convert the governing equations of the wave motion to the second-order ordinary differential equations with variable coefficients. By further a variables transform technique, these equations are transformed to the Whittaker equations so that analytic solutions can be obtained. For the spherical shell case, by satisfying the traction-free boundary conditions on both the inner and outer surfaces of the shell, we obtain the dispersion equations, which show that both the SH and Lamb-type wave modes are generated in the structure. The calculated dispersion curves in the functionally graded shell demonstrate a clear influence of the gradient coefficient as compared to those of the homogeneous shell, with the Lamb-type waves more sensitive to the gradient coefficient. The mode shapes and the distributions of stresses in the shells for various gradient coefficients are also presented to illustrate their dependence on the gradient coefficient.     

Explicit and implicit gradient series in damage mechanics

Gradient enhancement series are studied in the context of damage mechanics. Distinction is made between so-called explicit series and implicit series, both of which can be derived from a nonlocal damage model. The paper focuses on the difference between second-order and fourth-order truncations for either series. Dispersion analysis and numerical simulations are used to compare the various models. It is shown that for the explicit series the fourth-order term has a detrimental influence on the response, while for the implicit series the fourth-order term leads to a slightly closer approximation of the nonlocal model. The role of the critical wave length as it emerges from the dispersion analysis is shown to be decisive. When the critical wave length acts as an upper bound, a stable response is obtained and the critical wave length equals the width of the damaging zone. On the other hand, when the critical wave length acts as a lower bound, oscillations may appear of which the periodicity is set by this critical wave length.     

A higher-order nonlocal elasticity and strain gradient theory and its applications in wave propagation

In recent years there have been many papers that considered the effects of material length scales in the study of mechanics of solids at micro- and/or nano-scales. There are a number of approaches and, among them, one set of papers deals with Eringen's differential nonlocal model and another deals with the strain gradient theories. The modified couple stress theory, which also accounts for a material length scale, is a form of a strain gradient theory. The large body of literature that has come into existence in the last several years has created significant confusion among researchers about the length scales that these various theories contain. The present paper has the objective of establishing the fact that the length scales present in nonlocal elasticity and strain gradient theory describe two entirely different physical characteristics of materials and structures at nanoscale. By using two principle kernel functions, the paper further presents a theory with application examples which relates the classical nonlocal elasticity and strain gradient theory and it results in a higher-order nonlocal strain gradient theory. In this theory, a higher-order nonlocal strain gradient elasticity system which considers higher-order stress gradients and strain gradient nonlocality is proposed. It is based on the nonlocal effects of the strain field and first gradient strain field. This theory intends to generalize the classical nonlocal elasticity theory by introducing a higher-order strain tensor with nonlocality into the stored energy function. The theory is distinctive because the classical nonlocal stress theory does not include nonlocality of higher-order stresses while the common strain gradient theory only considers local higher-order strain gradients without nonlocal effects in a global sense. By establishing the constitutive relation within the thermodynamic framework, the governing equations of equilibrium and all boundary conditions are derived via the variational approach. Two additional kinds of parameters, the higher-order nonlocal parameters and the nonlocal gradient length coefficients are introduced to account for the size-dependent characteristics of nonlocal gradient materials at nanoscale. To illustrate its application values, the theory is applied for wave propagation in a nonlocal strain gradient system and the new dispersion relations derived are presented through examples for wave propagating in Euler–Bernoulli and Timoshenko nanobeams. The numerical results based on the new nonlocal strain gradient theory reveal some new findings with respect to lattice dynamics and wave propagation experiment that could not be matched by both the classical nonlocal stress model and the contemporary strain gradient theory. Thus, this higher-order nonlocal strain gradient model provides an explanation to some observations in the classical and nonlocal stress theories as well as the strain gradient theory in these aspects.     

Nonlocal damage model using the phase field method: Theory and applications

A new nonlocal, gradient based damage model is proposed for isotropic elastic damage using the phase field method in order to show the evolution of damage in brittle materials. The general framework of the phase field model (PFM) is discussed and the order parameter is related to the damage variable in continuum damage mechanics (CDM). The time dependent Ginzburg–Landau equation which is also termed the Allen–Cahn equation is used to describe the damage evolution process. Specific length scale which addresses the interface region in which the process of changing undamaged solid to fully damaged material (microcracks) occurs is defined in order to capture the effect of the damaged localization zone. A new implicit damage variable is proposed through the phase field theory. Details of the different aspects and regularization capabilities are illustrated by means of numerical examples and the validity and usefulness of the phase field modeling approach is demonstrated.     

Implicit and explicit secular equations for Rayleigh waves in two-dimensional anisotropic media

This paper is concerned with the derivation of implicit and explicit secular equations for Rayleigh waves polarized in a plane of symmetry of an anisotropic linear elastic medium. It has been confirmed, in accord with Ting’s paper [2], that the Rayleigh waves propagate with no geometric dispersion. Numerical evaluations of both the implicit and explicit equations give the same values of Rayleigh wave velocities. In the case of orthotropic material (thin composites) it has been found that Rayleigh wave velocity depends significantly (as with bulk waves) on the directions of principal material axes. For the same material model the analytical solutions, based on implicit and explicit secular equations, were compared against the finite element and experimental data that had been published by Cerv et al. [4] in 2010. It emerged that the theory was in accordance with the experiment.     

Analysis of plate in second strain gradient elasticity

The bending analysis of a thin rectangular plate is carried out in the framework of the second gradient elasticity. In contrast to the classical plate theory, the gradient elasticity can capture the size effects by introducing internal length. In second gradient elasticity model, two internal lengths are present, and the potential energy function is assumed to be quadratic function in terms of strain, first- and second-order gradient strain. Second gradient theory captures the size effects of a structure with high strain gradients more effectively rather than first strain gradient elasticity. Adopting the Kirchhoff’s theory of plate, the plane stress dimension reduction is applied to the stress field, and the governing equation and possible boundary conditions are derived in a variational approach. The governing partial differential equation can be simplified to the first gradient or classical elasticity by setting first or both internal lengths equal to zero, respectively. The clamped and simply supported boundary conditions are derived from the variational equations. As an example, static, stability and free vibration analyses of a simply supported rectangular plate are presented analytically.     

PROGRESSIVE FRACTURE MODELING OF THE FAILURE WAVE IN IMPACTED GLASS

The failure wave has been observed propagating in glass under impact loading since 1991. It is a continuous fracture zone which may be associated with the damage accumulation process during the propagation of shock waves. A progressive fracture model was proposed to describe the failure wave formation and propagation in shocked glass considering its heterogeneous meso-structures. The original and nucleated microcracks will expand along the pores and other defects with concomitant dilation when shock loading is below the Hugoniot Elastic Limit. The governing equation of the failure wave is characterized by inelastic bulk strain with material damage and fracture. And the inelastic bulk strain consists of dilatant strain from nucleation and expansion of microcracks and condensed strain from the collapse of the original pores. Numerical simulation of the free surface velocity was performed and found in good agreement with planar impact experiments on K9 glass at China Academy of Engineering Physics. And the longitudinal, lateral and shear stress histories upon the arrival of the failure wave were predicted, which present the diminished shear strength and lost spall strength in the failed layer.