PM2.5原理介绍
一、产品原理
该款产品主要是通过电气原理,利用脉冲、振荡电器将低电压升至直流负高压，利用碳毛刷尖端直流高压产生高电晕，高速地放出大量的电子（e-），而电子并无法长久存在于空气中（存在的电子寿命只有nS级），立刻会被空气中的中性分子如氧分子（O2）捕捉，形成负离子，再通过净化器内部的风机将负离子吹出来，达到净化空气和对人体保健作用。
二、产品优势
1、负离子输出浓度高。有环境部门进行过检测，在城市中，每立方厘米含有的负离子在室外也只有几百个，室内就更少了，而该产品输出的负离子浓度高达6.4×106pcs/cm3。
2、人能够快速大量呼吸到负离子。因其小巧，悬挂在胸骨前，负离子直接对人输出，人能尽快呼吸到清新的负离子，而且随时呼吸到的健康安全的空气，防止各种病毒传染。
3、如果在室内或车内使用松池净化器，可以不密闭门窗。其他厂家的净化器产品，几乎都不要求房间一定密闭净化，如果人在一个完全密闭的空间生活，氧气的数量逐渐减少，即使空气净化得很干净，但
是缺氧的空气对人来说还算好吗？开启门窗，进入室内或车内的空气有氧气和其他气体、尘埃等，松池净化器直接对人输出的负离子迅速净化人周围的小环境，使人呼吸到清新洁净的空气，同时车内补充了大量的氧气。
4、安全、高效、环保、节能。松池净化器插上电源，开启电源开关便能使用，操作简单，功率小、噪音小，无副作用。符合国家低碳、环保、节能要求。
5、外观精致小巧。
6、使用成本低。无需支付频繁更换滤网的费用。空气中的可吸粒物多了，用过滤网过滤达到PM2.5，其孔必须极其小，那么就很容易被堵住，导致净化不起作用，只能经常更换滤网，而滤网成本较高。
7、对灭菌和除尘PM2.5，有较好作用。在医学界，负离子被确认是具有杀灭病菌及净化空气的有效手段。使用松池产品几天后，就会发现在产品输出口和负离子发射头上有一些黑色的污物附在上面，用手去摸，就能感觉到那是很多极其细小的尘埃被负离子吸附过来而集聚在一起形成的，这就明显证明了负离子除尘的功效。
三、产品功效
1.改善空气质量
能有效的杀菌、除尘（PM2.5）、除臭、除异味，消除“二手烟”污染。切断病毒空气传播，防止交叉感染。
2.清除车内污染
能快速清除车内污染：甲醛、三苯系物、TOVC等一系列对人体有强烈危害的有毒有害物质。3.补充负氧离子
负离子被誉为“环境警察”、“空气维生素”、“长寿素”。
高浓度负离子发射头，源源不断输出负离子，保证其空气中的负离子含量远远大于日常空气中负离子的含量，让你置身于有松有池的美妙感受之中，不容易疲劳。
4.促进身体健康
增强肌体免疫抗病能力;
促使红细胞、血红蛋白、血钙的增加，降低血糖、血脂和血液粘稠度；促进机体新陈代谢和基础代谢；
改善心脏功能和心肌营养不良状况；
有效预防呼吸道疾病和改善肺功能等;
改善内分泌功能，如提高性腺功能和改善甲状腺功能等
有降低血压的辅助治疗作用;
有效预防“空调病”;
消除疲劳、改善睡眠、提振精神；
中和电视、电脑的静电，形成负离子保护层，有效预防近视；
对颈椎病、腰间盘突出、类风湿关节炎风湿关节炎有辅助治疗作用。
四、负离子对部分疾病的辅助治疗和预防原理
1、负离子对高血压的功效
负离子可有效扩张冠状动脉，调整心率，达到降低血压的作用。通过心电图X线发现，负离子可有效改善心脏功能和心肌营养不良状况，对高血压和高血脂等心脑血管疾病有良好作用。
2、负离子治疗失眠的功效
负离子通过促进单胺氧化酶（MAO）的氧化脱氨基作用，降低脑及组织内的5-HT（5-羟色胺）水平，调节大脑皮层的兴奋和抑止趋于正常，从而提高脑力活动和改善睡眠。
3、负离子对哮喘的功效
负氧离子能有效加强气管粘膜上皮的纤毛运动，影响上皮绒毛内呼吸酶的活性，促使哮喘病人的肺功能及肺泡的分泌功能恢复正常。
4、负离子对糖尿病的功效
负氧离子使酸性化的生物体变成弱碱性，促进细胞组织的胶原蛋白生成，并分解、燃烧糖，降低血糖，并改善糖尿病的神经，血管机能障碍等并发症，而且未见明显副作用。
5、负离子对抗癌症防衰老的功效
自由基破坏人体细胞及DNA，导致人体衰老或致癌，负离子通过调节体内的酸碱失衡及氧化还原状况失衡，消除自由基，达到延缓衰老，减轻，消除化疗的不良负作用，抑止癌细胞的扩散。负离子能有效将家装污染中的甲醛、笨等致癌物质分解成无毒无味的二氧化碳和水。
6、负离子对增强免疫力、提高智力的功效
可增加神经元突触数目，加强连接部的强度，扩大有效神经网络，达到提高儿童智能的效果；可调节大脑的PH值，Ph值越高，人的智商就越高，人的免疫力也越强。
7、负离子对孕婴养护的功效
改善孕妇呼吸环境，促进胎儿生长和智力发育；缓解孕妇疲劳；改善孕妇睡眠，有利孕妇身心健康；增强孕妇及胎儿的体质；给孕妇创造一个健康的生活环境。
8、负离子对不孕不育有辅助治疗功效
五、检测负离子存在的方法
方法1：在一个开放的空间里，将负离子发生器接通电源后，在把风机关上的情况下，距其碳毛刷正前端0-30cm处，利用空气负离子浓度测试仪测得其浓度。
方法2：用电笔靠近净化器输出窗口，电笔发亮，表明有负离子办出．
方法3：用万用表打在MV档或最低的Ｖ档，一手紧握接地黑表笔笔头（金属部分），另一只手拿红表笔的绝缘部分，将其表笔头靠近负离子输出口，表显示为负数，则有负离子输出。
方法4：取一个有盖的透明小容器，将烟尘置于其中，再将负离子发生器放入其中，盖上通电1-8秒后，您将会发现容器内烟尘消失。
方法5：将负离子发生器置于有异味的环境中，一段时间后，您将会发现异味明显被消除。
第一性原理计算原理和方法
第二章 计算方法及其基本原理介绍 化学反应的本质就是旧键的断裂与新建的形成,参与成键原子的电子壳层重新组合就是导致生成稳定多原子化学键的明显特征。因此阐述化学键的理论应当描写电子壳层的相互作用与重排,借助求解满足适当的Schrodinger 方程的波函数描写分子中电子分布的量子力学,为解决这一问题提供了一般的方法,然而,对于一些实际的体系,不引入一些近似,就不可能求解其Schrodinger 方程。这些近似使一般量子力学方程简化为现代电子计算机可以求解的方程。这些近似与关于分子波函数的方程形成计算量子化学的数学基础。 2、1 SCF-MO 方法的基本原理 分子轨道的自洽场计算方法 (SCF-MO)就是各种计算方法的理论基础与核心部分,因此在介绍本文计算工作所用方法之前,有必要对其关键的部分作一简要阐述。 2、1、1 Schrodinger 方程及一些基本近似 为了后面介绍各种具体在自洽场分子轨道(SCF MO)方法方便,这里将主要阐明用于本文量子化学计算的一些重要的基本近似,给出SCF MO 方法的一些基本方程,并对这些方程作简略说明,因为在大量的文献与教材中对这些方程已有系统的推导与阐述[1-5]。 确定任何一个分子的可能稳定状态的电子结构与性质,在非相对论近似下,须求解 R AB =R 图2-1分子体系的坐标
定态Schrodinger 方程 ''12121212122ψψT p B A q p A p pA A pq AB B A p A A A E R Z r R Z Z M =??????? ?-++?-?-∑∑∑∑∑∑≠≠ (2、1) 其中分子波函数依赖于电子与原子核的坐标,Hamilton 算符包含了电子p 的动能与电子p 与q 的静电排斥算符, ∑∑≠+?-=p q p pq p e r H 12121?2 (2、2) 以及原子核的动能 ∑?-=A A A N M H 2121? (2、3) 与电子与核的相互作用及核排斥能 ∑∑≠+-=p A B A AB B A pA A eN R Z Z r Z H ,21? (2、4) 式中Z A 与M A 就是原子核A 的电荷与质量,r pq =|r p -r q |,r pA =|r p -R A |与R AB =|R A -R B |分别就是电子p 与q 、核A 与电子p 及核A 与B 间的距离(均以原子单位表示之)。上述分子坐标系如图2、1所示。可以用V(R,r)代表(2、2)－(2、4)式中所有位能项之与 ∑∑∑-+=≠≠p A pA A B A q p pq AB B A r Z r R Z Z r R V ,1 2121),( (2、5) 原子单位 上述的Schrodinger 方程与Hamilton 算符就是以原子单位表示的,这样表示的优点在于简化书写型式与避免不必要的常数重复计算。在原子单位的表示中,长度的原子单位就是Bohr 半径
有限元法的基本思想及计算 步骤
有限元法的基本思想及计算步骤 有限元法是把要分析的连续体假想地分割成有限个单元所组成的组合体，简称离散化。这些单元仅在顶角处相互联接，称这些联接点为结点。离散化的组合体与真实弹性体的区别在于：组合体中单元与单元之间的联接除了结点之外再无任何关联。但是这种联接要满足变形协调条件，即不能出现裂缝，也不允许发生重叠。显然，单元之间只能通过结点来传递内力。通过结点来传递的内力称为结点力，作用在结点上的荷载称为结点荷载。当连续体受到外力作用发生变形时，组成它的各个单元也将发生变形，因而各个结点要产生不同程度的位移，这种位移称为结点位移。在有限元中，常以结点位移作为基本未知量。并对每个单元根据分块近似的思想，假设一个简单的函数近似地表示单元内位移的分布规律，再利用力学理论中的变分原理或其他方法，建立结点力与位移之间的力学特性关系，得到一组以结点位移为未知量的代数方程，从而求解结点的位移分量。然后利用插值函数确定单元集合体上的场函数。显然，如果单元满足问题的收敛性要求，那么随着缩小单元的尺寸，增加求解区域内单元的数目，解的近似程度将不断改进，近似解最终将收敛于精确解。 用有限元法求解问题的计算步骤比较繁多，其中最主要的计算步骤为： 1）连续体离散化。首先，应根据连续体的形状选择最能完满地描述连续体形状的单元。常见的单元有：杆单元，梁单元，三角形单元，矩形单元，四边形单元，曲边四边形单元，四面体单元，六面体单元以及曲面六面体单元等等。其次，进行单元划分，单元划分完毕后，要将全部单元和结点按一定顺序编号，每个单元所受的荷载均按静力等效原理移植到结点上，并在位移受约束的结点上根据实际情况设置约束条件。 2）单元分析。所谓单元分析，就是建立各个单元的结点位移和结点力之间的关系式。现以三角形单元为例说明单元分析的过程。如图1所示，三角形有三个结点i，j，m。在平面问题中每个结点有两个位移分量u，v和两个结点力分量F x，F y。三个结点共六个结点位移分量可用列
第一性原理简介
第一性原理是什么 第一性原理怎么用 1什么是第一性原理 根据原子核和电子互相作用的原理及其基本运动规律，运用，从具体要求出发，经过一些近似处理后直接求解的算法，称为第一性原理。广义 的第一原理包括两大类，以Hartree-Fock自洽场计算为基础的从头算和 （DFT计算。 从定义可以看出第一性原理涉及到量子力学、、Hartree-Fock自洽场、等许多对我来说很陌生的物理化学定义。因此我通过向师兄请教和上网查资料一点点的了解并学习这些知识。 2第一性原理的作用 以密度泛函理论（DFT）为基础以及在此基础上发展起来的简单而具有一定精度的局域密度近似（LDA）和广义梯度近似（GGA）的第一性原理电子结构计算方法，与传统的解析方法一样，不但能够给出描述体系微观电子特性的物理量如波函数、态密度、费米面、电子间互作用势等，以及在此基础上所得到的体现体系宏观物理特性的参量如结合能、电离能、比热、电导、光电子谱、穆斯堡尔谱等等，而且它还可以帮助人们预言许多新的
物理现象和物理规律。密度泛函计算的一些结果能够与实验直接进行比较一些应用程序的发展乃至商业软件的发布，导致了基于密度泛函理论的第 一原理计算方法的广泛应用。 密度泛函理论（DFT）为第一性原理中的一类，在物理系、化学、材料科学以及其他工程领域中，密度泛函理论（DFT及其计算已经快速发展成 为材料建模模拟的一种“标准工具”。 密度泛函理论可以计算预测固体的晶体结构、晶格参数、能带结构、态密度（DOS、光学性能、磁性能以及原子集合的总能等等。 3第一性原理怎么用 目前我所学到的利用第一性原理的软件为Material Studio 、VASP软件。其中Materials Studio （简称MS是专门为材料科学领域研究者幵发的一款可运行在PC上的模拟软件。使化学及材料科学的研究者们能更方便地建立三维结构模型，并对各种晶体、无定型以及高分子材料的性质及相关过程进行深入的研究。模拟的内容包括了催化剂、聚合物、固体及表面、晶体与衍射、化学反应等材料和化学研究领域的主要课题。 模块简介 Materials Studio 采用了大家非常熟悉的Microsoft标准用户界面， 允许用户通过各种控制面板直接对计算参数和计算结果进行设置和分析。 目前，Materials Studio 软件包括如下功能模块： Materials Visualizer: 提供了搭建分子、晶体及高分子材料结构模型所需要的所有工具，可以操作、观察及分析结构模型，处理图表、表格或文本等形式的数据，并提供软件的基本环境和分析工具以及支持Materials Studio 的其他产品。是Materials Studio 产品系列的核心模块。 Discover: Materials Studio 的分子力学计算引擎。使用多种分子力学和动力学 方法，以仔细推导的力场作为基础，可准确地计算出最低能量构型、分子体系的结构和动力学轨迹等。
有限元分析理论基础
有限元分析概念 有限元法：把求解区域看作由许多小的在节点处相互连接的单元（子域）所构成，其模型给出基本方程的分片（子域）近似解，由于单元（子域）可以被分割成各种形状和大小不同的尺寸，所以它能很好地适应复杂的几何形状、复杂的材料特性和复杂的边界条件 有限元模型：它是真实系统理想化的数学抽象。由一些简单形状的单元组成，单元之间通过节点连接，并承受一定载荷。 有限元分析：是利用数学近似的方法对真实物理系统（几何和载荷工况）进行模拟。并利用简单而又相互作用的元素，即单元，就可以用有限数量的未知量去逼近无限未知量的真实系统。 线弹性有限元是以理想弹性体为研究对象的，所考虑的变形建立在小变形假设的基础上。在这类问题中，材料的应力与应变呈线性关系，满足广义胡克定律；应力与应变也是线性关系，线弹性问题可归结为求解线性方程问题，所以只需要较少的计算时间。如果采用高效的代数方程组求解方法，也有助于降低有限元分析的时间。 线弹性有限元一般包括线弹性静力学分析与线弹性动力学分析两方面。 非线性问题与线弹性问题的区别： 1）非线性问题的方程是非线性的，一般需要迭代求解； 2）非线性问题不能采用叠加原理； 3）非线性问题不总有一致解，有时甚至没有解。 有限元求解非线性问题可分为以下三类：
1）材料非线性问题 材料的应力和应变是非线性的，但应力与应变却很微小，此时应变与位移呈线性关系，这类问题属于材料的非线性问题。由于从理论上还不能提供能普遍接受的本构关系，所以，一般材料的应力与应变之间的非线性关系要基于试验数据，有时非线性材料特性可用数学模型进行模拟，尽管这些模型总有他们的局限性。在工程实际中较为重要的材料非线性问题有：非线性弹性（包括分段线弹性）、弹塑性、粘塑性及蠕变等。 2）几何非线性问题 几何非线性问题是由于位移之间存在非线性关系引起的。 当物体的位移较大时，应变与位移的关系是非线性关系。研究这类问题一般都是假定材料的应力和应变呈线性关系。它包括大位移大应变及大位移小应变问题。如结构的弹性屈曲问题属于大位移小应变问题，橡胶部件形成过程为大应变问题。 3）非线性边界问题 在加工、密封、撞击等问题中，接触和摩擦的作用不可忽视，接触边界属于高度非线性边界。 平时遇到的一些接触问题，如齿轮传动、冲压成型、轧制成型、橡胶减振器、紧配合装配等，当一个结构与另一个结构或外部边界相接触时通常要考虑非线性边界条件。 实际的非线性可能同时出现上述两种或三种非线性问题。
有限元原理与步骤
2.1.1 有限元法基本原理（Basic Theory of FEM） 有限元法的基本思想是离散的概念，它是指假设把弹性连续体分割成数目有限的单元，并认为相邻单元之间仅在节点处相连。根据物体的几何形状特征、载荷特征、边界约束特征等，选择合适的单元类型。这样组成有限的单元集合体并引进等效节点力及节点约束条件，由于节点数目有限，就成为具有有限自由度的有限元计算模型，它替代了原来具有无限多自由度的连续体[24][25]。 有限元法从选择基本未知量的角度来看，可分为三类：位移法、力法和混合法。以节点位移为基本未知量的求解方法称为位移法；以节点力为基本未知量的求解方法称为力法；一部分以节点位移，另一部分以节点力作为基本未知量的求解方法称为混合法。由于位移法通用性强，计算机程序处理简单、方便，成为应用最广泛的一种方法[26]。 有限元法的求解过程简单、方法成熟、计算工作量大，特别适合于计算机计算。再加上它有成熟的大型软件系统支持，避免了人工在连续体上求分析解的数学困难，使其成为一种非常受欢迎的、应用极广泛的数值计算方法[27]。 2.1.2 有限元法基本步骤（Basic Process of FEM） 有限元法求解各种问题一般遵循以下的分析过程和步骤[28][29]： 1. 结构的离散化 结构的离散化是进行有限元法分析的第一步，它是有限元法计算的基础。将结构近似为具有不同有限大小和形状且彼此相连的有限个单元组成的计算模型，习惯上称为有限元网格划分。离散后单元与单元之间利用单元的节点相互连接起来，而单元节点的设置、性质、数目等应视问题的性质、描述变形形态的需要和计算精度而定。所以有限元法分析的结构已不是原有的物体或结构物，而是同种材料的由众多单元以一定方式连接成的离散物体。这样，用有限元分析计算所获得的结果是近似的。显然，单元越小（网格越密）则离散域的近似程度越好，计算结果也越精确，但计算量将增大，因此结构的离散化是有限元法的核心技术之一。有限元离散过程中又一重要环节是单元类型的选择，这应根据被分析结构的几何形状特点、载荷、约束等因素全面考虑。 2. 位移模式的选择 位移模式是表示单元内任意点的位移随位置变化的函数，位移模式的选择是有限元特性分析的第一步。由于多项式的数学运算比较简单、易于处理，所以通常是选用多项式作为位移函数。选择合适的位移函数是有限元分析的关键，它将决定有限元解的性质与近似程度。位移函数的选择一般遵循以下原则（有限元解的收敛条件）：
第一性原理计算原理和方法精编
第一性原理计算原理和 方法精编 Document number：WTT-LKK-GBB-08921-EIGG-22986
第二章 计算方法及其基本原理介绍 化学反应的本质是旧键的断裂和新建的形成，参与成键原子的电子壳层重新组合是导致生成稳定多原子化学键的明显特征。因此阐述化学键的理论应当描写电子壳层的相互作用与重排，借助求解满足适当的Schrodinger 方程的波函数描写分子中电子分布的量子力学，为解决这一问题提供了一般的方法，然而，对于一些实际的体系，不引入一些近似，就不可能求解其Schrodinger 方程。这些近似使一般量子力学方程简化为现代电子计算机可以求解的方程。这些近似和关于分子波函数的方程形成计算量子化学的数学基础。 SCF-MO 方法的基本原理 分子轨道的自洽场计算方 法(SCF-MO)是各种计算方法的理论基础和核心部分，因此在介绍本文计算工作所用方法之 前，有必要对其关键的部分作 一简要阐述。 Schrodinger 方程及一些基本近似 为了后面介绍各种具体在自洽场分子轨道(SCF MO)方法方便，这里将主要阐明用于本文量子化学计算的一些重要的基本 R AB =R 图2-1分子体系的坐标
近似，给出SCF MO 方法的一些基本方程，并对这些方程作简略说明，因为在大量的文献和教材中对这些方程已有系统的推导和阐述[1-5]。 确定任何一个分子的可能稳定状态的电子结构和性质，在非相对论近似下，须求解定态Schrodinger 方程 ''12121212122ψψT p B A q p A p pA A pq AB B A p A A A E R Z r R Z Z M =??????? ?-++?-?-∑∑∑∑∑∑≠≠ （） 其中分子波函数依赖于电子和原子核的坐标，Hamilton 算符包含了电子p 的动能和电子p 与q 的静电排斥算符， ∑∑≠+?-=p q p pq p e r H 12121?2 以及原子核的动能 ∑?-=A A A N M H 2121? 和电子与核的相互作用及核排斥能 ∑∑≠+-=p A B A AB B A pA A eN R Z Z r Z H ,21? 式中Z A 和M A 是原子核A 的电荷和质量，r pq =|r p -r q |，r pA =|r p -R A |和R AB =|R A -R B |分别是电子p 和q 、核A 和电子p 及核A 和B 间的距离（均以原子单位表示之）。上述分子坐标系如图所示。可以用V(R,r)代表－式中所有位能项之和 ∑∑∑-+=≠≠p A pA A B A q p pq AB B A r Z r R Z Z r R V ,12121),( 原子单位
第一性原理
第二章 第一性原理计算方法与软件介绍 19世纪末，科学家们发现经典力学和经典电动力学在描述物质的微观系统时存在明显不足，对实验中的许多现象不能做出真正合理的解释。鉴于此，20世纪初物理学家们在旧量子论的基础上建立了量子力学，主要研究原子、分子、凝聚态物质等内部微观粒子的结构、运动规律等性质，目前已广泛应用于物理、化学、材料等学科领域。随着量子力学理论的不断完善，并结合日趋成熟的计算机技术，量子计算模拟成为了现代科学中必不可少的研究手段之一。第一性原理计算(First-principles calculation)，亦称为从头算(Ab-initio calculation)。该计算方法可根据量子力学基本原理，基于密度泛函理论对材料微观体系的状态和性质进行理论上的预测，且计算过程中不需要使用任何经验参数，只需要一些基本物理量(电子电荷质量e 、电子静止质量m 0、光速c 、普朗克常数h 、波尔兹曼常数k B )。本工作所选用的计算程序为Materials Studio 软件中的CASTEP 量子力学模块，该模块是基于密度泛函理论的从头算量子力学程序。本章节将简要的介绍密度泛函理论和CASTEP 计算模块。 2.1密度泛函理论概述 第一性原理主要的研究对象是多原子体系。它依据量子力学原理，且在无任何实验参数引入的情况下，将多原子体系当作由自由电子和原子核组成的多粒子体系进行处理。然而，关于量子力学中多粒子体系处理的出发点则为著名的薛定谔方程(Schr?dinger Equation)。Schr?dinger 方程是量子力学的一个基本方程，也是第一性原理计算方法的核心，它是由奥地利物理学家薛定谔(Schr?dinger)于1926年提出的。该方程可用于描述微观粒子的运动规律，故亦被称为薛定谔波动方程(Schr?dinger Wave Equation)，其定态方程描述如下： 2 2[()]()(,)2V r r,t i r t t ψψμ?-?+=? (2-1) 式中?为约化普朗克(Plank)常数；μ和V(r)分别表示粒子质量和势场；r 和t 则为体系中所有电子与原子核的位置坐标；Ψ(r,t)是系统波函数，即运动的微观粒子
有限元法基本原理与应用
有限元法基本原理与应用 班级机械2081 姓名方志平 指导老师钟相强 摘要：有限元法的基础是变分原理和加权余量法，其基本求解思想是把计算域划分为有限个互不重叠的单元，在每个单元内，选择一些合适的节点作为求解函数的插值点，将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式，借助于变分原理或加权余量法，将微分方程离散求解。采用不同的权函数和插值函数形式，便构成不同的有限元方法。 关键词：有限元法；变分原理；加权余量法；函数。 Abstract：Finite element method is based on the variational principle and the weighted residual method, the basic idea is to solve the computational domain is divided into a finite number of non-overlapping units, each unit, select some appropriate function for solving the interpolation node points as , the differential variables rewritten or its derivative by the variable value of the selected node interpolation functions consisting of linear expressions, by means of variational principle or weighted residual method, the discrete differential equations to solve. Different forms of weight functions and interpolation functions, it constitutes a different finite element method. Keywords：Finite element method; variational principle; weighted residual method; function。 引言 有限元方法最早应用于结构力学，后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中，把计算域离散剖分为有限个互不重叠且相互连接的单元，在每个单元内选择基函数，用单元基函数的线形组合来逼近单元中的真解，整个计算域上总体的基函数可以看为由每个单元基函数组成的，则整个计算域内的解可以看作是由所有单元上的近似解构成。在河道数值模拟中，常见的有限元计算方法是由变分法和加权余量法发展而来的里兹法和伽辽金法、最小二乘法等。根据所采用的权函数和插值函数的不同，有限元方法也分为多种计算格式。从权函数的选择来说，有配置法、矩量法、最小二乘法和伽辽金法，从计算单元网格的形状来划分，有三角形网格、四边形网格和多边形网格，从插值函数的精度来划分，又分为线性插值函数和高次插值函数等。不同的组合同样构成不同的有限元计算格式。对于权函数，伽辽金(Galerkin)法是将权函数取为逼近函数中的基函数；最小二乘法是令权函数等于余量本身，而内积的极小值则为对代求系数的平方误差最小；在配置法中，先在计
第一性原理计算原理和方法
第二章 计算方法及其基本原理介绍 化学反应的本质是旧键的断裂和新建的形成，参与成键原子的电子壳层重新组合是导致生成稳定多原子化学键的明显特征。因此阐述化学键的理论应当描写电子壳层的相互作用与重排，借助求解满足适当的Schrodinger 方程的波函数描写分子中电子分布的量子力学，为解决这一问题提供了一般的方法，然而，对于一些实际的体系，不引入一些近似， 确定任何一个分子的可能稳定状态的电子结构和性质，在非相对论近似下，须求解定态Schrodinger 方程 ''12121212122 ψψT p B A q p A p pA A pq AB B A p A A A E R Z r R Z Z M =??? ?????-++?-?-∑∑∑∑∑∑≠≠ （2.1） 其中分子波函数依赖于电子和原子核的坐标，Hamilton 算符包含了电子p 的动能和电子p
与q 的静电排斥算符， ∑∑≠+?-=p q p pq p e r H 12121?2 (2.2) 以及原子核的动能 ∑?-=A A A M H 2? (2.3) 和电子与核的相互作用及核排斥能 ∑∑≠+-=p A B A AB B A pA A eN R Z Z r Z H ,21? (2.4) 式中Z A 和M A 是原子核A 的电荷和质量，r pq =|r p -r q |，r pA =|r p -R A |和R AB =|R A -R B |分别是电子p 和q 、核A 和电子p 及核A 和B 间的距离（均以原子单位表示之）。上述分子坐标系如图2.1所示。可以用V(R,r)代表(2.2)－(2.4)式中所有位能项之和 ∑∑∑-+= ≠≠p A pA A B A q p pq AB B A r Z r R Z Z r R V ,1 2121),( (2.5) 原子单位 上述的Schrodinger 方程和Hamilton 算符是以原子单位表示的，这样表示的优点在于简化书写型式和避免不必要的常数重复计算。在原子单位的表示中，长度的原子单位是Bohr 半径 能量是以Hartree 为单位，它定义为相距1Bohr 的两个电子间的库仑排斥作用能 质量则以电子制单位表示之，即定义m e =1 。
第一性原理简介
第一性原理是什么？ 第一性原理有什么用？ 第一性原理怎么用？ 怎样将第一性原理与实 践结合起来？ 什么是第一性原理？1原理，量子力学根据原子核和电子互相作用的原理及 其基本运动规律，运用第一性称为经过一些近似处理后直接求解薛定谔方程的算法，从具体要求出发，计算为基础的从头算。广义的第一原理包括两大类，以
Hartree-Fock自洽场原理DFT）计算。密度泛函理论和（自从定义可以看出第一性原理涉及到量子力学、薛定谔方程、Hartree-Fock因此我通过向师兄密度泛函理论等许多对我来说很陌生的物理化学定义。洽场、请教和上网查资料一点点 的了解并学习这些知识。 2第一性原理的作用为基础以及在此基础上发展起 来的简单而具有一定精(DFT)以密度泛函理论,的第一性原理电子结构计算方法 和广义梯度近似(GGA)度的局域密度近似(LDA)不但能够给出描述体系微观电子特性的物理量如波函与传统的解析方法一样,以及在此基础上所得到的体现体系宏,数、态密度、费米面、电子间互作用势等,穆斯堡尔谱等等比热、电导、观物理特性的参量如结合能、电离能、光电子谱、密度泛函计算的一些而且它还可以帮助人们预言许多新的物理现象和物理规律。． 导致了,结果能够与实验直接进行比较,一些应用程序的发展乃至商业软件的发布基于密度泛函理论的第一原理计算方法的广泛应用。为第一性原理中的一类，在物理系、化学、材料科学以(DFT)密度泛函理论）及其计算已经快速发展成为材料建模DFT及其他工程领域中，密度泛函理论（模拟的一种“标准工具”。密度泛函理论可以计算预测固体的晶体结构、晶格参数、能带结构、态密度（DOS）、 光学性能、磁性能以及原子集合的总能等等。 3第一性原理怎么用？其中ASP、软件。V目前我所学到的利用第一性原理的软件为Material Studio）是专门为材料科学领域研究者开发的一款可运行在MSMaterials Studio（简称使化学及材料科学的研究者们能更方便地建立三维结构模型，上的模拟软件。PC模拟无定型以及高分子材料的性质及相关过程进行深入的研究。并对各种晶体、的内容包括了催化剂、聚合物、固体及表面、晶体与衍射、化学反应等材料和化学研究领域的主要课题。模块简介Materials Studio采用了大家非常熟悉的Microsoft标准用户界面，允许用户通过各种控制面板直接对计算参数和计算结果进行设置和分析。目前，Materials Studio软件包括如下功能模块： Materials Visualizer: 提供了搭建分子、晶体及高分子材料结构模型所需要的所有工具，可以操作、观察及分析结构模型，处理图表、表格或文本等形式的数据，并提供软件的基本环境和分析工具以及支持Materials Studio的其他产品。是Materials Studio产品系列的核心模块。 Discover: Materials Studio的分子力学计算引擎。使用多种分子力学和动力学方法，以仔细推导的力场作为基础，可准确地计算出最低能量构型、分子体系的结构和动力学轨迹等。． COMPASS: 支持对凝聚态材料进行原子水平模拟的功能强大的力场。是第一个由凝聚态性质以及孤立分子的各种从头算和经验数据等参数化并经验证的从头算力场。可以在很大的温度、压力范围内精确地预测孤立体系或凝聚态体系中各种分子的结构、构象、振动以及热物理性质。 Amorphous Cell: 允许对复杂的无定型系统建立有代表性的模型，并对主要性质进行预测。通过观察系统结构和性质之间的关系，可以对分子的一些重要性质有更深入的了解，从
第一节第一性原理计算方法综述
第一性原理计算的理论方法 随着科技的发展，计算机性能也得到了飞速的提高，人们对物理理论的认识也更加的深入，利用计算机模拟对材料进行设计已经成为现代科学研究不可缺少的研究手段。这主要是因为在许多情况下计算机模拟比实验更快、更省，还得意于计算机模拟可以预测一些当前实验水平难以达到的情况。然而在众多的模拟方法中，第一性原理计算凭借其独特的精度和无需经验参数而得到众多研究人员的青睐，成为计算材料学的重要基础和核心计算。本章将介绍第一性原理计算的理论基础，研究方法和ABINIT软件包。 1.1 第一性原理 第一性原理计算( 简称从头计算，the abinitio calculation) ，指 从所要研究的材料的原子组分出发，运用量子力学及其它物理规律，通过自洽计算来确定指定材料的几何结构、电子结构、热力学性质和光学性质等材料物性的方法。基本思想是将多原子构成的实际体系理解成为只有电子和原子核组成的多粒子系统，运用量子力学等最基本的物理原理最大限度的对问题进行”非经验”处理。【1】第一性原理计算就只需要用到五个最基本的物理常量即( m o.e.h.c.k b ) 和元素周期表中各组分元素的电子结构，就可以合理地预测材料的许多物理性质。用第一性原理计算的晶胞大小和实验值相比误差只有几个百分点，其他性质也和实验结果比较吻合，体现了该理论的正确性。
第一性原理计算按照如下三个基本假设把问题简化： 1．利用Born-Oppenheimer 绝热近似把包含原子核和电子的多粒子问题转化为多电子问题。 2．利用密度泛函理论的单电子近似把多电子薛定谔方程简化为比较容易求解的单电子方程。 3．利用自洽迭代法求解单电子方程得到系统基态和其他性质。以下我将简单介绍这些第一性原理计算的理论基础和实现方法：绝热近似、密度泛函理论、局域密度近似(LDA)和广义梯度近似(GGA)、平面波及赝势方法、密度泛函的微扰理论、热力学计算方法和第一性原理计算程序包ABINIT。 1．2量子力学与Born-Oppenheimer 近似固体是由原子核和核外的电子组成的，在原子核与电子之间，电子与电子之间，原子核与原子核之间都存在着相互作用。从物理学的角度来看，固体是一个多体的量子力学体系【2】，相应的体系哈密顿量可以写成如下形式： H (r,R) E H(r ,R) (1-1) 其中r,R 分别代表所有电子坐标的集合、所有原子核坐标的集合。在不计外场作用下，体系的哈密顿量日包括体系所有粒子( 原子核和电子) 的动能和粒子之间的相互作用能，即 H H e H N H e N (1-2) 其中，以是电子部分的哈密顿量，形式为： 22 1 e2 H e(r) r2i 1 e(1-3)
有限元方法理论及其应用
1 课程论文：弹性力学有限元位移法原理（30分） 撰写一篇论文，对有限元位移法的原理作一般性概括和论述。要求论文论及但不限于下列内容：1）弹性力学有限元位移法的基本思想和数学、力学基础；2）有限元法求解的原理和过程，推导计算列式；对基本概念和矩阵符号进行解释和讨论；3）等参单元的概念、原理和应用。 1.1 对一维杆单元有限元形式的理解 我对此提出了几点疑问： 1)为什么边界条件u1=0，就要划去刚度矩阵[K]中对应的行列再解方程？ 2)为什么刚度矩阵[K]会奇异？ 3)为什么平衡方程本身是矛盾的，而加上边界条件u1=0之后就能解出一 个唯一的近似解？ 4)为什么刚度矩阵[K]是对称的？ 下面我谈谈自己的理解:节点平衡方程是在u1不定的前提下，假设单元内位移都是线性变化推导出来的，由此u1相当于一个不确定的定值约束，再加上中间两个节点的连续性要求，系统实际上只有三个独立的自由度（广义坐标）。 对于第一个问题，其实刚度矩阵[K]中的元素不是一成不变的，相反它是伴随边界条件动态变化的。当u1=0时由刚度矩阵的推导过程可以知道，刚度矩阵的第一行和第一列都会变为0，所以此时第一行和第一列对于求解方程是没有作用的。 对于第二个问题，由于系统自由度（广义坐标）只有三个，而我们的方程却列出
了四个，显然
这四个方程不可能线性无关，所以刚度矩阵奇异。 对于第三个问题，首先我们应该明确方程区别于等式，虽然左右两边都是用“=”连接，但是方程只在特殊条件下取得定解。由于平衡方程是在没有约束的条件下推导出来的，显然它不可能满足等式要求。宏观上看，系统在没有外部约束，而又施加有外力，显然系统会产生加速度而绝不会平衡。所以平衡方程本身是矛盾的。而加上边界条件之后，不但满足了平衡的前提，还改变了矩阵的结构和性质，所以有解。但是，由于我们提前假设了位移线性变化，相当于人为对单元施加了额外约束，让位移按照我们假设的规律变化，所以得到的解是过刚的近似解。但对于方程本身而言是精确解。 对于第四个问题，其力学的作用机理类似于作用力与反作用力，由于刚度矩阵不表征方向，所以其大小是相等的。 1.2 有限元法的思想 有限元法是求解连续介质力学问题的数值方法，更一般意义是一种分析结构问题和连续场数学物理问题的数值方法。 有限元法的基本思想是离散化和分片插值。 即把连续的几何机构离散成有限个单元，并在每一个单元中设定有限个节点，从而将连续体看作仅在节点处相连接的一组单元的集合体，同时选定场函数的节点值作为基本未知量并在每一单元中假设一个近似插值函数以表示单元中场函数的分布规律，再建立用于求解节点未知量的有限元方程组，从而将一个连续域中的无限自由度问题转化为离散域中的有限自由度问题。 求解得到节点值后就可以通过设定的插值函数确定单元上以至个集合体上的场函数。对每个单元，选取适当的插值函数，使得该函数在子域内部、在子域分界面上以及子域与外界面上都满足一定的条件。单元组合体在已知外载荷作用下处于平衡状态时，列出一系列以节点、位移为未知量的线性方程组，利用计算机解出节点位移后，再用弹性力学的有关公式，计算出各单元的应力、应变，当各单元小到一定程度，那么它就代表连续体各处的真实情况。
第一节第一性原理计算方法综述
第一性原理计算的理论方法 随着科技的发展，计算机性能也得到了飞速的提高，人们对物理理论的认识也更加的深入，利用计算机模拟对材料进行设计已经成为现代科学研究不可缺少的研究手段。这主要是因为在许多情况下计算机模拟比实验更快、更省，还得意于计算机模拟可以预测一些当前实验水平难以达到的情况。然而在众多的模拟方法中，第一性原理计算凭借其独特的精度和无需经验参数而得到众多研究人员的青睐，成为计算材料学的重要基础和核心计算。本章将介绍第一性原理计算的理论基础，研究方法和ABINIT 软件包。 1.1第一性原理 第一性原理计算(简称从头计算，the abinitio calculation)，指从所要研究的材料的原子组分出发，运用量子力学及其它物理规律，通过自洽计算来确定指定材料的几何结构、电子结构、热力学性质和光学性质等材料物性的方法。基本思想是将多原子构成的实际体系理解成为只有电子和原子核组成的多粒子系统，运用量子力学等最基本的物理原理最大限度的对问题进行”非经验”处理。【1】第一性原理计算就只需要用到五个最基本的物理常量即(b o k c h e m ....)和元素周期表中各组分元素的电子结构，就可以合理地预测材料的许多物理性质。用第一性原理计算的晶胞大小和实验值相比误差只有几个百分点，其他性质也和实验结果比较吻合，体现了该理论的正确性。
第一性原理计算按照如下三个基本假设把问题简化： 1．利用Born-Oppenheimer 绝热近似把包含原子核和电子的多粒子问题转化为多电子问题。 2．利用密度泛函理论的单电子近似把多电子薛定谔方程简化为比较容易求解的单电子方程。 3．利用自洽迭代法求解单电子方程得到系统基态和其他性质。 以下我将简单介绍这些第一性原理计算的理论基础和实现方法：绝热近似、密度泛函理论、局域密度近似(LDA)和广义梯度近似(GGA)、平面波及赝势方法、密度泛函的微扰理论、热力学计算方法和第一性原理计算程序包ABINIT 。 1．2量子力学与Born-Oppenheimer 近似 固体是由原子核和核外的电子组成的，在原子核与电子之间，电子与电子之间，原子核与原子核之间都存在着相互作用。从物理学的角度来看，固体是一个多体的量子力学体系【2】，相应的体系哈密顿量可以写成如下形式： ),(),(R r E R r H H ψψ= （1-1） 其中r,R 分别代表所有电子坐标的集合、所有原子核坐标的集合。在不计外场作用下，体系的哈密顿量日包括体系所有粒子(原子核和电子)的动能和粒子之间的相互作用能，即 N e N e H H H H -++= （1-2） 其中，以是电子部分的哈密顿量，形式为：
第一性原理
第一章引言 在21世纪的今天，全球都面对着资源的短缺和环境的污染这两大问题。氢能源的出现，不仅仅解决了能源短缺的问题（氢能源是二次能源），同时氢能源的使用对环境几乎是没有任何的污染（氢气和氧气的反应产物只有水）。因此，从上个世纪70年代就开始关注氢能源的研发。到21世纪的今天氢能源也逐渐开始走上舞台，但是对于氢能源在应用过程中会出现的问题也亟待解决。 本章内容之一将主要介绍氢能源应用中面临的一个严峻的问题——氢气的 。其二，存储，我们将会详细的论述最新出现的一种储氢材料：储氢合金——AlH 3 简述历年来在实验和理论上对于该材料在常压下的研究成果。同时，提出对于在高压条件下进行研究的必要性以及在现阶段的成果，指出我们理论研究的AlH 3 必要性。最后，将简单的介绍高压物理学在当今学科发展中的重要性以及高压物理的发展历史，当然我们将会简述由于现实实验条件上的限制，高压物理的理论研究对于材料性能的分析和高压物理未来发展方向上的重要性和必然性。 §1.1储氢材料简介 伴随着人类社会的发展和进步，人类赖以生存的环境却让全世界都开始担忧。环境的破坏的危机以及能源的短缺的意识，迫使人们一方面去寻求新的能源，另一方面又要考虑新能源对环境所造成的破坏问题。于是氢能源作为存储量丰富，无公害，无污染的新型能源而得到了全球的关注。在以氢作为能源媒介的氢能体系中，氢的存储和氢的运输成为氢能源的实际应用中的关键环节。近年来，人们注意到储氢合金由于其材料结构上的优势而成为一种新型的储氢功能材料。由于某些合金具备特殊的晶体结构，能够使氢原子很容易的进入晶格间隙中并且形成金属氢化物，由于这种氢与金属的结合力很弱，在加热的时候，氢就能从金属中释放出来。但是这些储氢合金的储氢量很大，可以存储比其自身体积要大上1000-1300倍的氢。目前，对于储氢合金的研究也进行的如火如荼。 1.1.1氢能 随着全球人口急增，人类的能源消耗大幅度的增长；而作为主要能源的煤炭和石油，它们又都是不可再生的能源，其储量极为有限。另外，大量矿物能源的燃烧，造成大气污染、"酸雨"和"温室效应"等环境问题。因此，从20世纪60年代以来，人类为了解决未来能源的供应和生存环境问题，高呼"能源革命"。"
一般有限元原理
一般有限元原理 一、基本理论 有限元单元法是数值计算方法中发展较早、应用最广的一种方法。利用有限元法，可以解决经典的传统的方法难以解决或无法求解的许多实际问题。其优点是部分地考虑边坡岩土体的非均质、不连续的介质特征，考虑岩土体的应力应变特征，可以避免将坡体视为刚体，过于简化边界条件的缺点，能够接近实际从应力应变的角度分析边坡的变形破坏机制。对了解边坡的应力分布及应变位移变化很有利。 有限单元法实质是变分法的一种特殊的有效形式，其基本思想是：把连续体离散化为一系列的连接单元，每个单元内可以任意指定各种不同的力学形态，从而可以在一定程度上更好地模拟地质体的实际情况，特殊的节理元，可以有效地模拟岩土体中的结构面。 在大多数情况下岩土体材料应采用非线形模型，其中包括岩体弹塑性、蠕变、不抗拉特性以及结构面性质的影响。下面简要叙述有限元法的求解过程和原理。 有限单元法的基本原理 1.有限单元法的实施步骤 有限元的重要步骤归纳起来，主要有以下几步： （1）建立离散化的计算模型，包括以一定型式的单元进行离散化，按照求解问题的具体条件确定荷载及边界条件； （2）建立单元的刚度矩阵； （3）由单元刚度矩阵组集总体刚度矩阵，并建立系统的整体方程组； （4）引入边界条件，解方程组，求得节点位移； （5）求各单元的应变、应力及主应力。 2位移模式与单元类型 在一般的有限单元法问题中，我们常以位移作为未知数，称为位移法。为保证解的收敛性，要求位移模式必须满足以下三条： （1）位移模式必须能包含单元的刚体位移。即当节点位移是由某个刚体位移所引起时，弹性体内不会有应变。 （2）位移模式必须能包含单元的常应变，即与位置坐标无关的那部分应变。
(完整版)有限元法的基本原理
第二章有限元法的基本原理 有限元法吸取了有限差分法中的离散处理内核，又继承了变分计算中选择试探函数并对区域积分的合理方法。有限元法的理论基础是加权余量法和变分原理，因此这里首先介绍加权余量法和变分原理。 2.1等效积分形式与加权余量法 加权余量法的原理是基于微分方程等效积分的提法，同时它也是求解线性和非线性微分方程近似解的一种有效方法。在有限元分析中，加权余量法可以被用于建立有限元方程，但加权余量法本身又是一种独立的数值求解方法。 2.1.1 微分方程的等效积分形式 工程或物理学中的许多问题，通常是以未知场函数应满足的微分方程和边界条件的形式提出来的，可以一般地表示为未知函数u 应满足微分方程组 12()()()0A A A ?? ?== ? ??? u u u M （在Ω内） （2-1） 域Ω可以是体积域、面积域等，如图2-1所示。同时未知函数u 还应满足边界条件 12()()()0B B B ?? ?== ? ??? u u u M （在Γ内） （2-2） 要求解的未知函数u 可以是标量场（例如压力或温度），也可以是几个变量组成的向量场（例如位移、应变、应力等）。A ，B 是表示对于独立变量（例如空间坐标、时间坐标等）的微分算子。微分方程数目应和未知场函数的数目相对应，因此，上述微分方程可以是单个的方程，也可以是一组方程。所以在以上两式中采用了矩阵形式。 以二维稳态的热传导方程为例，其控制方程和定解条件如下： ()()()0A k k q x x y y φφφ????=++=???? （在Ω内） （2-3）
0()0q B k q n φφφφφ?-=Γ?=??-=Γ???（在上）（在上） （2-4） 这里φ表示温度（在渗流问题中对应压力）；k 是流度或热传导系数（在渗流问题中对应流度/K μ）；φ和q 是边界上温度和热流的给定值（在渗流问题中分别对应边界上的压力和边界上的流速）；n 是有关边界Γ的外法线方向；q 是源密度（在渗流问题中对应井的产量）。 在上述问题中，若k 和q 只是空间位置的函数时，问题是线性的。若k 和q 是φ及其导数的函数时，问题则是非线性的。 由于微分方程组（2-1）在域Ω中每一点都必须为零，因此就有 1122()(()())0u d v A u v A u d ΩΩ Ω≡++Ω≡? ?T V A L （2-5） 其中 12v V v ?? ?= ? ??? M （2-6） 其中V 是函数向量，它是一组和微分方程个数相等的任意函数。 式（2-5）是与微分方程组（2-1）完全等效的积分形式。我们可以说，若积分方程对于任意的V 都能成立，则微分方程（2-1）必然在域内任一点都得到满足。同理，假如边界条件（2-2）亦同时在边界上每一点都得到满足，对于一组任意函数，下式应当成立 1122 ()(()())0u d v B u v B u d ΓΓΓ≡++Γ≡??VB L 因此积分形式 ()()0u d u d ΓΓ Ω+Γ=??T T V A V B 对于所有的V 和V 都成立是等效于满足微分方程（2-1）和边界条件（2-2）。我们把（2-7）式称为微分方程的等效积分形式。 2.1.2等效积分的“弱”形式 在一般情况下，对（2-7）式进行分部积分得到另一种形式： ()()()()0T T v d v d ΩΓ Ω+Γ=??C D u E F u （2-8） 其中C ，D ，E ，F 是微分算子，它们中所包含的导数的阶数较（2-7）式的低，这样对函数u 只需要求较低阶的连续性就可以了。在（2-8）式中降低连续性要求是以提高V 和V 的连续性要求为代价的，由于原来对V 和V （在（2-7）式中）并无连续性要求，但是适当提高对其连续性的要求并不困难，因为它们是可以选择的已知函数。这种降低对函数u 连续性要求的作法在近似计算中，尤其是在有限单元法中是十分重要的。（2-8）式称为微分方程