多孔介质弹塑性变形的积分求解与数值实现.pdfVIP

多孔介质弹塑性变形的积分求解与数值实现 黄林冲1，2 徐志胜3 1.中山大学工学院，广东广州410075 2.西北大学土木与环境工程系，伊利诺伊埃文斯顿60208 3.中南大学土木建筑学院，湖南长沙410075 摘要：针对多孔渗水介质，基于Borja和Andrade建立的临界状态塑性模型，进行了弹塑性变形的积分求解，编 制了大型有限元计算程序，实现了有限元数值模拟。从非线性的弹性本构方程出发，对其本构关系的率形式进行了积 分，得到了应力求解的迭代方程；在非线性准则下的弹塑性变形中，引入平均有效应力、剪应力以及一致性系数，设计 了牛顿迭代算法流程，得到了多孔介质弹塑性变形增量形式的积分求解；最后编制了大型有限元计算程序以及积分求 解核心算法程序，再现了饱和砂土边坡液化的形成。结果表明该积分算法以及算法程序能有效地求解多孔介质弹塑性 变形的力学行为；在边坡局部剪切带即将形成的区域，孔隙水压力明显增大，而且沿着剪切带的外法线方向向外流动； 在远离剪切带区域，水压力值分布较为均匀，且流动方向趋于水平；根据液化发生准则，得到了边坡发生液化的滑动 面，这对于预测边坡的发生具有重要的现实意义。 关键词 变形；多孔介质；积分求解；有限元；数值模拟 Integrationalgorithmandnumericalimplementation ofelasto-plasticdeformationinporemedium HUANGLin-chong1,2XUZhi-sheng3 1.SchoolofEngineering,SunYat-senUniversity,Guangzhou510275,China 2.DepartmentofCivilandEnvironmentalEngineering,NorthwesternUniversity, Evanston,IL60208,USA 3.SchoolofCivilandArchitecturalEngineering,CentralSouthUniversity, Changsha410075,China Abstract:BasedonthecriticalstateplasticitymodelbyBorjaandAndrade,anintegrationalgorithmwasutilizedtocalcu latetheelasto-plasticdeformationoftheporemedium,andafiniteelementprogramwascodedtoimplementthefiniteelement simulation.Theiterativefunctionsforstresscalculatingwereresolvedbyintegratingwiththerateformoftheelasticconstitutive relationships.Further,intheelasto-plasticdeformationunderthenonlinearlaw,themeaneffectivestress,theshearstress, andthecoincidencecoefficientwereintroduced,andanewNewtoniterativealgorithmwasdesignedtocalculatetheincremental formoftheelasto-plasticdeformation.Atlast,afiniteelementprogramandanintegrationalgorithmwerecodedbyFORTRAN tosimulatetheonsetoftheliquefactionofasandslope.Itwasshownthatthisintegrationalgorithmandalgorithmcodecancalcu latetheelasto-plasticdeformationoftheporemediumeffectively.Theporepressureaffectedth