1 引 言

当今世界，随着人类对清洁能源需求量的增加和防汛抗洪及灌溉用水的需要，越来越多的混凝土大坝正在建设或者已经建设完成。虽然和平与发展已经成为世界的主题，但是，局部冲突仍然存在，恐怖袭击时有发生。由于蓄水大坝具有很强的政治、军事和经济意义，极有可能成为现代战争中的潜在打击目标或恐怖袭击的重要对象。故对混凝土坝抗爆性能以及抗爆安全性的研究是工程设计中至关重要的一环[1]。

2)相比于还需要借助眼镜来实现的VR技术的AR，已经可以摆脱智能眼镜的捆绑，因此也获得了更大的发展空间，得到了很多科技企业和景区的青睐。

研究混凝土坝的抗爆安全性，最重要的是研究其在爆炸荷载下的动态响应。由于爆炸发生时的物理反应过程极其复杂，所以在爆炸荷载作用下混凝土坝的动态响应过程比由其他任何荷载所引起的动态响应过程都要复杂得多。近年来，国内外许多学者对混凝土坝的抗爆安全性做了一定的研究。李本平等[2]在综合考虑射弹的冲击侵彻作用、破碎弹壳爆炸能损耗和侵彻孔洞对爆轰产物的影响等因素的基础上，通过对制导炸弹水平侵彻混凝土重力坝爆炸全过程的三维数值仿真，研究了混凝土重力坝在射弹侵彻爆炸下的动态响应及破坏效应；张社荣等[3]对爆炸荷载下混凝土重力坝的抗爆性能和破坏模式进行了数值模拟，讨论了遭受水下爆炸情况下混凝土重力坝的可能失效模式，并且研究了大坝高度、爆炸距离以及上游水位对大坝抗爆性能的影响。徐俊祥等[4]建立了全耦合的模型，考虑了炸药的爆炸、爆炸波的水中传播以及水介质与坝体结构的相互作用，对混凝土重力坝的动力响应进行了模拟，得出了坝体的加速度响应以及位移的时程变化特点。王冰玲等[5]开发了一个适合于脆性材料破坏模拟的随机网格生成算法，利用该方法将介质用四面体单元离散化，并选择离散元法中的平行粘结模型模拟材料的强度特征；以某一混凝土坝为实例，利用开发的随机网格生成算法建立了大坝模型，对爆炸载荷下混凝土坝的溃坝过程进行了模拟仿真。余天堂[6]基于任意拉格朗日-欧拉算法，使用完全耦合模型，并利用LS -DYNA软件模拟受到接触爆炸载荷的混凝土坝的动态响应，提出了大坝混凝土损伤演变过程和大坝所承受的主要应力。Linsbauer[7]研究和论述了在冲击荷载条件下重力坝脆弱区潜在的危险，基于断裂力学分析了受影响坝体过度损坏的可能性。张启灵等[8]基于ABAQUS有限元程序，采用声学介质描述流体，考虑坝体-水流固耦合作用和混凝土的受拉及受压损伤，对某典型的重力坝坝段结构进行爆炸荷载作用下的损伤塑性时程分析。

由以上研究可以看出，采取必要措施提高坝体抗爆性能势在必行，而已有的研究结果表明，泡沫混凝土能够有效提高受保护物体的抗爆能力。杜玉兰等[9]采用室内小型平板模型试验，利用平面药包爆炸直接加载，系统研究了爆炸作用下复合防护结构的抗爆性能以及含泡沫混凝土吸能层结构中的泡沫层厚度和位置对抗爆性能的影响，结果表明采用含泡沫混凝土复合材料能大幅提高结构的抗爆性能。倪小军[10]研究了泡沫材料在水下爆炸冲击载荷作用下的吸能特性及其对结构动态响应的影响，以及泡沫材料爆炸冲击特性的表征。张景飞等[11]通过对自行设计的混凝土和泡沫混凝土复合防护结构在爆炸冲击波作用下的吸能特性进行了试验研究，就不同密度的泡沫混凝土在复合防护结构中的抗爆吸能作用进行了比较，结果表明复合防护结构中泡沫混凝土层具有明显的吸能作用，且其抗爆吸能是以牺牲外层及泡沫混凝土来保护下层结构为代价的。但是到目前为止，针对表面附着有泡沫混凝土保护层的情况下坝体抗爆性能的研究还相对较少，因此这是一个值得深入探讨的课题。

