针对CycLiq本构模型的自动校准方法，重点研究其液化行为

《Soil Dynamics and Earthquake Engineering》：Automatic calibration method for the CycLiq constitutive model with focus on liquefaction behaviour

自动参数校准方法ACM4CycLiq通过三阶段优化和Hausdorff距离目标函数实现CycLiq模型的客观校准，验证其在合成数据和离心机试验中的有效性，CMA-ES算法表现最佳，参数校准精度达98%。

　　
土壤本构模型参数标定是岩土工程领域的关键技术挑战，尤其在模拟饱和砂土循环液化行为时面临诸多难题。当前研究主要存在两大痛点：其一，传统手动标定方法依赖专家经验且效率低下，难以应对模型中参数数量增加带来的复杂性；其二，现有自动标定方法在算法选择、目标函数构建及多阶段标定流程等方面存在不足，导致参数敏感性差异显著，模型预测结果离散性较大。针对上述问题，中国建筑第八工程局有限公司研究团队在CycLiq本构模型标定领域取得突破性进展，提出了具有多阶段优化架构和新型目标函数的自动化标定方法ACM4CycLiq。

研究首先系统梳理了先进土壤本构模型的发展脉络。早期简单模型如Cam-Clay模型仅有少数物理意义明确的参数，可通过解析计算或单组试验数据直接标定。随着模型复杂度提升，如采用多曲面塑性理论的本构模型，参数数量激增且存在多重耦合关系，传统标定方法已无法满足需求。LEAP项目2022年的最新数据显示，即便使用相同本构模型和试验数据，不同研究者选择的参数组合会导致离心机试验模拟结果差异超过30%，这暴露了参数标定过程中主观因素影响过大的根本缺陷。

针对CycLiq模型这种具有明确物理意义的先进本构模型，研究团队创新性地构建了三阶段递进式标定框架。第一阶段采用直接回归法确定7个基础参数，这些参数具有明确的物理对应关系，如剪切模量、体积模量等，通过特征点匹配即可快速获取。第二阶段聚焦弹性与塑性剪切模量参数的优化，这些参数直接影响应力-应变曲线的形状特征。第三阶段则针对膨胀特性相关的4个关键参数进行标定，这部分参数对循环荷载下的孔隙水压力演化具有决定性影响。

在算法选择方面，研究团队对CMA-ES、GA、PSO、DE等主流优化算法进行了系统性对比。通过构建包含压缩、剪切、三轴循环等典型应力路径的合成数据集，评估不同算法在参数空间维度（最高达11维）、目标函数收敛速度（平均迭代次数）、标定精度（平均相对误差）等关键指标上的表现。实验结果表明，CMA-ES算法在参数维度超过8维时展现出显著优势，其收敛稳定性比传统启发式算法提升约40%，参数标定误差控制在3%以内，且对初始猜测不敏感的特性使其特别适合处理本构模型中参数间的非线性耦合问题。

实际应用验证阶段，研究团队选取Ottawa F-65砂作为基准材料。通过LEAP项目提供的离心机模型试验数据（包含11组不同工况测试），对比传统手动标定与ACM4CycLiq自动标定方法的预测精度。结果显示，在循环荷载作用下，自动标定模型对挡土墙位移的模拟平均误差仅为2.1%，最大误差不超过5%，且各工况下的预测离散系数低于12%。特别值得注意的是，在模拟7次往复循环的复合应力路径时，自动标定模型成功捕捉到膨胀参数随应力比值的非线性变化规律，而传统方法因参数数量多易出现过度拟合现象。

参数敏感性分析采用Sobol全局敏感性方法，量化了膨胀相关参数对预测结果的影响程度。研究发现，临界膨胀压力指数（λc）和膨胀模量（Kd）对循环孔隙水压力演化的敏感度最高（一阶指数贡献率超过65%），而静水压力系数（λp）和膨胀残余模量（Kdr）对长期累积变形的敏感性更为显著（方差贡献率约58%）。这种敏感性差异特征指导了研究团队在算法设计时采用分阶段优化的策略，既保证了基础参数的物理意义明确性，又通过分阶段优化有效降低了计算复杂度。

该方法的创新性体现在三个方面：首先，构建了包含应力比、累积塑性应变、孔隙水压力比等12项指标的复合目标函数，通过Hausdorff距离量化模拟曲线与试验数据的整体匹配度，相比传统单一指标（如R2系数）更能有效捕捉非平稳荷载下的动态响应特征。其次，采用参数解耦策略，将11个待标定参数划分为直接回归、弹性-塑性耦合优化和膨胀特性优化三个层次，既保留了专家经验的价值，又通过算法优化解决了高维参数空间搜索难题。最后，通过建立算法性能评价指标体系（涵盖收敛速度、鲁棒性、精度保持性等维度），实现了不同优化算法的横向比较，为后续研究提供了可复用的算法评估框架。

在工程应用层面，该方法显著提升了本构模型的应用门槛。以LEAP项目中的离心机挡土墙模型试验为例，传统标定流程需要3-5名工程师配合完成，平均耗时120小时。采用ACM4CycLiq后，仅需输入基础物性参数和试验数据，系统可在72小时内完成全参数标定，且不同工程师使用相同算法得到的结果差异系数从传统方法的38%降至7%以下。这种标准化流程使得本构模型在工程实践中的推广更具可行性，特别是在多团队协作的复杂项目中，可减少约70%的参数标定沟通成本。

研究团队通过构建标准化数据接口和算法工具箱，将ACM4CycLiq方法封装成模块化软件系统。该系统支持常见格式（如MAT、CSV、HDF5）的数据输入，提供可视化参数敏感性热力图和优化过程追踪功能。在LEAP 2022验证阶段，该系统成功处理了来自8个不同实验室的、包含327组试验数据的标定任务，平均处理时间从传统方法的45小时缩短至6.8小时，计算资源消耗降低约60%。

未来研究将重点拓展该方法的应用范围：在软土液化评估方面，计划将现有标定流程扩展至非饱和土体；在数值模拟优化方面，拟集成代理模型加速高维参数搜索；在工程应用方面，将开发基于移动端的前端工具，实现现场试验数据即时标定。该方法论的创新为先进本构模型的工程化应用提供了可复制的技术路径，特别在大型基础设施（如跨海通道、超长挡墙）的地震液化评估中展现出重要应用价值。研究团队正在与多个国际工程公司合作，将ACM4CycLiq集成到BIM平台和地震动响应分析系统中，预计可使复杂土工结构的液化评估效率提升3-5倍。