地基沉降大变形有限元分析的几何刚度效应.docx

《地基沉降大变形有限元分析的几何刚度效应.docx》由会员分享，可在线阅读，更多相关《地基沉降大变形有限元分析的几何刚度效应.docx（6页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、地基沉降大变形有限元分析的几何刚度效应 摘要：随着技术的发展，大变形有限元法在地基沉降分析中的应用也越来越广，几何刚度效应作为大变形有限元分析的基本概念。在本案，笔者简要阐释了大变形有限元法，并基于算例分析基础上，探析了地基沉降的几何刚度效应。 关键字：地基沉降 大变形有限元法 几何刚度效应 地基作为建筑工程的保障性工程之一，其施工质量将干脆影响到上部建筑物的稳定性及运用寿命等，但因众多建筑工程地基为土性较差的软土地基，且因此引发沉降大变形的平安事故时有发生。地基沉降大变形势必会对建筑物造成干脆性或间接性的影响，且修补工作难度系数相当大。在对地基沉降大变形分析方面，因连续介质大变形固结理论自身

2、存在肯定困难性，则较具可行性的分析方法就属以有限元法为主的数值方法。在本案，笔者将从各个角度对地基沉降大变形有限元分析的几何刚度效应进行阐释。 一、大变形有限元方程与几何刚度效应 应用较为普遍的大变形限制方程及大变形有限原方程均基于士力学符号推导所得，相对于其他大变形限制方程及大变形有限原方程，由士力学符号推导而来的方程式未简化位移-应变几何关系，这一特征在几何刚度效应中发挥着不容忽视的作用。地基沉降大变形增量FEM函数式可表达为： 上述方程式中，t时刻内总应力场等效节点力矢量；有限元系统刚度矩阵； 单元节点未知量矢量；t+t时刻内面力与体力载荷等效节点力矢量；其中，满意函数式：。 上述方程式

3、中，大变形应变转换矩阵；几何刚度矩阵；单元刚度矩阵；D本构矩阵； M应力矢量S相关重量建构的应力矩阵；G形函数梯度矩阵。 由以上分析结果可得，若选用困难性最低的线性次弹性本构关系所得大变形FEM增量方程组式亦为非线性，造成以上结果的缘由在于属节点位移增量隐函数。对于某些忽视内非线性项的做法，其造成的误差也不小，则可以确定的是将内非线性项忽视等效于将大位移造成的结构刚度改变忽视，且若将其应用于地基大变形沉降问题分析上更是不行取。相对于此，本案中提及的大变形限制方程及大变形有限原方程的精确性明显更高。 就力学意义而言，几何刚度矩阵于大变形FEM方程组内的作用体现为大变形条件下系统结构刚度受应力水平

4、的影响程度，即几何刚度效应。为了使计算过程更加的简化，笔者通过引入系数，并得到如下方程式： 若=1，则表示严格考虑几何刚度效应；若=0，则忽视几何刚度效应；若=-1，则几何刚度矩阵取相反数。若大变形系统刚度函数表达式由士力学符号推导所得，之前存在一个“-”；而固体力学大变形系统刚度矩阵函数表达式内，之前存在一个“+”。以上现象出现的缘由为应变度量及应力符号约定间存在差异性。 二、算例分析 荷载-变形曲线与几何刚度效应 下图为压力荷载作用下的荷载-沉降量关系图。由图可以得出一个规律，即压力荷载作用下地基沉降量s随着荷载p改变而改变。图内SS表示小变形法的结果，FS表示为大变形法的结果。若地基土体

5、变形小，SS与FS间差异较模糊，且值对大变形分析结果的影响程度也不明显，但不同值随着变形程度的加深而对大变形分析结果的影响不断上升。若将几何刚度矩阵=0忽视，则计算出的沉降量将比对几何刚度矩阵=1正确考虑的计算结果更小。针对某些学者提出将自修正迭代算法应用于大变形固结方程组内进行求解，及路基沉降量在1m2m范围内时，几何刚度矩阵受到的影响较小的论断，下图均一一对其进行了科学验证。下图计算结果显示：几何刚度矩阵应考虑全面，且几何刚度矩阵的影响程度会随着变形程度的加深而扩大；若将几何刚度矩阵忽视，则地基大变形系统刚度会随之增加，其造成的影响为实际计算结果与刚硬更加接近。 大变形计算本构关系与几何刚

