二维波动方程的动态有限元离散化的频散研究

本文利用二维波动方程，推导出波在单元内行进的动态形函数，在此基础上，得到了波在离散体内传播的频散关系，利用这一关系式，绘出了频谱曲线．分析表明：波在离散体内传播的频散程度不仅与波的传播方向有关而且与单元的质量分布、网格划分有密切的关系．     

有限元畸变单元离散模型中的SH波动

本文采用晶格动力学分析方法,研究了有限元畸变单元离散模型中的SH波的传播规律。结果表明单元畸变对波动规律影响显著,当单元畸变较小时,单元畸变对波动传播并不产生明显的不利影响,当单元畸变较大时,这一影响不可忽略。通过详细分析和对比集中质量、混合质量和一致质量模型中的波动规律,评价了不同有限元模型的抗畸变能力。文中还讨论了时域离散化与单元畸变联合作用对波动传播的影响。     

脆性材料动态破坏过程的数值模拟

根据连续介质力学基本原理和离散元法基本思想，针对脆性材料动态破坏过程问题，推导出基于圆盘单元的正交各向异性离散元二维计算模型，并利用此模型计算了各向异性薄板在冲击载荷下的应力波传播问题。通过将计算结果与LS－DYNA程序相比较，表明了本文提出的离散元模型的精确性。另外，对钢弹侵彻下滑凝土圆板破坏过程这个典型的发生连续介质向非连续介质转化的动力学问题进行了数值模拟和动画显示，并同侵彻实验中出现的几种破坏形式进行了分析和对比，从而证明了本计算方法是计算和模拟脆性材料出现从连续介质向非连续介质转变的动态破坏问题的有力工具。     

粘塑性靶板中冲击波的演化

讨论了粘塑性靶板中的一维应变波的传播规律。利用广义特征理论导出了应力波传播的特征线和特征关系。利用特征关系和冲击波阵面上的突跃条件，得出了冲击波在传播过程中的演化的规律，并以Ｂｏｄｎｅｒ－Ｐａｒｔｏｍ幂函数型粘塑性材料为例，计算和讨论了板中应力波传播规律的特点     

粘塑性靶板中冲击波的演化

讨论了粘塑性靶板中的一维应变波的传播规律。利用广义特征理论导出了应力波传播的特征线和特征关系。利用特征关系和冲击波阵面上的突跃条件，得出了冲击波在传播过程中的演化的规律，并以Ｂｏｄｎｅｒ－Ｐａｒｔｏｍ幂函数型粘塑性材料为例，计算和讨论了板中应力波传播规律的特点。     

滞后细观模型在岩石力学中的应用

对以砂岩为代表的所谓``NME材料'的力学行为研究方面的一些新的概念和模型进行了评介.首先介绍了一种基于所谓``滞后单元'的描述滞后现象的物理模型------Preisach-Mayergoyz(P-M)模型,然后详细阐述了P-M模型应用于模拟岩石的非线性滞后应力应变关系的过程和结果.这种唯象模型很好地描述了宏观上的滞后表现和``离散记忆'效应.接着本文对应变能耗散的力学机制进行了简单分析. 最后,介绍了一种描述弹性波在``NME材料'中传播规律的数学方法, 该方法从一般的弹性波传播规律出发,分析了``NME材料'特殊的力学性质给弹性波传播带来的影响,揭示了产生特殊的弹性波传播规律的原因.      

用二次形函数薄层法分析弹性层状地基中的动力问题

薄层法是分析和模拟弹性波在层状介质中传播的一种半解析半数值方法。本文在土层垂直方向离散中利用Galerkin加权残值法推导出二次形函数薄层元的计算公式。采用薄层单元模拟半空间上的层状场地,模型底面用阻尼器边界或傍轴边界代替半空间。利用点源简谐荷载作用下的土层反应与其它数值分析方法的对比讨论薄层模型的设置指标。并用一次和二次两种形函数离散方法计算了层状地基中的面波频散曲线、圆形均布简谐荷载作用下弹性半空间的位移反应和半无限地基中单桩的竖向阻抗函数,分别讨论其计算精度。     

Lamb波在复合材料板中传播的谱有限元模拟

