目前在实际工程中常用的桥面板形式按照材料划分主要有混凝土桥面板和正交异性钢桥面板.混凝土桥面板构造简单,造价低廉,在小跨径桥梁中被大量使用,但是自重较大,限制了其在中、大跨径桥梁中的应用.正交异性钢桥面板自重小,承载能力高,非常适合在中、大跨径桥梁中应用,但是疲劳裂纹和铺装问题一直是正交异性钢桥面板的“顽疾”[1-2].钢混凝土组合桥面板的出现为解决正交异性钢桥面板的疲劳和铺装问题提供了另外一种思路[3].混凝土板通过剪力连接件与正交异性钢桥面板连接在一起,混凝土板参与结构受力,这样降低了钢结构应力,提高了抗疲劳性能,同时也解决了桥面板的铺装问题.
当组合桥面板应用在连续梁中时,负弯矩区混凝土板开裂的问题不可避免.为了保证结构的耐久性,避免钢筋和下层钢结构的锈蚀,除了在铺装层中增设防水层等构造措施外,限制混凝土板的裂缝宽度也非常必要.目前混凝土构件受拉裂缝宽度的研究已经比较成熟,其裂缝宽度计算公式主要分为以下两种类型:一种是半经验半理论公式,以《混凝土结构设计规范》(GB 50010—2010) [4](后简称混凝土规范)为代表;另一种是经过数理统计得到的经验公式,以《公路钢筋混凝土及预应力混凝土桥涵设计规范》(JTG D62—2004) [5](后简称公预规)为代表.也有学者对组合梁负弯矩区裂缝宽度的计算进行了相关研究.聂建国等[6]提出组合梁负弯矩区混凝土裂缝宽度主要与配筋率、力比和栓钉连接件间距有关.余志武等[7]根据18根部分预应力混凝土组合梁试验得出,负弯矩区裂缝宽度主要与负弯矩区综合力比、剪力连接件间距和钢梁与混凝土的相对高度比有关.张彦玲等[8-9]和Ryu等[10]提出,实测的裂缝间距与横向钢筋间距有很强的相关性,也是影响裂缝宽度的重要因素.上述三位学者提出的裂缝宽度计算公式都是基于混凝土规范中的半经验半理论公式,形式上相同,只是对钢筋应力不均匀系数、裂缝间距的计算方法不同.《公路钢结构桥梁设计规范》(JTG D64—2015) [11](后简称钢桥规范)中规定,可以先按照开裂截面特性计算负弯矩区钢筋应力,再按照公预规中轴心受拉构件计算组合梁负弯矩区裂缝宽度.
对于组合桥面板,虽然与组合梁构造上相似,但是仍有其特殊之处.例如,组合桥面板混凝土板厚度与组合梁相比较薄,一般只布置一层钢筋.组合梁桥面板下部的钢结构为正交异性钢桥面板,其宽度与混凝土板相同,所以剪力滞效应没有组合梁明显.目前还未见有关组合桥面板中混凝土裂缝宽度计算的文献,因此有必要对组合桥面板负弯矩区混凝土裂缝宽度的计算方法进行研究.本文对5个简支组合桥面板试件进行试验研究,分析钢筋应力、连接件形式、钢筋直径、混凝土强度等级等因素对裂缝宽度的影响.最后,提出了适合组合桥面板的负弯矩区钢筋应力计算公式,并基于现有规范提出了适合组合桥面板的负弯矩区混凝土裂缝宽度计算方法.
1 试验研究 1.1 试验试件及加载方案共设计了5个组合桥面板试件,编号为S1~S5.每个试件的长度为4.3 m.以试件S1为例详细介绍试件的结构形式.试件S1截面形式如图 1a所示,主要由混凝土顶板、钢顶板、U形加劲肋以及连接件和钢筋组成.试件的混凝土强度等级为C60,钢材强度等级为Q345,宽度为1 800 mm,高度为406 mm.其中,混凝土板厚120 mm,加劲肋高280 mm,钢顶板厚6 mm,混凝土和钢顶板之间通过焊钉连接,焊钉在试件长度和宽度方向的间距为300 mm;试件宽度范围内有2个U肋,U肋板厚6 mm;混凝土中纵向钢筋直径为22 mm,钢筋间距为150 mm,横向钢筋直径为12 mm,钢筋间距为150 mm.以试件S1为基准试件,其他4个试件与试件S1的不同之处在于:试件S2混凝土纵向钢筋采用直径16 mm的细钢筋;试件S3混凝土采用高强钢纤维混凝土;试件S4加劲肋采用球扁钢,球扁钢的型号为280×11,同时试件宽度改为1 500 mm,截面形状如图 1b所示;试件S5剪力连接件采用开孔板连接件,开孔板连接件的厚度为12 mm,高度为120 mm,孔径为50 mm,孔中心间距为150 mm.5个试件的设计变量如表 1所示.试件的标准断面如图 1所示.
加载方式如图 2所示,试件的支承跨度均为4 m.在试件的两端约束竖向位移,在跨中下部设置一油压千斤顶施加向上的力使得跨中组合桥面板承受负弯矩作用.试件的现场加载如图 3所示.
各试验试件的荷载裂缝宽度曲线如图 4所示.由于S4试件的宽度与其他试件不同,这里为了方便与其他试件比较,将其对应的荷载放大1.2倍进行正则化处理.可以看出,达到相同的裂缝宽度时采用高强钢纤维混凝土的试件S3所能承受的荷载更大,采用φ16钢筋的试件S2比采用φ22钢筋的试件S1裂缝扩展的速度更快.采用球扁钢的试件S4比采用U肋的试件S1裂缝扩展速度略快.同时,试件S1与S5的裂缝扩展速率基本相同,说明剪力连接件的类型对负弯矩区裂缝的扩展影响不大.
各试件在最终破坏时的裂缝分布如图 5所示.
从图 5看出,对于试件S1、S2、S3和S5的混凝土裂缝主要集中在跨中L/2(2 m范围,L为跨度),而试件S4裂缝分布范围较大,在3L/4范围内.对于试件S3,由于采用纤维混凝土,其混凝土裂缝较少且稀疏,裂缝在试件宽度方向也没有形成通长裂缝,而其他试件均在试件宽度方向形成了通长裂缝.
2 混凝土裂缝宽度计算方法 2.1 钢筋应力目前混凝土裂缝宽度的计算公式中都将裂缝截面钢筋应变作为一个主要影响因素考虑在内,因此准确计算开裂截面钢筋应力是混凝土裂缝宽度计算分析的基础.钢桥规范中对于组合梁受弯构件钢筋应力的计算采用如下方法:不考虑混凝土作用的开裂截面模量,按照受弯构件的截面应力计算公式计算.在弯矩M作用下钢筋的应力
${{\sigma }_{\text{ss}}}=\frac{M{{y}_{\text{r}}}}{{{I}_{\text{cr}}}}$ | (1) |
式中:yr为混凝土中受拉钢筋距开裂截面中性轴的距离;Icr为开裂截面惯性矩.
将采用式(1) 计算的钢筋应力换算成应变,与试验测得的钢筋应变对比如图 6所示.对比发现,式(1) 的计算值明显偏小.上述试验试件的混凝土裂缝宽度达到0.2 mm时,将跨中开裂截面的钢筋、钢顶板、U肋(或球扁钢)下缘应变按照截面高度位置绘于图 7.可以看出,开裂截面应变已不再符合平截面假定,钢筋应变明显大于利用钢顶板、加劲肋底部应变按照平截面假定外推得到的应变.同时,对于常规混凝土试件,两者的差别较大,对于采用高强钢纤维混凝土的试件,两者的差别较小.可见,式(1) 中认为开裂截面仍然符合平截面假定,是导致计算结果偏小的原因.
基于上述试验结果和理论分析,本文提出了计算组合桥面板负弯矩区混凝土开裂截面钢筋应力的公式,其推导过程如下所示:
计算中假定:① 混凝土与钢顶板之间的剪力连接件完好,不会发生破坏,能够有效传递剪力;② 混凝土与钢筋之间的黏结完好,两者之间不发生滑移;③ 未开裂截面仍符合平截面假定.
取图 8的状态进行分析,此时,截面2位置的弯矩值刚好为开裂弯矩Mcr,跨中弯矩为M,在跨中截面1位置已经产生了一条贯通裂缝.取出截面2到截面1之间的组合桥面板段的混凝土板进行分析,其内力如图 9所示.
图 8中,Mcr=ftknscI0/yct,其中yct为混凝土板顶缘到组合截面形心的距离,nsc为钢与混凝土的弹性模量比,ftk为混凝土抗拉强度标准值,I0为换算截面惯性矩.
由图 9所示混凝土段x方向的合力为零可得
${{\sigma }_{\text{ss}}}{{A}_{\text{r}}}=\frac{1}{2}({{f}_{\text{tk}}}+{{f}_{\text{cb}}})\text{ }{{h}_{\text{c}}}b+b\int _{0}^{{{l}_{0}}}v\text{ }\left( x \right)\text{d}x+{{\sigma }_{\text{ss,cr}}}{{A}_{\text{r}}}$ | (2) |
式中:σss为开裂截面(截面1) 钢筋应力;Ar为钢筋截面积;fcb为开裂弯矩作用下混凝土截面下缘应力,fcb=Mcrycb/I0,其中ycb为混凝土下缘距组合截面中性轴的距离;hc为混凝土高度;b为混凝土板宽度;v (x)为混凝土与钢之间的剪应力;l0为截面1和2之间的距离;σss, cr为开裂弯矩作用下钢筋应力,σss, cr=Mcryr/I0.
对于承受跨中集中荷载的简支梁,截面1和2之间剪力分布
${{F}_{v}}\left( x \right)=\frac{\text{d}M\text{ }\left( x \right)}{\text{d}x}=\frac{M-{{M}_{\text{cr}}}}{{{l}_{0}}}=\frac{\Delta M}{{{l}_{0}}}$ | (3) |
由此可以得到从截面1到2钢与混凝土之间的剪应力分布
$v\text{ }\left( x \right)=\frac{{{F}_{\text{v}}}\left( x \right)\text{ }{{S}^{*}}}{{{I}_{0}}b}=\frac{\Delta M{{S}^{*}}}{{{I}_{0}}b{{l}_{0}}}$ | (4) |
其中,S*为混凝土板面积和钢筋面积对组合截面中性轴的静矩.
${{S}^{*}}=\frac{1}{{{n}_{\text{sc}}}}{{A}_{\text{c}}}{{y}_{\text{c}}}+{{A}_{\text{r}}}{{y}_{\text{r}}}$ | (5) |
其中,yc为混凝土板形心到组合截面形心的距离;Ac为混凝土板面积.
将式(4) 和(5) 代入式(2) 得
${{\sigma }_{\text{ss}}}=\frac{({{f}_{\text{tk}}}+{{f}_{\text{cb}}})\text{ }{{h}_{\text{c}}}b}{2{{A}_{r}}}+\frac{\Delta M\left( \frac{1}{{{n}_{\text{sc}}}}{{A}_{\text{c}}}{{y}_{\text{c}}}+{{A}_{\text{r}}}{{y}_{\text{r}}} \right)}{{{I}_{0}}{{A}_{\text{r}}}}+{{\sigma }_{\text{ss,cr}}}$ | (6) |
进行简化后得到钢筋应力计算公式为
${{\sigma }_{\text{ss}}}=\frac{({{f}_{\text{tk}}}+{{f}_{\text{cb}}})\text{ }{{h}_{\text{c}}}b}{2{{A}_{r}}}+\frac{\Delta M{{A}_{\text{c}}}{{y}_{\text{c}}}}{{{n}_{\text{sc}}}{{I}_{0}}{{A}_{\text{r}}}}+\frac{\Delta M{{y}_{\text{r}}}}{{{I}_{0}}}+{{\sigma }_{\text{ss,cr}}}$ | (7) |
式(7) 等号右边第一项和第四项的物理意义为混凝土截面开裂瞬间,原来混凝土所承担的力全部传递给钢筋承担,会使钢筋应力由σss, cr瞬间产生一个增量.此后随着荷载的增加,产生的钢筋应力增量由第二项和第三项计算得到.
此外,特别需要说明的是,即使在截面1(跨中开裂截面)和2(弯矩值等于开裂弯矩的截面)之间还有其他裂缝,也不会影响此公式计算的精度,因为只要混凝土和下部钢结构的剪力连接完好,混凝土和钢筋之间的黏结完好,就能保证剪力的顺畅传递,式(2) 就能成立,由此计算得到的钢筋应力就具有足够的精度.
将由式(7) 计算得到的钢筋应变和由钢桥规范推荐的式(1) 计算得到的钢筋应变与实测值进行比较并绘于图 10.
从图 10可以看出,由本文提出的式(7) 计算得到的钢筋应变与实测值吻合较好,钢桥规范中推荐的式(1) 计算出的结果明显小于实测值,而且变化趋势也与实测结果不符.可见平截面假定在开裂截面并不成立,基于此假定计算的钢筋应力、应变也会产生较大的偏差.
文献[12-15]中通过将红墨水或树脂注入裂缝观察混凝土开裂时的裂缝分布,得到了很多重要的结论,其中对混凝土开裂的特征作如下描述:① 裂缝表面是一个规则的曲面.混凝土裂缝宽度沿截面发生显著变化,在钢筋周界处的宽度最小,构件表面的混凝土裂缝宽度最大,两者相差3~7倍.② 钢筋周界处混凝土裂缝宽度很小,表明钢筋和混凝土的相对滑移小,即使是光圆钢筋,相对滑移也很小.③ 构件的受拉裂缝,除了表面垂直于钢筋轴线或间距和宽度都较大的主裂缝外,还有自钢筋表面横肋处向外延伸的内部裂缝(即次裂缝),如图 11所示.这些裂缝首先在裂缝截面附近产生,次裂缝数量较多,间距小,从钢筋向外延伸,但未到达构件表面[16].可以看出,在裂缝截面附近,混凝土与钢筋虽然黏结完好,但由于次裂缝的存在两者的应变并不协调,钢筋应变要比周围混凝土应变大.这个现象进一步验证了上述推导的假定.
组合桥面板混凝土裂缝宽度与许多因素有关.除了上面提到的开裂截面钢筋应力外,还与组合结构力比、横向钢筋间距、配筋率、剪力钉间距、保护层厚度等诸多因素有关,以公预规和混凝土规范为代表的两类混凝土裂缝宽度计算公式分别不同程度地考虑了上述因素的影响.本文采用四种方法计算了试验中试件的混凝土裂缝宽度,并将其与试验值进行对比.
方法一,采用式(7) 计算得到的钢筋应力,按照公预规中的混凝土裂缝宽度计算公式,采用偏心受拉构件的系数计算,计算公式如下:
${{w}_{\text{tk}}}={{C}_{1}}{{C}_{2}}{{C}_{3}}\frac{{{\sigma }_{\text{ss}}}}{{{E}_{\text{s}}}}\left( \frac{30+d}{0.28+10\rho } \right)$ | (8) |
式中:C1为钢筋表面形状系数,对带肋钢筋取1.0;C2为作用长期效应影响系数,由于采用试验试件,不考长期效应的影响,取1.5;C3为与构件受力性质有关的系数,这里取偏拉构件系数1.1;Es钢筋弹性模量;d钢筋直径;ρ为混凝土板配筋率.
由于组合桥面板在混凝土板下层还有正交异性钢板通过剪力连接件与混凝土连接,当混凝土开裂后,正交异性钢板的钢顶板也承受拉应力,和普通钢筋一样具有限制裂缝开展的作用.因此,在计算配筋率时也应该考虑这部分钢板的作用,经过笔者试算,发现不论普通钢筋如何配置,在考虑了钢顶板后,配筋率均会超过公预规公式的上限值0.02,所以这里5个试件的配筋率均取为0.02.
方法二,采用钢桥规范推荐的方法,即采用式(1) 计算得到的钢筋应力σss,按照公预规中的混凝土裂缝宽度计算公式,用轴心受拉构件的系数C3=1.2计算,其他系数取值同方法一.
方法三,采用文献[6]提出的组合梁负弯矩区混凝土裂缝宽度计算公式计算,计算公式如下:
${{w}_{\text{m}}}=1.45\psi \prime {{\varepsilon }_{\text{r}}}{{l}_{\text{cr}}}$ | (9) |
式中:
${{l}_{\text{cr}}}=1.1\left( 2.7c+\frac{0.11}{{{\rho }_{\text{ct}}}/d+0.25{{R}^{2}}/p} \right)\nu $ | (10) |
式中:c为最外层受拉区钢筋的保护层厚度;p为栓钉连接件间距;ν为与钢筋有关的系数,带肋钢筋为0.7,光圆钢筋为1.0.
方法四,采用文献[8]提出的组合梁负弯矩区混凝土裂缝宽度计算公式计算,计算公式如下:
${{w}_{\text{max}}}={{\alpha }_{\text{cr}}}\varphi ({{\sigma }_{\text{rk}}}/{{E}_{\text{r}}})\text{ }{{l}_{\text{cr}}}$ | (11) |
$\varphi =1.1-0.65\frac{{{f}_{\text{tk}}}\sqrt{{{R}_{\text{p}}}}}{{{\rho }_{\text{te}}}{{\sigma }_{\text{rk}}}}$ | (12) |
${{l}_{\text{cr}}}={{l}_{a}}(1-R_{p}^{3})$ | (13) |
式中:ρte为混凝土板配筋率; σrk为钢筋应力,采用式(1) 计算; Er为钢筋弹性模量;lcr为裂缝间距;φ为钢筋应力不均匀系数;Rp为考虑了预应力钢筋后的综合力比,对于组合桥面板,计算公式同方法三;αcr为构件的受力特征系数,取1.0;la为混凝土板横向钢筋间距.
需要说明的是,对于组合桥面板,由于混凝土板厚度较小,只配置一层钢筋,中性轴一般在混凝土和钢交界的位置,负弯矩区开裂时混凝土全截面承受拉应力,因此在计算截面配筋率时,不宜按照公预规中规定将轴心受拉构件受拉钢筋面积取一半,而应取所有受拉钢筋面积计算配筋率.
上述四种方法的计算结果与实测值的对比如图 12所示.
由图 12可以看出,对于试件S1,方法一的计算结果与实测值吻合较好,方法二、方法三和方法四的计算结果偏小.对于试件S2,方法一、方法三和方法四计算结果均与实测值吻合较好,其中方法一在混凝土裂缝宽度小于0.3 mm时计算结果较为准确,混凝土裂缝宽度大于0.3 mm后计算结果偏大;试件S3采用了钢纤维混凝土,方法一计算结果趋近实测值的上限,计算结果较为保守,方法二和方法四与实测值吻合较好,方法三结果偏小.对于试件S4和S5,方法一计算结果与实测值较为吻合,方法二、三和四的结果均偏小.总体而言方法一的计算结果与试验结果的吻合情况最好.
由上述分析可以看出,对于组合桥面板,在进行负弯矩区混凝土裂缝宽度验算时,采用本文提出的方法计算钢筋应力,并使用公预规提出的混凝土裂缝宽度计算公式按照偏心受拉构件计算,可以得到与实测值较为吻合的结果,且概念清晰,推荐在进行桥面板负弯矩区混凝土裂缝宽度验算时使用.对于采用钢纤维混凝土的组合桥面板,此方法计算结果略微保守.
3 结论(1) 根据5个试件的试验结果,可以观察到如下试验现象:混凝土标号、钢筋直径、加劲肋类型对负弯矩区裂缝的扩展速率有较为显著的影响,而剪力连接件形式的影响不大.其中,采用钢纤维混凝土的试件开裂荷载较高,采用细直径钢筋的试件裂缝扩展较快,采用球扁钢的试件在破坏时裂缝的分布范围较采用U肋的试件更大.由于试件数量有限,上述试验现象仍然需要大量试验试件进行验证.
(2) 开裂截面钢筋和正交异性钢板的应变分布并不符合平截面假定,实际钢筋的应变明显偏大,常规混凝土钢筋应力实测值与按平截面假定计算得到的结果差别较大,高强钢纤维混凝土两者的差别则较小.
(3) 推导了组合桥面板负弯矩区开裂截面钢筋应力的计算公式,计算结果与实测值吻合较好.
(4) 采用四种方法计算了5个试件的混凝土裂缝宽度,发现采用本文推荐钢筋应力公式,并使用公预规提出的混凝土裂缝宽度计算公式按照偏心受拉构件计算,可以得到与实测值吻合较好的结果.
[1] |
ZHANG Qinghua, CUI Chuang, BU Yizhi, et al. Fatigue tests and fatigue assessment approaches for rib-to-diaphragm in steel orthotropic decks[J]. Journal of Constructional Steel Research, 2015, 114: 110 DOI:10.1016/j.jcsr.2015.07.014 |
[2] |
XU W, ZHANG X N. Experimental and numerical simulation study on the crack of steel orthotropic bridge deck pavement[C]//6th RILEM International Conference on Cracking in Pavements. Chicago: CRC Press-Taylor & Francis Group, 2008.
|
[3] |
ALLAHYARI H, DEHESTANI M, BEYGI M H A, et al. Mechanical behavior of steel-concrete composite decks with perfobond shear connectors[J]. Steel and Composite Structures, 2014, 17(3): 339 DOI:10.12989/scs.2014.17.3.339 |
[4] |
中华人民共和国住房和城乡建设部. 混凝土结构设计规范: GB 50010—2010[S]. 北京: 中国建筑工业出版社, 2010. Ministry of Housing and Urban-Rural Development of the People's Republic of China. Code for design of concrete structures:GB 50010—2010[S]. Beijing: China Architecture & Building Press, 2010. |
[5] |
中华人民共和国交通运输部. 公路钢筋混凝土及预应力混凝土桥涵设计规范: JTG D62—2004[S]. 北京: 人民交通出版社, 2004. Ministry of Transport of the People's Republic of China. Code for design of highway reinforced concrete and prestressed concrete bridges and culverts: JTG D62—2004[S]. Beijing: China Communications Press, 2004. |
[6] |
聂建国, 张眉河. 钢混凝土组合梁负弯矩区板裂缝的研究[J]. 清华大学学报(自然科学版), 1997, 37(6): 97 NIE Jianguo, ZHANG Meihe. Study on the crack of concrete flange in tension of composite steel-concrete beams[J]. Journal of Tsinghua University(Science and Technology), 1997, 37(6): 97 |
[7] |
余志武, 郭风琪. 部分预应力钢混凝土连续组合梁负弯矩区裂缝宽度试验研究[J]. 建筑结构学报, 2004, 25(4): 55 YU Zhiwu, GUO Fengqi. Experimental study of crack width in negative bending region of partially prestressed continuous steel-concrete composite beams[J]. Journal of Building Structures, 2004, 25(4): 55 |
[8] |
张彦玲, 樊健生, 李运生. 连续组合梁桥裂缝发展规律分析及裂缝宽度计算[J]. 工程力学, 2011, 28(7): 84 ZHANG Yanling, FAN Jiansheng, LI Yunsheng. Law of crack development and calculation of crack width of continuous composite beams[J]. Engineering Mechanics, 2011, 28(7): 84 |
[9] |
张彦玲. 钢混凝土组合梁负弯矩区受力性能及开裂控制的试验及理论研究[D]. 北京: 北京交通大学, 2009. ZHANG Yanling. Theoretical analysis and experimental research on behavior and crack control of negative moment zone in steel-concrete composite beams [D]. Beijing: Beijing Jiaotong University, 2009. |
[10] |
RYU H-K, CHANG S-P, KIM Y-J, et al. Crack control of a steel and concrete composite plate girder with prefabricated slabs under hogging moments[J]. Engineering Structures, 2005, 27(11): 1613 DOI:10.1016/j.engstruct.2005.05.015 |
[11] |
中华人民共和国交通运输部. 公路钢结构桥梁设计规范: JTG D64—2015[S]. 北京: 人民交通出版社, 2015. Ministry of Transport of the People's Republic of China. Specifications for design of highway steel bridge: JTG D64—2015[S]. Beijing: China Communications Press, 2015. |
[12] |
BROMS B B. Crack width and crack spacing in reinforced concrete members[J]. ACI Structural Journal, 1965, 62(10): 1237 |
[13] |
BROMS B B, LUTZ L A. Effects of arrangement of reinforcement on crack width and spacing of reinforced concrete members[J]. ACI Journal Proceedings, 1965, 62(11): 1395 |
[14] |
BASE G D. Crack control in reinforced concrete-present position[C]//Syposium on Serviceability of Concrete. Melbourne: Institution of Engineers, 1975:1-10.
|
[15] |
REED J B, BASE G D, BEEBY A W, et al. An investigation of the crack control characteristics of various types of bar in reinforce concrete beams [R]. London: Cement and Concrete Association, 1966.
|
[16] |
丁嵬. 钢筋混凝土结构裂缝宽度计算方法研究[D]. 天津: 天津大学, 2007. DING Wei. The study on crack width calculation methods of reinforced concrete structure[D]. Tianjing: Tianjing University, 2007. |