QC法介绍.pptVIP

多尺度分析中QC法介绍 * 主要内容 多尺度分析研究背景 QC法的局部化计算方法 QC法的非局部化的计算方法 及两种计算方法的连接 * 研究背景 现实需求：随着在微米和纳米尺度领域内科学技术的发展，对于该领域内材料的数值模拟成为热点问题 局部原子运动分布上的不同将影响在该尺度下材料性能，连续介质力学无法进行模拟计算。 目前计算机的计算能力无法在原子理论下对大规模的原子进行模拟计算。 理论诉求： “还原论” 与“大统一”思想 * 钢在不同尺度下的计算模型[1] * 主要内容 多尺度分析研究背景 QC法的局部化计算方法 QC法的非局部化的计算方法 及两种计算方法的连接 * QC（Quasicontinuum）法有限元网格的划分 Form References [2] * QC法中局部化计算方法 应用Cauchy－Born法则计算有限元区域内原子间变形梯度，从而求出应变能密度。 Cauchy－Born法则：在均匀变形情形下, 连结初始未变形晶体上某两个原子间的晶键矢a 和变形后晶体上相应原子间的晶键矢量b 的关系为b = F ·a, 这里的F 为变形梯度张量[3]。该准则适用于空间块状晶体。 在CB法则下，晶体的应变能密度为变形梯度的函数，可根据相应的理论有已知的变形梯度直接求出。 * QC法局部化计算中有限元网格的自动划分 划分准则： 根据误差e的计算结果判定是否重新划分网格 * QC法局部化计算中存在的问题 单元内的应变能密度不能准确的反映单元边界上原子的应变能密度。 Form References [2] * 主要内容 多尺度分析研究背景 QC法的局部化计算方法 QC法的非局部化的计算方法 及两种计算方法的连接 * QC法中非局部计算方法 通过原子能量计算的理论（如嵌入原子势法等）计算代表原子的能量，乘以其相应的权重系数，从而求得总能量。 Form References [2] * 采用QC法中非局部计算方法的判断准则 判断准则： a,b为某代表原子为圆心，半径rnl内的任意有限元单元； rnl一般取原子计算中的截断半径rcut的两到三倍。 * QC法两种计算方法的交界面 Form References [2] * 鬼力修正技术 * * 修正后的能量函数： * 参考文献： [1] Liu W K, Karpov E G, Zhang S and Park H S 2004 《An introduction to computational nano mechanics and materials》 Comput. Methods Appl. Mech. Eng. 193 1529–78 [2] Miller R E and Tadmor E B 2002 《The Quasicontinuum Method: Overview, applications and current directions》 J. Computer-Aided Materials Design 9 203–239 [3] 郭旭 王晋宝 张洪武 2005 《基于高阶Cauchy-Born 准则的单壁碳纳米管本构模型》 计算力学学报 Vol.22 No.2 [4] Tadmor E B Ortiz M and Phillips R 1996《Quasicontinuum analysis of defects in solids》 Phil. Mag. A 73 1529–63 [5] 王建华 杨磊 沈为平 2000 《有限元后验误差估计方法的研究进展》 力学进展 Vol.30 No.2 * Thank you for your attention