高度非线性孤立波与弹性大板的耦合作用研究

一维颗粒链的一端受到一个有初速度颗粒的撞击,导致颗粒连中产生稳定传播的应力波——高度非线性孤立波,该应力波的波长、波速以及幅值都能保持很好的稳定性,且遇到边界才会反射. 孤立波是一种良好的信息载体,广泛应用于无损检测技术中. 基于孤立波的特性,研究高度非线性孤立波与弹性大板耦合作用,基于赫兹定律和板的内在非弹性理论,推导出晶体链与大板的耦合微分方程组. 用龙格库塔法求解该微分方程组,得到颗粒链中各颗粒的位移、速度曲线. 通过分析回弹波出现的时间、回弹波所携带的能量以及模量、厚度、重力等对孤立波的影响,发现反射孤立波对大板的弹性模量和厚度尤为敏感,此外,颗粒链的摆放对整个耦合过程也有影响. 研究的结果为孤立波对结构体的无损探伤提供了理论依据,该技术可实现对结构体的快速检查和可控性研究.     

高度非线性孤立波与Euler梁的耦合作用研究

无预压缩一维颗粒链中的机械撞击以高度非线性孤立波的形式传播,该孤立波是一种能量包,在较长距离内保持稳定波幅。本文研究高度非线性孤立波与Euler梁的耦合作用,对颗粒链采用离散元模型,基于Hertz定律和Euler梁振动方程并考虑Euler梁的横向振动,建立颗粒链与Euler梁耦合的微分方程组。采用Runge-Kutta法进行求解,得到晶体链中各颗粒的位移、速度解。分析回弹波出现的时间、波幅以及Euler梁的材料特性和几何特性参数等对反射孤立波的影响,结果表明,反射孤立波对Euler梁的弹性模量和几何尺寸在一定范围内是敏感的。本文的研究成果丰富和完善了孤立波探伤理论,为一维颗粒链的实际应用提供了广阔的应用前景。     

含电学边界的压电层合梁的非线性弯曲波

非线性科学己成为近代科学发展的一个重要标志, 特别是非线性动力学和非线性波的研究对于解决自然科学各领域中遇到的复杂现象和问题有着极其重要的意义. 本文研究了含电学边界条件的压电层合梁的非线性弯曲波传播特性.首先, 考虑几何非线性效应和压电耦合效应, 利用哈密顿原理建立了一维无限长矩形压电层合梁弯曲波的非线性方程.其次, 采用Jacobi椭圆函数展开法对非线性弯曲波方程进行求解, 得到了非线性弯曲波动方程在近似情况下对应的冲击波解和孤波解.最后, 利用约化摄动法得到了非线性薛定谔方程, 进一步得到了亮孤子和暗孤子解.基于两种方法具体研究了外加电压、压电层厚度等参数对冲击波和孤立波以及亮孤子和暗孤子特性的影响. 研究结果表明, 在波速较小时, 外加电压对冲击波的影响较大, 波速较大时, 外加电压对孤立波影响减弱.通过调整作用在压电层合梁上的电压发现了存在亮孤子和暗孤子, 分析结果表明随着外加电压值的增大, 亮孤子和暗孤子的振幅都增大.      

含界面相Voronoi单元的等效弹性常数的计算

采用含界面相Voronoi单元有限元法,根据广义胡克定律,计算了在给定边界条件下,颗粒增强复合材料的等效弹性常数。建立了含多个随机分布的椭圆形夹杂及界面相的VCFEM模型,分析了夹杂体分比,界面相厚度和界面相弹性模量等因素对颗粒增强复合材料等效弹性常数的影响,并利用普通有限元方法对比验证。结果表明,当界面相弹性模量小于基体与夹杂时,材料的等效弹性模量会随着界面相厚度的增大而减小,随着夹杂体分比的增大而减小,并且界面过薄时,材料的等效弹性模量会随着夹杂体分比的增大而增大;当界面相弹性模量大于基体或夹杂时,材料的等效弹性模量会随着夹杂体分比和界面相厚度的增大而增大。而界面相的厚度和弹性模量对材料的等效泊松比的影响较小,材料的等效泊松比主要受夹杂体分比的影响,与其呈反比关系。     