2 有限元模型

2.1 应变率模型

近些年来，各种各样用于模拟混凝土爆炸荷载响应的材料模型相继出现，在非线性有限元程序LS -DYNA中，材料编号为72R3，96，111和272的材料都能用于模拟混凝土爆炸荷载作用下的响应。在爆炸荷载下混凝土材料会表现出率相关性，其应变率会达到10 s-1～1000 s-1甚至是更高。混凝土材料的抗压强度和抗拉强度分别会提高100%和600%[12，13]。Malvar等[14]提出了用动态增长因数DIF来表示静态强度的动力放大率。对混凝土的抗压强度来说，动态增长因数可以表示为

 

(1)

式中表示应变率，介于30×10-6 s-1～300 s-1之间，是静态应变率，为混凝土的静态抗压强度。混凝土抗压强度的动力提高因子如图1(a)所示。混凝土抗拉强度的动力提高因子可以表示为

 

(2)

式中为应变率，介于10-6 s-1～160 s-1之间，为静态应变率，且混凝土抗拉强度的动力提高因子如图1(b)所示。

2.2 RHT混凝土本构模型

RHT混凝土动力损伤本构模型由Riedel等[15]发展而来。该本构模型通过引入弹性极限面、失效面和残余失效面这三个极限面来分别描述混凝土的初始屈服强度、失效强度以及残余强度的有关变化规律[16]。RHT本构模型既考虑了混凝土的大应变、高应变率和高压强的效应，又兼顾了混凝土应变的硬化和软化以及应力偏量第三不变量的影响。因为其能够描述混凝土从弹性阶段到损伤失效整个过程的变化状态，故广泛应用于模拟爆炸荷载下混凝土的动力响应以及损伤断裂过程。上述提到的三个不同的极限面如图2所示，所描述的混凝土性能沿着不同的拉压子午线衰减，并受到应变率效应的影响。

 

 

图1 混凝土抗压强度动力提高因数和抗拉强度动力提高因数[14]

Fig.1 DIF for compressive strength of concrete and for tensile strength of concrete

 

图2 RHT本构模型极限面

Fig.2 Limit surfaces of the RHT model

σy = [A+B(εp)n](1+C lnε*)[1-(Th)m]

[5] 王冰玲,刘 军.爆炸载荷下混凝土坝溃坝过程的连续仿真[J].系统仿真学报,2014,26(1):159-162.(WANG Bing-ling,LIU Jun.Numerical simulation of process of concrete dam-break under explosive loa-ding[J].Jouranal of System Simulation,2014,26(1):159-162.(in Chinese))

[2] 李本平,卢文波.制导炸弹水平侵彻爆炸作用下混凝土重力坝毁伤效应数值仿真[J].爆破,2007,24(1):1-5.(LI Ben-ping,LU Wen-bo.Numerical simulation for the damage of concrete gravity dam under horizontal penetration and explosion of GBU-28[J].Blasting,2007,24(1):1-5.(in Chinese))

对以废弃油脂和麻风树油脂为原料的两种生物柴油进行GC-MS分析。图1所示为废弃油脂生物柴油的谱图，与图库比对，其主要组分列于表 1。以麻风树油脂为原料的生物柴油的 GC-MS分析结果如图2所示，主要组分见表2。对比表1和表2，可以看出生物柴油的典型组分是棕榈酸甲酯（C17H34O2）、硬脂酸甲酯（C19H38O2）、油酸甲酯（C19H36O2）、亚油酸甲酯（C19H34O2）和亚麻酸甲酯（C19H32O2）。不难发现，原材料的差异会影响到组分构成以及组分含量，为进一步说明这一点，结合相关文献[9-10]，列举其他原料制备的生物柴油组分，汇总于表3。

