《7大晶系的弹性模量矩阵元》由会员分享，可在线阅读，更多相关《7大晶系的弹性模量矩阵元（6页珍藏版）》请在金锄头文库上搜索。
1、关于弹性常数的计算-弹性模量矩阵元2010-06-29 10:44:06| 分类：默认分类|举报|字号 订阅最近又开始辗转计算NH3BH3的弹性常数了。还是VASP软件，首先还是很感谢论坛上的各位前辈：)在计算弹性常数时，还是参考的是侯老师的指南，最近很迷恋这个。下面先转载一篇别人的日志：关键词：第一性原理计算(First-principles calculations ) 刚度矩阵(Stiffness matrix)晶系(Lattice system)1.弹性模量矩阵元对于足够小的变形，由胡克定律(Hookes Law )可知，应力与应变成正比，即 应力分量是应变分量的线性函数，用矩阵的形式
2、可以表示为：T1C11C12C13C14C15C16S1T2C12C22C23C24C25C26S2T3=C13C23C33C34C35C36?S3T4C14C24C34C44C45C46S4T5C15C25C35C45 C55C56S5T6C16C26C36C46C56C66S6式中Cj就是我们通常所说的弹性模量，可以证明，上述刚度矩阵为对称阵，Cj = Cji， 因此，弹性模量的独立张量元数目至多只有21个。晶系的对称性越高，独立的张量元数目就越少。需要指出的是，Cj的数目只与晶系有关，而与晶系中具体 的对称类型无关。下面分别讨论七种不同晶系的弹性模量矩阵元：1.1 三斜晶系(Tricli
3、nic system )三斜晶系是所有七大晶系中对称性最低的晶系，因此拥有最多的独立矩阵元，其形式为：C11C12C13C14C15C16C12C22C23C24C25C26C13C23C33C34C35C36C14C24C34C44C45C46C15C25C35C45C55C56C16C26C36C46C56C66共有21个独立的刚度矩阵元，求解过程也因此较为复杂考虑对称性后，单斜晶系有11个独立的矩阵单元1.2 单斜晶系(Monoclinic system)C11C12C1300C16C12C22C2300C26C13C23C3300C36000C44C450000C45C550C16C2
4、6C3600C66从上式可以看出，正交晶系拥有相当高的对成性，其独立刚度矩阵元的数目为8个。1.4 四方晶系(Tetragonal system )1.4.1 四方晶系(4，-4，4/m)对于具有4，-4，4/m对称操作的四方晶系，其弹性矩阵的形式为：C11C12C1300C16C12C22C2300-C16C13C23C33000000C44000000C440C16-C16000C66其独立刚度矩阵元的数目也为 8个。1.4.2 四方晶系(422, 4mm , -42m , 4/mmm )对于具有422 , 4mm , -42m , 4/mmm对称操作的四方晶系，其弹性矩阵的形式为：C11
5、C12C13000C12C11C13000C13C13C33000000C44000000C44000000C66独立刚度矩阵元的数目仅为6个。1.5 三角晶系(Trigonal system )1.5.1三角晶系(3，3)C11C12C13C1400C12C11C13-C1400C13C13C33000C14-C140C44000000C44C140000C14 (C11-C 12)/2三角晶系(32，3m，32/m )的独立刚度矩阵元的数目为6个1.6 六角晶系(Hexagonal system )C11C12C13000C12C11C13000C13C13C33000000C44000000C44000000(C11-C 12)/2六角晶系共有5个独立的刚度矩阵元1.7 立方晶系(Cubic system )C11C12C12000C12C11C12000C12C12C11000000C44000000C44000000C44立方晶系是所有晶系中对称度最高的晶系，其独立的刚度矩阵元数目仅为 3个 至此，我们列出了所有七大晶系的刚度矩阵元，只要求出各晶系对应的所有独立矩阵元，即可得到晶体的刚度矩阵。
