用快速多极边界元法模拟碳纳米管复合材料.pdfVIP

用快速多极边界元法模拟碳纳米管复合材料 1 2 1 1 姚振汉 ，王海涛 ，徐俊东 ，郑小平 1 清华大学工程力学系，100084 北京 2 清华大学核能与新能源研究院，100084 北京 Email ：demyzh@tsinghua.edu.cn 摘要：此文简单介绍快速多极边界元法在模拟碳纳米管复合材料方面的应用，包括弹性特 性的模拟，以及稳态导热和导电特性的模拟。碳纳米管被模拟为具有与基体材料不同物理 性质的纤维状连续体，和基体之间满足理想的边界条件。数值算例表明快速多极边界元法 能胜任碳纳米管复合材料模拟的大规模计算，但是要使结果和实验结果相符还需对计算模 型做较大的改进，需考虑非理想界面条件，以及纤维之间的搭接。 关键词：快速多极算法，边界元法，复合材料模拟，碳纳米管 一、 引言 近年来快速多极边界元法（FM-BEM ）受到许多研究者的重视。近期研 究成果表明，FM-BEM能解决一些有限元法和其他数值方法难以求解的大规 模问题。当务之急是将这种方法真正用于解决科学与工程中的实际问题。另 一方面，碳纳米管正由于其特有的优良性能而受到工程界和研究者越来越多 的重视。因此碳纳米管复合材料的数值模拟可能是FM-BEM 的有价值的实际 应用之一。 基于作者对颗粒与纤维增强复合材料的快速多极边界元模拟，在模拟碳 纳米管复合材料的弹性性能时，将每根碳纳米管看作等效的弹性纤维，而在 模拟其热与电性能时，看作具有很高导热与导电性能的直纤维。为了提高计 算效率采用了重复相似子域法，而对大规模计算采用并行算法。 文中给出典型算例来验证算法的可行性。计算得到的等效材料特性和文 献中的实验结果进行了比较。 二、 重复相似子域法模拟碳纳米管的弹性特性 对于图 1 所示的代表性体元，对于基体和第 k 根碳纳米管的边界积分方 程可以分别列出为 1 C ( x)u ( x) + T 0 ( x, y)u ( y)dS( y) U 0 ( x, y)t ( y)dS( y) (1) ij j ∫S0 +∑Sk ij j ∫S0 +∑Sk ij j C ( x)u ( x) + T k ( x, y)u ( y)dS( y) U k ( x, y)t ( y)dS( y) (2) ij j ∫Sk ij j ∫Sk ij j 其中S 为基体子域外边界，S , S , ..., S , ..., S 为各碳纳米管与基体的界面。 0 1 2 k n 图 1 ：含200根随机分布碳纳米管的代表性体元 经离散后方程(1)和(2)可以写成如下矩阵形式： ⎡A0 A0 A1 L Ak L An ⎤ U ⎡B1 L B k L B n ⎤ ⎡B 0 B 0 ⎤ ⎧ ⎫ ⎢ 11 12 13 13 13 ⎥ 0 ⎢ 13 13 13 ⎥ ⎢ 11 12 ⎥ 0 0 1 k n ⎪ ⎪ 1 k n 0 0 0 A A A L A L A T ⎧ ⎫ B L B L B T B B ⎢ 21