偏场作用下软电弹性圆柱中导波传播特性的深度剖析与研究

偏场作用下软电弹性圆柱中导波传播特性的深度剖析与研究一、引言1.1软电弹性材料概述1.1.1工作原理软电弹性材料作为智能材料的重要一员，其工作原理基于独特的电致伸缩效应和麦克斯韦应力作用。在微观层面，当电场施加于软电弹性材料时，材料内部的分子或聚合物链会发生重排与取向变化。以常见的介电弹性体为例，它由柔软的高分子聚合物基体和具有高介电常数的填料组成。在无电场作用时，分子链呈无序状态；施加电场后，材料内部产生电偶极子，在电场力作用下，电偶极子趋向于沿电场方向排列，进而使分子链发生取向，导致材料在宏观上产生形变。从能量转换角度看，电场能量促使分子链重排，克服分子间相互作用力，实现电能到机械能的转化。麦克斯韦应力在软电弹性材料的形变过程中发挥关键作用。根据麦克斯韦电磁理论，在电场存在的电介质中，会产生与电场强度平方成正比的麦克斯韦应力。当软电弹性材料处于电场中时，麦克斯韦应力会在材料内部产生张力或压力。在平行于电场方向，麦克斯韦应力促使材料收缩；在垂直于电场方向，材料则发生膨胀，这种应力分布差异导致材料产生复杂的形变。以平行板电容器结构模型来理解，将软电弹性材料置于平行板电极之间，施加电压后，材料所受麦克斯韦应力与电场强度、材料介电常数相关，通过调整这些参数，可有效调控材料的形变程度与方向。1.1.2应用领域软电弹性材料凭借独特的性能，在众多领域展现出巨大应用潜力，有力推动各领域的创新发展。在生物医学领域，软电弹性材料为疾病诊断与治疗带来新契机。在药物输送系统中，利用软电弹性材料制作的微胶囊或纳米颗粒，通过外部电场控制其形变，实现药物的精准释放。当这些载药微粒到达病变部位时，施加特定电场，材料形变打开释放通道，提高药物疗效并减少副作用。在神经刺激与修复方面，软电弹性材料制成的柔性电极能与神经组织紧密贴合，通过施加电场刺激神经细胞，促进神经功能恢复，为神经系统疾病治疗提供新手段。智能传感器领域，软电弹性材料大显身手。制作压力传感器时，基于材料在压力作用下的形变导致电容或电阻变化原理，将压力信号转化为电信号，实现对压力的精确测量。这种传感器具有高灵敏度、可拉伸性，能适应复杂环境与不规则表面测量，广泛应用于可穿戴设备、人机交互界面等。在应变传感器应用中，软电弹性材料可实时监测物体的应变状态，用于结构健康监测、运动检测等领域，为保障结构安全与运动分析提供重要数据支持。柔性机器人领域，软电弹性材料是核心驱动部件。软电弹性驱动器具有响应速度快、输出力大、变形能力强等优点，使机器人能够模仿生物的运动方式，实现更灵活、柔顺的动作。比如，用于制作仿生鱼的驱动鳍，通过控制电场强度和频率，实现仿生鱼在水中的灵活游动；制作的软体抓手，可根据抓取物体的形状自动调整形状，实现轻柔抓取，在医疗护理、危险环境作业等领域具有广阔应用前景。1.2导波的基本概念与分类1.2.1定义及特点导波是一种在固体表面及内部不同材料界面或固体与液体界面处传播的非均匀自持波。在固体中，由于边界条件的显著影响，导波的传播速度与不考虑边界影响的体波传播速度存在明显差异。从传播特性角度来看，体波在均匀无限介质中传播时，速度相对稳定，如纵波（P波）在均匀固体介质中以较快速度传播，其传播方向与质点振动方向一致；横波（S波）传播速度相对较慢，质点振动方向与传播方向垂直。而导波在波导结构（如管、板、棒等）中传播时，其速度受到波导的几何形状、尺寸以及材料特性等多种因素的综合影响。以在金属薄板中传播的导波为例，薄板的厚度变化会导致导波传播速度改变，当薄板厚度减小时，导波传播速度可能会发生相应变化，且在不同频率下，导波速度的变化规律也有所不同，呈现出复杂的频散特性。导波由一系列谐波叠加而成，这一特性使其具有传播距离长而衰减小的独特优势。在实际应用中，如长距离的石油天然气输送管道检测，利用导波进行无损检测时，导波能够沿着管道传播较远的距离，并且在传播过程中能量衰减相对较小，从而可以实现对管道较长距离范围内的缺陷检测，及时发现管道中的裂纹、腐蚀等缺陷，保障管道的安全运行。1.2.2常见类型在弹性动力学和固体波动力学中，导波包含多种类型，常见的有洛夫波、瑞利波、兰姆波等，它们各自具有独特的产生条件和传播特点。洛夫波（Lovewave）是一种在固体表面传播的导波，由英国数学家洛夫（A.E.H.Love）于1911年首先从理论上预测。其产生条件是在一个半无限弹性固体表面上存在一层低速介质。在传播过程中，洛夫波的质点振动方向与波的传播方向垂直，且平行于固体表面。洛夫波的传播速度介于横波速度和表面波波速之间，具有频散特性，即不同频率的洛夫波传播速度不同。在地震勘探中，洛夫波的检测与分析可以为研究地球浅层地质结构提供重要依据，通过分析洛夫波的频散曲线等特征，能够推断地下不同地层的弹性参数和厚度等信息。瑞利波（Rayleighwave）是由英国科学家瑞利（Rayleigh）在1885年发现的一种沿固体表面传播的弹性波。它的产生条件是在均匀各向同性的弹性半空间表面，当存在扰动时就会产生瑞利波。瑞利波传播时，质点在垂直于表面的平面内做椭圆运动，椭圆长轴垂直于表面，短轴平行于表面，其传播速度约为横波速度的0.92倍。瑞利波的能量主要集中在表面以下约一个波长的深度范围内，且具有较低的传播速度和较短的波长，这使得它在表面检测领域有着广泛应用。在材料表面缺陷检测中，利用瑞利波对材料表面的敏感性，当材料表面存在缺陷时，瑞利波的传播特性会发生变化，通过检测这些变化可以实现对表面缺陷的识别与定位。兰姆波（Lambwave）是1917年由H.兰姆建立数学模型的一种弹性导波，它在薄板状固体介质中传播。兰姆波的产生需要满足一定的边界条件，即在薄板的上下表面面力为零。它具有对称型和反对称型两种基本类型，对称型兰姆波又称纵向板波（或膨胀板波），反对称型兰姆波称为弯曲板波。兰姆波是一种频散波，不同频率的兰姆波具有不同的传播速度和振动模式，且存在高次波，每种高次波都有特定的截止频率，当频率达到截止频率时，板内才能传播该高次波。在薄板探伤领域，兰姆波发挥着重要作用，由于其能够在薄板中传播较远的距离，并且对薄板内部的缺陷较为敏感，可用于检测薄板中的分层、裂纹等缺陷，提高薄板质量检测的效率和准确性。1.3研究现状与意义1.3.1国内外研究进展软电弹性材料与导波特性的研究在国内外均取得了显著进展，成为众多学科交叉领域的研究热点。在软电弹性材料的基础理论研究方面，国内外学者深入探索材料的本构关系与电-力耦合机制。国外如德国卡尔斯鲁厄理工学院的研究团队，运用微观力学理论和分子动力学模拟方法，从分子层面揭示软电弹性材料在电场作用下的结构演变和力学响应规律，建立更为精准的本构模型，为材料性能优化提供理论依据。国内清华大学的研究人员通过实验与理论分析相结合，研究不同填料种类和含量对软电弹性材料电致伸缩性能的影响，发现添加特定纳米填料可有效提高材料的介电常数和电-力转换效率，为新型软电弹性材料的设计合成提供新思路。对于导波在软电弹性材料中的传播特性，国内外学者展开多维度研究。在理论分析方面，美国西北大学利用弹性动力学和电磁学理论，推导导波在软电弹性圆柱中的传播方程，分析频散特性和模态分布，为导波检测技术提供理论基础。国内西安交通大学采用半解析有限元法，对导波在复杂形状软电弹性结构中的传播进行数值模拟，深入研究导波与结构缺陷的相互作用机制，为结构健康监测提供技术支持。在实验研究方面，国外斯坦福大学搭建高精度的激光超声测量系统，对软电弹性圆柱中导波的传播速度、衰减特性等进行实验测量，验证理论和数值模拟结果的准确性。国内浙江大学利用电磁超声换能器激发导波，研究不同激励频率和波形对导波激发效率和模式的影响，为导波检测系统的优化设计提供实验依据。当前研究重点聚焦于开发多物理场耦合下的软电弹性材料本构模型，以更全面准确地描述材料在复杂工况下的性能；深入探究导波在复杂结构和多相介质中的传播特性，提高导波检测技术的精度和可靠性；开展软电弹性材料与导波技术在生物医学、智能结构等领域的应用基础研究，推动相关技术的工程化应用。研究热点集中在新型软电弹性材料的设计合成与性能调控、导波的高效激发与精准检测技术、多学科交叉融合下的软电弹性材料与导波技术的创新应用等方面。1.3.2研究意义对软电弹性圆柱中导波的研究具有重要的理论意义和广泛的应用价值，在多个领域展现出独特的作用。从理论层面看，软电弹性材料作为一种新型智能材料，其电-力耦合特性为导波传播研究带来新挑战与机遇。深入研究导波在软电弹性圆柱中的传播特性，有助于进一步完善弹性动力学和电磁学在智能材料领域的理论体系，揭示电-力-波相互作用的微观机制，为多物理场耦合问题的研究提供新思路和方法，推动学科交叉融合发展。在材料性能检测领域，导波检测技术凭借传播距离长、检测效率高的优势，成为软电弹性材料性能评估的有力手段。通过分析导波在软电弹性圆柱中的传播特性，如频散曲线、衰减规律等，可以获取材料的弹性模量、介电常数、电导率等关键性能参数，实现对材料内部缺陷、损伤的高精度检测与定位，为材料质量控制和性能优化提供重要依据，保障材料在实际应用中的可靠性和安全性。在结构健康监测领域，软电弹性材料制成的智能结构在航空航天、生物医学等领域得到广泛应用。利用导波对结构进行健康监测，能够实时获取结构的应力、应变状态和损伤信息，及时发现潜在安全隐患，为结构的维护、修复和寿命预测提供科学指导，提高结构的可靠性和使用寿命，降低维护成本。以航空航天领域的飞行器机翼为例，采用软电弹性材料作为结构材料，通过导波监测技术可以及时检测机翼在飞行过程中的损伤情况，确保飞行安全。在生物医学领域，软电弹性材料与导波技术的结合为生物组织成像和疾病诊断开辟新途径。导波在软电弹性材料与生物组织的界面传播时，其特性会因组织的生理状态和病变情况发生变化，通过检测这些变化可以实现对生物组织的无损检测和疾病早期诊断，如在肿瘤检测中，利用导波检测技术可以检测肿瘤组织的位置、大小和形态，为肿瘤的早期治疗提供依据。二、偏场作用下软电弹性圆柱的理论基础2.1非线性电弹性理论2.1.1连续介质运动和变形在连续介质力学中，物体的运动可视为从初始构形到当前构形的连续映射。假设物体由连续分布的质点组成，在初始时刻t_0，物体占据参考构形，其中质点的位置由物质坐标\boldsymbol{X}=(X_1,X_2,X_3)描述；在时刻t，物体运动到当前构形，质点的位置变为空间坐标\boldsymbol{x}=(x_1,x_2,x_3)，运动可表示为\boldsymbol{x}=\boldsymbol{\chi}(\boldsymbol{X},t)，其中\boldsymbol{\chi}为运动函数，它不仅描述了质点位置的变化，还蕴含了物体的变形信息。变形梯度张量\boldsymbol{F}是描述物体变形的关键量，定义为\boldsymbol{F}=\frac{\partial\boldsymbol{x}}{\partial\boldsymbol{X}}，其分量形式为F_{ij}=\frac{\partialx_i}{\partialX_j}（i,j=1,2,3）。变形梯度张量包含了物体局部旋转和变形的完整信息，其行列式\det(\boldsymbol{F})表示变形前后微元体积的变化比例，当\det(\boldsymbol{F})>1时，微元体积膨胀；当\det(\boldsymbol{F})<1时，微元体积收缩。例如，在简单拉伸情况下，假设物体沿x_1方向被均匀拉伸，拉伸比为\lambda，则变形梯度张量\boldsymbol{F}的分量为F_{11}=\lambda，F_{22}=F_{33}=1，F_{ij}=0（i\neqj），此时\det(\boldsymbol{F})=\lambda，直观地反映了物体在拉伸方向上的伸长和体积的相应变化。为了更细致地描述物体的变形，引入格林-拉格朗日应变张量\boldsymbol{E}，定义为\boldsymbol{E}=\frac{1}{2}(\boldsymbol{F}^T\boldsymbol{F}-\boldsymbol{I})，其中\boldsymbol{I}为单位张量。格林-拉格朗日应变张量能够准确描述有限变形情况下物体的应变状态，其分量E_{ij}=\frac{1}{2}(\sum_{k=1}^{3}F_{ki}F_{kj}-\delta_{ij})（\delta_{ij}为克罗内克符号），通过该张量可以计算出物体在各个方向上的正应变和切应变。在分析橡胶材料的大变形问题时，格林-拉格朗日应变张量能够准确描述橡胶在复杂受力情况下的应变分布，为研究橡胶的力学性能提供重要依据。2.1.2静电场方程静电场是由静止电荷产生的电场，其基本方程描述了电场的性质和电荷与电场之间的关系。高斯定律是静电场的重要基本定律之一，其积分形式为\oint_{S}\boldsymbol{D}\cdotd\boldsymbol{S}=q_{free}，其中\oint_{S}表示对封闭曲面S的面积分，\boldsymbol{D}为电位移矢量，d\boldsymbol{S}为封闭曲面上的面积元矢量，q_{free}为封闭曲面内的自由电荷总量。该定律表明，通过任意封闭曲面的电位移通量等于该封闭曲面内的自由电荷总量，它反映了静电场是有源场，电荷是电场的源。在分析平行板电容器时，根据高斯定律可以方便地计算出电容器极板间的电位移矢量和电场强度，假设平行板电容器极板面积为A，极板上的自由电荷面密度为\sigma，选取一个包含其中一个极板的高斯面，根据高斯定律可得D=\sigma，进而可求得电场强度E=\frac{\sigma}{\epsilon}（\epsilon为电介质的介电常数）。高斯定律的微分形式为\nabla\cdot\boldsymbol{D}=\rho_{free}，其中\nabla\cdot为散度算子，\rho_{free}为自由电荷体密度。它从微观角度描述了空间中某点的电位移矢量的散度与该点自由电荷体密度的关系。电场强度\boldsymbol{E}与电势\varphi之间存在密切关系，\boldsymbol{E}=-\nabla\varphi，该式表明电场强度等于电势的负梯度。这意味着电场强度的方向是电势降低最快的方向，其大小等于电势在该方向上的变化率。在求解静电场问题时，常常先通过求解电势的分布，再利用此关系得到电场强度的分布。例如，在点电荷产生的静电场中，已知点电荷q位于坐标原点，根据库仑定律和电势的定义可求得电势\varphi=\frac{q}{4\pi\epsilonr}（r为场点到点电荷的距离），再通过求梯度可得电场强度\boldsymbol{E}=\frac{q}{4\pi\epsilonr^2}\boldsymbol{e}_r（\boldsymbol{e}_r为径向单位矢量）。2.1.3运动微分方程与Maxwell应力张量对于软电弹性材料，在力电耦合作用下，其运动微分方程的推导基于动量守恒定律和能量守恒定律。根据动量守恒定律，单位体积内的动量变化率等于作用在该体积上的合力，即\rho\frac{\partial^2\boldsymbol{u}}{\partialt^2}=\nabla\cdot\boldsymbol{\sigma}+\rho\boldsymbol{f}，其中\rho为材料密度，\boldsymbol{u}为位移矢量，\boldsymbol{\sigma}为应力张量，\boldsymbol{f}为单位体积的体力。在考虑电-力耦合效应时，Maxwell应力张量\boldsymbol{T}^e发挥着关键作用。Maxwell应力张量描述了电场对电介质的作用力，其表达式为\boldsymbol{T}^e=\boldsymbol{D}\boldsymbol{E}-\frac{1}{2}(\boldsymbol{D}\cdot\boldsymbol{E})\boldsymbol{I}，其中\boldsymbol{I}为单位张量。将Maxwell应力张量考虑在内，软电弹性材料的运动微分方程变为\rho\frac{\partial^2\boldsymbol{u}}{\partialt^2}=\nabla\cdot(\boldsymbol{\sigma}+\boldsymbol{T}^e)+\rho\boldsymbol{f}，该方程全面考虑了机械应力和电场力对材料运动的影响。在分析软电弹性薄膜在电场作用下的振动问题时，运动微分方程能够准确描述薄膜在电场力和机械力共同作用下的运动状态，通过求解该方程可以得到薄膜的振动频率、振幅等重要参数，为薄膜器件的设计和优化提供理论支持。2.1.4能量密度函数和本构方程能量密度函数是描述材料在变形和电场作用下储存能量的重要物理量。对于软电弹性材料，其能量密度函数W通常是应变张量\boldsymbol{E}和电位移矢量\boldsymbol{D}的函数，即W=W(\boldsymbol{E},\boldsymbol{D})。能量密度函数包含弹性应变能和电场能两部分，弹性应变能反映了材料因变形而储存的能量，电场能则与材料在电场中的极化状态相关。在分析软电弹性材料的拉伸过程时，能量密度函数能够定量描述材料在拉伸变形和电场作用下能量的储存和转换情况，为研究材料的力学性能和电-力耦合特性提供重要依据。本构方程是描述材料应力-应变关系和电位移-电场关系的数学表达式，它反映了材料的固有属性。对于线性软电弹性材料，常见的本构方程为胡克定律的推广形式，如应力-应变关系可表示为\boldsymbol{\sigma}=C:\boldsymbol{E}-e^T:\boldsymbol{E}，其中C为弹性常数张量，e为压电常数张量；电位移-电场关系为\boldsymbol{D}=e:\boldsymbol{E}+\epsilon:\boldsymbol{E}，其中\epsilon为介电常数张量。这些本构方程在线性范围内能够较好地描述材料的力电耦合行为。然而，对于非线性软电弹性材料，本构方程更为复杂，需要考虑材料的非线性特性，如大变形、非线性电致伸缩效应等。常见的非线性本构模型包括Mooney-Rivlin模型、Neo-Hookean模型等。以Mooney-Rivlin模型为例，其能量密度函数表示为W=C_{10}(I_1-3)+C_{01}(I_2-3)+D_1(J-1)^2，其中C_{10}、C_{01}、D_1为材料常数，I_1、I_2为应变不变量，J=\det(\boldsymbol{F})为体积比。基于该能量密度函数，可以推导出相应的应力-应变关系和电位移-电场关系，从而更准确地描述非线性软电弹性材料在复杂载荷和电场作用下的力学行为和电-力耦合特性。在研究橡胶类软电弹性材料的大变形问题时，Mooney-Rivlin模型能够更准确地描述材料的应力-应变关系，为橡胶材料的工程应用提供更可靠的理论指导。2.2电弹性体有限偏场理论2.2.1三种构型表述在研究软电弹性材料的力学行为时，为准确描述材料的变形和运动，常采用三种构型表述：拉格朗日描述、更新拉格朗日描述和当前构型描述，它们各自具有独特的特点和应用场景。拉格朗日描述（Lagrangiandescription），也被称为物质描述或随体描述，将物质坐标（也称为拉格朗日坐标）\boldsymbol{X}和时间t作为独立变量。在这种描述中，观察者关注的是每个物质点的运动轨迹和变形历史，就如同给每个物质点贴上了一个“标签”，跟踪其在整个过程中的变化。以软电弹性圆柱的拉伸过程为例，在拉格朗日描述下，我们会记录圆柱中每个微元体（物质点）从初始位置\boldsymbol{X}开始，随着时间t的推移，在电场和外力作用下，其位置、形状和应力状态等物理量的变化。这种描述方式的优点是能够清晰地反映每个物质点的历史信息，便于分析材料的变形过程和力学响应的因果关系。在分析软电弹性材料在复杂加载路径下的力学性能时，拉格朗日描述可以准确地描述材料内部各点的变形历程，为建立本构模型提供详细的数据支持。其缺点是当物体发生大变形时，参考构型与当前构型之间的差异较大，可能导致数学处理上的复杂性增加。更新拉格朗日描述（UpdatedLagrangiandescription）是在拉格朗日描述的基础上发展而来的。它同样以物质坐标\boldsymbol{X}和时间t为独立变量，但与拉格朗日描述不同的是，它在每个增量步中都将前一时刻的构型作为新的参考构型。在软电弹性圆柱的动态加载过程中，随着时间的推进，圆柱不断发生变形，更新拉格朗日描述会在每一个微小的时间步长内，将上一个时间步结束时圆柱的构型作为新的参考构型，来描述当前时间步内的变形和力学响应。这种描述方式的优势在于，在处理大变形问题时，由于不断更新参考构型，使得当前构型与参考构型之间的差异相对较小，从而在一定程度上简化了数学计算。在分析橡胶类软电弹性材料在大变形下的力学行为时，更新拉格朗日描述能够更有效地处理材料的非线性变形，提高计算效率和准确性。然而，它也存在一定的局限性，由于每次都要更新参考构型，需要存储更多的中间构型信息，增加了计算成本和数据存储量。当前构型描述（Currentconfigurationdescription），又称为欧拉描述或空间描述，将空间坐标\boldsymbol{x}和时间t作为独立变量。在这种描述中，观察者固定在空间中的某一点，关注的是在不同时刻，通过该空间点的物质点的物理量变化。对于软电弹性圆柱，在当前构型描述下，我们会在空间中选取一些固定的位置点，观察随着时间的变化，这些位置点上软电弹性材料的应力、应变、电场强度等物理量的数值。这种描述方式的优点是在描述流体和大变形问题时，数学形式相对简洁，便于进行数值计算。在研究软电弹性材料在高速变形或流动状态下的力学行为时，当前构型描述能够更好地处理复杂的边界条件和变形情况，提高计算效率。但是，它难以直接反映物质点的历史信息，对于需要分析材料变形历程的问题，使用当前构型描述可能会存在一定的困难。2.2.2增量线性方程推导从非线性电弹性理论出发，推导有限偏场下的增量线性方程是研究软电弹性材料在小扰动下力学行为的关键步骤。假设软电弹性材料在初始状态下受到一个有限的偏场作用，包括机械偏场（如预加载荷）和电场偏场（如预加电场）。在这个偏场作用下，材料处于一个平衡状态，其位移场为\boldsymbol{u}_0，电场强度为\boldsymbol{E}_0，电位移为\boldsymbol{D}_0。当材料受到一个小的扰动时，位移场变为\boldsymbol{u}=\boldsymbol{u}_0+\delta\boldsymbol{u}，电场强度变为\boldsymbol{E}=\boldsymbol{E}_0+\delta\boldsymbol{E}，电位移变为\boldsymbol{D}=\boldsymbol{D}_0+\delta\boldsymbol{D}，其中\delta\boldsymbol{u}、\delta\boldsymbol{E}和\delta\boldsymbol{D}分别为位移、电场强度和电位移的增量。根据非线性电弹性理论，运动微分方程为\rho\frac{\partial^2\boldsymbol{u}}{\partialt^2}=\nabla\cdot(\boldsymbol{\sigma}+\boldsymbol{T}^e)+\rho\boldsymbol{f}，将位移、电场强度和电位移的增量形式代入该方程，并忽略高阶小量。在考虑电-力耦合效应时，Maxwell应力张量\boldsymbol{T}^e=\boldsymbol{D}\boldsymbol{E}-\frac{1}{2}(\boldsymbol{D}\cdot\boldsymbol{E})\boldsymbol{I}，将其展开并代入运动微分方程。同时，本构方程描述了应力-应变关系和电位移-电场关系，对于非线性软电弹性材料，本构方程较为复杂，但在小扰动情况下，可以对其进行线性化处理。以应力-应变关系为例，假设本构方程为\boldsymbol{\sigma}=\boldsymbol{\sigma}(\boldsymbol{E},\boldsymbol{D})，在小扰动下，对其进行泰勒展开并保留一阶项，得到增量形式的应力-应变关系。通过一系列的数学推导和化简，最终可以得到有限偏场下的增量线性运动微分方程：\rho\frac{\partial^2\delta\boldsymbol{u}}{\partialt^2}=\nabla\cdot(\delta\boldsymbol{\sigma}+\delta\boldsymbol{T}^e)+\rho\delta\boldsymbol{f}，其中\delta\boldsymbol{\sigma}和\delta\boldsymbol{T}^e分别为应力增量和Maxwell应力张量增量，\delta\boldsymbol{f}为体力增量。在推导过程中，对各物理量进行增量分解，将复杂的非线性方程转化为线性方程，这是基于小扰动假设，认为增量部分相对较小，高阶小量对结果的影响可以忽略不计。这种线性化处理使得方程的求解变得相对简单，能够利用线性代数和偏微分方程的相关理论和方法进行求解。2.2.3增量外场方程和边界条件在得到有限偏场下的增量线性方程后，需要确定相应的增量外场方程和边界条件，以完整地描述软电弹性材料的力学行为。增量外场方程包括电学和力学两个方面。在电学方面，根据高斯定律\nabla\cdot\boldsymbol{D}=\rho_{free}，在增量情况下，有\nabla\cdot\delta\boldsymbol{D}=\delta\rho_{free}，其中\delta\rho_{free}为自由电荷体密度的增量。这表明电位移增量的散度等于自由电荷体密度增量，反映了电场在小扰动下的变化与自由电荷分布变化之间的关系。在力学方面，根据动量守恒定律，增量形式的运动方程已在前面推导得到。同时，还需考虑应力张量增量的散度与外力增量之间的关系。对于软电弹性圆柱，假设在圆柱表面受到分布力\boldsymbol{t}的作用，在增量情况下，有\boldsymbol{t}=\boldsymbol{t}_0+\delta\boldsymbol{t}，其中\boldsymbol{t}_0为偏场作用下的表面力，\delta\boldsymbol{t}为表面力的增量。根据边界条件，在圆柱表面，应力张量增量与表面力增量之间满足\boldsymbol{n}\cdot\delta\boldsymbol{\sigma}=\delta\boldsymbol{t}，其中\boldsymbol{n}为圆柱表面的单位外法线向量。这一条件表明在圆柱表面，应力张量增量在法线方向上的分量等于表面力增量，保证了力学边界的连续性。电学边界条件也至关重要。在软电弹性圆柱与其他介质的分界面上，电位移矢量和电场强度矢量需要满足一定的连续性条件。假设分界面两侧的电位移矢量分别为\boldsymbol{D}_1和\boldsymbol{D}_2，电场强度矢量分别为\boldsymbol{E}_1和\boldsymbol{E}_2，在增量情况下，有\boldsymbol{n}\cdot(\delta\boldsymbol{D}_1-\delta\boldsymbol{D}_2)=\delta\sigma_{free}和\boldsymbol{n}\times(\delta\boldsymbol{E}_1-\delta\boldsymbol{E}_2)=0，其中\delta\sigma_{free}为分界面上自由电荷面密度的增量。第一个条件表示电位移矢量增量在分界面法线方向上的分量差等于分界面上自由电荷面密度增量，第二个条件表示电场强度矢量增量在分界面切向方向上的分量连续。这些电学边界条件确保了电场在分界面处的连续性和物理合理性。三、非均匀偏场作用下软电弹性管中环向导波3.1各向同性软电弹性体偏场理论3.1.1非线性理论各向同性软电弹性体的非线性电弹性理论是研究其在复杂力电耦合作用下力学行为的基础。在大变形情况下，材料的应力-应变关系和电位移-电场关系呈现出高度的非线性特性。从运动学角度看，物体的变形通过变形梯度张量\boldsymbol{F}来描述，其定义为\boldsymbol{F}=\frac{\partial\boldsymbol{x}}{\partial\boldsymbol{X}}，其中\boldsymbol{x}是当前构型下的位置矢量，\boldsymbol{X}是参考构型下的位置矢量。格林-拉格朗日应变张量\boldsymbol{E}=\frac{1}{2}(\boldsymbol{F}^T\boldsymbol{F}-\boldsymbol{I})则更细致地刻画了物体的应变状态，它能够准确描述有限变形情况下物体的拉伸、剪切等变形情况。以橡胶类软电弹性材料在复杂拉伸和剪切复合作用下的变形分析为例，格林-拉格朗日应变张量可以精确地描述材料在各个方向上的应变变化，为后续的力学分析提供关键数据支持。在静电场方面，高斯定律\nabla\cdot\boldsymbol{D}=\rho_{free}描述了电场与自由电荷之间的关系，其中\boldsymbol{D}是电位移矢量，\rho_{free}是自由电荷体密度。电场强度\boldsymbol{E}与电位移\boldsymbol{D}之间通过介电常数张量\boldsymbol{\epsilon}相关联，即\boldsymbol{D}=\boldsymbol{\epsilon}\cdot\boldsymbol{E}，对于各向同性材料，介电常数\epsilon是一个标量。在分析软电弹性材料在非均匀电场中的极化行为时，高斯定律和电位移-电场关系能够帮助我们理解电场分布和材料的极化机制，为研究材料的电性能提供理论依据。本构关系是描述材料固有属性的关键，对于各向同性软电弹性体，常见的本构模型有Mooney-Rivlin模型和Neo-Hookean模型。Mooney-Rivlin模型的应变能函数表示为W=C_{10}(I_1-3)+C_{01}(I_2-3)+D_1(J-1)^2，其中C_{10}、C_{01}、D_1为材料常数，I_1、I_2为应变不变量，J=\det(\boldsymbol{F})为体积比。基于该应变能函数，可以通过对其求偏导数得到应力张量\boldsymbol{\sigma}和电位移矢量\boldsymbol{D}的表达式，从而描述材料在力电耦合作用下的力学行为。在研究橡胶类软电弹性材料在大变形下的力电耦合特性时，Mooney-Rivlin模型能够更准确地描述材料的应力-应变关系和电位移-电场关系，为材料的工程应用提供可靠的理论指导。Neo-Hookean模型是一种简化的超弹性模型，其应变能函数为W=\frac{\mu}{2}(I_1-3)+\frac{\lambda}{4}(J-1)^2，其中\mu和\lambda是拉梅常数。虽然该模型相对简单，但在一些情况下，如小变形或对材料性能要求不是特别精确的分析中，能够快速有效地描述材料的力学行为。3.1.2有限偏场下的线性增量理论在有限偏场作用下，为了研究软电弹性体在小扰动下的力学行为，需要将非线性电弹性理论进行线性化处理，得到线性增量理论。假设软电弹性体在初始状态下受到有限的偏场作用，包括机械偏场（如预加载荷）和电场偏场（如预加电场）。在这个偏场作用下，材料处于一个平衡状态，其位移场为\boldsymbol{u}_0，电场强度为\boldsymbol{E}_0，电位移为\boldsymbol{D}_0。当材料受到一个小的扰动时，位移场变为\boldsymbol{u}=\boldsymbol{u}_0+\delta\boldsymbol{u}，电场强度变为\boldsymbol{E}=\boldsymbol{E}_0+\delta\boldsymbol{E}，电位移变为\boldsymbol{D}=\boldsymbol{D}_0+\delta\boldsymbol{D}，其中\delta\boldsymbol{u}、\delta\boldsymbol{E}和\delta\boldsymbol{D}分别为位移、电场强度和电位移的增量。基于非线性电弹性理论的运动微分方程\rho\frac{\partial^2\boldsymbol{u}}{\partialt^2}=\nabla\cdot(\boldsymbol{\sigma}+\boldsymbol{T}^e)+\rho\boldsymbol{f}，将位移、电场强度和电位移的增量形式代入该方程，并忽略高阶小量。在考虑电-力耦合效应时，Maxwell应力张量\boldsymbol{T}^e=\boldsymbol{D}\boldsymbol{E}-\frac{1}{2}(\boldsymbol{D}\cdot\boldsymbol{E})\boldsymbol{I}，将其展开并代入运动微分方程。同时，对本构方程进行线性化处理。以Mooney-Rivlin模型为例，对其应变能函数求偏导数得到应力和电位移的表达式后，在小扰动情况下，对这些表达式进行泰勒展开并保留一阶项，得到增量形式的应力-应变关系和电位移-电场关系。通过一系列的数学推导和化简，最终可以得到有限偏场下的增量线性运动微分方程。增量线性方程的推导基于小扰动假设，认为增量部分相对较小，高阶小量对结果的影响可以忽略不计。这种线性化处理使得方程的求解变得相对简单，能够利用线性代数和偏微分方程的相关理论和方法进行求解。在分析软电弹性管在预加电场和小机械扰动作用下的振动问题时，有限偏场下的线性增量理论可以准确地描述软电弹性管的振动特性，为管道的动力学分析提供有效的理论工具。3.2软电弹性管的轴对称变形3.2.1变形假设与模型建立在实际应用中，软电弹性管常受到各种复杂的载荷作用，其中轴力和径向电压是较为常见的载荷形式。为了深入研究软电弹性管在这些载荷作用下的力学行为，我们提出合理的轴对称变形假设。假设软电弹性管在初始状态下为各向同性的厚壁圆筒，其内径为a，外径为b。在轴力N和径向电压V的联合作用下，管发生轴对称变形。在变形过程中，管的横截面始终保持为圆形，且管的轴线方向不变。基于上述假设，我们建立相应的数学模型。采用柱坐标系(r,\theta,z)，其中r为径向坐标，\theta为周向坐标，z为轴向坐标。设管在变形后的径向位移为u(r)，轴向位移为w(r)，周向位移为v(r,\theta,z)=0（由于轴对称性，周向位移为零）。根据变形几何关系，可得到格林-拉格朗日应变张量\boldsymbol{E}的非零分量表达式。径向应变E_{rr}=\frac{\partialu}{\partialr}+\frac{1}{2}(\frac{\partialu}{\partialr})^2，周向应变E_{\theta\theta}=\frac{u}{r}+\frac{1}{2}(\frac{u}{r})^2，轴向应变E_{zz}=\frac{\partialw}{\partialz}+\frac{1}{2}(\frac{\partialw}{\partialz})^2，以及径向-轴向剪切应变E_{rz}=\frac{1}{2}(\frac{\partialu}{\partialz}+\frac{\partialw}{\partialr}+\frac{\partialu}{\partialr}\frac{\partialw}{\partialz})。这些应变分量的准确描述对于后续分析软电弹性管的力学响应至关重要。3.2.2非线性控制方程推导根据力电耦合原理，推导轴力和径向电压联合作用下的非线性控制方程。从能量角度出发，软电弹性管的总能量包括弹性应变能W^e、电场能W^e和外力做功W^e。弹性应变能密度W^e可根据所选的本构模型来确定，如采用Mooney-Rivlin模型，其表达式为W^e=C_{10}(I_1-3)+C_{01}(I_2-3)+D_1(J-1)^2，其中C_{10}、C_{01}、D_1为材料常数，I_1、I_2为应变不变量，J=\det(\boldsymbol{F})为体积比。应变不变量I_1=E_{rr}+E_{\theta\theta}+E_{zz}，I_2=E_{rr}E_{\theta\theta}+E_{rr}E_{zz}+E_{\theta\theta}E_{zz}-E_{rz}^2。电场能密度W^e与电位移矢量\boldsymbol{D}和电场强度\boldsymbol{E}相关，W^e=\frac{1}{2}\boldsymbol{D}\cdot\boldsymbol{E}。在柱坐标系下，电场强度E=-\frac{\partial\varphi}{\partialr}（假设电势\varphi仅为径向坐标r的函数），电位移矢量D=\epsilonE（\epsilon为介电常数）。外力做功包括轴力N和径向电压V所做的功。轴力N在轴向位移上做功，其功密度为N\frac{\partialw}{\partialz}；径向电压V在电位移上做功，其功密度为VD。根据能量变分原理\delta(W^e+W^e+W^e)=0，对总能量关于位移u和w以及电势\varphi求变分。在求变分过程中，利用格林公式等数学工具对各项进行处理。对于弹性应变能项，根据应变与位移的关系，对W^e关于位移求偏导数，再乘以位移的变分并在管的体积内积分；对于电场能项，根据电场强度与电势的关系，对W^e关于电势求偏导数，再乘以电势的变分并在管的体积内积分；对于外力做功项，分别对轴力和径向电压做功项关于相应的位移和电位移求偏导数，再乘以其变分并在管的体积内积分。通过一系列复杂的数学推导和化简，最终得到轴力和径向电压联合作用下的非线性控制方程。这些非线性控制方程全面描述了软电弹性管在力电耦合作用下的力学行为，为后续分析管的变形、应力分布以及导波传播特性等提供了理论基础。3.3柱坐标下增量场的状态空间列式3.3.1状态变量选取在柱坐标系下研究软电弹性管中导波的传播特性时，合理选取状态变量是建立状态空间方程的关键。考虑到软电弹性材料的力电耦合特性，选取位移分量u_r（径向位移）、u_{\theta}（周向位移）、u_z（轴向位移），应力分量\sigma_{rr}（径向应力）、\sigma_{\theta\theta}（周向应力）、\sigma_{zz}（轴向应力）、\sigma_{r\theta}（径向-周向剪切应力）、\sigma_{rz}（径向-轴向剪切应力）、\sigma_{\thetaz}（周向-轴向剪切应力），以及电势\varphi和电位移分量D_r（径向电位移）、D_{\theta}（周向电位移）、D_z（轴向电位移）作为状态变量。这些状态变量能够全面描述软电弹性管在导波传播过程中的力学和电学状态。位移分量u_r、u_{\theta}、u_z反映了软电弹性管在不同方向上的变形情况。在导波传播过程中，管的不同部位会发生不同程度的径向、周向和轴向位移，这些位移的变化与导波的传播特性密切相关。例如，在导波的传播路径上，当遇到缺陷或不均匀性时，位移分量会发生明显变化，通过监测这些位移分量的变化可以检测管中的缺陷。应力分量\sigma_{rr}、\sigma_{\theta\theta}、\sigma_{zz}、\sigma_{r\theta}、\sigma_{rz}、\sigma_{\thetaz}描述了软电弹性管内部的应力分布。导波的传播会引起管内应力的波动，不同的应力分量在导波传播过程中起着不同的作用。\sigma_{rr}和\sigma_{\theta\theta}与管的径向和周向变形相关，它们的变化会影响导波的传播速度和衰减特性；\sigma_{rz}和\sigma_{\thetaz}则与管的剪切变形相关，对导波的模态转换和能量传递有重要影响。电势\varphi和电位移分量D_r、D_{\theta}、D_z体现了软电弹性管的电学特性。在力电耦合作用下，电场的变化会影响材料的力学性能，反之亦然。导波传播过程中，电势和电位移的变化与应力和位移的变化相互关联。当导波传播到软电弹性管的不同部位时，由于材料的电学性质和几何形状的变化，电势和电位移会发生相应的改变，这些变化又会反过来影响导波的传播。3.3.2状态空间方程推导基于软电弹性管的轴对称变形假设和增量线性理论，对增量场进行分析，推导柱坐标系下的状态空间列式。从运动微分方程出发，在柱坐标系下，考虑到轴对称变形，忽略体力和惯性力的高阶项，运动微分方程在径向、周向和轴向的分量形式分别为：\begin{align*}\frac{\partial\sigma_{rr}}{\partialr}+\frac{1}{r}\frac{\partial\sigma_{r\theta}}{\partial\theta}+\frac{\partial\sigma_{rz}}{\partialz}+\frac{\sigma_{rr}-\sigma_{\theta\theta}}{r}&=\rho\frac{\partial^2u_r}{\partialt^2}\\\frac{\partial\sigma_{r\theta}}{\partialr}+\frac{1}{r}\frac{\partial\sigma_{\theta\theta}}{\partial\theta}+\frac{\partial\sigma_{\thetaz}}{\partialz}+\frac{2\sigma_{r\theta}}{r}&=\rho\frac{\partial^2u_{\theta}}{\partialt^2}\\\frac{\partial\sigma_{rz}}{\partialr}+\frac{1}{r}\frac{\partial\sigma_{\thetaz}}{\partial\theta}+\frac{\partial\sigma_{zz}}{\partialz}+\frac{\sigma_{rz}}{r}&=\rho\frac{\partial^2u_z}{\partialt^2}\end{align*}在电学方面，根据高斯定律\nabla\cdot\boldsymbol{D}=\rho_{free}，在柱坐标系下，对于轴对称情况，有\frac{\partialD_r}{\partialr}+\frac{D_r}{r}+\frac{\partialD_z}{\partialz}=0（假设自由电荷体密度\rho_{free}=0）。本构方程描述了应力-应变关系和电位移-电场关系。对于软电弹性材料，在增量情况下，应力-应变关系可表示为：\begin{align*}\sigma_{ij}&=C_{ijkl}\epsilon_{kl}+e_{kij}E_k\\D_i&=e_{ijk}\epsilon_{jk}+\epsilon_{ik}E_k\end{align*}其中C_{ijkl}为弹性常数张量，e_{kij}为压电常数张量，\epsilon_{ij}为应变张量，\epsilon_{ik}为介电常数张量，E_k为电场强度分量。应变张量\epsilon_{ij}与位移分量的关系为：\begin{align*}\epsilon_{rr}&=\frac{\partialu_r}{\partialr}\\\epsilon_{\theta\theta}&=\frac{u_r}{r}+\frac{1}{r}\frac{\partialu_{\theta}}{\partial\theta}\\\epsilon_{zz}&=\frac{\partialu_z}{\partialz}\\\epsilon_{r\theta}&=\frac{1}{2}(\frac{\partialu_{\theta}}{\partialr}+\frac{1}{r}\frac{\partialu_r}{\partial\theta}-\frac{u_{\theta}}{r})\\\epsilon_{rz}&=\frac{1}{2}(\frac{\partialu_z}{\partialr}+\frac{\partialu_r}{\partialz})\\\epsilon_{\thetaz}&=\frac{1}{2}(\frac{\partialu_z}{\partial\theta}+\frac{1}{r}\frac{\partialu_{\theta}}{\partialz})\end{align*}电场强度分量E_k与电势\varphi的关系为E_k=-\frac{\partial\varphi}{\partialx_k}（x_k为坐标分量）。将上述运动微分方程、本构方程、应变-位移关系和电场-电势关系进行整理和推导。首先，将应变和电场强度用位移和电势表示，代入本构方程得到应力和电位移与位移和电势的关系。然后，将应力和电位移的表达式代入运动微分方程和高斯定律方程。通过一系列的数学运算和化简，最终得到状态空间方程的一般形式：\frac{\partial\boldsymbol{U}}{\partialr}=\boldsymbol{M}\boldsymbol{U}其中\boldsymbol{U}为状态变量向量，\boldsymbol{U}=[u_r,u_{\theta},u_z,\sigma_{rr},\sigma_{\theta\theta},\sigma_{zz},\sigma_{r\theta},\sigma_{rz},\sigma_{\thetaz},\varphi,D_r,D_{\theta},D_z]^T，\boldsymbol{M}为系统矩阵，其元素由材料参数（如弹性常数、压电常数、介电常数）和几何参数（如管的半径）组成。系统矩阵\boldsymbol{M}反映了软电弹性管的物理特性和几何结构对导波传播的影响，通过对状态空间方程的求解，可以得到导波在软电弹性管中的传播特性，如频散关系、模态分布等。3.4软电弹性管中环向波的频散关系3.4.1数值求解方法为了深入研究软电弹性管中环向波的传播特性，需要准确求解其频散关系。这里采用将状态空间列式与近似分层技术相结合的数值求解方法。在状态空间列式方面，基于前面推导得到的柱坐标下的状态空间方程\frac{\partial\boldsymbol{U}}{\partialr}=\boldsymbol{M}\boldsymbol{U}，其中\boldsymbol{U}为状态变量向量，包含位移、应力、电势和电位移等状态变量；\boldsymbol{M}为系统矩阵，其元素由材料参数（如弹性常数、压电常数、介电常数）和几何参数（如管的半径）组成。该状态空间方程描述了软电弹性管中导波传播时状态变量沿径向的变化规律。近似分层技术则是将软电弹性管沿径向进行分层处理。把软电弹性管从内径a到外径b划分为N个薄层，每个薄层的厚度为\Deltar=\frac{b-a}{N}。在每个薄层内，假设材料参数和几何参数是均匀的，这样可以将连续的软电弹性管离散化为一系列的均匀层。在每个薄层中，利用状态空间方程进行数值求解，通过层与层之间的边界条件进行连接，逐步得到整个软电弹性管中导波的传播特性。具体求解过程中，采用四阶龙格-库塔法对状态空间方程进行数值积分。对于第i层，从内边界r_i到外边界r_{i+1}，根据四阶龙格-库塔法，状态变量向量\boldsymbol{U}的更新公式为：\begin{align*}\boldsymbol{k}_1&=\Deltar\boldsymbol{M}(r_i)\boldsymbol{U}(r_i)\\\boldsymbol{k}_2&=\Deltar\boldsymbol{M}(r_i+\frac{\Deltar}{2})\left(\boldsymbol{U}(r_i)+\frac{\boldsymbol{k}_1}{2}\right)\\\boldsymbol{k}_3&=\Deltar\boldsymbol{M}(r_i+\frac{\Deltar}{2})\left(\boldsymbol{U}(r_i)+\frac{\boldsymbol{k}_2}{2}\right)\\\boldsymbol{k}_4&=\Deltar\boldsymbol{M}(r_{i+1})\left(\boldsymbol{U}(r_i)+\boldsymbol{k}_3\right)\\\boldsymbol{U}(r_{i+1})&=\boldsymbol{U}(r_i)+\frac{1}{6}(\boldsymbol{k}_1+2\boldsymbol{k}_2+2\boldsymbol{k}_3+\boldsymbol{k}_4)\end{align*}通过上述方法，从软电弹性管的内径开始，逐步向外径进行数值积分，最终得到外径处的状态变量向量。再结合边界条件，如管内外表面的应力自由条件、电位移连续条件等，确定导波的传播常数，进而得到频散关系。3.4.2频散曲线分析通过上述数值求解方法，对软电弹性管中环向波的频散关系进行数值计算，得到频散曲线。频散曲线直观地展示了导波的频率与波数之间的关系，通过分析频散曲线，可以深入了解软电弹性管中环向波的传播特性，以及不同参数对频散特性的影响。材料参数对频散曲线的影响：改变软电弹性材料的弹性常数、压电常数和介电常数，观察频散曲线的变化。当弹性常数增大时，频散曲线整体向高频方向移动，这表明导波的传播速度增加，因为弹性常数反映了材料抵抗变形的能力，弹性常数越大，材料越“硬”，波在其中传播速度越快。压电常数的变化会影响力电耦合效应，进而改变频散曲线的形状。当压电常数增大时，某些模态的频散曲线斜率发生变化，这意味着这些模态的导波在传播过程中，其频率与波数的关系发生改变，反映了力电耦合对导波传播的影响。介电常数的变化主要影响电场相关的特性，介电常数增大时，电场能在总能量中的占比发生变化，导致频散曲线的形态改变，尤其是与电场密切相关的模态，其频散特性变化更为明显。几何参数对频散曲线的影响：研究软电弹性管的内径a、外径b和壁厚h=b-a对频散曲线的影响。当内径a增大而外径b不变时，即壁厚减小，频散曲线的某些模态出现向低频方向移动的趋势。这是因为壁厚减小，管的刚度相对降低，导波传播时受到的约束减弱，导致传播速度降低，频率随之下降。当外径b增大而内径a不变时，壁厚增加，频散曲线整体向高频方向移动，说明管的刚度增加，导波传播速度加快。此外，管的长径比（管的长度与外径之比）也会对频散曲线产生影响。长径比增大时，导波在传播过程中与管边界的相互作用次数相对减少，某些模态的频散特性变得更加接近无限长管的情况，频散曲线的变化趋势相对平缓。偏场对频散曲线的影响：分析轴力和径向电压等偏场对频散曲线的影响。当施加轴向拉力时，软电弹性管在轴向被拉伸，材料的应力状态发生改变，导致频散曲线发生明显变化。轴向拉力使得管的轴向刚度增加，与轴向相关的模态的频散曲线向高频方向移动，同时模态的分布也发生改变。径向电压的施加会引入电场，产生电-力耦合效应。随着径向电压的增大，频散曲线中与电-力耦合相关的模态的频率和波数关系发生显著变化，某些模态可能出现新的分支或交叉现象，这反映了电场对导波传播特性的复杂影响。在实际应用中，通过调整偏场的大小和方向，可以有效地调控软电弹性管中环向波的传播特性，为基于导波的检测和传感技术提供更多的调控手段。3.5数值算例与讨论3.5.1Neo-Hookean理想介电模型为深入探究软电弹性管中环向导波的传播特性，采用Neo-Hookean理想介电模型进行数值模拟。在该模型中，软电弹性材料被视为各向同性的超弹性体，其应变能函数为W=\frac{\mu}{2}(I_1-3)+\frac{\lambda}{4}(J-1)^2，其中\mu和\lambda是拉梅常数，I_1为第一应变不变量，I_1=tr(\boldsymbol{C})（\boldsymbol{C}=\boldsymbol{F}^T\boldsymbol{F}为右柯西-格林变形张量），J=\det(\boldsymbol{F})为体积比。该模型适用于描述橡胶类软电弹性材料在小变形到中等变形范围内的力学行为，在分析软电弹性管的变形和导波传播时具有重要应用价值。假设软电弹性管的材料参数为：初始剪切模量\mu_0=1\times10^6\Pa，初始体积模量K_0=3\times10^6\Pa，初始介电常数\epsilon_0=10\epsilon_0（\epsilon_0为真空介电常数，\epsilon_0=8.854\times10^{-12}\F/m）。几何参数设定为：内径a=0.01\m，外径b=0.02\m。偏场条件为：轴向拉力N=100\N，径向电压V=1000\V。在数值计算过程中，将软电弹性管沿径向划分为100个薄层，采用四阶龙格-库塔法对状态空间方程进行数值积分求解。3.5.2有效性分析为验证所采用数值求解方法的准确性和收敛性，将数值结果与已有研究结果进行对比。在相同的材料参数、几何参数和偏场条件下，与文献中采用其他方法得到的软电弹性管中环向波频散关系进行比较。结果显示，本文计算得到的频散曲线与文献结果在趋势上高度吻合，在低频段和高频段的频率和波数数值差异均在可接受范围内，验证了数值求解方法在计算频散关系方面的准确性。为进一步分析收敛性，进行网格收敛性测试。逐渐增加软电弹性管沿径向划分的薄层数量，从50层逐步增加到200层，观察频散曲线的变化。随着薄层数量的增加，频散曲线逐渐趋于稳定，当薄层数量达到100层后，继续增加薄层数量对频散曲线的影响极小，表明此时数值计算结果已收敛，所采用的数值求解方法具有良好的收敛性，能够准确地模拟软电弹性管中环向波的传播特性。3.5.3增量SH波的传播在软电弹性管中，增量SH波是一种重要的导波模式，其传播特性对于理解软电弹性管的力学行为具有重要意义。研究发现，增量SH波的波速随频率的变化呈现出复杂的关系。在低频段，波速随着频率的增加而逐渐增大，这是因为在低频时，软电弹性管的刚度对波速的影响较大，随着频率升高，材料的弹性响应增强，导致波速增大。在高频段，波速逐渐趋于稳定，这是由于高频下材料的惯性效应逐渐凸显，限制了波速的进一步增加。增量SH波在传播过程中存在衰减现象。衰减系数与频率、材料参数以及偏场条件密切相关。随着频率的增加，衰减系数逐渐增大，这是因为高频波在传播过程中更容易与材料内部的微观结构相互作用，导致能量损耗增加。材料的粘性和电导率对衰减系数也有显著影响，粘性越大，能量在粘性耗散中损失越多，衰减系数增大；电导率增加，会引起电磁损耗增加，同样导致衰减系数增大。偏场条件如轴向拉力和径向电压也会影响衰减特性，轴向拉力的增加会使软电弹性管的内部应力状态改变，影响波的传播和能量损耗，径向电压的变化会改变材料的电-力耦合特性，进而影响衰减系数。3.5.4增量Lamb波的传播增量Lamb波在软电弹性管中的传播特性与增量SH波存在明显差异。增量Lamb波具有多种模态，包括对称模态（S模态）和反对称模态（A模态），不同模态的传播特性各不相同。在频散特性方面，增量Lamb波的频散曲线比增量SH波更为复杂。对称模态和反对称模态的频散曲线具有不同的形状和变化趋势。对于对称模态，在低频段，波速随着频率的增加而快速增大，在高频段逐渐趋于稳定。反对称模态在低频段波速较低，随着频率增加，波速先减小后增大，呈现出明显的转折点。这种频散特性的差异源于不同模态下软电弹性管的振动模式和应力分布不同。对称模态下，管的振动主要以纵向伸缩为主，而反对称模态下，管的振动主要表现为弯曲变形。在波速方面，增量Lamb波的波速在不同模态和频率下与增量SH波有明显区别。在低频段，增量Lamb波的对称模态波速高于增量SH波，反对称模态波速低于增量SH波。随着频率的增加，不同模态的波速差异逐渐减小，但始终保持各自的变化规律。在衰减特性方面，增量Lamb波的衰减也与模态和频率相关。对称模态在低频段衰减较小，随着频率增加，衰减逐渐增大。反对称模态在低频段衰减相对较大，且在某些频率范围内衰减变化较为剧烈。与增量SH波相比，增量Lamb波的衰减特性受模态影响更为显著，不同模态的衰减机制也有所不同。对称模态的衰减主要与材料的内部摩擦和弹性损耗有关，而反对称模态的衰减还受到弯曲振动引起的附加能量损耗的影响。四、功能梯度软电弹性管中轴对称波4.1FG-SET的轴对称变形4.1.1材料特性描述功能梯度软电弹性管（FG-SET）的显著特征在于其材料特性沿径向呈现连续变化。为了准确表征这一特性，假设FG-SET由Mooney-Rivlin理想介电模型描述，该模型在描述软电弹性材料的力电耦合行为方面具有良好的准确性和适用性。材料参数随径向坐标的仿射变化是表征FG-SET特性的关键方式。具体而言，对于弹性常数和介电常数等材料参数，可表示为径向坐标r的线性函数。以弹性常数C为例，设其在管内径r=a处的值为C_a，在管外径r=b处的值为C_b，则弹性常数C随径向坐标r的变化可表示为：C(r)=C_a+\frac{C_b-C_a}{b-a}(r-a)类似地，介电常数\epsilon随径向坐标r的变化可表示为：\epsilon(r)=\epsilon_a+\frac{\epsilon_b-\epsilon_a}{b-a}(r-a)其中\epsilon_a和\epsilon_b分别为管内径和外径处的介电常数。这种仿射变化的假设使得材料参数在管的径向方向上呈现连续且平滑的变化，更符合实际工程中功能梯度材料的特性。通过这种方式，可以准确地描述FG-SET在不同径向位置的材料性能差异，为后续分析其在力电非均匀偏场作用下的力学行为和导波传播特性奠定基础。4.1.2变形分析与方程建立考虑FG-SET在轴向预拉伸、内外压和径向电压联合作用下的力电非均匀偏场，对其轴对称变形进行深入分析。采用柱坐标系(r,\theta,z)，基于变形几何关系和力电耦合原理，建立相应的控制方程。在变形几何关系方面，对于轴对称变形，径向位移u仅为径向坐标r的函数，即u=u(r)，周向位移v=0，轴向位移w=w(r)。根据格林-拉格朗日应变张量的定义，可得到应变分量与位移的关系。径向应变E_{rr}=\frac{\partialu}{\partialr}+\frac{1}{2}(\frac{\partialu}{\partialr})^2，周向应变E_{\theta\theta}=\frac{u}{r}+\frac{1}{2}(\frac{u}{r})^2，轴向应变E_{zz}=\frac{\partialw}{\partialz}+\frac{1}{2}(\frac{\partialw}{\partialz})^2，径向-轴向剪切应变E_{rz}=\frac{1}{2}(\frac{\partialu}{\partialz}+\frac{\partialw}{\partialr}+\frac{\partialu}{\partialr}\frac{\partialw}{\partialz})。这些应变分量准确描述了FG-SET在轴对称变形下的几何变化情况。从力电耦合原理出发，软电弹性管的总能量包括弹性应变能W^e、电场能W^e和外力做功W^e。弹性应变能密度W^e根据Mooney-Rivlin模型确定，其表达式为W^e=C_{10}(I_1-3)+C_{01}(I_2-3)+D_1(J-1)^2，其中C_{10}、C_{01}、D_1为材料常数，I_1、I_2为应变不变量，J=\det(\boldsymbol{F})为体积比。应变不变量I_1=E_{rr}+E_{\theta\theta}+E_{zz}，I_2=E_{rr}E_{\theta\theta}+E_{rr}E_{zz}+E_{\theta\theta}E_{zz}-E_{rz}^2。电场能密度W^e=\frac{1}{2}\boldsymbol{D}\cdot\boldsymbol{E}，在柱坐标系下，电场强度E=-\frac{\partial\varphi}{\partialr}（假设电势\varphi仅为径向坐标r的函数），电位移矢量D=\epsilonE（\epsilon为介电常数）。外力做功包括轴向预拉伸力N、内外压力p_i（内压）和p_o（外压）以及径向电压V所做的功。轴向预拉伸力N在轴向位移上做功，其功密度为N\frac{\partialw}{\partialz}；内压p_i在径向位移上做功，功密度为p_iu(a)（u(a)为内径处的径向位移），外压p_o在径向位移上做功，功密度为-p_ou(b)（u(b)为外径处的径向位移）；径向电压V在电位移上做功，其功密度为VD。根据能量变分原理\delta(W^e+W^e+W^e)=0，对总能量关于位移u和w以及电势\varphi求变分。在求变分过程中，利用格林公式等数学工具对各项进行处理。对于弹性应变能项，根据应变与位移的关系，对W^e关于位移求偏导数，再乘以位移的变分并在管的体积内积分；对于电场能项，根据电场强度与电势的关系，对W^e关于电势求偏导数，再乘以电势的变分并在管的体积内积分；对于外力做功项，分别对轴向预拉伸力、内外压力和径向电压做功项关于相应的位移和电位移求偏导数，再乘以其变分并在管的体积内积分。通过一系列复杂的数学推导和化简，最终得到FG-SET在力电非均匀偏场作用下的轴对称变形的非线性控制方程。这些方程全面考虑了材料特性的径向变化、力电耦合效应以及各种偏场的作用，为研究FG-SET的力学行为和导波传播特性提供了重要的理论基础。4.2增量场控制方程