这个例子展示了如何解决弹性-静电耦合问题。
压电材料在外加电压下会变形。相反,压电材料的变形会产生电压。因此,压电部件的分析需要求解一组以挠度和电势为因变量的耦合偏微分方程。
在这个例子中,模型是一个两层悬臂梁,两层都由相同的聚偏氟乙烯(PVDF)材料制成。极化方向向下(负y-direction)在顶层,并在底层指向上面。典型的长厚比是100。当你在光束的上下表面之间施加电压时,光束就会发生偏转y-方向,因为一层变短,另一层变长。
平衡方程描述了固体的弹性行为:
在这里, 是应力张量,和 是物体的力向量。高斯定律描述了固体的静电行为:
是电位移,和 是分配的免费费用。将这两个PDE系统合并为一个系统:
对于二维分析, 有组件 而且 , 有组件 而且 .
材料的本构方程用应变张量和电场来定义应力张量和电位移矢量。对于平面应力条件下正交各向异性压电材料的二维分析,可以将这些方程写成
弹性系数, 电介电常数和 为压电应力系数。这些方程中的压电应力系数符合压电材料的常规符号z-direction(第三个方向)与材料的“极化”方向对齐。对于2-D分析,将“极化”方向与y设在。把应变向量写成x位移 而且y位移 :
把电场写成电势的形式 :
你可以把应变-位移方程和电场方程代入本构方程,得到一个用位移和电势导数表示的应力和电位移方程组。将得到的方程代入PDE方程组,得到一个包含位移导数和电势导数散度的方程组。作为下一步,排列这些方程以匹配工具箱所需的形式。
偏微分方程工具箱™需要一个椭圆方程系统以矢量形式表示:
或者是张量形式:
这里的重复指标意味着求和。对于本例中的二维压电系统,系统矢量 是
这是一个 系统。的梯度 是
有关以工具箱所需的格式指定系数的详细信息,请参阅:
这个例子中的c系数是一个张量。你可以用2 × 2块的3 × 3矩阵表示它:
为了将本构方程的项映射到工具箱所需的形式,将c张量和解的梯度写成这样的形式:
从这个方程,你可以将传统的本构系数映射到c矩阵所需的形式。电场方程中的负号被并入到c矩阵中,以匹配工具箱的惯例。
创建PDE模型。线弹性方程有三个分量,因此模型必须有三个方程。
模型= createpde(3);
创建几何图形并将其包含在模型中。
L = 100e-3;%光束长度,单位为米H = 1e-3;%梁的总高度H2 = h /2;%每层高度,单位为米topLayer = [3 4 0 L L 0 0 0 H2 H2];bottomLayer = [3 4 0 L L 0 -H2 -H2 0 0];gdm = [topLayer;bottomLayer]';G = decsg(gdm,R1 + R2的, (R1的;R2的]“);geometryFromEdges(模型中,g);
用面和边标签绘制几何图形。
图pdegplot(模型,“EdgeLabels”,“上”,...“FaceLabels”,“上”)包含(“x坐标,米”) ylabel (“坐标,米”)轴([- 1*L,1.1*L,-4*H2,4*H2]广场
指定梁层的材料属性。两层材料都是聚偏氟乙烯(PVDF),一种具有压电性能的热塑性聚合物。
E = 2.0e9;%弹性模量,N/m^2Nu = 0.29;泊松比G = 0.775e9;%剪切模量,N/m^2D31 = 2.2e-11;%压电应变系数,C/ND33 = -3.0e-11;
指定材料在恒定应力下的相对电介电常数。
relPermittivity = 12;
指定真空的电介电常数。
permittivityFreeSpace = 8.854187817620e-12;% F / mC11 = e /(1 - ^2);C12 = nu * c11;c2d = [C11 C12 0;C12 c11 0;0 0 g];pzeD = [0 d31;0 d33;0 0];
指定压电应力系数。
pzeE = c2d*pzeD;D_const_stress = [relPermittivity 0;0 relPermittivity] * permittivityFreeSpace;
将介电矩阵从恒应力转换为恒应变。
D_const_strain = D_const_stress - pzeD'*pzeE;
你可以把这21个系数看成是2 × 2块的3 × 3矩阵。的cij矩阵是这个矩阵的上三角形中的2 × 2块。
C11 = [c2d(1,1) c2d(1,3) c2d(3,3)];C12 = [c2d(1,3) c2d(1,2);c2d(3)汇集(2、3)];C22 = [c2d(3,3) c2d(2,3) c2d(2,2)];c13 = [pzeE(1,1) pzeE(1,2);pzeE (3,1) pzeE (2)];c23 = [pzeE(3,1) pzeE(3,2);pzeE (2, 1) pzeE (2, 2)];c33 = [D_const_strain(1,1) D_const_strain(2,1) D_const_strain(2,2)];Ctop = [c11(:); c12(:); c22(:); -c13(:); -c23(:); -c33(:)]; cbot = [c11(:); c12(:); c22(:); c13(:); c23(:); -c33(:)]; f = [0 0 0]'; specifyCoefficients(model,“米”,0,' d ',0,“c”ctop,“一个”,0,“f”f“脸”2);specifyCoefficients(模型,“米”,0,' d ',0,“c”芝加哥期货交易所,“一个”,0,“f”f“脸”1);
设置光束顶部(边缘1)的电压(溶液组件3)为100伏。
voltTop =应用边界条件(模型,“混合”,...“边缘”, 1...“u”, 100,...“EquationIndex”3);
通过将电压设置为0,指定光束的底部(边缘2)接地。
volbot =应用边界条件(模型,“混合”,...“边缘”2,...“u”,0,...“EquationIndex”3);
属性指定左边(边6和边7)是夹紧的x- - -y-位移(溶液组分1和2)到0。
clpleft =应用边界条件(模型,“混合”,...“边缘”者,...“u”[0 0],...“EquationIndex”1:2);
梁右侧的应力和电荷为零。因此,对边3和边4使用默认边界条件。
生成网格并求解模型。
msh = generateMesh(模型,“Hmax”5的军医);Result = solvepde(模型)
结果= StationaryResults with properties: nodesolution: [3605x3 double] XGradients: [3605x3 double] YGradients: [3605x3 double] ZGradients: [0x3 double] Mesh: [1x1 FEMesh]
在节点位置访问解决方案。第一列包含x偏转。第二列包含y偏转。第三列是电势。
rs = result. nodesolution;
求最小值y偏转。
feTipDeflection = min(rs(:,2));流('有限元尖端挠度为:%12.4e\n', feTipDeflection);
有限元尖端挠度为:-3.2900e-05
将这个结果与已知的解析解进行比较。
tipDeflection = -3*d31*100*L^2/(8*H2^2);流('分析尖端挠度为:%12.4e\n', tipDeflection);
解析尖端偏转为:-3.3000e-05
画出偏转分量和电势。
varsToPlot = char(“X-Deflection,米”,...“y偏转,米”,...电势,伏特);为i = 1:size(varsToPlot,1) figure;pdeplot(模型,“XYData”rs(:我),“轮廓”,“上”)标题(varsToPlot(我:))%缩放轴,使它更容易查看轮廓xlabel([0, L, -4*H2, 4*H2])“x坐标,米”) ylabel (“坐标,米”)轴广场结束
黄禹锡,朴玄哲。压电传感器和执行器的有限元建模张仁杂志31日,5号(1993年5月):930-937。
Pieford, V。压电活性结构的有限元建模博士羞辱。,Universite Libre De Bruxelles, 2001.