[发明专利]一种钢桥面铺装力学控制指标动态计算方法

权利要求书

1.一种钢桥面铺装力学控制指标动态计算方法，其特征在于，包括如下步骤：

(1)获取不同温度和不同荷载频率下的沥青混合料动态模量数据；

(2)根据时温等效原理采用遗传基因算法对试验结果分析拟合，得出沥青混合料的动态模量主曲线公式及不同温度间的移位因子计算公式；

(3)在有限元软件ABAQUS中，利用Fortran语言对沥青混合料的动态模量主曲线公式及不同温度间的移位因子计算公式进行二次开发，将沥青混合料的动态模量主曲线公式及不同温度间的移位因子计算公式作为ABAQUS软件的用户子程序；

(4)当用ABAQUS软件求解力学控制指标时，调用动态模量用户子程序，输入温度与荷载频率相关参数，得到相对应的沥青混合料的动态模量，为后续的模型计算提供准确的模量值。

2.如权利要求1所述的钢桥面铺装力学控制指标动态计算方法，其特征在于，步骤(1)中的试验温度分别为10℃、20℃和30℃，荷载频率分别为0.01Hz、0.02Hz、0.05Hz、0.1Hz、0.2Hz、0.5Hz、1Hz、2Hz、5Hz、10Hz和15Hz。

3.如权利要求1所述的钢桥面铺装力学控制指标动态计算方法，其特征在于，步骤(2)中的动态模量主曲线公式对应于20℃参考温度下为

$\ln E=\ln E_0+A*\left\{1-e^{\left(\frac{-ln\mathrm{\ω}x}{\mathrm{\β}}\right)}*\[1+\frac{ln\mathrm{\ω}x}{\mathrm{\β}}+\frac{{\left(ln\mathrm{\ω}x\right)}^2}{2{\mathrm{\β}}^2}\]\right\}$

其中，ln E₀＝6.788，A＝2.999，β＝3.651，lnωx＝ln(freq)+9.882；

不同温度下沥青混合料动态模量对应于20℃的移位因子计算公式为

${\mathrm{ln\α}}_T=-19.732*(1-e^{-\frac{T-20}{44.257}})$

其中，α_T为移位因子，T为试验温度。

4.如权利要求1所述的钢桥面铺装力学控制指标动态计算方法，其特征在于，步骤(2)中的遗传基因算法包括如下步骤：

(1)定义回归函数形式和参数，回归函数为其中RSS(residual sum of squares)为差值平方和，y_i为试验测出的动态模量，为拟合出的动态模量；

(2)生成N组基因起始序列，N为基因组数且为偶数；在定义过程中，最初的p个参数(t1,t2,t3,…,tp)被挑选作为一组起始基因序列，t1，t2...，tp为不同温度下材料的模量值，并确定每个参数的数值取值范围，而后，通过均匀分布取值方式对每个参数在特定数值范围内选取此参数数值，组成由p个参数数值组成的基因组Λ_i；其中为水平移位因子，Λ_i代表第i组的基因组；最后，由N个基因组形成的基因组序列{Λ₁,Λ₂,…,Λ_N}被确定：

基因1：

基因2：

…

基因N：

(3)利用回归函数∏(Λ_i)对每组基因序列的适合度进行计算，此处回归函数采用最小二乘法函数形式；

(4)按照计算出的适合度对原基因序列进行等级排序形成新序列：

(5)将排序后最靠近的两组序列进行两两配对，利用线性组合函数产生新的后代序列Θ_i：

${\mathrm{\Θ}}_1={\mathrm{\Φ}}^{\left(1\right)}{\mathrm{\Λ}}_1^{*}+\[1-{\mathrm{\Φ}}^{\left(1\right)}\]{\mathrm{\Λ}}_2^{*}$

${\mathrm{\Θ}}_2={\mathrm{\Φ}}^{\left(2\right)}{\mathrm{\Λ}}_1^{*}+\[1-{\mathrm{\Φ}}^{\left(2\right)}\]{\mathrm{\Λ}}_2^{*}$

${\mathrm{\Θ}}_{N-1}={\mathrm{\Φ}}^{(N-1)}{\mathrm{\Λ}}_{N-1}^{*}+\[1-{\mathrm{\Φ}}^{(N-1)}\]{\mathrm{\Λ}}_N^{*}$

${\mathrm{\Θ}}_N={\mathrm{\Φ}}^{\left(N\right)}{\mathrm{\Λ}}_{N-1}^{*}+\[1-{\mathrm{\Φ}}^{\left(N\right)}\]{\mathrm{\Λ}}_N^{*}$

其中，Θ_i为经过配对后的新后代序列；Φ⁽ⁱ⁾＝Aφ⁽ⁱ⁾，φ⁽ⁱ⁾为[0,1]区间中按均匀分布方式选取的一个随机数值，

(6)检验取值范围：其中为左约束条件，为右约束条件，如果变量值超出约束范围，将其取为相应的约束值；

(7)弃除排序后底部不良的M组基因；

(8)将余下N－M组优良基因序列再加上M组新基因序列组成新的基因序列，重复上述步骤进行计算至规定次数后结束。