CN112464338A - 一种近场爆炸下简支钢筋混凝土梁弯曲响应预测方法 - Google Patents

一种近场爆炸下简支钢筋混凝土梁弯曲响应预测方法 Download PDF

Info

Publication number
CN112464338A
CN112464338A CN202011313265.9A CN202011313265A CN112464338A CN 112464338 A CN112464338 A CN 112464338A CN 202011313265 A CN202011313265 A CN 202011313265A CN 112464338 A CN112464338 A CN 112464338A
Authority
CN
China
Prior art keywords
calculating
reinforced concrete
concrete beam
load
predicting
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
CN202011313265.9A
Other languages
English (en)
Other versions
CN112464338B (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.)
Xian Modern Chemistry Research Institute
Original Assignee
Xian Modern Chemistry Research Institute
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 Xian Modern Chemistry Research Institute filed Critical Xian Modern Chemistry Research Institute
Priority to CN202011313265.9A priority Critical patent/CN112464338B/zh
Publication of CN112464338A publication Critical patent/CN112464338A/zh
Application granted granted Critical
Publication of CN112464338B publication Critical patent/CN112464338B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/08—Construction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • 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)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Human Resources & Organizations (AREA)
  • Computer Hardware Design (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Marketing (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • General Business, Economics & Management (AREA)
  • Computational Mathematics (AREA)
  • Structural Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Game Theory and Decision Science (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Civil Engineering (AREA)
  • Architecture (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种近场爆炸下简支钢筋混凝土梁弯曲响应预测方法，本方法可根据炸药TNT当量、炸药与目标的相对位置关系以及目标特征量预测结构的在近场爆炸下的弯曲响应历程。本发明首先根据炸药TNT当量及位置建立结构表面的冲击波压力峰值、比冲量分布模型，然后根据载荷特征计算弹性、塑性阶段的变形函数及等效转换系数，并建立结构在近场爆炸下的单自由度动力学方程，最后通过对方程差分求解可预测出结构的弯曲响应。

Description

一种近场爆炸下简支钢筋混凝土梁弯曲响应预测方法
技术领域
本发明属于毁伤评估领域，主要涉及一种近场爆炸下简支钢筋混凝土梁弯曲响应预测方法。
背景技术
等效单自由度方法是一种能够基于能量守恒原理将均布荷载下的连续梁式结构简化为集中载荷作用下的单自由度动力学系统的方法，该方法因计算方便、快捷，被广泛应用于钢筋混凝土梁在均布爆炸载荷下的弯曲响应预测。
Biggs教授于1964年出版的《结构动力学》一书中详细介绍了均布爆炸荷载下钢筋混凝土梁的等效单自由度模型，并给出了模型中诸多参数的计算方法。后来，美军联合出版的 TM5-1300手册中将该等效单自由度方法用于建筑抗核爆防护设计。在核爆冲击波作用下，建筑外表面的载荷往往是均匀的，且脉宽长、衰减慢。然而，当常规炸药距离建筑结构较近时，结构表面的冲击波压力往往是非均匀的，直接采用该单自由度方法预测结构响应往往会带来很大的误差。
为了能将单自由度方法应用于爆炸近场下的钢筋混凝土梁结构响应预测，国内外很多学者都开展了相关工作，如：汪维采用了载荷等效的方法，首先基于虚功原理将非均布载荷等效为均布载荷，然而再采用均布荷载的单自由度模型预测结构响应(Wang W,ZhangD,Lu F Y,et al.A new SDOF method of one-way reinforced concrete slab undernon-uniform blast loading. Structural Engineering&Mechanics,2013,46(5):595-613.)；Nagata通过引入载荷形状系数，假设冲击波载荷在结构表面指数分布，并对传统单自由度方法进行了修正，建立了载荷非均匀分布特征的单自由度方法(Nagata M,Beppu M,Ichino H,et al.Method for evaluating the displacement response of RC beamssubjected to close-in explosion using modified SDOF model. EngineeringStructures,2018,157:105-118.)。
然而目前采用的方法均存在一些不足：
第一：非均布载荷作用下，结构的挠曲线变形函数会发生改变，变形函数的差异会不仅会改变结构进入塑性状态的临界点，同时会影响到结构的抗力；而以上方法均未考虑到非均匀载荷下结构变形函数的差异；
第二：非均布载荷对结构弯曲抗力的影响同时与变形函数与载荷非均匀性相关，以上模型中均没有充分考虑；
第三：在单自由度模型迭代计算中，结构抗力与变形速率密切相关，动态变形速率会引起动态材料应变率，从而使得材料强度动态变化进而改变弯曲抗力，然而以上模型中均采用恒定的应变率强化系数；
第四：以上模型在计算质量转换系数、载荷转换系数时，未能充分考虑形函孔和非均布载荷特征。
发明内容
为了解决现有技术存在的问题，本发明提供了种近场爆炸下简支钢筋混凝土梁弯曲响应预测方法，包括如下步骤：
一种近场爆炸下简支钢筋混凝土梁弯曲响应预测方法，该方法包括如下步骤：
步骤1：计算炸药到结构表面的比例距离；
步骤2：拟合爆炸冲击波压力、比冲量分布函数；
步骤3：计算弹性、塑性阶段的结构变形函数；
步骤4：计算弹性、塑性阶段的载荷转换系数和质量转换系数；
步骤5：计算等效载荷时间历程；
步骤6：建立简支钢筋混凝土梁单自由度动力学方程；
步骤7：预测结构弯曲响应。
进一步地，步骤1：通过炸药爆炸当量及炸药爆心距确定炸药中心到钢筋混凝土梁结构表面的比例距离r(m/kg1/3)，具体采用如下公式：
r＝R1/W1/3
其中，W为炸药爆炸当量kg，R1为炸药爆心距m。
进一步地，步骤2采用经验公式或数值仿真，计算钢筋混凝土梁迎爆面的冲击波压力峰值分布P(x)以及比冲量分布I(x)，所述的步骤2采用经验公式或数值仿真，计算钢筋混凝土梁迎爆面的冲击波压力峰值分布P(x)以及比冲量分布I(x)，其中
Figure RE-GDA0002893486150000021
Figure RE-GDA0002893486150000022
或
Figure RE-GDA0002893486150000023
或其他形式，Pm为跨中位置的峰值反射压力，αp为压力分布系数，x为梁表面沿跨度方向的坐标，坐标原点位于跨中，L为梁的长度；
Figure RE-GDA0002893486150000024
或
Figure RE-GDA0002893486150000025
或其他形式，其中Im为跨中位置的峰值比冲量，αi为比冲量分布系数。
特别地步骤1，进一步计算压力载荷P(x)作用下结构的弹性挠曲线方程ye(x)、曲率方程 k(x)、弯矩方程M(x)，并根据结构截面尺寸及材料属性确定屈服曲率kc，并计算跨中屈服挠度yc：
Figure RE-GDA0002893486150000031
所述的步骤3:计算弹性变形函数
Figure RE-GDA0002893486150000032
和塑性变形函数
Figure RE-GDA0002893486150000033
Figure RE-GDA0002893486150000034
Figure RE-GDA0002893486150000035
x为梁表面沿跨度方向的坐标，坐标原点位于跨中，梁的长度为L。
所述的步骤4:计算弹性、塑性阶段的载荷转换系数和质量转换系数，其中弹性阶段载荷转换系数KLE和质量转换系数KME：
Figure RE-GDA0002893486150000036
Figure RE-GDA0002893486150000037
塑性阶段载荷转换系数KLP和质量转换系数KMP为：
Figure RE-GDA0002893486150000038
Figure RE-GDA0002893486150000039
所述的步骤5计算等效载荷时间历程，根据压力和比冲量算结构表面的总冲量Itotal，峰值载荷Fmax，载荷等效作用时间td，并确定等效载荷时间历程F(t)：
Figure RE-GDA00028934861500000310
Figure RE-GDA0002893486150000041
Figure RE-GDA0002893486150000042
Figure RE-GDA0002893486150000043
其中其中W为为炸药爆炸当量kg，w为结构宽度m。
所述的步骤6建立简支钢筋混凝土梁单自由度动力学方程如下：
Figure RE-GDA0002893486150000044
其中，KM为质量转换系数，KL为载荷转换系数，M为结构总质量，y、
Figure RE-GDA0002893486150000045
和
Figure RE-GDA0002893486150000046
分别为结构跨中挠度、速度和加速度，R为弯曲抗力。
所述的步骤7预测结构弯曲响应，即确定足够小的求解时间步长Δt，并初始化时间序列 ti＝iΔt，i＝0,1,2,3,…n；对加速度项进行中心差分得到迭代公式；
Figure RE-GDA0002893486150000047
构造
Figure RE-GDA0002893486150000048
y(t0)＝0；然后依次迭代计算y(ti+1)，i＝0,1,2,3,…n，直到 |y(ti+1)|＜|y(ti)|。
进一步地，本发明的步骤7迭代计算包括如下步骤：
第一步：根据y(ti)和y(ti-1)计算弯曲抗力R，包括以下步骤：
(1)计算
Figure RE-GDA0002893486150000049
Figure RE-GDA00028934861500000410
(2)计算跨中曲率k(ti)和曲率变化率
Figure RE-GDA00028934861500000411
若y(ti)<yc，结构处于弹性运动状态；
Figure RE-GDA00028934861500000412
其中A为弹性阶段的线性比例常数：
Figure RE-GDA0002893486150000051
否则，结构处于塑性运动状态；
Figure RE-GDA0002893486150000052
其中lp为塑性铰长度；
(3)根据截面纤维分层方法，通过跨中曲率k(ti)和曲率率
Figure RE-GDA0002893486150000053
计算截面弯矩M(ti)，
(4)计算抗力R(ti)；
Figure RE-GDA0002893486150000054
第二步：计算载荷转换系数KL和质量转换系数KM：
若结构处于弹性段，即y(ti)＜yc，
Figure RE-GDA0002893486150000055
若结构处于塑性段，即y(ti)≥yc，
Figure RE-GDA0002893486150000056
第三步：计算
Figure RE-GDA0002893486150000057
第四步：计算
Figure RE-GDA0002893486150000058
本发明的有益效果主要体现在以下几个方面。
1、本发明专利建立了近场非均布爆炸冲击波载荷作用下的简支钢筋混凝土梁弯曲动抗力计算方法，并且在抗力计算过程中考虑了材料应变率对抗力的增强效果；
2、本发明专利建立了一种简支钢筋混凝土梁在近场爆炸作用下的等效单自由度动力学模型，模型中考虑了材料应变率效应以及近场非均布爆炸载荷特征；
3、本发明专利能用于快速预测简支钢筋混凝土梁在近场爆炸下的结构弯曲响应。
附图说明
图1是本发明中炸药和结构的位置关系；
图2是本发明的计算流程；
图3是本发明实施例1中计算出的转换系数与载荷均匀性关系；其中图3a为KLE-αp关系图3(b)为KME-αp关系、图3(c)为KLP-αp关系、图3(d)KMP-αp关系。
图4是本发明实施例2的试验工况；
图5是本发明实施例2中计算出弹性、塑性阶段形状函数；
图6是本发明实施例2中计算出的结构静抗力曲线；
图7是文献(US Department of Defense.Structures to resist the effectsof accidental explosions(UFC 3-340-02)[R].Washington DC,2008.中给出的均布荷载及集中荷载下质量和载荷转换系数取值。
以下结合附图及具体实施方式对本发明的具体内容作进一步的详细说明。
具体实施方式
如图1、图2所示，炸药装药位于钢筋混凝土梁跨中正上方，结构长、宽、高分别为L、w、h (m)，炸药TNT当量为W(kg)，炸药几何中心到结构表面的垂直距离爆心距为R1(m)。具体采用如下公式：
r＝R1/W1/3
采用经验公式、数值仿真或实验，获取钢筋混凝土梁迎爆面的跨中方向多个位置处的冲击波压力峰值，并采用数学方法拟合出指数衰减压力峰值分布函数。通常，在近场对称爆炸条件下，结构跨中位置冲击波压力峰值最高，且压力峰值向支座方向快速衰减，比冲量也满足此特征。
因此结构表面的压力峰值形式为：
Figure RE-GDA0002893486150000061
或
Figure RE-GDA0002893486150000062
或其他形式，其中Pm为跨中位置的峰值反射压力，αp为压力分布系数，x为梁表面沿跨度方向的坐标，坐标原点位于跨中，L为梁的长度；
结构表面比冲量形式为：
Figure RE-GDA0002893486150000063
或
Figure RE-GDA0002893486150000064
或其他形式，其中Im为跨中位置的峰值比冲量，αi为比冲量分布系数；
在此处，我们假设
Figure RE-GDA0002893486150000065
来阐述整个计算过程，主要原因为该函数形式在计算过程中方便解析表达。计算过程如下：
根据弹性静力学理论，可计算在压力载荷P(x)的作用下结构的弹性静力学挠曲线方程ye(x)、曲率方程k(x)、弯矩分布M(x)，并根据结构截面尺寸及材料属性确定屈服曲率kc，进而可计算跨中屈服挠度yc：
Figure RE-GDA0002893486150000071
对于简支梁，在塑性阶段，通常在跨中区域形成塑性铰，塑性铰左右两侧结构做近似刚性转动。因此，可分别计算处弹性变形函数
Figure RE-GDA0002893486150000072
和塑性变形函数
Figure RE-GDA0002893486150000073
Figure RE-GDA0002893486150000074
Figure RE-GDA0002893486150000075
进一步的，分别计算弹性阶段载荷转换系数KLE和质量转换系数KME：
Figure RE-GDA0002893486150000076
Figure RE-GDA0002893486150000077
并分别计算塑性阶段载荷转换系数KLP和质量转换系数KMP：
Figure RE-GDA0002893486150000078
Figure RE-GDA0002893486150000079
根据压力和比冲量算结构表面的总冲量Itotal，峰值载荷Fmax，载荷等效作用时间td，并确定等效载荷时间历程F(t)：
Figure RE-GDA00028934861500000710
Figure RE-GDA00028934861500000711
Figure RE-GDA00028934861500000712
Figure RE-GDA0002893486150000081
其中W为为炸药爆炸当量kg，w为结构宽度m。
然后可以建立处结构的等效单自由度动力学模型：
Figure RE-GDA0002893486150000082
其中，KM为质量转换系数，KL为载荷转换系数，M为结构总质量，y、
Figure RE-GDA0002893486150000083
和
Figure RE-GDA0002893486150000084
分别为结构跨中挠度、速度和加速度，R为弯曲抗力；
为了求解该动力学方程，确定确定足够小的求解时间步长Δt，并初始化时间序列ti＝iΔt， i＝0,1,2,3,…n；进一步的，通过中心差分，可得到迭代计算公式；
Figure RE-GDA0002893486150000085
为了进行求解，需要构造
Figure RE-GDA0002893486150000086
y(t0)＝0；然后依次迭代计算y(ti+1)， i＝0,1,2,3,…n，直到|y(ti+1)|＜|y(ti)|；在每一步迭代计算中，包括以下四步：
第一步：根据y(ti)和y(ti-1)计算弯曲抗力R，包括以下步骤：
(1)计算
Figure RE-GDA0002893486150000087
Figure RE-GDA0002893486150000088
(2)计算跨中曲率k(ti)和曲率变化率
Figure RE-GDA0002893486150000089
若y(ti)<yc，结构处于弹性运动状态；
Figure RE-GDA00028934861500000810
其中A为弹性阶段的线性比例常数：
Figure RE-GDA00028934861500000811
否则，结构处于塑性运动状态：
Figure RE-GDA00028934861500000812
其中lp为塑性铰长度；
(3)根据截面纤维分层方法，通过跨中曲率k(ti)和曲率率
Figure RE-GDA0002893486150000091
计算截面弯矩M(ti)，
(4)计算抗力R(ti)；
Figure RE-GDA0002893486150000092
第二步：计算载荷转换系数KL和质量转换系数KM：
若结构处于弹性段，即y(ti)＜yc，
Figure RE-GDA0002893486150000093
若结构处于塑性段，即y(ti)≥yc，
Figure RE-GDA0002893486150000094
第三步：计算
Figure RE-GDA0002893486150000095
第四步：计算
Figure RE-GDA0002893486150000096
实施例1：
针对任一钢筋混凝土简支梁结构，长为L，当炸药放置于该结构跨中正上方时，结构表面载荷形式为
Figure RE-GDA0002893486150000097
其中Pm为跨中位置的峰值反射压力，αp为压力分布系数，x为梁表面沿跨度方向的坐标，坐标原点位于跨中；
根据上文实施方式，可得分别计算弹性阶段载荷转换系数KLE和质量转换系数KME：
Figure RE-GDA0002893486150000098
Figure RE-GDA0002893486150000099
并分别计算塑性阶段载荷转换系数KLP和质量转换系数KMP：
Figure RE-GDA00028934861500000910
Figure RE-GDA00028934861500000911
表1中列举了文献(US Department of Defense.Structures to resist theeffects of accidental explosions(UFC 3-340-02)[R].Washington DC,2008.)中给了的简支梁分别在集中和均布荷载下载荷转换系数和质量转换系数的值。本专利计算了弹性阶段载荷转换系数KLE、质量转换系数KME以及塑性阶段载荷转换系数KLP和质量转换系数KMP，见图3(a)～(d)。当压力分布系数αp接近0.01时，近似为均布荷载；当压力分布系数αp接近100时，近似集中荷载。从图中可见，随着压力分布系数的变化，以上四个参数的取值范围在表1给出的范围之间。且本专利的计算结果与表1中的数据吻合较好。
实施例2
如图4所示，在某次试验中，将3kg柱形TNT悬挂于简支钢筋混凝土梁跨中正上方0.72m 处。钢筋混凝土梁的有效尺寸为，长(L)1.4m，宽(w)0.13m，高(h)0.13m，包覆层厚度(b)为20mm，配筋情况图5所示。
紧邻钢筋混凝土梁安装一钢梁用于安装压力传感器，钢梁上表面与钢筋混凝土梁上表面平齐。压力传感器安装于钢梁中轴线上，共5个压力测点(P1，P2，P3，P4，P5)，测点间距为100mm，最左侧测点(P1)位于钢梁跨中位置。试验中，钢筋混凝土梁宽130mm，钢梁宽100mm，因此，炸药中心在钢筋混凝土梁上表面的投影与测点P1、P2、P3、P4、P5的距离约为115mm，152mm，231mm，321mm和416mm。
试验中压力传感器选用PCB公司的113B22和113B24传感器，量程分别为34MPa(可超量程至68MPa)和6.8MPa(可超量程至13.6MPa)。测试时，以炸药起爆作为数据采集仪器触发信号，数据采样率为1M点/s。
根据试验结果，结构表面载荷函数符合下面公式中的形式：
Figure RE-GDA0002893486150000101
Figure RE-GDA0002893486150000102
其中Pm、Im分别为跨中位置的峰值反射压力及比冲量，αp、αi分别为压力分布系数和冲量分布系数，x为梁表面沿跨度方向的坐标，坐标原点位于跨中，L为梁的长度；进而，可采用数学方法拟合出指数衰减压力峰值分布函数为：
Figure RE-GDA0002893486150000103
Figure RE-GDA0002893486150000104
根据弹性静力学理论，可计算在压力载荷P(x)的作用下结构的弹性静力学挠曲线方程 ye(x)、曲率方程k(x)、弯矩分布M(x)，并根据结构截面尺寸及材料属性确定屈服曲率kc为 0.028m-1，进而可计算跨中屈服挠度yc：
Figure RE-GDA0002893486150000111
对于简支梁，在塑性阶段，通常在跨中区域形成塑性铰，塑性铰左右两侧结构做近似刚性转动。因此，可分别计算处弹性变形函数
Figure RE-GDA0002893486150000112
和塑性变形函数
Figure RE-GDA0002893486150000113
如图5。
Figure RE-GDA0002893486150000114
Figure RE-GDA0002893486150000115
进一步的，分别计算弹性阶段载荷转换系数KLE和质量转换系数KME：
Figure RE-GDA0002893486150000116
Figure RE-GDA0002893486150000117
并分别计算塑性阶段载荷转换系数KLP和质量转换系数KMP：
Figure RE-GDA0002893486150000118
Figure RE-GDA0002893486150000119
根据压力和比冲量算结构表面的总冲量Itotal，峰值载荷Fmax，载荷等效作用时间td，并确定等效载荷时间历程F(t)：
Figure RE-GDA00028934861500001110
Figure RE-GDA00028934861500001111
Figure RE-GDA00028934861500001112
Figure RE-GDA00028934861500001113
然后可以建立处结构的等效单自由度动力学模型：
Figure RE-GDA0002893486150000121
其中，KM为质量转换系数，KL为载荷转换系数，M为结构总质量，y、
Figure RE-GDA0002893486150000122
和
Figure RE-GDA0002893486150000123
分别为结构跨中挠度、速度和加速度，R为弯曲抗力；
为了求解该动力学方程，确定确定足够小的求解时间步长Δt，并初始化时间序列ti＝iΔt， i＝0,1,2,3,…n；进一步的，通过中心差分，可得到迭代计算公式；
Figure RE-GDA0002893486150000124
为了进行求解，需要构造
Figure RE-GDA0002893486150000125
y(t0)＝0；然后依次迭代计算 y(ti+1)，i＝0,1,2,3,…n，直到|y(ti+1)|＜|y(ti)|；在每一步迭代计算中，包括以下四步：
第一步：根据y(ti)和y(ti-1)计算弯曲抗力R(ti)，包括如下步骤：
(1)计算跨中速度，
Figure RE-GDA0002893486150000126
(2)计算跨中曲率k(ti)和曲率变化率
Figure RE-GDA0002893486150000127
若y(ti)<yc，结构处于弹性运动状态；
Figure RE-GDA0002893486150000128
其中A为弹性阶段的线性比例常数：
Figure RE-GDA0002893486150000129
否则，结构处于塑性运动状态；
Figure RE-GDA00028934861500001210
其中lp为塑性铰长度；
(3)采用截面分层分析方法计算跨中弯矩M(ti)，在计算过程中考虑应变率对材料强度的提高以及箍筋包覆区域混凝土强度的提高；
(4)计算抗力：
Figure RE-GDA0002893486150000131
不考虑速度时，采用上述方法可计算出结构抗力位移曲线如图6所示。
第二步：计算载荷转换系数KL和质量转换系数KM：
若结构处于弹性段，即y(ti)＜yc，
Figure RE-GDA0002893486150000132
若结构处于塑性段，即y(ti)≥yc，
Figure RE-GDA0002893486150000133
第三步：计算
Figure RE-GDA0002893486150000134
第四步：计算
Figure RE-GDA0002893486150000135
通过计算，可得跨中峰值挠度为22.1mm，试验测量结果为25.06mm，误差约11.8％。以上计算表明，该方法能用于预测近场爆炸下简支钢筋混凝土梁的弯曲响应。

Claims (10)

1.一种近场爆炸下简支钢筋混凝土梁弯曲响应预测方法，其特征在于：该方法包括如下步骤：
步骤1：计算炸药到结构表面的比例距离；
步骤2：拟合爆炸冲击波压力、比冲量分布函数；
步骤3：计算弹性、塑性阶段的结构变形函数；
步骤4：计算弹性、塑性阶段的载荷转换系数和质量转换系数；
步骤5：计算等效载荷时间历程；
步骤6：建立简支钢筋混凝土梁单自由度动力学方程；
步骤7：预测结构弯曲响应。
2.如权利要求1所述的近场爆炸下简支钢筋混凝土梁弯曲响应预测方法，其特征在于：所述的步骤1：通过炸药爆炸当量及炸药爆心距确定炸药中心到钢筋混凝土梁结构表面的比例距离r(m/kg1/3)，具体采用如下公式：
r＝R1/W1/3
其中，W为炸药爆炸当量kg，R1为炸药爆心距m。
3.如权利要求1所述的近场爆炸下简支钢筋混凝土梁弯曲响应预测方法，其特征在于：所述的步骤2采用经验公式或数值仿真，计算钢筋混凝土梁迎爆面的冲击波压力峰值分布P(x)以及比冲量分布I(x)，其中
Figure FDA0002790499460000011
或
Figure FDA0002790499460000012
Figure FDA0002790499460000013
或其他形式，Pm为跨中位置的峰值反射压力，αp为压力分布系数，x为梁表面沿跨度方向的坐标，坐标原点位于跨中，L为梁的长度；
Figure FDA0002790499460000014
或
Figure FDA0002790499460000015
或其他形式，其中Im为跨中位置的峰值比冲量，αi为比冲量分布系数。
4.如权利要求3所述的近场爆炸下简支钢筋混凝土梁弯曲响应预测方法，其特征在于：所述的步骤1，进一步计算压力载荷P(x)作用下结构的弹性挠曲线方程ye(x)、曲率方程k(x)、弯矩方程M(x)，并根据结构截面尺寸及材料属性确定屈服曲率kc，并计算跨中屈服挠度yc：
Figure FDA0002790499460000016
5.如权利要求1所述的近场爆炸下简支钢筋混凝土梁弯曲响应预测方法，其特征在于：所述的步骤3:计算弹性变形函数
Figure FDA0002790499460000021
和塑性变形函数
Figure FDA0002790499460000022
Figure FDA0002790499460000023
Figure FDA0002790499460000024
x为梁表面沿跨度方向的坐标，坐标原点位于跨中，梁的长度为L。
6.如权利要求1所述近场爆炸下简支钢筋混凝土梁弯曲响应预测方法，其特征在于：所述的步骤4:计算弹性、塑性阶段的载荷转换系数和质量转换系数，其中弹性阶段载荷转换系数KLE和质量转换系数KME：
Figure FDA0002790499460000025
Figure FDA0002790499460000026
塑性阶段载荷转换系数KLP和质量转换系数KMP为：
Figure FDA0002790499460000027
Figure FDA0002790499460000028
7.如权利要求1所述近场爆炸下简支钢筋混凝土梁弯曲响应预测方法，其特征在于：所述的步骤5计算等效载荷时间历程，根据压力和比冲量算结构表面的总冲量Itotal，峰值载荷Fmax，载荷等效作用时间td，并确定等效载荷时间历程F(t)：
Figure FDA0002790499460000029
Figure FDA00027904994600000210
Figure FDA00027904994600000211
Figure FDA0002790499460000031
其中W为为炸药爆炸当量kg，w为结构宽度m。
8.如权利要求1所述近场爆炸下简支钢筋混凝土梁弯曲响应预测方法，其特征在于：
所述的步骤6建立简支钢筋混凝土梁单自由度动力学方程如下：
Figure FDA0002790499460000032
其中，KM为质量转换系数，KL为载荷转换系数，M为结构总质量，y、
Figure FDA0002790499460000033
和
Figure FDA0002790499460000034
分别为结构跨中挠度、速度和加速度，R为弯曲抗力。
9.如权利要求1所述近场爆炸下简支钢筋混凝土梁弯曲响应预测方法，其特征在于：所述的步骤7预测结构弯曲响应，即确定足够小的求解时间步长Δt，并初始化时间序列ti＝iΔt，i＝0,1,2,3,…n；对加速度项进行中心差分得到迭代公式；
Figure FDA0002790499460000035
构造
Figure FDA0002790499460000036
y(t0)＝0；然后依次迭代计算y(ti+1)，i＝0,1,2,3,…n，直到|y(ti+1)|＜|y(ti)|。
10.如权利要求9所述近场爆炸下简支钢筋混凝土梁弯曲响应预测方法，其特征在于：所述的步骤7迭代计算包括如下步骤：
第一步：根据y(ti)和y(ti-1)计算弯曲抗力R(ti)，包括以下步骤：
(1)计算
Figure FDA0002790499460000037
Figure FDA0002790499460000038
(2)计算跨中曲率k(ti)和曲率变化率
Figure FDA0002790499460000039
若y(ti)<yc，结构处于弹性运动状态；
Figure FDA00027904994600000310
其中A为弹性阶段的线性比例常数：
Figure FDA0002790499460000041
否则，结构处于塑性运动状态；
Figure FDA0002790499460000042
其中lp为塑性铰长度；
(3)根据截面纤维分层方法，通过跨中曲率k(ti)和曲率率
Figure FDA0002790499460000043
计算截面弯矩M(ti)，
(4)计算抗力R(ti)；
Figure FDA0002790499460000044
第二步：计算载荷转换系数KL和质量转换系数KM：
若结构处于弹性段，即y(ti)＜yc，
Figure FDA0002790499460000045
若结构处于塑性段，即y(ti)≥yc，
Figure FDA0002790499460000046
第三步：计算
Figure FDA0002790499460000047
第四步：计算
Figure FDA0002790499460000048
CN202011313265.9A 2020-11-20 2020-11-20 一种近场爆炸下简支钢筋混凝土梁弯曲响应预测方法 Active CN112464338B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011313265.9A CN112464338B (zh) 2020-11-20 2020-11-20 一种近场爆炸下简支钢筋混凝土梁弯曲响应预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011313265.9A CN112464338B (zh) 2020-11-20 2020-11-20 一种近场爆炸下简支钢筋混凝土梁弯曲响应预测方法

Publications (2)

Publication Number Publication Date
CN112464338A true CN112464338A (zh) 2021-03-09
CN112464338B CN112464338B (zh) 2023-09-12

Family

ID=74799061

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011313265.9A Active CN112464338B (zh) 2020-11-20 2020-11-20 一种近场爆炸下简支钢筋混凝土梁弯曲响应预测方法

Country Status (1)

Country Link
CN (1) CN112464338B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100256957A1 (en) * 2009-04-01 2010-10-07 Livermore Software Technology Corporation Explosion Simulation in Finite Element Analysis
JP2011257237A (ja) * 2010-06-08 2011-12-22 Ohbayashi Corp 地震被害予測方法
CN103512815A (zh) * 2013-09-25 2014-01-15 天津大学 爆炸荷载下钢筋混凝土柱动态响应判断方法
JP2014149222A (ja) * 2013-02-01 2014-08-21 Taisei Corp コンクリートの爆裂推測方法
CN109580193A (zh) * 2018-11-14 2019-04-05 中北大学 一种在爆炸荷载作用下等截面简支梁荷载系数的计算方法
CN109740300A (zh) * 2019-03-13 2019-05-10 中国人民解放军军事科学院国防工程研究院 基于等效单自由度方法的钢筋混凝土结构构件评估方法
CN109765025A (zh) * 2018-12-25 2019-05-17 哈尔滨理工大学 爆炸荷载作用下基于p-i曲线的rpc板损伤评估方法
CN109933903A (zh) * 2019-03-13 2019-06-25 中国人民解放军军事科学院国防工程研究院 一种构件表面爆炸载荷分布规律的分析方法
CN109933904A (zh) * 2019-03-13 2019-06-25 中国人民解放军军事科学院国防工程研究院 基于等效单自由度的钢筋混凝土结构构件p-i曲线评估方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20100256957A1 (en) * 2009-04-01 2010-10-07 Livermore Software Technology Corporation Explosion Simulation in Finite Element Analysis
JP2011257237A (ja) * 2010-06-08 2011-12-22 Ohbayashi Corp 地震被害予測方法
JP2014149222A (ja) * 2013-02-01 2014-08-21 Taisei Corp コンクリートの爆裂推測方法
CN103512815A (zh) * 2013-09-25 2014-01-15 天津大学 爆炸荷载下钢筋混凝土柱动态响应判断方法
CN109580193A (zh) * 2018-11-14 2019-04-05 中北大学 一种在爆炸荷载作用下等截面简支梁荷载系数的计算方法
CN109765025A (zh) * 2018-12-25 2019-05-17 哈尔滨理工大学 爆炸荷载作用下基于p-i曲线的rpc板损伤评估方法
CN109740300A (zh) * 2019-03-13 2019-05-10 中国人民解放军军事科学院国防工程研究院 基于等效单自由度方法的钢筋混凝土结构构件评估方法
CN109933903A (zh) * 2019-03-13 2019-06-25 中国人民解放军军事科学院国防工程研究院 一种构件表面爆炸载荷分布规律的分析方法
CN109933904A (zh) * 2019-03-13 2019-06-25 中国人民解放军军事科学院国防工程研究院 基于等效单自由度的钢筋混凝土结构构件p-i曲线评估方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
张秀华;段忠东;张春巍;: "爆炸荷载作用下钢筋混凝土梁的动力响应和破坏过程分析", 东北林业大学学报, no. 04, pages 50 - 53 *
方秦;陈力;张亚栋;柳锦春;: "爆炸荷载作用下钢筋混凝土结构的动态响应与破坏模式的数值分析", 工程力学, no. 2, pages 135 - 144 *
汪维;刘光昆;赵强;刘瑞朝;吴飚;周松柏;: "近爆作用下方形板表面爆炸载荷分布函数研究", 中国科学:物理学 力学 天文学, no. 02, pages 1 - 9 *
赵向军: "爆炸冲击载荷作用下钢筋混凝土介质裂纹扩展速度的实验研究", 火炸药学报, vol. 36, no. 1, pages 55 - 58 *

Also Published As

Publication number Publication date
CN112464338B (zh) 2023-09-12

Similar Documents

Publication Publication Date Title
Chaise et al. Modelling of multiple impacts for the prediction of distortions and residual stresses induced by ultrasonic shot peening (USP)
Gariépy et al. On the effect of the peening trajectory in shot peen forming
CN112464337B (zh) 一种近场爆炸下简支钢筋混凝土梁弯曲动抗力计算方法
Liao et al. Study on explosion resistance performance experiment and damage assessment model of high-strength reinforcement concrete beams
Mohammadzadeh et al. An analytical and numerical investigation on the dynamic responses of steel plates considering the blast loads
Yan et al. Compression behaviour of stainless-steel stub square tubular columns at cold-region low temperatures
Cybulski et al. Local buckling of cold-formed elements used in arched building with geometrical imperfections
Mirzaei et al. Analytical and numerical modeling of the transient elasto-dynamic response of a cylindrical tube to internal gaseous detonation
Machado Non-linear buckling and postbuckling behavior of thin-walled beams considering shear deformation
Quayyum et al. Initial residual stresses in hot-rolled wide-flange shapes: A computational technique and influence on structural performances
Wei et al. Modification of SDOF model for reinforced concrete beams under close-in explosion
Prochazka et al. Assessment of laminated cylindrical arch loaded by a shock wave
CN112464338A (zh) 一种近场爆炸下简支钢筋混凝土梁弯曲响应预测方法
Walentyński et al. Stability analyses and experimental investigation of doubly corrugated steel arch panels
Forcellese et al. Computer aided engineering of the sheet bending process
Kong et al. Research on the dynamic buckling of a typical deck grillage structure subjected to in-plane impact Load
Walentyński et al. Numerical stability analyses and preliminary experimental investigation of doubly corrugated steel arch panels
Lee Recent finite element applications of the incremental endochronic plasticity
Huang et al. Finite deformation analysis of the elastic circular plates under pressure loading
Manjuprasad et al. Non-linear dynamic response of a reinforced concrete secondary containment shell subjected to seismic load
Zhang et al. Eccentric compressive distortional buckling and design of non-symmetric cold-formed angular column with complex edges
Ikushima et al. Numerical analysis of residual stress distribution on peening process
Cho et al. Dynamic response and damage of composite shell under impact
Walentyński et al. Analysis of local stability of doubly corrugated cold-framed profiles
Luo et al. The mechanical behavior analysis of CFCC with overall anisotropic damage by the micro-macro scale method

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