• 主题
高级搜索  >
历史搜索 清空历史
首页> 中文学位 >弹性材料及超弹性材料中波传播问题的间断Galerkin有限元方法
【6h】

弹性材料及超弹性材料中波传播问题的间断Galerkin有限元方法

代理获取
页面导航
  • 摘要
  • 著录项
  • 相似文献
  • 相关主题

摘要

本文主要针对材料学中三种材料:棒状非线性弹性复合材料、细柱状非凸弹性材料以及杆状可压缩超弹性材料中波在传播过程中追赶碰撞问题以及波的相变问题进行研究。利用（局部）间断Galerkin方法及自适应Runge-Kutta间断Galerkin(RKDG)方法进行数值模拟。本文内容大致分成四个部分。
　　第一部分针对棒状非线性弹性复合材料，其模型方程可以看成是守恒律方程组，但是由于不同材料的（两层）复合，模型方程中流通量项出现间断，且在模型方程的右端源项出现δ函数，这些都给数值格式的设计增加了困难。由于流通量项的不连续性，因此在处理数值流通量的时候需要注意。另外由于物理界面的存在，波发生反射和透射，因此方程的解会包含多个激波，从而也增加了问题的难度。我们首先从物理的角度，引入耗散率和物理容许解的概念，证明间断Galerkin(DG)格式的能量稳定性和DG解是物理容许解。然后借助于耗散率，我们提出一种dissipation-rate reserving DG格式，并通过数值算例比较说明该格式的优点。最后我们利用dissipation-rate reserving DG格式数值模拟了复合材料中波在传播过程中发生碰撞的细节，包括拉伸激波、压缩激波和稀疏波之间的碰撞，并给出了详细的碰撞时间和位置;通过选取四组不同参数比较碰撞之后拉伸波的幅值减小率来说明参数的影响，结果也验证了材料的非线性性能够使得材料中具有破坏性的拉伸波降低幅值达400％。这也表明这样的非线性弹性材料具有很好的抗冲击能力。最后我们考虑应力参数ε以及右边界对波的影响，还有四层复合材料中材料排列方式对材料保护的影响。
　　第二部分针对细柱状非凸弹性材料，由于应力-应变关系的非凸性，将会导致解的不适定性，因此在材料中考虑柱体径向变形的影响得到新的模型方程，正是这种影响使得新的模型方程中包含高阶时间导数项，从而增加求解的难度。为了便于使用局部间断Galerkin(LDG)方法进行求解，我们引入新的时间辅助变量从而将模型方程转化为只具有一阶的时间导数的方程组，针对这一方程组设计了LDG方法。其次证明了LDG格式的稳定性，并通过数值算例给出LDG格式的精度:在L1模和L∞模意义下应变γ具有k+1阶精度。最后利用LDG格式进行了一系列数值模拟，由于色散和材料非凸性的影响，我们观察到了一些有趣的波动模式（wave pattern），例如，transformation front与类似于孤立子的波组成的波动模式和稀疏波与类似于孤立子的波组成的波动模式等等;另外我们给出了稀疏波与transformation fronts碰撞的详细过程。
　　第三部分针对杆状可压缩超弹性材料，其模型方程（可压缩超弹性杆波动方程）与Camassa-Holm方程形式很一致，但是在非线性项多一个常数系数γ。正是由于常数系数γ的存在，使得模型方程不再具有好的性质（例如可积性）。模型方程所具有的多种类型波（例如紧孤立波，孤立子波波）是我们重点考察对象，我们利用LDG格式数值模拟这些类型波之间的相互碰撞，发现一些有趣的现象:紧孤立波、孤立子波在Gauss扰动下仍然是稳定的;两个相向的紧孤立波之间碰撞后以及两个右行孤立子波之间碰撞后均有类似于cuspon波形出现;两个紧孤立波碰撞后会出现两个孤立子波;三个紧孤立波碰撞后会出现三个大小不同的孤立子波;一个右行紧孤立波与一个孤立子波碰撞后会出现两个大小不同的孤立子波。
　　第四部分同样针对棒状非线性弹性复合材料波的传播问题，从第一部分的数值计算结果可以看出应变γ包含多个激波，因此采用较大数目的网格进行模拟来抓住激波的信息，但是网格数目增大同时也大大增加了计算量。因为DG方法具有易于实现自适应的优点，我们采用自适应RKDG方法并结合KXRCF指示子(indicator)和TVB指示子这两种不同的指示子进行模拟，对是否需要TVB重构、参数M、有无耗散DG格式以及Pk元（所有次数不超过k的多项式）等多方面的影响进行考虑。总体而言，TVB指示子的计算结果要比KXRCF指示子好，因为存在的振荡少些。但是激波处出现少许光滑化的现象，因此还需通过选取合适的M值以使得结果更为精确。

著录项

相似文献

  • 中文文献
  • 外文文献
  • 专利
代理获取

客服邮箱：kefu@zhangqiaokeyan.com

京公网安备：11010802029741号 ICP备案号：京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有

  • 客服微信

  • 服务号