An Axisymmetric Contact Problem of Piezoelectric Materials Based on the Couple Stress Theory
-
摘要: 基于偶应力理论,研究了刚性绝缘球压头与横观各向同性压电半空间尺度依赖的轴对称接触问题. 利用Hankel积分变换和积分最小二乘法,获得了接触压力. 讨论了特征材料长度对接触压力分布、接触半径和压痕深度的影响. 结果表明,基于偶应力理论获得的接触压力结果明显大于经典结果.Abstract: Based on the couple stress theory, the axisymmetric contact problem between a rigid insulating spherical punch and a transversely isotropic piezoelectric half-space was studied. With the Hankel integral transform and the integral least squares approach, the analytical solutions of the contact pressure were obtained. The effects of the characteristic material length on the contact pressure distribution, the contact radius and the indentation depth were discussed. The results indicate that, the contact pressure obtained based on the couple-stress theory is significantly greater than the classical results.
-
Key words:
- axisymmetric contact /
- piezoelectric material /
- size effect /
- couple stress theory
-
0. 引言
经典连续介质力学是一种基于物质连续分布假设的近似分析方法,这一理论为分析宏观尺度上材料的力学行为提供了合理的依据. 在经典连续介质力学中,小变形假设下旋转运动不造成形变,所以位移梯度全部用对称的应力张量表示是合理的. 然而实验表明,当材料的尺度减小到微米量级时,材料内禀性能及应用规律和原理表现出与宏观尺度明显的不同[1-3]. 为了建立更完整的连续介质理论,需要引入新的形变度量来描述介质内部存在的微观结构(裂纹、孔洞、位错或晶格缺陷)所导致的力学不对称,例如曲率张量. 在这样的理论框架下就引入了偶应力这一概念. 偶应力弹性理论是最基本的梯度型广义连续介质力学理论,其他理论方法分支多是在偶应力理论的基础上建立起来的.
尺度效应是微/纳米材料的一个普遍特征,它导致了各种与尺寸有关的行为、特性和现象. 随着接触尺寸的减小,材料内部的特征长度对宏观响应的影响变得越来越明显. 基于偶应力理论,Zisis课题组[4-8]系统地研究了材料的接触行为,他们研究的结果表明,经典接触理论预测的结果与在偶应力理论框架下得到的结果有很大的差异. Karuriya和Bhandakkar[9]研究了有限厚度弹性层与不同形状压头的尺度依赖的压痕问题. 基于偶应力理论,Song等[10-11]考虑了均匀弹性半平面的滑动摩擦和均匀涂层半平面的无摩擦接触问题. Wang等[12]分析了尺度依赖的部分滑移接触问题. Çömez和El-Borgi[13]探究了尺度依赖的层状结构的滑动摩擦接触问题. 总而言之,偶应力理论在接触问题中得到了广泛应用[14-18].
随着智能器件的日益小型化,压电材料被广泛应用于各种微型设备中,如微型发电机、制动器和传感器等[19-21]. 这些压电微型设备在实际使用中通常会遭受集中的力电载荷,因而会不可避免地产生接触、摩擦和磨损. 在偶应力理论框架下,关于压电材料中波的传播[22]、弯曲[23]、自由振动[24]和断裂[25]等问题得到了较多关注. 但是关于压电材料尺度依赖接触问题的研究非常有限,目前仅见Sun等[26]对压电材料的二维尺度依赖接触问题进行了讨论. 轴对称接触问题在工业生产活动中广泛存在,同时轴对称接触可以消除只考虑二维近似时带来的边界效应,然而目前还没有发现关于压电材料尺度依赖轴对称接触问题的相关文献报道.
接触问题是一个混合边值问题,目前涌现出几种具有代表性的处理方法:第一种方法是利用积分变换,将混合边值问题转化为奇异积分方程,再将积分方程离散化为代数方程后进行数值求解[27-28]. 第二种方法是用近似的边界条件代替精确的边界条件,然后使用基函数迭代求解弹性接触问题[29-30]. 这两种方法都可以获得满意的结果[27-30]. 第二种方法能够给出接触压力的显式表达式,通常需要数量较少的基函数就能获得较为精确的结果,因此本文选择第二种方法来研究接触问题.
本文利用偶应力理论,对压电半空间与球压头尺度依赖的无摩擦接触问题进行了理论分析. 利用材料的特征材料长度来描述尺度效应. 利用Hankel积分变换和积分最小二乘法,得到了接触压力的解析解. 最后,讨论了特征材料长度对接触压力分布、接触半径和压痕深度的影响.
本文研究工作的创新点如下:1) 基于偶应力理论,建立了压电材料尺度依赖轴对称无摩擦接触模型,系统地分析了压电材料尺度依赖轴对称无摩擦接触行为. 2) 研究发现,利用偶应力理论可以有效地描述微尺度压电材料的尺度依赖接触力学行为,研究结果为改善微型智能压电器件的力电接触损伤以及压电材料的优化设计和工程应用提供了理论依据.
1. 偶应力理论下的压电材料接触问题
图 1描述了横观各向同性压电半空间与刚性绝缘球压头的尺度依赖的轴对称接触模型. 球压头的半径为R,球压头上表面承受集中载荷F,圆形接触区域的半径为a. 压电材料的极化方向沿着z轴方向.
1.1 基本方程
在偶应力弹性理论中,考虑轴对称问题,几何方程表示为[8]
$$\varepsilon_{r r}=\frac{\partial u_{r}}{\partial r}, \varepsilon_{\theta \theta}=\frac{u_{r}}{r}, \varepsilon_{z z}=\frac{\partial u_{z}}{\partial z}, \varepsilon_{r z}=\frac{1}{2}\left(\frac{\partial u_{z}}{\partial r}+\frac{\partial u_{r}}{\partial z}\right), $$ (1) $$D_{r}=-\frac{\partial \phi}{\partial r}, D_{z}=-\frac{\partial \phi}{\partial z} , $$ (2) $$\omega_{\theta}=\frac{1}{2}\left(\frac{\partial u_{r}}{\partial z}-\frac{\partial u_{z}}{\partial r}\right), $$ (3) $$k_{\theta r}=-\frac{\omega_{\theta}}{r}, k_{r \theta}=\frac{\partial \omega_{\theta}}{\partial r}, k_{z \theta}=\frac{\partial \omega_{\theta}}{\partial z}, $$ (4) 其中,εrr,εθθ,εzz和εrz为应变分量;ur和uz为位移分量;Dr和Dz为电位移分量;ϕ为电势;ωθ为旋转分量;kθr,krθ和kzθ为曲率张量的分量.
$$\frac{\partial \sigma_{r r}}{\partial r}+\frac{\partial \sigma_{z r}}{\partial z}+\frac{\sigma_{r r}-\sigma_{\theta \theta}}{r}=0 , $$ (5) $$\frac{\partial \sigma_{r z}}{\partial r}+\frac{\partial \sigma_{z z}}{\partial z}+\frac{\sigma_{r z}}{r}=0 , $$ (6) $$ \frac{\partial D_{r}}{\partial r}+\frac{\partial D_{z}}{\partial z}+\frac{D_{r}}{r}=0 , $$ (7) $$\frac{\partial m_{r \theta}}{\partial r}+\frac{\partial m_{z \theta}}{\partial z}+\frac{m_{r \theta}+m_{\theta r}}{r}+\sigma_{z r}-\sigma_{r z}=0 , $$ (8) 其中,mrθ,mθr,mzθ和mθz为偶应力分量. 在偶应力理论中,应力分量σzr和σrz是非对称的,可以写作[17]
$$\sigma_{z r}=\tau_{z r}+\alpha_{z r}, \sigma_{r z}=\tau_{r z}+\alpha_{r z}, $$ (9) 其中
$$\tau_{z r}=\tau_{r z}, \alpha_{z r}=-\alpha_{r z}, $$ (10) 这里,τrz和τzr是应力的对称部分;αrz和αzr是应力的反对称部分.
将方程(9)、(10)代入方程(8),可以得到
$$\tau_{z r}=-\tau_{r z}=-\frac{1}{2}\left(\frac{\partial m_{r \theta}}{\partial r}+\frac{m_{r \theta}+m_{\theta r}}{r}+\frac{\partial m_{z \theta}}{\partial z}\right) .$$ (11) Wang等[32]提出了具有偶应力的横观各向同性压电材料的本构方程. 对于横观各向同性材料,特征材料长度和剪切模量在不同方向上是不同的[33-34]. 在横观各向同性压电材料的尺度依赖接触问题的研究中应考虑这一效应. 然而,由于难以通过实验方法测量不同方向上的特征材料长度和剪切模量,目前还没有发现相关的报道. 为了得到解析解,我们假设不同方向上的特征材料长度和剪切模量相同. 压电半空间的本构方程表示为[8, 31]
$$\sigma_{z z}=C_{13}\left(\frac{\partial u_{r}}{\partial r}+\frac{u_{r}}{r}\right)+C_{33} \frac{\partial u_{z}}{\partial z}+e_{33} \frac{\partial \phi}{\partial z}, $$ (12) $$ \sigma_{r z}+\sigma_{z r}=2 C_{44}\left(\frac{\partial u_{r}}{\partial z}+\frac{\partial u_{z}}{\partial r}\right)+2 e_{15} \frac{\partial \phi}{\partial r} , $$ (13) $$D_{z}=e_{31}\left(\frac{\partial u_{r}}{\partial r}+\frac{u_{r}}{r}\right)+e_{33} \frac{\partial u_{z}}{\partial z}-\varepsilon_{33} \frac{\partial \phi}{\partial z} , $$ (14) $$m_{r \theta}=4 G l^{2}\left(\frac{\partial \omega_{\theta}}{\partial r}-\frac{\beta}{r} \omega_{\theta}\right), m_{\theta r}=4 G l^{2}\left(\beta \frac{\partial \omega_{\theta}}{\partial r}-\frac{1}{r} \omega_{\theta}\right), $$ (15) $$m_{z \theta}=4 G l^{2} \frac{\partial \omega_{\theta}}{\partial z}, m_{\theta z}=\beta m_{z \theta}, $$ (16) 其中,Cjk为弹性常数;ejk为压电常数;εjk为介电常数;l为特征材料长度,与材料的微观结构有关,包括微裂缝、空隙、位错或晶格缺陷等[10];这里,G=C66=(C11-C12)/2[26];β是一个无量纲参数[8],且-1<β<1.
把本构方程(12)—(16)代入到平衡方程(5)—(8),并利用几何方程(1)—(4),得到如下的位移控制方程:
$$ \begin{align*} & C_{11}( \left(\frac{\partial^{2} u_{r}}{\partial r^{2}}+\frac{1}{r} \frac{\partial u_{r}}{\partial r}-\frac{u_{r}}{r^{2}}\right)+C_{44} \frac{\partial^{2} u_{r}}{\partial z^{2}}+\left(C_{13}+C_{44}\right) \frac{\partial^{2} u_{z}}{\partial r \partial z}+ \\ &\ \ \ \ \ \ \ \ \left(e_{31}+e_{15}\right) \frac{\partial^{2} \phi}{\partial r \partial z}+G l^{2}\left[\frac{1}{r^{2}}\left(\frac{\partial^{2} u_{r}}{\partial z^{2}}-\frac{\partial^{2} u_{z}}{\partial r \partial z}\right)+\right. \\ &\ \ \ \ \ \ \ \ \left.\frac{1}{r}\left(\frac{\partial^{3} u_{z}}{\partial r^{2} \partial z}-\frac{\partial^{3} u_{r}}{\partial r \partial z^{2}}\right)-\frac{\partial^{4} u_{r}}{\partial r^{2} \partial z^{2}}-\frac{\partial^{4} u_{r}}{\partial z^{4}}+\frac{\partial^{4} u_{z}}{\partial r \partial z^{3}}+\frac{\partial^{4} u_{z}}{\partial r^{3} \partial z}\right]=0, \end{align*}$$ (17) $$ \begin{align*}& C_{44}\left(\frac{\partial^{2} u_{z}}{\partial r^{2}}+\frac{1}{r} \frac{\partial u_{z}}{\partial r}\right)+C_{33} \frac{\partial^{2} u_{z}}{\partial z^{2}}+\left(C_{13}+C_{44}\right)\left(\frac{\partial^{2} u_{r}}{\partial r \partial z}+\frac{1}{r} \frac{\partial u_{r}}{\partial z}\right)+ \\ &\ \ \ \ \ \ \ \ e_{15}\left(\frac{\partial^{2} \phi}{\partial r^{2}}+\frac{1}{r} \frac{\partial \phi}{\partial r}\right)+e_{33} \frac{\partial^{2} \phi}{\partial z^{2}}+G l^{2}\left[\frac{1}{r^{3}}\left(\frac{\partial u_{r}}{\partial z}-\frac{\partial u_{z}}{\partial r}\right)+\right. \\ &\ \ \ \ \ \ \ \ \frac{1}{r^{2}}\left(\frac{\partial^{2} u_{z}}{\partial r^{2}}-\frac{\partial^{2} u_{r}}{\partial r \partial z}\right)+\frac{1}{r}\left(\frac{\partial^{3} u_{r}}{\partial z^{3}}-\frac{\partial^{3} u_{z}}{\partial r \partial z^{2}}+2 \frac{\partial^{3} u_{r}}{\partial r^{2} \partial z}-2 \frac{\partial^{3} u_{z}}{\partial r^{3}}\right)+ \\ &\ \ \ \ \ \ \ \ \left.\frac{\partial^{4} u_{r}}{\partial r \partial z^{3}}+\frac{\partial^{4} u_{r}}{\partial r^{3} \partial z}-\frac{\partial^{4} u_{z}}{\partial r^{2} \partial z^{2}}-\frac{\partial^{4} u_{z}}{\partial r^{4}}\right]=0, \phi{\partial z^{2}}=0 . \end{align*} $$ (18) $$ \begin{align*}& e_{15}\left(\frac{\partial^{2} u_{z}}{\partial r^{2}}+\frac{1}{r} \frac{\partial u_{z}}{\partial r}\right)+e_{33} \frac{\partial^{2} u_{z}}{\partial z^{2}}+\left(e_{15}+e_{31}\right)\left(\frac{\partial^{2} u_{r}}{\partial r \partial z}+\frac{1}{r} \frac{\partial u_{r}}{\partial z}\right)- \\ & \ \ \ \ \ \ \ \ \varepsilon_{11}\left(\frac{\partial^{2} \phi}{\partial r^{2}}+\frac{1}{r} \frac{\partial \phi}{\partial r}\right)-\varepsilon_{33} \frac{\partial^{2} \phi}{\partial z^{2}}=0 . \end{align*}$$ (19) 注意,当l=0时,方程(17)—(19)可以退化为不考虑尺度效应的经典轴对称接触问题的控制方程.
定义Hankel积分变换为
$$\hat{H}(s)=\int_{0}^{\infty} H(r) r \mathrm{J}_{m}(s r) \mathrm{d} r, $$ (20) 则其逆变换为
$$H(r)=\int_{0}^{\infty} \hat{H}(s) s \mathrm{J}_{m}(s r) \mathrm{d} s ,$$ (21) 其中,Jm(·)是m阶第一类Bessel函数;s是积分变换的变量.
对方程(17)—(19)应用Hankel积分变换,方程(17)—(19)变为
$$ \begin{align*} & C_{44} \frac{\mathrm{d}^{2} \hat{u}_{r}}{\mathrm{d} z^{2}}-s^{2} C_{11} \hat{u}_{r}-s\left(C_{13}+C_{44}\right) \frac{\mathrm{d} \hat{u}_{z}}{\mathrm{d} z}-s\left(e_{31}+e_{15}\right) \frac{\mathrm{d} \hat{\phi}}{\mathrm{d} z}+ \\ & \ \ \ \ \ \ \ G l^{2}\left(s^{2} \frac{\mathrm{d}^{2} \hat{u}_{r}}{\mathrm{d} z^{2}}+s^{3} \frac{\mathrm{d} \hat{u}_{z}}{\mathrm{d} z}-\frac{\mathrm{d}^{4} \hat{u}_{r}}{\mathrm{d} z^{4}}-s \frac{\mathrm{d}^{3} \hat{u}_{z}}{\mathrm{d} z^{3}}\right)=0, \end{align*}$$ (22) $$ \begin{align*}& s\left(C_{13}+C_{44}\right) \frac{\mathrm{d} \hat{u}_{r}}{\mathrm{d} z}+C_{33} \frac{\mathrm{d}^{2} \hat{u}_{z}}{\mathrm{d} z^{2}}-s^{2} C_{44} \hat{u}_{z}+e_{33} \frac{\mathrm{d}^{2} \hat{\phi}}{\mathrm{d} z^{2}}-s^{2} e_{15} \hat{\phi}+ \\ & \ \ \ \ \ \ \ G l^{2}\left(s \frac{\mathrm{d}^{3} \hat{u}_{r}}{\mathrm{d} z^{3}}-s^{4} \hat{u}_{z}+s^{2} \frac{\mathrm{d}^{2} \hat{u}_{z}}{\mathrm{d} z^{2}}-s^{3} \frac{\mathrm{d} \hat{u}_{r}}{\mathrm{d} z}\right)=0, \end{align*} $$ (23) $$ \begin{align*}& s\left(e_{15}+e_{31}\right) \frac{\mathrm{d} \hat{u}_{r}}{\mathrm{d} z}+e_{33} \frac{\mathrm{d}^{2} \hat{u}_{z}}{\mathrm{d} z^{2}}-s^{2} e_{15} \hat{u}_{z}-\varepsilon_{33} \frac{\mathrm{d}^{2} \hat{\phi}}{\mathrm{d} z^{2}}+s^{2} \varepsilon_{11} \hat{\phi}=0 . \end{align*}$$ (24) 方程(22)—(24)的一般解为
$$\left[\hat{u}_{r}, \hat{u}_{z}, \hat{\phi}\right]=\sum\limits_{i=1}^{4}\left[1, a_{i}(s), b_{i}(s)\right] A_{i}(s) \mathrm{e}^{n_{i} z}, $$ (25) 这里Ai(s)(i=1, 2, 3, 4)是未知的常数,ni是下列特征方程具有正实部的特征根
$$\operatorname{det}\left[\bar{h}_{i j}\left(s, n_{i}\right)\right]=0, $$ (26) 其中
$$ \left[\bar{h}_{i j}\right]=\left[\begin{array}{ccc} C_{44} n^{2}-C_{11} s^{2}+G l^{2} s^{2} n^{2}-G l^{2} n^{4} & -s n\left(C_{13}+C_{44}\right)+G l^{2} s^{3} n-G l^{2} s n^{3} & -s\left(e_{31}+e_{15}\right) n \\ s n\left(C_{13}+C_{44}\right)+G l^{2} s n^{3}-G l^{2} s^{3} n & C_{33} n^{2}-C_{44} s^{2}-G l^{2} s^{4}+G l^{2} s^{2} n^{2} & e_{33} n^{2}-e_{15} s^{2} \\ s n\left(e_{31}+e_{15}\right) & e_{33} n^{2}-e_{15} s^{2} & \varepsilon_{11} s^{2}-\varepsilon_{33} n^{2} \end{array}\right], $$ (27) 这里,ai(s)和bi(s)的表达式为
$$a_{i}(s)=\frac{\bar{h}_{21} \bar{h}_{13}-\bar{h}_{11} \bar{h}_{23}}{\bar{h}_{12} \bar{h}_{23}-\bar{h}_{13} \bar{h}_{22}},\ b_{i}(s)=\frac{\bar{h}_{21} \bar{h}_{12}-\bar{h}_{11} \bar{h}_{22}}{\bar{h}_{13} \bar{h}_{22}-\bar{h}_{12} \bar{h}_{23}}, $$ (28) 其中,hij(i=1, 2;j=1, 2, 3)代表矩阵[hij]第i行j列的元素.
1.2 接触压力分布
本小节中,我们将给出压电半空间与球压头尺度依赖的轴对称无摩擦接触问题的接触压力分布. 假设p(r)是真实的接触压力分布,那么表面法向位移满足如下的方程:
$$u_{z}(r, 0)-u_{z}(0, 0)-f(r)=0, $$ (29) 其中,f(r)是变形前压头的表面轮廓函数. 对于球压头来说,当接触区远小于压头半径时,f(r)可以用抛物面近似表示,即
$$f(r)=\frac{r^{2}}{2 R} . $$ (30) 如果p(r)是近似的接触压力分布,那么
$$u_{z}(r, 0)-u_{z}(0, 0)-f(r) \neq 0.$$ (31) 由于近似带来的误差记作[30]
$$e=\int_{0}^{a}\left[u_{z}(r, 0)-u_{z}(0, 0)-f(r)\right]^{2} r \mathrm{d} r.$$ (32) 接触压力分布的近似表达式可以写成级数的形式,即
$$p(r)=\sum\limits_{i=1}^{M} q_{i} p_{i}(r), $$ (33) 其中qi是权重;M是基函数的个数; pi(r)是基函数,可以表示为
$$p_{1}(r)=\sqrt{1-r^{2} / a^{2}}, $$ (34) $$p_{i}(r)=\mathrm{J}_{0}\left(\eta_{i-1} r\right), \quad i>1, $$ (35) 这里,J0(·)是第一类零阶Bessel函数;ηi是J0(ηa)=0的第i个零点.
将pi(r)单独作用导致的相对法向表面位移记作δVi(r)=uzi(r, 0)-uzi(0, 0),根据式(33),表面相对位移为
$${\rm{ \mathsf{ δ}}} V(r)=\sum\limits_{i=1}^{M} q_{i} {\rm{ \mathsf{ δ}}} V_{i}(r).$$ (36) 未知数qi需要通过积分最小二乘误差法确定,即
$$\frac{\partial e}{\partial q_{i}}=0, \quad i=1, 2, \cdots, M.$$ (37) 将式(37)表示为矩阵的形式,即
$$ \left[\begin{array}{ccc} K_{11} & \cdots & K_{1 M} \\ \vdots & & \vdots \\ K_{M 1} & \cdots & K_{M M} \end{array}\right]\left[\begin{array}{c} q_{1} \\ \vdots \\ q_{M} \end{array}\right]=\left[\begin{array}{c} f_{1} \\ \vdots \\ f_{M} \end{array}\right], $$ (38) 其中
$$K_{i j}=\int_{0}^{a} {\rm{ \mathsf{ δ}}} V_{i} {\rm{ \mathsf{ δ}}} V_{j} r \mathrm{d} r, \quad i=1, 2, \cdots, M ; j=1, 2, \cdots, M, $$ (39) $$f_{i}=\int_{0}^{a} {\rm{ \mathsf{ δ}}} V_{i} f(r) r \mathrm{d} r . $$ (40) 注意到,接触半径a和权重qi由上面的迭代过程确定,接触半径的初始值取经典弹性接触下的结果,即$ a_{\mathrm{c}}=$\sqrt[1 / 3]{3 F R f_{21} / 8}$,f21是与压电材料常数相关的系数. 相对误差取作0.001,即收敛准则为εerr≤0.001,其中
$$\varepsilon_{\text {err }}=\left|F-P^{\prime}\right| / F, $$ (41) 这里
$$P^{\prime}=\int_{0}^{a} 2 {\rm{ \mathsf{ π}}} r p(r) \mathrm{d} r .$$ (42) 至此,接触压力分布的式(33)中包含的系数Ai(s)仍然未知,所以需要通过边界条件进一步确定.
1.3 边界条件
在变换域中,接触表面的边界条件为
$$\hat{\sigma}_{z z}=-\hat{p}(s), \hat{\sigma}_{z r}=0, \hat{D}_{z}=0, \hat{m}_{z \theta}=0, $$ (43) 其中
$$\hat{p}(s)=\hat{p}_{1}(s)+\sum\limits_{j=2}^{m} \hat{p}_{j}(s), $$ (44) 并且
$$ \hat{p}_{1}(s)=q_{1} \frac{\sin (a s)-a s \cos (a s)}{a s^{3}}, \hat{p}_{j}(s)=q_{j} \frac{a\left[\eta_{j} \mathrm{J}_{0}(a s) \mathrm{J}_{1}(a s)-s \mathrm{J}_{0}(a s) \mathrm{J}_{1}(a s)\right]}{\eta_{j}^{2}-s^{2}} . $$ 1.4 接触问题的解
根据式(12)—(16)和(25),在变换域中,我们有
$$\left[\hat{\sigma}_{z z}, \hat{\sigma}_{r z}, \hat{D}_{z}, \hat{m}_{z \theta}\right]^{\mathrm{T}}=\boldsymbol{T}(s, z) \boldsymbol{A}(s), $$ (45) 其中
$$\boldsymbol{A}(s)=\left[A_{1}(s), A_{2}(s), A_{3}(s), A_{4}(s)\right]^{\mathrm{T}}, $$ (46) $$ \boldsymbol{T}(s, z)=\left[T_{1 i}(s, z), T_{2 i}(s, z), T_{3 i}(s, z), T_{4 i}(s, z)\right]^{\mathrm{T}}, \quad i=1, 2, 3, 4, $$ (47) 这里
$$ \begin{aligned} & T_{1 i}=\left(C_{13} s+C_{33} n_{i} a_{i}+e_{33} n_{i} b_{i}\right) \mathrm{e}^{n_{i z}}, \\ & T_{2 i}=\left[C_{44} n_{i}-C_{44} s-e_{15} s-G l^{2}\left(-s^{2} n_{i}+n_{i}^{3}-s^{3} a_{i}+s n_{i}^{2} a_{i}\right)\right] \mathrm{e}^{n_{i} z}, \\ & T_{3 i}=\left(e_{31} s+e_{33} n_{i} a_{i}-\varepsilon_{33} n_{i} b_{i}\right) \mathrm{e}^{n_{i} z}, T_{4 i}=\left(2 G n_{i}^{2} l^{2}+2 G s a_{i} n_{i} l^{2}\right) \mathrm{e}^{n_{i} z} . \end{aligned} $$ 根据边界条件(43),有
$$\boldsymbol{T}(s, h) \boldsymbol{A}(s)=\boldsymbol{L}, $$ (48) 其中
$$ \boldsymbol{L}=[-\hat{p}(s), 0, 0, 0]^{\mathrm{T}}, $$ (49) 那么,未知系数的表达式为
$$ \boldsymbol{A}=\boldsymbol{T}^{-1} \boldsymbol{L}.$$ (50) 因此,根据式(50)并对式(45)进行Hankel积分逆变换,我们就可以得到该接触问题应力场和偶应力场的解.
2. 结果与讨论
本文研究了刚性绝缘球压头和横观各向同性压电半空间尺度依赖的轴对称接触问题. 选择的压电材料的材料属性列于表 1中[31]. 除非特别说明,材料选用PZT-4,法向加载力、球压头半径和无量纲特征材料长度分别取F=10 N,R=500 μm,l/a=1.
C11/GPa C12/GPa C13/GPa C33/GPa C44/GPa e13/(C/m2) e33/(C/m2) e15/(C/m2) PZT-4 139 77.8 74.3 115 25.6 -5.2 15.1 12.7 BaTiO3 150 6.6 6.6 14.6 4.4 -17.3 17.5 11.4 ε11/(10-10·C/(V·m)) ε33/(10-10·C/(V·m)) f21/(10-11·Pa-1) PZT-4 64.61 56.2 1.772 8 BaTiO3 98.7 112 1.376 6 2.1 收敛性分析
表 2给出了接触半径a的收敛性分析. 在规定的误差范围内(εerr≤0.001),当M=3时,接触半径a的值不再发生变化,故取基函数的个数为3.
表 2 接触半径a的收敛性分析Table 2. Convergence analysis of contact radius aM 2 3 4 a/μm 23.671 23.665 23.665 2.2 对比算例
为了验证结果的有效性,首先将压电材料简化为各向同性弹性材料,当前的接触问题可以化简为一个刚性球压头和弹性半空间之间尺度依赖的轴对称无摩擦接触问题,Gourgiotis等[8]对该接触问题进行了研究. 取F=10 N,R=100 μm,μ=27.3 GPa和ν=0.3时,图 2给出了该问题的接触压力分布. 通过对比发现,利用本文方法得到的结果与Gourgiotis等[8]的研究结果一致,说明本文的分析方法是正确的.
2.3 特征材料长度的影响
图 3给出了特征材料长度对接触压力分布的影响. 从图中可以看出,在偶应力理论的框架内,球压头作用下的接触压力分布是光滑的,无奇异性出现,其最大值出现在接触区域的中心. 考虑了偶应力的接触压力结果明显大于经典弹性理论的结果(l/a=0). 并且随着特征材料长度的增加,接触压力显著增大,但是增加的速率减小.
图 4给出了不同特征材料长度下,法向加载F对最大接触压力的影响. 对于给定的最大接触压力,随着特征材料长度的增加,所对应的法向加载减小. 图 5描述了不同特征材料长度下,法向加载F对压痕深度的影响,可以看到力-位移曲线是非线性的. 在微纳米压痕实验中,通过分析加载和卸载过程中的力-位移曲线,可以确定材料的微机械特性,如硬度和弹性模量. 这些数据对于了解材料的抗压强度、耐磨性和其他重要特性至关重要. 从图中可以看出,对于较大的特征材料长度,需要更大的法向载荷才能产生相同的压入深度.
图 6描述了球压头作用下接触区内不同位置上,特征材料长度对法向接触压力比值p(r)/pc(r)的影响. 其中,pc(r)表示没有尺度效应的经典无摩擦接触问题的接触压力. 不同接触点处这一比值与1的偏差程度可以说明这一接触点处尺度效应对接触压力的影响. 从图中我们发现比值p(r)/pc(r)始终大于1. 也就是说,在接触区内,经典弹性理论下的法向接触压力始终小于考虑了尺度效应的结果. 随着特征材料长度的增加,p(r)/pc(r)显著增大,最终趋近于一个常数,这个常数只与材料参数有关. 对于弹性材料尺度依赖的接触问题,也得到了类似的结论[4].
PZT-4和BaTiO3是典型的压电陶瓷,其压电参数列于表 1中. PZT-4的压电效应显著;BaTiO3具有高介电常数,既有压电性也有铁电性. 选用这两种材料,当F=10 N,R=500 μm时,图 7讨论了特征材料长度对接触半径a/ac、最大接触压力p0/p0c和压痕深度δ/δc的影响. 其中,ac,p0c和δc是经典弹性理论中的接触半径、最大接触压力和压痕深度. 从图中可以看到,经典弹性理论的接触半径和最大压入深度始终大于考虑了尺度效应的解,而经典弹性理论的最大接触压力始终小于考虑了尺度效应的解. 当法向载荷不变时,随着特征材料长度的增加,接触半径和压痕深度减小,最大接触压力增加. 接触半径、最大接触压力和压痕深度开始随特征材料长度的增加而发生明显变化,随着特征材料长度增加到一定程度时,接触半径、最大接触压力和压痕深度的变化趋于平缓. 当特征材料长度一定时,PZT-4所对应的a/ac和δ/δc的值要小于BaTiO3,但是PZT-4所对应的p0/p0c要大于BaTiO3.
3. 结论
利用偶应力理论,本文研究了压电半空间和刚性绝缘球压头的尺度依赖的轴对称无摩擦接触问题. 利用Hankel积分变换和积分最小二乘法,得到了接触压力分布,获得了如下的结论:
1) 经典弹性理论下的接触压力始终小于考虑了尺度效应的解,而经典弹性理论下的接触半径和压痕深度始终大于考虑了尺度效应的解;
2) 随着特征材料长度的增加,法向接触压力分布明显增加,但增加速率减缓;
3) 法向载荷一定时,随着特征材料长度的增加,最大法向接触压力显著增加,接触半径和压痕深度显著减小.
-
C11/GPa C12/GPa C13/GPa C33/GPa C44/GPa e13/(C/m2) e33/(C/m2) e15/(C/m2) PZT-4 139 77.8 74.3 115 25.6 -5.2 15.1 12.7 BaTiO3 150 6.6 6.6 14.6 4.4 -17.3 17.5 11.4 ε11/(10-10·C/(V·m)) ε33/(10-10·C/(V·m)) f21/(10-11·Pa-1) PZT-4 64.61 56.2 1.772 8 BaTiO3 98.7 112 1.376 6 表 2 接触半径a的收敛性分析
Table 2. Convergence analysis of contact radius a
M 2 3 4 a/μm 23.671 23.665 23.665 -
[1] STELMASHENKO N A, WALLS M G, BROWN L M, et al. Microindentations on W and Mo oriented single crystals: an STM study[J]. Acta Metallurgica et Materialia, 1993, 41(10): 2855-2865. doi: 10.1016/0956-7151(93)90100-7 [2] NIX W D, GAO H. Indentation size effects in crystalline materials: a law for strain gradient plasticity[J]. Journal of the Mechanics and Physics of Solids, 1998, 46(3): 411-425. doi: 10.1016/S0022-5096(97)00086-0 [3] MILLER R E, SHENOY V B. Size-dependent elastic properties of nanosized structural elements[J]. Nanotechnology, 2000, 11(3): 139-147. doi: 10.1088/0957-4484/11/3/301 [4] ZISIS T, GOURGIOTIS P A, BAXEVANAKIS K P, et al. Some basic contact problems in couple stress elasticity[J]. International Journal of Solids and Structures, 2014, 51(11/12): 2084-2095. [5] GOURGIOTIS P, ZISIS T. Two-dimensional indentation of microstructured solids characterized by couple-stress elasticity[J]. The Journal of Strain Analysis for Engineering Design, 2016, 51(4): 318-331. doi: 10.1177/0309324715611524 [6] ZISIS T. Anti-plane loading of microstructured materials in the context of couple stress theory of elasticity: half-planes and layers[J]. Archive of Applied Mechanics, 2018, 88(1): 97-110. [7] ZISIS T. Burmister's problem extended to a microstructuredlayer[J]. Journal of Mechanics of Materials and Structures, 2018, 13(2): 203-223. doi: 10.2140/jomms.2018.13.203 [8] GOURGIOTIS P A, ZISIS T, GIANNAKOPOULOS A E, et al. The Hertz contact problem in couple-stress elasticity[J]. International Journal of Solids and Structures, 2019, 168: 228-237. doi: 10.1016/j.ijsolstr.2019.03.032 [9] KARURIYA A N, BHANDAKKAR T K. Plane strain indentation on finite thickness bonded layer in couple stress elasticity[J]. International Journal of Solids and Structures, 2017, 108: 275-288. doi: 10.1016/j.ijsolstr.2016.12.027 [10] SONG H X, KE L L, WANG Y S. Sliding frictional contact analysis of an elastic solid with couple stresses[J]. International Journal of Mechanical Sciences, 2017, 133: 804-816. doi: 10.1016/j.ijmecsci.2017.09.037 [11] SONG H X, KE L, WANG Y, et al. Two-dimensional frictionless contact of a coated half-plane based on couple stress theory[J]. International Journal of Applied Mechanics, 2018, 10(5): 1850049. doi: 10.1142/S1758825118500497 [12] WANG Y, SHEN H, ZHANG X, et al. Semi-analytical study of microscopic two-dimensional partial slip contact problem within the framework of couple stress elasticity: cylindrical indenter[J]. International Journal of Solids and Structures, 2018, 138: 76-86. doi: 10.1016/j.ijsolstr.2017.12.030 [13] ÇÖMEZ I, EL-BORGI S. Sliding frictional contact problem of a layer indented by a rigid punch in couple stress elasticity[J]. Mathematics and Mechanics of Solids, 2023, 28(3): 730-747. doi: 10.1177/10812865221080551 [14] NABHANI M, EL KHLIFI M, GBEHE O S T, et al. Coupled couple stress and surface roughness effects on elasto-hydrodynamic contact[J]. Lubrication Science, 2014, 26(4): 251-271. doi: 10.1002/ls.1246 [15] WANG Y X, ZHANG X, SHEN H, et al. Three-dimensional contact analysis with couple stress elasticity[J]. International Journal of Mechanical Sciences, 2019, 153: 369-379. [16] WANG Y X, ZHANG X, SHEN H, et al. Couple stress-based 3D contact of elastic films[J]. International Journal of Solids and Structures, 2020, 191: 449-463. [17] LI P X, LIU T J. Axisymmetric adhesive contact of multi-layer couple-stress elastic structures involving graded nanostructured materials[J]. Applied Mathematical Modelling, 2022, 111: 501-520. doi: 10.1016/j.apm.2022.06.044 [18] ZHOU Y T, TIAN X J, DING S H. Microstructure size-dependent contact behavior of a thermoelectric film bonded to an elastic substrate with couple stress theory[J]. International Journal of Solids and Structures, 2022, 256: 111982. doi: 10.1016/j.ijsolstr.2022.111982 [19] LU R, LI M H, YANG Y, et al. Accurate extraction of large electromechanical coupling in piezoelectric MEMS resonators[J]. Journal of Microelectromechanical Systems, 2019, 28(2): 209-218. doi: 10.1109/JMEMS.2019.2892708 [20] WANG Z L, SONG J. Piezoelectric nanogenerators based on zinc oxide nanowire arrays[J]. Science, 2006, 312(5771): 242-246. doi: 10.1126/science.1124005 [21] MIRZAEI A, LEE J H, MAJHI S M, et al. Resistive gas sensors based on metal-oxide nanowires[J]. Journal of Applied Physics, 2019, 126(24): 241102. doi: 10.1063/1.5118805 [22] SAHU S A, SINGHAL A, CHAUDHARY S. Surface wave propagation in functionally graded piezoelectric material: an analytical solution[J]. Journal of Intelligent Material Systems and Structures, 2018, 29(3): 423-437. doi: 10.1177/1045389X17708047 [23] LI Y S, PAN E. Static bending and free vibration of a functionally graded piezoelectric microplate based on the modified couple-stress theory[J]. International Journal of Engineering Science, 2015, 97: 40-59. doi: 10.1016/j.ijengsci.2015.08.009 [24] RAZAVI H, BABADI A F, BENI Y T. Free vibration analysis of functionally graded piezoelectric cylindrical nanoshell based on consistent couple stress theory[J]. Composite Structures, 2017, 160: 1299-1309. doi: 10.1016/j.compstruct.2016.10.056 [25] WANG X, PAN E, FENG W J. Anti-plane Green's functions and cracks for piezoelectric material with couple stress and electric field gradient effects[J]. European Journal of Mechanics A: Solids, 2008, 27(3): 478-486. doi: 10.1016/j.euromechsol.2007.09.005 [26] SUN Y Y, SU J, SONG H X, et al. The size-dependent frictionless contact of piezoelectric materials[J]. International Journal of Mechanical Sciences, 2024, 261: 108685. doi: 10.1016/j.ijmecsci.2023.108685 [27] KE L L, YANG J, KITIPORNCHAI S, et al. Electro-mechanical frictionless contact behavior of a functionally graded piezoelectric layered half-plane under a rigid punch[J]. International Journal of Solids and Structures, 2008, 45(11/12): 3313-3333. [28] SU J, KE L L, WANG Y S. Axisymmetric frictionless contact of a functionally graded piezoelectric layered half-space under a conducting punch[J]. International Journal of Solids and Structures, 2016, 90: 45-59. doi: 10.1016/j.ijsolstr.2016.04.011 [29] CHEN W T, ENGEL P A. Impact and contact stress analysis in multilayer media[J]. International Journal of Solids and Structures, 1972, 8(11): 1257-1281. doi: 10.1016/0020-7683(72)90079-0 [30] O'SULLIVAN T C, KING R B. Sliding contact stress field due to a spherical indenter on a layered elastic half-space[J]. Journal of Tribology, 1988, 110(2): 235-240. doi: 10.1115/1.3261591 [31] LV X, KE L L, SU J, et al. Axisymmetric contact vibration analysis of a rigid spherical punch on a piezoelectric half-space[J]. International Journal of Solids and Structures, 2021, 210: 224-236. [32] WANG G F, YU S W, FENG X Q. A piezoelectric constitutive theory with rotation gradient effects[J]. European Journal of Mechanics A: Solids, 2004, 23(3): 455-466. doi: 10.1016/j.euromechsol.2003.12.005 [33] MINDLIN R D. Influence of couple-stresses on stress concentrations[J]. Experimental Mechanics, 1963, 3(1): 1-7. doi: 10.1007/BF02327219 [34] LIU C, YU J, ZHANG B, et al. Reflection and transmission of elastic waves in the multilayered orthotropic couple-stressed plates sandwiched between two elastic half-spaces[J]. Applied Mathematical Modelling, 2019, 75: 52-72. doi: 10.1016/j.apm.2019.05.023 -