矩形板中应力波主波和次波的传播与反射

弹性应力波的修正理论指出关于体积应变的波动方程与现有理论一致,但发展了一组关于体积应变和偏应变的弱耦合波动方程。针对矩形板受侧向集中载荷冲击下应力波的波动问题,建立了应力波传播的两组控制方程以及加载面和自由面的波动边界条件。采用有限差分方法求解波动方程,数值分析了应力波关于主波和次波的传播以及自由面上斜入射波的反射过程。偏应变在传播过程中分裂为两部分,一部分与体积应变共同传播组成主波,另一部分传播较慢形成次波。数值模拟结果显示与冲击载荷下纳钙玻璃板中应力波的传播图像是完全符合的。     

双相介质中储层波场特征数值模拟研究

根据双相介质弹性波方程,并对弹性波方程进行分离得到了膨胀波方程。利用高阶差分对膨胀波方程进行数值模拟,边界处理采用完全匹配层的吸收边界条件,算法实现了任意高阶差分和自动加载边界条件的纵波方程正演模拟。计算结果表明,在双相各向同性介质中存在快纵波、慢纵波,两种纵波有明显的区别,快纵波的速度远大于慢纵波的速度。在分界面上,快纵波要产生透射快纵波、透射转换慢纵波,反射快纵波、反射转换慢纵波。慢纵波具有很强的衰减性,耗散系数越小,慢纵波越明显,反之亦然。该算法的高阶差分形式可以显著地降低数值频散,有效提高地震波正演计算的精度,适应性好,操作简便灵活。     

源造波法在三维完全非线性数值波浪水槽中的应用

为消除波浪在水槽出流边界的反射及波浪遇到结构物后在入射边界形成的二次反射,基于时域高阶边界元方法建立三维完全非线性数值波浪水槽模型:利用源造波法产生入射波浪,采用人工阻尼层技术消除波浪反射及二次反射;建立水槽的格林函数并应用于整个计算域,消除水槽两侧壁和底面的积分,减少了计算量.利用所建模型分别对规则波和不规则波的数值实验表明,数值结果与理论解吻合良好,且在出流边界和入射边界均无反射现象.     

岩石在磨料射流作用下破坏机理

描述了岩石材料在磨料射流的冲击下破坏的主要原因，分析了应力波就是岩石局部受到拉伸和剪切井最终使岩石发生破坏的根本因素，在这一假设的前提下采用波动理论建立了岩石受应力波作用的力学模型即射流波动的基本方程，在已知的边界条件下推导出了由波动引起的最大正应力的计算公式。通过该公式对具体实例进行了理论计算，并把该计算结果与在河南油田所进行的地面磨料射流割缝的现场实验结果进行了比较，证明两种结果吻合的较好．说明这种理论分析方法是可行的。     

含垂直穿透裂纹板拉伸断裂的有限元模拟

对含垂直穿透裂纹板受拉伸载荷作用的断裂过程进行了数值模拟.选用增量型弹塑性本构关系,采用有限元方法求解虚功原理方程.首先对弹塑性杆的动态响应进行了数值计算.参照LS_DYNA中处理断裂问题的单元失效方法,对受拉伸载荷作用的垂直穿透裂纹板断裂过程进行了计算.结果表明:一维杆的弹性波和塑性波的计算结果和解析解相符,说明了该方法和计算程序计算结构动态响应的可靠性.受冲击拉伸载荷作用,应力波在板中传播并在裂纹表面反射后出现畸变,裂纹尖点出现应力集中,裂纹沿初始方向扩展,当裂纹扩展长度为初始长度2倍时,裂纹在板边界处出现分叉,裂纹扩展过程和理论分析结果相符.     

应力波界面效应试验

利用动光弹法模拟研究在冲击作用下界面上的应力波传播过程,利用应力波理论分析了界面上应力波的作用机理,证实了界面上存在着能量累积和阻滞能量传递的效应。并结合弹性波动理论和应力—光学定律导出了波动方程与条纹级数之间的关系式。     

基于反周期扩展边界方法的单程波正演模拟

