含裂紋結構的區(qū)間B樣條小波模糊有限元分析

彭惠芬，孟廣偉，周立明，李 鋒

（1.吉林大學機械科學與工程學院，130022 長春，phfdaqing@163.com;2.東北石油大學機械科學與工程學院，163318 大慶黑龍江）

含裂紋結構的區(qū)間B樣條小波模糊有限元分析

彭惠芬1，2，孟廣偉1，周立明1，李 鋒1

（1.吉林大學機械科學與工程學院，130022 長春，phfdaqing@163.com;2.東北石油大學機械科學與工程學院，163318 大慶黑龍江）

為了解決工程問題中材料和載荷的不確定性給含裂紋結構數(shù)值分析帶來的困難，提高數(shù)值計算的精度和效率，將模糊理論與小波有限元相結合，提出了基于區(qū)間B樣條小波(B-spline wavelet on the interval，BSWI)模糊有限元分析法.該方法將啞節(jié)點斷裂單元鑲嵌到含裂紋結構的小波有限元模型中，建立了含裂紋結構的小波有限元模型，推導了小波有限元模糊平衡方程，并采用λ水平截集及分解定理求解模糊平衡方程;在此基礎上，利用虛擬裂紋閉合法計算了不同裂紋長度下應力強度因子隸屬函數(shù)值，并將計算結果與解析解進行了比較.結果分析表明該方法可用較少單元更為真實、準確地反映結構響應的變化情況，為工程實際中實現(xiàn)含不確定參數(shù)的斷裂數(shù)值分析提供了一種新途徑.

區(qū)間B樣條小波;裂紋;模糊數(shù);應力強度因子;有限元

現(xiàn)有的斷裂力學數(shù)值分析方法大多都是基于傳統(tǒng)有限元模型的.傳統(tǒng)有限元理論對裂紋等奇異性問題，存在計算精度低和效率下降的弊端.近年來，人們提出了一些新的有限元法，如有限差分法［1］、邊界元法［2－3］和無網(wǎng)格法［4－5］等，但由于缺少相應軟件的支持，使這些數(shù)值方法的發(fā)展受到了限制.

小波有限元法(Wavelet Finite Element Method，WFEM)是近幾年發(fā)展起來的一種新的數(shù)值分析方法.以小波函數(shù)或小波尺度函數(shù)作為插值函數(shù)，充分利用小波尺度函數(shù)的多分辨特性，可根據(jù)實際需要任意改變分析尺度，具有算法穩(wěn)定性好、計算精度和效率高的優(yōu)點.因此，在裂紋等奇異性問題數(shù)值計算方面具有誘人的優(yōu)越性.目前，小波有限元法研究對象基本上是確定性問題.然而，在實際工程結構中，由于各種因素的影響，使得結構的物理參數(shù)、幾何參數(shù)及載荷等具有不確定性，從而導致結構的響應也具有不確定性，這種不確定性需要用模糊理論來進行研究［6－7］.本文針對工程實際中幾何參數(shù)及載荷等不確定性裂紋結構的有限元分析問題，提出了基于BSWI模糊有限元分析法，建立了BSWI有限元模糊斷裂分析模型，采用區(qū)間數(shù)分解方法求解模糊平衡方程，基于輸入材料和載荷模糊數(shù)的隸屬函數(shù)，利用虛擬裂紋閉合法，計算并分析了含裂紋結構響應量的可能性分布.數(shù)值算例證明了該方法切實可行，并具有較高的計算精度.

1 BSWI單元構造

采用BSWI尺度函數(shù)的張量積插值構造單元，位移函數(shù)表示為

式中:Φ1，Φ2分別為m階j尺度空間的n(n=2j+m －1)個一維BSWI尺度函數(shù)［8］;p、q分別為待求小波系數(shù)列向量;ε，η分別為單元局部坐標.局部坐標與整體坐標關系為

式中:x1，y1分別為整體坐標下單元起始坐標;Lex為單元長度;Ley為寬度.

單元節(jié)點位移列陣為

將各節(jié)點坐標和位移代入式(1)和式(2)，即得:

由彈性力學知，平面應力問題的勢能泛函為

式中:Ωe為單元求解域;t為單元厚度;F={fx，fy}T為體力向量;p={px，py}T為面力向量;R={u，v}T為位移場向量;D為彈性矩陣;ε為應變矩陣，其中

式中:E，μ分別為材料的彈性模量和泊松比.

將式(3)、式(5)和式(6)聯(lián)立，由變分原理，令δΠp=0，可得單元剛度矩陣為

2 BSWI斷裂單元構造

BSWI啞節(jié)點斷裂單元(構造如圖1所示)，是用來從小波有限元分析結果中提取相關信息，進而計算裂紋尖端后面的張開位移和裂紋尖端前面的虛擬裂紋擴展量［9－10］.