为了考察东亚夏季风对两类El Nio事件的响应及对我国夏季降水的影响,对比850 hPa风场在两类El Nio事件发生时的空间分布,比较两类El Nio事件对应的大尺度环流异常的差异。

 

(3)

式中 Y c(p*)为压缩子午线方程，并且

 

(4)

式中 p*为标准化静水压力且p* = p/fc,p表示实际的静水压力表示动力放大因数为应变率;fc为准静态单轴抗压强度为准静态单轴抗拉强度；r(θ)是偏平面上的角隅函数，为第二偏应力不变量J2和第三偏应力不变量J3及拉压子午线比率Ψ的函数，可以表示为

 

 

(5)

Lode角θ可以由下式得到。

 

(6)

2.2.2 弹性极限面

馆服的穿着能起到对馆员行为规范的约束作用，会给馆员一种心理暗示，自觉遵守学校和图书馆纪律，注意自己的言行举止，维护图书馆形象。无论是工作和休息期间，馆服会时刻提醒你是一名图书馆工作人员，是知识的传播者，身上理应体现文化人的要求，对自己的行为举止在一定程度上要做到自律，遵守学校的规章制度，避免和减少不良行为的发生。

[6] Yu T T.Dynamical response simulation of concrete dam subjected to underwater contact explosion load[A].2009 World Congress on Computer Science and Information Engineering[C].2009.

Y e = Y f×ω×Fcap(p)

(7)

式中 ω是沿径向的弹性强度与失效强度的比值，鉴于混凝土材料在其最大抗压强度的30%及最大抗拉强度的50%～80%时表现出非线性行为，故ω呈线性变化并与混凝土的单轴抗拉强度和抗压强度密切相关。Fcap(p)的取值范围是(0，1)，用于限制静水压力作用下的弹性力偏应力。

电气工程及其自动化作为一门新兴的朝阳学科，是科学技术与生产实践在新时代紧密结合的必然产物，在国民经济生产中占据着重要地位，也是提升我国工业水平、提高人民生活质量的重要技术。因此对科学合理的应用，是当今必不可少的探索过程。

2.2.3 残余失效面

残余失效面Y r用于描述混凝土材料在完全破碎时的强度，其可以表述为

Y r = B×(p*)M

(8)

式中 B和M分别是残余失效面常数和残余失效面指数。

3 模型实验

本文用Johnson-Cook模型来模拟钢筋材料，Johnson-Cook模型尤其适用于钢筋在大变形及高应变情况下的性能分析[17]，其屈服应力σy可以表示为

二是冲突双方事后都显示了理性。利夫西将军事后派人去当事的中餐馆结账，还给当时的华裔送餐员留下了小费，说明这位将军知错能改，并不是一味地任性，死要面子。考虑到利夫西将军与华裔送餐员地位悬殊，利夫西能这样做，姿态还是不错的。餐馆老板并不想一味把事情搞大，也有意撤诉。

单片机采用AT89C52，时钟晶体选择AT切型的11.059 2MHz，传感器SHT11温湿度传感器的DATA管脚接P1.1，SCK接P1.0。这两个管脚要接4.7kΩ的上拉电阻。仿真软件下不显示SHT11的电源和地端，实际接线时，在电源和地间接入100nF的去耦滤波电容。因为数据要上传到PC，为了实现单片机和PC间的电平转换，实际电路要采用RS232/RS485转换接口，但在仿真软件元件库中无此器件，仿真时采用RS485标准的器件MAX487的输出端接模拟PC串口的虚拟终端。电路连接如图4所示。

2.2.1 失效面

(9)

