CN115455712A - 一种基于响应曲面法的土体压实度检测方法 - Google Patents

一种基于响应曲面法的土体压实度检测方法 Download PDF

Info

Publication number
CN115455712A
CN115455712A CN202211152176.XA CN202211152176A CN115455712A CN 115455712 A CN115455712 A CN 115455712A CN 202211152176 A CN202211152176 A CN 202211152176A CN 115455712 A CN115455712 A CN 115455712A
Authority
CN
China
Prior art keywords
conductivity
dielectric constant
soil body
soil
penetration resistance
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.)
Pending
Application number
CN202211152176.XA
Other languages
English (en)
Inventor
王柳江
臧耀辉
毛航宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hohai University HHU
Original Assignee
Hohai University HHU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hohai University HHU filed Critical Hohai University HHU
Priority to CN202211152176.XA priority Critical patent/CN115455712A/zh
Publication of CN115455712A publication Critical patent/CN115455712A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/22—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance
    • G01N27/221—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating capacitance by investigating the dielectric properties
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Geometry (AREA)
  • Molecular Biology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Electrochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Computing Systems (AREA)
  • Computational Linguistics (AREA)
  • Evolutionary Biology (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Artificial Intelligence (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

本发明涉及一种基于响应曲面法的土体压实度检测方法，先取现场测试土体开展室内标定试验，建立土体的介电常数、电导率和最大贯入阻力随质量含水率和干密度变化的二阶响应曲面模型，其次，研究土体的电导率和介电常数之间的相关性，建立土体电导率的修正模型并计算电导率修正值；然后，将实测介电常数、实测最大贯入阻力以及电导率修正值代入二阶响应曲面模型，以实测值和计算值均方根误差最小建立目标函数；最后，采用粒子群优化算法进行反演，得到土体质量含水率和干密度，计算出土体压实度，本发明的测试方法能够对现场碾压后的土体压实度进行快速精准的检测，测得土体质量含水率和干密度，适用于不同的工程填筑土体。

Description

一种基于响应曲面法的土体压实度检测方法
技术领域
本发明涉及填筑土体含水率和干密度快速检测技术领域，具体为一种基于响应曲面法的土体压实度检测方法。
背景技术
土料在水利、交通、市政等基础设施建设领域应用十分广泛，如在水利工程中常作为土坝、心墙、堤防等挡水结构物的填筑材料，压实度是填方工程施工管理最为关键的指标之一，直接关系到填方工程的强度、变形及稳定性，其中质量含水率和干密度是确定压实度的两个关键参数，传统检测方法为：对于土体质量含水率，采用烘干法，而对于土体干密度，通常采用环刀法、灌砂法、灌水法等与烘干法联合测定；尽管这些传统检测方法测试结果准确可靠，但属于有损检验，且存在抽样点少、检测速度慢、费用高、对施工干扰大等缺点，难以满足当前快速机械化施工要求；因此，工程界一直在寻找快速无损、精度高的土体质量含水率和干密度快速测试方法，以提高土体压实度的检测效率。
目前，快速无损检测土体压实度的方法有基于时域反射技术（TDR）或频域反射技术（FDR）的一步检测法（陈赟, 陈伟, 陈仁朋, 等. TDR 联合监测土体含水率 和干密度的传感器的设计及应用[J]. 岩石力学与工程学报, 2011, 30(2): 418-426.），该方法主要是将土体介电常数、电导率与含水率和干密度进行一次线性拟合，进而推导出含水率和干密度的计算公式。
由于该一步检测法未考虑含水率和干密度对土体介电常数和电导率影响的二次效应和交互效应，且测试值介电常数难以全面反映含水率和干密度的影响，导致测得的土体质量含水率和干密度精度偏低，使其无法适用于不同的工程填筑土体，难以对现场碾压后的土体压实度进行准确评价，为此，我们提出一种基于响应曲面法的土体压实度检测方法。
发明内容
为解决现有技术中存在的问题，本发明提供一种基于响应曲面法的土体压实度检测方法。
为达到上述目的，本发明采用的方案如下：
一种基于响应曲面法的土体压实度检测方法，包括以下步骤：
S1：取部分现场土体做室内标定实验，制备不同含水率和干密度的试样，利用FDR传感器测试土体的介电常数和电导率，利用贯入仪测试土体最大贯入阻力；
S2：以土体的介电常数、电导率和最大贯入阻力作为响应目标，土体含水率和干密度作为影响因素，考虑线性效应、二次效应和交互效应，建立土体的介电常数、电导率和最大贯入阻力随质量含水率和干密度变化的二阶响应曲面模型；
S3：根据土体的电导率和介电常数之间的相关性，建立土体电导率的修正模型；
S4：将FDR传感器、贯入仪与电动反力升降装置连接，形成一套便携式快速检测装置，通过电机控制升降，将FDR传感器和贯入仪的探针***压实土体，分别测得压实土体的实测介电常数和实测最大贯入阻力，将实测介电常数代入土体电导率的修正模型，计算得到土体的电导率修正值；
S5：将实测介电常数、电导率修正值和实测最大贯入阻力作为已知量，代入对应的二阶响应曲面模型，以三个目标值的实测值和模型计算值均方根误差最小为目标函数，基于粒子群优化算法反演得到最优的一组计算参数，即为现场压实土体的质量含水率和干密度。
作为优选的技术方案：
如上所述的一种基于响应曲面法的土体压实度检测方法，步骤S2中，土体的介电常数、电导率和最大贯入阻力随质量含水率和干密度变化的二阶响应曲面模型，其表达式如下：
Figure DEST_PATH_IMAGE002
式中，
Figure DEST_PATH_IMAGE004
、
Figure DEST_PATH_IMAGE006
、
Figure DEST_PATH_IMAGE008
分别为土体的介电常数、电导率和最大贯入阻力的响应曲面函数；
Figure DEST_PATH_IMAGE010
、
Figure DEST_PATH_IMAGE012
、
Figure DEST_PATH_IMAGE014
为截断系数，
Figure DEST_PATH_IMAGE016
、
Figure DEST_PATH_IMAGE018
、
Figure DEST_PATH_IMAGE020
、
Figure DEST_PATH_IMAGE022
、
Figure DEST_PATH_IMAGE024
、
Figure DEST_PATH_IMAGE026
为线性项系数，
Figure DEST_PATH_IMAGE028
、
Figure DEST_PATH_IMAGE030
、
Figure DEST_PATH_IMAGE032
为交互项系数，
Figure DEST_PATH_IMAGE034
、
Figure DEST_PATH_IMAGE036
、
Figure DEST_PATH_IMAGE038
、
Figure DEST_PATH_IMAGE040
、
Figure DEST_PATH_IMAGE042
、
Figure DEST_PATH_IMAGE044
为二次项系数，均为通过最小二乘法拟合后获得的参数，由率定试验拟合确定；
Figure DEST_PATH_IMAGE046
为土体质量含水率，
Figure DEST_PATH_IMAGE048
为干密度。
如上所述的一种基于响应曲面法的土体压实度检测方法，步骤S3中的土体电导率修正模型采用如下函数：
Figure DEST_PATH_IMAGE050
式中，
Figure DEST_PATH_IMAGE052
,
Figure DEST_PATH_IMAGE054
,
Figure DEST_PATH_IMAGE056
,
Figure DEST_PATH_IMAGE058
为模型参数，
Figure DEST_PATH_IMAGE060
为电导率修正值，
Figure DEST_PATH_IMAGE062
为介电常数。
如上所述的一种基于响应曲面法的土体压实度检测方法，步骤S5中，选取实测介电常数、电导率修正值和实测最大贯入阻力作为已知量，采用等权重法构造目标函数：
Figure DEST_PATH_IMAGE064
式中，
Figure DEST_PATH_IMAGE066
为实测介电常数，
Figure DEST_PATH_IMAGE060A
为电导率修正值，P为实测最大贯入阻力；
Figure DEST_PATH_IMAGE004A
、
Figure DEST_PATH_IMAGE006A
和
Figure DEST_PATH_IMAGE008A
为土体的介电常数、电导率和最大贯入阻力对应的响应曲面函数。
如上所述的一种基于响应曲面法的土体压实度检测方法，步骤S4中，FDR传感器和贯入仪与电动反力升降装置的传力杆底部采用可拆卸式的螺纹连接，现场检测时，先连接贯入仪测最大贯入阻力，再连接FDR传感器测介电常数，数值均记录在采集仪中；测试时，通过电机控制探针贯入速度；
其中，FDR传感器包括电路板、设置在电路板上方的FDR传感器外螺纹柱、设置在电路板下方的FDR探针、以及与电路板连接的FDR传感器数据传导线；贯入仪包括荷载传感器、设置在荷载传感器上方的贯入仪外螺纹柱、设置在荷载传感器下方的微型探针，以及与荷载传感器连接的贯入仪数据传导线；电动反力升降装置包括底板和与底板平行的支撑板，底板的两侧各设置有一个脚踏板，底板和支撑板之间通过4根支撑螺杆固定连接，底板中间开设有可使FDR传感器和贯入仪穿过的圆形通孔，支撑板垂直穿设有传力杆，传力杆内部轴向上开设有与FDR传感器外螺纹柱和贯入仪外螺纹柱适配的螺纹孔，传力杆远离圆形通孔的一端设置有升降传动装置， 升降传动装置固定在支撑板上，且一侧设置有电机，电机通过采集仪数据传导线与采集仪连接，电机连接有电源。
本发明的原理为：
首先，取现场测试土体开展室内标定试验，研究土体的介电常数、电导率和最大贯入阻力与质量含水率和干密度的关系，考虑线性效应、二次效应和交互效应，建立土体的介电常数、电导率和最大贯入阻力随质量含水率和干密度变化的二阶响应曲面模型，并通过最小二乘法拟合确定相应的模型参数，同时，研究土体的电导率和介电常数之间的相关性，建立土体电导率的修正模型；其次，将FDR传感器和贯入仪与电动反力升降装置通过螺纹连接固定形成一套便携式的快速检测设备，测试现场压实土体的实测介电常数和实测最大贯入阻力，并利用实测介电常数计算得到电导率修正值；然后，将实测介电常数、实测最大贯入阻力以及电导率修正值代入二阶响应曲面模型，以实测值和计算值均方根误差最小建立目标函数；最后，采用粒子群优化算法进行反演，得到土体质量含水率和干密度，计算出土体压实度。
本发明提供了一种基于响应曲面法的土体压实度检测方法，具备以下有益效果：
（1）首先取现场测试土体开展室内率定试验，研究土体介电常数、电导率和最大贯入阻力与质量含水率和干密度的关系，建立土体介电常数、电导率和最大贯入阻力随质量含水率和干密度变化的二阶响应曲面模型，并确定响应的模型参数；同时，研究土体电导率和介电常数之间的相关性，建立土体电导率的修正模型，克服了现场检测土体孔隙水的电导率与室内率定试验孔隙水的电导率不一致的问题；
（2）将实测介电常数、实测最大贯入阻力以及电导率修正值代入二阶响应曲面模型，以实测值和计算值均方根误差最小建立目标函数，采用粒子群优化算法进行反演，得到土体质量含水率和干密度，并给出压实度；该方法在原有一步法的基础上增加了最大贯入阻力的检测，并考虑了含水率和干密度对介电常数、电导率以及最大贯入阻力影响的二次效应和交互效应，建立的二阶响应曲面模型基本能够应用于所有类型的土体，预测得到的土体质量含水率和干密度的精度大幅度提升；
（3）土体实测介电常数和实测最大贯入阻力可在电动反力升降装置的辅助下可快速测得，可以有效提高土体压实度的检测效率，从而缓解填筑施工质量检测与快速机械化施工之间的现实矛盾；
（4）本发明有效克服了压实度传统有损检测方法中存在的抽样点少、检测速度慢、费用高、对施工干扰大等缺点；同时，克服了基于介电常数和电导率测试的一步法在不同类型土体中适用性差，测试精度低的问题，大幅度提高了土体压实度检测效率及检测精度。
附图说明
图1为本发明的基于响应曲面法的土体压实度检测方法流程图；
图2为本发明中电动反力升降装置示意图；
图3为本发明的FDR传感器示意图；
图4为本发明的贯入仪示意图；
图5为本发明的介电常数响应曲面图；
图6为本发明的电导率响应曲面图；
图7为本发明的贯入阻力响应曲面图；
图中：2-1、脚踏板；2-2、底板；2-3、支撑螺杆；2-4、传力杆；2-5、支承板；2-6、升降传动装置；2-7、电机；2-8、采集仪数据传导线；2-9、采集仪；2-10、内螺纹；3-1、FDR探针；3-2、电路板；3-3、FDR传感器外螺纹柱；3-4、FDR传感器数据传导线；3-5、微型探针；3-6、荷载传感器；3-7、贯入仪数据传导线；3-8、贯入仪外螺纹柱。
具体实施方式
下面结合具体实施方式，进一步阐述本发明。应理解，这些实施例仅用于说明本发明而不用于限制本发明的范围。此外应理解，在阅读了本发明讲授的内容之后，本领域技术人员可以对本发明作各种改动或修改，这些等价形式同样落于本申请所附权利要求书所限定的范围。
实施例1
一种基于响应曲面法的土体压实度检测方法，如图1所示，具体步骤如下：
S1：取部分现场土体做室内标定实验，制备不同含水率和干密度的试样，利用FDR传感器测试土体介电常数和电导率，利用贯入仪测试土体最大贯入阻力；
S2：以土体介电常数、电导率和最大贯入阻力作为响应目标，土体含水率和干密度作为影响因素，考虑线性效应、二次效应和交互效应，建立土体介电常数、电导率和最大贯入阻力随质量含水率和干密度变化的二阶响应曲面模型；
其中，土体的介电常数、电导率和最大贯入阻力随质量含水率和干密度变化的二阶响应曲面模型，其表达式如下：
Figure DEST_PATH_IMAGE072
式中，
Figure DEST_PATH_IMAGE004AA
、
Figure DEST_PATH_IMAGE006AA
、
Figure DEST_PATH_IMAGE008AA
分别为土体的介电常数、电导率和最大贯入阻力的响应曲面函数；
Figure DEST_PATH_IMAGE010A
、
Figure DEST_PATH_IMAGE012A
、
Figure DEST_PATH_IMAGE014A
为截断系数，
Figure DEST_PATH_IMAGE016A
、
Figure DEST_PATH_IMAGE018A
、
Figure DEST_PATH_IMAGE020A
、
Figure DEST_PATH_IMAGE022A
、
Figure DEST_PATH_IMAGE024A
、
Figure DEST_PATH_IMAGE026A
为线性项系数，
Figure DEST_PATH_IMAGE028A
、
Figure DEST_PATH_IMAGE030A
、
Figure DEST_PATH_IMAGE032A
为交互项系数，
Figure DEST_PATH_IMAGE034A
、
Figure DEST_PATH_IMAGE036A
、
Figure DEST_PATH_IMAGE038A
、
Figure DEST_PATH_IMAGE040A
、
Figure DEST_PATH_IMAGE042A
、
Figure DEST_PATH_IMAGE044A
为二次项系数，均为通过最小二乘法拟合后获得的参数，由率定试验拟合确定；
Figure DEST_PATH_IMAGE046A
为土体质量含水率，
Figure DEST_PATH_IMAGE048A
为干密度；
S3：根据土体电导率和介电常数之间的相关性，建立土体电导率的修正模型；
其中，土体电导率修正模型采用如下函数：
Figure DEST_PATH_IMAGE050A
式中，
Figure DEST_PATH_IMAGE052A
,
Figure DEST_PATH_IMAGE054A
,
Figure DEST_PATH_IMAGE056A
,
Figure DEST_PATH_IMAGE058A
为模型参数，
Figure DEST_PATH_IMAGE060AA
为电导率修正值，
Figure DEST_PATH_IMAGE062A
为介电常数；
S4：将FDR传感器、贯入仪与电动反力升降装置连接，形成一套便携式快速检测装置，通过电机控制升降，将FDR传感器和贯入仪的探针***压实土体，分别测得压实土体的介电常数和最大贯入阻力，将介电常数代入电导率修正模型，计算得到土体电导率的修正值；
其中，如图2、图3、图4所示，FDR传感器包括电路板3-2、设置在电路板3-2上方的FDR传感器外螺纹柱、设置在电路板3-2下方的FDR探针3-1、以及与电路板3-2连接的FDR传感器数据传导线3-4；贯入仪包括荷载传感器3-6、设置在荷载传感器3-6上方的贯入仪外螺纹柱3-8、设置在荷载传感器3-6下方的微型探针3-5，以及与荷载传感器3-6连接的贯入仪数据传导线3-7；电动反力升降装置包括底板2-2和与底板2-2平行的支撑板2-5，底板的两侧各设置有一个脚踏板2-1，底板2-2和支撑板2-5之间通过4根支撑螺杆2-3固定连接，底板2-2中间开设有可使FDR传感器和贯入仪穿过的圆形通孔，支撑板2-5垂直穿设有传力杆2-4，传力杆2-4内部轴向上开设有与FDR传感器外螺纹柱3-3和贯入仪外螺纹柱3-8适配的螺纹孔2-10，传力杆2-4远离圆形通孔的一端设置有升降传动装置2-6， 升降传动装置2-6固定在支撑板2-5上，且一侧设置有电机2-7，电机2-7通过采集仪数据传导线2-8与采集仪2-9连接，电机2-7连接有电源；
FDR传感器、贯入仪与电动反力升降装置连接是指FDR传感器和贯入仪与电动反力升降装置的传力杆底部采用可拆卸式的螺纹连接，现场检测时，先连接贯入仪测最大贯入阻力，再连接FDR传感器测介电常数，数值均记录在采集仪中；测试时，通过电机控制探针贯入速度；
S5：将实测介电常数、电导率修正值和实测最大贯入阻力作为已知量，代入对应的二阶响应曲面模型，以三个目标值的实测值和模型计算值均方根误差最小为目标函数，基于粒子群优化算法反演得到最优的一组计算参数，即为现场压实土体的质量含水率和干密度；
其中，目标函数采用等权重法构造，目标函数为：
Figure DEST_PATH_IMAGE064A
式中，
Figure DEST_PATH_IMAGE066A
为实测介电常数，
Figure DEST_PATH_IMAGE060AAA
为电导率修正值，P为实测最大贯入阻力；
Figure DEST_PATH_IMAGE004AAA
、
Figure DEST_PATH_IMAGE006AAA
和
Figure DEST_PATH_IMAGE008AAA
为土体的介电常数、电导率和最大贯入阻力对应的响应曲面函数。
实施例2
一种基于响应曲面法的土体压实度检测方法，具体步骤如下：
步骤1：将现场取回土料烘干并过2mm筛，过筛后土与自来水混合，制备质量含水率为14%，16%，18%，20%的土样并密封72小时，在有机玻璃桶（直径100mm，高度150mm）内按目标干密度（1.3g/cm3、1.4g/cm3、1.5g/cm3、1.6g/cm3）分五层击实为不同含水率和干密度土柱，每层击实完毕后，用刀具对表面拉毛，以减小与后续击实层之间分层现象；
步骤2：土柱制备完成后，先利用贯入仪测试土体最大贯入阻力，然后利用FDR传感器测试土体的介电常数和电导率，测试时应保持探针匀速***，并改变探针***位置重复3~5次测量，取平均值；
其中，贯入仪包括荷载传感器和一根微型探针，探针总长为66mm，探头呈圆锥形且夹角为60°，探头上端直径为3.4mm，探杆直径为3.0mm，略小于探头，可消除贯入过程中探杆与侧壁的摩擦作用；FDR传感器包括主体和三根探针组成，探针长度为70mm，直径为3mm，间距13mm，可同时测得土体介电常数和电导率，两者的测试数据均记录在采集仪中；
步骤3：以土体的介电常数、电导率和最大贯入阻力作为响应目标，土体含水率和干密度作为影响因素，考虑线性效应、二次效应和交互效应，建立土体的介电常数、电导率和最大贯入阻力随质量含水率和干密度变化的二阶响应曲面模型：
Figure DEST_PATH_IMAGE072A
式中，
Figure DEST_PATH_IMAGE004AAAA
、
Figure DEST_PATH_IMAGE006AAAA
、
Figure DEST_PATH_IMAGE008AAAA
分别为土体的介电常数、电导率和最大贯入阻力的响应曲面函数；
Figure DEST_PATH_IMAGE010AA
、
Figure DEST_PATH_IMAGE012AA
、
Figure DEST_PATH_IMAGE014AA
为截断系数，
Figure DEST_PATH_IMAGE016AA
、
Figure DEST_PATH_IMAGE018AA
、
Figure DEST_PATH_IMAGE020AA
、
Figure DEST_PATH_IMAGE022AA
、
Figure DEST_PATH_IMAGE024AA
、
Figure DEST_PATH_IMAGE026AA
为线性项系数，
Figure DEST_PATH_IMAGE028AA
、
Figure DEST_PATH_IMAGE030AA
、
Figure DEST_PATH_IMAGE032AA
为交互项系数，
Figure DEST_PATH_IMAGE034AA
、
Figure DEST_PATH_IMAGE036AA
、
Figure DEST_PATH_IMAGE038AA
、
Figure DEST_PATH_IMAGE040AA
、
Figure DEST_PATH_IMAGE042AA
、
Figure DEST_PATH_IMAGE044AA
为二次项系数，均为通过最小二乘法拟合后获得的参数，由率定试验拟合确定；
Figure DEST_PATH_IMAGE046AA
为土体质量含水率，
Figure DEST_PATH_IMAGE048AA
为干密度，介电常数、电导率和最大贯入阻力的响应曲面图如图5、图6、图7所示；
步骤4：根据土体电导率和介电常数之间的相关性，建立土体电导率的修正模型：
Figure DEST_PATH_IMAGE050AA
式中，
Figure DEST_PATH_IMAGE052AA
,
Figure DEST_PATH_IMAGE054AA
,
Figure DEST_PATH_IMAGE056AA
,
Figure DEST_PATH_IMAGE058AA
为模型参数，
Figure DEST_PATH_IMAGE060AAAA
为电导率修正值，
Figure DEST_PATH_IMAGE062AA
为介电常数；
步骤5：将FDR传感器、贯入仪与电动反力升降装置连接，形成一套便携式快速检测装置，通过电机控制升降，将FDR传感器和贯入仪的探针***现场压实土体，分别测得现场压实土体的实测介电常数和实测最大贯入阻力，将实测介电常数代入电导率修正模型，计算得到现场压实土体的电导率修正值；
其中，FDR传感器和贯入仪与电动反力升降装置的传力杆底部采用可拆卸式的螺纹连接。现场测试时，先连接贯入仪测最大贯入阻力，再连接FDR传感器测介电常数和电导率，数值均记录在采集仪中，测试时，通过电机控制探针贯入速度；
步骤6：将实测介电常数、电导率修正值和实测最大贯入阻力作为已知量，代入对应的二阶响应曲面模型，以三个目标值的实测值和计算值均方根误差最小为原则，采用等权重法构造目标函数：
Figure DEST_PATH_IMAGE064AA
式中，
Figure DEST_PATH_IMAGE066AA
为实测介电常数，
Figure DEST_PATH_IMAGE060_5A
为电导率修正值，P为实测最大贯入阻力；
Figure DEST_PATH_IMAGE004_5A
、
Figure DEST_PATH_IMAGE006_5A
和
Figure DEST_PATH_IMAGE008_5A
为土体的介电常数、电导率和最大贯入阻力对应的响应曲面函数，然后，基于粒子群优化算法反演得到最优的一组计算参数，即为现场压实土体的质量含水率和干密度；
其中，步骤5、6中：现场数据采集完毕后，可在探针***部位进行环刀取样，结合环刀法和烘干法测得土体实际质量含水率和干密度，即可将利用本文提出的方法反演得到的质量含水率和干密度与环刀法结合烘干法测得的质量含水率和干密度进行对比；
对每块土样和土柱均进行两次检测，得到的测值数据如表1、表2所示，从测试结果看出，本发明基于响应曲面法的土体压实度检测方法与传统的环刀法和烘干法相比，对土体压实度测试结果较为精准，适用于不同的工程填筑土体。
表1
Figure DEST_PATH_IMAGE145
表2
Figure DEST_PATH_IMAGE147

Claims (5)

1.一种基于响应曲面法的土体压实度检测方法，其特征在于，包括以下步骤：
S1：取部分现场土体做室内标定实验，制备不同含水率和干密度的试样，利用FDR传感器测试土体的介电常数和电导率，利用贯入仪测试土体最大贯入阻力；
S2：以土体的介电常数、电导率和最大贯入阻力作为响应目标，土体含水率和干密度作为影响因素，考虑线性效应、二次效应和交互效应，建立土体的介电常数、电导率和最大贯入阻力随质量含水率和干密度变化的二阶响应曲面模型；
S3：根据土体的电导率和介电常数之间的相关性，建立土体电导率的修正模型；
S4：将FDR传感器、贯入仪与电动反力升降装置连接，形成一套便携式快速检测装置，通过电机控制升降，将FDR传感器和贯入仪的探针***压实土体，分别测得压实土体的实测介电常数和实测最大贯入阻力，将实测介电常数代入土体电导率的修正模型，计算得到土体的电导率修正值；
S5：将实测介电常数、电导率修正值和实测最大贯入阻力作为已知量，代入对应的二阶响应曲面模型，以三个目标值的实测值和模型计算值均方根误差最小为目标函数，基于粒子群优化算法反演得到最优的一组计算参数，即为现场压实土体的质量含水率和干密度。
2.根据权利要求1所述的一种基于响应曲面法的土体压实度检测方法，其特征在于，步骤S2中，土体的介电常数、电导率和最大贯入阻力随质量含水率和干密度变化的二阶响应曲面模型，其表达式如下：
Figure 590893DEST_PATH_IMAGE002
式中，
Figure 321083DEST_PATH_IMAGE004
、
Figure 163137DEST_PATH_IMAGE006
、
Figure 861141DEST_PATH_IMAGE008
分别为土体的介电常数、电导率和最大贯入阻力的响应曲面函数；
Figure 3410DEST_PATH_IMAGE010
、
Figure 107763DEST_PATH_IMAGE012
、
Figure 640376DEST_PATH_IMAGE014
为截断系数，
Figure 901593DEST_PATH_IMAGE016
、
Figure 350897DEST_PATH_IMAGE018
、
Figure 672157DEST_PATH_IMAGE020
、
Figure 692066DEST_PATH_IMAGE022
、
Figure 242127DEST_PATH_IMAGE024
、
Figure 296671DEST_PATH_IMAGE026
为线性项系数，
Figure 41029DEST_PATH_IMAGE028
、
Figure 813813DEST_PATH_IMAGE030
、
Figure 620095DEST_PATH_IMAGE032
为交互项系数，
Figure 342195DEST_PATH_IMAGE034
、
Figure 208520DEST_PATH_IMAGE036
、
Figure 999758DEST_PATH_IMAGE038
、
Figure 593419DEST_PATH_IMAGE040
、
Figure 419293DEST_PATH_IMAGE042
、
Figure 722098DEST_PATH_IMAGE044
为二次项系数，各系数均为模型最小二乘法拟合参数，由率定试验拟合确定；
Figure 689048DEST_PATH_IMAGE046
为土体质量含水率，
Figure 899450DEST_PATH_IMAGE048
为干密度。
3.根据权利要求1所述的一种基于响应曲面法的土体压实度检测方法，其特征在于，步骤S3中的土体电导率修正模型采用如下函数：
Figure 783092DEST_PATH_IMAGE050
式中，
Figure 40155DEST_PATH_IMAGE052
,
Figure 212510DEST_PATH_IMAGE054
,
Figure 961023DEST_PATH_IMAGE056
,
Figure 449904DEST_PATH_IMAGE058
为模型参数，
Figure 891250DEST_PATH_IMAGE060
为电导率修正值，
Figure 82060DEST_PATH_IMAGE062
为介电常数。
4.根据权利要求1所述的一种基于响应曲面法的土体干压实度预测模型及测试方法，其特征在于，步骤S5中，选取实测介电常数、电导率修正值和实测最大贯入阻力作为已知量，采用等权重法构造目标函数：
Figure 352373DEST_PATH_IMAGE064
式中，
Figure 679449DEST_PATH_IMAGE066
为实测介电常数，
Figure 494959DEST_PATH_IMAGE060
为电导率修正值，P为实测最大贯入阻力；
Figure DEST_PATH_IMAGE067
、
Figure DEST_PATH_IMAGE068
和
Figure 48431DEST_PATH_IMAGE008
为土体的介电常数、电导率和最大贯入阻力对应的响应曲面函数。
5.根据权利要求1所述的一种基于响应曲面法的土体压实度检测方法，其特征在于，步骤S4中，FDR传感器和贯入仪与电动反力升降装置的传力杆底部采用可拆卸式的螺纹连接，现场检测时，先连接贯入仪测最大贯入阻力，再连接FDR传感器测介电常数，数值均记录在采集仪中；测试时，通过电机控制探针贯入速度。
CN202211152176.XA 2022-09-21 2022-09-21 一种基于响应曲面法的土体压实度检测方法 Pending CN115455712A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211152176.XA CN115455712A (zh) 2022-09-21 2022-09-21 一种基于响应曲面法的土体压实度检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211152176.XA CN115455712A (zh) 2022-09-21 2022-09-21 一种基于响应曲面法的土体压实度检测方法

Publications (1)

Publication Number Publication Date
CN115455712A true CN115455712A (zh) 2022-12-09

Family

ID=84305257

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211152176.XA Pending CN115455712A (zh) 2022-09-21 2022-09-21 一种基于响应曲面法的土体压实度检测方法

Country Status (1)

Country Link
CN (1) CN115455712A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117074180A (zh) * 2023-10-17 2023-11-17 天津风霖物联网科技有限公司 一种建筑物地下土体压力变化测量方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117074180A (zh) * 2023-10-17 2023-11-17 天津风霖物联网科技有限公司 一种建筑物地下土体压力变化测量方法
CN117074180B (zh) * 2023-10-17 2023-12-12 天津风霖物联网科技有限公司 一种建筑物地下土体压力变化测量方法

Similar Documents

Publication Publication Date Title
Malicki et al. Application of time-domain reflectometry (TDR) soil moisture miniprobe for the determination of unsaturated soil water characteristics from undisturbed soil cores
US20040201385A1 (en) Method and apparatus for measuring dry density and water content of soil
CN104749205A (zh) 土体冻结过程水热力综合试验***及试验方法
CN108387619B (zh) 一种可控制应力状态的冻土未冻结含水量测试仪器
CN115455712A (zh) 一种基于响应曲面法的土体压实度检测方法
CN111650085A (zh) 一种便携式土体原位密度测试仪
CN110297017A (zh) 冻融循环下土体电参数时空分布室内测试***及工作方法
CN109470580A (zh) 一种评价我国不同海域粘性土强度参数的方法
CN113008935A (zh) 一种土壤含水量、热特性和电导率的原位测定装置及方法
Drnevich et al. Water content and density of soil insitu by the purdue TDR method
CN204536237U (zh) 一种土体冻结过程水热力综合试验***
CN109799140B (zh) 非饱和土强度特性的现场试验测定方法
CN106442603A (zh) 一种模拟不同热负荷下土体热物性参数的测试方法
CN211453441U (zh) 一种室内原状土样拉伸状态下电阻率测试试验装置
CN110455712B (zh) 一种倒置制样测量特殊曲面与土界面黏附力的制样装置与测试方法
CN216247627U (zh) 一种电阻率法测定原位土石混合料孔隙率的装置
CN212646386U (zh) 一种便携式土体原位密度测试仪
CN210626226U (zh) 软粘土不排水抗剪强度以及孔隙水压力的测试装置
Tian et al. Combined penetrometer and standing wave ratio probe to measure compactness and moisture content of soils.
CN2622713Y (zh) 交流低频土电阻率测试仪
CN212206982U (zh) 一种循环测定土体温湿度和贯入阻力的曲柄滑块装置
CN110411854B (zh) 软粘土不排水抗剪强度以及孔隙水压力的测试装置
CN202247834U (zh) 一种便携式自动贯入仪
CN110132679A (zh) 一种用于土工试验的试样制备器
CN210775028U (zh) 一种承载比测试仪器

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