CN111006941A - 一种真三轴试验力转变形及峰后扰动自动化控制方法 - Google Patents

一种真三轴试验力转变形及峰后扰动自动化控制方法 Download PDF

Info

Publication number
CN111006941A
CN111006941A CN201911280699.0A CN201911280699A CN111006941A CN 111006941 A CN111006941 A CN 111006941A CN 201911280699 A CN201911280699 A CN 201911280699A CN 111006941 A CN111006941 A CN 111006941A
Authority
CN
China
Prior art keywords
stress
control
deformation
disturbance
force
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.)
Granted
Application number
CN201911280699.0A
Other languages
English (en)
Other versions
CN111006941B (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.)
Northeastern University China
Original Assignee
Northeastern University 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 Northeastern University China filed Critical Northeastern University China
Priority to CN201911280699.0A priority Critical patent/CN111006941B/zh
Publication of CN111006941A publication Critical patent/CN111006941A/zh
Application granted granted Critical
Publication of CN111006941B publication Critical patent/CN111006941B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/02Details
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058Kind of property studied
    • G01N2203/0069Fatigue, creep, strain-stress relations or elastic constants
    • G01N2203/0075Strain-stress relations or elastic constants
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/0202Control of the test
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02Details not specific for a particular testing method
    • G01N2203/025Geometry of the test
    • G01N2203/0256Triaxial, i.e. the forces being applied along three normal axes of the specimen

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明涉及真三轴岩石力学试验机控制技术领域,提供一种真三轴试验力转变形及峰后扰动自动化控制方法。首先编写控制对话框;然后输入力加载速率、力转变形控制速率、扰动时机应力值占峰值强度的百分比,并根据力加载速率对岩石试样进行加载;接着根据最小主应力方向变形数据Di及采样时间数据Ti计算最小主应力方向变形速率vi,在vi≥a时对岩石试样进行力转最小主应力方向变形控制;最后,根据最大主应力方向应力数据Sj计算最大主应力方向应力变化值detaj,在detaj<0时计算ratioj,在ratioj≥1‑b时对岩石试样施加峰后扰动。本发明能够提高力转变形及峰后扰动控制的自动化程度、精度和效率。

Description

