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

来源 :上海大学 | 被引量 : 0次 | 上传用户:fanjun6699
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文基于能量平衡原理研究纳米材料的力学性能。微观层次应用经典分子动力学模拟纳米结构薄膜。设置初始边界条件以及内部原子形态分布下的晶体结构,研究其在小变形条件下的能量改变。建立了一种平衡单元法模型,应用此模型计算材料在与微观原子位移等效的宏观变形下能量的变化量,进而由能量平衡求解材料的弹性模量。 分子动力学模拟方法通过牛顿运动方程来模拟计算材料内部原子的运动情况,以期求得不同原子在任意时刻,任一坐标点的速度,建立速度和待求量之间的相互关系。模拟过程就是寻求合适的势函数,根据势函数求得原子之间的相互作用力,通过牛顿运动方程建立求解方程,应用差分法求得各个原子的速度。根据哈密顿原理而得到的势函数表明了原子之间的相互作用,并由此可得到原子之间的相互作用力Fi,由薛定谔微正则方程组可知Fi力是势函数的梯度,其数值取决于原子的坐标位置,故而模拟实现的关键一步在于势函数的选取,以及用其来计算材料在不同的边界条件下内在性质的改变。分子动力学模拟方法包含两个重要部分:势函数的选取是关键,通过牛顿运动方程建立方程组,应用差分法求解是核心。本文应用经典对势分子动力学模拟了纳米材料微观的力学性能,对初始边界条件下的晶体结构,研究了小变形情况下材料内部能量的改变。模拟了较多纳米原子组成的薄膜,材料长度单位0A(10-10m)、能量单位ev、温度单位K。 近年来,已有很多学者研究和建立了多种描述纳米材料结构的连续模型,OdgardGM应用平衡束模型建立了计算化学和固体力学之间的联系,经过宏观有限元得到材料的力学性质。Sun将原子间的相互作用简化为弹簧连接,由此建立了一种平板连续模型。Dow应用一种连续模型来研究航空材料的力学性质。所有这些不同方法都包含一些自身的假定,适合于不同的具体问题,而且微观层析能量的计算进行了较大的简化,有的模型直接由单个原子作为研究的对象,完全忽略了量子效应,从而导致能量计算的误差较大。本文建立了一种平衡单元法模型,这种模型以微观——宏观能量平衡为基础,微观层次下应用分子动力学模拟纳米材料在边界条件下的能量变化,避免以往学者建立连续模型时能量计算误差较大的缺陷,并且在分子动力学模拟中本文建立的初始原子模型允许单个原子在平衡位置可以微小波动,这样更接近于原子的实际运动情况。同时,本文建立了一种微观——宏观边界条件的关系,这是平衡单元法连接微观——宏观的枢纽。借鉴whitomb边界条件的处理方法,研究材料在宏观小变形情况下的变形能以及在等效于宏观变形下的微观原子位移状态下分子动力学模拟所得材料能量变化情况,进而建立宏观——微观能量平衡方程。 平衡单元法以较多原子团簇为研究对象,建立材料求解分析的宏观单元模型。单元内部必须包含足够数目的原予,这样即保证了单元内部原子非空的缺陷,也可减小量子效应的作用程度。在分子动力学模拟中势函数的作用范围在截断距离Rc之内,本文宏观单元的边长采用高于Rc一个量级的方法来保证单元内部包含足够的原子。单元之间假定为连续情况,宏观边界条件施加到纳米薄膜之上。 总之,应用分子动力学模拟材料在小变形情况下内部能量的变化,建立了一套平衡单元法模型,规定了单元的选取原则以及宏观边界条件与微观原子边界条件的等效关系。微观结构材料的能量应用分子动力学进行模拟,这样计算所得微观结构能量更接近实际情况,通过求解宏观——微观能量平衡方程得到纳米薄膜材料的弹性模量E,并模拟了不同纳米材料在不同工况下的力学性能。
其他文献
框架剪力墙结构体系和薄壁异型柱结构体系在多高层建筑结构中得到了广泛的应用 但是迄今为止 相关研究主要集中在薄壁异型柱或剪力墙本身的受力性能和计算分析方法上 而对宽
CTOD(Crack Tip Opening Displacement)和J积分是弹塑性断裂力学研究方法中的两个重要参数.CTOD不是一个直接而严密的应力、应变场参量,其确切定义和直接测量还存在困难,但由
CAN由于其突出的可靠性和灵活性,已成为当代汽车应用最广泛的现场总线.但是标准CAN协议没有提供足够的安全措施,易遭受窃听、重放、泛洪、拒绝服务攻击.为了有效检测CAN总线
多孔介质在实际生活中是广泛存在的。多孔介质的力学行为是生物力学和岩土工程等领域关注的重要方面,具有广泛的研究和应用前景。  饱和多孔介质在实际生活中广泛存在,其动力
高压直流输电(HVDC)线路在正常运行状态下是工作于双极运行方式的。直流输电电流通过两极的输电线构成回路。但在系统调试或发生故障的情况下,会采用单极或双极以大地为回路的运行方式。此时,巨大的直流电流经直流接地极流入大地,并在很大范围内造成地电位的明显变化。这种地电位的变化,对于同一地区的交流系统可能产生不利影响。因此,有必要研究直流单极运行时对交流系统的影响。为此,首先需要研究的是直流经大地回路运
多稳态结构具有多个稳定构型,通过不同构型间的转换可实现结构变体,从而成为变体结构重要的可选构型。变体结构稳态间的转换是一个短时动态过程,因此,变体结构形状控制需要准
通过对近十多年的破坏性地震所造成的人员伤亡及经济损失的分析,提出对地震断层错动所引起的地表破裂进行评价的必要性,并对地震地表破裂的类型进行了阐述。整理了拉梅常数不
随着电力系统的不断发展,系统容量的不断提高,对高压电器行业提出了更高的要求。对断路器灭弧室内电场的计算机辅助计算和分析,不仅可以缩短产品设计周期、减少试验的费用和盲目性,也为设计、研制和开发新产品创造了良好的条件。因此,如何采用先进的计算机技术和科学的计算方法对断路器的电场进行精确的计算成为断路器研究中的主要内容之一。本文采用有限元法对高压SF_6断路器灭弧室内电场分布进行数值仿真求解,计算分析了
风电资源的合理利用和开发在改善生态环境、优化能源结构、提升社会效益、促进能源和社会可持续发展等方面有着显著的作用,但由于风电系统的复杂性、时变性以及运行过程中的强干扰等因素,对风电系统进行精确的控制仍然是当前研究的一个难题。因此,研究风电技术尤其是控制技术对促进风电产业高效发展具有重要意义。本课题将大型直驱式风电机组变桨距控制技术作为研究重点,通过将支持向量机平台和机器学习方法应用到风电机组统一和