CN111598448B - 一种基于混凝土t型梁火灾后损伤评估方法 - Google Patents
一种基于混凝土t型梁火灾后损伤评估方法 Download PDFInfo
- Publication number
- CN111598448B CN111598448B CN202010413146.4A CN202010413146A CN111598448B CN 111598448 B CN111598448 B CN 111598448B CN 202010413146 A CN202010413146 A CN 202010413146A CN 111598448 B CN111598448 B CN 111598448B
- Authority
- CN
- China
- Prior art keywords
- matrix
- formula
- ideal solution
- weight
- evaluation
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 121
- 230000006378 damage Effects 0.000 title claims abstract description 48
- 239000004567 concrete Substances 0.000 title claims abstract description 17
- 238000011156 evaluation Methods 0.000 claims abstract description 125
- 230000008569 process Effects 0.000 claims abstract description 28
- 238000005457 optimization Methods 0.000 claims abstract description 16
- 239000006185 dispersion Substances 0.000 claims abstract description 12
- 239000011159 matrix material Substances 0.000 claims description 126
- 239000013598 vector Substances 0.000 claims description 35
- 238000010606 normalization Methods 0.000 claims description 19
- 238000013519 translation Methods 0.000 claims description 11
- 238000012545 processing Methods 0.000 claims description 9
- 230000009467 reduction Effects 0.000 claims description 8
- 230000008859 change Effects 0.000 claims description 5
- 230000009286 beneficial effect Effects 0.000 claims description 4
- 241000764238 Isis Species 0.000 claims description 3
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000000903 blocking effect Effects 0.000 claims 1
- 230000006872 improvement Effects 0.000 abstract description 4
- 238000012937 correction Methods 0.000 abstract description 3
- 230000007547 defect Effects 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 14
- 238000005259 measurement Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000008439 repair process Effects 0.000 description 3
- 230000002787 reinforcement Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- XLVOWGDFXAVSFE-UHFFFAOYSA-N (4-amino-2-methylsulfanylpyrimidin-5-yl)methanol Chemical compound CSC1=NC=C(CO)C(N)=N1 XLVOWGDFXAVSFE-UHFFFAOYSA-N 0.000 description 1
- 230000009172 bursting Effects 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000002950 deficient Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011241 protective layer Substances 0.000 description 1
- 239000011150 reinforced concrete Substances 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/08—Construction
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Educational Administration (AREA)
- Development Economics (AREA)
- Data Mining & Analysis (AREA)
- Entrepreneurship & Innovation (AREA)
- Mathematical Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Algebra (AREA)
- Health & Medical Sciences (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Quality & Reliability (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
Abstract
一种基于混凝土T型梁火灾后损伤评估方法,涉及火灾后建筑物损伤等级评估方法技术领域,所述的评估方法步骤如下:步骤1、确定参数指标;步骤2、根据所确定的参数指标,利用多种评估方法得到各评价对象的综合贴近度;所述的多种评估方法包括改进层次分析法、熵权法、逼近理想解法、离差平方和最优化法、灰色关联度法、正交投影法;步骤3、根据所述的各评价对象的综合相对贴近度,最终确定建筑物的损伤等级。本发明建立了一套基于最优组合权重的灰色关联度修正改进逼近理想解(MTOPSIS‑GRA)法的综合评估体系,使评估结果更加客观,避免了现有技术的评估方法的技术缺陷,适合于工程应用推广。
Description
技术领域
本发明涉及火灾后建筑物损伤等级评估方法技术领域,具体涉及一种基于混凝土T型梁火灾后损伤评估方法。
背景技术
近年来,建筑物火灾问题比较常见,为减少损失,防止火灾的发生,有关专家作出了相关的研究,但对于火灾后建筑物的损伤等级评定及灾后的修复工作研究较少,且具有参考价值的相关书籍、文献比较匮乏。
《火灾后建筑结构鉴定标准》(CECS252-2009)是现在使用最广、认可度最高的国家建议性标准。然而该标准基本考虑的都是灾后建筑物的外观现象,评价标准主观性较强,通过外观确定损伤等级模糊性较大,往往取决于人的主观意识,因此,评估结果的误差也会偏大。
《火灾后损伤评价建议准则》是通过对灾后承载能力的折减程度确定损伤等级,此时确定结果较准确,可信度较高,可是由于灾后承载力的获取难免会对构件有一定的损伤甚至是破坏,且仅靠一个指标来确定会有一定的偶然性。
也有人通过火灾后的承载力、基频、刚度和受火时间等客观参数给出了一套简单的各因素折减评价体系,但没有具体的求解方法,对于指标的确定具有一定的复杂度。
另外,现有的数学评估方法中,如层次分析法存在着主观性太强的缺陷,而逼近理想解法(TOPSIS)则由于权重信息是事先给定,结果有一定的主观性,且在应用中由于新增加方案而容易产生逆序问题。若直接将上述评估方法应用到火灾后建筑物损伤评估中,将难以得到良好的评估结果。
发明内容
为解决现有技术中存在的问题,本发明提供了一种基于混凝土T型梁火灾后损伤评估方法,该方法基于最优组合权重的灰色关联度修正改进逼近理想解(MTOPSIS-GRA)法,对钢筋混凝土构件灾后损伤进行综合评价,并对评价结果进行量化,是一种针对混凝土T形梁火灾后损伤评估方法。
为实现上述目的,本发明的技术方案是:
一种基于混凝土T型梁火灾后损伤评估方法,所述的评估方法步骤如下:
步骤1、确定参数指标;
所述的参数指标依据火灾后建筑物的如下评价指标确定:表观现象X1、火灾后建筑物的表面最高温度X2、承载力折减X3、基频折减X4、刚度折减X5、受火时间X6,所述的表观现象X1包括但不限于爆裂因素、裂缝深度和面积、及裂缝宽度;
步骤2、根据所确定的参数指标,利用多种评估方法得到各评价对象的综合相对贴近度;所述的多种评估方法包括改进层次分析法、熵权法、逼近理想解法、离差平方和最优化法、灰色关联度法、正交投影法;
步骤3、根据所述的各评价对象的综合相对贴近度,最终确定建筑物的损伤等级。
优选的,所述的步骤2中,确定各评价对象的综合相对贴近度包括如下具体步骤:
(1)、通过所述的改进层次分析法确定主观权重,通过熵权法确定客观权重,用所述的离差平方和最优化法对主观权重和客观权重进行组合赋权,得到组合权重;
(2)、通过所述的逼近理想解法结合组合权重获得正、负理想解,然后通过正交投影法计算评价对象与正负理想解的垂直距离,通过正交投影法对逼近理想解法的改进,从空间的垂直距离出发,获得准确度更高的正、负理想解,最终结果用Pi表示;通过灰色关联度法对逼近理想解法进行修正,并结合所述的垂直距离确定各评价对象的综合相对贴近度。
优选的,所述的改进层次分析法采用专家评分的方式进行操作,其获得主观权重的过程具体步骤包括:依据改进层次分析法获得决策矩阵的步骤;依据决策矩阵获得主观权重的步骤;
所述的依据改进层次分析法获得决策矩阵的步骤包括:
(1)、根据三标度法对各参数指标进行打分确立比较矩阵;参数指标两两相互对比确定比较矩阵,A方阵为:
i,j为所选取的指标数量,两两比较确定矩阵A;
(2)、将比较矩阵转化为判断矩阵:
所述的转化运用极差法,公式为:
式中:R=rmax-rmin,称为极差;rmax=max{r1,r2,...,rn};rmin=min{r1,r2,...,rn};cb是指指标元素对的相对重要度;
所述的依据决策矩阵获得主观权重的步骤所采用的公式为:
优选的,所述的熵权法确定客观权重的步骤包括:依据熵权法获得决策矩阵的步骤;依据决策矩阵获得客观权重的步骤;
所述的依据熵权法获得决策矩阵的步骤包括:
(1)、根据评估对象和参数指标确定原始矩阵,公式为:
(2)、矩阵归一化处理,归一化处理的原则为,根据指标所对应优劣的不同划分,分为经济型指标和成本型指标,成本型指标数值越小越优,经济型指标反之,值越大越优,具体公式如下:
经济型指标归一化公式:
成本型指标归一化公式:
式中:max{r′ij}—矩阵R'中第i行的最大值;min{r′ij}—矩阵R'中第i行的最小值;
归一化矩阵表示为:
所述的依据决策矩阵获得客观权重的步骤包括计算信息熵及熵权重的步骤,公式分别为:
第i个指标的信息熵计算公式为:
第i个指标的熵权重计算公式为:
优选的,所述的用离差平方和最优化法对主观权重和客观权重进行组合赋权,得到组合权重的步骤的具体步骤如下:
(1)、确定第i个评价对象与其他对象的组合评价值的离差平方和,公式为;
(2)、构造目标函数,使M个评价对象总的离差平方和达到最大,公式为:
(3)、将目标函数转化为矩阵函数,以利于下一步参数的获取和使用:
令矩阵A1为:
目标函数H(Wc)可表示为:
(4)、得出最优化公式,求解方式为:
式中,W为主客观权系数向量组成的分块矩阵,是线性表示出的瑞利商,且WTA1W为对称矩阵,根据瑞利商的理论可知,存在的最大值和最小值就是对称矩阵WTA1W所对应的最大特征值和最小特征值,则WTA1W矩阵的最大特征根即为的最大值,所对应的单位化特征向量就是所求的最优解,且根据矩阵WTA1W为对称非负的,根据Perron-Frobenius定理,可知特征向量全部为正实数;
(5)、确定最优组合权重:
对组合权重向量作归一化处理:
(1)、建立归一化矩阵,对原始矩阵无量纲化,公式为:
(2)、计算加权矩阵,将组合权重与无量纲化矩阵相乘,公式为:
rij=wj×vij,(i=1,2,…,m;j=1,2,…,n)
(3)、简化后计算正、负理想解;
无量纲化后,对于经济型指标,每个指标中的最大值就是评价对象的最优指标;对于成本型指标,对应最小值则为所有评价对象中的最优指标,
正理想解可表示为:
负理想解可表示为:
优选的,所述的通过正交投影法计算评价对象与正负理想解的垂直距离的具体步骤为:
经过平移后的理想解为0,则距离可表示为:
式中:Pi为评价对象与理想解的垂直距离;Ti为平移后矩阵的行向量;距离越近结果越优,可知Pi越小,评价对象越优。
优选的,所述的灰色关联度法对逼近理想解法进行修正,并结合所述的垂直距离确定各评价对象的综合相对贴近度的步骤包括:
(1)、确定最优指标集及评判矩阵;
(2)、对评判矩阵规范化处理;
(3)、计算加权标准化矩阵;
(4)、以正负理想解为参考数列计算灰色关联度;
(5)、确定各评价对象的综合相对贴近度;
所述的以正负理想解为参考数列计算灰色关联度的具体步骤为:
(1)、确定关联度矩阵,并获得灰色关联度:
将原始灰色评价体系确定关联度公式中的最优解更换为逼近理想解法得到的正负理想解,得到相应关联度系数为:
式中:ξ取0.5,代表分辨系数;
则正理想解的灰色关联度矩阵为:
从而得到正理想解的灰色关联度为:
同理可得由负理想解确定的关联度系数为:
则负理想解的灰色关联度矩阵为:
从而得到负理想解的灰色关联度为:
(2)、获得综合关联度,并确定最优评价对象:
将各个参数指标相对于理想解的垂直距离Pi与关联度合并;对结果进行无量纲化处理,确定新的关联度,公式为:
构造综合关联度,公式为:
式中:综合关联度Fi +表示与理想方案的贴近程度,数值越大表示评价对象越优;综合关联度Fi -表示远离程度,数值越大则表示评价对象越劣;β表示偏好系数,表示位置和趋势相近性的偏好程度,取值范围是(0,1);
所述确定各评价对象的综合相对贴近度的公式为:
式中:F1 *为综合相对贴近度,其中,综合相对贴近度越大,所得评价对象评价等级越优越,反之则评价等级越恶劣。所述的损伤等级根据评价对象的评价等级优劣来确定。
本发明一种基于混凝土T型梁火灾后损伤评估方法具有如下有益效果:
1、本发明为解决火灾后混凝土梁的损伤提供较为客观、可行的评价方法,为损伤后的加固与修复提供更精确的技术参考资料,可以节约加固修复成本,并保证火灾后混凝土梁的承载力。
2、本发明基于逼近理想法(TOPSIS法)的不足,通过正交投影法对其进行改进,从空间的垂直距离出发,获得准确度更高的正负理想解。
3、本发明建立了一套基于最优组合权重的灰色关联度修正改进逼近理想解(MTOPSIS-GRA)法的综合评估体系,使评估结果更加客观。
4、本发明通过采用改进层次分析法(MAHP)确定主观权重,熵权法确定客观权重,用离差平方和最优化法对主客观进行组合赋权,得到最后的最优组合权重,大大减少了人主观性的影响。
5、本发明通过灰色关联度法(GRA)对改进逼近理想解法(MTOPSIS)进行修正,用灰色关联度法结合改进逼近理想法(MTOPSIS-GRA)确定最优评价对象,保证指标变化趋势的相近性,将GRA融入到MTOPSIS中,可以解决MTOPSIS不能反映指标变化趋势相近性的问题,使结果可靠度更高。
附图说明
图1、为本发明的流程框图;
图2、规范中的火灾后损伤评价指标三阶框图;
图3、为本发明火灾后评价指标二阶框图;
具体实施方式
以下所述,是以阶梯递进的方式对本发明的实施方式详细说明,该说明仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
本发明的描述中,需要说明的是,术语“上”“下”“左”“右”“顶”“底”“内”“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以及特定的方位构造和操作,因此不能理解为对本发明的限制。
如图1所示,一种基于混凝土T型梁火灾后损伤评估方法,所述的评估方法步骤如下:
步骤1、确定参数指标;
如图3所示,所述的参数指标依据火灾后建筑物的如下评价指标确定:表观现象X1、火灾后建筑物的表面最高温度X2、承载力折减X3、基频折减X4、刚度折减X5、受火时间X6,所述的表观现象X1包括但不限于爆裂因素、裂缝深度和面积、及裂缝宽度;
如图2所示,通常火灾后损伤评价指标分类方式很多,然而细分到第三阶时变得很复杂,从而影响到火灾后损伤评价的实用性,因此,考虑到评价方式的可实现性,本发明将评价指标调整为二阶体系,如图3所示;
由于火灾后确定表观参数等因素主观性较大,本发明将图2中的指标定量转化主观性较高的灾后梁体表观参数及灾后梁体损伤特征指标综合考虑,由于爆裂因素统计复杂,且不规律性太强,裂缝深度和面积测量较为复杂,裂缝宽度等指标也受到火灾后发生爆裂的环境因素影响,因此最终确定为一个综合指标:表观现象,另外,受火时间、受火最高温度、受火面积三个指标间相对关联度较高,且受火面积范围划分较复杂,选取表面最高温度和受火时间为评价指标;灾后梁体力学参数是对残余性能较高的评价,指标选取综合考虑到每一个指标的重要性和相互影响作用,且理论计算及实际损伤都较易获取,因此选取基频折减、刚度折减和承载力折减作为评估指标;
步骤2、根据所确定的参数指标,利用多种评估方法得到各评价对象的综合相对贴近度;所述的多种评估方法包括改进层次分析法、熵权法、逼近理想解法、离差平方和最优化法、灰色关联度法、正交投影法;
步骤3、根据所述的各评价对象的综合相对贴近度,最终确定建筑物的损伤等级;
所述的步骤2中,确定各评价对象的综合相对贴近度包括如下具体步骤:
(1)、通过所述的改进层次分析法确定主观权重,通过熵权法确定客观权重,用所述的离差平方和最优化法对主观权重和客观权重进行组合赋权,得到组合权重;离差平方和最优化是计算组合权重中较为合理的方法,它取的是指标间的差的平方和最大化来确定综合权重,此方法对各种赋权方法考虑较为全面,结果合理,能够明显得出不同对象的相对大小,且准确性高;
(2)、通过所述的逼近理想解法结合组合权重获得正、负理想解,然后通过正交投影法计算评价对象与正负理想解的垂直距离;通过灰色关联度法对逼近理想解法进行修正,并结合所述的垂直距离确定各评价对象的综合相对贴近度;现有技术中的逼近理想解法由于权重信息事先给定,导致结果具有一定的主观性,另外逼近理想法在应用中由于新增加的方案容易产生逆序问题,故本发明通过正交投影法对其进行改进,通过从空间的垂直距离出发,可以获得精度更高的正负理想解;本发明通过灰色关联度法对经过改进的逼近理想解法进行进一步修正,可以确定最优评价对象,保证指标变化趋势的相近性,将灰色关联度法融入到改进的逼近理想解法中,可以解决改进的逼近理想解法不能反映指标变化趋势相近性的问题,使评价结果的可靠性更高;
所述的改进层次分析法采用专家评分的方式进行操作,其获得主观权重的过程具体步骤包括:依据改进层次分析法获得决策矩阵的步骤;依据决策矩阵获得主观权重的步骤,本发明所述的决策矩阵包括了比较矩阵和判断矩阵,具体论述如下:
所述的依据改进层次分析法获得决策矩阵的步骤包括:
(1)、根据三标度法对各参数指标进行打分确立比较矩阵;参数指标两两相互对比确定比较矩阵,A方阵为:
i,j为所选取的指标数量,两两比较确定矩阵A;
(2)、将比较矩阵转化为判断矩阵:
所述的转化运用极差法,公式为:
式中:R=rmax-rmin,称为极差;rmax=max{r1,r2,...,rn};rmin=min{r1,r2,...,rn};cb是指指标元素对的相对重要度;
所述的依据决策矩阵获得主观权重的步骤所采用的公式为:
所述的熵权法确定客观权重的步骤包括:依据熵权法获得决策矩阵的步骤;依据决策矩阵获得客观权重的步骤;
所述的依据熵权法获得决策矩阵的步骤包括:
(1)、根据评估对象和参数指标确定原始矩阵,公式为:
(2)、矩阵归一化处理,归一化处理的原则为,根据指标所对应优劣的不同划分,分为经济型指标和成本型指标,成本型指标数值越小越优,经济型指标反之,值越大越优,具体公式如下:
经济型指标归一化公式:
成本型指标归一化公式:
式中:max{r′ij}—矩阵R'中第i行的最大值;min{r′ij}—矩阵R'中第i行的最小值;
归一化矩阵表示为:
所述的依据决策矩阵获得客观权重的步骤包括计算信息熵及熵权重的步骤,公式分别为:
第i个指标的信息熵计算公式为:
第i个指标的熵权重计算公式为:
所述的用离差平方和最优化法对主观权重和客观权重进行组合赋权,得到组合权重的步骤的具体步骤如下:
(1)、确定第i个评价对象与其他对象的组合评价值的离差平方和,公式为;
(2)、构造目标函数,使M个评价对象总的离差平方和达到最大,公式为:
(3)、将目标函数转化为矩阵函数,以利于下一步参数的获取和使用:
令矩阵A1为:
目标函数H(Wc)可表示为:
(4)、得出最优化公式,求解方式为:
式中,W为主客观权系数向量组成的分块矩阵,是线性表示出的瑞利商,且WTA1W为对称矩阵,根据瑞利商的理论可知,存在的最大值和最小值就是对称矩阵WTA1W所对应的最大特征值和最小特征值,则WTA1W矩阵的最大特征根即为的最大值,所对应的单位化特征向量就是所求的最优解,且根据矩阵WTA1W为对称非负的,根据Perron-Frobenius定理,可知特征向量全部为正实数;
(5)、确定最优组合权重:
对组合权重向量作归一化处理:
(1)、建立归一化矩阵,对原始矩阵无量纲化,公式为:
(2)、计算加权矩阵,将组合权重与无量纲化矩阵相乘,公式为:
rij=wj×vij,(i=1,2,…,m;j=1,2,…,n)
(3)、简化后计算正、负理想解;
无量纲化后,对于经济型指标,每个指标中的最大值就是评价对象的最优指标;对于成本型指标,对应最小值则为所有评价对象中的最优指标,
正理想解可表示为:
负理想解可表示为:
优选的,所述的通过正交投影法计算评价对象与正负理想解的垂直距离的具体步骤为:
经过平移后的理想解为0,则距离可表示为:
式中:Pi为评价对象与理想解的垂直距离;Ti为平移后矩阵的行向量;距离越近结果越优,可知Pi越小,评价对象越优。
所述的灰色关联度法对逼近理想解法进行修正,并结合所述的垂直距离确定各评价对象的综合相对贴近度的步骤包括:
(1)、确定最优指标集及评判矩阵;
(2)、对评判矩阵规范化处理;
(3)、计算加权标准化矩阵;
(4)、以正负理想解为参考数列计算灰色关联度;
(5)、确定各评价对象的综合相对贴近度;
所述的以正负理想解为参考数列计算灰色关联度的具体步骤为:
(1)、确定关联度矩阵,并得到灰色关联度:
将原始灰色评价体系确定关联度公式中的最优解更换为逼近理想解法得到的正负理想解,得到相应关联度系数为:
式中:ξ取0.5,代表分辨系数;
则正理想解的灰色关联度矩阵为:
从而得到正理想解的灰色关联度为:
同理可得由负理想解确定的关联度系数为:
则负理想解的灰色关联度矩阵为:
从而得到负理想解的灰色关联度为:
(2)、获得综合关联度,并确定最优评价对象:
将参数指标相对于理想解的垂直距离Pi与关联度合并;对结果进行无量纲化处理,确定新的关联度,公式为:
构造综合关联度,公式为:
式中:综合关联度Fi +表示与理想方案的贴近程度,数值越大表示评价对象越优;综合关联度Fi -表示远离程度,数值越大则表示评价对象越劣;β表示偏好系数,表示位置和趋势相近性的偏好程度,取值范围是(0,1);
所述确定各评价对象的综合相对贴近度的公式为:
式中:F1 *为综合相对贴近度,其中,综合相对贴近度越大,所得评价对象评价等级越优越,反之则评价等级越恶劣。
实验例:
设计10根T形梁长度均为3000mm,净跨为2800mm,截面及配筋情况完全相同,保护层厚度为30mm,1根对比梁只做静载试验,另外9根T形梁进行不同工况的横载升温试验,试验分组见下表:
试验分组情况
注:B-T06L0解释:B代表梁,T06表示受火60min,L0表示荷载比为0;N代表自然状态;
表示开裂荷载,使梁在高温前处于未开裂状态;常遇荷载为恒载+0.4活荷载,0.4为常遇系数,此代表正常使用状态下的荷载。
B-L0、B-L12试验梁损伤评定
主观指标权重运用改进层次分析法,此项权重各指标的确定是根据专家打分来确定的,跟具体试验所得结果没有关系,因此主观权重为:
W1=[0.0684,0.1282,0.4498,0.1754,0.1282,0.0500]T
火灾后各指标损伤程度根据量测值及熵权法计算如下表所示
B-L0、B-L12试验梁灾后各指标量测结果
根据离差平方和最优化公式,得到非负定矩阵:
由主客观权重组合,构造分块矩阵:
从而得到最大特征根λmax=20.3594,
对主客观权重进行线性组合得到最优化组合权重:
归一化处理后得到各指标的组合权重为:
确定指标的多属性综合评价值后,根据MTOPSIS-GRA求解每个评价对象的损伤等级,得到修正后的贴近度值为:
F* 1=0.6655,F* 2=0.5694,F* 3=0.4071,F* 4=0.3727,F* 5=0.2324,F* 6=0.2169
可知损伤程度从小到大排序为:F* 1<F* 2<F* 3<F* 4<F* 5<F* 6
根据所得结果,确定火灾后B-L0、B-L12试验梁损伤等级及评分见下表
B-L0、B-L12试验梁灾后损伤结果
B-L44试验梁的损伤评定
B-L44试验梁灾后各指标量测结果
主观权重:W1=[0.0684,0.1282,0.4498,0.1754,0.1282,0.0500]T
客观权重:W2=[0.1755,0.1453,0.1502,0.2206,0.1510,0.1574]T
归一化后最优组合权重:Wc *=(0.1220,0.1368,0.2999,0.1980,0.1396,0.1037)T
通过计算,得到B-L44三根试验梁相对贴近度分别为:
F* 1=0.5046,F* 2=0.3494,F* 3=0.2064
对应损伤等级见下表
B-L44试验梁灾后损伤结果
由此可见,运用本发明综合评价体系求解得到的损伤等级与初步评级的损伤等级基本一致,说明了本发明所使用的评价体系的合理性与适应性,B-T09L12试验梁评级有所不同,但通过评分可以看出,贴近度值在等级的界限处,因此导致结果有所差异。从损伤等级可以看出,在受火前没有裂缝的情况下,受火时间越长损伤程度越严重,综合考虑多指标对评价结果的可信度高,考虑的客观指标颇多,对最终结果所占比重大,使结果更准确。
Claims (4)
1.一种基于混凝土T型梁火灾后损伤评估方法,其特征为:所述的评估方法步骤如下:
步骤1、确定参数指标;
所述的参数指标由火灾后建筑物的如下评价指标构成:表观现象X1、火灾后建筑物的表面最高温度X2、承载力折减X3、基频折减X4、刚度折减X5、受火时间X6,所述的表观现象X1包括爆裂因素、裂缝深度和面积、及裂缝宽度;
步骤2、根据所确定的参数指标,利用多种评估方法得到各评价对象的综合相对贴近度;所述的多种评估方法包括改进层次分析法、熵权法、逼近理想解法、离差平方和最优化法、灰色关联度法、正交投影法;
步骤3、根据所述的各评价对象的综合相对贴近度,最终确定建筑物的损伤等级;
所述的步骤2中,确定各评价对象的综合相对贴近度包括如下具体步骤:
(1)、通过所述的改进层次分析法确定主观权重,通过熵权法确定客观权重,用所述的离差平方和最优化法对主观权重和客观权重进行组合赋权,得到组合权重;
(2)、通过所述的逼近理想解法结合组合权重获得正、负理想解,然后通过正交投影法计算评价对象与正负理想解的垂直距离,通过正交投影法对逼近理想解法的改进,从空间的垂直距离出发,获得准确度更高的正、负理想解,最终结果用表示;通过灰色关联度法对逼近理想解法进行修正,并结合所述的垂直距离确定各评价对象的综合相对贴近度;
所述的改进层次分析法采用专家评分的方式进行操作,其获得主观权重的过程具体步骤包括:依据改进层次分析法获得决策矩阵的步骤;依据决策矩阵获得主观权重的步骤;
所述的依据改进层次分析法获得决策矩阵的步骤包括:
(1)、根据三标度法对各参数指标进行打分确立比较矩阵;参数指标两两相互对比确定比较矩阵,A方阵为:
i,j为所选取的指标数量,两两比较确定矩阵A;
(2)、将比较矩阵转化为判断矩阵:
所述的转化运用极差法,公式为:
式中:R=rmax-rmin,称为极差;rmax=max{r1,r2,...,rn};rmin=min{r1,r2,...,rn};cb是指指标元素对的相对重要度;
所述的依据决策矩阵获得主观权重的步骤所采用的公式为:
所述的熵权法确定客观权重的步骤包括:依据熵权法获得决策矩阵的步骤;依据决策矩阵获得客观权重的步骤;
所述的依据熵权法获得决策矩阵的步骤包括:
(1)、根据评估对象和参数指标确定原始矩阵,公式为:
(2)、矩阵归一化处理,归一化处理的原则为,根据指标所对应优劣的不同划分,分为经济型指标和成本型指标,成本型指标数值越小越优,经济型指标反之,值越大越优,具体公式如下:
经济型指标归一化公式:
成本型指标归一化公式:
式中:max{r′ij}—矩阵R'中第i行的最大值;min{r′ij}—矩阵R'中第i行的最小值;
归一化矩阵表示为:
所述的依据决策矩阵获得客观权重的步骤包括计算信息熵及熵权重的步骤,公式分别为:
第i个指标的信息熵计算公式为:
第i个指标的熵权重计算公式为:
所述的用离差平方和最优化法对主观权重和客观权重进行组合赋权,得到组合权重的步骤的具体步骤如下:
(1)、确定第i个评价对象与其他对象的组合评价值的离差平方和,公式为;
(2)、构造目标函数,使M个评价对象总的离差平方和达到最大,公式为:
(3)、将目标函数转化为矩阵函数,以利于下一步参数的获取和使用:
令矩阵A1为:
目标函数H(Wc)可表示为:
(4)、得出最优化公式,求解方式为:
式中,W为主客观权系数向量组成的分块矩阵,是线性表示出的瑞利商,且WTA1W为对称矩阵,根据瑞利商的理论可知,存在的最大值和最小值就是对称矩阵WTA1W所对应的最大特征值和最小特征值,则WTA1W矩阵的最大特征根即为的最大值,所对应的单位化特征向量就是所求的最优解,且根据矩阵WTA1W为对称非负的,根据Perron-Frobenius定理,可知特征向量全部为正实数;
(5)、确定最优组合权重:
对组合权重向量作归一化处理,公式为:
2.如权利要求1所述的一种基于混凝土T型梁火灾后损伤评估方法,其特征为:通过所述的逼近理想解法结合组合权重获得正、负理想解的具体步骤为:
(1)、建立归一化矩阵,对原始矩阵无量纲化,公式为:
(2)、计算加权矩阵,将组合权重与无量纲化矩阵相乘,公式为:
rij=wj×vij,(i=1,2,…,m;j=1,2,…,n)
(3)、简化后计算正、负理想解;
无量纲化后,对于经济型指标,每个指标中的最大值就是评价对象的最优指标;对于成本型指标,对应最小值则为所有评价对象中的最优指标,
正理想解可表示为:
负理想解可表示为:
4.如权利要求3所述的一种基于混凝土T型梁火灾后损伤评估方法,其特征为:所述的灰色关联度法对逼近理想解法进行修正,并结合所述的垂直距离确定各评价对象的综合相对贴近度的步骤包括:
(1)、确定最优指标集及评判矩阵;
(2)、对评判矩阵规范化处理;
(3)、计算加权标准化矩阵;
(4)、以正负理想解为参考数列计算灰色关联度;
(5)、确定各评价对象的综合相对贴近度;
所述的以正负理想解为参考数列计算灰色关联度的具体步骤为:
(1)、确定关联度矩阵,并获得灰色关联度:
将原始灰色评价体系确定关联度公式中的最优解更换为逼近理想解法得到的正负理想解,得到相应关联度系数为:
式中:ξ取0.5,代表分辨系数;
则正理想解的灰色关联度矩阵为:
从而得到正理想解的灰色关联度为:
同理可得由负理想解确定的关联度系数为:
则负理想解的灰色关联度矩阵为:
从而得到负理想解的灰色关联度为:
(2)、获得综合关联度,并确定最优评价对象:
将各个参数指标相对于理想解的垂直距离Pi与关联度合并;对结果进行无量纲化处理,确定新的关联度,公式为:
构造综合关联度,公式为:
式中:综合关联度Fi +表示与理想方案的贴近程度,数值越大表示评价对象越优;综合关联度Fi -表示远离程度,数值越大则表示评价对象越劣;β表示偏好系数,表示位置和趋势相近性的偏好程度,取值范围是(0,1);
所述确定各评价对象的综合相对贴近度的公式为:
式中:Fi*为综合相对贴近度,其中,综合相对贴近度越大,所得评价对象评价等级越优越,反之则评价等级越恶劣。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010413146.4A CN111598448B (zh) | 2020-05-15 | 2020-05-15 | 一种基于混凝土t型梁火灾后损伤评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010413146.4A CN111598448B (zh) | 2020-05-15 | 2020-05-15 | 一种基于混凝土t型梁火灾后损伤评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111598448A CN111598448A (zh) | 2020-08-28 |
CN111598448B true CN111598448B (zh) | 2021-05-18 |
Family
ID=72192179
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010413146.4A Active CN111598448B (zh) | 2020-05-15 | 2020-05-15 | 一种基于混凝土t型梁火灾后损伤评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111598448B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112488463A (zh) * | 2020-11-13 | 2021-03-12 | 山东大学 | 一种基于组合赋权法与优劣解距离法的滑坡评价方法 |
CN112633631B (zh) * | 2020-11-26 | 2023-03-24 | 上海交通大学 | 多电源系统互补性的评价方法 |
CN112634078B (zh) * | 2020-12-18 | 2021-09-21 | 南京工程学院 | 一种基于多维指标融合的大工业负荷中断优先级评估方法 |
CN112668078B (zh) * | 2020-12-24 | 2022-05-17 | 青岛理工大学 | 一种锈蚀钢筋混凝土梁火灾后损伤识别方法 |
CN112699553B (zh) * | 2020-12-29 | 2021-11-02 | 昆明理工大学 | 一种岩爆烈度等级智能预测系统方法 |
CN112949115B (zh) * | 2021-02-01 | 2022-07-01 | 青岛理工大学 | 一种t形叠合梁火灾后残余承载力计算方法 |
CN113361923B (zh) * | 2021-03-29 | 2024-03-26 | 沈阳工业大学 | 一种基于spa和topsis耦合的工业园区负荷灵活度量化方法 |
CN113450333B (zh) * | 2021-06-30 | 2022-01-28 | 哈尔滨工业大学 | 基于机器学习的钢筋混凝土柱地震损伤程度评估方法 |
CN114565249A (zh) * | 2022-02-18 | 2022-05-31 | 西安建筑科技大学 | 基于改进熵权-critic法的社区防灾韧性评估方法 |
CN114693123B (zh) * | 2022-03-30 | 2024-08-06 | 中南大学 | 一种既有路基挡土墙安全评估方法 |
CN115018452A (zh) * | 2022-05-20 | 2022-09-06 | 中国长江三峡集团有限公司 | 一种工程建设项目进度管理指标体系构建和综合评价方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104217367A (zh) * | 2014-09-16 | 2014-12-17 | 北京交通大学 | 一种轨道交通线路动态安全风险评价方法 |
CN104361240A (zh) * | 2014-11-19 | 2015-02-18 | 北京市建筑工程研究院有限责任公司 | 基于移动无线技术的火灾后钢结构损伤评估系统 |
CN107742040A (zh) * | 2017-10-31 | 2018-02-27 | 广东电网有限责任公司惠州供电局 | 一种基于topsis和最优组合权重的输电线路综合风险评估方法 |
CN108629064A (zh) * | 2017-03-22 | 2018-10-09 | 安徽省交通控股集团有限公司 | 耦合数值模拟技术的火灾后混凝土桥梁安全状况评定方法 |
CN109740863A (zh) * | 2018-12-13 | 2019-05-10 | 国网山东省电力公司经济技术研究院 | 基于大电源接入系统的综合评价方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103164748A (zh) * | 2011-12-16 | 2013-06-19 | 安徽省电力科学研究院 | 一种电力变电站火灾风险评估方法 |
US10839334B2 (en) * | 2015-07-23 | 2020-11-17 | Tata Consultancy Services Limited | Evaluating performance of organizational units using human capital values |
CN105825320A (zh) * | 2016-03-07 | 2016-08-03 | 天津工业大学 | 一种石油化工企业承灾体脆弱性评估方法 |
CN109377117A (zh) * | 2018-12-21 | 2019-02-22 | 中国水利水电科学研究院 | 一种基于改进topsis法的应急预案评价方法 |
-
2020
- 2020-05-15 CN CN202010413146.4A patent/CN111598448B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104217367A (zh) * | 2014-09-16 | 2014-12-17 | 北京交通大学 | 一种轨道交通线路动态安全风险评价方法 |
CN104361240A (zh) * | 2014-11-19 | 2015-02-18 | 北京市建筑工程研究院有限责任公司 | 基于移动无线技术的火灾后钢结构损伤评估系统 |
CN108629064A (zh) * | 2017-03-22 | 2018-10-09 | 安徽省交通控股集团有限公司 | 耦合数值模拟技术的火灾后混凝土桥梁安全状况评定方法 |
CN107742040A (zh) * | 2017-10-31 | 2018-02-27 | 广东电网有限责任公司惠州供电局 | 一种基于topsis和最优组合权重的输电线路综合风险评估方法 |
CN109740863A (zh) * | 2018-12-13 | 2019-05-10 | 国网山东省电力公司经济技术研究院 | 基于大电源接入系统的综合评价方法 |
Non-Patent Citations (5)
Title |
---|
"AHP+熵权法"的CW-TOPSIS煤矿内因火灾评价模型;秦忠诚 等;《西安科技大学学报》;20180331;第38卷(第2期);第193-201页,引言,1.1节,1.2节,1.3节,第2节,第3.1节,第3.3节,图1、图2,表1 * |
Creep behavior of high-strength concrete with polyp ropylene fibers at elevated temperatures;Wu B,etal;《ACI Materials Journal》;20101231;第107卷(第2期);第176-184页 * |
三标度模糊层次分析法在桥梁安全评估中的应用;常健 等;《中国安全生产科学技术》;20140830;第10卷(第8期);第156-161页 * |
基于GRA-TOPSIS的城市轨道交通PPP项目风险分担;有维宝 等;《土木工程与管理学报》;20180627;第35卷(第3期);第15-21,27页 * |
基于熵权法的煤矿应急救援能力评价;杨力 等;《软科学研究成果与动态》;20131128(第11期);第185-192页,第四节 * |
Also Published As
Publication number | Publication date |
---|---|
CN111598448A (zh) | 2020-08-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111598448B (zh) | 一种基于混凝土t型梁火灾后损伤评估方法 | |
CN111898839A (zh) | 电力用户的重要程度分类方法及装置 | |
CN113935535A (zh) | 一种面向中长期预测模型的主成分分析方法 | |
CN106599927B (zh) | 基于模糊art划分的目标分群方法 | |
CN111105136A (zh) | 一种基于归一化样本的灰色层次评估方法 | |
CN110705856A (zh) | 基于层次分析法的核电设备锻件材料制造质量评估方法 | |
CN106442122A (zh) | 基于图像分割和辨识的钢材料落锤撕裂试验断口韧性断面百分比检测方法 | |
CN104156878B (zh) | 用于农网改造升级工程评估指标权重的确定方法 | |
CN109284574B (zh) | 一种串联桁架结构体系非概率可靠性分析方法 | |
CN113947280A (zh) | 一种基于反馈调整权值的组合评价方法 | |
CN113743662A (zh) | 一种基于机器学习的船舶关键技术参数预测方法及系统 | |
CN112001606A (zh) | 环保工作评价方法及系统 | |
CN108304649B (zh) | 一种高层建筑形变预测方法 | |
CN108051676B (zh) | 一种雷电流幅值累积概率分布曲线拟合计算方法 | |
CN114358612A (zh) | 汽车主观评价方法 | |
CN115496342A (zh) | 一种基于主客观动态权重的冲击地压预警方法和装置 | |
CN110516750B (zh) | 一种基于cart决策树的钢板板形质量异常检测方法 | |
Pryimak et al. | Fuzzy technologies modeling the level of welfare of the population in the system of effective management | |
CN105931133A (zh) | 配电变压器更换优先级评定方法和装置 | |
CN114781829B (zh) | 基于多矩阵ahp的无人机自主作战能力灰色评估方法 | |
Velychko et al. | A comparative analysis of the assessment results of the competence of technical experts by methods of analytic hierarchy process and with using the Rasch model | |
CN111062548A (zh) | 一种企业决策辅助实时数据采集及评价方法 | |
Lee et al. | A Machine Learning Program for Impact Fracture Analysis | |
CN110298590A (zh) | 一种质量基础设施发展水平评估方法 | |
Duan et al. | Selection of dimensionless methods in the back-to-back brand valuation |
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 |