不可压流体饱和多孔弹性梁的变分原理及有限元方法.pdfVIP

第 30卷 第 1期 固 体 力 学 学 报 V01．30NO．1 2009年 2月 CHINESEJOURNALOFSOLID MECHANICS Feb2009 不可压流体饱和多孔弹性梁的变分原理及有限元方法 杨 骁 王佩菁 (上海大学土木工程系，上海 ，200072)(上海大学 ，上海市应用数学和力学研究所 ，上海 ，200072) 摘 要 基于不可压饱和多孔弹性梁动力弯曲的数学模型，建立了以多孔弹性梁挠度和孔隙流体压力等效力 偶为宗量的Gurtin型变分原理，并给 出了特殊边界条件下解耦时的仅 以挠度为宗量 的变分原理．同时，作为动力响 应 的退化情形 ，讨论 了拟静态情形下的相应变分原理．根据所建立 的变分原理 ，导出了一个有限元离散公式．由于 Gurtin型变分原理是关于时间的卷积型的泛函，空间的有限元离散导致一个关于时间的对称微分～积分方程组， 此方程组可进一步转化为常微分方程组 ．利用隐式 Euler法 ，给出了时间区域的计算格式．作为一个数值例子 ，分析 了饱和