单程波动方程的正演模拟相对于基于双程波方程的正演模拟,有很多的优势。但是基于波动方程波场延拓中大都会存在边界反射吸收问题,常规的衰减型边界吸收条件在单程波正演模拟中难以达到较好的效果。而本文采用了反周期扩展边界条件可以很好的解决这个问题。经过模型试算对比可以说明,该方法具有更好的效果。     

旋转长杆弹侵彻混凝土靶应力波传播特性研究

利用动量定理和动量矩定理,建立了旋转长杆弹垂直侵彻混凝土时的运动微分方程,并对此微分方程进行了数值求解,讨论了撞击力和侵彻位移的时程特性。在此基础上,利用波动理论建立了纵波的波动方程和扭转波的波动方程。利用有限差分法对波动方程进行了数值求解,讨论了弹体内部应力波的传播特性,重点研究了初始旋转角速度对长杆弹内应力波传播特性的影响。     

特征线坐标系下柱对称超音速流膨胀波计算方法研究

在静压圆盘止推气体轴承柱对称超音速流膨胀波系的计算中,借鉴Liepmann的特征线坐标系与自然坐标系变换理论和Liepmann引入两个黎曼不变量的方法,使得边界条件的处理更为直观和简便。采用特征线坐标系下的特征线法,推导计算该膨胀波系所用到的基本计算公式,将其应用于P-M流动并与P-M流动理论解进行比较。结果表明,数值计算结果与P-M流动理论解吻合,该计算方法可行。     

耦合横向和微旋转波的界面反射与折射

基于广义微极磁热弹模型讨论了两微极弹性固体介质平面界面上耦合横向和微旋转波的反射与折射问题.利用两微极弹性固体平面分界面上位移和应力连续条件,建立了反射与折射波振幅的线性系统方程.数值计算给出各类反射、折射波振幅随耦合横向和微旋转波入射角的变化关系曲线.研究了热松弛时间对反射和折射热波、纵向位移波、耦合横向和微旋转波的振幅比的影响.结果表明,热松弛时间对热波、纵向位移波影响显著,而对耦合横向和微旋转波影响很小.     

非均匀介质中应力波反射诱发层裂过程的数值模拟

采用基于细观损伤力学基础上开发的动态版RFPA2。数值模拟软件，对不同冲击栽荷作用下非均匀介质中应力波反射谤发层裂过程进行数值模拟。通过对不同应力波延续时间、不同应力波峰值和不同应力波波形下应力波诱发层裂过程的数值分析，对应力波反射诱发层裂的规律和机制进行探讨。研究结果表明：应力波延续时间和应力幅值对岩石的层裂过程影响较大，不同形状的应力波实际上是在相同时段内应力波具有不同的应力幅值和不同的应力加栽速率，不同形状的应力波对介质的作用结果是这些因素同时作用的综合反映;应力波反射谤发层裂过程的发生，是各个时刻不同加载条件下微破裂累积的宏观体现。     

波浪在变深区的向岸传播

基于非线性长波方程,除在开边界对幅射波进行了线化处理外,还对局部波数(波长)、局部相位以及入射波速进行了修正。以圆柱绕射、斜坡底折射和波浪向岸传播为实例进行了计算。结果表明:对于波向变化、破波位置及波幅变化诸方面其计算与经验是一致的。因此本文的工作拓宽了非线性长波方程的适用范围,而且为应用Boussinesq方程研究浅水波变形打下了基础。     

面向海洋工程应用的数值波浪水池

