论文题目: 快速多极边界元法研究及其在复合材料模拟中的应用
论文类型: 博士论文
论文专业: 力学
作者: 王海涛
导师: 姚振汉
关键词: 边界元,快速多极,大规模,复合材料
文献来源: 清华大学
发表年度: 2005
论文摘要: 边界元法是伴随着有限元法发展起来的一种有效的数值方法。边界离散的优点使边界元法很适合模拟具有复杂边界或者界面的结构,如颗粒增强或纤维增强复合材料。然而,边界元法形成的系数矩阵通常是非对称的满阵,常规求解技术效率低下,使得边界元法不能有效处理大规模问题。为了克服上述困难,快速多极算法逐步用于加速边界元法的数值求解。在固体力学领域,由于基本解的复杂性,寻找基本解合适的展开和传递格式仍然是目前快速多极边界元法的研究重点。另外,新型算法的研究工作还很有限。本文提出了一种适用于二维弹性力学问题的初始快速多极边界元格式,将边界元法的存储和计算复杂度都降到O(N)的量级,从而极大提高了计算效率。数值算例表明,在二维问题上,快速多极边界元法在保证精度的前提下,在计算量和存储量上都比常规求解技术有数量级上的减少,从而能够在普通的个人电脑上有效求解大规模问题。在初始格式的基础上,进一步提出适用于二维弹性力学问题的新型快速多极边界元格式。从理论分析和数值实验上比较了初始和新型两种格式的存储与计算效率的差别。与初始格式相比,新型格式在略微增加存储量的前提下,计算速度又有一定程度的提高。将二维问题扩展到三维,建立三维弹性力学问题基本解的统一展开格式,从而能够有效处理固体力学中常见的混合边界条件。与二维问题不同,对于一般的求解规模,三维问题的初始算法的求解效率相比常规求解技术改善不明显,而新型算法的求解效率仍然有大幅度的提高。复合材料在工业领域得到了越来越广泛应用。本文利用快速多极边界元法对二维含有上千个随机分布夹杂的复合材料和三维含有上百个随机分布颗粒或纤维的复合材料进行大规模计算模拟,研究了复合材料的微观结构对宏观属性的影响,与部分经典理论近似解法比较,得到一些有重要参考价值的结论。总之,快速多极边界元法在大规模计算上的优势使边界元法有了更好的发展空间,在复合材料的模拟研究领域具有广泛的应用前景。
论文目录:
第1章 绪论
1.1 边界元法概述
1.2 快速多极边界元法的研究状况
1.2.1 快速多极算法的发展历史
1.2.2 快速多极边界元法的发展历史
1.2.3 用于边界元法的预处理技术概述
1.3 复合材料的边界元法数值模拟
1.4 本文研究的主要内容
1.4.1 用于二维弹性力学的初始快速多极边界元法
1.4.2 用于二维弹性力学的新型快速多极边界元法
1.4.3 用于三维弹性力学的初始和新型快速多极边界元法
1.4.4 复合材料的快速多极边界元模拟
第2章 用于二维弹性力学的初始快速多极边界元法
2.1 引言
2.2 二维弹性力学问题的边界元法
2.3 边界元的迭代求解算法
2.4 初始快速多极边界元法基本原理
2.4.1 基本解的多极展开
2.4.2 多极展开系数的传递
2.4.3 基本解的局部展开
2.4.4 局部展开系数的传递
2.4.5 自适应树结构
2.5 初始快速多极边界元法实现步骤
2.5.1 边界单元离散和自适应树结构的生成
2.5.2 上行遍历计算多极展开系数
2.5.3 下行遍历计算局部展开系数
2.5.4 利用树结构计算积分
2.5.5 迭代求解
2.6 快速多极边界元法的预处理技术
2.7 数值算例
2.7.1 初始快速多极边界元法的精度验证
2.7.2 初始快速多极边界元法的大规模计算精度验证
2.7.3 初始快速多极边界元法与高斯消去法的效率比较
2.7.4 初始快速多极边界元法不同展开截断阶次的效率比较
2.7.5 初始快速多极边界元法的迭代收敛性验证
2.8 本章小结
第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.3.4 下行遍历传递指数展开系数
3.3.5 下行遍历计算局部展开系数
3.3.6 利用树结构计算积分
3.3.7 迭代求解
3.4 初始算法与新型算法的效率比较分析
3.4.1 计算量比较
3.4.2 存储量比较
3.5 数值算例
3.5.1 展开阶数的选取
3.5.2 新型快速多极边界元法的精度验证
3.5.3 新型快速多极边界元法的大规模计算精度验证
3.5.4 初始、新型快速多极边界元法与高斯消去法的效率比较
3.6 本章小结
第4章 用于三维弹性力学的初始和新型快速多极边界元法
4.1 引言
4.2 三维弹性力学问题的边界元法
4.3 三维初始快速多极边界元法基本原理
4.3.1 基本解的多极展开
4.3.2 多极展开系数的传递
4.3.3 基本解的局部展开
4.3.4 局部展开系数的传递
4.3.5 自适应树结构
4.4 三维新型快速多极边界元法基本原理
4.5 三维初始算法与新型算法的效率比较分析
4.5.1 计算量比较
4.5.2 存储量比较
4.6 数值算例
4.6.1 三维展开阶数的选取
4.6.2 三维快速多极边界元法的精度验证
4.6.3 三维快速多极边界元法的大规模计算精度验证
4.6.4 三维快速多极边界元法与高斯消去法的效率比较
4.7 本章小结
第5章 复合材料的快速多极边界元大规模模拟
5.1 引言
5.2 边界元重复相似子域法
5.3 用于含夹杂复合材料结构的预处理技术
5.4 数值算例
5.4.1 二维含大量随机分布的圆形夹杂的复合材料大规模计算
5.4.2 三维含大量随机分布的球形颗粒的复合材料大规模计算
5.4.3 三维短纤维增强复合材料大规模计算
5.4.4 纤维长度和分布方向对复合材料的影响计算
5.5 本章小结
结论
参考文献
致谢与声明
附录
个人简历、在学期间发表的学术论文与研究成果
发布时间: 2006-06-29
参考文献
- [1].三维弹性问题边界元法并行计算及其工程应用[D]. 尹欣.清华大学2000
- [2].二维含缺陷弹性体移动和滚动接触的边界元法[D]. 蒲军平.清华大学2000
- [3].Galerkin边界元法及其在结构软化分析中的应用研究[D]. 刘清珺.清华大学1995
- [4].三维弹塑性分析的有限元—边界元耦合方法[D]. 岑章志.清华大学1984
- [5].若干弹性结构边界元优化问题及其显示系统[D]. 宋国书.清华大学1989
- [6].应用边界元法的弹性结构灵敏度分析及其形状优化[D]. 霍同如.清华大学1989
- [7].多裂纹体分析的边界元方法[D]. 秦飞.清华大学1996
- [8].弹塑性边界元法的若干基础性研究及在接触问题上的应用[D]. 董春迎.清华大学1992
- [9].快速多极边界元并行算法的研究与工程应用[D]. 雷霆.清华大学2006
- [10].涂层结构和V形切口界面强度的边界元法分析研究[D]. 程长征.合肥工业大学2007
相关论文
- [1].三维弹性体移动接触问题的边界元法研究[D]. 刘永健.清华大学2003
- [2].快速多极子边界元方法在完全非线性水波问题中的应用[D]. 宁德志.大连理工大学2006
- [3].边界型数学规划非线性多极边界元法[D]. 于春肖.燕山大学2006
- [4].平面正交各向异性体材料参数识别的边界元法[D]. 黄立新.清华大学2005
- [5].碳纳米管及其复合材料的力学性能研究[D]. 施冬莉.清华大学2005
- [6].颗粒增强复合材料界面脱层和热机疲劳的数值模拟[D]. 郭然.清华大学2003
- [7].三维弹性问题边界元法并行计算及其工程应用[D]. 尹欣.清华大学2000
- [8].一种新的边界类型无网格法[D]. 张见明.清华大学2002