一种考虑含石量的土石混合料本构模型建立方法技术

一种考虑含石量的土石混合料本构模型建立方法。其包括建立含石量与粘聚力及内摩擦角的关系式；基于扩展SMP准则推导粘聚力不为0时的变换应力；建立三维化变换应力空间中考虑含石量土石混合料本构模型的屈服函数；计算塑性应变增量；计算塑性标量因子；推导弹塑性刚度矩阵；推导对称的弹塑性刚度矩阵并确定出考虑含石量的土石混合料本构模型等步骤。本发明专利技术在本构模型中考虑含石量，将含石量和粘聚力及内摩擦角关系引入进来，建立了考虑含石量的土石混合料本构模型，能够更加合理地反映土石混合料的抗剪强度和应力应变特性，使得对实际高填方工程问题的分析更加合理。所建立的土石混合料本构模型参数简单易测，便于在有限元中应用。

【技术实现步骤摘要】
一种考虑含石量的土石混合料本构模型建立方法
本专利技术属于岩土本构模型
，特别涉及一种考虑含石量的土石混合料本构模型建立方法。
技术介绍
随着航空业的发展，机场建设规模越来越大，尤其是西部机场建设也发展迅速，西部机场建设多为高填方工程，其特点为填方量大，填方高度高，而且填料多为土石混合料，土石混合料的组成主要以砾石、块石与砂土、粘土等为主，与一般土体相比，其力学性质复杂。在机场高填方工程建设中，工程的工后沉降会引起许多沉降病害和稳定性问题，因此机场高填方工程的沉降变形及稳定性问题成为了机场建设的核心问题。目前，岩土数值分析是岩土工程问题的常用分析方法，而本构模型是数值分析的核心技术，能否对机场高填方工程变形及稳定性进行合理分析，关键取决于对高填方工程土石混合料的力学特性的准确模拟，也就是需要建立一种能准确反映土石混合料应力应变特性的本构模型，这样就能够对机场高填方工程问题进行合理的分析。目前在数值分析中广泛应用的MCC等弹塑性模型，主要针对的是细粒土，只能反映正常固结土或弱超固结土的应力应变特性，对超固结土或土石混合料的应力应变特性还不能够合理反映。所以迫切需要建立一种针对土石混合料的、能合理反映土石混合料应力应变特性的本构模型，以便该模型能方便应用于数值分析中，从而能够对机场高填方工程地基变形及稳定性进行合理分析，为实际工程施工提供指导依据。含石量是影响土石混合料抗剪强度的重要因素，影响着内摩擦角和粘聚力，所以如果在本构模型中考虑含石量，将含石量和粘聚力及内摩擦角的关系引入进来，则能够更加合理地反映土石混合料的抗剪强度和应力应变特性，使得对实际高填方工程问题的分析更加合理。
技术实现思路
为了解决上述问题，本专利技术的目的在于提供一种考虑含石量的土石混合料本构模型建立方法。为了达到上述目的，本专利技术提供的考虑含石量的土石混合料本构模型建立方法包括按顺序进行的下列步骤：1)根据土石混合料的含石量与粘聚力及内摩擦角增量的关系图，获得土石混合料的含石量与内摩擦角关系式以及含石量与土的粘聚力关系式，由此得到土石混合料的内摩擦角及土的粘聚力；2)根据土的强度准则中扩展松冈-中井准则，并参考统一硬化模型中变换应力公式，结合步骤1)得到的土石混合料的内摩擦角及土的粘聚力，获得土的粘聚力不为0时的变换应力；3)基于上述获得的土的粘聚力不为0时的变换应力，考虑含石量和内摩擦角及土的粘聚力的关系，参考统一硬化模型的屈服函数公式，在三维化的变换应力空间中，得到考虑含石量的土石混合料本构模型的屈服函数；4)利用塑性应变增量表达式计算出塑性应变增量；5)对上述步骤3)获得的屈服函数求微分并代入步骤4)获得的塑性应变增量，得到塑性标量因子；6)将步骤5)获得的塑性标量因子代入步骤4)中的塑性应变增量表达式，然后再代入应力增量和应变增量的关系式，得到弹塑性刚度矩阵；7)对上述步骤6)获得的弹塑性刚度矩阵进行对称化处理，得到对称的弹塑性刚度矩阵，最后基于对称的弹塑性刚度矩阵确定出考虑含石量的土石混合料本构模型。在步骤1)中，所述根据土石混合料的含石量与粘聚力及内摩擦角增量的关系图，获得土石混合料的含石量与内摩擦角关系式以及含石量与土的粘聚力关系式，由此得到土石混合料的内摩擦角及土的粘聚力的具体方法如下：土石混合料的内摩擦角与含石量的关系式为：其中，式中，为纯土的内摩擦角；为土石混合料的内摩擦角；PR为土石混合料的含石量；为当含石量为PR时土石混合料的内摩擦角较相应纯土的内摩擦角增量；为当含石量PR超过70％时土石混合料的内摩擦角较相应纯土的内摩擦角增量；土的粘聚力与含石量的关系式为：c＝ccg(PR)(3)其中，cc为纯土的粘聚力；c为土的粘聚力。在步骤2)中，所述根据土的强度准则中扩展松冈-中井准则，并参考统一硬化模型中变换应力公式，结合步骤1)得到的土石混合料的内摩擦角及土的粘聚力，获得土的粘聚力不为0时的变换应力的具体方法如下：土的粘聚力不为0时的变换应力为：其中，其中，σ为正应力；为土的粘聚力c＝0情况下用扩展SMP准则表示的变换应力；为土的粘聚力c≠0情况下用扩展SMP准则表示的平均主应力；为土的粘聚力c≠0情况下用扩展SMP准则表示的三轴压缩时的偏应力；为土的粘聚力c≠0情况下用扩展SMP准则表示的应力不变量；为扩展SMP面上的三个主应力。在步骤3)中，所述基于上述获得的土的粘聚力不为0时的变换应力，考虑含石量和内摩擦角及土的粘聚力的关系，参考统一硬化模型的屈服函数公式，在三维化的变换应力空间中，得到考虑含石量的土石混合料本构模型的屈服函数的具体方法如下：考虑含石量的土石混合料本构模型的屈服函数为：式中，为变换应力空间当前应力点上的平均主应力；为该点上的广义剪应力；为相应的初始平均主应力；cp＝(λ-κ)/1+e0，λ为等向压缩线的斜率，κ为等向回弹线的斜率，e0为初始孔隙比；M为特征状态和临界状态下的应力比，为硬化参数，可表示为：式中，为变换应力空间的潜在强度，为变换应力空间的超固结参数，为前期固结压力，为塑性体积应变，为变换应力空间的应力比。在步骤4)中，所述利用塑性应变增量表达式计算出塑性应变增量的具体方法如下：根据塑性位势理论获得塑性应变增量的表达式为：式中，dεp为塑性应变增量；Λ为塑性标量因子。在步骤5)中，所述对上述步骤3)获得的屈服函数求微分并代入步骤4)获得的塑性应变增量，得到塑性标量因子的具体方法如下：塑性标量因子为：式中，在步骤6)中，所述将步骤5)获得的塑性标量因子代入步骤4)中的塑性应变增量表达式，然后再代入应力增量和应变增量的关系式，得到弹塑性刚度矩阵的具体方法如下：将公式(15)代入式(14)，得到：又知，应力增量和应变增量的关系式为：dσ＝De(dε-dεp)(19)将公式(18)代入上述应力增量和应变增量的关系式(19)得到：式(20)可以写为：dσ＝Depdε(21)从而得到弹塑性刚度矩阵Dep为：在步骤7)中，所述对上述步骤6)获得的弹塑性刚度矩阵进行对称化处理，得到对称的弹塑性刚度矩阵，最后基于对称的弹塑性刚度矩阵确定出考虑含石量的土石混合料本构模型的具体方法如下：对弹塑性刚度矩阵进行对称化处理，得到对称的弹塑性刚度矩阵为：其中，基于对称的弹塑性刚度矩阵获得考虑含石量的土石混合料本构模型为：dσ＝Depdε(25)。本专利技术的效果在于：1)在土石混合料本构模型中考虑含石量，将含石量和土的粘聚力及内摩擦角的关系引入进来，则能够更加合理地反映土石混合料本文档来自技高网...

