CN113297665A - 一种考虑含石量的土石混合料本构模型建立方法 - Google Patents

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

Info

Publication number
CN113297665A
CN113297665A CN202110735907.2A CN202110735907A CN113297665A CN 113297665 A CN113297665 A CN 113297665A CN 202110735907 A CN202110735907 A CN 202110735907A CN 113297665 A CN113297665 A CN 113297665A
Authority
CN
China
Prior art keywords
soil
stone
stress
content
stone content
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
CN202110735907.2A
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.)
Civil Aviation University of China
Original Assignee
Civil Aviation University of China
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 Civil Aviation University of China filed Critical Civil Aviation University of China
Priority to CN202110735907.2A priority Critical patent/CN113297665A/zh
Publication of CN113297665A publication Critical patent/CN113297665A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Computing Systems (AREA)
  • Evolutionary Computation (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

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

Description

一种考虑含石量的土石混合料本构模型建立方法
技术领域
本发明属于岩土本构模型技术领域,特别涉及一种考虑含石量的土石混合料本构模型建立方法。
背景技术
随着航空业的发展,机场建设规模越来越大,尤其是西部机场建设也发展迅速,西部机场建设多为高填方工程,其特点为填方量大,填方高度高,而且填料多为土石混合料,土石混合料的组成主要以砾石、块石与砂土、粘土等为主,与一般土体相比,其力学性质复杂。
在机场高填方工程建设中,工程的工后沉降会引起许多沉降病害和稳定性问题,因此机场高填方工程的沉降变形及稳定性问题成为了机场建设的核心问题。目前,岩土数值分析是岩土工程问题的常用分析方法,而本构模型是数值分析的核心技术,能否对机场高填方工程变形及稳定性进行合理分析,关键取决于对高填方工程土石混合料的力学特性的准确模拟,也就是需要建立一种能准确反映土石混合料应力应变特性的本构模型,这样就能够对机场高填方工程问题进行合理的分析。
目前在数值分析中广泛应用的MCC等弹塑性模型,主要针对的是细粒土,只能反映正常固结土或弱超固结土的应力应变特性,对超固结土或土石混合料的应力应变特性还不能够合理反映。所以迫切需要建立一种针对土石混合料的、能合理反映土石混合料应力应变特性的本构模型,以便该模型能方便应用于数值分析中,从而能够对机场高填方工程地基变形及稳定性进行合理分析,为实际工程施工提供指导依据。
含石量是影响土石混合料抗剪强度的重要因素,影响着内摩擦角和粘聚力,所以如果在本构模型中考虑含石量,将含石量和粘聚力及内摩擦角的关系引入进来,则能够更加合理地反映土石混合料的抗剪强度和应力应变特性,使得对实际高填方工程问题的分析更加合理。
发明内容
为了解决上述问题,本发明的目的在于提供一种考虑含石量的土石混合料本构模型建立方法。
为了达到上述目的,本发明提供的考虑含石量的土石混合料本构模型建立方法包括按顺序进行的下列步骤:
1)根据土石混合料的含石量与粘聚力及内摩擦角增量的关系图,获得土石混合料的含石量与内摩擦角关系式以及含石量与土的粘聚力关系式,由此得到土石混合料的内摩擦角及土的粘聚力;
2)根据土的强度准则中扩展松冈-中井准则,并参考统一硬化模型中变换应力公式,结合步骤1)得到的土石混合料的内摩擦角及土的粘聚力,获得土的粘聚力不为0时的变换应力;
3)基于上述获得的土的粘聚力不为0时的变换应力,考虑含石量和内摩擦角及土的粘聚力的关系,参考统一硬化模型的屈服函数公式,在三维化的变换应力空间中,得到考虑含石量的土石混合料本构模型的屈服函数;
4)利用塑性应变增量表达式计算出塑性应变增量;
5)对上述步骤3)获得的屈服函数求微分并代入步骤4)获得的塑性应变增量,得到塑性标量因子;
6)将步骤5)获得的塑性标量因子代入步骤4)中的塑性应变增量表达式,然后再代入应力增量和应变增量的关系式,得到弹塑性刚度矩阵;
7)对上述步骤6)获得的弹塑性刚度矩阵进行对称化处理,得到对称的弹塑性刚度矩阵,最后基于对称的弹塑性刚度矩阵确定出考虑含石量的土石混合料本构模型。
在步骤1)中,所述根据土石混合料的含石量与粘聚力及内摩擦角增量的关系图,获得土石混合料的含石量与内摩擦角关系式以及含石量与土的粘聚力关系式,由此得到土石混合料的内摩擦角及土的粘聚力的具体方法如下:
土石混合料的内摩擦角与含石量的关系式为:
Figure BDA0003141687750000031
其中,
Figure BDA0003141687750000032
式中,
Figure BDA0003141687750000033
为纯土的内摩擦角;
Figure BDA0003141687750000034
为土石混合料的内摩擦角;PR为土石混合料的含石量;
Figure BDA0003141687750000035
为当含石量为PR时土石混合料的内摩擦角较相应纯土的内摩擦角增量;
Figure BDA0003141687750000036
为当含石量PR超过70%时土石混合料的内摩擦角较相应纯土的内摩擦角增量;
土的粘聚力与含石量的关系式为:
c=ccg(PR) (3)
Figure BDA0003141687750000037
其中,cc为纯土的粘聚力;c为土的粘聚力。
在步骤2)中,所述根据土的强度准则中扩展松冈-中井准则,并参考统一硬化模型中变换应力公式,结合步骤1)得到的土石混合料的内摩擦角及土的粘聚力,获得土的粘聚力不为0时的变换应力的具体方法如下:
土的粘聚力不为0时的变换应力为:
Figure BDA0003141687750000041
其中,
Figure BDA0003141687750000042
Figure BDA0003141687750000043
Figure BDA0003141687750000044
Figure BDA0003141687750000045
Figure BDA0003141687750000046
其中,σ为正应力;
Figure BDA0003141687750000047
为土的粘聚力c=0情况下用扩展SMP准则表示的变换应力;
Figure BDA0003141687750000048
为土的粘聚力c≠0情况下用扩展SMP准则表示的平均主应力;
Figure BDA0003141687750000049
为土的粘聚力c≠0情况下用扩展SMP准则表示的三轴压缩时的偏应力;
Figure BDA00031416877500000410
为土的粘聚力c≠0情况下用扩展SMP准则表示的应力不变量;
Figure BDA00031416877500000411
为扩展SMP面上的三个主应力。
在步骤3)中,所述基于上述获得的土的粘聚力不为0时的变换应力,考虑含石量和内摩擦角及土的粘聚力的关系,参考统一硬化模型的屈服函数公式,在三维化的变换应力空间中,得到考虑含石量的土石混合料本构模型的屈服函数的具体方法如下:
考虑含石量的土石混合料本构模型的屈服函数为:
Figure BDA0003141687750000051
式中,
Figure BDA0003141687750000052
为变换应力空间当前应力点上的平均主应力;
Figure BDA0003141687750000053
为该点上的广义剪应力;
Figure BDA0003141687750000054
为相应的初始平均主应力;cp=(λ-κ)/1+e0,λ为等向压缩线的斜率,κ为等向回弹线的斜率,e0为初始孔隙比;M为特征状态和临界状态下的应力比,
Figure BDA0003141687750000055
为硬化参数,可表示为:
Figure BDA0003141687750000056
Figure BDA0003141687750000057
式中,
Figure BDA0003141687750000058
为变换应力空间的潜在强度,
Figure BDA0003141687750000059
为变换应力空间的超固结参数,
Figure BDA00031416877500000510
为前期固结压力,
Figure BDA00031416877500000511
为塑性体积应变,
Figure BDA00031416877500000512
为变换应力空间的应力比。
在步骤4)中,所述利用塑性应变增量表达式计算出塑性应变增量的具体方法如下:
根据塑性位势理论获得塑性应变增量的表达式为:
Figure BDA00031416877500000513
式中,dεp为塑性应变增量;Λ为塑性标量因子。
在步骤5)中,所述对上述步骤3)获得的屈服函数求微分并代入步骤4)获得的塑性应变增量,得到塑性标量因子的具体方法如下:
塑性标量因子为:
Figure BDA0003141687750000061
式中,
Figure BDA0003141687750000062
Figure BDA0003141687750000063
在步骤6)中,所述将步骤5)获得的塑性标量因子代入步骤4)中的塑性应变增量表达式,然后再代入应力增量和应变增量的关系式,得到弹塑性刚度矩阵的具体方法如下:
将公式(15)代入式(14),得到:
Figure BDA0003141687750000064
又知,应力增量和应变增量的关系式为:
dσ=De(dε-dεp) (19)
将公式(18)代入上述应力增量和应变增量的关系式(19)得到:
Figure BDA0003141687750000065
式(20)可以写为:
dσ=Depdε (21)
从而得到弹塑性刚度矩阵Dep为:
Figure BDA0003141687750000071
在步骤7)中,所述对上述步骤6)获得的弹塑性刚度矩阵进行对称化处理,得到对称的弹塑性刚度矩阵,最后基于对称的弹塑性刚度矩阵确定出考虑含石量的土石混合料本构模型的具体方法如下:
对弹塑性刚度矩阵进行对称化处理,得到对称的弹塑性刚度矩阵为:
Figure BDA0003141687750000072
其中,
Figure BDA0003141687750000073
基于对称的弹塑性刚度矩阵获得考虑含石量的土石混合料本构模型为:
dσ=Depdε (25)。
本发明的效果在于:
1)在土石混合料本构模型中考虑含石量,将含石量和土的粘聚力及内摩擦角的关系引入进来,则能够更加合理地反映土石混合料的抗剪强度和应力应变特性,使得对实际高填方工程问题的分析更加合理。
2)所建立的土石混合料本构模型,弥补了传统的岩土本构模型仅能描述细粒土且为正常固结或若超固结土的力学特性,本模型能够考虑土石混合料的含石量和超固结性,合理反映土石混合料的力学特性。
3)所建立的土石混合料本构模型参数简单易测,便于在有限元中应用。
附图说明
图1为本发明中扩展SMP面上的剪应力和正应力。
图2为本发明提供的考虑含石量的土石混合料本构模型建立方法流程图。
图3为本发明中应用考虑含石量的土石混合料本构模型对某机场高填方地基进行数值分析所得到的沉降监测点的竖向位移随时间变化曲线。
具体实施方式
下面结合附图和具体实施例对本发明提供的考虑含石量的土石混合料本构模型建立方法进行详细说明。
如图2所示,本发明提供的考虑含石量的土石混合料本构模型建立方法包括按顺序进行的下列步骤:
1)根据土石混合料的含石量与粘聚力及内摩擦角增量的关系图,获得土石混合料的含石量与内摩擦角关系式以及含石量与土的粘聚力关系式,由此得到土石混合料的内摩擦角及土的粘聚力;
土石混合料的内摩擦角与含石量的关系式为:
Figure BDA0003141687750000081
其中,
Figure BDA0003141687750000082
式中,
Figure BDA0003141687750000083
为纯土的内摩擦角;
Figure BDA0003141687750000084
为土石混合料的内摩擦角;PR为土石混合料的含石量;
Figure BDA0003141687750000085
为当含石量为PR时土石混合料的内摩擦角较相应纯土的内摩擦角增量;
Figure BDA0003141687750000086
为当含石量PR超过70%时土石混合料的内摩擦角较相应纯土的内摩擦角增量。
土的粘聚力与含石量的关系式为:
c=ccg(PR) (3)
Figure BDA0003141687750000087
其中,cc为纯土的粘聚力;c为土的粘聚力。
2)根据土的强度准则中扩展松冈-中井(SMP)准则,如图1所示,并参考现有技术的统一硬化模型(UH模型)中变换应力公式,结合步骤1)得到的土石混合料的内摩擦角及土的粘聚力,获得土的粘聚力不为0时的变换应力;
土的粘聚力不为0时的变换应力为:
Figure BDA0003141687750000091
其中,
Figure BDA0003141687750000092
Figure BDA0003141687750000093
Figure BDA0003141687750000094
Figure BDA0003141687750000095
Figure BDA0003141687750000096
其中,σ为正应力;
Figure BDA0003141687750000097
为土的粘聚力c=0情况下用扩展SMP准则表示的变换应力;
Figure BDA0003141687750000098
为土的粘聚力c≠0情况下用扩展SMP准则表示的平均主应力;
Figure BDA0003141687750000099
为土的粘聚力c≠0情况下用扩展SMP准则表示的三轴压缩时的偏应力;
Figure BDA00031416877500000910
为土的粘聚力c≠0情况下用扩展SMP准则表示的应力不变量;
Figure BDA00031416877500000911
为扩展SMP面上的三个主应力。
3)基于上述获得的土的粘聚力不为0时的变换应力,考虑含石量和内摩擦角及土的粘聚力的关系,参考统一硬化模型的屈服函数公式,在三维化的变换应力空间中,得到考虑含石量的土石混合料本构模型的屈服函数;
考虑含石量的土石混合料本构模型的屈服函数为:
Figure BDA0003141687750000101
式中,
Figure BDA0003141687750000102
为变换应力空间当前应力点上的平均主应力;
Figure BDA0003141687750000103
为该点上的广义剪应力;
Figure BDA0003141687750000104
为相应的初始平均主应力;cp=(λ-κ)/1+e0,λ为等向压缩线的斜率,κ为等向回弹线的斜率,e0为初始孔隙比;M为特征状态和临界状态下的应力比,
Figure BDA0003141687750000105
为硬化参数,可表示为:
Figure BDA0003141687750000106
Figure BDA0003141687750000107
式中,
Figure BDA0003141687750000108
为变换应力空间的潜在强度,
Figure BDA0003141687750000109
为变换应力空间的超固结参数,
Figure BDA00031416877500001010
为前期固结压力,
Figure BDA00031416877500001011
为塑性体积应变,
Figure BDA00031416877500001012
为变换应力空间的应力比。
4)利用塑性应变增量表达式计算出塑性应变增量;
根据塑性位势理论获得塑性应变增量的表达式为:
Figure BDA00031416877500001013
式中,dεp为塑性应变增量;Λ为塑性标量因子。
5)对上述步骤3)获得的屈服函数求微分并代入步骤4)获得的塑性应变增量,得到塑性标量因子;
塑性标量因子为:
Figure BDA0003141687750000111
式中,
Figure BDA0003141687750000112
Figure BDA0003141687750000113
6)将步骤5)获得的塑性标量因子代入步骤4)中的塑性应变增量表达式,然后再代入应力增量和应变增量的关系式,得到弹塑性刚度矩阵;
将公式(15)代入式(14),得到:
Figure BDA0003141687750000114
又知,应力增量和应变增量的关系式为:
dσ=De(dε-dεp) (19)
将公式(18)代入上述应力增量和应变增量的关系式(19)得到:
Figure BDA0003141687750000115
式(20)可以写为:
dσ=Depdε (21)
从而得到弹塑性刚度矩阵Dep为:
Figure BDA0003141687750000121
7)对上述步骤6)获得的弹塑性刚度矩阵进行对称化处理,得到对称的弹塑性刚度矩阵,最后基于对称的弹塑性刚度矩阵确定出考虑含石量的土石混合料本构模型。
对弹塑性刚度矩阵进行对称化处理,得到对称的弹塑性刚度矩阵为:
Figure BDA0003141687750000122
其中,
Figure BDA0003141687750000123
基于对称的弹塑性刚度矩阵获得考虑含石量的土石混合料本构模型为:
dσ=Depdε (25)
图3为利用本发明提供的考虑含石量的土石混合料本构模型对机场高填方地基进行数值分析所得到的沉降监测点的竖向位移随时间变化曲线与实测数据的对比,由图3可知,两条曲线比较吻合,说明采用本发明提供的考虑含石量的土石混合料本构模型能够比较合理地分析机场高填方地基变形问题。

