CN103942387A - 一种基于变分数阶导数建立岩石蠕变本构模型的新方法 - Google Patents
一种基于变分数阶导数建立岩石蠕变本构模型的新方法 Download PDFInfo
- Publication number
- CN103942387A CN103942387A CN201410157685.0A CN201410157685A CN103942387A CN 103942387 A CN103942387 A CN 103942387A CN 201410157685 A CN201410157685 A CN 201410157685A CN 103942387 A CN103942387 A CN 103942387A
- Authority
- CN
- China
- Prior art keywords
- creep
- rock
- alpha
- fractional
- derivative
- 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
Links
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明公开了一种基于变分数阶导数建立岩石蠕变本构模型的新方法。本方法它包括以下步骤:建立变分数阶导数蠕变模型;获取实验数据;参数拟合。基于本方法建立的模型可以描述全过程蠕变曲线,尤其是非线性加速蠕变特征。本发明中的变分数阶导数可以充分体现岩石在蠕变过程中性质会随时间改变这一特征,并且应用变分数阶导数建立的岩石蠕变本构模型更加简单有效,可以准确的判断和预测岩体的加速变形,做出合理的防护措施,减少工程损失,具有巨大的工程应用价值和经济效益。
Description
技术领域:
本发明属于岩石工程技术领域,具体涉及一种基于变分数阶导数建立岩石蠕变本构模型的新方法。
背景技术:
岩石蠕变是指在常值应力持续作用下,岩体变形随时间而持续增长发展的过程。
岩石的蠕变性能是决定工程围岩长期稳定性的重要因素。如软岩流变常常会引起如岩体滑坡、地基失稳、坝基开裂以及矿柱岩爆等自然灾害。而通过岩石蠕变本构模型可准确判断和预测岩体的变形,做出合理的防护措施。因此,对岩石蠕变本构模型的研究是大型岩石工程必不可缺的重要内容,其重要性不言而喻。
由于分数阶导数算子具有记忆效应,近年来不少研究者提出了基于分数阶导数的岩石蠕变本构模型。但是常见的分数阶模型往往都很难反映岩石的全过程蠕变曲线,尤其是非线性加速蠕变阶段。
发明内容:
本发明的目的在于克服现有技术的不足,提供一种基于变分数阶导数建立岩石蠕变本构模型的新方法。
为了解决背景技术所存在的问题,本发明采用以下技术方案:
一种基于变分数阶导数建立岩石蠕变本构模型的新方法,它包括以下步骤:
步骤一:建立变分数阶导数蠕变模型:假设分数阶蠕变模型中的分数阶阶数是一个关于时间的函数,建立起变分数阶导数蠕变模型;
步骤二:获取实验数据:通过室内试验获得岩石出现加速蠕变特征的蠕变数据,进而根据蠕变实验数据绘出蠕变率曲线,从曲线上找出岩石进入加速蠕变阶段的时刻;
步骤三:参数拟合:根据岩石进入加速蠕变阶段的时刻,对变分数阶蠕变模型进行分段拟合,确定出蠕变本构模型的参数值。
进一步的,在步骤一中,分数阶阶数是一个关于时间的函数,可按以下步骤确定函数的形式:
(1)通过室内试验等技术手段得到岩石的蠕变实验曲线;
(2)再根据岩石蠕变实验曲线来确定函数的具体形式:假设分数阶阶数在出现加速蠕变的时刻之前是一个特定的值,出现加速蠕变的时刻之后变成另一个特定的值,这样就确定了函数的形式为分段阶跃函数;
基于Caputo导数定义,分数阶Maxwell由胡克体和Abel粘壶串联而成,
Abel粘壶(A)的应力应变关系为
其中表示对ε(t)的β阶分数阶导数,定义如下
式中:n为大于β的最小正整数;f(n)(τ)为函数f(τ)的n阶导数;
则式中σ为Abel粘壶的应力,即总应力;η0为其粘性系数;
认为Abel粘壶的分数阶阶数是一个关于时间的函数,即β=α(t),那么式中α(t)中为关于时间t的分数阶函数,ηα(t)为对应的粘性系数;
胡克体(H)的应力应变关系为式中σ为胡克体的应力,即总应力;E0为胡克体中弹簧的弹性模量;
综合考虑两部分应变,则变分数阶导数蠕变模型的本构方程可表示为:
进一步的,在步骤二中,假设α(t)在tc点以前是个特定的值,在tc点以后变成另一个特定的值,即:
式中tc为岩石进入加速蠕变的时刻;
其蠕变本构方程为: 其中
进一步的,在步骤三中,利用非线性最小二乘法拟合出蠕变本构方程:中的E0、ηα1、ηα2、α1、α25个参数。
本发明对比现有技术,有如下的有益效果:基于本方法建立的模型可以描述全过程蠕变曲线,尤其是非线性加速蠕变特征。本发明中的变分数阶导数可以充分体现岩石在蠕变过程中性质会随时间改变这一特征,并且应用变分数阶导数建立的岩石蠕变本构模型更加简单有效,可以准确的判断和预测岩体的加速变形,做出合理的防护措施,减少工程损失,具有巨大的工程应用价值和经济效益。
附图说明:
图1是为分数阶Maxwell模型图。
图2是岩石的典型蠕变曲线。
图3是通过室内试验得到的盐岩蠕变曲线。
图4是根据蠕变数据得到的蠕变率曲线。
图5是利用变分数阶导数蠕变模型对试验数据的拟合结果。
图6是利用分数阶导数蠕变模型对试验数据的拟合结果。
图7是本发明方法流程图。
具体实施方式:
下面结合附图和具体实施方式对本发明作进一步描述:
图7是本发明方法流程图。一种基于变分数阶导数建立岩石蠕变本构模型的新方法,它包括以下步骤:
步骤一:建立变分数阶导数蠕变模型:假设分数阶蠕变模型中的分数阶阶数是一个关于时间的函数,建立起变分数阶导数蠕变模型;
步骤二:获取实验数据:通过室内试验获得岩石出现加速蠕变特征的蠕变数据,进而根据蠕变实验数据绘出蠕变率曲线,从曲线上找出岩石进入加速蠕变阶段的时刻;
步骤三:参数拟合:根据岩石进入加速蠕变阶段的时刻,对变分数阶蠕变模型进行分段拟合,确定出蠕变本构模型的参数值。
方法原理:
由于常见的分数阶导数蠕变模型能很好的描述减速蠕变和稳定蠕变阶段特征,却不能反映加速蠕变特征。那么我们就可以假设分数阶阶数在出现加速蠕变的时刻之前是一个特定的值,出现加速蠕变的时刻之后变成另一个特定的值,这样就确定了函数的形式为分段阶跃函数。
确定求导阶数的函数形式关键在于判断岩石进入加速蠕变的时间点。因为在稳定蠕变阶段其蠕变率是几乎不变的,到了加速阶段其蠕变率会突然增大,所以可根据蠕变实验数据绘出蠕变率曲线,从曲线上找出岩石进入加速蠕变阶段的时刻。
进一步的,在步骤一中,分数阶阶数是一个关于时间的函数,可按以下步骤确定函数的形式:
(1)通过室内试验等技术手段得到岩石的蠕变实验曲线;
(2)再根据岩石蠕变实验曲线来确定函数的具体形式:假设分数阶阶数在出现加速蠕变的时刻之前是一个特定的值,出现加速蠕变的时刻之后变成另一个特定的值,这样就确定了函数的形式为分段阶跃函数;
基于Caputo导数定义,分数阶Maxwell由胡克体和Abel粘壶串联而成,
Abel粘壶(A)的应力应变关系为
其中表示对ε(t)的β阶分数阶导数,定义如下
式中:n为大于β的最小正整数;f(n)(τ)为函数f(τ)的n阶导数;
则式中σ为Abel粘壶的应力,即总应力;η0为其粘性系数;
认为Abel粘壶的分数阶阶数是一个关于时间的函数,即β=α(t),那么式中α(t)中为关于时间t的分数阶函数,ηα(t)为对应的粘性系数;
胡克体(H)的应力应变关系为式中σ为胡克体的应力,即总应力;E0为胡克体中弹簧的弹性模量;
综合考虑两部分应变,则变分数阶导数蠕变模型的本构方程可表示为:
进一步的,在步骤二中,假设α(t)在tc点以前是个特定的值,在tc点以后变成另一个特定的值,即:
式中tc为岩石进入加速蠕变的时刻;
其蠕变本构方程为: 其中
进一步的,在步骤三中,利用非线性最小二乘法拟合出蠕变本构方程:中的E0、ηα1、ηα2、α1、α25个参数。
图1是为分数阶Maxwell模型图。将整数阶Maxwell模型中的牛顿粘壶变为Abel粘壶,就是分数阶Maxwell模型。
图2是岩石的典型蠕变曲线。岩石的典型蠕变曲线一般分为三阶段:1、减速蠕变阶段;2、稳定蠕变阶段;3、加速蠕变阶段
图3是通过室内试验得到的盐岩蠕变曲线。图3是通过室内试验获得盐岩蠕变数据,并绘制成而成的曲线。
图4是根据蠕变数据得到的蠕变率曲线。蠕变数据进行差分求导得到蠕变率曲线。
图5是利用变分数阶导数蠕变模型对试验数据的拟合结果,是根据表1中得到的参数值代入蠕变本构方程中绘出的蠕变理论曲线。从图5可以看出,由于采用了本文提出的变分数阶导数建立的蠕变模型能很好的反映盐岩加速蠕变特征。
表1是基于盐岩蠕变实验的参数拟合结果
(表1中σ表示实验中的应力大小、E0表示试件的弹性模量、ηα1表示蠕变第一阶段Abel粘壶的粘性系数、ηα2表示第二阶段Abel粘壶的粘性系数、α1表示第一阶段的分数阶导数阶数、α2表示第二阶段的分数阶导数阶数。)
图6是利用分数阶导数蠕变模型对试验数据的拟合结果。从图6可以看出,常用的分数阶蠕变模型是不能描述加速蠕变特征的。
本发明首次提出基于变分数阶导数建立岩石蠕变本构模型的新方法,并利用室内试验等技术手段找到岩石进入加速蠕变的时刻,确定了分数阶阶数的分段函数形式,从而建立了变分数阶导数蠕变模型,和常见的分数阶导数蠕变模型相比,基于变分数阶导数建立的蠕变模型能够很好的反映蠕变过程中的加速阶段;本发明中的变分数阶导数可以充分体现岩石在蠕变过程中性质会随时间改变这一特征;本发明提出基于变分数阶导数建立岩石蠕变本构模型的新方法,使得所建模型简单,参数确定容易,可以更准确地判断和预测岩体的变形,做出合理的防护措施,减少工程损失。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种基于变分数阶导数建立岩石蠕变本构模型的新方法,其特征在于,它包括以下步骤:
步骤一:建立变分数阶导数蠕变模型:假设分数阶蠕变模型中的分数阶阶数是一个关于时间的函数,建立起变分数阶导数蠕变模型;
步骤二:获取实验数据:通过室内试验获得岩石出现加速蠕变特征的蠕变数据,进而根据蠕变实验数据绘出蠕变率曲线,从曲线上找出岩石进入加速蠕变阶段的时刻;
步骤三:参数拟合:根据岩石进入加速蠕变阶段的时刻,对变分数阶蠕变模型进行分段拟合,确定出蠕变本构模型的参数值。
2.根据权利要求1所述的一种基于变分数阶导数建立岩石蠕变本构模型的新方法,其特征在于,在步骤一中,分数阶阶数是一个关于时间的函数,可按以下步骤确定函数的形式:
(1)通过室内试验等技术手段得到岩石的蠕变实验曲线;
(2)再根据岩石蠕变实验曲线来确定函数的具体形式:假设分数阶阶数在出现加速蠕变的时刻之前是一个特定的值,出现加速蠕变的时刻之后变成另一个特定的值,这样就确定了函数的形式为分段阶跃函数;
基于Caputo导数定义,分数阶Maxwell由胡克体和Abel粘壶串联而成,
Abel粘壶(A)的应力应变关系为
其中表示对ε(t)的β阶分数阶导数,定义如下
式中:n为大于β的最小正整数;f(n)(τ)为函数f(τ)的n阶导数;
则式中σ为Abel粘壶的应力,即总应力;η0为其粘性系数;
认为Abel粘壶的分数阶阶数是一个关于时间的函数,即β=α(t),那么式中α(t)中为关于时间t的分数阶函数,ηα(t)为对应的粘性系数;
胡克体(H)的应力应变关系为式中σ为胡克体的应力,即总应力;E0为胡克体中弹簧的弹性模量;
综合考虑两部分应变,则变分数阶导数蠕变模型的本构方程可表示为:
3.根据权利要求1所述的一种基于变分数阶导数建立岩石蠕变本构模型的新方法,其特征在于,在步骤二中,
假设α(t)在tc点以前是个特定的值,在tc点以后变成另一个特定的值,即:
式中tc为岩石进入加速蠕变的时刻;
其蠕变本构方程为: 其中
4.根据权利要求1所述的一种基于变分数阶导数建立岩石蠕变本构模型的新方法,其特征在于,在步骤三中,
利用非线性最小二乘法拟合出蠕变本构方程:中的E0、ηα1、ηα2、α1、α25个参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410157685.0A CN103942387A (zh) | 2014-04-16 | 2014-04-16 | 一种基于变分数阶导数建立岩石蠕变本构模型的新方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410157685.0A CN103942387A (zh) | 2014-04-16 | 2014-04-16 | 一种基于变分数阶导数建立岩石蠕变本构模型的新方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103942387A true CN103942387A (zh) | 2014-07-23 |
Family
ID=51190055
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410157685.0A Pending CN103942387A (zh) | 2014-04-16 | 2014-04-16 | 一种基于变分数阶导数建立岩石蠕变本构模型的新方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103942387A (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105528491A (zh) * | 2015-12-24 | 2016-04-27 | 北京农业信息技术研究中心 | 一种植物动态虚拟仿真方法及*** |
CN106339545A (zh) * | 2016-08-26 | 2017-01-18 | 中国水利水电科学研究院 | 一种基于非连续变形分析的岩体流变模拟方法 |
CN108009311A (zh) * | 2017-10-30 | 2018-05-08 | 合肥通用机械研究院 | 一种针对蠕变试验的蠕变本构模型参数识别方法 |
CN108982818A (zh) * | 2018-08-30 | 2018-12-11 | 东北石油大学 | 一种描述外载作用岩体微细观破裂与微裂隙形成的方法 |
CN109142083A (zh) * | 2018-09-05 | 2019-01-04 | 南京航空航天大学 | 一种变载荷历程下的蠕变损伤计算方法及模型 |
CN109271742A (zh) * | 2018-10-29 | 2019-01-25 | 成都师范学院 | 磁控忆阶元 |
CN110631908A (zh) * | 2019-09-17 | 2019-12-31 | 西安建筑科技大学 | 一种岩石单轴压缩全过程蠕变损伤模型的构建方法 |
CN111044351A (zh) * | 2019-12-16 | 2020-04-21 | 合肥通用机械研究院有限公司 | 一种基于dic技术的焊接接头蠕变变形预测方法 |
CN111967120A (zh) * | 2020-02-21 | 2020-11-20 | 河海大学常州校区 | 基于变阶数分数阶导数的纳米银烧结体剪切变形建模方法 |
CN112329224A (zh) * | 2020-10-29 | 2021-02-05 | 四川大学 | 一种岩石时效变形的预测方法及其应用 |
CN112485114A (zh) * | 2020-11-23 | 2021-03-12 | 河海大学 | 一种预测混凝土特慢蠕变的方法 |
CN112613212A (zh) * | 2020-12-22 | 2021-04-06 | 中国人民解放军陆军装甲兵学院 | 修正分数阶Maxwell锰铜基阻尼合金本构模型构建方法及装置 |
CN114065499A (zh) * | 2021-11-10 | 2022-02-18 | 西安建筑科技大学 | 岩石单轴蠕变全过程改进Maxwell模型及构建方法 |
CN114550834A (zh) * | 2022-01-26 | 2022-05-27 | 河海大学常州校区 | 一种基于变阶数分数阶导数的高聚物变形的模型构建方法 |
CN115640683A (zh) * | 2022-10-20 | 2023-01-24 | 青岛理工大学 | 一种岩土体蠕变预测方法及*** |
CN116011191A (zh) * | 2022-12-13 | 2023-04-25 | 广西大学 | 一种表征真三轴下岩石蠕变启动和加速的模型构建方法 |
CN116665820A (zh) * | 2023-06-06 | 2023-08-29 | 河海大学常州校区 | 考虑温度影响的非晶玻璃态聚合物的本构模型构建方法 |
CN117910278A (zh) * | 2024-03-15 | 2024-04-19 | 太原科技大学 | 氢燃料电池双极板超薄不锈钢变分数阶本构模型构建方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130090902A1 (en) * | 2010-06-28 | 2013-04-11 | Yao Yao | Method and System for Modeling Fractures in Ductile Rock |
-
2014
- 2014-04-16 CN CN201410157685.0A patent/CN103942387A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20130090902A1 (en) * | 2010-06-28 | 2013-04-11 | Yao Yao | Method and System for Modeling Fractures in Ductile Rock |
Non-Patent Citations (2)
Title |
---|
李彦青 等: "变分数阶微积分在描述材料力学性质演变方面的应用", 《三峡大学学报(自然科学版)》 * |
段志强 等: "含Caputo 分数阶导数的流变模型研究", 《北京力学会第十六届学术年会论文集》 * |
Cited By (30)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105528491A (zh) * | 2015-12-24 | 2016-04-27 | 北京农业信息技术研究中心 | 一种植物动态虚拟仿真方法及*** |
CN106339545A (zh) * | 2016-08-26 | 2017-01-18 | 中国水利水电科学研究院 | 一种基于非连续变形分析的岩体流变模拟方法 |
CN106339545B (zh) * | 2016-08-26 | 2019-05-10 | 中国水利水电科学研究院 | 一种基于非连续变形分析的岩体流变模拟方法 |
CN108009311A (zh) * | 2017-10-30 | 2018-05-08 | 合肥通用机械研究院 | 一种针对蠕变试验的蠕变本构模型参数识别方法 |
CN108009311B (zh) * | 2017-10-30 | 2020-07-07 | 合肥通用机械研究院有限公司 | 一种针对蠕变试验的蠕变本构模型参数识别方法 |
CN108982818B (zh) * | 2018-08-30 | 2021-04-06 | 东北石油大学 | 一种描述外载作用岩体微细观破裂与微裂隙形成的方法 |
CN108982818A (zh) * | 2018-08-30 | 2018-12-11 | 东北石油大学 | 一种描述外载作用岩体微细观破裂与微裂隙形成的方法 |
CN109142083A (zh) * | 2018-09-05 | 2019-01-04 | 南京航空航天大学 | 一种变载荷历程下的蠕变损伤计算方法及模型 |
CN109271742A (zh) * | 2018-10-29 | 2019-01-25 | 成都师范学院 | 磁控忆阶元 |
CN110631908A (zh) * | 2019-09-17 | 2019-12-31 | 西安建筑科技大学 | 一种岩石单轴压缩全过程蠕变损伤模型的构建方法 |
CN110631908B (zh) * | 2019-09-17 | 2021-10-22 | 西安建筑科技大学 | 一种岩石单轴压缩全过程蠕变损伤模型的构建方法 |
CN111044351A (zh) * | 2019-12-16 | 2020-04-21 | 合肥通用机械研究院有限公司 | 一种基于dic技术的焊接接头蠕变变形预测方法 |
CN111044351B (zh) * | 2019-12-16 | 2021-04-20 | 合肥通用机械研究院有限公司 | 一种基于dic技术的焊接接头蠕变变形预测方法 |
CN111967120A (zh) * | 2020-02-21 | 2020-11-20 | 河海大学常州校区 | 基于变阶数分数阶导数的纳米银烧结体剪切变形建模方法 |
CN111967120B (zh) * | 2020-02-21 | 2022-07-22 | 河海大学常州校区 | 基于变阶数分数阶导数的纳米银烧结体剪切变形建模方法 |
CN112329224A (zh) * | 2020-10-29 | 2021-02-05 | 四川大学 | 一种岩石时效变形的预测方法及其应用 |
CN112485114A (zh) * | 2020-11-23 | 2021-03-12 | 河海大学 | 一种预测混凝土特慢蠕变的方法 |
CN112485114B (zh) * | 2020-11-23 | 2021-09-07 | 河海大学 | 一种预测混凝土特慢蠕变的方法 |
CN112613212A (zh) * | 2020-12-22 | 2021-04-06 | 中国人民解放军陆军装甲兵学院 | 修正分数阶Maxwell锰铜基阻尼合金本构模型构建方法及装置 |
CN114065499A (zh) * | 2021-11-10 | 2022-02-18 | 西安建筑科技大学 | 岩石单轴蠕变全过程改进Maxwell模型及构建方法 |
CN114065499B (zh) * | 2021-11-10 | 2024-06-28 | 西安建筑科技大学 | 岩石单轴蠕变全过程改进Maxwell模型及构建方法 |
CN114550834A (zh) * | 2022-01-26 | 2022-05-27 | 河海大学常州校区 | 一种基于变阶数分数阶导数的高聚物变形的模型构建方法 |
CN114550834B (zh) * | 2022-01-26 | 2022-12-13 | 河海大学常州校区 | 一种基于变阶数分数阶导数的高聚物变形的模型构建方法 |
CN115640683A (zh) * | 2022-10-20 | 2023-01-24 | 青岛理工大学 | 一种岩土体蠕变预测方法及*** |
CN115640683B (zh) * | 2022-10-20 | 2024-02-23 | 青岛理工大学 | 一种岩土体蠕变预测方法及*** |
CN116011191A (zh) * | 2022-12-13 | 2023-04-25 | 广西大学 | 一种表征真三轴下岩石蠕变启动和加速的模型构建方法 |
CN116011191B (zh) * | 2022-12-13 | 2024-05-10 | 广西大学 | 一种表征真三轴下岩石蠕变启动和加速的模型构建方法 |
CN116665820A (zh) * | 2023-06-06 | 2023-08-29 | 河海大学常州校区 | 考虑温度影响的非晶玻璃态聚合物的本构模型构建方法 |
CN117910278A (zh) * | 2024-03-15 | 2024-04-19 | 太原科技大学 | 氢燃料电池双极板超薄不锈钢变分数阶本构模型构建方法 |
CN117910278B (zh) * | 2024-03-15 | 2024-05-14 | 太原科技大学 | 氢燃料电池双极板超薄不锈钢变分数阶本构模型构建方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103942387A (zh) | 一种基于变分数阶导数建立岩石蠕变本构模型的新方法 | |
CN115527635B (zh) | 一种适用于谱载荷下的疲劳裂纹扩展寿命的预测方法 | |
Shaat et al. | Size-dependent analysis of functionally graded ultra-thin films | |
CN115640683B (zh) | 一种岩土体蠕变预测方法及*** | |
CN111680434A (zh) | 包含裂纹尖端非奇异应力项的裂缝扩展路径数值计算方法 | |
Jacklin et al. | Finite-element modeling of a light-framed wood roof structure | |
Mac et al. | A constitutive model for time-dependent behavior of clay | |
CN104359771B (zh) | 一种岩体结构面峰值摩擦角的分维数确定方法 | |
Ayas et al. | Energy release rate during the cracking of composite materials | |
Hunnell et al. | Numerical simulation of fatigue crack growth behavior by crack-tip blunting | |
CN104318014A (zh) | 一种二元法模拟钢和混凝土组合构件的计算方法 | |
CN104765898B (zh) | 一种锥环与密封环的组合结构及其设计方法 | |
CN207987266U (zh) | 一种拉链窗口折边机模具 | |
CN207380958U (zh) | 一种物理模拟中剪应力施加装置 | |
Baojian et al. | Study on Damage Factor of ABAQUS Concrete Plastic Damage Model | |
CN207380956U (zh) | 一种等体积拉张应力模拟装置 | |
Wang | State of the art-Reanalysis: recent development review | |
Lu et al. | Study on compressing principle of the fresh municipal solid waste | |
Dlouhý et al. | Micromechanical aspects of transgranular and intergranular failure competition | |
Zhou et al. | Creep Damage Analytical Model of RCCD | |
Wang et al. | Research on Shape Optimization of Spillway Tunnel | |
Frydrýšek et al. | Hard rock disintegration process and its numerical solution | |
Yang | Shape optimum design of shear panel damper made of low yield steel | |
Yu et al. | Numerical simulation of metal cutting process based on ANSYS/LS-DYNA | |
Shen et al. | A flow stress model for AZ31B considering recrystallization kinetics |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20140723 |
|
WD01 | Invention patent application deemed withdrawn after publication |