Loading [MathJax]/jax/output/SVG/jax.js

PVDF压电传感器及敏感单元设计关键技术

张旭 覃双 杨舒棋 彭文杨 赵锋 于君 钟斌

王志鹏, 韩志军, 王龙飞. 基于Galerkin法研究应力波作用下复合材料板的动力学失稳[J]. 高压物理学报, 2021, 35(5): 054204. doi: 10.11858/gywlxb.20210705
引用本文: 张旭, 覃双, 杨舒棋, 彭文杨, 赵锋, 于君, 钟斌. PVDF压电传感器及敏感单元设计关键技术[J]. 高压物理学报, 2020, 34(2): 023401. doi: 10.11858/gywlxb.20190796
WANG Zhipeng, HAN Zhijun, WANG Longfei. Dynamic Instability of Composite Plate under Stress Wave Based on Galerkin Method[J]. Chinese Journal of High Pressure Physics, 2021, 35(5): 054204. doi: 10.11858/gywlxb.20210705
Citation: ZHANG Xu, QIN Shuang, YANG Shuqi, PENG Wenyang, ZHAO Feng, YU Jun, ZHONG Bin. Key Design Techniques for PVDF Sensitive Element Used in Dynamic Compression Experiments[J]. Chinese Journal of High Pressure Physics, 2020, 34(2): 023401. doi: 10.11858/gywlxb.20190796

PVDF压电传感器及敏感单元设计关键技术

doi: 10.11858/gywlxb.20190796
基金项目: 国防科工局技术基础项目(JSZL2016212C001);科学挑战计划(TZ20118001)
详细信息
    作者简介:

    张 旭(1972-),男,研究员,博士生导师,主要从事流体动力学研究. E-mail:caepzx@sohu.com

  • 中图分类号: O521.3; TJ55