Claims (8)

1.一种考虑含石量的土石混合料本构模型建立方法,其特征在于:所述考虑含石量的土石混合料本构模型建立方法包括按顺序进行的下列步骤:
1)根据土石混合料的含石量与粘聚力及内摩擦角增量的关系图,获得土石混合料的含石量与内摩擦角关系式以及含石量与土的粘聚力关系式,由此得到土石混合料的内摩擦角及土的粘聚力;
2)根据土的强度准则中扩展松冈-中井准则,并参考统一硬化模型中变换应力公式,结合步骤1)得到的土石混合料的内摩擦角及土的粘聚力,获得土的粘聚力不为0时的变换应力;
3)基于上述获得的土的粘聚力不为0时的变换应力,考虑含石量和内摩擦角及土的粘聚力的关系,参考统一硬化模型的屈服函数公式,在三维化的变换应力空间中,得到考虑含石量的土石混合料本构模型的屈服函数;
4)利用塑性应变增量表达式计算出塑性应变增量;
5)对上述步骤3)获得的屈服函数求微分并代入步骤4)获得的塑性应变增量,得到塑性标量因子;
6)将步骤5)获得的塑性标量因子代入步骤4)中的塑性应变增量表达式,然后再代入应力增量和应变增量的关系式,得到弹塑性刚度矩阵;
7)对上述步骤6)获得的弹塑性刚度矩阵进行对称化处理,得到对称的弹塑性刚度矩阵,最后基于对称的弹塑性刚度矩阵确定出考虑含石量的土石混合料本构模型。
2.根据权利要求1所述考虑含石量的土石混合料本构模型建立方法,其特征在于:在步骤1)中,所述根据土石混合料的含石量与粘聚力及内摩擦角增量的关系图,获得土石混合料的含石量与内摩擦角关系式以及含石量与土的粘聚力关系式,由此得到土石混合料的内摩擦角及土的粘聚力的具体方法如下:
土石混合料的内摩擦角与含石量的关系式为:
Figure FDA0003141687740000021
其中,
Figure FDA0003141687740000022
式中,
Figure FDA0003141687740000023
为纯土的内摩擦角;
Figure FDA0003141687740000024
为土石混合料的内摩擦角;PR为土石混合料的含石量;
Figure FDA0003141687740000025
为当含石量为PR时土石混合料的内摩擦角较相应纯土的内摩擦角增量;
Figure FDA0003141687740000026
为当含石量PR超过70%时土石混合料的内摩擦角较相应纯土的内摩擦角增量;
土的粘聚力与含石量的关系式为:
c=ccg(PR) (3)
Figure FDA0003141687740000027
其中,cc为纯土的粘聚力;c为土的粘聚力。
3.根据权利要求1所述考虑含石量的土石混合料本构模型建立方法,其特征在于:在步骤2)中,所述根据土的强度准则中扩展松冈-中井准则,并参考统一硬化模型中变换应力公式,结合步骤1)得到的土石混合料的内摩擦角及土的粘聚力,获得土的粘聚力不为0时的变换应力的具体方法如下:
土的粘聚力不为0时的变换应力为:
Figure FDA0003141687740000031
其中,
Figure FDA0003141687740000032
Figure FDA0003141687740000033
Figure FDA0003141687740000034
Figure FDA0003141687740000035
Figure FDA0003141687740000036
其中,σ为正应力;
Figure FDA0003141687740000037
为土的粘聚力c=0情况下用扩展SMP准则表示的变换应力;
Figure FDA0003141687740000038
为土的粘聚力c≠0情况下用扩展SMP准则表示的平均主应力;
Figure FDA0003141687740000039
为土的粘聚力c≠0情况下用扩展SMP准则表示的三轴压缩时的偏应力;
Figure FDA00031416877400000310
为土的粘聚力c≠0情况下用扩展SMP准则表示的应力不变量;
Figure FDA00031416877400000311
为扩展SMP面上的三个主应力。
4.根据权利要求1所述考虑含石量的土石混合料本构模型建立方法,其特征在于:在步骤3)中,所述基于上述获得的土的粘聚力不为0时的变换应力,考虑含石量和内摩擦角及土的粘聚力的关系,参考统一硬化模型的屈服函数公式,在三维化的变换应力空间中,得到考虑含石量的土石混合料本构模型的屈服函数的具体方法如下:
考虑含石量的土石混合料本构模型的屈服函数为:
Figure FDA0003141687740000041
式中,
Figure FDA0003141687740000042
为变换应力空间当前应力点上的平均主应力;
Figure FDA0003141687740000043
为该点上的广义剪应力;
Figure FDA0003141687740000044
为相应的初始平均主应力;cp=(λ-κ)/1+e0,λ为等向压缩线的斜率,κ为等向回弹线的斜率,e0为初始孔隙比;M为特征状态和临界状态下的应力比,
Figure FDA0003141687740000045
为硬化参数,可表示为:
Figure FDA0003141687740000046
Figure FDA0003141687740000047
式中,
Figure FDA0003141687740000048
为变换应力空间的潜在强度,
Figure FDA0003141687740000049
为变换应力空间的超固结参数,
Figure FDA00031416877400000410
为前期固结压力,
Figure FDA00031416877400000411
为塑性体积应变,
Figure FDA00031416877400000412
为变换应力空间的应力比。
5.根据权利要求1所述考虑含石量的土石混合料本构模型建立方法,其特征在于:在步骤4)中,所述利用塑性应变增量表达式计算出塑性应变增量的具体方法如下:
根据塑性位势理论获得塑性应变增量的表达式为:
Figure FDA00031416877400000413
式中,dεp为塑性应变增量;Λ为塑性标量因子。
6.根据权利要求1所述考虑含石量的土石混合料本构模型建立方法,其特征在于:在步骤5)中,所述对上述步骤3)获得的屈服函数求微分并代入步骤4)获得的塑性应变增量,得到塑性标量因子的具体方法如下:
塑性标量因子为:
Figure FDA0003141687740000051
式中,
Figure FDA0003141687740000052
Figure FDA0003141687740000053
7.根据权利要求1所述考虑含石量的土石混合料本构模型建立方法,其特征在于:在步骤6)中,所述将步骤5)获得的塑性标量因子代入步骤4)中的塑性应变增量表达式,然后再代入应力增量和应变增量的关系式,得到弹塑性刚度矩阵的具体方法如下:
将公式(15)代入式(14),得到:
Figure FDA0003141687740000054
又知,应力增量和应变增量的关系式为:
dσ=De(dε-dεp) (19)
将公式(18)代入上述应力增量和应变增量的关系式(19)得到:
Figure FDA0003141687740000061
式(20)可以写为:
dσ=Depdε (21)
从而得到弹塑性刚度矩阵Dep为:
Figure FDA0003141687740000062
8.根据权利要求1所述考虑含石量的土石混合料本构模型建立方法,其特征在于:在步骤7)中,所述对上述步骤6)获得的弹塑性刚度矩阵进行对称化处理,得到对称的弹塑性刚度矩阵,最后基于对称的弹塑性刚度矩阵确定出考虑含石量的土石混合料本构模型的具体方法如下:
对弹塑性刚度矩阵进行对称化处理,得到对称的弹塑性刚度矩阵为:
Figure FDA0003141687740000063
其中,
Figure FDA0003141687740000064
基于对称的弹塑性刚度矩阵获得考虑含石量的土石混合料本构模型为:
dσ=Depdε (25)。
CN202110735907.2A 2021-06-30 2021-06-30 一种考虑含石量的土石混合料本构模型建立方法 Pending CN113297665A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110735907.2A CN113297665A (zh) 2021-06-30 2021-06-30 一种考虑含石量的土石混合料本构模型建立方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110735907.2A CN113297665A (zh) 2021-06-30 2021-06-30 一种考虑含石量的土石混合料本构模型建立方法