一种真三轴试验力转变形及峰后扰动自动化控制方法
技术领域
本发明涉及真三轴岩石力学试验机控制技术领域,特别是涉及一种真三轴试验力转变形及峰后扰动自动化控制方法。
背景技术
在进行常规三轴试验或者真三轴岩石力学试验时,目前的加载控制方法主要有力控制、变形控制、力转变形控制三种。岩石力学压缩试验全过程采用力控制方法无法获得峰后的应力应变曲线。若全过程采用变形控制方法则试验时间过长。为了获得岩石峰后曲线同时缩短试验时长,一般采用力转变形控制方法,在峰前弹性段采用力控制等到塑性段转为变形控制。目前大部分力转变形控制的转变时机都是根据最小主应力方向的变形速率人为判断,即对最小主应变进行观察,等应变时间曲线斜率到达一定数值时按照该斜率进行变形控制。这种方法需要试验操作人员不停的对曲线进行求斜率操作,自动化程度较低。
同时随着岩石力学试验设备的发展,现有设备希望更多的模拟现场工程实际状况,对压缩到峰后的岩石进行不同波形的扰动。然而峰后在什么时机进行扰动还没有成形的方法,若采用人工观察峰后曲线判断扰动时机则误差较大、很难严格的按照我们预期的时机施加扰动。
发明内容
针对现有技术存在的问题,本发明提供一种真三轴试验力转变形及峰后扰动自动化控制方法,能够提高力转变形及峰后扰动控制的自动化程度、精度和效率。
本发明的技术方案为:
一种真三轴试验力转变形及峰后扰动自动化控制方法,其特征在于,包括下述步骤:
步骤1:在真三轴岩石力学试验机中编写控制对话框,所述控制对话框可输入的参数包括力加载速率、力转变形控制速率、扰动时机应力值占峰值强度的百分比;
步骤2:在控制对话框中输入力加载速率、力转变形控制速率、扰动时机应力值占峰值强度的百分比,控制器将力加载速率编辑成力加载指令发送给作动器、将力转变形控制速率赋值给参数a、将扰动时机应力值占峰值强度的百分比赋值给参数b;
步骤3:作动器根据接收到的力加载指令对岩石试样进行加载;
步骤4:进行力转变形控制:
步骤4.1:令i=1;
步骤4.2:利用变形传感器监测岩石试样在各个主方向的变形值并记录监测时间,提取最小主应力方向变形数据Di及采样时间数据Ti;其中,Di为第i次变形采样时的最小主应力方向变形值,Ti为第i次变形采样对应的时间;
步骤4.3:令i=i+1,重复步骤4.2;
步骤4.4:计算第i次变形采样时的最小主应力方向变形速率
Figure BDA0002316660920000021
若vi<a则转至步骤4.3,若vi≥a则控制器将最小主应力方向变形速率vi编辑成变形控制指令来控制最大主应力方向的作动器对岩石试样进行力转最小主应力方向变形控制,进入步骤5;
步骤5:进行峰后扰动控制:
步骤5.1:令j=1;
步骤5.2:利用应力传感器监测岩石试样在各个主方向的应力值并记录监测时间,提取最大主应力方向应力数据Sj;其中,Si为第j次应力采样时的最大主应力方向应力值;
步骤5.3:令Smax=Sj
步骤5.4:令j=j+1,利用应力传感器监测岩石试样在各个主方向的应力值并记录监测时间,提取最大主应力方向应力数据Sj
步骤5.5:计算第j次应力采样时的最大主应力方向应力变化值detaj=Sj-Smax,若detaj≥0,则令Smax=Sj,转至步骤5.4;若detaj<0,则:令
Figure BDA0002316660920000022
若ratioj≥1-b则控制器将Smax作为峰值强度并发送峰后扰动指令给作动器、作动器根据接收到的峰后扰动指令对岩石试样施加峰后扰动,若ratioj<1-b则转至步骤5.4。
本发明的有益效果为:
本发明根据实时采集的最小主应力方向变形数据计算实时的最小主应力方向变形速率,在最小主应力方向变形速率大于或等于预设的力转变形控制速率时对岩石试样进行力转最小主应力方向变形控制,实现了力转变形的自动控制。本发明根据实时采集的最大主应力方向应力数据计算实时的最大主应力方向应力变化值,并在变化值为负时计算变化值绝对值占历史最大应力的比值,根据预设的扰动时机应力值占峰值强度的百分比来判断扰动时机,能够自动识别峰值强度且在一次试验中自动、准确、定量地把握扰动的施加时机。
附图说明
图1为本发明的真三轴试验力转变形及峰后扰动自动化控制方法的流程图;
图2为本发明的真三轴试验力转变形及峰后扰动自动化控制方法中力转变形控制的流程图。
图3为本发明的真三轴试验力转变形及峰后扰动自动化控制方法中峰后扰动控制的流程图。
具体实施方式
下面将结合附图和具体实施方式,对本发明作进一步描述。
如图1所示,本发明的真三轴试验力转变形及峰后扰动自动化控制方法,包括下述步骤:
步骤1:在真三轴岩石力学试验机中编写控制对话框,所述控制对话框可输入的参数包括力加载速率、力转变形控制速率、扰动时机应力值占峰值强度的百分比。
本实施例中,真三轴岩石力学试验机为具备整体锻造的高刚度框架以及动静分离式框架结构的能够实现峰后扰动试验的试验机。真三轴岩石力学试验机包括控制器、作动器、变形传感器、应力传感器。
步骤2:在控制对话框中输入力加载速率、力转变形控制速率、扰动时机应力值占峰值强度的百分比,点击对话框上的确定按钮,控制器将力加载速率编辑成力加载指令发送给作动器、将力转变形控制速率赋值给参数a、将扰动时机应力值占峰值强度的百分比赋值给参数b。
步骤3:作动器根据接收到的力加载指令对岩石试样进行加载;
步骤4:如图2所示,进行力转变形控制:
步骤4.1:令i=1;
步骤4.2:利用变形传感器监测岩石试样在各个主方向的变形值并记录监测时间,提取最小主应力方向变形数据Di及采样时间数据Ti;其中,Di为第i次变形采样时的最小主应力方向变形值,Ti为第i次变形采样对应的时间;
步骤4.3:令i=i+1,重复步骤4.2;
步骤4.4:计算第i次变形采样时的最小主应力方向变形速率
Figure BDA0002316660920000031
若vi<a则转至步骤4.3,若vi≥a则控制器将最小主应力方向变形速率vi编辑成变形控制指令来控制最大主应力方向的作动器对岩石试样进行力转最小主应力方向变形控制,进入步骤5。本实施例中,力转变形控制速率为0.015mm/min即0.9mm/h。
步骤5:如图3所示,进行峰后扰动控制:
步骤5.1:令j=1;
步骤5.2:利用应力传感器监测岩石试样在各个主方向的应力值并记录监测时间,提取最大主应力方向应力数据Sj;其中,Si为第j次应力采样时的最大主应力方向应力值;
步骤5.3:令Smax=Sj
步骤5.4:令j=j+1,利用应力传感器监测岩石试样在各个主方向的应力值并记录监测时间,提取最大主应力方向应力数据Sj
步骤5.5:计算第j次应力采样时的最大主应力方向应力变化值detaj=Sj-Smax,若detaj≥0,则令Smax=Sj,转至步骤5.4;若detaj<0,则:令
Figure BDA0002316660920000041
若ratioj≥1-b则控制器将Smax作为峰值强度并发送峰后扰动指令给作动器、作动器根据接收到的峰后扰动指令对岩石试样施加峰后扰动,若ratioj<1-b则转至步骤5.4。
显然,上述实施例仅仅是本发明的一部分实施例,而不是全部的实施例。上述实施例仅用于解释本发明,并不构成对本发明保护范围的限定。基于上述实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,也即凡在本申请的精神和原理之内所作的所有修改、等同替换和改进等,均落在本发明要求的保护范围内。

