1. 赛题回顾与核心挑战拆解2023年五一杯数学建模竞赛的C题题目是“深部岩石力学与安全开采”。这个题目一出来当时就在参赛圈里引起了不小的讨论。它不像一些纯数据预测或者优化类题目有明确的数据集和指标。这道题更像是一个“半开放”的综合应用题背景是煤矿深部开采中的巷道支护问题要求我们建立数学模型来分析围岩的稳定性并设计支护方案。很多同学拿到题的第一感觉是“懵”因为涉及到的知识面很广从岩石力学、弹塑性理论到有限元数值模拟再到多目标决策感觉每个部分都能写一篇论文。但恰恰是这种“综合性”才是数学建模竞赛考察的核心——如何将一个复杂的工程问题抽象、简化为可计算的数学模型并给出有指导意义的结论。这道题的核心价值在于它逼着参赛者去进行“跨学科”的思考。你不需要是岩石力学专家但你需要快速理解“地应力”、“支护阻力”、“塑性区”这些关键概念并用数学语言描述它们之间的关系。最终的模型和方案其“物理意义”和“工程可行性”往往比复杂的数学公式堆砌更重要。评委想看的是你解决问题的逻辑链条是否完整从问题分析、假设提出、模型建立、求解到结果分析每一步是否清晰合理。所以这篇解析不会给你一个“标准答案”事实上也不存在而是带你复盘一个高完成度的解题思路框架重点讲清楚“为什么这么想”以及“如何落地”希望能为你应对类似综合性赛题提供一套可迁移的方法论。2. 问题一围岩应力分布与塑性区半径模型问题一通常是整个赛题的基石要求建立巷道围岩的应力分布模型并计算塑性区范围。这是典型的岩石力学中的“圆形巷道”问题在很多教材和文献中都有经典解。但直接套公式是不够的关键在于理解公式的推导逻辑和适用条件。2.1 模型选择与假设的合理性首先面对一个实际的地下巷道我们必须进行简化。最经典的简化就是将其视为无限大岩体中的一个圆形孔洞。这里就引入了第一个关键假设围岩是连续、均质、各向同性的线弹性-理想塑性材料。这个假设是后续所有经典解析解的基础。为什么是“线弹性-理想塑性”岩石在受力初期应力应变关系近似线性弹性阶段当应力超过其屈服强度后会发生不可逆的塑性变形但假设其承载力不再增加理想塑性。这比纯弹性模型更能反映岩石破坏的实际情况又比复杂的硬化/软化模型更易于求解。“各向同性”意味着什么意味着岩石在各个方向上的力学性质相同。这显然是对实际地层的一种简化因为岩层往往有层理、节理。但在初步分析中这个假设可以大大降低模型复杂度。在论文中需要明确指出这一点并讨论其局限性。基于以上假设我们就可以引用著名的Kastner方程或Fenner公式来描述塑性区半径。其核心思想是在巷道周边径向应力最小切向应力集中。当切向应力超过岩石的抗压强度时该处岩石进入塑性状态。从巷道壁面向岩体深处应力状态逐渐恢复存在一个弹性区和塑性区的分界这个分界的半径就是塑性区半径 ( R_p )。其经典表达式为 [ R_p a \cdot \left[ \frac{(P_0 c \cdot \cot \varphi)(1 - \sin \varphi)}{P_i c \cdot \cot \varphi} \right]^{\frac{1 - \sin \varphi}{2 \sin \varphi}} ] 其中( a )巷道半径已知( P_0 )原岩应力关键参数通常等于上覆岩层重量 ( \gamma H )( P_i )支护阻力问题二的核心此处可先设为0或一个较小值进行分析( c )岩石粘聚力( \varphi )岩石内摩擦角2.2 参数获取与敏感性分析公式有了但参数 ( c ) 和 ( \varphi ) 从哪里来题目通常不会直接给。这里就体现了数学建模中“数据获取与处理”的能力。有几种思路查阅文献根据题目描述的“深部”、“煤矿”、“砂质泥岩”等关键词去中英文数据库如知网、Web of Science检索类似地质条件下的岩石力学参数实验值给出一个合理的取值范围。反向标定如果题目给出了某些间接信息比如巷道变形量可以尝试建立参数与观测值之间的关系进行反演分析。这在论文中是加分项。参数敏感性分析这是必须做的一步。即使你只能给出参数的估计范围也要系统地分析 ( c )、( \varphi )、( P_0 ) 的变化会如何影响塑性区半径 ( R_p )。你可以设计一个正交实验或简单的单变量分析。注意在计算原岩应力 ( P_0 \gamma H ) 时容重 ( \gamma ) 的取值也要有依据。煤矿岩石的容重大约在25 kN/m³左右。深度 ( H ) 是题目给定的核心变量之一。通过敏感性分析你可以得出一些有工程指导意义的结论例如“内摩擦角 ( \varphi ) 对塑性区范围的影响最为敏感提高岩石的内摩擦角例如通过注浆加固比提高粘聚力更能有效控制塑性区发展。” 这样的结论就比单纯算出一个 ( R_p ) 值更有价值。3. 问题二支护方案设计与多目标优化模型问题一告诉我们围岩在哪里、有多大范围可能失稳。问题二则要求我们主动干预设计支护方案主要是确定支护阻力 ( P_i )。这本质上是一个优化问题在一定的成本和安全约束下寻找最佳的支护参数。3.1 目标函数与决策变量的定义首先明确我们要优化什么。对于巷道支护核心目标通常有两个安全性最大化表现为塑性区半径 ( R_p ) 最小化或者围岩变形量 ( u ) 最小化。经济性最优化表现为支护成本 ( C ) 最小化。支护成本与支护阻力 ( P_i ) 正相关阻力越大需要的支护材料越多、强度越高、施工越复杂。决策变量很明确就是支护阻力 ( P_i )。但 ( P_i ) 不是一个孤立的数值它对应着具体的支护方式和参数。例如采用“锚杆锚索喷射混凝土”联合支护时锚杆的间排距、长度、直径锚索的预紧力喷层的厚度等共同决定了最终作用于围岩的等效支护阻力 ( P_i )。你需要建立一个从支护设计参数到等效 ( P_i ) 的换算模型。这可能需要查阅《煤矿巷道支护技术规范》或相关研究论文中的经验公式。因此目标函数可以初步构建为目标1安全Minimize ( R_p(P_i) ) 或 Minimize ( u(P_i) ) ( u ) 可通过弹性区或塑性区理论计算位移公式求得目标2经济Minimize ( C(P_i) )3.2 约束条件的建立优化不能天马行空必须满足实际约束技术约束支护阻力 ( P_i ) 有上下限。下限通常大于0至少要抵消部分地压上限受限于当前支护技术的最大能力如锚杆的破断力、混凝土的强度。安全约束塑性区半径 ( R_p ) 必须小于某个临界值 ( [R_{p}^{cr}] )或者巷道顶底板移近量必须小于安全规程允许值。这个临界值需要根据巷道用途回采巷道、准备巷道来设定。几何约束支护参数如锚杆长度不能超过巷道尺寸。3.3 多目标优化求解策略面对“安全”和“经济”这两个相互冲突的目标我们处理的是一个多目标优化问题。常用解法有权重系数法将两个目标通过权重 ( w_1 )、( w_2 ) 线性加权为一个总目标Minimize ( w_1 * R_p(P_i) w_2 * C(P_i) )。难点在于权重的确定具有主观性。可以在论文中设置几组不同的权重如更侧重安全 ( w_10.8, w_20.2 )均衡考虑 ( w_10.5, w_20.5 )对比分析结果并说明权重的工程意义。约束法将一个目标转化为约束。例如以经济成本最小化为目标同时约束塑性区半径 ( R_p \leq [R_{p}^{cr}] )。这样更符合工程思维——“在满足安全底线的前提下寻找最省钱的方案”。帕累托前沿Pareto Front求解这是更高级、也更受评委青睐的方法。通过算法如NSGA-II多目标遗传算法求解出一系列非支配解这些解构成了帕累托前沿。在这个前沿上任何一个安全性的提升都必须以经济性的下降为代价。决策者可以根据实际情况在这个前沿上选择最合适的方案。在论文中如果你采用了帕累托前沿一定要将求解结果可视化散点图并对其上的几个典型解进行解释“A点代表高安全储备方案成本较高适用于服务年限长的主巷道B点代表经济型方案安全性在可接受范围内适用于服务年限短的采准巷道。”4. 问题三基于可靠度理论的稳定性综合评价问题三往往是赛题的升华要求评价巷道围岩的“稳定状态”。如果仅仅用问题一算出的一个确定的 ( R_p ) 值来判断“稳定”或“不稳定”就显得过于武断和理想化。因为我们的模型输入参数( c, \varphi, P_0 ) 等本身存在不确定性变异性。因此引入可靠度理论或敏感性分析来评价稳定性是体现模型深度和实用性的关键。4.1 从确定性分析到不确定性分析在问题一中我们假设参数是确定的得到的是一个确定的 ( R_p )。但实际上岩石参数来自实验或经验存在统计分布特征。例如粘聚力 ( c ) 可能服从正态分布 ( N(\mu_c, \sigma_c^2) )。那么计算出的 ( R_p ) 也将是一个随机变量。4.2 建立极限状态函数可靠度分析的核心是定义“失效”的界限。对于巷道我们可以定义极限状态函数 ( Z ) [ Z R_{p}^{cr} - R_p(c, \varphi, P_0, ...) ] 其中 ( R_{p}^{cr} ) 是前面提到的临界塑性区半径可视为确定值。那么当 ( Z 0 ) 时巷道处于安全状态。当 ( Z 0 ) 时巷道处于极限平衡状态。当 ( Z 0 ) 时巷道失效塑性区过大。由于 ( R_p ) 是随机变量的函数因此 ( Z ) 也是一个随机变量。4.3 可靠度指标 ( \beta ) 的计算与解释巷道稳定的可靠度 ( P_s ) 定义为 ( P(Z 0) )。直接求这个概率比较困难通常引入可靠度指标 ( \beta )。最常用的是一次二阶矩法FORM。将非正态分布的随机变量当量正态化。将极限状态函数在验算点最可能失效点线性化。计算 ( \beta \mu_Z / \sigma_Z )其中 ( \mu_Z ) 和 ( \sigma_Z ) 是 ( Z ) 的均值和标准差。( \beta ) 的物理意义非常直观它表示标准正态空间下均值点到失效边界( Z0 )的距离。( \beta ) 值越大系统越可靠。工程上常有对应关系例如 ( \beta 3.0 ) 对应失效概率 ( P_f \approx 0.0013 )。在论文中你需要合理假设关键参数如 ( c, \varphi )的概率分布类型和变异系数标准差/均值。通过蒙特卡洛模拟Monte Carlo Simulation或FORM法计算不同支护方案不同 ( P_i ) 下的可靠度指标 ( \beta )。对比分析展示随着支护阻力 ( P_i ) 增加( \beta ) 如何变化。这实际上是将问题二的“确定性优化”推进到了“可靠性优化”的层面。你可以提出一个新的优化模型在满足目标可靠度 ( \beta \geq [\beta_{target}] ) 的前提下最小化支护成本 ( C )。4.4 稳定性等级划分最后基于可靠度指标 ( \beta )你可以建立一个简单的稳定性等级评价体系。例如( \beta \geq 4.0 )稳定。塑性区发展可能性极低常规支护即可。( 3.0 \leq \beta 4.0 )基本稳定。存在一定风险需加强监测和支护。( 2.0 \leq \beta 3.0 )稳定性较差。塑性区可能持续发展需采取强力联合支护措施。( \beta 2.0 )不稳定。有较大失稳风险需重新设计支护方案或采取特殊工法。这样的评价结果比单纯说“塑性区半径是3.5米”要丰富和深刻得多它给出了一个考虑不确定性的、概率意义上的安全评价。5. 模型求解、论文写作与避坑指南有了清晰的思路接下来就是具体的求解实现和论文撰写。这部分往往决定了论文的最终档次。5.1 求解工具与实现要点解析计算问题一的公式、问题二的简单优化完全可以用MATLAB或PythonNumPy, SciPy快速实现。重点在于编写清晰的函数并做好参数的全局管理方便后续调用和敏感性分析。数值模拟如果想做得更出彩可以引入数值方法。例如使用COMSOL Multiphysics或ABAQUS建立巷道围岩的二维轴对称有限元模型直接模拟开挖和支护过程得到应力、位移云图并与解析解相互验证。这能极大提升论文的“颜值”和说服力但对软件操作和计算资源要求较高。优化算法对于问题二的多目标优化如果采用智能算法求帕累托前沿MATLAB的Global Optimization Toolbox或Python的PyGMO、DEAP库都内置了NSGA-II算法可以直接调用。关键是要正确编写目标函数和约束函数。5.2 论文写作的核心结构一篇好的数模论文结构比文采更重要。摘要重中之重用一段话精炼概括每个问题用了什么方法、建立了什么模型、得到了什么关键结论。避免出现公式和细节但要体现模型名称如“基于Kastner方程的塑性区模型”、“多目标优化模型”、“可靠度评价模型”和核心结果如“得到了支护阻力与成本的帕累托前沿”、“当支护阻力达到X MPa时可靠度指标可提升至X”。问题重述与分析不要照抄题目要用自己的话梳理问题的背景、条件和待解决的点并简要分析解题思路。这部分展现你对题目的理解。模型假设与符号说明假设要合理且必要符号表格要清晰。模型的建立与求解这是主体。严格按照问题一、二、三来组织。每一部分都应包含分析 - 模型建立公式推导- 参数说明 - 求解过程可附关键代码截图或流程图- 结果分析与讨论。讨论部分尤其重要要解释你的结果说明了什么有何工程意义。模型的评价与推广客观评价自己模型的优点如引入可靠度理论使评价更科学和缺点如假设岩石各向同性未考虑节理影响。提出可能的改进方向如考虑流变、动态扰动等。参考文献规范引用特别是岩石力学经典教材、支护规范以及你所用算法或公式的出处。5.3 常见“坑”与应对策略坑一盲目追求模型复杂度。有的队伍一上来就想用复杂的本构模型如D-P准则、Hoek-Brown准则和三维有限元。结果时间都花在建模和调试上忽略了最核心的分析和优化。策略先基于经典弹性-理想塑性模型把整个解题逻辑跑通得到完整的结果。如果有余力再考虑用更复杂的模型进行对比作为模型的改进部分写在论文里。坑二参数处理不当。直接给参数赋值不做任何来源说明和敏感性分析。策略务必说明关键参数的取值依据文献、规范或合理假设并必须进行敏感性分析用图表展示结果这能体现模型的鲁棒性。坑三优化问题求解草率。对于多目标优化只是简单试算几个点没有系统性地求解帕累托前沿。策略即使只用权重法也要系统性地变化权重展示结果谱系。如果能用算法求出前沿哪怕参数设置简单些也是巨大的亮点。坑四论文像实验报告。通篇都是“我们做了A然后做了B结果如图C”缺乏分析、讨论和洞察。策略在每一个结果图、表后面紧跟一段“结果分析”解释这个图/表揭示了什么规律为什么会出现这样的规律其物理或工程意义是什么。坑五忽视可视化。大段的文字和公式堆砌。策略多使用图表。例如塑性区半径随各参数变化的曲线图、多目标优化的帕累托前沿散点图、可靠度指标随支护阻力变化的关系图、敏感性分析的雷达图或柱状图。一图胜千言。这道“深部岩石力学与安全开采”的题目本质上是一个“物理机理建模 参数优化决策 不确定性评价”的经典数模框架。掌握这个框架并能够根据具体问题填充合适的学科知识如本题的岩石力学就能应对很多类似的综合性赛题。在实际比赛中团队分工要明确一人主攻模型推导与建立一人主攻算法实现与求解一人主攻论文写作与整合。保持沟通确保思路一致最终呈现出一份逻辑严密、论据充分、表达清晰的作品。