第二章材料科学中常用的数值

第二章材料科学研究中常用的数值分析方法有限差分法和有限元法第一节有限差分法一、概述以有限差分代替无限微分、以差分代数方程代替微分方程、以数值计算代替数值推导的过程，从而将连续函数离散化，以有限的、离散的数值代替连续的函数分布。有限差分法的主要步骤如下：1）构成差分格式。首先选择网格布局、差分形式和步长；其次，以有限差分代替无限微分，以差商代替微商，以差分方程代替微分方程及边界条件。2）求解差分方程：直接法和间接法3）对所得到的数值解进行精度与收敛性分析和检验。二、差分方程的建立导出差分方程的途径有两种：一是从微分方程出发，以泰勒级数截断，从有限差分的数学含义去建立有限差分和差分方程；另一种则是从由网格所划分的单元体的能量平衡出发，由积分方程去建立差分方程，该方法又称单元体平衡法。两种方法各具特色，但无论采取哪种方法，在建立差分方程前，均需对所论区域进行离散化。1、合理选择网格布局及步长

区域离散化离散化网格的选择有两种方法：物理划分法和以几何区域形状为依据去划分。

2、将微分方程转化为差分方程将微分方程转化为差分方程实际上就是以差分代替微分、以差商代替微商，是以有限小量代替无限微量的近似化过程。（1）差分所谓差分，就是某物理量的有限增量。1）向前差分2）向后差分3）中心差分（2）差商差商为函数的差分与自变量差分之比。

3.有限差分法的求解步骤第二节有限元法finiteelementmethod

1用有限个单元将连续体离散化，通过对有限个单元作分片插值求解各种力学、物理问题的一种数值方法。有限元法把连续体离散成有限个单元：杆系结构的单元是每一个杆件；连续体的单元是各种形状（如三角形、四边形、六面体等）的单元体。每个单元的场函数是只包含有限个待定节点参量的简单场函数，这些单元场函数的集合就能近似代表整个连续体的场函数。根据能量方程或加权残量方程可建立有限个待定参量的代数方程组，求解此离散方程组就得到有限元法的数值解。有限元法已被用于求解线性和非线性问题，并建立了各种有限元模型，如协调、不协调、混合、杂交、拟协调元等。有限元法十分有效、通用性强、应用广泛，已有许多大型或专用程序系统供工程设计使用。结合计算机辅助设计技术，有限元法也被用于计算机辅助制造中。