Claims (1)

1.一种真三轴试验力转变形及峰后扰动自动化控制方法,其特征在于,包括下述步骤:
步骤1:在真三轴岩石力学试验机中编写控制对话框,所述控制对话框可输入的参数包括力加载速率、力转变形控制速率、扰动时机应力值占峰值强度的百分比;
步骤2:在控制对话框中输入力加载速率、力转变形控制速率、扰动时机应力值占峰值强度的百分比,控制器将力加载速率编辑成力加载指令发送给作动器、将力转变形控制速率赋值给参数a、将扰动时机应力值占峰值强度的百分比赋值给参数b;
步骤3:作动器根据接收到的力加载指令对岩石试样进行加载;
步骤4:进行力转变形控制:
步骤4.1:令i=1;
步骤4.2:利用变形传感器监测岩石试样在各个主方向的变形值并记录监测时间,提取最小主应力方向变形数据Di及采样时间数据Ti;其中,Di为第i次变形采样时的最小主应力方向变形值,Ti为第i次变形采样对应的时间;
步骤4.3:令i=i+1,重复步骤4.2;
步骤4.4:计算第i次变形采样时的最小主应力方向变形速率
Figure FDA0002316660910000011
若vi<a则转至步骤4.3,若vi≥a则控制器将最小主应力方向变形速率vi编辑成变形控制指令来控制最大主应力方向的作动器对岩石试样进行力转最小主应力方向变形控制,进入步骤5;
步骤5:进行峰后扰动控制:
步骤5.1:令j=1;
步骤5.2:利用应力传感器监测岩石试样在各个主方向的应力值并记录监测时间,提取最大主应力方向应力数据Sj;其中,Si为第j次应力采样时的最大主应力方向应力值;
步骤5.3:令Smax=Sj
步骤5.4:令j=j+1,利用应力传感器监测岩石试样在各个主方向的应力值并记录监测时间,提取最大主应力方向应力数据Sj
步骤5.5:计算第j次应力采样时的最大主应力方向应力变化值detaj=Sj-Smax,若detaj≥0,则令Smax=Sj,转至步骤5.4;若detaj<0,则:令
Figure FDA0002316660910000012
若ratioj≥1-b则控制器将Smax作为峰值强度并发送峰后扰动指令给作动器、作动器根据接收到的峰后扰动指令对岩石试样施加峰后扰动,若ratioj<1-b则转至步骤5.4。
CN201911280699.0A 2019-12-13 2019-12-13 一种真三轴试验力转变形及峰后扰动自动化控制方法 Active CN111006941B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911280699.0A CN111006941B (zh) 2019-12-13 2019-12-13 一种真三轴试验力转变形及峰后扰动自动化控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911280699.0A CN111006941B (zh) 2019-12-13 2019-12-13 一种真三轴试验力转变形及峰后扰动自动化控制方法

Publications (2)

Publication Number Publication Date
CN111006941A true CN111006941A (zh) 2020-04-14
CN111006941B CN111006941B (zh) 2021-04-13

Family

ID=70114436

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911280699.0A Active CN111006941B (zh) 2019-12-13 2019-12-13 一种真三轴试验力转变形及峰后扰动自动化控制方法

Country Status (1)

