详情

原文传递基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法

专利名称：	基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法
摘要：	基于最优松弛时间域的沥青混合料Prony级数模型参数的确定方法，本发明属于沥青混合料线性黏弹性能的评价方法，它为了解决目前缺乏有效、准确确定Prony级数模型参数的方法。本发明以GS模型为储能模量主曲线模型，以近似Kramers‑Kroning关系函数为损耗模量主曲线模型，以WLF方程为时间‑温度位移方程，得到参考温度下的主曲线，确定松弛强度与松弛时间之间的关系，按照松弛时间域的中心点的不同设置A、B两个试算组，计算含有不同项数的模型的误差，由试算组的特征项数和中心点计算得到特征松弛时间范围，从而确定的最优松弛时间域。本发明在最优松弛时间域内基函数充分发挥作用，得到准确的Prony级数表达式。
专利类型：	发明专利
国家地区组织代码：	黑龙江;23
申请人：	哈尔滨工业大学
发明人：	谭忆秋;吕慧杰;孙志棋;孟安鑫;曲元魁
专利状态：	有效
申请日期：	2019-02-27T00:00:00+0800
发布日期：	2019-06-14T00:00:00+0800
申请号：	CN201910155265.1
公开号：	CN109884288A
代理机构：	哈尔滨市松花江专利商标事务所
代理人：	岳泉清
分类号：	G01N33/42(2006.01);G;G01;G01N;G01N33
申请人地址：	150001 黑龙江省哈尔滨市南岗区西大直街92号
主权项：	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项的模型的误差。
所属类别：	发明专利