网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

超大断面隧道变截面衬砌台车应力特征  PDF

  • 王浩 1,2
  • 杨新安 1,2
  • 罗驰 1,2
  • 刘臻武 3
1. 同济大学 道路与交通工程教育部重点实验室,上海 201804; 2. 同济大学 交通运输工程学院,上海 201804; 3. 中铁五局集团第四工程有限责任公司,广东 韶关 512031

中图分类号: U25

最近更新:2021-03-17

DOI:10.11908/j.issn.0253-374x.20232

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
目录contents

摘要

以京张高铁八达岭长城站为工程依托,对超大断面隧道变截面衬砌台车的动态、静态应力进行研究,分析了截面变化对台车不同结构应力的影响规律。通过数值模拟分析台车结构的静态应力特征,利用现场监测分析台车结构的动态应力特征以及变截面对台车结构应力变化的影响规律。结果表明:台车结构静态应力和位移存在不均匀分布现象,各结构呈多变特征,部分结构出现明显的应力集中;浇筑过程中,应力呈现外侧大内部小的规律;台车浇筑区域位于支撑结构上方时应力最大,混凝土凝固过程中台车结构应力基本不变;台车截面由小至大变化时,各结构应力增大,其中拱肩丝杆是应力变化最大的构件。

随着铁路隧道复合式衬砌的普及,二次衬砌施工已成为隧道建设的重要环节,其施工方式直接影响了隧道的结构稳定和外观质量。从20世纪80年代衬砌台车进入中国以来,已成为衬砌施工的主要方式,从而有效保证了衬砌的施工质量。然而,对于超大断面变截面隧道,若每个断面单独制作台车则造成大量的资源浪费,因此变截面衬砌台车成为这类工程的首选。

随着衬砌台车的大规模应用,在应力特征方面取得了一些研究成果。目前,理论分析多是将台车拆解成单独的杆或板构件进行计算。曹小

1、苏睿2指出,台车底部约束框架承受荷载最大,该位置的约束反力是影响台车稳定的主要因素。蔡明3认为,台车门架受力最大,是控制台车位移的关键结构。超大断面隧道衬砌台车结构复杂,先拆解再单独计算的方式由于计算量大、精度不足等原因并不适用于超大断面隧道衬砌台车,而数值模拟可以实现台车的整体计算。姬海东4通过三维数值模拟分析指出,台车应力最大位置在纵向连接梁的中部。徐爱英5通过有限元计算,得到门架立柱是台车应力最大构件的结论。郝良果6对台车模板的受力进行了分析,认为台车模板应力最大位置在拱顶处。张凤7通过Ansys分析表明,边墙模板是台车位移的最大部位。以上模拟中,多是先建模后单次施加全部荷载进行计算,因此未能反映台车在工作过程中的动态应力变化。除此之外,在工程实践方面,但淑英8通过对石牙山隧道衬砌台车分析后认为,门架横梁是台车应力最大的结构;陈9对胶州湾隧道衬砌台车的有限元分析结果表明,门架是台车应力最大结构。

衬砌台车存在应力集中现象,但目前的研究成果多是对普通断面隧道衬砌台车静态应力特征进行分析,对于超大断面隧道衬砌台车的静态应力特征和工作时的动态应力特征的研究则很少。依托京张高铁八达岭长城站大跨过渡段衬砌施工,基于数值模拟的方式,对超大断面隧道衬砌台车的静态应力特征进行分析。使用振弦式钢筋应力计对台车在工作过程中的应力进行动态监测,探究台车工作过程中各结构应力的变化规律和关键值,并对比截面变化时台车各结构的应力变化。

1 工程概况

京张高铁是国家规划实施的重点建设项目,也是2022年北京张家口联合举办冬季奥运会的配套工程。京张高铁八达岭长城站位于八达岭长城景区下方新八达岭隧道之中,是目前国内埋深最大的暗挖高铁地下车站。八达岭长城站的大跨过渡段是连接车站与隧道正线之间的过渡结构,张家口方向过渡段里程为DK68+285―DK68+448,全长163 m,埋深90 m,过渡段设计为超大断面、变截面型式,由车站向隧道正线分5个断面进行过渡,如图1所示。北京方向由正线向车站过渡也采用同样的对称设计。

