• ISSN 1001-1455  CN 51-1148/O3
  • EI Compendex、CA收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊
高级检索 E-mail Alert

激波管聚酯膜片变形过程分析

何起光 张伟 陈小伟 徐剑鹏

引用本文:
Citation:

激波管聚酯膜片变形过程分析

    作者简介: 何起光(1994- ),男,硕士,478440913@qq.com;
    通讯作者: 张伟, zhdawei@hit.edu.cn
  • 中图分类号: O347.5; O354

Analysis on the deformation process of PET shock tube diaphragm

    Corresponding author: ZHANG Wei, zhdawei@hit.edu.cn ;
  • CLC number: O347.5; O354

  • 摘要: 使用聚酯薄膜作为激波管膜片,通过施加不同压力的激波管实验,获得了膜片厚度及多张膜片的组合方式对膜片所能承受最大压力的影响。利用高速相机对激波管膜片从开始变形到破坏的全过程进行拍摄,使用三维DIC软件获得膜片在变形过程中的位移场。实验发现了膜片会出现圆弧反翘并快速破坏的特别现象,并以此为特征将变形过程分为2个阶段。给出任意厚度膜片第1阶段圆弧变形的数学规律及第2阶段圆弧反翘的形状特征,以及全过程中膜片厚度变化的数学规律。
  • 图 1  实验装置示意图及照片

    Figure 1.  Schematic diagram and photo of experimental equipment

    图 2  膜片的应力应变曲线

    Figure 2.  Stress-strain curves of the diaphragm

    图 3  不同厚度及组合方式的膜片的破坏压力

    Figure 3.  Maximum bearing pressure of diaphragm with different thickness and combination

    图 4  三种厚度的膜片在不同压力条件下的变形

    Figure 4.  Deformation of diaphragms with three different thicknesses under different pressures

    图 5  DIC处理后所得不同压力条件下0.2 mm厚膜片变形的典型结果

    Figure 5.  Deformation of diaphragm with 0.2 mm in thickness by 3D-DIC analysis

    图 6  三种不同厚度膜片的变形轮廓曲线

    Figure 6.  Profiles of deformation of diaphragms with different thicknesses

    图 7  膜片中心挠度与压力曲线

    Figure 7.  Relation between maximum deflection of diaphragm and pressure

    图 8  膜片厚度减小量的分布曲线

    Figure 8.  Distribution curves of thickness reduction

    图 9  膜片中心的厚度减小量与压力曲线

    Figure 9.  Relation between maximum thickness reduction of diaphragm and pressure

    图 10  模型计算所得曲线与实验结果的对比

    Figure 10.  Contrast of deflection and radius between experiment and calculation

    图 11  模型计算所得膜片厚度变化与实验结果的对比

    Figure 11.  Contrast of thickness between experiment and calculation

    表 1  不同膜片组合及厚度的最大破坏压力

    Table 1.  The maximum bearing pressure of different diaphragms thickness and combinations

    膜片使用数量总厚度/mm最大破坏压力/MPa膜片使用数量总厚度/mm最大破坏压力/MPa
    0.05 mm0.1 mm0.2 mm0.05 mm0.1 mm0.2 mm
    1000.050.3114210.63.446
    2000.10.6260030.63.344
    0100.10.5660040.84.411
    4000.21.18140415.568
    0200.21.21402415.514
    0010.21.08900515.524
    4010.42.2444241.26.632
    0210.42.2920061.26.659
    0400.42.417
    下载: 导出CSV

    表 2  三种不同厚度膜片临界点的数据

    Table 2.  Data of critical point in diaphragm with three different thickness

    临界点厚度/mm压力/MPa挠度/mm平均挠度/mmk'/(GPa·m−1)\small$\overline {k{\rm{'}}} $/(GPa·m−1)
    Q10.050.27122.8122.845.425.440
    Q20.10.54522.955.45
    Q30.21.08 22.755.45
    Q1*0.050.27826.0226.475.565.563
    Q2*0.10.55726.015.56
    Q3*0.21.11427.435.57
    下载: 导出CSV
  • [1] 韩肇元, 王震球, 尹协振, 等. 用于研究双波迎面相互作用的电控双驱动激波管实验技术 [J]. 力学学报, 1982, 18(4): 86–92 doi: 10.6052/0459-1879-1982-4-1982-047
    HAN Zhaoyuan, WANG Zhengqiu, YIN Xiezhen, et al. An experimental technique of electrical controlled double driver shock tube for studying head-on moving shock-bow shock interaction [J]. Chinese Journal of Theoretical and Applied Mechanics, 1982, 18(4): 86–92 doi: 10.6052/0459-1879-1982-4-1982-047
    [2] 黄文斌, 谭多望, 闫石, 等. 电破膜式激波管的破膜装置: CN201859050U [P]. 2011-02-16.
    [3] NGUYEN T T N, WILGEROTH J M, PROUD W G. Controlling blast wave generation in a shock tube for biological applications [C]// Journal of Physics: Conference Series, 2014. DOI: 10.1088/1742-6596/500/14/142025.
    [4] GHARABABAEI H, DARVIZEH A. Experimental and analytical investigation of large deformation of thin circular plates subjected to localized and uniform impulsive loading [J]. Mechanics Based Design of Structures & Machines, 2010, 38(2): 171–189. doi: 10.1080/15397730903554633
    [5] BRADLEY J N, LINDSAY R B. Shock waves in chemistry and physics [J]. Physics Today, 1963, 16(8): 60. doi: 10.1063/1.3051077
    [6] BARSOUM I, FALESKOG J. Rupture mechanisms in combined tension and shear—experiments [J]. International Journal of Solids & Structures, 2007, 44(6): 1768–1786. doi: 10.1016/j.ijsolstr.2006.09.031
    [7] ROTHKOPF E M. Diaphragm opening process in shock tubes [J]. Physics of Fluids, 1974, 17(6): 1169–1173. doi: 10.1063/1.1694860
    [8] CAMPBELL G A, KIMBER G M, NAPIER D H. Bursting of diaphragms as related to the operation of shock tubes [J]. Journal of Scientific Instruments, 1965, 42(6): 381. doi: 10.1088/0950-7671/42/6/303
    [9] 佟富强, 刘宝璋. 聚酯(PET)薄膜形变性能的研究 [J]. 包装工程, 1997(1): 26–28
    TONG Fuqiang, LIU Baozhang, ZHANG Jinlong, et al. Study on the deforming properties of the PET film [J]. Packaging Engineering, 1997(1): 26–28
    [10] 黄官强, 鞠开阳, 朱强, 等. PET膜片压力胀型失效试验与有限元分析 [J]. 现代制造工程, 2015(10): 88–92 doi: 10.3969/j.issn.1671-3133.2015.10.020
    HUANG Guanqiang, JU Kaiyang, ZHU Qiang, et al. Experiment and finite element analysis of pressure bulging for PET film [J]. Modern Manufacturing Engineering, 2015(10): 88–92 doi: 10.3969/j.issn.1671-3133.2015.10.020
  • [1] 黄熙龙廖深飞邹立勇刘金宏曹仁义 . 激波与椭圆形重气柱相互作用的PLIF实验. 爆炸与冲击, 2017, 37(5): 829-836. doi: 10.11883/1001-1455(2017)05-0829-08
    [2] 刘金宏邹立勇柏劲松谭多望黄文斌郭文灿 . 激波冲击下air/SF6界面的Richtmyer-Meshkov不稳定性. 爆炸与冲击, 2011, 31(2): 135-140. doi: 10.11883/1001-1455(2011)02-0135-06
    [3] 丁永红尤文斌马铁华 . 舰用动爆冲击波记录系统的设计与应用. 爆炸与冲击, 2013, 33(2): 194-199. doi: 10.11883/1001-1455(2013)02-0194-06
    [4] 彭常贤 , . 声振耦合分析在冲击波测量中的应用. 爆炸与冲击, 2008, 28(5): 427-432. doi: 10.11883/1001-1455(2008)05-0427-06
    [5] 焦清介魏继锋周钢严楠蔡瑞娇 . 爆轰波拐角传播三维数值模拟. 爆炸与冲击, 2003, 23(6): 534-538.
    [6] 王正国孙立英杨志焕冷华光俞鸿儒谷家华李振华李仲发 . 系列生物激波管的研制与应用. 爆炸与冲击, 1993, 13(1): 77-83.
    [7] 易翔宇朱雨建杨基明 . 激波诱导高速气流中液滴的初期变形. 爆炸与冲击, 2017, 37(5): 853-862. doi: 10.11883/1001-1455(2017)05-0853-10
    [8] 何长江于志鲁冯其京 . 高速碰撞的三维欧拉数值模拟方法. 爆炸与冲击, 1999, 19(3): 216-221.
    [9] 初文华朱东俊梁德利封峰韦斯俊 . 基于耦合算法的三维复杂结构冲击动力学特性. 爆炸与冲击, 2018, 38(4): 725-734. doi: 10.11883/bzycj-2016-0283
    [10] 李磊沈兆武李学岭倪小军 . SPH 方法在聚能装药射流三维数值模拟中的应用. 爆炸与冲击, 2012, 32(3): 316-322. doi: 10.11883/1001-1455(2012)03-0316-07
    [11] 唐志平 . 激光辐照下充压柱壳失效的三维离散元模拟. 爆炸与冲击, 2001, 21(1): 1-7.
    [12] 杨宏伟范宝春李鸿志 . 障碍物和管壁导致火焰加速的三维数值模拟. 爆炸与冲击, 2001, 21(4): 259-264.
    [13] 岳鹏涛徐胜利彭金华 . FAE爆炸波对地面目标作用的三维数值研究. 爆炸与冲击, 2000, 20(2): 97-102.
    [14] 邹定祥 . 计算露天矿台阶爆破块度分布的三维数学模型. 爆炸与冲击, 1984, 4(3): 48-59.
    [15] 张伟胡德安韩旭谭柱华 . 陶瓷/金属复合装甲冲击响应的三维SPH 法分析. 爆炸与冲击, 2011, 31(4): 373-379. doi: 10.11883/1001-1455(2011)04-0373-07
    [16] 余育苗王肖均李永池王志海 . 三维正交机织复合材料的弹道性能和破坏模式. 爆炸与冲击, 2009, 29(5): 523-528. doi: 10.11883/1001-1455(2009)05-0523-06
    [17] 张阿漫姚熊亮 . 近边界三维水下爆炸气泡动态特性研究. 爆炸与冲击, 2008, 28(2): 124-130. doi: 10.11883/1001-1455(2008)02-0124-07
    [18] 宋顺成才鸿年王富耻 . 关于弹体冲击和贯穿混凝土的三维数值模拟. 爆炸与冲击, 2006, 26(1): 1-6. doi: 10.11883/1001-1455(2006)01-0001-06
    [19] 桂毓林于川刘仓理孙承纬 . 带尾翼的翻转型爆炸成形弹丸的三维数值模拟. 爆炸与冲击, 2005, 25(4): 313-318. doi: 10.11883/1001-1455(2005)04-0313-06
    [20] 刘凯欣张凯 . 特征线法在三维弹/粘塑性波数值解析中的应用. 爆炸与冲击, 1996, 16(2): 137-142.
  • 加载中
