CN114580160B 一种用于含裂纹功能梯度板传热问题的非局部方法 （河海大学）.docxVIP

(19)国家知识产权局

(12)发明专利

(10)授权公告号CN114580160B(45)授权公告日2025.01.10

(21)申请号202210174201.8

(22)申请日2022.02.24

(65)同一申请的已公布的文献号申请公布号CN114580160A

(43)申请公布日2022.06.03

(73)专利权人河海大学

地址210098江苏省南京市鼓楼区西康路1

号

(72)发明人黄丹周保良李志远

(74)专利代理机构南京苏高专利商标事务所(普通合伙)32204

专利代理师柏尚春

(51)Int.CI.

GO6F30/20(2020.01)

G16C60/00(2019.01)

GO6F17/13(2006.01)

GO6F17/16(2006.01)

G06F119/08(2020.01)

(56)对比文件

CN112116128A,2020.12.22

CN113486543A,2021.10.08审查员徐生芹

权利要求书3页说明书6页附图4页

(54)发明名称

一种用于含裂纹功能梯度板传热问题的非局部方法

(57)摘要

CN114580160B本发明公开了一种用于含裂纹功能梯度板传热问题的非局部方法，包括以下步骤：确定功能梯度材料的参数以及板的尺寸；确定初始裂纹位置；以合适的网格间距离散实体模型；确定功能梯度板温度场的控制方程并设置模型的边界条件；基于近场动力学微分算子理论将控制方程与边界条件转化为非局部形式，离散并构建求解方程组；求解方程组，得到每个物质点的温度值。对于含裂纹的传热工程构件，通过本发明对裂纹尖端处热通量奇异性的预测，可以有效阻止并减

CN114580160B

Tz

TzW

W

CN114580160B权利要求书1/3页

2

1.一种用于含裂纹功能梯度板传热问题的非局部方法，其特征在于，包括以下步骤：

S1、确定功能梯度材料的参数以及板的尺寸；

S2、确定初始裂纹位置；

S3、以合适的网格间距离散实体模型；

S4、确定功能梯度板温度场的控制方程并设置模型的边界条件；

S5、基于近场动力学微分算子理论将控制方程与边界条件转化为非局部形式，离散并构建求解方程组；

S6、求解方程组，得到每个物质点的温度值；

其中，S4中温度场的控制方程推导过程为：

稳态导热方程为：

▽·(k▽T)+Q=0

式中，T为温度，k为导热系数，Q为热源；

对于一个二维问题，假设功能梯度板的材料性能随x变化，即：

k=k(x)

在没有热源的情况下，功能梯度板的温度场控制方程为：

热传导中有以下三类边界条件：

Dirichlet边界条件：T(x,y)=T(x,y)

Neumann边界条件：-kVT·n=q(x,y)

第三类边界条件：-k▽T·n=h(T-T)

式中，T(x,y)为已知温度，q(x,y)为已知热通量，h为对流换热系数，T?为环境温度；

S5中控制方程离散的具体计算方法为：

对于一个二维空间中的场函数f(x)=f(x,y)在点x处进行N阶泰勒展开：

其中，R(N,x)指可忽略的余数，n?和n?表示对x和y的微分阶数，ξ=ξ1e?+5?e2,e?和e?表示x和y方向上的单位向量，ξ?和?表示两物质点之间长度在x和y方向上的投影；

利用正交函数的性质，用PD函数将偏微分转化为积分形式

其中1≤p?+p?≤N,H表示x点的近场范围，PD函数为：

CN114580160B权利要求书2/3页

3

其中是权函数，是未知系数矩阵，可由下式得到：

其中A(nn?)(992)是形状矩阵，是已知的系数矩阵；

S4中控制方程及涉及微分的Neumann边界条件和第三类边界条件可由其微分形式改写为积分形式如下：

其中T′=T(x′,y′),gPP2=gPIP2(ξ,52),将上式离散化得：

其中，含下标(k)和(j)的符号分别表示对