CN110954443B - 采用插板法计算沥青材料动态接触角的方法 - Google Patents

采用插板法计算沥青材料动态接触角的方法 Download PDF

Info

Publication number
CN110954443B
CN110954443B CN201911247367.2A CN201911247367A CN110954443B CN 110954443 B CN110954443 B CN 110954443B CN 201911247367 A CN201911247367 A CN 201911247367A CN 110954443 B CN110954443 B CN 110954443B
Authority
CN
China
Prior art keywords
slide
contact angle
equation
formula
liquid
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
CN201911247367.2A
Other languages
English (en)
Other versions
CN110954443A (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 CN201911247367.2A priority Critical patent/CN110954443B/zh
Publication of CN110954443A publication Critical patent/CN110954443A/zh
Application granted granted Critical
Publication of CN110954443B publication Critical patent/CN110954443B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N13/00Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种采用插板法计算沥青材料动态接触角的方法,包括以下步骤:S1、采用插板法测试沥青材料不同浸入深度h对应天平受力变化值ΔF;S2、采用线性方程拟合根据不同浸入深度h对应天平受力变化值ΔF数据,并求解模型参数;S3、根据黏滞力作用下动态接触角的求解方程计算动态接触角。本发明所采用的计算方法考虑了黏滞力对动态接触角的影响,使得计算结果更加准确,为准确获取材料的表面能参数奠定了基础。

Description

采用插板法计算沥青材料动态接触角的方法
技术领域
本发明涉及道路工程领域,具体涉及一种采用插板法计算沥青材料动态接触角的方法。
背景技术
沥青路面是我国公路与城市道路的主要结构之一。而沥青与集料的黏附性能,直接关系到沥青混合料的疲劳寿命、自愈合能力及水稳定性等路用性能。目前国际上常常采用表面自由能法确定沥青–集料界面的黏附结合能(简称黏附能),以此作为量化沥青与集料黏附性能的指标。确定黏附能指标的核心是准确地测定沥青和集料的表面能参数,一般认为物质表面自由能由3个分量组成:(a)非极性分量
Figure GDA0002694865450000011
(b)极性分量酸分量
Figure GDA0002694865450000012
和(c)极性碱分量
Figure GDA0002694865450000013
对于沥青这种表面能较低的固体物质,表面自由能分量一般通过测定液体在其表面形成的接触角θ,并根据Young-Dupre方程计算。
根据Young-Dupre方程可知,求解沥青3个表面自由能,需要至少3种已知表面自由能的试剂与沥青的接触角θ。目前,常用的沥青接触角测量方法有静滴法和插板法:静滴法通过图像分析技术计算得到沥青的接触角,但此方法不能克服液滴自身重力等对于接触角结果的影响,同时试验结果直接受到设备精度的影响;而插板法是另一种能够快速的表征沥青接触角的有效方法,此方法通过记录接触线浸入溶液前后力的变化求解接触角。因为插板法能有效减少静滴法测试过程中人为观察产生的误差,被广泛用于测量物质表面能参数。
但插板法相比于静滴法最大的区别在于引入了浸润速度这一影响因素。插板法测试过程是一个动态过程,而图像法测试过程是一个静态过程,所以一般将插板法测试结果称为动态接触角θD,而光学法测试接触角称为静态接触角θ0。当沥青玻片匀速浸入液体时,液体在沥青玻片表面不断施加黏滞力,黏滞力使得沥青与液体之间的接触角发生变化。然而在计算动态接触角时,很少有人考虑黏滞力对接触角的影响。
发明内容
本发明要解决的技术问题在于针对现有技术中插板法计算沥青动态接触角时为考虑黏滞力影响的局限性缺陷,提供一种能够准确测量沥青动态接触角的方法。
为实现上述目的,本发明提供的采用插板法计算沥青材料动态接触角的方法,其特征在于:包括以下步骤:
S1、采用插板法测试沥青材料不同浸入深度h对应天平受力变化值ΔF;
S2、采用线性方程拟合根据不同浸入深度h对应天平受力变化值ΔF数据,并求解模型参数;
S3、根据黏滞力作用下动态接触角的求解方程计算动态接触角:
计算接触角>90°时,<90°时,=90°时在黏滞力作用下的动态接触角大小,具体方法为:
1)接触角>90°时:
当玻片未浸入试剂时,玻片在Z方向上受力平衡;玻片受到拉力=玻片重力:
F=W (1);
当玻片匀速向下浸入试剂时,若接触角>90°,液体曲面形状低于水平液面;由于液面内外存在压力差使得玻片浸润深度为h=h′-Δh0;其中Δh0可以使用公式(2)计算:
Figure GDA0002694865450000031
此时,玻片在Z方向上受力平衡;
F(t)=W-FbLPtcosθD-Fv (3);
其中,F(t)表示玻片受到天平施加的拉力;Fb表示液体对玻片施加浮力;W表示玻片自身重力;γL表示液体表面张力;Fv表示液体对玻片施加黏滞力;
当玻片底部截面的长宽a、b已知时,液体对玻片的浮力可以表示为:
Figure GDA0002694865450000032
液体表面张力γL在Z方向上的分量可以采用公式(5)计算:
γLPtcosθ=2(a+b)γLcosθD (5);
液体在玻片表面作用一个黏滞力,黏滞力的计算公式如下所示:
Figure GDA0002694865450000033
其中,μ是测试试剂的粘度,
Figure GDA0002694865450000034
是x方向(水平)上的速度梯度;
因此,将公式(4)—(6)带入公式(3)中,可以得到:
Figure GDA0002694865450000035
所以,天平施加拉力的变化量ΔF(t)=F(t)-F可以表示为:
Figure GDA0002694865450000041
由于玻片匀速浸入试剂,此时可以认为
Figure GDA0002694865450000042
保持不变为一个常数;ΔF可以进一步化简为:
ΔF(t)=m+kh(t) (9);
其中m、k分别用公式(10)、公式(11)表示:
Figure GDA0002694865450000043
Figure GDA0002694865450000044
在公式(10)和公式(11)中所有参数均为定值,说明ΔF(t)和h(t)之间满足线性关系;将公式(11)带入公式(10),可进一步化简公式:
Figure GDA0002694865450000045
通过试验记录不同时刻ΔF(t)和h(t)的数据,并采用线性方程进行数据拟合,可以得到参数m和k;采用规化求解的方法可以计算得到θD
2)接触角<90°时:
当接触角<90°时,液体曲面高于水平液面,由于液面内外存在压力差使得玻片浸润深度为h=h′+Δh0
F(t)=W-FbLcosθ-Fv (13);
类似的可以得到:
Figure GDA0002694865450000051
由于玻片匀速浸入试剂,因此可以认为
Figure GDA0002694865450000052
保持不变为一个常数;ΔF可以进一步化简为:
ΔF(t)=m+kh(t) (15);
其中m、k分别用公式(16)、公式(17)表示:
Figure GDA0002694865450000053
Figure GDA0002694865450000054
在公式(16)和公式(17)中所有参数均为定值,说明ΔF(t)和h(t)之间满足线性关系;将公式(17)带入公式(16),可进一步化简公式:
Figure GDA0002694865450000055
通过试验记录不同时刻ΔF(t)和h(t)的数据,并采用线性方程进行数据拟合,可以得到参数m和k;此时,公式(18)中只含有一个未知数θD,采用规化求解的方法可以计算得到θD
3)接触角=90°时:
当接触=90°时,此时玻片浸润深度为h=h′,cosθD=0,同理可以得到:
Figure GDA0002694865450000061
本发明的优点及有益效果如下:
本发明的插板法计算沥青动态接触角方法,针对现有计算方法未考虑黏滞力作用的局限性,提出了改进的方法。该法基于沥青玻片浸入测试试剂前后力的平衡方程,提出动态接触角计算方法。为准确测量沥青接触角提供理论基础和计算方法。
附图说明
图1是本发明中接触角>90°时沥青玻片受力图。
图2是本发明中接触角<90°时沥青玻片受力图。
图3是本发明的流程图。
图4是本发明70#基质沥青与化学试剂实验数据图;
其中:(a)70#基质沥青与蒸馏水的结果;(b)70#基质沥青与甲酰胺的结果;(c)70# 基质沥青与乙二醇的结果;(d)70#基质沥青与丙三醇的结果。
图5是本发明SBS改性沥青与化学试剂实验数据图;
图中:(a)SBS改性沥青与蒸馏水的结果;(b)SBS改性沥青与甲酰胺的结果;(c)SBS 改性沥青与乙二醇的结果;(d)SBS改性沥青与丙三醇的结果。
具体实施方式
以下将结合附图及具体实施例对本发明作进一步详细说明。
首先制备70#基质沥青与SBS改性沥青玻片各50个。并使用游标卡尺测量沥青玻片底部截面的长宽a、b。在不同浸润速度下接触角,同一个速度进行3次平行试验取。沥青玻片底部界面的长度为24.50±0.5mm,宽度0.65±0.5mm,具体尺寸如表1和表2所示:
表1 70#基质沥青玻片底部截面尺寸
Figure GDA0002694865450000071
表2 SBS改性沥青玻片底部截面尺寸
Figure GDA0002694865450000081
通过恒温循环浴保证测试试剂的温度为20±0.5℃,然后将沥青玻片放入动态接触角仪中,记录沥青玻片在不同浸入深度h时天平受力变化值ΔF。
根据公式(8)(14)(19)可知不论接触角是>90°、<90°、=90°,ΔF(t)和h(t)之间均满足线性关系。因此,采用线性拟合上述试验数据,得到直线方程的截距m和斜率k。求解结果如下表4所示:
表3 70#基质沥青试验结果
Figure GDA0002694865450000091
表4 SBS改性沥青沥青试验结果
Figure GDA0002694865450000101
上表中所有方程的拟合度R2均>0.99,然后将拟合得到的参数m和k,带入公式(8)(14)(19)中,再根据EXCEL中规划求解功能,计算得到黏滞力作用下的动态接触角θD。计算结果如下表5和表6所示:
表5不同速率测得70#沥青接触角
Figure GDA0002694865450000111
表6不同速率测得SBS改性沥青接触角
Figure GDA0002694865450000112
当动态接触角>90°时,考虑浸润速度的动态接触角>不考虑浸润速度的接触角;当动态接触角<90°时,考虑浸润速度的动态接触角<不考虑浸润速度的接触角。
接触角在90°附近时,浸润速度对接触角的影响较小,一般在1°以内;但是接触角>90°或者<90°时,浸润速度对接触角的影响较大,接触角的变化在5°左右。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。

