弹性矩阵课件.ppt

上传人:夺命阿水 文档编号:85482 上传时间:2025-07-10 格式:PPT 页数:115 大小:588KB
下载 相关 举报
弹性矩阵课件.ppt_第1页
第1页 / 共115页
弹性矩阵课件.ppt_第2页
第2页 / 共115页
弹性矩阵课件.ppt_第3页
第3页 / 共115页
弹性矩阵课件.ppt_第4页
第4页 / 共115页
弹性矩阵课件.ppt_第5页
第5页 / 共115页
点击查看更多>>
资源描述

1、晶体的弹性性质晶体的弹性性质应力、应变张量,虎克定律弹性常数与对称性弹性波在晶体中传播压电铁电晶体是电介质,它具有介电性质;同时压电铁电晶体又是弹性介质,它又具有弹性性质,而压电效应就是反映了它的介电性质和弹性性质之间的耦合作用。不同晶体结构的压电铁电晶体,各向异性程度不一样,或者说独立的弹性常数的数目与晶体的对称性有关。形变形变 deformationdeformation在外力作用下,物体的大小和形状都要发生变化,通常称为形变。当讨论物体的转动和移动时,形变对运动物体的影响很小,是次要问题,一般可以忽略不计。当讨论振动的传播或压电效应等问题时,形变就成了重要问题,需要进行深入的讨论和研究。

2、塑性和弹性Plastic and Elastic如果外力撤消后,物体不能恢复原状,这种性质就称为物体的塑性;如果外力撤消后,物体能恢复原状,这种性质就称为物体的弹性。自然界不存在完全的弹性体,也不存在完全的塑性体;只存在既有弹性又有塑性的物体。当外力较小时,形变也小,外力撤消后,形变消失,物体恢复原状;当外力较大时,外力撤消后,物体不能恢复原状。可见物体的弹性有一定的限度,超过这个限度就成为塑性。与压电有关的问题,都属于弹性范围内的问题,应力、应变应力、应变 应变张量:strain tensor晶体中任一点的位置可以用所选定的坐标系的位置矢量来描述,它的三个分量为x1、x2、x3。当晶体发生形

3、变时,其中每一点的位置均会发生改变。设形变前的某一点的位置矢量为r r,形变后为rr(其分量为x1、x2、x3),由于形变这一点的位移可以用位置矢衣来表示:当晶体形变时,晶体内任意两点间的距离都会发生变化,设最近邻的两点形变前的距离为dl l(分量为dxi),形变后的距离为dll(分量为dxi),因为dxi=dxi+di,而 于是:利用以下关系:于是有:最后可得到形变前后距离的变化为:其中张量eik由下式给出:该式给出了在物体形变时,它的长度单元的改变。例如(i/xk),当i=k时,代表伸缩应变(纵向应变),而当ik时,代表切应变(横向应变)。一般称eik为应变张量元。从上式直接可以看出eik

4、eki,即应变张量是对称的。在大多数情况下,应变是很小的,所以上式右方的第三项可以略去,于是应变张量元为:应变张量元的矩阵形式二级对称的张量,有六个独立元素 如果用x、y、z代表位置矢量r r的三个分量;u、v、w代表位移矢量 的三个分量;那么这六个张量元可写成为:应变张量元的几何意义应变张量元的几何意义 正应变体积元的体积改变量:由形变引起的体积相对增量称为体膨胀为:切应变切应变 shearshear由于发生切应变,原来的正方形变成了菱形,它的边长不改变 切应变exy=(v/x+u/y)/2的几何意义 由于切变AA,BB,CC,DD,图中u、v代表A点位移的分量,令AD=AD=x,AB=A

5、B=y,则:所以exy=(v/x+u/y)/2=(1+2)/2 由于应变张量是个对称的二阶张量,只有六个独立的元素,因此常被写成一个纵列矩阵,用S代表张量元,用一个新的足标=1、2、6来代替原来的足标,其对应关系如下:S1S2S3S4S5S应力张量应力张量 stress tensorstress tensor在没有形变的固体中,分子的排列是处于在没有形变的固体中,分子的排列是处于热平衡状态,作用在固体中任意一部分的热平衡状态,作用在固体中任意一部分的合力都等于零。如果固体有形变,那么它合力都等于零。如果固体有形变,那么它就不再处于原来的平衡态,而会受到力的就不再处于原来的平衡态,而会受到力的作

