平衡单元法——米力学多尺度数值模拟方法研究.pdf

摘要(ABSTRACT) 摘要 本文基于能量平衡原理研究纳米材料的力学性能．微观层次应用经典分子动 力学模拟纳米结构薄膜。设置初始边界条件以及内部原子形态分布下的晶体结 构，研究其在小变形条件下的能量改变。建立了一种平衡单元法模型，应用此模 型计算材料在与微观原子位移等效的宏观变形下能量的变化量，进而由能量平衡 求解材料的弹性模量。 分子动力学模拟方法通过牛顿运动方程来模拟计算材料内部原子的运动情 况。以期求得不同原子在任意时刻，任一坐标点的速度，建立速度和待求置之润 的相互关系。模拟过程就是寻求合适的势函数，根据势函数求得原子之间的相互 作用力，通过牛顿运动方程建立求解方程，应用差分法求得各个原子的速度。根 据啥密顿原理而得到的势函数表明了原子之间的相互作用，并由此可得到原子之 间的棍互作用力F，，由薛定谔微正则方程组可知F。力是势函数的梯度，其数懂取 决于原子的坐标位置，故而模拟实现的关键一步在于势函数的选取。以及用其来 计算材料在不同的边界条件下内在性质的改变。分子动力学模拟方法包含两个重 要部分：势竭数的选取是关键，通过牛顿运动方程建立方程组，应用差分法求解 是核心。本文应用经典对势分子动力学模拟了纳米材料微观的力学性能，对视始 边界条件下的晶体结构，研究了小交形情况下材料内部能量的改变。模拟了较多 纳米原子组成的薄膜，材料长度单位j(1010m)、能量单位ev、温度单位k 近年来，已有很多学者研究和建立了多种描述纳米材料结构的连续模型， Odgard翎应用平衡柬模型建立了计算化学和固体力学之间的联系，经过宏观有 限元得到材料的力学性质。Sun将原子间的相互作用简化为弹簧连接，由此建立 了一种平板连续模型。Oow应用一种连续模型来研究航空材料的力学性质。所有 这些不同方法都包含一些自身的假定，适合于不同的具体问题，而且微观层析能 量的计算进行了较大的简化。有的模型直接由单个原子作为研究的对象，完全忽 略了量子效应，从而导致能量计算的误差较大。本文建立了一种平衡单元法模型， 这种模型以微观——宏躐能量平衡为基础，微观层次下应用分予动力学模拟纳米 材料在边界条件下的能量变化，避免以往学者建立连续模型时能量计算误差较大 的缺陷，并且在分子动力学模拟中本文建立的初始原子模型允许单个原子在平衡 撬要(ABSTRACT) 位置可以微小波动，这样更接近于原子的实际运动情况。同时，本文建立了一种 微观——宏观边界条件的关系，这是平衡单元法连接微观——宏观的枢纽。借鉴 Whitomb边界条件的处理方法，研究材料在宏观小变形情况下的变形能以及在等 效于宏观变形下的微观原子位移状态下分子动力学模拟所得材料能薰变化情况， 进而建立宏观——微观能量平衡方程。 平衡单元法以较多原子团簇为研究对象，建立材料求解分析的宏观单元模 型。单元内部必须包含足够数目的原予，这样即保证了单元内部原子非空的缺陷， 也可减小量子效应的作用程度。在分子动力学模拟中势函数的作用范围在截断距 离Rc之内，本文宏鹰单元的边长采用高于一个量级的方法来保证单元内部包 含足够的原子。单元之间假定为连续情况，宏观边界条件施加到纳米薄膜之上。 总之，应用分子动力学模拟材料在小变形情况下内部能量的变化，建立了一 套平衡单元法模型，规定了单元的选取原则以及宏观边界条件与微观原予边界条 件的等效关系。微观结构材料的能量应用分子动力学进行模拟，这样计算所得微 观结构能量更接近实际情况，通过求解宏观——微观能量平衡方程得到纳米薄膜 材料的弹性模量E，并模拟了不同纳米材料在不同工况下的力学性能≯ 关键词：平衡单元法，多尺度纳米力学，分子动力学模拟，六边形单元 Abstract findthe ofnano Basedon methodto mechanic energy properties equivalent filmMDmethods