一种边坡渐进破坏全过程计算方法
专利汇可以提供一种边坡渐进破坏全过程计算方法专利检索，专利查询，专利分析的服务。并且本 发明 提出了一种边坡渐进破坏全过程计算方法，在假设地质材料全过程主应 力 -主应变满足相关本构关系条件下，计算材料全过程主 应力 -主应变的任意应力和应变，利用坐标转换与可逆应变和不可逆应变无关的特点，获得任意应力状态下的应力和应变值；在滑面破坏满足某一峰值应力准则条件下，由于破坏区摩阻力的下降，产生剩余下滑力，增加滑体前面单元的法向应力和沿滑面的驱动剪应力，致使边坡一步一步向前发生破坏；通过计算滑面上滑体单元应变，决定滑体与滑床的相对应变值，利用该应变值计算滑面摩阻力。本发明直接计算滑面剪应力和剪应变，不存在条 块 之间 变形 协调假设，也可以计算滑体的应力和应变。，下面是一种边坡渐进破坏全过程计算方法专利的具体信息内容。
1.一种边坡渐进破坏全过程计算方法，其特征在于，包括以下步骤：
当滑体和滑面为同种材料时，
步骤1：对材料实施主应力-应变全过程曲线试验，试验获得峰值主应力、应变和全过程曲线特征，计算滑体材料的物理力学特征方程，包括主应力-主应变方程、峰值应力方程、临界主应变方程、弹性模量变化方程、软化系数方程；
步骤1.1：主应力-主应变方程
主应力-主应变为四参数本构方程：
式中，σi、εi分别为主应力和主应变，Eii为初始弹性模量，mi、Si、ρi为主应力方向的常系数，σi、Eii的单位为MPa或kPa或Pa，mi、Si、ρi为无单位参数；
其行为描述为：对于具有软化特征的材料行为，则有-1<ρi≤0和1+miSi≠0；临界应变空间满足关系式：
式中，εipeak为临界主应力对应的主应变；
步骤1.2：假设临界主应力满足摩尔库伦准则，即
式中，C为凝聚力，C和σi的单位为MPa或kPa或Pa，φ为材料摩擦角；
步骤1.3：设临界主应变空间仅相关于第三主应力，临界主应变εipeak采用关系式：
或
式中，ai,1,ai,2,ai,3,ζi,N, 为常系数；ai,1,ai,2单位为MPa或kPa或Pa，ai,3,ζi,N为无量纲系数，或 的量纲为1/MPa,1/MPa2或1/kPa,1/kPa2或1/Pa,1/Pa2；
步骤1.4：初始弹性模量方程：
Eii＝Eii,0+bi,1σ3+bi,2σ32
式中，Eii,0为第三主应力σ3为零值的Eii值，bi,1,bi,2为常系数，单位分别为无量纲和1/MPa或1/kPa或1/Pa；
步骤1.5：对无量纲参数ρi，软化系数演化方程表示为：
式中，ρi,0为第三主应力σ3为零值的ρi值，ρi,c为σi,3等于σi,3c时的ρi值，为常系数；
步骤1.6：计算滑体主应变在滑面单元上产生的应变εij'表达式为：
主应变ε1的方向与x轴间所夹角度 按式子
和 求得，当γxy为正直，且εxx>εyy时， 角为负值；计算滑面上的单元的应力，计算公式为：
步骤2：在峰值主应力满足摩尔库伦准则条件下，其破坏面上剪应力与最小主应力的夹角为45°+φ/2，φ为材料摩擦角，单元剪应力摩尔库伦准则为：
τpeak＝C+σn tanφ
其中， τpeak,σn,C分别为破坏面
最大剪应力、法向应力和凝聚力，单位为MPa或kPa或Pa；
步骤3：滑面上的剪应力和剪应变满足的本构方程，具体如下：
步骤3.1：剪应力-剪应变为四参数本构方程：
τ＝Gγ[1+γq/p]ξ
式中，τ、γ分别为剪应力和剪应变，G为剪切模量，p、q、ξ为在不同法向应力下的常系数，τ、G的单位为MPa或kPa或Pa，p、q、ξ为无单位参数；
其行为描述为：对于具有软化特征的材料行为，则有-1<ξ≤0和1+qξ≠0；临界应变空间满足关系式：
p+(1+qξ)γqpeak＝0
式中，γpeak为临界应力对应的应变；
步骤3.2：滑体主应变在滑面上产生的法向应变εn'和剪应变γ'表达式：
γ'＝(ε1-ε3)sin2α
式中，α为滑面与最小主应力方向的夹角；
步骤3.3：初始剪切模量方程：
式中，ν为泊松比；
步骤3.4：对于无量纲参数ξ，软化系数演化方程表示为：
ξ＝ξ0/(1+(ξ0/ξc-1)(σn/σnc)ψ)
ξ0为法向应力σn为零值的ξ值，ξc为σn等于σnc时的ξ值，ψ为常系数；
对于同一种地质材料，直接利用单元在滑面上求解的驱动剪应力、剪应变和法向应力求解软化系数：求解方法为：
利用破坏面上所得到的法向应力，对应地求解出主应力σ1,σ3，方程为：
滑面上的法向应力：
临界剪应力方程：
求解得临界主应力和围压：
利用方程
获得滑面上临界剪应变、临界剪应力和对应的法向应力，计算常系数：
p＝-(γpeak)q(1+qξ)
再利用滑面上的实际计算剪应力，以方程τ＝Gγ[1+γq/p]ξ计算软化系数ξ；
若滑面材料与滑体不一致，则滑面计算方法为：
步骤4：对于不同于滑体材料滑面的剪应力-剪应变四参数本构方程、临界应力方程、临界剪应变方程和软化系数演化方程均与方程：
τ＝Gγ[1+γq/p]ξ
τpeak＝C+σntanφ
或 二者之一
ξ＝ξ0/(1+(ξ0/ξc-1)(σn/σnc)ψ)一致，其初始剪切模量方程为：
G＝G0+b1σn+b2σn2
式中，G0为法向应力σn为零值的G值，b1,b2为常系数，单位为无量纲和1/MPa或1/kPa或1/Pa；
步骤5：沿滑面的变形分布特征
当M条块或单元为临界状态时，沿滑面的力分布特征为：在M+1至N条块间，下滑力等于摩阻力，压应力等于反压应力，在这个阶段，沿滑面的剪应变γ以下滑应力等于抗滑剪应力加以计算；在M条块或单元下滑应力等于抗滑临界摩阻应力，以此计算该条块的临界剪应变γ；
在第1至M-1条块或单元间，下滑应力大于摩阻应力，其任一条块或单元L存在两个位移：第一种位移为第L条块或单元本身产生的位移，即：由第L-1的右边和第L+1的左边位移之差加以计算，第二种位移为第L条块或单元本身产生的位移加上第L条块或单元的底边中心相对该点沿滑床原始位置滑动的位移；第一种位移用于计算滑体的驱动下滑剪应力，第二种位移用于计算滑床提供给滑体的摩阻应力；
L条块或单元的第一种底边位移SL,j1计算：
1 L L
SL,j＝εjl，j∈(x,y,z)
L条块或单元的第二种底边位移SL,j2计算：
式中，εji,li分别为沿滑面第i条块或单元在x、y、z轴方向的应变和底边长度；在滑面方向的位移为三者的矢量和，滑面剪应变γ为上述第一、二种位移除以滑面长度；
步骤6：硬化特征描述
当滑体仅表现出硬化特征时，其计算采用简化的模型，将本构方程
取ξ＝-1，q＝1，则a’＝1/(Ga”)，b'＝1/(Gp)，在主应力条件下其方程形式与邓肯—张模型一致，此时描述材料的硬化行为特征：
式中，ai’、bi’为常系数；
在单元主应力满足方程 条件下，相对应的单元应变
定义为峰值应变，满足方程 或
在峰值应力和应变下，方程 变为：
定义单元割线模量 则
对方程 求导，相应的导数为切线模量，在任一应力状态条件下，切向
模量Ei表示为：
得到初始弹性模量为1/ai'，利用方程 在主应力满足方程
时的切线模量定义为Et,i，则有：
Et,i＝ai’Kscant,i2
在峰值应力时，研究试验曲线的切向模量为Et,i，设其具有如下特征：
σ3≠0
式中，αi,ki为常系数，pa为大气压；
方程 的常系数利用三个不同法向应力试验加以确定；
步骤7：针对条分法，边坡渐进破坏计算方法包括如下步骤：
步骤7.1：利用现行方法进行条块划分；
步骤7.2：竖向应力以比重与高度之积加以计算，水平应力和剪应力以剩余推力为水平方向和垂直于水平方向力的矢量和，水平方向力和垂直于水平方向力满足应力分布条件；
计算相对应的应力，同时计算主应力和旋转角，根据主应力计算对应主应变，通过坐标旋转计算滑动方向和其它方向的应变和位移；
步骤7.3：根据现场实测或假设确定第一条块水平和竖直方向的位移，该位移加上步骤
7.2计算的位移，所得位移除以滑面长度即为滑面剪应变或视剪应变，利用该剪应变计算条块底边摩阻应力；
步骤7.4：重复步骤7.2和步骤7.3，直至现场判断的临界状态条块为止；
步骤8：实施边坡渐进破坏的数值计算。
一种边坡渐进破坏全过程计算方法
技术领域
背景技术
发明内容
的应力，计算公式为：
系数：
方程形式与邓肯—张模型一致，此时描述材料的硬化行为特征：
在峰值应力和应变下，方程 变为：
具体实施方式
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种陶瓷基复合材料界面参数识别方法 | 2020-05-11 | 115 |
| 光缆元件 | 2020-05-17 | 541 |
| 吹胀薄膜 | 2020-05-20 | 97 |
| 用于滚塑和注塑制品的以茂金属催化剂制备的聚合物 | 2020-05-18 | 548 |
| 具有金属化非取向薄膜的阻隔包装金属薄片 | 2020-05-20 | 448 |
| 具有高强度、高起始模量和低缩率的聚酯拉伸丝 | 2020-05-12 | 655 |
| 用于滚塑和注塑制品的以茂金属催化剂制备的聚合物 | 2020-05-15 | 829 |
| 具有提高的机械性能和抗渗性的聚乙烯膜及其制造方法 | 2020-05-26 | 579 |
| 用于重型车辆的轮胎 | 2020-05-17 | 406 |
| 包含涂覆的水溶胀性聚合物的水溶胀性材料 | 2020-05-25 | 852 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
