下限原理有教材限单元法及其在土工问题中的应用.pdfVIP

　第 19 卷 　第 5 期 岩 　土 　工 　程 　学 　报 V ol. 19 　No . 5 　 　1997 年 　　 9 月 Chinese Journal of Geotechnical Engineering Sept . , 　1997 　 下限原理有限单元法及其 在土工问题中的应用 李国英 　　沈珠江 (南京水利科学研究院土工所, 210024) 1 　前　　言 50 年代, Drucker 和 Prager 把静力场和运动场结合起来并提出极值原理, 建立了土体极限 ( ) 分析理论, 为土体极限平衡课题的求解开辟了新的途径。W. F. Chen 1975 发表的专著《极限 分析和土的塑性力学》, 进一步阐明了极限分析理论在土工问题中的应用。从理论上讲, 由平 衡条件、屈服条件、流动法则以及相应的边界条件, 足以确定应力场、速度场和破坏荷载, 但由 于实际问题的复杂性, 要求全面满足静力方程、运动方程以及相应边界条件的解答几无可能。 为此, 人们致力于开辟另一种求解途径[ 1 , 2] , 即引入变分原理, 通过求极值的方法求问题的解, 这对上限原理的研究和应用起到了较大的促进作用。下限原理的求解有赖于一个好的应力场 的构造, 然而, 实际问题中应力场的构造十分困难, 故下限原理的应用一直难以推广。随着现 代计算技术的发展和有限单元法的推广应用, 近年来已有人探索运用有限单元法来求解极限 平衡课题的下限值[ 3 , 4] 。运用有限单元法求解的基本思路是, 对土工结构进行有限元离散, 并 构造相应的单元应力模式, 根据下限原理的理论, 建立满足平衡方程、屈服条件和应力边界条 件的数学模型, 引入数学规划的方法寻求问题的下限解。计算所得的最后结果不仅有问题的 下限解, 还有结构内部相应的应力场。 2 　基本方程 本文采用三角形单元对计算区域进行离散, 单元中应力模式采用如下线性形式 : σ ψ σe ( ) { } = [ ]{ } 1 T 其中 σ σ σ τ { } = [ 　 　 ] ; x y xy e i i i j j j k k k T σ σ σ τ σ σ τ σ σ τ { } = [ 　 　 　 　 　 　 　 　 ] x y xy x y xy x y xy ψ [ ] = [ N I 　N I 　N I ] i j k 这里 i , j , k 为三角形单元的三个