6、用,该力会使物体具有恢复到平衡的趋作用,该力会使物体具有恢复到平衡的趋向。向。这种在固体形变时,作用在固体中单位面这种在固体形变时,作用在固体中单位面积上的力称为应力。应力是一个二级张量,积上的力称为应力。应力是一个二级张量,其各个分量为其各个分量为 xxxx、yyyy、zzzz、yzyz、zyzy、xyxy、yxyx、xzxz、zxzx。为此我们也把应力为此我们也把应力称为应力张量。张量元的前一个足标代表称为应力张量。张量元的前一个足标代表应力的方向,后一个足标代表应力所作用应力的方向,后一个足标代表应力所作用面的法线方向。面的法线方向。作用在立方体上的应力张量元 例如,作用在垂直于x轴的单

7、位面积上沿x方向的应力是xx,这类应力是垂直于表面的,代表张力或压力;作用在垂直于x轴的单位面积上沿y方向的应力是yx,这类应力是沿着表面的,即平行于表面的切向,代表切应力。内应力作用在物体上的总力矩等于零,因此,存在如下关系:应力张量:角标xxyy zzyzzx xy双脚标1122 332331 12单角标123456 单、双脚标之间的对应关系应力张量也是对称的二级张量,它只有六个独立的张量元素。所以对于晶体,也常常把应力张量写成一个纵列矩阵,以T(=1、2、6)来表示,其对应关系为:T1T2T3T4T5T6 xx yy zz yz zx xy应力张量和应变张量的情况有一点不同,当=4、5、

8、6时,T=ij,而S=2cij(ij)。作用在体积元上xyz的力与应力张量元ij之间的关系。如图所示,沿x方向力的分量有三个:所以作用在体积元xyz上力的x分量为:作用在单位体积上力的x量为:同理:以上公式在建立描述固体中弹性波传播的方程时将会用到。胡克定律胡克定律 Hooks lawHooks law对于足够小的形变,应变与应力成正比,因此应变分量是应力分量的线性函数,这一规律称为胡克定律,写成矩阵形式为:弹性柔顺常数弹性柔顺常数 compliancecompliance弹性柔顺常数的物理意义:s11=(S1/T1)Tk,当其它应力分量Tk(k1)为常数(或Tk=0)时,由于沿x方向的伸缩应

9、力T1的改变引起x方向伸缩应变S1的改变,与伸缩应力T1的改变成正比。可见s11只与x方向的伸缩应力T1和伸缩应变S1有关。s12=(S1/T2)Tk,当其它应力分量Tk(k2)为常数(或Tk=0)时,由于沿y方向的伸缩应力T2的改变引起x方向伸缩应变S1的改变,与伸缩应力T2的改变成正比。或s12=(S2/T1)Tk,当其它应力分量Tk(k1)为常数(或Tk=0)时,由于沿x方向的伸缩应力T1的改变引起y方向伸缩应变S2的改变,与伸缩应力T1的改变成正比。可见s12 为与y方向的伸缩应力T2和x方向的伸缩应变S1有关的弹性柔顺常数;或者为与x方向的伸缩应力T1和y方向的伸缩应变S2有关的弹性

10、柔顺常数。s14=(S1/T4)Tk,当其它应力分量Tk(k4)为常数(或Tk=0)时,由于x面上的切应力T4的改变引起x方向伸缩应变S1的改变与切应力T4的改变之比(x面即yz平面)。或s14=(S4/T1)Tk,当其它应力分量Tk(k1)为常数(或Tk=0)时,由于x面上的伸缩力T1的改变引起x面上伸切应变S4的改变与伸缩应力T1的改变之比。可见s14 为与x面上的切应力T4和x方向的伸缩应变S1有关的弹性柔顺常数;或者为与x方向的伸缩应力T1和x面上的切应变S4有关的弹性柔顺常数。s44=(S4/T4)Tk,当其它应力分量Tk(k4)为常数(或Tk=0)时,由于x面上的切应力T4的改变引