式中 εp是等效塑性应变，ε* 是规范化的塑性应变率，温度比值Th是室温T r和溶解温度T m的函数，即T h = (T-T r)/(T m-T r)。室温T r为300 K，熔化温度T m =1793 K，比热为477 J/kg·K，相对密度=7.83 g/cm3，剪切模量G=81.8 GPa，屈服强度A =792 MPa，硬化常数B=510 MPa，硬化指数n=0.25，应变速率常数C=0.014，热软化指数m=1.03。

炸药材料爆炸产生的压力P采用JWL状态方程来描述，其表达式为

 

 

(10)

式中 V是相对体积，E=7.0×109 J/m3是每单位体积炸药的初始内能，A=3.712×1011 Pa，B=3.231×109 Pa，R1=4.15，R2=0.95，ω=0.3，以上取值都是基于相关实验得到的炸药材料特征常数。

由于爆炸会引起周围空气的强烈压缩，在周围产生以炸点为核心的冲击波，此时，理想气体的线性多项式状态方程可以表示为

无论是从普希金的“人的影响短暂而微弱，书的影响则广泛而深远”一语，还是从歌德的“读一本好书，就是和许多高尚的人谈话”，都可以悟出读书对于人的重要性。

 

(11)

式中 是空气的压缩程度，且 0-1,是当前的空气密度，0是初始的空气密度。对于理想气体，C1,C2和C3都为0，C4 = C5 = 0.401，密度为1.293 kg/m3，初始内能E=2.5×105 J/m3。

游戏为学生创造了一番愉快的天地，它是培养学生良好行为习惯和优秀品德的重要形式之一。德国著名学前教育家福禄培尔特别强调游戏对儿童教育的重要性。他认为游戏是人类童年时代生活中最为快乐的一种现象，是儿童本能的表现，富有浓厚的生活气息、充满儿童情趣的游戏是启人心智的力量。在活动设计时，教师不妨从儿童喜闻乐见的游戏中找找灵感，或许能使课堂趣味盎然。

为了验证RHT混凝土本构模型对爆炸荷载下模拟混凝土动态响应的适应性和准确性，用适当规模的数值模拟与现场的试验结果进行对比是很有必要的。本文采用1000 mm×1000 mm×40 mm的钢筋混凝土板，几何模型的尺寸如图3所示[18]。钢筋的直径为6 mm，横向和纵向的间隔均为75 mm，则其配筋率为1.43%。钢筋的屈服强度为600 MPa，弹性模量为200 GPa。混凝土的抗压强度为39.5 MPa，抗拉强度为4.2 MPa，弹性模量为28.3 GPa。现场试验中，TNT炸药的重量分别为0.31 kg和0.46 kg，起爆距离为400 mm。RHT本构模型的相关参数列入表1。

 

图3 钢筋混凝土板尺寸及炸药布置

Fig.3 Dimensions of the reinforced concrete slab and layout of explosive

表1 RHT本构模型参数取值[1]

Tab.1 Parameters of the RHT model

  

参数 取值 参数 取值密度/(kg·m-3)2.314×103断裂压缩应变率Ec3.0×1022剪切模量/Pa1.67×1010断裂拉伸应变率Et3.0×1022B0 1.22压缩应变率因子 0.032B1 1.22拉伸应变率因子 0.036T1/Pa3.527×1010抗压屈服面参数 0.53T2 0.0抗压屈服面参数 0.70失效面常数A 1.6剪切模量换算系数 0.50失效面常数N 0.61损伤参数D1 0.04抗压强度fc/Pa3.5×107损伤参数D2 1.0相对抗剪强度 0.18失效面常数Af 1.6相对抗拉强度 0.10失效面指数Nf 0.61拉压子午线比值 0.6805A1/Pa3.53×1010脆柔转变系数B 0.0105A2/Pa3.96×1010参考压缩应变率3.0×10-8A3/Pa9.04×109参考拉伸应变率3.0×10-9压实应力/Pa6.00×109

 

图4 1/4有限元模型

Fig.4 A quarter finite element model

