CN113408103B - 基于扰动状态概念的结构面剪切本构模型的构建方法 - Google Patents

基于扰动状态概念的结构面剪切本构模型的构建方法 Download PDF

Info

Publication number
CN113408103B
CN113408103B CN202110557607.XA CN202110557607A CN113408103B CN 113408103 B CN113408103 B CN 113408103B CN 202110557607 A CN202110557607 A CN 202110557607A CN 113408103 B CN113408103 B CN 113408103B
Authority
CN
China
Prior art keywords
shear
response
state
rock
disturbance
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.)
Active
Application number
CN202110557607.XA
Other languages
English (en)
Other versions
CN113408103A (zh
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.)
Dragon Totem Technology Hefei Co ltd
Original Assignee
Central South University
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 Central South University filed Critical Central South University
Priority to CN202110557607.XA priority Critical patent/CN113408103B/zh
Publication of CN113408103A publication Critical patent/CN113408103A/zh
Application granted granted Critical
Publication of CN113408103B publication Critical patent/CN113408103B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • 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
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Health & Medical Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明公开了一种基于扰动状态概念的结构面剪切本构模型的构建方法，其包括如下步骤：基于扰动状态概念理论，定义扰动函数DF，建立RI响应及FA响应相对于观测响应的贡献；结合典型结构面剪切变形本构关系，分别获取初始响应及临界响应；通过剪切试验观测量表征扰动函数DF，建立基于扰动状态概念的结构面剪切变形本构模型；确定基于扰动状态概念的结构面剪切变形本构模型中的DSC参数；将结构面剪切试验数据代入基于扰动状态概念的结构面剪切本构模型中，验证得到的结构面剪切本构模型能够准确模拟结果面剪切变形行为。本发明通过采用扰动状态概念理论，建立节理剪切本构模型，模拟各剪切变形阶段的非线性特性，提升模型的实用性。

Description

