CN106529016A - 一种爆热弹爆炸容腔半径的计算方法 - Google Patents

一种爆热弹爆炸容腔半径的计算方法 Download PDF

Info

Publication number
CN106529016A
CN106529016A CN201610976784.0A CN201610976784A CN106529016A CN 106529016 A CN106529016 A CN 106529016A CN 201610976784 A CN201610976784 A CN 201610976784A CN 106529016 A CN106529016 A CN 106529016A
Authority
CN
China
Prior art keywords
quick
radius
fried
thickness
thermoelastic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201610976784.0A
Other languages
English (en)
Other versions
CN106529016B (zh
Inventor
张建明
董明
王翔
谭凯元
赵维
樊星
范玉德
文尚刚
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
SICHUAN HAOKANG TECHNOLOGY Co.,Ltd.
Original Assignee
Institute of Chemical Material of CAEP
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Institute of Chemical Material of CAEP filed Critical Institute of Chemical Material of CAEP
Priority to CN201610976784.0A priority Critical patent/CN106529016B/zh
Publication of CN106529016A publication Critical patent/CN106529016A/zh
Application granted granted Critical
Publication of CN106529016B publication Critical patent/CN106529016B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种爆热弹爆炸容腔半径的计算方法，所设计的爆热弹在工作时，弹体最大应力都小于屈服强度，确保弹体处于弹性变形范围内，不发生屈服变形。即使因工作超压，弹体也能因屈服硬化不至于发生破片。合理的安全系数能够确保弹体处于安全工作状态。爆热弹外形规则，厚度均匀，传热速度较快。设计的爆热弹在量热系统中比热容最小，测量精度高。

Description

一种爆热弹爆炸容腔半径的计算方法
技术领域
本发明涉及一种爆热弹爆炸容腔半径的计算方法，属于爆热弹制作领域。
背景技术
目前，爆热弹设计方法鲜有公开文献报道，特别是如何确定爆炸容腔半径，即使有相关文献也侧重于爆热弹抗爆炸冲击性能的设计，或者用等当距离比来确定爆炸容腔半径。如何设计合适的内径，在保证爆热弹安全的前提下降低壁厚和爆热弹重量，降低系统比热容目前尚未有公开文献报道。
发明内容
本发明的发明目的是：提供一种爆热弹爆炸容腔半径的计算方法，不仅能够使爆热弹安全可靠的工作，同时能够最大程度的降低爆热弹重量、减薄壁厚和减小系统比热容，以克服现有技术的不足。
本发明是这样实现的：一种爆热弹爆炸容腔半径的计算方法，该设计方法通过如下步骤进行：
步骤一：首先建立系统的比热容模型。假设爆热弹的内腔半径r，内腔高度h，弹底厚度h1，弹盖厚度h2，壁厚δ，水套厚度t。并通过这些数据建立系统的比热容C；
C＝C水m水+C钢m钢 (1)
m水＝ρ水V水＝ρ水π[(r+δ+t)2(h+h1+h2+2t)-(r+δ)2(h+h1+h2)] (2)
m钢＝ρ钢V钢＝ρ钢π[(r+δ)2(h+h1+h2)-r2h] (3)
步骤二：假设内腔高度h为内腔半径r的2.4倍；
步骤三：根据爆炸当量选择确定合适的螺距P，100gTNT当量以下可以选择6或8，100gTNT当量及100gTNT当量以上选择8或12；选定螺距后可确定弹盖厚度h2；弹盖厚度h2为八倍螺距p加50mm；
步骤四：通过使用设计爆炸容器的动力系数法确定壁厚δ和内腔半径r关系；先确定爆热弹的制造材料，小当量宜选择低合金高强度钢，大当量宜选择超高强度钢；依据选择所用材料常温状态下的屈服强度除以安全系数作为材料许用应力，确保弹体工作在弹性变形范围,安全系数选3～5；其次根据K-G(Kinney and Graham)公式计算出爆炸容腔内反射超压ΔPm，根据反射超压ΔPm计算瞬间载荷Pr，再采用动力系数法计算等效静压Pe；最后按压力容器设计准则，依据等效静压Pe和材料许用应力计算出弹体壁厚δ；
步骤五:假设弹底壁厚h1为1.5倍壁厚δ；
步骤六:根据所设计的爆热测试仪假设水套厚度t，t为20～25mm；
通过以上假设，其中公式(1)、(2)、(3)中只有内腔半径r未知，因此系统比热容C只与内腔半径r相关，计算出系统比热容C最小时的内腔半径r；
步骤七：对计算出的内腔半径r进行圆整和验证，所述验证主要包括验证所取的螺距是否满足强度要求，验证最大应力是否小于材料许用应力或验证计算等效静压Pe应不大于100MPa。
假设爆炸容腔半径r，设计当量为Q，则当距比为
1)根据K-G(Kinney and Graham)计算出爆炸容腔内反射超压ΔPm；
其中P0为当地大气压。
2)根据反射超压计算瞬间载荷Pr为：
3)反射冲击波超压的作用时间τ和容器的自振角频率ω。
其中Q0为TNT热值，η为经验系数，爆热弹为柱对称取值为0.5，E和ρ分别为为材料的弹性模量和密度。
4)动力系数Cd由下式表达：
5)根据瞬间载荷采用动力系数法计算等效静压Pe
Pe＝PrCd (10)
所述弹体壁厚δ为：
其中di爆炸容腔的直径，[σ]t为材料许用应力，φ为焊接系数，爆热弹为整体加工，取值为1。
据绝热式爆热量热法的工作原理，在释放同等热量的情况下，系统的比热容越小，温升越高，系统的比热容越大，温升越低，考虑测温元件的不确定性，系统温升越高相对测量误差越小。因此通过合理设计爆热弹降低系统比热容比选择高精度温度传感器更容易提升测量精度。
爆热弹采用圆柱形结构，弹盖采用单螺母固定承压形式，整体外形规则均匀。根据爆炸冲击理论爆炸容腔内径越大，爆炸弹所需的壁厚越小，爆热弹外形质量越大，但量热过程包覆爆热弹所需的蒸馏水质量越大，系统比热容也越大。爆炸容腔越小，反射超压越高，弹体所需壁厚越厚，传热速度慢，弹体质量增加，系统比热容也会相应增加，因此合理的爆炸容腔半径是设计一个爆热弹关键。
爆热弹放置在内桶内，四周均布满蒸馏水，假设水套厚度为t，爆炸容腔半径为r，弹体壁厚为δ。测试样品主要为圆柱体，因此，爆炸容腔高度应略大于直径，一般为2.4r；弹底壁厚h1为壁厚δ的1.5倍；弹盖厚度h2为8倍螺距P加50mm。为此根据爆热弹初步外形建立系统比热容模型：
C＝C水m水+C钢m钢 (1)
弹体质量m钢由下式可得到
m钢＝ρ钢V钢＝ρ钢π[(r+δ)2(2.4r+1.5δ+8P+50)-2.4r3]
水套质量m水由下式可得到
m水＝ρ水V水＝ρ水π[(r+δ+t)2(1.5δ+2.4r+8P+50+2t)-(r+δ)2(1.5δ+2.4r+8P+50)]
螺距根据爆炸当量确定的，100gTNT当量以上的爆热弹螺距一般为8或12mm，小当量可以选更小的螺距，如6或8mm，在确定最终爆炸半径后需对螺纹强度进行校核。弹底壁厚δ和爆炸容腔半径r相关。水套厚度一般是确定值。因此系统比热容大小只与爆炸容腔半径相关。最终设计出的爆热弹能够在安全工作前提下系统比热容达到最低。
本发明与现有技术相比，具有以下的有益效果：设计的爆热弹在工作时，弹体最大应力都小于屈服强度，确保弹体处于弹性变形范围内，不发生屈服变形。即使因工作超压，弹体也能因屈服硬化不至于发生破片。合理的安全系数能够确保弹体处于安全工作状态。爆热弹外形规则，厚度均匀，传热速度较快。设计的爆热弹在量热系统中比热容最小，测量精度高。
优选的，所述爆热弹关键尺寸包括爆热弹的内腔半径r，内腔高度h，弹底厚度h1，弹盖厚度h2，壁厚δ。
附图说明
图1为爆热弹初步外形建立系统比热容模型示意图；
具体实施方式
本发明的实施例：
实施例1：
以100gTNT当量爆热弹为例，首先拟定水套厚度t为25mm，选择制造加工爆热弹的材料，本实施例选用的具体材料是30CrMnSiA.确定材料屈服强度为835MPa，安全系数取4；根据经验，螺纹螺距选择12mm，因此可以确定弹盖厚度为146mm；
爆热弹的设计步骤如下:
确定爆热弹制造材料，根据材料手册上该材料的屈服强度除以安全系数作为材料许用应力，确保弹体工作在弹性变形范围；根据安屈服强度设计，则最大工作应力不会超过屈服变形，认为在弹性变形范围内
假设爆炸容腔半径r，则当距比为
根据K-G(Kinney and Graham)计算出爆炸容腔内反射超压ΔPm；
根据反射超压计算瞬间载荷Pr为：
反射冲击波超压的作用时间τ和容器的自振角频率ω。
动力系数Cd由下式表达：
根据瞬间载荷采用动力系数法计算等效静压Pe
Pe＝PrCd
根据等效静压和材料许用应力按《钢制压力容器》设计准则计算弹体壁厚δ为
联立上述方程,则壁厚δ可转换为与半径r相关的方程，并带入系统比热容方程中。
以系统比热容最小为优化目标，计算出r值为90。
对爆炸容腔半径r进行圆整和校核，最终确定爆热容腔半径和爆热弹其他相关尺寸，其他相关尺寸的具体值是100。
实施例2：
以50gTNT当量爆热弹为例，首先拟定水套厚度t为25mm，选择制造加工爆热弹的材料，本实施例选用的具体材料是30CrMnSiA.确定材料屈服强度为835MPa，安全系数取4；根据经验，螺纹螺距选择8mm，因此可以确定弹盖厚度为114mm；
爆热弹的设计步骤如下:
确定爆热弹制造材料，根据材料手册上该材料的屈服强度除以安全系数作为材料许用应力，确保弹体工作在弹性变形范围；根据安屈服强度设计，则最大工作应力不会超过屈服变形，认为在弹性变形范围内
假设爆炸容腔半径r，则当距比为
根据K-G(Kinney and Graham)计算出爆炸容腔内反射超压ΔPm；
根据反射超压计算瞬间载荷Pr为：
反射冲击波超压的作用时间τ和容器的自振角频率ω。
动力系数Cd由下式表达：
根据瞬间载荷采用动力系数法计算等效静压Pe
Pe＝PrCd
根据等效静压和材料许用应力按《钢制压力容器》设计准则计算弹体壁厚δ为
联立上述方程,则壁厚δ可转换为与半径r相关的方程，并带入系统比热容方程中。
以系统比热容最小为优化目标，计算出r值为56。
对爆炸容腔半径r进行圆整和校核，最终确定爆热容腔半径和爆热弹其他相关尺寸，其他相关尺寸的具体值是60。
实施例3：
以300gTNT当量爆热弹为例，首先拟定水套厚度t为25mm，选择制造加工爆热弹的材料，本实施例选用的具体材料是AF1410，确定材料屈服强度为1480MPa，安全系数取4，则材料许用应力为370MPa；根据经验，螺纹螺距选择12mm，因此可以确定弹盖厚度为146mm；
爆热弹的设计步骤如下:
确定爆热弹制造材料，根据材料手册上该材料的屈服强度除以安全系数作为材料许用应力，确保弹体工作在弹性变形范围；根据安屈服强度设计，则最大工作应力不会超过屈服变形，认为在弹性变形范围内
假设爆炸容腔半径r，则当距比为
根据K-G(Kinney and Graham)计算出爆炸容腔内反射超压ΔPm；
根据反射超压计算瞬间载荷Pr为：
反射冲击波超压的作用时间τ和容器的自振角频率ω。
动力系数Cd由下式表达：
根据瞬间载荷采用动力系数法计算等效静压Pe
Pe＝PrCd
根据等效静压和材料许用应力按《钢制压力容器》设计准则计算弹体壁厚δ为
联立上述方程,则壁厚δ可转换为与半径r相关的方程，并带入系统比热容方程中。
以系统比热容最小为优化目标，计算出r值为134。
对爆炸容腔半径r进行圆整和校核，最终确定爆热容腔半径和爆热弹其他相关尺寸，其他相关尺寸的具体值是140。