沿裂紋方向，啞節(jié)點斷裂單元從兩側的小波單元中共提取5個節(jié)點:節(jié)點1和節(jié)點2對應于裂紋尖端，并用特殊剛度的彈簧連接;節(jié)點3和節(jié)點4在裂紋尖端的后面;節(jié)點5在裂紋尖端的前面.在實際應用時，節(jié)點1和節(jié)點2，節(jié)點3和節(jié)點4的坐標分別是重合的.

圖1 啞節(jié)點斷裂單元及節(jié)點位移矢量

設Kx和Ky分別為x和y方向的彈簧剛度，由圖1可知，裂紋尖端垂直節(jié)點力為

裂紋尖端后面的張開垂直位移為

裂紋前面的虛擬擴展為

3 BSWI模糊有限元分析

式中:Em為的主值，其值與不考慮彈性模量模糊性時的E值相同;EL，ER分別為的左、右展形.

根據(jù)式(7)，相應地，小波模糊單元剛度矩陣為

將式(14)代入式(12)得:

將式(15)代入式(13)得:

式中［Km］e與式(7)普通 BSWI單元剛度矩陣相同.

式中:Fm為的主值，F(xiàn)L，F(xiàn)R分別為的左、右展形.

假設彈簧連接裂尖節(jié)點i和i+1，組集啞節(jié)點斷裂單元和各BSWI模糊單元的剛度矩陣可得含裂紋結構總模糊剛度矩陣為

L－R型模糊結構平衡方程為

對模糊結構小波有限元平衡方程作λ水平截集，可得區(qū)間方程:

根據(jù)區(qū)間數(shù)運算法則和區(qū)間數(shù)分解形式的唯一性，λ水平截集下結構區(qū)間方程的解［11］為

4 BSWI虛擬裂紋閉合法

虛擬裂紋閉合法是用于二維斷裂參數(shù)問題計算的一步分析法，具有精度高、方法簡單的優(yōu)點.本文針對I型裂紋斷裂問題進行討論，但思路和方法適用于其他斷裂模式.如圖2所示，將裂紋從α擴展到α+Δα所需要的功與將裂紋從α+Δα閉合到α所需要的功是相等的.

圖2 虛擬裂紋閉合法

在虛擬裂紋線上裂紋尖端節(jié)點力在節(jié)點位移上做功為

假設虛擬裂紋尖端后面的張開位移和初始裂紋尖端后面的張開位移近似相等，式(22)改寫為

數(shù)值實驗表明虛擬裂紋閉合法對有限元網(wǎng)格尺寸并不敏感，因此，I型裂紋能量釋放率G1可近似表達為

式中B為裂紋體的厚度.

對于各向同性均勻線彈性材料而言，應力強度因子KI和能量釋放率GI的關系為

當λ取遍［0，1］中的一切值后，根據(jù)模糊分解定理可得結構模糊應力強度因子為

若只取有限個λ值進行計算，則:

5 算例分析

如圖3所示，共線雙邊裂紋單向拉伸平板.其中:L=2 m，W=2 m，裂紋長度a，板厚t=0.01 m，彈性模量 E=(210，0.04，0.04)LRGPa，泊松比μ =0.3，σ =(0.2，0.01，0.01)LRMPa，計算時，彈性模量、載荷均為L－R型模糊數(shù)，其隸屬函數(shù)取為線性函數(shù)為

假設為平面應力問題，試用基于BSWI小波有限元的模糊虛擬裂紋閉合法計算其應力強度因子.

圖3 裂紋板單向拉伸

為更好地利用BSWI有限元虛擬裂紋閉合法，考慮結構幾何和載荷的對稱性，取該結構1/2作為計算模型，在其底部施加對稱邊界條件，采用2個二階三尺度區(qū)間B樣條小波(BSWI23)和1個啞節(jié)點斷裂單元求解.

表1給出了當隸屬度為1.0時，采用2個BSWI23單元的模糊虛擬裂紋閉合法計算不同裂紋長度下應力強度因子與解析解的相對誤差.從表1中可看出最大相對誤差值為1.583%，說明本文提出的計算方法切實可行的，并可用較少單元獲得較高計算精度.

表1 隸屬度為1.0時模糊算法與解析解的相對誤差

表2給出了不同裂紋長度α下的應力強度因子KI隸屬函數(shù)值.由表2看出不同隸屬度下應力強度因子隨裂紋長度變化的關系.

圖4，5分別為裂紋面上距裂尖r=0.125，0.250 m處垂直位移v的隸屬函數(shù)圖.由于選取左、右基準函數(shù)L(x)，R(x)相同，垂直位移v分布均為中心對稱分布;由圖4，5可見，對于不同r值，當隸屬度較低時，垂直位移v的變化區(qū)間較大，當隸屬度較高時，垂直位移v變化區(qū)間較小，當隸屬度為1.0時，左右端點值均相等，分別為v0.125=4.569 6×10－5m和v0.250=6.446 6 ×10－5m，結果與確定性小波有限元結果一致.由此看出，利用本文建立的區(qū)間B樣條小波有限元模糊虛擬裂紋閉合法對含裂紋結構進行斷裂分析比確定性小波有限元分析更能全面反映結構響應的真實變化情況.