【技术保护点】
1.一种考虑含石量的土石混合料本构模型建立方法，其特征在于：所述考虑含石量的土石混合料本构模型建立方法包括按顺序进行的下列步骤：/n1)根据土石混合料的含石量与粘聚力及内摩擦角增量的关系图，获得土石混合料的含石量与内摩擦角关系式以及含石量与土的粘聚力关系式，由此得到土石混合料的内摩擦角及土的粘聚力；/n2)根据土的强度准则中扩展松冈-中井准则，并参考统一硬化模型中变换应力公式，结合步骤1)得到的土石混合料的内摩擦角及土的粘聚力，获得土的粘聚力不为0时的变换应力；/n3)基于上述获得的土的粘聚力不为0时的变换应力，考虑含石量和内摩擦角及土的粘聚力的关系，参考统一硬化模型的屈服函数公式，在三维化的变换应力空间中，得到考虑含石量的土石混合料本构模型的屈服函数；/n4)利用塑性应变增量表达式计算出塑性应变增量；/n5)对上述步骤3)获得的屈服函数求微分并代入步骤4)获得的塑性应变增量，得到塑性标量因子；/n6)将步骤5)获得的塑性标量因子代入步骤4)中的塑性应变增量表达式，然后再代入应力增量和应变增量的关系式，得到弹塑性刚度矩阵；/n7)对上述步骤6)获得的弹塑性刚度矩阵进行对称化处理，得到对称的弹塑性刚度矩阵，最后基于对称的弹塑性刚度矩阵确定出考虑含石量的土石混合料本构模型。/n...

