电致伸缩材料力学问题的位移函数解法.docVIP

电致伸缩材料力学问题的位移函数解法 蒋 泉1* 匡震邦2 高存法1 1南京航空航天大学，航空宇航学院，江苏南京（210016） 2上海交通大学， 船舶海洋与建筑工程学院，上海 （200240） 摘要：本文在Knops、Smith和Warren等人在电致伸缩问题研究的基础上，对电致伸缩问题基本方程作了一些修正。将Stratton、Landau和Lifshitz导出的电致伸缩体积力引入电致伸缩力学问题中，得出了较为合理的电致伸缩基本方程。并通过构造位移函数，建立了电致伸缩力学问题的位移解法。 关键词：电致伸缩材料 应力场 位移函数 1、引言 在机－电耦合场中，电致伸缩是一个关于电场二次项的非线性本构关系，有着一定的复杂性。目前流行的文献中，一般应用Stratton[1]和Landau、Lifshitz[2]导出的电介质在电场中受到的应力张量作为本构关系，或者在此应力张量的基础上加上一些跟体积力有关的电场力（例如表征体积力的应力张量）。Knops[3]在1963年论文中给出了本构关系包含有前面所讲的表征体积力张量。随后Smith和Warren[4]所用的本构关系也是照搬Knops的本构关系。另外其他的很多作者也给出了一些别的形式的本构方程。例如Yang和Suo[5]等人在论文中提出的低电场中的二次电场强度以及高电场中电致伸缩达到饱和状态的阶跃式电致伸缩本构关系。同时他们对电致伸缩材料中的传导缺陷进行了一定的研究。匡[7]通过引入热力学内变量提出了较为一般的力电耦合场中材料的本构方程。从而通过一定的简化，忽略电场的一次项以及与温度等条件有关的项，直接导出了较为正确的关于电致伸缩材料的本构关系。修正了前面的一些错误，重构了在力电耦合场中的边界条件。本论文的一些理论和计算都是基于这些结果得出的。 电致伸缩材料力电方面的具体解法是Knops提出的，其应用力学基本方程通过构造应力函数，并发展了Myskhelisvili[6]对弹性体的解法，利用复变函数方法给出了标准的电致伸缩力学问题的解法。最后利用两个应力函数给出了平面应变情况下电致伸缩问题中的应力、位移以及边界条件的具体表达式。随后Smith基于Knops的工作，给出了更为详尽的电致伸缩的解法，并对平面应变和平面应力问题作出了系数上不同的区分，并通过含有狭长椭圆的无限大电致伸缩材料做了详细的分析，其中包含第一类边界条件（孔口自由）和第二类边界条件（刚性夹杂）的情况。在解的过程中作出了电学上的近似，即根据电场对应力、应变的贡献是关于电场强度的二次项这一事实，将力、电关系解耦。先不考虑电致伸缩效应的影响，单独通过静电场基本方程解出电场后，再利用解出的电场方程代入Knops的关于复变函数方法的基本方程得出了最后的结果。随后利用复变函数的保角变换，把含有椭圆缺陷的无限大体变换到含有单位圆的无限大体的变换域中，再对变换域中的孔口边界条件进行Cauchy积分，求出两个应力函数，最后即可求出变换域中的应力场和位移场。通过保角变换的逆变换即可解出物理平面上的应力和位移。Smith等人还对狭长椭圆蜕化成裂纹的情况进行了渐近展开，指出在裂纹尖端应力具有的奇异性。 2、电致伸缩材料的基本方程 从热力学理论可以得到弹性电介质本构方程如下[7] (1) 上式中为应力，电场电位移，应变和电场强度；和分别为等温等电场条件下四阶和六阶弹性系数；和为等温等应力二阶及三阶介电常数；， 和为等温条件下的压电系数，电弹性系数和电致伸缩系数。另外，式中下标小写拉丁字母取并且满足张量求和规则。对于没有压电效应的各向同性电致伸缩材料，其等温本构关系可以写为如下形式 (2) 其中，为等温条件下的两个相互独立的电致伸缩系数；为无应变情况下的介电常数；为Kronecker delta； 和为 常数，可以用杨氏弹性模量和泊松比表示为 另外，由电场产生的体积力为 (3) 这里， (4) 其中，为电场中电致伸缩应力，如果对非电致伸缩材料，上式可以简化为Maxwell应力。 将式()代入平衡方程，得到 (5) 我们称为伪总应力，利用式()以及式()，伪总应力可以取如下形式 (6) 图1 不同介质分界面示意图 对不均匀电介质，在不同介质分界面上，这个界面看作为内边界，即在所研究的对象内部（见图），力学量应该满足 (7) 上式中，和为介质I和介质II在分界面上的位移；，和分别为总面力，外加机械面力和电致伸缩面力。上标“＋”、“－”分别表示边界绕行方