基于扰动状态概念的结构面剪切本构模型的构建方法
技术领域
本发明涉及岩石结构面本构模型技术领域，尤其是涉及一种基于扰动状态概念的结构面剪切本构模型的构建方法。
背景技术
结构面纵横穿插在节理岩体中，由于结构面的承载能力、变形能力和稳定性弱于完整岩体，因而控制了岩体结构的变形和稳定。节理岩体在受到外部荷载作用或地下工程开挖、地震等环境剧烈变化时，可能会沿结构面滑移。因此，结构面的剪切力学行为对评价岩质基础、矿山、隧道等岩土结构的变形能力和安全性具有重要影响。对于结构面剪切变形破坏机理已有大量的理论建模和预测研究。 Goodman首先在有限元方法的基础上建立了岩石结构面模型，开创了数值模拟在结构面建模中的应用；考虑到结构面剪切变形的峰后软化特性，Simon建立了 CSDS(completestress-displacementsurface)本构模型，模拟了全剪切应力-剪切位移曲线；Lietal.考虑了初始张开度对结构面剪切力学性能的影响，建立了新的可以纳入数值软件模拟的剪切本构模型；PouyaandBemaniYazdi针对岩石结构面弹塑性变形损伤过程提出了一种新的结构面完整本构-塑性模型；另一方面，岩土工程从业专家和工程师利用弹塑性理论对岩石结构面剪切特性进行了广泛的研究，也极大地提高了地质工程工作者对结构面剪切变形行为的认识。
尽管上述建立的大部分本构模型都能简单地描述结构面剪切变形的阶段特性，但仍存在一些缺陷。例如，有的剪切本构模型忽略了结构面剪切过程中的非线性变形特性，将剪切变形视为线性变形，这与实际情况有很大的不同；有些模型的参数没有物理和力学意义；有些模型过于复杂，在数值计算中难以迭代，进而降低本构模型的实用性，不利于本构模型的推广及研究。
扰动状态概念(DSC)理论最早由FrantziskonisandDesai[1]创立，为地质材料和界面提供了一种独特的本构模拟方法。该方法从微观角度阐述了地质材料变形过程的非线性行为，为将细观单元的力学响应与宏观本构关系[2]相结合架起了桥梁，从扰动状态概念(DSC)的角度研究地质材料的破坏机理和本构关系已成为近年来的研究热点[3，4]。
[1]Frantziskonis G,Desai CS.Elastoplastic model with damage forstrain softening geomaterials.Acta Mechanica.1987；68:151-70.
[2]Sane SM,Desai CS,Jenson JW,Contractor DN,Carlson AE,Clark PU.Disturbed State constitutive modeling of two Pleistocene tills.QuaternaryScience Reviews.2008；27:267-83.
[3]Xiao Y,Desai CS.Constitutive Modeling for Overconsolidated ClaysBased on Disturbed State Concept.I:Theory.International Journal OfGeomechanics. 2019；19:18.
[4]Fan R-D,Liu M,Du Y-J,Horpibulsuk S.Estimating the compressionbehaviour of metal-rich clays via a Disturbed State Concept(DSC)model.AppliedClay Science.2016；132-133:50-8.
发明内容
基于此，本发明的目的在于提供一种基于扰动状态概念的结构面剪切本构模型的构建方法，采用扰动状态概念理论，建立节理剪切本构模型，模拟各剪切变形阶段的非线性特性，提升模型的实用性。
为解决上述技术问题，本发明采用以下技术方案：
本发明提供了一种基于扰动状态概念的结构面剪切本构模型的构建方法，其包括如下步骤：
基于扰动状态概念理论，定义扰动函数DF，建立RI响应及FA响应相对于观测响应的贡献；其中，FA响应相对于观测响应的贡献通过DF来进行表示， RI响应相对于观测响应的贡献通过1-DF来进行表示，0≤DF≤1；RI响应指代为处于RI状态的材料剪切力学响应，即初始响应；FA响应指代为处于FA状态的材料剪切力学响应，即临界响应；观测响应指代为观测到的岩石结构面材料整体的剪切力学响应；
结合典型结构面剪切变形本构关系，分别获取初始响应及临界响应；
通过剪切试验观测量表征扰动函数DF，建立基于扰动状态概念的结构面剪切变形本构模型；
确定基于扰动状态概念的结构面剪切变形本构模型中的DSC参数；
将结构面剪切试验数据代入基于扰动状态概念的结构面剪切本构模型中，验证得到的结构面剪切本构模型能够准确模拟结果面剪切变形行为。
在其中一个实施例中，所述步骤基于扰动状态概念理论，定义扰动函数DF，建立RI响应及FA响应相对于观测响应的贡献之前，还包括：
定义岩石结构面材料整体由无数个岩石单元组成，其中，岩石单元总数为l， l→∞，单个岩石单元的面积是A0，剪切面积A可以表示为A＝lA0，任意时刻处于RI状态的岩石单元数量为m，处于FA状态的岩石单元数量为n，l＝m+n；
获取力学平衡方程τA＝τRIARI+τFAAFA，其中，τ和A代表剪切应力和对应的剪切面积，τRI指的是处于RI状态下岩石单元所承受的剪切应力，ARI是处于RI状态下岩石单元的总面积，τFA表示处于FA状态下岩石单元所承受的剪切应力，AFA是处于FA状态下岩石单元的总面积；
获取剪切面积A、处于RI状态下岩石单元的总面积ARI、处于FA状态下岩石单元的总面积AFA：
Figure BDA0003077759360000041
将
Figure BDA0003077759360000042
代入到力学平衡方程τA＝τRIARI+τFAAFA中，并同时除以lA0后，得到剪切应力τ：
Figure BDA0003077759360000043
获取扰动函数
Figure BDA0003077759360000044
将表达式
Figure BDA0003077759360000045
及表达式l＝m+n代入到表达式
Figure BDA0003077759360000046
中，剪切应力τ的表达式变为：
τ＝τRI(1-DF)+τFADF。
在其中一个实施例中，所述步骤结合典型结构面剪切变形本构关系，分别获取初始响应及临界响应的方法，具体操作过程包括：
结合典型结构面剪切变形本构关系曲线剪切应力τ-剪切位移δ，分别获取初始响应τ＝ksδ及临界响应τ＝τr，其中，ks指代为剪切刚度，τr指代为残余剪切强度。
在其中一个实施例中，所述步骤通过宏观观测量表征扰动函数DF，建立基于扰动状态概念的结构面剪切变形本构模型的方法，包括如下步骤，
随着剪切变形过程的推进如剪切位移δ的增加，由RI状态变为FA状态的岩石单元数量也在增加，处于FA状态的岩石单元数量n与剪切位移δ的关系如下：
Figure BDA0003077759360000047
其中，r为比例参数，表示岩石单元的增长速率；
将表达式
Figure BDA0003077759360000048
两边求导得到
Figure BDA0003077759360000049
后，结合等式
Figure BDA00030777593600000410
及
Figure BDA0003077759360000051
获取dDF/dδ，
Figure 1
对
Figure 2
进行积分后可得到扰动函数DF的表达式为：
Figure BDA0003077759360000054
其中，C1和C2是
Figure 3
积分过程中产生的积分参数；
定义参数
Figure BDA0003077759360000056
将表达式
Figure BDA0003077759360000057
代入到τ＝ksδ(1-DF)+τrDF中后，获得基于扰动状态概念的结构面剪切变形本构模型
Figure BDA0003077759360000058
其中，ks、τr为剪切力学参数，可通过剪切试验直接获取，r、η为DSC参数。
在其中一个实施例中，所述步骤确定基于扰动状态概念的结构面剪切变形本构模型中的DSC参数的方法，具体操作包括，
根据峰值强度点(δp，τp)对应的剪切应力-剪切位移关系满足等式
Figure BDA0003077759360000059
并且在峰值强度点(δp，τp)处剪切应力τ相对于剪切位移δ的导数为0，获得如下公式：
Figure BDA00030777593600000510
Figure BDA00030777593600000511
联立上述两公式后，获得基于扰动状态概念的结构面剪切变形本构模型中的DSC参数η和r的计算公式：
Figure BDA0003077759360000061
在其中一个实施例中，所述基于扰动状态概念理论，定义扰动函数DF，建立RI响应及FA响应相对于观测响应的贡献的方法，包括如下步骤：
在DSC理论框架中，将变形岩石单元设定为RI状态岩石单元及FA状态岩石单元的组合；
定义扰动函数DF，获取RI状态岩石材料单元的响应及FA状态岩石材料单元的响应分别相对于观测到的岩石结构面材料整体的剪切力学响应的贡献。
综上所述，本发明提供的基于扰动状态概念的结构面剪切本构模型的构建方法采用扰动状态概念理论，建立节理剪切本构模型，模拟各剪切变形阶段的非线性特性，提升模型的实用性。
附图说明
图1为本发明实施例提供的基于扰动状态概念的结构面剪切本构模型的构建方法的流程示意图；
图2为本发明实施例提供的扰动状态概念下岩石结构面剪切过程示意图；
图3为本发明实施例提供的结构面剪切力学响应曲线示意图；
图4为本发明实施例提供的本构模型曲线与基于石英岩结构面的剪切试验数据对比图；
图5为本发明实施例提供的本构模型曲线与其中四种结构面的剪切试验数据对比图；
图6为本发明实施例提供的本构模型曲线与另外其中四种结构面的剪切试验数据对比图。
具体实施方式
为能进一步了解本发明的特征、技术手段以及所达到的具体目的、功能，下面结合附图与具体实施方式对本发明作进一步详细描述。
图1是本发明实施例提供的第一种基于扰动状态概念的结构面剪切本构模型的构建方法的流程示意图，如图1所示，该基于扰动状态概念的结构面剪切本构模型的构建方法，具体包括如下步骤：
步骤S1、基于扰动状态概念理论，定义扰动函数DF，建立RI响应及FA响应相对于观测响应的贡献；其中，FA响应相对于观测响应的贡献通过DF来进行表示，RI响应相对于观测响应的贡献通过1-DF来进行表示，0≤DF≤1；RI 响应指代为处于RI状态(相对完整状态)的材料剪切力学响应，即初始响应； FA响应指代为处于FA状态(完全调整状态)的材料剪切力学响应，即临界响应；观测响应指代为观测到的岩石结构面材料整体的剪切力学响应。
所述步骤S1、基于扰动状态概念理论，定义扰动函数DF，建立RI响应及 FA响应相对于观测响应的贡献的方法，包括如下步骤：
步骤S11、在DSC理论框架中，将变形岩石单元设定为RI状态岩石单元及 FA状态岩石单元的组合；
步骤S12、定义扰动函数DF，获取RI状态岩石材料单元的响应及FA状态岩石材料单元的响应分别相对于观测到的岩石结构面材料整体的剪切力学响应的贡献；其中，结构面切向荷载施加后，观测到的岩石结构面材料整体的剪切力学响应包括RI状态岩石材料单元的响应和/或FA状态岩石材料单元的响应，FA 状态岩石材料单元的响应相对于观测到的岩石结构面材料整体的剪切力学响应的贡献为DF，RI状态岩石材料单元的响应相对于观测到的岩石结构面材料整体的剪切力学响应的贡献为1-DF。
请参考图2，在其中一个实施例中，所述步骤S1、基于扰动状态概念理论，定义扰动函数DF，建立RI响应及FA响应相对于观测响应的贡献之前，还包括
步骤S1-1、定义岩石结构面材料整体由无数个岩石单元组成，其中，岩石单元总数为l，l→∞，单个岩石单元的面积是A0，剪切面积A可以表示为A＝lA0，任意时刻处于RI状态的岩石单元数量为m，处于FA状态的岩石单元数量为n， l＝m+n。
步骤S1-2、获取力学平衡方程τA＝τRIARI+τFAAFA，其中，τ和A代表剪切应力和对应的剪切面积，τRI指的是处于RI状态下岩石单元所承受的剪切应力，ARI是处于RI状态下岩石单元的总面积，τFA表示处于FA状态下岩石单元所承受的剪切应力，AFA是处于FA状态下岩石单元的总面积；在岩石结构面剪切破坏全过程的各个阶段，外力(包括法向和剪切方向)的施加是一个准静态过程。因此，沿岩石结构面的剪切方向应处于准平衡力状态，考虑剪切应力加载方向的力学平衡，从而可以得到所述力学平衡方程τA＝τRIARI+τFAAFA。
在剪切过程发生前，各岩石单元均处于RI状态；在整个岩石结构面材料整体剪切过程中，在剪切作用下，由RI状态转向FA状态的岩石单元逐渐增多，这个过程是快速且不可逆转的；在剪切过程进入残余阶段时，各岩石单元全部进入FA状态。
步骤S1-3、获取剪切面积A、处于RI状态下岩石单元的总面积ARI、处于 FA状态下岩石单元的总面积AFA：
Figure BDA0003077759360000091
步骤S1-4、将步骤S1-3中的
Figure BDA0003077759360000092
代入到步骤S1-2的力学平衡方程τA＝τRIARI+τFAAFA中，并同时除以lA0后，得到剪切应力τ的表达式：
Figure BDA0003077759360000093
步骤S1-5、获取扰动函数
Figure BDA0003077759360000094
FA状态下岩石单元数量及其增加速率从另一个方面表征了岩石材料的破坏程度，此为已知技术，在此不必赘述，因此，通过用FA状态下的岩石单元数量n与岩石单元总数l之比来表示扰动函数DF。
步骤S1-6、将步骤S1-5的表达式
Figure BDA0003077759360000095
及步骤S1-1的表达式l＝m+n代入到步骤S1-4的表达式
Figure BDA0003077759360000096
中，步骤S1-4中剪切应力τ的表达式变为：
τ＝τRI(1-DF)+τFADF。
步骤S2、结合典型结构面剪切变形本构关系，分别获取初始响应及临界响应。
请参考图3，所述步骤S2、结合典型结构面剪切变形本构关系，分别获取初始响应及临界响应的方法，具体操作过程包括：
结合典型结构面剪切变形本构关系曲线剪切应力τ-剪切位移δ，分别获取初始响应τ＝ksδ及临界响应τ＝τr，其中，ks指代为剪切刚度，τr指代为残余剪切强度。
典型结构面剪切全过程本构曲线τ-δ及岩石结构面材料的力学响应可以分为 4个阶段，具体为准弹性阶段(oa区域)、峰前非线性阶段(ab区域)、峰后软化阶段(bc区域)、残余强度阶段(cd区域)；RI响应表示准弹性阶段，此时剪切应力τ和剪切位移δ近似为线弹性关系；FA响应代表了残余强度阶段，剪切应力τ已经达到一个稳定值称为残余剪切强度τr；在结构面剪切试验中观察到的剪切力学响应是RI响应与FA响应的反应，在准弹性阶段，扰动很小，近似为无扰动，岩石单元可以近似认为均处于RI状态，即n＝0，步骤1-5中
Figure BDA0003077759360000101
及步骤1-6中τ＝τRI(1-DF)+τFADF的扰动函数DF＝0，此时RI状态岩石单元的变形特征近似处于准弹性阶段，因此，由等式τ＝τRI(1-DF)+τFADF可以得到表达式τ＝τRI＝ksδ，其中，ks指代为剪切刚度，通过准弹性阶段典型结构面剪切全过程本构模型曲线τ-δ的斜率进行表示；当扰动达到完全扰动阶段时，等式τ＝τRI(1-DF)+τFADF 中DF＝1，图3中所有剪切力学响应均为FA响应，即剪切应力进入残余强度阶段，因此由等式τ＝τRI(1-DF)+τFADF可以得到表达式τ＝τFA＝τr。
具体地，对典型结构面剪切本构关系分为四个阶段作进一步的描述：
(1)准弹性阶段(图3所示oa区域)；剪切刚发生时(起始点o)，结构面表面微凸体在较小剪切位移范围内变形不明显。在该阶段，剪切本构曲线τ-δ近似呈线性形状，其中τ-δ曲线的斜率为剪切刚度ks。在这一过程中，由RI状态向FA状态转变的岩石单元数量很少，因此认为所有岩石单元均处于“相对完整”状态 (RI)；
(2)峰前非线性阶段(图3所示ab区域)，此时结构面剪切变形的非线性特征开始明显。与前一个阶段(oa区域)相比，随着剪切位移δ的增加，剪切应力呈现非线性增加，而且增加的速度明显慢于前一个阶段(oa区域)。此时，转换为“完全调整”状态(FA)的岩石单位数量大幅增加；
(3)峰后软化阶段(图3所示bc区域)。在这一阶段，随着剪切位移δ的增加，剪切应力τ急剧下降，结构面表面的凸体被大量破坏。“完全调整”状态(FA)的岩石单元数量大大增加；
(4)残余强度阶段(图3所示cd区域)。在这一阶段，随着剪切位移δ的增加，剪切应力几乎没有变化。岩石结构面完全失去了粘结力。所有的粗糙表面都经过了磨损，仅保留了稳定的残余摩擦，即残余剪切强度τr。岩石单元几乎都处于“完全调整”状态(FA)。
步骤S3、通过剪切试验观测量表征扰动函数DF，建立基于扰动状态概念的结构面剪切变形本构模型。
具体地，所述步骤S3、通过剪切试验观测量表征扰动函数DF，建立基于扰动状态概念的结构面剪切变形本构模型的方法，包括如下步骤：
步骤S3-1、随着剪切变形过程的推进如剪切位移δ的增加，由RI状态变为 FA状态的岩石单元数量也在增加，某一剪切位移δ下处于FA状态的岩石单元数量n与该剪切位移δ的关系如下：
Figure BDA0003077759360000111
其中，r为比例参数，表示岩石单元的增长速率。
步骤S3-2、将步骤S3-1中的表达式
Figure BDA0003077759360000112
两边求导得到
Figure BDA0003077759360000113
后，结合等式
Figure BDA0003077759360000114
及
Figure BDA0003077759360000115
获取dDF/dδ，
Figure 4
其中，Liuetal.和Liuetal.发现岩石材料破坏规律与生物科学中的种群增长数学模型相似，此为已知技术，本发明将其应用于岩石力学，以获得一个合适的扰动函数DF来描述扰动与岩石结构面剪切特性之间的关系dDF/dδ。
步骤S3-3、对
Figure 5
进行积分后可得到扰动函数 DF的表达式为：
Figure BDA0003077759360000118
其中，C1和C2是
Figure 6
积分过程中产生的积分参数。
步骤S3-4、定义参数
Figure BDA0003077759360000121
将步骤S3-3中的表达式
Figure BDA0003077759360000122
代入到τ＝ksδ(1-DF)+τrDF中后，获得基于扰动状态概念的结构面剪切变形本构模型
Figure BDA0003077759360000123
其中，ks、τr为剪切力学参数，可通过剪切试验直接获取；r、η为DSC参数。
如表1所示，将本发明中基于扰动状态概念的结构面剪切变形本构模型与其余基于弹塑性理论建立的模型、经验模型等进行对比，本发明基于扰动状态概念的结构面剪切变形本构模型形式简单，模型参数少，易于求解且均具有物理意义。
表1常用剪切本构模型及其特点
Figure BDA0003077759360000124
Figure BDA0003077759360000131
步骤S4、确定基于扰动状态概念的结构面剪切变形本构模型中的DSC参数 r、η；其中，根据结构面剪切变形峰值强度点(δp，τp)获取基于扰动状态概念的结构面剪切变形本构模型的DSC参数r和η。
请参阅图3，具体地，所述步骤S4、确定基于扰动状态概念的结构面剪切变形本构模型中的DSC参数r、η的方法，具体操作包括：
根据峰值强度点(δp，τp)对应的剪切应力-剪切位移关系满足等式
Figure BDA0003077759360000132
并且在峰值强度点(δp，τp)处剪切应力τ相对于剪切位移δ的导数为0，获得如下公式：
Figure BDA0003077759360000133
Figure BDA0003077759360000134
联立上述两公式后，获得基于扰动状态概念的结构面剪切变形本构模型中的DSC参数η和r的计算公式：
Figure BDA0003077759360000135
步骤S5、将结构面剪切试验数据代入基于扰动状态概念的结构面剪切本构模型中，验证得到的结构面剪切本构模型能够准确模拟结果面剪切变形行为，并能充分反映出剪切过程各个阶段的剪切变形特征。
具体地，所述步骤S5、将结构面剪切试验数据代入基于扰动状态概念的结构面剪切本构模型中，验证得到的结构面剪切本构模型能够准确模拟结果面剪切变形行为，并能充分反映出剪切过程各个阶段的剪切变形特征，其具体步骤包括：
基于石英岩结构面的剪切试验数据验证
为了比较真实天然岩石结构面与其复制品的剪切特性差异，SinghandBasu 对石英岩结构面进行了多次正应力为0.22～0.71MPa的剪切试验。
本次验证以正应力σn＝0.23MPa和σn＝0.44MPa的试验数据为例，对基于扰动状态概念的结构面剪切变形本构模型进行了详细的验证。σn＝0.23MPa和σn＝0.44MPa的试验数据如表2所示。
表2 SinghandBasu剪切试验结果
σ<sub>n</sub>(MPa) k<sub>s</sub>(MPa/mm) τ<sub>p</sub>(MPa) δ<sub>p</sub>(mm) τ<sub>r</sub>(MPa)
0.23 0.237 0.34 1.50 0.208
0.44 0.279 0.55 2.00 0.418
将表2中所有试验参数如峰值剪切应力τp、峰值剪切位移δp等代入等式
Figure BDA0003077759360000141
中，得到的参数r和η如下:
Figure BDA0003077759360000142
Figure BDA0003077759360000143
将表2中的剪切试验结果和计算所得模型参数r和η分别代入等式
Figure BDA0003077759360000144
得到最终的基于扰动状态概念的结构面剪切变形本构模型表达式如下:
Figure BDA0003077759360000151
请参考图4，图4展示了将基于扰动状态概念的结构面剪切本构模型计算得到的剪切变形曲线与实际剪切变形结果进行对比的示意图。实际剪切本构关系与基于扰动状态概念的结构面剪切本构模型预测曲线吻合较好。从结果来看，两者具有相当的一致性，其中，相关系数R2分别为0.901和0.959，基于扰动状态概念的结构面剪切本构模型能够充分模拟结构面剪切变形中的线性与非线性的行为。
四种天然结构面和人工结构面的直剪试验数据验证
为了直观方便地验证所提出的基于扰动状态概念的结构面剪切本构模型的适应性和合理性，以8个不同类型天然结构面和人工结构面的剪切试验为例进行了验证。这些被选择的天然结构面和人工结构面的试验数据如表3所示。
表3四种天然结构面和人工结构面试验数据及其基于扰动状态概念的结构面剪切本构模型
Figure BDA0003077759360000152
Figure BDA0003077759360000161
如等式
Figure BDA0003077759360000162
所示，分别计算表3中对应试验数据的基于扰动状态概念的结构面剪切本构模型参数、本构模型表达式和相关系数R2等。模型曲线与试验结果的本构关系(τ-δ)曲线对比如图5～6所示。从图5～6的对比结果可以看出，虽然模型曲线与试验数据之间的确存在某些细微差异，但模型曲线能充分反映试验数据呈现出的各个阶段的剪切变形特征。图5～图6的相关系数R2均大于0.88，也证明了模型具有很好的描述性。
综上所述，本发明基于扰动状态概念的结构面剪切变形本构模型，相对于现有剪切本构模型，模型公式简洁、参数易求解且具有明确物理意义，模型曲线表达式与结构面试验数据吻合良好，证明所建立的模型合理，提升模型的实用性。
以上所述实施例仅表达了本发明的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对本发明范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保护范围。因此，本发明的保护范围应以所附权利要求为准。

