学位论文 > 优秀研究生学位论文题录展示

非均质材料力学性能与失效分析的多尺度有限元法研究

作　者: 吴敬凯
导　师: 张洪武
学　校: 大连理工大学
专　业: 工程力学
关键词: 非均质材料多尺度计算力学性能应变局部化裂纹扩展
分类号: O344.3
类　型: 博士论文
年　份: 2013年
下　载: 107次
引　用: 0次
阅　读: 论文下载

内容摘要

自然界中存在的、以及工程应用中的大部分材料都具有多尺度特征,且当考察尺度小到一定程度后,都将表现为非均质性。对这类材料而言,其微观非均质特征对宏观行为影响非常大,实际上,材料在宏观尺度所表现出来的大部分物理力学现象大多是由微观非均质结构决定的。以航空航天工业中广泛使用的复合材料为例,其力学性能与各组分相的材料属性、尺寸、形状、空间分布及界面连接特性等细观特征密切相关。另一方面,材料的宏观失效通常是由微观缺陷、微观裂纹等引起并发展所致。所以,研究微观现象和宏观行为的相互影响已成为工程应用中的一个迫切需要解决的问题。当采用传统数值方法(如有限元法)直接求解这类多尺度特征问题时,常常因为需要大量的计算机存储和计算时间而遇到求解困难。多尺度计算方法作为求解这类问题的一种有效方法,已成为近年来的一个研究热点。本文首先提出了非均质材料弹性力学性能分析的扩展多尺度有限元法(EMsFEM)。基于多尺度有限元法(MsFEM)的思想,数值构造了能反映粗网格单元内部小尺度非均质信息的矢量场多尺度基函数。和标量场问题不同的是,在考虑固体变形的矢量场问题中材料内部会发生体积膨胀/收缩效应(泊松效应),从而导致材料各个方向的变形是相互耦合的。为了处理这种耦合关系,我们提出了一种新型的多尺度基函数构造方式,并在多尺度基函数中添加了耦合附加项。对不同构造方法下的多尺度基函数性能进行比较,结果表明,新型多尺度基函数能正确地反映粗网格单元内部的耦合变形,在相同的边界条件下可得到和常规精细有限元解完全一致的微观变形场,从而很大程度地提高了多尺度方法的计算精度。比较了不同边界条件下多尺度基函数的构造对多尺度方法精度的影响,并分析了尺度效应和误差产生的原因。所发展的方法实施方便,具有较好的数值精度。与传统均匀化方法相比,本文方法无尺度分离和周期性假设的限制,且易于进行降尺度计算以获得微观尺度的真实应力应变信息。其次,基于扩展多尺度有限元法基本思想,提出了一类新的平面矩形单元和空间长方体单元(广义单元)。通过弹性力学伽辽金方程,理论推导出单元内部位移场和节点值之间的函数关系,并应用于新单元形函数的构造中。新单元能更合理地考虑单元内部的耦合变形,在不增加计算自由度的情况下得到比传统单元更精确的结果,为构造精度更高的新型有限单元提供新的思路。进一步,发展了非均质材料小变形弹塑性分析的扩展多尺度有限元法。对于弹塑性材料的多尺度分析,随着微观物理变量的演化,微观尺度节点上将产生不平衡节点力。为了处理这些不平衡节点力产生的响应,我们提出了一种位移分解技术,将微观不平衡节点力响应等效成宏观等效力和微观扰动力的共同作用。其中宏观等效力用来驱动结构的宏观位移响应,而微观扰动力用来获得局部的微观位移扰动响应。在位移分解技术的基础上,进一步提出了两尺度协同迭代算法。在两尺度计算框架下,对宏微观非线性方程组进行协同增量迭代计算。最后,对方法的计算复杂度进行分析并与传统有限元法进行比较。数值结果表明,所发展的方法在模拟非均质材料弹塑性问题时具有很好的精度,且能显著地减少计算机存储需求和计算时间。此外,发展了周期性桁架材料应变局部化分析的自适应多尺度方法(AMM)。局部化现象是引起工程结构破坏的主要原因之一,传统基于体积平均思想的均匀化方法由于受到局部周期性假设和尺度分离假设的限制而不能很好地用于模拟这类问题。在计算过程中,AMM自适应地将整体计算域分解为两个部分,即细网格计算域和粗网格计算域,EMsFEM提供了多尺度框架。这样,变形剧烈的局部化域可采用有限元法在细尺度网格上直接求解,保证了算法的精度；而对于其它变形较缓和的域,则仍采用EMsFEM在粗网格尺度上进行非线性计算。同时,尺度连接通过建立在尺度界面上的主从节点约束方程完成。提出一种基于宏观单元位移梯度的网格识别器来自动更新细网格计算域。数值结果表明AMM具有较高的计算效率,能准确预测局部化问题的剪切带位置,且避免了宏观单元网格依赖性问题。最后,发展了一种模拟准静态裂纹扩展问题的自适应多尺度方法。采用扩展有限元法(XFEM)和水平集法(LSM)的组合对细尺度模型的不连续问题(包括强不连续和弱不连续问题)进行模拟,其中,XFEM用来计算不连续模型的应力应变场,而LSM用来对固定不连续界面(材料界面)和移动不连续界面(裂纹面)进行几何描述。这样,有限元网格可独立于不连续界面,模拟裂纹扩展问题时不需要对有限元网格重新剖分。而对于不需要精细分析的其它子域,采用EMsFEM在粗尺度网格上进行求解。提出一种基于距离参数的自适应网格识别技术来对求解域进行自适应分解,其距离参数可由水平集函数的演化方程自动获取。方法具有很好的精度,能准确地预测裂纹扩展路径,有效地减少裂纹扩展模拟所需的计算自由度和计算时间。

