弹性半无限域问题边界元法研究

论文摘要

边界元法作为一种数值方法,因其计算量小、精度高而得到了发展。边界元法能够自动满足无穷远处边界条件,不用进行人工截断,无需施加人工边界,处理无限域问题具备先天优势。有限域和无限域的边界元法已经取得较为完善的研究成果,为处理工程中众多半无限域问题,本课题基于现有弹性静力学及动力学边界元法的理论及计算程序,通过引进无限边界单元,处理角点问题,编制相应的Matlab程序处理半无限域的静力学及动力学波动问题,并计算交通荷载下边坡的动力响应。本文进行了弹性静力学半无限域边界元法的理论推导。针对半无限域问题,提出将其分为有限域部分和无限域部分。有限域部分基于已有的边界元理论,着重研究了角点问题,提出边界边坡角点属于混合角点,应分别对其进行几何角点和物理角点的处理。无限域中引进无限边界单元,推导出坐标插值函数以及相应的位移衰减函数。因Kelvin基本解含有刚体位移,对其进行修正,并嵌套进边界元理论中进行计算,选取算例进行验证。弹性静力学的方法为动力学提供理论基础和支撑。在静力学之后,进行了弹性动力学半无限域边界元法的理论推导。基于现在较成熟的时域边界元法,采用与静力学相同的解决思路,将空间离散分为常规单元和无限边界单元。对三维基本解的第三坐标进行积分得到二维时域基本解,经过一系列的积分运算得到平面半无限域弹性动力学积分方程。进行数值计算,处理奇异元素的奇异性,对积分方程进行求解。求解交通荷载下边坡的动力响应。选取爆炸双指数和突加荷载对半无限域动力时域边界元法的理论和程序进行验证。从算例结果可以看出精度和计算效率都较可靠,最后计算交通荷载下边坡的动力响应。

论文目录

• 摘要
• Abstract
• 第1章 绪论
•   1.1 课题背景
•   1.2 课题研究的目的和意义
•   1.3 边界元法的发展
•   1.4 时域边界元的历史与现状
•   1.5 主要研究内容
•   1.6 技术路线
• 第2章 弹性静力学半无限域问题边界元法
•   2.1 引言
•   2.2 弹性静力学半无限域边界积分方程的建立
•     2.2.1 弹性静力学的基本解
•     2.2.2 性静力学半无限域边界积分方程的建立
•   2.3 边界积分方程的数值计算
•     2.3.1 边界积分方程的离散
•     2.3.2 有限域的H、G矩阵元素的求解方法
•   2.4 角点问题的处理
•   2.5 半无限域边界积分方程的数值计算
•     2.5.1 无限边界单元插值函数的构造
•     2.5.2 无限边界单元的理论准备
•     2.5.3 半无限域边界积分方程的离散
•     2.5.4 弹性静力学半无限域H′、G′矩阵元素求解方法
•   2.6 算例验证
•     2.6.1 算例介绍
•     2.6.2 算例解析解
•     2.6.3 计算结果和分析
•   2.7 本章小结
• 第3章 动力学半无限域问题时域边界元法
•   3.1 引言
•   3.2 弹性动力学边界积分方程
•     3.2.1 弹性动力学基本方程
•     3.2.2 弹性动力学问题的基本解
•     3.2.3 边界积分方程的建立
•     3.2.4 半无限域问题的边界积分方程
•   3.3 边界积分方程的数值处理
•     3.3.1 边界积分方程的离散
•     3.3.2 边界积分方程的组装
•     3.3.3 处理非奇异性积分与奇异性积分
•   3.4 弹性动力学边界积分方程求解
•     3.4.1 H、G矩阵的形成
•     3.4.2 计算H、G矩阵中非对角元素
•     3.4.3 计算H、G矩阵中对角元素
•   3.5 本章小结
• 第4章 交通荷载下边坡的动力响应
•   4.1 引言
•   4.2 半无限域验证
•     4.2.1 算例描述
•     4.2.2 算法流程
•     4.2.3 结果和分析
•   4.3 计算交通荷载下边坡动力响应
•     4.3.1 交通荷载介绍
•     4.3.2 计算结果
•   4.4 本章小结
• 结论
• 参考文献
• 致谢

文章来源

类型: 硕士论文

作者: 许先林

导师: 雷卫东

关键词: 时域边界元法,弹性动力学,半无限域,无限边界单元

来源: 哈尔滨工业大学

年度: 2019

分类: 基础科学

专业: 数学,数学

单位: 哈尔滨工业大学

分类号: O241.8

DOI: 10.27061/d.cnki.ghgdu.2019.002547

总页数: 82

文件大小: 3744K

下载量: 67