Claims (4)

1.一种基于扰动状态概念的结构面剪切本构模型的构建方法，其特征在于，包括如下步骤：
基于扰动状态概念理论，定义扰动函数DF，建立RI响应及FA响应相对于观测响应的贡献；其中，FA响应相对于观测响应的贡献通过DF来进行表示，RI响应相对于观测响应的贡献通过1-DF来进行表示，0≤DF≤1；RI响应指代为处于RI状态的材料剪切力学响应，即初始响应；FA响应指代为处于FA状态的材料剪切力学响应，即临界响应；观测响应指代为观测到的岩石结构面材料整体的剪切力学响应；
结合典型结构面剪切变形本构关系，分别获取初始响应及临界响应；
通过剪切试验观测量表征扰动函数DF，建立基于扰动状态概念的结构面剪切变形本构模型；
确定基于扰动状态概念的结构面剪切变形本构模型中的DSC参数；
将结构面剪切试验数据代入基于扰动状态概念的结构面剪切本构模型中，验证得到的结构面剪切本构模型能够准确模拟结果面剪切变形行为；
其中，所述步骤基于扰动状态概念理论，定义扰动函数DF，建立RI响应及FA响应相对于观测响应的贡献之前，还包括，
定义岩石结构面材料整体由无数个岩石单元组成，其中，岩石单元总数为l，l→∞，单个岩石单元的面积是A0，剪切面积A可以表示为A＝lA0，任意时刻处于RI状态的岩石单元数量为m，处于FA状态的岩石单元数量为n，l＝m+n；
获取力学平衡方程τA＝τRIARI+τFAAFA，其中，τ和A代表剪切应力和对应的剪切面积，τRI指的是处于RI状态下岩石单元所承受的剪切应力，ARI是处于RI状态下岩石单元的总面积，τFA表示处于FA状态下岩石单元所承受的剪切应力，AFA是处于FA状态下岩石单元的总面积；
获取剪切面积A、处于RI状态下岩石单元的总面积ARI、处于FA状态下岩石单元的总面积AFA：
Figure FDA0003605021510000021
将
Figure FDA0003605021510000022
代入到力学平衡方程τA＝τRIARI+τFAAFA中，并同时除以lA0后，得到剪切应力τ：
Figure FDA0003605021510000023
获取扰动函数
Figure FDA0003605021510000024
将表达式
Figure FDA0003605021510000025
及表达式l＝m+n代入到表达式
Figure FDA0003605021510000026
中，剪切应力τ的表达式变为τ＝τRI(1-DF)+τFADF；
所述步骤通过宏观观测量表征扰动函数DF，建立基于扰动状态概念的结构面剪切变形本构模型的方法，包括如下步骤，
随着剪切变形过程的推进，由RI状态变为FA状态的岩石单元数量也在增加，处于FA状态的岩石单元数量n与剪切位移δ的关系如下：
Figure FDA0003605021510000027
其中，r为比例参数，也为DSC参数，表示岩石单元的增长速率；
将表达式
Figure FDA0003605021510000028
两边求导得到
Figure FDA0003605021510000029
后，结合等式
Figure FDA00036050215100000210
及
Figure FDA00036050215100000211
获取dDF/dδ，
Figure FDA00036050215100000212
对
Figure FDA00036050215100000213
进行积分后可得到扰动函数DF的表达式为：
Figure FDA00036050215100000214
其中，C1和C2是
Figure FDA00036050215100000215
积分过程中产生的积分参数；
定义参数
Figure FDA0003605021510000031
将表达式
Figure FDA0003605021510000032
代入到τ＝ksδ(1-DF)+τrDF中后，获得基于扰动状态概念的结构面剪切变形本构模型
Figure FDA0003605021510000033
其中，ks、τr为剪切力学参数，可通过剪切试验直接获取，r、η为DSC参数。
2.根据权利要求1所述的基于扰动状态概念的结构面剪切本构模型的构建方法，其特征在于，所述步骤结合典型结构面剪切变形本构关系，分别获取初始响应及临界响应的方法，具体操作过程包括：
结合典型结构面剪切变形本构关系曲线剪切应力τ-剪切位移δ，分别获取初始响应τ＝ksδ及临界响应τ＝τr，其中，ks指代为剪切刚度，τr指代为残余剪切强度。
3.根据权利要求1所述的基于扰动状态概念的结构面剪切本构模型的构建方法，其特征在于：所述步骤确定基于扰动状态概念的结构面剪切变形本构模型中的DSC参数的方法，具体操作包括，
根据峰值强度点(δp，τp)对应的剪切应力-剪切位移关系满足等式
Figure FDA0003605021510000034
并且在峰值强度点(δp，τp)处剪切应力τ相对于剪切位移δ的导数为0，获得如下公式：
Figure FDA0003605021510000035
Figure FDA0003605021510000036
联立上述两公式
Figure FDA0003605021510000037
及
Figure FDA0003605021510000041
后，获得基于扰动状态概念的结构面剪切变形本构模型中的DSC参数η和r的计算公式：
Figure FDA0003605021510000042
4.根据权利要求1所述的基于扰动状态概念的结构面剪切本构模型的构建方法，其特征在于，所述基于扰动状态概念理论，定义扰动函数DF，建立RI响应及FA响应相对于观测响应的贡献的方法，包括如下步骤：
在DSC理论框架中，将变形岩石单元设定为RI状态岩石单元及FA状态岩石单元的组合；
定义扰动函数DF，获取RI状态岩石材料单元的响应及FA状态岩石材料单元的响应分别相对于观测到的岩石结构面材料整体的剪切力学响应的贡献。
CN202110557607.XA 2021-05-21 2021-05-21 基于扰动状态概念的结构面剪切本构模型的构建方法 Active CN113408103B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110557607.XA CN113408103B (zh) 2021-05-21 2021-05-21 基于扰动状态概念的结构面剪切本构模型的构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110557607.XA CN113408103B (zh) 2021-05-21 2021-05-21 基于扰动状态概念的结构面剪切本构模型的构建方法

