发明名称 一种爆破振动多参数安全评价方法
摘要 本发明一种爆破振动多参数安全评价方法,采用小波包分解和重构方法将爆破振动信号分解为一系列简谐荷载,然后计算这些简谐荷载的能量并进行归一化处理,再计算简谐荷载作用下结构体系的动力响应,同时将归一化后的能量比例融入结构振动响应中,而后根据结构振动响应容许值可给出结构爆破振动安全的判据,进而评判结构的安全性。该判据是一个反映结构特性、爆破荷载幅值、多个优势频率、持续时间及能量比例等参数的多参数爆破振动安全判据,比现行的爆破振动安全允许标准更细致、合理和科学。
申请公布号 CN104182612A 申请公布日期 2014.12.03
申请号 CN201410355221.0 申请日期 2014.07.24
申请人 华侨大学 发明人 陈士海;常方强;张子华
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 泉州市文华专利代理有限公司 35205 代理人 陈智海
主权项 一种爆破振动多参数安全评价方法,其特征在于包括如下步骤:步骤1、通过小波包分解和重构技术将采集到的爆破振动信号表示为由不同振幅和不同频率的简谐波叠加,获得简谐荷载作用下结构体系的动态响应表达式:若单支重构信号V<sub>ij</sub>所在频段的幅值为<img file="FDA0000543170740000011.GIF" wi="98" he="101" />经过傅立叶变换后求得的该频段的主频为<img file="FDA0000543170740000012.GIF" wi="97" he="97" />则圆频率<img file="FDA0000543170740000013.GIF" wi="247" he="99" />将该支信号表示为简谐正弦波形式:<img file="FDA0000543170740000014.GIF" wi="359" he="96" />其中t为振动持续时间,i为小波包分解的层数,取正整数,j为小波包分解后的第i层的第j个节点,则爆破振动信号V可表示为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>V</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><msup><mn>2</mn><mi>i</mi></msup><mo>-</mo><mn>1</mn></mrow></munderover><msubsup><mi>V</mi><mi>ij</mi><mo>*</mo></msubsup><mi>sin</mi><mrow><mo>(</mo><msubsup><mi>&omega;</mi><mi>ij</mi><mo>*</mo></msubsup><mi>t</mi><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000543170740000015.GIF" wi="522" he="171" /></maths>结构体系的动态响应表达式可近似为:<img file="FDA0000543170740000016.GIF" wi="709" he="169" />式中η<sub>ij</sub>、<img file="FDA0000543170740000017.GIF" wi="63" he="85" />分别表示单支重构信号V<sub>ij</sub>作用下引起的结构动态响应的速度因子和相位角,这些值是小波分解时计算得到的<img file="FDA0000543170740000018.GIF" wi="94" he="96" />ω<sup>*</sup><sub>ij</sub>为单支重构信号V<sub>ij</sub>所在频段的幅值、主圆频率,β<sub>ij</sub>为稳态反应放大系数;步骤2、计算简谐载荷的振动能量并进行归一化处理:设单支重构信号V<sub>ij</sub>对应的振动能量为E<sub>ij</sub>,则有<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>E</mi><mi>ij</mi></msub><mo>=</mo><mo>&Integral;</mo><msup><mrow><mo>|</mo><msub><mi>V</mi><mi>ij</mi></msub><mrow><mo>(</mo><mi>t</mi><mo>)</mo></mrow><mo>|</mo></mrow><mn>2</mn></msup><mi>dt</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>|</mo><msub><mi>v</mi><mi>jk</mi></msub><mo>|</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000543170740000019.GIF" wi="525" he="138" /></maths>式中:v<sub>jk</sub>(k=1,2,...,n)为重构信号V<sub>ij</sub>离散点的幅值;归一化后的能量比例en<sub>ij</sub>表示为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>en</mi><mi>ij</mi></msub><mo>=</mo><mfrac><msub><mi>E</mi><mi>ij</mi></msub><mi>E</mi></mfrac><mo>=</mo><mfrac><msub><mi>E</mi><mi>ij</mi></msub><msup><mrow><mo>(</mo><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><msup><mn>2</mn><mi>i</mi></msup><mo>-</mo><mn>1</mn></mrow></munderover><msup><mrow><mo>|</mo><msub><mi>E</mi><mi>ij</mi></msub><mo>|</mo></mrow><mn>2</mn></msup><mo>)</mo></mrow><mfrac><mn>1</mn><mn>2</mn></mfrac></msup></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000543170740000021.GIF" wi="533" he="283" /></maths>步骤3、将归一化后的能量比例en<sub>ij</sub>融入结构振动响应中:将归一化后的能量比例en<sub>ij</sub>融入爆破振动信号V作用下单自由度结构的速度响应幅值中,获得归一化能量后的结构动态响应幅值<img file="FDA0000543170740000022.GIF" wi="106" he="96" /><maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msubsup><mi>V</mi><mi>s</mi><mo>*</mo></msubsup><mo>=</mo><mi>&lambda;</mi><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><msup><mn>2</mn><mi>i</mi></msup><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>en</mi><mi>ij</mi></msub><msub><mi>&eta;</mi><mi>ij</mi></msub><msub><mi>&beta;</mi><mi>ij</mi></msub><msubsup><mi>V</mi><mi>ij</mi><mo>*</mo></msubsup></mrow>]]></math><img file="FDA0000543170740000023.GIF" wi="442" he="156" /></maths>其中λ为一常系数,基本范围为1.2~1.6;步骤4、根据结构振动响应容许值<img file="FDA0000543170740000028.GIF" wi="102" he="80" />计算出表示爆破地面振动速度的容许幅值[V<sub>g</sub>],以此作为结构爆破振动安全的判据,根据标准或规程来评判结构的安全性:若V<sub>g</sub>为结构体系所在处爆破地震波原始信号引起的爆破地面振动速度幅值,则有:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msubsup><mi>V</mi><mi>s</mi><mo>*</mo></msubsup><mo>=</mo><mi>&lambda;</mi><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><msup><mn>2</mn><mi>i</mi></msup><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>en</mi><mi>ij</mi></msub><msub><mi>&eta;</mi><mi>ij</mi></msub><msub><mi>&beta;</mi><mi>ij</mi></msub><mrow><mo>(</mo><mfrac><msubsup><mi>V</mi><mi>ij</mi><mo>*</mo></msubsup><msub><mi>V</mi><mi>g</mi></msub></mfrac><mo>)</mo></mrow><mo>&times;</mo><msub><mi>V</mi><mi>g</mi></msub></mrow>]]></math><img file="FDA0000543170740000024.GIF" wi="593" he="178" /></maths>令<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msub><mi>D</mi><mi>ij</mi></msub><mo>=</mo><msub><mi>en</mi><mi>ij</mi></msub><msub><mi>&eta;</mi><mi>ij</mi></msub><msub><mi>&beta;</mi><mi>ij</mi></msub><mrow><mo>(</mo><mfrac><msubsup><mi>V</mi><mi>ij</mi><mo>*</mo></msubsup><msub><mi>V</mi><mi>g</mi></msub></mfrac><mo>)</mo></mrow><mo>,</mo></mrow>]]></math><img file="FDA0000543170740000025.GIF" wi="442" he="190" /></maths>为一无量纲系数,则<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><msubsup><mi>V</mi><mi>s</mi><mo>*</mo></msubsup><mo>=</mo><mi>&lambda;</mi><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><msup><mn>2</mn><mi>i</mi></msup><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>D</mi><mi>ij</mi></msub><msub><mi>V</mi><mi>g</mi></msub></mrow>]]></math><img file="FDA0000543170740000026.GIF" wi="325" he="154" /></maths>若<img file="FDA0000543170740000027.GIF" wi="100" he="79" />表示结构体系速度响应的容许幅值,则[V<sub>g</sub>]表示爆破地面振动速度的容许幅值,有<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><mo>[</mo><msub><mi>V</mi><mi>g</mi></msub><mo>]</mo><mo>=</mo><mfrac><mrow><mo>[</mo><msubsup><mi>V</mi><mi>s</mi><mo>*</mo></msubsup><mo>]</mo></mrow><mrow><mi>&lambda;</mi><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>0</mn></mrow><mrow><msup><mn>2</mn><mi>i</mi></msup><mo>-</mo><mn>1</mn></mrow></munderover><msub><mi>D</mi><mi>ij</mi></msub></mrow></mfrac></mrow>]]></math><img file="FDA0000543170740000031.GIF" wi="319" he="229" /></maths>式中,<img file="FDA0000543170740000032.GIF" wi="102" he="83" />可通过<img file="FDA0000543170740000033.GIF" wi="238" he="148" />求出,其中[σ]为建筑物结构介质的容许应力,ρ为介质密度;c<sub>p</sub>为纵波波速;步骤5、由此得到可以作为爆破振动安全评价的标准值:爆破地面振动速度幅值V<sub>g</sub>,该值隐含了振动频率<img file="FDA0000543170740000034.GIF" wi="96" he="95" />振动持续时间t、振动能量E<sub>ij</sub>、以及体现结构特性的结构动态响应的速度因子η<sub>ij</sub>。
地址 362000 福建省泉州市丰泽区城东华侨大学