高阶弹性常数，研究报道的最高弹性常数为6阶弹性常数，大部分超硬材料的高阶弹性常数仍然是未知的。

• 中文名称 高阶弹性常数
• 外文名称 Highorder elastic constants
• 最高 6阶弹性常数
• 现状 大部分超硬材料的仍然是未知

导论

来自　　众所周知的是二阶弹性常数是一个二阶四秩张量(cijkl)，三阶弹性常数是三阶六秩张量(Cijklmn);更高阶的弹城肉皮盐言细音微性常数还包括四阶弹性常数(Cijklmnpq)，五阶弹性常数(Cijklmnpqrs)和六阶弹性常数(Cijklmnpqrsuv)，分别是八秩，十秩和1输激款刻2秩张量。研究报道的最高弹性常数为6阶弹性常数。二阶弹性常数的计算由于张量分量较少，同时在计算过程中主要是采用了线性hooker定理，所施加的应变量很小，因此在材料力学性能表征中得到了广泛的应用，三阶弹性常数描述了非线性Hooker定理或者非线性力作用下的材料的力学响应问题360百科。

　　三阶弹性常数矩阵形式蒸品夫粉杀校干顶左十分复杂，即使对于立方晶体结构也有6个独立分量，对于对称性更低的晶体结构则分实室总占量更多，如果采用能量-应变的方法来计算各个独立分量，则计算量相当可观。如立方晶体(点群O装什吧甲势阻事信娘图规h，O，Td)有六个独立分量，则需要六个不同的应变模式，得到六个多元一次方程组，联立求解得到各个分量数值。若采用应力-应变(Strain-Stress relations)则可以明显减少应变模式的数量，但主要问题在于这要求第一原理计算软件具有计算晶体Cauchy应力张量的能力，然而，广泛采用的DFT室合名圆施没级计算软件，如VASP，Wine2K等等是不能直接得到Cauchy应力张伯还补乙县乡还量的，只能计算特定应变下的应变能数值。Materials Studio软件中的CASTEP模块是为数不多的具有直接计算应变结构Cauchy应力张量的软件，因此CASTEP模块在计算材料的二阶弹性常数方面十分的方便。此外，对于三阶弹性常数，没有盟火收较式软件可以直接计算，需要研究者自行设计方法进行计算。

　　三阶弹性常数可以描述材料在高压下的力学响应情况，鉴于高压物理学维货现货州伤，行星结构科学等相关领域的飞速发展，对于超硬材料的非线性力学常数受到了越来越多的关注。采用超声腔共振法可以方便的测量材料的二阶弹性常数，但对于非线性弹性常数，试验方面进展十分的缓慢，时至今日，大部分超硬材料的高阶弹黄做执苗玉氢伯良析双性常数仍然是未知的。

　　器数高阶弹性常数 - 1 广义Hoo训算运许探位福体曾集盟ker定理以及晶体的Neumann原理

