CN111460545A - 一种高效计算超高层结构温度应变的方法和*** - Google Patents
一种高效计算超高层结构温度应变的方法和*** Download PDFInfo
- Publication number
- CN111460545A CN111460545A CN202010168266.2A CN202010168266A CN111460545A CN 111460545 A CN111460545 A CN 111460545A CN 202010168266 A CN202010168266 A CN 202010168266A CN 111460545 A CN111460545 A CN 111460545A
- Authority
- CN
- China
- Prior art keywords
- strain
- temperature
- sun
- super high
- rise structure
- 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
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000005855 radiation Effects 0.000 claims abstract description 58
- 239000006096 absorbing agent Substances 0.000 claims abstract description 14
- 238000004088 simulation Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000005259 measurement Methods 0.000 claims description 5
- 238000011156 evaluation Methods 0.000 claims description 3
- 239000012530 fluid Substances 0.000 claims description 3
- 230000002093 peripheral effect Effects 0.000 claims description 3
- 238000012546 transfer Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 abstract description 5
- 238000004458 analytical method Methods 0.000 description 9
- 230000004044 response Effects 0.000 description 9
- 238000005452 bending Methods 0.000 description 8
- 238000006073 displacement reaction Methods 0.000 description 8
- 238000012544 monitoring process Methods 0.000 description 5
- 229910000831 Steel Inorganic materials 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000036541 health Effects 0.000 description 4
- 238000013507 mapping Methods 0.000 description 4
- 238000007619 statistical method Methods 0.000 description 4
- 239000010959 steel Substances 0.000 description 4
- 230000009471 action Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012916 structural analysis Methods 0.000 description 2
- 238000002076 thermal analysis method Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
Abstract
本发明公开了一种高效计算超高层结构温度应变的方法和***,其中方法包括:根据太阳相对位置建立虚拟太阳,根据超高层结构建筑的几何形状建立有限元模型;对有限元模型进行网格划分得到多个结构单元,将虚拟太阳设定为辐射源,结构单元设定为吸收体,计算辐射源和每个吸收体之间的辐射系数,利用辐射系数判断结构单元为阴面单元或阳面单元;在阴面单元上施加第三类热边界条件时,太阳辐射强度取0,在阴面单元和阳面单元上施加第三类热边界条件后,对有限元模型进行仿真模拟,得到模拟温度场,基于温度应变与模拟温度场成正比的原则计算温度应变。本发明可以高效准确地计算温度应变、进而得到结构构件的真实应变变化、准确评估结构的安全状态。
Description
技术领域
本发明属于建筑结构分析计算技术领域,更具体地,涉及一种高效计算超高层结构温度应变的方法和***。
背景技术
结构健康监测技术已被广泛应用于大型土木结构的安全状态评估。从结构健康监测***中收集的数据不仅测量了结构的真实响应,更重要的是,还测量了环境引起的结构响应。温度应变作为最重要的环境响应,它的变化可能掩盖结构损伤引起的变化,从而影响安全评估的准确性。近年来,大量超高层结构安装了健康监测***,但是,关于如何计算温度应变的报道还较少。
目前,计算温度应变的方法主要有热传导有限元法和统计分析法。热传导有限元法首先通过建立热传导分析有限元模型,计算得到结构整体的温度场,然后建立结构分析有限元模型,将温度场转化为温度荷载,利用热-结构耦合分析来计算结构响应。在建立热传导分析有限元模型时,由于混凝土中钢筋对整体温度场的贡献可以忽略不计,因此可以不用考虑钢筋的影响,直接按照结构几何形状建立纯混凝土实体模型即可获得较精确的温度场。但是在建立结构分析有限元模型时,混凝土中的钢筋对结构整体刚度贡献较大,在建立有限元模型时必须考虑钢筋的影响,这对于体型巨大的超高层结构来说是非常困难的,不仅建模过程耗时耗力,最终的有限元模型也会因为单元数量巨大而无法完成计算。统计分析通常采用回归分析方法建立测点温度与结构响应之间的映射关系。统计分析的缺点是计算精度依赖于测点数据量,并且无法解释模型系数的物理意义。此外,统计分析是通过以往的响应数据来预测结构的响应,这依赖于过往响应数据的真实性。这种响应-响应映射容易出现***偏差。例如,在构建温度到应变的映射关系时,如果先前的应变数据包含结构损伤,那么所建立的温度应变定量关系也包含结构损伤,这不能反映温度应变之间原有的力学关系。目前亟需提出一种从温度场到温度应变的映射方法。
由此可见,现有技术存在不能准确计算温度应变、无法得到结构构件的真实应变变化、无法准确评估结构的安全状态的技术问题。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供了一种高效计算超高层结构温度应变的方法和***,由此解决现有技术存在不能准确计算温度应变、无法得到结构构件的真实应变变化、无法准确评估结构的安全状态的技术问题。
为实现上述目的,按照本发明的一个方面,提供了一种高效计算超高层结构温度应变的方法,包括以下步骤:
(1)根据太阳相对位置建立虚拟太阳,根据超高层结构建筑的几何形状建立有限元模型;
(2)对有限元模型进行网格划分得到多个结构单元,将虚拟太阳设定为辐射源,结构单元设定为吸收体,计算辐射源和每个吸收体之间的辐射系数,利用辐射系数判断结构单元为阴面单元或阳面单元;
(3)在阴面单元上施加第三类热边界条件时,太阳辐射强度取0,在阴面单元和阳面单元上施加第三类热边界条件后,对有限元模型进行仿真模拟,得到模拟温度场,基于温度应变与模拟温度场成正比的原则计算温度应变。
进一步地,所述太阳相对位置包括太阳高度角β、太阳方位角α和日地距离d,所述日地距离d为常数,所述太阳高度角β和太阳方位角α为:
进一步地,所述虚拟太阳的具体建立方式为:
根据太阳高度角β和太阳方位角α计算太阳在任意时刻的位置,建立虚拟太阳。
进一步地,所述步骤(2)中网格划分时,对超高层结构建筑中的柱、核心筒、梁和板采用不同的网格尺寸进行划分。
进一步地,所述步骤(2)中利用辐射系数判断的具体实现方式为:
若辐射系数为0,则结构单元被判定为阴面单元;若辐射系数不等于0,则结构单元被判定为阳面单元。
进一步地,所述步骤(3)还包括:
将气象参数转化为第三类热边界条件,所述气象参数包括环境温度Ta、风速v和太阳辐射强度I,所述第三类热边界条件包括表面传热系数hf和周围流体温度Tf,
hf=9.8+3.8v
其中,αs为热辐射吸收系数,0<αs<1,在阴面单元上施加第三类热边界条件时,太阳辐射强度I为0。
进一步地,所述步骤(3)还包括:
根据实测应变与模拟温度场之间的比例关系得到比例系数,对模拟温度场进行积分后乘以比例系数得到温度应变。
进一步地,步骤(3)还包括:
利用模拟温度场和温度应变公式计算温度应变,所述温度应变公式为:
其中,ε为待计算截面的温度应变,A表示待计算截面的截面面积,T(x,y)表示待计算截面上坐标为(x,y)处的模拟温度场,E代表弹性模量,αT表示热膨胀系数,H为待计算截面所处的高度,k1为待计算截面因周边构件约束而产生的轴向弹簧刚度。
进一步地,所述方法计算得到的温度应变应用于评估超高层结构建筑的安全状态,具体地:
从超高层结构建筑的实测应变中剔除温度应变,得到超高层结构建筑的真实应变,利用真实应变评估超高层结构建筑的安全状态。
按照本发明的另一方面,提供了一种高效计算超高层结构温度应变的***,包括以下模块:
模型建立模块,用于根据太阳相对位置建立虚拟太阳,根据超高层结构建筑的几何形状建立有限元模型;
阴面判断模块,用于对有限元模型进行网格划分得到多个结构单元,将虚拟太阳设定为辐射源,结构单元设定为吸收体,计算辐射源和每个吸收体之间的辐射系数,利用辐射系数判断结构单元为阴面单元或阳面单元;
模拟计算模块,用于在阴面单元上施加第三类热边界条件时,太阳辐射强度取0,在阴面单元和阳面单元上施加第三类热边界条件后,对有限元模型进行仿真模拟,得到模拟温度场,基于温度应变与模拟温度场成正比的原则计算温度应变。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,能够取得下列有益效果:
(1)世界高层建筑与都市人居学会(CTBUH)将高度超过300米的建筑定义为超高层结构建筑。超高层结构建筑在进行热传导分析时非常困难,而本发明根据太阳相对位置建立虚拟太阳,根据超高层结构建筑的几何形状建立有限元模型;建立模型后自动进行阴面单元或阳面单元的判断,在阴面单元和阳面单元上施加第三类边界条件后,进行仿真模拟,由此模拟得到的模拟温度场更加高效、准确,现有技术中仅仅只有单点温度与应变之间的关系,本申请创造性的利用温度应变与模拟温度场成正比这一原则计算温度应变,可以高效准确地计算温度应变、进而得到结构构件的真实应变变化、准确评估结构的安全状态。
(2)本发明在网格划分时,对超高层结构建筑中的柱、核心筒、梁和板采用不同的网格尺寸进行划分,以实现计算精度和计算效率之间的平衡。
(3)如果辐射系数为0,则表明辐射源和吸收体单元之间有遮挡,该结构单元被判定为阴面单元;如果辐射系数不等于0,则表明辐射源和吸收体单元之间没有遮挡,该结构单元被判定为阳面单元,在阴面单元上施加第三类热边界条件时,太阳辐射强度取0,在阴面单元和阳面单元上施加第三类热边界条件,可以提高仿真模拟结果的准确率。
(4)在监测过程中布置的传感器,一般只能测得单个点的温度和应变,但是通过观测测量数据发现,实测的温度和应变并不是线性相关,因而很难通过测量数据获取温度应变。以往的文献也没有涉及截面温度场与应变关系的公式,只有单点应变-单点温度公式。本发明创造性发现温度应变与模拟温度场成正比,基于此可以通过比例系数和公式两种方式计算温度应变,根据实测应变与模拟温度场之间的比例关系得到比例系数进而计算温度应变较为方便、快捷,根据公式计算温度应变较为准确
附图说明
图1是本发明实施例提供的一种高效计算超高层结构温度应变的方法的流程图;
图2是本发明实施例1提供的武汉长江航运中心三维实体模型;
图3(a)是本发明实施例1提供的柱的网格划分示意图;
图3(b)是本发明实施例1提供的核心筒的网格划分示意图;
图3(c)是本发明实施例1提供的梁和板的网格划分示意图;
图4(a)是本发明实施例1提供的环境温度图;
图4(b)是本发明实施例1提供的风速图;
图4(c)是本发明实施例1提供的太阳辐射强度图;
图5是本发明实施例1提供的温度应变的计算与比较图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
如图1所示,一种高效计算超高层结构温度应变的方法,包括以下步骤:
(1)根据太阳相对位置建立虚拟太阳,根据超高层结构建筑的几何形状建立有限元模型;
(2)对有限元模型进行网格划分得到多个结构单元,将虚拟太阳设定为辐射源,结构单元设定为吸收体,计算辐射源和每个吸收体之间的辐射系数,利用辐射系数判断结构单元为阴面单元或阳面单元;
(3)在阴面单元上施加第三类热边界条件时,太阳辐射强度取0,在阴面单元和阳面单元上施加第三类热边界条件后,对有限元模型进行仿真模拟,得到模拟温度场,基于温度应变与模拟温度场成正比的原则计算温度应变。
步骤(1)包括:
所述太阳相对位置包括太阳高度角β、太阳方位角α和日地距离d,所述日地距离d为常数,所述太阳高度角β和太阳方位角α为:
其中,表示超高层结构建筑所处的纬度,τ为辐射历计时角,单位为度,τ以年为周期变化,以春分日(3月21日)0:00:00时刻为计时零点,τ每小时增加15度,τ∈[0,131400);k为整数,k=0,1,2,...,364。赤纬角为0的时候就是春分。这是因为太阳辐射历时角是在赤纬角的基础上改进的。
根据太阳高度角β和太阳方位角α计算太阳在任意时刻的位置,建立虚拟太阳。在ANSYS平台上建立虚拟太阳和有限元模型。
步骤(2)包括:
网格划分时,对超高层结构建筑中的柱、核心筒、梁和板采用不同的网格尺寸进行划分。
若辐射系数为0,则结构单元被判定为阴面单元;若辐射系数不等于0,则结构单元被判定为阳面单元。
步骤(3)还包括:
将气象参数转化为第三类热边界条件,所述气象参数包括环境温度Ta、风速v和太阳辐射强度I,所述第三类热边界条件包括表面传热系数hf和周围流体温度Tf,
hf=9.8+3.8v
其中,上式中,v为风速,单位:m/s;Ta为环境温度,单位:℃;αs为热辐射吸收系数,0<αs<1;I为太阳辐射强度,单位:W/m2。
在阴面单元上施加第三类热边界条件时,太阳辐射强度I取为0。在ANSYS平台上定义分析类型为瞬态热分析,进行温度场求解。
推导悬臂柱模型的温度-应变理论公式。假定超高层结构在日照作用下为弯曲变形,将结构简化为悬臂柱模型,该模型描述为,柱底面被完全约束,柱高为H,悬臂柱仅受温度作用,截面内的温度场为T(x,y)。现推导温度和应变之间的理论关系。当柱顶被完全约束住时,截面内任一点(x,y)的应力为:
σ(x,y)=EαTT(x,y)
上式中E代表弹性模量;αT表示热膨胀系数。
此时,需要在顶部施加的轴力为:
FN=∫σ(x,y)dA=∫EαTT(x,y)dA
除了轴力外,由不均匀分布的应力还会产生弯矩。微面dA对x轴的弯矩为EαTT(x,y)dA,沿整个截面积分可得对x轴的总弯矩为:
Mx=∫EαTT(x,y)ydA
类似的,对y轴的总弯矩为:
My=∫EαTT(x,y)xdA
现在释放顶部的约束,为了保证悬臂柱的变形模态和初始状态相同,需要在顶部施加外力FN,Mx和My。假设截面变形满足平截面假定,则由轴力和弯矩引起的任意点(x0,y0)处的弯矩为:
上式中Ix和Iy分别表示截面对x和y轴的惯性矩。
推导广义悬臂柱模型的温度-应变理论公式。广义悬臂柱模型描述为在悬臂柱的基础上,顶部截面还被轴向弹簧、x和y方向的转动弹簧约束。设轴向弹簧刚度、x和y方向的转动弹簧刚度为k1,k2和k3。在截面温度场T(x,y)的作用下,假设释放顶部的弹簧约束,则由温度场T(x,y)造成的顶部等效力为:
FN=∫EαTT(x,y)dA
Mx=∫EαTT(x,y)ydA
My=∫EαTT(x,y)xdA
然而,顶部弹簧会约束部分顶部位移,约束的位移将会产生柱截面内力。根据内力和外力的平衡关系,可得:
Mx-k2θx=EIxδ″x
My-k3θy=EIyδ″y
上式中u0为顶部线位移,θx和θy为顶部转角;δx和δy分别表示x和y方向的顶部水平位移。
现在考虑顶部横向位移和转角之间的关系。在弯曲型变形模式中,侧向位移的一阶导数是转角,转角的一阶导数是弯矩。因此,转角是高度的一次函数,侧向位移是高度的二次函数。我们知道,在该模型中底部和顶部的x向旋转角分别为0和θx,则任意高度z处的转角方程表示为:
对上式两侧求导,可得:
由上式可知δ″x在高度方向保持为常数。将上式带入内外力平衡等式得:
类似的可以计算k3。因此,k1,k2和k3的表达式为:
k1,k2和k3为刚度系数,它们均为常数。根据上面三个等式可以推导位移和温度之间的关系:
由平截面假定,有:
将u0,θx和θy带入上式得:
因为FN,Mx和My为温度场T(x,y)的函数,所以任一点(x0,y0)处的应变为:
如果令,
则有,
ε(x0,y0)=αT(TA+Tex·y0+Tey·x0)
上式中,TA,Tex和Tey分别称为面平均温度、x和y方向的梯度温度,其中梯度温度表征了温度的不均匀程度。对比悬臂柱模型和广义悬臂柱模型的温度-应变公式,可发现两公式的表达形式很相似。如果令k1=k2=k3=0,广义悬臂柱模型公式退化成悬臂柱模型公式,因此广义悬臂柱温度-应变公式更通用,后面将应用该公式进行温度应变的计算。
利用模拟温度场和温度应变公式计算温度应变。在实际高层结构中,线弹簧和转动弹簧刚度一般难以求解,可将温度应变公式简化为:
上式表明温度应变大约与温度场成正比。温度应变的计算方法为,首先利用有限元热传导分析计算结构整体的温度场,然后计算待研究截面的面平均温度TA,最后根据上述公式计算温度应变。本发明利用推导的温度应变理论公式来计算温度应变,避免建立复杂的结构分析有限元模型,极大的提高了温度应变的计算效率,在结构健康监测领域具有较高的工程应用价值。
因此,本发明基于温度应变与模拟温度场成正比的原则计算温度应变,具体地:
根据实测应变与模拟温度场之间的比例关系得到比例系数,对模拟温度场进行积分后乘以比例系数得到温度应变。
或者,利用模拟温度场和温度应变公式计算温度应变,所述温度应变公式为:
其中,ε为待计算截面的温度应变,A表示待计算截面的截面面积,T(x,y)表示待计算截面上坐标为(x,y)处的模拟温度场,E代表弹性模量,αT表示热膨胀系数,H为待计算截面所处的高度,k1为待计算截面因周边构件约束而产生的轴向弹簧刚度。
实施例1
如图2所示,以335米高的武汉长江航运中心(所处纬度约为北纬31度)为对象,在ANSYS平台上根据武汉长江航运中心的结构施工图建立三维实体有限元模型,单位类型选择为Solid90,该单元为20节点六面体单元,每个节点有温度一个自由度。如图3(a)-3(c)所示,采用多尺度的网格划分方式,柱、核心筒、梁和板的单元尺寸可分别取为0.12m、0.20m、0.48m和0.60m。材料密度取为2635kg/m3,弹性模量取为3.07e10N/m2,热膨胀系数取为10e-6/℃。
本发明实施例1选择2017年11月18号至2017年12月16号为分析时段,计算该时段内的温度场。2017年11月18号00:00的太阳位置确定如下:计算该时刻的辐射历计时角为87120°,将维度和辐射历计时角τ=87120°带入下式中计算得太阳高度角β=-79.1°,α=120.0°。
当β<0时可认为太阳位于地平面以下。利用上述过程依次计算其他时刻的太阳相对位置。
确定太阳相对结构的相对位置后,在此位置上构建一个正六面体虚拟太阳,可取正六面体的长、宽、高分别0.3m×0.3m×0.03m,太阳距结构中心的距离为500m。
将虚拟太阳划分为一个单位,并与结构单元构成辐射对,利用ANSYS的热辐射模块计算每个辐射对的辐射系数,据此判断阴、阳面。
如图4(a)-4(c)所示,在气象站查询2017年11月18号至2017年12月16号之间的大气温度、太阳辐射强度和风速等气象参数,在阴面和阳面单元上施加第三类热边界条件(阴面单元上的太阳辐射强度取值为0)。
定义分析类型为瞬态热分析,定义初始时刻2017年11月18号00:00为均匀温度场,然后每小时施加一次热边界条件,求解可得研究时段内任一时刻的温度场。
随后进行温度应变的计算,例如要计算结构第18层某柱内一测点的温度应变,可由下式计算:
上式中的T(x,y)为有限元模拟所得的柱截面温度场。由上式知温度应变与面平均温度TA成正比,即与∫T(x,y)dA成正比。若取2017年11月18号0:00时刻的应变为0,则相对于该时刻的温度应变可通过计算∫T(x,y)dA得到,如图5所示虚线为计算所得的温度应变,实线为实测值,与实测值对比发现两者之间的相对均方根误差小于7με。本发明提供的计算方法可以用于剔除实测应变中的温度应变部分,从而得到结构构件的真实应变变化,这对评估结构的安全状态具有重要意义。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种高效计算超高层结构温度应变的方法,其特征在于,包括以下步骤:
(1)根据太阳相对位置建立虚拟太阳,根据超高层结构建筑的几何形状建立有限元模型;
(2)对有限元模型进行网格划分得到多个结构单元,将虚拟太阳设定为辐射源,结构单元设定为吸收体,计算辐射源和每个吸收体之间的辐射系数,利用辐射系数判断结构单元为阴面单元或阳面单元;
(3)在阴面单元上施加第三类热边界条件时,太阳辐射强度取0,在阴面单元和阳面单元上施加第三类热边界条件后,对有限元模型进行仿真模拟,得到模拟温度场,基于温度应变与模拟温度场成正比的原则计算温度应变。
3.如权利要求2所述的一种高效计算超高层结构温度应变的方法,其特征在于,所述虚拟太阳的具体建立方式为:
根据太阳高度角β和太阳方位角α计算太阳在任意时刻的位置,建立虚拟太阳。
4.如权利要求1-3任一所述的一种高效计算超高层结构温度应变的方法,其特征在于,所述步骤(2)中网格划分时,对超高层结构建筑中的柱、核心筒、梁和板采用不同的网格尺寸进行划分。
5.如权利要求1-3任一所述的一种高效计算超高层结构温度应变的方法,其特征在于,所述步骤(2)中利用辐射系数判断的具体实现方式为:
若辐射系数为0,则结构单元被判定为阴面单元;若辐射系数不等于0,则结构单元被判定为阳面单元。
7.如权利要求1-3任一所述的一种高效计算超高层结构温度应变的方法,其特征在于,所述步骤(3)还包括:
根据实测应变与模拟温度场之间的比例关系得到比例系数,对模拟温度场进行积分后乘以比例系数得到温度应变。
9.如权利要求1-3任一所述的一种高效计算超高层结构温度应变的方法,其特征在于,所述方法计算得到的温度应变应用于评估超高层结构建筑的安全状态,具体地:
从超高层结构建筑的实测应变中剔除温度应变,得到超高层结构建筑的真实应变,利用真实应变评估超高层结构建筑的安全状态。
10.一种高效计算超高层结构温度应变的***,其特征在于,包括以下模块:
模型建立模块,用于根据太阳相对位置建立虚拟太阳,根据超高层结构建筑的几何形状建立有限元模型;
阴面判断模块,用于对有限元模型进行网格划分得到多个结构单元,将虚拟太阳设定为辐射源,结构单元设定为吸收体,计算辐射源和每个吸收体之间的辐射系数,利用辐射系数判断结构单元为阴面单元或阳面单元;
模拟计算模块,用于在阴面单元上施加第三类热边界条件时,太阳辐射强度取0,在阴面单元和阳面单元上施加第三类热边界条件后,对有限元模型进行仿真模拟,得到模拟温度场,基于温度应变与模拟温度场成正比的原则计算温度应变。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010168266.2A CN111460545B (zh) | 2020-03-11 | 2020-03-11 | 一种高效计算超高层结构温度应变的方法和*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010168266.2A CN111460545B (zh) | 2020-03-11 | 2020-03-11 | 一种高效计算超高层结构温度应变的方法和*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111460545A true CN111460545A (zh) | 2020-07-28 |
CN111460545B CN111460545B (zh) | 2022-05-31 |
Family
ID=71684233
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010168266.2A Active CN111460545B (zh) | 2020-03-11 | 2020-03-11 | 一种高效计算超高层结构温度应变的方法和*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111460545B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112380595A (zh) * | 2020-10-27 | 2021-02-19 | 华中科技大学 | 超高层结构的温度变形预测模型建立方法及预测方法 |
CN113094795A (zh) * | 2021-04-21 | 2021-07-09 | 中南大学 | 一种桥梁温度测量方法、装置、设备及可读存储介质 |
CN118243267A (zh) * | 2024-05-27 | 2024-06-25 | 中建三局城建有限公司 | 一种用于超高层结构应力分布的光纤监测方法及*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5645773A (en) * | 1994-11-30 | 1997-07-08 | E.R.C. Co., Ltd. | Method for placing concrete for construction of a master concrete structure |
CN101251865A (zh) * | 2008-03-27 | 2008-08-27 | 上海交通大学 | 基于整体的大型钢结构防火的***化设计方法 |
CN102607644A (zh) * | 2012-02-23 | 2012-07-25 | 中建钢构有限公司 | 一种建筑施工监测***及其监测方法 |
CN106092628A (zh) * | 2016-06-06 | 2016-11-09 | 武汉理工大学 | 一种土木工程结构太阳辐射温度效应分析方法及*** |
-
2020
- 2020-03-11 CN CN202010168266.2A patent/CN111460545B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5645773A (en) * | 1994-11-30 | 1997-07-08 | E.R.C. Co., Ltd. | Method for placing concrete for construction of a master concrete structure |
CN101251865A (zh) * | 2008-03-27 | 2008-08-27 | 上海交通大学 | 基于整体的大型钢结构防火的***化设计方法 |
CN102607644A (zh) * | 2012-02-23 | 2012-07-25 | 中建钢构有限公司 | 一种建筑施工监测***及其监测方法 |
CN106092628A (zh) * | 2016-06-06 | 2016-11-09 | 武汉理工大学 | 一种土木工程结构太阳辐射温度效应分析方法及*** |
Non-Patent Citations (7)
Title |
---|
FEIGAO 等: ""Efficient calculation and monitoring of temperature actions on supertall"", 《ENGINEERING STRUCTURES》 * |
JIA ZHANSU: ""Field monitoring and numerical simulation of the thermal actions of a supertall structure"", 《STRUCTURE CONTROL AND HEALTH MONITORING》 * |
康为江: ""钢筋混凝土箱梁日照温度效应研究"", 《中国优秀硕士学位论文全文期刊库工程科技Ⅱ辑》 * |
谢鹏 等: ""基于ANSYS软件的混凝土薄壁箱梁日照温度场有限元分析"", 《公路工程》 * |
陈潘: ""日照非均匀温度场对施工期内超高层结构的影响分析"", 《日照非均匀温度场对施工期内超高层结构的影响分析》 * |
顾斌 等: ""大跨桥梁结构三维日照温度场计算方法"", 《东南大学学报》 * |
高飞 等: ""非均匀日照条件下结构的三维温度场分析"", 《土木工程与管理学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112380595A (zh) * | 2020-10-27 | 2021-02-19 | 华中科技大学 | 超高层结构的温度变形预测模型建立方法及预测方法 |
CN112380595B (zh) * | 2020-10-27 | 2024-04-19 | 华中科技大学 | 超高层结构的温度变形预测模型建立方法及预测方法 |
CN113094795A (zh) * | 2021-04-21 | 2021-07-09 | 中南大学 | 一种桥梁温度测量方法、装置、设备及可读存储介质 |
CN118243267A (zh) * | 2024-05-27 | 2024-06-25 | 中建三局城建有限公司 | 一种用于超高层结构应力分布的光纤监测方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN111460545B (zh) | 2022-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111460545B (zh) | 一种高效计算超高层结构温度应变的方法和*** | |
CN109992900B (zh) | 一种大体积混凝土多场实时在线协同智能仿真方法和*** | |
JP4404220B2 (ja) | 気体状況予測装置、方法、プログラム、および拡散状況予測システム | |
Tome et al. | Structural response of a concrete cable-stayed bridge under thermal loads | |
CN111723509B (zh) | 一种桥梁结构温度场监测方法 | |
CN114357571B (zh) | 建成建筑环境下大气边界层风场特性的反演方法和*** | |
CN105956216A (zh) | 大跨钢桥基于均匀温度响应监测值的有限元模型修正方法 | |
CN105930571A (zh) | 基于单位温度响应监测值的大跨钢桥有限元模型修正方法 | |
Taylor et al. | Hybrid pressure integration and buffeting analysis for multi-row wind loading in an array of single-axis trackers | |
CN112580241B (zh) | 一种基于结构降阶模型的非线性气动弹性动响应分析方法 | |
Larsson et al. | Modelling of climatic thermal actions in hollow concrete box cross-sections | |
Gao et al. | Efficient calculation and monitoring of temperature actions on supertall structures | |
CN105824987A (zh) | 一种基于遗传算法的风场特征统计分布模型建立方法 | |
Aly et al. | Wind design of solar panels for resilient and green communities: CFD with machine learning | |
CN115470731A (zh) | 一种基于微气象的风场风电功率预测方法及*** | |
Li et al. | Global temperature behavior monitoring and analysis of a three-tower cable-stayed bridge | |
CN104091003A (zh) | 一种基础运动时柔性壳结构大变形响应的有限元建模方法 | |
Aßmus et al. | Consideration of non-uniform and non-orthogonal mechanical loads for structural analysis of photovoltaic composite structures | |
Reda et al. | Large eddy simulation of wind flow through an urban environment in its full-scale wind tunnel models | |
Abid et al. | Three-dimensional thermal modeling of temperature variation in concrete box-girders using COMSOL Multiphysics® | |
Leitão et al. | Heat transfer analysis of infrastructures subjected to environmental actions: A finite element solver PAT | |
CN109252441B (zh) | 一种变截面箱梁剪力滞效应的分析方法 | |
Prilliman et al. | Technoeconomic analysis of changing PV array convective cooling through changing array spacing | |
CN115859715A (zh) | 一种周期性多孔结构材料等效性能的计算方法 | |
CN113176054B (zh) | 桥梁钢管拱肋变形位移监控*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |