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

力学环境约束下剪切销可靠性分析及优化设计

王若冰 赵志军 肖和业 陈家文 蒋晓磊

引用本文:
Citation:

力学环境约束下剪切销可靠性分析及优化设计

    作者简介: 王若冰(1986- ),男,博士,高级工程师,robbin203@163.com;
  • 中图分类号: O347; TJ45; TB114.3

Reliability analysis and design optimization of a shear pin constrained by mechanical boundaries

  • CLC number: O347; TJ45; TB114.3

  • 摘要: 剪切销是火工装置关键部件,其可靠性不仅表现为点火作用下可靠剪断,还表现为受力学环境激励不发生断裂。本文中以多项式混沌展开方法为基础,建立了力学环境约束下的剪切销分析模型,结合序贯优化与可靠性分析方法,提出了剪切销可靠性优化设计的思路。以某型火工作动装置为应用实例,依据实用的力学环境,进行了剪切销可靠性分析及优化设计,揭示了设计参数与力学环境之间的关系,并获得了影响可靠性的关键参数。最后,开展了优化后的火工作动装置实验测试,结果佐证了优化设计的有效性。
  • 图 1  SORA-PMA方法流程示意图

    Figure 1.  Flowchart of the SORA-PMA method

    图 2  火工作动装置内部结构图

    Figure 2.  Internal structure of the explosive actuated device

    图 3  容差、目标函数和归一化约束函数值的迭代收敛过程

    Figure 3.  Convergence history of tolerance, objective function and normalized constraints

    图 4  火工作动装置工作后的照片

    Figure 4.  A photo of an explosive-actuated device after work

    表 1  火工作动装置工作环境

    Table 1.  Work situation of the pyrotechnic pin

    序号工作状态过载类型加速度作用时间/ms要求
    1点火工作静压力剪切销可靠剪断
    2起落冲击后峰锯齿波20g11剪切销可靠锁止
    3发射冲击半正弦波3 000g2剪切销可靠锁止
    4载机加速恒加速度10g剪切销可靠锁止
    5挂载飞行随机振动参见国家军用标准GJB 150.16A—2009附件中
    表C.4[22]42和图C.10[22]53中直升机挂飞振动条件
    剪切销可靠锁止
    下载: 导出CSV

    表 2  影响火工作动装置可靠性的因素

    Table 2.  Influence factors on the pyrotechnic pin reliability

    序号参数符号参数名称参数均值参数标准差
    1ppow火药工作产生压强10 MPa2.5 MPa
    2S1火药压力作用面积500 mm25 mm2
    3d剪切销截面直径2.5 mm0.1 mm
    4E7075铝合金弹性模量70 GPa0.005 GPa
    5σb7075铝合金极限强度429 MPa5 MPa
    6σfatigue7075铝合金真实断裂强度627 MPa5 MPa
    7εfatigue7075铝合金真实断裂延性0.290.000 5
    8μ泊松比0.330.000 5
    9l剪切面长度3.5 mm0.1 mm
    10ξ损耗因子0.0050.001
    11mcylinderpin圆柱销质量15 g0.1 g
    12acylinderpin圆柱销加速度200 m/s210 m/s2
    13tvib振动时长10 h1 h
    14σmean平均预紧力400 MPa5 MPa
    15aimpact起落冲击(发射冲击)加速度幅值200(30 000) m/s2100 m/s2
    16timpact起落冲击(发射冲击)时长11(3) ms1 ms
    下载: 导出CSV

    表 3  变量对不同工作环境的不确定性灵敏度

    Table 3.  Uncertainty sensitivity of variables in different mechanic environments

    参数点火工作产品发射载机加速载机起落挂载飞行综合影响系数
    τ1σ1σ2τ2τ3σ3σ4τ4σ5
    ppow0.2600.0000.0000.0000.0000.0000.0000.0000.0000.029
    S10.3020.0000.0000.0000.0000.0000.0000.0000.0000.034
    d0.2890.2260.1990.1780.3780.2330.2360.2580.2540.250
    E0.0000.1060.0920.1280.0000.1100.0640.0720.0630.070
    σb0.0000.0000.0000.0000.0000.0000.0000.0000.0620.007
    σfatigue0.0000.0000.0000.0000.0000.0000.0000.0000.0460.005
    εfatigue0.0000.0000.0000.0000.0000.0000.0000.0000.0570.006
    μ0.0000.0000.0000.1070.0000.0000.0000.0710.0000.020
    l0.0000.0000.1260.0990.0000.0000.1350.0430.0000.045
    ξ0.0000.0190.0120.0100.0000.0160.0190.0170.0340.014
    mcylinderpin0.0000.1120.1150.0900.2020.1210.0650.0650.0600.092
    acylinderpin0.0000.0000.0000.0000.2030.0000.0000.0000.0000.023
    tvib0.0000.0000.0000.0000.0000.0000.0000.0000.1210.014
    σmean0.0000.0000.0000.0000.0000.0000.0000.0000.1360.015
    aimpact0.0000.1730.1560.1300.0000.1580.1730.1590.0000.106
    timpact0.0000.1970.1570.1330.0000.1780.1520.1410.0000.107
     注:τ1τ2τ3τ4为剪切应力,σ1σ3为弯曲应力,σ2σ4为拉伸应力,σ5为疲劳破坏应力。
    下载: 导出CSV

    表 4  设计变量初值、取值范围和优化值

    Table 4.  Initial values, ranges and optimal values of design variables

    设计变量单位初值优化值取值范围
    ppowMPa10.0013.51[8, 15]
    dmm2.501.49[1, 2]
    lmm3.502.62[1, 3]
    σmeanMPa400300[300, 400]
    下载: 导出CSV
  • [1] 荣吉利, 宋乾强, 张涛. 一种预测航天火工装置可靠性的小样本方法 [J]. 宇航学报, 2015, 36(3): 360–364. DOI: 10.3873/j.issn.1000-1328.2015.03.016.
    RONG Jili, SONG Qianqiang, ZHANG Tao. A small sample method for predicting reliability of space pyrotechnic devices [J]. Journal of Astronautics, 2015, 36(3): 360–364. DOI: 10.3873/j.issn.1000-1328.2015.03.016.
    [2] 荣吉利, 宋乾强, 张涛, 等. 冗余航天火工装置可靠性评估方法及应用 [J]. 宇航学报, 2013, 34(7): 1021–1026. DOI: 10.3873/j.issn.1000-1328.2013.07.019.
    RONG Jili, SONG Qianqiang, ZHANG Tao. Reliability assessment method for redundant aerospace pyrotechnic devices and its application [J]. Journal of Astronautics, 2013, 34(7): 1021–1026. DOI: 10.3873/j.issn.1000-1328.2013.07.019.
    [3] 董海平, 蔡瑞娇, 穆慧娜. 火工品可靠性计量计数评估方法的有效性研究 [J]. 含能材料, 2008, 16(5): 553–555. DOI: 10.3969/j.issn.1006-9941.2008.05.020.
    DONG Haiping, CAI Ruijiao, MU Huina. Validity of variables attributes assessment method for reliability of initiating explosive devices [J]. Chinese Journal of Energetic Materials, 2008, 16(5): 553–555. DOI: 10.3969/j.issn.1006-9941.2008.05.020.
    [4] 董海平, 赵霞, 蔡瑞娇. 基于信息等值的火工品可靠性评估小样本方法 [J]. 兵工学报, 2011, 32(5): 554–558.
    DONG Haiping, ZHAO Xia, CAI Ruijiao. An assessment method of reliability of initiating devices based on information measure equivalency [J]. Acta Armamentarii, 2011, 32(5): 554–558.
    [5] 董海平, 蔡瑞娇, 穆慧娜, 等. 火工品可靠性的感度参数设计法 [J]. 爆炸与冲击, 2009, 29(6): 613–616. DOI: 10.11883/1001-1455(2009)06-0613-04 .
    DONG Haiping, CAI Ruijiao, MU Huina, et al. A sensitivity parameter design method for reliability of explosive initiators [J]. Explosion and Shock Waves, 2009, 29(6): 613–616. DOI: 10.11883/1001-1455(2009)06-0613-04 .
    [6] 伊枭剑, 董海平, 翟志强, 等. 基于应力-强度干涉模型的火工品可靠性设计方法 [J]. 北京理工大学学报, 2014, 34(10): 1007–1011.
    YI Xiaojian, DONG Haiping, ZHAI Zhiqiang, et al. Reliability design for initiating devices based on stress-strength interference model [J]. Transactions of Beijing Institute of Technology, 2014, 34(10): 1007–1011.
    [7] OUCH R, UKRITCHON B, PIPATPONGSA T, et al. Finite element analyses of the stability of a soil block reinforced by shear pins [J]. Geomechanics and Engineering, 2017, 12(6): 1021–1046.
    [8] 孙洁, 郭崇星, 吴瑞德, 等. 拔销器作用过程的仿真研究 [J]. 火工品, 2017(1): 10–13. doi: 10.3969/j.issn.1003-1480.2017.01.004
    SUN Jie, GUO Chongxing, WU Ruide, et al. Simulation research on the function process of the pin puller [J]. Initiators and Pyrotechnics, 2017(1): 10–13. doi: 10.3969/j.issn.1003-1480.2017.01.004
    [9] 高滨. 火工作动装置设计参数的敏感性分析 [J]. 航天返回与遥感, 2006, 27(3): 57–60. doi: 10.3969/j.issn.1009-8518.2006.03.011
    GAO Bin. Design parameter sensitivity analysis of pyrotechnically actuated devices [J]. Space Craft Recovery and Remote Sensing, 2006, 27(3): 57–60. doi: 10.3969/j.issn.1009-8518.2006.03.011
    [10] COLANTUONO T, GROSSONI I, ALLEN P, et al. The optimal strength of shear pins: requirements in point run-throughs [J]. Journal of Rail and Rapid Transit, 2016, 232(2): 2–19. DOI: 10.1177/0954409716682657.
    [11] SALTELLI A, RATTO M, ANDRES T, et al. Global sensitivity analysis: the primer [M]. Chichester: John Wiley and Sons, 2008.
    [12] HELTON J C, JOHNSON J D, SALLABERRY C J, et al. Survey of sampling-based methods for uncertainty and sensitivity analysis [J]. Reliability Engineering and System Safety, 2006, 91(10/11): 1175–1209.
    [13] IMAN R L, HELTON J C. An investigation of uncertainty and sensitivity analysis techniques for computer models [J]. Risk Analysis, 2010, 8(1): 71–90.
    [14] WIENER N. The homogeneous chaos [J]. American Journal of Mathematics, 1938, 60(4): 897–936. doi: 10.2307/2371268
    [15] DU Xiaoping, GUO Jia, BEERAM H. Sequential optimization and reliability assessment for multidisciplinary systems design [J]. Structural and Multidisciplinary Optimization, 2008, 35(2): 117–130. DOI: 10.1007/s00158-007-0121-7.
    [16] 崔卫东, 关云翔, 惠蕾. 拔销器的小型化设计 [J]. 火工品, 2006(2): 17–19. doi: 10.3969/j.issn.1003-1480.2006.02.005
    CUI Weidong, GUAN Yunxiang, HUI Lei. Design of the small-sized piston actuator [J]. Initiators and Pyrotechnics, 2006(2): 17–19. doi: 10.3969/j.issn.1003-1480.2006.02.005
    [17] 汪凤泉. 电子设备振动与冲击手册 [M]. 北京: 科学出版社, 1998.
    [18] 苟文选. 材料力学 [M]. 北京: 科学出版社, 2005.
    [19] 徐鹏. 金属材料应变寿命曲线估算的新方法[D]. 南京: 南京航空航天大学, 2012.
    XU Peng. A new method for the estimation of strain-life curve of metals [D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2012.
    [20] 陈小前, 姚雯, 欧阳琦. 飞行器不确定性多学科设计优化理论与应用 [M]. 北京: 科学出版社, 2013.
    [21] 王若冰, 谷良贤, 龚春林. 随机-区间混合不确定性分层序列化多学科可靠性分析方法 [J]. 西北工业大学学报, 2016(1): 139–146. DOI: 10.3969/j.issn.1000-2758.2016.01.021.
    WANG Ruobing, GU Liangxian, GONG Chunlin. A stratified sequencing multi-disciplinary reliability analysis method under random and interval uncertainty [J]. Journal of Northwestern Polytechnical University, 2016(1): 139–146. DOI: 10.3969/j.issn.1000-2758.2016.01.021.
    [22] 施荣明、朱广荣、吴飒, 等. 军用装备实验室环境试验方法: 第16部分: 振动试验: GJB 150. 16A-2009 [S]. 北京: 总装备部军标出版发行部, 2009: 42; 53.
  • [1] 樊正复姚尧孟冲 . 计算接力式爆破网路可靠性的一种方法. 爆炸与冲击, 1994, 14(3): 283-288.
    [2] 董海平蔡瑞娇穆慧娜曹建华 . 火工品可靠性的感度参数设计法. 爆炸与冲击, 2009, 29(6): 613-616. doi: 10.11883/1001-1455(2009)06-0613-04
    [3] 王丹于海江 . 试验设计在爆炸网络可靠性研究中的应用. 爆炸与冲击, 2015, 35(2): 184-190. doi: 10.11883/1001-1455(2015)02-0184-07
    [4] 曹建华蔡瑞娇董海平李意起 . 电火工品桥丝与药剂升温模型及在可靠性设计中的应用. 爆炸与冲击, 2004, 24(1): 90-95.
    [5] 秦栋泽范宁军 . 自毁装置的安全性和可靠性. 爆炸与冲击, 2014, 34(1): 111-114. doi: 10.11883/1001-1455(2014)01-0111-04
    [6] 王军波高敏 . 引信安全系统过程识别理论及其可靠性. 爆炸与冲击, 1999, 19(1): 44-49.
    [7] 刘明安伟光宋向华 , . 随机参数超空泡射弹的动力稳定性和可靠性. 爆炸与冲击, 2013, 33(5): 525-530. doi: 10.11883/1001-1455(2013)05-0525-06
    [8] 洪东跑赵宇温玉全 . 基于虚拟样本的火工品可靠性评估. 爆炸与冲击, 2009, 29(6): 669-672. doi: 10.11883/1001-1455(2009)06-0669-04
    [9] 温玉全洪东跑王玮 . 基于试验熵的火工品可靠性评估理论与方法研究. 爆炸与冲击, 2007, 27(6): 553-556. doi: 10.11883/1001-1455(2007)06-0553-04
    [10] 殷之平李玉龙黄其青 . 含诱导缺陷薄壁圆管耐撞性优化设计. 爆炸与冲击, 2011, 31(4): 418-422. doi: 10.11883/1001-1455(2011)04-0418-05
    [11] 谭丽辉徐涛崔晓梅张炜赵世佳 . 带有圆弧形凹槽金属薄壁圆管抗撞性优化设计. 爆炸与冲击, 2014, 34(5): 547-553. doi: 10.11883/1001-1455(2014)05-0547-07
    [12] 王瑞利刘全温万治 . 非嵌入式多项式混沌法在爆轰产物JWL参数评估中的应用. 爆炸与冲击, 2015, 35(1): 9-15. doi: 10.11883/1001-1455(2015)01-0009-07
    [13] 侯海量朱锡古美邦 . 爆炸载荷作用下加筋板的失效模式分析及结构优化设计. 爆炸与冲击, 2007, 27(1): 26-33. doi: 10.11883/1001-1455(2007)01-0026-08
    [14] 林俊德 . 仪器车减震台的设计优化. 爆炸与冲击, 1993, 13(2): 145-152.
    [15] 王继海 . 多项式形式Mie-Grneisen物态方程及其等熵线. 爆炸与冲击, 1992, 12(1): 1-10.
    [16] 樊成飞王耀华王强 . 爆燃动力装置内弹道的优化设计. 爆炸与冲击, 2015, 35(2): 267-272. doi: 10.11883/1001-1455(2015)02-0267-06
    [17] 卢维娴陆在庆 . 钛合金应力波铆接中绝热剪切的显微分析. 爆炸与冲击, 1985, 5(1): 67-72.
    [18] 李磊马宏昊沈兆武 . 基于正交设计方法的双锥罩结构优化设计. 爆炸与冲击, 2013, 33(6): 567-573. doi: 10.11883/1001-1455(2013)06-0567-07
    [19] 徐全军白帆李裕春龚自正张庆明 . 微声聚能穿孔装置优化设计及IED 安全销毁实验. 爆炸与冲击, 2012, 32(1): 79-83. doi: 10.11883/1001-1455(2012)01-0079-05
    [20] 顾强张世豪安晓红张亚 . 基于灰色理论的两点爆炸起爆参数优化设计. 爆炸与冲击, 2015, 35(3): 359-365. doi: 10.11883/1001-1455(2015)03-0359-07
  • 加载中
图(4)表(4)
计量
  • 文章访问数:  35
  • HTML全文浏览量:  64
  • PDF下载量:  5
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-04-27
  • 录用日期:  2018-08-09
  • 网络出版日期:  2019-06-25
  • 刊出日期:  2019-07-01

力学环境约束下剪切销可靠性分析及优化设计

    作者简介:王若冰(1986- ),男,博士,高级工程师,robbin203@163.com
  • 西安现代控制技术研究所,陕西 西安 710065

摘要: 剪切销是火工装置关键部件,其可靠性不仅表现为点火作用下可靠剪断,还表现为受力学环境激励不发生断裂。本文中以多项式混沌展开方法为基础,建立了力学环境约束下的剪切销分析模型,结合序贯优化与可靠性分析方法,提出了剪切销可靠性优化设计的思路。以某型火工作动装置为应用实例,依据实用的力学环境,进行了剪切销可靠性分析及优化设计,揭示了设计参数与力学环境之间的关系,并获得了影响可靠性的关键参数。最后,开展了优化后的火工作动装置实验测试,结果佐证了优化设计的有效性。

English Abstract

  • 火工装置在制导武器、航天等领域有极其重要的应用[1-2],它能够在相当短的时间内释放出相当大的能量,破坏剪切销以完成预定的动作。火工装置作为航天、武器等产品中的子产品,也会受到振动、冲击、恒加速过载等力学环境激励,如果剪切销在此情况下失效,火工装置不能正常工作,会导致整个飞行任务的失败,甚至危及人员和设备的安全。为了保证火工装置的可靠性,不仅需要关注火工装药感度、点火可靠性,还需要对剪切销进行可靠性设计并进行评估分析。

    针对火工装置的可靠性设计、分析、评估,学者们已开展了大量的研究。荣吉利等[1-2]对航天火工装置可靠性的小样本预测方法以及可靠性方法在冗余火工装置中的应用、评估开展了研究。董海平等[3-5]、伊枭剑等[6]研究了计量计数法、火工品可靠性小样本评估方法及感度参数设计法。以上研究主要将火工品作为一个整体对象进行可靠性方面的分析,未针对剪切销具体结构进行可靠性分析。在火工装置剪切销的力学性能研究中,学者们主要关注剪切销在火工作用下失效过程仿真分析[7-8]、优化设计[9-10],并未考虑力学环境激励对剪切销性能及可靠性的影响。综上所述,针对力学环境下的剪切销可靠性分析、设计还没有相应的理论方法。因此,开展力学环境约束下剪切销可靠性分析及优化设计的研究,不仅能够为火工装置可靠性设计及分析方法提供支撑,还能用于指导火工装置中剪切销的设计,提升火工装置的环境适应性。

    由于力学环境约束下剪切销的可靠性分析是多参数解耦的可靠性问题,采用微分和基于抽样仿真的方法已不适用,需借鉴其他可靠性分析设计领域的灵敏度分析方法[11-13]进行突破性的研究。因此,本文中基于多项式混沌展开(polynomial chaos expansions, PCE)方法[14],建立力学环境约束下剪切销可靠性分析模型。然后,结合SORA (sequential optimization and reliability assessment)优化算法[15],以火药工作压强最小为目标函数,提出力学环境约束下剪切销可靠性优化设计思路。以某型火工装置为实例,进行力学环境下剪切销可靠性的灵敏度分析,并对剪切销的参数开展可靠性的优化设计。最后,通过抽样实验佐证优化设计结果的正确性,验证本文中剪切销可靠性分析及优化设计的正确性。

    • 在剪切销的可靠性分析模型中,作以下简化:(1)在保证火工品安全性的条件下,将火工品的工作发火的可靠性设定为0.999 9;(2)剪切销及其它结构件不考虑缺陷、裂纹情况;(3)不考虑温度、外部气压等自然环境变化对结构的影响;(4)结构件之间的摩擦系数为常数。

    • 剪切销力学环境可以分为点火工作及环境激励两类情况。点火工作状态指火工品产生压力剪断剪切销失效的工况;环境激励状态指剪切销随产品一起承受外部力学载荷激励的工况,如恒加速、稳态振动、冲击等条件。根据火工作动装置的工作环境,建立剪切销在不同载荷下的分析模型。

    • 以摩擦力因素作为剪切力的修正项[16],可得点火工作时剪切销所受到的剪切力:

      $ {\tau _{{\rm{ignit}}}} = \frac{{{p_{{\rm{pow}}}}{S_1}-{F_{{\rm{friction}}}}}}{{{S_2}}} $

      式中:ppow为火药工作产生的压强;S1为火药压力作用面积;S2为剪切销截面积,当剪切销截面为直径d的圆形时,S2 = πd2/4;Ffriction为密封环与壁面的摩擦力。

    • 剪切销和被固定结构在动力学环境下为近似线性的单振子系统,则剪切销受冲击时的最大响应加速度[17-18]:

      $ {\ddot Z_{\left( {\rm{B,L,S}} \right)\max }} = \frac{{J{\omega _{\rm{n}}}_{\left( {{\rm{B,L,S}}} \right)}{{\rm{e}}^{ - \xi {\omega _{\rm{n}}}_{\left( {{\rm{B,L,S}}} \right)}{t_{\max }}}}}}{{\sqrt {1 - {\xi ^2}} }}\sin \left( {{\omega _{{\rm{n}}\left( {\rm{B,L,S}} \right)}}\sqrt {1 - {\xi ^2}} {t_{\max }} + \varphi '} \right),\quad{\omega _{{{\rm{n}}_{\left( {{\rm{B,L,S}}} \right)}}}} = {n_{{\rm{shearpin}}}}\sqrt {\frac{{{K_{\left( {{\rm{B,L,S}}} \right)}}}}{{{m_{{\rm{cylinderpin}}}}}}} $

      式中:下标B、L、S分别表示弯曲、拉伸、剪切状态;J为加速度幅值为aimpt的阶跃函数在冲击时长timpt结束时的积分值;${\omega _{\rm{n}}}$为固有频率;ξ为损耗因子;$\varphi '$为相位角;tmax为冲击响应时间;${n_{{\rm{shearpin}}}}$为剪切销个数;${m_{{\rm{cylinderpin}}}}$为圆柱销质量;KBKLKS为剪切销的弯曲、拉伸、剪切刚度,且${K_{\rm{B}}} = E{W_{\rm{y}}}$${K_{\rm{L}}} = E\displaystyle\frac{{{S_2}}}{l}$${K_{\rm{S}}} = G\displaystyle\frac{{{S_2}}}{l} = \displaystyle\frac{{E{S_2}}}{{2l\left( {1 + \mu } \right)}}$E为弹性模量,${W_{\rm{y}}}$为截面惯性矩,l为剪切面长度,G为剪切模量,μ为泊松比。

      剪切销在加速度${\ddot Z_{\left( {\rm S} \right)\max }}$${\ddot Z_{\left( {\rm L} \right)\max }}$${\ddot Z_{\left( {\rm B} \right)\max }}$下所受应力:

      $ {\tau _{{\rm{impact}}\left( {\rm{S}} \right)}} = \frac{{{m_{{\rm{cylinderpin}}}}{{\ddot Z}_{\left( {\rm{S}} \right)\max }}}}{{{S_2}}},\quad{\sigma _{{\rm{impact}}\left( {\rm{L}} \right)}} = \frac{{{m_{{\rm{cylinderpin}}}}{{\ddot Z}_{\left( {\rm{L}} \right)\max }}}}{{{S_2}}},\quad{\sigma _{{\rm{impact}}\left( {\rm{B}} \right)}} = \frac{{{m_{{\rm{cylinderpin}}}}{{\ddot Z}_{\left( {\rm{B}} \right)\max }}}}{{{S_2}}} $

    • 恒加速度环境下剪切销所受切应力:

      ${\tau _{{\rm{accelerate}}}} = \frac{{{a_{{\rm{cylinderpin}}}}{m_{{\rm{cylinderpin}}}}}}{{{S_2}}}$

      式中:${a_{{\rm{cylinderpin}}}}$为圆柱销加速度。

    • 疲劳线性累积系数为:

      ${K_{{\rm{fatigue}}}} = \sum\limits_i^{} {\frac{{{n_i}}}{{{N_i}}} = } \sum\limits_i^{} {\frac{{{t_{{\rm{vib}}}}{f_i}}}{{{N_i}\left( {{\varepsilon _a}\left( {{f_i}} \right)} \right)}}} $

      式中:ni为振动时长tvib内的振动次数;Ni为应变εa(fi)对应的破坏次数,${f_i}$为振动频率。本文采用Manson-Coffin方法[19]求得εa对应的Ni,,且有:

      ${\varepsilon _a} = \frac{{{{\sigma '}_{{{\rm{fatigue}}}}} - {\sigma _{{\rm{mean}}}}}}{E}{\left( {2{N_i}} \right)^b} + {\varepsilon '_{{{\rm{fatigue}}}}}{\left( {2{N_i}} \right)^c}$

      式中:σmean为平均预紧力,参数${\sigma '_{{{\rm{fatigue}}}}}$${\varepsilon '_{{{\rm{fatigue}}}}}$bc采用四点关联法[19]求得。

      又有:

      ${\varepsilon _a} = \frac{{{k_{\rm{c}}}A\left( {{f_i}} \right){m_{{\rm{cylinderpin}}}}}}{{{S_2}E}}$

      式中:A(·)为激励加速度关于频率的函数;kc为线性弹簧振子的放大系数。设z=fi/ωn,则[17]:

      ${k_{\rm{c}}} = \frac{{\sqrt {1 + {{\left( {2\xi {\textit z}} \right)}^2}} }}{{\sqrt {{{\left( {1 - {{\textit z}^2}} \right)}^2} + {{\left( {2\xi {\textit z}} \right)}^2}} }}$

      联立式(6)和(7)可得fi对应的破坏次数Ni

      在以上力学环境条件下,如果剪切销的应力超过材料的许用强度就会发生破坏。在点火工作状态下剪切销破坏表明工作可靠,而在环境激励状态下剪切销不发生破坏表明使用可靠。因此,剪切销的可靠性是两种状态可靠性的综合,设计时需要同时考虑并进行灵敏度分析,找出可靠性设计中的关键参数。

    • 对于随机变量y,PCE方法将其近似展开成如下级数形式[20]

      $y = {c_0}{\varGamma _0} + \sum\limits_{{i_1} = 1}^\infty {{c_{{i_1}}}{\varGamma _1}({x_{{i_1}}}) + \sum\limits_{{i_1} = 1}^\infty {\sum\limits_{{i_2} = 1}^{{i_1}} {{c_{{i_1}{i_2}}}{\varGamma _2}({x_{{i_1}}},{x_{{i_2}}})} } } + \sum\limits_{{i_1} = 1}^\infty {\sum\limits_{{i_2} = 1}^{{i_1}} {\sum\limits_{{i_3} = 1}^{{i_2}} {{c_{{i_1}{i_2}{i_3}}}{\varGamma _3}({x_{{i_1}}},{x_{{i_2}}},{x_{{i_3}}})} } } + \cdots $

      式中:y为PCE模型响应值,x为随机变量;c为确定性系数;Γpp阶正交多项式。本文中采用Hermite多项式作为正交多项式,适用于系统响应是高斯随机过程、变量是高斯分布的情况。依据各随机设计变量产生的样本点建立二阶PCE模型,以最小二乘法确定PCE模型的系数:

      ${C_i} = {c_i} + 2\sqrt {{c_{ii}}} + \sum\limits_{j = 1}^{i - 1} {\sqrt {{c_{ij}}} } $

      并根据式(10)求得各变量对不同工作环境下应力值的不确定性灵敏度。

      本文中假设各可靠性约束的权值相同(也可根据实际情况确定各约束的权值),则各系数归一化后,以算数平均计算综合影响系数。

    • 本文中以火药工作压强最小为目标函数,建立火工作动装置可靠性设计优化模型如下:

      $\left\{ {\begin{array}{*{20}{l}} {{{{X}}_{{\rm{opt}}}} = \mathop {\min }\limits_{X \in D} F({{X}}) } \\ {D = {G_1} \cap {G_2} \cap \cdots \cap {G_9}} \end{array}} \right.$

      式中:${{{X}}_{{\rm{opt}}}}$为优化解向量,${{X}}$为设计变量ppowdlσmean组成的向量;$F({{X}}) $为目标函数,本文中选取火药工作产生压强ppow为目标函数;D为设计变量的取值范围,是约束条件集合G1G9的交集:

      $\left\{ {\begin{array}{*{20}{l}} {{G_1} = \left\{ {{{X}}|P\left( {{\sigma _{{\rm{fatigue}}}}-{\tau _1}\left( {{X}} \right) \text{≤}0} \right) \text{≥} R} \right\}} \\ {{G_2} = \left\{ {{{X}}|P\left( {{\tau _2}\left( {{X}} \right)-{\sigma _{\rm{b}}}\text{≤}0} \right) \text{≥} R} \right\}} \\ {{G_3} = \left\{ {{{X}}|P\left( {{\tau _3}\left( {{X}} \right)-{\sigma _{\rm{b}}} \text{≤}0} \right) \text{≥} R} \right\}} \\ {{G_4} = \left\{ {{{X}}|P\left( {{\tau _4}\left( {{X}} \right)-{\sigma _{\rm{b}}} \text{≤} 0} \right) \text{≥} R} \right\}} \\ {{G_5} = \left\{ {{{X}}|P\left( {{\sigma _1}\left( {{X}} \right)-{\sigma _{\rm{b}}} \text{≤} 0} \right) \text{≥} R} \right\}} \\ {{G_6} = \left\{ {{{X}}|P\left( {{\sigma _2}\left( {{X}} \right)-{\sigma _{\rm{b}}} \text{≤} 0} \right) \text{≥} R} \right\}} \\ {{G_7} = \left\{ {{{X}}|P\left( {{\sigma _3}\left( {{X}} \right)-{\sigma _{\rm{b}}} \text{≤} 0} \right) \text{≥}R} \right\}} \\ {{G_8} = \left\{ {{{X}}|P\left( {{\sigma _4}\left( {{X}} \right)-{\sigma _{\rm{b}}} \text{≤} 0} \right) \text{≥}R} \right\}} \\ {{G_9} = \left\{ {{{X}}|P\left( {{K_{{\rm{fatigue}}}}\left( {{X}} \right) - 1 \text{≤} 0} \right) \text{≥} R} \right\}} \end{array}} \right.$

      式中:P(·)为满足括号内条件的概率;R为设计可靠度;τ1为火药工作环境下剪切销的剪应力;${\sigma _1}$${\sigma _2}$${\tau _2}$分别为起降冲击载荷下的弯曲应力、拉伸应力和切应力;${\tau _3}$分别载机恒加速环境下的切应力;${\sigma _3}$${\sigma _4}$${\tau _4}$分别为发射环境中冲击载荷下的弯曲应力、拉伸应力和切应力;${K_{{\rm{fatigue}}}}$为挂载飞行环境下剪切销累积疲劳系数。

    • 本文中采用SORA方法[20]的框架,以性能测度方法(performance measurement approach,PMA)[21]求解约束函数的可靠性分析问题,在可靠性分析中引入PCE模型以提高求解效率,其求解流程如图1所示,原可靠性优化设计问题被分解为确定性优化和可靠性分析两部分迭代求解的问题[21]

      图  1  SORA-PMA方法流程示意图

      Figure 1.  Flowchart of the SORA-PMA method

    • 本文中以某机载火工作动装置为研究对象,来验证分析方法的有效性。该主要由电起爆器、本体、锁销、剪切销、压螺等组成。工作前,锁销在剪切销作用下保持伸出本体状态。工作时,电起爆器产生燃气,推动锁销向上运动,剪断剪切销,直至锁销收缩至本体内,完成预定动作。其内部结构如图2 所示。

      图  2  火工作动装置内部结构图

      Figure 2.  Internal structure of the explosive actuated device

    • 火工作动装置用于机载武器系统,需满足:(1)在载机起降、产品发射、载机恒加速、长时间挂载飞行等工作环境下可靠锁止;(2)在电起爆器工作压力为8~15 MPa时,剪切销被可靠剪断。根据国家军用标准GJB 150.16A—2009[22],火工作动装置使用环境与要求如表1所示。将力学环境中的条件量化为影响火工作动装置可靠性的因素,并根据工程经验和实验数据,定义各参数服从正态分布,如表2所示。发射冲击与起落冲击的都可通过式(2)~(3)计算获得,仅幅值和时间不同。

      序号工作状态过载类型加速度作用时间/ms要求
      1点火工作静压力剪切销可靠剪断
      2起落冲击后峰锯齿波20g11剪切销可靠锁止
      3发射冲击半正弦波3 000g2剪切销可靠锁止
      4载机加速恒加速度10g剪切销可靠锁止
      5挂载飞行随机振动参见国家军用标准GJB 150.16A—2009附件中
      表C.4[22]42和图C.10[22]53中直升机挂飞振动条件
      剪切销可靠锁止

      表 1  火工作动装置工作环境

      Table 1.  Work situation of the pyrotechnic pin

      序号参数符号参数名称参数均值参数标准差
      1ppow火药工作产生压强10 MPa2.5 MPa
      2S1火药压力作用面积500 mm25 mm2
      3d剪切销截面直径2.5 mm0.1 mm
      4E7075铝合金弹性模量70 GPa0.005 GPa
      5σb7075铝合金极限强度429 MPa5 MPa
      6σfatigue7075铝合金真实断裂强度627 MPa5 MPa
      7εfatigue7075铝合金真实断裂延性0.290.000 5
      8μ泊松比0.330.000 5
      9l剪切面长度3.5 mm0.1 mm
      10ξ损耗因子0.0050.001
      11mcylinderpin圆柱销质量15 g0.1 g
      12acylinderpin圆柱销加速度200 m/s210 m/s2
      13tvib振动时长10 h1 h
      14σmean平均预紧力400 MPa5 MPa
      15aimpact起落冲击(发射冲击)加速度幅值200(30 000) m/s2100 m/s2
      16timpact起落冲击(发射冲击)时长11(3) ms1 ms

      表 2  影响火工作动装置可靠性的因素

      Table 2.  Influence factors on the pyrotechnic pin reliability

      表2中各设计参数的分布参数产生10 000样本点,并根据式(10)求得各变量对不同工作环境下应力值的不确定性灵敏度。各变量不确定性灵敏度如表3所示:

      参数点火工作产品发射载机加速载机起落挂载飞行综合影响系数
      τ1σ1σ2τ2τ3σ3σ4τ4σ5
      ppow0.2600.0000.0000.0000.0000.0000.0000.0000.0000.029
      S10.3020.0000.0000.0000.0000.0000.0000.0000.0000.034
      d0.2890.2260.1990.1780.3780.2330.2360.2580.2540.250
      E0.0000.1060.0920.1280.0000.1100.0640.0720.0630.070
      σb0.0000.0000.0000.0000.0000.0000.0000.0000.0620.007
      σfatigue0.0000.0000.0000.0000.0000.0000.0000.0000.0460.005
      εfatigue0.0000.0000.0000.0000.0000.0000.0000.0000.0570.006
      μ0.0000.0000.0000.1070.0000.0000.0000.0710.0000.020
      l0.0000.0000.1260.0990.0000.0000.1350.0430.0000.045
      ξ0.0000.0190.0120.0100.0000.0160.0190.0170.0340.014
      mcylinderpin0.0000.1120.1150.0900.2020.1210.0650.0650.0600.092
      acylinderpin0.0000.0000.0000.0000.2030.0000.0000.0000.0000.023
      tvib0.0000.0000.0000.0000.0000.0000.0000.0000.1210.014
      σmean0.0000.0000.0000.0000.0000.0000.0000.0000.1360.015
      aimpact0.0000.1730.1560.1300.0000.1580.1730.1590.0000.106
      timpact0.0000.1970.1570.1330.0000.1780.1520.1410.0000.107
       注:τ1τ2τ3τ4为剪切应力,σ1σ3为弯曲应力,σ2σ4为拉伸应力,σ5为疲劳破坏应力。

      表 3  变量对不同工作环境的不确定性灵敏度

      Table 3.  Uncertainty sensitivity of variables in different mechanic environments

      (1)对于点火工作环境,参数对σ1灵敏度排序为ppowS1d,量级相当。

      (2)对于产品发射和载机起落的冲击过载环境,各参数对弯曲应力、拉伸应力和剪切应力的灵敏度排序具有一致性,几何参数dl的灵敏度最高,对冲击载荷参数aimpacttimpact的灵敏度次之,对材料参数Eμ${m_{{\rm{cylinderpin}}}}$的灵敏度较小。

      (3)对于载机恒加速环境,剪切销所受剪切应力对参数灵敏度排序为d${a_{{\rm{cylinderpin}}}}$${m_{{\rm{cylinderpin}}}}$,量级相当。

      (4)对于挂载飞行环境,剪切销的疲劳破坏强度对几何参数d的灵敏度最高,对平均预紧力σmean和挂飞时长tvib的灵敏度次之,对材料参数Eσbm${\varepsilon _{{\rm{fatigue}}}}$${\sigma _{{\rm{fatigue}}}}$ξ的灵敏度较小。

      (5)综合各工作环境的灵敏度可知,几何参数的灵敏度最高,输入载荷参数的灵敏度次之,材料参数的灵敏度较小。

      将上述参数灵敏度分析结果与其对应的实际工程意见进行联系,综合考虑可知:

      (1)参数可靠度灵敏性分析可以将约束的可靠度转化为各设计变量和参数的标准差信息,这对设计中公差范围的选取和可靠度指标的分配具有指导作用。

      (2)灵敏度分析结果可以反映不同参数对约束条件的影响程度大小,因此在加工中结合工艺可行性分析,通过控制高灵敏度参数的公差,例如提高结构尺寸的加工精度等,可以有效控制产品性能和质量。

    • 设计可靠度设为R=99.99%,设置初始设计变量初值和取值范围如表4所示,收敛判定条件为当前最优点与上一步最优点之间的距离小于容差$1 \times {10^{-6}}$,且各约束函数的极限状态值小于等于0,满足以上条件时优化收敛并结束。本实例的优化问题经6次迭代收敛,优化结果如表4所示,设计变量初值为(10.00,2.50,3.50,400),优化值为(13.51,1.49,2.62,300)。容差的迭代收敛过程如图3(a)所示,由于采用SORA-PMA方法,确定性优化与可靠性分析序贯执行,迭代过程表现为振荡收敛。而目标函数的迭代收敛过程,在优化中出现先减小、后逐渐增大趋于收敛的现象,具体如图3(b)所示。这是由于第二步开始引入可靠性分析,考虑了变量和参数的不确定性,增大的火药压强才能满足可靠性约束条件,使得目标函数值变大然后收敛。归一化后的约束函数收敛过程如图3(c)所示,由图3(c)可知:(1)在迭代收敛结束时,所有约束的极限状态值均小于0,说明所得设计变量优化解满足所有约束条件,求解正确;(2)在优化过程中,约束G1G9接近于0,而其他约束远小于0,说明在本次取值空间内,G1G9为有效约束,最优解收敛于G1G9构成的曲面边界内。综上所述,本次优化过程收敛性好,结果正确、有效。

      设计变量单位初值优化值取值范围
      ppowMPa10.0013.51[8, 15]
      dmm2.501.49[1, 2]
      lmm3.502.62[1, 3]
      σmeanMPa400300[300, 400]

      表 4  设计变量初值、取值范围和优化值

      Table 4.  Initial values, ranges and optimal values of design variables

      图  3  容差、目标函数和归一化约束函数值的迭代收敛过程

      Figure 3.  Convergence history of tolerance, objective function and normalized constraints

    • 采用本文优化结果试制样品,剪切销直径取1.5 mm,伸出长度2.6 mm,电起爆器装药150 mg,产生压力8~15 MPa。每批次产品按照15%抽样依次进行随机振动、功能性冲击、强冲击、恒加速度等实验考核,然后实验样品工作点火,其他样品随武器系统进行实验考核。点火后的作动装置图片如图4所示,电起爆器点火破坏剪切销有效推动锁销缩入本体结构内。同批次产品通过力学实验及随武器系统实验考核的样本,其可靠度及其置信度评估公式如下:

      图  4  火工作动装置工作后的照片

      Figure 4.  A photo of an explosive-actuated device after work

      $P\left( {{N_{{\rm{fatigue}}}} \text{≤} {N_r}} \right) = \sum\limits_{i = 0}^r {C_q^iR_{\rm{L}}^{\left( {q - i} \right)}{{\left( {1 - {R_{\rm{L}}}} \right)}^i}} = 1 - \alpha $

      式中:q为火工作动装置实验共使用的样本量,Nfatigue为实验失败数;$P\left( {{N_{{\rm{fatigue}}}} \text{≤} {N_r}} \right)$表示当失效数Nfatigue不超过Nr时,在置信度$\alpha $的条件下,火工作动装置的可靠度为${R_{\rm{L}}}$。本次产品样本量为120件,抽样18件,力学实验点火及随导弹实验均可靠工作,失效数为0件。则根据式(13)可知,其可靠度为99%,置信度为70%。由于验证可靠度99.99%的产品,在80%置信度条件下的样本量需要达到16 000以上,难以通过实验验证,本文中仅能以小样本初步验证设计结果。

    • 为了提升火工装置对力学环境的适应性,解决力学环境下剪切销可靠性分析及优化设计的问题,本文中基于多项式混沌展开(polynomial chaos expansions, PCE)方法建立了力学环境约束下剪切销可靠性分析模型,结合SORA (sequential optimization and reliability assessment)-PMA(performance measurement approach)方法提出了剪切销可靠性优化设计思路。然后,以某型火工装置的剪切销作为实例,进行了力学环境下剪切销可靠性灵敏度分析,并以最小点火压力为目标函数,开展了力学环境约束下的可靠型优化设计。最后,按照优化后参数进行了剪切销的加工试制,组装成产品后经历抽样力学环境实验,再工作点火。实验结果佐证了优化设计结果的有效性,也表明本文中提出的可靠性分析及优化设计思路可行、有效。总结实例中的分析结果可知,对剪切销进行可靠性分析及优化设计还具有以下益处:

      (1)在剪切销设计中应用可靠性分析方法,可以将约束的可靠度转化为各设计变量和参数的标准差信息,这对设计中公差范围的选取和可靠度指标的分配具有指导作用。

      (2)可靠性灵敏度分析可以发掘剪切销设计参数与力学环境条件的内在联系,反映不同参数受力学环境的影响程度大小。同时结合加工中的制造工艺可行性分析,通过控制高灵敏度参数的公差,可以提高火工装置对力学环境的适应性。

      (3)可靠性设计优化结果可以给出关于装配参数的允许范围,如火工作动装置可靠性设计中对平均预紧力的要求,可以为装配工作提出更细致和具有可行性的装配工艺要求。

参考文献 (22)

目录

    /

    返回文章
    返回