傅里叶变换在材料中的应用

The Application of Fourier Transform in Materials

摘要: 从现代数学的眼光来看，傅里叶变换是一种特殊的积分变换。它能将满足一定条件的某个函数表示成正弦基函数的线性组合或者积分。在不同的研究领域，傅里叶变换具有多种不同的变体形式，如连续傅里叶变换和离散傅里叶变换。本文主要介绍了傅里叶变换在处理压电材料与热电材料的性能、接触、断裂等问题中的应用。

Abstract: From the perspective of modern mathematics, the Fourier transform is a special integral transform. It can express a certain function that satisfies certain conditions as a linear combination or integration of sine basis functions. In different research fields, Fourier transforms have a variety of different variations, such as a continuous Fourier transform and a discrete Fourier transform. This article mainly introduces the application of Fourier transform in dealing with the performance, contact and fracture of piezoelectric materials and thermoelectric materials.

1. 引言

Kelvin说，傅里叶理论不仅是现代分析中最美妙的结果之一，而且可以说，它为现代物理中每一个深奥问题的处理提供了一件必不可少的工具。随着近代物理的飞速发展，越来越多的实际问题需要用数学理论来解决。傅里叶变换是一类重要的积分变换，而积分变换能够将分析运算转换为代数运算，正是由于积分变换这一特性，在微分方程、偏微分方程、积分方程的求解中成为重要的方法之一。本文主要综述了傅里叶变换在压电材料、热电材料中的应用。

2. 傅里叶变换

2.1. 傅里叶变换的提出

傅里叶是一位法国数学家和物理学家，英文原名是Jean Baptiste Joseph Fourier (1768~1830)，Fourier对热传递很感兴趣，于1807年在法国科学学会上发表了一篇论文，运用正弦曲线来描述温度分布，论文里有个在当时具有争议性的决断：任何连续周期信号可以由一组适当的正弦曲线组合而成。但论文审查者之一的拉格朗日坚决反对这一观点，坚持认为傅里叶的方法无法表示带有棱角的信号，如在方波中出现非连续变化斜率。孰对孰错？拉格朗日是对的：正弦曲线无法组合成一个带有棱角的信号。但是，我们可以用正弦曲线来非常逼近地表示它，逼近到两种方法不存在能量差别，基于此，傅里叶是对的。

最初傅里叶分析是作为热过程的解析分析的工具被提出的。由于不少人对其进行不停地研究，所以在不同的研究领域，傅里叶变换具有多种不同的变体形式，如连续傅里叶变换和离散傅里叶变换。

2.2. 傅里叶变换的数学意义

傅里叶变换是一种分析信号的方法，它可以分析信号的成分，也可用这些成分合成信号。它表示能将满足一定条件的某个函数表示成三角函数(正弦或余弦函数)或者它们的积分的线性组合。

$f\left(t\right)$ 是t的周期函数，如果t满足狄利克雷条件：在一个以2T为周期内 $f\left(X\right)$ 连续或只有有限个第一类间断点，附 $f\left(X\right)$ 单调或可划分成有限个单调区间，则 $f\left(X\right)$ 以2T为周期的傅里叶级数收敛，和函数 $S\left(X\right)$ 也是以2T为周期的周期函数，且在这些间断点上，函数是有限值；在一个周期内具有有限个极值点；绝对可积 [1] 。则有(1)式成立

$F\left(\omega \right)=F\left[f\left(t\right)\right]={\displaystyle \underset{-\infty }{\overset{\infty }{\int }}f\left(t\right){\text{e}}^{-jwt}\text{d}t}$ (1)

上式称为积分运算 $f\left(t\right)$ 的傅里叶变换。

式(2)称为傅里叶逆变换

$f\left(t\right)=F^{-1}\left[F\left(\omega \right)\right]=\frac{1}{\text{2π}}{\displaystyle \underset{-\infty }{\overset{\infty }{\int }}F\left(\omega \right){\text{e}}^{jwt}\text{d}\omega }$ (2)

3. 傅里叶变换在材料中的应用

3.1. 傅里叶变换在压电材料中的应用

3.1.1. 由傅里叶变换得到奇异积分方程