爆炸波作用下地下防护结构与围岩的非线性动力相互作用分析

针对实际地下工程中普遍存在的材料非线性以及半无限介质域的处理问题,给出了基于时间有关基本解的时域边界元法与非线性动力有限元法的耦合方法,应用该耦合方法计算了一马蹄形截面地下防护结构与围岩受爆炸冲击波作用下非线性相互作用的时间历程,并与线弹性情况进行了比较分析。结果表明:本文的方法具有较高精度,真实地再现了波在弹性层中传播以及反射的全过程。     

非稳态纯滚动接触弹塑性分析

建立了二维弹塑性非稳态循环纯滚动接触有限元模型.材料本构采用一种较好的循环塑性模型,并通过材料用户子程序在通用有限元软件ABAQUS中自定义该本构模型.通过在弹塑性无限半空间表面上重复移动随时间按简谐规律变化的赫兹法向载荷来模拟非稳态循环纯滚动接触过程.通过数值模拟,得到接触表面附近的残余累积变形、应变和残余应力.不同的最大赫兹接触压力对残余应力和残余应变影响较大.在简谐变化的法向接触载荷作用下接触表面的变形呈波浪形,随着滚动次数的增加,该波状表面沿载荷移动相反方向逐渐移动,但移动速率要衰减.波状表面波谷处的残余应力、应变和变形大于波峰处.随滚动次数的增加,残余应力增大但很快趋于稳定,残余应变也增大但增大速率衰减.     

球链中冲击扰动的传播

采用Hertz接触理论处理牛顿摆中球与球之间的相互接触,建立球链碰撞模型,观察到碰撞过程中类似孤立波的传播现象.为了便于处理"n对m"碰撞的计算,编写了带有GUI(graphic user interface)界面的计算程序,可以研究不同撞击球个数对撞击过程的影响,包括球心位置、球心速度以及相邻小球间的变形量与相互作用力随时间的变化规律.进一步比较了一维球链模型和三维有限元模型的结果,并实验观察了碰撞过程的细节,得到一致性结论.     

微结构固体中钟型与扭结孤立波的演化

根据Mindlin理论和Murnaghan模型,首先建立了描述耗散、频散及非线性微结构固体中一维纵波传播的一种简单模型.然后利用有限差分方法,数值模拟了微结构效应对钟型与扭结孤立波演化的影响. 结果表明,随着微结构效应的减弱,钟型孤立波的幅度衰减以及非对称特征变得越来越明显;随着微结构效应的增强,扭结孤立波顶部出现的“帽子”状变化以及由此产生的非对称特征变得越来越明显.      

基于FEM-DEM的三体摩擦界面中接触行为与应力的多尺度分析

为了研究三体摩擦界面中第一体变形与第三体状态的相互影响,利用耦合有限元法和离散元法的多尺度法模拟了平行板剪切颗粒第三体的过程。整个模型分为两个区域：有限元区域（上板）和离散元区域（第三体和下板）,上板在一定的外载荷压应力下挤压颗粒第三体,下板以恒定的速度剪切颗粒第三体。为实现两个子区域间的相互联系,建立了子区域间应力应变的传递机制。实现了三体摩擦界面的多尺度分析,模拟了平行板剪切颗粒第三体的过程。模拟结果表明：当外载荷压应力低于10MPa时,颗粒间的碰撞增多使得第三体内量纲归一化平均应力增大,宏观摩擦系数也随之增大;在剪切过程中,第三体内部颗粒间的接触随接触角度的分布呈现出一定的规律性,0~90°各区间内的强接触较多,尤其在54°~72°之间;颗粒接触随接触力大小的分布也具有一定的规律性,接触力与第三体颗粒平均接触力的比值在0~0.6之间内的接触较多,随后接触力越大,接触数越少。同时,第一体内压应力分布与第三体内力链的分布相对应,力链越强则与其接触的第一体的压应力越大。     

PROPAGATION OF IMPACT PULSE IN A CHAIN OF ELASTIC BALLS