Country Link
CN (1) CN111006941B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112903456A (zh) * 2021-01-20 2021-06-04 东北大学 一种主应力可变换的真三轴加卸载试验方法
CN113567246A (zh) * 2021-07-26 2021-10-29 东北大学 一种快速测定岩石真三轴强度的方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102221504A (zh) * 2011-04-14 2011-10-19 中国矿业大学 一种测定破裂后岩石三轴强度的加载路径方法
CN102435503A (zh) * 2011-11-24 2012-05-02 长江水利委员会长江科学院 一种大型真三轴试验测试方法及测试设备
CN106526131A (zh) * 2016-11-29 2017-03-22 广西大学 一种模拟岩墙轴向微扰动触发型岩爆的实验方法
CN108661553A (zh) * 2018-05-11 2018-10-16 东北大学 一种工程岩体致裂用微波功率自适应调控系统及方法
CN110411842A (zh) * 2019-08-01 2019-11-05 东北大学 可获取脆岩峰后曲线和残留强度的真三轴实验装置及方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102221504A (zh) * 2011-04-14 2011-10-19 中国矿业大学 一种测定破裂后岩石三轴强度的加载路径方法
CN102435503A (zh) * 2011-11-24 2012-05-02 长江水利委员会长江科学院 一种大型真三轴试验测试方法及测试设备
CN106526131A (zh) * 2016-11-29 2017-03-22 广西大学 一种模拟岩墙轴向微扰动触发型岩爆的实验方法
CN108661553A (zh) * 2018-05-11 2018-10-16 东北大学 一种工程岩体致裂用微波功率自适应调控系统及方法
CN110411842A (zh) * 2019-08-01 2019-11-05 东北大学 可获取脆岩峰后曲线和残留强度的真三轴实验装置及方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HOUQIANG YANG 等: "Stability Control of a Goaf-Side Roadway under the Mining Disturbance of an Adjacent Coal Working Face in an Underground Mine", 《SUSTAINABILITY》 *
PING YANG 等: "Research on Unloading Deformation Characteristics and The Law of Mechanical Parameters Based on Unloading Experiments of Marble", 《2014 7TH INTERNATIONAL CONFERENCE ON INTELLIGENT COMPUTATION TECHNOLOGY AND AUTOMATION》 *
邵生俊 等: "一室四腔刚–柔加载机构真三轴仪的改进与强度试验—西安理工大学真三轴仪", 《岩土工程学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112903456A (zh) * 2021-01-20 2021-06-04 东北大学 一种主应力可变换的真三轴加卸载试验方法
CN113567246A (zh) * 2021-07-26 2021-10-29 东北大学 一种快速测定岩石真三轴强度的方法
CN113567246B (zh) * 2021-07-26 2022-09-09 东北大学 一种快速测定岩石真三轴强度的方法

Also Published As

Publication number Publication date
CN111006941B (zh) 2021-04-13

Similar Documents

Publication Publication Date Title
CN111006941B (zh) 一种真三轴试验力转变形及峰后扰动自动化控制方法
JP6836577B2 (ja) 研磨工具摩耗量予測装置、機械学習装置及びシステム
CN113988202B (zh) 一种基于深度学习的机械臂异常振动检测方法
CN107272412B (zh) 一种暂冲式风洞流场控制的辩识方法
CN114619292B (zh) 一种基于小波降噪和注意力机制融合gru网络的铣削加工刀具磨损监测方法
CN114326596A (zh) 使用自学习系统模拟在机床上进行机加工的方法和设备
CN111832624A (zh) 一种基于对抗迁移学习的刀具剩余寿命预测方法
CN111639749A (zh) 一种基于深度学习的工业机器人摩擦力辨识方法
CN109815855B (zh) 一种基于机器学习的电子设备自动测试方法及系统
CN110631849B (zh) 一种摩擦系统磨损状态在线识别方法及系统
CN115628930B (zh) 一种基于rbf神经网络的掘进机井下截割工况预测方法
CN110889091A (zh) 基于温度敏感区间分段的机床热误差预测方法与系统
CN114398049A (zh) 一种离散制造车间数字孪生模型自适应动态更新方法
CN113534741A (zh) 一种铣削加工薄壁件的控制方法及系统
WO2020175418A1 (ja) 予測制御開発装置、予測制御開発方法、及び予測制御開発プログラム
CN110631792A (zh) 基于卷积神经网络的抗震混合试验模型更新方法
CN114253219A (zh) 一种基于端面磨削的磨削力自适应控制方法及系统
CN114936428A (zh) 一种基于adams的转向机构多体动力学建模及参数确认方法
CN116599857B (zh) 一种适用于物联网多场景的数字孪生应用系统
CN109739178A (zh) 一种监测数据驱动的工件加工变形在线预测方法
CN113326774A (zh) 一种基于AlexNet网络的机床能耗状态识别方法及系统
CN112475410A (zh) 一种铣削温度与多元影响因子的关联分析系统及方法
CN111045415A (zh) 一种基于局部概率密度双子空间的多模态过程故障检测方法
KR20210060147A (ko) 구조물 성능기반 재료의 역설계 방법
CN109886420B (zh) 一种自适应的采煤机截割高度智能预测系统

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