由于实验模型的对称性，只需建立1/4模型即可完成对实验的数值模拟，其1/4模型如图4所示。在LS -DYNA中，采用Eulerian和 Lagran-gian 耦合的方法建立爆炸冲击下的钢筋、混凝土和空气的全耦合模型，其中钢筋和混凝土采用 Lag-rangian 网格进行模拟，网格大小为2.5 mm，空气和炸药采用Eulerian网格进行模拟，网格大小为10 mm。周围空气包络整个混凝土板和炸药，计算域外围距离混凝土板四边各50 mm，距离混凝土板上表面500 mm，下表面50 mm，在空气周围采用无反射的透射边界，使得空气中的能量可在截断处向无限域散发。

在0.31 kg TNT爆炸冲击下，将数值模拟的结果与现场试验结果进行对比。现场试验结果中，混凝土板的中心挠度最大值达到15 mm，在钢筋混凝土板的迎爆面，如图5(a)所示，出现了平行于固定约束边界的裂纹以及环向和向四周散射的裂纹，平行于固定约束边界的裂纹贯穿于整个混凝土板，向四周散射的裂纹有一部分延伸到板的边界，环向裂纹的半径约为280 mm；背爆面主要是由于反射形成的拉力超过了混凝土的抗拉强度而形成的开裂破坏，在中心形成了半径约为90 mm的混凝土脱落区，如图5(c)所示。数值模拟结果中，混凝土板的中心挠度最大值达到14 mm，图5(b)所示的迎爆面裂纹走向以及环向裂纹的直径与现场试验中大致相同，但是环向裂纹的半径约为290 mm；图5(d)所示的背爆面中心形成半径约为100 mm的混凝土完全破碎区。在0.46 kg TNT爆炸冲击下，现场试验结果中混凝土板的中心挠度最大值达到35 mm，钢筋混凝土板的迎爆面和背爆面的破坏范围明显加大，破坏程度明显加深，向四周散射的裂纹几乎全部延伸到板的边界，环向裂纹的半径约为300 mm，如图6(a)所示；背爆面混凝土脱落区域半径约为120 mm，如图6(c)所示。数值模拟得到的混凝土板中心挠度最大值达到32 mm，环向裂纹的半径约为310 mm，如图6(b)所示；背爆面中心形成半径约为140 mm的混凝土完全破碎区，如图6(d)所示。

 

 

 

 

图5 0.31 kg TNT爆炸冲击下钢筋混凝土板的破坏

Fig.5 Experimental and numerical results for reinforced concrete slab under explosive loads of 0.31 kg TNT

 

 

 

 

图6 0.46 kg TNT爆炸冲击下钢筋混凝土板的破坏

Fig.6 Experimental and numerical results for reinforced concrete slab under explosive loads of 0.46 kg TNT

总体来说，数值模拟得到的混凝土板的损伤比现场试验得到的结果稍微偏大。这主要是因为现场试验中混凝土板是采用塞木板的方式进行固定，故混凝土板约束端并未按照试验描述进行全约束，木板的振动以及木板本身的材料属性会消耗和吸收一部分能量；而数值模拟中混凝土板采用全约束的方式进行固定，混凝土板没有与除空气外的任何物质相接触，故混凝土板内能量的耗散少了一条途径，混凝土板的损伤也必然会有所增加。但是综合以上描述可以看出，数值模拟得到的混凝土板迎爆面和背爆面损伤区域的分布以及板的中心挠度最大值与试验结果大致相符。说明采用数值方法能够有效模拟爆炸冲击荷载下混凝土板的损伤分布和裂纹扩展，故采用RHT模型进行爆炸荷载下混凝土板动态响应的数值模拟是可行的。

4 坝体抗爆反应

4.1 泡沫混凝土

泡沫混凝土材料能够部分吸收爆炸产生的冲击能量，并且能起到减振的效果。对于水工大坝，可以在表面覆盖适当厚度的泡沫混凝土以减小坝体在爆炸冲击作用下的动力响应，起到很好的保护效果。泡沫混凝土材料采用可压扁的Mat_Crushable_Foam泡沫模型，该模型受应变率效应影响明显，其弹性模量为常数且应力为弹性响应，具体表达式为

 

