Mechanical Properties and Water Holding Characteristics of Initially Isotropic Soils and Transversely Isotropic Soils
-
摘要: 为了探讨不同初始状态下非饱和黄土的力学性质和水量变化规律之间的差异,提出了制备各向同性及横观各向同性试样的方法,用改进的非饱和土三轴仪对制备的各向同性、横观各向同性及常规重塑试样进行了72个三轴固结排水剪切试验. 试验结果表明:对不同削样角度的各向同性试样进行试验验证,发现其近似各向同性;3种不同初始状态下各向同性非饱和试样的破坏应力最大,横观各向同性非饱和试样次之;横观各向同性试样在剪切阶段的初始切线模量最大,其他两组试样的初始切线模量较小且近似相等;各向同性试样及横观各向同性试样在剪切的初始阶段均有轻微的剪胀现象,常规重塑试样均处于剪缩状态;不同初始状态下的各非饱和试样在剪切阶段吸力变化对排水的影响并不显著,其中横观各向同性试样的偏应力-含水率关系曲线的斜率绝对值最小,其他两组试验的斜率几乎相同,各向同性、横观各向同性及常规的重塑试样的力学特性及持水特性不尽相同. 本文的研究成果可为不同初始状态下各模型的建立及工程应用提供有益的参考.Abstract: In order to investigate the difference between the mechanical properties and water quantity variation of unsaturated loess under different initial conditions, a method for preparing isotropic and transversely isotropic specimens was suggested. An improved triaxial apparatus for unsaturated soils was used to carry out 72 consolidated drained triaxial tests on prepared isotropic, transversely isotropic and regular specimens. Isotropic samples with different cutting angles were studied experimentally and found to be approximately isotropic. For three different initial states, the failure stress of the isotropic unsaturated specimen is the largest, followed by that of the transversely isotropic unsaturated specimen; the initial tangent modulus of the transversely isotropic specimen is the largest at the shear stage, while those of the other two groups are smaller and approximately equal. Both the isotropic and transversely isotropic specimens have slight dilatancy at the initial stage of shear, while all ordinary remoulded specimens are in the state of shear contraction. Under different initial conditions, the change of the suction of unsaturated samples has little effect on the drainage, the absolute slope value of the deviator stress-water content curve of the transversely isotropic sample is the lowest, and those of the other two groups are almost the same. The mechanical properties and water holding characteristics of the isotropic, transversely isotropic and ordinary remoulded specimens are different. The research results obtained may provide a useful reference for the establishment of different models of soil under different initial conditions and for various engineering applications.
-
随着我国城市化进程不断加快,隧道数量和复杂性都在不断增加. 目前,我国已成为世界上隧道数量最多、条件最复杂、技术发展最快的国家[1-2]. 为解决交通便捷性需求,城市隧道结构逐渐丰富,各类分岔隧道数量也急剧增加.
火灾作为隧道中危害性最大的灾害,一旦发生,将造成严重的人员财产损失. 相较于普通地面建筑,隧道中火灾的烟气蔓延情况及控制手段都更为复杂,使得隧道发生火灾时人员的安全疏散及灭火救援难度高于普通地面建筑[3-4].
针对于单直隧道火灾,国内外进行了大量的研究. Alpert[5]首次提出了预测顶棚射流最高羽流温度的经验模型;Kurioka等[6]通过改变隧道的横截面以及火源的热释放速率进行实验,提出隧道顶棚下的最高温度模型;Li等[7]考虑连续火焰长度低于顶棚高度的情况下,通过理论分析及实验研究,探讨2种不同尺寸隧道中,不同热释放速率、纵向通风风速及隧道几何形状对顶棚下烟气最高温度的影响;王钟宽等[8]以美国Memorial隧道为原型,应用Froude准则建立1∶20缩尺隧道模型,探究不同坡度的隧道火灾自熄特性;姜学鹏等[9]提出侧部点式排烟模式隧道火灾临界风速的无量纲计算式.
随着地下空间的开发与应用,隧道结构逐渐丰富. 因此,针对于分岔隧道火灾的研究也快速展开. Yang等[10]采用理论分析和数值模拟相结合的方法,研究在具备机械排烟系统的倾斜分岔隧道中烟气流动的多重行为;Huang等[11]通过小尺寸分岔隧道火灾实验,研究隧道顶棚下气体最大超限温度,提出一个分为两部分的数学模型来预测分岔隧道的最大上限温度;Zhang等[12]研究T型分岔隧道结构火灾的顶棚最大温升,并提出在目前的模型结构规模和热释放率的研究范围内T型分岔隧道的最大顶烟温度随着无量纲支管宽度的增加而略有下降.
目前,空气幕也被广泛应用于隧道火灾排烟控制中,其工作原理是让高速喷出的空气流形成一道墙,以起到隔烟隔热的效果. 学者们研究了空气幕对火灾烟气蔓延的控制作用. Shu等[13]利用一系列典型的空气动力学性能曲线,建立一种在无风条件下计算空气幕效率系数的新模型,对分岔隧道火灾空气幕控烟有借鉴作用. Yu等[14]提出当动量比约为10∶1时,吹风角度垂直向下的空气幕封烟效果最佳,且其效果随缝宽和喷射角的增大而提升;当火源位于空气幕一侧时,空气幕的最佳喷射角为朝向火源倾斜30°. Li等[15]提出空气幕射流的均匀性受气幕单元的横截面影响,随着射流距离的增加,空气幕的射流速度先急剧下降,然后缓慢下降,直至为0.
隧道火灾中,单独使用纵向通风或空气幕对分岔隧道进行烟气的控制效果仍不尽如人意[16-17]. 而针对纵向通风与空气幕协同作用下的分岔隧道烟气控制及最高温度的研究还较少,故本文根据57组缩尺寸实验数据,对分岔隧道多向气流耦合下的烟气运动和最高温度变化进行研究.
1. 实验布置
如图1所示,在缩小比例(1∶10)分岔隧道模型中进行实验,隧道长10.0 m、宽1.0 m、高0.6 m,分岔隧道长5.0 m、宽1.0 m、高0.6 m,夹角45°. 以液化石油气为燃料,火源位于主隧道中间.
分岔隧道的岔道口处安装有空气幕. 该空气幕采用铁制栅格通过软性连接与风机相连,并通过底部支架固定于隧道内,从而保证气流的稳定性,且各个位置的风速分布均匀. 同时,可通过侧面螺丝及手柄调节空气幕的角度A在 −45°~45° (面向主隧道为正,分岔隧道口为负),空气幕喷射宽度(厚度T)为0~0.15 m. 通过变频器也可随意调节空气幕风速大小,使空气幕射流速度$ {v_{{\mathrm{air}}}} $在0~4.0 m/s. 空气幕的整体结构如图1(b)、(c)所示. 在分岔隧道的左侧安装变频风机,用于控制纵向通风. 为监测隧道中最高温度的分布,将127个K型热电偶沿着主隧道和分岔隧道的中心轴放置,热电偶设置细节如图2所示. 纪杰等[18]研究狭长空间内的烟气流动特性发现,在小尺寸隧道中进行火灾实验,雷诺数Re (式(1))和弗洛德数Fr (式(2))对流体运动具有决定性影响.
Re=vlμ, (1) Fr=glv2, (2) 式中:$ v $为纵向通风速度,m/s;$ l $为长度尺寸,m;$ \mu $为流体的运动黏度.
基于几何相似原理,热量的相似关系可表示为
Qm=Qf(LmLf)5/2, (3) 式中:$ {Q_{\mathrm{m}}} $为模型尺寸热释放速率,kW;$ {Q_{\mathrm{f}}} $为全尺寸热释放速率,kW;$ {L_{\mathrm{m}}} $为模型隧道的长度,m;$ {L_{\mathrm{f}}} $为全尺寸隧道的长度,m.
基于该相似关系,当全尺寸隧道中热释放速率 (HRR)为5、15、25 MW时,缩尺寸隧道模型中分别对应为15.9、47.9、77.7 kW.
本文总共进行了57组实验,变量包括热释放速率、纵向通风速度、空气幕角度、厚度和风速,拟探究纵向通风与空气幕之间的耦合效应. 详细实验工况见表1.
表 1 实验工况Table 1. Experimental conditions工况 HRR/kW 纵向通风/(m•s−1) 空气幕风速/(m•s−1) 角度/(°) 厚度/m 1~15 15.9,47.9,77.7 0.4,0.8,1.2,1.6,2.0 1.0 0 0.15 16~33 47.9 0.4,0.8,1.2 1.5,2.0,2.5,3.0,3.5,4.0 0 0.15 34~51 47.9 0.4,0.8,1.2 3.5 −45,−30,−15,15,30,45 0.15 52~57 47.9 0.4,0.8,1.2 3.5 30 0.05,0.10 2. 不同纵向通风下的隧道温度分布
2.1 主线隧道
当空气幕射流速度较小且恒定($ {v_{{\mathrm{air}}}} $=1.0 m/s)时,3种不同热释放速率下纵向通风大小对主隧道内温度分布影响如图3所示. 图中,横坐标x0为距离主隧道火源的距离. 由图可知,当纵向风速较小时($ v $=0.4 m/s),其对烟气回流层的阻隔作用较小,从而导致隧道上游出现温度波动的情况,而火源热释放速率越大,该现象越明显. 当纵向风速为0.8 m/s及以上时,对主线隧道内烟气的降温效果由纵向通风主导,能在隧道上游形成空气屏障,有效缩短烟气回流层长度. 烟气的蔓延过程为非定常的三维湍流运动[19],而空气幕的作用方向与烟气在顶棚形成的射流方向垂直,导致空气幕通风会对烟气的一维水平蔓延形成剪切作用,在此基础上再添加纵向通风,会对火灾烟气在主线隧道下游蔓延的湍流运动产生较大扰乱,且火源功率越小影响越大. 在3种热释放速率下,随着纵向风速的增加,不仅隧道内最高气温降低,而且烟气最高温度的位置也向隧道下游有一定偏移.
由图3(a)可知:当热释放速率为15.9 kW时,纵向风速从0.4 m/s增加到2.0 m/s,隧道顶棚最高温度从150.1 ℃降低到49.5 ℃;纵向风速达到2.0 m/s时,烟气温度已无明显峰值,火羽流在纵向通风的作用下向火源下游倾斜;当纵向通风增大时,火羽流的倾斜程度增加,对隧道顶棚的直接冲击作用减弱,因此,在隧道顶棚下方形成的高温区域逐渐减小,并且向火源下游偏移. 由图3(b)可知:当热释放速率为47.9 kW时,纵向通风从0.4 m/s增加到2.0 m/s,顶棚最高温度从381.4 ℃降低到83.6 ℃;当纵向风速为2.0 m/s时,隧道顶棚最高温度位置根据火源位置向下游偏移1.0 m. 由图3(c)可知:当热释放速率为77.7 kW时,纵向风速从0.4 m/s增加到2.0 m/s,最高温度从512.3 ℃降低到149.2 ℃,主线隧道下游沿程温度变化随纵向风速的增加而趋于稳定. 综上所述,当空气幕射流速度不变时,随着热释放速率的增大,主线隧道的纵向通风降温效果会达到阈值,因此,有必要探究纵向通风与空气幕协同作用时的最佳风速组合.
2.2 分岔隧道
当$ {v_{{\mathrm{air}}}} $=1.0 m/s时,不同热释放速率下纵向风速大小对分岔隧道温度分布影响如图4所示. 分岔隧道的最高温度随着热释放速率的增大而显著升高,因空气幕良好的隔热阻烟效果,温度在x0=0.8 m处呈断崖式下降. 当纵向风速不同时,分岔隧道中x0=−0.5 m处的温度差异较大. 纵向风速过小(如0.4、0.8 m/s)会导致高温烟气在火源附近滞留,并向分岔隧道蔓延. 然而,加大纵向风速虽能一定程度上避免烟气向分岔隧道蔓延,但会加剧主线隧道下游烟气的沉降效果.
综上,在城市隧道分岔段采用空气幕进行隔热阻烟具有一定的应用价值,其对火灾烟气的具体影响将在下文进行论述. 由于当纵向风速超过1.2 m/s时其降温效果不再显著变化,因此,后续实验选择在纵向风速为0.4、0.8、1.2 m/s条件下进行.
3. 空气幕对隧道最高温升的影响
3.1 空气幕射流速度
为探究不同纵向风速和空气幕射流速度条件下分岔隧道的最大温升,基于Chen等[20]的研究,建立式(4)~(7)所示关系.
ΔTmax=f(Q,Cp,g,v,ρ∞,T∞,H), (4) ΔTmaxT∞=γ(Q2/30Fr1/3)ε, (5) Q0=Qρ∞CPT∞g1/2H5/2, (6) ΔTmaxT∞∼Q0Fr, (7) 式中:$ \gamma 、\varepsilon $为常数,当$ \dfrac{{\mathop Q\nolimits_{\text{0}}^{{2 \mathord{\left/ {\vphantom {2 3}} \right. } 3}} }}{{F{r^{{1 \mathord{\left/ {\vphantom {1 3}} \right. } 3}}}}} < 1.35 $时,$ \gamma =1.77,\varepsilon = \text{1}\text{.2} $,当$ \dfrac{{\mathop Q\nolimits_{\text{0}}^{{2 \mathord{\left/ {\vphantom {2 3}} \right. } 3}} }}{{F{r^{{1 \mathord{\left/ {\vphantom {1 3}} \right. } 3}}}}} \geqslant 1.35 $时,$ \gamma =\text{2}\text{.54},\varepsilon =\text{0} $;$ \Delta {T_{\max }} $为隧道顶棚最高温升,K;$ Q $为热释放速率,kW;Cp为恒压比热容,kJ/(kg•K);$ {\rho _\infty } $为环境空气密度,kg/m3;$ {T_\infty } $为环境温度,K;H为燃烧器表面到隧道顶棚的距离,m;$ \mathop Q\nolimits_{\text{0}} $为无量纲热释放速率.
根据陈贻来[21]的研究,单个汽车燃烧时的热释放速率为3~5 MW,当隧道中单个汽车发生火灾并引燃两侧汽车时,火灾热释放速率最高可达15 MW. 而夏永旭等[22]研究表明,大巴车在隧道中平均热释放速率为8 MW. 因此,为更好地与实际相结合,探究隧道火灾极端情况下不同纵向风速与空气幕对隧道顶棚温度的影响,后续实验选取热释放速率为15 MW,即在缩尺寸隧道中热释放速率为47.9 kW时进行. 因此,在该情况下,$ {Q_0} $可以作为一个常数而忽略,故而式(7)可以写为
ΔTmaxT∞∼Fr∼v2gH. (8) 图5显示在不同空气幕射流速度下$ {{{v^2}} / {gH}} $与$ {{\Delta {T_{\max }}} / {{T_\infty }}} $之间的关系. 图中,x、y分别为拟合曲线横、纵坐标参量,余图同. 由此可以得出,在不同空气幕射流速度的影响下,分岔隧道的最大温升与纵向通风之间的关系如式(9)所示.
ΔTmaxT∞=m(v2gH)−n, (9) m=−0.2vair√gH+1.75, (10) n=−0.08vair√gH+0.18, (11) 式中:m、n为系数.
拟合的详细描述见图6.
将m和n代入式(9),便可以得到Q = 47.9 kW,空气幕角度为0° 时,纵向通风与空气幕协同作用下分岔隧道最大温升的理论模型,如式(12)所示.
ΔTmaxT∞=(−0.2vair√gH+1.75)(v2gH)0.08vair√gH−0.18. (12) 将不同空气幕射流速度下主隧道最大温升实验值与理论模型(式(12))中的预测值进行比较,如图7所示,其结果吻合度较高.
3.2 空气幕角度和厚度
根据图5可以看出,在3种纵向通风下,当空气幕射流速度达到3.5 m/s时,空气幕射流速度对主隧道最高温影响减弱,此时,空气幕和纵向通风产生的气流达到动态平衡. 因此,为探索空气幕的角度A和厚度T对隧道温度分布的影响,将空气幕的射流速度$ {v_{{\mathrm{air}}}} $设置为固定值3.5 m/s,以排除其影响.
不同空气幕角度下,$ {v^2}/(gH) $与$ {{\Delta {T_{\max }}} / {{T_\infty }}} $之间的关系如图8所示. 由此可以得出,在不同空气幕角度(−45° 除外)的影响下,分岔隧道最大温升与纵向通风的关系如图9所示,其拟合式如式(13)所示.
图 8 不同空气幕角度下$ {{{v^2}} \mathord{\left/ {\vphantom {{{v^2}} {(gH)}}} \right. } {(gH)}} $与$ {{\Delta {T_{\max }}} \mathord{\left/ {\vphantom {{\Delta {T_{\max }}} {{T_\infty }}}} \right. } {{T_\infty }}} $的关系Figure 8. Relationship between $ v^{2} / ({gH}) $and $ \Delta T_{\max } / T_{\infty} $at different angles of air curtainΔTmaxT∞=αv2gH+β, (13) α = −0.9sin2A+0.8sinA−0.3, (14) β = −0.61sinA+1.57, (15) 式中:α、β为拟合系数.
将系数$ \alpha $和$ \beta $代入式(13)中,得到热释放速率为47.9 kW,$ {v_{{\mathrm{air}}}} $=3.5 m/s时,纵向通风与空气幕协同作用下空气幕角度对分岔隧道最大温升影响的理论模型,如式(16)所示.
ΔTmaxT∞=(−0.9sin2A+0.8sinA−0.3)v2gH−0.61sinA+1.57. (16) 如图10所示,当火源功率47.9 kW,$ {v_{{\mathrm{air}}}} $=3.5 m/s 时,不同空气幕角度下,隧道最大温升拟合得到的理论公式较好与实验结果基本一致.
如图11所示,$ {v_{{\mathrm{air}}}} $=3.5 m/s,A=30° 条件下进行了3种不同空气幕厚度的实验,当T=0.10 m时分岔隧道降温效果最佳. 由图11可以得出,当vair=3.5 m/s,A=30° 时,在不同空气幕厚度的影响下,分叉隧道最大温升与纵向通风之间的关系如式(17)所示.
ΔTmaxT∞=av2gH+bv√gH+c, (17) 式中:a、b、c为拟合系数,拟合的详细描述见图12.
将系数a、b和c代入式(17)中,便可以得到热释放速率为47.9 kW时,纵向通风和空气幕协同作用下空气幕厚度对分岔隧道最大温升影响的理论模型,如式(18)所示.
ΔTmaxT∞=av2gH−bv√gH+24T2H2−13TH+2.7, (18) 式中:a = −26.88T2/H2−12T/H−0.44,b = −44.8T2/H2 +22T/H−1.3.
将不同空气幕喷射宽度下主隧道最大温升实验值与理论模型(式(18))中的预测值进行比较,如图13所示,具有较好的吻合性.
4. 结 论
1) 空气幕能有效协助纵向通风降低主线隧道温度,防止烟气进入分岔隧道.
2) 当空气幕射流速度较小时($ {v_{{\mathrm{air}}}} $=1.0 m/s),纵向风速的增加能有效防止烟气在分岔口积聚,提高空气幕对烟气的阻隔效率,降低隧道分岔点处的烟气温度.
3) 根据隧道火灾顶棚最大温升的无量纲经验相关公式,提出在固定热释放速率47.9 kW条件下,隧道最大温升与纵向通风和空气幕射流速度之间的关系,该模型的预测结果与实验结果吻合度较高.
4) 在固定热释放速率47.9 kW,空气幕射流速度3.5 m/s条件下,提出分岔隧道最大温升与纵向风速和空气幕角度之间的关系,其预测结果与实验结果吻合度较高.
5) 在固定热释放速率47.9 kW,空气幕射流速度3.5 m/s及空气幕角度30° 条件下,提出分岔隧道最大温升与纵向风速和空气幕厚度之间的关系,其预测结果与实验结果吻合度较高.
-
表 1 土样的基本物理指标
Table 1. Physical parameters of soil samples
相对密度ds 塑限wp/% 液限wL/% 最大干密度
${\rho _{{\rm d}\max }}$/(g•cm−1)最优含水率wop/% 2.71 17.3 31.1 1.91 12.5 表 2 试验研究方案
Table 2. Experimental research programs
试验组别 试验简称 控制竖向压力/kPa 削样角度/(°) 控制净围压/kPa 控制基质吸力/kPa 1 各向同性试验 0 100、200、300 50、100、200 45 60 90 2 K0CD试验 100 100、200、300 0、50、100、200 200 3 CD试验 100、200、300 0、50、100、200 表 3 土样的强度参数
Table 3. Strength parameters
吸力/kPa 净围压p/kPa qf/kPa pf/kPa φ/(°) c/kPa 各向
同性K0CD
试验CD
试验各向
同性K0CD
试验CD
试验各向
同性K0CD
试验CD
试验各向
同性K0CD
试验CD
试验0 100 367.7 201.3 226.6 167.1 31.39 26.51 41.63 16.11 200 596.2 400.9 398.7 333.7 300 831.0 521.5 593.7 473.8 50 100 341.5 312.1 276.0 213.8 204.0 192.0 26.61 28.93 27.13 58.08 37.25 30.98 200 524.1 503.8 424.5 374.7 367.9 341.5 300 665.4 686.9 610.8 521.8 529.0 503.6 100 100 399.0 333.0 325.5 229.7 211.0 208.3 28.44 30.04 28.21 67.07 42.07 47.90 200 564.3 557.7 525.1 388.1 385.9 375.0 300 768.4 733.0 683.8 556.1 544.3 527.9 200 100 447.4 433.3 409.3 249.1 244.4 236.4 31.35 31.61 29.2 70.56 62.61 67.50 200 704.9 666.0 605.4 435.0 422.0 410.7 300 892 873.8 790.3 602.7 591.3 563.4 表 4 剪切过程中土样的q-w关系曲线的斜率
Table 4. Slope of the q-w relation curve of soil samples in the shear process
试样 吸力/kPa 净围压/kPa $\alpha$/(× 10–5) 各吸力下$\alpha$的平均值/(× 10–5) 平均值/(× 10–5) 各向同性 50 100 –3.27 –2.68 –2.62 200 –2.63 300 –2.15 100 100 –3.09 –2.61 200 –2.48 300 –2.26 200 100 –3.30 –2.56 200 –2.39 300 –2.00 横观各向同性 0 100 –4.98 –4.43 –1.55 200 –4.06 300 –3.72 50 100 –2.16 –1.55 200 –1.33 300 –1.16 100 100 –1.44 –1.52 200 –1.84 300 –1.27 200 100 –1.63 –1.59 200 –1.69 300 –1.45 重塑 0 100 –6.82 –5.46 –2.59 200 –5.05 300 –4.50 50 100 –3.02 –2.47 200 –2.49 300 –1.89 100 100 –2.95 –2.48 200 –2.49 300 –1.99 200 100 –3.20 –2.83 200 –2.78 300 –2.52 -
曾国熙,龚晓南,盛进源. 正常固结黏土K0固结剪切试验研究[J]. 浙江大学学报,1987,21(2): 1-9.ZENG Guoxi, GONG Xiaonan, SHENG Jinyuan. Research on normally consolidated clay by K0 consolidated shear test[J]. Journal of Zhejiang University, 1987, 21(2): 1-9. 黄茂松,宋晓宇,秦会来. K0固结粘土基坑抗隆起稳定性上限分析[J]. 岩土工程学报,2008,30(2): 250-255. doi: 10.3321/j.issn:1000-4548.2008.02.016HUANG Maosong, SONG Xiaoyu, QIN Huilai. Basal stability of braced excavations in K0-consolidated soft clay by upper bound method[J]. Chinese Journal of Geotechnical Engineering, 2008, 30(2): 250-255. doi: 10.3321/j.issn:1000-4548.2008.02.016 李校兵,郭林,蔡袁强,等. K0固结饱和软黏土的三轴应力路径试验研究[J]. 中南大学学报(自然科学版),2015,46(5): 1820-1825.LI Xiaobing, GUO Lin, CAI Yuanqiang, et al. Stress path triaxial tests on K0-consolidated saturated soft clay[J]. Journal of Central South University (Science and Technology), 2015, 46(5): 1820-1825. 王立忠,叶盛华,沈恺伦,等. K0固结软土不排水抗剪强度[J]. 岩土工程学报,2006,28(8): 970-977. doi: 10.3321/j.issn:1000-4548.2006.08.008WANG Lizhong, YE Shenghua, SHEN Kailun, et al. Undrained shear strength of K0 consolidated soft clays[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(8): 970-977. doi: 10.3321/j.issn:1000-4548.2006.08.008 褚福永,朱俊高,王平,等. K0固结条件下粗粒土变形及强度特性研究[J]. 岩土力学,2012,33(6): 1625-1630. doi: 10.3969/j.issn.1000-7598.2012.06.005CHU Fuyong, ZHU Jungao, WANG Ping, et al. Study of deformation and strength characteristics of coarse grained soil under K0-consolidation condition[J]. Rock and Soil Mechanics, 2012, 33(6): 1625-1630. doi: 10.3969/j.issn.1000-7598.2012.06.005 吕玺琳,钱建固,黄茂松. 不排水加载条件下K0固结饱和砂土失稳预测[J]. 岩土工程学报,2015,37(6): 1010-1015.LÜ Xilin, QIAN Jiangu, HUANG Maosong. Prediction of instability of K0-consolidated saturated sands under undrained loading conditions[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(6): 1010-1015. 刘清秉,吴云刚,项伟,等. K0及三轴应力状态下压实膨胀土膨胀模型研究[J]. 岩土力学,2016,37(10): 2795-2802.LIU Qingbing, WU Yungang, XIANG Wei, et al. Swelling model study of expansive soil at K0 and triaxial stress state[J]. Rock and Soil Mechanics, 2016, 37(10): 2795-2802. 程海涛, 刘保健, 谢永利. 压实黄土连续加载K0固结特性[J]. 岩石力学与工程学报, 2007, 26(增刊1): 3203-3208CHENG Haitao, LIU Baojian, XIE Yongli. K0 consolidation characteristics of compacted loess under continuous loading[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(S1): 3203-3208 侯伟,姚仰平,崔文杰. K0超固结土的不排水抗剪强度[J]. 力学学报,2008,40(6): 795-803. doi: 10.3321/j.issn:0459-1879.2008.06.010HOU Wei, YAO Yangping, CUI Wenjie. Undrained shear strength of K0 over consolidated soil[J]. China Journal of Theoretical and Applied Meachanics, 2008, 40(6): 795-803. doi: 10.3321/j.issn:0459-1879.2008.06.010 刘俊卿,李倩,王保实. 考虑横观各向同性疲劳损伤的 沥青路面数值模型[J]. 西南交通大学学报,2016,51(4): 1-9.LIU Junqing, LI Qian, WANG Baoshi. Numerical model of asphalt pavement considering transversely isotropic fatigue damage[J]. Journal of Southwest Jiaotong University, 2016, 51(4): 1-9. 殷宗泽. 土体本构模型剖析[J]. 岩土工程学报,1996,18(4): 95-97. doi: 10.3321/j.issn:1000-4548.1996.04.016YIN Zongze. Analysis of soil constitutive model[J]. Chinese Journal of Geotechnical Engineering, 1996, 18(4): 95-97. doi: 10.3321/j.issn:1000-4548.1996.04.016 陈正汉. 重塑非饱和黄土的变形、强度、屈服和水量变化特性[J]. 岩土工程学报,1999,21(1): 82-90. doi: 10.3321/j.issn:1000-4548.1999.01.018CHEN Zhenghan. Deformation,strength,yield and moismre change of a remolded unsaturated loess[J]. Chinese Journal of Geotechnical Engineering, 1999, 21(1): 82-90. doi: 10.3321/j.issn:1000-4548.1999.01.018 殷宗泽,周建,赵仲辉,等. 非饱和土本构关系及变形计算[J]. 岩土工程学报,2006,28(2): 137-146. doi: 10.3321/j.issn:1000-4548.2006.02.003YIN Zongze, ZHOU Jian, ZHAO Zhonghui, et al. Constitutive relations and deformation calculation for unsaturated soils[J]. Chinese Journal of Geotechnical Engineering, 2006, 28(2): 137-146. doi: 10.3321/j.issn:1000-4548.2006.02.003 ALONSO E E, GENS A, JOSA A. A constitutive model for partially saturated soils[J]. Géotechnique, 1990, 40(3): 405-430. doi: 10.1680/geot.1990.40.3.405 SHENG D, FREDLUND D G, GENS A. A new modelling approach for unsaturated soils using independent stress variables[J]. Canadian Geotechnical Journal, 2008, 45(4): 511-534. doi: 10.1139/T07-112 卢再华,陈正汉,蒲毅斌. 原状膨胀土损伤演化的三轴CT试验研究[J]. 水利学报,2002(6): 106-112. doi: 10.3321/j.issn:0559-9350.2002.06.019LU Zaihua, CHEN Zhenghan, PU Yibin. Study on damage evolution of natural expansive soil with computerized tomography during triaxial shear test[J]. Shuili Xuebao, 2002(6): 106-112. doi: 10.3321/j.issn:0559-9350.2002.06.019 朱元青,陈正汉. 原状Q3黄土在加载和湿陷过程中细观结构动态演化的CT-三轴试验研究[J]. 岩土工程学报,2009,31(8): 1219-1228. doi: 10.3321/j.issn:1000-4548.2009.08.011ZHU Yuanqing, CHEN Zhenghan. Experimental study on dynamic evolution of meso-structure of intact Q3 loess during loading and collapse using CT and triaxial apparatus[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(8): 1219-1228. doi: 10.3321/j.issn:1000-4548.2009.08.011 -