H.S. Nan，B.L. Wang [2] 通过考虑残余表面应力对裂纹表面的影响，研究压电纳米材料的纳米级裂纹问题时，根据解决压电介质裂纹问题的方法，利用傅里叶积分变换和特征方程导出位移和势函数，通过引入满足边界条件的辅助函数，最终得到了奇异积分方程组，得到了裂纹尖端场。Jie Su，Liao-Liang Ke [3] 等首次尝试对刚性圆柱形冲头下压电半平面微动接触进行理论研究时，利用Fourier积分变换技术和叠加定理，将问题简化为一组耦合的Cauchy奇异积分方程。Zengtao Chen [4] 在研究尺寸对冲击载荷作用下压电带裂纹扭结的影响时，应用傅里叶变换和拉普拉斯变换将混合模式裂纹问题(模式I和II)减少为奇异积分方程组。Yue-Ting Zhou [5] 在研究三角形和圆柱形冲头对单斜压电材料摩擦滑动接触问题时，通过傅里叶变换将混合边值问题简化为第二类Cauchy型奇异积分方程。基于简化奇异积分方程的精确解，得到了各种表面应力和电位移的闭合表达式。钟献词 [6] 在积分变换方法求解压电带形中静止或运动的共线裂纹问题时，在电渗透裂纹和非电渗透型的裂纹假设下，利用Fourier积分变换把相应的混合边值问题转化为求解奇异积分方程。陈建 [7] 研究压电功能梯度材料层中平行于边界的动态反平面裂纹问题时，采用拉普拉斯、傅里叶积分变换和位错函数法将问题简化为Cauchy奇异积分方程，最后给出数值结果。贾洪波 [8] 在研究正交各向异性压电材料中III型裂纹对SH波的散射问题时，使用Fourier变换并引入位错密度函数，将混合边值问题转化为奇异积分方程。魏苓彦 [9] 在研究正交各向异性压电材料中III型裂纹受冲击力电载荷作用时的瞬态断裂问题，基于电可渗透裂纹模型，采用积分变换技术和数值求解奇异积分方程，最后反演得到了动态应力强度因子的表达式。B. Gu [10] 等在研究机械和电气作用下不同压电层间绝缘裂纹的瞬态响应问题时，利用拉普拉斯、傅里叶积分变换和Cauchy奇异积分方程方法，在电绝缘的边界条件下，解决了两个有限宽度的不同压电层受到机电冲击的界面裂纹面内问题。刘俊俏 [11] 用奇异积分方程方法研究SH波和压电压磁材料中周期裂纹的相互作用。根据电磁弹性材料的控制微分方程，通过傅立叶积分变换，将问题转化为Hilbert核的奇异积分方程组。李永东 [12] 在研究压电界面电极的应力强度因子时，建立了压电界面电极的力学模型，利用有限傅里叶变换方法将相应的力电耦合问题描述为奇异积分方程，并得到了界面应力及其强度因子的解析解。常志艳 [13] 在研究压电材料全平面含平行线裂纹问题时，通过引入势函数将原来获得的压电方程整理为较简单的偏微分方程组，再利用算子作用和Fourier变换等数学工具将偏微分方程转化为常微分方程求解。赵明皞 [14] 等运用广义不连续位移边界元法研究二维压电半导体材料在多场耦合环境下的断裂问题，利用傅里叶变换，给出二维压电半导体的广义不连续位移基本解。邢时超 [15] 研究了含半绝缘裂纹的压电条上、下表面在受到不同温度与电载荷作用时的断裂问题。基于压电材料的本构方程、热控制方程和边界条件，并结合Fourier变换和奇异积分公式法，获得了压电板条内温度场和应力场的解。杨昌锦，李尧臣 [16] 在进行有限厚压电层表面金属电极脱层的屈曲研究时，用Fourier积分变换，界面连续条件和上、下电极表面的边界条件将问题归为第二类Cauchy型奇异积分方程组。乐金朝 [17] 等使用Fourier积分变换，对非对称载荷作用下的Griffith裂纹问题进行了研究，利用裂纹两侧的位移和应力联结条件，将问题归结为一组带Cauchy核的奇异积分方程，通过Cauchy反演求得了问题的精确解。彭达仁 [18] 等在研究层状介质垂直界面有限裂纹问题时，运用傅里叶变换及引用