Key Design Techniques for PVDF Sensitive Element Used in Dynamic Compression Experiments

  • 摘要: 研究了聚偏二氟乙烯(PVDF)柔性压电传感器敏感单元设计关键技术和工艺。通过传感器敏感单元设计、电极制作、极化和封装研究实现功能,结合动态冲击压缩实验,对自主研制的PVDF传感器在不同压力段进行标定和A类不确定度评价,实验数据的A类不确定度优于10%,说明所研制的PVDF压电传感器精度高、重复性好,能够满足0.3~10.0 GPa动态冲击压力测量应用需求。目前产品设计水平已达到6级,后续将进一步开展更高压段有效性的应用设计、温度补偿效应研究和工程化产品设计。

     

  • 复合材料板具有较高的强度质量比、良好的耐腐蚀性和优异的可设计性,被广泛应用于航空航天和工业制造等领域[1]。在实际使用中复合材料板经常受到不同形式的冲击荷载,从而产生振动和屈曲问题,因此冲击载荷作用下复合材料板的动力稳定性问题备受关注。

    近年来关于复合材料板的研究越来越多,尤其是冲击荷载作用下复合材料板的动力稳定性问题研究[2],对工程部件结构设计和使用具有重要的意义。Sun等[3]研究了在加热环境中应力波对功能梯度圆柱壳轴向冲击屈曲的影响;毛柳伟等[4]对弹性直杆在应力波作用下的动力分叉屈曲进行了分析与探讨,提出了求解应力波作用下直杆动力屈曲的数值方法;Lepik[5]讨论了在应力波影响下轴向压缩的弹塑性梁的屈曲;Zhang等[6]分析了不确定初始几何缺陷对薄板屈曲的影响;Abdelaziz等[7]利用双曲线剪切变形理论,分析了在各种边界条件下复合材料板的弯曲变形和屈曲;Kouchakzadeh等[8]采用线性和旋转弹簧的均匀分布来模拟边界条件,对矩形层压复合板的屈曲进行了分析;Czapski等[9]通过数值和实验方法,研究了残余应力对压缩至破坏期间薄壁层压板屈曲性能的影响。

    在实际工程中,复合材料板多应用于振动环境,其在动力响应下的动态特性和振动分析是必不可少的,因而对该类材料的振动屈曲研究至关重要。Kuo[10]研究了两种非均匀分布纤维复合材料板的振动屈曲问题,Villarreal等[11]对典型正交异性板的本征频率和振动屈曲进行了理论分析,Eftekhari等[12]提出了通过组合应用有限元方法和微分正交方法求解矩形板的振动屈曲问题,Rehman等[13]探讨了壳体结构的缺陷和损坏对结构振动屈曲的影响,Sayyad等[14]将三角剪切变形理论应用于复合板的变形和振动屈曲研究。

    关于复合材料板的振动屈曲问题已开展了较多的研究,但大多未考虑应力波效应对振动屈曲的影响,而动力屈曲一般与应力波相联系且具有局部发生的特点,研究含初始缺陷的复合材料板能更好地揭示实际工程中复合材料板在不同工况下发生动力屈曲的机理。基于此,本研究利用Kirchhoff薄板理论和Hamilton原理,考虑应力波效应,建立含初始几何缺陷的四边简支复合材料板的振动控制方程,得到板的屈曲临界荷载表达式,在此基础上通过数值计算讨论初始几何缺陷、振型函数初相位、铺层角度、屈曲模态阶数和铺层层数对复合材料板振动屈曲临界荷载的影响,为工程实际提供理论依据。

    复合材料板在x = La处为固定边界条件,其余3边为简支边界条件,在z = 0的中性面上受x方向的面内阶跃荷载N作用,如图1所示。板在z方向上含初始几何缺陷w1,且xyz方向的位移分别为uvw。复合材料板的长度、宽度和厚度分别为LaLbh,由n层单层板组成,θ为单层板的铺层角度,即纤维材料铺设方向与x方向的夹角(见图1)。

    图  1  复合材料板结构示意图
    Figure  1.  Schematic diagram of composite plate structure

    根据Kirchhoff薄板理论及经典弹性理论,复合材料板的位移与应变、弹性曲面的曲率和扭率的表达式为

    {w=w1+w0u=u0z(ww1)xv=v0z(ww1)y (1)
    {ε0x=u0xε0y=v0yγ0xy=u0y+v0x (2)
    {κx=2(ww1)x2κy=2(ww1)y2κxy=22(ww1)xy (3)
    {εx=ε0x+zκxεy=ε0y+zκyγxy=γ0xy+zκxy (4)

    式中:u0v0w0分别为复合材料板在xyz方向上的中面位移,ε0xε0yγ0xy分别为复合材料板中面应变分量,κxκyκxy为中面的曲率和扭率,εxεyγxy为复合材料板任意一点的应变。

    板的内力(NxNyNxy)与内力矩(MxMyMxy)之间的关系为

    [NxNyNxyMxMyMxy]=[A11A12A16B11B12B16A12A22A26B12B22B26A16A26A66B16B26B66B11B12B16D11D12D16B12B22B26D12D22D26B16B26B66D16D26D66][ε0xε0yγ0xyκxκyκxy] (5)

    式中:AijBijDij(i,j=1,2,6)分别表示板的拉伸刚度、耦合刚度和弯曲刚度系数[15],表达式如下

    {Aij=nk=1¯Qkij(hkhk1)Bij=12nk=1¯Qkij(h2kh2k1)Dij=13nk=1¯Qkij(h3kh3k1) (6)

    式中:¯Qkij为复合材料板第k层的偏轴刚度系数

    ¯Q=P1Q(P1)T (7)

    式中:P为坐标转换矩阵,Q为刚度矩阵。

    P=[cos2θsin2θ2sinθcosθsin2θcos2θ2sinθcosθsinθcosθsinθcosθcos2θsin2θ] (8)
    Q=[Q11Q120Q12Q22000Q66] (9)

    考虑材料为正交各向异性材料,设E1E2G12μ12μ21分别为板材料xy方向的拉压弹性模量、剪切弹性模量、主泊松比和副泊松比,则有

    Q11=E11μ12μ21,Q22=E21μ12μ21,Q12=μ12E21μ12μ21=μ21E11μ12μ21,Q66=G12 (10)

    复合材料板在左端受面内冲击荷载N作用(见图1)时,应力波沿x方向在板内传播,其应力变化如图2所示。

    图  2  应力波传播示意图
    Figure  2.  Schematic diagram of stress wave propagation

    当应力波传播至波阵面位置Lcr(临界长度)时,板发生振动屈曲,板的内力Nt和应力波波速c分别表示为

    Nt={N0xLcr0x>Lcr,c=A11ρh (11)

    板发生振动屈曲时的变形能可以表示为

    U=12Lcr0Lb0(Nxε0x+Nyε0y+Nxyγ0xy+Mxκx+Myκy+Mxyκxy)dxdy=12Lcr0Lb0[Nxu0x+Nyv0y+Nxy(u0x+v0y)+Mx(2w0x2)+My(2w0y2)+Mxy(22w0xy)]dxdy (12)

    发生振动屈曲时的动能(考虑转动惯量)可以表示为

    T=12h/2h/2Lcr0Lb0ρ(k)[(ut)2+(vt)2+(w0t)2]dxdydz=12Lcr0Lb0[I0(u0t)2+I0(v0t)2+I0(w0t)22I1u0t2w0xt2I1v0t2w0yt+I2(2w0xt)2+I2(2w0yt)2]dxdy (13)

    式中:(I0,I1,I2)=Nkk=1hkhk1ρ(k)(1,z,z2)dzρ(k)为第k层材料的密度。

    板发生振动屈曲时的外力功可以表示为

    W=12Lcr0Lb0Nt(wx)2dxdy (14)

    考虑Hamilton变分原理,即

    δt2t1(TU+W)dt=0 (15)

    将式(12)~式(14)代入式(15)并进行变分计算,可得

    Nxx+NxyxI02u0t2+I13w0xt2=0 (16)
    Nyy+NxyyI02v0t2+I13w0yt2=0 (17)
    2Mxx2+2Myy222MxyxyNt(2w0x2+2w1x2)=I02w0t2+I1(3u0xt2+3v0xt2)I2(4w0x2t2+4w0y2t2) (18)

    对于正交各向异性正规对称正交铺设的复合材料板,其刚度矩阵满足[16]

    D16=D26=0,A16=A26=0,Bij=0 (19)

    根据 Kirchhoff 薄板理论及经典弹性理论,薄板中面在变形过程中没有伸长变形,将板的本构关系代入式(16)~式(18)中,略去含u0v0的项,得到复合材料板在面向阶跃荷载激励下的控制方程

    D114w0x4+(2D12+4D66)4w0x2y2+D224w0y4+Nt(2w0x2+2w1x2)I02w0t2+I2(4w0x2t2+4w0y2t2)=0 (20)

    设三边简支和应力波传播到Lcrw0(Lcr,y,t)=w0(Lcr,y,t)=0)时的振型函数[17]

    w0(x,y,t)=Rij[siniπxLcr+ii+1sin(i+1)πxLcr]sinjπyLbsin(ωt+φ) (21)

    由屈曲模态确定的缺陷分布形式是板结构最有可能发生的屈曲形式,能够很好地确定结构的缺陷敏感性[18]。对于复合材料板在制造过程中出现的初始几何缺陷,引入屈曲模态的ε倍变形作为初始几何缺陷[19],可以表示为

    w1(x,y)=εRij[siniπxLcr+ii+1sin(i+1)πxLcr]sinjπyLb (22)

    式中:ij为屈曲模态阶数,i,j=1,2,3,Rij为板的第(i, j)阶模态幅值;ε表示初始几何缺陷系数。

    根据式(21)和式(22),利用棣莫弗公式对控制方程式(20)中的各项进行求导并化简,得到

    D11(πLcr)4[i4+(i+1)4]sin(ωt+φ)+(2D12+4D66)(πLcr)2[i2+(i+1)2](jπLb)2sin(ωt+φ)+2D22(jπLb)4sin(ωt+φ)Nt(πLcr)2[i2+(i+1)2][sin(ωt+φ)+ε]+2I0ω2sin(ωt+φ)+I2ω2{(πLcr)2[i2+(i+1)2]+2(jπLb)2}sin(ωt+φ)=0 (23)

    根据式(23)可以得到N的表达式

    N={D11(πLcr)4[i4+(i+1)4]+(2D12+4D66)(πLcr)2[i2+(i+1)2](jπLb)2+2D22(jπLb)4+2I0ω2+I2ω2{(πLcr)2[i2+(i+1)2]+2(jπLb)2}}sin(ωt+φ)(πLcr)2[i2+(i+1)2][sin(ωt+φ)+ε] (24)

    板发生屈曲时,临界条件为ω=0[20],代入式(24)可得振动屈曲临界荷载为

    Ncr=D11(πLcr)2[i4+(i+1)4]sinφ[i2+(i+1)2](sinφ+ε)+(2D12+4D66)(jπLb)2sinφ(sinφ+ε)+2D22(j2πL2b)2L2crsinφ[i2+(i+1)2](sinφ+ε) (25)

    利用MATLAB数值分析应力波未反射时初始几何缺陷、初相位、铺层角度、屈曲模态阶数、铺设厚度以及铺层层数对复合材料板振动屈曲临界荷载的影响,使用的材料参数见表1[21]

    表  1  复合材料板参数[21]
    Table  1.  Material parameters of composite plate[21]
    E1/GPaE2/GPaG12/GPaμ12La/mLb/m
    140.08.65.00.350.600.50
    下载: 导出CSV 
    | 显示表格

    图3显示了复合材料板的模态阶数i为1、2、3,j为1时的屈曲模态。当板的x方向模态增大时,x方向的屈曲模态第一峰值增大且波数增加,而y方向的屈曲模态呈正对称分布。模态阶数的增加使板振动的屈曲模态变得更复杂。

    图  3  x方向模态取值增大时板的屈曲模态
    Figure  3.  Buckling mode of composite plate with increasing mode value in x direction

    设置7组算例,分别以初始几何缺陷、初相位、铺层角度、xy两个方向屈曲模态阶数、铺层层数及铺设厚度为变量进行算例分析,研究以上因素对板振动屈曲临界荷载的影响,算例参数见表2

    表  2  算例分析参数表
    Table  2.  Example analysis parameter table
    GroupInitial defect
    coefficient
    Order of modeLaying angle/(°)Initial phaseNumber of
    layers laid
    Thickness
    of the plate/m
    x directiony direction
    AVariablei = 1j = 1[0, 0, 0, 0, 0]π/250.01
    B0.1Variablej = 1[0, 0, 0, 0, 0]π/250.01
    C0.1i =1Variable[0, 0, 0, 0, 0]π/250.01
    D0.1i =1j = 1Variableπ/250.01
    E0.1i =1j = 1[0, 0, 0, 0, 0]Variable50.01
    F0.1i =1j = 1[0, 0, 0, 0, 0]π/2Variable0.01
    G0.1i =1j = 1[0, 0, 0, 0, 0]π/25Variable
    下载: 导出CSV 
    | 显示表格

    将A组数据代入式(25)中,可以得到不同初始缺陷系数对复合材料板振动屈曲的影响,如图4所示。由Lcr-Ncr曲线可知:在应力波传播过程中,Ncr呈指数下降,分为敏感区和非敏感区。应力波在Lcr < 0.4 m区域传播时,Ncr的变化较陡峭,该区域为敏感区;应力波在Lcr > 0.4 m区域传播时,Ncr的变化趋于平缓,该区域为非敏感区,因此敏感分界点为0.4 m。当选取的初始缺陷系数增大时,临界荷载Ncr也随之增大。在敏感区,初始缺陷系数对临界荷载Ncr的影响较大,且随应力波传播呈减小趋势。此外,初始缺陷系数对非敏感区的影响较小。图4表明,初始几何缺陷系数越大,板越容易发生屈曲。

    图  4  不同初始缺陷系数条件下NcrLcr的关系曲线
    Figure  4.  Relationship between Ncr and Lcr under different initial defect coefficients

    将B组数据代入式(25)中,可以得到不同x方向模态阶数对复合材料板振动屈曲的影响,如图5所示。由Lcr-Ncr曲线可知:当选取的x方向模态阶数增大时,临界荷载Ncr随之明显增大。在敏感区,x方向模态阶数对临界荷载Ncr的影响很大,并随应力波的传播不断减小,到达非敏感区之后影响较小并趋于稳定。图5表明,x方向模态阶数的增加会显著增大板的屈曲临界荷载。

    图  5  x方向模态阶数不同时NcrLcr的关系曲线
    Figure  5.  Relationship between Ncr and Lcr with different order of modes in x direction

    将C组数据代入式(25)中,可以得到不同的y方向模态阶数对复合材料板振动屈曲的影响,如图6所示。由Lcr-Ncr曲线可知:当选取的y方向模态阶数增大时,临界荷载Ncr也随之变大。在应力波传播过程中,在敏感区y方向模态阶数对临界载荷基本没有影响,而在非敏感区有极小的影响。图6表明,y方向模态阶数的变化对板屈曲临界荷载基本没有影响。

    图  6  y方向模态阶数不同时NcrLcr的关系曲线
    Figure  6.  Relationship between Ncr and Lcr with different order of modes in y direction

    将D组数据代入式(25)中,得到不同铺层角度对复合材料板振动屈曲的影响,如图7所示。由Lcr-Ncr曲线可知:在敏感区,不同的铺设角度对临界荷载Ncr的影响较大,且随应力波的传播不断减小,到达非敏感区后趋于稳定。图7表明,铺层角度小的单层板的层数越多,板的临界屈曲荷载越大,说明复合材料板的铺设角度直接影响板的屈曲临界荷载。

    图  7  不同铺层角度条件下NcrLcr的关系曲线
    Figure  7.  Relationship between Ncr and Lcr under different laying angles

    将E组数据代入式(25)中,得到不同初相位对复合材料板振动屈曲的影响,如图8所示。由Lcr-Ncr曲线可知:振型函数的初相位越大,对应的临界荷载越大。在敏感区,振型函数的初相位对临界荷载Ncr的影响较小,且随应力波的传播不断减小;到达非敏感区之后,影响趋于平缓。图8表明,振型函数的初相位越大,板的屈曲临界荷载越大。

    图  8  不同初相位条件下NcrLcr的关系曲线
    Figure  8.  Relationship between Ncr and Lcr under the condition of initial phase of different mode functions

    将F组数据代入式(25)中,得到不同铺层层数对复合材料板振动屈曲的影响,如图9所示。由Lcr-Ncr曲线可知:当按照不同层数铺设时,敏感区的临界荷载Ncr的变化较大,且随应力波的传播不断减小;到达非敏感区之后变化较小并趋于平缓。图9表明,对于厚度固定且对称铺设的板,当铺设层数达到7时,其屈曲荷载随层数增加趋于稳定。

    图  9  不同铺层层数下NcrLcr的关系曲线
    Figure  9.  Relationship between Ncr and Lcr under different laying modes

    将G组数据代入式(25)中,得到不同铺设厚度对复合材料板振动屈曲的影响,如图10所示。由Lcr-Ncr曲线可知:板的铺设厚度越大,对应的临界荷载越大。在敏感区,不同的板厚对临界荷载Ncr的影响很大,且随应力波的传播不断减小,到达非敏感区后趋于稳定。图10表明,复合材料板的铺设厚度将直接决定板的屈曲临界荷载。

    图  10  不同铺设厚度下NcrLcr的关系曲线
    Figure  10.  Relationship between Ncr and Lcr under different thicknesses

    基于Kirchhoff薄板理论和Hamilton变分原理,建立了具有初始几何缺陷的四边简支复合材料板的振动控制方程。采用伽辽金法,选取符合边界条件的振型函数求解控制方程,得到屈曲临界载荷表达式。数值计算结果表明:应力波在未发生反射前的传播过程中,复合材料板的振动屈曲临界载荷随着临界长度的增大、铺设厚度的减小、初始几何缺陷系数的增大、振型函数初相位的减小而减小;复合材料板的各层铺层角度与荷载作用方向的夹角越小,屈曲临界载荷越大,当对称铺设层数达7层时,临界荷载趋于稳定。研究结果可为工程中复合材料板的结构设计与应用提供一定的参考。

  • 图  敏感单元形状设计

    Figure  1.  Sensitive element design

    图  PVDF薄膜电极结构和热极化装置

    Figure  2.  Schematic diagram of PVDF film electrode and thermal polarization

    图  20μm厚PVDF薄膜的热极化

    Figure  3.  Effect of polarized electric field on 20μm thick PVDF film

    图  PVDF传感器结构

    Figure  4.  PVDF sensor structure

    图  压电薄膜传感器电荷模式等效电路

    Figure  5.  The equivalent circuit of electric charge mode of PVDF film sensor

    图  测试系统测量电路

    Figure  6.  Measurement circuit of PVDF film sensor test system

    图  对称碰撞实验示意图

    Figure  7.  Experimental scheme of symmetrical impact

    图  30μm-PVDF传感器低压段标定曲线

    Figure  8.  Calibration curve of 30μm-PVDF in low pressure zone

    图  30μm-PVDF传感器高压段标定曲线

    Figure  9.  Calibration curve of 30μm-PVDF in high pressure zone

    表  1  Al和45钢的Hugoniot参数

    Table  1.   Hugoniot parameters of Al and 45 steel

    Materialρ0/(g·cm–3)c0/(km·s–1)λ
    Al2.7855.3281.338
    45 steel7.8504.4831.332
    下载: 导出CSV

    表  2  电荷法实验结果

    Table  2.   Experimental results by charge method

    Exp.No.Material of flyeru0/(m·s−1)p/GPaUmax/V(Q/A)max/(μC·cm−2d33/(pC·N−1)Pr/(μC·cm−2
    0830-1LY12 Al48.4340.361 20.598 670.665 212.08.279 7
    0830-245 steel83.3841.473 01.097 901.219 910.87.528 9
    0830-345 steel109.5571.942 81.367 271.519 210.57.878 1
    0830-445 steel129.5362.303 81.433 011.592 215.57.987 1
    0831-1LY12 Al126.2210.950 60.813 580.904 015.210.301 7
    1128-1LY12 Al47.1370.351 50.459 010.510 015.17.806 1
    1030LY12 Al628.5315.395 01.636 011.817 815.66.934 0
    1026PMMA1 004.7042.426 11.120 071.244 515.37.832 2
    1107LY12 Al970.2758.069 31.883 432.092 711.47.216 3
    1108LY12 Al1 250.44110.725 52.357 622.619 612.88.287 2
    下载: 导出CSV
  • [1] KAWAI H. The piezoelectricity of poly (vinylidene fluoride) [J]. Japanese Journal of Applied Physics, 1969, 8(7): 975–976. doi: 10.1143/JJAP.8.975
    [2] BAUER F. Method and device for polarizing ferroelectric materials: 4611260 [P]. 1986–09–09.
    [3] BAUER F, LICHTENBERGER A. Use of PVF2 shock gauges for stress measurements in Hopkinson bar [C]// SCHMIDT S C, HOLMES N C. Shock Waves in Condensed Matter. California, 1987: 751–755.
    [4] OBARA T, BOURNE N K, MEBAR Y. The construction and calibration of an inexpensive PVDF stress gauge for fast pressure measurements [J]. Measurement Science and Technology, 1995, 6(4): 345. doi: 10.1088/0957-0233/6/4/001
    [5] SALISBURY D A, WINTER R E, TAYLOR P, et al. The response of foams to shock compression [J]. AIP Conference Proceedings, 2000, 505(1): 197–200.
    [6] FIELD J E, WALLEY S M, PROUD W G, et al. Review of experimental techniques for high rate deformation and shock studies [J]. International Journal of Impact Engineering, 2004, 30(7): 725–775. doi: 10.1016/j.ijimpeng.2004.03.005
    [7] GAMA B A, LOPATNIKOV S L, GILLESPIE Jr J W. Hopkinson bar experimental technique: a critical review [J]. Applied Mechanics Reviews, 2004, 57(4): 223–250. doi: 10.1115/1.1704626
    [8] YANG L M, SHIM V P W. An analysis of stress uniformity in split Hopkinson bar test specimens [J]. International Journal of Impact Engineering, 2005, 31(2): 129–150. doi: 10.1016/j.ijimpeng.2003.09.002
    [9] 席道瑛, 郑永来. PVDF压电计在动态应力测量中的应用 [J]. 爆炸与冲击, 1995, 15(2): 174–179.

    XI D Y, ZHENG Y L. Application of PVDF gauges to dynamical stress measurements [J]. Explosion and Shock Waves, 1995, 15(2): 174–179.
  • 加载中
图(9) / 表(2)
计量
  • 文章访问数:  9588
  • HTML全文浏览量:  3485
  • PDF下载量:  53
出版历程
  • 收稿日期:  2019-06-26
  • 修回日期:  2019-07-14
  • 刊出日期:  2020-02-25

目录

/

返回文章
返回