当前位置: 首页 > 专利查询>中南大学专利>正文

界面裂纹的热力耦合应力强度因子计算方法及系统技术方案

技术编号:39653871 阅读:12 留言:0更新日期:2023-12-09 11:22
本发明专利技术公开了一种界面裂纹的热力耦合应力强度因子计算方法,包括获取待分析的含界面裂纹双材料板的材料参数数据信息;根据裂纹面的类型计算得到对应的应力场和位移场;根据边界条件和连续性条件对应力场和位移场中的参数进行二次计算;根据功的互等定理完成界面裂纹的热力耦合应力强度因子的计算

【技术实现步骤摘要】
界面裂纹的热力耦合应力强度因子计算方法及系统


[0001]本专利技术属于土木工程
,具体涉及一种界面裂纹的热力耦合应力强度因子计算方法及系统


技术介绍

[0002]随着经济技术的发展,风力发电

水利发电等可再生能源不断被开发

为了解决可再生能源利用的间歇性问题,近年来,压气储能技术得到了较快的发展,也建立了一定数量的地下储气库

[0003]在建设地下储气库的过程中,密封层

衬砌层及围岩层受高温高压气体的影响,其间界面易产生裂纹并不断向外扩展,对地下储气库的安全性及稳定性产生严重影响

因此,有必要对热

力耦合作用下,界面裂纹的应力强度因子计算开展研究,从而为后续的热

力耦合下界面裂纹起裂判据的建立提供基础数据

[0004]目前,已经有学者提出了一些界面裂纹的应力强度因子计算方法

但是,该类方法一般都没有考虑到热力耦合因素,这使得目前的方案的可靠性不高,而且精确性较差


技术实现思路

[0005]本专利技术的目的之一在于提供一种可靠性高且精确性好的界面裂纹的热力耦合应力强度因子计算方法

[0006]本专利技术的目的之二在于提供一种实现所述界面裂纹的热力耦合应力强度因子计算方法的系统

[0007]本专利技术提供的这种界面裂纹的热力耦合应力强度因子计算方法,包括如下步骤:
S1. 获取待分析的含界面裂纹双材料板的材料参数数据信息;
S2. 根据步骤
S1
获取的数据信息,根据裂纹面的类型,计算得到对应的应力场和位移场;
S3. 根据边界条件和连续性条件,对步骤
S2
计算得到的应力场和位移场中的参数,进行二次计算;
S4. 根据功的互等定理,基于步骤
S2
和步骤
S3
的计算结果,完成界面裂纹的热力耦合应力强度因子的计算

[0008]所述的步骤
S2
,具体包括如下步骤:根据步骤
S1
获取的数据信息,将裂纹面分为绝热裂纹面和非绝热裂纹面;根据材料的剪切模量和材料的导热系数,计算绝热裂纹面的应力场和位移场;根据材料的剪切模量,计算非绝热裂纹面的应力场和位移场

[0009]所述的步骤
S2
,具体包括如下内容:根据步骤
S1
获取的数据信息,将裂纹面分为绝热裂纹面和非绝热裂纹面;针对绝热裂纹面,采用如下算式计算得到应力场和位移场:
;针对非绝热裂纹面,采用如下算式计算得到应力场和位移场:针对非绝热裂纹面,采用如下算式计算得到应力场和位移场:针对非绝热裂纹面,采用如下算式计算得到应力场和位移场:针对非绝热裂纹面,采用如下算式计算得到应力场和位移场:针对非绝热裂纹面,采用如下算式计算得到应力场和位移场:式中为第
j
种材料的径向应力;
Re
表示复变函数的实数部分;为第一奇异性指数;
r
为极坐标系的级径;
j
取值为1或2,分别对应
于第一种材料和第二种材料;为第
j
种材料的第一不定权函数;为极坐标系的级角;为第
j
种材料的第二不定权函数;为第
j
种材料的第三不定权函数;为第
j
种材料的第四不定权函数;为第
j
种材料的剪切模量;为第二奇异性指数;为第一热力耦合系数;为第
j
种材料的热膨胀系数;为第
j
种材料的第一温度不定权函数;为第
j
种材料的第二温度不定权函数;为第二热力耦合系数;为界面处温度;为第
j
种材料的环向应力;为第
j
种材料的剪切应力;为第
j
种材料的径向位移;为第
j
种材料由泊松比计算得到的系数;为第
j
种材料的环向位移;的计算公式为;的计算公式为;的计算公式为,其中为第
j
种材料的材料泊松比

[0010]所述的步骤
S3
,具体包括如下步骤:根据边界条件和连续性条件,对步骤
S2
中的第一奇异性指数

第二奇异性指数和不定权函数进行求解;针对第一奇异性指数和第二奇异性指数的求解:根据边界条件及连续性条件,以步骤
S2
中构建的应力场及位移场均有唯一解

且构成的行列式的值为0为目标进行求解,得到第一奇异性指数和第二奇异性指数;针对不定权函数的求解:将边界条件及连续性条件导入步骤
S2
中构建的应力场及位移场表达式,得到不定权函数的关系式;将不定全函数的关系式再代入步骤
S2
中构建的应力场及位移场表达式,并结合应力强度因子的定义表达式,从而对不定权函数进行求解

[0011]所述的步骤
S3
,具体包括如下内容:根据边界条件和连续性条件,对步骤
S2
中的第一奇异性指数

