Spatial Development Characteristics of Post-Fire Debris Flow in Bajiaolou Town
-
摘要: 作为一种与林火密切相关的特殊泥石流，火后泥石流的形成具有明显的空间差异性，对其空间发育特征及影响因素的分析有利于此类泥石流的防灾减灾. 运用现场调查、室内外试验和遥感解译等手段，对雅江县八角楼乡火烧迹地15条沟道的火后泥石流发育特征及地形、林火、植被特征、降雨条件进行了多元分析. 结果表明：火烈度越强，土壤的物理力学性质劣化越严重，斥水性越强，渗透性越低，泥石流越易形成；流域过火面积越大，植被覆盖率越低，泥石流暴发频率越高；相对于大流域，小而陡的流域反而更易暴发火后泥石流. 此外，研究区内暴发的5次火后泥石流中有4次泥石流其激发雨强重复周期均小于5 a，且随着时间的增长激发雨强也逐渐升高.Abstract: As a special type of debris flows closely relevant to forest fire, the generation of post-fire debris flow often shows a spatial variation. It is of great significance to analyze the spatial development characteristics and its influencing factors for hazard mitigation and reduction. By means of field investigation, in-situ and laboratory tests, and remote sensing interpretation, multi-factor analysis was conducted on the 15 channels in the burned area in Bajiaolou town of Yajiang county in terms of their spatial development characteristics, topography, forest fire, vegetation and rainfall conditions. The results show that the higher the fire intensity, the stronger the soil water repellency and the lower the permeability; moreover, the physical and mechanical properties of soil deteriorate more seriously, which contribute to the generation of post-fire debris flow. The occurrence frequency of post-fire debris flow rises when the burned area become larger and the vegetation cover coefficient lower. In contrast to the large basins, the smaller basins show more vulnerability to post-fire debris flow for their short channel lengths and high longitudinal slopes. Besides, the recurrence time of the triggering rainfall intensity that corresponds to four of five debris-flow events in the study areas is less than 5 years, and the triggering rainfall intensity also increases with time.
-
与普通混凝土箱梁相比，波形钢腹板组合箱梁（composite box girders with corrugated steel webs, CBGCSWs）具有自重轻、预应力导入效率高、完全避免腹板开裂等诸多优势[1-3]，因此在我国得以迅速推广. 近十几年来，国内外学者对CBGCSWs的抗弯[4-5]、抗剪[6-7]等方面的受力性能进行了广泛研究. 然而，聚焦其纯扭性能方面的研究却较为匮乏. 事实上，用轻薄的波形钢腹板取代传统的混凝土腹板后，CBGCSWs的抗扭刚度会被大幅削弱（大约仅为传统混凝土箱梁的30%~40%）[8]. 因此，在扭转效应突出的情形下（如曲线梁桥、偏心荷载等），该结构的抗扭问题值得重点关注和深入研究.
针对波形钢腹板PC组合箱梁（prestressed concrete composite box girders with corrugated steel webs, PCCBGCSWs）纯扭作用下的受力问题，诸多学者均提出过相应的分析模型. Mo等[9]将变角软化桁架理论（rotating angle softened truss model for torsion, RA-STMT）应用于PCCBGCSWs，提出了一种预测该结构纯扭作用下力学性能的理论分析模型. 基于RA-STMT，同时考虑结构的受力特点，聂建国等[10]建立了PCCBGCSWs在纯扭作用下的非线性分析模型. 丁勇等[11]以扭转理论及固角软化桁架模型为基础，建立了PCCBGCSWs 在纯扭作用下的抗扭承载力计算模型. Ko等[12]提出了预测PCCBGCSWs纯扭受力行为的改进分析模型，该模型考虑了混凝土抗拉钢化效应，并对混凝土翼缘板内的剪力流有效厚度进行了修正. 但作者仅将该模型用于混凝土开裂后阶段的分析，而针对开裂前阶段的预测则基于Bredt薄壁构件扭转理论[13]. 沈孔健等[14]将PCCBGCSWs的全过程扭矩-扭率曲线分成混凝土开裂前及开裂后两个阶段，并针对上述两个阶段分别采用扭转刚度修正和RA-STMT进行分析. 基于RA-STMT，同时将混凝土抗拉强度考虑在内，Shen等[15]提出了针对单箱多室PCCBGCSWs纯扭全过程的分析模型. Shen等[16]以软化薄膜元模型（softened membrane model for torsion, SMMT）[17]及RA-STMT为基础，建立了PCCBGCSWs纯扭全过程分析模型. 之后，Shen等[16]对该模型进行了改进，在模型中考虑了预应力效应对初始应力及应变的影响[18].
从以上文献可以看出：已有的PCCBGCSWs纯扭分析模型大都基于RA-STMT. 然而，由于RA-STMT忽略了混凝土的抗拉强度和抗拉刚化效应，也不能考虑开裂混凝土提供的抗剪强度，因此基于该理论提出的分析模型无法预测结构在混凝土开裂前的扭矩-扭率曲线，对混凝土开裂后纯扭力学行为的预测也不够精确. 另外，有部分模型（如Ko模型[12]和沈模型[14]）分别采用不同的计算理论对PCCBGCSWs开裂前、后两个阶段进行分析，这样的处理方式显然不具有理论一致性. 除此之外，已有理论模型的求解程序均包含多个迭代循环，因此导致这些模型的求解效率比较低.
本文基于SMMT提出了用于分析PCCBGCSWs纯扭全过程受力行为的改进软化薄膜元模型（improved softened membrane model for torsion，ISMMT）[19-20]. 该分析模型基于单一软化薄膜元理论，考虑了混凝土抗拉强度及开裂混凝土提供的抗剪强度，且对剪力流有效厚度进行了合理修正. 文献[19]将国内外8根PCCBGCSWs纯扭试件的扭矩-扭率曲线与ISMMT预测的理论曲线进行了比较，初步验证了ISMMT的适用性. 事实上，ISMMT除了能够预测扭矩-扭率曲线，还能对各构件在加载全过程的应变历程进行模拟. 因此，该理论模型的准确性尚需进一步的验证.
本文首先对ISMMT的平衡、变形协调、材料本构方程以及通用求解程序进行介绍；之后，针对ISMMT和同类型分析模型中求解程序过于复杂且耗时的问题，提出当波形钢腹板、预应力及普通钢筋均处于弹性阶段时的简化求解程序框图；最后，完成1根PCCBGCSWs试件的纯扭模型试验，并将模型试验结果与采用ISMMT计算的理论结果进行对比，来进一步验证该理论模型的适用性与准确性.
1. 软化薄膜元模型（SMMT）简介
通过考虑受扭构件中混凝土斜压杆的受压软化效应，Hsu等[13]对空间桁架模型进行了改进，提出了经典的软化桁架模型（softened truss model, STM）. 该模型能准确预测纯扭构件在极限状态下的扭矩与扭转角，但由于其未考虑混凝土的抗拉强度，无法对结构在开裂前的扭转性能进行预测，因而具有局限性[21].
为了克服上述不足，Jeng等[17]将用于预测RC受剪构件受力行为的软化薄膜元模型（softened membrane model, SMM）进行了改进，提出了针对RC受扭构件的分析模型（softened membrane model for torsion, SMMT）. 该模型的创新点体现在：通过修正混凝土本构关系来考虑斜压杆的应变梯度效应，并首次将混凝土的抗拉强度考虑在内. 通过与已有文献中的纯扭试验数据进行对比，表明SMMT能对RC实心纯扭构件的全过程受力行为作出精准预测.
2. 改进软化薄膜元模型（ISMMT）
ISMMT是以软化薄膜元理论为基础，并将PCCBGCSWs的结构和力学特点考虑在内，进而得出的一种预测该结构纯扭全过程受力行为的理论分析模型. 以下对ISMMT的平衡、变形协调及材料本构方程以及通用求解程序进行简要介绍，其详细推导过程参见文献[19-20].
2.1 平衡、变形协调及材料本构方程
2.1.1 平衡方程
对于承受外扭矩的PCCBGCSWs，其截面内会形成连续闭合的剪力流q来平衡外部扭矩，如图1所示. 图中：
A0 为组合箱梁中剪力流中心线所包含的截面面积，A0 = b(h−td) ，b、h分别为组合箱梁梁宽和梁高；td为剪力流有效厚度；T 为箱梁承受的扭矩；τlt 为l-t坐标系中的平均剪应力；s 为箍筋间距.
图 1 纯扭作用下波形钢腹板PC组合箱梁Figure 1. Prestressed concrete composite box girder with corrugated steel webs subjected to pure torsion根据组合箱梁截面平衡状态可得
T={2A0τwtw,τw<τwy,2A0τwytw+2×2As(τlttd−τwtw),τw=τwy, (1) 式中：
As 为在混凝土翼缘板内单独闭合的剪力流q2 所包含的截面面积；τw 、τwy 分别为钢腹板的剪应力和剪切屈服强度；tw为钢腹板的厚度.在混凝土翼缘板内剪力流区域取微元A，该微元受纯剪作用. 根据微元A面内平衡状态可得
12(σ2c+σ1c)+τ21c−σci+ρlfl+ρpfp−ρlifli−ρpifpi=0, (2) 12(σ2c+σ1c)−τ21c+ρtft=0, (3) τlt=12(−σ2c+σ1c), (4) 式中：
σ1c 、σ2c 分别为1与2方向的混凝土平均正应力；σ21c 为1-2坐标系中混凝土的平均剪应力；fl 、ft 、fp 分别为纵筋、箍筋及预应力筋的应力；ρl 、ρt 、ρp 分别为纵筋、箍筋及预应力筋的含筋率，如式（5）；σci 、fli 、fpi 分别为由预应力引起的混凝土（式（6））、纵筋（式（7））及预应力筋的初应力；ρli 、ρpi 分别为扭转前纵筋及预应力筋的含筋率，如式（8）.{ρl=−Albtd,ρt=−Atstd,ρp=−Apbtd, (5) σci=−Ecˉεli=−EcApfpiEcAcc+EsAl, (6) fli=−Esˉεli=−EsApfpiEcAcc+EsAl, (7) {ρli=AlwAcc,ρpi=ApwAcc, (8) 式（5）~（8）中：
Al 、Ap 分别为剪力流区间内的纵筋和预应力筋面积（单块混凝土翼缘板内）；At 为单根箍筋面积；Ec 、Es 分别为混凝土和钢筋的弹性模量；ˉεli 为预应力引起的l方向初应变；Acc=Ac−Alw−Apw ，Ac 为混凝土翼缘板截面面积，Alw 、Apw 分别为全部纵筋和预应力筋的面积.2.1.2 变形协调方程
根据Bredt薄壁构件扭转理论[13]，可得组合箱梁截面变形协调方程为
∮γds=2(γltb+γwhw)=2θA0, (9) 式中：
γlt 为l-t坐标系中的平均剪应变；γw 、hw 分别为钢腹板的剪应变和高度；θ 为组合箱梁扭率.基于莫尔应变圆相关理论，可推导出微元A的3个变形协调方程为
{εl=12(ε2+ε1+γ21),εt=12(ε2+ε1−γ21),γlt=(−ε2+ε1), (10) 式中：
εl 、εt 分别为沿l、t方向钢筋的双轴应变；ε1 、ε2 分别为1、2方向的混凝土平均双轴正应变；γ21 为1-2坐标系中混凝土的平均剪应变.假定混凝土薄膜元中应变是线性分布的，可得变形协调方程为
{ϕ=θsin2α2=θ,td=−ˉε2sϕ, (11) 式中：
ϕ 为混凝土薄膜元应变梯度曲率；α2 为混凝土主压应力方向与纵筋的夹角，取45°；ˉε2s 为2方向的混凝土表面单轴正应变.微元A处于双向受力状态，双轴应变与单轴应变之间的关系为
{ˉε1=ε1+ε2ν12,（12a）ˉε2=ε2,（12b）ˉεl=12(ˉε2+ˉε1+γ21),（12c）ˉεt=12(ˉε2+ˉε1−γ21),（12d） 式中：
ˉε1 、ˉε2 、ˉεl 、ˉεt 分别为对应于ε1 、ε2 、εl 、εt 的单轴应变值；ν12 为纯扭作用下的Hsu/Zhu比，如式(13).ν12={0.8(0.2+850εsf),εsf⩽εy,0.8×1.9,εsf>εy, (13) 式中：
εsf 为纵筋与箍筋中先屈服者的应变；εy 为钢筋屈服应变.2.1.3 材料本构方程
受压混凝土的本构方程如下：
σ2c=−ηk1cζfc, (14) k1c=(−Ecˉε1i2fc+1ζfcˉε2s∫ˉε2s0σ2c(¯ε2,ζ)d¯ε2), (15) ζ={0.9,¯ε1,con=¯ε1+¯ε1i⩽0,5.8√fc1√1+400¯ε1,con(1−|β|24°)⩽0.9且5.8√fc⩽0.9,¯ε1,con=¯ε1+¯ε1i>0, (16) β=12tan−1[γ21/(ε2−ε1)]180π, (17) 式（14）~（17）中：
η 为修正系数，按式（19）计算；k1c 为平均压应力系数：ζ 为软化系数；fc 为混凝土轴心抗压强度标准值；σ2c(ˉε,ζ) 按文献[19]分段选取；ˉε1i=0.625ˉεli ，为预应力引起的1方向混凝土初应变；β 为裂缝旋转角度.受拉混凝土的本构方程如下：
σ1c={12η(ˉEcε1s+σci),ˉˉε1s⩽ˉεcx,η[14(ˉεcx−ˉε1i)ˉˉε1sσci+12ˉˉEc(ˉε1s+ˉε1i−ˉεcx)(1−ˉεcx−ˉε1iˉε1s)],ˉεcx<ˉˉε1s⩽εcr,η{σci4(ˉεcx−ˉε1i)ˉε1s+fcr2(εcr−ˉεcx)ˉε1s+53fcrε0.4crˉε1s[(ˉε1s+ˉε1i)0.6−ε0.6cr]},ˉˉε1s>εcr, (18) 式中：
ˉˉε1s=ˉε1s+ˉε1i ，ˉε1s 为1方向的混凝土表面单轴正应变；εcr=0.00008μ ，为混凝土开裂应变；fcr=λμ(0.31)√fc ，为混凝土开裂应力，λ 、μ 分别为受扭混凝土相较于受剪混凝土的弹性模量、开裂应变放大系数，如式（19）；ˉεcx=ˉε1i−σci2ˉEc 为混凝土消压应变；ˉEc=2fcε0 ；ˉˉEc=fcrεcr−ˉεcx .{λ=μ=0.93,η=0.033√fc+0.73,tc⩽0.91td,λ=μ=1.20,η=0.094√fc+0.43,tc>0.91td且fc⩽47.85MPa,λ=μ=1.129,η=8.45/√fc+0.017,tc>0.91td且fc>47.85MPa. (19) 普通钢筋的本构方程如下：
fs={Esˉεs,ˉεs⩽ˉεn,fy[(0.91−2B)+(0.02+0.25B)ˉεsεy],ˉεs>ˉεn, (20) {ˉεn=εy(0.93−2B),B=1ρ(fcrfy)1.5, (21) 式中：
fs 为嵌入普通钢筋的应力；fy 为裸筋的屈服应力；ˉεs 、ˉεn 为嵌入钢筋的单轴应变和单轴屈服应变.对于嵌入纵筋，
ˉεls=ˉεl+ˉεli .预应力钢筋的本构方程如下：
fp={Epεp,εp⩽εpro,Epεp[1+(Epεpfpu)4.38]−14.38,εp>εpro, (22) 式中：fp为预应力筋的应力；
Ep 、εp 、εpro 、fpu 分别为预应力筋的弹性模量、应变、比例极限应变和极限强度，εp=fpi/Ep+ˉεl .波形钢腹板的本构方程如下：
{τw={Geγw,τw⩽τwy,τy,τw>τwy,Ge=Gs(aw+bw)/(aw+cw), (23) 式中：
Ge 为波形钢腹板的有效剪切模量；Gs 为钢材的剪切模量；aw 、bw 和cw 分别为波形钢腹板的直腹板段、斜腹板段投影和斜腹板段长度.混凝土剪应力与剪应变之间的关系为
τ21c=12(σ1c−σ2c)ε1−ε2γ21+σci2. (24) 需要指出的是，一些已有的PCCBGCSWs纯扭分析模型（如文献[9-10]）忽略了混凝土的抗拉强度，因此仅能预测结构在混凝土开裂后的抗扭行为. 而ISMMT成功地将混凝土抗拉强度考虑在内，从而能够预测PCCBGCSWs包括开裂前上升段在内的整个扭矩-扭率曲线.
2.2 求解程序框图
2.2.1 通用求解程序
联立式（9）、（11）可求得
td 的表达式，但可能出现计算所得td 大于混凝土翼缘板厚度tc 的不实际情况. 因此，对td 进行判定并修正，如式（25）.td={−2ˉε2bh(γltb+γwhw−2ˉε2b),td⩽tc,tc,td>tc. (25) 式（25）的推导是为了避免采用“试错法”迭代计算td，从而提高求解效率. ISMMT的通用求解程序见图2（a），式（26）~（28）为迭代计算的3个收敛判别准则.
ρlfl+ρpfp+ρtft−ρlifli−ρpifpi=−(σ2c+σ1c−σci), (26) ρlfl+ρpfp−ρtft−ρlifli−ρpifpi=−(2τ21c−σci), (27) {τlttd=τwtw,τw<τwy,γw=γlt,τw=τwy. (28) 2.2.2 简化求解程序
Jeng等[17]给出的SMMT求解程序适用于PC构件的纯扭全过程分析，但其中包含2层嵌套迭代循环，较为复杂，求解耗时也较长. 为提高求解效率，Jeng等在此基础上提出了当普通钢筋处于弹性阶段时的简化求解程序，该简化程序仅包含1层迭代循环[22]. 同样地，图2（a）所示的通用求解程序适用于PCCBGCSWs纯扭全过程分析，但包含3层嵌套迭代循环，求解更为复杂与耗时. 因此，本文在充分考虑PCCBGCSWs结构特点的基础上，参考Jeng等[17]提出的简化模型，提出了当波形钢腹板、预应力及普通钢筋均处于弹性阶段时的简化程序框图（图2（b））. 当纵筋、箍筋及预应力钢筋均处于弹性阶段时，其本构关系可分别表示为
{fl=Es(ˉεl+ˉεli),ft=Esˉεt,fp=Ep(ˉεl+εpi). (29) 将式（12c）、（12b）、（24）、（29）代入式（26）和式（27），可得
AlEs2btd(ˉεli+12(ˉε2+ˉε1+γ21))+AtEs2std(ˉε2+ˉε1−γ21)+ApEp2btd(εpi+12(ˉε2+ˉε1+γ21))−ρlifli−ρpifpi=−(σ2c+σ1c−σci), (30) AlEs2btd(ˉεli+12(ˉε2+ˉε1+γ21))−AtEs2std(ˉε2+ˉε1−γ21)+ApEp2btd(εpi+12(ˉε2+ˉε1+γ21))−ρlifli−ρpifpi=(σ2c−σ1c)ε1−ε2γ21. (31) 显然，
γ21 可通过求解式（30）、（31）得出，求解γ21 的目的是为了消除与γ21 相关联的迭代循环. 此外，由于波形钢腹板也处于弹性阶段，与γw 相关联的迭代循环可通过代入γw=τlttd/τlttdGestwGestw （根据式（23）和q=τlttd=τwtw 求出）予以消除. 至此，通过上述简化，求解程序仅包含1层迭代循环，如图2（b）所示. 需要注意的是，图2（a）所示的通用求解程序适用于纯扭全过程分析，而图2（b）所示的简化求解程序仅适用于当波形钢腹板、预应力及普通钢筋均处于弹性阶段时的纯扭分析.与通用求解程序相比，该简化求解程序的计算效率有显著提高，主要原因为：1） 通用求解程序与简化求解程序分别包含3层嵌套迭代循环和1层迭代循环，其计算复杂度可分别用O(
n3 )和O(n)来表示（n为问题的维度）. 因此，后者的求解效率要显著高于前者；2） 通用求解程序采用“试错法”来假定未知变量γ21 、ε1 和γw 的值，具有盲目性，求解效率低. 而在简化求解程序中，当收敛判别式不满足后，程序可自动生成下一个迭代循环所需的初始值ˉε2 ，从而避免采用“试错法”盲目寻找，可有效提高求解效率. 同时，与本文通用求解程序类似，已有的同类型分析模型（文献[9-12、14-16、18]中模型）均包含3层嵌套迭代循环. 因此，本文简化求解程序的求解效率同样远高于已有同类型分析模型.根据图2所示求解程序框图，本文采用MATLAB编制了PCCBGCSWs纯扭全过程分析模型，运行程序可得到一组扭矩与扭率的离散点，最终绘成扭矩-扭率图.
2.3 ISMMT与已有分析模型比较
与文献中已有的同类型分析模型相比，ISMMT有着较为显著的优势. 首先，由于ISMMT基于单一软化薄膜元理论提出，因而具有理论一致性. 同时，SMMT成功地将混凝土抗拉强度、抗拉刚化效应以及开裂混凝土提供的抗剪强度考虑在内. 本文ISMMT基于SMMT提出，同样能够考虑上述效应，因而能准确模拟结构包括开裂前阶段在内的全过程扭转行为. 此外，与已有分析模型[9-12,14-16,18]的求解程序相比，本文提出的简化求解程序效率更高、求解耗时更少. 因此，ISMMT为纯扭作用下PCCBGCSWs的受力分析提供了有效途径.
3. 模型试验
3.1 试件概况
设计并制造了一根等截面PCCBGCSWs试件，如图3所示. 图3（a）为试件的立面图，梁的长度及高度分别为2.8 m和0.55 m， 为方便固定与加载，在梁的两端各布置一道0.4 m厚的横隔板. 图3（b）为组合箱梁横截面尺寸和预应力及普通钢筋布置示意图，普通钢筋型号均为HRB335，规格均为Φ10.0，箍筋间距为100 mm，混凝土保护层厚度均为20 mm. 预应力钢束采用单根直径15.2 mm、抗拉强度标准值1 860 MPa的低松弛钢绞线，在加载前施加初应力. 波形钢腹板采用Q235钢材，图3（c）为其尺寸示意图. 为了确保波形钢腹板与混凝土顶、底板连接处在加载过程中不提前发生破坏，剪力连接件采用刚度较大的双PBL形式. 此外，模型试件的材料特性列于表1.
表 1 试件材料参数Table 1. Material properties of the test beam混凝土棱柱体抗压强度/MPa 波形钢腹板 预应力钢筋初应力/MPa 普通钢筋 厚度/mm 屈服强度/MPa 直径/mm 屈服强度/MPa 24.9 3.7 235 800 10 357 3.2 加载装置
自行设计了一套纯扭加载装置，如图4所示. 试件的一端通过千斤顶进行固定，另一端通过加载梁夹紧后作为转动端，在转动端底部放置铰支座. 在压力传感器与反力梁之间放置球铰，采用液压千斤顶在加载梁上施加竖向偏心荷载，可使得试件在加载过程中处于近乎纯扭的受力状态. 加载梁的悬臂段长1.65 m，通过压力传感器对竖向荷载的大小进行控制. 由于加载梁自重较大，考虑加载梁偏心对试件产生的初始扭矩.
3.3 测试方案
3.3.1 扭矩与扭转角
扭矩可近似按竖向偏心荷载与加载梁力臂的乘积计算，该近似值与考虑扭转角的精确值相比相差较小[23]. 此外，在试件的1/2及3/4截面分别布置位移千分表来测量试件在加载过程中的扭转竖向位移. 按式（32）换算得到该截面的扭转角，图5为计算示意图.
[L2c+(δc−h2)2]sin2αc−Lchsinαc−δ2c+δch=0, (32) 式中：
Lc 为梁中心至千分表的距离；δc 为千分表测得的扭转竖向位移；αc 为试件的扭转角度.3.3.2 测点布置
试件的测点布置情况为：在跨中截面左、右两侧波形钢腹板上各布置6组（共12组）应变花，编号分别为LW-1~LW-6（左侧）和RW-1~RW-6（右侧）；跨中截面混凝土顶、底板沿横向分别均匀布置5组（共10组）应变花，编号分别为TC-1~TC-5（顶板）和BC-1~BC-5（底板）；预应力钢束采用后张法进行张拉（未灌浆，属无黏结预应力），沿对角选取2根钢束在其张拉端布置预应力传感器，编号分别为P-ZS及P-YX（如图4所示）；在混凝土顶、底板内各选取8个普通钢筋（包括纵筋和箍筋）测点（共16个）布置应变片，编号分别为TG-1~TG-8（顶板）及BG-1~BG-8（底板）. 上述主要测点的详细布置情况如图6所示.
4. 结果对比
为验证ISMMT能否准确预测PCCBGCSWs中各构件在纯扭作用下的应变历程，将ISMMT预测的各项理论结果与对应的试验结果进行对比，如图7所示.
4.1 扭矩-扭率曲线
将扭矩-扭率曲线的理论和试验结果列于图7（a），其中，扭率的试验结果取所有4个竖向位移测点按式（27）换算得到扭率值的平均值. 从图7（a）可以看出：当试件处于弹性阶段时，扭率试验结果要略大于对应的理论结果，造成该现象的原因可能是该阶段扭率值相对较小，千分表读数易受周围环境的干扰. 但由于本文仅完成了1根模型梁的纯扭试验，扭率的试验与理论结果存在偏差的原因还有待进一步研究. 在混凝土开裂后，理论与试验曲线遵循着相同的发展趋势，且吻合良好. 此外，将理论和试验曲线在纯扭加载过程中关键节点处的对应结果列于表2~4. 表中：T 和 θ 分别为扭矩和扭率，下标第1个数（1、2、3）代表状态（开裂状态、屈服状态和极限状态），下标第2个数1、2分别代表理论值与试验值. 表中结果表明：除开裂状态下的扭率试验值与对应的理论结果相差较大外，其他结果均十分接近. 由此可知：ISMMT能准确预测PCCBGCSWs在纯扭状态下的全过程扭矩-扭率曲线.
4.2 波形钢腹板剪应变
图7（b）为由ISMMT和试验给出的波形钢腹板剪应变结果（绝对值）. 从图中可以看出：理论与试验曲线吻合良好，且由分析模型给出波形钢腹板平均剪应变位于各分散试验曲线的区间内；左侧与右侧波形钢腹板对应位置处的剪应变值很接近，且在凹进（LW-1~LW4、RW-1~RW4）与凸出（LW-5~LW6、RW-5~RW6）直钢板处的剪应变值也相差很小. 在钢腹板屈服前，波形钢腹板上的剪应变随着外扭矩的增加而近乎线性增加. 而当钢腹板达到其屈服强度时，由于模型梁已基本达到其抗扭极限承载力，因此尽管波形钢腹板上的剪应变迅速增加，施加在试件上的外扭矩仍基本保持不变.
4.3 混凝土翼缘板剪应变
图7（c）给出了混凝土翼缘板剪应变（绝对值）的理论和试验结果对比情况. 图中结果表明：当模型梁处于弹性阶段时，混凝土翼缘板上的剪应变随着外扭矩的增加而近乎线性增加，但增长幅度很小；在该阶段内，理论结果与试验结果高度吻合；在混凝土开裂后，混凝土顶、底板上的剪应变迅速增加，尽管此时试验曲线较为分散，但由ISMMT给出混凝土翼缘板平均剪应变仍处于各分散试验曲线的区间内，且二者遵循着同样的变化趋势. 造成精度下降的主要原因在于混凝土裂缝出现位置及开裂角度的不确定性.
图 7 ISMMT预测结果与试验结果对比Figure 7. Comparison of the results obtained from the ISMMT and experiment表 2 ISMMT预测结果与试验结果对比（开裂状态）Table 2. Comparison of the predicted torques and twists from the ISMMT and the experiment (cracking state)参数 T11/
(kN•m)T12/
(kN•m)T11T12 θ11/
((°)•m−1)θ12/
((°)•m−1)θ11θ12 取值 160.6 174.5 0.92 0.073 0.161 0.45 表 3 ISMMT预测结果与试验结果对比（屈服状态）Table 3. Comparison of the predicted torques and twists from the ISMMT and the experiment (yield state)参数 T21/
(kN•m)T22/
(kN•m)T21T22 θ21/
((°)•m−1)θ22/
((°)•m−1)θ21θ22 取值 351.9 337.9 1.04 1.043 0.973 1.07 表 4 ISMMT预测结果与试验结果对比（极限状态）Table 4. Comparison of the predicted torques and twists from the ISMMT and the experiment (limit state)参数 T31/
(kN•m)T32/
(kN•m)T31T32 θ31/
((°)•m−1)θ32/
((°)•m−1)θ31θ32 取值 354.9 339.5 1.04 1.547 1.590 0.97 4.4 预应力钢束应变
将由ISMMT和试验给出的预应力钢束应变结果列于图7（d）. 图中结果表明：理论与试验曲线吻合较好且遵循着同样的发展规律，在混凝土开裂前，预应力钢束的应变基本保持不变；当所施加的外扭矩达到开裂扭矩后，混凝土顶底板上斜裂缝的数量和宽度随着外扭矩的增加而不断扩张，从而导致预应力钢束应变迅速增加.
4.5 纵筋及箍筋应变
图7（e）、（f）分别给出了纵筋和箍筋应变的理论和试验结果对比. 从图中可以看出：由于试验误差等原因，纵、箍筋各测点的应变试验结果较为离散，但总体保持着一致的发展趋势，且理论与试验结果也较为吻合；预应力的存在会使得模型梁沿桥轴向产生初始应变. 因此，在正式施加纯扭荷载前，纵筋处于受压状态（如图7（e））；图7（d） 中预应力钢束应变的变化趋势和幅值与图7（e）中纵筋应变一致，该试验现象与理论模型中的假设吻合，从而进一步证明了理论模型的准确性.
4.6 求解效率
当采用ISMMT预测试件的纯扭全过程受力行为时，在波形钢腹板、预应力及普通钢筋均处于弹性阶段时分别采用简化求解程序和通用求解程序进行计算，并对二者的求解耗时进行对比. 其中，两种求解程序的收敛误差均控制在0.1%以内. 结果表明，对于某指定的
ε2 ，采用通用求解程序求得γ21 、ε1 和γw 真实值所需的总迭代次数最少为4.9 × 105次，最多达7.9 × 106次（需要说明的是，该迭代次数与γ21 、ε1 和γw 所假定的迭代步长有关）. 而采用简化求解程序进行计算时，对于某指定的ˉε1 ，求得所有未知变量的最少迭代次数为5次，最多也仅为193次. 由此可见，采用本文提出的简化求解程序可以极大地提高求解效率.5. 结 论
在软化薄膜元理论的基础上提出了改进软化薄膜元模型（ISMMT）来预测PCCBGCSWs在纯扭作用下的全过程受力行为. 针对ISMMT以及已有同类型理论模型中求解程序所存在的迭代循环多、求解效率低等问题，提出了当波形钢腹板、预应力及普通钢筋均处于弹性阶段时的简化求解程序框图. 此外，还完成了纯扭模型试验对ISMMT的适用性与准确性进行验证. 主要结论如下：
1） 本文所提出的简化求解程序框图仅有一层迭代循环，其求解效率要远高于国内外已有同类型分析模型.
2） 由模型试验得到的试件扭矩-扭率曲线、波形钢腹板和混凝土翼缘板剪应变、预应力及普通钢筋应变与ISMMT预测的理论结果吻合良好，表明ISMMT除了能准确预测PCCBGCSWs的全过程扭矩-扭率曲线外，还能模拟混凝土翼缘板、波形钢腹板、预应力和普通钢筋等构件的整个应变发展历程. 该理论模型为更全面地了解PCCBGCSWs的纯扭力学性能提供了有效途径.
-
图 1 研究区火烧迹地及火后泥石流发育工程地质平面
Figure 1. Burnt area and geological development plane of post-fire debris flow
图 4 各火烈度区土壤内聚力和内摩擦角
Figure 4. Cohesion and internal friction of soil in areas with different fire intensities
图 5 各火烈度区坡表不同深度处土壤渗透性
Figure 5. Soil permeability at different depths in areas with different fire intensities
图 6 火烧区面积占比与泥石流暴发频次的关系
Figure 6. Relationship between proportion of burned area and frequency of post-fire debris flow
图 9 植被覆盖率与火后泥石流暴发频次的关系
Figure 9. Relationship between vegetation coverage rate and frequency of post-fire debris flow
图 10 地形因子及火后泥石流爆发频次与流域面积的关系
Figure 10. Relationship between terrain factors，post-fire debris flow frequency and basin area
表 1 2018年火烧迹地泥石流暴发情况统计
Table 1. Statistics of post-fire debris flow occurred in burnt area in 2018
流域 是否发生泥石流和洪水 是否过火 暴发频次/次 易发性 6月13日 6月24日 6月30日 7月5日 7月12日 D1 × × × √ √ 是 2 L D2 √ × √ √ √ 是 4 H D3 √ √ √ √ √ 是 5 H D4 × × √ √ √ 是 3 L D5 √ √ √ √ √ 是 5 H D6 × × √ √ √ 是 3 L D7 √ √ √ √ √ 是 5 H D8 × × √ √ √ 是 3 L D9 √ √ √ √ √ 是 5 H D10 √ √ √ √ √ 是 5 H F1 × × × × × 是 0 N F2 × × × × × 是 0 N F3 × × × × × 是 0 N F4 × × × × × 否 0 N F5 × × × × × 否 0 N 注：“√”表示发生了泥石流，“×”表示只发生了洪水；“L”表示易发，“H”表示高易发，“N”表示不易发. 
下载: 导出CSV
表 2 火烧迹地坡表土壤物理性质试验结果
Table 2. Tested physical properties of soil on slope in study area
火烈度 密度/（g•cm−3） 含水率/% 比重 干密度/（g•cm−3） 孔隙度/% 未火烧（对照组） 1.51 ± 0.08 （a） 4.77 ± 0.30 （a） 2.62 ± 0.01 （a,b） 1.44 ± 0.08 （a） 45.14 ± 3.08 （b） 轻度 1.30 ± 0.03 （a,b） 3.89 ± 0.44 （a,b） 2.62 ± 0.01 （a,b） 1.25 ± 0.03 （a,b） 52.28 ± 1.28 （a,b） 中度 1.20 ± 0.08 （b） 2.69 ± 0.34 （b） 2.63 ± 0.01（a） 1.17 ± 0.08 （b） 55.66 ± 3.17 （a） 重度 1.14 ± 0.04 （b） 1.14 ± 0.26 （c） 2.59 ± 0.01 （b） 1.12 ± 0.04 （b） 56.66 ± 1.40 （a） 注：表中数据为均值 ± 标准误差；a、b、c反映单因素方差分析结果，同列数据不同字母表示均值差异性显著（P值小于0.05）.
下载: 导出CSV
表 3 各火烈度区坡表不同深度处不同斥水土壤样本含量占比
Table 3. Proportion of hydrophobic soil samples at different depths in areas with different fire intensities
深度/cm 火烈度 样本含量占比/% 斥水性土壤
占比合计/%N L M H 0 W 72.48 23.91 2.45 1.16 27.52 Q 42.22 26.93 7.81 23.04 57.78 Z 29.79 28.61 14.62 26.98 70.21 Y 25.42 28.59 12.05 33.94 74.58 1 W 93.04 4.16 1.53 1.27 6.96 Q 93.53 2.57 0.86 3.04 6.47 Z 51.07 32.01 7.64 9.28 48.93 Y 52.24 33.81 5.97 7.98 47.76 2 W 97.05 1.23 0.44 1.28 2.95 Q 97.12 1.54 0.53 0.81 2.88 Z 88.21 6.53 2.64 2.62 11.79 Y 87.53 7.62 1.56 3.29 12.47 3 W 97.13 0.82 1.21 0.84 2.87 Q 96.43 1.68 0.87 1.02 3.57 Z 96.49 1.07 0.99 1.45 3.51 Y 94.35 1.86 1.98 1.81 5.65 注：“W”为未火烧，“Q”为轻度火烧，“Z”为中度火烧，“Y”为重度火烧；“N”为亲水，“L”为轻度斥水，“M”为中度斥水，“H”为严重斥水.
下载: 导出CSV
表 4 各流域基本地形参数
Table 4. Basic topographic parameters of each basin
流域编号 面积/km2 流域高
差/m主沟长/km 主沟纵比降/‰ 沟道密度/
（km•km−2）D1 3.19 1430 3.88 336.04 1.87 D2 0.85 825 2.00 370.61 3.58 D3 0.25 650 1.28 501.70 8.80 D4 2.20 1445 3.35 395.28 2.67 D5 0.20 624 0.95 717.63 9.93 D6 0.31 761 1.17 722.32 6.42 D7 0.64 870 1.48 555.21 4.95 D8 2.82 1077 2.72 335.14 3.13 D9 0.31 592 0.97 590.58 6.29 D10 0.24 492 0.83 549.63 8.21 F1 4.13 1449 5.06 286.82 1.80 F2 6.98 1501 6.45 206.52 2.54 F3 5.15 1443 4.34 289.46 1.48 F4 0.44 697 1.13 581.34 6.18 F5 1.38 1062 2.81 348.70 3.01
下载: 导出CSV
-
NYMAN P, SMITH H G, SHERWIN C B, et al. Predicting sediment delivery from debris flows after wildfire[J]. Geomorphology, 2015, 250: 173-186. doi: 10.1016/j.geomorph.2015.08.023 NYMAN P, SHERIDAN G J, SMITH H G, et al. Evidence of debris flow occurrence after wildfire in upland catchments of south-east Australia[J]. Geomorphology, 2011, 125(3): 383-401. doi: 10.1016/j.geomorph.2010.10.016 GARTNER J E, CANNON S H, BIGIO E R, et al. Compilation of data relating to the erosive response of 606 recently-burned basins in the Western U. S.[R]. Washington D. C.: U. S. Geological Survey, 2005. GABET E J, BOOKTER A. A morphometric analysis of gullies scoured by post-fire progressively bulked debris flows in southwest Montana,USA[J]. Geomorphology, 2008, 96(3/4): 298-309. CANNON S H, GARTNER J E, JAKOB M, et al. Wildfire-related debris flow from a hazards perspective[M]. Berlin: Springer, 2005: 363-385. HYDE K D, WILCOX A C, JENCSO K, et al. Effects of vegetation disturbance by fire on channel initiation thresholds[J]. Geomorphology, 2014, 214: 84-96. doi: 10.1016/j.geomorph.2014.03.013 KEAN J W, STALEY D M, CANNON S H. In situ measurements of post-fire debris flows in southern California:comparisons of the timing and magnitude of 24 debris-flow events with rainfall and soil moisture conditions[J]. Journal of Geophysical Research: Earth Surface, 2011, 116(4): F4019.1-F4019.21. 胡卸文,王严,杨瀛. 火后泥石流成灾特点及研究现状[J]. 工程地质学报,2018,26(6): 1562-1573.HU Xiewen, WANG Yan, YANG Ying. Research actuality and evolution mechanism of post-fire debris flow[J]. Journal of Engineering Geology, 2018, 26(6): 1562-1573. 胡卸文,韩玫,梁敬轩,等. 汶川地震灾区泥石流若干关键问题[J]. 西南交通大学学报,2016,51(2): 331-340. doi: 10.3969/j.issn.0258-2724.2016.02.012HU Xiewen, HAN Mei, LIANG Jingxuan, et al. Some key problems on debris flow in Wenchuan earthquake area[J]. Journal of Southwest Jiaotong University, 2016, 51(2): 331-340. doi: 10.3969/j.issn.0258-2724.2016.02.012 邸雪颖,陶玉柱. 火后泥石流研究进展[J]. 应用生态学报,2013,24(8): 2383-2392.DI Xueying, TAO Yuzhu. Research progress in post-fire debris flow[J]. Chinese Journal of Applied Ecology, 2013, 24(8): 2383-2392. 陈子龙. 仁额拥沟火后泥石流发育特征及成因机制研究[D]. 成都: 成都理工大学, 2016. 任云,胡卸文,王严,等. 四川省九龙县色脚沟火后泥石流成灾机理[J]. 水文地质工程地质,2018,45(6): 150-156.REN Yun, HU Xiewen, WANG Yan, et al. Disaster mechanism of the Sejiao post-fire debris flow in Jiulong county of Sichuan[J]. Hydrogeology and Engineering Geology, 2018, 45(6): 150-156. KEY C H, BENSON N C. Landscape assessment (LA) sampling and analysis methods[M]. Washington D. C.: USDA Forest Service, 2006: 25-26. 王晓莉,王文娟,常禹,等. 基于NBR指数分析大兴安岭呼中森林过火区的林火烈度[J]. 应用生态学报,2013,24(4): 967-974.WANG Xiaoli, WANG Wenjuan, CHANG Yu, et al. Fire severity of burnt area in Huzhong forest region of Great Xing’an Mountains,Northeast China based on normalized burn ration analysis[J]. Journal of Applied Ecology, 2013, 24(4): 967-974. 中华人民共和国水利部. 土工试验方法标准: GB/T 50123—1999[S]. 北京: 中国计划出版社, 1999. BADÍA-VILLAS D, GONZÁLEZ-PÉREZ J A, AZNAR J M, et al. Changes in water repellency,aggregation and organic matter of a mollic horizon burned in laboratory:Soil depth affected by fire[J]. Geoderma, 2014, 213: 400-407. doi: 10.1016/j.geoderma.2013.08.038 ROBICHAUD P R, WAGENBRENNER J W, PIERSON F B, et al. Infiltration and interrill erosion rates after a wildfire in western Montana,USA[J]. Catena, 2016, 142: 77-88. doi: 10.1016/j.catena.2016.01.027 张颖,牛健植,谢宝元,等. 森林植被对坡面土壤水蚀作用的动力学机理[J]. 生态学报,2008,28(10): 5084-5094. doi: 10.3321/j.issn:1000-0933.2008.10.056ZHANG Ying, NIU Jianzhi, XIE Baoyuan, et al. Dynamics mechanism of the effect of forest vegetation on hill-slop water erosion[J]. Journal of Applied Ecology, 2008, 28(10): 5084-5094. doi: 10.3321/j.issn:1000-0933.2008.10.056 陈晋,陈云浩,何春阳,等. 基于土地覆盖分类的植被覆盖率估算亚像元模型与应用[J]. 遥感学报,2001,5(6): 416-422. doi: 10.11834/jrs.20010603CHEN Jin, CHEN Yunhao, HE Chunyang, et al. Sub-pixel model for vegetation fraction estimation based on land cover classification[J]. Journal of Remote Sensing, 2001, 5(6): 416-422. doi: 10.11834/jrs.20010603 NGUYEN P, SHEARER E J, TRAN H, et al. The CHRS data portal,an easily accessible public repository for Persiann global satellite precipitation data[J]. Nature Scientific Data, 2019, 6: 180296.1-180296.10. doi: 10.1038/sdata.2018.296 期刊类型引用(8)
1. 曹希超，胡卸文，金涛，周瑞宸，周永豪，何坤，胡亚运. 玉溪市2023年“4·11”森林火灾火烧迹地坡面侵蚀及火后泥石流动储量预测. 工程地质学报. 2025(01): 258-268 . 
百度学术2. 林尧，李涛，章浩淳，王志恒，郑伟皓，唐智慧. 基于复杂网络的火后泥石流成灾机理的建模方法. 防灾科技学院学报. 2024(02): 75-86 . 百度学术3. 胡卸文，周永豪，何坤，曹希超，胡亚运. 不同岩性区火后泥石流成灾机理. 山地学报. 2024(04): 535-545 . 百度学术4. 周瑞宸，胡卸文，金涛，曹希超，周永豪，张瑜. 重庆市2022年8月森林火灾火烧区特点及火后泥石流易发性评价. 水文地质工程地质. 2024(05): 150-160 . 百度学术5. 魏凌浩，张可能，胡丰云，李斌田，王剑，伍联合. 薄土层高植被盖度区域低频泥石流成因分析. 四川建筑. 2023(03): 139-141 . 百度学术6. 黄健，胡卸文，金涛，曹希超，杨相斌. 四川西昌“3·30”火烧区响水沟火后泥石流成灾机理. 中国地质灾害与防治学报. 2022(03): 15-22 . 百度学术7. 吉日伍呷，田宏岭，刘雅琴. 喜德县“6.4”山火泥石流灾害特征与防治建议. 科学技术与工程. 2022(17): 6880-6887 . 百度学术8. 卢洪沛，胡卸文. 森林火灾后泥石流的成因与特征分析. 四川水力发电. 2021(05): 1-5+127 . 百度学术其他类型引用(4)
-
下载:
下载:
百度学术
点击查看大图
计量
- 文章访问数: 457
- HTML全文浏览量: 229
- PDF下载量: 18
- 被引次数: 12