位错密度函数，导出了反映裂纹尖端奇异性的奇异积分方程组。S.Ueda [19] 研究了热电荷载作用下压电材料平行裂纹的热机电断裂行为。通过傅立叶变换，热机械问题和机电问题分别被简化为一个奇异积分方程组，并进行数值求解。S.Ueda [20] 在研究热电荷载作用下裂纹压电板的瞬态响应时，通过使用拉普拉斯变换和傅里叶变换，热机械问题和机电问题被简化为一个奇异积分方程组，并被数值求解。S.Ueda [21] 在研究压电带两个共面裂纹的热应力强度因子时，应用傅里叶变换将混合边值问题转换为奇异积分方程。S.Ueda [22] 在研究热载作用下功能梯度压电材料裂纹表面电导对强度因子的影响时，通过傅里叶变换，热机电问题分别被简化为奇异积分方程和奇异积分方程组，并进行了数值求解。S.Ueda [23] 在研究热电加载条件下压电带中两个平行裂纹之间的热机电相互作用时，通过傅里叶变换，热机械问题和机电问题分别被简化为奇异积分方程。S.Ueda & H.Hatano [24] 研究考虑了含有T形裂纹的无限压电材料的热电机械断裂问题。 通过使用傅里叶变换，热和机电问题分别被简化为一个奇异积分方程组。S.Ueda & H.Kondo [25] 在热冲击负载条件下研究了功能梯度压电材料板中的裂纹。通过使用拉普拉斯变换和傅立叶变换，热机械问题和机电问题被简化为一个奇异积分方程和一个奇异积分方程组，它们被数值求解。S.Ueda [26] 在研究功能梯度压电层的热诱导时，傅里叶变换技术用于奇异积分方程来解决问题。S.Ueda & Y.Ashida [27] 分析了功能梯度压电材料带中平行裂纹的无限行问题，在静态力学和瞬态热载荷条件下进行了分析。通过拉普拉斯变换和傅里叶变换，热机电问题被简化为一个奇异积分方程。S.Ueda & A.Lshii [28] 在研究热冲击载荷下压电材料带中两个不同长度的平行裂纹的热机电断裂行为，通过使用拉普拉斯变换和傅里叶变换，热机械问题和机电问题分别被简化为两个奇异积分方程组，并被数值求解。

以钟献词 [6] 的论文为例，我们来看一下傅里叶积分变换在压电材料中的具体运用。

考虑一高为h₂的带形压电材料,内部含有两条平行于边界表面的长为(a-c)的共线裂纹.即裂纹位于 $c\le \left|x\right|\le a$ 上，且裂纹距上下表面的距离分别为场 $h_I,h_{II}$ ( $h_I+h_{II}=2h$ )如图1。

这里xoy平面为各向同性面，与之垂直的方向(即z方向)为压电介质的极化方向。当带形表面受反平面机械荷载和平面内电场作用时，压电介质内仅有非零的反平面位移分量 $\omega \left(x,y\right)$ 与平面内电势 $\varphi \left(x,y\right)$ 在线电弹性理论框架下(忽略体力与自由电荷)，该问题的控制微分方程为

$c_{44}{\nabla }^2\omega +e_{15}{\nabla }^2\varphi =0,e_{15}{\nabla }^2\omega -{\epsilon }_{11}{\nabla }^2\varphi =0$ (3)

其中 $c_{44}$ 为压电材料的弹性常数， $e_{15}$ 为压电常数，为介电常数， ${\nabla }^2=\frac{{\partial }^2}{x^2}+\frac{{\partial }^2}{y^2}$ 为二维拉普拉斯算子。对一般的压电陶瓷，有 $c_{44}+\frac{e_{15}^2}{{\epsilon }_{11}}\ne 0$ ，所以，方程(3.1)可以简化为两个独立的调和方程：

${\nabla }^2\omega =0,{\nabla }^2\varphi =0$ (4)

