弹性动力学问题的常应力组合杂交有限元方法.pdfVIP

第 35卷 第 2期 西 华 师 范 大 学 学 报 (自然 科 学 版 ) 20l4年 6月 Vo1．35 No．2 JournalofChinaWestNormalUniversity(NaturalSciences) Jun．2014 文章编号 ：1673-5072(2014)02-0103-05 弹性动力学问题的常应力组合杂交有限元方法 余正琴 ，陈豫眉 (1．四川大学数学学院，四川 成都 610064；2．西华师范大学数学与信息学院。四川 南充 637009) 摘 要 ：本文将组合杂交有限元法应用于求解弹性动力学问题．位移选取标准的双线性元，应力采用分片常数 时间方向上采用中心差分格式．数值算例表明，组合杂交方法具有较好的数值精度． 关键词：平面弹性动力学问题；双线性元；组合杂交有限元 中图分类号：0242．21 文献标识码：A 1 引 言 ( ( ( E 关于平面弹性振动问题的求解方法讨论已经有很多，例如，Joly和Tsogka把二阶系统 (1)改写成c一；阶 ∈ ∈ ∈ 系统 ¨ ，并为新系统构造新的四边形和立方体混合有限元．Boulaajine，Farhloul和 Paquet在求解线弹性动 Q C；r × × × 力学问题时采用 了一种新 的对偶混合有 限元方法 ，并且在时间离散上采取 了显式和隐[ 式[ 差分[ 格式． 0 0 0 Douglas和 Gupta用混合有限元方法讨论了在平面区域中波动方程的超收敛性 ，其中，在得到Sobolev空间 明 中非正指数的误差估计时使用 了拟投影分析方法．为了解决平面弹性振动问题 ，Lai，Huang和 Chen建立 了 一 个新 C。一连续时间步位移型有限元方法 ．显式的差分格式在适当的CFL条件下稳定，而隐式的差分格 式是无条件稳定的．成和谢在处理振动问题的时候采取了时空非协调有限元方法 ，在空间上采取了非协 调不完全二次矩形元和 c—R元两种方法 ，在时间上采用 了连续的分片二次多项式插值． 本文将组合杂交有限元方法 用于如下平面弹性振动问题的求解 ： M 一divq= ，t) or=￡8()=2 (u)+AdivuI = 0 u(x，0)= o()，u(，0)= () ／L 其 中，nci表示有界开 区域 ，边界记为 F，T为正常数 ，=(u。， )为位移向量 ，ro=(ro)：：为对称应力张 量，8(Ⅱ)=(Vu+(Vu))／2为应变张量，，A为 Lam6常数 ，对于平面应变问题取 A= ’ = ≠_l，对于平面应力问题取A= ，=≠ ，这里。。．5为泊松比，E为杨氏模 量 ∈(L( ×[0，T]))表示体力密度，咖。()，。()∈( (Q))分别为初始位移和初始速度． 2 组合杂交有限元格式 定义应力和位移两个空间如下 ： ∑：=L2( ；i2一x2)，：=( ( )) 收稿 日期 ：2014—