多孔功能梯度压电纳米壳中波传播特性.pdf

应用数学和力学编委会，：多孔功能梯度压电纳米壳中波传播特性王鑫特，刘 娟，胡 彪，张 波，沈火明（西南交通大学 力学与航空航天学院，成都）摘要：基于非局部应变梯度理论，探究了含孔隙的功能梯度压电陶瓷纳米壳中波传播特性 利用 原理和一阶剪切理论推导了控制方程 结合非局部应变梯度理论和谐波解得到了尺度依赖的特征方程 数值讨论了尺度参数、波数、梯度指数、壳厚、孔隙率及电压对波传播特性的影响 研究表明：非局部参数和应变梯度参数对波传播频率的影响与波数密切相关，在一定范围内波数越大，尺度参数对频率的影响越大；另外，孔隙和梯度指数对频率具有耦合作用 关 键 词：压电纳米壳；波动特性；非局部应变梯度理论；功能梯度中图分类号：；文献标志码：，（，）：，：；引 言随着新型智能器件的高速发展，开发各类多功能材料已成为一种趋势 其中，功能梯度压电材料（应用数学和力学 卷 期 年 月 ，收稿日期：；修订日期：基金项目：国家自然科学基金（）作者简介：王鑫特（），男，硕士（：）；刘娟（），女，副教授，博士（通讯作者：）引用格式：王鑫特，刘娟，胡彪，张波，沈火明 多孔功能梯度压电纳米壳中波传播特性 应用数学和力学，（）：，）是采用先进的材料复合技术将两种或多种不同材料耦合而成的非均质材料 的构成要素（组份、组织、显微气孔率等）随空间位置呈梯度性连续变化，能充分减少和克服组份材料结合部位界面性能的不匹配因素，具有良好的抗变形能力、机电转换性及耐腐蚀性 在 制造过程中可通过调控两种或两种以上材料组分的体积分数、孔隙率等参数来设计材料性能，以满足工程结构特定的功能要求 随着智能器件逐渐向微型化方向发展，的应用范围逐渐延伸至纳机电系统领域 纳机电系统中承载构件可以简化为特征尺度处于纳米量级的梁、板、壳等结构，典型应用为纳米谐振器、纳米传感器、纳米执行器、纳米开关等 随着构件特征尺寸减小至纳米尺度，原子之间长程作用力会显著增加，尺度效应变得不可忽略 分子动力学模拟表明，经典连续介质力学理论所预测的纳米结构振动频率与实测结果相差 因此，寻找和发展适应于纳米力学研究的新途径和新方法，是当前连续力学的研究热点之一 由 提出的非局部弹性理论在纳米力学研究中扮演着重要角色，非局部连续力学将分子间作用力是长程力的思想引入到传统连续介质力学中，认为连续体内某一点的应力不仅取决于该点的应变，还是连续体内所有点的应变及应变梯度的函数 非局部理论被诸多学者证明能成功地预测软化效应 然而 等研究了含流体的双壁碳纳米管的非线性振动，表明非局部弹性模型无法预测材料可能存在的刚度硬化效应 而且在不同金属材料的微纳米压痕实验中，同样观察到了非局部理论无法解释的微尺度下材料强度比常规尺度下材料强度显著提高的现象 而应变梯度假设则是将连续体中的每一个物质点看作含有高阶应变的胞元，据此引入长度尺度参数来表征其对结构力学性能的影响，该参数可合理地预测结构中的刚度硬化效应 等基于非局部理论并结合应变梯度假设，提出了新的非局部应变梯度理论 该理论考虑了应变梯度的非局部效应，也考虑了材料物质点处高阶应力梯度的非局部效应 等将该理论应用于纳米梁和碳纳米管中的波传播分析中，揭示了关于晶格动力学和波传播实验的一些新发现 等通过两种非局部应变梯度壳理论，研究了磁电弹性纳米壳中的波传播特性，总结了电磁机械荷载与截止波数的关系 和他的同事利用非局部应变梯度理论，对功能分级纳米板在弯曲、屈曲和轴向运动时的振动特性和稳定性进行了系统研究，并且在这些机械运动中都观察到了刚度软化和硬化机制 此外，不少学者也利用非局部应变梯度理论探究了 纳米结构的动态特性 然而，以上研究对孔隙作用下 纳米壳的波动特性关注较少 本文采用基于非局部应变梯度理论的一阶剪切壳模型研究了多孔 纳米材料的波传播特性 根据 原理推导了结构的非局部运动微分方程，进而通过数值分析研究了尺度参数、几何参数、内部结构参数、外部电压等参数对波传播特性的影响 多孔 纳米壳的动力学建模 非局部应变梯度理论由 等建立的非局部应变梯度理论综合考虑了应力场与高阶应力场的非局部性 其应力场可表示为 （）（），（）其中，（）是经典非局部应力，（）是高阶非局部应力，是 算子 该理论通过引入两个尺度参数，同时涉及了应变和应变梯度的非局部效应，式（）中的两个应力可由尺度参数及对应应变表示为（）（，）（），（）（）（，），（），（）式中，是表征高阶应变梯度效应的材料特征长度参数，与 是非局部核函数，和 是表征非局部效应的非局部参数，是弹性系数，和，分别是应变和应变梯度 将式（）、（）代入式（）中，引入 算子，可得到一种简化应力和非局部应力的关系：（）（）（）（）（）经典非局部应变梯度理论假设 ，并忽略高阶项（），即可得到基于该理论的本构方程为 （），（）应 用 数 学 和 力 学 年 第 卷其中，（），在电场的影响下，纳米壳结构的矩阵形式本构方程可以定义为（）（），（）（）（），（）其中，为材料的弹性系数、压电系数及介电系数，为纳米壳任意一点处 方向的电位移 多孔 纳米壳的几何描述本文采用由压电陶瓷 与 耦合组成的多孔 纳米壳模型，如图 所示 壳长度为，半径为，厚度为，图中（，），（，）和（，）是壳上任意一点处在，三个方向上的位移 壳内部含有的孔隙分布形式可看做是均匀分布，纳米壳的材料常数可表示为（）（）（），（）其中，为纳米壳中任意一点到中表面的垂直距离，参数，）表示幂律指数，即功能梯度指数，参数 为孔隙率 为了满足 方程，纳米壳外部沿厚度方向分布的电势定义为（，）（）（，），（）其中，为一个线性常数，（，）表示壳平面上分布的电势，代表初始的外加电压 图 多孔 纳米壳的模型示意图 利用一阶剪切变形壳模型方程分析位移场，得到壳中任意一点处的位移分量为（，）（，）（，）（，）（，）（，）（，）（，）（）纳米壳任意一点处的电位移为 ，（），（），（），（），（）（）第 期 王鑫特，等：多孔功能梯度压电纳米壳中波传播特性 多孔 纳米壳的运动微分方程 纳米壳上各个点处的动能 为 （）模型应变能 为（）（）轴向和环向电场力对 纳米壳做功为 （）将式（）（）代入 原理（），（）可得如下控制微分方程：，（），（），（），（），（）（）（）（），（）其中（，）（，），（）（，）（，），（）（，）（，），（）（，）（，）（）将非局部应变梯度本构方程（）沿纳米壳径向进行积分，并将所得力和力矩表达式代入控制方程（）（）中，可得位移形式的结构运动微分方程如下：，（），（）应 用 数 学 和 力 学 年 第 卷 ，（）（），（），（）（），（）其中，（），（）多孔 纳米壳波传播问题的求解针对 纳米壳内每一点处波传播响应的模拟问题，引入壳结构波传播通解：（，）（，）（，）（，）（，）（，）（），（）其中，分别表示纳米壳各个点处波传播的位移与转角的幅值，表示电势幅值，为虚数单位，和 表示波传播时沿着 方向产生的纵向波数和沿着转角 方向产生的周向波数 将通解方程（）代入运动微分方程（）（）中，即可得到 纳米壳的波传播特征方程为（），（）其中，与 分别为刚度系数矩阵以及质量系数矩阵，角频率 的值可由系数矩阵的行列式等于零求得，而波传播频率 的值可根据 与 的关系式得到：（）第 期 王鑫特，等：多孔功能梯度压电纳米壳中波传播特性 数值结果及分析本节通过求解特征方程（），研究了波数、尺度参数、几何参数、外部电势、孔隙率以及功能梯度指数对图 所示 圆柱形纳米壳波传播特性的影响，两种材料参数如表 所示 表 的材料特性数值 ，（），有效性验证为了检验理论模型的正确性，将 纳米壳的结果和磁电弹纳米壳的结果进行验证 基于非局部理论将本文得到的波传播频散关系与 等的结果（文献中的图（）进行对比（图）可以看出，本研究的结果与 等的结果良好吻合，验证了本研究的有效性 图 退化验证（）（）注 为了解释图中的颜色，读者可以参考本文的电子网页版本，后同 数值讨论表 表示前三阶固有频率（，）与纵向波数和周向波数的关系 从表中可以看出三阶频率随着纵向及周向波数的增加而增加，且两方向波数对频率的增加影响呈现互补效果 本文后续工作主要分析探究纳米壳中波传播的纵向波数 对基频 的影响 图 显示了在不同波数 下，两个尺度参数间比值 的变化对 纳米壳波传播频率 的影响 图中 ，从图中可以看出，整体上，当 时，非局部应变梯度理论对应的频率都小于通过经典弹性理论得到的频率，这是因为非局部效应占主导地位时，会使得壳的刚度减弱，呈现软化效应 相反地，当 时，非局部应变梯度理论获得的频率都大于通过经典弹性理论得到的频率，此时占主导地位的应变梯度效应使得壳的刚度增强，呈现硬化效应 当 时，非局部应变梯度理论和经典弹性理论得到的频率是等价的，说明两种类型的尺寸依赖效应会相互抵消 此外，固定长度尺度参数 不变，在不同波数下同等程度的改变非局部参数，发现当波数较大时导致的频移大于波应 用 数 学 和 力 学 年 第 卷数较小时对应的频移 同理，固定 不变时，对 的影响也随着 的增大而增大 这表明由非局部参数引起的软化效应和由应变梯度参数引起的硬化效应在波数较大时更加明显，在波数较小时较为微弱 所以对于纳米壳尺度效应的频散特性研究不仅需要考虑尺度效应的作用，也需要关注波数的影响 表 前三阶固有频率（，）与纵向波数 和周向波数 的关系 （）（）（）图 基于不同波数，变化的 比值对 纳米壳波传播频率的影响 图 为基于非局部应变梯度理论，在不同波数 下，变化的功能梯度指数 对 纳米壳波传播频率的影响 图中 ，在图 中也可以观察到图 中的类似现象 这是由于当 和 时，非局部参数的软化效应和应变梯度参数的硬化效应得到了体现 特别地，当两个参数取值不同但比值相同时，其对应的 曲线在波数较小时间距很大，而在波数较大时近似重合 这也证明第 期 王鑫特，等：多孔功能梯度压电纳米壳中波传播特性了纳米结构中的波传播行为研究需要考虑尺度效应和波数的共同作用 此外，还可以看出，波传播频率 随着功能梯度指数 的增加而减小 这是因为在由 与 耦合组成的 纳米壳中，当 增大时，材料的整体性能 会更接近于，而 的弹性模量小于，从而整个纳米壳的弹性模量会减小，与之相对应的频率也会减小 且通过对比发现，频率 受梯度指数 的影响程度随着波数的增加而增加，也即纳米壳频率受到功能梯度指数和波数的耦合作用（）（）（）图 基于不同波数，变化的尺度参数和功能梯度指数对 纳米壳波传播频率的影响 图 探讨了对于不同的壳厚，纳米壳的频率 与功能梯度指数 间的关系 图中 ，从图中可以看出，波传播频率 随着壳厚 的增加而降低，且在壳厚 越小时，频率对壳厚 变化越敏感 这是因为壳厚变大会导致结构刚度和质量的同时增加，而对于由 和 复合而成的 纳米壳，在一定范围内，由质量增加导致波传播频率的降低程度比由刚度增加导致频率的增大程度更大 因此，更大的厚度会导致更低的频率，而且这种影响程度随着壳厚的减小而更加显著 图 中给出了孔隙率 和功能梯度指数 对 纳米壳波传播频率 的影响 图中 ，由图可知，孔隙率 和功能梯度指数 对波传播频率有耦合影响，当 从 开始小范围内增加时，纳米壳的频率 会随着孔隙率 的增加而增加 当 时，规律则会相反，即频率 随着孔隙率 的增加而减少 这是因为孔隙率对材料刚度的影响与此时的功能梯度指数有关 当 时，减小，材料的刚度在增加，而 的增大会减小材料的刚度，但此时 较小，刚度较大的材料 占比很大，即在这一耦合作用中，由 减小导致增加的材料刚度大于由 增大而减小的材料刚度，所以频率会增大 当 时，与 的增加同时导致材料刚度的减少，从而导致频率快速减小 这说明对于两种材料耦合组成的多孔，可以通过改变结构的孔隙率和功能梯度指数来控制波在该应 用 数 学 和 力 学 年 第 卷结构中的传播 图 为在不同外加电势作用下，纳米壳的频率与功能梯度指数的函数关系图 图中 ，由图可知，波传播频率会随着正向电压增大而减小，而随反向电压的增大而增大 这是由于在对模型结构施加正 反向电势时，会在结构的轴向和周向产生压 拉力，使得结构的刚度发生软 硬化效应，从而减少 增加波传播频率 特别地，从图中可看出，频率对反向电压变化的敏感度大于对正向电压的敏感度 此外，还可观察到频率受正向与反向电压的影响均随着功能梯度指数的增加而变得明显 这是因为 越大，材料 占比越多，而 的压电系数大于，与之相对应的电压对频率的影响程度也会增加 图 不同壳厚下 纳米壳的频率与功能梯度指数的关系 图 不同孔隙率下 纳米壳的频率图 不同外电压作用下 纳米壳的与功能梯度指数的关系频率与功能梯度指数的关系 结 论本文基于非局部应变梯度理论建立了多孔 纳米壳的尺度依赖模型 数值分析了尺度参数、梯度指数、孔隙率、壳厚和电压对波传播频率的影响 主要结论如下：）当尺度参数间比值 时，会使材料刚度呈现软化效应，而 时，则呈现硬化效应，另外，增大波数或减小壳厚会加剧尺度效应；）增大功能梯度指数或纳米壳厚度会减小波传播频率，其中，波数越大，梯度指数的影响程度越大；）孔隙和功能梯度指数对频率具有耦合作用，功能梯度指数较小时，功能梯度指数和孔隙对频率的作用总体表现为增强，反之表现为降低；第 期 王鑫特，等：多孔功能梯度压电纳米壳中波传播特性）对结构施加正向电压，波传播频率会减小，施加反向电压，频率会增大，且电压对频率的影响程度随着功能梯度指数的增加而增加 参考文献（）：夏巍，冯浩成 热过屈曲功能梯度壁板的气动弹性颤振 力学学报，（）：（，（）：（），（）：，：陈明飞，刘坤鹏，靳国永，等 面内功能梯度三角形板等几何面内振动分析 应用数学和力学，（）：（，（）：（），：王平远，李成，姚林泉 基于非局部应变梯度理论功能梯度纳米板的弯曲和屈曲研究 应用数学和力学，（）：（，（）：（），（）：，（）：，（）：，（）：，（）：，（）：，（）：徐晓建，邓子辰 基于简化的应变梯度理论下 板模型边值问题的提法及其应用 应用数学和力学，（）：（，（）：（），：，：，应 用 数 学 和 力 学 年 第 卷，（）：，（）：，：，：（），（）：，（）：，（），（）：，（）：，：，（）：，：，：，（）：，第 期 王鑫特，等：多孔功能梯度压电纳米壳中波传播特性