一旦反平面位移分量 $\omega \left(x,y\right)$ 与平面内电势 $\varphi \left(x,y\right)$ 被确定，反平面剪应力和平面内电位移可由下列本构方程给出：

${\sigma }_{xz}=c_{44}{\gamma }_{xz}-e_{15}{E}_x,{\sigma }_{yz}=c_{44}{\gamma }_{yz}-e_{15}{E}_y$ (5)

$D_x=e_{15}{\gamma }_{xz}+{\epsilon }_{11}{E}_x,D_y=e_{15}{\gamma }_{yz}+{\epsilon }_{11}{E}_y$ (6)

${\gamma }_{xz}=\frac{\partial \omega }{\partial x},{\gamma }_{yz}=\frac{\partial \omega }{\partial y}$ (7)

$E_x=\frac{\partial \phi }{\partial x},D_y=e_{15}{\gamma }_{yz}+{\epsilon }_{11}{E}_y$ (8)

Figure 1. Piezoelectric material with two collinear cracks

图1. 内含两共线裂纹的带形压电材料

基于Dunn，Zhang和Tong的精确分析，电渗透型裂纹面上边界条件为：

$$ (9)

$D_y^I\left(x,0^{+}\right)=D_y^{II}\left(x,0^{-}\right)c<\left|x\right|<a,$ (10)

${\sigma }_{yz}^I\left(x,0^{+}\right)={\sigma }_{yz}^{II}\left(x,0^{-}\right)=0c<\left|x\right|<a.$ (11)

而在y = 0面上的其他部分电弹性场为连续的，即：

${\phi }^I\left(x,0^{+}\right)={\phi }^{II}\left(x,0^{-}\right)\left|x\right|<c,\left|x\right|>a$ (12)

$D_y^I\left(x,0^{+}\right)=D_y^{II}\left(x,0^{-}\right)\left|x\right|<c,\left|x\right|>a$ (13)

${\sigma }_{yz}^I\left(x,0^{+}\right)={\sigma }_{yz}^{II}\left(x,0^{-}\right)\left|x\right|<c,\left|x\right|>a$ (14)

${\omega }^I\left(x,0^{+}\right)={\omega }^{II}\left(x,0^{-}\right)\left|x\right|<c,\left|x\right|>a$ (15)

其中上标I和II分别表示 $0<y<h_I$ 和 $-h_{II}<y<0$ 区域内对应的量。

当压电介质在边界表面上受反平面剪应力和平面内电场作用时，即：

$\begin{array}{l}{\sigma }_{yz}^I\left(x,h_I\right)={\sigma }_{yz}^{II}\left(x,-h_{II}\right)={\sigma }_0,-\infty <x<+\infty ,\\ E_y^I\left(x,h_I\right)=E_y^{II}\left(x,-h_{II}\right)=E_0,-\infty <x<+\infty .\end{array}$ (16)

求解方法

从断裂力学的观点来看，感兴趣的是裂纹尖端的电弹性场。因此，在以后的解答中我们将着眼于解决裂纹尖端的奇异场，并给出相应的强度因子计算式和能量释放率。首先我们求出调和方程(4)的一般解。其次对所求问题予以讨论(由于对称性，只考虑压电介质的右半部分)。

由所考虑的间题关于y轴具有对称性，采用Fourier余弦变换及导数的Fourier变换定理，得：

${\displaystyle {\int }_0^{\infty }\left({\nabla }^2\omega \right)}\cos \left(\xi x\right)\text{d}x=\left(\frac{{\text{d}}^{\text{2}}}{\text{d}{y}^2}-{\xi }^2\right)\bar{\omega }\left(\xi ,y\right)$ (17)

其中是位移分量 $\omega \left(x,y\right)$ 的Fourier余弦变换，并且有反演公式：

$\omega \left(x,y\right)=\frac{2}{\text{π}}{\displaystyle {\int }_0^{\infty }\bar{\omega }\left(\xi ,y\right)\cos \left(\xi x\right)\text{d}\xi }$ (18)

由方程 ${\nabla }^2\omega =0$ 可得：

$\bar{\omega }\left(\xi ,y\right)=A\left(\xi \right)\exp \left(-\xi y\right)+B\left(\xi \right)\exp \left(\xi y\right)$ (19)

