改进时域间断Galerkin有限元方法在弹塑性波传播数值模拟中的应用.pdfVIP

第34卷第3期 固体力学学报 V01．34NO．3 OFSOLIDMECHANICS CHINESEJOURNAL 2013 2013年6月 June 改进时域间断Galerkin有限元方法 在弹塑性波传播数值模拟中的应用+ 郭 攀1 武文华p。 吴志刚2 (1大连理工大学运载工程与力学学部力学系，工业装备结构分析国家重点实验室，大连，116024) (2大连理工大学运载工程与力学学部航空航天学院，工业装备结构分析国家重点实验室，大连，116024) 摘要对于高频、强脉动荷载作用下的结构动力学波传播分析，对比于传统的时域算法，时域间断Galerkin 方法能捕捉到波阵面的间断，有效地避免了由于间断引起的数值振荡．但时域间断方法却带来了波前面的虚假数 值振荡．论文针对上述波前数值振荡的现象进行研究，通过引入人工阻尼的方法对时域间断Galerkin有限元方法 进行进一步改进．数值结果表明，所发展的方法能够有效的滤掉强动荷载产生的波前数值振荡现象，同时降低了时 域间断Galerkin方法的网格依赖性． 关键词 改进时域间断Galerkin方法(mDGFEM)，人工阻尼，波传播，强动荷载 0简介 刻和终点时刻的节点基本未知向量．这样的处理 以Newmark方法等为代表性的时域有限元法 过程使得算法计算效率较低．李锡夔等凹’1叩在传统 通过在空间域采用有限元离散，时域采用有限差分 DGFEM方法的基础上做了进一步的发展．仅使位 法，在低频响应为主的结构动力学问题分析中通常 移基本向量的时间导数向量存在间断，节点位移则 能得到满意的数值求解结果[1]．但对于强动荷载作 在离散的时间段之间自动保证连续．具体的在时域 用下，具有双曲型特征的结构动力方程进行波传播 对位移进行三次Hermite插值，对位移导数进行时 的数值分析时，上述算法常呈现出虚假的数值振荡 间的线性插值．这样大大提高了计算效率．但在无阻 现象．尤其的，不能有效捕捉在波阵面上的解的间断 尼的强动荷载作用下波传播问题的数值模拟中，虽 K 行为．为此，TammaK等[21提出采用人工阻尼的然波后相比Newmark方法较好的滤去数值振荡的 方法消减强动荷载作用下动力响应中称为第二声的 情况下，波前的计算结果却往往会出现虚假的数值 数值振荡现象．国内方秦等通过隐式和显式New— 振荡现象．同时，随着网格的加粗，波前的数值振荡 mark数值阻尼的技巧较好的消除了这种数值振荡进一步加剧． 本文在李锡夔等所发展DGFEM的基础上，提 的现象[3“]．众多研究表明[511]，间断Galerkin有限 出通过引入刚度比例人工阻尼的新的DGFEM．数 元法(DGFEM)具备自动引入数值耗散和滤去虚假 的高阶模式和数值振荡效应的能力，进而得到了广 值结果显示出所提出的改进DGFEM方法有效的 泛的注意和研究[5．6]．DGFEM的重要特征是在空间滤去了这种波前数值振荡的现象，同时降低了 域和时间域同时采用有限元离散，对问题的半离散 DGFEM对网格较强的依赖性． (空间离散)控制方程中节点基本未知数向量及其时 1 含人工阻尼的改进时域间断伽辽金有限元方法 间导数向量在时间域中独立分片插值，并允许它们 简介 在离散的时间段之间存在间断，其间断值通过变分 双曲型动力