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

Hopkinson拉杆平板挂钩试件结构智能协同优化

黄德东 王清华 邢亮亮 徐丰 吴斌

引用本文:
Citation:

Hopkinson拉杆平板挂钩试件结构智能协同优化

    作者简介: 黄德东(1982- ),男,博士,助理研究员,huangdedong@nwpu.edu.cn;
    通讯作者: 徐丰, xufeng@nwpu.edu.cn
  • 中图分类号: O383; TB302.3

Intelligent collaborative optimization of structural parameters for hook-sheet specimens used in split Hopkinson tensile bar

    Corresponding author: XU Feng, xufeng@nwpu.edu.cn ;
  • CLC number: O383; TB302.3

  • 摘要: 与霍普金森拉杆装置中常用的螺纹、胶粘等固定连接方式相比,平板挂钩试件具有连接形式简单、可实现快速组装等优势。针对平板挂钩试件在拉伸过程中因结构几何效应引起的数据测量误差问题,基于影响拉伸试件测量精度的指标:应力平衡达到时间、变形均匀程度、过渡段相对变形以及非轴向力水平,采用正交试验设计、反向传播(back propagation,BP)神经网络与遗传算法相结合的多目标智能协同优化算法对平板挂钩试件的结构参数进行优化,得到了平板挂钩试件最优的结构参数组合,有限元模拟和实验验证了最优结构参数的有效性。该研究结果可为基于平板挂钩试件的霍普金森拉伸实验的数据可靠性分析提供参考。
  • 图 1  分离式霍普金森拉杆装置示意图

    Figure 1.  Schematic of a split Hopkinson tensile bar device

    图 2  与试件连接的杆端的结构及尺寸

    Figure 2.  Structure and dimensions of the tensile bar end connected to the specimen

    图 3  1 500 s−1应变率下AA5182真实应力-真实应变实验曲线

    Figure 3.  True stress-true strain curves of AA5182 at the strain rate of 1 500 s−1

    图 4  试件结构及尺寸

    Figure 4.  Structure and dimensions of the specimen

    图 5  试件及其连接区域的细化网格

    Figure 5.  Refined mesh of the specimen and its connected zone

    图 6  标距段范围内沿试件中心线选取的路径

    Figure 6.  The path taken along the centerline of the specimen in the central section

    图 7  试件的轴向应力偏差随时间变化

    Figure 7.  Relative deviation of axial stress in the specimen with time

    图 8  轴向应力与非轴向应力沿路径的分布

    Figure 8.  Axial stress and non-axial stress along the path

    图 9  非轴向应力/轴向应力比值沿路径的分布

    Figure 9.  Ratio of non-axial stress/axial stress along the path

    图 10  单波加载后试件中的轴向应变分布

    Figure 10.  Distribution of axial strain in specimen after single wave

    图 11  智能协同优化方案流程图

    Figure 11.  The flow chart of intelligent collaborative optimization

    图 12  最小目标值随遗传代数的变化

    Figure 12.  The change of minimum objective function in each generation

    图 13  优化后试件的结构

    Figure 13.  The structure and dimensions of the optimized specimen

    图 14  优化试件的1/4模型网格

    Figure 14.  The 1/4 meshed model of the optimized specimen

    图 15  优化试件轴向应变云图

    Figure 15.  Distribution of axial strain in the optimized specimen

    图 16  加载前的最优结构试件

    Figure 16.  The specimen with optimal structure before loading

    图 17  试件轴向应变分布云图

    Figure 17.  Distributed cloud map of axial strain of the specimen

    图 18  轴向应变分布数值模拟与实验结果的比较

    Figure 18.  Comparison of axial strain distribution between simulation and experiment

    图 19  试件轴向变形分布云图

    Figure 19.  Contour of axial deformation of the specimen

    表 1  试件各段变形量及过渡段相对变形

    Table 1.  The deformation of each section of the specimen and the relative deformation of the transition zone

    前过渡段变形量标距段变形量后过渡段变形量过渡段相对变形
    0.20 mm1.62 mm0.21 mm20.20%
    下载: 导出CSV

    表 2  结构参数正交试验设计表

    Table 2.  Orthogonal test table of structural parameters

    试验编号结构参数(单位:mm)应力平衡达到时间E/μs应变方差V/10−3过渡段相对
    变形D
    非轴向应力
    水平N
    L1W1L2W2RT
    01621230.50.618.000.470 50.051 70.038 1
    026313410.918.000.442 90.136 00.042 1
    036414521.218.610.566 70.249 90.049 7
    04661562.51.518.850.571 30.309 10.056 6
    056816731.820.800.525 60.351 40.081 5
    06721352.51.821.450.139 60.221 90.021 4
    077314630.621.000.195 40.283 10.023 8
    08741570.50.921.000.879 60.014 90.075 0
    097616311.221.500.950 60.114 00.082 9
    107812421.522.000.805 50.230 20.081 0
    118214711.522.500.228 90.074 00.032 9
    128315321.823.000.235 80.196 70.030 0
    13841642.50.622.790.281 70.244 50.032 9
    148612530.922.810.371 70.287 40.046 2
    15881360.51.222.811.593 20.074 60.125 3
    169215431.224.000.085 90.196 90.011 0
    17931650.51.523.500.774 20.034 60.048 7
    189412611.824.000.317 80.061 00.044 9
    199613720.623.410.472 10.194 60.050 4
    20981432.50.924.000.607 00.230 40.060 7
    2110216620.925.200.154 70.138 70.013 5
    221031272.51.225.200.192 60.191 70.019 4
    2310413331.525.000.211 60.224 40.022 7
    241061440.51.825.011.442 40.027 00.086 6
    2510815510.624.001.117 40.096 90.088 8
    下载: 导出CSV

    表 3  测试样本数据

    Table 3.  The data of test samples

    测试试验结构参数(单位:mm)应力平衡达到时间E/μs应变方差V/10−3过渡段相对
    变形D
    非轴向应力
    水平N
    L1W1L2W2RT
    016314611.518.770.285 40.119 20.047 7
    028212521.222.680.203 50.162 40.018 7
    031041672.50.925.120.276 90.205 30.029 0
    下载: 导出CSV

    表 4  各项指标及目标值网络预测与实际情况的比较

    Table 4.  The comparison of predicted and actual values of indicators and objective function

    指标/目标值123
    实际值预测值误差/%实际值预测值误差/%实际值预测值误差/%
    应力平达到衡时间 E/μs18.7717.89 4.722.6820.1511.225.1227.459.2
    应变方差 V/10−30.285 40.263 3 7.70.203 50.185 7 8.70.276 90.299 48.1
    过渡段相对变形 D0.119 20.104 012.80.162 40.163 8 0.90.205 30.210 72.6
    非轴向应力水平 N0.047 70.050 1 5.00.018 70.020 811.20.029 00.029 51.7
    目标函数值 Obj3.009 92.887 7 4.12.716 72.620 6 3.53.417 33.604 65.5
    下载: 导出CSV

    表 5  优化试件各段变形量及过渡段相对变形

    Table 5.  The deformation of each section of the optimized specimen and the relative deformation of the transition zone

    前过渡段变形量标距段变形量后过渡段变形量过渡段相对变形
    0.16 mm2.03 mm0.16 mm13.62%
    下载: 导出CSV

    表 6  优化前后各项指标的比较

    Table 6.  The comparison of various indicators before and after the optimization

    指标优化前优化后增大(↑)或降低(↓)
    应力平衡达到时间E/μs22.8124.758.51%↑
    应变方差V/10−30.598 80.135 177.44%↓
    过渡段相对变形 D0.202 00.136 232.57%↓
    非轴向应力水平 N0.056 30.018 866.61%↓
    下载: 导出CSV

    表 7  各段变形值实验与计算的比较

    Table 7.  Comparison of deformation of transition zones between simulations and experiments

    变形段计算值实验值相对偏差
    前过渡段0.16 mm0.13 mm18.75%
    后过渡段0.16 mm0.15 mm6.25%
    标距段2.03 mm2.12 mm4.43%
    相对变形13.62%11.67%14.32%
    下载: 导出CSV

    表 8  参考试件与优化试件数值模拟和实验结果的对比

    Table 8.  The comparison of simulated and experimental results between reference specimen and optimized specimen

    指标参考试件优化试件
    计算实验相对偏差计算实验相对偏差
    过渡段相对变形/%20.2022.3110.45%13.6211.6714.32%
    标距段变形均匀度/10−30.598 80.544 69.05%0.135 10.145 57.70%
    下载: 导出CSV
  • [1] KOLSKY H. An investigation of the mechanical properties of material at very high rates of loading [J]. Proceedings of the Physical Society: Section B, 1949, 62(1): 676–700.
    [2] HARDING J, WOOD E O, CAMPBELL J D. Tensile testing of materials at impact rates of strain [J]. Journal of Mechanical Engineering Science, 1960, 2(2): 88–96. DOI: 10.1243/JMES_JOUR_1960_002_016_02.
    [3] OGAWA K. Impact-tension compression test by using a split Hopkinson bar [J]. Experimental Mechanics, 1984, 24(2): 81–86. DOI: 10.1007/BF02324987.
    [4] 申海艇, 蒋招绣, 王贝壳, 等. 基于超高速相机的数字图像相关性全场应变分析在SHTB实验中的应用 [J]. 爆炸与冲击, 2017, 37(1): 15–20. DOI: 10.11883/1001-1455(2017)01-0015-06.
    SHEN Haiting, JIANG Zhaoxiu, WANG Beike, et al. Full field strain measurement in split Hopkinson tension bar experiments by using ultra-high-speed camera with digital image correlation [J]. Explosion and Shock Waves, 2017, 37(1): 15–20. DOI: 10.11883/1001-1455(2017)01-0015-06.
    [5] GILAT A, GOLDBERG R K, ROBERTS G D. Experimental study of strain-rate-dependent behavior of carbon/epoxy composite [J]. Composites Science and Technology, 2002, 62(10/11): 1469–1479. DOI: 10.1016/S0266-3538(02)00100-8.
    [6] CHEN X, LI Y, ZHI Z, et al. The compressive and tensile behavior of a 0/90 C fiber woven composite at high strain rates [J]. Carbon, 2013, 61: 97–104. DOI: 10.1016/j.carbon.2013.04.073.
    [7] OWENS A T, TIPPUR H V. A tensile split Hopkinson bar for testing particulate polymer composites under elevated rates of loading [J]. Experimental Mechanics, 2008, 49(6): 799–811. DOI: 10.1007/s11340-008-9192-7.
    [8] 朱耀, 庞宝君, 盖秉政. 一种用于动态拉伸试验装置的新型试件装卡方式 [J]. 实验力学, 2009, 24(5): 433–438. DOI: 1001-4888(2009)05-0433-06.
    ZHU Yao, PANG Baojun, GAI Bingzheng. A new specimen fastener for dynamic tensile testing apparatus [J]. Journal of Experimental Mechanics, 2009, 24(5): 433–438. DOI: 1001-4888(2009)05-0433-06.
    [9] TAN X, GUO W, GAO X, et al. A new technique for conducting split Hopkinson tensile bar test at elevated temperatures [J]. Experimental Techniques, 2017, 41(2): 191–201. DOI: 10.1007/s40799-017-0167-4.
    [10] VERLEYSEN P, DEGRIECK J. Non-homogeneous and multi-axial stress distribution in concrete specimens during split Hopkinson tensile tests [J]. Computers and Structures, 2000, 77(6): 669–676. DOI: 10.1016/S0045-7949(00)00022-5.
    [11] CHEN W, LU F, CHENG M. Tension and compression tests of two polymers under quasi-static and dynamic loading [J]. Polymer Testing, 2002, 21(2): 113–121. DOI: 10.1016/S0142-9418(01)00055-1.
    [12] PRABOWO Q A, KARIEM M A, GUNAWAN L. The effect of specimen dimension on the results of the split Hopkinson tension bar testing [J]. Procedia Engineering, 2017, 173: 608–614. DOI: 10.1016/j.proeng.2016.12.114.
    [13] VERLEYSEN P, DEGRIECK J, VERSTRAETE T, et al. Influence of specimen geometry on split Hopkinson tensile bar test on sheet materials [J]. Experimental Mechanics, 2008, 48(5): 587–598. DOI: 10.1007/s11340-008-9149-x.
    [14] NGUYEN K H, KIM H C, SHIN H, et al. Numerical investigation into the stress wave transmitting characteristics of threads in the split Hopkinson tensile bar test [J]. International Journal of Impact Engineering, 2017, 109: 253–263. DOI: 10.1016/j.ijimpeng.2017.07.004.
    [15] 王礼立. 应力波基础[M]. 第2版. 北京: 国防工业出版社, 2005: 52−60.
    [16] SMERD R, WINKLER S, SALISBURY C, et al. High strain rate tensile testing of automotive aluminum alloy sheet [J]. International Journal of Impact Engineering, 2005, 32(1): 541–560. DOI: 10.1016/j.ijimpeng.2005.04.013.
    [17] SONG B, CHEN W. Dynamic stress equilibrium in split Hopkinson pressure bar tests on soft materials [J]. Experimental Mechanics, 2004, 44: 300–312. DOI: 10.1007/BF02427897.
    [18] 施建俊, 李庆亚, 张琪, 等. 基于Matlab和BP神经网络的爆破振动预测系统 [J]. 爆炸与冲击, 2017, 37(6): 1087–1092. DOI: 10.11883/1001-1455(2017)06-1087-06.
    SHI Jianjun, LI Qingya, ZHANG Qi. Forecast system for blasting vibration velocity peak based on Matlab and BP neural network [J]. Explosion and Shock Waves, 2017, 37(6): 1087–1092. DOI: 10.11883/1001-1455(2017)06-1087-06.
    [19] WEN J, LIU C H, YAO H, et al. A nonlinear dynamic model and parameters identification method for predicting the shock pulse of rubber waveform generator [J]. International Journal of Impact Engineering, 2018, 120: 1–15. DOI: 10.1016/j.ijimpeng.2018.05.009.
    [20] 李守巨, 刘迎曦, 刘玉静, 等. 基于遗传算法的爆炸冲击荷载参数识别方法 [J]. 爆炸与冲击, 2002, 22(4): 295–300.
    LI Shouju, LIU Yingxi, LIU Yujing, et al. Parameter identification procedures of explosion shock loading based on genetic algorithm [J]. Explosion and Shock Waves, 2002, 22(4): 295–300.
  • [1] 朱秀云潘蓉林皋李亮 . 基于ANSYS/LS-DYNA的钢板混凝土墙冲击实验的有限元分析. 爆炸与冲击, 2015, 35(2): 222-228. doi: 10.11883/1001-1455(2015)02-0222-07
    [2] 李志强刘晓明赵勇刚赵隆茂 . 微爆索线性切割航空有机玻璃的实验研究和有限元分析. 爆炸与冲击, 2007, 27(5): 385-389. doi: 10.11883/1001-1455(2007)05-0385-05
    [3] 莫根林金永喜李忠新吴志林 . 多目标优化算法在弹头侵彻明胶运动模拟中的应用. 爆炸与冲击, 2018, 38(6): 1404-1411. doi: 10.11883/bzycj-2017-0160
    [4] 曹丽娜韩秀清张永光曹宇欣 . 射孔弹爆轰波间干扰的有限元分析. 爆炸与冲击, 2011, 31(2): 220-224. doi: 10.11883/1001-1455(2011)02-0220-05
    [5] 张永康李玉龙汪海青 . 典型梁-缘结构鸟撞破坏的有限元分析. 爆炸与冲击, 2008, 28(3): 236-242. doi: 10.11883/1001-1455(2008)03-0236-07
    [6] 侯海量朱锡古美邦 . 爆炸载荷作用下加筋板的失效模式分析及结构优化设计. 爆炸与冲击, 2007, 27(1): 26-33. doi: 10.11883/1001-1455(2007)01-0026-08
    [7] 米双山张锡恩陶贵明 . 钨球侵彻LY-12铝合金靶板的有限元分析. 爆炸与冲击, 2005, 25(5): 477-480. doi: 10.11883/1001-1455(2005)05-0477-04
    [8] 李典朱锡侯海量仲强 . 高速杆式弹体侵彻下蓄液结构载荷特性的有限元分析. 爆炸与冲击, 2016, 36(1): 1-8. doi: 10.11883/1001-1455(2016)01-0001-08
    [9] 宋延泽李志强赵隆茂 . 基于十四面体模型的闭孔泡沫材料动态力学性能的有限元分析. 爆炸与冲击, 2009, 29(1): 49-55. doi: 10.11883/1001-1455(2009)01-0049-07
    [10] 张磊胡时胜 . 混凝土层裂强度测量的新方法. 爆炸与冲击, 2006, 26(6): 537-542. doi: 10.11883/1001-1455(2006)06-0537-06
    [11] 毋玲郭英男李玉龙 . 蜂窝夹芯雷达罩结构的鸟撞数值分析. 爆炸与冲击, 2009, 29(6): 642-647. doi: 10.11883/1001-1455(2009)06-0642-06
    [12] 侯海量朱锡梅志远 . 舱内爆炸载荷及舱室板架结构的失效模式分析. 爆炸与冲击, 2007, 27(2): 151-158. doi: 10.11883/1001-1455(2007)02-0151-08
    [13] 刘烨王蕊李志刚 . CFRP-混凝土-钢管组合结构在低速侧向撞击下的动力响应. 爆炸与冲击, 2018, 38(4): 759-767. doi: 10.11883/bzycj-2016-0349
    [14] 徐挺缪馥星周风华杨黎明 . 多网格防阻块护栏系统的冲击响应. 爆炸与冲击, 2018, 38(4): 820-826. doi: 10.11883/bzycj-2016-0354
    [15] 魏晓凯李杰郑涛张樨冯凯强钱海宁 . 被动式半捷联稳定平台抗高过载结构设计. 爆炸与冲击, 2019, 39(7): 075102-1-075102-10. doi: 10.11883/bzycj-2018-0064
    [16] 钱七虎黄小平唐德高李志成 . 冲击荷载作用下有限元方法求解波动过程的精度研究. 爆炸与冲击, 1995, 15(1): 1-10.
    [17] 任波王乘 . 并行冲击/侵彻有限元数值模拟技术. 爆炸与冲击, 2005, 25(3): 260-264. doi: 10.11883/1001-1455(2005)03-0260-05
    [18] 李妍妍郑志军虞吉林王长峰 . 闭孔泡沫金属变形模式的有限元分析. 爆炸与冲击, 2014, 34(4): 464-470. doi: 10.11883/1001-1455(2014)04-0464-07
    [19] 卫丰张光 . 药筒发射应力和抽壳力的有限元分析. 爆炸与冲击, 2001, 21(1): 72-75.
    [20] 谢宗蕻卞文杰昂海松乔新 . 峰窝夹芯结构雷达罩鸟撞有限元分析与模拟. 爆炸与冲击, 1999, 19(3): 235-242.
  • 加载中