众所周知Lamb波在复合材料中的传播呈各向异性的特点,经典有限元法模拟这类问题效率不高,所以,本文采用谱有限元法进行研究。先建立了一种新的谱有限板单元,该单元以Gauss-Lobatto-Legendre点作为节点,使质量矩阵是对角矩阵;另外,该单元采用了扩展的位移场,能够较好地模拟板结构的三维特性。然后,对复合材料板结构中Lamb波在对称模式与反对称模式下的传播速度进行了求解,将计算结果与Mindlin板谱单元的结果以及三维弹性理论的结果进行了比较,并讨论了Lamb波在反对称层合板中的传播特点。最后,模拟了Lamb波在含和不含损伤复合材料层合板中的传播,数值结果表明所建立的谱有限板单元可以较好地模拟出Lamb波在复合材料板结构中的传播特性。     

研究材料动态本构特性中的重要作用

在材料动态本构关系的研究中，不论是由波传播信息反求材料本构关系，即所谓解第二类反问题，还是利用应力波效应和应变率效应解耦的方法（如SHPB技术），应力波传播实际上都起着关键作用。在一般性讨论的基础上，就SHPB试验技术分析了应力波传播如何影响材料动态本构特性的有效确定。对于应力/应变沿试件长度均匀分布假定以及一维应力波假定，着重进行了分析。     

细观结构固体层中孤立波的传播及相互作用特性

利用直接微扰方法.确定了孤立波的放大或衰减与孤立波的初始幅度以及介质的结构参数之间的关系.然后利用线性化技术构造出一种二阶精度的稳定差分格式,并对孤立波在细观结构固体层中传播特性进行了数值模拟,特别对细观结构固体层中传播的不同幅度的孤立波的相互作用进行了详细的数值模拟,从而得到在适当条件下细观结构固体层中孤立波传播时即可以衰减、放大又可以稳定传播,且相互作用不影响这种传播特性.     

弹性半空间二维出平面自由波场的一维化时域算法

提出了一种计算出平面SH波斜入射时弹性半空间自由波场时域计算的一维化有限元方法。首先利用Snell定律确定平面波沿水平方向的传播规律,在用有限元法对弹性半空间进行离散化时,竖向单元尺寸根据波动有限元模拟精度要求确定,而水平向有限元网格尺寸根据水平向波的传播规律和采用的离散时间步长确定,使得有限元离散模型中任意节点的运动可以用水平向相邻节点的运动表示,从而将二维有限元节点运动方程组化为一维的形式。求解此一维方程组,可得到弹性半空间中一列节点的运动,再根据行波的传播规律,可确定全空间自由波场。理论分析和数值算例表明,该方法具有较高的精度和良好的稳定性。     

原发性肺冲击伤生物力学研究概述（二）

原发性肺冲击伤生物力学研究概述（二）赵敏，王正国（重庆第三军医大学野战外科研究所重庆６３００４２）（续第１４卷第１期第９６页）４激波作用时体内应力波的传播生物体是个非常复杂的结构。激波作用下体内应力波的传播一直是研究者们非常关心的问题。但应力波在体内     

离散元法在求解三维冲击动力学问题中的应用

提出了三维连结型离散模型,建立了可实现连结型模型(用于连续介质)-接触型模型(用于非连续介质)转化的三维离散元计算程序,用来模拟连续介质转变为非连续介质的力学过程．利用该计算程序对冲击载荷下混凝土块体内(连续体情况下)的应力波传播过程进行了数值模拟．将计算结果的数值与LS-DYNA程序计算的结果进行比较,验证了该计算程序的计算精度．在此基础上,模拟了混凝土块体的动态破坏(连续介质向非连续介质转化)过程．其计算结果可用动画显示,得到的破坏形式与由实验得到的破坏形式相近．两个算例说明该离散元模型及其计算程序是模拟计算伴随有连续介质向非连续介质转变的动态破坏问题的有力工具．     

三维离散网格中的弹性波动

本文采用形式解法研究三维有限元离散模型中波动的传播规律。文中由波动方程得到它的特征方程，进而给出了三维离散网格中波的频散关系，分析了三维网格中平面波、非匀匀波、寄生振荡及其它可能存在的运动形式的基本特点；详细讨论了平面波的频散、截止频率及偏振漂移等重要的不同于连续介质中的波动传播性质的特点，提示了其中不同于低维网格中 波动的现象。     