同时，为了使得 $\omega \left(x,y\right)$ 满足边界条件，取：

$\omega \left(x,y\right)=\frac{2}{\text{π}}{\displaystyle {\int }_0^{\infty }\left[A\left(\xi \right)\exp \left(-\xi y\right)+B\left(\xi \right)\exp \left(\xi y\right)\right]\cos \left(\xi x\right)\text{d}\xi }+A_{0}y$ (20)

同样的方法，可得：

$\phi \left(x,y\right)=\frac{2}{\text{π}}{\displaystyle {\int }_0^{\infty }\left[C\left(\xi \right)\exp \left(-\xi y\right)+D\left(\xi \right)\exp \left(\xi y\right)\right]\cos \left(\xi x\right)\text{d}\xi }+C_{0}y$ (21)

并且有应力，电位移可按下列关系表示

$\begin{array}{c}{\sigma }_{yz}\left(x,y\right)=\frac{2c_{44}}{\text{π}}{\displaystyle {\int }_0^{\infty }\left(-\xi \right)}\left[A\left(\xi \right)\exp \left(-\xi y\right)-B\left(\xi \right)\exp \left(\xi y\right)\right]\cos \left(\xi x\right)\text{d}\xi +A_0{c}_{44}+C_0{e}_{15}\\ +\frac{2e_{15}}{\text{π}}{\displaystyle {\int }_0^{\infty }\left(-\xi \right)}\left[C\left(\xi \right)\exp \left(-\xi y\right)-D\left(\xi \right)\exp \left(\xi y\right)\right]\cos \left(\xi x\right)\text{d}\xi ,\end{array}$ (22)

$\begin{array}{c}{D}_y\left(x,y\right)=\frac{2e_{15}}{\text{π}}{\displaystyle {\int }_0^{\infty }\left(-\xi \right)}\left[A\left(\xi \right)\exp \left(-\xi y\right)-B\left(\xi \right)\exp \left(\xi y\right)\right]\cos \left(\xi x\right)\text{d}\xi +A_0{e}_{15}-C_0{\epsilon }_{11}\\ -\frac{2{\epsilon }_{11}}{\text{π}}{\displaystyle {\int }_0^{\infty }\left(-\xi \right)}\left[C\left(\xi \right)\exp \left(-\xi y\right)-D\left(\xi \right)\exp \left(\xi y\right)\right]\cos \left(\xi x\right)\text{d}\xi ,\end{array}$ (23)

由边界条件及一系列代数运算，该问题可被化简为(具体过程可见论文原文)

$\frac{4}{\text{π}}{\displaystyle {\int }_c^a\Delta {{\omega }^{\prime }}_1}\left(t\right)\text{d}t{\displaystyle {\int }_0^{\infty }\frac{\tanh \left(h_I\xi \right)\tanh \left(h_{II}\xi \right)}{\tanh \left(h_I\xi \right)+\tanh \left(h_{II}\xi \right)}}\cos \left(x\xi \right)\sin \left(t\xi \right)\text{d}\xi =-\frac{{\sigma }_0}{c_{44}}c<x<a$ (24)

进一步，由已知公式

$\frac{2}{\text{π}}{\displaystyle {\int }_0^{\infty }\cos }\left(x\xi \right)\sin \left(t\xi \right)\text{d}\xi =\frac{1}{\text{π}}\left(\frac{\text{1}}{\text{t}+x}+\frac{1}{t-x}\right)$

方程可被化为求解下列奇异积分方程：

$\frac{2}{\text{π}}{\displaystyle {\int }_c^a\frac{\Delta {{\omega }^{\prime }}_1\left(t\right)}{t-x}\text{d}t}+\frac{2}{\text{π}}{\displaystyle {\int }_c^a\Delta {{\omega }^{\prime }}_1\left(t\right)k^{\left(1\right)}\left(t,x,h_I,h_{II}\right)\text{d}t}=-\frac{{\tau }^{\left(1\right)}}{c_{44}}$ (25)

其中 ${\tau }^{\left(1\right)}={\sigma }_0$ ，且