(12)

式中 εi j为应变率，E为弹性模量，t为时间，在准静态时其应力应变曲线如图7所示，常用的泡沫混凝土参数[19]如下，密度799 kg/m3，泊松比0.1，杨氏模量342.2 MPa，抗拉强度0.20 MPa，阻尼系数0.10。

4.2 坝体抗爆

本文所采用的模型是我国云南澜沧江上游某重力坝的某一坝段，其高程为159 m，正常蓄水位为149 m，坝段宽度为15 m，顺河向坝底宽度为122 m，坝顶宽度为16 m，其有限元网格模型如图8所示，计算结构模型如图9所示。

假设爆炸点离坝体上游面的水平距离为10 m，分别模拟其在水下10 m,20 m和30 m处在2 t TNT炸药爆炸下的动力响应，并计算在无泡沫混凝土保护层和1 m和2 m泡沫混凝土保护层情况下坝体的动力响应。炸药在水下爆炸会产生高压冲击波、压力波以及气泡脉动。水下爆炸产生的初始冲击波作用在坝体上游面，其中有一部分冲击波会发生反射并传向库水区域，另一部分则将入射到坝体内部继续传播。

 

图7 泡沫混凝土的应力应变曲线

Fig.7 Stress-strain curve of foamed concrete

 

图8 坝体有限元模型混凝土

Fig.8 Finite element model of the dam

 

图9 大坝结构模型示意图

Fig.9 Structure model of the dam

从图10～图12的坝体损伤示意图可以看出，在没有泡沫混凝土保护层时，坝体上游面主要损伤区域的位置总是随着起爆深度的增加而向坝体底部移动，总体上在坝高方向的宽度为25 m左右，深度为4.5 m左右；下游面的损伤面积则整体向逐渐变小的趋势发展。由于缺少泡沫混凝土保护层，大量的冲击波从坝体上游面入射并传播到下游面，这些入射波的一部分将在下游面发生反射而形成强拉伸波，故下游面的混凝土会发生较大面积的拉伸损伤破坏，其损伤破坏最大深度达到3.7 m。正对爆破点的坝体损伤横跨整个坝体上游面和下游面。整个坝头部位有较大面积的损伤，由此可见坝头是坝体抗爆的薄弱部位。当上游表面附着有1 m厚的泡沫混凝土保护层，坝体上游表面的损伤明显变小，下游面的损伤较之无保护层情况也明显减小，这是因为上游面覆盖的泡沫混凝土具有消波吸能减振的作用，泡沫混凝土的大塑性变形吸收了大量的冲击和振动能量，使冲击振动能量在泡沫混凝土中急剧衰减，从而有力减少了坝体的振动效应。下游面的损伤明显大于上游面，这仍然是由于一部分入射冲击波在下游面的反射造成的强拉伸波大于混凝土的抗拉强度而产生的拉伸破坏。当上游表面附着有2 m厚的泡沫混凝土保护层时，则坝体上游面已经没有明显的损伤，与附着有1 m厚保护层的情况相比，下游面的损伤区域大大减小，随着爆破深度的增加，减小的程度越发明显。这是因为2 m厚泡沫混凝土保护层能够吸收相对更多的冲击和振动能量，更进一步地减少了坝体的振动效应。坝体下游面的损伤仍然是部分入射波造成的拉伸破坏。

图13～图15是坝顶中心点的位移时程曲线，可以看出，无泡沫混凝土保护层时，在爆炸荷载冲击下坝顶位移明显比有混凝土保护层时大得多，而且随着混凝土保护层厚度的增加，坝顶位移明显减小，位移增加速度趋于平缓，这说明泡沫混凝土能够有效保护坝体，减小因爆炸荷载引起的坝顶位移，有效提高坝体的抗爆性能。

 

 

 

