趙新娟,趙吉義

(1.河南工業(yè)大學(xué),鄭州450052;2.平原工業(yè)濾器有限公司,河南新鄉(xiāng)453000)

位勢問題的雜交有限元算法研究

趙新娟1,趙吉義2

(1.河南工業(yè)大學(xué),鄭州450052;2.平原工業(yè)濾器有限公司,河南新鄉(xiāng)453000)

利用各向異性位勢問題的基本解,推導(dǎo)出了一種用于求解二維各向異性位勢問題的雜交有限元方法.在該方法中,用基本解的線性組合近似表示的單元域內(nèi)的位勢場分布可以解析滿足控制方程,而采用普通的形函數(shù)插值表示的單元邊界位勢場被用于保證單元之間的協(xié)調(diào)性要求;結(jié)合新構(gòu)造的積分雜交泛函,可以得到只涉及單元邊界積分的計算公式.算例結(jié)果表明了該算法的精確性和有效性,且可以推廣到其他問題的求解.

雜交有限元法;基本解;各向異性;位勢問題

作為一種基于基本解(格林函數(shù))的邊界積分型算法,邊界元法(邊界積分方程方法)存在邊界點及近邊界點奇異和幾乎奇異積分的計算難點,其準確性很大程度上依賴于這些奇異積分的計算精度,且當源點為域內(nèi)近邊界點時,其物理參量的計算通常由于幾乎奇異積分導(dǎo)致失效,這種現(xiàn)象稱為邊界層效應(yīng)[1-2].除了邊界元法之外,基于基本解插值的雜交有限元法是另外一種新興的計算方法.該方法融合了有限元法[3]和邊界元法的計算優(yōu)勢,可以同時消除邊界元法中的奇異積分計算和邊界層效應(yīng),提高近邊界位勢梯度計算精度.該算法利用相互獨立的單元域內(nèi)和單元邊界插值方案,可以構(gòu)造出滿足需要的單元類型和保持單元邊界積分特性.為了保證控制方程的解析滿足,基本解的線性組合被用于構(gòu)造單元域內(nèi)位勢場,而普通有限元插值則被用來構(gòu)造單元邊界位勢場,以保證相鄰單元之間的協(xié)調(diào)性要求.文獻[4]導(dǎo)出了二維位勢問題的雜交有限元計算公式,并討論了不同的插值方案對計算結(jié)果的影響;文獻[5]研究了正交各向異性材料中彈性響應(yīng)的雜交有限元分析;文獻[6]研究了復(fù)雜生物結(jié)構(gòu)眼睛中生物熱的雜交有限元分析.除此之外,國內(nèi)外尚無這種算法研究的公開報道.本文將推導(dǎo)各向異性位勢問題的雜交有限元法實現(xiàn)過程,并采用此算法計算各向異性介質(zhì)內(nèi)的位勢場分布.

1 雜交有限元法分析過程

二維各向異性材料位勢問題的控制方程為[7]:

式中:u是待求的勢函數(shù);kij(i=1,2;j=1,2)是各向異性材料特性系數(shù).對該問題,其對應(yīng)的雜交泛函可以為:

式中:,表示法向勢流;是給定的邊界條件;u和u分別表示相互獨立的單元域內(nèi)插值場和邊界插值場.

在基于基本解的雜交有限元算法中,單元域內(nèi)位勢場被假設(shè)為基本解的線性組合,例如:

而單元邊界位勢場則被假設(shè)為普通的形函數(shù)插值:

式中,{ce}和{de}分別為插值系數(shù)和未知的節(jié)點位勢矢量.其基本解為:

(x1,x2)為場點坐標,(xs1,xs2)為源點坐標.

對式(2),用高斯散度定理可以消除域積分,從而得到如下形式的單元邊界積分方程:

把式(3)和式(4)代入式(6)得到

其中:

為了建立{ce}和{de}之間的轉(zhuǎn)換關(guān)系并得到剛度方程,給出泛函Πe的極值條件:

進而可以得到:

至此,我們已經(jīng)推導(dǎo)出了各向異性材料位勢問題的雜交有限元分析公式.從中我們發(fā)現(xiàn),獨立的插值方案和完全的邊界積分允許自行構(gòu)造單元,而不必局限于普通有限元法中的單元類型,這為構(gòu)造滿足需要的特殊單元提供了保障.下面考察各向異性正方形區(qū)域內(nèi)的位勢分布.

假設(shè)材料的各向異性參數(shù)為k11=1,k22=5,k12=2,正方形邊長為6 m,邊界條件如圖1所示.該問題的精確解為T(x1,x2)=3x21-x22+x1x2.采用49個8節(jié)點單元,方形區(qū)域內(nèi)的位勢分布如圖2所示.可以看到,采用本文所提出的雜交有限元算法得到的計算結(jié)果和ABAQUS結(jié)果吻合較好,說明該算法具有較好的計算精度.

2 結(jié) 語

本文導(dǎo)出了一種用于求解二維各向異性位勢問題的雜交有限元方法.該方法采用基本解的線性組合近似表示單元域內(nèi)的位勢場分布,單元邊界位勢場采用普通的形函數(shù)插值,結(jié)合新構(gòu)造的積分雜交泛函,可以得到只涉及單元邊界積分的計算公式.算例結(jié)果表明了該算法的精確性和有效性,且可以推廣到其他問題的求解.

[1]Brebbia C A,Telles J C F,W robel L C.Boundary Element Techniques[M].Berlin:Sp ringer-Verlag,1984.

[2]周煥林,牛忠榮,王秀喜,等.正交各向異性位勢問題邊界元法中幾乎奇異積分的解析算法[J].應(yīng)用力學(xué)學(xué)報,2005,22(2):193-198.

[3]Chandrupatla T R,Belegundu A D.Introduction to Finite Elements in Engineering(3rd ed)[M].New Jersey:Prentice Hall,2002.

[4]Wang H,Qin Q H.Hybrid FEM w ith Fundamental Solutions as Trial Functions for Heat Conduction Simulation[J].Acta Mechanica Solida Sinica,2009,22(5):487-498.

[5]Wang H,Qin Q H.Fundamental-solution-based Finite Element Model for Plane Orthotropic Elastic Bodies[J].European Journal of Mechanics-A/Solids,2010,29(5):801-809.

[6]Wang H,Qin Q H.FEApp roach w ith Green’s Function as Internal Trial Function fo r Simulating Bioheat Transfer in the Human Eye[J].A rchievesof Mechanics,2010,62(6):493-510.

[7]Carslaw H S,Jaeger J C.Conduction of Heat in Solids(2nd ed)[M].Oxfo rd:Oxfo rd University Press,1959.

Potential Problems in Anisotropic Solids Using Hybrid Finite Element Model

ZHAO Xin-juan1,ZHAO Ji-yi2
(1.Henan University of Technology,Zhengzhou 450052;2.Pingyuan Industry Filter Co.,L td,Xinxiang 453000,China)

The hybrid finite element formulation is developed in the paper for p lane anisotropic potential analysis.The linear com bination of fundamental solution is used to app roximate the element interio r field,w hile the app roximation of shape function is used to give element boundary field.The hybrid variational is constructed to link these two fields and p roduce the final system of equations involving boundary integrals only.Numerical results show that the p roposed app roach has good accuracy,and can be extended to solve other p roblem s.

hybrid finite element;fundamental solution;anisotropy;potential p roblem s

O343.1

A

10.3969/j.issn.1671-6906.2011.01.016

1671-6906(2011)01-0059-03

2011-01-15

趙新娟(1976-),女,河南新鄉(xiāng)人,講師.