沥青混合料剪胀力学行为分析及数值模拟解析方案【附数据】

✨ 长期致力于沥青混合料、剪胀、局部三轴试验、三轴试验、X-ray CT扫描技术、图像处理、离散元模拟、本构模型、细观结构研究工作擅长数据搜集与处理、建模仿真、程序编写、仿真设计。✅ 专业定制毕设、代码✅如需沟通交流点击《获取方式》1沥青混合料剪胀行为的局部三轴试验与机理分析利用UTM-25多功能材料系统配备高精度环向变形传感器量程±5mm分辨率0.001mm进行沥青混合料AC-13的局部三轴压缩试验。试件直径100mm高150mm试验温度40℃围压50kPa、100kPa、200kPa加载速率1mm/min。在试件中部30mm高度范围内贴附三个LVDT测量径向变形轴向变形通过两个对称LVDT测量。结果表明在偏应力达到峰值后体积应变从压缩转为膨胀剪胀起始点对应轴向应变的0.5%-0.8%。围压越大剪胀越不明显围压200kPa时体积膨胀最大值为0.3%而50kPa时达到0.9%。剪胀角ψ随塑性应变增加而增长从初始5度增至20度。通过Reynolds剪胀理论分析剪胀率dε_v/dε_q与应力比q/p呈线性关系斜率即为剪胀系数。实验数据拟合得到剪胀系数β0.32。利用X-ray CT扫描剪胀前后的试件空隙率从初始4.2%增加到剪切带区域的8.5%证实了剪胀源于集料颗粒的重排和互锁作用释放。2基于离散元方法的三维细观数值模拟及剪胀机制采用PFC3D建立沥青混合料的离散元模型包含粗集料粒径2.36mm、沥青胶浆通过平行粘结模型模拟和空隙三相。集料颗粒采用随机多面体生成算法颗粒数约12000个胶浆用半径为0.3mm的小颗粒填充并赋予平行粘结法向刚度1e9 N/m切向刚度4e8 N/m粘结强度法向3MPa切向2MPa。采用Burgers接触模型模拟沥青的粘弹性Maxwell粘度η1.2e5 Pa·sKelvin刚度5e7 Pa粘度6e4 Pa·s。在三轴压缩模拟中围压100kPa轴向应变达到2%时体积应变开始转为正膨胀与室内试验一致。通过追踪颗粒位移场发现剪胀主要发生在局部剪切带内带内颗粒的法向接触力减小40%配位数从4.2降到2.8。模拟得到剪胀角ψ18度与试验值21度接近。空隙率沿剪切带分布的不均匀性导致剪胀行为呈现应变软化特征。3考虑剪胀的粘弹塑性本构模型与参数标定基于Drucker-Prager屈服准则和非关联流动法则构建一个包含剪胀效应的粘弹塑性本构模型。模型中屈服函数 F α I1 sqrt(J2) - k其中α和k与内摩擦角φ和粘聚力c相关。剪胀势函数采用与屈服面不同的斜率剪胀角ψ作为状态变量随等效塑性应变演化。粘弹性部分采用Burgers模型塑性部分采用Chaboche各向同性硬化。通过三轴蠕变和循环加载试验标定模型参数Burgers参数E11200MPaη11500MPa·sE2800MPaη2600MPa·s塑性参数φ38度c0.35MPa初始ψ5度最终ψ22度。将模型嵌入ABAQUS UMAT对沥青路面车辙进行模拟计算结果与现场实测车辙深度加载10万次后12mm误差小于1.5mm。考虑剪胀的模型预测车辙比忽略剪胀的模型深8%说明剪胀导致材料软化加速车辙发展。该本构模型为沥青路面设计提供了更准确的力学依据。import numpy as np from scipy.optimize import curve_fit def local_triaxial_dilatancy(axial_strain, volumetric_strain): # 剪胀角计算根据Reynolds公式 # 剪胀率 dε_v / dε_q d_eps_v np.gradient(volumetric_strain, axial_strain) d_eps_q np.gradient(axial_strain - volumetric_strain/3, axial_strain) # 偏应变 psi np.arcsin(d_eps_v / (2 * d_eps_q - d_eps_v)) return np.degrees(psi) def dem_contact_model(particle_radius, kn1e9, ks4e8): # 离散元线性接触模型参数PFC3D风格 # 返回接触刚度和阻尼系数 return dict(knkn, ksks, friction0.5, damping0.7) def viscous_plastic_model(stress, strain, params): # 粘弹塑性本构模型简化一维 E1, eta1, E2, eta2, phi, c, psi0, psi_inf params # Burgers 粘弹性应变计算 eps_viscous stress/E1 stress/eta1 * time stress/E2 * (1 - np.exp(-E2/eta2 * time)) # 塑性屈服判断Drucker-Prager I1 stress[0] stress[1] stress[2] J2 0.5 * ((stress[0]-stress[1])**2 (stress[1]-stress[2])**2 (stress[2]-stress[0])**2) 3*stress[3]**2 alpha 2 * np.sin(np.radians(phi)) / (np.sqrt(3)*(3-np.sin(np.radians(phi)))) k 6 * c * np.cos(np.radians(phi)) / (np.sqrt(3)*(3-np.sin(np.radians(phi)))) F alpha * I1 np.sqrt(J2) - k if F 0: # 塑性流动非关联 psi psi0 (psi_inf - psi0) * (1 - np.exp(-eps_plastic)) # 简化塑性应变增量 deps_plastic 0.001 eps_plastic deps_plastic return eps_viscous eps_plastic # 示例数据模拟 axial_strain_exp np.linspace(0, 0.03, 30) vol_strain_exp -0.001 * axial_strain_exp 0.003 * axial_strain_exp**2 psi_angle local_triaxial_dilatancy(axial_strain_exp, vol_strain_exp) print(剪胀角 (度) 峰值:, np.max(psi_angle))