Publications (1)

Publication Number Publication Date
CN113297665A true CN113297665A (zh) 2021-08-24

Family

ID=77329999

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110735907.2A Pending CN113297665A (zh) 2021-06-30 2021-06-30 一种考虑含石量的土石混合料本构模型建立方法

Country Status (1)

Country Link
CN (1) CN113297665A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114112685A (zh) * 2021-12-29 2022-03-01 国网新源控股有限公司 一种现场压实土石混合料前期固结应力确定方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110926944A (zh) * 2019-12-05 2020-03-27 中国科学院武汉岩土力学研究所 一种适用于循环荷载下岩石损伤本构模型的构建方法
US20200394347A1 (en) * 2019-06-12 2020-12-17 Sichuan University Method for assessing fatigue damage and fatigue life based on abaqus
CN113032955A (zh) * 2021-02-05 2021-06-25 中国科学院武汉岩土力学研究所 一种适用于地震荷载下岩石动态本构模型的构建方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200394347A1 (en) * 2019-06-12 2020-12-17 Sichuan University Method for assessing fatigue damage and fatigue life based on abaqus
CN110926944A (zh) * 2019-12-05 2020-03-27 中国科学院武汉岩土力学研究所 一种适用于循环荷载下岩石损伤本构模型的构建方法
CN113032955A (zh) * 2021-02-05 2021-06-25 中国科学院武汉岩土力学研究所 一种适用于地震荷载下岩石动态本构模型的构建方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
冯兴 等: "山区机场高填方地基变形分析", 《北京航空航天大学学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114112685A (zh) * 2021-12-29 2022-03-01 国网新源控股有限公司 一种现场压实土石混合料前期固结应力确定方法

