基于zig-zag假设的压电阻尼层合结构机电耦合建模与分析

2022-05-11王雄高英山张顺琦窦伟元

西北工业大学学报 2022年2期

王雄, 高英山, 张顺琦,, 窦伟元

1.榆林学院能源工程学院, 陕西榆林 719000; 2.上海大学机电工程与自动化学院, 上海 200072;3.北京交通大学机械与电子控制工程学院, 北京 100044

薄壁结构由于质量轻,广泛应用于航空、航天领域,如太阳能帆板、飞机的机翼和舱体等。单一功能的薄壁结构极易发生振动与变形,且结构阻尼低、不稳定[1]。压电阻尼层合结构具有被动阻尼和主动致动性能,也称为主被动结构,具有良好的力学特性,广泛应用于形状控制、振动控制、噪声控制等方面[2]。黏弹性阻尼材料特殊的物理性能,使得主被动阻尼层合结构拥有复杂的层间物理特性,增加了主被动阻尼层合结构动力学建模难度,且不易实现。

最早Kerwin等[3]对黏弹性材料的阻尼参数进行了定量分析,并提出了阻尼材料的损耗因子。王金朝等[4]采用多输入多输出锤击法对约束阻尼板进行模态实验并与有限元结果互相验证。陈威和夏利娟[5]基于模态应变能法推导出模态损耗因子的修正公式,结果表明阻尼加筋板的阻尼弹性模量频变效应对模态损耗因子影响明显,阻尼材料的弹性模量可取为定值。黄微波等[6]研究了阻尼层合结构的层间厚度对约束阻尼结构振动性能带来的影响,结果表明约束层厚度与基层厚度相等时结构复合损耗因子最大,并且随着基层厚度与约束层厚度之比逐渐减小,厚度比对结构损耗因子的影响也变小。Jin等[7]对黏弹性阻尼夹层复合矩形板的振动和阻尼进行了分析,基于Reddy的高阶剪切理论,并且考虑了弯曲-拉伸、弯曲-扭转、拉伸-扭转以及泊松效应等不同情况下的材料耦合情况。

阻尼层合圆柱状壳体结构在工程中的应用也很广泛。Mokhtari等[8]利用Rayleigh-Ritz方法提出了夹层圆柱壳的一个计算公式,讨论了关键参数对动力性能的影响。Yang等[9]研究了中间夹层为黏弹性材料的复合材料板的超音速颤振阻尼效应。Huang等[10]建立了2种基于剪切和压缩阻尼机制的夹层结构分析模型,并说明了这2种不同阻尼机制的适用范围。Ebrahimi和Barati[11]建立了一个非局部应变梯度板模型,研究了黏弹性、非局部参数、电压等因素对压电阻尼层合板振动特性的影响。

通过文献调研可知,多数文献集中于阻尼层合板结构的研究,对于在工程中应用前景很广的压电阻尼层合壳体结构研究较少。本文拟构建基于zig-zag板壳假设的“机-电”耦合压电阻尼层合结构有限元模型。用文献数据验证模型的正确性,并进一步研究材料增强角度、阻尼层厚度和结构曲率对压电阻尼层合结构频率和损耗系数的影响,可以为压电阻尼层合结构的减振降噪及结构设计提供必要的理论基础。

1 压电阻尼层合结构

图1 结构变形与自由度表示

则板壳结构上任意点的位移向量u可以用7个运动学参数表示为

u=Zv=

(1)

式中:v为位移向量;Z为变换矩阵,其中zv和ze值取决于计算的结构位置,如表1所示。

表1 zv和ze取值

2 应变位移关系

假设结构在厚度方向上不可压缩,这时Green-Lagrange应变张量主要由面内应变εαβ和横向剪切应变εα3组成,可以写为[12]

(2)式与(3)式中各应变分量与位移场之间的关系可以表示为:

1) 对于计算结构位置Θ3≤z1与Θ3≥z2

(4)

(5)

(6)

(7)

(8)

计算Θ3≤z1时d=z1,计算Θ3≥z2时d=z2。

2) 对于计算结构位置z1<Θ3

(9)

(10)

(11)

(12)

(13)

式中

(14)

3 动力学有限元模型

3.1 本构方程

压电阻尼层合结构机电耦合有限元本构方程的矩阵形式为

式中，σ,ε,D和E分别表示应力向量、应变向量、电位移向量和电场向量。此外,c为弹性常数矩阵,其与压电常数矩阵d和e之间的关系可以用矩阵形式表示为e=dc,χ为介电常数矩阵。

图2 局部坐标系和材料坐标系

(17)

3.2 动力学有限元方程

采用Hamilton原理来推导压电阻尼层合结构的动力学方程,可以表述为从t1到t2的虚功积分等于零,即

(18)

式中，δT为压电阻尼层合结构的动能微分,可以表示为

(19)

(20)

外力功的变分可以表示为

(21)

将(19)～(21)式代入(18)式可以得到压电阻尼层合结构的机电耦合动力学有限元方程