Claims (5)

1.一种爆热弹爆炸容腔半径的计算方法，其特征在于：该方法通过如下步骤进行：
步骤一：首先建立系统的比热容模型，假设爆热弹的内腔半径r，内腔高度h，弹底厚度h1，弹盖厚度h2，壁厚δ，水套厚度t，并通过上述数据建立系统的比热容C；
C＝C水m水+C钢m钢 (1)
m水＝ρ水V水＝ρ水π[(r+δ+t)2(h+h1+h2+2t)-(r+δ)2(h+h1+h2)] (2)
m钢＝ρ钢V钢＝ρ钢π[(r+δ)2(h+h1+h2)-r2h] (3)
步骤二：假设内腔高度h为内腔半径r的2.4倍；
步骤三：根据爆炸当量选择合适的螺距P，100gTNT当量以下选择6或8，100gTNT当量及100gTNT当量以上选择8或12；选定螺距后可确定弹盖厚度h2；弹盖厚度h2为八倍螺距P加50mm；
步骤四：通过使用设计爆炸容器的动力系数法确定壁厚δ和内腔半径r关系；先确定爆热弹的制造材料，小当量选择低合金高强度钢，大当量选择超高强度钢；依据选择所用材料常温状态下的屈服强度除以安全系数作为材料许用应力，确保弹体工作在弹性变形范围,安全系数选3～5；其次根据K-G(Kinney and Graham)公式计算出爆炸容腔内反射超压ΔPm，根据反射超压ΔPm计算瞬间载荷Pr，再采用动力系数法计算等效静压Pe；最后按压力容器设计准则，依据等效静压Pe和材料许用应力计算出弹体壁厚δ与r的关系；
步骤五:假设弹底壁厚h1为1.5倍壁厚δ；
步骤六:根据所设计的爆热测试仪假设水套厚度t，t为20～25mm；
通过以上假设，其中公式(1)、(2)、(3)中只有内腔半径r未知，因此系统比热容C只与内腔半径r相关，计算出系统比热容C最小时的内腔半径r；
步骤七：对计算出的内腔半径r进行圆整和验证，所述验证主要包括验证所取的螺距是否满足强度要求，验证弹体最大应力是否小于材料许用应力且验证计算等效静压Pe应不大于100MPa。
2.根据权利要求1所述的爆热弹爆炸容腔半径的计算方法，其特征在于：所述步骤一中，先建立系统比热容模型，所述比热容模型与爆热弹关键尺寸和水套厚度t相关。
3.根据权利要求2所述的爆热弹爆炸容腔半径的计算方法，其特征在于：所述步骤二、三、四、五、六中，通过对爆热弹关键尺寸进行假设和转换，最终确定系统的比热容仅仅与内腔半径r相关。
4.根据权利要求1所述的爆热弹爆炸容腔半径的计算方法，其特征在于：所述步骤四中，所述材料许用应力通过常温下屈服强度除安全系数确定。
5.根据权利要求1所述的爆热弹爆炸容腔半径的计算方法，其特征在于，所述圆整过程中，内腔半径r为5mm的倍数。
CN201610976784.0A 2016-11-07 2016-11-07 一种爆热弹爆炸容腔半径的计算方法 Active CN106529016B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610976784.0A CN106529016B (zh) 2016-11-07 2016-11-07 一种爆热弹爆炸容腔半径的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610976784.0A CN106529016B (zh) 2016-11-07 2016-11-07 一种爆热弹爆炸容腔半径的计算方法