第二奇异性指数和不定权函数
、、、、
和进行求解;针对第一奇异性指数和第二奇异性指数的求解:根据边界条件及连续性条件,以步骤
S2
中构建的应力场有唯一解

步骤
S2
中构建的位移场有唯一解

且构成的行列式的值为0为目标进行求解,得到第一奇异性指数和第二奇异性指数;针对不定权函数的求解:所述的步骤
S3
,具体包括如下步骤:根据边界条件和连续性条件,对步骤
S2
中的第一奇异性指数

第二奇异性指数和不定权函数进行求解;针对第一奇异性指数和第二奇异性指数的求解:根据边界条件及连续性条件,以步骤
S2
中构建的应力场及位移场均有唯一解

且构成的行列式的值为0为目标进行求解,得到第一奇异性指数和第二奇异性指数;针对不定权函数的求解:将边界条件及连续性条件导入步骤
S2
中构建的应力场及
位移场表达式,得到不定权函数的关系式;将不定全函数的关系式再代入步骤
S2
中构建的应力场及位移场表达式,并结合应力强度因子的定义表达式,从而对不定权函数进行求解

[0012]所述的步骤
S3
,具体包括如下内容:根据边界条件和连续性条件,对步骤
S2
中的第一奇异性指数

第二奇异性指数和不定权函数
、、、、
和进行求解;针对第一奇异性指数和第二奇异性指数的求解:根据边界条件及连续性条件,以步骤
S2
中构建的应力场有本文档来自技高网
...

【技术保护点】

【技术特征摘要】
1.
一种界面裂纹的热力耦合应力强度因子计算方法,其特征在于包括如下步骤:
S1. 获取待分析的含界面裂纹双材料板的材料参数数据信息;
S2. 根据步骤
S1
获取的数据信息,根据裂纹面的类型,计算得到对应的应力场和位移场;
S3. 根据边界条件和连续性条件,对步骤
S2
计算得到的应力场和位移场中的参数,进行二次计算;
S4. 根据功的互等定理,基于步骤
S2
和步骤
S3
的计算结果,完成界面裂纹的热力耦合应力强度因子的计算
。2.
根据权利要求1所述的界面裂纹的热力耦合应力强度因子计算方法,其特征在于所述的步骤
S2
,具体包括如下步骤:根据步骤
S1
获取的数据信息,将裂纹面分为绝热裂纹面和非绝热裂纹面;根据材料的剪切模量和材料的导热系数,计算绝热裂纹面的应力场和位移场;根据材料的剪切模量,计算非绝热裂纹面的应力场和位移场
。3.
根据权利要求2所述的界面裂纹的热力耦合应力强度因子计算方法,其特征在于所述的步骤
S2
,具体包括如下内容:根据步骤
S1
获取的数据信息,将裂纹面分为绝热裂纹面和非绝热裂纹面;针对绝热裂纹面,采用如下算式计算得到应力场和位移场:针对绝热裂纹面,采用如下算式计算得到应力场和位移场:针对绝热裂纹面,采用如下算式计算得到应力场和位移场:针对绝热裂纹面,采用如下算式计算得到应力场和位移场:针对绝热裂纹面,采用如下算式计算得到应力场和位移场:;针对非绝热裂纹面,采用如下算式计算得到应力场和位移场:
式中为第
j
种材料的径向应力;
Re
表示复变函数的实数部分;为第一奇异性指数;
r
为极坐标系的级径;
j
取值为1或2,分别对应于第一种材料和第二种材料;为第
j
种材料的第一不定权函数;为极坐标系的级角;为第
j
种材料的第二不定权函数;为第
j
种材料的第三不定权函数;为第
j
种材料的第四不定权函数;为第
j
种材料的剪切模量;为第二奇异性指数;为第一热力耦合系数;为第
j
种材料的热膨胀系数;为第
j
种材料的第一温度不定权函数;为第
j
种材料的第二温度不定权函数;为第二热力耦合系数;为界面处温度;为第
j
种材料的环向应力;为第
j
种材料的剪切应力;为第
j
种材料的径向位移;为第
j
种材料由泊松比计算得到的系数;为第
j
种材料的环向位移;的计算公式为;的计算公式为;的计算公式为,其中为第
j
种材料的材料泊松比
。4.
根据权利要求3所述的界面裂纹的热力耦合应力强度因子计算方法,其特征在于所述的步骤
S3
,具体包括如下步骤:根据边界条件和连续性条件,对步骤
S2
中的第一奇异性指数

第二奇异性指数和不定权函数进行求解;针对第一奇异性指数和第二奇异性指数的求解:根据边界条件及连续性条件,以步骤
S2
中构建的应力场及位移场均有唯一解

且构成的行列式的值为0为目标进行求解,得到
第一奇异性指数和第二奇异性指数;针对不定权函数的求解:将边界条件及连续性条件导入步骤
S2
中构建的应力场及位移场表达式,得到不定权函数的关系式;将不定全函数的关系式再代入步骤
S2
中构建的应力场及位移场表达式,并结合应力强度因子的定义表达式,从而对不定权函数进行求解
。5.
根据权利要求4所述的界面裂纹的热力耦合应力强度因子计算方法,其特征在于所述的步骤
S3
,具体包括如下内容:根据边界条件和连续性条件,对步骤
S2
中的第一奇异性指数

第二奇异性指数和不定权函数
、...

【专利技术属性】
技术研发人员:饶秋华黄典一沈恺马燕易威
申请(专利权)人:中南大学
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1