CN112836336A - 根据灰色关联度分析水泥稳定碎石微裂程度中因素的方法 - Google Patents
根据灰色关联度分析水泥稳定碎石微裂程度中因素的方法 Download PDFInfo
- Publication number
- CN112836336A CN112836336A CN202011464055.XA CN202011464055A CN112836336A CN 112836336 A CN112836336 A CN 112836336A CN 202011464055 A CN202011464055 A CN 202011464055A CN 112836336 A CN112836336 A CN 112836336A
- Authority
- CN
- China
- Prior art keywords
- degree
- correlation
- data
- grey
- microcracking
- 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
- 239000004568 cement Substances 0.000 title claims abstract description 45
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000012360 testing method Methods 0.000 claims abstract description 37
- 230000006835 compression Effects 0.000 claims abstract description 23
- 238000007906 compression Methods 0.000 claims abstract description 23
- 238000004458 analytical method Methods 0.000 claims abstract description 21
- 239000011159 matrix material Substances 0.000 claims abstract description 11
- 230000008901 benefit Effects 0.000 claims abstract description 9
- 238000012935 Averaging Methods 0.000 claims abstract 3
- 230000009466 transformation Effects 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims description 2
- 238000011160 research Methods 0.000 abstract description 6
- 230000001419 dependent effect Effects 0.000 abstract description 3
- 238000004451 qualitative analysis Methods 0.000 abstract description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 230000000875 corresponding effect Effects 0.000 description 6
- 238000012163 sequencing technique Methods 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 4
- 238000010219 correlation analysis Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000005336 cracking Methods 0.000 description 3
- 230000001276 controlling effect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000035876 healing Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- BPQQTUXANYXVAA-UHFFFAOYSA-N Orthosilicate Chemical compound [O-][Si]([O-])([O-])[O-] BPQQTUXANYXVAA-UHFFFAOYSA-N 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 239000010426 asphalt Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000011423 initialization method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/38—Concrete; Lime; Mortar; Gypsum; Bricks; Ceramics; Glass
- G01N33/383—Concrete or cement
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- General Physics & Mathematics (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Medicinal Chemistry (AREA)
- Food Science & Technology (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Ceramic Engineering (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明涉及一种根据灰色关联度分析水泥稳定碎石微裂程度中因素的方法,属于土木工程技术领域,主要包括以下步骤:步骤一:通过研究目的确定相应的系统特征变量Yi(k)和相关因素变量Xj(k);步骤二:依据试验数据和相关理论确定是否需要进行正负相关性的转换,并利用算子D0的作用进行相应转换;步骤三:对数据进行无量纲化,通过均值化算子D1作用将其转化为数量级大致相近的无量纲数据,并求各变量序列的均值象;步骤四:根据GRA原理进行灰色关联度的计算;步骤五:针对灰色关联度的数据结果,得到灰色关联矩阵E:E=[ε11 ε12 ε13],根据E中灰色关联度的大小来判断相关因素对于抗压回弹模量的影响程度的大小顺序,并进行优势分析。此方法的优点在于思路明晰,可以在很大程度上减少由于信息不对称带来的损失,并且对数据要求较低,工作量较少,其结果与定性分析结果也比较吻合。
Description
技术领域
本发明涉及一种对于水泥稳定碎石微裂程度中影响影响因素的分析方法,尤其涉及一种基于系列数据以及因素之间发展趋势的相似或相异程度(亦即“灰色关联度”)、通过对数据较为科学的组合变换以寻求微裂程度影响因素的排列顺序。
背景技术
近些年来,我国公路交通事业发展迅速。我国的公路总里程自1996年至2019年,由118.58万公里增长至484.65万公里。其中半刚性基层沥青路面是我国路面结构的主要形式。水泥稳定碎石基层也成为我国高等级公路路面结构的主要基层形式,但这种路面在使用初期极易产生收缩裂缝,而后会扩展导致反射性裂缝的出现,使路面加速破坏。其关键之处在于早期的收缩应力。为此,许多专家进行对级配、水泥用量、外加剂含量等进行研究,而通过微裂技术来降低早期收缩应力的研究不多。微裂技术是基于工程措施来减小收缩应力的有效措施,应用前景广阔。通过研究发现,微裂后水泥稳定碎石基层材料力学强度的愈合过程是一个复杂的过程,影响其力学强度恢复或微裂缝自愈合能力的因素主要有水泥稳定碎石材料结构类型、水泥剂量、实施微裂的时间、微裂荷载以及愈合时间等,在进行相关研究后需对混合料的水泥剂量、微裂时间以及微裂荷载这三个主要影响因素进行有效分析。为了得到不同因素的影响程度,通常对实验数据进行分析,常见的分析方法有层次分析法、模糊综合评价以及利用SPSS显著性分析等。
为了得到不同因素的影响程度,通常对实验数据进行SPSS显著性分析,河北工业大学李娜、魏连雨在《基于微观结构分析的水泥稳定碎石微裂自愈性能》中(李娜,魏连雨.基于微观结构分析的水泥稳定碎石微裂自愈性能.《硅酸盐通报》2018年第2期)通过SPSS显著性分析方法确定水泥稳定碎石微裂程度的主要影响因素。但在其分析过程中,仅仅可以表明微裂时间和微裂荷载对试验结果影响显著,并未在两者之间确定一个最优因素。因此,对于该研究,SPSS显著性分析有一定的局限性。
而灰色关联分析是按照发展趋势做分析,通过计算得出的关联度对相关因素变量进行排序。因此对样本量的多少没有过多要求,而且计算量比较小,其结果与定性分析结果也比较吻合。因此,灰色关联分析是一种比较简单而又可靠的方法。
发明内容
本发明的目的是提供一种可以分析水泥稳定碎石微裂影响因素程度大小的方法。该方法是基于灰色关联分析法原理(GRA),并根据研究目的确定系统特征变量和相关因素变量后,通过对数据的正负相关转换、无量纲化,进而求得相应的灰色关联度,最后按照灰色关联度进行排序,得到水泥稳定碎石微裂程度影响因素的排序。
为实现上述目的,本发明的技术方案是:
一种根据灰色关联度分析水泥稳定碎石微裂程度中因素的方法,该方法包括以下步骤:
步骤一、通过研究目的确定相应的系统特征变量Y1(k)和相关因素变量X1(k)、X2(k)、X3(k)。其中Y1(k)表示系统特征量抗压回弹模量在第k次试验时的测试数据,X1(k)表示水泥用量在第k次试验时的测试数据,X2(k)表示水泥用量在第k次试验时的测试数据,X3(k)表示水泥用量在第k次试验时的测试数据。那么试验中所测得的相应数据就可以组成系统特征变量数据序列和相关因素数据数列,即:
Yi=[yi(1),yi(2),…,yi(k),…,yi(27)]
i=1
k=1,2,…,27
Xj=[xj(1),xj(2),…,xj(k),…,xj(27)]
j=1,2,3
k=1,2,…,27
式中:i表示第i个系统特征变量,j表示第j个相关因素变量,k表示第k次试验。
步骤二、依据试验数据和相关理论确定微裂荷载和抗压回弹模量呈负相关关系,则需要进行正负相关性的转换:
若X3(k)和Y1(k)呈负相关关系,通过算子D0的作用,使得X3(k)的对称化像X3D0与Y1(k)具有正相关关系。
X3D0=[x3(1)d0,x3(2)d0,…,x3(k)d0]
式中:x3(k)d0=2max(x3)-x3(k),k=1,2,3,…,27,则称D0为对称化算子,X3D0为相应对称化像。
步骤三、对数据进行无量纲化,通过均值化算子D1作用将水泥用量、微裂时间、微裂荷载对应试验数据转化为数量级大致相近的无量纲数据,并求各变量序列的均值象。
k=1,2,…,27
i=1
j=1,2,3
步骤四、进行相应的灰色关联度的计算,根据GRA原理进行如下运算:
j=1,2,3
k=1,2,…,27
根据灰色关联度定义,
步骤五、进行优势分析。
针对灰色关联度的数据结果,得到灰色关联矩阵E:
E=[ε11 ε12 ε13]
根据E中灰色关联度的大小来判断相关因素对于抗压回弹模量的影响程度的大小顺序,并进行优势分析。
进一步地,在步骤一中,系统特征变量在本实例中仅有抗压回弹模量Y1(k),相关因素变量有三个,分别是水泥用量X1(k)、微裂时间X2(k)、微裂荷载X3(k)。本实验的目的在于探寻影响水泥微裂的三个相关因素对于反应微裂程度的抗压回弹模量的影响程度。
进一步地,在步骤二中,首先是根据试验结果分析出微裂荷载X3(k)与抗压回弹模量Y1(k)具有负相关关系,然后对其进行转化。
进一步地,在步骤三中,进行无量纲化的常用方法有初值化法,均值化法和区间值化法等,在此次分析中,运用均值化法对27组试验数据进行处理。即先算出每个变量相对应的27组试验数据的算术平均值,然后用各个值除以平均值,最后得出均值象矩阵。
进一步地,在步骤四中,进行灰色关联度计算时,因只有一个系统特征变量,所以只有一组两极最大差和最小差。
在步骤五中,按照计算得出的灰色关联度来分析水泥用量、微裂时间和微裂荷载对于抗压回弹模量的影响程度的大小。
与现有技术相比,本发明的有益效果在于:
本发明基于灰色关联法对水泥微裂实验繁杂的试验数据进行一系列处理,最终整理为灰色关联度矩阵,可以明显看出微裂程度相关因素间的关联性,进而得出相关因素影响程度的大小排序。
在SPSS显著性分析下,虽然能分析出微裂中某些相关因素的显著性相同,但是P值均为0的情况下,很难在多个极为显著的相关因素中分析出最为重要的优势因素,也不能给出一个微裂影响因素程度大小的排序。本发明基于灰色关联理论进一步得出水泥用量、微裂时间和微裂荷载对于抗压回弹模量的排序,突破了不能比较各因素影响程度的局限。在分析过程中对样本量的多少没有过多要求,而且计算量比较小,其结果与定性分析结果也比较吻合。通过灰色关联分析的结果可以得到影响程度大小的排序,进而可以在实践中控制其最优或准优因素来显著提高实际工程的效率。
附图说明
图1是本发明的流程图。
具体实施方式
下面以基于灰度分析原理分析水泥稳定碎石微裂程度因素为例,对本发明作进一步说明。
首先整合试验数据,进行数据集的建立。对于该平行实验而言,共有三组,所以采用抗压回弹模量的算术平均值,并确定抗压回弹模量为本次试验中唯一的系统特征变量Y1。实验采用结构类型为骨架密实结构的混合料,进一步分析确定水泥用量、微裂时间和微裂荷载为相关因素变量X1、X2、X3,水泥用量分别为3%、4%、5%。微裂时间分别为1d、2d和3d,微裂荷载主要取决于振动时间,分别为90s、120s、150s。
根据各变量定义,以水泥稳定碎石抗压回弹模量为研究对象进行室内试验,将试验数据整理为表格,如表1所示:
表1
通过观察试验数据,发现微裂荷载X3和抗压回弹模量Y1呈负相关,而水泥用量和微裂时间均与抗压回弹模量呈正相关关系,此时需要利用对称化算子D0的作用,使微裂荷载X3与抗压回弹模量Y1也呈正相关。
X3D0=[4.543.54.543.54.543.5…]
将X3D0代替相关因素变量矩阵中的X3后,采用均值化法对相关因素变量矩阵和系统特征变量矩阵进行无量纲化,得到均值象矩阵X1和Y1,将其用表格表示,见表2。
表2
表3
通过分析表格3数据,
取分辨系数ρ=0.5,根据公式
表4
继而,根据公式:
表5
根据所得GRA结果,对灰色关联度进行排序和优势分析,从大到小优势因素依次为:微裂时间、微裂荷载和水泥用量。并且从排序来看,最优因素为微裂时间,可以通过主控微裂时间来增强其抗压强度模量。相比于SPSS显著性分析只能得出微裂时间和微裂荷载的影响都比较显著,本发明更加区分出微裂时间才是主控因素,也是增加抗压回弹模量的关键之处。
综上所述,运用本发明可以得到关于水泥稳定碎石微裂程度影响因素的排序从大到小为:微裂时间、微裂荷载和水泥用量,并从中分析出最优因素为微裂时间。因此可以通过控制微裂时间来有效改善微裂效果,从而提升水泥稳定碎石的强度,这为工程中的优化提供了重要的依据和方向。
尽管已经示出和描述了本发明的实施例,对于本领域的普通技术人员而言,可以理解在不脱离本发明的原理和精神的情况下可以对这些实施例进行多种变化、修改、替换和变型,本发明的范围由所附权利要求及其等同物限定。
Claims (5)
1.一种根据灰色关联度分析水泥稳定碎石微裂程度中因素的方法,其特征在于:包括以下步骤:
步骤一、确定相关因素变量X1(k)、X2(k)、X3(k)和系统特征变量Y1(k),Y1(k)表示系统特征量抗压回弹模量在第k次试验时的测试数据,系统特征变量为抗压回弹模量,相关因素变量有三个,分别为水泥用量、微裂时间、微裂荷载,X1(k)表示水泥用量在第k次试验时的测试数据,X2(k)表示水泥用量在第k次试验时的测试数据,X3(k)表示水泥用量在第k次试验时的测试数据,那么试验中所测得的相应数据组成系统特征变量数据序列和相关因素数据数列,即:
Yi=[yi(1),yi(2),…,yi(k),…,yi(27)]
i=1
k=1,2,…,27
Xj=[xj(1),xj(2),…,xj(k),…,xj(27)]
j=1,2,3
k=1,2,…,27
式中:i表示第i个系统特征变量,j表示第j个相关因素变量,k表示第k次试验;
步骤二、依据试验数据和理论确定微裂荷载和抗压回弹模量呈负相关关系,进行正负相关性的转换:
若X3(k)和Y1(k)呈负相关关系,通过算子D0的作用,使得X3(k)的对称化像X3D0与Y1(k)具有正相关关系;
X3D0=[x3(1)d0,x3(2)d0,…,x3(k)d0]
式中:x3(k)d0=2max(x3)-x3(k),k=1,2,3,…,27,称D0为对称化算子,X3D0为相应对称化像;
步骤三、对数据进行无量纲化,通过均值化算子D1作用将水泥用量、微裂时间、微裂荷载对应试验数据转化为无量纲数据,并求各变量序列的均值象;
k=1,2,…,27
i=1
j=1,2,3;
步骤四、进行相应的灰色关联度的计算,根据GRA原理进行如下运算:
j=1,2,3
k=1,2,…,27
根据灰色关联度定义,
步骤五、进行优势分析。
针对灰色关联度的数据结果,得到灰色关联矩阵E:
E=[ε11 ε12 ε13]
根据E中灰色关联度的大小来判断相关因素对于抗压回弹模量的影响程度的大小顺序,并进行优势分析。
2.根据权利要求1所述的一种根据灰色关联度分析水泥稳定碎石微裂程度中因素的方法,其特征在于:所述步骤二中先根据试验结果分析出微裂荷载与抗压回弹模量具有负相关关系,然后进行转化。
3.根据权利要求2所述的一种根据灰色关联度分析水泥稳定碎石微裂程度中因素的方法,其特征在于:所述步骤三中进行无量纲化的常用方法有初值化法,均值化法和区间值化法,运用均值化法对27组试验数据进行处理,先算出每个变量相对应的27组试验数据的算术平均值,后用各个值除以平均值,最后得出均值象矩阵。
4.根据权利要求3所述的一种根据灰色关联度分析水泥稳定碎石微裂程度中因素的方法,其特征在于:对所述步骤四进行灰色关联度计算时,系统特征变量为1个,两极最大差和最小差为一组。
5.根据权利要求1或3或4所述的一种根据灰色关联度分析水泥稳定碎石微裂程度中因素的方法,其特征在于:对所述步骤五按照计算得出的灰色关联度来分析水泥用量、微裂时间和微裂荷载对于抗压回弹模量的影响程度的大小。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011464055.XA CN112836336B (zh) | 2020-12-10 | 2020-12-10 | 根据灰色关联度分析水泥稳定碎石微裂程度中因素的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011464055.XA CN112836336B (zh) | 2020-12-10 | 2020-12-10 | 根据灰色关联度分析水泥稳定碎石微裂程度中因素的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112836336A true CN112836336A (zh) | 2021-05-25 |
CN112836336B CN112836336B (zh) | 2022-03-25 |
Family
ID=75923533
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011464055.XA Active CN112836336B (zh) | 2020-12-10 | 2020-12-10 | 根据灰色关联度分析水泥稳定碎石微裂程度中因素的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112836336B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113506599A (zh) * | 2021-07-19 | 2021-10-15 | 中国水利水电科学研究院 | 一种胶结坝数字化拌和与智能化动态调控系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108197403A (zh) * | 2018-01-19 | 2018-06-22 | 扬州大学 | 基于灰色关联分析法的沥青路面最佳面层厚度组合确定方法 |
-
2020
- 2020-12-10 CN CN202011464055.XA patent/CN112836336B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108197403A (zh) * | 2018-01-19 | 2018-06-22 | 扬州大学 | 基于灰色关联分析法的沥青路面最佳面层厚度组合确定方法 |
Non-Patent Citations (2)
Title |
---|
SUN Z H等: "Grey Relational Analysis of Fatigue Performance of Semi-rigid Pavement Structure", 《APPLIED MECHANICS AND MATERIALS》 * |
马士宾等: "微裂后水泥粉煤灰稳定碎石力学性能变化规律研究", 《硅酸盐通报》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113506599A (zh) * | 2021-07-19 | 2021-10-15 | 中国水利水电科学研究院 | 一种胶结坝数字化拌和与智能化动态调控系统 |
CN113506599B (zh) * | 2021-07-19 | 2024-03-22 | 中国水利水电科学研究院 | 一种胶结坝数字化拌和与智能化动态调控系统 |
Also Published As
Publication number | Publication date |
---|---|
CN112836336B (zh) | 2022-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108960651B (zh) | 一种致密油气藏多级压裂水平井完井效率的综合评价方法 | |
Ghosh et al. | Particle fracture simulation in non-uniform microstructures of metal–matrix composites | |
CN103630450A (zh) | 考虑疲劳-蠕变交互损伤作用的沥青混合料寿命预测方法 | |
CN112836336B (zh) | 根据灰色关联度分析水泥稳定碎石微裂程度中因素的方法 | |
Ren et al. | Recognition of shear and tension signals based on acoustic emission parameters and waveform using machine learning methods | |
CN110658086A (zh) | 一种考虑拉压模量差异性的沥青路面荷载响应分析方法 | |
CN112070788B (zh) | 基于分块梯度分割的快速统计形变孪晶的图像处理方法 | |
Chen et al. | Identification of fracture damage characteristics in ultra-high performance cement-based composite using digital image correlation and acoustic emission techniques | |
Ma et al. | Experimental study on the influence of height and dip angle of asperity on the mechanical properties of rock joints | |
CN111062078B (zh) | 适用于混合型再生粗骨料的品质分级方法 | |
CN107923830A (zh) | 材料的断裂表面的开裂机制的表征方法 | |
CN109284574B (zh) | 一种串联桁架结构体系非概率可靠性分析方法 | |
Yan | Microdefect interacting with a finite main crack | |
CN116486960A (zh) | 一种纳米增强水泥基材料智能化设计方法及装置 | |
CN108333335A (zh) | 一种混凝土梁桥抗力折减系数确定方法 | |
CN118347882B (zh) | 一种沥青混合料疲劳寿命预测方法 | |
Fang | Crack repair model of ancient ceramics based on digital image | |
CN111339497A (zh) | 一种评价再生混凝土强度指标的概率预测模型 | |
CN1687935A (zh) | 一种极少新钢种数据样本下的质量设计方法 | |
Özen | Investigation of relationship between aggregate shape parameters and concrete strength using imaging techniques | |
CN112485253A (zh) | 玄武岩纤维沥青混合料中纤维最佳长度组合比的确定方法 | |
CN105528657A (zh) | 基于北斗和向量机的建筑物震害预测方法 | |
Yong et al. | A real-time detection algorithm based on nanodet for pavement cracks by incorporating attention mechanism | |
CN110457865A (zh) | 基于数字散斑方法的离散元图像建模方法 | |
Wang et al. | Crack pattern identification in cementitious materials based on acoustic emission and machine learning |
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 |