Loading [MathJax]/jax/output/HTML-CSS/jax.js
  • ISSN 1001-1455  CN 51-1148/O3
  • EI、Scopus、CA、JST收录
  • 力学类中文核心期刊
  • 中国科技核心期刊、CSCD统计源期刊

基于微裂纹界面摩擦生热的点火模型

楼建锋 张延耿 洪滔 周婷婷 郭少冬

李刚, 王晓放, 殷硕, 李文亚. 粒子入射角度对冷喷涂涂层形成的影响[J]. 爆炸与冲击, 2007, 27(5): 477-480. doi: 10.11883/1001-1455(2007)05-0477-04
引用本文: 楼建锋, 张延耿, 洪滔, 周婷婷, 郭少冬. 基于微裂纹界面摩擦生热的点火模型[J]. 爆炸与冲击, 2015, 35(6): 807-811. doi: 10.11883/1001-1455(2015)06-0807-05
LI Gang, WANG Xiao-fang, YIN Shuo, LI Wen-ya. Study of incidence angle of particle on its coating formation in cold spraying[J]. Explosion And Shock Waves, 2007, 27(5): 477-480. doi: 10.11883/1001-1455(2007)05-0477-04
Citation: Lou Jian-feng, Zhang Yan-geng, Hong Tao, Zhou Ting-ting, Guo Shao-dong. Study on the model of hot-spot ignition based on friction generated heat on the microcrack face[J]. Explosion And Shock Waves, 2015, 35(6): 807-811. doi: 10.11883/1001-1455(2015)06-0807-05

基于微裂纹界面摩擦生热的点火模型

doi: 10.11883/1001-1455(2015)06-0807-05
基金项目: 国家自然科学基金项目(11302031, 11372053, 11402031);中国工程物理研究院科学技术发展基金项目(2012A0101004, 2014B0101014, 2014A0201008)
详细信息
    作者简介:

    楼建锋(1980—), 男, 博士, 副研究员

    通讯作者:

    张延耿, zhang_yangeng@iapcm.ac.cn

  • 中图分类号: O389