多洞室对矢量波散射引起半空间表面位移的边界元解

采用边界元方法研究了半空间中近表面多洞室对矢量波的散射问题,给出了以全空间格林函数为基本解且半空间表面离散的边界积分方程,在这一边界各分方程中,较好的消除了主值积分,在半空间表面进行离散时,采用无限单元与有限单元相结合的方法,大大减少了计算量,提出了精度。     

爆轰波在扩张管内临界传播特性的研究

由于稀疏波的作用，爆轰波在变截面扩张管内传播时会发生衰减。临界条件下，爆轰波恰能在扩张管内维持自持传播的特性，称为爆轰波的临界传播特性。针对Ｃ＿２Ｈ＿２－Ｏ＿２－Ａｒ不同的爆轰气体成份，在不同扩张半角的方形扩张管内，作者研究了临界初压ｐ＿（ｃｒ）、胞格尺寸Ｓ与扩张半角的关系，用Ｗｈｉｔｈａｍ理论与重点火准则分析了爆轰波在扩张管的传播，得到胞格尺寸的计算值与实验值基本相符。     

非均质变截面杆中弹性波的反射和透射

基于一维弹性波理论，本文对应力波在非均质变截面杆中传播问题进行了一维简单波分析，并把分析结果与二维轴对称有限元分析结果进行了比较，表明一维简单波分析是非常有效和实用的。利用一维简单波分析方法，本文还揭示了应力波在非均质变截面杆中的传播规律，特别对含有内部交界面的非均质变截面杆（带有连接段）进行了一维等效简化分析，研究了连接段对应力波传播的影响。     

激波与小扰动波的相互作用(激波稳定性问题)

本文讨论朗道和栗弗席兹意义下的激波稳定性问题,也就是求解激波与小扰动波的相互作用问题。在无耗散介质的定常、平面激波两边产生了小扰动。朗道等假定小扰动波传播方向是垂直于激波的,得到的稳定条件是M_>1,M_2<1。我们则把小扰动波取为二维的,即推广为小扰动波可以斜着激波传播的情形,这是更常见的。结论是,无论按群速度方向来区分趋向波和离散波,还是像朗道等那样,按相速度方向来区分趋向波和离散波,既便满足条件M_1>1,M_2<1,对于某些纵向扰动波长与频率范围的声波而言,激波都不稳定。接着提出了几个实验,根据可能的实验结果讨论了趋向波与离散波的定义应按相速度方向还是群速度方向,说明了其意义。 对于管道内的定常激波有类似结果。     

阶梯压电层合梁的波动动力学特性

采用行波理论系统地研究了压电阶梯梁的自由振动分析以及强迫响应的分析方法. 基于分布 参数理论研究了压电阶梯梁的波传播特性，忽略柔性梁横向剪切和转动惯量的影响，给出了 梁的轴向和横向的简谐波解. 将压电阶梯梁离散化为单元，考虑压电片的刚度和质量的影响， 建立了节点散射模型. 应用位移连续和力平衡条件，推导了节点的波反射和波传递矩阵，在 此基础上，引入波循环矩阵的概念，给出波循环矩阵、波传递系数矩阵的确定方法. 应用波 循环矩阵可以有效地计算结构的固有频率. 另外，应用波传递系数研究了压电陶瓷作动器位 置对其驱动能力的影响. 得出两个主要结论：1)作动器靠近悬臂梁固定端将有较强的驱动 能力，悬臂梁边界反射行波产生弯曲消失波有利于增大压电波的模态传递系数；2)模态传 递系数与固有频率的灵敏度密切相关，波传递系数越大, 对应该处固有频率变化灵敏度越大. 另外，数值算例表明了行波方法比有限元方法具有更高的计算精度.     

非均匀水流中非线性波传播的数值模拟

以一种考虑波流相互作用的新型{Boussinesq}型方程为控制方程组， 采用五阶{Runge}-{Kutta}-{England}格式离散时间积分，采用七点 差分格式离散空间导数，并通过采用恰当的出流边界条件，从而建立了非均匀水流中非线性 波传播的数值模拟模型. 通过对均匀水流与水深水域内和潜堤地形上存在弱流或强流时波浪 传播的数值模拟，说明模型能有效地反映水流对波浪传播的影响.