发明名称 基于等效单自由度体系的钢筋混凝土框架结构整体地震损伤水平评估方法
摘要 基于等效单自由度体系的钢筋混凝土框架结构整体地震损伤水平评估方法,涉及地震损伤评估技术领域,本发明为了实现在地震发生后利用结构获得的强震记录能够快速、准确的评估结构整体地震损伤水平。将多自由度体系结构等效为单自由度体系结构,将具有N层的多自由度体系等效为一个单自由度体系,单自由度体系质量M<sub>e</sub>为多自由度体系总质量;求解等效单自由度体系在相同地震动作用下的最大加速度反应;求解等效单自由度体系在相同地震动作用下延性系数;损坏指标计算;结构损伤评估。本发明克服了基于模态参数变化损伤评估存在较大误差以及基于物理参数耗时耗力且结果容易发散等缺点,实现简单、快速、有效的定量评估结构地震损伤水平。
申请公布号 CN103678937B 申请公布日期 2016.08.17
申请号 CN201310738272.7 申请日期 2013.12.29
申请人 中国地震局工程力学研究所 发明人 公茂盛;孙静;谢礼立
分类号 G06F19/00(2011.01)I 主分类号 G06F19/00(2011.01)I
代理机构 哈尔滨市松花江专利商标事务所 23109 代理人 杨立超
主权项 一种基于等效单自由度体系的钢筋混凝土框架结构整体地震损伤水平评估方法,其特征在于,所述方法的实现过程为:步骤一、将多自由度体系结构等效为单自由度体系结构:针对多自由度体系结构,对其进行等效简化,将具有N层的多自由度体系等效为一个单自由度体系,单自由度体系质量M<sub>e</sub>为多自由度体系总质量,等效高度为h<sub>e</sub>;<maths num="0001"><math><![CDATA[<mrow><msub><mi>M</mi><mi>e</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>m</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000942795820000011.GIF" wi="1052" he="127" /></maths>其中,m<sub>i</sub>为多自由度体系第i层的质量,M<sub>e</sub>为等效单自由度体系的质量,N为多自由度体系结构层数;步骤二、求解等效单自由度体系在相同地震动作用下的最大加速度反应:假定多自由度体系结构在地震中实际受到地震力分布如公式(2)所示:<maths num="0002"><math><![CDATA[<mrow><msub><mi>A</mi><mi>i</mi></msub><mo>=</mo><mfrac><mrow><msub><mi>G</mi><mi>i</mi></msub><msub><mi>H</mi><mi>i</mi></msub></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>G</mi><mi>j</mi></msub><msub><mi>H</mi><mi>j</mi></msub></mrow></mfrac><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>N</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000942795820000012.GIF" wi="1050" he="175" /></maths>其中,G<sub>i</sub>、G<sub>j</sub>分别为第i、j层重力,根据第i、j层结构集中质量计算;H<sub>i</sub>、H<sub>j</sub>分别为第i、j层距离地面高度;N为结构层数;假定地震中结构底层和顶层获得了强震记录,地震中顶层遭受的最大地震力用公式(3)估计:F<sub>N</sub>=m<sub>N</sub>a<sub>N max</sub>   (3)m<sub>N</sub>为顶层的质量,a<sub>N max</sub>为顶层最大加速度;根据公式(2)推导出,第i层的最大地震力F<sub>i</sub>:<maths num="0003"><math><![CDATA[<mrow><msub><mi>F</mi><mi>i</mi></msub><mo>=</mo><mfrac><msub><mi>A</mi><mi>i</mi></msub><msub><mi>A</mi><mi>N</mi></msub></mfrac><msub><mi>F</mi><mi>N</mi></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>...</mo><mi>N</mi><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000942795820000013.GIF" wi="1052" he="118" /></maths>假定多自由度体系和其等效单自由度体系的在相同地震动作用下的基底剪力和倾覆弯矩都相同,得到作用在单自由度体系的等效地震力,从V<sub>bS</sub>=V<sub>bM</sub>,可知等效地震力用公式(5)计算:<maths num="0004"><math><![CDATA[<mrow><msub><mi>F</mi><mi>e</mi></msub><mo>=</mo><msub><mi>V</mi><mrow><mi>b</mi><mi>S</mi></mrow></msub><mo>=</mo><msub><mi>V</mi><mrow><mi>b</mi><mi>M</mi></mrow></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>F</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000942795820000021.GIF" wi="1052" he="118" /></maths>其中,F<sub>e</sub>为相同地震动作用下等效单自由度体系ESDOF的地震力,V<sub>bS</sub>为ESDOF基底剪力,V<sub>bM</sub>为多自由度体系MDOF的基底剪力;根据公式(5)计算得到的等效地震力,等效单自由度体系ESDOF在同等地震动作用下最大加速度反应利用公式(7)计算:a<sub>Smax</sub>=F<sub>e</sub>/M<sub>e</sub>    (7)a<sub>Smax</sub>为相同地震动作用下等效单自由度体系ESDOF的最大加速度反应,a<sub>Smax</sub>在下一步中用来求解结构的等效延性系数;步骤三、求解等效单自由度体系在相同地震动作用下延性系数,利用非弹性等延性反应谱通过插值计算得出延性系数的具体过程为:首先根据结构基底地震动记录计算等延性反应谱,如公式(8)所示:S<sub>a</sub>=S<sub>a</sub>(T,ξ,μ)=|a(t,T,ξ,μ)|<sub>max</sub>(8)一条地震动的非弹性等延性谱S<sub>a</sub>为周期T,阻尼比ξ,及延性系数μ的函数;已知单自由度体系最大地震反应,通过结构的周期和阻尼比,可推得延性系数;通过计算得到基底地震动的两条等延性反应谱S<sub>aμ1</sub>(T,ξ,μ)和S<sub>aμ2</sub>(T,ξ,μ);S<sub>aμ1</sub>表示延性系数为μ<sub>1</sub>时的等延性反应谱,S<sub>aμ2</sub>表示延性系数为μ<sub>2</sub>时的等延性反应谱;延性系数μ<sub>1</sub>和μ<sub>2</sub>已知,当单自由度体系周期T已知时,可分别计算当单自由度体系延性达到μ<sub>1</sub>和μ<sub>2</sub>时的反应值,该值是单自由度体系在地震动作用下的最大反应,其计算表达式如公式(9)、(10)所示:<maths num="0005"><math><![CDATA[<mrow><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><mi>T</mi><mo>-</mo><msub><mi>T</mi><mn>1</mn></msub></mrow><mrow><msub><mi>T</mi><mn>2</mn></msub><mo>-</mo><msub><mi>T</mi><mn>1</mn></msub></mrow></mfrac><mo>(</mo><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><msub><mi>T</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>1</mn></mrow></msub><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo><mo>)</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000942795820000022.GIF" wi="1228" he="109" /></maths><maths num="0006"><math><![CDATA[<mrow><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>+</mo><mfrac><mrow><mi>T</mi><mo>-</mo><msub><mi>T</mi><mn>1</mn></msub></mrow><mrow><msub><mi>T</mi><mn>2</mn></msub><mo>-</mo><msub><mi>T</mi><mn>1</mn></msub></mrow></mfrac><mrow><mo>(</mo><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>2</mn></mrow></msub><mo>(</mo><msub><mi>T</mi><mn>2</mn></msub><mo>)</mo><mo>-</mo><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>2</mn></mrow></msub><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000942795820000031.GIF" wi="1252" he="110" /></maths>根据以上方程,当T=T<sub>1</sub>,简化为:S<sub>aμ1</sub>(T)=S<sub>aμ1</sub>(T<sub>1</sub>)     (11)S<sub>aμ2</sub>(T)=S<sub>aμ2</sub>(T<sub>1</sub>)     (12)或当T=T<sub>2</sub>时,简化为:S<sub>aμ1</sub>(T)=S<sub>aμ1</sub>(T<sub>2</sub>)     (13)S<sub>aμ2</sub>(T)=S<sub>aμ2</sub>(T<sub>2</sub>)        (14)根据方程(8)中非弹性反应谱的定义,周期为T的单自由度体系在某地震动作用下的最大反应,等于该条地震动延性系数为μ的等延性谱上周期T处对应的谱值:S<sub>aμ</sub>(T)=a<sub>Smax</sub>       (15)单自由度体系的最大反应a<sub>Smax</sub>已求得,单自由度体系的延性系数可以表示为插值公式(16)或(17):<maths num="0007"><math><![CDATA[<mrow><mi>&mu;</mi><mo>=</mo><msub><mi>&mu;</mi><mn>1</mn></msub><mo>-</mo><mfrac><mrow><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow></mfrac><mrow><mo>(</mo><msub><mi>&mu;</mi><mn>2</mn></msub><mo>-</mo><msub><mi>&mu;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>&mu;</mi><mn>1</mn></msub><mo>+</mo><mfrac><mrow><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>a</mi><mrow><mi>S</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub></mrow><mrow><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow></mfrac><mrow><mo>(</mo><msub><mi>&mu;</mi><mn>2</mn></msub><mo>-</mo><msub><mi>&mu;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>16</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000942795820000032.GIF" wi="1437" he="127" /></maths><maths num="0008"><math><![CDATA[<mrow><mi>&mu;</mi><mo>=</mo><msub><mi>&mu;</mi><mn>2</mn></msub><mo>-</mo><mfrac><mrow><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow></mfrac><mrow><mo>(</mo><msub><mi>&mu;</mi><mn>2</mn></msub><mo>-</mo><msub><mi>&mu;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>&mu;</mi><mn>2</mn></msub><mo>-</mo><mfrac><mrow><msub><mi>a</mi><mrow><mi>S</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow></msub><mo>-</mo><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow><mrow><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>1</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow><mo>-</mo><msub><mi>S</mi><mrow><mi>a</mi><mi>&mu;</mi><mn>2</mn></mrow></msub><mrow><mo>(</mo><mi>T</mi><mo>)</mo></mrow></mrow></mfrac><mrow><mo>(</mo><msub><mi>&mu;</mi><mn>2</mn></msub><mo>-</mo><msub><mi>&mu;</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>17</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000942795820000033.GIF" wi="1436" he="127" /></maths>在一条地震动作用下,当单自由度体系延性系数为μ时,其最大反应为a<sub>Smax</sub>,反过来讲,对于同一个固定的单自由度体系,当其加速度最大反应为a<sub>Smax</sub>时,其延性系数必然为μ;至此,求出了等效单自由度体系在基底地震动作用下的延性系数μ;步骤四、损坏指标计算过程:钢筋混凝土结构损伤指标采用Park&amp;Ang损伤指标模型,该模型为如公式(18)所示的双参数损伤指标:<img file="FDA0000942795820000034.GIF" wi="1078" he="134" />其中:DI<sub>P&amp;A</sub>为Park&amp;Ang损伤指数;μ<sub>m</sub>为地震动激励下单自由度体系延性系数,用最大反应位移除以屈服位移得到;μ<sub>u</sub>为单调荷载下极限延性系数;E<sub>h</sub>为地震动作用过程中体系滞回耗能;F<sub>y</sub>为体系的屈服强度;δ<sub>y</sub>为地震动作用下的体系屈服位移;β为一个无量纲常数;多自由度体系结构等效为单自由度体系结构后,利用公式(18)计算等效单自由度体系在基底地震动作用下达到第三步中估计的延性系数时的损伤指数;或利用公式(18)计算基底地震动的等延性损伤谱,根据周期T读出相应的损伤指数;步骤五、结构损伤评估:得到等效单自由度体系损伤指数DI<sub>P&amp;A</sub>后,根据表1中Park&amp;Ang给出的损伤指数与结构损伤水平之间的关系,可评估原多自由度体系结构的损伤水平,并给出结构可否修复的建议;表1 整体损伤指标和损伤水平对应关系<tables num="0001" wi="165"><table><tgroup cols="4"><colspec colname="c001" colwidth="16%" /><colspec colname="c002" colwidth="46%" /><colspec colname="c003" colwidth="21%" /><colspec colname="c004" colwidth="17%" /><tbody><row><entry morerows="1">损伤程度</entry><entry morerows="1">物理描述和表现</entry><entry morerows="1">损伤指数DI<sub>P&amp;A</sub></entry><entry morerows="1">结构状态</entry></row><row><entry morerows="1">倒塌</entry><entry morerows="1">建筑部分或全部倒塌</entry><entry morerows="1">≥1.0</entry><entry morerows="1">失效</entry></row><row><entry morerows="1">严重</entry><entry morerows="1">混凝土大量裂缝,钢筋露出、屈服</entry><entry morerows="1">0.4≤D&lt;1.0</entry><entry morerows="1">不能修复</entry></row><row><entry morerows="1">中等</entry><entry morerows="1">大量较大裂缝,薄弱构件混凝土脱落</entry><entry morerows="1">0.25≤D&lt;0.4</entry><entry morerows="1">可以修复</entry></row><row><entry morerows="1">轻微</entry><entry morerows="1">细小裂缝,部分柱混凝土脱落</entry><entry morerows="1">0.1≤D&lt;0.25</entry><entry morerows="1">简单修复</entry></row><row><entry morerows="1">完好</entry><entry morerows="1">零星的小裂缝</entry><entry morerows="1">D&lt;0.1</entry><entry morerows="1">无需修复</entry></row></tbody></tgroup></table></tables>
地址 150080 黑龙江省哈尔滨市南岗区学府路29号