式中：q为节点自由度向量；Muu,Kuu,Kuφ和Kφφ分别为总质量矩阵、刚度矩阵、耦合刚度矩阵和压电矩阵;Fue和Gφe为外力和电场向量。

3.3 阻尼材料模型与复刚度矩阵

通常认为阻尼材料的弹性模量是复弹性模量,表示为

Gv=GR+iGI=GR(1+iηv)

(24)

弹性模量实部GR称为储能模量,弹性模量虚部GI称为损耗模量,ηv表示阻尼材料的损耗因子。则结构的总刚度矩阵Kuu可以写为

(25)

压电阻尼层合结构的自然频率f和损耗系数η可表达为

(26)

式中，λn是第n阶复特征角频率。弹性层的阻尼矩阵用CS表示,根据Moita等[14]使用的方法来求解阻尼材料的阻尼矩阵CV,可以表示为

(27)

根据(25)式和(27)式,(22)式可以整理为

(28)

4 有限元仿真及数据分析

4.1 模型验证

采用如图3所示的对称阻尼层合板与文献[15]的固有频率和损耗系数结果进行对比。板的尺寸为348 mm×304.8 mm×0.762 mm,中间阻尼层厚度为0.254 mm。阻尼材料的损耗系数ηv为常数0.5,阻尼层与弹性层的材料参数如表2所示。

图3 阻尼层合结构示意图

表2 阻尼结构材料参数

采用12×10网格划分,八节点矩形单元,两短边简支(simply supported:S),两长边固支(clamped:C),即CSCS。且认为压电材料在无外电场情况下为闭环状态。则结构前5阶固有频率f和结构损耗系数η数据对比如表3所示,仿真结果与文献结果吻合较好。

表3 结果数据对比验证

4.2 压电阻尼层合壳体参数化研究

压电阻尼层合壳体结构如图4所示,半径用R表示,压电层厚度为1 mm,复合弹性层(T300/976)的厚度为1.5 mm,宽度W=160 mm,中性面弧长为314 mm。复合材料(T300/976)和压电层的材料参数在表4中列出,且认为它们为弹性层。阻尼材料参数如表2所示。在Θ1和Θ2方向上分别划分了8个和10个网格,采用七自由度八节点的矩形单元。

图4 压电阻尼层合壳示意图

表4 压电阻尼层合壳结构材料参数

壳体两直边固支,其余简支,即CSCS。令结构半径R=100 mm,阻尼层的厚度为1 mm,复合弹性层(T300/976)的增强角度β变化范围为0°～90°,结果如表5所示。当增强角度增加到55°时,第一阶结构频率f1逐渐增加,随着增强角度继续增加,频率开始缓缓降低,而第一阶损耗因子η1一直处于下降状态;角度从0°逐渐变化到90°时,第二阶频率f2一直增加,第二阶损耗因子η2整体呈降低趋势。

表5 增强角度β对结构频率f和损耗因子η的影响

图5 曲率Ω对结构频率f和损耗因子η的影响

同时固定压电阻尼层合壳体的2条直角边,即边界条件为CFCF。结构的半径R为100 mm,增强角度β=0°保持不变,而阻尼层厚度从0.1 mm逐渐增加至5 mm。结果如图6所示,当阻尼层厚度从0.1 mm增加到约1 mm时,频率降低非常快,而当阻尼层厚度继续增加时,频率开始缓慢增加;而损耗系数在阻尼层厚度从0.1～5 mm时一直减小,但减小的速度越来越慢。即阻尼层厚度对结构特性的影响并不是单调递增或单调递减的。

图6 阻尼层厚度对结构频率f和损耗因子η的影响

4.3 静力学和动力学仿真

本节研究边界条件为CFFF的压电阻尼层合壳体(见图4)的静力和动力学特性。令壳的半径R=100 mm,阻尼层的厚度为0.5 mm和增强角度β=0°。其他尺寸、材料、单元类型和网格划分与4.2节相同。压电片黏贴面的电势为零,且整个上表面都是等电位的。对壳体结构施加150 V电压,壳体结构在Θ2方向上中心线位移如图7静力学响应所示,壳体距离固定端越远,节点位移就越大。

图7 压电阻尼层合壳体的静力和动力学响应

对结构施加电压后,结构会在一个范围内来回振动。由图5可知壳体的频率为35.8 Hz,取0.1 s即3.58个振动周期来研究该壳体的动力学特性。不考虑阻尼情况下自由端中间节点的环向振动响应如图7所示。现实中结构不会一直振动下去,结构会逐渐趋于平稳,这个过程是阻尼效应作用的结果。不同阻尼效应情况下壳体的自由端中间节点的环向振动,结果如图8所示。图例中CS表示只考虑了结构自身的阻尼效应;CS+CV则表示既考虑了结构自身的阻尼效应又考虑了黏弹性结构的阻尼效应。可以看出黏弹层阻尼可以使结构的振动衰减很快,对抑制结构振动有很好的效果。