靳静：压电材料平面问题的基本解在复变函数中的应用 3基本公式推导 4无限介质平面问题复势函数基本解 考虑一横观各向同性压电体，设其各向同性平面平行 如图1所示，坪面内%点受任意集中力X+Y和 于一平面，且沿轴方向处于平面应变状态，则对其x- 集中点电荷Q昀为单位厚度上的载单位厚度上的载荷) 平面内的二维问题，其面内应力场、位移场和电位移场可 作用，此时复势函数可表示为： 表示为三个复势函数Pk(,(k=123多下同)的线性组 Pk()=A四-面) (15) 合，其中Pk()由下列边界条件确定 其中=音十“kA为常复数，确定如下： 2Re Pk()=- 3 2Re 9k()三 ds (8) 入kPk() D ds 3 2RR()=u 图1无限压电介质含集中荷载 49k()=v (9) 把式(15分别代入式(12).(13)，并计算绕名点旋转 一周时相关量的增量。由位移单值条件知，位移增量为零： 2Re K9k()= 由力平衡条件知，外力主矢增量为一(X+Y):由高斯定 其中“R"表示取实部：X和y分别为边界上的外力和 理知 中=Q由电场环流定理知。 年d-0:而ra- 位移的直角坐标分量：D和E分别为电位移的法向分量和 Z0)的增量为2x由此可得确定A的代数方程： 电场的切向分量：为边界上的弧长：其余为已知常数 现引入下列矩阵： 《-[(=(多 (16) 11 R B (-[-(0 [= [= 其中 kk -(44月 (》=(y-x 同时式(16又可展开为： 4A=1月 Ak号 (17) 2π台 当名点作用逆时针转向集中力偶M时，可证： 9k()=M(-6) (18) 其中M=M(AM+iMe)(1-ek), 其中表示矩阵转置，则式(8，(9汉可表示为： [政9+[(9=〈 (10) 参考文就 [以9+比(g=( 【刂丁皓江，王国庆，陈伟球.压电材料平面问题的基本解[手 (11) 中国科学，199727(3：224-228 可证当k不相等时，[日，[均可逆，故由式(10，(11) 【2]金吞根，等.压电材料平面问题的一般解及其在Te法 中的应用[」.自然科学进展，200212(9,713一719 [3引丁皓江，王国庆，梁剑.压电介质平面问题的一般解和基本 (12) 解[」.力学学报.19968(4,441-448 9+K9=[以 13) [4丁皓江，等.压电介质半平面的Gren晒数解[.中国科 其中[小=[∧h=[-,[d=[G闭=[- 学，199727(5,392-396 [B=[0-[a.[i=[-[ 【5高寸法崔得密.压电介质平面问题的基本解[」.应用力 学学报。199916(1)：140-143 式(12)又可展开为 【6克劳斯.电磁学1M.北京：人民邮电出版社.1979 ()+PK=月 [7刀M uskhelish ivilN数学弹性力学的几个基本向题IM. (14) 北京：科学出版社1953 ?1994-2016 China Academic Journal Electronic Publishing House.All rights reserved. http://www.dki.net3 基本公式推导 考虑一横观各向同性压电体 , 设其各向同性平面平行 于 x-z平面 , 且沿 z轴方向处于平面应变状态 , 则对其 x￾y平面内的二维问题 ,其面内应力场 、位移场和电位移场可 表示为三个复势函数 φk(zk), (k=1, 2, 3;下同 )的线性组 合 , 其中 φk(zk)由下列边界条件确定 2Re∑ 3 k=1 φk(zk)=- ∫ s 0 Ynds 2Re∑ 3 k=1 μkφk(zk)=- ∫ s 0 Xnds 2Re∑ 3 k=1 λkφk(zk)=- ∫ s 0 Dnds (8) 2Re∑ 3 k=1 pk(zk)=u 2Re∑ 3 k=1 qkφk(zk)=v 2Re∑ 3 k=1 Kkφk(zk)= ∫ s 0 Efds (9) 其中 “Re”表示取实部 ;Xn, Yn和 u, v分别为边界上的外力和 位移的直角坐标分量 ;Dn和 Ef分别为电位移的法向分量和 电场的切向分量 ;s为边界上的弧长 ;其余为已知常数 。 现引入下列矩阵 : H = 1 1 1 μ1 μ2 μ3 λ1 λ2 λ3 U = p1 p2 p3 q1 q2 q3 k1 k2 k3 fa = - ∫ s 0 Ynds, ∫ s 0 Xnds, ∫ s 0 Dnds T fb = u, v, ∫ s 0 Eids T φ = φ1 , φ2 , φ3 T 其中 T表示矩阵转置 ,则式 (8), (9)又可表示为 : H φ + H φ = fa (10) U φ + U φ = fb (11) 可证当 μk不相等时 , H , U均可逆 , 故由式 (10), (11) 得 φ + B φ = ∧ fa (12) φ + F φ = G fa (13) 其中 ∧ = ∧ kj = H -1 , G = Gkj = U -1 B = H -1 H , F = U -1 U 式 (12)又可展开为 φk(zk)+∑ 3 j=1 Bkjφj(zj)=∑ 3 j=1 ∧ kjfαj (14) 4 无限介质平面问题复势函数基本解 如图 1所示 , x-y平面内 Z0点受任意集中力 X+iY和 集中点电荷 Q(均为单位厚度上的载单位厚度上的载荷 ) 作用 ,此时复势函数可表示为 : φk(zk)=Akln(zk -zk0) (15) 其中 zk0 =x0 +μky0;Ak为常复数 , 确定如下 : 图 1 无限压电介质含集中荷载 把式 (15)分别代入式 (12)、(13),并计算绕 Z0 点旋转 一周时相关量的增量 。由位移单值条件知 , 位移增量为零 ; 由力平衡条件知 , 外力主矢增量为 -(X+iY);由高斯定 理知 n∮ds=Q;由电场环流定理知 , ∮Etds=0 ;而 ln(Zk - Zk0)的增量为 2πi, 由此可得确定 Ak的代数方程 : A - B A = 1 2πi ∧ f 0 j A - F A = 0 (16) 其中 A = A1 , A2 , A3 T fj 0 = Y, -X, Q T 同时式 (16)又可展开为 : Ak -∑ 3 j=1 BkjAj = 1 2πi∑ 3 j=1 ∧ kjfj 0 (17) 当 Z0 点作用逆时针转向集中力偶 M时 ,可证 : φk(zk)=cMk (zk -zk0) -1 (18) 其中 cMk =iM(∧ k1 +i∧ k2 )(1 -iμk)/8π。 参考文献: [ 1] 丁皓江, 王国庆, 陈伟球.压电材料平面问题的基本解 [ J] . 中国科学, 1997, 27(3):224-228. [ 2]金吾根, 等.压电材料平面问题的一般解及其在 Trefftz法 中的应用[ J] .自然科学进展, 2002, 12(9):713 -719. [ 3]丁皓江, 王国庆, 梁剑.压电介质平面问题的一般解和基本 解[ J] .力学学报, 1996, 28(4):441 -448. [ 4]丁皓江, 等.压电介质半平面的 Green函数解[ J] .中国科 学, 1997, 27(5), 392-396. [ 5]高寸法, 崔得密.压电介质平面问题的基本解 [ J] .应用力 学学报, 1999, 16(1):140-143. [ 6]克劳斯.电磁学[ M] .北京:人民邮电出版社, 1979. [ 7] MuskhelishiviliNI.数学弹性力学的几个基本问题[ M] . 北京:科学出版社, 1953. · 117· 靳 静:压电材料平面问题的基本解在复变函数中的应用