Claims (1)

1.一种采用插板法计算沥青材料动态接触角的方法,其特征在于:包括以下步骤:
S1、采用插板法测试沥青材料不同浸入深度h对应天平受力变化值ΔF;
S2、采用线性方程拟合根据不同浸入深度h对应天平受力变化值ΔF数据,并求解模型参数;
S3、根据黏滞力作用下动态接触角的求解方程计算动态接触角:
计算接触角>90°时,<90°时,=90°时在黏滞力作用下的动态接触角大小,具体方法为:
1)接触角>90°时:
当玻片未浸入试剂时,玻片在Z方向上受力平衡;玻片受到拉力=玻片重力:
F=W (1);
当玻片匀速向下浸入试剂时,若接触角>90°,液体曲面形状低于水平液面;由于液面内外存在压力差使得玻片浸润深度为h=h(t)-h0;其中h0可以使用公式(2)计算:
Figure FDA0002694865440000011
此时,玻片在Z方向上受力平衡;
F(t)=W-FbLPtcosθD-Fv (3);
其中,F(t)表示玻片受到天平施加的拉力;Fb表示液体对玻片施加浮力;W表示玻片自身重力;γL表示液体表面张力;Fv表示液体对玻片施加黏滞力;
当玻片底部截面的长宽a、b已知时,液体对玻片的浮力可以表示为:
Figure FDA0002694865440000012
液体表面张力γL在Z方向上的分量可以采用公式(5)计算:
γLPtcosθ=2(a+b)γLcosθD (5);
液体在玻片表面作用一个黏滞力,黏滞力的计算公式如下所示:
Figure FDA0002694865440000021
其中,μ是测试试剂的粘度,
Figure FDA0002694865440000022
是x方向即水平上的速度梯度;
因此,将公式(4)—(6)带入公式(3)中,可以得到:
Figure FDA0002694865440000023
所以,天平施加拉力的变化量ΔF(t)=F(t)-F可以表示为:
Figure FDA0002694865440000024
由于玻片匀速浸入试剂,此时可以认为
Figure FDA0002694865440000025
保持不变为一个常数;ΔF(t)可以进一步化简为:
ΔF(t)=m+kh(t) (9);
其中m、k分别用公式(10)、公式(11)表示:
Figure FDA0002694865440000026
Figure FDA0002694865440000027
在公式(10)和公式(11)中所有参数均为定值,说明ΔF(t)和h(t)之间满足线性关系;将公式(11)带入公式(10),可进一步化简公式:
Figure FDA0002694865440000031
通过试验记录不同时刻ΔF(t)和h(t)的数据,并采用线性方程进行数据拟合,可以得到参数m和k;采用规化求解的方法可以计算得到θD
2)接触角<90°时:
当接触角<90°时,液体曲面高于水平液面,由于液面内外存在压力差使得玻片浸润深度为h=h(t)+h0
F(t)=W-FbLPtcosθD-Fv (13);
类似的可以得到:
Figure FDA0002694865440000032
由于玻片匀速浸入试剂,因此可以认为
Figure FDA0002694865440000035
保持不变为一个常数;ΔF(t)可以进一步化简为:
ΔF(t)=m+kh(t) (15);
其中m、k分别用公式(16)、公式(17)表示:
Figure FDA0002694865440000033
Figure FDA0002694865440000034
在公式(16)和公式(17)中所有参数均为定值,说明ΔF(t)和h(t)之间满足线性关系;将公式(17)带入公式(16),可进一步化简公式:
Figure FDA0002694865440000041
通过试验记录不同时刻ΔF(t)和h(t)的数据,并采用线性方程进行数据拟合,可以得到参数m和k;此时,公式(18)中只含有一个未知数θD,采用规化求解的方法可以计算得到θD
3)接触角=90°时:
当接触角 =90°时,此时玻片浸润深度为h=h(t),cosθD=0,同理可以得到:
Figure FDA0002694865440000042
CN201911247367.2A 2019-12-09 2019-12-09 采用插板法计算沥青材料动态接触角的方法 Active CN110954443B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911247367.2A CN110954443B (zh) 2019-12-09 2019-12-09 采用插板法计算沥青材料动态接触角的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911247367.2A CN110954443B (zh) 2019-12-09 2019-12-09 采用插板法计算沥青材料动态接触角的方法