$k^{\left(1\right)}\left(t,x,h_I,h_{II}\right)={\displaystyle {\int }_0^{\infty }2\left[\frac{\tanh \left(h_I\xi \right)\tanh \left(h_{II}\xi \right)}{\tanh \left(h_I\xi \right)+\tanh \left(h_{II}\xi \right)}-1\right]\cos \left(x\xi \right)\sin \left(t\xi \right)\text{d}\xi }+\frac{1}{t+x}$

以上就是利用Fourier积分变换把相应的混合边值问题转化为求解奇异积分方程的过程。

3.1.2. 由傅里叶变换得到对偶积分方程组

Keqiang Hu [29] 等研究Dugdale型裂纹沿着两种不同的压电材料的界面移动的问题时，应用傅里叶变换将裂纹的混合边值问题转化为对偶积分方程，得到了精确求解。Yi-Ze Wang [30] 在讨论残余极化和周围介质中麦克斯韦应力对两种不同压电材料之间的反平面界面裂纹的影响时，通过傅里叶变换和它的逆变换，它被转化成两对双积分方程。推导了应力强度因子和电位移的闭合解。Jeong Woo Shin [31] 采用积分变换方法，分析了两种不同功能梯度压电材料界面裂纹在反平面剪切冲击载荷作用下的瞬态响应，用拉普拉斯变换和傅里叶变换将问题简化为两组对偶积分方程，然后将其表示为第二类的Fredholm积分方程。R. Bagheri [32] 等在研究功能梯度压电条中多个运动裂纹的解析解时，利用复数傅里叶变换获得了变换的位移和应力场。马旭 [33] 通过对裂纹的无限长功能梯度/压电材料层合的SH波散射问题的研究，用Fourier积分变换把混合边值问题的求解转化为对偶积分方程，然后利用Copson方法将得到的对偶积分方程转化为Fredholm积分方程再进行数值求解。郭玉彬 [34] 在忽略界面上裂纹尖端裂纹面相互叠入的前提下，讨论了功能梯度材料与压电材料拼接界面上的反平面运动裂纹问题。通过Fourier积分变换，将混合边值问题转化为对偶积分方程，并利用Copson-Sih方法将对偶积分方程转化为第二类Fredholm积分方程进行求解，给出了反平面位移、电势及应力分量的解析表达式。杜勇峰、王建国 [35] 采用Schmidt方法研究条形压电材料和弹性材料的界面III型裂纹问题，主要探讨条形压电材料和弹性材料受反平面剪应力和平面电位移作用的情形；通过Fourier积分变换和界面裂纹位移差的车贝雪夫多项式假设，并利用Schmidt方法得到数值解。边文凤 [36] 在研究压电材料I型裂纹动态问题的对偶方程组及其求解的情况时，采用Laplace变换、半无限对称Fourier 正弦变换和Fourier余弦变换，对用势函数表示的压电材料基本微分方程进行变换和初步求解；然后通过Fourier反演和引入边界条件，建立了二维压电材料动态裂纹问题的对偶方程组，曾云 [37] 用含有电场梯度效应的电弹性体理论分析了压电体中的反平面裂纹问题。利用Fourier积分变换方法，将相应的复合边值问题转化为对偶积分方程组。求解这些方程组，获得了裂纹尖端的强度因子和能量释放率。刘海涛 [38] 在研究压电材料中矩形裂纹材料的三维断裂问题时，通过二维Fourier变换技术将此类问题的求解转化为三对以矩形裂纹上下面的位移间断函数为未知变量的对偶积分方程，并采用Schmidt方法进行数值求解。胡克强 [39] 在研究矩形压电体中反平面裂纹的电弹性场的问题时，利用Fourier积分变换和Fourier正弦级数将电绝缘型裂纹问题化为对偶积分方程，并进一步归结为易于求解的第二类Fredholm积分方程。哈元军 [40] 在讨论无限长功能梯度压电板中双裂纹尖端热响应分析时，用Fourier积分变换把混合边值问题转化为以裂纹面上位移为未知量的对偶积分方程，利用Schmidt方法对其求解。杨娟 [41] 利用积分变换及奇异积分方程技术研究电磁复合材料底层处裂纹对SH波散射问题。假定裂纹面的边界条件为电渗透性，通过Fourier余弦变化将问题转化为对偶积分方程，并利用Copson方法将对偶积分方程转化为第二类Fredholm积分方程。冯文杰 [42] 等在研究加层压电条界面裂纹的稳态扩展时，应用Fourier积分变换将问题化为以第二类Fredholm积分方程表示的对偶积分方程，导出了相应的动应力强度因子表达式。陆万顺 [43] 等分析了功能梯度压电材料中裂纹尖端的热应力，针对考虑的问题，通过Fourier积分变换把混合边值问题的求解转化为裂纹面上位移间断为未知量的对偶积分方程，然后利用Schmidt方法来求解。胡克强 [44] 等基于三维弹性理论和压电理论，研究了功能梯度压电板条中的电渗透型运动裂纹问题。利用Fourier积分变换方法，将混合边值问题化为对偶积分方程，并进一步归结为易于数值求解的第二类Fredholm积分方程。曾云 [45] 等用含有电场梯度效应的电弹性体理论分析了压电体中的III型裂纹问题，利用傅里叶积分变换方法，将相应的混合边值问题转化为对偶积分方程组，求解这些方程组，得到了裂纹尖端处的电场强度因子。

