CN112270032B - Method for designing dynamic coefficient of anti-exponential explosion load of critical structure with damping flexible rigidity - Google Patents

Method for designing dynamic coefficient of anti-exponential explosion load of critical structure with damping flexible rigidity Download PDF

Info

Publication number
CN112270032B
CN112270032B CN202011304175.3A CN202011304175A CN112270032B CN 112270032 B CN112270032 B CN 112270032B CN 202011304175 A CN202011304175 A CN 202011304175A CN 112270032 B CN112270032 B CN 112270032B
Authority
CN
China
Prior art keywords
load
stage
vibration
explosion
structural
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.)
Expired - Fee Related
Application number
CN202011304175.3A
Other languages
Chinese (zh)
Other versions
CN112270032A (en
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.)
North University of China
Original Assignee
North University of China
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 North University of China filed Critical North University of China
Priority to CN202011304175.3A priority Critical patent/CN112270032B/en
Publication of CN112270032A publication Critical patent/CN112270032A/en
Application granted granted Critical
Publication of CN112270032B publication Critical patent/CN112270032B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13Differential equations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Geometry (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Data Mining & Analysis (AREA)
  • Civil Engineering (AREA)
  • Architecture (AREA)
  • Structural Engineering (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Buildings Adapted To Withstand Abnormal External Influences (AREA)

Abstract

The invention relates to a dynamic coefficient method for designing an anti-exponential explosion load of a damping-containing flexible rigid critical structure, which belongs to the technical field of antiknock design, and specifically relates to a dynamic coefficient method for designing an anti-exponential explosion load of a damping-containing flexible rigid critical structure, which comprises the following steps: exponential blast load action real time t i The value is just that the structural vibration reaches the maximum value y of the elastic vibration T At a certain moment, after the explosive load is unloaded, the structure starts to perform plastic vibration by means of inertia force, and at a certain moment t m The maximum value y of structural elastoplastic displacement is reached m The method comprises the steps of carrying out a first treatment on the surface of the According to the action process of explosion on the building structure, dividing the process into two stages of forced vibration in the explosion load action elastic stage and free vibration in the explosion load unloading plastic stage; through the calculation of the different stages, the influence of the structural type, the accurate description of the explosion load and the structural damping on the explosion load dynamic coefficient is fully considered, so that the designed building structure is more practical as much as possible, the cost is reduced, and the antiknock requirement is met.

Description

含阻尼柔性刚性临界结构抗指数型爆炸荷载设计动力系数 方法Design dynamic coefficients of damped flexible rigid critical structures against exponential blast loads method

技术领域technical field

本发明涉及一种含阻尼柔性刚性临界结构抗指数型爆炸荷载设计动力系数方法,属于抗爆设计技术领域。The invention relates to a method for designing dynamic coefficients of a flexible rigid critical structure with damping against exponential explosion loads, and belongs to the technical field of explosion-resistant design.

背景技术Background technique

目前,现有的建筑物在进行抗爆设计时,均需要考虑抗爆结构构件的抗爆能力。常规爆炸施加在结构上的爆炸荷载作用时长ti很短,我国及国外人防结构进行抗爆设计时,按等冲量线性荷载处理,进一步将该爆炸动荷载超压峰值Δpm与结构弹塑性阶段抗力动力系数kh相乘后,作为静载进行结构抗爆设计数值。其中规范给出的动力系数kh公式,未考虑结构阻尼和抗爆设计结构类型,对实际结构精准化设计缺乏有力支撑;规范简化的等冲量线性衰减荷载等效作用时长也小于真实爆炸作用时长。如:耿少波,李洪,葛培杰.考虑跃迁的指数型炸药空爆荷载等效静载动力系数[J].爆炸与冲击和2019,39(03):33-41;BAKER WE.Explosion hazards and evaluation[M].Amsterdam:Elsevier ScientificPub.Co.1983.两部文献中,分别介绍了真实的空气爆炸冲击波衰减模式采用指数型函数描述将会使得计算结果更加准确。因此现有规范在进行结构抗爆设计时,设计偏保守,在部分应用环境中,偏保守的设计,将造成建设成本的增加。At present, when the anti-blast design of existing buildings is carried out, the anti-blast capability of the anti-blast structural components needs to be considered. The duration t i of the explosive load imposed on the structure by conventional explosion is very short. When the anti-explosion design of civil air defense structures in China and abroad is carried out, it is treated as an equal-impulse linear load, and the peak value of the overpressure of the explosive dynamic load After multiplied by the resistance dynamic coefficient k h , it is used as the static load to carry out the anti-blast design value of the structure. Among them, the dynamic coefficient k h formula given in the code does not consider the structural damping and anti-blast design structure type, which lacks strong support for the precise design of the actual structure; the equivalent duration of the equal-impulse linear attenuation load simplified by the code is also shorter than the real explosion duration . Such as: Geng Shaobo, Li Hong, Ge Peijie. Equivalent static load dynamic coefficient of exponential explosive air-burst load considering transition[J]. Explosion and Shock and 2019,39(03):33-41; BAKER WE.Explosion hazards and evaluation[ M].Amsterdam:Elsevier ScientificPub.Co.1983. In the two documents, it is introduced that the real air explosion shock wave attenuation mode is described by an exponential function, which will make the calculation result more accurate. Therefore, when the existing codes carry out the anti-blast design of the structure, the design is conservative. In some application environments, the conservative design will increase the construction cost.

发明内容Contents of the invention

为解决现有技术存在的技术问题,本发明提供了一种含阻尼柔性刚性临界结构抗指数型爆炸荷载设计动力系数方法。In order to solve the technical problems existing in the prior art, the present invention provides a method for designing dynamic coefficients of a flexible rigid critical structure with damping against exponential explosion loads.

为实现上述目的,本发明所采用的技术方案为含阻尼柔性刚性临界结构抗指数型爆炸荷载设计动力系数方法,所述的柔性刚性临界结构指的是:爆炸荷载作用真实时长ti数值恰为结构振动达到弹性振动最大值yT时刻,爆炸荷载卸载后,结构依靠惯性力开始进行塑性振动,在某一时刻tm,达到了结构弹塑性位移最大值ymIn order to achieve the above-mentioned purpose, the technical solution adopted in the present invention is a method for designing dynamic coefficients of a flexible rigid critical structure with resistance to exponential explosion loads, and the flexible rigid critical structure refers to: the actual duration t of the explosive load is just When the structural vibration reaches the maximum value of elastic vibration y T , after the explosion load is unloaded, the structure starts to undergo plastic vibration relying on the inertial force, and at a certain time t m , the maximum elastic-plastic displacement of the structure y m is reached;

根据爆炸对建筑结构的作用过程,将该过程分为爆炸荷载作用弹性阶段强迫振动和爆炸荷载卸载塑性阶段自由振动两个阶段;According to the action process of the explosion on the building structure, the process is divided into two stages: the forced vibration in the elastic stage of the explosion load and the free vibration in the plastic stage of the explosion load unloading;

a、爆炸荷载作用弹性阶段强迫振动a. Forced vibration in the elastic stage of explosion load

在弹性阶段且在荷载作用时长范围0<t≤ti内,结构等效体系的动力方程为In the elastic stage and within the range of load duration 0<t≤t i , the dynamic equation of the structural equivalent system is

Figure BDA0002787793000000021
Figure BDA0002787793000000021

其中,t为结构振动时间参数,ti为爆炸荷载持续时长,Me为弹性阶段等效结构质量,Ce为弹性阶段等效结构阻尼,Ke为弹性阶段等效结构刚度,

Figure BDA0002787793000000022
为结构等效体系振动加速度,/>
Figure BDA0002787793000000023
为结构等效体系振动速度,y为结构等效体系振动位移,ΔPe(t)为结构承受的随时间t变化的爆炸动荷载,等效结构系数计算公式分别为:Among them, t is the structural vibration time parameter, t i is the duration of the explosion load, M e is the equivalent structural mass in the elastic stage, C e is the equivalent structural damping in the elastic stage, K e is the equivalent structural stiffness in the elastic stage,
Figure BDA0002787793000000022
is the vibration acceleration of the structural equivalent system, />
Figure BDA0002787793000000023
is the vibration velocity of the structural equivalent system, y is the vibration displacement of the structural equivalent system, ΔP e (t) is the explosive dynamic load on the structure that changes with time t, and the calculation formulas of the equivalent structural coefficient are:

Figure BDA0002787793000000024
Figure BDA0002787793000000024

其中,m为真实结构每延米质量,l为真实结构跨长,ξ为真实结构阻尼比,K为真实结构刚度,kM为弹性阶段质量变换系数,kL为弹性阶段荷载变换系数。指数型爆炸动荷载为:Among them, m is the mass per linear meter of the real structure, l is the span length of the real structure, ξ is the damping ratio of the real structure, K is the stiffness of the real structure, k M is the mass conversion coefficient in the elastic stage, and k L is the load conversion coefficient in the elastic stage. The exponential explosion dynamic load is:

Figure BDA0002787793000000025
Figure BDA0002787793000000025

其中,ti为爆炸荷载作用时长,Δpm为爆炸荷载超压峰值,a为指数型荷载衰减系数,由无阻尼自振频率ω、含阻尼自振频率ωd、阻尼调整系数γ、爆炸荷载峰值作为静载时对应的静位移yst、弹性阶段质量变换系数与荷载变换系数之比值kM-L,各参数计算如下:Among them, t i is the duration of the blast load, Δp m is the overpressure peak value of the blast load, a is the exponential load attenuation coefficient, which is composed of the undamped natural frequency ω, the damped natural frequency ω d , the damping adjustment coefficient γ, and the explosive load The peak value is used as the static displacement y st corresponding to the static load, and the ratio k ML of the mass conversion coefficient in the elastic stage to the load conversion coefficient k ML . The parameters are calculated as follows:

Figure BDA0002787793000000031
Figure BDA0002787793000000031

可知,结构承受爆炸荷载之前的初始位移及速度均为0,求解该微分方程后,可确定此阶段位移和速度表达式为:It can be seen that the initial displacement and velocity of the structure before bearing the explosion load are both 0. After solving the differential equation, the expression of displacement and velocity at this stage can be determined as:

Figure BDA0002787793000000032
Figure BDA0002787793000000032

Figure BDA0002787793000000041
Figure BDA0002787793000000041

在爆炸荷载作用结束卸载后,完成弹性阶段的振动时,对应ti时刻位移yi和速度vi为:After the explosion load is unloaded and the vibration in the elastic stage is completed, the displacement y i and velocity v i corresponding to time t i are:

Figure BDA0002787793000000042
Figure BDA0002787793000000042

Figure BDA0002787793000000051
Figure BDA0002787793000000051

b、爆炸荷载卸载塑性阶段自由振动b. Free vibration in the plastic stage of explosion load unloading

当结构振动时刻大于ti时刻,爆炸荷载消失,且结构刚进入塑性振动状态,此时结构为以yi及vi为初始条件的含阻尼塑性阶段的自由振动,即当ti<t<tm时,结构等效体系动力方程为When the structural vibration moment is greater than t i , the explosion load disappears, and the structure just enters the plastic vibration state. At this time, the structure is free vibration in the damped plastic stage with y i and v i as the initial conditions, that is, when t i <t< At t m , the dynamic equation of the structural equivalent system is

Figure BDA0002787793000000052
Figure BDA0002787793000000052

其中,me为塑性阶段等效结构质量,ce为塑性阶段等效结构阻尼,qm为结构塑性阶段的最大抗力,其计算公式为:Among them, m e is the equivalent structural mass in the plastic stage, c e is the equivalent structural damping in the plastic stage, q m is the maximum resistance force in the plastic stage of the structure, and its calculation formula is:

Figure BDA0002787793000000053
Figure BDA0002787793000000053

其中,km为塑性阶段质量变换系数,kl为塑性阶段荷载变换系数,解此运动微分方程,求出此阶段位移和速度解为:Among them, k m is the mass transformation coefficient in the plastic stage, and k l is the load transformation coefficient in the plastic stage. Solve this differential equation of motion to obtain the solution of displacement and velocity in this stage:

Figure BDA0002787793000000061
Figure BDA0002787793000000061

Figure BDA0002787793000000062
Figure BDA0002787793000000062

c、弹塑性阶段基于动力系数的延性比c. Ductility ratio based on dynamic coefficient in elastic-plastic stage

当结构振动至最大位移ym时,对应的时刻为tm,此时速度vm=0,代入(12)式,则:When the structure vibrates to the maximum displacement y m , the corresponding moment is t m , and the velocity v m =0 at this time, substituting into formula (12), then:

Figure BDA0002787793000000063
Figure BDA0002787793000000063

将tm带入到(11)中得出结构弹塑性振动最大位移为:Substituting t m into (11), the maximum displacement of the structural elastic-plastic vibration is:

Figure BDA0002787793000000064
Figure BDA0002787793000000064

令km-l为塑性阶段质量变换系数与荷载变换系数之比值,即km-l=km/kl,则Let k ml be the ratio of the mass conversion coefficient to the load conversion coefficient in the plastic stage, that is, k ml =k m /k l , then

Figure BDA0002787793000000065
Figure BDA0002787793000000065

由弹塑性阶段的抗爆设计,可知对于柔性结构、刚性结构临界区分结构,抗力动力系数kh和延性β分别为:From the anti-knock design in the elastoplastic stage, it can be known that for the flexible structure and the rigid structure critically divided structure, the resistance dynamic coefficient k h and the ductility β are respectively:

Figure BDA0002787793000000066
Figure BDA0002787793000000066

Figure BDA0002787793000000067
Figure BDA0002787793000000067

将(15)带入延性比公式(17)后,则得到After substituting (15) into the ductility ratio formula (17), we get

Figure BDA0002787793000000071
Figure BDA0002787793000000071

其中,yi与vi的表达式为(7)、(8)。Among them, the expressions of y i and vi are (7), (8).

与现有技术相比,本发明具有以下技术效果:本发明根据实际情况,充分考虑结构的类型及结构阻尼对爆炸荷载动力系数的影响,尽量使设计的建筑结构更符合实际,在降低成本的同时,满足抗爆要求。并且通过该方法能够实现实际结构的精准化设计,也为抗爆设计奠定基础。Compared with the prior art, the present invention has the following technical effects: according to the actual situation, the present invention fully considers the type of structure and the influence of structural damping on the explosive load dynamic coefficient, and tries to make the designed building structure more realistic, while reducing costs. At the same time, it meets the anti-knock requirements. And through this method, the precise design of the actual structure can be realized, and it also lays the foundation for the anti-knock design.

具体实施方式Detailed ways

为了使本发明所要解决的技术问题、技术方案及有益效果更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。In order to make the technical problems, technical solutions and beneficial effects to be solved by the present invention clearer, the present invention will be further described in detail below in conjunction with the embodiments. It should be understood that the specific embodiments described here are only used to explain the present invention, not to limit the present invention.

含阻尼柔性刚性临界结构抗指数型爆炸荷载设计动力系数方法,所述的柔性刚性临界结构指的是:爆炸荷载作用真实时长ti数值恰为结构振动达到弹性振动最大值yT时刻,爆炸荷载卸载后,结构依靠惯性力开始进行塑性振动,在某一时刻tm,达到了结构弹塑性位移最大值ymThe dynamic coefficient method for designing a flexible rigid critical structure with damping against exponential explosion loads, the flexible rigid critical structure refers to: the real duration t i value of the explosion load is just the moment when the structural vibration reaches the maximum value of elastic vibration y T , and the explosion load After unloading, the structure begins to vibrate plastically depending on the inertial force, and at a certain time t m , the maximum elastic-plastic displacement of the structure y m is reached;

根据爆炸对建筑结构的作用过程,将该过程分为爆炸荷载作用弹性阶段强迫振动和爆炸荷载卸载塑性阶段自由振动两个阶段;According to the action process of the explosion on the building structure, the process is divided into two stages: the forced vibration in the elastic stage of the explosion load and the free vibration in the plastic stage of the explosion load unloading;

a、爆炸荷载作用弹性阶段强迫振动a. Forced vibration in the elastic stage of explosion load

在弹性阶段且在荷载作用时长范围0<t≤ti内,结构等效体系的动力方程为In the elastic stage and within the range of load duration 0<t≤t i , the dynamic equation of the structural equivalent system is

Figure BDA0002787793000000081
Figure BDA0002787793000000081

其中,t为结构振动时间参数,ti为爆炸荷载持续时长,Me为弹性阶段等效结构质量,Ce为弹性阶段等效结构阻尼,Ke为弹性阶段等效结构刚度,

Figure BDA0002787793000000082
为结构等效体系振动加速度,/>
Figure BDA0002787793000000083
为结构等效体系振动速度,y为结构等效体系振动位移,ΔPe(t)为结构承受的随时间t变化的爆炸动荷载,等效结构系数计算公式分别为:Among them, t is the structural vibration time parameter, t i is the duration of the explosion load, M e is the equivalent structural mass in the elastic stage, C e is the equivalent structural damping in the elastic stage, K e is the equivalent structural stiffness in the elastic stage,
Figure BDA0002787793000000082
is the vibration acceleration of the structural equivalent system, />
Figure BDA0002787793000000083
is the vibration velocity of the structural equivalent system, y is the vibration displacement of the structural equivalent system, ΔP e (t) is the explosive dynamic load on the structure that changes with time t, and the calculation formulas of the equivalent structural coefficient are:

Figure BDA0002787793000000084
Figure BDA0002787793000000084

其中,m为真实结构每延米质量,l为真实结构跨长,ξ为真实结构阻尼比,K为真实结构刚度,kM为弹性阶段质量变换系数,kL为弹性阶段荷载变换系数。指数型爆炸动荷载为:Among them, m is the mass per linear meter of the real structure, l is the span length of the real structure, ξ is the damping ratio of the real structure, K is the stiffness of the real structure, k M is the mass conversion coefficient in the elastic stage, and k L is the load conversion coefficient in the elastic stage. The exponential explosion dynamic load is:

Figure BDA0002787793000000085
Figure BDA0002787793000000085

其中,ti为爆炸荷载作用时长,Δpm为爆炸荷载超压峰值,a为指数型荷载衰减系数,由无阻尼自振频率ω、含阻尼自振频率ωd、阻尼调整系数γ、爆炸荷载峰值作为静载时对应的静位移yst、弹性阶段质量变换系数与荷载变换系数之比值kM-L,各参数计算如下:Among them, t i is the duration of the blast load, Δp m is the overpressure peak value of the blast load, a is the exponential load attenuation coefficient, which consists of the undamped natural frequency ω, the damped natural frequency ω d , the damping adjustment coefficient γ, and the explosive load The peak value is used as the static displacement y st corresponding to the static load, and the ratio k ML of the mass conversion coefficient in the elastic stage to the load conversion coefficient k ML . The parameters are calculated as follows:

Figure BDA0002787793000000086
Figure BDA0002787793000000086

可知,结构承受爆炸荷载之前的初始位移及速度均为0,求解该微分方程后,可确定此阶段位移和速度表达式为:It can be seen that the initial displacement and velocity of the structure before bearing the explosion load are both 0. After solving the differential equation, the expression of displacement and velocity at this stage can be determined as:

Figure BDA0002787793000000091
Figure BDA0002787793000000091

Figure BDA0002787793000000092
Figure BDA0002787793000000092

在爆炸荷载作用结束卸载后,完成弹性阶段的振动时,对应ti时刻位移yi和速度vi为:After the explosion load is unloaded and the vibration in the elastic stage is completed, the displacement y i and velocity v i corresponding to time t i are:

Figure BDA0002787793000000101
Figure BDA0002787793000000101

Figure BDA0002787793000000102
Figure BDA0002787793000000102

b、爆炸荷载卸载塑性阶段自由振动b. Free vibration in the plastic stage of explosion load unloading

当结构振动时刻大于ti时刻,爆炸荷载消失,且结构刚进入塑性振动状态,此时结构为以yi及vi为初始条件的含阻尼塑性阶段的自由振动,即当ti<t<tm时,结构等效体系动力方程为When the structural vibration moment is greater than t i , the explosion load disappears, and the structure just enters the plastic vibration state. At this time, the structure is free vibration in the damped plastic stage with y i and v i as the initial conditions, that is, when t i <t< At t m , the dynamic equation of the structural equivalent system is

Figure BDA0002787793000000111
Figure BDA0002787793000000111

其中,me为塑性阶段等效结构质量,ce为塑性阶段等效结构阻尼,qm为结构塑性阶段的最大抗力,其计算公式为:Among them, m e is the equivalent structural mass in the plastic stage, c e is the equivalent structural damping in the plastic stage, q m is the maximum resistance force in the plastic stage of the structure, and its calculation formula is:

Figure BDA0002787793000000112
Figure BDA0002787793000000112

其中,km为塑性阶段质量变换系数,kl为塑性阶段荷载变换系数,解此运动微分方程,求出此阶段位移和速度解为:Among them, k m is the mass transformation coefficient in the plastic stage, and k l is the load transformation coefficient in the plastic stage. Solve this differential equation of motion to obtain the solution of displacement and velocity in this stage:

Figure BDA0002787793000000113
Figure BDA0002787793000000113

Figure BDA0002787793000000114
Figure BDA0002787793000000114

c、弹塑性阶段基于动力系数的延性比c. Ductility ratio based on dynamic coefficient in elastic-plastic stage

当结构振动至最大位移ym时,对应的时刻为tm,此时速度vm=0,代入(12)式,则:When the structure vibrates to the maximum displacement y m , the corresponding moment is t m , and the velocity v m =0 at this time, substituting into formula (12), then:

Figure BDA0002787793000000115
Figure BDA0002787793000000115

将tm带入到(11)中得出结构弹塑性振动最大位移为:Substituting t m into (11), the maximum displacement of the structural elastic-plastic vibration is:

Figure BDA0002787793000000121
Figure BDA0002787793000000121

令km-l为塑性阶段质量变换系数与荷载变换系数之比值,即km-l=km/kl,则Let k ml be the ratio of the mass conversion coefficient to the load conversion coefficient in the plastic stage, that is, k ml =k m /k l , then

Figure BDA0002787793000000122
Figure BDA0002787793000000122

由弹塑性阶段的抗爆设计,可知对于柔性结构、刚性结构临界区分结构,抗力动力系数kh和延性β分别为:From the anti-knock design in the elastoplastic stage, it can be known that for the flexible structure and the rigid structure critically divided structure, the resistance dynamic coefficient k h and the ductility β are respectively:

Figure BDA0002787793000000123
Figure BDA0002787793000000123

Figure BDA0002787793000000124
Figure BDA0002787793000000124

将(15)带入延性比公式(17)后,则得到After substituting (15) into the ductility ratio formula (17), we get

Figure BDA0002787793000000125
Figure BDA0002787793000000125

其中,yi与vi的表达式为(7)、(8)。Among them, the expressions of y i and vi are (7), (8).

通过上述方法,结构实际抗爆设计进行如下举例。Through the above method, the actual anti-knock design of the structure is given as an example as follows.

1、一种防护结构进行抗爆设计时,要求结构设计为柔性刚性临界过渡性结构,该结构自振圆频率ω与爆炸动荷载作用时长ti乘积ωti为3.2,延性比β为1.61,阻尼比ξ为1%时,指数型荷载衰减系数a为1.27,km-l与km-l取值分别为0.66、0.78(《地下防护结构》方秦,柳锦春编著,ISBN 9787508470009中国水利水电出版社,2010年出版),其它参数由其它公式分别计算得出,代入式(7)、(8)、(18),由迭代方法求出抗力动力系数kh为0.574。1. When a protective structure is designed for explosion resistance, it is required that the structure be designed as a flexible rigid critical transitional structure. The product ωt i of the circular frequency ω of the structure and the duration t i of the explosive dynamic load is 3.2, and the ductility ratio β is 1.61. When the damping ratio ξ is 1%, the exponential load attenuation coefficient a is 1.27, and the values of k ml and k ml are 0.66 and 0.78 respectively ("Underground Protective Structure", edited by Fang Qin and Liu Jinchun, ISBN 9787508470009 China Water Resources and Hydropower Press, 2010 Published in 2010), other parameters were calculated from other formulas, and substituted into formulas (7), (8) and (18), and the resistance dynamic coefficient k h was obtained as 0.574 by the iterative method.

2、一种防护结构进行抗爆设计时,要求结构设计为柔性刚性临界过渡性结构,该结构自振圆频率ω与爆炸动荷载作用时长ti乘积ωti为2.93,延性比β为1.26,阻尼比ξ为5%时,指数型荷载衰减系数a为1.27,km-l与km-l取值分别为0.66、0.78(《地下防护结构》方秦,柳锦春编著,ISBN 9787508470009中国水利水电出版社,2010年出版),其它参数由其它公式分别计算得出,代入式(7)、(8)、(18),由迭代方法求出抗力动力系数kh为0.616。2. When a protective structure is designed for explosion resistance, it is required that the structure be designed as a flexible rigid critical transition structure. The product ωt i of the circular frequency ω of the structure and the duration t i of the explosive dynamic load is 2.93, and the ductility ratio β is 1.26. When the damping ratio ξ is 5%, the exponential load attenuation coefficient a is 1.27, and the values of k ml and k ml are 0.66 and 0.78 respectively ("Underground Protective Structure", edited by Fang Qin and Liu Jinchun, ISBN 9787508470009 China Water Resources and Hydropower Press, 2010 Published in 2010), other parameters were calculated by other formulas, and substituted into formulas (7), (8) and (18), and the dynamic coefficient k h of resistance was obtained by iterative method, which was 0.616.

3、一种防护结构进行抗爆设计时,要求结构设计为柔性刚性临界过渡性结构,该结构自振圆频率ω与爆炸动荷载作用时长ti乘积ωti为1.46,延性比β为1.10,阻尼比ξ为10%时,指数型荷载衰减系数a为1.27,km-l与km-l取值分别为0.66、0.78(《地下防护结构》方秦,柳锦春编著,ISBN 9787508470009中国水利水电出版社,2010年出版),其它参数由其它公式分别计算得出,代入式(7)、(8)、(18),由迭代方法求出抗力动力系数kh为0.375。3. When a protective structure is designed for explosion resistance, it is required that the structure be designed as a flexible rigid critical transitional structure. The product ωt i of the structure's natural vibration circular frequency ω and the duration of the explosion dynamic load t i is 1.46, and the ductility ratio β is 1.10. When the damping ratio ξ is 10%, the exponential load attenuation coefficient a is 1.27, and the values of k ml and k ml are 0.66 and 0.78 respectively ("Underground Protective Structure", edited by Fang Qin and Liu Jinchun, ISBN 9787508470009 China Water Resources and Hydropower Press, 2010 Published in 2010), other parameters were calculated from other formulas, and substituted into formulas (7), (8) and (18), and the dynamic coefficient k h of resistance was obtained by iterative method, which was 0.375.

以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包在本发明范围内。The above descriptions are only preferred embodiments of the present invention, and are not intended to limit the present invention. Any modifications, equivalent replacements and improvements made within the spirit and principles of the present invention should be included within the scope of the present invention .

Claims (1)

1.含阻尼柔性刚性临界结构抗指数型爆炸荷载设计动力系数方法,其特征在于:所述的柔性刚性临界结构指的是:爆炸荷载作用真实时长ti数值恰为结构振动达到弹性振动最大值yT时刻,爆炸荷载卸载后,结构依靠惯性力开始进行塑性振动,在某一时刻tm,达到了结构弹塑性位移最大值ym1. The dynamic coefficient method for designing a flexible rigid critical structure with damping against exponential explosion loads is characterized in that: the flexible rigid critical structure refers to: the real duration t of the explosive load is just the value of the structural vibration reaching the maximum value of elastic vibration At time y T , after the explosion load is unloaded, the structure begins to vibrate plastically depending on the inertial force, and at a certain time t m , the maximum elastic-plastic displacement of the structure y m is reached; 根据爆炸对建筑结构的作用过程,将该过程分为爆炸荷载作用弹性阶段强迫振动和爆炸荷载卸载塑性阶段自由振动两个阶段;According to the action process of the explosion on the building structure, the process is divided into two stages: the forced vibration in the elastic stage of the explosion load and the free vibration in the plastic stage of the explosion load unloading; a、爆炸荷载作用弹性阶段强迫振动a. Forced vibration in the elastic stage of explosion load 在弹性阶段且在荷载作用时长范围0<t≤ti内,结构等效体系的动力方程为In the elastic stage and within the range of load duration 0<t≤t i , the dynamic equation of the structural equivalent system is
Figure QLYQS_1
Figure QLYQS_1
其中,t为结构振动时间参数,ti为爆炸荷载持续时长,Me为弹性阶段等效结构质量,Ce为弹性阶段等效结构阻尼,Ke为弹性阶段等效结构刚度,
Figure QLYQS_2
为结构等效体系振动加速度,/>
Figure QLYQS_3
为结构等效体系振动速度,y为结构等效体系振动位移,ΔPe(t)为结构承受的随时间t变化的爆炸动荷载,等效结构系数计算公式分别为:
Among them, t is the structural vibration time parameter, t i is the duration of the explosion load, M e is the equivalent structural mass in the elastic stage, C e is the equivalent structural damping in the elastic stage, K e is the equivalent structural stiffness in the elastic stage,
Figure QLYQS_2
is the vibration acceleration of the structural equivalent system, />
Figure QLYQS_3
is the vibration velocity of the structural equivalent system, y is the vibration displacement of the structural equivalent system, ΔP e (t) is the explosive dynamic load on the structure that changes with time t, and the calculation formulas of the equivalent structural coefficient are:
Figure QLYQS_4
Figure QLYQS_4
其中,m为真实结构每延米质量,l为真实结构跨长,ξ为真实结构阻尼比,K为真实结构刚度,kM为弹性阶段质量变换系数,kL为弹性阶段荷载变换系数,指数型爆炸动荷载为:Among them, m is the mass per linear meter of the real structure, l is the span length of the real structure, ξ is the damping ratio of the real structure, K is the stiffness of the real structure, k M is the mass conversion coefficient in the elastic stage, k L is the load conversion coefficient in the elastic stage, and the exponent The type explosion dynamic load is:
Figure QLYQS_5
Figure QLYQS_5
其中,ti为爆炸荷载作用时长,Δpm为爆炸荷载超压峰值,a为指数型荷载衰减系数,由无阻尼自振频率ω、含阻尼自振频率ωd、阻尼调整系数γ、爆炸荷载峰值作为静载时对应的静位移yst、弹性阶段质量变换系数与荷载变换系数之比值kM-L,各参数计算如下:Among them, t i is the duration of the blast load, Δp m is the overpressure peak value of the blast load, a is the exponential load attenuation coefficient, which consists of the undamped natural frequency ω, the damped natural frequency ω d , the damping adjustment coefficient γ, and the explosive load The peak value is used as the static displacement y st corresponding to the static load, and the ratio k ML of the mass conversion coefficient in the elastic stage to the load conversion coefficient k ML . The parameters are calculated as follows:
Figure QLYQS_6
Figure QLYQS_6
可知,结构承受爆炸荷载之前的初始位移及速度均为0,求解该微分方程后,可确定此阶段位移和速度表达式为:It can be seen that the initial displacement and velocity of the structure before bearing the explosion load are both 0. After solving the differential equation, the expression of displacement and velocity at this stage can be determined as:
Figure QLYQS_7
Figure QLYQS_7
Figure QLYQS_8
Figure QLYQS_8
在爆炸荷载作用结束卸载后,完成弹性阶段的振动时,对应ti时刻位移yi和速度vi为:After the explosion load is unloaded and the vibration in the elastic stage is completed, the displacement y i and velocity v i corresponding to time t i are:
Figure QLYQS_9
Figure QLYQS_9
Figure QLYQS_10
Figure QLYQS_10
b、爆炸荷载卸载塑性阶段自由振动b. Free vibration in the plastic stage of explosion load unloading 当结构振动时刻大于ti时刻,爆炸荷载消失,且结构刚进入塑性振动状态,此时结构为以yi及vi为初始条件的含阻尼塑性阶段的自由振动,即当ti<t<tm时,结构等效体系动力方程为When the structural vibration moment is greater than t i , the explosion load disappears, and the structure just enters the plastic vibration state. At this time, the structure is free vibration in the damped plastic stage with y i and v i as the initial conditions, that is, when t i <t< At t m , the dynamic equation of the structural equivalent system is
Figure QLYQS_11
Figure QLYQS_11
其中,me为塑性阶段等效结构质量,ce为塑性阶段等效结构阻尼,qm为结构塑性阶段的最大抗力,其计算公式为:Among them, m e is the equivalent structural mass in the plastic stage, c e is the equivalent structural damping in the plastic stage, q m is the maximum resistance force in the plastic stage of the structure, and its calculation formula is:
Figure QLYQS_12
Figure QLYQS_12
其中,km为塑性阶段质量变换系数,kl为塑性阶段荷载变换系数,解此运动微分方程,求出此阶段位移和速度解为:Among them, k m is the mass transformation coefficient in the plastic stage, and k l is the load transformation coefficient in the plastic stage. Solve this differential equation of motion to obtain the solution of displacement and velocity in this stage:
Figure QLYQS_13
Figure QLYQS_13
Figure QLYQS_14
Figure QLYQS_14
c、弹塑性阶段基于动力系数的延性比c. Ductility ratio based on dynamic coefficient in elastic-plastic stage 当结构振动至最大位移ym时,对应的时刻为tm,此时速度vm=0,代入(12)式,则:When the structure vibrates to the maximum displacement y m , the corresponding moment is t m , and the velocity v m =0 at this time, substituting into formula (12), then:
Figure QLYQS_15
Figure QLYQS_15
将tm带入到(11)中得出结构弹塑性振动最大位移为:Substituting t m into (11), the maximum displacement of the structural elastic-plastic vibration is:
Figure QLYQS_16
Figure QLYQS_16
令km-l为塑性阶段质量变换系数与荷载变换系数之比值,即km-l=km/kl,则Let k ml be the ratio of the mass conversion coefficient to the load conversion coefficient in the plastic stage, that is, k ml =k m /k l , then
Figure QLYQS_17
Figure QLYQS_17
由弹塑性阶段的抗爆设计,可知对于柔性结构、刚性结构临界区分结构,抗力动力系数kh和延性β分别为:From the anti-knock design in the elastoplastic stage, it can be known that for the flexible structure and the rigid structure critically divided structure, the resistance dynamic coefficient k h and the ductility β are respectively:
Figure QLYQS_18
Figure QLYQS_18
Figure QLYQS_19
Figure QLYQS_19
将(15)带入延性比公式(17)后,则得到After substituting (15) into the ductility ratio formula (17), we get
Figure QLYQS_20
Figure QLYQS_20
其中,yi与vi的表达式为(7)、(8)。Among them, the expressions of y i and vi are (7), (8).
CN202011304175.3A 2020-11-19 2020-11-19 Method for designing dynamic coefficient of anti-exponential explosion load of critical structure with damping flexible rigidity Expired - Fee Related CN112270032B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011304175.3A CN112270032B (en) 2020-11-19 2020-11-19 Method for designing dynamic coefficient of anti-exponential explosion load of critical structure with damping flexible rigidity

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011304175.3A CN112270032B (en) 2020-11-19 2020-11-19 Method for designing dynamic coefficient of anti-exponential explosion load of critical structure with damping flexible rigidity

Publications (2)

Publication Number Publication Date
CN112270032A CN112270032A (en) 2021-01-26
CN112270032B true CN112270032B (en) 2023-07-07

Family

ID=74340273

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011304175.3A Expired - Fee Related CN112270032B (en) 2020-11-19 2020-11-19 Method for designing dynamic coefficient of anti-exponential explosion load of critical structure with damping flexible rigidity

Country Status (1)

Country Link
CN (1) CN112270032B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20180009888A (en) * 2016-07-20 2018-01-30 연세대학교 산학협력단 Apparatus and system for measuring deformation of concrete structure under internal blast loading
CN109765025A (en) * 2018-12-25 2019-05-17 哈尔滨理工大学 Damage assessment method of RPC slab based on P-I curve under blast load
CN111753472A (en) * 2020-06-16 2020-10-09 中国人民解放军海军工程大学 Prediction method of ship motion response under underwater explosion considering damping effect

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20180009888A (en) * 2016-07-20 2018-01-30 연세대학교 산학협력단 Apparatus and system for measuring deformation of concrete structure under internal blast loading
CN109765025A (en) * 2018-12-25 2019-05-17 哈尔滨理工大学 Damage assessment method of RPC slab based on P-I curve under blast load
CN111753472A (en) * 2020-06-16 2020-10-09 中国人民解放军海军工程大学 Prediction method of ship motion response under underwater explosion considering damping effect

Also Published As

Publication number Publication date
CN112270032A (en) 2021-01-26

Similar Documents

Publication Publication Date Title
CN109765025B (en) Damage assessment method of RPC plate based on P-I curve under blast load
CN108416092B (en) Determination method for equivalent static load of RC beam with wave reinforcement
CN111753472B (en) Ship motion response prediction method under underwater explosion effect considering damping effect
CN110795827A (en) Vortex-induced vibration simulation method of elastically supported rigid cylindrical structures under the action of nonlinear energy wells
CN105183025A (en) Quasi-zero stiffness vibration isolation system and nonlinear feedback control method thereof
CN109826334B (en) Cylindrical tuned liquid column damper
CN104763765B (en) A piecewise linear vibration isolator with high static and low dynamic stiffness and its working method
CN112395678B (en) Dynamic coefficient method for linear explosion resistance design of flexible rigid critical structures with damping
CN112364425B (en) Dynamic coefficient method for design of rigid structures with damping against linear blast loads
CN110502787A (en) Optimal Design Method of Quasi-Zero Stiffness Shock Absorber
CN112270032B (en) Method for designing dynamic coefficient of anti-exponential explosion load of critical structure with damping flexible rigidity
CN106777813A (en) A kind of method for improving and controlling Multi-halfwaves Filters energy absorption efficiency
CN112417559B (en) Method for designing dynamic coefficient of anti-exponential explosion load of damping-containing flexible structure
CN106706190B (en) A Method of Determining the Equivalent Static Load of Shock Load
CN112417561B (en) Dynamic coefficient method for design of rigid structures with damping against exponential blast loads
CN106021630A (en) Structure/damping integrated design method suitable for novel spacecraft
CN114329736B (en) Solving method for residual deformation of low-damping rigid beam member under action of explosion load
CN106446386B (en) In mode energy method between mode stiffness of coupling a kind of confining method
CN114329735A (en) Solving method for residual deformation of mid-damping flexible beam members under blast load
Luo et al. Use of a novel rate-independent damper for performance improvement of a seismically isolated building structure
CN112883477B (en) Wall body reinforcing performance evaluation method and device based on bidirectional stress model
CN112395676A (en) Method for designing power coefficient of damping-containing flexible structure for resisting linear explosive load
Zhao et al. Studies on influence of geometrical parameters of an Acoustic Black Hole
CN117993159A (en) Numerical calculation method and system for structural fragmentation process under shock wave
CN103216566B (en) Vibration attenuating and impact resisting parameter design method adaptable to inertial navigation in supersonic flight environment

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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20230707