matlab层间位移角的程序,地震作用下空间框架结构层间相对位移响应梯度和海赛矩阵的计算...

1引言结构动力响应梯度和海赛矩阵的计算是结构动力响应优化设计的关键技术之一。现有计算方法主要是有限差分法、解析法和半解析法[1]。根据各国规范的规定、震害经验和实验研究结果及工程实例分析,在地震作用下采用房屋建筑结构层间相对位移作为衡量结构变形能力从而判断结构是否满足建筑功能要求的指标是合理的。因此,在框架结构抗震优化设计中,框架结构的层间相对位移在结构优化设计的数学模型中常作为优化的目标函数或约束函数。然而,在地震作用下框架结构层间相对位移响应梯度和海赛矩阵的计算是一件棘手的事,可以采用有限差分法进行计算[2],但有限差分法通常用数值差商近似代替微商,导致结构重分析次数成倍的增加,降低了求解的效率,并且求解的精度与差分步长有关,差分步长的确定比较困难。因此,框架结构的抗震优化设计一般采用不需要求层间相对位移响应梯度和海赛矩阵的优化设计方法,如优化准则法[3,4]、进化算法[5]等。本文提出一种在地震作用下空间框架结构层间相对位移响应梯度和海赛矩阵的计算方法。在纽马克-法的基础上推导出在地震作用下空间框架结构层间相对位移响应梯度和海赛矩阵的计算公式,用Matlab语言编制了层间相对位移响应梯度和海赛矩阵的计算程序,实现了空间框架结构层间相对位移响应梯度和海赛矩阵的精确计算。2空间框架总的刚度矩阵、总的质量矩阵和总的阻尼矩阵对设计变量的一阶和二阶导数的计算21空间梁单元刚度矩阵对设计变量的一阶和二阶导数的计算空间梁单元如图1所示,单元编号为e,局部坐标系-xo-y的正向由结点i指向结点j,局部坐标系下单元结点位移向量-ex=[-iu-vi-wii-xi-yi-z-uj-vj-wjj-x-jy-jz]T,整体坐标系xoy下结点位移向量xe=[uiviwiixiyizujvjwjjxjyjz]T,其中u,v和w为线位移,为角位移。单元长度为le,单元截面为矩形截面(也可以是其他截面形式),矩形截面的宽be和高he称为单元e的设计变量。图1空间梁单元及单元的设计变量Fig.1Spatialbeamelementanddesignvariables在局部坐标系下单元的刚度矩阵为-Ke=Ebehele00000-Ebehele000000Ehebe3l3e000Ehebe32l2e0-Ehebe3l3e000Ehebe32le200Ebehe33le0-Ebehe32le2000-Ebehe33le0-Ebehe32le20000GIple00000-GIple0000-Ebehe32le20Ebehe33le000Ebehe32le20Ebehe36le00Ehebe32le2000Ehebe33le0-Ehebe32le2000Ehebe36le-Ebehele00000Ebehele000000-Ehebe33le000-Ehebe32le20Ehebe33le000-Ehebe32le200-Ebehe3l3e0Ebehe32l2e000Ebehe3l3e0Ebehe32le20000-GIple00000GIple0000-Ebehe32l2e0Ebehe36le000Ebehe32l2e0Ebehe33le00Ehebe32l2e000Ehebe36le0-Ehebe32l2e000Ehebe33le(1)式中Ip=heb3e12+beh3e12。等式(1)两边分别对单元设计变量求导一次、二次获得在局部坐标系下单元的刚度矩阵对单元设计变量的一阶导数-Kebe,-Kehe和二阶导数2-Kebe2,2-Kehe2,2-Kebehe和2-Kehebe。