Publications (2)

Publication Number Publication Date
CN106529016A true CN106529016A (zh) 2017-03-22
CN106529016B CN106529016B (zh) 2019-11-05

Family

ID=58350667

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610976784.0A Active CN106529016B (zh) 2016-11-07 2016-11-07 一种爆热弹爆炸容腔半径的计算方法

Country Status (1)

Country Link
CN (1) CN106529016B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106949824A (zh) * 2017-04-17 2017-07-14 中水东北勘测设计研究有限责任公司 水下淤泥爆破空腔半径电极阵列测试法
CN110287635A (zh) * 2019-07-03 2019-09-27 合肥工业大学 一种同步配重盖式防爆门设计方法
CN112364425A (zh) * 2020-11-19 2021-02-12 中北大学 含阻尼刚性结构抗线性爆炸荷载设计动力系数方法
CN113158489A (zh) * 2021-05-13 2021-07-23 中国石油大学(华东) 一种基于等效载荷的抗爆管道壁厚的计算方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104931334A (zh) * 2015-06-01 2015-09-23 武汉大学 一种基于爆炸自毁作用的加卸载装置
CN105160144A (zh) * 2015-06-18 2015-12-16 武汉科技大学 一种适用于软介质爆破鼓包与爆腔运动的监测方法
CN105202362A (zh) * 2015-09-24 2015-12-30 中国工程物理研究院化工材料研究所 爆热弹密封结构

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104931334A (zh) * 2015-06-01 2015-09-23 武汉大学 一种基于爆炸自毁作用的加卸载装置
CN105160144A (zh) * 2015-06-18 2015-12-16 武汉科技大学 一种适用于软介质爆破鼓包与爆腔运动的监测方法
CN105202362A (zh) * 2015-09-24 2015-12-30 中国工程物理研究院化工材料研究所 爆热弹密封结构

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
冯晓军 等: "200g TNT当量爆热弹体的设计", 《万方数据库》 *
龙建华 等: "1OOg(TNT) 当量真空密封爆炸容器的设计", 《设计与研究》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106949824A (zh) * 2017-04-17 2017-07-14 中水东北勘测设计研究有限责任公司 水下淤泥爆破空腔半径电极阵列测试法
CN110287635A (zh) * 2019-07-03 2019-09-27 合肥工业大学 一种同步配重盖式防爆门设计方法
CN110287635B (zh) * 2019-07-03 2022-10-18 合肥工业大学 一种同步配重盖式防爆门设计方法
CN112364425A (zh) * 2020-11-19 2021-02-12 中北大学 含阻尼刚性结构抗线性爆炸荷载设计动力系数方法
CN112364425B (zh) * 2020-11-19 2023-07-07 中北大学 含阻尼刚性结构抗线性爆炸荷载设计动力系数方法
CN113158489A (zh) * 2021-05-13 2021-07-23 中国石油大学(华东) 一种基于等效载荷的抗爆管道壁厚的计算方法

