CN112699529B - 一种沥青混合料水气运动参数计算方法、存储介质及系统 - Google Patents

一种沥青混合料水气运动参数计算方法、存储介质及系统 Download PDF

Info

Publication number
CN112699529B
CN112699529B CN202011434387.3A CN202011434387A CN112699529B CN 112699529 B CN112699529 B CN 112699529B CN 202011434387 A CN202011434387 A CN 202011434387A CN 112699529 B CN112699529 B CN 112699529B
Authority
CN
China
Prior art keywords
asphalt mixture
water vapor
water
test piece
mass
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
Application number
CN202011434387.3A
Other languages
English (en)
Other versions
CN112699529A (zh
Inventor
黄婷婷
罗蓉
涂崇志
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University of Technology WUT
Original Assignee
Wuhan University of Technology WUT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Wuhan University of Technology WUT filed Critical Wuhan University of Technology WUT
Priority to CN202011434387.3A priority Critical patent/CN112699529B/zh
Publication of CN112699529A publication Critical patent/CN112699529A/zh
Application granted granted Critical
Publication of CN112699529B publication Critical patent/CN112699529B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • 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
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • General Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Sampling And Sample Adjustment (AREA)

Abstract

本发明涉及一种沥青混合料水气运动参数计算方法，其包括步骤，制备并成型圆柱体的待测沥青混合料试件；对沥青混合料试件进行干燥处理，并置于磁悬浮重量平衡系统中开展积聚型水气运动试验，获取沥青混合料试件吸收水气质量随时间的变化数据；基于获取的沥青混合料试件吸收水气质量随时间的变化数据，建立描述沥青混合料内部水气扩散和吸附的三维水气运动模型。本发明还提供一种存储介质及沥青混合料水气运动参数计算系统，本发明提供的沥青混合料水气运动参数计算方法、存储介质及系统可更加准确的描述沥青混合料内水气运动特性的理论计算方法，从而为研究沥青混合料水损害机理提供理论基础。

Description

一种沥青混合料水气运动参数计算方法、存储介质及系统
技术领域
本发明涉及道路工程领域，尤其涉及一种沥青混合料水气运动参数计算方法、存储介质及系统。
背景技术
水损害作为沥青路面早期病害形式之一，严重影响道路的服役寿命，造成社会和经济损失，因此研究沥青路面水损害形成原因及评价指标，对指导沥青路面设计与施工具有十分重要的意义。现阶段沥青路面水损害研究主要集中在液态水引发的各种破坏上，国内外沥青混合料设计方法中以沥青混合料浸水或者冻融前后的力学指标比作为评价其水稳定性的指标。但是大量干旱地区的工程实例标明，水气也是引发沥青路面的水损害重要原因。干旱地区沥青路面周围环境和路基、土基内部水分以气态形式在沥青混合料中运动并积聚是引发水损害的关键因素，水分子在外界湿度、温度、混合料内部结构等因素的影响下气化并扩散，在此过程中水气不断携带水分子进入到沥青混合料内部，附着在沥青膜上或穿透沥青膜进入到沥青与集料界面间，从而降低沥青膜自身的内聚结合能或沥青与集料界面之间的粘附结合能，引起沥青膜的内聚破坏或沥青与集料之间的粘聚破坏，进而导致沥青路面产生水损害。
已有相关研究基于理想的菲克扩散，建立了严谨的三维水气运动模型，提出采用扩散系数和保水度作为评价沥青混合料水气运动特性的指标，但是这一体系的建立前提是沥青混合料内部的水气运动符合理想菲克扩散，即不发生水分子吸附与聚集，这与沥青混合料中的水气存在形式存在一定的差异。
为此，现阶段急需一种更加准确的描述沥青混合料内水气运动特性的理论计算方法，从而为研究沥青混合料水损害机理提供理论基础。
发明内容
有鉴于此，本发明提供一种沥青混合料水气运动参数计算方法、存储介质及系统解决准确的描述沥青混合料内水气运动特性的理论的问题。
为了达到上述目的，本发明解决技术问题的技术方案是提供一种沥青混合料水气运动参数计算方法，其包括步骤：制备并成型圆柱体的待测沥青混合料试件；对沥青混合料试件进行干燥处理，并置于磁悬浮重量平衡系统中开展积聚型水气运动试验，获取沥青混合料试件吸收水气质量随时间的变化数据；基于获取的沥青混合料试件吸收水气质量随时间的变化数据，建立描述沥青混合料内部水气扩散和吸附的三维水气运动模型。
进一步，所述对沥青混合料试件进行干燥处理，并置于磁悬浮重量平衡系统中开展积聚型水气运动试验，获取沥青混合料试件吸收水气质量随时间的变化数据包括步骤：进行试验准备；进行试验；对获取的数据进行处理。
进一步，所述进行试验包括步骤：称取真空环境下样品桶的质量；将沥青混合料试件放置在样品桶里，对反应腔进行抽真空处理，并不断称量样品桶和试件质量，直至质量保持恒定；设置试验温度和水气压力。
进一步，所述基于获取的沥青混合料试件吸收水气质量随时间的变化数据，建立描述沥青混合料内部水气扩散和吸附的三维水气运动模型包括步骤：通过移动态水分子满足菲克扩散定律得到三维水气扩散基本模型；确定初始条件和边界条件；进行拉普拉斯变换和逆变换；建立三维水气扩散模型。
进一步，所述建立三维水气扩散模型包括步骤：分离变量；施加边界条件和初始条件；积分求解水气质量。
进一步，所述对获取的数据进行处理指通过对试验数据进行处理，得到沥青混合料中的水气扩散、积聚等参数。
本发明还提供一种存储介质，所述存储介质中存储有计算机程序，其中，所述计算机程序被设置为运行时执行沥青混合料水气运动参数计算方法。
本发明还提供一种沥青混合料水气运动参数计算系统，其包括处理器和存储器,所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行时，实现沥青混合料水气运动参数计算方法。
与现有技术相比，本发明所提供的沥青混合料水气运动参数计算方法、存储介质及系统具有以下有益效果：
通过基于获取的沥青混合料试件吸收水气质量随时间的变化数据，建立描述沥青混合料内部水气扩散和吸附的三维水气运动模型，从而更加准确的描述沥青混合料内水气运动特性的理论计算方法，从而为研究沥青混合料水损害机理提供理论基础。
以上所述本发明的具体实施方式，并不构成对本发明保护范围的限定。任何根据本发明的技术构思所做出的各种其他相应的改变与变形，均应包含在本发明的保护范围内。
附图说明
图1为本发明第一实施例提供的一种沥青混合料水气运动参数计算方法的步骤流程示意图；
图2为图1中步骤S2的子步骤流程图；
图3为图2中步骤S22的子步骤流程图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
请参阅图1，本发明提供的一种沥青混合料水气运动参数计算方法，其包括步骤：
S1，制备并成型圆柱体的待测沥青混合料试件；
具体的，基于道路沥青混合料设计方法，设计、制备并成型圆柱体的待测沥青混合料试件，记试件截面半径为a，高度为H。
S2，对沥青混合料试件进行干燥处理，并置于磁悬浮重量平衡系统中开展积聚型水气运动试验，获取沥青混合料试件吸收水气质量随时间的变化数据；
具体的，对成型的沥青混合料样品进行干燥处理，直至样品质量无变化，此时认定试件初始相对湿RH1度为0；将沥青混合料样品置于磁悬浮重量平衡系统开展积聚型水气运动试验，在试验装置内通入设定大小的水气压（对应相对湿度为RH2），水气进入样品后样品的质量由试验装置的高精度磁悬浮天平称量得到，称量频率为5次/分钟，最终获得沥青混合料试件吸收水气质量随时间的变化数据。
S3，基于获取的沥青混合料试件吸收水气质量随时间的变化数据，建立描述沥青混合料内部水气扩散和吸附的三维水气运动模型；
具体的，建立的模型为：
Figure SMS_1
(1)
或者，
Figure SMS_2
(2)
式中：M(t)=时刻t时沥青混合料样品内部积聚的水气质量；
Figure SMS_3
，表示在温度T和纯水气压P环境下沥青混合料样品能够容纳的最大水气质量，其中C0为水气运动试验中沥青混合料样品表面恒定的水气浓度；γ=扩散进入沥青混合料试件中的水气分子从移动态变成束缚态的概率，s-1；β =同一时刻水气分子从束缚态变成移动态的概率，s-1；/>
Figure SMS_4
=第一类零阶贝塞尔函数的第m个零点；D=移动态水气扩散系数；t=扩散时间。
在本实施例中，推导方法如下：
假设沥青混合料试件中的移动态水气分子的扩散满足菲克扩散定律，试验采用圆柱体试件，故在柱坐标系下，沥青混合料试件中的移动态水分子浓度n (r, z, t)和束缚态水分子浓度N (r, z, t)满足下面的三维水气扩散基本模型：
Figure SMS_5
(3)
Figure SMS_6
(4)
Figure SMS_7
(5)
式中：γ、β、t、r、z、D含义同前；n = t时刻沥青混合料试件中移动态水分子的浓度，g/mm3；N = t时刻沥青混合料试件中束缚态水分子的浓度，g/mm3；n∞、N∞ =时间无穷大时移动态和束缚态水分子的浓度，g/mm3。
基于本发明开展的积聚型水气运动试验方案，可以确定初始条件和边界条件：
（1）初始条件：
Figure SMS_8
，/>
Figure SMS_9
(6)
Figure SMS_10
，/>
Figure SMS_11
(7)
（2）边界条件：
Figure SMS_12
，/>
Figure SMS_13
(8)
Figure SMS_14
，/>
Figure SMS_15
(9)
Figure SMS_16
，/>
Figure SMS_17
(10)
式中：a=沥青混合料试件的半径，mm；H=沥青混合料试件的高度，mm；
Figure SMS_18
=任意时间t沥青混合料试件侧表面浓度，g/mm3；/>
Figure SMS_19
=任意时间t时，沥青混合料试件下表面的水气浓度，g/mm3；/>
Figure SMS_20
=任意时间t时，沥青混合料试件上表面的水气浓度，g/mm3；/>
Figure SMS_21
=水气运动试验中沥青混合料试件周边的恒定水气浓度，g/mm3。
拉普拉斯变换和逆变换：
结合公式4将公式3符号统一、进行简化，可得如下扩散模型：
Figure SMS_22
(11)
则初始条件变化成：
Figure SMS_23
，/>
Figure SMS_24
(12)
Figure SMS_25
，/>
Figure SMS_26
(13)
同理简化边界条件成为：
Figure SMS_27
(14)
Figure SMS_28
(15)
Figure SMS_29
(16)
首先采用拉普拉斯变换对边界条件进行转换：
Figure SMS_30
，/>
Figure SMS_31
(17)
Figure SMS_32
，/>
Figure SMS_33
(18)
Figure SMS_34
，/>
Figure SMS_35
(19)
再采用拉普拉斯逆变换，可得边界条件如下：
Figure SMS_36
；/>
Figure SMS_37
(20)
Figure SMS_38
，/>
Figure SMS_39
(21)
Figure SMS_40
，/>
Figure SMS_41
(22)
基于上述步骤，详细的模型求解过程如下：
（1）分离变量：
为了求解模型公式11，首先对边界条件式20~22进行齐次化，即令：
Figure SMS_42
(23)
公式11变成：
Figure SMS_43
(24)
边界条件成为：
Figure SMS_44
；/>
Figure SMS_45
(25)
Figure SMS_46
，/>
Figure SMS_47
(26)
Figure SMS_48
，/>
Figure SMS_49
(27)
初始条件成为：
Figure SMS_50
(28)
Figure SMS_51
(29)
接下来进行求解，首先采用分离变量法：
Figure SMS_52
(30)
则公式24成为：
Figure SMS_53
(31)
上式两边仅在同时等于一个常数下成立，假设这个常数为
Figure SMS_54
，则有：
Figure SMS_55
(32)
Figure SMS_56
(33)
首先求解T (t)的表达式，公式32重新整理后有：
Figure SMS_57
(34)
公式34的特征方程为：
Figure SMS_58
(35)
公式35的解为：
Figure SMS_59
(36)
所以公式34的解表示如下：
Figure SMS_60
(37)
接着继续求解Z (z)和R (r)，公式33重新整理如下：
Figure SMS_61
(38)
上式38仅在等号两边同时等于一个常数时成立，假设该常数为
Figure SMS_62
，则有
Figure SMS_63
(39)
Figure SMS_64
(40)
式39的通解为：
Figure SMS_65
(41)
式40的通解为：
Figure SMS_66
(42)
因此
Figure SMS_67
的通解可以表示如下：
Figure SMS_68
(43)
式中：
Figure SMS_69
,/>
Figure SMS_70
，/>
Figure SMS_71
Figure SMS_72
(44)
Figure SMS_73
(45)
（2）施加边界条件和初始条件
①边界条件1（公式25）：
Figure SMS_74
Figure SMS_75
(46)
由于对任意的z、t上式都成立，因此有：
Figure SMS_76
，因此/>
Figure SMS_77
是/>
Figure SMS_78
的解。
②边界条件2（公式26）：
Figure SMS_79
Figure SMS_80
(47)
因为对任意的r、t都成立，因此有：
Figure SMS_81
③边界条件3（公式27）：
Figure SMS_82
Figure SMS_83
(48)
因为对任意r、t都成立，因此有：
Figure SMS_84
，即/>
Figure SMS_85
，因此
Figure SMS_86
Figure SMS_87
(49)
将前面确定的
Figure SMS_88
，/>
Figure SMS_89
代入到/>
Figure SMS_90
中，可以得到：
Figure SMS_91
(50)
式中：
Figure SMS_92
，/>
Figure SMS_93
④初始条件1（公式28）：
Figure SMS_94
Figure SMS_95
(51)
⑤初始条件式2（公式29）：
Figure SMS_96
Figure SMS_97
(52)
因为上式对任意r、z都成立，因此有：
Figure SMS_98
(53)
由公式51可知，等号两侧对任意r和z都成立，因此假设等号右侧的常数为r、z的函数：
Figure SMS_99
(54)
式中：
Figure SMS_100
是r的函数；/>
Figure SMS_101
是z的函数。对这两个函数进行贝塞尔和傅里叶展开：
Figure SMS_102
(55)
Figure SMS_103
(56)
式中的
Figure SMS_104
、/>
Figure SMS_105
为未知参数，求解过程如下：
①公式55两侧同乘以
Figure SMS_106
，并对r从0~a积分，结合贝塞尔函数性质得到：
Figure SMS_107
(57)
②公式56两侧同时乘以
Figure SMS_108
，并对z从0~H积分，结合三角函数性质得到:
Figure SMS_109
(58)
结合
Figure SMS_110
、/>
Figure SMS_111
及/>
Figure SMS_112
和/>
Figure SMS_113
表达式，可以确定/>
Figure SMS_114
表示如下：
Figure SMS_115
(59)
将公式59与公式51进行对比，可得：
Figure SMS_116
(60)
考虑贝塞尔函数及三角函数的如下性质：
Figure SMS_117
，/>
Figure SMS_118
因此公式60可以简化：
Figure SMS_119
(61)
结合式53，可以得到
Figure SMS_120
和/>
Figure SMS_121
的表达式如下：
Figure SMS_122
(62)
Figure SMS_123
(63)
因此
Figure SMS_124
的表达式如下：
Figure SMS_125
(64)
式中：
Figure SMS_126
,/>
Figure SMS_127
。
因此根据公式23，可以得到
Figure SMS_128
的表达式：
Figure SMS_129
(65)
式中：
Figure SMS_130
,/>
Figure SMS_131
。
将
Figure SMS_132
表达式代入，可得：
Figure SMS_133
(66)
Figure SMS_134
(67)
式中：
Figure SMS_135
,/>
Figure SMS_136
。
（3）积分求解水气质量
束缚态水分子质量
Figure SMS_137
可以通过对/>
Figure SMS_138
积分得到：
Figure SMS_139
(68)
式中：
Figure SMS_140
,/>
Figure SMS_141
。
上式中当n=偶数时
Figure SMS_142
，n=奇数时/>
Figure SMS_143
，因此公式68可以简化为：
Figure SMS_144
(69)
式中：
Figure SMS_145
,/>
Figure SMS_146
。
对
Figure SMS_147
进行积分，可以得到移动态水分子的质量/>
Figure SMS_148
：
Figure SMS_149
(70)
式中：
Figure SMS_150
,/>
Figure SMS_151
。
因此进入沥青混合料试件中的总质量
Figure SMS_152
为：
Figure SMS_153
(71)
式中的
Figure SMS_154
和/>
Figure SMS_155
按照公式44、45计算，这两个参数是公式35的特征解，满足下列等式：
Figure SMS_156
(72)
Figure SMS_157
(73)
因为
Figure SMS_159
，且/>
Figure SMS_162
，因此/>
Figure SMS_164
和/>
Figure SMS_160
均为负值，所以当/>
Figure SMS_163
，/>
Figure SMS_165
和/>
Figure SMS_166
Figure SMS_158
，因此当/>
Figure SMS_161
时，
Figure SMS_167
(74)
则任意时间t时，
Figure SMS_168
如下：
Figure SMS_169
(75)/>
假设γ和β在数值上与扩散系数相近，考虑到扩散系数的数量级，则γβ可以认为远小于
Figure SMS_170
，因此公式44、45可以重新写成如下形式：
Figure SMS_171
(76)
Figure SMS_172
(77)
Figure SMS_173
(78)
因此，将公式76~78代入公式75中，可以得到：
Figure SMS_174
(79)
因为
Figure SMS_175
，/>
Figure SMS_176
，/>
Figure SMS_177
，随着m、n增大，如下等式成立：
Figure SMS_178
(80)
Figure SMS_179
(81)
所以公式79可以简化为下式：
Figure SMS_180
(82)
因为
Figure SMS_181
是/>
Figure SMS_182
的根，前文已经定义/>
Figure SMS_183
满足/>
Figure SMS_184
，因此将
Figure SMS_185
代入上式中，可以得到：
Figure SMS_186
(83)
当时间t = 0时，可以认为
Figure SMS_187
=0，即上式为0：
Figure SMS_188
(84)
即有：
Figure SMS_189
(85)
易知：
Figure SMS_190
，/>
Figure SMS_191
，所以有：
Figure SMS_192
(86)
将公式86代入公式83中，同时将
Figure SMS_193
、/>
Figure SMS_194
、/>
Figure SMS_195
（
Figure SMS_196
）代入，可得沥青混合料内部的水气运动与积聚模型：
Figure SMS_197
(87)
也可以表述成下式：
Figure SMS_198
(88)
根据上述步骤，完成沥青混合料内部的水气运动与积聚模型的创建。
请参阅图2，步骤S2包括子步骤：
S21，进行试验准备；
具体的，试验准备主要包括沥青混合料试验试件准备，进行沥青混合料级配设计，成型、切割沥青混合料试件，得到直径12mm、高度20mm的沥青混合料圆柱体试件，并进行低温烘干，直至质量恒定；制备纯蒸馏水，在蒸馏水小瓶内注入10ml蒸馏水，利用液氮对蒸馏水小瓶,冷冻至冰点以下，真空泵对蒸馏水小瓶内存在的其他气体抽真空，获得纯蒸馏水，便于后续得到纯水蒸气。
S22，进行试验；
具体的，设置试验温度为20℃，称取真空环境下样品桶的质量M1；将沥青混合料试件放置在样品桶里，对反应腔进行抽真空处理，并不断称量样品桶和试件质量，直至质量保持恒定，记此时干燥沥青混合料和样品桶的质量为M2，则干燥沥青混合料试件的质量M0=M2-M1，初始干燥沥青混合料试件内部相对湿度RH1为0%；设置试验温度和水气压力（对应相对湿度RH2），试验方案如表1所示，试验过程中维持反应腔内温度和水气压的大小保持不变，水分子不断进入沥青混合料试件，磁悬浮天平对沥青混合料试件和样品桶进行称量，获得各时间点t下沥青混合料试件和样品桶质量为M3，则沥青混合料试件吸收的水气质量为M(t)=M3-M0，绘制沥青混合料吸收水气质量随时间的变化曲线，以17.17%相对湿度下的吸收水气质量随时间变化数据为例。
表1 水气运动试验方案
Figure SMS_199
注：20℃水的饱和蒸气压为23.3mbar，相对湿度定义为试验温度下施加水气压与饱和蒸气压的比值。
S23，对获取的数据进行处理；
具体的，沥青混合料试件吸收的水气质量M(t)满足公式86所示的模型，采用该模型对试验所得数据进行处理，得到沥青混合料中的水气扩散、积聚等参数，如表2所示。
表2 各相对湿度下沥青混合料水气运动参数汇总表
Figure SMS_200
请参阅图3，步骤S22还包括子步骤：
S221，称取真空环境下样品桶的质量；
具体的，设置试验温度为20℃，称取真空环境下样品桶的质量M1。
S222，将沥青混合料试件放置在样品桶里，对反应腔进行抽真空处理，并不断称量样品桶和试件质量，直至质量保持恒定；
具体的，将沥青混合料试件放置在样品桶里，对反应腔进行抽真空处理，并不断称量样品桶和试件质量，直至质量保持恒定，记此时干燥沥青混合料和样品桶的质量为M2，则干燥沥青混合料试件的质量M0=M2-M1，初始干燥沥青混合料试件内部相对湿度RH1为0%。
S223，设置试验温度和水气压力；
具体的，设置试验温度和水气压力（对应相对湿度RH2），试验过程中维持反应腔内温度和水气压的大小保持不变，水分子不断进入沥青混合料试件，磁悬浮天平对沥青混合料试件和样品桶进行称量，获得各时间点t下沥青混合料试件和样品桶质量为M3，则沥青混合料试件吸收的水气质量为M(t)=M3-M0，绘制沥青混合料吸收水气质量随时间的变化曲线，以17.17%相对湿度下的吸收水气质量随时间变化数据为例。
本发明还提供一种存储介质，该存储介质中存储有计算机程序，其中，该计算机程序被设置为运行时执行上述方法步骤。存储介质可以包括如软盘、光盘、DVD、硬盘、闪存、U盘、CF卡、SD卡、MMC卡、SM卡、记忆棒（MemoryStick）、XD卡等。
计算机软件产品存储在存储介质中，包括若干指令用以使得一台或多台计算机设备（可以是个人计算机设备、服务器或其他网络设备等）用以执行本发明方法的全部或部分步骤。
本发明还提供一种沥青混合料水气运动参数计算系统，该沥青混合料水气运动参数计算系统包括处理器和存储器,存储器上存储有计算机程序，计算机程序被处理器执行时，实现沥青混合料水气运动参数计算方法。
与现有技术相比，本发明所提供的沥青混合料水气运动参数计算方法、存储介质及系统具有以下有益效果：
通过基于获取的沥青混合料试件吸收水气质量随时间的变化数据，建立描述沥青混合料内部水气扩散和吸附的三维水气运动模型，从而更加准确的描述沥青混合料内水气运动特性的理论计算方法，从而为研究沥青混合料水损害机理提供理论基础。
以上所述本发明的具体实施方式，并不构成对本发明保护范围的限定。任何根据本发明的技术构思所做出的各种其他相应的改变与变形，均应包含在本发明的保护范围内。

Claims (7)

1.一种沥青混合料水气运动参数计算方法，其特征在于，包括步骤：
制备并成型圆柱体的待测沥青混合料试件；
对沥青混合料试件进行干燥处理，并置于磁悬浮重量平衡系统中开展积聚型水气运动试验，获取沥青混合料试件吸收水气质量随时间的变化数据；
基于获取的沥青混合料试件吸收水气质量随时间的变化数据，建立描述沥青混合料内部水气扩散和吸附的三维水气运动模型，具体包括：
通过移动态水分子满足菲克扩散定律得到三维水气扩散基本模型；建立三维水气扩散模型所述三维水气运动模型为：
Figure QLYQS_1
或者，
Figure QLYQS_2
，
其中，M(t)表示时刻t时沥青混合料样品内部积聚的水气质量；
Figure QLYQS_3
，表示在温度T和纯水气压P环境下沥青混合料样品能够容纳的最大水气质量；C0表示水气运动试验中沥青混合料样品表面恒定的水气浓度；γ表示扩散进入沥青混合料试件中的水气分子从移动态变成束缚态的概率，s-1；β 表示同一时刻水气分子从束缚态变成移动态的概率，s-1；/>
Figure QLYQS_4
表示第一类零阶贝塞尔函数的第m个零点；D表示移动态水气扩散系数；t表示扩散时间；a为沥青混合料试件的半径，H为沥青混合料试件的高度；
确定初始条件和边界条件；
进行拉普拉斯变换和逆变换；
其中，所述初始条件为：
Figure QLYQS_5
，/>
Figure QLYQS_6
；
Figure QLYQS_7
，/>
Figure QLYQS_8
；
所述边界条件为：
Figure QLYQS_9
，/>
Figure QLYQS_10
；
Figure QLYQS_11
，/>
Figure QLYQS_12
；
Figure QLYQS_13
，/>
Figure QLYQS_14
；
Figure QLYQS_15
为扩散时间t时，沥青混合料试件侧表面浓度，/>
Figure QLYQS_16
为扩散时间t时，沥青混合料试件下表面的水气浓度，/>
Figure QLYQS_17
为扩散时间t时，沥青混合料试件上表面的水气浓度，/>
Figure QLYQS_18
为水气运动试验中沥青混合料试件周边的恒定水气浓度；
进行拉普拉斯变换和逆变换，包括：
基于简化后的扩散模型
Figure QLYQS_19
对初始条件进行转换，则初始条件变化成：
Figure QLYQS_20
，/>
Figure QLYQS_21
；
Figure QLYQS_22
，/>
Figure QLYQS_23
；
同理简化边界条件成为：
Figure QLYQS_24
；
Figure QLYQS_25
；
Figure QLYQS_26
；
采用拉普拉斯变换对边界条件进行转换：
Figure QLYQS_27
，/>
Figure QLYQS_28
；
Figure QLYQS_29
，/>
Figure QLYQS_30
；
Figure QLYQS_31
，/>
Figure QLYQS_32
；
采用拉普拉斯逆变换，可得边界条件如下：
Figure QLYQS_33
；/>
Figure QLYQS_34
；
Figure QLYQS_35
，/>
Figure QLYQS_36
；
Figure QLYQS_37
，/>
Figure QLYQS_38
。
2.如权利要求1所述的一种沥青混合料水气运动参数计算方法，其特征在于，所述对沥青混合料试件进行干燥处理，并置于磁悬浮重量平衡系统中开展积聚型水气运动试验，获取沥青混合料试件吸收水气质量随时间的变化数据包括步骤：
进行试验准备；
进行试验；
对获取的数据进行处理。
3.如权利要求2所述的一种沥青混合料水气运动参数计算方法，其特征在于，所述进行试验包括步骤：
称取真空环境下样品桶的质量；
将沥青混合料试件放置在样品桶里，对反应腔进行抽真空处理，并不断称量样品桶和试件质量，直至质量保持恒定；
设置试验温度和水气压力。
4.如权利要求1所述的一种沥青混合料水气运动参数计算方法，其特征在于，所述建立三维水气扩散模型包括步骤：
分离变量；
施加边界条件和初始条件；
积分求解水气质量。
5.如权利要求2所述的一种沥青混合料水气运动参数计算方法，其特征在于：
所述对获取的数据进行处理指通过对试验数据进行处理，得到沥青混合料中的水气扩散、积聚参数。
6.一种存储介质，其特征在于：
所述存储介质中存储有计算机程序，其中，所述计算机程序被设置为运行时执行所述权利要求1-5中任一项中所述的沥青混合料水气运动参数计算方法。
7.一种沥青混合料水气运动参数计算系统，其特征在于：
所述沥青混合料水气运动参数计算系统包括处理器及存储器，所述存储器上存储有计算机程序，所述计算机程序被所述处理器执行，实现如权利要求1-5任一项所述的沥青混合料水气运动参数计算方法。
CN202011434387.3A 2020-12-10 2020-12-10 一种沥青混合料水气运动参数计算方法、存储介质及系统 Active CN112699529B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011434387.3A CN112699529B (zh) 2020-12-10 2020-12-10 一种沥青混合料水气运动参数计算方法、存储介质及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011434387.3A CN112699529B (zh) 2020-12-10 2020-12-10 一种沥青混合料水气运动参数计算方法、存储介质及系统

Publications (2)

Publication Number Publication Date
CN112699529A CN112699529A (zh) 2021-04-23
CN112699529B true CN112699529B (zh) 2023-06-23

Family

ID=75506974

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011434387.3A Active CN112699529B (zh) 2020-12-10 2020-12-10 一种沥青混合料水气运动参数计算方法、存储介质及系统

Country Status (1)

Country Link
CN (1) CN112699529B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114323831B (zh) * 2021-11-25 2023-09-26 河海大学 一种泡沫沥青混合料残存水分迁移路径的评价方法
CN114235640B (zh) * 2021-12-17 2024-04-19 武汉理工大学 沥青混合料内部水气扩散路径长度的计算方法及装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110455665A (zh) * 2019-08-16 2019-11-15 武汉理工大学 沥青混合料积聚型水气运动参数检测装置及方法
CN110595935A (zh) * 2019-09-20 2019-12-20 武汉理工大学 一种道路沥青混合料水气运动参数的确定方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105928814A (zh) * 2016-04-15 2016-09-07 武汉理工大学 沥青混合料穿透型水气扩散系数检测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110455665A (zh) * 2019-08-16 2019-11-15 武汉理工大学 沥青混合料积聚型水气运动参数检测装置及方法
CN110595935A (zh) * 2019-09-20 2019-12-20 武汉理工大学 一种道路沥青混合料水气运动参数的确定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
冻融循环对沥青混合料内水气扩散的影响;罗蓉;柳子尧;黄婷婷;涂崇志;冯光乐;;中国公路学报(第09期);24-30 *
基质吸力对非饱和离子型稀土抗剪强度的影响;洪本根;罗嗣海;胡世丽;王观石;姚康;;岩土力学(第06期);281-288 *

Also Published As

Publication number Publication date
CN112699529A (zh) 2021-04-23

Similar Documents

Publication Publication Date Title
CN112699529B (zh) 一种沥青混合料水气运动参数计算方法、存储介质及系统
Gelb et al. Characterization of porous glasses: Simulation models, adsorption isotherms, and the Brunauer− Emmett− Teller analysis method
Myers et al. Adsorption in porous materials at high pressure: theory and experiment
Uchida et al. Dissociation condition measurements of methane hydrate in confined small pores of porous glass
Kim et al. Adsorption equilibria of water vapor on alumina, zeolite 13X, and a zeolite X/activated carbon composite
Allardice et al. The-brown coal/water system: Part 2. Water sorption isotherms on bed-moist Yallourn brown coal
Clifford Properties of water in capillaries and thin films
Miyahara et al. Modeling capillary condensation in cylindrical nanopores: a molecular dynamics study
Awschalom et al. Liquid-film instabilities in confined geometries
Kishimoto et al. Diffusion coefficients for amorphous polymer and water systems
Nguyen et al. Adsorption of supercritical gases in porous media: determination of micropore size distribution
Scaife et al. Characterization of porous materials by gas adsorption: do different molecular probes give different pore structures?
Daiko et al. Proton conduction and pore structure in sol− gel glasses
CN110595935A (zh) 一种道路沥青混合料水气运动参数的确定方法
Jagiełło et al. Characterization of microporous carbons using adsorption at near ambient temperatures
Liu et al. A water retention curve model describing the effect of temperature
Tang et al. In situ study of diffusion and interaction of water and mono-or divalent anions in a positively charged membrane using two-dimensional correlation FT-IR/attenuated total reflection spectroscopy
Day The moisture potential of soils
Si et al. Gas-liquid competitive adsorption characteristics and coal wetting mechanism under different pre-adsorbed gas conditions
WO2022227385A1 (zh) 一种半刚性基层水气扩散系数的检测方法
Cai et al. Evaluation of methods for determining the pore size distribution and pore-network connectivity of porous carbons
Shimizu et al. Stochastic pore blocking and gating in PDMS–glass nanopores from vapor–liquid phase transitions
Wang et al. Experimental study on adsorption promotion and desorption restrain of gas-containing coal in low-temperature environment
Abdrashitov et al. Modification of the proton conducting membranes of MF-4SC with a carbon phase via gamma-initiated graft polymerization of vinyliden chloride
Yuan et al. Research on the Effect of Freeze–Thaw Cycles at Different Temperatures on the Pore Structure of Water-Saturated Coal Samples

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