余能原理基面力元法在重力壩應(yīng)力分析中的應(yīng)用

張昭鋒，彭一江，郭 慶，單巖巖

（北京工業(yè)大學建筑工程學院，北京 100124）

近年來，隨著計算機技術(shù)的高速發(fā)展和廣泛應(yīng)用，有限元法有了很大的進展，已經(jīng)成為處理力學、物理、工程等問題的有效方法之一［1－4］。但是現(xiàn)有的有限元法采用的都是傳統(tǒng)二階應(yīng)力張量的描述體系，這種理論框架讓建立和表述數(shù)學模型較為繁瑣，難以尋求更大的突破。

對此，GAO 提出了一個遠較傳統(tǒng)的應(yīng)力張量簡單的新概念——“基面力”來描述應(yīng)力狀態(tài)，并較系統(tǒng)地描述了基面力的理論體系［5－6］。彭一江等人在基于基面力概念的余能原理有限元法方面進行了研究，推導出了適用于大位移、大轉(zhuǎn)動問題的有限元列式，具有邊中節(jié)點的四邊形基面力元的幾何非線性有限元列式以及任意多面體基面力元的幾何非線性余能原理有限元列式，采用Matlab語言編制相應(yīng)的余能原理基面力元分析軟件，研究了該方法的計算性能，在適應(yīng)任意網(wǎng)格、提高計算精度和減少網(wǎng)格依賴方面取得了一些成果，并將該方法簡稱為“基面力元法”（base force element method，BFEM）［7－14］。目前，基于余能原理的基面力元法已經(jīng)得到了一些國際同行的關(guān)注［15］。

本文在介紹描述平面應(yīng)力狀態(tài)的基線力概念和具有邊中節(jié)點的平面4 節(jié)點基面力元模型的基礎(chǔ)上，探索了相應(yīng)的基面力元法在重力壩應(yīng)力分析中的應(yīng)用。

1 二維問題的余能原理基面力元法

1.1 二維受力狀態(tài)的基面力表征

考慮二維彈性體區(qū)域，P，Q分別表示一物質(zhì)點變形前后的位置矢量，xα（α＝1，2）表示該點的Lagrange坐標，則變形前后的基矢為

在向量dx1Q1，dx2Q2上作一個平行四邊形微元，以描述Q點附近的應(yīng)力狀態(tài)，將基線力［9］定義為

式中：約定dx3＝dx1。

1.2 單元應(yīng)力的表征

如圖1所示，一個具有邊中節(jié)點的基面力元。

當單元足夠小時，可以假設(shè)應(yīng)力均勻地分布在每一邊上，則單元的平均應(yīng)力的表達式為

式中：A為單元的面積；PI為邊中節(jié)點的徑矢。

單元應(yīng)力的展開表達式可寫為

式中：e1，e2表示單位基矢量。

1.3 單元柔度矩陣顯式表達式

單元柔度矩陣CIJ的展開表達式為［9］

單元柔度矩陣CIJ的矩陣形式為

式（6）為平面應(yīng)力狀態(tài)下的表達式。對于平面應(yīng)變問題，只需將式中的E換為E／（1－v2），v換為v／（1－v）。

1.4 余能原理基面力元法控制方程

余能原理基面力元法的控制方程是如下泛函的約束極值問題：

式中：WeC為單元的余能為位移邊界Γu上的給定位移；T（Γu）為Γu上的面力。

利用Lagrange乘子法，放松平衡條件約束，可得單元的修正泛函為

式中：λ，λ3為Lagrange乘子。

系統(tǒng)的修正泛函為

泛函的駐值條件為

由式（10）得到下列線性方程組

求解該線性方程組，得到各個單元的面力。

1.5 節(jié)點位移的表征

節(jié)點位移可由各單元的支配方程求得

式中：ε為置換張量，在直角坐標系下的表達式為

節(jié)點位移的展開表達式為

2 應(yīng)用算例

一混凝土重力壩，壩高65m，底寬49m，混凝土彈性模量E＝15GPa，泊松比取v＝0.2，水容重1t／m3，混凝土容重2.45t／m3。汛期水位為60m，考慮水壓力和混凝土重力共同作用下的應(yīng)力分析。計算時按平面應(yīng)變問題考慮，采用20×40的4邊中節(jié)點四邊形單元，單元剖分網(wǎng)格如圖2所示。

圖2 重力壩及基面力元網(wǎng)格Fig.2 A gravity dam and the meshes of BFEM

計算所得擋水面上各單元的應(yīng)力及節(jié)點位移，取大壩擋水面上單元的應(yīng)力、位移繪于圖3。

圖3 大壩擋水面的應(yīng)力和位移Fig.3 Stresses and displacements at the upstream face of the dam

計算結(jié)果表明：二維線彈性的余能原理基面力元法能夠用于重力壩應(yīng)力分析，與ABAQUS有限元分析結(jié)果吻合較好，計算精度較高。

