发明名称 一种框架结构损伤识别方法
摘要 一种框架结构损伤识别方法。首先,对损伤前的框架结构测试,对框架结构离散化初始模型进行模型修正,得到精确有限元理论模型,计算结构模态参数(ωA0,ΦA0);然后,对损伤后的结构进行测试,并进行模态分析,得到模态参数(ωt2,Φt2);建立直接解析法损伤识别方程,利用结构损伤前后的前n(n≤N)阶频率、振型变化值Δωr,ΔΦr(r=1,2,…,n)来求解单元损伤系数Di(i=1,2,…,n),达到同时识别整个结构存在的损伤的数量、位置和大小的目的。该方法适用于建筑框架结构的损伤识别,检测所需数据量少,效率高,具有较高的精度。
申请公布号 CN102043019A 申请公布日期 2011.05.04
申请号 CN201010513214.0 申请日期 2010.10.21
申请人 重庆大学 发明人 杨佑发;梁文广;徐典;李帅
分类号 G01N29/44(2006.01)I 主分类号 G01N29/44(2006.01)I
代理机构 重庆大学专利中心 50201 代理人 郭吉安
主权项 一种框架结构损伤识别方法,其步骤包括:(1)在损伤前框架结构的每层安装加速度传感器,在外激励作用下测试结构的时程响应,提取前8阶频率和振型数据,对框架结构离散化初始模型进行模型修正,得到精确有限元理论模型,计算结构模态参数(ωA0,ΦA0);(2)对损伤后的结构每层安装加速度传感器进行测试,在外激励作用下测试结构的时程响应并进行模态分析,得到模态参数(ωt2,Φt2);(3)求得固有频率的一阶灵敏度为: <mrow> <mfrac> <mrow> <mo>&PartialD;</mo> <msub> <mi>&omega;</mi> <mi>r</mi> </msub> </mrow> <mrow> <mo>&PartialD;</mo> <msub> <mi>D</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <msub> <mi>&omega;</mi> <mi>r</mi> </msub> </mrow> </mfrac> <mrow> <mo>(</mo> <msup> <msub> <mi>&Phi;</mi> <mi>r</mi> </msub> <mi>T</mi> </msup> <mfrac> <mrow> <mo>&PartialD;</mo> <msub> <mi>K</mi> <mi>d</mi> </msub> </mrow> <mrow> <mo>&PartialD;</mo> <msub> <mi>D</mi> <mi>i</mi> </msub> </mrow> </mfrac> <msub> <mi>&Phi;</mi> <mi>r</mi> </msub> <mo>-</mo> <msup> <msub> <mi>&Phi;</mi> <mi>r</mi> </msub> <mi>T</mi> </msup> <msubsup> <mi>&omega;</mi> <mi>r</mi> <mn>2</mn> </msubsup> <mfrac> <mrow> <mo>&PartialD;</mo> <mi>M</mi> </mrow> <mrow> <mo>&PartialD;</mo> <msub> <mi>D</mi> <mi>i</mi> </msub> </mrow> </mfrac> <msub> <mi>&Phi;</mi> <mi>r</mi> </msub> <mo>)</mo> </mrow> <mo>;</mo> </mrow>(4)求得振型对损伤系数的一阶偏导数为: <mrow> <mfrac> <mrow> <mo>&PartialD;</mo> <msub> <mi>&Phi;</mi> <mi>r</mi> </msub> </mrow> <mrow> <mo>&PartialD;</mo> <msub> <mi>D</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>=</mo> <munderover> <mi>&Sigma;</mi> <mrow> <mi>s</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <msub> <mi>&alpha;</mi> <mi>s</mi> </msub> <msub> <mi>&Phi;</mi> <mi>s</mi> </msub> </mrow>其中: <mrow> <msub> <mi>&alpha;</mi> <mi>s</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msubsup> <mi>&Phi;</mi> <mi>r</mi> <mi>T</mi> </msubsup> <mfrac> <mrow> <mo>&PartialD;</mo> <mi>M</mi> </mrow> <mrow> <mo>&PartialD;</mo> <msub> <mi>D</mi> <mi>i</mi> </msub> </mrow> </mfrac> <msub> <mi>&Phi;</mi> <mi>r</mi> </msub> </mtd> <mtd> <mi>s</mi> <mo>=</mo> <mi>r</mi> </mtd> </mtr> <mtr> <mtd> <munderover> <mi>&Sigma;</mi> <mrow> <mi>s</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>N</mi> </munderover> <mfrac> <mrow> <mn>2</mn> <msubsup> <mi>&omega;</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>&omega;</mi> <mi>r</mi> <mn>2</mn> </msubsup> </mrow> <mrow> <msubsup> <mi>&omega;</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <msubsup> <mi>&omega;</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msubsup> <mi>&omega;</mi> <mi>r</mi> <mn>2</mn> </msubsup> <mo>)</mo> </mrow> </mrow> </mfrac> <msup> <msub> <mi>&Phi;</mi> <mi>s</mi> </msub> <mi>T</mi> </msup> <mrow> <mo>(</mo> <msubsup> <mi>&omega;</mi> <mi>r</mi> <mn>2</mn> </msubsup> <mfrac> <mrow> <mo>&PartialD;</mo> <mi>M</mi> </mrow> <mrow> <mo>&PartialD;</mo> <msub> <mi>D</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>-</mo> <mfrac> <mrow> <mo>&PartialD;</mo> <msub> <mi>K</mi> <mi>d</mi> </msub> </mrow> <mrow> <mo>&PartialD;</mo> <msub> <mi>D</mi> <mi>i</mi> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> <msub> <mi>&Phi;</mi> <mi>r</mi> </msub> </mtd> <mtd> <mi>s</mi> <mo>&NotEqual;</mo> <mi>r</mi> </mtd> </mtr> </mtable> </mfenced> </mrow>(5)建立直接解析法损伤识别方程[S]{D}={Δf},其中:{Δf}为频率与振型差列阵,{D}为待求的单元损伤参数构成的列阵,[S]为频率及振型对单元损伤参数一阶偏导数构成的矩阵;(6)求解超定方程组[S]{D}={Δf},一次性得到全部损伤参数值,利用参数值{D}确定损伤的程度,通过查看损伤构件的编号,获得损伤位置,这样,达到同时识别整个结构存在的损伤的数量、位置和大小的目的。
地址 400044 重庆市沙坪坝区沙坪坝正街174号