表2 不同裂紋長度下應力強度因子隸屬函數(shù)值 105Pa·m1/2

圖4 裂紋面上距裂尖r=0.125 m垂直位移v隸屬函數(shù)

圖5 裂紋面上距裂尖r=0.250 m垂直位移v隸屬函數(shù)

6 結 論

1)將本文方法的計算結果與解析解進行了比較，結果表明:當隸屬度為1.0時，不同裂紋長度下，應力強度因子隸屬函數(shù)值最大相對誤差為1.583%.驗證了該方法在解決含有不確定性參數(shù)的裂紋結構有限元分析方面的可靠性.

2)利用本方法計算了不同裂紋長度和隸屬度下的應力強度因子KI隸屬函數(shù)值.結果表明:對于不同r值，當隸屬度較低時，垂直位移v的變化區(qū)間較大，當隸屬度較高時，垂直位移v變化區(qū)間較小，當隸屬度為1.0時，左右端點值均相等.由此看出，本文建立的計算方法對含裂紋結構進行斷裂分析比確定性有限元分析法更能全面反映結構響應的真實變化情況.

3)小波有限元與傳統(tǒng)有限元相比，在計算裂紋等奇異性問題方面可用較少單元和自由度數(shù)，獲得較高計算精度和效率，因此，在數(shù)值計算方法上為求解具有材料和載荷不確定性的裂尖奇異性問題提供了一種可行的途徑.

［1］ZHANG Yang.A finite difference method for fractional partial differential equation［J］.Applied Mathematics and Computation，2009，215(2):524－529.

［2］程長征，牛忠榮，葉建喬.邊界元法計算淺表面裂紋應力強度因子［J］.合肥工業(yè)大學學報(自然科學版)，2009，32(4):503－507.

［3］SUN Yuzhou，WANG Jinyan，XU Junfeng.Boundary element method for the plane elastic problem with cracks［J］.Journal of Zhong yuan Institute of Technology，2004，15(5):46－54.

［4］張雄，劉巖，馬上.無網(wǎng)格法的理論及應用［J］.力學進展，2009，39(1):1－36.

［5］RAO B N，RAHMAN S.Probabilistic fracture mechanics by Galerkin meshless methods-Part I:Rates of stress intensity factors［J］.Computational Mechanics，2002，28(5):351－364.

［6］蘇靜波，卲國建.基于區(qū)間分析的工程結構不確定性研究現(xiàn)狀與展望［J］.力學進展，2005，35(3):338－344.

［7］RASHID M K，AL-ARAIMI S A.Fuzzy algorithm and structural stiffness in error attenuation of intelligent toolpost［J］.Journal of Intelligent Manufacturing，2005(16):277－286.

［8］何正嘉，陳雪峰，李兵，等.小波有限元理論及其工程應用［M］.北京:科學出版社，2006:40－42.

［9］XIE D.Damage progression in tailored laminated panels with a cutout and delamination growth in sandwich panels with tailored face sheets［D］.Clemson:Clemson U-niversity，2002.

［10］XIE D，Jr BIGGERS S B.Progressive crack growth analysis using interface element based on the virtual crack closure technique［J］.Finite Elements in Analysis and Design，2006，42(11):977－984.

［11］肖盛燮，王平義，呂恩琳.模糊數(shù)學在土木與水利工程中的應用［M］.北京:人民交通出版社，2004:349－368.

Fuzzy finite element analysis of B-spline wavelet on the interval of structure with cracks

PENG Hui-fen1，2，MENG Guang-wei1，ZHOU Li-ming1，LI Feng1

(1.College of Mechanical Science and Engineering，Jilin University，130022 Changchun，China，phfdaqing@163.com;2.College of Mechanical Science and Engineering，Northeast Petroleum University，163318 Daqing，China)

To overcome the difficults caused by uncertainties of material and load in numerical analysis of structure with cracks and to improve the accuracy and efficiency of numerical calculation，a fuzzy finite element analysis method of B-spline wavelet on the interval(BSWI)is put forward by combing the fuzzy theories with wavelet finite element.Fracture elements of dummy nodes are embedded into the wavelet finite element model with cracks，fuzzy equilibrium equations on the basis of finite element of BSWI are established，and the fuzzy equilibrium equations are solved by using λ level set and decomposition theorem.Subsequently，the membership function values of stress intensity factor with different crack lengths are calculated by using virtual crack closure technique，and the calculated results are compared with analytical solution.The analysis results show that the proposed method can accurately reflect changes in structural response with fewer elements and can be a new way for engineering fracture analysis of complex structures with uncertainties.

BSWI;crack;fuzzy numbers;stress intensive factor;finite element method

O242;TB115

A

0367－6234(2011)09－0134－05

2010－12－20.

高等學校博士學科專項科研基金資助項目(20060183063);吉林省科學技術廳基金資助項目(20090540);吉林大學“985工程”資助項目.

彭惠芬(1969—)，女，博士研究生，副教授;

孟廣偉(1959—)，男，教授，博士生導師.

(編輯 張 紅)