装呀歌表误要　　(The gen针歌额打片娘纸手给杂eralized Hooker's Law and Neumann Principle)

　　1.1 广义Hooker定理

　　(The Generalized Hooker's Law)

　　与传统线弹性力学中广泛采用的Hooker定理相比，广义Hooker定理可以认为是包含了高阶神己殖负王武的零非线性应力和应变关系项的Hooker定理的Taylor级质吧称史增象屋我数展开形式，因此从数学意义上来讲这种应力对应变的阶数可以无限制的进行，正如前文所说，广泛景岁采用的弹性常数为应力对应变展开的线性项，展开系数即为Cijkl，Cijkl就是弹性常数，由于根据张量运算法则可知Cijkl是一格机右个四秩张量，描述两个二秩张量应力(stress)和应变(Strain)之间的关系。Cijkl矩阵元素有3^4个，考虑到Lagrange应变为对称矩阵，同时应变自由能与应变路径无关，可以用一个6×6的矩阵来描述，有36个矩阵元素。进一步的矩阵元素化简来自于晶体结构点群对称性对物理学性质的限制，即Neumann原理。例如对于三斜晶体(Triclinic Crystal Class)，独立矩阵元素为21个，如果是正交晶体(Orthorhombic Crystal Class)则减少到9个，对称性最高的立方晶系(Cubic Crystal Class)，仅为三个(C11， C12和C44)。

　　简单归纳一下Cijklmnpq....的性质:

　　Cijkl: Second order 4th rank tesnor, 6×6 matrix with 36 matrix elements and reduces to 21 independent elements for Triclinic crystal class (The matrix form see Fig);

　　Cijklmn: Third order 6th rank tensor, 21×6 matrix with 126 elements and for triclinic crystal class the independent numer is 56;

　　Cijklmnpq: Fourth order 8th rank tensor, 56×6 matrix with 336 elements and for cubic crystal class the independent nuber is 11;

　　Using Voigt notations, we have Cijkl-Cij， Cijklmn-Cijk and Cijklmnpq-Cijkl;

　　Voigt Notations:

　　examples:C1111-C11， C111111-C111 and C11111111-C1111;

　　C1122-C12;C112233-C123; C11112233-C1123;C2332=C2323=C44 and etc. etc.

　　1122 3323(32) 13(31) 12(21)

　　12 3 45 6

　　1.2 Neumann 原理

　　(Neumann Principle)

　　Neumann原理指出，任何晶体结构的物理性质所具有的对称性不低于晶体点群的对称性，这表明张量分身最少具有晶体点群对称性，将晶体点群对称操作作用到各个张量分量Cijkl上，得到新的张量Cmnpq ，则Cijkl=Cmnpq ，即张量各个分量相等。对于三阶和四阶弹性常数则可以表示为

　　对于高秩张量，上述计算必须在计算机上编程进行，否则计算量是十分巨大的。对于晶体结构，如果张量采用了Voigt标记法之后，张量在晶体点群对称操作下的变化相当于对张量下标进行变换，对Voigt标记数字进行对称变换得到的下标分量的对应关系与对张量本身变换得到的分量关系是一致的。存在的问题是晶体点群包含的对称操作的个数很多，如果将每个点群元素对张量进行操作，则不能有效的减少运算量，并且这样的操作大多数是重复的。根据点群的性质，对于每个点群总能找到几个代表性的操作元素(Representative Operations)，点群中所有的对称操作元素可以通过代表性元素的幂乘积得到，因此代表性元素称为点群的生成操作(generators)。因此只要找到点群的生成操作，将生成操作对张量进行变换即可得到点群对称性限制下独立矩阵元素。对于32个空间群点群的生成操作在大部分群论方面的书籍中均可以找到，当然也可以借助群论分析软件Isobyu获得。

　　简化后的二阶，三阶以及四阶弹性常数矩阵形式依然十分复杂，这里给出几个例子:

　　examples: Orthorhombic crystal class: Space Group Pmm2; Point group mm2

　　Cijkl (Cij)

　　Cijklmn (Cijk)

　　Cijklmnpq (Cijkl)

　　Cubic Crystal Class: Space Group P-43m; Point group -43m

　　Cij

　　Cijk

　　Cijkl

　　(下文未注明时弹性常数均为Voigt notation， Cijkl-Cijklmnpq)

　　高阶弹性常数 - 2 应变自由能与弹性常数的关系

　　(The Relationship between Strain energy and Elastic Constants)

　　应变自由能(Strain energy)定义为应变晶胞与初始平衡结构总能量的差值，应变自由能与弹性常数不能直接对应。实际上如果将应变自由能除以初始晶胞的晶体得到新的参数称为应变自由能密度 (Strain energy density)，该参数对Lagrange应变的Taylor级数展开直接对应弹性常数的某个组合或者弹性常数本身。在应变形式比较复杂的情况下，应变自由能密度展开系数一般为弹性常数的线性组合形式，因此如果要通过应变自由能密度来计算弹性常数则需要求解一组多元一次方程组，方程组的个数即为独立弹性常数的个数。如对于四阶弹性常数Cijkl对于正交晶体有42个独立变量，则至少应该有42个多元一次方程组，同理立方晶体至少为11个。

常数计算

　　(The calculations for Higher order elastic constants)

实验测量

　　(Experimental determination)

　　高阶弹性常数的试验测量主要是采用超声腔共振法(Ultrasonic Resonance Method)，如对于二阶弹性常数Cij，可以直接采用该方法，设置不同波长的纵波(P Wave)和横波(S Wave)，其中纵波对应晶体结构中的疏密变化，即拉伸模量数值，如C11，C22，C33等;而S波则对应剪切变形即C44，C55，C66等。在施加混合波形情况下可以得到C12，C13等混合应变模式下的弹性常数数值。对于C11而言可以施加一个沿[100]方向传播的P波得到，具体关系为Cijkl=pv^2(p: density; v: Velocity)。对于三阶弹性常数，可以通过测量二阶弹性常数与温度之间的关系得到，同时注意到Cijk是Cij对应变一阶导数，因此也可以通过对晶体预先施加特定的应变(应力模式)，然后测定声速的变化来对曲线进行拟合得到。更高阶的弹性常数如Cijkl则需要通过计算Cijk在特定应变下的波速变化计算，由于高阶弹性常数独立变量很多，因此试验测定高阶弹性常数遇到的极大的困难，到位置大部分晶体的三阶以上弹性常数都是未知的。

理论计算

　　(Theoretical Calculations)

　　由于计算机的飞速发展和固体理论，凝聚态物理理论方面的巨大进度。高阶弹性常数可以通过理论计算获得，主要是基于密度泛函理论的第一原理计算(The ab-initio calculations based on density functional theory).正如前文提到的，计算方法有两个:应力-应变关系和应变能-应变关系。前者基于广义Hooker定理，后者基于应变能密度的Taylor级数展开。两个方法各有优势，Stress-Strain方法最大的特点是所需要计算的应变模式很少就能计算出全部的弹性常数，如Cubic晶体，只需要一个应变模式即可得到C11，C12和C44，但若采用应变能密度-应变关系则需要设计三个应变模式分别计算获得三元一次方程组，求解得到全部常数，应变-应力计算主要需要计算出晶体的Cauchy应力张量，但由于晶体Cauchy应力张量对边界条件敏感，因此如果精度计算方法对拟合精度都有严重影响;应变能密度-应变关系特点是计算量很大，但数值精度稳定，计算结果比较可靠。

参考文献

　　(Refe来自rences)

　　1. Hiroshi Soma, Yosio Hiki, Journal of the Physical Society of Japan, 37 (2), 544 (1974)

　　2. Shigeru Mori, Yosio Hiki, Journal of the Physical Soc360百科iety of Japan, 45 (5),1449 (1978)

　　3. Yosio Hiki, Annual Review of Materials Science, 11, 51-73 (1981)

　　4. Xanthi基下善生满七歌错黄ppi Markenscoff, Journal of Applied Physics, 50 (3), 1325 (1979) (The 4th order elastic constants for hexagonal and trigonal crystal classes are discusssed)

　　5. P.B. Ghate, Journal of Applied Physics, 35 (2), 337 (1964)

　　6. T.S.G. Krishnamurty, Acta Crystallgraphy,16, 839 (1963)

　　7.REMI策肉评江晚句蒸Brendel, Acta Crystallgraphy, 唱关群陆A35, 525-533 (1979) (The independent fourth order elastic constants for all crystal classes are listed)

　　8. David Y. Chung and Yuan Li, Acta复互干血. Cryst, A30, 1 (1974巴百) (The 4th, 5th and 6th order tensors are discussed in this paper and the matrix forms of them 酸良单续粉概声才粉套席are also tabulated, very good article!)

　　9. Lars Fa任st and et.al., Physical Revi钱ew B, 51 (24), 17431 (1995) (The Second order elastic constants for Hexagonal crystals were studied)

　　10. S.K.R. Patil and et.al., Physical Review B, 73, 104118 (2006纸居意呀香员未汽脸细) (The second order elastic constants for Orthorho弦得回引宗往卷众穿mbic and Tet等供ragonal, Cubic crystals)

　　11. O.H. Nielsen, Physical Review B, 34 (8), 5808 (1986) (both of the Sec战吸器然ond and thir满d order elastic 引血宗模胶向皮被流constants of Diamond were calculated)

　　12. Jijun Zhao and el.al., Physical Review B, 75半确弦啊谓伯油东众, 094105 (2007) (The general method for the calculation 响轮没照调他of higher ord补给回候混er elastic constants using density functional the