图10 水下10 m爆炸坝体损伤

Fig.10 Dam damage of 10 m under water explosion

 

 

 

图11 水下20 m爆炸坝体损伤

Fig.11 Dam damage of 20 m under water explosion

 

 

 

图12 水下30 m爆炸坝体损伤

Fig.12 Dam damage of 30 m under water explosion

 

图13 水下10 m爆炸坝顶位移时程

Fig.13 Displacement and time curve of 10 m under water explosion

 

图14 水下20 m爆炸坝顶位移时程

Fig.14 Displacement and time curve of 20 m under water explosion

 

图15 水下30 m爆炸坝顶位移时程

Fig.15 Displacement and time curve of 30 m under water explosion

5 结 论

RHT本构模型能够有效模拟混凝土在爆炸荷载作用下的动力反应。在没有泡沫混凝土保护层时，坝体上游面主要损伤区域的位置总是随着起爆深度的增加而向坝体底部移动，且坝头是坝体抗爆的薄弱部位。当上游表面附着泡沫混凝土保护层时，坝体上游表面的损伤明显变小，下游面的损伤较无保护层情况也明显减小，说明泡沫混凝土能够有效减小混凝土大坝在爆炸荷载作用下的损伤，提高混凝土大坝的抗爆性能，所以对混凝土重力坝来说，在其表面附着一层泡沫混凝土，是一项良好的抗爆措施。通过数值模拟结果可得，当混凝土坝上游表面附着有2 m厚的泡沫混凝土保护层时，能够对坝体起到很好的保护作用，在爆炸冲击作用下，能够保证坝体的安全性。

参考文献(References):

[1] 王高辉,张社荣,卢文波,等.水下爆炸冲击荷载下混凝土重力坝的破坏效应[J].水利学报,2015,46(6):723-731.(WANG Gao -hui,ZHANG She -rong,LU Wen-bo,et al.Damage effects of concrete gravity dams subjected to underwater explosion[J].Journal of Hydraulic Engineering,2015,46(6):723-731.(in Chinese))

失效面Y f的表达式为

[3] Zhang S,Wang G,Wang C,et al.Numerical simulation of failure modes of concrete gravity dams sub-jected to underwater explosion[J].Engineering Fai-lure Analysis,2014,36:49-64.

随着全球化的进程，各学科之间相互共生，促进发展。一部经典的中国文化著作涵盖不只是文学知识，还包括语言学、文学、艺术、心里学、社会学、哲学、政治学等各个文化层面。当前单一的语言翻译教学体系显然无法满足文化“走出去”的深度和广度需求。为了提高翻译专业学生的综合素质和翻译教学质量，北外、上外、广外三所试点院校以及其他翻译院系都提出，应该在翻译本科阶段开设古代汉语、现代汉语、中国经典作品研读、中国现当代文学、高级汉语写作、中西方文化、中西方哲学、中西思想史和逻辑学等课程，打下比较坚实的国学基础，培养学生的思考能力，提高其认识水平和文化素养。

[4] 徐俊祥,刘西拉.水中爆炸冲击下混凝土坝动力响应的全耦合分析[J].上海交通大学学报,2008,42(6)：1001-1004.(XU Jun-xiang,LIU Xi-la.Full coupled simulation of concrete dams subjected to underwater explosion[J].Journal of Shanghai Jiaotong University,2008,42(6):1001-1004.(in Chinese))

数字化转型推动着当今工业市场的发展。汉诺威工业博览会将举办90多场会议，涵盖工业领域的数字化转型以及其他重大趋势和话题。

在RHT本构模型中，弹性极限面Y e用来考虑材料的应变硬化效应，并由失效面确定，其可以表达为

[7] Linsbauer H.Hazard potential of zones of weakness in gravity dams under impact loading conditions[J].Frontiers of Architecture and Civil Engineering in China,2011,5(1):90-97.