采用Hertz 接触理论处理牛顿摆中球与球之间的相互接触,建立球链碰撞模型,观察到碰撞过程中类似孤立波的传播现象。为了便于处理“n 对m”碰撞的计算,编写了带有GUI (graphic user interface) 界面的计算程序,可以研究不同撞击球个数对撞击过程的影响,包括球心位置、球心速度以及相邻小球间的变形量与相互作用力随时间的变化规律。进一步比较了一维球链模型和三维有限元模型的结果,并实验观察了碰撞过程的细节,得到一致性结论。     

High-amplitude elastic solitary wave propagation in 1-D granular chains with preconditioned beads: Experiments and theoretical analysis

Elastic solitary waves resulting from Hertzian contact in one-dimensional (1-D) granular chains have demonstrated promising properties for wave tailoring such as amplitude-dependent wave speed and acoustic band gap zones. However, as load increases, plasticity or other material nonlinearities significantly affect the contact behavior between particles and hence alter the elastic solitary wave formation. This restricts the possible exploitation of solitary wave properties to relatively low load levels (up to a few hundred Newtons). In this work, a method, which we term preconditioning, based on contact pre-yielding is implemented to increase the contact force elastic limit of metallic beads in contact and consequently enhance the ability of 1-D granular chains to sustain high-amplitude elastic solitary waves. Theoretical analyses of single particle deformation and of wave propagation in a 1-D chain under different preconditioning levels are presented, while a complementary experimental setup was developed to demonstrate such behavior in practice. The experimental results show that 1-D granular chains with preconditioned beads can sustain high amplitude (up to several kN peak force) solitary waves. The solitary wave speed is affected by both the wave amplitude and the preconditioning level, while the wave spatial wavelength is still close to 5 times the preconditioned bead size. Comparison between the theoretical and experimental results shows that the current theory can capture the effect of preconditioning level on the solitary wave speed.     

Directional Wave Propagation in a Highly Nonlinear Square Packing of Spheres

We studied the dynamic response of a two-dimensional square packing of uncompressed stainless steel spheres excited by impulsive loadings. We developed a new experimental measurement technique, employing miniature tri-axial accelerometers, to determine the stress wave properties in the array resulting from both an in-plane and out-of-plane impact. Results from our numerical simulations, based on a discrete particle model, were in good agreement with the experimental results. We observed that the impulsive excitations were resolved into solitary waves traveling only through initially excited chains. The observed solitary waves were determined to have similar (Hertzian) properties to the extensively studied solitary waves supported by an uncompressed, uniform, one-dimensional chain of spheres. The highly directional response of this system could be used as a basis to design granular crystals with predetermined wave propagation paths capable of mitigating stress wave energy.     

On the coupling mechanism between nonlinear solitary waves and slender beams

This paper describes the coupling mechanism between highly nonlinear solitary waves propagating along a granular system and slender beams in contact with the granular medium. Nonlinear solitary waves are compact non-dispersive waves that can form and travel in nonlinear systems such as one-dimensional chains of particles, where they are conventionally generated by the mechanical impact of a striker. These waves have a constant spatial wavelength and their speed, amplitude, and duration can be tuned by modifying the particles’ material or size, or the velocity of the striker. In the study presented in this article we investigated numerically the interaction between solitary waves propagating along a chain of granular particles and a slender beam. Some of the numerical findings were validated experimentally. We found that the geometric and mechanical properties of the beam or thermal stress applied to the beam alter certain features of the solitary waves. In the future, these findings may be used to develop a novel sensing system for the nondestructive assessment of beams.     

Solitary wave interactions and reshaping in coupled systems

The interaction of the components of composite solitary waves governed by nonlinear coupled equations is studied numerically. It is shown how predictions of the known exact traveling wave solutions may help in understanding and explaining the process of reshaping seen as head-on and take-over collisions of individual solitary waves. The most interesting results concern the switch in the sign or the periodic modulation of the amplitude of the solitary wave and the direction of its propagation due to collisions.     

Propagation of solitary waves over underwater ridges

