CN110110382B - 一种内部双爆源烟花厂房结构爆炸荷载的确定方法 - Google Patents
一种内部双爆源烟花厂房结构爆炸荷载的确定方法 Download PDFInfo
- Publication number
- CN110110382B CN110110382B CN201910289348.XA CN201910289348A CN110110382B CN 110110382 B CN110110382 B CN 110110382B CN 201910289348 A CN201910289348 A CN 201910289348A CN 110110382 B CN110110382 B CN 110110382B
- Authority
- CN
- China
- Prior art keywords
- explosion
- wall surface
- factory building
- areas
- wall
- 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/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Architecture (AREA)
- Civil Engineering (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Description
技术领域
本发明是一种烟花厂房结构爆炸荷载的确定方法,尤其是涉及一种内部双爆源烟花厂房结构爆炸荷载的确定方法。
背景技术
爆炸荷载常采用的确定方法有两种,一种是参照规范进行确定;另一种是采用显式动力分析软件进行数值模拟确定。上述两种方法在应用于内部双爆源烟花厂房的爆炸荷载确定时,存在诸多不足。根据相关规范,单爆源室内爆炸在确定爆炸荷载时,假定爆炸荷载沿结构各墙面均匀分布,然而双爆源室内爆炸的爆炸荷载与单爆源对比,沿墙面分布极不均匀,因此在等效为单爆源进行荷载确定时,存在很大的不准确性,在实际爆炸荷载的局部峰值区域存在很大的安全隐患。而采用显式动力分析软件可以准确的计算出其爆炸荷载,但流固耦合算法耗时极长,软件操作对非专业人员而言要求高。因此以上两种方法在确定内部双爆源爆炸荷载时缺乏针对性,难以应用于工程实际。
发明内容
本发明的目的就是为了克服上述现有技术的缺陷,提供一种可以准确、快速的确定内部双爆源烟花厂房结构的爆炸荷载的内部双爆源烟花厂房结构爆炸荷载的确定方法。
本发明所采用的技术方案是:
一种内部双爆源烟花厂房结构爆炸荷载的确定方法,包括如下步骤:
(1)烟花厂房结构为典型的三面墙和一面屋顶的钢筋混凝土结构,三面墙分别为两面侧墙和一面后墙,将厂房各墙面进行如下划分:
对于左右两面侧墙:将左右两面侧墙等分成九个相等的区域;
对于后墙和屋顶:以后墙和屋顶的中轴线为对称轴将后墙和屋顶分为左右两部分区域,将左右两部分区域分别分成九个相等的区域;
(2)根据规范UFC 3-340-02计算出厂房内两个间隔距离为20cm-240cm的爆源分别在不同间距下同时爆炸分别作用在烟花厂房结构内部的爆炸荷载,所述的爆炸荷载包括厂房内各墙面的平均反射压力和平均反射冲量/>然后将两个爆源同时作用在一个墙面上的平均反射压力/>相加得到该墙面的总平均反射压力,将两个爆源同时作用在一个墙面上的平均反射冲量/>分别相加得到该墙面的总平均反射冲量,两爆源均位于结构中心竖轴一上一下并列排布,关于结构中心点对称;
(3)将厂房在某一爆源间隔距离下各墙面的总平均反射压力和总平均反射冲量分别乘以该间隔距离下各墙面各区域的爆炸荷载增强系数λ,得到烟花厂房结构内部双爆源在该间隔距离下同时爆炸作用下各墙面各区域调整后的爆炸荷载;
所述的爆炸荷载增强系数λ是通过数值模拟得到,具体步骤如下:
第一步,采用LS-DYNA显式动力分析软件对内部双爆源烟花厂房的爆炸进行数值模拟,建立炸药和空气模型;采用单点多物质ALE算法进行模拟,计算出内部双爆源间隔距离为20cm时烟花厂房结构各墙面的反射压力和反射冲量;
第二步,采用多点积分方法,在LS-DYNA烟花厂房结构模型的左右两面侧墙的九个区域中以及后墙和屋顶对称轴一侧的九个区域中分别均布选取25个测点,将每个区域中的25个测点的反射压力和反射冲量分别进行多点积分,求出各墙面各区域的反射压力和反射冲量,其中后墙和屋顶对称轴另一侧的九个区域的反射压力和反射冲量与后墙和屋顶对称轴一侧的九个区域反射压力和反射冲量关于对称轴对称;
第三步,将计算出的各墙面各区域反射压力Pr和反射冲量ir,分别除以各墙面的总平均反射压力和总平均反射冲量,得到各墙面九个区域的爆炸荷载增强系数λ,当λ<1时,取λ=1;
第四步,重复第一步、第二步、第三步,在重复的过程中只改变第一步中构建LS-DYNA有限元模型时两爆源的间隔距离,将两爆源的间隔距离由20cm开始递增至240cm,每次递增10cm,求出每个爆源间隔距离下的后墙中间水平排列的六个区域中对称轴一侧的三个区域以及侧墙墙面中间水平排列的三个区域的爆炸荷载增强系数,每个爆源间隔距离下的后墙中间水平排列的六个区域中对称轴另一侧的三个区域的爆炸荷载增强系数关于对称轴对称,绘制后墙和侧墙的增强系数λ随两爆源的间隔距离变化的关系图。
本发明的一种内部双爆源烟花厂房结构爆炸荷载的确定方法,解决了传统方法难以应用于内部双爆源烟花厂房结构爆炸荷载确定的缺陷,具有以下优点:
1、考虑了爆炸荷载的分布不均匀性,可以准确的确定内部双爆源烟花厂房结构的爆炸荷载;
2、步骤简单,不需要使用高性能计算机,对使用人员的技术要求低,可以快速的确定内部双爆源烟花厂房结构的爆炸荷载。
附图说明
图1-1为本发明中烟花厂房结构侧墙九个区域分区图;
图1-2为本发明中烟花厂房结构后墙十八个区域分区图;
图1-3为本发明中烟花厂房结构屋顶十八个区域分区图;
图2-1为本发明中所公开烟花厂房结构爆源间隔20cm后墙各区域增强系数;
图2-2为本发明中所公开烟花厂房结构爆源间隔20cm屋顶各区域增强系数;
图2-3为本发明中所公开烟花厂房结构爆源间隔20cm侧墙各区域增强系数;
图3-1为本发明中后墙反射压力的增强系数随两爆源间隔距离变化图(间隔距离变化范围为20cm-240cm);
图3-2为本发明中后墙反射冲量的增强系数随两爆源间隔距离变化图(间隔距离变化范围为20cm-240cm);
图4-1为本发明中侧墙反射压力的增强系数随两爆源间隔距离变化图(间隔距离变化范围为20cm-240cm);
图4-2为本发明中侧墙反射冲量的增强系数随两爆源间隔距离变化图(间隔距离变化范围为20cm-240cm);
图5为实施例中通过本发明的方法计算的烟花厂房结构后墙各区域的爆炸荷载数值(反射压力Pr和作用时间T0)。
具体实施方式
下面结合实施例和附图对本发明做出详细说明。
本发明的一种内部双爆源烟花厂房结构爆炸荷载的确定方法,包括如下步骤:
(1)烟花厂房结构为典型的三面墙和一面屋顶的钢筋混凝土结构,三面墙分别为两面侧墙和一面后墙。将厂房各墙面进行如下划分:
对于左右两面侧墙:将左右两面侧墙等分成九个相等的区域(如图所示用罗马数字表示为Ⅰ、Ⅱ、Ⅲ、、、Ⅶ、Ⅷ、Ⅸ九个区域)。
对于后墙和屋顶:以后墙和屋顶的中轴线为对称轴将后墙和屋顶分为左右两部分区域,将左右两部分区域分别分成九个相等的区域(如图所示用罗马数字表示为Ⅰ、Ⅱ、Ⅲ、、、Ⅶ、Ⅷ、Ⅸ九个区域)。
(2)根据规范UFC 3-340-02计算出厂房内两个间隔距离为20cm-240cm的爆源同时爆炸分别作用在烟花厂房结构内部的爆炸荷载,所述的爆炸荷载包括厂房内各墙面的平均反射压力和平均反射冲量/>然后将两个爆源同时作用在一个墙面上的平均反射压力/>相加得到该墙面的总平均反射压力,将两个爆源同时作用在一个墙面上的平均反射冲量/>分别相加得到该墙面的总平均反射冲量,两爆源均位于结构中心竖轴一上一下并列排布,关于结构中心点对称;
(3)将厂房在某一爆源间隔距离下各墙面的总平均反射压力和总平均反射冲量分别乘以该间隔距离下各墙面各区域的爆炸荷载增强系数λ,得到烟花厂房结构内部双爆源在该间隔距离下同时爆炸作用下各墙面各区域调整后的爆炸荷载。
所述的爆炸荷载增强系数λ是通过数值模拟得到,具体步骤如下:
第一步,采用LS-DYNA显式动力分析软件对内部双爆源烟花厂房的爆炸进行数值模拟,建立炸药和空气模型;采用单点多物质ALE算法进行模拟,计算出内部双爆源间隔距离为20cm时烟花厂房结构各墙面的反射压力和反射冲量;
第二步,采用多点积分方法,在LS-DYNA烟花厂房结构模型的左右两面侧墙的九个区域中以及后墙和屋顶对称轴一侧的九个区域中分别均布选取25个测点,将每个区域中的25个测点的反射压力和反射冲量分别进行多点积分,求出各墙面各区域的反射压力和反射冲量,其中后墙和屋顶对称轴另一侧的九个区域的反射压力和反射冲量与后墙和屋顶对称轴一侧的九个区域反射压力和反射冲量关于对称轴对称;
第三步,将计算出的各墙面各区域反射压力Pr和反射冲量ir,分别除以参照规范(UFC 3-340-02)中计算各墙面的总平均反射压力和总平均反射冲量,得到各墙面九个区域的爆炸荷载增强系数λ,即本发明所公开的结构各墙面各区域的爆炸荷载增强系数λ,当λ<1时,取λ=1。
第四步,重复第一步、第二步、第三步,在重复的过程中只改变第一步中构建LS-DYNA有限元模型时两爆源的间隔距离,将两爆源的间隔距离由20cm开始递增至240cm,每次递增10cm,求出每个爆源间隔距离下的后墙中间水平排列的六个区域中对称轴一侧的三个区域(对称轴另一侧的三个区域的爆炸荷载增强系数关于对称轴对称)以及侧墙墙面中间水平排列的三个区域(即Ⅳ、Ⅴ、Ⅵ区域)的爆炸荷载增强系数,绘制后墙和侧墙的增强系数λ随两爆源的间隔距离变化的关系图。
两个爆源之间的间隔距离发生变化时,后墙中间水平排列的六个区域和侧墙墙面中间水平排列的三个区域(即Ⅳ、Ⅴ、Ⅵ区域)的增强系数λ会受此影响发生变化,其他区域受间隔距离变化影响小,增强系数不变。
本发明的一种内部双爆源烟花厂房结构爆炸荷载的确定方法,解决了传统方法难以应用于内部双爆源烟花厂房结构爆炸荷载确定的缺陷,具有以下优点:
1、考虑了爆炸荷载的分布不均匀性,可以准确的确定内部双爆源烟花厂房结构的爆炸荷载;
2、步骤简单,不需要使用高性能计算机,对使用人员的技术要求低,可以快速的确定内部双爆源烟花厂房结构的爆炸荷载。
实施例1
烟花厂房结构尺寸为6000mm*3000mm*3000mm,两爆源均位于结构中心轴一上一下并列排布,关于结构中心点对称,间隔距离为150cm,两爆源同时起爆。
1)参照规范UFC 3-340-02,确定烟花厂房结构所有关键参数如表1、表2。
表1烟花厂房参数值
表2烟花厂房关键参数值
其中N为反射面的个数;h为爆源至底面的距离;H为结构边墙高度;l为爆源至边墙的距离;L为结构底面长度;W为炸药重量;RA为爆源至后墙的距离;表中的ZA表示比例距离ZA=RA/W1/3;表中的ft表示所采用的长度单位英尺,表中的lb表示所采用的重量单位磅;表中的ft/lb1/3即表示英尺/磅1/3;表2为表1中各参数的比值和比例距离ZA;
爆源O1爆炸时作用在后墙墙面的爆炸荷载:
由于结构为对称结构,所以爆源O2爆炸时作用在后墙墙面的爆炸荷载与爆源O1相等。
2)将爆源O1和爆源O2爆炸后同时作用在后墙墙面的平均反射压力和平均反射冲量/>分别相加后,得到后墙墙面的总平均反射压力和总平均反射冲量,根据图2-1确定后墙墙面相应区域的增强系数λ。由于间隔距离为150cm时,后墙中间水平排列的六个区域的增强系数发生改变,该六个区域的增强系数需要通过图3-1、图3-2进行确定,其他区域的增强系数不变(按照图2-1取值)。将后墙墙面的总平均反射压力和总平均反射冲量分别乘以后墙墙面各区域的增强系数λ,得到内部双爆源烟花厂房结构后墙墙面各区域的爆炸荷载(反射压力/>和反射冲量/>为了更直观的反应爆炸荷载的作用情况,通过反射冲量可计算得出反射压力的作用时间T0),如图5。
以后墙墙面Ⅴ区域的爆炸荷载计算为例:
Pr=1.834×2×λpr=1.834×2×1.000=3.668MPa
ir=1.1×2×λir=1.1×2×1.438=3.1636MPa·ms
上述对实施例的详细描述是为便于该技术领域的普通技术人员能理解和使用发明。熟悉本领域的人员可以对这些实施例做出各种修改,并把在此说明的一般原理应用到其他实施例中而不必经过创造性的劳动。因此,本发明不限于上述实施例,本领域技术人员在本发明的范畴内所作出的改进和修改都应该在本发明的保护范围之内。
Claims (1)
1.一种内部双爆源烟花厂房结构爆炸荷载的确定方法,包括如下步骤:
(1)烟花厂房结构为典型的三面墙和一面屋顶的钢筋混凝土结构,三面墙分别为两面侧墙和一面后墙,将厂房各墙面进行如下划分:
对于左右两面侧墙:将左右两面侧墙等分成九个相等的区域;
对于后墙和屋顶:以后墙和屋顶的中轴线为对称轴将后墙和屋顶分为左右两部分区域,将左右两部分区域分别分成九个相等的区域;
(2)根据规范UFC3-340-02计算出厂房内两个间隔距离为20cm-240cm的爆源分别在不同间距下同时爆炸分别作用在烟花厂房结构内部的爆炸荷载,所述的爆炸荷载包括厂房内各墙面的平均反射压力和平均反射冲量/>然后将两个爆源同时作用在一个墙面上的平均反射压力/>相加得到该墙面的总平均反射压力,将两个爆源同时作用在一个墙面上的平均反射冲量/>分别相加得到该墙面的总平均反射冲量,两爆源均位于结构中心竖轴一上一下并列排布,关于结构中心点对称;
(3)将厂房在某一爆源间隔距离下各墙面在的总平均反射压力和总平均反射冲量分别乘以该间隔距离下各墙面各区域的爆炸荷载增强系数λ,得到烟花厂房结构内部双爆源在该间隔距离下同时爆炸作用下各墙面各区域调整后的爆炸荷载;
所述的爆炸荷载增强系数λ是通过数值模拟得到,具体步骤如下:
第一步,采用LS-DYNA显式动力分析软件对内部双爆源烟花厂房的爆炸进行数值模拟,建立炸药和空气模型;采用单点多物质ALE算法进行模拟,计算出内部双爆源间隔距离为20cm时烟花厂房结构各墙面的反射压力和反射冲量;
第二步,采用多点积分方法,在LS-DYNA烟花厂房结构模型的左右两面侧墙的九个区域中以及后墙和屋顶对称轴一侧的九个区域中分别均布选取25个测点,将每个区域中的25个测点的反射压力和反射冲量分别进行多点积分,求出各墙面各区域的反射压力和反射冲量,其中后墙和屋顶对称轴另一侧的九个区域的反射压力和反射冲量与后墙和屋顶对称轴一侧的九个区域反射压力和反射冲量关于对称轴对称;
第三步,将计算出的各墙面各区域反射压力Pr和反射冲量ir,分别除以各墙面的总平均反射压力和总平均反射冲量,得到各墙面九个区域的爆炸荷载增强系数λ,当λ<1时,取λ=1;
第四步,重复第一步、第二步、第三步,在重复的过程中只改变第一步中构建LS-DYNA有限元模型时两爆源的间隔距离,将两爆源的间隔距离由20cm开始递增至240cm,每次递增10cm,求出每个爆源间隔距离下的后墙中间水平排列的六个区域中对称轴一侧的三个区域以及侧墙墙面中间水平排列的三个区域的爆炸荷载增强系数,每个爆源间隔距离下的后墙中间水平排列的六个区域中对称轴另一侧的三个区域的爆炸荷载增强系数关于对称轴对称,绘制后墙和侧墙的增强系数λ随两爆源的间隔距离变化的关系图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910289348.XA CN110110382B (zh) | 2019-04-11 | 2019-04-11 | 一种内部双爆源烟花厂房结构爆炸荷载的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910289348.XA CN110110382B (zh) | 2019-04-11 | 2019-04-11 | 一种内部双爆源烟花厂房结构爆炸荷载的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110110382A CN110110382A (zh) | 2019-08-09 |
CN110110382B true CN110110382B (zh) | 2023-06-09 |
Family
ID=67484083
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910289348.XA Active CN110110382B (zh) | 2019-04-11 | 2019-04-11 | 一种内部双爆源烟花厂房结构爆炸荷载的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110110382B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106055832A (zh) * | 2016-06-22 | 2016-10-26 | 国网上海市电力公司 | 一种变电站防爆安全性评估方法 |
CN109376459A (zh) * | 2018-11-14 | 2019-02-22 | 中北大学 | 一种在爆炸荷载作用下等截面固端梁荷载系数的计算方法 |
CN109580193A (zh) * | 2018-11-14 | 2019-04-05 | 中北大学 | 一种在爆炸荷载作用下等截面简支梁荷载系数的计算方法 |
-
2019
- 2019-04-11 CN CN201910289348.XA patent/CN110110382B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106055832A (zh) * | 2016-06-22 | 2016-10-26 | 国网上海市电力公司 | 一种变电站防爆安全性评估方法 |
CN109376459A (zh) * | 2018-11-14 | 2019-02-22 | 中北大学 | 一种在爆炸荷载作用下等截面固端梁荷载系数的计算方法 |
CN109580193A (zh) * | 2018-11-14 | 2019-04-05 | 中北大学 | 一种在爆炸荷载作用下等截面简支梁荷载系数的计算方法 |
Non-Patent Citations (2)
Title |
---|
cui j.Study of concrete damage mechanism under hydrostatic pressure numerical simulations.《web of science》.2018,全文. * |
丁阳 ; 陈晔 ; 师燕超 ; .室内爆炸超压荷载简化模型.工程力学.2015,(第03期),全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN110110382A (zh) | 2019-08-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Montuori et al. | Geometrical patterns for diagrid buildings: Exploring alternative design strategies from the structural point of view | |
CN107391804B (zh) | 基于综合造价法的高层建筑结构抗震性能优化方法 | |
CN104242306B (zh) | 一种基于主成分分析方法的电力系统自适应分区方法 | |
JPH102004A (ja) | 住宅の水平ブレース配置方法 | |
CN112528540B (zh) | 一种基于力热耦合计算的结构优化方法 | |
CN109918712A (zh) | 一种基于p型有限元法和围线积分法求解复合型应力强度因子的计算方法 | |
Panton et al. | Correlation of pressure fluctuations in turbulent wall layers | |
CN110110382B (zh) | 一种内部双爆源烟花厂房结构爆炸荷载的确定方法 | |
CN106202666A (zh) | 一种船舶轴系轴承位移调整的计算方法 | |
WO2018014454A1 (zh) | 一种多点吊装计算方法 | |
Acosta et al. | Towards an analysis of the performance of lightwell skylights under overcast sky conditions | |
CN109580193A (zh) | 一种在爆炸荷载作用下等截面简支梁荷载系数的计算方法 | |
CN110378011A (zh) | 一种金属波纹管液压胀形工艺稳健设计方法 | |
CN107655650B (zh) | 用于结构试验的杠杆加载系统碰撞检测方法 | |
CN104947838B (zh) | 一种用于寻找索杆结构可行预应力的方法 | |
CN106202759A (zh) | 一种可编程仿生负泊松比结构及设计方法 | |
Mohammadnejad et al. | Dynamic response analysis of a combined system of framed tubed, shear core and outrigger-belt truss | |
Feng et al. | Multi-objective morphology optimization of free-form cable-braced grid shells | |
CN104050342B (zh) | 变电站软导线装配仿真优化设计系统 | |
CN104008253B (zh) | 一种伸缩吊臂动态模型的集成优化方法 | |
CN109164852A (zh) | 一种大体积混凝土自适应温度与应力控制方法 | |
CN112989475B (zh) | 一种适用于中国抗震设计规范的位移设计谱方法 | |
CN107103126B (zh) | 一种铆接结构干涉量快速预测方法 | |
CN108920787B (zh) | 一种基于自适应配点的结构模糊不确定性分析方法 | |
Tarabieh et al. | Parametric investigation of three types of brick bonds for thermal performance in a hot arid climate zone |
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 |