采用FLUENT软件中的用户自定义函数功能和动网格功能进行造波和消波,采用VOF方法捕捉自由液面,开发了一个数值波浪水池,该水池可以生成规则波浪、按照波谱生成不规则波浪以及生成既定的不规则波波列.文中开展了网格和时间步长试验,力图找到经济精确的网格划分方案和高效的时间步长方案.文中按照Jonswap谱进行了不规则波浪的数值模拟,采用了等分能量法分解波谱的频率区间来得到组成波的圆频率,模拟结果显示计算波谱与目标波谱一致,波浪的统计特征得到了很好的体现,模拟不规则波波谱的缺点是耗时过长.由于真实海况下对海洋浮式结构物造成损坏的往往不是整个波浪时间序列,而是少数的几个大波或者波列中的波群.因此采用CFD方法准确模拟出既定的波列将会对了解极端海况下海洋结构物的损坏事故提供极大帮助.采用本文提出的方法成功模拟出了既定的极大波,这个方法的关键在于确定造波板的运动规律.类似按照波谱造波,文中将造波板处的波浪时历进行了分解,得到了若干个小振幅波,通过傅里叶级数法推导出了每个小振幅波的特征参数.同样采用水力传递函数确定每个小振幅波对应的推板振幅,之后将这些组成波对应的推板运动规律进行叠加得到了模拟不规则波波列时需要的推板运动规...     

爆炸应力波作用下的凿岩爆破破岩模拟

爆破破岩是爆炸应力波和爆生气体载荷共同作用的结果,本文采用改进的DDA方法模拟爆炸应力波作用下的爆破破岩过程. 采用Hopkinson层裂算例验证了DDA方法对应力波传播和动态断裂破坏问题模拟的准确性. 通过在炮孔壁上施加三角冲击载荷,模拟再现了凿岩爆破中爆炸应力波作用下的孔周破坏、自由表面剥落、碎裂漏斗区域形成、较远区裂纹的扩展,以及岩体的抛掷,结合应力波传播模拟结果对破岩过程与机理进行了分析. 研究为进一步采用DDA模拟两种形式爆破载荷共同作用下的凿岩爆破创造了条件.     

水波传播的无网格数值模拟

在水波传播的数值模拟中采用了一种基于配点和径向基函数的无网格方法.采用Laplace方程的基本解作为径向基函数,将源点布置在模拟波浪场之外,沿边界布置配点而不是划分网格,从而自动满足控制方程,且不存在奇点,不需要求解积分方程.数值造波采用给定入射波面和速度势的方法,数值消波综合采用阻尼层消波和Sommerfeld辐射条件,非线性自由面的演化追踪采用二阶Taylor级数展开式.对线性规则波和非线性三阶Stokes波的模拟显示,数值结果与理论解吻合良好.表明无网格方法不但形式简单、计算速度快,而且稳定性和准确性令人满意,有望成为水波模拟问题的一种有效的数值方法.     

岩石的损伤软化对应力波传播的影响

针对不含水干岩弹塑性变形的基本特征，考虑岩石的剪胀效应和软化效应，并根据岩石不同变形阶段的特点，给出了不同的屈服面形式和损伤演化方程，建立了岩石的弹塑性连续损伤本构模型．在此基础上研究了一维应变波在岩石介质中的传播规律，揭示了应力波在岩石类脆性损伤软化材料的传播过程中，损伤软化使得材料分为损伤演化区和完全破坏区，应力平台在两区交界面处产生凹陷．这一不稳定的凹陷在传播过程中逐渐消失，因而应力波传播在损伤的作用下经历了一个从不稳定到稳定的演化过程．     

弹性方柱中波的传播规律Ⅲ. 无限长方柱的瞬态响应

提出行波追踪方法,在以前工作中得到解析形式的频散方程、群速度方程和稳态响应的基础上,完全得到无限长方柱在任意外力作用下解析形式的瞬态响应。绘制了无限长方柱在外力作用下的瞬态响应,并同有限长方柱中的行波进行了比较。     

完全匹配层吸收边界条件的研究

边界反射是数值模拟中的一个关键问题。构建完全匹配层(PML)吸收边界条件是消除边界反射的理想方法之一。此方法是在研究区域的边界上加入吸收层,使边界上传入吸收层的波随传播距离按指数规律衰减,不产生任何反射,达到消除边界反射的目的。首先给出了二维声波方程完全匹配层的公式,然后通过地面地震模型和井间地震模型数值算例与常规吸收边界条件进行了比较。结果表明,PML边界条件具有吸收更干净、相位对应关系好且能够吸收各种角度的边界反射等优点,从而验证了完全匹配层吸收效果的优越性。