本文研究热荷载下正交各向异性带型中的裂纹问题。利用傅里叶变换,将热边值问题转化为奇异积分方程。采用Lobatto-Chebyshev公式得到线性代数方程系统。给出了裂纹诱导的温度场的数值解,通过数值计算分析了裂缝表面温度差的变化规律,采用图型表明了裂纹位置和裂纹尺寸对温度场的影响。 This paper investigates the problem of a crack embedded in an orthotropic strip under a thermal loading. Using the Fourier transform, the thermal boundary value problem is transformed into a singular integral equation. The Lobatto-Chebyshev formula is used to further derive the linear system of algebraic equations. Through the numerical calculation, the temperature jump across the crack surface is obtained. Finally, the influence of crack position and crack size on the temperature field is reflected in the form of a graph.
廖珊莉,吴远波,张利花
广西大学数学与信息科学学院,广西 南宁
收稿日期:2018年10月7日;录用日期:2018年10月22日;发布日期:2018年10月29日
本文研究热荷载下正交各向异性带型中的裂纹问题。利用傅里叶变换,将热边值问题转化为奇异积分方程。采用Lobatto-Chebyshev公式得到线性代数方程系统。给出了裂纹诱导的温度场的数值解,通过数值计算分析了裂缝表面温度差的变化规律,采用图型表明了裂纹位置和裂纹尺寸对温度场的影响。
关键词 :正交各向异性带型,裂纹,傅里叶变换,奇异积分方程,Lobatto-Chebyshev公式
Copyright © 2018 by authors and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY).
http://creativecommons.org/licenses/by/4.0/
断裂力学是研究带有裂纹的物体在载荷的作用下裂纹扩展规律的一门学科 [
本文侧重研究裂纹位置、裂纹尺寸对正交各向异性带型裂纹温度场的影响,基于热断裂物理边值问题转化为奇异积分方程的思想,通过数值计算求解奇异积分方程,获得了与实际生活的断裂现象较吻合的数值结果,表明了该方法的有效性。本文的研究不仅给出了裂纹位置、裂纹尺寸对温度场的影响,还为热弹性场的分析提供了研究基础,具有一定的实际意义。
考虑如图1所示嵌入在正交各向异性带型中的Griffith裂纹问题,该裂纹与带型面平行。建立笛卡尔坐标系
q x = − λ x ∂ θ ∂ x , q y = − λ y ∂ θ ∂ y ,
其平衡方程为
∂ q x ∂ x + ∂ q y ∂ y = 0
从而得到以下的控制偏微分方程
图1. 正交异性带型中的裂纹
λ 2 ∂ 2 θ ∂ x 2 + ∂ 2 θ ∂ y 2 = 0 , (1)
其中 λ 2 = λ x / λ y 。
其热边界条件为:
θ ( I ) ( x , h I ) = θ 0 , θ ( I I ) ( x , − h I I ) = − θ 0 , x > 0 , (2)
θ ( I ) ( x , 0 ) = θ ( I I ) ( x , 0 ) , x > a , (3)
q y ( I ) ( x , 0 ) = q y ( I I ) ( x , 0 ) , x > a , (4)
q y ( I ) ( x , 0 ) = q y ( I I ) ( x , 0 ) = q c , 0 < x < a . (5)
我们将研究偏微分方程(1)在以上热物理边值条件下的数值解。综合考虑带型边界初始温度、裂纹位置、裂纹尺寸对温度场的影响。
傅里叶变换在求解微分方程中有着重要的应用。其解题基本思想是:通过傅里叶变换将复杂的方程化成简单的代数方程,然后求解代数方程,再对代数方程的解作反变换,就得到了原方程的解 [
在求解之前,我们先引进傅里叶余弦变换 [
f ¯ ( ξ ) = 2 π ∫ 0 + ∞ f ( x ) cos ( ξ x ) d x ,
f ( x ) = ∫ 0 + ∞ f ¯ ( ξ ) cos ( ξ x ) d ξ .
以及傅里叶变换的性质之一
f ( n ) ( x ) = ( − i ξ ) n f ¯ ( ξ ) .
同时,介绍如下求解常微分方程的特征值法。求解常微分方程
d 2 f d y 2 + p d f d y + q f = 0 , (6)
其特征方程为
m 2 + p m + q = 0.
当特征方程存在两个不同实数根,设为 m 1 , m 2 ,则该微分方程的通解可表为
f = c 1 e m 1 y + c 2 e m 2 y , (7)
当特征方程存在两个相同的实数根,即
f = c 1 e m 1 y + c 2 y e m 1 y , (8)
其中 c 1 , c 2 为任意常数。又由于指数函数与双曲函数之间的转换关系
e x = sinh ( x ) + cosh ( x ) ,
所以(7)、(8)可分别表示为
f = ∑ i = 1 2 c i [ cosh ( m i y ) + sinh ( m i y ) ] , (9)
f = ( c 1 + c 2 y ) [ cosh ( m 1 y ) + sinh ( m 1 y ) ] . (10)
微分方程(6)的通解(9)或(10)将运用在温度场的求解过程中。
假设 I 和 I I 分别表示在坐标系中 y > 0 和 y < 0 的区域。由于该裂纹关于 y 轴对称,所以本文只取 y 轴右侧的带型进行研究。为了满足温度场的物理边值条件,温度场的解由两部分组成 [
θ 1 ( j ) ( x , y ) = 2 π ∫ 0 ∞ [ A ( j ) ( ξ ) cosh ( λ ξ y ) + B ( j ) ( ξ ) sinh ( λ ξ y ) ] cos ( ξ x ) d ξ ,
其中 j = I , I I 。此外, θ 2 ( y ) 满足如下3个方程
∂ 2 θ 2 ( y ) ∂ y = 0 , (11)
θ 2 ( h I ) = θ 0 , (12)
θ 2 ( − h I I ) = − θ (13)
用特征方程法求解常微分方程(11),获得 θ 2 ( y ) 的通解
θ 2 ( y ) = A 0 y + B 0 ,
又根据(12)、(13),求得
A 0 = 2 θ 0 h I + h I I , B 0 = θ 0 ( h I I − h I ) h I + h I I .
所以,温度场的解为
θ ( j ) ( x , y ) = 2 π ∫ 0 ∞ [ A ( j ) ( ξ ) cosh ( λ ξ y ) + B ( j ) ( ξ ) sinh ( λ ξ y ) ] cos ( ξ x ) d ξ + A 0 y + B 0 . (14)
其中 A ( j ) ( ξ ) , B ( j ) ( ξ ) 为待求函数。将(14)应用于热边界条件(2),得 A ( j ) ( ξ ) 与 B ( j ) ( ξ ) 的关系
B ( j ) ( ξ ) = δ { j } coth ( λ ξ h i ) A ( j ) ( ξ ) ,
将上式代入(14)中,温度场的解可表示为
θ ( j ) ( x , y ) = 2 π ∫ 0 ∞ sinh [ λ ξ ( h j + δ { j } y ) ] sinh ( λ ξ h j ) A ( j ) ( ξ ) cos ( ξ x ) d ξ + A 0 y + B 0 , (15)
综上,
q y ( j ) ( x , y ) = − λ y δ { j } 2 π ∫ 0 ∞ cosh [ λ ξ ( h j + δ { j } y ) ] sinh ( λ ξ h j ) A ( j ) ( ξ ) cos ( ξ x ) d ξ − λ y A 0 . (16)
其中, δ { I } = − 1 , δ { I I } = 1 。
为了求解 A ( j ) ( ξ ) ,首先引进辅助函数
φ 1 ( x ) = ∂ [ θ ( I ) ( x , 0 ) − θ ( I I ) ( x , 0 ) ] ∂ x , (17)
根据温度场的解(15)得
θ ( I ) ( x , 0 ) = 2 π ∫ 0 ∞ A ( I ) ( ξ ) cos ( ξ x ) d ξ + B 0 , (18)
θ ( I I ) ( x , 0 ) = 2 π ∫ 0 ∞ A ( I I ) ( ξ ) cos ( ξ x ) d ξ + B 0 , (19)
将方程(18)、(19)代入(17)且由边值条件(3),可得
A ( I ) ( ξ ) − A ( I I ) ( ξ ) = − 1 ξ ∫ 0 a φ 1 ( s ) sin ( ξ s ) d s . (20)
然后,由条件(4)得
A ( I I ) ( ξ ) = − coth ( λ ξ h I ) coth ( λ ξ h I I ) A ( I ) ( ξ ) , (21)
联立(20)、(21)有
A ( I ) ( ξ ) = − 1 ξ ∫ 0 a coth ( λ ξ h I I ) coth ( λ ξ h I ) + coth ( λ ξ h I I ) φ 1 ( s ) sin ( ξ s ) d s , (22)
A ( I I ) ( ξ ) = 1 ξ ∫ 0 a coth ( λ ξ h I ) coth ( λ ξ h I ) + coth ( λ ξ h I I ) φ 1 ( s ) sin ( ξ s ) d s , (23)
将方程(22)、(23)分别代入方程(18)、(19),可得裂纹面的温度分布公式
θ ( I ) ( x , 0 ) = − 2 π ∫ 0 a [ ∫ 0 ∞ 1 ξ coth ( λ ξ h I I ) coth ( λ ξ h I ) + coth ( λ ξ h I I ) sin ( ξ s ) cos ( ξ x ) d ξ ] φ 1 ( s ) d s + B 0 ,
将(22)代入方程(16),并根据边界条件(5)得到以下积分方程
2 π ∫ 0 a [ ∫ 0 ∞ 1 ξ λ y λ ξ coth ( λ ξ h I ) coth ( λ ξ h I I ) coth ( λ ξ h I ) + coth ( λ ξ h I I ) sin ( ξ s ) cos ( ξ x ) d ξ ] φ 1 ( s ) d s = − ( q c + λ y A 0 ) . (24)
奇异积分方程在断裂力学方面有着重要的应用。对奇异积分方程数值解的研究,不少学者们已经取得了相对成熟的成果。例如,配位法、基于Fourier展开的Galekrin-Petrvo方法、样条方法等 [
为了方便求解奇异积分方程,假设
M ( ξ ) = λ y λ ξ coth ( λ ξ h I ) coth ( λ ξ h I I ) coth ( λ ξ h I ) + coth ( λ ξ h I I ) ,
所以,方程(24)可以写为
2 π ∫ 0 a [ ∫ 0 ∞ M ( ξ ) ξ sin ( ξ s ) cos ( ξ x ) d ξ ] φ 1 ( s ) d s = τ 0 . (25)
根据文献 [
2 π ∫ 0 ∞ sin ( ξ s ) cos ( ξ x ) d ξ = 1 π ( 1 s − x + 1 s + x ) ,
且假设
Λ = lim x → + ∞ M ( ξ ) ξ ,
Λ 是由材料性质决定的常数,可以通过数值计算解得。所以方程(25)可化为以下奇异积分方程
1 π ∫ − a a Λ s − x φ 1 ( s ) d s + 1 π ∫ − a a k ( s , x ) φ 1 ( s ) d s = τ 0 , (26)
其中,
k ( s , x ) = ∫ 0 ∞ [ M ( ξ ) ξ − Λ ] sin ( ξ s ) cos ( ξ x ) d ξ .
引进以下变量变换
x = a x ¯ , s = a s ¯ , φ 1 ( s ) = Y ( s ¯ ) 1 − s ¯ 2 ,
方程(26)可表示为
1 π ∫ − 1 1 Λ s ¯ − x ¯ Y ( s ¯ ) 1 − s ¯ 2 d s ¯ + 1 π ∫ − 1 1 k ¯ ( s ¯ , x ¯ ) Y ( s ¯ ) 1 − s ¯ 2 d s ¯ = τ 0 , (27)
其中,
k ¯ ( s ¯ , x ¯ ) = a ∫ 0 ∞ [ M ( ξ ) ξ − Λ ] sin ( ξ a s ¯ ) cos ( ξ a x ¯ ) d ξ .
由单值性条件(3)得
∫ − 1 1 φ 1 ( x ) d x = 0. (28)
一般的,奇异积分方程(27)可以通过数值计算方法求解。本文对方程(27)、(28)应用
Lobatto-Chebyshev积分公式,获得以下代数方程组
1 L ∑ l = 0 L λ l Λ s ¯ l − x ¯ m Y ( s ¯ l ) + 1 L ∑ l = 0 L λ l k ¯ ( s ¯ , x ¯ ) Y ( s ¯ l ) = τ 0 , (29)
∑ l = 0 L λ l Y ( s ¯ l ) = 0 , (30)
其中
x ¯ m = cos [ ( 2 m − 1 ) π / ( 2 L ) ] , m = 1 , 2 , ⋯ , L ,
s ¯ l = cos ( l π / L ) , l = 0 , 1 , 2 , ⋯ , L ,
λ 0 = λ L = 1 / 2 , λ 1 = ⋯ = λ L − 1 = 1.
在上式中的 L 取偶数。
只要代数方程(29)、(30)得以求解,那么裂纹上下表面的温度差就可以确定了。比如,
∫ x a φ 1 ( s ) d s = ∫ x a ∂ [ θ ( I ) ( s , 0 ) − θ ( I I ) ( s , 0 ) ] ∂ s d s = [ θ ( I ) ( s , 0 ) − θ ( I I ) ( s , 0 ) ] | a x = − [ θ ( I ) ( x , 0 ) − θ ( I I ) ( x , 0 ) ] ,
所以,
θ ( I ) ( x , 0 ) − θ ( I I ) ( x , 0 ) = − a ∫ x ¯ 1 Y ( s ¯ ) 1 − s ¯ 2 d s ¯ .
特别地,当 x ¯ = 0 ,
θ ( I ) ( x , 0 ) − θ ( I I ) ( x , 0 ) ≈ − a ∑ l = 0 L / 2 λ l Y ( s ¯ l ) . (31)
综上,我们便获得了裂纹上下表面温度差的数值结果。
例:为了体现本文方法的有效性,取陶瓷材料Tyrannohex [
根据公式(29)、(30)、(31),我们可以得到以下数值结果。当我们令裂纹面的热流密度 q c = 0 ,则该带型裂纹就变成热不渗透裂纹,此时只考虑裂纹位置、裂纹尺寸、带型边界初始温度对温度场的影响。图2中的 a / h 可表示裂纹相对于带型的尺寸, h I / h 表示裂纹相对于带型中心的位置,其中 Δ θ = θ ( I ) ( x , 0 ) − θ ( I I ) ( x , 0 ) 。通过观察图2,我们不难发现,若 h I / h 固定,当裂纹尺寸相对于带型宽度的比值越大,裂纹上下表面的温度差越大;若 a / h 固定,当裂纹逐渐偏离带型边界,越向带型中间位置靠拢,裂纹上下表面的温度差越大。此外,裂纹上下表面的温度差也与带型边界初始温度成正比。
在本文中,通过傅里叶余弦变换和数值计算,给出了裂纹面的温度分布公式以及裂纹上下表面温度
图2. 裂纹位置以及裂纹尺寸对 Δ θ / θ 0 的影响曲线
差的数值结果。并且通过数值结果验证了裂纹上下表面的温度差与裂纹位置、裂纹尺寸、带型边界初始温度有关。本文不仅为裂纹温度场的求解提供了新的思路,采用的方法能够广泛的应用于各种裂纹问题的求解,所得到的结果还为热弹性场的分析奠定了研究基础。
广西自然科学基金(2016GXNSFAA380261)。
廖珊莉,吴远波,张利花. 正交各向异性带型中裂纹诱导的温度场分析Analysis of the Temperature Field Induced by a Crack in an Orthotropic Strip[J]. 应用物理, 2018, 08(10): 439-447. https://doi.org/10.12677/APP.2018.810056