[发明专利]动模量、阻尼比计算方法及装置

权利要求书

1.一种动模量计算方法，其特征在于，应用于计算设备，所述计算设备中存储有在动三轴试验过程中采集的应力数据点信息及应变数据点信息，所述方法包括：

根据相关函数及采集的应力数据点信息得到应力时程的第一自相关函数，其中，所述应力数据点信息包括时刻及对应的应力值，所述第一自相关函数为：

其中，R_σ(iΔτ)表示第iΔτ时刻应力时程的自相关函数值，σ_k表示第kΔτ时刻的应力值，σ_k+i表示第(k+i)Δτ时刻的应力值；

根据相关函数及采集的应变数据点信息得到应变时程的第二自相关函数，其中，所述应变数据点信息包括时刻及对应的应变值，所述第二自相关函数为：

其中，R_ε(iΔτ)表示第iΔτ时刻应变时程的自相关函数值，Δτ表示数据采集的时间间隔，N表示循环周数，N₀表示每周期采集的应力数据点信息、应变数据点信息的个数，ε_k表示第kΔτ时刻的应变值，ε_k+i表示第(k+i)Δτ时刻的应变值，应力时程及应变时程均为弦波形式；

通过所述第一自相关函数计算得到应力时程的平均幅值；

通过所述第二自相关函数计算得到应变时程的平均幅值；

根据应力时程的平均幅值、应变时程的平均幅值及动模量计算公式得到用于评估土体的动力稳定性的动模量，其中，动模量计算公式为：

E＝σ_a/ε_a

其中，E表示动模量，σ_a表示应力时程的平均幅值，ε_a表示应变时程的平均幅值。

2.根据权利要求1所述的方法，其特征在于，所述通过所述第一自相关函数计算得到应力时程的平均幅值，包括：

计算所述第一自相关函数的极大值及极小值；

将所述第一自相关函数中的多个极大值进行线性拟合，计算得到与所述第一自相关函数中的多个极大值对应的拟合直线的第一截距；

将所述第一自相关函数中的多个极小值进行线性拟合，计算得到与所述第一自相关函数中的多个极小值对应的拟合直线的第二截距；

根据所述第一截距、第二截距及应力时程的平均幅值计算公式计算出应力时程的平均幅值，其中，所述第一截距及第二截距均为纵截距，所述应力时程的平均幅值计算公式为：

其中，σ_a表示应力时程的平均幅值，R_σ,max(0)表示第一截距，R_σ,min(0)表示第二截距。

3.根据权利要求1所述的方法，其特征在于，所述通过所述第二自相关函数计算得到应变时程的平均幅值，包括：

计算所述第二自相关函数的极大值及极小值；

将所述第二自相关函数中的多个极大值进行线性拟合，计算得到与所述第二自相关函数中的多个极大值对应的拟合直线的第三截距；

将所述第二自相关函数中的多个极小值进行线性拟合，计算得到与所述第二自相关函数中的多个极小值对应的拟合直线的第四截距；

根据所述第三截距、第四截距及应变时程的平均幅值计算公式计算出应变时程的平均幅值，其中，所述第三截距及第四截距均为纵截距，所述应变时程的平均幅值计算公式为：

其中，ε_a表示应变时程的平均幅值，R_ε,max(0)表示第三截距，R_ε,min(0)表示第四截距。