冻土损伤的本构模型和耦合问题的数值模拟概念.docVIP

冻土损伤本构模型和耦合问题的数值模拟 作者：朱志武①， ， ， 610031； ， 100081； ， 730000 收稿日期： 2010-01-06； ： 2010-02-23 ，。，，。 ，，，，， 。，。 ，， ， 。 ， ， ， ， 1 引言 冻土是一种对温度敏感和易变的特殊低温地质体， 23%， 21。5%。 ，，， 。， ， ， ，， 。 ，，，，，。-土的相互作用已经有了比较深刻的认识和研究。 Fremond[1]； Gary[2]； Konard[3]； Dennis[4]；[5]先后对冻土的水分迁移与热质迁移， ，。 ， (温度场与水分场)的耦合作用， ，， 。， 。Foriero [6]用有限元模拟了冻土边坡的蠕变；Su [7]对多年冻土区通风路基进行了数值分析研究；Dempsey[8]； Rai [9]导出了流体、热流和土体变形完全耦合的控制微分方程；Masters[10]， ；Zhang [11]对寒区碎石路基在冻结过程中边界条件的影响进行了深入研究。， 。， ， ，，， 。，， ，。 2 2。1 从细观力学的角度出发， ，，， ， 。，。 Es 和Ei， Gs 和Gi， s 和νi。 s 和i分别表示冻土中土和冰的各个分量。K 和等效剪切模量G 为 式中：；。： K =E 3(1 2ν )， G = E 2(1+ )。 E和等效泊松比ν 的