The propagation of solitary waves is investigated on the basis of a nonlinear system of equations of hyperbolic type describing the motion of the crest of a solitary wave over the surface of a liquid of variable depth [1]. The existence of solutions with discontinuities, the boundary conditions at which are introduced on the basis of [2, 3], is assumed. In the case of an infinite cylindrical ridge both solitary and periodic captured waves are found. Depending upon the height of the ridge and the parameters of the wave, the encounter between a uniform wave and a semi-infinite ridge yields qualitatively different solutions — continuous and discontinuous, where the primary wave is broken down by the ridge into several solitary waves. The amplitude of the wave may either increase or decrease over the ridge.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 36–93, January–February, 1985.The author is grateful to A. G. Kulikovskii and A. A. Barmin for their interest in his work, useful discussions and valuable comments offered during the preparation of the article for the press.     

Formation of solitary waves on gas-sheared liquid layers

The origin of solitary waves on gas-liquid sheared layers is studied by comparing the behavior of the wave field at sufficiently low liquid Reynolds number, R_L, where solitary waves are observed to form, to measurements at higher R_L where solitary waves do not occur. Observations of the wave field with high-speed video imaging suggest that solitary waves, which appear as a secondary transition of the stratified gas-liquid interface, emanate from existing dominant waves, but that not all dominant waves are transformed. From measurements of interface tracings it is found that for low R_L, waves which have amplitude/substrate depth (a/h) ratios of 0.5–1 occur while for higher R_L, no such waves are observed. A comparison of amplitude/wavelength ratios shows no distinction for different R_L. Consequently, it is conjectured that solitary waves originate from waves with sufficiently large a/h ratios; this change of form being similar to wave breaking. The dimensionless wavenumber is found to be smaller at low R_L, where solitary waves are observed. This suggests that perhaps, larger precursor (to solitary wave) waves are possible because the degree of dispersion, which acts to break waves into separate modes, is lower.     

Frontal collision of internal solitary waves of first mode

The dynamics and energetics of a frontal collision of internal solitary waves (ISW) of first mode in a fluid with two homogeneous layers separated by a thin interfacial layer are studied numerically within the framework of the Navier–Stokes equations for stratified fluid. It was shown that the head-on collision of internal solitary waves of small and moderate amplitude results in a small phase shift and in the generation of dispersive wave train travelling behind the transmitted solitary wave. The phase shift grows as amplitudes of the interacting waves increase. The maximum run-up amplitude during the wave collision reaches a value larger than the sum of the amplitudes of the incident solitary waves. The excess of the maximum run-up amplitude over the sum of the amplitudes of the colliding waves grows with the increasing amplitude of interacting waves of small and moderate amplitudes whereas it decreases for colliding waves of large amplitude. Unlike the waves of small and moderate amplitudes collision of ISWs of large amplitude was accompanied by shear instability and the formation of Kelvin–Helmholtz (KH) vortices in the interface layer, however, subsequently waves again become stable. The loss of energy due to the KH instability does not exceed 5%–6%. An interaction of large amplitude ISW with even small amplitude ISW can trigger instability of larger wave and development of KH billows in larger wave. When smaller wave amplitude increases the wave interaction was accompanied by KH instability of both waves.     

梯度介质半空间Rayleigh面波传播性能研究

对剪切弹性模量沿深度以指数函数变化的非均质半空间，本文用摄动法得到了Rayleigh面波的波函数解答及相速度方程。以不同金属与陶瓷复合而成的几种梯度材料为例，用数值方法求解了相速度方程，给出了相应的波的弥散曲线，结果表明，梯度介质半空间自由表面附近的Rayleigh波通常有两种不同的弥散形式，即正常弥散和非正常弥散。     

相变复合波在伪弹性TiNi合金薄壁管中的传播

相变可以改变材料的性质,从而严重影响波在介质中传播的结构。采用考虑静水压力和偏应力联合作用的增量型相变本构模型,研究了在拉(压)-扭联合作用下半无限长TiNi合金薄壁管内相变复合波的传播规律。基于广义特征理论分析了相变复合波的特征波速及简单波解的基本性质。利用数值方法研究了两种典型情况下管内相变耦合波传播的规律,管内传播的应力路径和波的结构与初始状态及加载幅值有关,展现出和普通弹塑性材料截然不同的性质。