3.2. 傅里叶积分变换在热电材料中的应用

热电材料是一种将热能和电能相互转换的功能材料。其应用主要包括：温差发电、半导体制冷，以及作为传感器和温度控制器在微电子器件和MEMS中的应用等 [46] 傅里叶积分变换在处理热电问题时，也起到了很大的作用。

Magdy A. Ezzat [47] 在修改傅里叶定律的热电MHD时，构建了一个具有分数阶导数和热电性质的热方程的数学模型。在边界层框架内的控制耦合方程应用于斯托克斯在存在垂直于板的横向磁场的情况下热源分布的无限板上的导电流的第一个问题，拉普拉斯变换和状态空间方法技术被用来获得解决方案。拉普拉斯变换是利用变换的复数反演公式和傅里叶展开技术得到的。Magdy A. Ezzat和Hamd M. Youssef [48] 在研究斯托克斯对于具有热电性质的导电微极流体第一个问题时，将边界层模型框架内的控制耦合方程应用于斯托克斯热源的第一个问题。拉普拉斯变换和傅立叶变换技术被用来获得解决方案。傅里叶变换的倒数是分析获得的。拉普拉斯变换是利用变换的复数反演公式和傅里叶展开技术得到的。该方法也在Magdy A. Ezzat [49] 在研究具有分数导数热传递的热电MHD非牛顿流体中得到应用。B.L.Wang [50] [51] [52] 等在研究裂纹对热电转换效率的影响时，利用傅里叶变换技术和奇异积分方程方法，得到了电势和温度场分布的解析表达式，得到了裂纹尖端的电流密度和强度因子，有助于确定热电转换效率。以B.L.Wang文章为例，在研究裂纹对热电转换效率的影响时，为了使计算过程简洁容易，引进一总的势函数F，F表示总的电势，通过对原有的本构方程进行代入和化简，可得到T和温度场T的表达式

${\nabla }^{2}F=0$ (26)

$k{\nabla }^{2}T+\sigma {\left(\nabla F\right)}^2=0$ (27)

应用傅里叶变换可得(3)式的解为

$F\left(x,y\right)\frac{F_{10}-F_{20}}{2h}y+\frac{F_{10}+F_{20}}{2}=+{\displaystyle \underset{n=1}{\overset{\infty }{\sum }}{\displaystyle \underset{m=1}{\overset{2}{\sum }}\sin \left(s_{n}y\right)}}\exp \left(s_{n}x{\lambda }_m\right)A_{mn}$ $+\frac{2}{\text{π}}{\displaystyle {\int }_0^{\infty }{\displaystyle \underset{m=1}{\overset{2}{\sum }}\cos \left(sx\right)}}\exp \left(sy{\gamma }_m\right)B_{m}ds$ (28)

温度场分布的解析表达式也可由傅里叶变换得到，进而得到了裂纹尖端的电流密度和强度因子，确定热电转换效率。