6、度效应 下图为小变形法与大变形法计算本构关系图： 备注：小变形法应变度量为工程应变、应力为Cauchy应力；大变形法应变度量为Green应变E11、应力为Kirchhoff应力S11。 由上图可知，简洁的线性次弹性本构模型所对应的大变形应变与全量应变间的关系亦为非线性。通过对以上问题进行深化探讨可得，本构关系曲线形态受到值的影响程度较小，且值仅对p-s曲线最终走向及应变值和应力最终大小造成影响，缘由是变形量因条件不同而有所差异。 下表显示了以有限元分析结果为依据所得的小变形法原点各线压缩模量理论值及计算值。 由上表可得，压缩模量相对误差肯定值均0.02%，则本案中笔者选用的数值计算方法及计算结

7、果均存在可靠性。 不同率型本构关系与几何刚度效应 前文中笔者提及的大变形率型本构关系均属于Jaumann应力率，从某种程度上讲，其在实际应用中存在肯定局限性。随着技术的发展，Truesdell应力率应运而生，探讨结果表明，后者的客观性略强于强者。探讨结果表明，若将领域变形忽视，Jaumann应力率可转化为Truesdell应力率。在本案，笔者将基于前人就 Jaumann 应力率与Truesdell应力率本构关系的探讨基础上，探析Jaumann应力率及Truesdell应力率与几何刚度效应间的关系。 下图分别为Truesdell率型本构关系于不同值条件下的计算本构关系曲线图及大变形荷载-沉降变形

8、曲线图。 由图可得，忽视几何刚度效应所引起的沉降计算结果将随着变形程度的加深而持续增加，且其所引起的变形结果刚硬度更强；若几何刚度效应于=-1条件下，其引起的计算偏差将略显著于其他条件下的计算偏差，此时的沉降值也最小。上述结论均在Jaumann率型几何刚度效应分析结果中得到了印证。 因几何刚度效应对大变形条件下系统结构刚度的影响程度较为明显，则在值不同、荷载水平相同的条件下，结构刚度沉降变形计算值将存在肯定差异。 三、总结 就荷载-沉降曲线关系而言，若将几何刚度效应忽视或为妥当处理好几何刚度效应，其所引起的地基沉降计算结果将会略小于精确值，且计算结果刚硬度更强，其为实际工程中重要的不平安因素，

9、则应得到更多的重视。Truesdell率型大变形理论所预料的最终沉降量符合Gibson理论完整解析解。Jaumann率型大变形法地基沉降预料结果比小变形法预料结果略大。须要留意的是，若没有妥当处理好几何刚度效应，大变形分析结构性质将因此受到很大的影响。 参考文献： 1 丁洲祥.地基沉降大变形计算方法的关键问题探讨J.岩土力学，2022，30：1920-1926. 2 丁洲祥，朱合华，丁文其等.地基沉降大变形有限元分析的几何刚度效应J.岩土力学，2022，30：1275-1280. 3 黎莉，赵明华，刘晓明等.高填石路堤施工阶段地基沉降分析方法初探J.马路，2002，：64-66. 4 吴健.饱和软土困难非线性大变形固结特性及应用探讨D.浙江高校建筑工程学院，2022. 5 刘祚秋，周翠英.软粘土地基非线性一维大变形固结的有限差分法分析J.中山高校学报，2022，44：25-28，41. 第6页 共6页第 6 页 共 6 页第 6 页 共 6 页第 6 页 共 6 页第 6 页 共 6 页第 6 页 共 6 页第 6 页 共 6 页第 6 页 共 6 页第 6 页 共 6 页第 6 页 共 6 页第 6 页 共 6 页