[8] 张启灵,李端有,李 波.水下爆炸冲击作用下重力坝的损伤发展及破坏模式[J].爆炸与冲击,2012,32(6):609-615.(ZHANG Qi-ling,LI Duan-you,LI Bo.Damage propagation and failure mode of gravity dam subjected to underwater explosion[J].Explosion and Shock Waves,2012,32(6):609-615.(in Chinese))

[9] 杜玉兰,王代华,刘殿书,等.含泡沫混凝土层复合结构抗爆性能试验研究[A].首届全国水工抗震防灾学术会议[C].南京，2006.(DU Yu-lan,WANG Dai-hua,LIU Dian-shui,et al.Experimental research on the characteristics of anti-blast compound structures including foam concrete[A].1s t National Conference on Hydraulic Seismic Disaster Prevention[C].Nanjing,2006.(in Chinese))

[10] 倪小军.泡沫材料对水下爆炸能量吸收及其冲击特性的研究[D].中国科学技术大学,2012.(NI Xiao -jun.Effect and Description of Foam Materials on the Energy Absorption to UNDEX Shock[D].University of Science and Technology of China,2012.(in Chinese))

[11] 张景飞,冯明德,陈金刚.泡沫混凝土抗爆性能的试验研究[J].混凝土,2010(10):10-12.(ZHANG Jing-fei,FENG Ming-de,CHEN Jin-gang.Study on the knock characteristic of foam concrete[J].Concrete,2010(10)：10-12.(in Chinese))

先说“太幼稚”。正如网友所说，无论“上下班坚决不迟到早退，若临时请假必须事前写好请假条”，还是“杜绝开发商请客吃饭与办证速度挂钩问题，杜绝一切有业务往来的饭局”，都属于“幼儿园认知级别”的问题，本就是基本常识，根本就不需要“承诺”。的确，如果不迟到早退也要进行公开承诺，那就说明要么该部门平时的确连纪律底线也没有守住，要么就是敷衍了事，用原本不需要承诺的问题滥竽充数、走过场。

[12] Bischoff P H,Perry S H.Compressive behaviour of concrete at high strain rates[J].Materials and structures,1991,24(6)：425-450.

[13] Malvar L J,Ross C A.Review of strain rate effects for concrete in tension[J].Materials Journal,1998,95(6)：735-739.

[14] Malvar L J,Crawford J E,Morrill K B.K&C Concrete Material Model,Release III-Automated Generation of Material Model Input.Karagozian and Case structural Engineers[R].Technical report 2000 TR-99-24.3,2000.

[15] Riedel W,Thoma K,Hiermaier S,et al.Penetration of reinforced concrete by BETA-B-500 numerical analysis using a new macroscopic concrete model for hydrocodes[A].9t h International Symposium on the Effects of Munitions with Structures [C].1999.

[16] Tu Z,Lu Y.Evaluation of typical concrete material models used in hydrocodes for high dynamic response simulations[J].International Journal of Impact Engineering,2009,36(1):132-146.

[17] Johnson G R,Cook W H.A constitutive model and data for metals subjected to large strains,high strain rates and high temperatures[A]. 7t h International Symposium on Ballistics[C].1983.

[18] Wang W,Zhang D,Lu F,et al.Experimental study on scaling the explosion resistance of a one -way square reinforced concrete slab under a close-in blast loading[J].International Journal of Impact Engineering,2012,49:158-164.

天晓得，这些事别人除过自己的工作会不会干，一天下来，人都累成个瘫人儿，连老婆把洗脚水端到跟前都懒得把脚伸到盆子里。

[19] 唐德高,王昆明,贺虎成,等.泡沫混凝土回填层在坑道中的耗能作用[J].解放军理工大学学报(自然科学版),2006,7(4):365-370.(TANG De -gao,WANG Kun-ming,HE Hu-cheng,et al.Energy dissipation mechanism of foamed concrete backfill layers in underground tunnels[J].Journal of PLA University of Science and Technology(Natural Science Edition),2006,7(4):365-370.(in Chinese))