CN109884288A - 基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法 - Google Patents
基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法 Download PDFInfo
- Publication number
- CN109884288A CN109884288A CN201910155265.1A CN201910155265A CN109884288A CN 109884288 A CN109884288 A CN 109884288A CN 201910155265 A CN201910155265 A CN 201910155265A CN 109884288 A CN109884288 A CN 109884288A
- Authority
- CN
- China
- Prior art keywords
- equation
- model
- slack time
- time domain
- modulus
- 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
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法,本发明属于沥青混合料线性黏弹性能的评价方法,它为了解决目前缺乏有效、准确确定Prony级数模型参数的方法。本发明以GS模型为储能模量主曲线模型,以近似Kramers‑Kroning关系函数为损耗模量主曲线模型,以WLF方程为时间‑温度位移方程,得到参考温度下的主曲线,确定松弛强度与松弛时间之间的关系,按照松弛时间域的中心点的不同设置A、B两个试算组,计算含有不同项数的模型的误差,由试算组的特征项数和中心点计算得到特征松弛时间范围,从而确定的最优松弛时间域。本发明在最优松弛时间域内基函数充分发挥作用,得到准确的Prony级数表达式。
Description
技术领域
本发明属于沥青混合料线性粘弹性能的评价方法,具体涉及一种沥青混合料Prony级数模型参数的确定方法。
背景技术
沥青混合料是一种由沥青、集料和外掺剂等按照一定比例组合而成的复合型粘弹性材料,在小应变条件下,可采用线性粘弹性理论来表征沥青混合料的力学行为。沥青混合料线性粘弹性力学参数被广泛应用于沥青混合料的损伤识别和性能评价,准确测定线粘弹性力学参数对沥青路面的研究具有非常重要的意义。广义Maxwell模型是一种物理模型,被广泛应用于沥青混合料线性粘弹性力学行为的表征,具有计算高效且模拟准确的特点。广义Maxwell模型所对应的数学表达式为Prony级数模型。Prony级数模型中参数的准确确定是利用该模型准确表征沥青混合料线性粘弹性力学行为的基础。
目前,用于确定Prony级数模型参数的方法可以分为两类,一类是直接法,另一类是连续松弛时间谱法。直接法是直接通过数据拟合的方法确定模型参数,该方法操作方便,但是存在以下几点不足:(1)模型参数的解不唯一,(2)模型参数存在负值。连续松弛时间谱法是指利用连续松弛时间谱与Prony级数参数之间的关系确定模型参数,该方法所确定的模型参数均为正值,有效解决了直接法中出现的问题。但是在利用连续松弛时间谱法确定模型参数的过程中,需要预先设置一系列的松弛时间点,这些松弛时间点的选择直接影响到Prony级数模型的优劣性。目前,松弛时间点的设置主要是基于研究者的经验判断,缺乏统一的方法理论,求解结果的精度依赖于操作者的主观判断。针对目前确定Prony级数模型参数方法的不足,本发明提出一种基于最优松弛时间域的沥青混合料Prony级数模型粘弹参数的确定方法。
发明内容
本发明的目的是为了解决目前缺乏有效、准确确定Prony级数模型参数的方法,而提出一种基于最优松弛时间域的沥青混合料Prony级数模型粘弹参数的确定方法,用于研究沥青混合料的线性粘弹性行为。
本发明基于最优松弛时间域的沥青混合料Prony级数模型粘弹参数的确定方法按照以下步骤实现:
一、对沥青混合料试件进行动态模量试验,得到各温度和加载频率下的储能模量E'和损耗模量E″,然后以加载频率为横坐标,以存储模量为纵坐标绘制储能模量曲线,再以加载频率为横坐标,以损耗模量为纵坐标绘制损耗模量曲线;
二、以方程(1)广义西格摩德模型(GS模型)为储能模量主曲线模型,以方程(2)基于近似Kramers-Kroning关系确定的函数为损耗模量主曲线模型,以方程(3)WLF方程为时间—温度位移方程,选取参考温度,以方程(4)为拟合目标误差方程,采用EXCEL规划求解功能模块确定储能模量主曲线模型、损耗模量主曲线模型和WLF方程的参数,从而得到参考温度下的储能模量主曲线和损耗模量主曲线,并确定模型参数;
式中:δ=储能模量最小值的对数值,MPa;α=储能模量最大值与最小值之差的对数值,MPa;β,γ,λ=储能模量主曲线形状参数;αT=时温等效因子;T=试验温度,℃;Tr=参考温度,℃;C1=(第一)模型参数;C2=(第二)模型参数,℃;N=测量数据点的个数;E′CS,k=第k个利用GS模型计算得到的储能模量,MPa;E′m,k=第k个试验得到的储能模量,MPa;E″CS,k=第k个利用GS模型计算得到的损耗模量,MPa;E″m,k=第k个试验得到的损耗模量,MPa;
三、利用方程(5)所示的储能模量与连续松弛时间谱之间的关系确定松弛时间谱函数表达式(6),将步骤二中确定的模型参数代入连续松弛时间谱函数中得到连续松弛时间谱表达式,并绘制连续松弛时间谱曲线图;
其中:
式中:H(τ)=连续松弛时间谱,MPa;Im=函数的虚部;
四、利用连续松弛时间谱与离散松弛时间谱之间满足的方程(8)和方程(9)近似关系,从而确定松弛强度与松弛时间之间如方程(10)的关系;
Ei=H(τi)·△lnτi=H(τi)·(ln10)·△lgτi(10)
式中:Ei=松弛强度,MPa;τi=松弛时间,s;△lgτi=以10为底的对数刻度上松弛时间之间的距离;
五、对连续松弛时间谱函数求导数,确定连续松弛时间谱曲线峰值点对应的松弛时间点τpeak,计算方法如式(11)所示;
六、按照松弛时间域的中心点的不同设置A、B两个试算组,试算组A的中心点表达式如方程(12-1)所示,试算组B的中心点表达式如方程(12-2)所示:
其中,=试算组A的中心点;=试算组B的中心点;
模型误差和阈值误差的计算如下:
模型误差通过如式(13)所示的方程计算得到:
其中:
式中:N=测量数据点的个数;E′CP,k=第k个利用Prony级数模型计算得到的储能模量,如方程(14-1)所示,MPa;E′m,k=第k个试验得到的储能模量,MPa;E″CP,k=第k个利用Prony级数模型计算得到的损耗模量,如方程(14-2)所示,MPa;E″m,k=第k个试验得到的损耗模量,MPa;n=Prony级数模型中所含项数值;
在每个试算组中,利用方程(13)、方程(14-1)和方程(14-2)计算含有不同项数的模型的误差,当模型误差不大于阈值误差时,计算结束,得到的模型为特征模型,特征模型中所含项数值为特征项数,记为nC;
七、由试算组的特征项数和中心点计算得到特征松弛时间范围,计算方法如下:
其中:
八、计算步骤七中两个典型松弛时间范围的交集,采用方程(17)确定最优松弛时间域;
九、基于确定的最优松弛时间域,利用方程(6)和方程(10)确定Prony级数中的参数。
本发明利用最优松弛时间域法所确定的Prony级数模型,既保证了模型准确性,又简化了模型应用时的计算量。在确定Prony级数模型参数时,若所设定的松弛时间域未能覆盖最优松弛时间域,所确定的Prony级数模型误差偏大。
本发明作为一种沥青混合料线性粘弹性物理模型参数的确定方法,对于研究沥青混合料的粘弹性力学行为具有十分重要的意义。Prony级数模型是一种物理模型,具有明确的物理意义,被广泛应用于路面力学分析。该模型表达式为一系列指数基函数的和的形式,适于进行与线性粘弹性本构方程相关的遗传积分计算,是大多数商用有限元软件(如ANSYA和ABAQUS)中用于力学分析的线性粘弹性模型。另外,在路面结构的损伤分析中,Prony级数模型通常作为粘弹性本构方程用于计算材料的损伤程度,并用于路面使用寿命的预测。Prony级数模型广泛应用的前提是Prony级数模型参数的准确确定。
本发明提出了一种准确确定最优松弛时间域的方法,并利用最优松弛时间域确定了储能模量和损耗模量的Prony级数模型。利用该方法能够准确确定最优松弛时间域的原理为:储能模量和损耗模量Prony级数模型的基函数分别为和当ωk一定时,Prony级数模型为这些基函数之和,当基函数值非常小时,基函数的值对Prony级数贡献非常小,可以忽略不计,将这样的基函数舍弃可以简化计算。因此可以根据实际工程和研究需要,确定最优松弛时间域。在最优松弛时间域内基函数充分发挥作用,得到准确的Prony级数表达式。
附图说明
图1是实施例步骤一得到的储能模量曲线图,其中◆代表5℃,■代表20℃,▲代表35℃,●代表55℃;
图2是实施例步骤一得到的损耗模量曲线图,其中◆代表5℃,■代表20℃,▲代表35℃,●代表55℃;
图3是实施例步骤二得到的参考温度为20℃的储能模量主曲线图,其中◆代表5℃,■代表20℃,▲代表35℃,●代表55℃;
图4是实施例步骤二得到的参考温度为20℃的损耗模量主曲线图,其中◆代表5℃,■代表20℃,▲代表35℃,●代表55℃;
图5是实施例步骤三得到的连续松弛时间谱曲线图;
图6是实施例步骤六得到的模型误差随项数增加的变化情况图,其中◆代表试算组A,■代表试算组B;
图7是实施例步骤八得到的最优松弛时间域图。
具体实施方式
具体实施方式一:本实施方式基于最优松弛时间域的沥青混合料Prony级数模型粘弹参数的确定方法按照以下步骤实施:
一、对沥青混合料试件进行动态模量试验,得到各温度和加载频率下的储能模量E'和损耗模量E″,然后以加载频率为横坐标,以储能模量为纵坐标绘制储能模量曲线,再以加载频率为横坐标,以损耗模量为纵坐标绘制损耗模量曲线;
二、以方程(1)广义西格摩德模型(GS模型)为储能模量主曲线模型,以方程(2)基于近似Kramers-Kroning关系确定的函数为损耗模量主曲线模型,以方程(3)WLF方程为时间—温度位移方程,选取参考温度,以方程(4)为拟合目标误差方程,采用EXCEL规划求解功能模块确定储能模量主曲线模型、损耗模量主曲线模型和WLF方程的参数,从而得到参考温度下的储能模量主曲线和损耗模量主曲线;
式中:δ=储能模量最小值的对数值,MPa;α=储能模量最大值与最小值之差的对数值,MPa;β,γ,λ=储能模量主曲线形状参数;αT=时温等效因子;T=试验温度,℃;Tr=参考温度,℃;C1=模型参数;C2=模型参数,℃;N=测量数据点的个数;E′CS,k=第k个利用GS模型计算得到的储能模量,MPa;E′m,k=第k个试验得到的储能模量,MPa;E″CS,k=第k个利用GS模型计算得到的损耗模量,MPa;E″m,k=第k个试验得到的损耗模量,MPa;
三、利用方程(5)所示的储能模量与连续松弛时间谱之间的关系确定松弛时间谱函数表达式(6),将步骤二中确定的模型参数代入连续松弛时间谱函数中得到连续松弛时间谱表达式,并绘制连续松弛时间谱曲线图;
其中:
式中:H(τ)=连续松弛时间谱,MPa;Im=函数的虚部;
四、利用连续松弛时间谱与离散松弛时间谱之间满足的方程(8)和方程(9)近似关系,从而确定松弛强度与松弛时间之间如方程(10)的关系;
Ei=H(τi)·△lnτi=H(τi)·(ln10)·△lgτi(10)
式中:Ei=松弛强度,MPa;τi=松弛时间,s;△lgτi=以10为底的对数刻度上松弛时间之间的距离;
五、对连续松弛时间谱函数求导数,确定连续松弛时间谱曲线峰值点对应的松弛时间点τpeak,计算方法如式(11)所示;
六、按照松弛时间域的中心点的不同设置A、B两个试算组,试算组A的中心点表达式如方程(12-1)所示,试算组B的中心点表达式如方程(12-2)所示:
其中,=试算组A的中心点;=试算组B的中心点;
模型误差和阈值误差的计算如下:
模型误差通过如式(13)所示的方程计算得到:
其中:
式中:N=测量数据点的个数;E′CP,k=第k个利用Prony级数模型计算得到的储能模量,如方程(14-1)所示,MPa;E′m,k=第k个试验得到的储能模量,MPa;E″CP,k=第k个利用Prony级数模型计算得到的损耗模量,如方程(14-2)所示,MPa;E″m,k=第k个试验得到的损耗模量,MPa;n=Prony级数模型中所含项数值;
在每个试算组中,利用方程(13)、方程(14-1)和方程(14-2)计算含有不同项数的模型的误差,当模型误差不大于阈值误差时,计算结束,得到的模型为特征模型,特征模型中所含项数值为特征项数,记为nC;
七、由试算组的特征项数和中心点可以计算得到特征松弛时间范围,计算方法如下:
其中:
八、计算步骤七中两个典型松弛时间范围的交集,采用方程(17)确定最优松弛时间域;
九、基于确定的最优松弛时间域,利用方程(6)和方程(10)确定Prony级数中的参数。
具体实施方式二:本实施方式与具体实施方式一不同的是步骤一采用DTS多功能力学试验仪对对沥青混合料试件进行动态模量试验。
具体实施方式三:本实施方式与具体实施方式一或二不同的是步骤一动态模量试验中设置的试验温度为5℃~55℃。
具体实施方式四:本实施方式与具体实施方式三不同的是步骤一动态模量试验中设置的加载频率为0.1Hz~25Hz。
具体实施方式五:本实施方式与具体实施方式一至四之一不同的是步骤一沥青混合料试件为圆柱体。
具体实施方式六:本实施方式与具体实施方式一至五之一不同的是步骤三中选取20℃作为参考温度。
具体实施方式七:本实施方式与具体实施方式一至六之一不同的是步骤四中△lgτi=0.5。
具体实施方式八:本实施方式与具体实施方式一至七之一不同的是步骤六中方程(12-1)和方程(12-2)中JA=1,JB=20~40。
本实施方式JB的取值需要保证试算组B是有效试算组,即保证试算组的特征项数小于100项,建议取值为20~40之间。
具体实施方式九:本实施方式与具体实施方式八不同的是步骤六中方程(12-1)和方程(12-2)中JA=1,JB=31。
具体实施方式十:本实施方式与具体实施方式一至九之一不同的是步骤六中所述的阈值误差是指含有项数为100项的模型的误差。
本实施方式在每个试算组中,模型的误差值随着模型项数的增大而减小,当项数大于一定数量时,模型误差基本保持不变,该模型误差变为阈值误差,项数为阈值项数。为了实际应用,将含有项数为100项的模型的误差定义为阈值误差,特征项数小于100项的即为有效试算组。
实施例:本实施例基于最优松弛时间域的沥青混合料Prony级数模型粘弹参数的确定方法按照以下步骤实施:
一、采用DTS多功能力学试验仪对AC-13C型圆柱体沥青混合料试件(尺寸:直径100mm,高150mm)进行动态模量试验,试验在4种温度(5℃,20℃,35℃和55℃)和6个频率(0.1Hz,0.5Hz,1Hz,5Hz,10Hz,25Hz)下进行,得到各温度和加载频率下的储能模量E'和损耗模量E″,然后以加载频率为横坐标,以储能模量为纵坐标绘制储能模量曲线,如图1所示;再以加载频率为横坐标,以损耗模量为纵坐标绘制损耗模量曲线,如图2所示;
二、以方程(1)广义西格摩德模型(GS模型)为储能模量主曲线模型,以方程(2)基于近似Kramers-Kroning关系确定的函数为损耗模量主曲线模型,以方程(3)WLF方程为时间—温度位移方程,选取参考温度,以方程(4)为拟合目标误差方程,采用EXCEL规划求解功能模块确定储能模量主曲线模型、损耗模量主曲线模型和WLF方程的参数,从而得到参考温度20℃下的储能模量主曲线和损耗模量主曲线;
式中:δ=储能模量最小值的对数值,MPa;α=储能模量最大值与最小值之差的对数值,MPa;β,γ,λ=储能模量主曲线形状参数;αT=时温等效因子;T=试验温度,℃;Tr=参考温度,℃;C1=模型参数;C2=模型参数,℃;N=测量数据点的个数;E′CS,k=第k个利用GS模型计算得到的储能模量,MPa;E′m,k=第k个试验得到的储能模量,MPa;E″CS,k=第k个利用GS模型计算得到的损耗模量,MPa;E″m,k=第k个试验得到的损耗模量,MPa;
表1 参数拟合结果
参数 | δ | α | β | γ | λ | C<sub>1</sub> | C<sub>2</sub> |
拟合值 | 2.02 | 2.55 | –0.98 | –0.46 | 0.35 | 20.58 | 168.59 |
三、利用方程(5)所示的储能模量与连续松弛时间谱之间的关系确定松弛时间谱函数表达式(6),将步骤二中确定的参数代入连续松弛时间谱函数中得到连续松弛时间谱表达式,连续松弛时间谱曲线图如图5所示;;
其中:
式中:H(τ)=连续松弛时间谱,MPa;Im=函数的虚部;
四、利用连续松弛时间谱与离散松弛时间谱之间满足的方程(8)和方程(9)近似关系,从而确定松弛强度与松弛时间之间如方程(10)的关系;
Ei=H(τi)·△lnτi=H(τi)·(ln10)·△lgτi(10)
式中:Ei=松弛强度,MPa;τi=松弛时间,s;△lgτi=以10为底的对数刻度上松弛时间之间的距离,通常取0.5;
五、对连续松弛时间谱函数求导数,确定连续松弛时间谱曲线峰值点对应的松弛时间点τpeak,计算方法如式(11)所示;
利用方程(11)计算得到沥青混合料的τpeak,AC-13型沥青混合料的τpeak为10-2.54s;
六、按照松弛时间域的中心点的不同设置A、B两个试算组,试算组A的中心点表达式如方程(12-1)所示,试算组B的中心点表达式如方程(12-2)所示:
其中,=试算组A的中心点;=试算组B的中心点,JA=1;JB=31;
模型误差和阈值误差的计算如下:
模型误差通过如式(13)所示的方程计算得到:
其中:
式中:N=测量数据点的个数;E′CP,k=第k个利用Prony级数模型计算得到的储能模量,如方程(14-1)所示,MPa;E′m,k=第k个试验得到的储能模量,MPa;E″CP,k=第k个利用Prony级数模型计算得到的损耗模量,如方程(14-2)所示,MPa;E″m,k=第k个试验得到的损耗模量,MPa;n=Prony级数模型中所含项数值;Ei=(ln10)H(τi)△lgτi;
计算发现,在每个试算组中,模型的误差值随着模型项数的增大而减小,当项数大于一定数量时,模型误差基本保持不变,该模型误差变为阈值误差,项数为阈值项数,将含有项数为100项的模型的误差定义为阈值误差,特征项数小于100项的即为有效试算组;
在每个试算组中,模型项数n由1到100依次递增取值,利用方程(13)、方程(14-1)和方程(14-2)计算含有不同项数的模型的误差,当模型误差不大于阈值误差时,计算结束,得到的模型为特征模型,特征模型中所含项数值为特征项数,记为nC;
七、由试算组的特征项数和中心点可以计算得到特征松弛时间范围,计算方法如下:
其中:
A试算组和B试算组的典型松弛时间域计算结果如下式(17)和式(18)所示,图7为两个试算组的典型松弛时间范围;
八、计算步骤七中两个典型松弛时间范围的交集,采用方程(19)确定最优松弛时间域;
AC-13型沥青混合料最终确定的最优松弛时间域lgτi∈[-5.54,8.96];
九、基于确定的最优松弛时间域,利用方程(6)和方程(10)确定Prony级数中的参数。
本实施例确定Prony级数中的参数,AC-13型沥青混合料的Prony级数参数如表2所示。
表2 拟合参数列表
Claims (10)
1.基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法,其特征在于该Prony级数模型参数的确定方法按下列步骤实现:
一、对沥青混合料试件进行动态模量试验,得到各温度和加载频率下的储能模量E'和损耗模量E″,然后以加载频率为横坐标,以储能模量为纵坐标绘制储能模量曲线,再以加载频率为横坐标,以损耗模量为纵坐标绘制损耗模量曲线;
二、以方程(1)广义西格摩德模型为储能模量主曲线模型,以方程(2)基于近似Kramers-Kronig关系确定的函数为损耗模量主曲线模型,以方程(3)WLF方程为时间—温度位移方程,选取参考温度,以方程(4)为拟合目标误差方程,采用EXCEL规划求解功能模块确定储能模量主曲线模型、损耗模量主曲线模型和WLF方程的参数,从而得到参考温度下的储能模量主曲线和损耗模量主曲线,并确定模型参数;
式中:δ=储能模量最小值的对数值,MPa;α=储能模量最大值与最小值之差的对数值,MPa;β,γ,λ=储能模量主曲线形状参数;αT=时温等效因子;T=试验温度,℃;Tr=参考温度,℃;C1=模型参数;C2=模型参数,℃;N=测量数据点的个数;E′CS,k=第k个利用GS模型计算得到的储能模量,MPa;E′m,k=第k个试验得到的储能模量,MPa;E″CS,k=第k个利用GS模型计算得到的损耗模量,MPa;E″m,k=第k个试验得到的损耗模量,MPa;
三、利用方程(5)所示的储能模量与连续松弛时间谱之间的关系确定松弛时间谱函数表达式(6),将步骤二中确定的模型参数代入连续松弛时间谱函数中得到连续松弛时间谱表达式,并绘制连续松弛时间谱曲线图;
其中:
A=ln10δ;B=αln10;C=β-γlg2π;D=γ/ln10;
式中:H(τ)=连续松弛时间谱,MPa;Im=函数的虚部;
四、利用连续松弛时间谱与离散松弛时间谱之间满足的方程(8)和方程(9)近似关系,从而确定松弛强度与松弛时间之间如方程(10)的关系;
Ei=H(τi)·△lnτi=H(τi)·(ln10)·△lgτi(10)
式中:Ei=松弛强度,MPa;τi=松弛时间,s;△lgτi=以10为底的对数刻度上松弛时间之间的距离;
五、对连续松弛时间谱函数求导数,确定连续松弛时间谱曲线峰值点对应的松弛时间点τpeak,计算方法如式(11)所示;
六、按照松弛时间域的中心点的不同设置A、B两个试算组,试算组A的中心点表达式如方程(12-1)所示,试算组B的中心点表达式如方程(12-2)所示:
其中,
模型误差和阈值误差的计算如下:
模型误差通过如式(13)所示的方程计算得到:
其中:
式中:N=测量数据点的个数;E′CP,k=第k个利用Prony级数模型计算得到的储能模量,如方程(14-1)所示,MPa;E′m,k=第k个试验得到的储能模量,MPa;E″CP,k=第k个利用Prony级数模型计算得到的损耗模量,如方程(14-2)所示,MPa;E″m,k=第k个试验得到的损耗模量,MPa;n=Prony级数模型中所含项数值;
在每个试算组中,利用方程(13)、方程(14-1)和方程(14-2)计算含有不同项数的模型的误差,当模型误差不大于阈值误差时,计算结束,得到的模型为特征模型,特征模型中所含项数值为特征项数,记为nC;
七、由试算组的特征项数和中心点计算得到特征松弛时间范围,计算方法如下:
其中:
八、计算步骤七中两个典型松弛时间范围的交集,采用方程(17)确定最优松弛时间域;
九、基于确定的最优松弛时间域,利用方程(6)和方程(10)确定Prony级数中的参数。
2.根据权利要求1所述的基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法,其特征在于步骤一采用DTS多功能力学试验仪对对沥青混合料试件进行动态模量试验。
3.根据权利要求1所述的基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法,其特征在于步骤一动态模量试验中设置的试验温度为5℃~55℃。
4.根据权利要求1所述的基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法,其特征在于步骤一动态模量试验中设置的加载频率为0.1Hz~25Hz。
5.根据权利要求1所述的基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法,其特征在于步骤一沥青混合料试件为圆柱体。
6.根据权利要求1所述的基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法,其特征在于步骤三中选取20℃作为参考温度。
7.根据权利要求1所述的基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法,其特征在于步骤四中△lgτi=0.5。
8.根据权利要求1所述的基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法,其特征在于步骤六中方程(12-1)和方程(12-2)中JA=1,JB=20~40。
9.根据权利要求1所述的基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法,其特征在于步骤六中方程(12-1)和方程(12-2)中JA=1,JB=31。
10.根据权利要求1所述的基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法,其特征在于步骤六中所述的阈值误差是指含有项数为100项的模型的误差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910155265.1A CN109884288B (zh) | 2019-02-27 | 2019-02-27 | 基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910155265.1A CN109884288B (zh) | 2019-02-27 | 2019-02-27 | 基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109884288A true CN109884288A (zh) | 2019-06-14 |
CN109884288B CN109884288B (zh) | 2021-05-07 |
Family
ID=66930233
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910155265.1A Active CN109884288B (zh) | 2019-02-27 | 2019-02-27 | 基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109884288B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110555268A (zh) * | 2019-09-02 | 2019-12-10 | 内蒙动力机械研究所 | 一种固体推进剂松弛模量主曲线试验数据的分析软件 |
CN110658079A (zh) * | 2019-09-18 | 2020-01-07 | 浙江大学 | 沥青面层在多重环境梯度耦合下的室内表征方法 |
CN112507287A (zh) * | 2020-11-10 | 2021-03-16 | 山东省交通科学研究院 | 聚合物改性沥青相容性设计方法 |
CN112966430A (zh) * | 2021-01-13 | 2021-06-15 | 天津城建大学 | 一种考虑应力松弛的沥青面层累积温度收缩应力计算方法 |
CN113533133A (zh) * | 2021-08-24 | 2021-10-22 | 哈尔滨工业大学 | 一种基于机器学习的沥青动态黏弹特性预测方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0466060A2 (en) * | 1990-07-10 | 1992-01-15 | The B.F. Goodrich Company | Method and apparatus for determining the fundamental viscoelastic properties of a material |
JPH10326271A (ja) * | 1997-05-23 | 1998-12-08 | Ube Ind Ltd | 溶融粘弾性推算装置及びプログラム記録媒体 |
CN107132114A (zh) * | 2017-07-04 | 2017-09-05 | 长安大学 | 一种沥青类材料蠕变柔量参数测试方法 |
CN107391801A (zh) * | 2017-06-23 | 2017-11-24 | 中国人民解放军国防科学技术大学 | 含脱湿的推进剂动态热粘弹性本构模型的构建和应用方法 |
CN108846143A (zh) * | 2018-03-26 | 2018-11-20 | 哈尔滨工业大学 | 一种基于综合数据优化准则的沥青黏弹模型参数辨识方法 |
-
2019
- 2019-02-27 CN CN201910155265.1A patent/CN109884288B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0466060A2 (en) * | 1990-07-10 | 1992-01-15 | The B.F. Goodrich Company | Method and apparatus for determining the fundamental viscoelastic properties of a material |
JPH10326271A (ja) * | 1997-05-23 | 1998-12-08 | Ube Ind Ltd | 溶融粘弾性推算装置及びプログラム記録媒体 |
CN107391801A (zh) * | 2017-06-23 | 2017-11-24 | 中国人民解放军国防科学技术大学 | 含脱湿的推进剂动态热粘弹性本构模型的构建和应用方法 |
CN107132114A (zh) * | 2017-07-04 | 2017-09-05 | 长安大学 | 一种沥青类材料蠕变柔量参数测试方法 |
CN108846143A (zh) * | 2018-03-26 | 2018-11-20 | 哈尔滨工业大学 | 一种基于综合数据优化准则的沥青黏弹模型参数辨识方法 |
Non-Patent Citations (1)
Title |
---|
陈静云等: "沥青混合料动态粘弹性行为分析的模拟方法", 《中国公路学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110555268A (zh) * | 2019-09-02 | 2019-12-10 | 内蒙动力机械研究所 | 一种固体推进剂松弛模量主曲线试验数据的分析软件 |
CN110658079A (zh) * | 2019-09-18 | 2020-01-07 | 浙江大学 | 沥青面层在多重环境梯度耦合下的室内表征方法 |
CN110658079B (zh) * | 2019-09-18 | 2024-04-05 | 浙江大学 | 沥青面层在多重环境梯度耦合下的室内表征方法 |
CN112507287A (zh) * | 2020-11-10 | 2021-03-16 | 山东省交通科学研究院 | 聚合物改性沥青相容性设计方法 |
CN112507287B (zh) * | 2020-11-10 | 2022-11-29 | 山东省交通科学研究院 | 聚合物改性沥青相容性设计方法 |
CN112966430A (zh) * | 2021-01-13 | 2021-06-15 | 天津城建大学 | 一种考虑应力松弛的沥青面层累积温度收缩应力计算方法 |
CN113533133A (zh) * | 2021-08-24 | 2021-10-22 | 哈尔滨工业大学 | 一种基于机器学习的沥青动态黏弹特性预测方法 |
CN113533133B (zh) * | 2021-08-24 | 2022-02-15 | 哈尔滨工业大学 | 一种基于机器学习的沥青动态黏弹特性预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109884288B (zh) | 2021-05-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109884288A (zh) | 基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法 | |
Singh et al. | Evaluation of predictive models for estimating dynamic modulus of hot-mix asphalt in Oklahoma | |
CN103630450A (zh) | 考虑疲劳-蠕变交互损伤作用的沥青混合料寿命预测方法 | |
CN109933925B (zh) | 一种金属板材的冲压成形性能预测方法 | |
CN101746995A (zh) | 基于最紧密状态的沥青混合料配合比设计方法 | |
CN101339150B (zh) | 基于介电谱技术测定汽油辛烷值的方法 | |
CN103246821A (zh) | 一种基于仿真的多应力小样本加速寿命试验方案设计优化方法 | |
CN107341332A (zh) | 一种车用电机驱动系统的评价指标权重的确定方法 | |
CN105486845B (zh) | 基于轴载谱的沥青混合料多级加载高温蠕变曲线分析方法 | |
CN103345544B (zh) | 采用逻辑回归方法预测有机化学品生物降解性 | |
CN114235890B (zh) | 一种基于耗散能变化的沥青时温敏感性测试与评价方法 | |
CN101839836B (zh) | 一种测定高粘度改性沥青零剪切粘度的简化方法及其应用 | |
Kassem et al. | Quantification of the inherent uncertainty in the relaxation modulus and creep compliance of asphalt mixes | |
CN105004745A (zh) | 一种由核磁共振谱预测原油粘度的方法 | |
CN110334473B (zh) | 一种基于模量衰变的沥青路面结构相对弯沉的计算方法 | |
CN104483654B (zh) | 一种智能电能表测量正负偏差的综合评价方法及其系统 | |
CN109342310A (zh) | 一种基于老化损伤的橡胶密封制品寿命预测方法 | |
CN103134764B (zh) | 由透射红外光谱预测原油实沸点蒸馏曲线的方法 | |
CN1727877A (zh) | 由近红外光谱测定汽油性质数据的方法 | |
CN105181526A (zh) | 上浆剂生产过程中粘度的检测方法 | |
CN107144683A (zh) | 一种基于西格摩德函数的沥青混合料主曲线的构建方法 | |
CN109858699B (zh) | 水质定量模拟方法、装置、电子设备及存储介质 | |
CN101858877A (zh) | 一种生物柴油混合燃料混合比检测方法 | |
Dick | Practical Rubber Rheology and Dynamic Properties | |
Lee et al. | Performance-based mix design of unmodified and lime-modified hot mix asphalt |
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 |