Publications (2)

Publication Number Publication Date
CN110954443A CN110954443A (zh) 2020-04-03
CN110954443B true CN110954443B (zh) 2020-11-17

Family

ID=69980273

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911247367.2A Active CN110954443B (zh) 2019-12-09 2019-12-09 采用插板法计算沥青材料动态接触角的方法

Country Status (1)

Country Link
CN (1) CN110954443B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112816366B (zh) * 2021-01-20 2024-09-06 武汉理工大学 一种沥青材料表面能测试中化学试剂的选取方法及系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2782502B2 (ja) * 1995-09-25 1998-08-06 工業技術院長 生体分子吸着層の動的界面張力測定法とこれに用いる測定装置
CN103558124A (zh) * 2013-11-01 2014-02-05 同济大学 一种沥青材料表面能的测试方法
CN104406892A (zh) * 2014-12-04 2015-03-11 上海梭伦信息科技有限公司 一种基于分析天平的界面张力和接触角测试装置和方法
CN109374479B (zh) * 2018-10-10 2021-03-30 太原理工大学 基于力学测试确定变温条件下固体材料表面能的方法

Also Published As

Publication number Publication date
CN110954443A (zh) 2020-04-03

Similar Documents

Publication Publication Date Title
Schindler et al. The evaporation method: Extending the measurement range of soil hydraulic properties using the air‐entry pressure of the ceramic cup
CN101477012B (zh) 一种新型的蜡封测定沥青混合料毛体积密度的方法
CN111398572A (zh) 一种评价沥青混合料沥青用量的试验方法
CN112816366B (zh) 一种沥青材料表面能测试中化学试剂的选取方法及系统
CN110954443B (zh) 采用插板法计算沥青材料动态接触角的方法
Ranieri et al. Measurement of hydraulic conductivity in porous mixes
CN110082250A (zh) 一种快速检测非黏性土含水率的方法
CN105527190A (zh) 一种定压力的土壤水分特征曲线测定装置及测定方法
CN109374478B (zh) 基于力学测试确定变温条件下固定材料表面能测试方法
Martins Reis et al. Determination of the soil-water retention curve and the hydraulic conductivity function using a small centrifuge
CN104502227A (zh) 一种纤维密度的测定方法
CN105319233A (zh) 一种液体浸润法测试推进剂装药绝热层材料体膨胀系数的方法
CN110736675A (zh) 一种砂石骨料饱和面干吸水率测定装置及方法
CN105938084A (zh) 一种化学渗吸剂渗透性能评价方法
CN107860685B (zh) 一种砂质土含水率的检测方法
Tu et al. Effect of water vapor concentration on adhesion between asphalt and aggregate
Bhattacharjee et al. An alternative approach for the determination of bulk specific gravity and permeability of hot mix asphalt (HMA)
CN205333459U (zh) 一种定压力的土壤水分特征曲线测定装置
Rice Maximum Specific Gravity of Bituminous Mixtures by Vacuum Saturation Procedure
CN111625984A (zh) 一种基于有限元的沥青路面水损伤数值模拟方法
CN116067822A (zh) 一种测定混凝土粗骨料饱和面干吸水率的方法
CN113834840B (zh) 一种测试岩心渗吸效率的方法
Tu et al. Influence of infiltration velocity on the measurement of the surface energy components of asphalt binders using the Wilhelmy plate method
CN1540311A (zh) 一种压实沥青混合料的密度测定方法
Valenza et al. Measuring Permeability of Rigid Materials by a Beam‐Bending Method: V, Isotropic Rectangular Plates of Cement Paste

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