发明名称 单层网壳结构在地震作用下的动力破坏成因判断方法
摘要 本发明公开单层网壳结构在地震作用下的动力破坏成因判断方法,它包括以下步骤:(1)完成增量步平衡迭代;(2)将每一个时间增量步结束时得到的杆件杆端力代入ISO圆钢管的失稳判别方程,以判断结构杆件是否发生失稳;(3)如果判断结构杆件没有发生失稳,继续按照步骤1的方法模拟该杆件计算杆件的两端受力,如果判断结构杆件发生了失稳,按照Marshall模型模拟该杆件;(4)对杆件失稳类型进行判断;(5)当步骤(1)中的结构非线性动力平衡方程组经过15次迭代仍不收敛时,可认为结构的动态承载力与地震作用不能维持平衡而发生动力破坏。采用本方法能够模拟地震作用下杆件可能反复经历的失稳-拉直与塑性铰形成-消失过程。
申请公布号 CN102127932A 申请公布日期 2011.07.20
申请号 CN201110022422.5 申请日期 2011.01.19
申请人 天津大学 发明人 丁阳;齐麟;李忠献
分类号 E04B1/32(2006.01)I 主分类号 E04B1/32(2006.01)I
代理机构 天津市北洋有限责任专利代理事务所 12201 代理人 王丽英
主权项 1.单层网壳结构在地震作用下的动力破坏成因判断方法,其特征在于它包括以下步骤:(1)完成增量步平衡迭代;求解过程如下:步骤101:采用Newmark时间积分法求解结构非线性动力平衡方程组,每次迭代得到杆端的位移增量Δu;步骤102:设杆件两端编号为I、J,I端平衡方程为<maths num="0001"><![CDATA[<math><mrow><msub><mi>F</mi><mi>iI</mi></msub><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>J</mi><mo>=</mo><mn>1</mn></mrow><mn>2</mn></munderover><munderover><mi>&Sigma;</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mn>6</mn></munderover><msubsup><mi>K</mi><mi>iI jJ</mi><mi>e</mi></msubsup><mrow><mo>(</mo><msub><mi>u</mi><mi>jJ</mi></msub><mo>-</mo><msubsup><mi>u</mi><mi>jJ</mi><mi>p</mi></msubsup><mo>)</mo></mrow></mrow></math>]]></maths>式中,F<sub>iI</sub>为I端截面力分量;<img file="FDA0000044474300000012.GIF" wi="79" he="56" />为单元弹性刚度矩阵中的分块矩阵,计算杆件弹性刚度矩阵时杆件的横向位移采用4次多项式插值函数表示,转动位移为横向位移插值函数对长度坐标的导函数,轴向位移采用2次多项式插值函数表示,扭转位移采用线性插值函数表示;u<sub>jJ</sub>为由步骤101中的位移增量累加得到的J端节点位移分量,<img file="FDA0000044474300000013.GIF" wi="54" he="56" />为前一个迭代结束时J端节点的塑性位移分量;步骤103:计算I端的屈服面函数Φ<sub>I</sub>为<maths num="0002"><![CDATA[<math><mrow><msub><mi>&Phi;</mi><mi>I</mi></msub><mo>=</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>N</mi><mi>xI</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>NxI</mi></msub></mrow><msub><mi>N</mi><mi>xu</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>M</mi><mi>xI</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>MxI</mi></msub></mrow><msub><mi>M</mi><mi>xu</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>M</mi><mi>yI</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>MyI</mi></msub></mrow><msub><mi>M</mi><mi>yu</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mfrac><mrow><msub><mi>M</mi><mi>zI</mi></msub><mo>-</mo><msub><mi>&alpha;</mi><mi>MzI</mi></msub></mrow><msub><mi>M</mi><mi>zu</mi></msub></mfrac><mo>)</mo></mrow><mn>2</mn></msup><mo>-</mo><mn>1</mn></mrow></math>]]></maths>式中,N<sub>xu</sub>为杆件仅受轴力作用时任意截面全截面屈服时的轴力值;M<sub>xu</sub>为杆件仅受扭矩作用时任意截面全截面屈服时的扭矩值;M<sub>yu</sub>为杆件仅受绕y轴的弯矩作用时任意截面全截面屈服时的弯矩值;M<sub>zu</sub>为杆件仅受绕z轴的弯矩作用时任意截面全截面屈服时的弯矩值;N<sub>xI</sub>为轴力;M<sub>yI</sub>、M<sub>zI</sub>分别为I端截面绕2个主轴的弯矩;M<sub>xI</sub>为I端截面的扭矩;α<sub>NxI</sub>、α<sub>MxI</sub>、α<sub>MyI</sub>、α<sub>MzI</sub>为背应力分量,Φ<sub>I</sub>≥0表示I端全截面屈服,形成塑性铰;步骤104:计算当前迭代步的杆端塑性位移增量<maths num="0003"><![CDATA[<math><mrow><mi>&Delta;</mi><msubsup><mi>u</mi><mi>I</mi><mi>p</mi></msubsup><mo>=</mo><mi>&Delta;</mi><msub><mi>&lambda;</mi><mi>I</mi></msub><mfrac><msub><mrow><mo>&PartialD;</mo><mi>&Phi;</mi></mrow><mi>I</mi></msub><mrow><mo>&PartialD;</mo><msub><mi>S</mi><mi>I</mi></msub></mrow></mfrac></mrow></math>]]></maths>式中,Δλ<sub>I</sub>为比例系数;S<sub>I</sub>=F<sub>I</sub>-α<sub>I</sub>,α<sub>I</sub>为I端截面背应力向量;F<sub>I</sub>为I端杆端力向量;步骤105:检查结构非线性动力平衡方程组迭代是否收敛,如果收敛,则本时间增量步计算完成;如果迭代不收敛,则将时间增量步减小一半重复进行步骤101-105;(2)将每一个时间增量步结束时由步骤102得到的杆件杆端力代入ISO圆钢管的失稳判别方程,以判断结构杆件是否发生失稳;(3)如果判断结构杆件没有发生失稳,则继续按照步骤1的方法模拟该杆件计算杆件的两端受力,如果判断结构杆件发生了失稳,则按照Marshall模型模拟该杆件;(4)对杆件失稳类型进行判断:若满足圆钢管的失稳判别方程时杆端屈服面方程Φ≥0,则杆件为第一种失稳类型,所述的第一种失稳类型为杆端截面产生塑性区或形成塑性铰后杆件长细比变大导致杆件失稳;若满足圆钢管的失稳判别方程时杆端屈服面方程Φ<0,则杆件为第二种失稳类型,所述的第二种失稳类型为杆端截面尚未产生塑性区或形成塑性铰,杆件因承受较大轴力作用达到失稳临界条件而发生失稳;(5)当步骤(1)中的结构非线性动力平衡方程组经过15次迭代仍不收敛时,则判定结构的动态承载力与地震作用不能维持平衡而发生动力破坏,对单层网壳结构在地震作用下的动力破坏成因进行判断:若破坏时结构中失稳杆件的类型为第一种,则单层网壳结构的动力破坏由杆端塑性铰集中出现导致结构局部变为机构引起;若破坏时结构中失稳杆件的类型为第二种,则单层网壳结构的动力破坏由杆件失稳导致的结构承载力下降引起;若破坏时两种失稳类型的杆件在结构中都存在,则单层网壳结构的动力破坏由失稳杆件引起结构承载力降低与杆端塑性铰使结构局部形成机构共同引起。
地址 300072 天津市南开区卫津路92号