冻胀问题建模.ppt

冻胀问题建模;土壤冻胀问题;冻胀过程中的物理原理：;为简单计，假定水平方向上土壤质地是各向同性的，故只考虑

在水平方向上有变化的一维模型，其中x轴的正向垂直向上。;在任何层区域中，总体积为Ｖ的块区域中有三种物质成份：;由质量守恒定律，可导出：;由热量守恒定律，可导出：;这是因为在⊿t时段内，单位体积块区域(V＝１)中第i号物质;令⊿t→0,⊿x→0,即得偏微分方程：;注意到;;所以在任何区域中，温度函数T(x,t)均服从以下状态方程：;而在冰、水融混〔mush〕区域R1中，冰与土壤颗粒无迁移，

故有;为了在区域R1内进一步简化温度函数所服从的状态方程，假定

此区域内水的块密度与水的温度有关，即可表达为;???冰、水融混层区域R1内温度函数所服从的状态方程可简化为：;在完全冻结的冰土层区域R2中，显然有;综合之，在三个区域里描述温度T(x,t)所服从的状态方程模型

如下：;在任何两个区域的交界面x=ni(t),(i=0,1)上，温度函数

T(x,t)均是连续的，故有：;根据交界面x=ni(t),(i=0,1)两侧热量传递时的守恒关系，

还有：;单位体积元中三种物质在物相变化过程中，需要的热量数为：;现在，在湿土层区域R0和冰、水融混层区域R1的交界面

n0(t)上,由于温度的连续性和冻结温度为零度，那么有：;在冰、水融混层区域R1和冻结层区域R2的交界面n1(t)上,;再根据冰、水在跨越交界面n1(t)时的质量守恒关系;因为在交界面n0(t)上,水流量是连续的，故有：;归纳之，在交界面n1(t)处，有：;再根据热传导过程中的Fourier定律：;为了简化问题，假定，注意到就有;最后，我们建立了土壤冻胀的数学模型如下：;土壤冻胀的数学模型的图示如下：