基于能量原理的束筒结构简化计算.docVIP

基于能量原理的束筒结构简化计算 　　摘要：根据连续化原理，把束筒等效连续化为由各向异形板和角柱围成的等效实腹薄壁筒，记其剪切变形和纵向翘曲位移，假定了纵向位移函数，引入相应的对偶变量，建立了考虑剪力滞后影响的束筒结构受弯分析的哈密顿对偶求解体系，能方便地通过精细积分法求出高精度数值解。该方法是哈密顿力学在薄壁杆件结构弯曲分析中的应用，数学推导简单，且有成熟高效的数值算法，思路清晰、精度高、易于接受，对问题的求解提供了新的思路。计算结果表明，本方法具有较高的精度和适用性，可方便地用于束筒结构的初步计算。 　　关键词：束筒结构；单向弯曲；剪力滞后；哈密顿对偶求解体系；精细积分法； 　　 中图分类号： TU311文献标识码： A 文章编号： 　　束筒结构[1]是将几个筒体（通常为框筒）排放在一起，由楼盖将它们约束起来共同抵抗水平作用的结构体系。由于这种结构侧向刚度很大，从而可允许建筑物达到很高的高度。从结构上看，束筒结构的承载能力强、侧向刚度大，具有良好的延性和抗震性能；从建筑上看，束筒结构可以提供较大的使用空间，使建筑布置灵活，这是这些优点，束筒结构体系常被高层建筑[2]所采用。所以，对束筒结构的研究不仅具有一定的科学理论意义，而且也有广泛地工程应用价值。 　　1.基本假定和计算模型 　　目前束筒结构最有代表性的建筑是美国芝加哥的Sears大厦，这种结构具有空间刚度大、抗侧性能好、平面布置灵活等优点，是超高层建筑的最佳结构形式之一。本文采用连续化的数学模型，将九孔束筒的八榀框架模拟为图1中的八片等效均质正交板。从任一榀框架中取界于四个相邻反弯点之间的一标准十字元，并将它视为等效的矩形元。等效薄壁板的折算厚度为及折算剪切模量为，以上等效参数可由等效连续化法[3]求出。 　　在束筒结构单项弯曲分析中，我们可以采用如下基本假定： 　　（1）假定九孔束筒结构中个个模数筒之间是通过刚性楼板连接在一起，即同一标高处的沿方向的平动位移和绕轴的截面转角都是相同的。 　　（2）假定九孔束筒结构截面侧向内框架的纵向位移仍然采用材料力学里的平截面假定。 　　（3）假定九孔束筒结构截面法向外框架的纵向位移也采用材料力学里的平截面假定，只不过要乘以一个增大系数，这个增大系数可以通过内框架的纵向位移表示出。 　　（4）假定九孔束筒截面法向外框架的纵向位移中间为二次抛物线，两边为字形线性变化。 　　（5）假定九孔束筒截面法向内框架的纵向位移为线性变化且平行于断面，可以通过侧向内框架的纵向位移表示出。 　　（6）放弃了符拉索夫关于沿薄壁中心线剪应力为零的假定。 　　（7）忽略了环向应力和法向应力。 　　在束筒结构单项弯曲分析中，我们可以采用如下计算模型： 　　首先按照连续化原理将束筒结构等效连续化，如右图所示：可以求得侧向内框架的等效板厚，侧向外框架的等效板厚，法向内框架的等效板厚，法向外框架的等效板厚。可以求得侧向内框架的等效弹性模量和剪切模量、，侧向外框架的等效弹性模量和剪切模量、，法向内框架的等效弹性模量和剪切模量、，法向外框架的等效弹性模量和剪切模量、。角柱的等效弹性模量和截面面积分别为、。 　　2.束筒结构弯曲问题的位移场函数和拉格朗日函数 　　弯曲问题的位移场函数可表示为： 　　 　　 　　四个角柱的纵向位移函数 　　、、为法向外框架的纵向位移差函数。 　　通过能量原理得出体系的总的能量，可导出体系的拉格朗日函数[4]： 　　 　　3.束筒结构的哈密顿对偶求解体系及求解 　　将拉格朗日函数写成标准的矩阵形式 　　 　　式中：；； 　　 　　 　　其中的、、、、、、、、为简化矩阵而表示的。具体代表的数据如下所示： 　　 　　 　　此后即可按照基于两端边值问题精细积分法编制的程序求解结构体系的内力和位移。 　　4.计算实例 　　如图3-7所20层九孔束筒结构，层高,柱距，中柱截面为，角柱截面为，柱距为,梁截面均为楼板厚度取；钢筋混凝土弹性模量混凝土为。 　　表1侧移对比 　　Table1 Comparison of displacements 　　 　　 　　表2内力比较 　　Table2 Comparison of node forces 　　 　　 　　5.结论 　　通过算例对比可看出，本文的计算结果与有限元方法结果很接近。此方法以增加变量为代价，将问题降阶为一个一阶常微分方程组，从而回避了求解高阶微分方程不易得到解答的缺点，可通过精细积分法求出高精度的数值解。本文提出的束筒结构弯曲分析[5]的新方法，还适用于框筒、筒中筒等高层建筑结构的弯曲分析，具有一定的通用性，便于对此类实际工程的设计提供依据。 　　 　　参考文献 　　龚耀清，李珂.考虑地基