11、起x面上切应变S4的改变与切应力T4的改变之比。可见s44 只与x面上的切应力T4和切应变S4有关的弹性柔顺常数。其它弹性常数所代表的意义与s11、s12、s14和s44类似。矩阵形式的胡克定律还可写为 由于应力、应变都是二级张量,所以弹性柔顺常数sij共有36个,其中独立的弹性柔顺常数共21个(如果应力、应变都是9个独立分量,则弹性柔顺常数将为81个)。弹性柔顺常数是一个四级张量。弹性刚度常数弹性刚度常数 stiffnessstiffness胡克定律也可表示为其中cij为弹性刚度常数,共有36个,独立的弹性刚度常数共21个,弹性刚度常数也是一个四级张量。矩阵形式:对于各向同性的多晶体,胡克定

12、律的形式为 式中s44=2(s11-s12),有两个独立的弹性常数;写成分量形式为:从胡克定律看出,各向同性物体的的应力T1、T2、T3只对伸缩形变有贡献,而与切应变无关;切应力T4、T5、T6只对切应变有贡献,而与伸缩形变无关。并且切应力T4只对切应变S4有贡献,与切应变S5、S6无关;T5、T6情况与此相同。实验上常采用杨氏模量Y,泊松比及切变模量G来代替弹性柔顺常数s11、s12、s44。它们之间的关系为:Yangs modulusPoission ratioShear 用杨氏模量Y,泊松比及切变模量表示的各向同性介质的虎克定律:弹性常数与对称性弹性常数与对称性和介电常数一样,材料的弹性

13、常数也与对称性有关。描写各各向向同同性性弹弹性性介介质质的独立弹性常数只有两个;描写完完全全各各向向异异性性介介质质弹性介质的独立弹性常数共有21个;而介于各向同性与完全各向异性之间的介质,它们的独立弹性常数个数则介于221之间。例如,属于三斜晶系1点群的压电晶体是完全各向异性的,独立的弹性常数共有21个。属于立方晶系的23点群和3m点群的压电晶体,是对称性最高的晶体,它接近于完全各向同性。独立的弹性常数只有三个。属于四方晶系4mm点群的BaTiO3晶体,独立弹性常数共有六个。属于六角晶系32点群的-石英晶体和点群的LiNbO3,独立弹性常数都是六个。属于正交晶系mm2点群的铌酸钡钠(Ba2N

14、aNb5O15)晶体和222点群的酒石酸钾钠(NaKC4H4O64H2O)晶体,独立弹性常数有9个。不同对称性晶类的独立弹性常数不同对称性晶类的独立弹性常数根据Neumann原则,晶体的对称性不仅表现在结构上,也表现在它的物理特性上,因此晶体的弹性常数必然和晶体的对称性密切相关。通常是晶体的对称性愈高,其独立的弹性常数分量数目愈少。为了确定晶体具有的独立弹性常数,通常有两种方法:一种是一种是脚标代换法脚标代换法;另一种是另一种是坐标变换法坐标变换法。由于坐标变换法具有普适性因此我们首先讨论弹性常数张量的坐标变换。此外对于各向异性晶体,其弹性常数的数值都是对于正常晶体坐标系给出的,而实际使用的晶

15、片往往是旋转切割的,其坐标选取与正常的晶体坐标系不同,为此必须将弹性常数张量从晶体坐标系变换到实际采用的坐标系中下面首先讨论应力和应变张量的坐标变换。应力张量应力张量X X和应变张量和应变张量x x坐标变换坐标变换设新旧坐标系分别为O-x1,x2,x3和O-x1,x2,x3,如右图示。新旧坐标系的坐标轴的方向余弦为:x1x2x3X1a11a12a13X2a21a22a23X上述新旧坐标系的方向余弦的9个数构成一个正交矩阵因为应力T和应变S是二阶张量,所以它们的坐标变换遵从二阶张量的变换规则。首先考虑应力张量的坐标变换,设T和T分别为坐标变换前后的应力张量,则根据二阶张量的变换法则有 T=AA