图(19)表(8)
计量
  • 文章访问数:  51
  • HTML全文浏览量:  75
  • PDF下载量:  4
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-09-26
  • 录用日期:  2019-04-01
  • 网络出版日期:  2019-09-25
  • 刊出日期:  2019-10-01

Hopkinson拉杆平板挂钩试件结构智能协同优化

    作者简介:黄德东(1982- ),男,博士,助理研究员,huangdedong@nwpu.edu.cn
    通讯作者: 徐丰, xufeng@nwpu.edu.cn
  • 1. 西北工业大学航天学院,陕西 西安 710072
  • 2. 北京电子工程总体研究所,北京 100854

摘要: 与霍普金森拉杆装置中常用的螺纹、胶粘等固定连接方式相比,平板挂钩试件具有连接形式简单、可实现快速组装等优势。针对平板挂钩试件在拉伸过程中因结构几何效应引起的数据测量误差问题,基于影响拉伸试件测量精度的指标:应力平衡达到时间、变形均匀程度、过渡段相对变形以及非轴向力水平,采用正交试验设计、反向传播(back propagation,BP)神经网络与遗传算法相结合的多目标智能协同优化算法对平板挂钩试件的结构参数进行优化,得到了平板挂钩试件最优的结构参数组合,有限元模拟和实验验证了最优结构参数的有效性。该研究结果可为基于平板挂钩试件的霍普金森拉伸实验的数据可靠性分析提供参考。