Study on the model of hot-spot ignition based on friction generated heat on the microcrack face

  • 摘要: 开展了基于微裂纹界面摩擦生热的细观点火模型研究,采用有限元方法对包含化学反应放热和摩擦生热的热传递方程进行了离散求解,计算模型中考虑了炸药颗粒熔化对升温过程的影响。着重分析了点火模型中主要参数(热点尺度、应变率和界面压力)对炸药点火的影响规律。数值研究表明,随着热点尺度的增大,热点的温度上升越快,越容易发生点火;应变率越大或者界面压力越高,热量积累越快,炸药越容易点火。
  • 薄壁金属构件由于其低成本和高吸能性而在汽车结构中作为一种性能良好的缓冲吸能元件得到广泛的应用。为了确保碰撞过程中乘员和汽车主要部件的安全, 薄壁构件作为受撞时主要的承载和吸能元件[1], 在耐撞性设计中, 薄壁构件应能够以稳定可控的变形方式吸收尽可能多的能量; 同时最大限度的降低碰撞过程中薄壁构件所承受的压溃载荷, 较高的压溃载荷会使汽车在碰撞过程中产生较大的加速度, 使车内乘员受到严重的伤害, 故在碰撞过程最大压溃载荷应尽可能的小[2]

    诱导结构是薄壁构件应力集中的地方, 可以控制薄壁构件变形形式和降低碰撞时的压溃载荷, 从而使乘员免受高压溃力, 同时也降低了其它部件所承受的压溃载荷。Y.B.Cho等[3]采用数值方法对具有孔洞型诱导结构的前梁装置进行了抗撞性模拟研究, 计算表明设定合理的诱导结构能提高能量吸收同时降低压溃载荷; 张涛等[4]研究了薄壁组合结构及其设计缺陷结构的吸能特性, 提出设置一定的诱导缺陷可以降低碰撞过程中的峰值载荷。

    本文中提出在薄壁圆管上增加圆弧形凹槽的诱导变形原则, 研究其对抗撞性的影响。在碰撞分析中大多采用最大峰值压溃力为评价指标, 为了提高结构的抗撞性, 同时满足轻量化的设计要求, 综合考虑比吸能、最大峰值压溃力及压溃力的平均值, 建立多目标薄壁构件优化任务, 并结合径向基函数法构造优化目标的近似函数, 采用理想点法求解多目标优化问题, 详细分析凹槽数量与其半径对薄壁构件抗撞性的影响, 以期为吸能构件抗撞性设计提供依据。

    (1) 最大峰值压溃力Fmax:在压溃过程中沿轴向产生的压溃载荷的最大值[5]。最大峰值压溃力会引起较高的加速度, 是碰撞过程中需要考虑的重要参数。为了减小对乘员的伤害, 故希望在整个压溃过程中最大峰值压溃力尽可能的小。

    (2) 平均压溃力Fm:薄壁构件在碰撞过程中的平均压溃力为:

    Fm=δmax0F dδδmax (1)

    式中:F为压溃力, δ为压溃位移, δmax最大压溃位移。在碰撞过程中压溃力一般围绕平均压溃力波动。

    (3) 压溃力效率η, 即平均压溃力与最大压溃力之比, 其表达式为:

    η=FmFmax (2)

    压溃力效率η综合考虑了平均压溃力与最大压溃力对薄壁构件抗撞性的影响, η越大说明最大压溃力越小越接近平均压溃力, 压溃力曲线平稳, 理想的吸能构件η应为1。

    (4) 比吸能χ, 单位质量的薄壁结构所吸收的能量[6]为:

    χ=EW (3)

    式中:E为薄壁结构所吸收的总能量, W为薄壁结构的总质量。为了提高结构的抗撞性, 同时满足轻量化的设计要求, 应使薄壁构件的比吸能最大化。

    以比吸能和压溃力效率为优化设计指标, 在优化设计过程中, 有2个设计参数即凹槽数量m和凹槽半径r, 它们在一定范围内变化。该多目标优化问题可表示为

    {maxη(m,r)maxχ(m,r) s. t. mminmmmaxrminrrmax (4)

    事实上, 目标函数很难用数学方法精确的表达, 本文中采用径向基函数模型构造其近似表达式, 并采用理想点法寻找多目标问题最优解。

    建立薄壁金属圆管模型, 如图 1所示, 圆管下端固定, 上端自由并承受质量25kg, 冲击速度v=15.49m/s的刚性体冲击作用, 圆管厚度t=0.8mm, 长度L=90mm, 直径d=31mm。

    图  1  薄壁圆管模型
    Figure  1.  Analysis model of the thin-walled cylinder

    薄壁构件的材料选为高强度钢, 密度ρ=7.82×103 kg/m3, 弹性模量E=207.2GPa, 泊松比ν=0.3, 屈服应力σ0=446MPa。高强度钢材料的动态变形受材料应变率的影响较大, 需要在材料模型中考虑应变率的影响。一般采用Cowper-Symonds塑性材料模型[7]:

    σy=σ0(1+ε/c)1/p (5)

    式中:σy为考虑应变率之后的动态屈服应力, σ0为静态屈服应力, ε为应变率, c=40s-1p=5为该模型中与应变率相关的参数, 静态应力应变变化曲线如图 2所示。

    图  2  静态应力应变关系曲线
    Figure  2.  Relation between static stress and strain

    为了验证有限元分析的精确性, 与文献[7]中的实验数据进行比较, 表 1给出了薄壁构件原模型由LS-DYNA计算得到的最大峰值压溃力、总吸能、最大压溃位移及其与文献[7]中对应的实验值。通过对比可以看出有限元结果和实验值之间的相对误差很小, 数值模拟的结果与实验值能够较好地吻合。

    表  1  有限元模型结果与实验结果[7]对比
    Table  1.  Comparison of experimental result[7] and finife element model
    Fmax/kN E/kJ δmax/mm
    实验[7] 有限元 ε/% 实验[7] 有限元 ε/% 实验[7] 有限元 ε/%
    120.30 120.09 0.2 2.95 2.94 0.3 64.10 58.40 8.9
    下载: 导出CSV 
    | 显示表格

    诱导结构是薄壁构件应力集中的地方, 使结构的屈曲变形稳定。为了提高原模型在碰撞过程中的抗撞性, 本文对原模型结构做了改进, 如图 3所示, 沿光滑表面的薄壁结构(原模型结构尺寸不变)添加均布的圆弧形凹槽。现分别施加1~4个半径为0.8mm诱导槽, 在前文的工况下对4个模型进行实验及有限元数值模拟计算得到其碰撞响应, 如表 2所示。由表 2可知, 改进后的模型与原模型相比η值增大, 最大峰值压溃力显著降低了, 平均压溃力却没有明显减少, 压溃力曲线趋于平稳; 而吸能几乎没有受到影响, 由于增加凹槽导致结构质量的增加, 因而比吸能χ有所减小, 可见添加圆弧形凹槽, 显著提高了薄壁结构的抗撞性。

    图  3  改进结构模型
    Figure  3.  The improved structural model
    表  2  改进模型的有限元分析结果及实验值对比
    Table  2.  Comparison of the improved model between experiment and simulation
    m Fmax/kN Fm/kN η E/kJ χ/(kJ·kg-1)
    实验 有限元 实验 有限元
    0 120.30 120.09 50.36 0.420 0 2.95 2.94 55.60
    1 80.26 77.82 46.14 0.592 9 2.93 2.91 54.29
    2 91.63 88.39 42.48 0.480 6 2.87 2.90 53.90
    3 94.32 89.26 42.32 0.474 1 2.90 2.93 54.26
    4 90.53 87.12 40.92 0.469 7 2.91 2.94 54.24
    下载: 导出CSV 
    | 显示表格

    为了进一步研究凹槽结构对抗撞性的影响, 取凹槽数m=1~4, 半径r=0.8~3.2mm。采用全因子实验设计, 在mr的设计域内, 选用了28个样本点研究其抗撞性, 并对凹槽结构进行优化设计。表 2所示为原模型和4个改进模型的有限元分析结果与实验值对比, 表 2m=0表示未加凹槽, 为原模型, m=1~4表示加入0.8mm凹槽的数量, 也即是4个改进模型, 从数据对比看有限元分析结果与实验值较接近, 从而验证了改进模型有限元分析的可靠性。为了节约资源, 对余下的24个设计样本点只进行有限元分析计算从而得到其碰撞响应。

    该多目标优化问题数学表达式可进一步表示为:

    {minη(m,r)minχ(m,r) s. t. 1m40.8 mmr3.2 mm (6)

    径向函数[8-11]是以预测点与样本点间的欧氏距离为自变量的一族函数。以径向函数为基函数, 通过线性组合构造出来的近似模型即为径向基函数模型。

    在设计空间中, 关于设计变量响应函数的近似表达式可以定义为:

    ˜y(x)=ni=1λiφ(xxi) (7)

    式中:λT=(λ1, λ2, …, λn)为权系数, ‖x-xi‖为预测点和样本点间的欧氏距离, φ(‖x-xi‖)为径向函数。将n个样本点及每个样本点对应的响应函数值代入到式(7), 可以得到方程组矩阵形式为:

    Y=Φλ (8)

    式中: Y=[˜y(x1),˜y(x2),,˜y(xn)]T,Φ=[Φij]=[φ(xixj)](i,j=1,2,,n)

    在得到响应函数的径向基函数模型后, 需要对近似模型精度进行验证, 通过额外测试点(不包括样本点)决定系数R2和误差平方和方均根 ¯R2来验证径向基函数模型对设计变量的拟合情况。决定系数R2和误差平方和方均根 ¯R2分别定义为:

    R2=1ni=1(yi˜yi)2ni=1(yiˉyi)2 (9)
    ¯R2=ni=1(yi˜yi)2n (10)

    式中:n是测试点的个数, 其中yii个测试点的有限元分析结果yi的平均值, ˜yi是第i个测试点的径向基函数近似值。

    多目标优化问题在数学上一般可表示为

    {minf1(x1,,xn)minfp(x1,,xn) s. t. gi(x)0,i=1,2,,mhi(x)=0,i=1,2,,l (11)

    先求解p个单目标问题

    minxDfj(x),j=1,2,,p (12)

    设其最优值为fj*, 称f*=(f1*, …, fp*)T为值域中的一个理想点, 因为一般很难达到, 故在期望的某种度量下, 寻求距离f*最近的f作为近似值。最直接的方法是构造评价函数并极小化, 既求解:

    minxDφ[f(x)]=pi=1[fi(x)fi]2 (13)

    并将它的最优解x*作为式(8)在这种意义下的最优解。

    于是多目标优化问题可进一步表示为

    {minf1(x)=χ(x)minf2(x)=η(x) s. t. 1m40.8 mmr3.2 mm (14)

    对于带有圆弧形凹槽的薄壁构件, 凹槽数量m和半径r是影响结构抗撞性能的2个重要参数, 通过有限元计算得到了28个样本点的碰撞响应, 从而得到ηχ关于设计变量mr的响应函数和响应面, 如图 4~5所示。从图中可见凹槽数量半径对结构的比吸能和压溃力效率有很大影响, 呈现出一定的非线性关系, 不具备单调性。

    图  4  比吸能与变量mr的变化关系
    Figure  4.  Specific energy absorption varied with m and r
    图  5  压溃力效率与变量mr的变化关系
    Figure  5.  Crushing force efficiency varied with m and r

    表 3给出了χη这2个单目标函数的优化值, 可见2个单目标函数最优解不同, 2个目标函数很难同时达到最优。采用理想点法对式(14)的多目标优化问题进行求解, 得到带有圆弧形凹槽薄壁构件的最优值为m=3和r=1.224mm, χ=53.723kJ/kg, η=0.673。图 6分别给出了采用有限元方法计算得到的最大的比吸能和压溃力效率的结构, 以及采用理想点法计算得到多目标问题的优化结构。

    图  6  优化后薄壁结构模型
    Figure  6.  Optimized thin-walled structure models
    表  3  单目标函数优化结果
    Table  3.  Optimums of single objective functions
    单目标函数 m r/mm χ/(kJ·kg-1) η
    max χ 1 0.800 54.29 0.592 9
    max η 4 1.198 51.94 0.682 5
    下载: 导出CSV 
    | 显示表格

    优化后构件的η=0.673, 比表 2m=0的原模型的η=0.420 0提高了52.90%, 而χ=53.723kJ/kg, 与55.06kJ/kg相比仅降低了4.05%。图 7~8分别给出了该多目标问题达到最优时, 在碰撞过程中压溃力和吸收能量随压溃位移和压溃时间变化情况, 并将其与相应条件的原模型进行比较。从压溃载荷图 7可以看出, 与原模型相比优化后构件最大峰值压溃力降低了52.1%, 最大峰值压溃力更接近平均压溃力, 压溃力曲线平稳, 可见以η为优化设计指标是合理的, η综合考虑了最大峰值压溃力和压溃力曲线的平稳性。图 8所示为总吸能与压溃时间的关系, 可以看出, 优化后构件所吸收的能量并没有很大的减少, 在碰撞时间为10ms时, 优化后构件所吸收的能量与同时间原模型吸收的能量比较相近了, 可见优化后构件的抗撞性得到了显著提高。

    图  7  压溃力随压溃位移变化关系
    Figure  7.  Crushing force varied with crushing displacement
    图  8  总吸能随压溃时间变化关系
    Figure  8.  Total energy absorption varied with crushing time

    优化后构件在碰撞过程中发生叠缩变形如图 9所示, 这种变形模式有利于结构在变形过程中吸收更多的碰撞动能。带有圆弧形诱导凹槽结构的薄壁构件变形首先在凹槽处开始, 最大的塑性弯矩和塑性铰发生在凹槽处, 薄壁结构在该处容易压溃。凹槽能控制塑性变形, 塑性应变在该处达到最大值, 这也就意味着沿着变形的凹槽处能消散大量的压溃力, 大部分的冲击能在该处以塑性变形的方式被吸收。可见增加圆弧形凹槽诱导结构后, 可大大提高薄壁构件的抗撞性能。

    图  9  优化后构件叠缩变形
    Figure  9.  Progressive folding deformation of optimized component

    对薄壁圆管结构进行改进引入圆弧形诱导凹槽, 以比吸能和压溃力效率为评价指标, 建立了兼顾最大峰值压溃力、平均压溃力及比吸能的多目标优化问题的数学模型, 研究凹槽数量和半径对抗撞性能的影响。数值分析的结果表明:引入圆弧形诱导凹槽结构在保证不减少薄壁构件吸能能力的情况下达到了减小最大峰值压溃力、使压溃力曲线平稳的目的。优化结果表明最优构件的最大峰值压溃力比原模型减小了52.1%, 压溃力曲线与原模型相比更趋于平稳, 从而有效的保证了碰撞过程中乘员的安全性; 从吸能角度来看, 最优构件所吸能的能量并没有很大减少, 在碰撞时间为10ms时, 最优构件所吸收的能量达到了与同一时间原模型吸收的能量; 而最优构件的比吸能只比原模型减少了4.05%, 可见引入圆弧形诱导凹槽提高了薄壁结构的抗撞性, 为进一步研究吸能构件的抗撞性设计提供了依据。

  • 图  1  微裂纹面热点模型示意图

    Figure  1.  Configuration of hot-spot model showing friction generated heat

    图  2  一维热传递单元

    Figure  2.  The one dimensional heat transfer element

    图  3  不同热点尺度情况下热点中心温度随时间的变化

    Figure  3.  Temperature at center of hot-spot vs. time for the cases of different hot-spot sizes

    图  4  不同时刻距离裂纹面不同位置的温度分布

    Figure  4.  Temperature distribution vs. distance from crack face at various times

    图  5  不同应变率下点火时间的比较

    Figure  5.  Effect of strain rate on ignition time for a constant pressure

    图  6  不同界面压力下点火时间的比较

    Figure  6.  Effect of pressure on ignition time for a constant strain rate

  • [1] Asay Blaine W. Non-shock initiation of explosives[M]. Heidelberg: Springer-Verlag, 2010: 15-18.
    [2] Bowden F P, Yoffe A D. Initiation and growth of explosives in liquids and solids[M]. Cambridge: Cambridge University Press, 1952.
    [3] Bowden F P, Yoffe A D. Hot spots on rubbing surfaces and the detonation of explosives by friction[J]. Proceedings of the Royal Society of London, Series A: Mathematical & Physical Sciences, 1947, 188(10): 329-349. doi: 10.1098/rspa.1947.0012
    [4] Amosov A P, Bostandzhiyan S A, Kozlov V S. Ignition of solid explosives by the heat of dry friction[J]. Fizika Goreniya i Vzryva, 1972, 8(3): 362-368. http://www.onacademic.com/detail/journal_1000034907003910_0866.html
    [5] Amosov A P, Bostandzhiyan S A, Kozlov V S, et al. Mechanism of heating up and ignition of solid explosives due to external friction as a result of mechanical stimulations[J]. Fizika Goreniya i Vzryva, 1976, 12(5): 699-703. doi: 10.1007/BF00743166
    [6] Wiegand D A, Redingius B, Ellis K, et al. Pressure and friction dependent mechanical strength-cracks and plastic flow[J]. International Journal of Solids and Structures, 2011, 48(11/12): 1617-1629. http://www.sciencedirect.com/science/article/pii/S0020768311000424
    [7] Wiegand D A, Redingius B, Ellis K, et al. Evidence for fricgtion between crack surfaces during deformation of dcomposite plastic bonded explosives[C]∥Elert M L, Buttle W T, Furnish M D, et al. Proceedings of Shock Compression of Condensed Matter-2009. Nashville, Tennessee, 2009: 349-352.
    [8] Wiegand D A, Redingius B. The role of friction in the mechanical failure properties of a polymer particulate composite[C]∥APS March Meeting. New Orleans, 2008.
    [9] 陈文.高速侵彻条件下战斗部装药安全性研究[D].北京: 北京理工大学, 2009.
    [10] Boyle V, Frey R, Blake O. Combined pressure shear ignition of explosive[C]∥Lee E L, Short J M. Proceedings of the 9th International Detonation Symposium. Oregon, Portland, 1989: 3-17.
    [11] Frey R B. The initiation of explosive charges by rapid shear[C]∥Proceedings of the 7th International Detonation Symposium. Annapolis, Maryland, 1981: 36-42.
    [12] Dienes J K. A unified theory of flow, hot spots, and fragmentation with an application to explosive sensitivity[M]. New York: Springer, 1996: 366-398.
    [13] Dienes J K, Kershner J D. Multiple-shock initiation via statistical crack mechanics[C]∥Short J M, Kennedy J E. Proceeding of the 11th International Detonation Symposium. Snowmass, Colorado, 1998: 717-724.
    [14] Dienes J K, Kershner J D. Crack dynamics and explosive burn via generalized coordinates[J]. Journal of Computer-Aided Materials Design, 2000, 7(3): 217-237. doi: 10.1023/A:1011874909560
    [15] Bennett J G, Haberman K S, Johnson J N, et al. A constitutive model for the non-shock ignition and mechanical response of high explosives[J]. Journal of Mechanical Physical Solids, 1998, 46(12): 2303-2322. doi: 10.1016/S0022-5096(98)00011-8
    [16] Linan A, Williams F A. Theory of ignition of a reactive solid by constant energy flux[J]. Combustion Science and Technology, 1971, 3(2): 91-98. doi: 10.1080/00102207108952276
  • 期刊类型引用(7)

    1. 雷成,李福胜,缪得祥,朱涛. 基于回归分析法的矩形槽圆管吸能特性研究. 机械强度. 2024(02): 292-300 . 百度学术
    2. 王春华,杨雨泽,安达,唐治. 液压支架凹角圆管式薄壁构件吸能特性分析与优化. 机械设计与制造. 2024(10): 266-271+276 . 百度学术
    3. 邓敏杰,刘志芳. 仿马尾草薄壁结构的设计与耐撞性研究. 高压物理学报. 2022(03): 111-120 . 百度学术
    4. 张欣玥,惠旭龙,葛宇静,舒挽,白春玉,刘小川. 中低速压缩加载下不同截面构型复合材料薄壁结构吸能特性及失效分析. 爆炸与冲击. 2022(06): 36-49 . 本站查看
    5. 霍鹏,许述财,范晓文,李建平,杨欣,黄晗. 鹿角骨单位仿生薄壁管斜向冲击耐撞性研究. 爆炸与冲击. 2020(11): 127-138 . 本站查看
    6. 谭丽辉,谭洪武,鲁帅. 诱导缺陷结构在薄壁构件耐撞性中的应用. 河南科技. 2016(21): 74-75 . 百度学术
    7. 李玉如,柳忠彬,肖守讷,王欢. 预制螺旋槽薄壁管的吸能研究. 四川理工学院学报(自然科学版). 2015(02): 21-25 . 百度学术

    其他类型引用(16)

  • 加载中
图(6)
计量
  • 文章访问数:  4280
  • HTML全文浏览量:  436
  • PDF下载量:  806
  • 被引次数: 23
出版历程
  • 收稿日期:  2014-04-22
  • 修回日期:  2014-07-24
  • 刊出日期:  2015-12-10

目录

/

返回文章
返回