16、1 或 T=ATA-1采用爱因斯坦脚标重复自动求和规则,变换前后应力的分量可写成:Tij=aimajnTmn。将上式展开,首先考虑到i=j=1时应力分量T11,有:将上式完全展开,得整理,得:应力矩阵元采用缩写下标,上式可写为:同理可以得出变换后应力张量的其余五个分量上面六组联立代数方程组的矩阵形式式中:应力张量的坐标变换关系矩阵M叫做应力张量的变换矩阵按照完全相同的方法可以得出式中:应变张量S的坐标变换关系矩阵N称为应变张量的变换矩阵注意矩阵注意矩阵M M、N N的的不同!不同!使用完全类似的方法,还可以求出应力张量的逆变换矩阵M-1和应变张量的逆变换矩阵N-1,即T=M-1T,S=N-1S

17、并得出如下关系式:M-1=Nt,N-1=Mt式中,Mt为M转置矩阵,Nt为N的转置矩阵弹性常数张量的坐标变换弹性常数张量的坐标变换设胡克定律在原是坐标中表示为S=cT,T=sS 在坐标变换后新坐标系中表示为S=cT,T=sS 根据新旧坐标变换关系式有S=MS,T=N-1T 故有S=McT=McN-1T将上式和(1)比较得出弹性刚度常数在新旧坐标系中的变换关系为c=McMt (2)同理,可得到弹性柔顺常数在新旧坐标系中的变换关系s=NsNt (3)由式(2)和(3)可知,只要知道了应力和应变张量的坐标变换矩阵M和N,就可以求出弹性常数张量的坐标变换。为了确定晶体独立弹性常数,必须根据晶体的对称性

18、并应用Neumann原则来完成,现在以三角系3m点群晶体为例子来进行讨论。对于三角晶系3m点群的晶系,x=0的面是对称面,z轴为三阶转轴,根据Neumann原则,晶体的弹性常数张量经上述对称性操作,其值不应改变。对于x1=0的对称面,新旧坐标选取如下图,新旧坐标系之间的方向余弦矩阵为:将上式代入应变张量矩阵的坐标变换矩阵N为:将坐标变换矩阵代入弹性柔顺常数在新旧坐标系中的变换式:s=NsNt得,由于x=0面为对称面,新旧坐标系的弹性柔顺常数矩阵应该相等,即sij=sij,为此只有下式成立时才能满足s15=s16=s25=s26=s35=s36=s45=s46=0所以弹性柔顺常数矩阵变成如下:

19、由于z(x3)轴为三阶转轴,新旧坐标系选取如图示,对此新旧坐标系的变换矩阵为:将上式代入坐标变换矩阵N的转置矩阵是将N和Nt代入s=NsNt,再令s=s,得到s11=s22,s13=s23,s14=-s24,s34=0,s44=s55,s56=2s14,s66=2(s11-s12)又因为:由此可见,独立的弹性柔顺常数只有s11,s12,s13,s14,s33,s44共六个所以足标代换法 以222点群为例。222点群表示有三个二次旋转轴,分别沿x,y,z方向。先考察沿z轴的二次旋转轴。因为z轴是二阶轴,当晶体绕z轴转180后,晶体坐标变换为,或由于弹性柔顺常数是一个四阶对称张量,完善的写法应有四

20、个足标。例如,s1111、s1122、s1123、s1212等等,通常为了方便,常用二个足标(缩写下标)代替四个足标。四足标与双足标之间的关系为112233231312123456例如:s1111s11,s1122s12,s1123s14,s1212s66等。当晶体绕z轴转180后,四足标中的变换为,1111,2222,3333,23-23,13-13,1212;双足标中的变换为11,22,33,4-4,5-5,66。11,22,33,4-4,5-5,绕z轴转180后:绕x、y轴转180后:因为z轴是二阶轴,当晶体绕z轴转180后,弹性柔顺常数应保持不变,这就要求sij=sij,i,j=1,2