图1 大跨过渡段平面图

Fig.1 Plan of large-span transition section

过渡段各断面轮廓如图2所示。断面1为最小断面,断面5为最大断面,最大开挖断面面积为494.4 m2。围岩等级主要为Ⅲ、Ⅳ级,局部为Ⅴ级,是目前世界上开挖跨度最大、开挖断面面积最大的交通隧道。大跨过渡段二次衬砌施工是施工难点,国内外尚无相似工程参考。

图2 大跨过渡段断面图

Fig.2 Cross-sectional view of large-span transition section

大跨过渡段二次衬砌采用长寿命混凝土浇筑,断面5衬砌厚度为60 cm,断面4和断面3衬砌厚度为55 cm,断面2和断面1衬砌厚度为45 cm。过渡段断面变化5次,衬砌厚度变化3次,若每个断面单独制作台车则造成严重的资源浪费,因此衬砌台车设计要满足断面可变、多变的要求。

2 衬砌台车结构设计

衬砌台车的组成系统主要包括门架系统、支撑系统、模板系统、行走系统、电气化操作系统等。变截面方式主要包括改变门架系统、调整支撑系统结构、增减和旋转模板等。

(1)门架系统

台车门架系统是主要承重系统,由四榀门架组成,如图3所示。每榀门架包括5根立柱、4根短横梁、1根长横梁,各榀之间间距为3 m,分别在门架的上、中、下位置布置纵梁将各榀门架连接在一起。立柱、斜拉杆、短横梁和纵梁设计为工字钢,长横梁设计为厚度20 mm与16 mm的钢板围成的箱型梁,高度为1 200 mm,宽度为384 mm。根据各断面跨度在长横梁对应位置布置定位孔,台车跨度变化时,两侧立柱移动,通过定位孔定位。

图3 门架系统结构

Fig.3 Structure of gantry system

(2)支撑系统

支撑系统将门架系统和模板系统连在一起,将模板系统荷载传递至门架系统,如图4图5所示。主要构件为28号可拆卸工字钢和可拉伸丝杆,丝杆由Φ114×12无缝钢管与Φ80圆钢焊接。边墙支撑设计为可拉伸丝杆,与门架立柱相铰接,拱部支撑下层为可拆卸工字钢,上层为可拉伸丝杆。支撑系统共七榀,每榀之间的间距为1.5 m,各榀之间通过工字钢或槽钢刚接。

图4 边墙支撑结构

Fig.4 Structure of sidewall support

图5 拱部支撑结构

Fig.5 Structure of arch support

(3)模板系统

模板系统决定了衬砌混凝土表面完好程度,由多块小模板组成,纵向用螺栓连接成需要的长度,环向采用铰接连接,模板可以绕着铰接转动,通过伸缩丝杆、增加小模板将台车调整至隧道设计轮廓线。

大跨过渡段各断面台车门架结构如图6所示。台车使用Q235钢,最大跨度为28.95 m,最大高度为13.71 m,纵向长度为9.00 m。

图6 大跨过渡段衬砌台车结构

Fig.6 Structure of lining trolley for long-span transition section

3 台车静态应力特征分析

图6各断面台车结构的对比可以看出,断面5是截面面积最大的断面,也是结构最复杂的断面,因此断面5衬砌台车的静态应力特征分析对整个过渡段衬砌台车受力特征的研究有着重要意义。利用有限元软件MIDAS GTS‒NX建立断面5衬砌台车三维数值模型,分析各组成系统的轴力、轴应力和位移分布特征以及最大值所在位置,为现场监测提供参考。

衬砌台车的模拟过程分为建模、输入荷载、结构分析和输出结果等4个步骤,即首先定义台车模型材料和各结构截面,建立台车数值模型,然后根据混凝土的浇筑状态确定台车模板承受的荷载,运行结构分析,最终输出分析结果。

3.1 数值模型建立

台车模型如图7所示。门架系统和支撑系统选择梁单元和板单元,模板系统选择板单元。各结构模型输入参数如表1所示。台车竖向由门架立柱与油缸共同支撑,中心立柱与油缸放置于钢轨上,故仅施加竖向约束,不限制中心立柱的左右转动;剩余立柱为固定端,故施加三向约束,同时限制剩余立柱在各方向的转动。

图7 断面5台车数值模型