Also Published As

Publication number Publication date
CN106529016B (zh) 2019-11-05

Similar Documents

Publication Publication Date Title
CN106529016A (zh) 一种爆热弹爆炸容腔半径的计算方法
Turner et al. Underwater implosion of cylindrical metal tubes
CN108280268B (zh) 用于检验水下爆破器材性能的压力容器参数设计方法
Cheng et al. Detonation characteristics of emulsion explosives sensitized by MgH 2
US8707764B1 (en) Assembly and method for standardized insensitive munitions testing
Xiang et al. Effect of Al/O ratio on the detonation performance and underwater explosion of HMX‐based aluminized explosives
CN107632043B (zh) 一种炸药弱冲击点火阈值测试装置和方法
Ren et al. Experimental study on dynamic buckling of submerged grid-stiffened cylindrical shells under intermediate-velocity impact
CN115221814A (zh) 一种氢气管道化学爆炸冲击波强度预测方法
CN109916586A (zh) 一种推进剂贮箱爆炸碎片初始速度的计算方法
Wang et al. Detonation characteristics of emulsion explosives sensitized by hydrogen‐storage glass microballoons
CN111191382A (zh) 一种计算内爆炸下金属管道前后向裂纹渐进扩展长度的方法
CN108426498B (zh) 用于模拟离心机高g值试验的爆炸装置起爆可靠性检验装置
Nowak et al. Small-scale underwater explosion in shallow-water tank
CN109086538A (zh) 一种水下用高冲击波能炸药能量设计方法
CN106248732A (zh) 一种炸药温度减威度的测量装置及表征方法
Shimamura et al. Study of water entry of high-špeed projectile
Iyama et al. Numerical simulation on manufacturing of pressure vessel for shock food processing using explosive forming
Chen et al. A numerical study on the blast wave distribution and propagation characteristics of cylindrical explosive in motion
Cai et al. Simulation study on the influence of liquid level on anti-jet penetration ability of single-cell composite structure
YUE et al. Experimental and simulation study of afterburning effect for blast load in confined cabin
ZHANG et al. Numerical simulation of dynamic response of laminated glass subjected to blast load
Schleyer et al. Experimental investigation of large plastic deformation and fracture in explosively loaded open-ended steel cylinders
FAN et al. Simulation and experimental study of tantalum liner to form dual-mode damage element by detonation
Dhouibi et al. Modeling and simulation of the engraving process in different life stages of small caliber guns

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20210927

Address after: No.1, kekecheng Avenue, Youxian District, Mianyang City, Sichuan Province

Patentee after: SICHUAN HAOKANG TECHNOLOGY Co.,Ltd.

Address before: Mianyang City, Sichuan Province, 621000 Mountain Road No. 64

Patentee before: INSTITUTE OF CHEMICAL MATERIALS OF CHINA ACADEMY OF ENGINEERING AND PHYSICS