【技术特征摘要】
1.一种考虑含石量的土石混合料本构模型建立方法，其特征在于：所述考虑含石量的土石混合料本构模型建立方法包括按顺序进行的下列步骤：
1)根据土石混合料的含石量与粘聚力及内摩擦角增量的关系图，获得土石混合料的含石量与内摩擦角关系式以及含石量与土的粘聚力关系式，由此得到土石混合料的内摩擦角及土的粘聚力；
2)根据土的强度准则中扩展松冈-中井准则，并参考统一硬化模型中变换应力公式，结合步骤1)得到的土石混合料的内摩擦角及土的粘聚力，获得土的粘聚力不为0时的变换应力；
3)基于上述获得的土的粘聚力不为0时的变换应力，考虑含石量和内摩擦角及土的粘聚力的关系，参考统一硬化模型的屈服函数公式，在三维化的变换应力空间中，得到考虑含石量的土石混合料本构模型的屈服函数；
4)利用塑性应变增量表达式计算出塑性应变增量；
5)对上述步骤3)获得的屈服函数求微分并代入步骤4)获得的塑性应变增量，得到塑性标量因子；
6)将步骤5)获得的塑性标量因子代入步骤4)中的塑性应变增量表达式，然后再代入应力增量和应变增量的关系式，得到弹塑性刚度矩阵；
7)对上述步骤6)获得的弹塑性刚度矩阵进行对称化处理，得到对称的弹塑性刚度矩阵，最后基于对称的弹塑性刚度矩阵确定出考虑含石量的土石混合料本构模型。


2.根据权利要求1所述考虑含石量的土石混合料本构模型建立方法，其特征在于：在步骤1)中，所述根据土石混合料的含石量与粘聚力及内摩擦角增量的关系图，获得土石混合料的含石量与内摩擦角关系式以及含石量与土的粘聚力关系式，由此得到土石混合料的内摩擦角及土的粘聚力的具体方法如下：
土石混合料的内摩擦角与含石量的关系式为：



其中，



式中，为纯土的内摩擦角；为土石混合料的内摩擦角；PR为土石混合料的含石量；为当含石量为PR时土石混合料的内摩擦角较相应纯土的内摩擦角增量；为当含石量PR超过70％时土石混合料的内摩擦角较相应纯土的内摩擦角增量；
土的粘聚力与含石量的关系式为：
c＝ccg(PR)(3)



其中，cc为纯土的粘聚力；c为土的粘聚力。


3.根据权利要求1所述考虑含石量的土石混合料本构模型建立方法，其特征在于：在步骤2)中，所述根据土的强度准则中扩展松冈-中井准则，并参考统一硬化模型中变换应力公式，结合步骤1)得到的土石混合料的内摩擦角及土的粘聚力，获得土的粘聚力不为0时的变换应力的具体方法如下：
土的粘聚力不为0时的变换应力为：



其中，















其中，σ为正应力；为土的粘聚力c＝0情况下用扩展SMP准则表示的变换应力；为土的粘聚力c≠0情况下用扩展SMP准则表示的平均主应力；为土的粘聚力c≠0情况下用扩展SMP准则表示的三轴压缩时的偏应力；为土的粘聚力c≠0情况下用扩展SMP准则表...

【专利技术属性】
技术研发人员：冯兴，张献民，
申请(专利权)人：中国民航大学，
类型：发明
国别省市：天津;12