CN109470553A - 应用Laplace变换实现沥青混合料材料参数快速转换的方法 - Google Patents
应用Laplace变换实现沥青混合料材料参数快速转换的方法 Download PDFInfo
- Publication number
- CN109470553A CN109470553A CN201811384847.9A CN201811384847A CN109470553A CN 109470553 A CN109470553 A CN 109470553A CN 201811384847 A CN201811384847 A CN 201811384847A CN 109470553 A CN109470553 A CN 109470553A
- Authority
- CN
- China
- Prior art keywords
- laplace
- asphalt
- equation
- relaxation modulus
- domain
- 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
- 238000000034 method Methods 0.000 title claims abstract description 77
- 239000010426 asphalt Substances 0.000 title claims abstract description 74
- 239000000203 mixture Substances 0.000 title claims abstract description 29
- 239000000463 material Substances 0.000 title claims abstract description 28
- 230000009466 transformation Effects 0.000 claims abstract description 24
- 238000006243 chemical reaction Methods 0.000 claims abstract description 12
- 239000013013 elastic material Substances 0.000 claims abstract description 8
- 238000012360 testing method Methods 0.000 claims description 11
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000003068 static effect Effects 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 abstract description 4
- 230000010354 integration Effects 0.000 description 6
- 235000008331 Pinus X rigitaeda Nutrition 0.000 description 4
- 235000011613 Pinus brutia Nutrition 0.000 description 4
- 241000018646 Pinus brutia Species 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- LFEUVBZXUFMACD-UHFFFAOYSA-H lead(2+);trioxido(oxo)-$l^{5}-arsane Chemical compound [Pb+2].[Pb+2].[Pb+2].[O-][As]([O-])([O-])=O.[O-][As]([O-])([O-])=O LFEUVBZXUFMACD-UHFFFAOYSA-H 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010835 comparative analysis Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 235000015220 hamburgers Nutrition 0.000 description 1
- 238000004154 testing of material Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0001—Type of application of the stress
- G01N2203/0003—Steady
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/0058—Kind of property studied
- G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
- G01N2203/0071—Creep
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/0202—Control of the test
- G01N2203/0212—Theories, calculations
- G01N2203/0218—Calculations based on experimental data
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
应用Laplace变换实现沥青混合料材料参数快速转换的方法,本发明属于沥青混合料材料参数换算领域,它为了解决现有计算沥青混合料的松弛模量的方法精度难以保证,算法复杂的问题。参数转换方法:一、计算沥青混合料的蠕变柔量,采用Prony级数模型对沥青混合料的蠕变柔量进行拟合;二、对沥青混合料Prony级数模型施加Laplace正变换,通过黏弹性材料中松弛模量与蠕变柔量的换算关系计算得到Laplace域松弛模量;三、采用FT法施加Laplace逆变换求解沥青混合料松弛模量时域解。本发明可快速、精确的获取松弛模量值,通过Laplace积分逆变换,可完成任意时刻下松弛模量的求解。
Description
技术领域
本发明属于沥青混合料材料参数换算领域,具体涉及一种基于Laplace变换实现沥青混合料材料参数快速转换的方法。
背景技术
沥青混合料的松弛模量是开展沥青路面结构力学计算的材料参数之一。然而,目前并没有直接测试沥青混合料松弛模量仪器,而是采用沥青混合料静态蠕变试验获取沥青混合料的蠕变柔量,然后利用黏弹性力学中松弛模量与蠕变柔量的卷积公式进行松弛模量的计算。该方法首先是将卷积公式离散成有限个无穷积分式,然后将区段中值时刻的模量值假定为该区段的模量值,即将模量假定为一个定值,这样的简化处理使得卷积公式变成一个简单的积分式,通过数值积分的方法获取每个区段的积分值,然后通过等式反算松弛模量值。但该方法存在如下不足:
(1)每一个区段的模量值假定为一个定值,这样的处理方法需要将区间划分的很小,才能保证计算的精度,而多小的区段才能保证计算精度并没有明确处理方法。
(2)该种方法只能计算有限个时间点下的松弛模量值,不能计算任意时刻的松弛模量值。
(3)该方法需要处理复杂的数值积分,需要编写复杂的计算程序,计算效率低。
发明内容
本发明的目的是为了解决现有计算沥青混合料的松弛模量的方法精度难以保证,算法复杂的问题,而提供一种应用Laplace变换实现沥青混合料材料参数快速转换的方法。
本发明应用Laplace变换实现沥青混合料材料参数快速转换的方法按下列步骤实现:
一、沥青混合料蠕变柔量表征
测试固定荷载下沥青混合料的蠕变变形量△l,然后采用方程(1)计算沥青混合料的蠕变柔量;
方程(1)中△l是加载前与加载结束后高度差,h是试件的高度,σ0是荷载大小;
采用Prony级数模型对沥青混合料的蠕变柔量进行拟合,Prony级数的表达式如方程(2)所示:
方程(2)中J0和Ji是Prony级数的参数(其中J0是t=0时的Prony参数,Ji通过Matlab数据拟合获得),τi为延迟时间,t是时间,然后曲线拟合得到沥青混合料Prony级数模型;
二、Laplace域松弛模量表征
对沥青混合料Prony级数模型施加Laplace正变换,Prony级数模型的Laplace正变换如方程(3)所示:
方程(3)中s是Laplace变量;
沥青混合料为(典型的)黏弹性材料,黏弹性材料中松弛模量与蠕变柔量的Laplace域换算关系如方程(4)所示:
方程(4)中E(s)是Laplace域的松弛模量;
根据方程(4)得到Laplace域的松弛模量计算公式如下:
方程(5)中L-1()是Laplace逆变换,E(t)是时域中的松弛模量;
将方程(3)代入到方程(5)中,解得:
即得到Laplace域松弛模量;
三、松弛模量时域值求解
采用FT法(Fixed Talbot算法)对方程(6)施加Laplace逆变换,完成沥青混合料松弛模量时域解的求解,从而完成沥青混合料材料的参数快速转换。
本发明先在室内开展沥青混合料蠕变试验,并采用Prony级数来表征沥青混合料的蠕变柔量获取模型参数,然后对Prony级数施加Laplace正变换;之后基于黏弹性力学当中松弛模量与蠕变柔量的Laplace转换公式,将蠕变柔量的Laplace表达式代入转换公式当中,完成Laplace域松弛模量表达式的求解;最后对表达式施加Laplace数值逆变换,完成松弛模量时域值的求解。
本发明应用Laplace变换实现沥青混合料材料参数快速转换的方法包含以下有益效果:
1、Laplace域松弛模量表达式的推导简单易行,概念清晰;
2、避免了复杂的数值积分运算及程序编写;
3、可快速、精确的获取松弛模量值;
4、通过Laplace积分逆变换,可完成任意时刻下松弛模量的求解。
附图说明
图1是实施例一中20℃下四种沥青混合料蠕变柔量,其中■代表type-A,●代表type-B,▲代表type-C,▼代表type-D;
图2是实施例一中广义Kelvin模型的结构示意图;
图3是实施例一中采用卷积方法和FT法计算得到的松弛模量值(R2=0.999),其中▼代表卷积方法,▲代表FT法;
图4是实施例一中M=5和M=8的沥青混合料松弛模量值(R2=1.00),其中▼代表M=5,▲代表M=8;
图5是实施例一中M=5和M=10的沥青混合料松弛模量值(R2=1.00),其中▼代表M=5,▲代表M=10;
图6是实施例一中M=5和M=15的沥青混合料松弛模量值(R2=1.00),其中▼代表M=5,▲代表M=15;
图7是本发明应用Laplace变换实现沥青混合料材料参数快速转换的方法的流程图。
具体实施方式
具体实施方式一:本实施方式应用Laplace变换实现沥青混合料材料参数快速转换的方法按下列步骤实施:
一、沥青混合料蠕变柔量表征
测试固定荷载下沥青混合料的蠕变变形量△l,然后采用方程(1)计算沥青混合料的蠕变柔量;
方程(1)中△l是加载前与加载结束后高度差,h是试件的高度,σ0是荷载大小;
采用Prony级数模型对沥青混合料的蠕变柔量进行拟合,Prony级数的表达式如方程(2)所示:
方程(2)中J0和Ji是Prony级数的参数(其中J0是t=0时的Prony参数,Ji通过Matlab数据拟合获得),τi为延迟时间,t是时间,然后曲线拟合得到沥青混合料Prony级数模型;
二、Laplace域松弛模量表征
对沥青混合料Prony级数模型施加Laplace正变换,Prony级数模型的Laplace正变换如方程(3)所示:
方程(3)中s是Laplace变量;
沥青混合料为(典型的)黏弹性材料,黏弹性材料中松弛模量与蠕变柔量的Laplace域换算关系如方程(4)所示:
方程(4)中E(s)是Laplace域的松弛模量;
根据方程(4)得到Laplace域的松弛模量计算公式如下:
方程(5)中L-1()是Laplace逆变换,E(t)是时域中的松弛模量;
将方程(3)代入到方程(5)中,解得:
即得到Laplace域松弛模量;
三、松弛模量时域值求解
采用FT法(Fixed Talbot算法)对方程(6)施加Laplace逆变换,完成沥青混合料松弛模量时域解的求解,从而完成沥青混合料材料的参数快速转换。
本实施方式步骤三中尝试采用Durbin法和Dubner法进行求解,其公式分别如下:
式中a,b,T——计算参数;
t——具体计算时间;
Re,Im——复数的实部和虚部。
在完成松弛模量时域值求解之前,需要选取计算效率高、精度高的Laplace数值逆变换方法,采用如下原函数形式,函数的参数如下表1。
表1 Prony级数系数
表2不同数值反演方法计算效率
本实施方式发现随着时间的增加,Durbin法和Dubner法反演的精度不断下降,而FT法能够较好满足反演的计算精度。从表2三者计算效率可以看出,FT法不仅计算精度要高于Durbin法和Dubner法,而且其计算效率非常高,逆变换50个数据点只需要0.04ms。FT在Prony级数Laplace反演中不仅保证了计算精度,而且能够很好的保证计算效率。FT法能够较好的应用到沥青混合料黏弹参数的换算当中。
FT法中的参数M是指将Laplace反演的路径划分的步数,这类似于数值积分方法中的积分步数。步长越小,步数越多,计算精度越高,但计算精度越低。Kim在研究层状黏弹体系理论中,提出FT法中任意的参数M值不会影响反演精度。但其研究的是短时间内荷载函数的Laplace逆变换,其研究结论不能应用于Prony级数的Laplace逆变换中。因此,研究参数M等于5、8、10、15时FT法对Prony级数反演的精度和效率。
表3不同M值下FT反演效率
从表3可以看出,无论M值选为多少,对FT法的计算精度和效率的影响都很小。所以,采用FT法研究沥青混合料黏弹参数换算。
因此,对方程(6)施加Laplace逆变换(FT法),即可完成松弛模量时域解的求解。
具体实施方式二:本实施方式与具体实施方式一不同的是步骤一中沥青混合料制成直径100mm,高度150mm的标准试件。
具体实施方式三:本实施方式与具体实施方式一不同的是步骤一中蠕变变形量的测试过程是在静态加载下,控制荷载大小为0.3MPa,加载时间为240s。
具体实施方式四:本实施方式与具体实施方式一至三之一不同的是步骤一中Prony级数的表达式中当i=1时,则τi=0.001;当i=2时,则τi=0.01;当i=3时,则τi=0.1;当i=4时,则τi=1;当i=5时,则τi=10;当i=6时,则τi=100;当i=7时,则τi=1000。
具体实施方式五:本实施方式与具体实施方式一至四之一不同的是步骤一中利用Maltab软件中曲线拟合得到沥青混合料Prony级数模型。
具体实施方式六:本实施方式与具体实施方式一至五之一不同的是步骤三中所述的FT法的算法公式如方程(9)所示,方程(6)中Laplace域松弛模量表达式E(s)为方程(9)中F()函数的表达式,即将F(r)和F(s(θk))中的自变量值r和s(θk)的值代入E(s)中的s变量,分别求取E(s)的值即为F(r)和F(s(θk))的值:
其中,θk=kπ/M,s(θ)=rθ(cotθ+i),σ(θ)=θ+(θcotθ-1)cotθ,r=2M/(5t);
式中t——计算时间,r,M——模型参数。
具体实施方式七:本实施方式与具体实施方式六不同的是M的取值为5~15。
具体实施方式八:本实施方式与具体实施方式七不同的是M的取值为5~8。
实施例一:本实施例应用Laplace变换实现沥青混合料材料参数快速转换的方法按下列步骤实施:
一、沥青混合料蠕变柔量表征
在实验室内成型4种直径为150mm,高度170mm的圆柱体沥青混合料试件,分别命名为type-A、type-B、type-C、type-D,然后钻芯切割成直径100mm,高度150mm的标准试件,采用多功能材料试验机(UTM-250)在静态加载下,荷载大小为0.3MPa,加载时间为240s,通过公式(1)计算沥青混合料的蠕变柔量,如图1所示;
方程(1)中△l是加载前与加载结束后高度差,h是试件的高度,σ0是荷载大小;
将蠕变柔的试验数据导入到matlab中,通过自带工具箱Curve fitting确定不同沥青混合料的Prony参数如下;
表4四种沥青混合料Prony级数参数值
二、Laplace域松弛模量表征
线黏弹性材料的应力-应变-时间关系,主要有微分和积分型两大类。恒定荷载下黏弹材料的力学行为通过时间相关的材料函数来描述,而对于交变荷载作用下的黏弹材料性能,采用频率相关的复函数来描述。该部分内容主要基于广义Kelvin模型微分型本构关系,介绍黏弹材料蠕变柔量、松弛模量及动态模量的计算过程;
经典的黏弹本构由弹簧和黏糊串/并联组合而成,由此组成不同的本构模型,如:Maxwell、Kelvin模型及Burgers模型等,广义Kelvin模型由n个Kelvin链组成,设Kelvin链的第i个单元的应变为εi,其弹簧弹性模量和阻尼器黏性系数分别为Ei和ηi,用微分算子表示应变:
由n个Kelvin单元组成的广义Kelvin模型的总应变:
将此展开经整理后,得到一般模型的本构方程:
写作
或
Pσ=Qε (22)
其中微分算子:
为求得式(22)的解,对公式施加Laplace积分变换,则得到代数方程:
或
为了求得本构方程的蠕变柔量函数,将固定荷载σ=σ0H(t)的Laplace变换式代入(24)中,求得蠕变柔量的像函数:
通过Laplace逆变换求得蠕变柔量的时域解:
同样,为了求得松弛模量的时域解,将应变函数的Laplace积分变换式代入(24)中,求得松弛模量的像空间函数如公式(27),然后通过Laplace逆变换完成时域求解如式(28)。
由(25)和(28)可求得蠕变柔量与松弛模量像空间解的数学关系:
通过Laplace逆变换,可求得蠕变柔量与松弛模量时域解的数学关系:
目前,针对沥青及沥青混合料黏弹参数的换算研究主要集中蠕变柔量和松弛模量的换算,首先通过试验测得沥青的蠕变柔量,然后将蠕变柔量在时域内离散,利用蠕变柔量与松弛模量的时域下的卷积关系如式(29),计算得到松弛模量。首先是将卷积公式离散成有限个积分公式如下式:
将松弛模量E(τ)看成区间[ti-1,ti]之间的中值,即方程(24)即简化为下式:
然后应用数值求积公式,如Gauss积分,求得每个区间的积分值,然后求解得到松弛模量值;
应用Laplace数值反演方法是首先得到蠕变柔量Prony级数的Laplace变换式,然后通过方程(4)完成Laplace域松弛模量表达式的求解,如下式:
其目标就是采用Laplace数值逆变换的方法,完成上式的求解;
三、松弛模量时域值求解
这里,采用type-A的蠕变柔量数据,对比分析Laplace方法(FT法)与卷积方法计算的精度和效率,卷积方法中的积分采用精度较高的Gauss积分法,计算效率采用计算总时长来评价,计算的精度采用平均误差来评价,平均误差的计算公式如下:
方程中,n是数据总量,yA、yB分别为两种方法的计算结果;
如图3所示,两种方法计算结果相似,通过计算分析,两种方法的计算结果平均误差误差小于10-5,在240个时间点的计算中,卷积方法的所用时长为0.204s,而FT方法的计算时长为0.008s,FT方法比卷积方法快近25倍,这为沥青路面层状黏弹体系力学长时间计算节约了大量的时间;
在FT方法中,参数M值会影响Laplace反演的精度和效率,普遍认为M值取5~15,这里将M的值取为5、8、10及15。如图4-6所示和表5所示,M值的大小对计算精度影响很小,对计算效率存在一定的影响,在计算过程可选取较小的M值。
表5不同M值下反算时长
Claims (8)
1.应用Laplace变换实现沥青混合料材料参数快速转换的方法,其特征在于该方法按下列步骤实现:
一、沥青混合料蠕变柔量表征
测试固定荷载下沥青混合料的蠕变变形量,然后采用方程(1)计算沥青混合料的蠕变柔量;
方程(1)中△l是加载前与加载结束后高度差,h是试件的高度,σ0是荷载大小;
采用Prony级数模型对沥青混合料的蠕变柔量进行拟合,Prony级数的表达式如方程(2)所示:
方程(2)中J0和Ji是Prony级数的参数,τi为延迟时间,t是时间,然后曲线拟合得到沥青混合料Prony级数模型;
二、Laplace域松弛模量表征
对沥青混合料Prony级数模型施加Laplace正变换,Prony级数模型的Laplace正变换如方程(3)所示:
方程(3)中s是Laplace变量;
沥青混合料为黏弹性材料,黏弹性材料中松弛模量与蠕变柔量的Laplace域换算关系如方程(4)所示:
方程(4)中E(s)是Laplace域的松弛模量;
根据方程(4)得到Laplace域的松弛模量计算公式如下:
方程(5)中L-1()是Laplace逆变换,E(t)是时域中的松弛模量;
将方程(3)代入到方程(5)中,解得:
即得到Laplace域松弛模量;
三、松弛模量时域值求解
采用FT法对方程(6)施加Laplace逆变换,完成沥青混合料松弛模量时域解的求解,从而完成沥青混合料材料的参数快速转换。
2.根据权利要求1所述的应用Laplace变换实现沥青混合料材料参数快速转换的方法,其特征在于步骤一中沥青混合料制成直径100mm,高度150mm的标准试件。
3.根据权利要求1所述的应用Laplace变换实现沥青混合料材料参数快速转换的方法,其特征在于步骤一中蠕变变形量的测试过程是在静态加载下,控制荷载大小为0.3MPa,加载时间为240s。
4.根据权利要求1所述的应用Laplace变换实现沥青混合料材料参数快速转换的方法,其特征在于步骤一中Prony级数的表达式中当i=1时,则τi=0.001;当i=2时,则τi=0.01;当i=3时,则τi=0.1;当i=4时,则τi=1;当i=5时,则τi=10;当i=6时,则τi=100;当i=7时,则τi=1000。
5.根据权利要求1所述的应用Laplace变换实现沥青混合料材料参数快速转换的方法,其特征在于步骤一中利用Maltab软件中曲线拟合得到沥青混合料Prony级数模型。
6.根据权利要求1所述的应用Laplace变换实现沥青混合料材料参数快速转换的方法,其特征在于步骤三中所述的FT法的算法公式如方程(9)所示,方程(6)中Laplace域松弛模量表达式E(s)为方程(9)中F()函数的表达式,即将F(r)和F(s(θk))中的自变量值r和s(θk)的值代入E(s)中的s变量,分别求取E(s)的值即为F(r)和F(s(θk))的值:
其中,θk=kπ/M,s(θ)=rθ(cotθ+i),σ(θ)=θ+(θcotθ-1)cotθ,r=2M/(5t);
式中t——计算时间,r,M——模型参数。
7.根据权利要求6所述的应用Laplace变换实现沥青混合料材料参数快速转换的方法,其特征在于M的取值为5~15。
8.根据权利要求7所述的应用Laplace变换实现沥青混合料材料参数快速转换的方法,其特征在于M的取值为5~8。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811384847.9A CN109470553B (zh) | 2018-11-20 | 2018-11-20 | 应用Laplace变换实现沥青混合料材料参数快速转换的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811384847.9A CN109470553B (zh) | 2018-11-20 | 2018-11-20 | 应用Laplace变换实现沥青混合料材料参数快速转换的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109470553A true CN109470553A (zh) | 2019-03-15 |
CN109470553B CN109470553B (zh) | 2021-01-12 |
Family
ID=65674008
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811384847.9A Expired - Fee Related CN109470553B (zh) | 2018-11-20 | 2018-11-20 | 应用Laplace变换实现沥青混合料材料参数快速转换的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109470553B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110031611A (zh) * | 2019-04-19 | 2019-07-19 | 哈尔滨工业大学 | 基于多种群遗传算法的分数阶粘弹性模型参数识别的方法 |
CN110658079A (zh) * | 2019-09-18 | 2020-01-07 | 浙江大学 | 沥青面层在多重环境梯度耦合下的室内表征方法 |
CN112881196A (zh) * | 2021-01-14 | 2021-06-01 | 清华大学深圳国际研究生院 | 一种玻璃材料高温粘弹性参数检测方法 |
CN112966430A (zh) * | 2021-01-13 | 2021-06-15 | 天津城建大学 | 一种考虑应力松弛的沥青面层累积温度收缩应力计算方法 |
CN113533133A (zh) * | 2021-08-24 | 2021-10-22 | 哈尔滨工业大学 | 一种基于机器学习的沥青动态黏弹特性预测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0348745A (ja) * | 1989-07-18 | 1991-03-01 | Shimizu Corp | コンクリート多目的試験装置 |
US5187987A (en) * | 1991-11-19 | 1993-02-23 | The Pennsylvania Research Corporation | Bending beam creep test device with piston having a gas bearing |
CN103105334A (zh) * | 2013-02-07 | 2013-05-15 | 湘潭大学 | 一种构建黏弹性材料力学性能主曲线的方法 |
CN105787167A (zh) * | 2016-02-24 | 2016-07-20 | 重庆大学 | 聚合物基复合材料有效应力松弛系数的预测方法 |
CN107132114A (zh) * | 2017-07-04 | 2017-09-05 | 长安大学 | 一种沥青类材料蠕变柔量参数测试方法 |
CN107480835A (zh) * | 2017-09-12 | 2017-12-15 | 南通大学 | 一种纤维沥青混凝土黏弹性预测模型的构建方法 |
-
2018
- 2018-11-20 CN CN201811384847.9A patent/CN109470553B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0348745A (ja) * | 1989-07-18 | 1991-03-01 | Shimizu Corp | コンクリート多目的試験装置 |
US5187987A (en) * | 1991-11-19 | 1993-02-23 | The Pennsylvania Research Corporation | Bending beam creep test device with piston having a gas bearing |
CN103105334A (zh) * | 2013-02-07 | 2013-05-15 | 湘潭大学 | 一种构建黏弹性材料力学性能主曲线的方法 |
CN105787167A (zh) * | 2016-02-24 | 2016-07-20 | 重庆大学 | 聚合物基复合材料有效应力松弛系数的预测方法 |
CN107132114A (zh) * | 2017-07-04 | 2017-09-05 | 长安大学 | 一种沥青类材料蠕变柔量参数测试方法 |
CN107480835A (zh) * | 2017-09-12 | 2017-12-15 | 南通大学 | 一种纤维沥青混凝土黏弹性预测模型的构建方法 |
Non-Patent Citations (2)
Title |
---|
S. CHEN ET AL.: "Implement the Laplace transform to convert viscoelastic functions of asphalt mixtures", 《CONSTRUCTION AND BUILDING MATERIALS》 * |
吕慧杰 等: "沥青混合料蠕变柔量转换松弛模量的新方法", 《公路交通科技》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110031611A (zh) * | 2019-04-19 | 2019-07-19 | 哈尔滨工业大学 | 基于多种群遗传算法的分数阶粘弹性模型参数识别的方法 |
CN110031611B (zh) * | 2019-04-19 | 2022-02-11 | 哈尔滨工业大学 | 基于多种群遗传算法的分数阶粘弹性模型参数识别的方法 |
CN110658079A (zh) * | 2019-09-18 | 2020-01-07 | 浙江大学 | 沥青面层在多重环境梯度耦合下的室内表征方法 |
CN110658079B (zh) * | 2019-09-18 | 2024-04-05 | 浙江大学 | 沥青面层在多重环境梯度耦合下的室内表征方法 |
CN112966430A (zh) * | 2021-01-13 | 2021-06-15 | 天津城建大学 | 一种考虑应力松弛的沥青面层累积温度收缩应力计算方法 |
CN112966430B (zh) * | 2021-01-13 | 2022-04-22 | 天津城建大学 | 一种考虑应力松弛的沥青面层累积温度收缩应力计算方法 |
CN112881196A (zh) * | 2021-01-14 | 2021-06-01 | 清华大学深圳国际研究生院 | 一种玻璃材料高温粘弹性参数检测方法 |
CN113533133A (zh) * | 2021-08-24 | 2021-10-22 | 哈尔滨工业大学 | 一种基于机器学习的沥青动态黏弹特性预测方法 |
CN113533133B (zh) * | 2021-08-24 | 2022-02-15 | 哈尔滨工业大学 | 一种基于机器学习的沥青动态黏弹特性预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109470553B (zh) | 2021-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109470553A (zh) | 应用Laplace变换实现沥青混合料材料参数快速转换的方法 | |
CN109408969B (zh) | 一种应用有限元软件辨识橡胶粘弹参数建立本构模型的方法 | |
CN106021928A (zh) | 一种综合应力加速试验方法 | |
CN110146279B (zh) | 一种基于向量证据推理的船舶轴系不平衡故障诊断方法 | |
CN109871625A (zh) | 基于Johnson变换的非高斯风压模拟方法 | |
Han et al. | Theory of control systems described by differential inclusions | |
CN106063132B (zh) | 嵌入式装置的芯片上模数转换器(adc)线性度测试 | |
CN103278387A (zh) | 一种联合有限元软件模拟与实验的柔性复合材料材料参数测定的方法 | |
CN104615840A (zh) | 一种数字仿真模型的修正方法及系统 | |
CN109508498A (zh) | 基于bp人工神经网络的橡胶减振器配方设计系统及方法 | |
CN107341332A (zh) | 一种车用电机驱动系统的评价指标权重的确定方法 | |
Ruiz | Methodologies of ai technologies in science and culture | |
CN104504197A (zh) | 阿基米德螺旋线平面螺纹测量模型与其偏心参量修正方法 | |
Tichy et al. | nlrpBENCH: a benchmark for natural language requirements processing | |
CN105787265A (zh) | 基于综合集成赋权法的原子自旋陀螺随机误差建模方法 | |
CN104008234B (zh) | 密集模态含阻尼结构模型修正方法 | |
CN112200252A (zh) | 基于概率盒全局灵敏度分析和主动子空间的联合降维方法 | |
CN115687854B (zh) | 高精度土样参数测量方法及其系统 | |
Özbek et al. | A novel interactive system identification and control toolbox dedicated to real-time identification and model reference adaptive control experiments | |
Ouyang et al. | A rank-based high-dimensional test for equality of mean vectors | |
JP2023161480A (ja) | 処理システム、処理装置、処理方法、処理プログラム | |
Yang et al. | Modelling concrete strength using support vector machines | |
CN110276045B (zh) | 解析装置 | |
CN105893698A (zh) | 一种求解结构工程的非概率可靠性指标的仿射算法 | |
Zhang et al. | Research on fault diagnosis of diesel engine based on PCA-RBF neural network |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210112 |