一种土石混合料的回弹模量确定方法与流程

本发明涉及路基工程
技术领域：
，尤其涉及一种土石混合料的回弹模量确定方法。
背景技术：
：目前，我国山区、丘陵地带公路修筑过程中产生的大量土石混合料，主要用于路基修筑，然而，由于土石混填路基土的组成复杂、颗粒差距大，材料不均匀等问题，导致其物理力学性质受含水量、土石性质影响外，还与其含石量、粒径分布、土石比例等密切相关，因此，土石混填路基的工程力学特性呈现复杂规律，易出现不均匀沉降等病害，因此，需要对土石混填路基土的力学性能进行有效预测。目前，现有的《公路土工试验规程》(jtge40-2007)规定采用击实法、振动台法或表面振动压实法来成型小于60mm路基土的力学试件，由此测试试件的回弹模量等力学参数。然而，土石混填路基土中含有大量的卵砾石或碎石，卵石粒径为60～200mm、砾石粒径为2～60mm；粒径大于40mm或60mm卵砾石的加入，大大改变了土石混合料的力学强度。显然，现有方法并能测试含有卵砾石的土石混合料的力学强度，导致室内成型试件的工程性质与现场芯样强度差异较大，无法揭示材料组成结构与性能之间客观规律，导致土石混填路基的沉降变形与失稳特性没有得到有效预测。技术实现要素：为了解决上述问题，本发明的目的是提出一种土石混合料的回弹模量确定方法，能够有效预测路基土中含粒径大于40mm的土石混合料的回弹模量，为土石混填路基的设计提供精确参数，提高了土石混填路基设计的可靠性；同时，对提升土石混填路基的强度与承载力，降低土石混填路基的病害发生率，延长道路的使用寿命都具有重要意义。为了达到上述目的，本发明采用以下技术方案予以解决。一种土石混合料的回弹模量确定方法，包括以下步骤：步骤1：取部分土石混合料作为试样，将试样分为a1组和a2组；其中，a1组的土石混合料的粒径为0～40mm，a2组的土石混合料的粒径大于40mm；包含以下子步骤：1.1，将a1组的土石混合料分为粒径范围不同的n组，分别为a11,…a1i,…,a1n，i＝1,2,…n，n≥2；其中，a1i的粒径范围为di-1～dimm，d0＝0，dn＝40mm；1.2，将a2组的土石混合料配制成粒径范围不同的m组土石混合料，分别为a21,…a2j,…,a2m，j＝1,2,…m，m≥3；其中，a2j的粒径范围为d40～djmm，40＜d1＜…＜dm，d40＝40mm；步骤2：确定a1组中的p组土石混合料的实测回弹模量ex，x＝1,2,….,p，p≥4；包含以下子步骤：2.1，将a11至a1n这n组土石混合料进行混合，分别混合p组；对于p组土石混合料中的第一组土石混合料进行取样，烘干至恒重，重复多次取样并烘干，测定每份试样的含水率和干密度；2.2，绘制干密度-含水率的关系曲线，通过干密度-含水率的关系曲线，得到第一组土石混合料的最大实测干密度ρ1,max，最大实测干密度ρ1,max对应的含水率为最佳含水率ω1；2.3，按照步骤2.2确定的最大实测干密度ρmax和最佳含水率ω，并结合预设压实度k，采用振动击实法，将第一组土石混合料成型为圆柱体试件；施加竖向荷载pi，读取每级荷载下的回弹变形li，得到第一组土石混合料的实测回弹模量其中，d是圆柱体试件的直径，μ0是土石混合料的泊松比；pi每级加载的荷载，li每级荷载对应的回弹变形量；2.4，按照步骤2.1至步骤2.3，测试第二组至第p组土石混合料的实测回弹模量e2,….,ep，并结合第一组土石混合料的实测回弹模量，得到p组土石混合料的实测回弹模量ex。步骤3：确定土石混合料的最佳微力学参数；步骤4：确定a2组与a1组混合后的土石混合料的数值试验回弹模量e2j,s；包含以下子步骤：4.1，将a21组土石混合料与a1组中的p组土石混合料分别进行混合，根据土石混合料的最佳微力学参数对混合后的p组土石混合料分别进行数值模拟试验，得到混合后土石混合料的数值试验回弹模量e21,s，4.2，将a22,…a2j,…,a2m组土石混合料分别与a1组中的p组土石混合料进行混合，对混合后的[(m-1)×p]组土石混合料分别进行数值模拟试验，得到混合后的[(m-1)×p]组土石混合料的数值试验回弹模量，分别为e22,s,…,e2j,s,...,e2m,s，结合e21,s，得到a2组与a1组混合后的土石混合料的数值试验回弹模量e2j,s，其中，步骤5，根据a1组土石混合料的实测回弹模量ex、a2组与a1组混合后的土石混合料的数值试验回弹模量e2j,s，通过第一预设公式：确定回归系数a和b；第一预设公式中，ex为a1组土石混合料中第x组的实测回弹模量，单位为mpa；e2j,s为a2组与a1组混合后的土石混合料的数值试验回弹模量，单位为mpa；ρx,max为ex为a1组土石混合料中第x组的实测最大干密度，单位为g/cm3；λ为a2组土石混合料的质量占比；为a1组土石混合料的质量，单位为g；为a2组土石混合料的质量，单位为g；η＝d2max/40，d2max为a2组土石混合料的最大粒径，单位为mm，η为颗粒特征参数；a，b分别为回归系数；步骤6：根据回归系数a和b，确定土石混合料的回弹模量回归方程，通过土石混合料的回弹模量回归方程，确定待测土石混合料的回弹模量e。本发明技术方案的特点和进一步的改进为：(1)子步骤2.1包含以下子步骤：2.1.1，对于每份试样，以水质量与土石混合料试样的质量比为2％～3％的比例加水，搅拌均匀，密封，得润湿试样；2.1.2，将润湿试样装入试模，采用振动击实仪进行振动击实，脱模，得脱模试件，测试脱模试件的质量m、脱模试件的高度h和脱模试件的直径d；2.1.3，对脱模试件进行破碎，得破碎料，称取多份相同质量的破碎料，烘干至恒重，分别测试多份破碎料的含水率，计算多份破碎料的含水率平均值，即为试样的含水率ω；2.1.4，根据脱模试件的质量m、脱模试件的高度h、脱模试件的直径d和试样的含水率ω，计算出试样的干密度(2)步骤3包含以下子步骤：3.1，对p组土石混合料中每组土石混合料分别进行回弹模量的数值模拟试验，分别得到每组土石混合料的数值试验回弹模量ex,s，x＝1,2,….,p；3.2，根据第二预设公式计算实测回弹模量ex与数值试验回弹模量ex,s的累计相对偏差δ；反复调整土石混合料的微力学参数，得到多个δ，当δ最小时对应的微力学参数为最佳微力学参数；第二预设公式中，ex是p组土石混合料中第x组土石混合料的实测回弹模量，ex,s是p组土石混合料中第x组土石混合料的数值试验回弹模量，δ是ex与ex,s的相对偏差累加值。(3)步骤3.1包含以下子步骤：3.1.1、生成贯入杆模型、试筒模型和底座模型；根据p组土石混合料的颗粒组成、几何形状及物理参数，分别生成p组土石混合料的颗粒模型；采用重力作用，分别对p组土石混合料的颗粒模型分别进行初步密实，得p个初始虚拟试件；采用带阻尼的线弹性模型，设定土石混合料的微力学参数的初始值，分别对p个初始虚拟试件进行振动击实，即得p个虚拟试件；3.1.2、对p个虚拟试件分别逐级加荷载，每级加载的荷载为pi，读取每级荷载下的回弹变形li，得到p组土石混合料的数值试验回弹模量为：其中，ex,s为p组土石混合料中第x组土石混合料的数值试验回弹模量；d为p个虚拟试件中每个试件的直径；μ0为土石混合料的泊松比；pi为每级竖向加载的荷载，li为每级荷载对应的回弹变形量。进一步地，子步骤3.1.1中，所述物理参数为p组土石混合料中每组土石混合料的密度、模量和泊松比。进一步地，子步骤3.2中，土石混合料的微力学参数包含：土颗粒的湿润摩擦系数、土颗粒的阻尼系数、石颗粒的湿润摩擦系数、石颗粒的阻尼系数。更进一步地，所述土颗粒的湿润摩擦系数的取值范围为：0.20～0.46，土颗粒的阻尼系数取值范围为0.30～0.40，石颗粒的湿润摩擦系数取值范围为0.42～0.46，石颗粒的阻尼系数取值范围为0.65～0.75。(4)步骤4.1包含以下子步骤：4.1.1、生成贯入杆模型、试筒模型和底座模型；根据混合后的p组土石混合料的颗粒组成、几何形状及物理参数，生成虚拟的p组土石混合料颗粒，其中，物理参数为混合后的p组土石混合料中每组土石混合料的密度、模量与泊松比；在重力作用下完成初步密实；选用带阻尼的线弹性模型，输入步骤3确定的土石混合料的最佳微力学参数，施加激振力，对混合后的p组土石混合料振动击实形成p个虚拟试件。4.1.2、对步骤4.1.1得到的p个虚拟试件分别加荷载，每级加载的压力pi，读取每级压力下的回弹变形li。每级加载后都运算至模型达到静定状态，并记录承载板位置，同时卸载使模型的变形恢复。运算至静定状态后，再次记录承载板位置，同时施加下一级荷载。重复以上步骤进行加载，直至最后一级荷载；按公式计算，得到p个虚拟试件的数值试验回弹模量e21,s，(5)子步骤4.2中，具体为：对a22,…a2j,…,a2m组土石混合料分别与a1组中的p组土石混合料进行混合，得到混合后的[(m-1)×p]组土石混合料按照子步骤4.1.1-4.1.2分别进行数值模拟试验，得到混合后的[(m-1)×p]组土石混合料的数值试验回弹模量，分别为e22,s,…,e2j,s,...,e2m,s，其中，(6)步骤6包含以下子步骤：步骤6包含以下子步骤：6.1，根据回归系数a和b，确定土石混合料的回弹模量回归方程为：其中，e40为颗粒粒径小于等于40mm的土石混合料的实测回弹模量，单位为mpa；e为待测土石混合料的回弹模量，单位为mpa；ρmax为颗粒粒径小于等于40mm的土石混合料的实测最大干密度，单位为g/cm3；λ为颗粒粒径大于40mm的土石混合料的质量占比；m1为颗粒粒径小于等于40mm的土石混合料的质量，单位为g；m2为颗粒粒径大于40mm的土石混合料的质量，单位为g；η＝dmax/40，dmax为土石混合料的最大粒径，单位为mm；，η为颗粒特征参数；a，b分别为常数；6.2，取待测土石混合料，将待测土石混合料分为b1组和b2组，其中，b1组土石混合料的颗粒粒径为0-40mm，b2组土石混合料的颗粒粒径大于40mm；测定b1组的土石混合料的实测最大干密度ρmax和实测回弹模量e40，计算出b2土石混合料的质量占比λ和待测土石混合料的最大粒径与40的比值η；将ρmax、e40、λ、η代入土石混合料的回弹模量回归方程，得到待测土石混合料的回弹模量e。与现有技术相比，本发明的有益效果为：克服了现有技术只能测试粒径不超过40mm土石混合料回弹模量的不足；通过实际测量与数值模拟相结合的方法，建立了颗粒粒径大于40mm的含大颗粒的土石混合料的回弹模量拟合方程，通过该方程实现了对待测土石混合料回弹模量的快速确定；能够有效预测路基土中含粒径大于40mm的土石混合料的回弹模量，为土石混填路基的设计提供精确参数，提高了土石混填路基设计的可靠性；同时，对提升土石混填路基的强度与承载力，降低土石混填路基的病害发生率，延长道路的使用寿命都具有重要意义。附图说明下面结合附图和具体实施例对本发明做进一步详细说明。图1是本发明的土石混合料的回弹模量确定方法的流程示意图。图2是本发明实施例中5组试样的实测回弹模量与数值试验回弹模量的对比图；其中，a为5组试样的实测回弹模量；b为5组试样的数值试验回弹模量。图3是本发明实施例中最大粒径60mm的土石混合料的数值试验回弹模量与最大粒径40mm的土石混合料的实测回弹模量对比图；其中，a为最大粒径40mm的土石混合料的实测回弹模量，b为最大粒径60mm且40-60mm颗粒质量占比为10％的土石混合料数值试验回弹模量，c为最大粒径60mm且40-60mm颗粒质量占比为20％的土石混合料数值试验回弹模量，d为最大粒径60mm且40-60mm颗粒质量占比为30％的土石混合料数值试验回弹模量。图4是本发明实施例中最大粒径为80mm的土石混合料的数值试验回弹模量与最大粒径40mm的土石混合料的实测回弹模量对比图；其中，a为最大粒径40mm的土石混合料的实测回弹模量，b为最大粒径80mm且40-80mm颗粒质量占比为10％的土石混合料数值试验回弹模量，c为最大粒径80mm且40-80mm颗粒质量占比为20％的土石混合料数值试验回弹模量，d为最大粒径80mm且40-80mm颗粒质量占比为30％的土石混合料数值试验回弹模量。图5是本发明实施例中最大粒径为100mm的土石混合料的数值试验回弹模量与最大粒径40mm的土石混合料的实测回弹模量对比图；其中，a为最大粒径40mm的土石混合料的实测回弹模量，b为最大粒径100mm且40-100mm颗粒质量占比为10％的土石混合料数值试验回弹模量，c为最大粒径100mm且40-100mm颗粒质量占比为20％的土石混合料数值试验回弹模量，d为最大粒径100mm且40-100mm颗粒质量占比为30％的土石混合料数值试验回弹模量。图6是本发明的实施例中步骤6中b1组土石混填路基土的干密度-含水率关系图。具体实施方式下面将结合实施例对本发明的实施方案进行详细描述，但是本领域的技术人员将会理解，以下实施例仅用于说明本发明，而不应视为限制本发明的范围。在实际工程应用中，比如铺路过程中，需要测定较长一段土石混填路基不同点处土石混合料的回弹模量，而针对含粒径大于40mm的土石混合料，目前传统的方法是无法测定的，因此，本发明提出以下方法对含粒径大于40mm的土石混合料回弹模量进行确定。以下给出本发明的土石混合料的回弹模量的具体确定方法，参考图1，包括以下步骤：步骤1：取待测定的某条土石混填路基上某一处的土石混合料作为试样，将试样分为a1组和a2组，其中a1组的混合料粒径为0～40mm，a2组的混合料的最小粒径大于40mm；包含以下子步骤：1.1，将a1组的土石混合料分为粒径范围不同的n组，分别为a11,…a1i,…,a1n，i＝1,2,…n，n≥2；其中，a1i的粒径范围为di-1～dimm，d0＝0，dn＝40mm。优选的，根据工程实际应用过程最佳筛分粒径，将a1组分为两组，a11组的土石混合料的粒径范围为0～5mm，a12组的土石混合料的粒径范围为5～40mm。1.2，将a2组的土石混合料配制成粒径范围不同的m组土石混合料，分别为a21,…a2j,…,a2m，j＝1,2,…m，m≥3；其中，a2j的粒径范围为d40～djmm，40＜d1＜…＜dm，d40＝40mm；具体地，可选dj-dj-1＝20mm。优选的，为了使后期拟合结果更加准确，将a2组分为至少三组，如对于最大粒径为100mm的混合料，粒径范围分别为：40～60mm、40～80mm、40～100mm。步骤2：测定a1组土石混合料的实测回弹模量ex，x＝1,2,….,p，p≥4。传统是直接用仪器测试的，不用分组，本申请进行分组是为了获得不同的数值用于后续拟合，具体包含以下子步骤：2.1，将a11至a1n这n组混合料按照一定质量比进行混合，分别混合p组，需要至少四个点才能保证后期拟合精度，每组的比例关系各异；如n＝2时，a11与a12的质量比可为70:30、60:40、50:50、40:60与30:70；对于5组土石混合料中的第一组土石混合料进行取样，烘干至恒重，重复多次取样并烘干，具体地为重复5-6次，测定每份试样的干密度。包含以下子步骤：2.1.1、对每份试样，以水质量与土石混合料试样的质量比为2％～3％的比例加水，搅拌均匀，装入塑料袋或塑料桶内密封一段时间，得润湿试样。2.1.2、将润湿试样装入试模，采用振动击实仪进行振动击实，脱模，得脱模试件，测试脱模试件的质量m、脱模试件的高度h和脱模试件的直径d。2.1.3、对脱模试件进行破碎，得破碎料，称取至少两份相同重量的破碎料，烘干至恒重，分别测试多份破碎料的含水率，计算多份破碎料的含水率平均值，即为试样的含水率ω。2.1.4、根据脱模试件的质量m、脱模试件的高度h、脱模试件的直径d和试样的含水率ω，计算出该份试样的干密度2.2，绘制干密度-含水率的关系曲线，通过干密度-含水率的关系曲线，得到第一组土石混合料的最大实测干密度ρ1,max，最大实测干密度ρ1,max对应的含水率为最佳含水率ω1。2.3，按照步骤2.2确定的最大实测干密度ρ1,max和最佳含水率ω1，并结合预设压实度k，采用振动击实法，将第一组土石混合料成型为圆柱体试件；施加竖向荷载pi，读取每级荷载下的回弹变形li，得到第一组土石混合料的实测回弹模量其中，d是圆柱体试件的直径，μ0是土石混合料的泊松比；pi每级加载的荷载，li每级荷载对应的回弹变形量。2.4，按照步骤2.1至步骤2.3，测试第二组至第p组土石混合料的实测回弹模量e2,….,ep，并结合第一组土石混合料的实测回弹模量，得到p组土石混合料的实测回弹模量ex。步骤3：确定土石混合料的最佳微力学参数；包含以下子步骤：3.1，对p组土石混合料中每组土石混合料分别进行回弹模量的数值模拟试验，分别得到每组土石混合料的数值试验回弹模量ex,s，x＝1,2,….,p；包含以下子步骤：3.1.1、生成贯入杆模型、试筒模型和底座模型；根据p组土石混合料模型的颗粒组成、几何形状及物理参数，生成p组土石混合料的颗粒模型；采用重力作用，对p组土石混合料的颗粒模型分别进行初步密实，得p个初始试件；采用带阻尼的线弹性模型，输入土石混合料中土颗粒与石颗粒的微力学参数的初始值，施加激振力的初始值，对p个初始试件进行振动击实，即得p个虚拟试件；其中振动击实时间根据压实度k选定，如表1所示。表1不同压实度要求下试件振动成型所需对应时间压实度k(％)969594939290振动时间t(s)6765646260583.1.2、对p个虚拟试件分别逐级加荷载，每级加载的荷载为pi，读取每级荷载下的回弹变形li；每级加载后都运算至模型达到静定状态，直至最后一级荷载。得到土石混合料的数值试验回弹模量其中，ex,s为p组土石混合料中的第x组土石混合料的数值试验回弹模量；d为p个虚拟试件中每个试件的直径；μ0为土石混合料的泊松比；pi每级竖向加载的荷载，li每级荷载对应的回弹变形量。3.2，根据第二预设公式计算实测回弹模量ex与数值试验回弹模量ex,s的累计相对偏差δ；反复调整土石混合料的微力学参数，得到多个δ，当δ最小时对应的微力学参数为最佳微力学参数；第二预设公式中，ex是p组土石混合料中第x组土石混合料的实测回弹模量，ex,s是p组土石混合料中第x组土石混合料的数值试验回弹模量，δ是ex与ex,s的相对偏差累加值。进一步地，子步骤3.1.1中，所述物理参数为p组土石混合料中每组土石混合料的密度、模量和泊松比。进一步地，子步骤3.2中，土石混合料的微力学参数包含：土颗粒的湿润摩擦系数、土颗粒的阻尼系数、石颗粒的湿润摩擦系数、石颗粒的阻尼系数。更进一步地，所述土颗粒的湿润摩擦系数的取值范围为：0.20～0.46，土颗粒的阻尼系数取值范围为0.30～0.40，石颗粒的湿润摩擦系数取值范围为0.42～0.46，石颗粒的阻尼系数取值范围为0.65～0.75。步骤4：确定a2组与a1组混合后的土石混合料的数值试验回弹模量e2j,s；包含以下子步骤：4.1，将a21组土石混合料与a1组中的p组土石混合料分别进行混合，根据土石混合料的最佳微力学参数对混合后的p组土石混合料分别进行数值模拟试验，得到混合后土石混合料的数值试验回弹模量e21,s，包含以下子步骤：4.1.1、生成贯入杆、试筒和底座的模型；根据混合后的p组土石混合料的颗粒组成、几何形状及物理参数，生成虚拟的p组土石混合料颗粒，其中物理参数为p组土石混合料中每组土石混合料的密度、模量与泊松比；在重力作用下完成初步密实；选用带阻尼的线弹性模型，输入步骤3确定的土石混合料的最佳微力学参数，施加激振力，对混合后的土石混合料振动击实形成试件，共有p个虚拟试件。4.1.2、对步骤4.1.1得到的p个虚拟试件分别加荷载，每级加载的压力pi，读取每级压力下的回弹变形li。每级加载后都运算至模型达到静定状态，并记录承载板位置，同时卸载使模型的变形恢复。运算至静定状态后，再次记录承载板位置，同时施加下一级荷载。重复以上步骤进行加载，直至最后一级荷载；按公式计算，得到p个虚拟试件数值试验回弹模量e21,s，4.2，将a22,…a2j,…,a2m组土石混合料分别与a1组中的p组土石混合料进行混合，对混合后的[(m-1)×p]组土石混合料分别进行数值模拟试验，得到混合后的[(m-1)×p]组土石混合料的数值试验回弹模量，分别为e22,s,…,e2j,s,...,e2m,s，结合e21,s，得到a2组与a1组混合后的土石混合料的数值试验回弹模量e2j,s，其中，具体为：对a22,…a2j,…,a2m组土石混合料分别与a1组中的p组土石混合料进行混合，得到的混合后的[(m-1)×p]组土石混合料按照子步骤4.1.1-4.1.2，分别进行数值模拟试验，得到混合后的[(m-1)×p]组土石混合料的数值试验回弹模量，分别为e22,s,…,e2j,s,...,e2m,s，其中，更进一步的，为保证后期拟合精度，可根据需要将a2j组和a1组的任一组的土石混合料再按照不同质量比细分为r组，细分后的土石混合料的数值回弹模量测试过程与步骤2相同，而此时得到的ρ2j,s的个数为p×r个。步骤5，根据a1组土石混合料的实测回弹模量ex和a2组与a1组混合后的土石混合料的数值试验回弹模量e2j,s，通过第一预设公式：确定回归系数a和b；第一预设公式中，ex为a1组土石混合料中第x组的实测回弹模量，单位为mpa；e2j,s为a2组与a1组混合后的土石混合料的数值试验回弹模量，单位为mpa；ρx,max为a1组土石混合料的实测最大干密度，单位为g/cm3；λ为a2组土石混合料的质量占比；为a1组土石混合料的质量，单位为g；为a2组土石混合料的质量，单位为g；η＝d2max/40，d2max为a2组土石混合料的最大粒径，单位为mm，η为颗粒特征参数；a，b分别为回归系数；步骤6：根据回归系数a和b，确定土石混合料的回弹模量回归方程，通过土石混合料的回弹模量回归方程，确定待测土石混合料的回弹模量e。包含以下子步骤：6.1，根据回归系数a和b，确定土石混合料的回弹模量回归方程为：其中，e40为颗粒粒径小于等于40mm的土石混合料的实测回弹模量，单位为mpa；e为待测土石混合料的回弹模量，单位为mpa；ρmax为颗粒粒径小于等于40mm的土石混合料的实测最大干密度，单位为g/cm3；λ为颗粒粒径大于40mm的土石混合料的质量占比；m1为颗粒粒径小于等于40mm的土石混合料的质量，单位为g；m2为颗粒粒径大于40mm的土石混合料的质量，单位为g；η＝dmax/40，dmax为土石混合料的最大粒径，单位为mm；η为颗粒特征参数；a，b分别为常数。6.2，取待测土石混合料，将待测土石混合料分为b1组和b2组，其中，b1组土石混合料的颗粒粒径为0-40mm，b2组土石混合料的颗粒粒径大于40mm；测定b1组的土石混合料的实测最大干密度ρmax和实测回弹模量e40，计算出b2土石混合料的质量占比λ和待测土石混合料的最大粒径与40的比值η；将ρmax、e40、λ、η代入土石混合料的回弹模量回归方程，得到待测土石混合料的回弹模量e。实施例将某种土石混合料用于铺设路基时需要测定该混合料的回弹模量，对于含粒径大于40mm的土石混合料，取足量的土石混合料试样，进行如下步骤：步骤1：采用40mm筛将试样分成a1(0～40mm)与a2(40～100mm)两部分；再采用5mm筛将a1分成a11(0～5mm)与a12(5～40mm)两部分；根据最大粒径，又将a2组分为三种情况：a21(40～60mm)、a22(40～80mm)、a23(40～100mm)。步骤2：通过室内振动击实试验测试粒径为0～40mm的土石混合料的实测回弹模量，具体为：首先，按以下比例将a11与a12混合为不同质量比的5组，a11:a12为70:30、60:40、50:50、40:60、30:70，并对这5组混合料编号为①、②、③、④、⑤；其次，对这5组混合料进行室内振动击实试验，振动压实仪参数设置为：振动频率为25hz、激振力5.3kn、振动击实75s、名义振幅1.3mm、上车质量107.08～115.01kg、下车质量170.59～179.33kg、偏心块夹角0°、静面压力154.0～163.2kpa。对于①号试样，取若干土石混填路基土试样置于烘箱中烘干至恒重，温度控制范围为105±5℃，时间不少于6小时；通过四分法配制5份相同的试样，混合均匀，以水质量与土石混填路基土干重比为2～3％向试样中加入水，拌合至均匀状态，装入塑料袋或塑料桶中封闭闷料一夜；取土样装入直径×高度＝150mm×230mm的试模中，并将试模固定到振动击实仪上，振动击实时间75s；用电动脱模机将试件从试模中完整脱出，量取试件高度h和直径d，称量试件质量m；将试件破碎从试件中部取约600g(平行两份)土石混填路基土放入小铁盒中，然后置于105±5℃烘箱中烘干约6小时，测试含水率ω(以百分数计，取两份样的平均值，两者误差不超过1％)，计算试样的干密度ρ1，计算公式(1)如下：重复上述步骤，得到①号混合料其他4份试样的含水率ω与干密度ρ1；绘制干密度-含水率关系曲线，其中横坐标为含水率，纵坐标为干密度，曲线峰值所对应的干密度为最大干密度ρ1,max。按上述同样步骤测定剩余②、③、④、⑤号土石混合料的最大干密度ρ2,max、ρ3,max、ρ4,max、ρ5,max。结果如表2所示。表2粒径小于40mm土石混合料的实测最大干密度和最佳含水率编号①②③④⑤土石比(a11:a12)70:3060:4050:5040:6030:70最大干密度(g/cm3)1.9631.9671.9801.9841.978最佳含水率(％)14.514.313.813.312.7采用振动压实方法成型0～40mm土石混合料的圆柱体试件，直径为152mm，高120mm；压实度k为95％，振动时间65s。将圆柱体试件浸水4天，测试①号至⑤号试样的实测回弹模量ex。具体地，对于①号试样，先进行预压：承载板上赋予最大单位压力p进行预压，p＝200kpa。预压一次；预压完成后，缓慢将最大单位压力卸载，重新记录承载板位置，准备试验。其次进行加载：将预定最大单位压力分成4份，每级50kpa，作为每级加载的压力pi，并读取每级压力下的回弹变形li。重复以上步骤进行逐级加载，直至最后一级荷载；按以下公式(2)计算，其中，μ0为土石混合料的泊松比0.25，结果如表3所示。表3粒径0～40mm土石混合料的5组试样的实测回弹模量编号①②③④⑤土石比(a11:a12)70:3060:4050:5040:6030:70ex(mpa)68.171.577.483.285.3步骤3：按照室内振动击实试验条件，对步骤2中的5组试样进行数值模拟试验，确定土石混合料的最佳微力学参数：首先，采用pfc软件，生成承载板、试筒和底座的模型，各种构件的模量取为30gpa(与钢板的基本一致)、泊松比为0.2；试筒内壁直径150mm、外壁160mm、高150mm，承载板直径149mm、高25mm，底座直径180mm、高25mm。然后，根据土石混填路基土的颗粒组成、几何形状以及物理参数，生成相应的土石混合料。土颗粒均简化为圆形，各种粒径的土颗粒数按照实际的颗粒组成计算；物理参数包括土颗粒(≦5mm)与石颗粒(5mm﹥)的模量，泊松比和实测容重，其中土颗粒的模量为100mpa，泊松比为0.30，石颗粒的模量为20000mpa，泊松比为0.25；在重力作用下土石混合料完成自由初步排布后，密实形成试件；施加振动频率为25hz与振幅1.3mm的激振力，激振力为5.3kn，振动击实时间根据压实度k选定，如表1所示；最后，进行回弹模量数值模拟试验，具体为：(1)预压：承载板上赋予最大单位压力p进行预压，p＝200kpa。预压一次，运算至模型达到静定状态。预压完成后，缓慢将最大单位压力卸载，并运算达到静定状态，重新记录承载板位置，准备试验。(2)加载：将预定最大单位压力分成4份，作为每级加载的压力pi，并读取每级压力下的回弹变形li。每级加载后都运算至模型达到静定状态，并记录承载板位置，同时卸载使模型的变形恢复。运算至静定状态后，再次记录承载板位置，同时施加下一级荷载。重复以上步骤进行加载，直至最后一级荷载；按以下公式(3)计算，其中，μ0为土石混合料的泊松比0.25，得到土石混合料的数值试验回弹模量ex,s。图2为5组试样的实测回弹模量与数值试验回弹模量的对比图，从图中可以看出，两者结果很接近。获得数值试验回弹模量ex,s后，分别按下式(4)计算5组试样的实测回弹模量ex与数值试验回弹模量ex,s的累计相对偏差δ；参考表4，以0.02为间隔，反复调整土颗粒与石颗粒的湿润摩擦系数与阻尼系数，得到多个δ；当土颗粒的湿润摩擦系数为0.31、阻尼系数为0.34，石颗粒的湿润摩擦系数为0.41、阻尼系数为0.65时，实测回弹模量与数值试验回弹模量的累积相对误差δ＝7.2％，达到最小值，即此时的土颗粒与石颗粒的微力学参数为土石混合料的最佳微力学参数。表4土颗粒与颗粒的微力学参数参考范围材料湿润摩擦系数阻尼系数土颗粒0.20～0.460.30～0.40石颗粒0.42～0.460.65～0.75步骤4：对于最大粒径为60mm的混合料，按以下质量比将a21与a1混合为不同质量比的三组，三组中颗粒粒径大于40mm的混合料的质量占比分别为10％、20％、30％，具体见表5；根据步骤3中确定的土石混合料的最佳微力学参数，具体地为：土颗粒的摩擦系数0.31、阻尼系数0.34，石颗粒的摩擦系数0.41、阻尼系数0.65，进行数值模拟试验：施加激振力，振动击实形成试件；对试件加荷载，记录每级加载的压力pi，并读取每级压力下的回弹变形li，按上述公式(3)计算，得到该组混合料的数值试验回弹模量e21,s，结果如图3所示，图3中，a为最大粒径40mm的土石混合料的实测回弹模量，b为最大粒径60mm且40-60mm颗粒质量占比为10％的土石混合料数值试验回弹模量，c为最大粒径60mm且40-60mm颗粒质量占比为20％的土石混合料数值试验回弹模量，d为最大粒径60mm且40-60mm颗粒质量占比为30％的土石混合料数值试验回弹模量。按上述步骤2中的最大干密度确定方法，确定最大粒径为80mm和100mm土石混合料的多组最大干密度ρ2j,s，按步骤4的数值模拟试验过程，分别得到最大粒径为80mm和100mm土石混合料的数值试验结果如图4和图5所示；图4中，a为最大粒径40mm的土石混合料的实测回弹模量，b为最大粒径80mm且40-80mm颗粒质量占比为10％的土石混合料数值试验回弹模量，c为最大粒径80mm且40-80mm颗粒质量占比为20％的土石混合料数值试验回弹模量，d为最大粒径80mm且40-80mm颗粒质量占比为30％的土石混合料数值试验回弹模量。图5中，a为最大粒径40mm的土石混合料的实测回弹模量，b为最大粒径100mm且40-100mm颗粒质量占比为10％的土石混合料数值试验回弹模量，c为最大粒径100mm且40-100mm颗粒质量占比为20％的土石混合料数值试验回弹模量，d为最大粒径100mm且40-100mm颗粒质量占比为30％的土石混合料数值试验回弹模量。表5数值模拟试验方案步骤5：根据步骤2得到的小于40mm的土石混合料的实测回弹模量e1、e2、e3、e4、e5以及步骤4得到的多个数值试验回弹模量e2j,s，通过matlab按下式(5)进行拟合，得到回归系数a＝1.56，b＝1.2。步骤6：将回归系数a＝1.56，b＝1.2代入土石混合料的回弹模量回归方程，即得本实施例中0～100mm土石混合料的回弹模量e的计算公式为：取需要测定路段的现场挖坑土样，按照步骤1将土样分为b1组和b2组，并对b1组(0～40mm)土石混填路基土进行振动击实试验，测定的干密度-含水率曲线如图6所示。由图6可得，0～40mm土石混合料的最大干密度ρmax＝1.903g/cm3，最佳含水率ω40＝14.8％，根据ρmax和ω40对土样进行密实测定0～40mm土石混合料的回弹模量e40＝73.2mpa。再进行颗粒分析，计算b2组土石混合料在待测定土石混合料中的质量占比λ和颗粒特征参数η。将得到的测试参数ρmax、e40、λ、η代入上述回弹模量e的计算公式中，计算得到0～100mm土石混合料的回弹模量e，具体测试结果和计算结果如表6所示，在该待测定路段上的不同纵向点和不同横向点取土样作为试样，计算结果如表6所示。表60～100mm土石混合料的回弹模量计算结果从表6可以看出，断面ii处的不同横向位置的土石混合料的回弹模量基本相同，说明断面ii处的土石颗粒在横向分布均匀；断面i和断面iii处的不同横向位置的土石混合料的回弹模量存在差异，说明断面i和断面iii处的土石颗粒的横向分布存在较大差异；而对比断面i、ii、iii的回弹模量，可以看出，三处存在较大差异，因此，对于同一路段的不同位置的回弹模量的预测具有非常重要的意义，本发明所述方法可有效的实现该段路的力学性能和使用寿命的预测，具有实际应用价值。图2为本发明实施例中5组土石混合料的实测回弹模量与数值试验回弹模量的对比图，从图2中可以看出，采用本发明的数值模拟试验方法得到的数值试验回弹模量与实际测量的回弹模量结果接近，说明本发明的数值模拟试验方法结果可靠。图3为最大粒径为60mm的土石混合料的数值试验回弹模量与最大粒径40mm土石混合料的实测回弹模量对比图。从图3可以看出，采用本发明中的数值模拟试验方法确定的最大粒径为60mm的土石混合料的数值试验回弹模量，与最大粒径为40mm的土石混合料的实测回弹模量，表现出基本相同的规律，都随着土石比的减小而逐渐增加。图4为最大粒径为80mm的土石混合料的回弹模量与最大粒径40mm土石混合料的实测回弹模量对比图。从图4可以看出，采用本发明中的数值模拟试验方法确定的最大粒径为80mm的土石混合料的数值试验回弹模量，与最大粒径为40mm的土石混合料的实测回弹模量，表现出基本相同的规律，都随着土石比的减小而逐渐增加；同时可以看出，最大粒径为80mm的土石混合料的回弹模量明显升高。图5为最大粒径为100mm的土石混合料的回弹模量与最大粒径40mm土石混合料的实测回弹模量对比图。从图5可以看出，采用本发明中的数值模拟试验方法确定的最大粒径为100mm的土石混合料的数值试验回弹模量，与最大粒径为40mm的土石混合料的实测回弹模量，表现出基本相同的规律，都随着土石比的减小而逐渐增加；同时可以看出，最大粒径为100mm的土石混合料的回弹模量明显升高。对比图3、图4和图5，可以看出，随着试样中土石混合料最大粒径的逐渐增大，土石混合料的回弹模量逐渐增大；随着巨粒颗粒质量占比的增大，土石混合料的回弹模量逐渐增大，且所有的试验结果与最大粒径40mm土石混合料的实测回弹模量表现出基本相同的规律，说明基于本发明中的数值模拟试验得到的土石混合料的回弹模量确定方法的结果可靠。综上所述，使用本发明的方法可以快速准确的得到铺设路基不同路段上的土石混合料的回弹模量，为土石混填路基的压实质量控制提供了关键参数，客观上提升了土石混填路基的施工质量与使用品质；同时，可有效预测路基土中含粒径大于40mm颗粒的土石混合料的回弹模量，为土石混填路基的设计提供精确参数，提高了土石混填路基设计的可靠性，对降低土石混填路基的病害发生率，延长道路的使用寿命都具有重要意义。本发明在实际工程应用中，可应用于铺路过程中，需要测定较长一段土石混填路基不同断面处土石混合料的回弹模量，用于预测该段路的力学性能和使用寿命，或者某种土石混合料用于铺设路基时需要确定该混合料的回弹模量，主要针对土石混合料中含粒径大于40mm颗粒的情况。本发明中的土石混合料为土颗粒和石颗粒组成的混合料。土颗粒为粒径小于等于5mm的颗粒，石颗粒为粒径大于5mm的粒径，说明书附图中的土石比为土颗粒与石颗粒的质量比。本发明中的土石混合料的微力学参数是指土颗粒的湿润摩擦系数和阻尼系数，石颗粒的湿润摩擦系数和阻尼系数。显然，本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些改动和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。当前第1页12