English Abstract

  • 材料在高应变率下的力学性能是材料科学领域研究的重要课题之一。Kolsky[1]改进的分离式Hopkinson压杆(split Hopkinson pressure bar,SHPB)自问世以来成为研究材料在高应变率下压缩性能的重要装置。Harding等[2]在压杆的基础上开展了分离式Hopkinson拉杆(split Hopkinson tensile bar,SHTB)的研究。随后,Ogawa[3]改进后的拉杆装置成为现在常见分离式Hopkinson拉杆的原型。

    试件与拉杆的连接问题是SHTB实验技术的关键问题之一,应用较广泛的连接方式有螺纹连接和片状试件的胶粘连接[3-6]。但以上两种连接方式均存在突出缺点:螺纹连接试件加工费料费时,且对精度要求较高[2];胶粘连接试件在每次实验前需将试件与杆端粘连、固化[5-6],实验效率低,易对杆端造成累积损伤等。鉴于此,Owens等[7]最早提出了平板挂钩连接方式,该方式具有形式简单、可实现快速组装等优点。朱耀等[8]通过数值模拟和实验证明了平板挂钩式连接的可行性,Tan等[9]将这一连接方式应用于分离式Hopkinson拉杆高温同步组装实验技术中。

    分离式Hopkinson拉杆实验中,影响测量精度的试件几何效应一直以来备受关注[3,10-12]。近年来,Verleysen等[13]通过实验证明了胶粘连接方式片状试件的过渡段半径和标距段长宽比对实验的测量精度存在明显的影响。Nguyen等[14]采用有限元模拟研究了螺纹连接方式中试件螺牙形状、螺纹内外径、螺距等结构参数对实验测量精度的影响。为了进一步提高SHTB实验试样测量的精确度,Verleysen等[13]、Nguyen等[14]还分别对胶粘连接和螺纹连接试样的结构参数进行了优化。针对SHTB实验中连接形式简单、可快速组装的平板挂钩试样,如何实现其结构最优化以提高实验的测量精度鲜有人研究[7-14];Verleysen等[13]、Nguyen等[14]所采用的试样结构优化方法均为控制变量式的研究方法,即通过改变参考组单一变量,控制其余参数不变进行对比寻优,该方法仅能获得局部次优解而不能进行参数取值范围内的全局寻优,为此,本文在分析平板挂钩连接试样结构参数对实验精度影响的基础上,采用有限元数值模拟与智能协同优化算法相结合的方法对平板挂钩试件的结构进行全局优化,建立了结构参数与测量精度指标之间的非线性映射关系,得出最优的结构设计方案,并通过有限元模拟和实验进行了验证。

    • 为分离式Hopkinson拉杆实验装置,如图1所示。实验时,撞击杆(子弹)由气室内的高压气体驱动,沿炮管向左撞击位于入射杆末端的凸台。撞击产生的拉伸波沿入射杆向右传播,当到达入射杆与试件的连接处时,一部分将以压缩波的形式反射回入射杆,另一部分仍以拉伸波的形式经过试件传入透射杆。分别通过入射杆与透射杆上的应变片采集反射应变εR和透射应变εT,则利用公式(1)可得到试件的应力、应变以及应变率[15]

      图  1  分离式霍普金森拉杆装置示意图

      Figure 1.  Schematic of a split Hopkinson tensile bar device

      $\left\{ {\begin{aligned} &{{\sigma _{\rm{s}}} = \frac{{AE}}{{{A_{\rm{s}}}}}{\varepsilon _{\rm{T}}}} \\ & {{\varepsilon _{\rm{s}}} = - \frac{{2{C_0}}}{{{L_{\rm{s}}}}}\int_0^t {{\varepsilon _{\rm{R}}}{\rm{d}}\tau } } \\ &{{{\dot \varepsilon }_{\rm{s}}} = - \frac{{2{C_{\rm{0}}}}}{{{L_{\rm{s}}}}}{\varepsilon _{\rm{R}}}} \end{aligned}} \right.$

      式中:AEC0分别为杆的截面积、弹性模量和弹性波速,AsLs分别为试件测量段的截面积与长度,${\sigma _{\rm{s}}}$${\varepsilon _{\rm{s}}}$${\dot \varepsilon _{\rm{s}}}$别为试件应力、应变和应变率。

      需要说明的是,式(1)是基于一维应力假设且以试件受力、变形均匀,试件前后两端处于应力平衡状态为前提推导得出的。利用式(1)进行材料力学性能分析时,测量精度须满足以上假设和前提。

    • 平板挂钩连接结构中与试件相连接的杆端的结构和参考尺寸[9](单位:mm)如图2所示,试件的结构和参考尺寸[9](单位:mm)如图3所示。由图3可以看出平板挂钩试件的结构主要分三部分:前后过渡段①、③,标距段②以及实现试件与杆端相连接的长方体挂钩。本文所关注的试件关键结构参数在图3中亦有所标识:标距段长L1、标距段宽W1、过渡段半径R、试件厚度T、挂钩与杆接触面长L2、挂钩与杆接触面宽W2

      图  2  与试件连接的杆端的结构及尺寸

      Figure 2.  Structure and dimensions of the tensile bar end connected to the specimen

      图  3  1 500 s−1应变率下AA5182真实应力-真实应变实验曲线

      Figure 3.  True stress-true strain curves of AA5182 at the strain rate of 1 500 s−1

    • 进行数值模拟之前,先采用3D建模软件SolidWorks建立SHTB实验装置的模型,由于SHTB装置整体上具有很好的轴对称性,为减少计算工作量,建模时简化为1/4模型。所模拟装置的几何尺寸如下:入射杆、透射杆直径均为19 mm,长度分别为2 800 mm、1 400 mm;子弹长500 mm,内径、外径分别为21 mm、28 mm,凸台厚10 mm,外径为28 mm;试件及其连接区域的尺寸如图2所示。将建立好的模型导入有限元软件ABAQUS(Explicit)模拟分离式Hopkinson拉伸实验,为撞击杆定义初始速度场代替气炮发射实现撞击凸台的过程。

    • 入射杆、透射杆以及子弹的材料采用45钢,材料模型为线弹性模型,密度为7 800 kg/m3,弹性模量为210 GPa,泊松比为0.28。试件材料采用AA5182型铝合金,材料模型为弹塑性模型,密度为2 700 kg/m3,弹性模量为70 GPa,泊松比为0.3;塑性段采用AA5182合金在应变率为1 500 s−1时由实验测得的真实应力-真实应变曲线[16],如图3中黑色曲线所示,红色标记点则表示输入ABAQUS(Explicit)用于有限元计算的数据点。为验证采用该材料模型的可行性与正确性,选择具有图4所示结构和尺寸的试件进行计算验证,得到平均应变率为1 500 s−1时的计算结果如图3蓝色曲线所示,可以看出有限元模拟的结果曲线虽有所波动但与实验曲线基本吻合,说明采用该材料模型具有正确性。

      图  4  试件结构及尺寸

      Figure 4.  Structure and dimensions of the specimen

      本文关注的重点是试件及其连接区域的应力应变状态,因此需对试件及其连接区域进行网格细化以提高计算精度,图5所示即为细化后的网格。杆端利用Tet单元实现从粗网格到细网格的过渡,单元类型为C3D4;其余网格均采用Hex单元,单元类型为C3D8R。入射杆、透射杆、子弹以及试件的单元数量分别为:6 949、4 674、715和1 955;其中,杆端过渡段区域的Tet单元的数量为889。

      图  5  试件及其连接区域的细化网格

      Figure 5.  Refined mesh of the specimen and its connected zone

    • 根据2.1节中式(1)成立的假设和前提,为提高测量精度,试件标距段需满足非轴向应力尽量小[10],达到应力平衡的时间尽量短[11],变形尽可能均匀[13]。此外,由于过渡段变形将导致实测应变偏大,应使过渡段变形尽量小[13]

    • 定义试件的应力平衡达到时间为应力波到达试件标距段末端(轴向应力开始大于1 MPa)时起,直至试件标距段前后两端轴向应力的相对偏差开始小于5%为试件达到应力平衡所需要的时间[17],并将其作为衡量试件达到应力平衡难易程度的指标,时间越长试件越不易达到平衡。图6所示是在标距段两端选取的前后端点,图7所示为图6前后端点轴向应力偏差随时间的变化。

      图  6  标距段范围内沿试件中心线选取的路径

      Figure 6.  The path taken along the centerline of the specimen in the central section

      图  7  试件的轴向应力偏差随时间变化

      Figure 7.  Relative deviation of axial stress in the specimen with time

      图7可以看出,轴向应力偏差随时间逐渐减小,当应力波到达试件标距段末端时应力偏差为26.35%、22.81 μs后降低至4.86%,应力平衡达到时间为22.81 μs。

    • 试件中点的轴向应力达到最大值时计算W1(见图3)方向上的应力与轴向应力的比值,将路径上比值的均值作为衡量非轴向应力水平的指标。图8所示为轴向应力与非轴向应力沿路径的分布,非轴向应力相对轴向应力明显较小且分布的均匀性较差,沿路径,非轴向应力中间区域相对两侧较小。图9为非轴向应力与轴向应力的比值沿路径的分布,可以看出分布并不均匀,中间区域的比值相对两侧较小。整个路径上应力比值的均值即非轴向应力水平为0.056 3。

      图  8  轴向应力与非轴向应力沿路径的分布

      Figure 8.  Axial stress and non-axial stress along the path

      图  9  非轴向应力/轴向应力比值沿路径的分布

      Figure 9.  Ratio of non-axial stress/axial stress along the path

    • 拉伸波对试件加载完成后试件轴向应变的分布如图10所示。可以看出,试件中点附近的应变最大,往两侧应变逐渐减小,体现了参考试件标距段变形的不均匀性。为进一步衡量试件标距段变形的不均匀程度,以路径上应变值的方差为指标,方差越大则变形越不均匀。图10所示应变值的方差为0.598 8×10−3

      图  10  单波加载后试件中的轴向应变分布

      Figure 10.  Distribution of axial strain in specimen after single wave

    • 图10可以看出,标距段发生变形的过程中,过渡段也发生了相对较小的变形,可以用前、后过渡段变形量占全段(过渡段和标距段)变形量的比重来衡量过渡段相对变形的大小,结果如表1所示。

      前过渡段变形量标距段变形量后过渡段变形量过渡段相对变形
      0.20 mm1.62 mm0.21 mm20.20%

      表 1  试件各段变形量及过渡段相对变形

      Table 1.  The deformation of each section of the specimen and the relative deformation of the transition zone

      从以上结果可以看出,式(1)成立的假设和前提在实际情况下都得不到满足,但可以根据以上结果提出衡量试件测量精度的指标:应力平衡达到时间、非轴向应力水平、变形均匀程度以及过渡段相对变形。很明显,应力平衡达到时间越短,非轴向应力水平越低,变形越均匀,过渡段变形越小则试件的测量精度越高。因此,兼顾以上4个指标使优化后的平板挂钩试件的指标水平总体处于最低是进行试件结构参数优化的目的所在。

    • 智能协同优化的目的旨在通过正交试验设计、BP神经网络以及遗传算法的相互补充,克服控制变量法无法进行全局寻优的缺陷,建立智能、准确的协同优化方案,获得分离式霍普金森拉杆平板挂钩试件结构参数在设定取值区间内的最优组合。

    • BP神经网络是一种信号向前传递、误差向后传播的多层前馈神经网络,具有强大的非线性建模能力和预测能力[18-19] 。正交试验设计所选取的数据具有均匀性、代表性,可以以最少的试验次数达到与全面试验等效的效果。因此,可通过正交试验设计建立样本,BP神经网络学习样本,建立结构参数与精确测量各指标间的非线性映射关系。

      遗传算法(genetic algorithm, GA)是一种具备全局寻优能力的优化算法[20] 。可将BP神经网络与GA相结合,通过BP神经网络预测个体目标值从而确定适应度值,GA则根据个体适应度值进行选择、进化得到最优个体即最优结构参数组合。智能协同优化方案流程如图11所示。

      图  11  智能协同优化方案流程图

      Figure 11.  The flow chart of intelligent collaborative optimization

    • 以分离式Hopkinson拉杆平板挂钩试件的关键结构参数L1W1L2W2RT为研究对象,设计六因素五水平正交试验建立BP神经网络的样本数据库,如表2所示:第1列为试验编号,六因素五水平正交试验仅需25组试验;第2~7列为平板挂钩试件结构参数;第8~11列为通过有限元模拟获得的试件测量精度的各指标,此处需要说明的是,由于数值模拟采用的材料模型应变率为1 500 s−1,为提升数值模拟的准确度,各组正交试验的仿真均是在调整撞击杆初始撞击速度使试件平均应变率在1 500 s−1左右(1 450 s−1~1 550 s−1范围内)的条件下进行的。

      试验编号结构参数(单位:mm)应力平衡达到时间E/μs应变方差V/10−3过渡段相对
      变形D
      非轴向应力
      水平N
      L1W1L2W2RT
      01621230.50.618.000.470 50.051 70.038 1
      026313410.918.000.442 90.136 00.042 1
      036414521.218.610.566 70.249 90.049 7
      04661562.51.518.850.571 30.309 10.056 6
      056816731.820.800.525 60.351 40.081 5
      06721352.51.821.450.139 60.221 90.021 4
      077314630.621.000.195 40.283 10.023 8
      08741570.50.921.000.879 60.014 90.075 0
      097616311.221.500.950 60.114 00.082 9
      107812421.522.000.805 50.230 20.081 0
      118214711.522.500.228 90.074 00.032 9
      128315321.823.000.235 80.196 70.030 0
      13841642.50.622.790.281 70.244 50.032 9
      148612530.922.810.371 70.287 40.046 2
      15881360.51.222.811.593 20.074 60.125 3
      169215431.224.000.085 90.196 90.011 0
      17931650.51.523.500.774 20.034 60.048 7
      189412611.824.000.317 80.061 00.044 9
      199613720.623.410.472 10.194 60.050 4
      20981432.50.924.000.607 00.230 40.060 7
      2110216620.925.200.154 70.138 70.013 5
      221031272.51.225.200.192 60.191 70.019 4
      2310413331.525.000.211 60.224 40.022 7
      241061440.51.825.011.442 40.027 00.086 6
      2510815510.624.001.117 40.096 90.088 8

      表 2  结构参数正交试验设计表

      Table 2.  Orthogonal test table of structural parameters

    • 采用单隐含层BP神经网络依据样本数据库进行训练,单隐含层BP神经网络总共分为三层:输入层、隐含层和输出层。根据需要拟合的非线性函数特点确定网络结构为6-25-4,即:输入层有6个节点,分别对应6个结构参数;输出层有4个节点,分别对应4个指标;隐含层有25个节点,连接输入层与输出层。设置隐含层间的传递函数为正切S型传递函数tansig,输出层间的传递函数为线性传递函数purelin,训练函数为基于Levenberg-Marquardt算法的trainlm。设置训练次数为100,学习速率为0.001,学习目标为1×10−5

      设输入矩阵X为结构参数矩阵,输出矩阵Y为精度指标矩阵:

      ${{ X}^{\rm{T}}} = \left[ {{{{{L}}}_1}},\;{{{{W}}_1}},\;{{{{L}}_2}}, \;{{{{W}}_2}},\;{{R}},\;{{T}} \right]$

      ${Y^{\rm{T}}} = \left[ {{E}},\;{{V}},\;{{D}},\;{{N}} \right]$

      式中:EVDN分别表示试件测量精度指标应力平衡达到时间、应变方差、过渡段相对变形以及非轴向应力水平,则在训练结束时,合格的神经网络将学会输入矩阵X与输出矩阵Y之间蕴含的内在关系,具备从输入矩阵X到输出矩阵Y之间的非线性映射能力及预测能力。

    • 利用BP神经网络拟合出的XY之间的非线性映射,兼顾3.3中的各指标,构建体现指标整体水平的目标函数:

      ${\rm{Obj}} = \frac{e}{{\sum\limits_{i = 1}^k {{E_i}\bigg/k} }} + \frac{v}{{\sum\limits_{i = 1}^k {{V_i}\bigg/k} }} + \frac{d}{{\sum\limits_{i = 1}^k {{D_i}\bigg/k} }} + \frac{n}{{\sum\limits_{i = 1}^k {{N_i}\bigg/k} }}$

      式中:(e, v, d, n)为任意结构参数组合(l1, w1, l2, w2, r, t)在BP神经网络映射下输出的各项指标值,EiViDiNi均取自样本库表2的后4列,k=25,Obj为目标函数值,不同指标对目标函数值的影响权重相同。

      利用MTALAB谢菲尔德遗传算法工具箱进行遗传算法的编程,采用rank函数根据目标值为个体分配适应度值,目标值越大所得适应度值越小,越容易被淘汰。

    • BP神经网络的学习或训练过程实际上是确定各层之间权值与阈值的过程,即使学习同样的样本这一过程也存在一定的随机性,因此在利用神经网络建立的模型确定个体目标值之前,需先验证网络模型的有效性和准确性。为此,设置表3所示的用于网络模型检测的测试样本数据。

      测试试验结构参数(单位:mm)应力平衡达到时间E/μs应变方差V/10−3过渡段相对
      变形D
      非轴向应力
      水平N
      L1W1L2W2RT
      016314611.518.770.285 40.119 20.047 7
      028212521.222.680.203 50.162 40.018 7
      031041672.50.925.120.276 90.205 30.029 0

      表 3  测试样本数据

      Table 3.  The data of test samples

      以测试样本中的结构参数作为输入矩阵Xt(6×3)输入到BP神经网络,训练好的BP神经网络会将矩阵Xt映射到输出矩阵Yt。比较输出矩阵Yt中的预测值与表3中有限元模拟获得的实际值即可验证网络模型的有效性和准确性,结果如表4所示。由表4可知三组测试试验各项指标的预测值与实际值之间最大偏差为12.8%、平均偏差为7.0%;目标值的预测值与实际值之间最大偏差为5.5%,平均偏差为4.4%。可以认为所训练的神经网络为合格的网络。

      指标/目标值123
      实际值预测值误差/%实际值预测值误差/%实际值预测值误差/%
      应力平达到衡时间 E/μs18.7717.89 4.722.6820.1511.225.1227.459.2
      应变方差 V/10−30.285 40.263 3 7.70.203 50.185 7 8.70.276 90.299 48.1
      过渡段相对变形 D0.119 20.104 012.80.162 40.163 8 0.90.205 30.210 72.6
      非轴向应力水平 N0.047 70.050 1 5.00.018 70.020 811.20.029 00.029 51.7
      目标函数值 Obj3.009 92.887 7 4.12.716 72.620 6 3.53.417 33.604 65.5

      表 4  各项指标及目标值网络预测与实际情况的比较

      Table 4.  The comparison of predicted and actual values of indicators and objective function

    • 表3设置遗传算法的各项参数,基于式(4)所示目标函数采用4.1中训练好的网络确定个体适应度,根据适应度选择或淘汰个体进行寻优。种群共进化50代,记录进化过程中每一代的最优个体(最小目标值)如图12所示。可以看出随着进化的进行更加优秀的个体不断被选择,直至20代之后进化过程趋于稳定,而40代之后每代中最优个体的目标值基本不再改变;因此可认为种群在末代进化出最优个体,利用神经网络协同的遗传算法寻优成功。基于给定目标函数(式(4)),分离式Hopkinson拉杆平板挂钩试件关键结构参数(单位:mm)的最优组合为:

      图  12  最小目标值随遗传代数的变化

      Figure 12.  The change of minimum objective function in each generation

      $ \left( {{L_1}},\;{{W_1}},\;{{L_2}},\; {{W_2}},\;R,\;T \right) = \left({9.6},\;{2.0},\;{14.3},\; {4.9},\;{1.8},\;{1.5} \right) $

    • 根据智能协同优化的结果,试件的最优结构尺寸如图13所示(单位:mm),用于有限元模拟验证的1/4模型的网格如图14所示。优化试件的分离式Hopkinson拉伸有限元模拟结果如下:

      图  13  优化后试件的结构

      Figure 13.  The structure and dimensions of the optimized specimen

      图  14  优化试件的1/4模型网格

      Figure 14.  The 1/4 meshed model of the optimized specimen

      (1) 应力平衡达到时间

      应力波到达试件标距段末端(轴向应力开始大于1 MPa)时起,直至试件标距段前后端(图14)轴向应力偏差小于5%,历时24.75 μs,即应力平衡达到时间为24.75 μs。

      (2) 非轴向应力水平

      试件中点的轴向应力达到最大值时,计算法向的应力与轴向应力的比值,得到在整个优化路径(见图14)上应力比值的均值为0.018 8。

      (3) 变形均匀程度

      单波加载完成后,试件轴向应变的分布如图15所示,可以看出优化后试件的变形较均匀地集中在标距段,标距段各点应变值的方差为0.135 1×10−3

      图  15  优化试件轴向应变云图

      Figure 15.  Distribution of axial strain in the optimized specimen

      (4) 过渡段变形

      拉伸波单波加载完成后优化试件前过渡段、标距段、后过渡段的变形量以及过渡段相对变形如表5所示。

      前过渡段变形量标距段变形量后过渡段变形量过渡段相对变形
      0.16 mm2.03 mm0.16 mm13.62%

      表 5  优化试件各段变形量及过渡段相对变形

      Table 5.  The deformation of each section of the optimized specimen and the relative deformation of the transition zone

      表6所示为智能协同优化前后试件各项指标的对比,除应力平衡达到时间有所变长以外,其余各项指标均有明显改善,试件测量精度的整体性能得以提升。

      指标优化前优化后增大(↑)或降低(↓)
      应力平衡达到时间E/μs22.8124.758.51%↑
      应变方差V/10−30.598 80.135 177.44%↓
      过渡段相对变形 D0.202 00.136 232.57%↓
      非轴向应力水平 N0.056 30.018 866.61%↓

      表 6  优化前后各项指标的比较

      Table 6.  The comparison of various indicators before and after the optimization

    • 为验证基于优化算法和有限元仿真所得最优试件结构的有效性,采用数字图像相关(digital image correlation, DIC)技术对其进行实验验证。配合高速摄像的DIC技术因其可以快速捕捉结构件的变形、振动等信号,可应用于冲击测量等高速、高应变率场合[4]

      图16所示为加载前的最优结构试件,在试件测量区域喷洒散斑以提高实验的测量精度。DIC技术难以测量试件内部的应力状态,但可以测量试件表面的变形和应变信号,因此可以对试件标距段变形的均匀度以及过渡段的相对变形进行验证。为保证试件在视场中的完整性,采用图16所示的视角进行拍摄测量。

      图  16  加载前的最优结构试件

      Figure 16.  The specimen with optimal structure before loading

      调整霍普金森拉杆设备气室气压,使得试件应变率为1 500 s−1,选取试件中间区域轴向应变最大的一帧图像做后处理,所得试件轴向应变的分布如图17所示,图中 AB分别为试件测量段和过渡段交界处在试件厚度方向上的中点,连接AB两点构成路径AB,提取路径AB上各点轴向应变值沿轴中心线的分布并与优化后试件轴向应变分布的计算结果进行对比,结果如图18所示,红色曲线是实验所得的的轴向应变曲线,黑色曲线是数值模拟所得轴向应变分布曲线。由图18可以看出,实验所得轴向应变沿路径AB的分布与计算结果基本吻合。轴向应变方差的实验结果为0.145 5×10−3,与计算结果相差仅7.70%。

      图  17  试件轴向应变分布云图

      Figure 17.  Distributed cloud map of axial strain of the specimen

      图  18  轴向应变分布数值模拟与实验结果的比较

      Figure 18.  Comparison of axial strain distribution between simulation and experiment

      图19所示是与图17相同时刻,试件轴向变形(单位:mm)的分布云图,从图中可以看出,试件各点的轴向变形值从左至右逐渐增大,图中点AB为前过渡段的分界线在试件厚度方向上的中点,点CD为后过渡段分界线在试件厚度方向的中点。实验时,追踪加载过程中AB以及CD四点的轴向位移,在图19所示时刻,B点的轴向位移差A点的轴向位移可得前过渡段的轴向变形,C点的轴向位移差D点的轴向位移可得后过渡段的轴向变形,D点轴向位移差B点轴向位移可得标距段的轴向变形。将实验所得各段的轴向变形以及过渡段相对全段的相对变形与计算结果进行对比,结果如表7所示。由表7可以看出,前过渡段的相对误差最大,为18.75%,其余各项相对偏差均小于15%;整体上,计算值与实验值之间的相对偏差处于较低水平,可以认为计算所得过渡段的变形及其相对全段的相对变形具有有效性和正确性。

      图  19  试件轴向变形分布云图

      Figure 19.  Contour of axial deformation of the specimen

      变形段计算值实验值相对偏差
      前过渡段0.16 mm0.13 mm18.75%
      后过渡段0.16 mm0.15 mm6.25%
      标距段2.03 mm2.12 mm4.43%
      相对变形13.62%11.67%14.32%

      表 7  各段变形值实验与计算的比较

      Table 7.  Comparison of deformation of transition zones between simulations and experiments

      表8所示,为参考试件与优化试件关于过渡段相对变形和标距段变形均匀度在仿真和实验结果上的对比。从表8可以看出,参考试件、优化试件的计算结果和实验结果均可较好地吻合,计算结果具有一定程度的有效性;对参考试件与优化试件的实验结果进行对比,优化后试件过渡段的相对变形减小了47.69%,标距段变形的均匀度则提升了73.28%,这一定程度上证实了优化后试件结构的优越性和本文工作的有效性。

      指标参考试件优化试件
      计算实验相对偏差计算实验相对偏差
      过渡段相对变形/%20.2022.3110.45%13.6211.6714.32%
      标距段变形均匀度/10−30.598 80.544 69.05%0.135 10.145 57.70%

      表 8  参考试件与优化试件数值模拟和实验结果的对比

      Table 8.  The comparison of simulated and experimental results between reference specimen and optimized specimen

    • 本文中针对平板挂钩试件在拉伸过程中因试件结构几何效应引起的数据测量误差问题,提出衡量平板挂钩试件测量精度的指标;通过正交试验设计、BP神经网络、遗传算法相结合的智能协同优化,可以得到以下结论。

      (1)分离式Hopkinson拉杆实验中关于实现精确测量的假设:试件处于应力平衡状态、只承受轴向应力、变形均匀以及过渡段不发生变形实际情况下均不能满足,但是可以将标距段应力平衡达到时间、应变方差、非轴向应力水平以及过渡段相对变形作为试件在测量精度方面的性能指标。

      (2)结合正交试验的样本数据库,BP神经网络能够实现平板挂钩试件结构参数与测量精度指标之间的非线性建模,建立的模型具有有效性和准确性;遗传算法也可以就试件的结构参数实现全局寻优。结果证明,采用正交试验设计、BP神经网络和遗传算法相结合的智能协同优化方法对铝合金平板挂钩试件的优化具有可行性;同时,该方法优化方案对其它不同材料、不同连接形式的试件优化具有可借鉴性。

      (3)采用智能协同优化方案对平板挂钩试件的结构参数进行协同优化,针对直径为19 mm的拉杆,试件采用AA5182型铝合金在平均应变率为1 500 s−1时得到试件结构参数的最优组合为:标距段长L1=9.6 mm,标距段宽W1=2.0 mm,挂钩与杆接触面长L2=14.3 mm,挂钩与杆接触面宽W2=4.9 mm,过渡段半径R=1.8 mm,试件厚度T=1.5 mm。

      (4)通过有限元模拟的验证分析,优化后的试件与参考试件相比应力平衡达到时间增加了8.51%,试件变形均匀程度指标提高了77.44%,过渡段相对的变形量降低了32.57%,非轴向应力水平指标降低了66.61%。总体上,除应力平衡达到时间变长之外,其余各指标的变化均有利于试件测量精度的提升。最后,经过实验验证,基于有限元模拟和优化算法所得平板挂钩试件的最优结构具有一定的有效性。

参考文献 (20)

目录

    /

    返回文章
    返回