图(11)表(2)
计量
  • 文章访问数:  22
  • HTML全文浏览量:  123
  • PDF下载量:  5
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-11-10
  • 录用日期:  2018-05-07
  • 网络出版日期:  2019-03-25
  • 刊出日期:  2019-03-01

激波管聚酯膜片变形过程分析

    作者简介:何起光(1994- ),男,硕士,478440913@qq.com
    通讯作者: 张伟, zhdawei@hit.edu.cn
  • 1. 哈尔滨工业大学航天学院高速撞击动力学实验室,黑龙江 哈尔滨 150001
  • 2. 北京理工大学前沿交叉科学研究院,北京 100081

摘要: 使用聚酯薄膜作为激波管膜片,通过施加不同压力的激波管实验,获得了膜片厚度及多张膜片的组合方式对膜片所能承受最大压力的影响。利用高速相机对激波管膜片从开始变形到破坏的全过程进行拍摄,使用三维DIC软件获得膜片在变形过程中的位移场。实验发现了膜片会出现圆弧反翘并快速破坏的特别现象,并以此为特征将变形过程分为2个阶段。给出任意厚度膜片第1阶段圆弧变形的数学规律及第2阶段圆弧反翘的形状特征,以及全过程中膜片厚度变化的数学规律。

English Abstract

  • 激波管已广泛用于航空、化学动力学、等离子体物理等领域。近年开始使用聚酯薄膜作为激波管的膜片以替代原来的金属膜片,并出现了多种破膜方法[1-2]。聚酯薄膜作为膜片具有价格低廉、便于安装、易于控制破膜压力等优点。膜片安装在激波管的高压段与低压段之间,当高压段内气体压力超过膜片承载时,膜片将破裂。高压气体向低压段方向传播并产生激波。激波的强度、速度和激波位置与高压段内压力有关,同时与膜片破坏形状和破坏形式有关。所以,获得膜片在高压段内充压过程中的变形规律对激波管激波分析是有意义的。

    Nguyen等[3]研究了膜片承载与膜片厚度的关系。Gharababaei等[4]给出了金属薄板在静态压力下变形的数学规律。Bradley等[5]使用理想的变形模型对激波管中的膜片变形进行过分析。Barsoum等[6]、Rothkopf等[7]和Campbell等[8]研究了膜片破坏时的速度与形状对激波的影响。佟富强等[9]和黄官强等[10]对聚酯薄膜的一维应力应变关系进行过深入研究,得到了单向拉伸的应力应变关系。以上研究都局限在弹性变形范围,假设薄膜为圆弧状变形,而对后续更严重变形没有涉及。

    本文中主要研究不同厚度的聚酯薄膜在不同压力下的变形全过程。发现接近破坏时,膜片会出现圆弧反翘的现象。以此为依据,将全过程分为2个阶段,给出膜片变形过程的数学规律。

    • 图1所示,实验装置由高压段、靶舱和高速摄像机组成。其中,高压段长128 mm,直径为66 mm,最大承受压力为20 MPa。通过压力传感器测量高压段的压力。在高压段的管口,使用夹具装有不同厚度的聚酯薄膜膜片。本实验采用常州新遥光绝缘材料有限公司生产的BE12型号的PET聚酯薄膜作为膜片。如图2所示,采用GB/T1040标准对聚酯薄膜进行常温下准静态拉伸,获得膜片纵向及横向的应力应变曲线。可以看出,薄膜两个方向的弹性模量近似相等,但在塑性阶段,受制作工艺影响,其拉伸强度和硬化模量存在较大差异。

      图  1  实验装置示意图及照片

      Figure 1.  Schematic diagram and photo of experimental equipment

      图  2  膜片的应力应变曲线

      Figure 2.  Stress-strain curves of the diaphragm

      三维DIC技术全称数字图像相关技术,使用两台摄像机对已有散斑的试件进行拍摄后,运用光学原理对某一时刻两个视角的照片进行处理后,可获得试件上各点当前时刻在空间中的坐标。本实验中,使用两台高速摄像机同步记录膜片的破坏过程后,用DIC分析软件ARAMIS对所获得的照片进行计算,获得膜片从充气到破坏的变形历程。在本实验中,压力增速约为12 kPa/s。

    • 在激波管实验中,可采用单膜或多膜结构的聚酯薄膜膜片。膜片的最大承载直接决定激波管高压段内压力大小,而这与聚酯薄膜膜片的厚度有直接关系。本文考虑0.05、0.1和0.2 mm三种不同厚度聚酯薄膜组合,以此获得不同厚度或同厚度下不同组合的薄膜承载能力。

      在上述实验装置中安装不同厚度与不同组合方式的膜片。在高压段中逐渐增加气体压力,由压力传感器得到高压段内的压力。表1给出膜片破坏压力,共有8种不同厚度的不同组合方式。由表1可知,在相同厚度的条件下,不同的膜片组合对膜片最大破坏压力影响甚小。可认为:同种材质、同种实验环境下的聚酯薄膜膜片的最大破坏压力仅与其厚度有关,与组合方式无关。

      膜片使用数量总厚度/mm最大破坏压力/MPa膜片使用数量总厚度/mm最大破坏压力/MPa
      0.05 mm0.1 mm0.2 mm0.05 mm0.1 mm0.2 mm
      1000.050.3114210.63.446
      2000.10.6260030.63.344
      0100.10.5660040.84.411
      4000.21.18140415.568
      0200.21.21402415.514
      0010.21.08900515.524
      4010.42.2444241.26.632
      0210.42.2920061.26.659
      0400.42.417

      表 1  不同膜片组合及厚度的最大破坏压力

      Table 1.  The maximum bearing pressure of different diaphragms thickness and combinations

      用最小二乘法拟合表1的数据,获得最大破坏压力pmax与膜片厚度的关系见图3,图中直线方程为:

      图  3  不同厚度及组合方式的膜片的破坏压力

      Figure 3.  Maximum bearing pressure of diaphragm with different thickness and combination

      $ {p_{\max }} = k\delta $

      式中:斜率k受材料及制作工艺的影响,本实验中k=5.56 GPa/m。各条件下的k可由同样方法得到。

    • 为研究聚酯薄膜膜片的破坏过程,需获得膜片在不同压力下对应的变形。在不同厚度的膜片表面绘制5 mm×5 mm的网格,使用高速相机记录膜片的变形过程。图4 给出了三种厚度膜片在不同压力下的变形过程,可以看出,膜片中心的网格变形最大,四周较小,网格变形基本呈轴对称分布。

      图  4  三种厚度的膜片在不同压力条件下的变形

      Figure 4.  Deformation of diaphragms with three different thicknesses under different pressures

      使用三维DIC技术进一步观察膜片的变形,将两台相机同步记录的图片导入ARAMIS软件进行计算后可得膜片的位移场。图5为0.2 mm厚的膜片在不同压力条件下的变形图。与图4相同,膜片位移场呈轴对称变化,中心处位移最大,因此可由轮廓曲线描述膜片完整形状。

      图  5  DIC处理后所得不同压力条件下0.2 mm厚膜片变形的典型结果

      Figure 5.  Deformation of diaphragm with 0.2 mm in thickness by 3D-DIC analysis

      图6给出了三种厚度的膜片在不同压力下的轮廓曲线。可以看出:随着压力增加,膜片轮廓曲线由圆弧形状逐渐变化为仅头顶部呈圆弧的锥形,最后锥形两侧的直线发生内凹现象,出现圆弧反翘。当轮廓曲线刚出现内凹时,以该锥形作为其极限轮廓,用两条虚线表示,可将变形过程分为2个阶段:第1阶段为前7条曲线,即轮廓曲线在极限轮廓内变化,其中前3条曲线呈圆弧状,各点曲率大致相等且逐渐增大,第4至第7条曲线达到极限轮廓的部分挠度近乎不变,基本呈现直线,而圆弧部分各点曲率仍大致相等且继续增大;第2阶段为后2条曲线,轮廓曲线开始超过极限轮廓,轮廓曲线中央部分出现反翘现象,曲线中点的挠度、曲率迅速增大直至膜片破坏。

      图  6  三种不同厚度膜片的变形轮廓曲线

      Figure 6.  Profiles of deformation of diaphragms with different thicknesses

      对于不同厚度的膜片,图6中虚线的斜率基本相等,平均值为0.73。令轮廓曲线即将超过极限轮廓时对应压力临界点为Qi,膜片破坏时对应压力临界点为Qi*,其中i=1,2,3,分别代表厚度为0.05、0.1和0.2 mm三种膜片,如图7所示。表2给出了这6个点的压力及挠度参数,其中:k′为各点压力与厚度之比,$\overline {k{\rm{'}}} $为不同厚度k′的平均值。在Qi*处,即有k′=k;根据表2,在Qi处,k′=0.974k,此时:

      临界点厚度/mm压力/MPa挠度/mm平均挠度/mmk'/(GPa·m−1)\small$\overline {k{\rm{'}}} $/(GPa·m−1)
      Q10.050.27122.8122.845.425.440
      Q20.10.54522.955.45
      Q30.21.08 22.755.45
      Q1*0.050.27826.0226.475.565.563
      Q2*0.10.55726.015.56
      Q3*0.21.11427.435.57

      表 2  三种不同厚度膜片临界点的数据

      Table 2.  Data of critical point in diaphragm with three different thickness

      图  7  膜片中心挠度与压力曲线

      Figure 7.  Relation between maximum deflection of diaphragm and pressure

      $ \left\{ {\begin{aligned} &{p = 0.974k\delta }\\ &{{Z_{\max }} = ak} \end{aligned}} \right.\quad\quad\quad\quad\quad\quad\quad\;\;(2)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad $

      式中:p为压力;Zmax为膜片中心挠度,即最大挠度;a≈4.100×10-12 m2/Pa。

      膜片的应变场可由位移场计算得到。显然,膜片的应变分布也是轴对称的,只需考虑轮廓曲线上的应变分布。由于Gharababaei等[4]对鼓包后金属薄板的应变计算方法是连续的,但DIC获得的数据是离散点的位移信息,所以对其离散化后得到:

      $ \left\{ \begin{array}{l} {\varepsilon _{r}} = \displaystyle\frac{{\sqrt {{{[{u^{\rm{'}}}\left( {a + 1,b} \right) - u'\left( {a - 1,b} \right)]}^2} + {{[{w^{\rm{'}}}\left( {a + 1,b} \right) - w'\left( {a - 1,b} \right)]}^2}} }}{{u\left( {a + 1,b} \right) - u\left( {a - 1,b} \right)}}\\ {\varepsilon _{\rm{\theta }}} = \displaystyle\frac{{v'\left( {a,b + 1} \right) - v'\left( {a,b - 1} \right)}}{{v\left( {a,b + 1} \right) - v\left( {a,b - 1} \right)}}\\ {\varepsilon _{\rm{t}}} = 1 - \displaystyle\frac{1}{{(1 + {\varepsilon _{{r}}})\left( {1 + {\varepsilon _{\rm{\theta }}}} \right)}} \end{array} \right. $

      式中:ab分别表示网格节点在xy方向的编号;uv分别表示节点在xy方向的初始坐标;u′、v′和w′分别表示在xyz方向变形后的坐标;εrεθ分别为径向应变和环向应变;εt为厚度减少的百分比。

      由式(3) 计算所得膜片厚度变化情况如图8所示,图中曲线与图6中曲线分别对应。膜片大变形时,边缘处与相机的夹角很小,会造成较大的计算误差,所以仅给出膜片中间部分的厚度变化分布。可以看出:处于第1阶段的膜片厚度变化较为均衡,中间处的厚度较边缘减少量略多,厚度略小,分布的形状近似于一段圆弧。随着压力的增加,膜片厚度的变化逐渐加快。在第2阶段中,分布呈锥形,中心处膜片的厚度最小且远小于膜片边缘的。膜片中心处厚度随着压力的增加快速减小,距离中心点越远,厚度的变化速度越慢。中心处也是膜片破裂位置。

      图  8  膜片厚度减小量的分布曲线

      Figure 8.  Distribution curves of thickness reduction

      图9给出三种初始厚度膜片的中心处厚度变化与压力的关系。与中心处挠度与压力的关系类似,三条曲线的趋势是相似的,仅斜率不同。受膜片初始厚度的影响,三条曲线的斜率与膜片的初始厚度也呈反比关系。

      图  9  膜片中心的厚度减小量与压力曲线

      Figure 9.  Relation between maximum thickness reduction of diaphragm and pressure

    • 根据表2Qi点的压力数据可知,对于厚度为δ的膜片,p$\text{≤} $0.974为第1阶段,p$\text{>} $0.974为第2阶段,p=0.974为膜片即将出现反翘现象的压力。

      对于第1阶段,使用最小二乘法对图7进行拟合,并代入激波管内径与膜片初始厚度的参数,在数学上使用两段线拟合第1阶段曲线,可以得到经验公式:

      $ {Z_{{\rm{max}}}} \!\!=\!\! \left\{ {\begin{aligned} &{\left[ {{{0.54p}}/{{(\delta k)}} + 2} \right]{R}/{{33}}\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;p \text{≤} 0.6k\delta }\\& {\left[ {{{\rm{e}}^{{{0.4p/}}{{(\delta k) }}}} /500\;00\!+\!{{0.64p}}/{{(\delta k)}}} \right]{R}/{{33}}\;\;p \text{>} 0.6k\delta } \end{aligned}} \right.\;(4)\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\;\;\;\;\;\; $

      式中:R为激波管半径。

      式(4)将第1阶段的膜片变形分为2个部分:第1部分,p$ \text{≤} $0.6,轮廓曲线近似为一段圆弧;第2部分,p$\text{>} $0.6,膜片轮廓曲线中间近似为一段圆弧,边侧部分近似为直线,其斜率即为极限轮廓的斜率,约为0.73。直线长度(L)、圆弧半径(r)与最大挠度的关系为:

      $ r = \left\{ {\begin{aligned} &{\displaystyle{{{Z_{{\rm{max}}}}}}/{2} + {{{R^2}}}/{{(2{Z_{{\rm{max}}})}}}\;\;\;\;\;\;\;\;\;\;p \text{≤} 0.6k\delta }\\& {3.067R - \displaystyle{{{Z_{{\rm{max}}}}}}/{{0.238}}\;\;\;\;\;\;\;\;p \text{>} 0.6k\delta } \end{aligned}} \right. $

      $ L = 0.59\left( {{Z_{{\rm{max}}}} - r} \right) + 0.81R $

      将式(4)~(6)的计算结果与实验结果对比,如图10所示。可以看出:由式(4)计算所得三种厚度膜片的最大挠度与压力的关系曲线与实验结果对比吻合较好(图10(a));对于圆弧半径(图10(b)),由于式(5)和(6)中使用了式(4)的计算结果,误差累计导致0.2 mm厚度的曲线左端误差较大,但总体而言,计算结果与实验结果吻合较好。由于圆弧的半径和最大挠度决定直线部分长度,后者误差取决于前两者误差,所以直线部分长度误差同样较小。

      图  10  模型计算所得曲线与实验结果的对比

      Figure 10.  Contrast of deflection and radius between experiment and calculation

      膜片厚度的变化直接取决于膜片的变形状态,而图9中曲线与图7类似,可将膜片中心处最大挠度与厚度变化相联系。利用式(4)结果,得到膜片中心处厚度变化与挠度之间的关系:

      $ {\varepsilon _{\rm{t}}} = \left\{ \begin{aligned} &1.415\;7{\left( {{{{Z_{{\rm{max}}}}}}/{R}} \right)^2} - 0.095\;7{{{Z_{{\rm{max}}}}}}/{R} + 0.006\;6\;\;\;\;\;\;p \text{≤} 0.6k\delta \\ &1.056{{{Z_{{\rm{max}}}}}}/{R} - 0.23\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;\;p \text{>} 0.6k\delta \end{aligned} \right. $

      将式(4)代入式(7),可计算膜片中心厚度减少量,如图11所示。

      图  11  模型计算所得膜片厚度变化与实验结果的对比

      Figure 11.  Contrast of thickness between experiment and calculation

    • 本文对激波管聚酯膜片的最大承载进行试验,得到膜片的最大承载与其组合方式无关,仅与其厚度相关且呈线性关系。使用DIC技术对膜片破坏全过程进行观察分析,获得膜片从静止到破坏全过程的变化规律,观察到膜片在破坏前出现圆弧反翘的特别现象。以出现反翘为判据,将整个过程分为2阶段,给出了任意厚度膜片在第一阶段圆弧变形过程中最大挠度与压力及膜片形状的关系,以及膜片在第二阶段圆弧反翘的变形特征。对数据进一步处理获得了膜片在变化过程中的厚度变化,及其与膜片最大挠度的关系。

参考文献 (10)

目录

    /

    返回文章
    返回