单层球面网壳结构抗震性能研究
第7卷第4期2001年12月
[文章编号]1006-6578(2001)04-0030-07
空 间 结 构SPATIALSTRUCTURES
Vol.7No.4Dec.2001
单层球面网壳结构抗震性能研究
陈军明, 陈应波, 吴代华
(武汉理工大学理学院结构工程与力学系, 湖北 武汉430070)
[摘 要] 本文研究了单层网壳结构地震响应的计算理论和计算方法。以K8型单层球面网壳为研究对象,得出了单层网壳在水平地震作用和竖向地震作用下的力学响应特性的一些重要结论。进一步针对地震响应的主要影响参数,如结构几何参数、边界刚度、阻尼比等,对单层球面网壳地震响应规律作了系统的研究,并提出了对单层网壳结构抗震设计有应用价值的结论和建议。
[关键词] 单层网壳;地震响应;计算理论;参数分析;抗震设计[中图分类号] TU311.3 [文献标识码] A
1 引 言
鉴于大跨建筑物的重要性,其地震作用效应引起了工程界的关注。近年来,许多学者对平板网架结构在地震作用下的反应及抗震计算方法进行了系统的研究,并制定了《网架结构设计与施工规范》。但国内外学者对单层网壳结构的研究主要集中在静力稳定性能研究,对其抗震性能研究进行较少。网壳结构的地震反应特征是否与网架结构相同?在各种结构参数影响下
网壳结构地震内力响应规律如何?对罕遇地震作用下结构的弹塑性反应如何计算?这些都是急待解决的问题。研究网壳结构的抗震性能是其在地震地区广泛应用的前提和基础。
2 运动平衡方程
2.1 运动方程
从单层网壳结构的合理传力方式来看,一般认为空间刚接是其最佳节点形式。故可将网壳
[收稿日期] 2001-06-30
[基金项目] 高等学校博士后流动站科研基金资助。
[作者简介] 陈军明(1966—),女,湖南人,博士,主要从事大跨结构抗震性能的研究。
2001年空 间 结 构第4期
简化为空间梁体系,并将荷载转化至结构的结点上。取空间梁单元,每个节点六个自由度,即U
T--={u v -w θx θy θz}。
地震作用下,结构的振动方程为
+Ku=-Mu +Cug(1)Mu
式中,M为结构质量矩阵;C为阻尼矩阵;K为结构刚度矩阵;u、 、 uu分别为相对于地面的位移矢量、速度矢量、加速度矢量; ug为地面地震运动加速度矢量。
将线位移和角位移分开排列,
则式(1)可改写成
Ms0
usCssCs usKssKsusMs0ugs
++=-(2) 0uCθsC
uKθsKu0ugθ,式中,us, θ,uθ分别为节点us, us分别是节点线位移矢量,线速度矢量和线加速度矢量;uu
角位移,角速度和角加速度矢量; ugs, ugθ分别是地面地震运动的线加速度和角加速度矢量。这
T--}T,u---里,us={u -v wθ={θx θy θz}。
进一步假定:(1)只考虑线位移加速度引起的惯性力,不考虑角加速度引起的惯性力;(2)
作用在质点上的阻尼力与节点相对线速度成正比,不考虑由角速度引起的阻尼力;(3)不考虑地面运动的转动分量。
根据以上基本假定,式(2)可写成
Ms us+Css us+Kssus+Ksθuθ=-Ms ugs
+KCθsusθsus+Kθθuθ=0由式(3)和式(4)消除uθ,并注意到Ksθ=Kθs,则有Ms us+Cs us+Ksus=-Ms ugs
-1-1
式中,s=Css-KsθKθθθs, Ks=Kss-KsθKθθKθsCC求解方程(5)可得到节点的线位移。再由式(4)可求得角位移计算公式
-1-1 s-Kuθ=-KθθCθsuθθKθsus
以上通过矩阵缩减和自由度凝聚,使6n×6n的矩阵降为3n×3n的矩阵。
2.2 质量矩阵和阻尼矩阵的计算
基于动力分析在运算上方便的考虑,采用团聚质量矩阵,荷载按静力等效原则作用于网壳节点,这些等效集中力转换为节点的等价集中质量,且在空间三个自由度方向具有相同惯性作用。采用近似Rayleigh阻尼考虑阻尼的影响,记阻尼矩阵C=α+βK,其中M和K分别为总M体质量矩阵和总体刚度矩阵,系数T、U的计算可根据结构第i、j振型的频率ki和kj及相应阻
i、aj确定。计算式为尼比a
2222
=2ikj(aikj-ajki)/(kj-ki), U=2(ajkj-aiki)/(kj-ki)(8)Tk
一般ki取结构的第一圆频率,kj取在结构抗震响应中的有效(贡献较大)最高频率,通常
(3)(4)(5)(6)(7)
取30Hz作为相应较高频率。本文中对各振型均取阻尼比a=0.02。
2.3 地震波的选用
地震波的波形对响应结果影响很大,其中对结构破坏有重要影响的因素为地震动强度、频谱特性和强震持续时间。故在选择强震记录时,场地条件应尽量接近;原则上应采用持续时间长的波;另外,当所选择的实际地震记录的加速度峰值与建筑地区设防烈度所对应的加速度峰值不一致时,可将实际地震记录的加速度按比例放大或缩小来加以修正。对应于不同设防烈度
2001年陈军明等:单层球面网壳结构抗震性能研究第4期
的多遇地震与罕遇地震的峰值加速度见建筑抗震设计规范[1],可采用小震时的峰值加速度作为地震输入以使结构处于弹性阶段工作;当要对结构作罕遇地震下的弹塑性分析时,则输入大震时峰值加速度进行计算。
2.4 动力方程求解
用来分析复杂的有限单元系统的直接积分法应该是无条件稳定的[2],否则,为了保持动力响应高阶分量的稳定性,所需的时间步长就太小了。Wilson-θ法在θ≥1.37时是无条件稳定的,但对结构的高阶振型具有较高的算法阻尼,对于频谱相当密集的网壳结构而言,用Wilson-θ法会产生较大的误差。当T≥0.5,U≥0.25(0.5+T)时,Newmark法是无条件稳定的,当T>
0.5时,算法引进了“算法阻尼”。据已有的理论分析结果和经验,本文建议:对于单层网壳结构,一种很好的求解技术是梯形法则,即T=0.5,U=0.25的Newmark法。运用Newmark法求解方程(5)可得到各时间步线位移,再用式(7)即可求得各时间步角位移。
2
3 单层球面网壳的地震响应[3]
文献[3]以跨度为40m,矢跨比为F/B=1/5,屋面荷载q=2.5kN/m,周边简支的K8型单层球面网壳为例(结构模型见图1),选取了适用于二类场地土的EL-Centro强震记录,按8度地震烈度,分别考虑水平向地震作用和竖向地震作用,进行多遇地震作用的弹性时程分析,得出了如下结论:
(1)水平地震作用下,主肋的动内力较小,而环向杆、斜向杆的动内力较大。竖向地震作用下,主肋、环向杆、斜向杆的动内力均较小。设防烈度为8度时,杆件竖向地震内力系数一般在1.1左右;而水平地震作用所产生的杆件内力已达静内力的20%——50%左右(主肋除外,10%左右)。总的来说,水平地震作用下的动响应较竖向地震作用下的动响应强烈得多,这与网架结构不一样(网架结构主要考虑竖向地震作用)[4]。因此,在8度、9度地区,单层球面网壳应考虑水平地震作用,9度地区还应该考虑竖向地震作用。
(2)采用地震内力系数的方法考虑地震作用,如平板网架结构的抗震计算,比较简单,但比较粗略。
2
图1 K8型单层球面网壳模型
由于网壳结构杆件的内力在地震作用下并不是按同一比例增加的,特别是对水平地震内力系数取统一值显然是不合适的,环杆与斜杆的地震内力系数值较大,而主肋杆的地震内力系数值
较小。因此,对于平面复杂或重要的大跨度网壳结构,建议采用时程法作专门的分析和验算。
4 单层球面网壳地震响应的主要影响参数
针对网壳的载荷、跨度、矢跨比、支座刚度以及阻尼等多种工况,仍以图1所示凯威特型网
2001年空 间 结 构第4期
壳为研究对象,重点研究水平地震作用,考虑8度设防。表1、2中均为最不利内力情况。(1)质量影响
屋面静荷载增大时,动内力迅速增大,近似呈线性关系,见图2。但杆件地震内力系数对荷载改变不敏感,见表1。
(2)跨度影响
表1 水平地震内力系数
F/B
B
q
C1
C2
C3
C4
C5
C6
C7
C8
C9
(m)(kN/m2)40
[***********][***********][***********]4040
1.01.52.02.51.01.01.01.01.01.01.01.52.02.51.01.01.01.01.01.01.01.52.02.51.01.01.01.01.01.0
主 肋
1/51/51/51/51/51/51/51/61/71/81/51/51/51/51/51/51/51/61/71/81/51/51/51/51/51/51/51/61/71/8
1.0111.0141.0131.0141.0091.0091.0071.0081.0071.005
1.0371.0431.0481.0511.0431.0361.0391.0461.0401.033
1.0541.0591.0551.0801.0491.0411.0601.0481.0441.038
1.0781.0701.0591.0881.0601.0651.0761.0781.0691.061
1.1241.1171.1081.1031.0891.0841.0871.1101.0851.074
1.1641.1511.1631.1861.1181.1191.0851.1691.1291.081
1.2171.1241.126
1.2901.141
1.300
环 杆1.1541.2471.3631.4911.5731.1811.2011.1931.1661.1561.1481.1871.172
1.2701.2051.3231.1671.1731.1871.1611.152
1.3021.3051.3451.2541.2161.2111.2421.197
1.4001.3861.4081.3031.2931.2811.3781.247
1.5281.5021.5341.5041.3531.3331.4351.273
1.5811.4931.376
1.5171.507
1.503
1.1451.1171.1761.2311.208
斜 杆
1.2211.2391.2961.3441.335
1.2811.2961.2701.2301.2161.2131.2721.2451.210
1.2771.2151.2801.1571.1731.1481.1531.1651.114
1.2611.1931.2761.1721.1931.1561.1901.1501.119
1.2891.2341.2701.2161.2221.1841.2471.1891.153
1.2961.2521.2431.2881.2441.2101.2601.1921.154
1.2971.2821.233
1.2861.276
1.287
2001年陈军明等:单层球面网壳结构抗震性能研究第4期
杆件地震内力系数对跨度改变不敏感,见表1。但杆件动内力随跨度增加而增加,见图3。(3)矢跨比影响
网壳环杆和斜杆水平地震内力系数随矢跨比增大而明显增大,对主肋影响不大,见表1。这表明随着网壳矢高的增大,其水平地震反应为主的特性更加突出
。
图2 设计控制杆件动内力随荷载的变化
图3 设计控制杆件动内力随跨度的变化
(4)支座影响
根据工程实际支承情况,应考虑网壳结构与下部支承结构协同工作。支承结构的弹性效应:考虑支承的弹性刚度,来模拟结构上下部共同作用的工况。根据实际网壳结构工程支承情况,一般采用周边切向约束和竖向约束,仅考虑结构法向弹性刚度。支承结构的惯性效应:上下两个结构体系之间的相互惯性作用比较复杂,为了近似地描述下部支承结构的惯性效应,引入支座的附加等价质量,并与弹性效应同时考虑。取网壳下部的柱和墙的一半质量,作为各相应支座的附加质量。处理方法是近似的,但在一定程度上能体现下部结构对上部结构的影响。
考虑三种分析方案:1.简支边界条件;2.弹性边界条件;3.考虑下部结构惯性效应。计算结果见表2。
表2 水平地震作用下,轴力响应随支座刚度的变化(8度,ξ=0.02,q=1kN/m2)
主 肋
杆件编号
方案1
C1C2C3C4C5C6
0.3341.2301.8482.6823.5632.104
方案20.5512.0954.3905.0284.3813.252
方案30.5851.9894.4295.1574.1012.544
方案15.2649.59713.47015.81014.060
方案23.3595.8266.5508.92714.320
方案34.0028.1709.79611.70015.040
3.5366.0568.81110.418.684
2.6835.1166.5456.9687.546
3.1576.2918.2328.7918.492
方案1
方案2
方案3
环 杆
斜 杆
注:方案2中,Ks=5.0×105N/m,方案3中,支座附加质量m=5kg。表中力的单位:kN
由表2可见,环杆和斜杆的水平地震内力随支承刚度减小而减小,而主肋内力有所增加。这主要是边缘环参与了工作,使内力变化趋于缓和。考虑弹性效应,并进一步考虑下部结构的惯性效应,环杆和斜杆动响应明显增大,主肋影响不大。结果表明:支座刚度的变化改变网壳地
2001年空 间 结 构第4期
震内力的分布,起控制作用的水平地震内力有明显改变。另外,网壳结构下部支承结构的附加惯性效应不容忽视。对于网壳这类大空间结构,其水平抗震设计宜考虑网壳与支承结构共同工作。
(5)阻尼影响
本文对阻尼比对结构动力响应的影响也作了研究,研究结果显示:阻尼比取值不同,对杆件中动内力有一定影响。阻尼比减小,地震反应有所增大,同样在El-Centro波作用下,当阻尼比为0.05时,主肋、环杆、斜杆的最不利动内力与静内力的比值一般分别在0.05,0.25,0.15左右;当阻尼比为0.02时,主肋、环杆、斜杆的最不利动内力与静内力的比值分别增大至0.10,0.35,0.30左右。动内力增长率随阻尼比的减小近似线性增大。在没有充分实测资料及理论分析之前,对于大跨网壳这种体系阻尼比取值不宜过高。
5 结 论
(1)设防8度地区应考虑水平地震,设防9度地区则水平和竖向地震均应考虑,且宜采用直接积分法进行地震响应计算。
(2)在单层球面网壳常用矢跨比1/5～1/8范围内,结构的水平地震响应较竖向地震响应大。矢跨比较大时,网壳以水平地震响应为主;随矢跨比减小,竖向地震响应呈增大的趋势。(3)工程设计中不应孤立地设计网壳,应整体地考虑网壳与下部结构的协同工作,尤其是水平抗震设计。另外,也可以充分利用边界影响来调整结构构件的内力,得出优化设计方案。
参 考 文 献
[1] GBJ11-89.建筑抗震设计规范[S].
[2] BatheK.J.FiniteElementProceduresinEngineeringAnalysis[M],Prentice-Hall,Inc.,Englewood
Cliffs,NewJersey,1982.
[3] 陈军明,陈应波,吴代华.单层球面网壳结构地震响应的动力时程分析[J].空间结构,1999,5(4):15-20.
[4] 陈军明,吴代华.样条配点法分析正交正放类网壳的动力响应[J].空间结构,1999,5(1):14-19.
STUDYONASEISMICBEHAVIOROFSINGLE-LAYERSPHERICALLATTICESHELLS
CHENJun-ming, CHENYing-bo, WUDai-hua
(Dept.ofStructuralEnginnering&Mechanics,WuhanUniversityofTechnology,Wuhan430070,China)[Abstract] Thisinvestigationisconcernedwiththeanalyticaltheoriesandcomputationalmethodsofseismicresponseforsingle-layersphericallatticeshells.BasedontheresultsofK8single-layersphericallatticeshells,
2001年陈军明等:单层球面网壳结构抗震性能研究第4期
seismicresponsesofsingle-layerlatticeshellsunderverticalandhorizontalearthquakeexcitationsareobtained.Furthermore,effectsofmajorstructuralparameters,suchasthegeometricalparameters,therigidityofsupportsanddampingratios,onthedynamicbehaviorofsingle-layersphericallatticeshellsareinvestigatedindetail.Theobtainedresultsarevaluableforthedesignofsingle-layerlatticeshellsinseismicregions.
[KeyWords] single-layerlatticeshell;seismicresponse;computationaltheories;parametricstudy;aseismicdesign