CN104866688A - 一种基于加速度敏感性分析的结构构件尺寸修正方法 - Google Patents
一种基于加速度敏感性分析的结构构件尺寸修正方法 Download PDFInfo
- Publication number
- CN104866688A CN104866688A CN201510314146.8A CN201510314146A CN104866688A CN 104866688 A CN104866688 A CN 104866688A CN 201510314146 A CN201510314146 A CN 201510314146A CN 104866688 A CN104866688 A CN 104866688A
- Authority
- CN
- China
- Prior art keywords
- phi
- structural elements
- wind
- formula
- unit
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Buildings Adapted To Withstand Abnormal External Influences (AREA)
Abstract
本发明涉及一种基于加速度敏感性分析的结构构件尺寸修正方法,用于对初始设计的超高层建筑的结构构件进行尺寸修正使其满足加速度响应限值,包括以下步骤:1)建立横向风加速度响应与结构构件尺寸的关系;2)获取初始设计的超高层建筑的结构参数,计算各结构构件的圆频率平方造价敏感性系数和质量归一化振型顶点分量造价敏感性系数;3)根据步骤2)的计算结果及风场环境参数计算各结构构件的横风向加速度响应造价敏感性系数;4)根据各结构构件的横风向加速度响应造价敏感性系数大小对结构构件进行相应修正。与现有技术相比,本发明可有效得到满足加速度限值的结构设计方案,且经济性好。
Description
技术领域
本发明涉及结构工程技术领域,尤其是涉及一种基于加速度敏感性分析的结构构件尺寸修正方法。
背景技术
现代超高层建筑不断往更高、更柔的方向发展,结构在强风作用下由于脉动风的影响产生的振动越来越显著。风荷载作用下,高层建筑结构振动加速度达到某一限值时,人们心理上就会出现不舒适的感觉。随着越来越多的对风荷载敏感的超高层建筑在各地的兴起,风振舒适度性能设计已经成为研究的热点问题,合理的抗风设计对结构安全性和经济性有重大的意义。
大量试验及实测数据表明,现代超高层建筑的横风向动力响应通常比顺风向的大,有时高达顺风向的数倍。超高层建筑横风向风致振动在风振舒适度设计中将起到控制作用。与顺风向风荷载相比,横风向风荷载形成机理要复杂得多,横风向风荷载不再符合准定常假定,横风向风荷载谱不能根据脉动风速谱直接给出。顾明、全涌对影响横风向一阶广义气动力谱因素进行分析,基于试验结果对横风向一阶广义气动力谱的拟合公式,并考虑结构气动阻尼影响推导出了一套高层建筑横风向风致响应的计算方法。这一拟合公式可准确给出横风向加速度响应与结构构件尺寸关系。
当结构振动不满足舒适度性能水平要求时,通常有两种方案降低加速度响应:动力修正方法修正结构刚度分布以改变结构动力特性(振型及周期)及布置阻尼器增加结构阻尼比。两种方案中布置阻尼器方案对于降低加速度较为高效但造价高,更适用于舒适度大幅度超限情况。而对于工程中大量存在的超高层小幅度超限情况(超限比例小于20%),结构动力修正(构件尺寸修正)方案由于不需增加装备费用,其增加的结构造价可能更低。
工程实际中增加阻尼器布置来降低风振的研究比较充分,缺少通过构件尺寸的修正来降低风振的对比研究。在加速度超限幅度不大时应对比两种方案的经济性,需要考虑在成本增加最少的情况下满足业主需求。另一方面,对于工程师来说,为使加速度响应降低到限值以下,可通过加速度响应与构件尺寸之间关系推导求解各构件成本敏感性的大小找到对加速度较为敏感的构件,最终根据成本敏感性大小修正构件尺寸,得到满足加速度限值的结构设计方案。工程实践中急需发展以顶点加速度为约束条件的基于加速度敏感性分析的结构构件尺寸修正方法,以降低结构顶点加速度响应到限值以下。
结构顶点横风向加速度响应除与结构一阶横风向自振频率相关外,还与结构质量归一化振型顶点分量有关。结构自振频率为整体指标,频率的修正代价相对较高。修正过程中可基于振型修正对横风向顶点加速度约束下的超高层建筑结构进行构件修正,在每一轮振型修正后,综合考虑结构自振周期变长对顶点加速度的不利影响。
发明内容
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种基于加速度敏感性分析的结构构件尺寸修正方法。
本发明的目的可以通过以下技术方案来实现:
一种基于加速度敏感性分析的结构构件尺寸修正方法,用于对初始设计的超高层建筑的结构构件进行尺寸修正使其满足加速度响应限值,包括以下步骤:
1)建立横向风加速度响应与结构构件尺寸的关系;
2)获取初始设计的超高层建筑的结构参数,计算各结构构件的圆频率平方造价敏感性系数和质量归一化振型顶点分量造价敏感性系数;
3)根据步骤2)的计算结果及风场环境参数计算各结构构件的横风向加速度响应造价敏感性系数;
4)根据各结构构件的横风向加速度响应造价敏感性系数大小对结构构件进行相应修正。
所述结构构件包括梁单元、柱单元和壳单元。
所述步骤2)中,圆频率平方造价敏感性系数的具体计算公式为:
所述结构构件为梁单元或柱单元时,第i个结构构件的圆频率平方造价敏感性系数SITi如下:
所述结构构件为壳单元时,第i个结构构件的圆频率平方造价敏感性系数SITi如下:
SITi=(ewT1LD/co2+3ewT2L3D3/co4)i
式中,E为材料弹性模量,G为材料剪切模量,co为构件成本,k为构件截面的高宽比,L为构件长度,D为构件高度,ebT1、ebT2、ebT3、ebT4、ebT5、ebT6为以下公式经积分获得的常量:
ewT1、ewT2为以下公式经高斯积分获得的常量:
式中,Lbi为第i个梁、柱单元的长度,Lwi为第i个壳单元的长度,Dwi为第i个壳单元的高度,FX、FY、FZ、MX、MY、MZ、F11、F22、F12、V13、V23、M11,M22和M12为modal工况作用下单元的内力,A为梁、柱单元的横截面积,AY和AZ分别为梁、柱单元的剪切面积,IX,IY和IZ分别为惯性力作用下梁、柱单元的扭转和弯曲惯性矩,B为壳单元的厚度,v为泊松比。
所述步骤2)中,质量归一化振型顶点分量造价敏感性系数的具体计算公式为:
所述结构构件为梁单元或柱单元时,第i个结构构件的质量归一化振型顶点分量造价敏感性系数SIφi如下:
所述结构构件为壳单元时,第i个结构构件的质量归一化振型顶点分量造价敏感性系数SIφi如下:
SIφi=(ewφ1LD/co2+3ewφ2L3D3/co4)i
式中,E为材料弹性模量,G为材料剪切模量,co为构件成本,k为构件截面的高宽比,L为构件长度,D为构件高度,ebφ1、ebφ2、ebφ3、ebφ4、ebφ5、ebφ6为以下公式经积分获得的常量:
ewφ1、ewφ2为以下公式经高斯积分获得的常量:
式中,Lbi为第i个梁、柱单元的长度,Lwi为第i个壳单元的长度,Dwi为第i个壳单元的高度,FX、FY、FZ、MX、MY、MZ、F11、F22、F12、V13、V23、M11,M22和M12为modal工况作用下单元的内力,fX、fY、fZ、mX、mY、mZ、f11、f22、f12、v13、v23,m11,m22和m12为虚拟工况作用下单元的内力,A为梁、柱单元的横截面积,AY和AZ分别为梁、柱单元的剪切面积,IX,IY和IZ分别为惯性力作用下梁、柱单元的扭转和弯曲惯性矩,B为壳单元的厚度,v为泊松比。
所述步骤3)中,横风向加速度响应造价敏感性系数表达式为:
式中,c由结构参数及风场环境参数确定,
G1(φ)=φ2
a1=β(B/UHfp)2-2(B/UHfp)2
a2=(B/UHfp)4
式中,H表示结构高度,B分别表示迎风面长度,f1表示模型的一阶自振周期,β为结构的振型修正系数,ζs为结构阻尼比,ζa为结构气动阻尼比,α为地面粗糙度系数,ωH为结构顶点处基本风压,UH为结构顶点处基本风速,gf为高斯因子,α、fp、Sp、β为曲线拟合得到的参数,f表示结构一阶自振频率,φ表示振型。
所述步骤4)中,进行结构构件修改时,增加横风向加速度响应造价敏感性系数大的结构构件的尺寸,减小横风向加速度响应造价敏感性系数小的结构构件的尺寸。
与现有技术相比,本发明基于能量法和虚功原理进行结构圆频率平方和质量归一化振型顶点分量关于构件造价的敏感性系数的求解,进而获得横风向加速度响应造价敏感性系数,根据该敏感性系数的大小进行结构构件尺寸修正,能够获得横风向加速度响应在限值内且造价经济的结构设计,方便可靠,且不需借助计算机优化算法,便于工程师理解。本发明适用于加速度响应超过规范要求的限值幅度不大(20%以下)的情况。
附图说明
图1为本发明实施例中结构抗侧力体系及伸臂桁架布置示意图;
图2为本发明实施例中结构剪力墙平面布置示意图;
图3为横风向加速度响应计算流程示意图;
图4为本发明横风向加速度响应关于构件造价敏感性计算流程示意图;
图5为实施例中结构各区伸臂及核心筒质量归一化振型顶点分量造价敏感性系数示意图;
图6为实施例中构各区伸臂及核心筒圆频率平方关于构件成本敏感性示意图;
图7为实施例中横风向加速度响应关于成本敏感性示意图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
如图1-图2所示,本实施例的背景工程为集商业、办公以及酒店为一体的综合性超高层建筑,塔楼的建筑高度为729m,结构高度为598m,共138层,整个建筑接近方形布置,建筑平面沿高度逐渐缩进,采用巨型框架—环带桁架—伸臂桁架—核心筒结构体系。塔楼沿高度设置9道环带桁架(1~9区分界处),以9道环带桁架为界,将塔楼沿高度分为9个区。塔楼50年重现期的风荷载,基本风压为0.45kN/m2,由于结构体型复杂,风荷载由风洞实验提供。结构参数表及风环境参数表见表1及表2。
表1结构参数表
H(m) | B(m) | D(m) | f1(Hz) | β | ζs(%) |
729 | 68 | 68 | 0.107 | 1.2 | 1.5 |
表2风环境参数表
风场类型 | ω0(kN/m2) | α | HT | ωH(kN/m2) | UH(m/s) |
C | 0.30 | 0.22 | 400 | 0.87 | 37.3 |
表1和表2中,H表示结构高度,B,D分别表示迎风面长度和宽度,f1表示案例模型的一阶自振周期,β为结构的振型修正系数,ζs为结构阻尼比,ω0为案例所在地基本风压,α为地面粗糙度系数,HT为梯度风高度,ωH为结构顶点处基本风压,UH为结构顶点处基本风速。
结构初始方案设计的剪力墙墙厚及伸臂桁架尺寸见表3及表4。
表3剪力墙初始墙厚
剪力墙 | W1E | W2E | W3E | W4E | W5E | W1F | W2F | W3F |
墙厚/mm | 900 | 800 | 750 | 550 | 450 | 1200 | 1100 | 1000 |
剪力墙 | W4F | W5F | W6F | W7F | W8F | W9F | W1I | W2I |
墙厚/mm | 900 | 750 | 750 | 600 | 550 | 500 | 900 | 800 |
剪力墙 | W3I | W4I | W5I | W6I | W7I | W8I | W9I | |
墙厚/mm | 800 | 700 | 650 | 650 | 550 | 450 | 400 |
注:W表示剪力墙,数字表示剪力墙所在区,E表示外墙,F表示翼墙,I表示内墙。
表4伸臂桁架初始截面尺寸
注:O表示伸臂桁架,数字表示伸臂桁架所在区,D表示斜腹杆,U表示上弦杆,L表示下弦杆。
超高层建筑结构为超静定结构,在外荷载作用下,修正变量的内力不仅与自身特性有关,也与修正变量之间的刚度比等有关,因此,一个约束条件常常与多个修正变量相关。结构的设计约束与多个修正变量的关系可用式(1)表示:
gi=fi(va,vb,vc,…) (1)
式中,gi表示第i个设计约束,va、vb和vc表示与设计约束gi相关的修正变量,fi表示第i个约束条件关于修正变量的函数。
根据布置位置、截面形式、材料类型和连接形式等的不同,各构件对约束条件的贡献程度、敏感性程度有很大差异。例如组成伸臂桁架的构件对增加结构刚度的贡献巨大,而组成环带桁架的构件对结构刚度作用甚微。因此,结构构件修正过程中,可根据约束条件关于修正变量的敏感程度对修正变量进行修正,增加对约束条件敏感构件的尺寸,减小对约束条件不敏感构件的尺寸,使修正变量间的材料得到合理的重分布。
约束条件关于结构构件造价的敏感程度可用敏感性系数表示,敏感性系数定义如下:
当Δcj趋向于0时,式(2)可表示为式(3):
SIi,j=dgi/dcj (3)
式中,SIi,j表示第i个约束条件关于修正变量j的敏感性系数,表示修正变量修正前结构第i个约束条件数值,表示修正变量修正后结构第i个约束条件数值,Δcj表示结构构件尺寸变化导致的结构构件成本变化。
如图4所示,本发明提供的一种基于加速度敏感性分析的结构构件尺寸修正方法,包括以下步骤:
1)建立横向风加速度响应与结构构件尺寸的关系;
2)获取初始设计的超高层建筑的结构参数,计算各结构构件的圆频率平方造价敏感性系数和质量归一化振型顶点分量造价敏感性系数;
3)根据步骤2)的计算结果及风场环境参数计算各结构构件的横风向加速度响应造价敏感性;
4)根据各结构构件的横风向加速度响应造价敏感性大小对结构构件进行相应修正。
1、质量归一化振型顶点分量造价敏感性系数
由虚功原理可知结构顶点质量归一化振型可由结构各构件modal工况下的节点力和虚拟工况(结构顶层施加单位力)下节点位移乘积总和来表示,第i个单元对结构质量归一化振型顶点分量贡献部分为
SITi=dφi/dV (5)
φ为结构的质量归一化振型顶点分量,N为单元数目,eφi为单元i在外荷载和虚拟荷载作用下的虚功,Fi为第i个单元在modal工况作用下的节点力向量,Δi为第i个单元在modal工况作用下的节点位移向量。
式中,E表示材料弹性模量,G表示材料剪切模量,v表示泊松比,dx表示沿梁、柱单元长度方向单位长度变化微分,dx1dx2表示沿壳单元长度和高度两个方向单位长度变化微分。
假设待修正模型为静定结构,修正变量尺寸发生变化时,单元内力保持不变,梁、柱单元和壳单元的内功eφi可分别表示为式(8)和式(9):
eφi=(ewφ1/B+ewφ2/(B)3)i(壳单元) (9)
式中,Lbi为第i个梁、柱单元的长度,Lwi为第i个壳单元的长度,Dwi为第i个壳单元的高度,FX,FY,FZ,MX,MY,MZ,F11,F22,F12,V13,V23,M11,M22和M12为modal工况作用下单元的内力,包括截面绕水平中性轴x,y及竖向中性轴z的扭转和弯曲惯性矩,A为梁、柱单元的横截面积,AY和AZ分别为梁、柱单元的剪切面积,IX,IY和IZ分别为惯性力作用下梁、柱单元的扭转和弯曲惯性矩,B为壳单元的厚度,ebφ和ewφ为对式(6)、(7)进行高斯积分得到的常量。
对于一般的建筑结构,框架单元通常具有矩形截面,如立柱和横梁。假设矩形截面的高宽比为k,截面属性AY,AZ,IX,IY和IZ可用单元体积vo(co表示构件成本,可由体积乘以成本系数得到)表示如下:
AY=AZ=5/6BH=5/6A=5vo/6L (10)
IX=0.2B3H=0.2A2L2/kL2=0.2vo2/kL2 (11)
IY=B3H/12=A2L2/12kL2=vo2/12L2k (12)
IZ=BH3/12=kA2/12=kvo2/12L2 (13)
式中,B为单元横截面的宽,H为单元横截面的高,L表示构件长度。
若修正模型为静定结构,结构约束条件的变化值等于第i个单元内能的变化值,可得到式(14):
dg/dcoi=dφ/dcoi=deφi/dcoi (14)
将式(15)对体积voi求导,质量归一化振型顶点分量造价敏感性系数SIφi可表示为式(16):
对于壳单元,假定长度L和高度D保持不变,则第i个壳单元的惯性力作用下的内能与单元成本co的关系可表示如下:
eφi=(ewφ1LD/co+ewφ2L3D3/co3)i (17)
将式(17)对单元成本coi求导,质量归一化振型顶点分量成本敏感性系数SIφi可表示为式(18):
SIφi=(ewφ1LD/co2+3ewφ2L3D3/co4)i (18)
2、圆频率平方造价敏感性系数
自振周期是建筑结构设计的重要参数,自振周期较大的建筑结构不仅会增加结构在地震作用和风荷载作用下强度破坏的可能性,还会导致正常使用状态下的变形和舒适性问题。为了降低长周期引起的有害作用,通常做法是在超高层建筑结构设计中对周期进行限制。对于体系和布置形式确定的结构,设计人员常通过修正单元尺寸减小结构周期。周期敏感性分析方法可区分周期对各单元的敏感性大小,高效地指导单元的材料重分布。
瑞利原理是计算振动系统固有频率的近似值,特别是最小固有频率(即基频)的上界的一个原理。对于一个在稳定平衡位置附近振动的保守系统,假设它以某一满足变形连续条件和位移边界条件的可能位移为振型作简谐振动,圆频率为ω。根据机械能守恒,系统最大势能与最大动能的关系可表示为式(19):
ω2φTMφ(最大动能)=φTKφ(最大动能) (19)
ω2=φTKφ/φTMφ=φTF/φTMφ=W/φTMφ (20)
将模态正规化,可以得到式(21):
ω2=4π2/T2=φTF/φTMφ=φTF=W (21)
式中,M为结构的质量矩阵,K为结构的刚度矩阵,φ为惯性力F作用下的第一阶振型,W为惯性力F所做的功,T是结构的第一阶自振周期。
对于一个无阻尼的混合结构,惯性力所做的功等于惯性力作用下所有构件的内功和,如式(22)所示:
假设待修正模型为静定结构,修正变量尺寸发生变化时,单元内力保持不变,梁、柱单元和壳单元的内功eTi可分别表示为式(25)和式(26):
eTi=(ewT1/B+ewT2/(B)3)i(壳单元) (26)
式中,FX,FY,FZ,MX,MY,MZ,F11,F22,F12,V13,V23,M11,M22和M12为惯性力作用下单元的内力,A为梁、柱单元的横截面积,AY和AZ为梁、柱单元的剪切面积,IX,IY和IZ为惯性力作用下梁、柱单元的扭转和弯曲惯性矩,B为壳单元的厚度,ebT和ewT为对式(23)、(24)进行高斯积分得到的常量。
若修正模型为静定结构,则在第i个单元尺寸发生变化时,仅影响第i个单元的内能,其它单元的内能保持不变。因此,结构约束条件的变化值等于第i个单元内能的变化值,可得到式(27):
dg/dcoi=dW/dcoi=deTi/dcoi (27)
由式(27)可知,圆频率平方关于单元i的成本的敏感性系数等于单元i在惯性力作用下的内功关于构件成本的敏感性系数。
将式(28)对单元成本coi求导,圆频率平方成本敏感性系数SITi可表示为式(29):
对于壳单元,假定长度L和高度D保持不变,则第i个壳单元的惯性力作用下的内能可用单元成本co表示如下:
eTi=(ewT1LD/co+ewT2L3D3/co3)i (30)
将式(30)对单元成本coi求导,圆频率平方造价敏感性系数SITi可表示为式(31):
SITi=(ewT1LD/co2+3ewT2L3D3/co4)i (31)
3、横风向加速度响应计算原理
在频域内求解结构运动方程,得到结构高度z处结构的加速度响应共振分量均方值近似计算公式:
该式的物理意义为:
表示力谱能量在频率上的分布,量纲为N2/Hz;
表示力谱能量,量纲为N2;
表示广义质量,量纲为kg;
表示结构顶点归一化振型值。
ζs1为结构阻尼比。
顾明、全涌对影响横风向一阶广义气动力谱因素做了分析,基于试验结果对横风向一阶广义气动力谱进行了拟合。并考虑结构气动阻尼影响推导出了一套高层建筑横风向风致响应的计算方法,对上式改进如下:
ζa1为结构气动阻尼比,用来考虑结构与空气之间的相互影响。
由此可知,结构顶点横风向加速度响应均方根由广义风荷载力谱的总能量与广义质量比值按结构振型值分配到顶点而得到。
通过对模拟风场中的方形及矩形截面高层建筑刚性模型的基底弯矩进行测量及拟合得到了横风向折算气动基底弯矩功率谱公式,与结构自振频率、高宽厚及风场环境相关。
通过数学转换,横风向一阶加速度响应可以表示如下:
U*=UH/(f1B)(42)
横风向加速度响应计算过程如图3所示,包括:
1)获取结构参数,包括单元高度、宽度、厚度、结构横向风基阶频率、振型函数、结构阻尼比、结构广义质量等;
2)获取风环境系统,包括基本风压、平均风速剖面指数、建筑顶部设计风压、梯度风高度、建筑顶部设计风速等;
3)计算横风向广义气动力的模态修正因子:
4)计算横风向折算气动基底弯矩功率谱
5)计算横风向气动阻尼值ζa;
6)计算结构的总加速度响应:
4、横风向加速度响应造价敏感性系数
假设结构顶层横风向加速度响应与结构一阶自振频率f及振型φ关系可简化如下式表示:
a(φ,f)=cG1(φ)G2(f) (45)
f表示结构自振频率,则顶点加速度关于构件成本敏感性系数可表示如下:
式中,c由结构参数及风场环境参数确定,对特定结构为常数,
G1(φ)=φ2(49)
a1=β(B/UHfp)2-2(B/UHfp)2 (51)
a2=(B/UHfp)4 (52)
式中,H表示结构高度,B分别表示迎风面长度,f1表示案例模型的一阶自振周期,β为结构的振型修正系数,ζs为结构阻尼比,ζa为结构气动阻尼比,α为地面粗糙度系数,ωH为结构顶点处基本风压,UH为结构顶点处基本风速,gf为高斯因子,α,fp,Sp,β为曲线拟合得到的参数。
可以看出,加速度响应关于构件成本敏感性由结构质量归一化振型关于构件成本敏感性和结构圆频率平方关于构件成本敏感性两部分构成,最终可计算得到加速度响应关于各构件成本敏感性的相对大小关系。
采用上述方法对本实施例超高层建筑进行结构构件尺寸修正,具体为:
1)横风向加速度响应计算
将结构参数表及风环境参数表中参数带入式(44)可知横风向加速度响应计算结果为0.26m/s2,超出规范限值0.25m/s24%,采用本发明提供的横风向加速度响应关于构件成本敏感性计算方法计算各构件敏感性大小,根据成本敏感性大小修正构件尺寸,得到使加速度响应降低到限值以下结构设计方案。
2)质量归一化振型顶点分量关于成本敏感性系数
结构各区伸臂及核心筒顶点质量归一化振型关于成本敏感性系数如图5所示。
3)圆频率平方关于成本敏感性系数
结构各区伸臂及核心筒圆频率平方关于成本敏感性系数如图6所示。
4)横风向加速度响应关于成本敏感性系数
结构各区伸臂及核心筒横风向加速度响应关于成本敏感性系数如图7所示。
5)根据成本敏感性系数的相对关系,修正该系数较大的剪力墙厚度,得到使加速度响应降低到限值以下结构设计方案。修正后的剪力墙厚度见表5。
表5修正后剪力墙厚度
剪力墙 | W1E | W2E | W3E | W4E | W5E | W1F | W2F | W3F |
墙厚/mm | 1100 | 1000 | 850 | 650 | 650 | 1400 | 1300 | 1200 |
剪力墙 | W4F | W5F | W6F | W7F | W8F | W9F | W1I | W2I |
墙厚/mm | 1000 | 800 | 800 | 800 | 750 | 700 | 750 | 700 |
剪力墙 | W3I | W4I | W5I | W6I | W7I | W8I | W9I | |
墙厚/mm | 700 | 650 | 650 | 650 | 650 | 650 | 600 |
自振周期从9.29s增加到9.41s导致加速度响应增加1.8%。而振型降低导致加速度响应降低5.8%,最终加速度响应从0.26降低到限值0.25m/s2(4%),满足规范的限值要求。修正后剪力墙体积增加了10418m3,需增加造价417万元,相比阻尼器方案增加造价更为经济。
Claims (6)
1.一种基于加速度敏感性分析的结构构件尺寸修正方法,用于对初始设计的超高层建筑的结构构件进行尺寸修正使其满足加速度响应限值,其特征在于,包括以下步骤:
1)建立横向风加速度响应与结构构件尺寸的关系;
2)获取初始设计的超高层建筑的结构参数,计算各结构构件的圆频率平方造价敏感性系数和质量归一化振型顶点分量造价敏感性系数;
3)根据步骤2)的计算结果及风场环境参数计算各结构构件的横风向加速度响应造价敏感性系数;
4)根据各结构构件的横风向加速度响应造价敏感性系数大小对结构构件进行相应修正。
2.根据权利要求1所述的基于加速度敏感性分析的结构构件尺寸修正方法,其特征在于,所述结构构件包括梁单元、柱单元和壳单元。
3.根据权利要求2所述的基于加速度敏感性分析的结构构件尺寸修正方法,其特征在于,所述步骤2)中,圆频率平方造价敏感性系数的具体计算公式为:
所述结构构件为梁单元或柱单元时,第i个结构构件的圆频率平方造价敏感性系数SITi如下:
所述结构构件为壳单元时,第i个结构构件的圆频率平方造价敏感性系数SITi如下:
SITi=(ewT1LD/co2+3ewT2L3D3/co4)i
式中,E为材料弹性模量,G为材料剪切模量,co为构件成本,k为构件截面的高宽比,L为构件长度,D为构件高度,ebT1、ebT2、ebT3、ebT4、ebT5、ebT6为以下公式经积分获得的常量:
ewT1、ewT2为以下公式经高斯积分获得的常量:
式中,Lbi为第i个梁、柱单元的长度,Lwi为第i个壳单元的长度,Dwi为第i个壳单元的高度,FX、FY、FZ、MX、MY、MZ、F11、F22、F12、V13、V23、M11,M22和M12为modal工况作用下单元的内力,A为梁、柱单元的横截面积,AY和AZ分别为梁、柱单元的剪切面积,IX,IY和IZ分别为惯性力作用下梁、柱单元的扭转和弯曲惯性矩,B为壳单元的厚度,v为泊松比。
4.根据权利要求2所述的基于加速度敏感性分析的结构构件尺寸修正方法,其特征在于,所述步骤2)中,质量归一化振型顶点分量造价敏感性系数的具体计算公式为:
所述结构构件为梁单元或柱单元时,第i个结构构件的质量归一化振型顶点分量造价敏感性系数SIφi如下:
所述结构构件为壳单元时,第i个结构构件的质量归一化振型顶点分量造价敏感性系数SIφi如下:
SIφi=(ewφ1LD/co2+3ewφ2L3D3/co4)i
式中,E为材料弹性模量,G为材料剪切模量,co为构件成本,k为构件截面的高宽比,L为构件长度,D为构件高度,ebφ1、ebφ2、ebφ3、ebφ4、ebφ5、ebφ6为以下公式经积分获得的常量:
ewφ1、ewφ2为以下公式经高斯积分获得的常量:
式中,Lbi为第i个梁、柱单元的长度,Lwi为第i个壳单元的长度,Dwi为第i个壳单元的高度,FX、FY、FZ、MX、MY、MZ、F11、F22、F12、V13、V23、M11,M22和M12为modal工况作用下单元的内力,fX、fY、fZ、mX、mY、mZ、f11、f22、f12、v13、v23,m11,m22和m12为虚拟工况作用下单元的内力,A为梁、柱单元的横截面积,AY和AZ分别为梁、柱单元的剪切面积,IX,IY和IZ分别为惯性力作用下梁、柱单元的扭转和弯曲惯性矩,B为壳单元的厚度,v为泊松比。
5.根据权利要求1所述的基于加速度敏感性分析的结构构件尺寸修正方法,其特征在于,所述步骤3)中,横风向加速度响应造价敏感性系数表达式为:
式中,c由结构参数及风场环境参数确定,
a1=β(B/UHfp)2-2(B/UHfp)2
a2=(B/UHfp)4
式中,H表示结构高度,B分别表示迎风面长度,f1表示模型的一阶自振周期,β为结构的振型修正系数,ζs为结构阻尼比,ζa为结构气动阻尼比,α为地面粗糙度系数,ωH为结构顶点处基本风压,UH为结构顶点处基本风速,gf为高斯因子,α、fp、Sp、β为曲线拟合得到的参数,f表示结构一阶自振频率,φ表示振型。
6.根据权利要求1所述的基于加速度敏感性分析的结构构件尺寸修正方法,其特征在于,所述步骤4)中,进行结构构件修改时,增加横风向加速度响应造价敏感性系数大的结构构件的尺寸,减小横风向加速度响应造价敏感性系数小的结构构件的尺寸。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510314146.8A CN104866688B (zh) | 2015-06-09 | 2015-06-09 | 一种基于加速度敏感性分析的结构构件尺寸修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510314146.8A CN104866688B (zh) | 2015-06-09 | 2015-06-09 | 一种基于加速度敏感性分析的结构构件尺寸修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104866688A true CN104866688A (zh) | 2015-08-26 |
CN104866688B CN104866688B (zh) | 2017-11-24 |
Family
ID=53912513
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510314146.8A Active CN104866688B (zh) | 2015-06-09 | 2015-06-09 | 一种基于加速度敏感性分析的结构构件尺寸修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104866688B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105404714A (zh) * | 2015-10-28 | 2016-03-16 | 北京海颐威工程技术有限公司 | 一种在建模软件中在位修改构件位置的方法和装置 |
CN106709460A (zh) * | 2016-12-28 | 2017-05-24 | 华南理工大学 | 一种高频底座力天平的动力校准方法 |
CN108051197A (zh) * | 2018-01-15 | 2018-05-18 | 江苏壹鼎崮机电科技有限公司 | 一种建筑抗震支吊架抗震性能异常预警方法 |
CN115017576A (zh) * | 2022-05-19 | 2022-09-06 | 中信建筑设计研究总院有限公司 | 考虑刚度影响的弯剪型高层结构刚重比计算的修正方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005047724A2 (en) * | 2003-11-14 | 2005-05-26 | University Of Maryland, Baltimore County | System and method for damping vibrations in elevator cables |
CN1851436A (zh) * | 2006-05-31 | 2006-10-25 | 汕头大学 | 大跨度屋盖和超高层建筑结构风振响应的检测计算方法 |
-
2015
- 2015-06-09 CN CN201510314146.8A patent/CN104866688B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005047724A2 (en) * | 2003-11-14 | 2005-05-26 | University Of Maryland, Baltimore County | System and method for damping vibrations in elevator cables |
CN1851436A (zh) * | 2006-05-31 | 2006-10-25 | 汕头大学 | 大跨度屋盖和超高层建筑结构风振响应的检测计算方法 |
Non-Patent Citations (5)
Title |
---|
CHUN-MAN CHAN 等: "Optimal lateral stiffness design of tall buildings of mixed steel and concrete construction", 《THE STRUCTURAL DESIGN OF TALL BUILDINGS》 * |
WINSTON CHAI 等: "Vibration control of super tall buildings subjected to wind loads", 《NON-LINEAR MECHANICS》 * |
全涌 等: "超高层建筑横风向气动力谱", 《同济大学学报》 * |
张建国 等: "典型超高层建筑横风向气动力谱的构成分析", 《北京工业大学学报》 * |
赵昕 等: "周期约束下超高层结构伸臂构件数值优化研究", 《建筑结构》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105404714A (zh) * | 2015-10-28 | 2016-03-16 | 北京海颐威工程技术有限公司 | 一种在建模软件中在位修改构件位置的方法和装置 |
CN105404714B (zh) * | 2015-10-28 | 2019-01-15 | 北京瀚方伟业工程技术有限公司 | 一种在建模软件中在位修改构件位置的方法和装置 |
CN106709460A (zh) * | 2016-12-28 | 2017-05-24 | 华南理工大学 | 一种高频底座力天平的动力校准方法 |
CN108051197A (zh) * | 2018-01-15 | 2018-05-18 | 江苏壹鼎崮机电科技有限公司 | 一种建筑抗震支吊架抗震性能异常预警方法 |
CN108051197B (zh) * | 2018-01-15 | 2019-11-26 | 江苏壹鼎崮机电科技有限公司 | 一种建筑抗震支吊架抗震性能异常预警方法 |
CN115017576A (zh) * | 2022-05-19 | 2022-09-06 | 中信建筑设计研究总院有限公司 | 考虑刚度影响的弯剪型高层结构刚重比计算的修正方法 |
CN115017576B (zh) * | 2022-05-19 | 2024-05-03 | 中信建筑设计研究总院有限公司 | 考虑刚度影响的弯剪型高层结构刚重比计算的修正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104866688B (zh) | 2017-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Lee et al. | Outrigger placement in tall buildings using topology optimization | |
JP5911077B2 (ja) | 空気と水と土骨格の連成計算装置および連成計算方法並びに連成計算プログラム | |
Tavakoli et al. | Seismic performance of outrigger-braced system based on finite element and component-mode synthesis methods | |
Venanzi et al. | Optimal design of an array of active tuned mass dampers for wind‐exposed high‐rise buildings | |
Li et al. | State-space approach for transverse vibration of double-beam systems | |
CN106650147A (zh) | 一种基于有界不确定性的连续体结构非概率拓扑优化方法 | |
Bobby et al. | Data-driven performance-based topology optimization of uncertain wind-excited tall buildings | |
CN105868477A (zh) | 基于层间位移角约束的结构抗震敏感性优化方法 | |
CN104866688A (zh) | 一种基于加速度敏感性分析的结构构件尺寸修正方法 | |
CN108595728A (zh) | 一种蜂窝材料的铺层等效有限元模型构建方法 | |
CN110717216B (zh) | 不规则波下带柔性气囊直升机横摇响应预报方法 | |
CN103615054A (zh) | 一种基于区格剪切变形的屈曲约束支撑布置方法 | |
CN104405054B (zh) | 一种设置黏滞消能器的结构设计方法 | |
CN103294857A (zh) | 以频响位移为优化目标的约束阻尼板的拓扑优化方法 | |
CN105447269A (zh) | 一种飞机在定常过载状态下的非结构质量配平方法 | |
Ruiz et al. | Eigenvector sensitivity when tracking modes with repeated eigenvalues | |
Li et al. | Linear statics and free vibration sensitivity analysis of the composite sandwich plates based on a layerwise/solid-element method | |
CN109858071B (zh) | 一种考虑剪力滞后作用的薄壁箱梁结构动力特性分析方法 | |
CN105676903A (zh) | 一种基于非概率可靠性优化的振动最优控制系统设计方法 | |
Aly et al. | Dynamics and Control of High‐Rise Buildings under Multidirectional Wind Loads | |
CN106372324A (zh) | 基于地震剪力系数约束的结构抗震敏感性优化方法 | |
Fu et al. | Wind resistant size optimization of geometrically nonlinear lattice structures using a modified optimality criterion method | |
Lee et al. | Investigation on material layouts of structural diagrid frames by using topology optimization | |
He et al. | New speedup algorithms for nonlinear dynamic time history analysis of supertall building structures under strong earthquakes | |
CN104820762A (zh) | 含有钢管混凝土柱的高层建筑框架结构优化设计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
EXSB | Decision made by sipo to initiate substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |