CN113158489B - 一种基于等效载荷的抗爆管道壁厚的计算方法 - Google Patents
一种基于等效载荷的抗爆管道壁厚的计算方法 Download PDFInfo
- Publication number
- CN113158489B CN113158489B CN202110522262.4A CN202110522262A CN113158489B CN 113158489 B CN113158489 B CN 113158489B CN 202110522262 A CN202110522262 A CN 202110522262A CN 113158489 B CN113158489 B CN 113158489B
- Authority
- CN
- China
- Prior art keywords
- pipeline
- gas
- detonation
- wall thickness
- calculating
- 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.)
- Active
Links
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
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)
- Pipeline Systems (AREA)
- Shaping Metal By Deep-Drawing, Or The Like (AREA)
Abstract
本发明公开一种基于等效载荷的抗爆管道壁厚的计算方法,该方法首先根据管道内可燃气体类型,以及可能混入的氧化气体种类,考虑最危险情况,即管道内可燃气体与氧化气体的混合比例为化学计量比时,计算管道内混合气体发生爆轰的爆轰压力,以爆轰压力的k倍作为等效静态压力,然后根据公式初步计算管道壁厚,再向上取整得到最终的管道壁厚。本发明的方法充分考虑了管道的服役安装条件以及内部可能出现的不同***场景,可更有针对性地满足管道设计需求,使所设计管道满足抗爆性能的同时具有更合适的壁厚,从而具有更好的经济性,且该方法实施方便、迅速,对人员专业程度要求低。
Description
技术领域
本发明涉及管道***预测领域,具体涉及一种基于等效载荷的抗爆管道壁厚的计算方法。
背景技术
在国防、化工、核电、公共安全等领域,某些圆柱壳结构或管道设计时需满足抗爆性能的要求,即内部发生***结构也不会失效,达到其本质安全的设计要求。抗爆容器或管道设计中的一个关键是壁厚的计算与确定,目前主要有两类方法。一类是理论或经验方法,如等效单自由度方法。该类方法的基本思路是计算***载荷的等效静态压力,之后按照常规容器的设计方法计算壁厚,但该类方法主要适用于***物是TNT的情况,对于常见的可燃气体管道,管道内可能发生爆燃、爆轰以及爆燃转爆轰等不同***场景,不同***场景对应的等效静态压力是不同的,而这是目前方法所没有考虑的。第二类是数值模拟类方法,该类方法需对所设计管道进行弹塑性有限元动力学分析,根据分析结果判断结构是否能抵抗相应的***载荷,若不满足,则增加壁厚重复进行分析。该类方法对实施人员的专业水平要求较高,且数值分析结果受本构模型选取、网格尺寸等影响,在具体实施前往往需要对模型进行验证,实施周期长,难度较高。
发明内容
针对现有技术的不足,本发明提供一种基于等效载荷的抗爆管道壁厚的计算方法,该方法更有针对性地满足管道设计需求,且计算快捷方便。
本发明的目的通过如下的技术方案来实现:
一种基于等效载荷的抗爆管道壁厚的计算方法,该方法具体包括如下步骤:
(1)根据管道内可燃气体类型,以及可能混入的氧化气体种类,考虑最危险情况,即管道内可燃气体与氧化气体的混合比例为化学计量比时,计算管道内混合气体发生爆轰的爆轰压力pCJ;
(2)取爆轰压力pCJ的k倍作为等效静态压力peq,即
peq=k·pCJ
其中,k在不同的场景下取值不同;
(3)根据下式初步计算管道壁厚δ
(4)将初步计算的管道壁厚δ叠加上管道腐蚀余量,并根据管道设计国标向上取值,得到最终的管道壁厚δn。
进一步地,所述爆轰压力pCJ通过以下两种方法得到:
(1)通过如下的理论公式计算得到
式中,Q为管道内单位质量混合气体发生化学反应所释放的能量,根据反应气体的燃烧热计算;γ为爆轰产物的绝热指数,可根据爆轰产物气体种类及体积百分数计算得到,γ取值范围为1.1~1.4;为管道内初始混合气体的平均摩尔质量;R为理想气体常数;T为管道内混合气体温度,可查询管道运行工艺参数得到,或取环境温度。
(2)采用程序CEARUN计算。
进一步地,所述k取值如下:
当所设计管道两端连接其他管道或开口,且仅承受气体稳态爆轰载荷时,k=0.70;
当所设计管道两端连接其他管道或开口,且可能发生气体爆燃转爆轰时,k=5.25;
当所设计管道一端封闭,且仅承受稳态气体爆轰载荷,k=1.68;
当所设计管道一端封闭,且可能发生气体爆燃转爆轰时,k=12.6。
本发明的有益效果如下:
(1)本发明充分考虑了管道的服役安装条件以及内部可能出现的不同***场景,对应给出了不同的等效静态压力计算缩放系数,可更有针对性地满足管道设计需求,使所设计管道满足抗爆性能的同时具有更合适的壁厚,从而具有更好的经济性;
(2)本发明中壁厚设计主要依据公式进行计算,且步骤清晰、参数明确,相比数值模拟类方法实施方便、迅速,对人员专业程度要求低。
附图说明
图1为本发明的基于等效载荷的抗爆管道壁厚的计算方法的流程图;
图2为实施例中采用程序CEARUN计算爆轰压力的示意图。
具体实施方式
下面根据附图和优选实施例详细描述本发明,本发明的目的和效果将变得更加明白,应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
如图1所示,本发明的基于等效载荷的抗爆管道壁厚的计算方法,具体包括如下四部分。
1.明确管道设计要求与基本参数
根据管道设计要求,明确管道内气体类型,获取管道设计内径Di,获取管道运行压力p,根据管道材质查询该材料的许用应力[σ]和屈服强度σy,通常[σ]=σy/n,n为安全系数。
2.计算管道内气体爆轰压力
根据管道内可燃气体类型,以及可能混入的氧化气体种类(如空气、氧气),考虑最危险情况,即管道内可燃气体与氧化气体的混合比例为化学计量比,计算管道内混合气体发生爆轰的爆轰压力pCJ,具体可采用以下两种方法计算:
①理论公式法
采用下式计算管道内混合气体的爆轰压力pCJ:
式中,Q为管道内单位质量混合气体发生化学反应所释放的能量,可根据反应气体的燃烧热计算;γ为爆轰产物(气体)的绝热指数,可根据爆轰产物气体种类及体积百分数计算得到,常见气体的取值范围为:1.1-1.4;为管道内初始混合气体的平均摩尔质量;R为理想气体常数;T为管道内混合气体温度,可查询管道运行工艺参数得到,或取环境温度。
②采用程序CEARUN计算
程序CEARUN可提供更为精确的气体爆轰压力计算结果,且操作简便,具有网页版用户操作界面,得到领域研究人员较广泛的使用,具体网址为:https://cearun.grc.nasa.gov/index.html。根据程序说明输入混合气体种类、体积百分数,以及气体初始压力、温度等参数,即可计算获得爆轰压力pCJ。
3.计算管道设计等效静态压力
定义等效静态压力为:管道在该压力作用下的塑性变形与相应气体***场景导致的管道塑性变形一致,大小等于k倍的气体爆轰压力,即:
peq=k·pCJ
k值按照以下原则确定:
①所设计管道两端通过法兰连接其他管道或开口,且仅承受气体稳态爆轰载荷:
该场景下,爆轰波不会在所设计管道内发生轴向反射,k取值0.70。
②所设计管道两端通过法兰连接其他管道或开口,且可能发生气体爆燃转爆轰(DDT,Deflagration to Detonation Transition):
该场景下,管道内最大动压由DDT产生,且可能叠加初始爆燃波压力,k取值5.25。
③所设计管道一端连接有法兰盲板或封闭,且仅承受稳态气体爆轰载荷:
该场景下,爆轰波在管道内轴向传播中可能发生反射,形成更大的动压力,此时k取值1.68。
④所设计管道一端连接有法兰盲板或封闭,且可能发生气体爆燃转爆轰:
考虑该场景下最危险情形,即管道内气体DDT与爆轰波轴向反射同时发生,并叠加了初始爆燃波压力,此时k取值12.6。
4.计算并确定管道名义壁厚
根据等效静态压力大小,选择合适的压力容器设计准则(如弹性失效准则、塑性失效准则、***失效准则等),计算管道壁厚δ:
①若peq<0.4[σ],则一般根据弹性失效设计准则和中径公式计算管道壁厚δ:
②若peq≥0.4[σ],则可根据塑性失效准则和拉美公式计算管道壁厚δ:
根据计算壁厚δ,叠加管道腐蚀余量并向上圆整确定管道最终名义厚度δn。
下面给出一个具体的实施例,进一步说明本发明的方法的计算过程和结果。
某一管道内径为40mm,材质为6061-T6铝合金,内部气体为按化学计量比混合的乙烯和氧气,气体初始压力为180kPa,环境温度为300K。该管道安装条件为一端通过法兰连接其他管道,另一端通过聚酯膜片密封。设计要求其能承受内部乙烯、氧气混合气体的稳态爆轰载荷,以便使用该管道开展相关实验研究。下面根据本发明对该抗爆管道壁厚进行设计:
步骤1:明确管道设计要求与基本参数
根据管道设计要求,管道内气体类型为乙烯(C2H4)和氧气(O2),管道设计内径Di=40mm,管道运行压力p=180kPa,查询6061-T6材料力学性能参数,其屈服强度σy=240MPa,取安全系数1.5,其许用应力[σ]=σy/1.5=160MPa。
步骤2:计算管道内气体爆轰压力
该管道内部本身为按化学计量比混合的乙烯(C2H4)和氧气(O2),即为最危险情况,气体温度为环境温度300K,初始压力为180kPa,采用以下两种方法计算气体爆轰压力:
①理论公式法
乙烯与氧气的化学反应方程式为:
C2H4+3O2=2CO2+2H2O
②采用程序CEARUN计算
打开网页https://cearun.grc.nasa.gov/index.html,选择问题类型为“det”(求解气体爆轰压力),输入气体温度为300K,初始压力为180kPa,选择燃烧物为乙烯C2H4,氧化物为氧气O2,设置两者的摩尔混合比为1:3,点击计算,程序求解的爆轰压力为6.10MPa,如图2所示。
步骤3:计算管道设计等效静态压力
该管道安装条件为一端通过法兰连接其他管道,另一端通过聚酯膜片密封。爆轰波作用下聚酯膜片会迅速破裂,不会导致爆轰波反射,即管道安装条件按照两端开口处理,设计要求其仅需承受气体稳态爆轰载荷,根据本发明中的规定,k取值0.70,选取CEARUN程序计算的气体爆轰压力(6.10MPa)为基准,则该管道的设计等效静态压力计算为:
peq=k·pCJ=0.70×6.10=4.27MPa
步骤4:计算并确定管道名义壁厚
该管道的设计等效静态压力peq小于0.4倍材料许用应力,则按照弹性失效准则计算管道壁厚为:
查询《GB/T 4436-2012铝及铝合金管材外形尺寸及允许偏差》国家标准,考虑介质腐蚀性极微,腐蚀余量取0,向上取值至0.75mm,则该管道的名义厚度最终设计为0.75mm。
对该管道开展气体爆轰加载实验,结果表明该管道未发生塑性变形,满足相应设计要求。
本领域普通技术人员可以理解,以上所述仅为发明的优选实例而已,并不用于限制发明,尽管参照前述实例对发明进行了详细的说明,对于本领域的技术人员来说,其依然可以对前述各实例记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在发明的精神和原则之内,所做的修改、等同替换等均应包含在发明的保护范围之内。
Claims (1)
1.一种基于等效载荷的抗爆管道壁厚的计算方法,其特征在于,该方法具体包括如下步骤:
(1)根据管道内可燃气体类型,以及可能混入的氧化气体种类,考虑最危险情况,即管道内可燃气体与氧化气体的混合比例为化学计量比时,计算管道内混合气体发生爆轰的爆轰压力pCJ;
(2)取爆轰压力pCJ的k倍作为等效静态压力peq,即
peq=k·pCJ
其中,k在不同的场景下取值不同;
(3)根据下式初步计算管道壁厚δ
其中,Di为管道设计内径;[σ]为管道材料的许用应力;σy为管道材料的屈服强度;
(4)将初步计算的管道壁厚δ叠加上管道腐蚀余量,并根据管道设计国标向上取值,得到最终的管道壁厚δn;
所述爆轰压力pCJ通过以下两种方法得到:
(1)通过如下的理论公式计算得到
式中,Q为管道内单位质量混合气体发生化学反应所释放的能量,根据反应气体的燃烧热计算;γ为爆轰产物的绝热指数,可根据爆轰产物气体种类及体积百分数计算得到,γ取值范围为1.1~1.4;为管道内初始混合气体的平均摩尔质量;R为理想气体常数;T为管道内混合气体温度,可查询管道运行工艺参数得到,或取环境温度;p为管道运行压力;
(2)采用程序CEARUN计算;
所述k取值如下:
当所设计管道两端连接其他管道或开口,且仅承受气体稳态爆轰载荷时,k=0.70;
当所设计管道两端连接其他管道或开口,且可能发生气体爆燃转爆轰时,k=5.25;
当所设计管道一端封闭,且仅承受稳态气体爆轰载荷,k=1.68;
当所设计管道一端封闭,且可能发生气体爆燃转爆轰时,k=12.6。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110522262.4A CN113158489B (zh) | 2021-05-13 | 2021-05-13 | 一种基于等效载荷的抗爆管道壁厚的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110522262.4A CN113158489B (zh) | 2021-05-13 | 2021-05-13 | 一种基于等效载荷的抗爆管道壁厚的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113158489A CN113158489A (zh) | 2021-07-23 |
CN113158489B true CN113158489B (zh) | 2022-08-26 |
Family
ID=76874849
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110522262.4A Active CN113158489B (zh) | 2021-05-13 | 2021-05-13 | 一种基于等效载荷的抗爆管道壁厚的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113158489B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115221814B (zh) * | 2022-08-03 | 2023-04-25 | 中国石油大学(华东) | 一种氢气管道化学***冲击波强度预测方法 |
CN115033942B (zh) * | 2022-08-10 | 2022-11-04 | 中国长江三峡集团有限公司 | 砂浆内衬壁厚设计方法、管道修复方法及装置 |
CN115994420A (zh) * | 2022-11-22 | 2023-04-21 | 中国电力工程顾问集团中南电力设计院有限公司 | 一种适用于超超临界机组的仪表导管壁厚计算方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106874547A (zh) * | 2017-01-08 | 2017-06-20 | 浙江大学 | 一种预测纤维全缠绕复合材料气瓶***压力的方法 |
CN108318654A (zh) * | 2017-12-29 | 2018-07-24 | 中国人民解放军陆军工程大学 | 用于模拟离心机高g值试验的***装置起爆状态检验方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI390145B (zh) * | 2008-03-31 | 2013-03-21 | Rohm & Haas | 用於爆燃壓力衰減之方法及設備 |
CN106529016B (zh) * | 2016-11-07 | 2019-11-05 | 中国工程物理研究院化工材料研究所 | 一种爆热弹***容腔半径的计算方法 |
CN108280268B (zh) * | 2017-12-29 | 2021-01-19 | 中国人民解放军陆军工程大学 | 用于检验水下***器材性能的压力容器参数设计方法 |
CN111191382B (zh) * | 2020-01-09 | 2024-05-24 | 中国石油大学(华东) | 一种计算内***下金属管道前后向裂纹渐进扩展长度的方法 |
CN111611680B (zh) * | 2020-03-30 | 2024-02-09 | 浙江大学 | 基于失效模式设计的内压钢制椭圆形封头制造方法 |
-
2021
- 2021-05-13 CN CN202110522262.4A patent/CN113158489B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106874547A (zh) * | 2017-01-08 | 2017-06-20 | 浙江大学 | 一种预测纤维全缠绕复合材料气瓶***压力的方法 |
CN108318654A (zh) * | 2017-12-29 | 2018-07-24 | 中国人民解放军陆军工程大学 | 用于模拟离心机高g值试验的***装置起爆状态检验方法 |
Also Published As
Publication number | Publication date |
---|---|
CN113158489A (zh) | 2021-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113158489B (zh) | 一种基于等效载荷的抗爆管道壁厚的计算方法 | |
CN108828008A (zh) | 一种基于冲击波超压和准静压的***内***威力评价方法 | |
Tang et al. | Effect of vent size on explosion overpressure and flame behavior during vented hydrogen–air mixture deflagrations | |
CN115221814B (zh) | 一种氢气管道化学***冲击波强度预测方法 | |
Xiang et al. | Underwater explosion performance of RDX/AP-based aluminized explosives | |
Du et al. | Dynamic fracture response of pre-flawed elbow pipe subjected to internal hydrogen-oxygen detonation | |
Jing et al. | Ignition sensitivity and explosion behaviors of micron-sized aluminum powder: Comparison between flake aluminum powder and spherical aluminum powder | |
Liu et al. | Effects of initial temperature on the deflagration characteristics and flame propagation behaviors of CH4 and its blends with C2H6, C2H4, CO, and H2 | |
Yu et al. | Autoignition and detonation development induced by temperature gradient in n-C7H16/air/H2O mixtures | |
Keshavarz et al. | A simple method for reliable estimation of the bubble energy in the underwater explosion | |
Stamps et al. | Hydrogen-air-diluent detonation study for nuclear reactor safety analyses | |
Li et al. | Fundamental combustion characteristics of lean and stoichiometric hydrogen laminar premixed flames diluted with nitrogen or carbon dioxide | |
Morii et al. | What connects ignition and deflagration?--On explosive transition of deflagration | |
Kuhl et al. | Numerical simulations of thermobaric explosions | |
Sochet et al. | Detonability of fuel-oxygen and fuel-air mixtures | |
Meybodi et al. | Gas detonation forming by a mixture of H2+ O2 detonation | |
Li et al. | Rational design of pin puller and its closed bomb test-based equivalent testing | |
Stewart | Transition to detonation in a model problem | |
Leishear | Fluid Transients Ignited the San Bruno Gas Pipeline Explosions | |
Shigarkanthi et al. | Experimental investigation and modeling of cycle by cycle variations in a gas fuelled SI engine | |
Fomin | Reduced chemical kinetic model of detonation combustion of one-and multi-fuel gaseous mixtures with air | |
Haghgoo et al. | Numerical Simulation of Triangular Plate Deformation Profile Under Gaseous Detonation Loading | |
CN118332231A (zh) | 一种战斗部冲击波超压分布预测及验证方法 | |
Kuznetsov et al. | Structural response of DN15-tubes under radiolysis gas detonation loads for BWR safety applications | |
Clutter et al. | Assessment of Gaseous Detonations in Piping Systems and Pipelines |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20210723 Assignee: Qingdao Zhongchuang Zhihui Electronic Technology Co.,Ltd. Assignor: CHINA University OF PETROLEUM (EAST CHINA) Contract record no.: X2022980024511 Denomination of invention: A Calculation Method of Anti explosion Pipe Wall Thickness Based on Equivalent Load Granted publication date: 20220826 License type: Common License Record date: 20221227 |