收藏本站
收藏 | 手机打开
二维码
手机客户端打开本文

内域波动数值模拟的显式方法

刘恒  
【摘要】: 如何更精确高效地模拟大型、复杂系统内域的波动是发展和完善近场波动数值模拟技术的一个重要研究课题。内域波动的数值模拟通常采用计算量较小的显式方法,但现有的时空解耦显式有限元方法的精度只有二阶;低精度不仅影响数值模拟的精度,而且制约着计算效率的提高。鉴于此,本文旨在探索内域波动数值模拟具有更高精度且稳定的显式方法。作者发展了内域波动数值模拟的现有显式有限元解耦技术,提出了一种高精度且稳定的显式数值模拟方法。 第一,依据波速有限的概念,从波动方程的精确解出发提出了一种新的显式数值模拟方法。此方法与现有有限元技术相似之处在于:适于处理非规则网格节点,且节点递推公式是具有显式和时空解耦特征。其不同之处在于:新方法可给出时空离散精度皆为2M阶的稳定格式,M为正整数。本文通过一维模型详细论述了这一方法的可行性:建立了非规则网格节点递推公式,详细分析了均匀网格标量波动数值模拟的精度和稳定性,提出了构建时空精度皆为2M阶( M为正整数)的稳定递推公式的技术途径,并以构建二阶( M = 1)和四阶( M = 2)公式为例予以说明。 第二,将一维情形研究结果推广到高维情形,分别建立了二维、三维非规则网格节点的递推公式。针对二维正方形均匀网格详细论述了时空离散精度皆为2M阶的稳定递推公式的构建方法,给出了二阶和四阶递推公式的具体结果;并对三维立方体均匀网格模型作了简要的讨论。 第三,本文通过一维、二维模型算例检验了波动方程数值模拟显式方法的精度和稳定性等理论结果,指出了高阶公式对提高计算效率的价值。 第四,从空间解耦有限元常微分方程组(结构动力学方程)出发,通过被积函数的拉格朗日多项式内插和分部积分导出了一组具有高阶精度的时域显式积分格式,并将其推广为适用于一般数学物理方程的显式积分格式。作者以一个简单的线性时不变系统为例,初步考察了此积分格式的稳定性。


知网文化
【相似文献】
中国重要会议论文全文数据库 前10条
1 陆春阳;苏庆田;黄生富;;有限元计算中混凝土徐变应力分析方法比较[A];第18届全国结构工程学术会议论文集第Ⅰ册[C];2009年
2 蒋跃文;张伟伟;叶正寅;;一种基于CFD技术的高效气动弹性时域求解方法[A];第十届全国空气弹性学术交流会会议论文集[C];2007年
3 田稳苓;邢超;赵欣;李宁;;无加劲钢板剪力墙结构滞回性能数值分析[A];庆祝刘锡良教授八十华诞暨第八届全国现代结构工程学术研讨会论文集[C];2008年
4 刘恒;廖振鹏;刘晶波;;三维波动方程的高精度显式格式[A];第八届全国地震工程学术会议论文集(Ⅰ)[C];2010年
5 范尚武;;波动方程快速三维深度偏移方法[A];1993年中国地球物理学会第九届学术年会论文集[C];1993年
6 李志成;卢红标;钱七虎;;抗爆结构计算中外推法应用研究[A];第五届全国结构工程学术会议论文集(第一卷)[C];1996年
7 吴静萍;郑伟涛;杨鹏;高玉玲;秦江涛;;自由面流动黏性数值模拟关键技术研究[A];第二十届全国水动力学研讨会文集[C];2007年
8 刘祥庆;刘晶波;丁桦;;八种结构动力反应积分算法的对比分析[A];第14届全国结构工程学术会议论文集(第一册)[C];2005年
9 叶正寅;王刚;;非结构混合网格下的复杂外形粘性流动计算研究[A];计算流体力学研究进展——第十二届全国计算流体力学会议论文集[C];2004年
10 黄乘明;;泛函微分与泛函方程的数值稳定性[A];第九届全国微分方程数值方法暨第六届全国仿真算法学术会议论文集[C];2004年
中国博士学位论文全文数据库 前5条
1 刘恒;内域波动数值模拟的显式方法[D];中国地震局工程力学研究所;2009年
2 张晓志;近场SH波动数值模拟若干问题研究[D];中国地震局工程力学研究所;2002年
3 王鹏;随机常微分方程数值分析中的若干方法[D];吉林大学;2008年
4 马克颖;发展方程的并行GALERKIN区域分解方法[D];山东大学;2009年
5 张荣培;一类非线性扩散方程的间断有限元方法研究[D];中国工程物理研究院;2012年
中国硕士学位论文全文数据库 前10条
1 冯光;一类分数阶微分方程初值问题的数值方法[D];湘潭大学;2008年
2 王波兰;二维非结构网格Euler方程的LU-SGS算法[D];南京航空航天大学;2003年
3 王德才;拟动力试验数值积分方法及子结构拟动力试验研究[D];合肥工业大学;2007年
4 李赳华;管材液压胀形与板材液压拉深数值模拟方法的研究[D];华中科技大学;2007年
5 常元章;基于质点—弹簧模型的布料仿真[D];天津大学;2008年
6 金永虎;三角指数拟合Runge-Kutta方法[D];华中科技大学;2011年
7 王晓波;碟形弹簧的力学性能研究[D];郑州大学;2007年
8 詹涌强;改进的Hamming预测—校正法[D];暨南大学;2007年
9 李晓红;常微分方程数值解法及其应用[D];东北师范大学;2008年
10 武继银;基于OSG的实时布料仿真[D];山东大学;2009年
中国知网广告投放
 快捷付款方式  订购知网充值卡  订购热线  帮助中心
  • 400-819-9993
  • 010-62982499
  • 010-62783978