Similar Documents

Publication Publication Date Title
Wang et al. Shear behaviour and acoustic emission characteristics of bolted rock joints with different roughnesses
Su et al. Influences of overburden pressure and soil dilation on soil nail pull-out resistance
Shan et al. Study on the triaxial unloading creep mechanical properties and damage constitutive model of red sandstone containing a single ice-filled flaw
Xu et al. Three-dimensional stability analysis of slope in unsaturated soils considering strength nonlinearity under water drawdown
CN110595887B (zh) 一种浅埋硬岩大型地下洞室成洞条件量化评价方法
CN113820750A (zh) 基于弹塑性力学定量预测泥岩构造裂缝的方法
Tu et al. Elastoplastic coupling analysis of high-strength concrete based on tests and the Mohr-Coulomb criterion
CN110057675A (zh) 一种土体受扰动后的e-logp’曲线的确定方法
Wang et al. Numerical investigation on anisotropy and shape effect of mechanical properties of columnar jointed basalts containing transverse joints
Zhu et al. Experimental study on the strength failure characteristics of columnar jointed rock masses under three-dimensional stress
Zhou et al. Study of the coupling effect of elliptical cavities and cracks on tunnel stability under dynamic loads
CN113297665A (zh) 一种考虑含石量的土石混合料本构模型建立方法
Wu et al. Modeling of compaction grouting considering the soil unloading effect
Oliaei et al. Evaluation of soil-nail pullout resistance using mesh-free method
Chen et al. Study on the influence of crack depth on the safety of tunnel lining structure
Wang et al. Model tests on dynamic responses of surrounding rock and support structure on underground tunnel under combined dynamic and static loading
Xia et al. Kernel broken smooth particle hydrodynamics method for crack propagation simulation applied in layered rock cells and tunnels
CN116244780A (zh) 高填方湿陷性黄土路基强夯填筑施工计算方法
Wen et al. The use of Hoek Brown failure criterion on determination of the geo-mechanical parameters of a grouting consolidation body
Ma et al. The calculation for ultimate surrounding earth pressure on deep-buried tunnels in aeolian sand stratum to prevent surface collapse
Liu et al. Experimental study on mechanical and fracture characteristics of inclined weak-filled rough joint rock-like specimens
Lu et al. Mechanical behavior of grouted rock joints
CHENG et al. Establishment of constitutive relation of shear deformation for irregular joints in sandstone
Wang et al. Experimental investigation of crack initiation and propagation in the unreinforced masonry specimen subjected to vertical settlement
Lyu et al. Characteristics of early-damaged cement stone and the creep behavior after damage

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20210824