3 結(jié) 論

1）算例結(jié)果表明，基于余能原理的基面力元法和大型程序ABAQUS計算結(jié)果吻合，可以應(yīng)用于重力壩應(yīng)力分析問題。

2）基于余能原理的基面力元法的特點是：以基線力為基本未知量，來表征物體的受力狀態(tài)，以其對偶量位移梯度表征物體的變形狀態(tài)；該方法不用引入單元插值函數(shù)；單元應(yīng)力及節(jié)點位移列式均為顯式，無需進行數(shù)值積分。

3）研究工作表明，基于余能原理的基面力元法對單元網(wǎng)格的長寬比變化不敏感，具有較好的計算性能，具有較廣闊的應(yīng)用前景。

／References：

［1］ 王勖成.有限單元法［M］.北京：清華大學出版社，2003.103－121.WANG Maocheng.Finite Element Method［M］.Beijing：Tsinghua University Press，2003.103－121.

［2］ 王勖成，邵 敏.有限單元法基本原理和數(shù)值方法［M］.第2版.北京：清華大學出版社，1997.88－95.WANG Maocheng，SHAO Min.Basic Principles of Finite Element Method and Numerical Method［M］.2nd ed.Beijing：Tsinghua University Press，1997.88－95.

［3］ 王 軍，杜 壯，王文娜，等.基于ANSYS的橋式起重機主梁三維有限元分析［J］.河北工業(yè)科技，2010，27（6）：361－363.WANG Jun，DU Zhuang，WANG Wenna，et al.3Dfinite element analysis of bridge crane girder based on ANSYS［J］.Hebei Journal of Industrial Science and Technology，2010，27（6）：361－363.

［4］ 孫 旸，劉金祥.6114柴油機氣缸蓋的有限元分析［J］.河北工業(yè)科技，2009，26（3）：147－150.SUN Yang，LIU Jinxiang.Method of finite element analysis for cylinder head of engine［J］.Hebei Journal of Industrial Science and Technology，2009，26（3）：147－150.

［5］ GAO Y C.Large deformation field near a crack tip in rubberlike material［J］.Theor Appl Fract Mech，1997，26：155－162.

［6］ GAO Y C.A new description of the stress state at a point with applications［J］.Archive Appl Mech，2003，73：171－183.

［7］ 彭一江，金 明.基于基面力概念的一種新型余能原理有限元方法［J］.應(yīng)用力學學報，2006，23（4）：649－652.PENG Yijiang，JIN Ming.A new finite element method of complementary energy principle based on the concept of the base forces［J］.Chinese Journal of Applied Mechanics，2006，23（4）：649－652.

［8］ 彭一江，金 明.基于基面力概念的余能原理任意網(wǎng)格有限元方法［J］.工程力學，2007，24（6）：41－45.PENG Yijiang，JIN Ming.Finite element method for arbitrary meshes based on complementary energy principle using base forces［J］.Engineering Mechanics，2007，24（6）：41－45.

［9］ LIU Y H，PENG Y J.Base force element method（BFEM）on omplementary energy principle for linear elasticity problem［J］.Science China：Physics，Mechanics ＆ Astronomy，2011，54（11）：2025－2032.

［10］ 彭一江，劉應(yīng)華.基面力概念在幾何非線性余能有限元中的應(yīng)用［J］.力學學報，2008，40（4）：496－501.PENG Yijiang，LIU Yinghua.The application of the concept of the base forces to geometrically nonlinear analysis［J］.Chinese Journal of Theoretical and Applied Mechanics，2008，40（4）：496－501.

［11］ 彭一江，劉應(yīng)華.一種基于基線力概念的平面幾何非線性余能原理有限元模型［J］.固體力學學報，2008，29（4）：365－372.PENG Yijiang，LIU Yinghua.A plane geometric nonlinear finite element model of omplementary energy principle based on the concept of the baseline forces［J］.Chinese Journal of Solid Mechanics，2008，29（4）：365－372.

［12］ 彭一江，劉應(yīng)華.基于余能原理的有限變形問題有限元列式［J］.計算力學學報，2009，26（4）：460－465.PENG Yijiang，LIU Yinghua.Finite element formulas of finite deflection programs based on complementary energy principle［J］.Chinese Journal of Computational Mechanics，2009，26（4）：460－465.

［13］ PENG Y J，LIU Y H.Base force element method（BFEM）of complementary energy principle for large rotation problems［J］.Acta Mechanica Sinica，2009，25（4）：507－515.

［14］ PENG Y J，DONG Z L，PENG B，et al.The application of 2Dbase force element method（BFEM）to geometrically nonlinear analysis［J］.International Journal of Non－Linear Mechanics，2012，47（3）：153－161.

［15］ SANTOS H A F A，PIMENTA P M，ALMEIDA J P M.A hybrid－mixed finite element formulation for the geometrically exact analysis of three－dimensional framed structures［J］.Computational Mechanics，2011，48（5）：591－613.