$j_{ey}\left(x,y\right)$ $=j_{e0}-\sigma {\displaystyle \underset{n=1}{\overset{\infty }{\sum }}{\displaystyle \underset{m=1}{\overset{2}{\sum }}{s}_n}}\cos \left(s_{n}y\right)\exp \left(s_{n}x{\lambda }_m\right)A_{mn}-\frac{2\sigma }{\text{π}}{\displaystyle \underset{m=1}{\overset{2}{\sum }}{\gamma }_m}{\displaystyle {\int }_0^{\infty }s\cos \left(sx\right)}\exp \left(sy{\gamma }_m\right)B_m\text{d}s$ (29)

$j_{ex}\left(x,y\right)=$ $-\sigma {\displaystyle \underset{n=1}{\overset{\infty }{\sum }}{\displaystyle \underset{m=1}{\overset{2}{\sum }}{s}_n}}{\lambda }_m\sin \left(s_{n}y\right)\exp \left(s_{n}x{\lambda }_m\right)A_{mn}+\frac{2\sigma }{\text{π}}{\displaystyle \underset{m=1}{\overset{2}{\sum }}{\displaystyle {\int }_0^{\infty }s\sin \left(sx\right)\exp \left(sy{\gamma }_m\right)B_m\text{d}s}}$ (30)

这儿 $j_{e0}=\frac{-\sigma \left(F_{10}-F_{20}\right)}{2h}$ ，未知系数 $A_{mn}$ 和 $B_m$ 可被新的函数 $f_0\left(x\right)$ 表示， $f_0\left(x\right)$ 被定义为

$f_0\left(x\right)=\{\begin{array}{l}\frac{\partial }{\partial x}\left[F\left(x,+0\right)-F\left(x,-0\right)\right],\left|x\right|\le a\\ 0,\left|x\right|>a\end{array}$ (31)

经过一系列的运算，未知系数 $B_m$ 可被表示为

$B_m=N_m\left(s\right){\displaystyle {\int }_0^a{f}_0\left(r\right)\sin \left(sr\right)\text{d}r}$ (32)

此时

$N_1\left(s\right)=-\frac{1}{s\left[1-\exp \left(2sh{\gamma }_1\right)\right]},N_2\left(s\right)=-\exp \left(2sh{\gamma }_1\right)N_1(s)$

根据对称性，可得到 $A_{mn}$ 的表达式

$A_{1n}=A_{2n}=-\frac{1}{2h}{\displaystyle {\int }_0^a{f}_0\left(r\right)v\left(s_n,r\right)\text{d}r}$ (33)

此时 $v\left(s_n,r\right)=\frac{\exp \left[-\left(c-r\right)s_n\right]-\exp \left[-\left(c+r\right)s_n\right]}{s_n\sinh \left(s_{n}c\right)}$ 。

经过运算以及边界条件，电流密度可用奇异积分方程表示出来

$\frac{1}{\text{π}}{\displaystyle {\int }_{-a}^a{f}_0\left(r\right)\left[\frac{1}{r-x}+M_0\left(r,x\right)\right]\text{d}r}=\frac{j_{e0}}{\sigma }$ (34)

其中 $f_0\left(x\right)$ 还满足条件

${\displaystyle {\int }_{-a}^a{f}_0\left(r\right)\text{d}r}=0$ (35)

一旦 $f_0\left(x\right)$ 被确定了， $F\left(x,y\right)$ 也被确定了。

4. 结束语

傅里叶变换在物理学、电子类学科、数论、组合数学、信号处理、概率论、统计学、声学、光学、海洋学等领域都有着广泛的应用。傅里叶变换在求解压电材料的断裂、接触等问题的过程中起到了重大作用。热电材料作为一种新型的材料，具有重要的学术价值和应用前景，但傅里叶变换在热电材料的很多领域都未涉及到，作者认为在研究热电材料的断裂、接触等问题上，应在借鉴处理压电材料问题的基础上，有所创新。

基金项目

国家自然科学基金项目(11762016)；宁夏自然科学研究项目(NZ17009)；宁夏高等学校优秀轻年培育项目(NGY2017002)。