全文目录

摘要  4-6
Abstract  6-12
1 绪论  12-27
  1.1 研究的背景及意义  12-15
  1.2 多尺度分析方法研究进展  15-22
    1.2.1 非耦合层级多尺度方法  16-18
    1.2.2 耦合层级多尺度方法  18-21
    1.2.3 嵌入式协同多尺度方法  21-22
  1.3 多尺度有限元法介绍  22-25
    1.3.1 多尺度有限元法的研究概况  22-24
    1.3.2 多尺度有限元法的基本思想  24-25
  1.4 本文主要研究内容  25-27
2 非均质材料弹性分析的扩展多尺度有限元法  27-58
  2.1 引言  27
  2.2 扩展多尺度有限元法基本原理  27-34
    2.2.1 多尺度基函数构造  29-33
    2.2.2 宏观有限元计算  33-34
    2.2.3 降尺度计算  34
  2.3 多尺度基函数性能比较  34-38
  2.4 其它边界条件  38-41
    2.4.1 基于超样本技术的振荡边界条件  38-40
    2.4.2 周期性边界条件  40-41
    2.4.3 基于超样本技术的周期性边界条件  41
  2.5 数值算例与分析  41-57
    2.5.1 周期性点阵类桁架材料力学性能分析  41-47
    2.5.2 连续体材料力学性能分析  47-57
  2.6 本章小结  57-58
3 一类新的平面矩形单元和空间长方体单元  58-71
  3.1 引言  58
  3.2 传统平面矩形单元(R4)简介  58-60
  3.3 广义平面矩形单元(GR4)与空间长方体单元(GC8)  60-65
    3.3.1 平面问题  60-63
    3.3.2 空间问题  63-65
  3.4 分片试验  65
  3.5 单元基本变形模态  65-67
  3.6 数值算例与分析  67-69
  3.7 本章小结  69-71
4 非均质材料弹塑性分析的扩展多尺度有限元法  71-99
  4.1 引言  71
  4.2 塑性力学的基本方程  71-72
  4.3 弹塑性分析的两尺度协同迭代算法  72-79
    4.3.1 微观节点力的处理  72-76
    4.3.2 两尺度协同迭代算法描述  76-79
  4.4 数值算例与分析  79-97
  4.5 计算复杂度比较  97-98
  4.6 本章小结  98-99
5 周期性桁架材料应变局部化分析的自适应多尺度方法  99-115
  5.1 引言  99
  5.2 自适应多尺度方法框架  99-101
  5.3 网格识别器  101-102
  5.4 尺度连接  102
  5.5 求解步骤  102-103
  5.6 数值算例与分析  103-113
  5.7 本章小结  113-115
6 裂纹扩展模拟的自适应多尺度方法  115-143
  6.1 引言  115-117
  6.2 线弹性断裂力学基本理论  117-123
    6.2.1 裂纹类型  117
    6.2.2 裂纹尖端附近的应力和位移场  117-120
    6.2.3 开裂准则  120-121
    6.2.4 应力强度因子计算  121-123
  6.3 小尺度不连续模型的建模  123-129
    6.3.1 扩展有限元法基本原理  124-126
    6.3.2 固定和移动不连续界面的几何描述  126-129
  6.4 不连续问题的自适应多尺度方法  129-133
    6.4.1 求解域的分解  129-131
    6.4.2 尺度连接  131-132
    6.4.3 求解步骤  132-133
  6.5 数值算例与分析  133-141
  6.6 本章小结  141-143
结论  143-146
参考文献  146-157
附录A 非均质材料多尺度分析程序介绍  157-161
创新点摘要  161-162
攻读博士学位期间科研项目及科研成果  162-163
致谢  163-164
作者简介  164-165

非均质材料力学性能与失效分析的多尺度有限元法研究

内容摘要

全文目录

相似论文