Fig.7 Numerical model of section 5 trolley

表1 数值模型输入参数
Tab.1 Numerical model input parameters
序号台车结构结构型式规格材质容重/(kN∙m-3
1 门架立柱 工字钢 320 mm×260 mm×19 mm Q235 78.5
2 门架短横梁 工字钢 320 mm×260 mm×19 mm Q235 78.5
3 门架长横梁 箱型梁 1 200 mm×384 mm×16 mm×20 mm Q235 78.5
4 门架斜拉杆 工字钢 280 mm×280 mm×15 mm Q235 78.5
5 可拉伸丝杆 圆钢管 114 mm×12 mm Q235 78.5
6 可拆卸工字钢 工字钢 280 mm×280 mm×15 mm Q235 78.5
7 门架纵梁 工字钢 280 mm×280 mm×15 mm Q235 78.5
8 支撑纵梁 槽钢 125.0 mm×53.0 mm×5.5 mm Q235 78.5
9 模板 曲面钢板 69.25 mm Q235 78.5

3.2 输入荷载确定

拱部模板与门架长横梁连接成一个穹形钢体,边墙模板与门架立柱连接成一个穹形钢体,因此只对模板施加荷载。模板的荷载方向主要沿着拱形构件的径向,用垂直方向荷载与水平方向荷载进行等效处理,同时对整个模型结构施加重力。

根据GB50068―2018《建筑结构可靠性设计统一标准》,台车荷载由永久荷载Pd和可变荷载Pv组成,根据断面5衬砌厚度和混凝土密度,垂直方向荷载P1的永久荷载Pd1=25 kN∙m-2,可变荷载Pv1=4.5 kN∙m-2。水平方向荷载P2的永久荷载Pd2根据《建筑手册》中“现浇混凝土结构模板的设计”侧压力公式计算,如下所示:

Pd2=0.22ρ0t0β1β2v12 (1)

式中:Pd2为新浇筑混凝土对模板的最大侧压力,kN∙m-2ρ0为混凝土密度,取25 kN∙m-3t0为新浇筑混凝土的初时间,取5 h;β1为外加剂影响修正参数,取1.2;β2为混凝土坍落度影响修正参数,取1.15;v为混凝土的浇筑速度,取2.5 m∙h-1。水平方向荷载的永久荷载Pd2=60 kN∙m-2

水平方向的可变荷载Pv2=6.0 kN∙m-2。根据荷载组合公式,水平方向荷载P2的计算式为

P2=1.3Pd2+1.5Pv2 (2)

计算得水平方向荷载P2=87 kN∙m-2

3.3 模拟结果分析

台车轴力模拟结果如图8所示。门架系统和支撑系统以承受压力为主,门架系统和边墙丝杆的轴力大于拱部支撑系统,压力最大值为1.15×106 N,位于门架中立柱,拉力最大值为1.75×105 N,位于边墙弧板。

图8 台车轴力

Fig.8 Axial force of the trolley

台车轴应力模拟结果图9所示。拱部支撑工字钢、门架横梁和两侧门架立柱轴应力约为-1.70 MPa,门架中立柱和拱顶丝杆轴应力为-15.06 MPa,拱部两侧丝杆轴应力为-28.43 MPa,边墙丝杆轴应力为-41.79 MPa。通过不同系统对比可以看出,与模板相连丝杆的轴应力最大,其次为拱部工字钢,门架系统最小。

图9 台车轴应力

Fig.9 Axial stress of the trolley

台车位移模拟结果如图10所示。门架系统位移较小,位移最大值约为0.57 mm,拱部工字钢位移为1.14~1.70 mm,拱顶模板位移为3.42 mm,边墙模板位移为6.84 mm。通过不同位置的位移对比可以看出,模板系统是整个台车位移最大的部分,其中边墙模板处位移最大。

图10 台车位移

Fig.10 Displacement of the trolley

从台车的组成可以看出,与门架系统相比,支撑结构的丝杆和工字钢间排距较小,因此当荷载确定时,支撑结构的轴力小于门架系统,位移大于门架系统。与其他大断面隧道衬砌台车相

10,大跨过渡段断面5衬砌台车的门架系统由于增加了中立柱,轴力向中间转移,因此门架中立柱成为轴力最大结构,轴力值约为门架其他立柱的2倍。由于边墙模板以水平荷载为主,拱部模板以竖向荷载为主,水平荷载大于竖向荷载,因此边墙丝杆的轴力和轴应力大于拱部丝杆。

4 台车动态应力特征分析

从二次衬砌开始浇筑至凝固拆模,台车各结构应力处于动态变化过程中,对台车结构的应力实施动态监测,不仅可以掌握应力在整个过程中的变化规律和关键值,还可以验证模拟计算的准确性。为此,选择断面5的衬砌台车作为分析对象,构建监测方案并进行结果分析。

4.1 台车应力动态监测方案

(1)测点布置

由数值模拟结果可知,边墙丝杆是轴应力最大位置,其次为拱部丝杆,台车门架中立柱是轴力最大位置,因此分别在以上位置布置监测点,记录台车在浇筑过程中的应力值。由于浇筑过程中台车对称受力,因此只在单侧布置监测点,共布置6个测点,各测点位置如图11所示。

图11 监测点布置

Fig.11 Layout of monitoring points

(2)仪器选择

为掌握应力变化规律和关键数据,必须对浇筑和凝固过程的应力值进行长时间、短间隔的采集,选用自动化数据采集仪器成为必

11-14。测量仪器选择振弦式钢筋应力计,采集仪器选择DT80G采集器和CEM20扩展模块,数据采集频率为15 s。

4.2 应力动态监测结果分析

断面5衬砌自下向上分4层浇筑,各层浇筑时间、浇筑高度和混凝土凝固时间如表2所示,浇筑及凝固过程应力‒时间关系如图12所示。

表2 断面5监测记录
Tab.2 Section 5 monitoring record
进度时间/h高度/m
第1层浇筑 2.80 6.38
第2层浇筑 3.70 4.01
第3层浇筑 4.10 3.08
第4层浇筑 5.80 0.89
拱顶压浆 4.55
衬砌凝固 101.05

图12 断面5台车应力‒时间曲线

Fig.12 Stress-time curve of section 5 trolley

(1)边墙丝杆应力变化分析

1号点应力值在第1层、第2层浇筑时减小至 -7.90 MPa,第3层浇筑时增大至-1.20 MPa,第4层浇筑、拱顶压浆时应力值减小至-4.00 MPa,凝固过程基本不变。2号点应力值在第3层浇筑时减小至-25.01 MPa,第4层浇筑时应力值增大至 -14.30 MPa,拱顶压浆、混凝土凝固时应力基本不变。1、2号点的应力值经历了前2层浇筑而减小,此后增大再稳定的过程,说明第1、2层浇筑对边墙丝杆的影响大于其他阶段。

(2)拱部丝杆应力变化分析

3号点应力值在前3层浇筑时减小至-24.18 MPa,第4层浇筑期时,应力值增大至-12.50 MPa,拱顶压浆时减小至-16.01 MPa。4号点应力值前3层浇筑时减小至-13.20 MPa,此后应力值不变,拱顶压浆时应力值增大至-8.70 MPa。5号点在第4层浇筑时压力值开始减小,压浆结束时应力值为 -11.10 MPa。从断面5拱部丝杆的压力显现顺序看出,浇筑先影响拱腰丝杆,后影响拱顶丝杆,混凝土凝固期间应力值基本不变。

(3)门架立柱应力变化分析

6号点应力值在第4层浇筑时突变至-4.83 MPa,拱顶压浆时减小至-17.88 MPa,混凝土凝固过程中,无明显变化,基本与浇筑完成时的应力值一致。从6号点的变化规律可以看出,门架中立柱受力主要受拱顶压浆的影响。

从应力变化曲线可以看出,衬砌开始浇筑后,边墙丝杆应力持续增大。随着浇筑区域向台车中上部转移,边墙丝杆的应力达到峰值后开始减小或小幅变化,拱部两侧丝杆开始受力并持续增大。拱顶浇筑时,拱部丝杆应力达到峰值后开始减小或小幅变化,拱顶丝杆和门架中立柱应力持续增大。拱顶压浆时拱顶丝杆和门架中立柱应力持续增大,其他位置应力持续减小或小幅变化,监测结果反映了台车在工作过程中的应力变化。

从数值上来看,拱肩2号点和3号点的应力监测最大值分别为-25.01 MPa和-24.18 MPa,该处应力模拟值约为-28.43 MPa,拱顶4号点、5号点的应力监测最大值分别为-13.20 MPa和-11.10 MPa,模拟值为-15.06 MPa,均较为吻合。从模拟计算的台车轴应力图来看,拱肩位置丝杆轴向压应力较大,拱顶其次,而拱脚位置压应力最小,甚至可能出现拉应力。从实测的丝杆轴应力可以看出,拱肩2号点、3号点轴应力最大,拱顶处的4号点、5号点轴应力其次,而1号点轴应力最小,甚至接近于零,这与数值模拟结果基本一致。

从台车断面的轮廓线和结构可以看出,台车模板的拱顶和边墙曲率较大,拱顶以竖向荷载为主,边墙以水平荷载为主,拱肩处曲率存在突变,在拱肩曲率突变处会产生较大的应力集中荷载,所以2号点和3号点的应力监测值最大,其次为4号点、5号点。

5 变截面对应力变化影响分析

过渡段由断面1向断面5分5个断面过渡,为体现截面变化对台车应力的影响,基于台车应力的动态监测,对不同截面台车的边墙丝杆、拱肩丝杆、拱顶丝杆和门架立柱应力峰值进行对比,分析变截面对衬砌台车应力的影响规律及特征。

台车不同结构峰值应力随截面的变化如图13所示。由于边墙模板的曲率变化较小,故边墙丝杆的应力变化最小,断面1至断面5变化时仅增加了2.60 MPa;台车由断面1变化至断面3时,拱顶曲率逐渐增大,拱肩曲率发生突变,因此拱肩丝杆的应力增量大于拱顶丝杆,拱肩丝杆应力增加了11.10 MPa,拱顶丝杆应力增加了3.20 MPa,断面4向断面5变化时,拱顶丝杆和拱肩丝杆的应力变化较小;台车由断面1向断面2变化时,门架系统增加了中立柱,荷载向中立柱转移,应力增加了3.60 MPa,其他截面变化时应力变化较小。

图13 台车结构峰值应力随截面的变化

Fig.13 Change of peak stress of trolley structure with section

综合边墙丝杆、拱肩丝杆、拱顶丝杆和门架立柱应力峰值的变化趋势可以看出,随着台车跨度和高度的增加,各结构的应力也相应增大。台车不同位置模板的曲率变化是影响台车受力变化的重要因素,曲率变化越大,应力变化越大。

6 结论

(1)变截面衬砌台车由门架系统、支撑系统和模板系统等组成,通过伸缩门架横梁和边墙丝杆实现跨度变化,增加拱部可拆卸工字钢实现高度变化,增加旋转模板改变截面轮廓。

(2)台车各系统的应力和位移呈多变特征,门架系统轴力大于其他系统,拱肩丝杆是轴应力最大位置,模板系统位移大于门架系统和支撑系统,边墙模板是位移最大位置。

(3)应力在拱部浇筑后出现峰值,边墙丝杆在拱部两侧浇筑时应力达到-25.01 MPa,拱肩丝杆在拱顶浇筑时应力达到-24.18 MPa,拱顶丝杆和门架中立柱在拱顶压浆时应力达到-11.10、-17.88 MPa,混凝土凝固过程中应力与浇筑结束时应力基本一致。

(4)浇筑过程中应力显现顺序依次为边墙、拱部两侧、拱顶、门架中立柱,总体呈“先边后中”、“先低后高”的趋势,台车外部丝杆应力大于内部门架系统,两侧丝杆大于中间丝杆,总体呈“外大内小”、“边大中小”的规律。

(5)台车前4个截面变化时应力变化较大,截面面积和模板曲率是影响各结构应力变化的重要因素。截面增大时,各结构应力也相应增大,模板曲率变化越大,应力变化越大。

作者贡献声明

王 浩:参与台车结构设计,建立台车数值模型,对数值模型进行综合分析,设计台车监测方案,对台车进行现场监测和数据分析,完成文章初稿,进行稿件修改及专家意见反馈工作。

杨新安:指导数值模型分析和监测方案,指导后续分析工作,指导文章撰写及修改工作,项目管理。

罗 驰:参与数值模拟分析和台车现场监测,参与文章初稿。

刘臻武:参与台车结构设计,参与台车现场监测。

参考文献

1

曹小林.隧道混凝土衬砌施工中衬砌台车受力分析[J].铁道工程学报20023):51. [百度学术

CAO Xiaolin. Loading analysis of lining jumbo during construction of concrete tunnel lining[J].Journal of Railway Engineering Society20023):51. [百度学术

2

苏睿李生宏孟祥义.大模板液压台车受力及行走稳定检算[J].铁道建筑技术20025):34. [百度学术

SU RuiLI ShenghongMENG Xiangyiet al. Checking of force and traveling stability of large formwork hydraulic trolley[J].Railway Construction Technoligy20025):34. [百度学术

3

蔡明慧.衬砌台车关键受力部位的确定及强度计算[J].石家庄铁路职业技术学院学报20071):65. [百度学术

CAI Minghui. Finite element structure simulation and optimization design of tunnel lining trolley[J]. Journal of Shijiazhuang Institute of Railway Technology20071):65. [百度学术

4

姬海东刘在政张海涛.新型带压浇筑隧道数字化衬砌台车研究与应用[J].隧道建设2018388):1384. [百度学术

JI HaidongLIU ZaizhengZHANG Haitao. Research and application of a new type of digital tunnel lining trolley with pressured casting[J]. Tunnel Construction2018388): 1384. [百度学术

5

徐爱英牛旭赵增耀.隧道衬砌台车有限元结构仿真及优化设计[J].塔里木大学学报2009214): 26. [百度学术

XU AiyingNIU XuZHAO Zengyao. Trolley tunnel lining structure of the finite element simulation and optimal design[J]. Journal of Tarim University2009214):26. [百度学术

6

郝良果陈守强康博.隧道衬砌台车模板的有限元分析[J].山西建筑20093531):323. [百度学术

HAO LiangguoCHEN ShouqiangKANG Bo. Finite analysis of the template of the trolley tunnel[J]. Shanxi Architecture20093531):323. [百度学术

7

张凤雷. 湘桂高铁刘家冲隧道模板台车的结构分析与振动器选择[D]. 成都西南交通大学2012. [百度学术

ZHANG Fenglei. The structure analysis and choice of vibrators for model board trolley used in Xiang-Gui High-Speed Rail Liujiachong Tunnel[D]. ChengduSouthwest Jiaotong University2012. [百度学术

8

但淑英刘涛.隧道双功能全液压伸缩全自动行走混凝土衬砌台车结构受力分析[J].公路20097):281. [百度学术

DAN ShuyingLIU Tao. Structural force analysis of tunnel dual-function full-hydraulic telescopic full-automatic walking concrete lining trolley[J]. Highway20097):281. [百度学术

9

陈鹰. 青岛胶州湾隧道超大断面变截面隧道二次衬砌台车结构受力分析[J].现代隧道技术2012增刊):249. [百度学术

CHEN Ying. Mechanical analysis of the secondary lining trolley structure of the Jiaozhou Bay Tunnel in Qingdao[J]. Modern Tunnelling Technology2012(Z):249. [百度学术

10

唐果良.特大断面隧道液压模板台车的研制与施工技术[J].现代隧道技术20064):77. [百度学术

TANG Guoliang. Development and fabrication of a hydraulic-drive tunnel shutter jumbo with super large cross-section[J]. Modern Tunnelling Technology20064):77. [百度学术

11

YIN HLEFTICARIU LWEI Jet al. In situ dynamic monitoring of stress revolution with time and space under coal seam floor during longwall mining[J]. Environmental Earth Sciences20167518):1249. [百度学术

12

ZHOU YangSHI JunfengYU Tianqiu. Steel truss girder bridge reinforcement removing phase analysis of stress dynamic monitoring[J]. Advanced Materials Research2013668650. [百度学术

13

HE JianDOU LinmingGONG Siyuanet al. Rock burst assessment and prediction by dynamic and static stress analysis based on micro-seismic monitoring[J]. International Journal of Rock Mechanics & Mining Sciences20179346. [百度学术

14

XU CIMAI YWANG Wet al. Real-time monitoring of dynamic stress concentration by mechanoluminescent sensing film[J]. Applied Mechanics & Materials200813247. [百度学术