21、3,4,5,6可见,只有当s14=s15=s16=s34=s35=s46=s56=0时两者才完全一致。再利用222点群晶体的轴是二阶轴,重复上述方法,可以得到s16=s26=s36=s45=0。最后得到属于正交晶系222点群晶体的弹性柔顺常数如下。其中独立弹性柔顺常数共9个。晶类 弹性柔顺常数矩阵 弹性刚度常数矩阵 三斜:1点群 c的矩阵与s的矩阵完全相似 单斜:2点群m点群 c的矩阵与s的矩阵完全相似 晶类弹性柔顺常数矩阵 弹性刚度常数矩阵 正交:222,2mmc的矩阵与s的矩阵完全相似 四方:4 c的矩阵与s的矩阵完全相似 四方:422,4mm,2m c的矩阵与s的矩阵完全相似 三方:3

22、点群 除去c46=c25,c56=c14,c66=(c11-c12)/2外,其余与s的矩阵相似 三方32点群3m点群:除去c56=c14,c66=(c11-c12)/2外,其余与s的矩阵相似 六方:6 点群 点群622 点群6mm 点群m2 点群 除去c66=(c11-c12)/2外,其余与s的矩阵相似 立方:23 点群3m 点群 c的矩阵与s的矩阵完全相似 晶类 弹性柔顺常数矩阵 弹性刚度常数矩阵 压电陶瓷 除去c66=(c11-c12)/2外,其余与s的矩阵相似 各向同性介质 除去c44=c55=c66=(c11-c12)/2外,其余与s的矩阵相似 压电陶瓷的弹性常数 先看弹性柔顺常数s1

23、1、s22、s33。因为s11反映压电陶瓷沿x方向伸缩形变的弹性柔顺常数,它只与x方向的弹性性质有关。同理,s22、s33是分别是反映压电陶瓷沿y方向和z 方向伸缩形变的弹性柔顺常数,它们分别只与y方向和z方向的弹性性质有关。而经过极化处理后的压电陶瓷(设轴z轴为极化轴),沿z方向的性质就与沿x、y轴方向的性质就不一样了。所以s11s33、s22s33,但是对于x轴与y轴之间则没有任何差别,即xy平面是各向同性面。这就要求弹性柔顺常数s11=s22,于是得到s11=s22s33的结论。s12、s13、s23等弹性柔顺常数。因为任何弹性介质,当纵向伸长时,横向都要产生收缩。而s12、s13、s2

24、3就是反映这一性质的弹性柔顺常数,所以s12、s13、s23都不可能等于零。由于压电陶瓷的xy平面是各向同性面,当其x轴与y轴之间互换时,弹性常数应保持不变,这就要求s13=s23。但是压电陶瓷的x轴与z轴之间不能互换,这就要求s12s23,于是得到s12s23=s23。此外还有s14=s24,s15=s25,s16=s26。s44、s55、s66等弹性柔顺常数又因为xy平面的各向同性,当x轴与y轴之间互换时,弹性柔顺常数应保持不变,即,因为:而x轴、y轴与z轴之间不能互换,固有s44=s55s66。s14、s15、s16、s24、s25、s26、s34、s35、s36、s45、s46、s56

25、等弹性柔顺常数,可用足标代换法证明它们皆等于零。最后得到压电陶瓷的弹性柔顺常数用矩阵表示为:其中独立弹性柔顺常数为s11、s12、s13、s33、s44五个。用上述相同的方法,可得到压电陶瓷的弹性刚度常数用矩阵表示为:其中独立弹性柔顺常数为c11、c12、c13、c33、c44五个。晶体中的弹性波晶体中的弹性波 以平面波为例,简单介绍描述压电体中声波传播的方法对于各向异性的晶体,它的弹性性质由胡克定律来描述。压电器件中有一类重要和应用广泛的领域:声表面波器件(Surface Acoustic Wave Devices,SAWD),电声电而应变张量元为:式中u、v、w代表晶体中质点位移沿x、y、