2有限元分析（FEA，FiniteElementAnalysis）的基本概念是用较简单的问题代替复杂问题后再求解。它将求解域看成是由许多称为有限元的小的互连子域组成，对每一单元假定一个合适的(较简单的）近似解，然后推导求解这个域总的满足条件(如结构的平衡条件），从而得到问题的解。这个解不是准确解，而是近似解，因为实际问题被较简单的问题所代替。由于大多数实际问题难以得到准确解，而有限元不仅计算精度高，而且能适应各种复杂形状，因而成为行之有效的工程分析手段。

3有限元是那些集合在一起能够表示实际连续域的离散单元。有限元的概念早在几个世纪前就已产生并得到了应用，例如用多边形（有限个直线单元）逼近圆来求得圆的周长，但作为一种方法而被提出，则是最近的事。有限元法最初被称为矩阵近似方法，应用于航空器的结构强度计算，并由于其方便性、实用性和有效性而引起从事力学研究的科学家的浓厚兴趣。经过短短数十年的努力，随着计算机技术的快速发展和普及，有限元方法迅速从结构工程强度分析计算扩展到几乎所有的科学技术领域，成为一种丰富多彩、应用广泛并且实用高效的数值分析方法。

直接字刚度耽法对于坡不同垃物理路性质虑和数稀学模谎型的归问题远，有关限元茶求解美法的珍基本摩步骤厉是相腹同的要，只络是具催体公角式推凤导和窗运算根求解延不同挽。有五限元熟求解归问题烧的基怠本步糕骤通诵常为防：第一互步：鸣问题缓及求酬解域郊定义坛：根碌据实刑际问屋题近醋似确津定求水解域待的物督理性印质和胀几何两区域肃。第二昂步：拌求解擦域离潜散化抓：将捎求解寇域近罗似为多具有葬不同锁有限览大小竹和形伐状且扁彼此爹相连呆的有摊限个昂单元错组成特的离描散域欧，习族惯上先称为诸有限秃元网至络划墙分。苍显然做单元输越小迹（网祸络越肯细）前则离紧散域浪的近反似程倍度越抽好，珠计算妄结果所也越招精确显，但匪计算循量及鹿误差破都将恰增大苍，因壤此求选解域苗的离难散化辩是有男限元答法的除核心抹技术积之一茧。第三箩步：详确定教状态敢变量商及控且制方笑法：处一个迈具体修的物洞理问亮题通睁常可午以用届一组顺包含懒问题镜状态惧变量其边界桌条件异的微遣分方串程式稠表示级，为仔适合装有限右元求昆解，碧通常秆将微程分方篮程化逆为等砌价的美泛函侨形式姿。第四丝式步：宝单元晴推导弊：对失单元俊构造霜一个待适合涌的近蛛似解匀，即乌推导握有限柄单元责的列期式，匀其中股包括画选择湖合理桨的单剩元坐雅标系使，建局立单机元试因函数版，以妄某种额方法武给出质单元熟各状著态变索量的玩离散梨关系手，从愁而形讨成单省元矩逮阵（芽结构潮力学纵中称安刚度派阵或岔柔度粥阵）巨。趁为保拌证问去题求惠解的愧收敛劲性，躬单元弦推导倍有许嘱多原国则要秧遵循拘。探对工钢程应陶用而叼言，叨重要碌的是士应注书意每亡一种侵单元更的解败题性跟能与校约束艇。例双如，出单元殊形状遣应以送规则肉为好观，畸辫形时驳不仅产精度缓低，序而且勒有缺牢秩的房诚危险舞，将惭导致摇无法贯求解恒。第五呢步：迎总装门求解涛：将肤单元火总装碎形成棒离散奶域的明总矩战阵方求程（迈联合冰方程松组）预，反汉映对椒近似爽求解沉域的苗离散位域的下要求惜，即宗单元绪函数妨的连嘴续性当要满攀足一隙定的沉连续弱条件辣。总往装是辟在相懒邻单悬元结越点进惨行，击状态辞变量日及其脆导数惩（可推能的螺话）扇连续水性建扶立在拦结点违处。第六很步：登联立议方程瞒组求悬解和甚结果碰解释烈：有慢限元扯法最迹终导紫致联器立方算程组睡。联谱立方火程组戏的求纽奉解可乱用直泛接法铁、选糠代法搜和随范机法话。求兽解结革果是普单元雄结点种处状沸态变粘量的贪近似熟值。补对于粮计算赤结果恨的质战量，针将通氏过与闯设计捎准则叔提供等的允觉许值霸比较押来评稀价并揉确定词是否怠需要凉重复勿计算粮。简言贱之，吨有限泥元分史析可兆分成槽三个饲阶段捞，前处兼理、钢处理窜和后践处理。前辛处理菜是建劫立有胳限元眨模型鹊，完归成单睛元网毫格划软分；骂后处何理则煤是采选集处认理分躺析结皮果，重使用减户能穷简便回提取溉信息蛋，了栽解计判算结息果。有限里元分雪析，宁表明疫上看谈起来脆简单哭，实盆际上甘是有喜较大区科技嗽含量倡的一形项辛兽苦工衡作。典在一涉些情路况下拣，找恋出分疗析对辆象的斩载荷炸和边荡界条绣件是似非常迈困难攻的，熔要对矩分析侄对象旋的使动用情愉况进宫行仔浊细的旬观察痕和分角析，锅最后遥还要叫对计询算结姿果进青行分草析和唤确认输，不斥是像膀有人窑说的报那样款做有僻限元霸分析敢很容臂易，孝所谓蔬一按持电钮物，计摄算结曾果就津出来昼了。罩看别丘人有滚限元玻分析鲁容易骂，自当己做从起来猴就不服一定借容易僚。有川限元划分析芳入门给容易含，但扫做好千分析鹿就不帆容易颠了。迁国内慨外的棵许多宫实例霉表明椅，有渣限元挠分析气人员帐责任锡重大额，做驼好有罗限元盛分析春就能希给公带司和造企业衣带来滴辉煌挽，做受不好就有限激元分目析就衔会误跳导设包计，够给公腔司和蜓企业卷带来叠损失膏。有刊限元涝建模灿分析迷，60％在吧于理术论，40％在巴于经贿验。乱有限径元分