Publications (2)

Publication Number Publication Date
CN113408103A CN113408103A (zh) 2021-09-17
CN113408103B true CN113408103B (zh) 2022-05-31

Family

ID=77679241

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110557607.XA Active CN113408103B (zh) 2021-05-21 2021-05-21 基于扰动状态概念的结构面剪切本构模型的构建方法

Country Status (1)

Country Link
CN (1) CN113408103B (zh)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7526385B2 (en) * 2007-06-22 2009-04-28 Schlumberger Technology Corporation Method, system and apparatus for determining rock strength using sonic logging
CN107036912A (zh) * 2017-05-27 2017-08-11 中国建筑东北设计研究院有限公司 一种粗颗粒土非线性弹性本构模型的构建方法
US11556683B2 (en) * 2018-06-14 2023-01-17 The Government Of The United States Of America, As Represented By The Secretary Of The Navy Fibrous composite failure criteria with material degradation for finite element solvers
CN110135113B (zh) * 2019-06-05 2023-07-07 中南大学 考虑尺寸效应的岩石结构面损伤统计本构模型的构建方法

Also Published As

Publication number Publication date
CN113408103A (zh) 2021-09-17

Similar Documents

Publication Publication Date Title
Jiang et al. Calibration of the continuous surface cap model for concrete
Banerjee et al. Historical origin and recent development on normal directional impact models for rigid body contact simulation: a critical review
Nikravesh et al. Improved state-based peridynamic lattice model including elasticity, plasticity and damage
WO2023005754A1 (zh) 一种季冻区层理岩石中岩土工程结构的稳定性分析方法
Ghazlan et al. An incremental formulation for the linear analysis of thin viscoelastic structures using generalized variables
Pandey et al. Numerical modelling of infilled clay brick masonry under blast loading
CN111967121B (zh) 一种基于区间法的粘弹性介电弹性体不确定性准静态和非线性动力学分析方法
Grassl et al. A damage-plasticity model for the dynamic failure of concrete
CN113408103B (zh) 基于扰动状态概念的结构面剪切本构模型的构建方法
Tian et al. Investigating time-dependent behavior of rocks using kinematic-constraint-inspired non-ordinary state-based peridynamics
Ukritchon et al. Soil parameter optimization of the NGI-ADP constitutive model for Bangkok soft clay
CN116629052A (zh) 一种粘弹性非线性介电弹性体本构模型的数值模拟方法
Chaimoon et al. Crack propagation due to time-dependent creep in quasi-brittle materials under sustained loading
Attia et al. Nonlinear analysis of frictional thermo-viscoelastic contact problems using FEM
Zbiciak et al. One-dimensional rheological models of asphalt-aggregate mixtures
Pruška Comparison of geotechnic softwares-Geo FEM, Plaxis, Z-Soil
CN116011191B (zh) 一种表征真三轴下岩石蠕变启动和加速的模型构建方法
CN116680832B (zh) 一种考虑结构面蠕变特性的隧道围岩稳定性评价方法
Kumar et al. Plasticity based interface model for failure modelling of unreinforced masonry under cyclic loading
CN106339545A (zh) 一种基于非连续变形分析的岩体流变模拟方法
Lu et al. Dynamic stability analysis of arch dam abutment based on strength reduction method
Kazerani et al. Application of numerical modelling for large-scale underground excavation in foliated rock mass
Kłosowski et al. Numerical treatment of elasto viscoplastic shells in the range of moderate and large rotations
Rašeta et al. Use of Finite Element Method for Simulation of RC Beam Nonlinear Behavior
Kim Comparisons of overstress theory with an empirical model in creep prediction for cohesive soils

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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20240123

Address after: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province

Patentee after: Dragon totem Technology (Hefei) Co.,Ltd.

Country or region after: China

Address before: Yuelu District City, Hunan province 410083 Changsha Lushan Road No. 932

Patentee before: CENTRAL SOUTH University

Country or region before: China

TR01 Transfer of patent right