26、z方向的分量。式中代表晶体密度,按(3-8)对应关系,上式可改写为质点的运动方程可写为质点的运动方程可写为 把(3-12)、(3-18)代入上式,原则上就得出波的传播方程式,从而解出位移三个分量u、v、w。u(t),v(t),w(t)显然,这样的运算比较繁杂,下面我们介绍另外的方法求解晶体中所传播的弹性波方程的方法。讨论的是各向异性介质,波在不同方向上的传播情况显然是不同的。例如,各个方向上的有效弹性常数不同,波的传播速度也不同。下面我们具体来考虑某一方向上的传播情况。设任意传播方向,它的方向余弦为l、m、n,在这个方向上的某点P(x、y、z)同原点的距离为 ,则显然,P的位移分量u、v、w依

27、赖于(即质点的位移)波的传播方向l、m、n,所以可以把(3-18)式改写为:应变张量元(3-18)式改写为:应力和应变之间由虎克定律联系应力和应变之间由虎克定律联系(3-21)式代入上式得到质点位移(u,v,w)与应力(T1-T6)之间的关系质点运动方程可写为:式中ij称为克利斯托夫(Christoffel)张量,它共有六个张量元,组成对称二级张量。克利斯托夫(Christoffel)张量ij的具体表示式为:-(3-23)设 表示沿传播的波在晶体中所引起的弹性位移矢量,它的分量为u、v、w。位移矢量 的方向余弦为p、q、r,那么,为位移矢量 的长度。把的表达式代入到质点运动方程(3-22)式中

28、得到,这就是 沿方向传播的弹性波方程,式中ceff代表沿方向晶体的有效弹性常数。晶体沿方向 的有效弹性常数ceff须满足下列方程组:因而得久期方程 由 此 可 知,一 般 情 况 下,ceff有 三 个 解ceff,i(i=1、2、3),它们分别对应于三个不同的波,其对应的传播速度为(ceff,i/)1/2。(3-26)对于这三个波,质点分别有相应的三个位移,把由久期方程(3-26)式求得的ceff,i代入到有效弹性常数(3-25)式,即可解出与ceff,i对应的位移的方向余弦pi、qi、ri。下面以立方晶系为例,对上面的讨论加以说明。立方晶系的弹性刚度常数立方晶系的弹性刚度常数立方晶系只有三

29、个独立的弹性常数c11、c12和c44。由(3-23)式可得到它的克利斯托夫张量为(3-27)把上式代入到久期方程(3-26)式,得到(3-28)当 沿100方向时,l=1、m=0、n=0,这时上式(3-28)变为 即可解得Ceff,1=C11,Ceff,2=C44,Ceff,3=C44。把这些式子代入到(3-24)式,可得到三个弹性波的波速和对应质点的位移方向。沿100方向三个弹性波的波速和对应质点的位移方向(1)即沿100方向;(2)即沿010方向;(3)即沿001方向;当 沿110方向时,l=m=1/、n=0,同理可知也有三个弹性波:(1)即沿110方向;(2)即沿001方向;(3)即沿

30、 方向;当 沿111方向时,l=m=n=1/,三个弹性波为:(1)即沿111方向;(2)和(3)即沿垂直于111方向平面上的任意方向;从以上讨论可以知道,弹性波在晶体中传播一般有三个不同速度,其中一个波的位移方向和波矢方向 相同,我们把它称为纵波;而另外两个波的位移方向垂直于波矢方向,则称为横波。最后,我们还应指出一点,对于各向异性介质(如单晶)来说,在一般的方向上,虽有三个传播速度,但是对应于某一传播速度,往往质点位移方向既不平行也不垂直于波矢方向,因此纵波和横波不能加以分辨,只能在某些特殊方向上,这种分辨才有可能。小结应力张量、应变张量,虎克定律,弹性柔顺、刚度常数对称性,弹性常数数目弹性波在固体中的传播

展开阅读全文
相关资源
猜你喜欢
相关搜索

当前位置:首页 > 高等教育 > 大学课件

宁ICP备18001539号-1