弹性纵波垂直通过非线性结构面的傅里叶级数分析法

利用傅里叶级数计算任意函数形式的弹性纵波垂直通过非线性结构面的透射波和反射波的速度波形。采用双曲模型描述结构面的变形特性,基于位移不连续模型,结合波振面处动量守恒定律,推导了弹性纵波垂直通过非线性结构面的基本方程。假设应力波在含结构面岩体中传播时,结构面的存在不改变应力波波形函数的最小正周期,运用傅里叶级数理论和周期延拓方法,得到了任意函数形式的弹性纵波垂直入射时透射波和反射波速度波形的傅里叶级数解,并验证了傅里叶级数解是合理的。利用傅里叶级数解,分析了单一频率的正弦谐波入射至结构面时,透射波中各阶谐波的振幅和相位与谐波阶数的关系。研究结果表明,各阶谐波的振幅与阶数呈负指数关系衰减,前7阶谐波振幅的衰减指数为谐波阶数的二次函数,当谐波阶数大于7时,衰减指数为谐波阶数的一次函数;各阶谐波的相位与谐波阶数呈线性关系。     

基于MPGA的复杂应力状态边坡稳定性分析

极限平衡法在边坡稳定性分析中处于主导地位,该算法通常建立在一定的假设之上,没有考虑边坡土体的非线性应力-应变关系,不能用于未达到极限状态以及加固后边坡的稳定性分析。为了解决以上问题,将有限元计算与多种群遗传算法（MPGA）相结合,建立一种基于MPGA的复杂应力状态边坡稳定性分析通用模型,通过数值应力场求解安全系数,并为多种群遗传算法构建适应度函数;再利用多种群遗传算法为安全系数的计算提供滑移面。为了保证分析的高效、合理,根据滑移面发展趋势,动态产生初始滑移面,并增加一个滑移面约束条件。最后,通过均质边坡和软弱夹层边坡两个典型算例分析验证了该方法的合理性;通过分析土钉加固的软弱夹层边坡,证明了该方法可用于加固等复杂应力状态边坡的稳定性分析。     

南北地震带壳幔结构成像研究

<正>对地壳和上地幔的速度结构进行成像是研究地球深部结构及动力学问题的一个有效途径。南北地震带地区是青藏高原与相邻块体的过渡带,这里的地质构造十分复杂,地震活动频繁。在南北地震带地区进行壳幔结构成像研究,有利于提高对青藏高原及其与相邻块体相     

2008-2015年北京高速公路道面结冰特征分析

道面结冰是北京地区高速公路交通安全的主要危害因素。利用北京地区2008-2015年28个ROSA交通气象站的结冰数据分析了北京地区高速公路道面结冰的基本特征。结果表明：（1）北京市高速公路道面结冰空间差异显著,不同高速公路结冰特征显著不同,同一高速公路不同路段结冰特征也明显不同,这与高速公路各路段局地气候存在差异有关。（2）道面结冰主要发生在北部和东部高速公路所在路段,以小汤山西桥站、顾家庄桥站、六道口桥站和丁各庄桥站道面结冰灾害最为严重。冬季道面结冰月变化显著,11月和3月,结冰次数较少,结冰持续时间短;12月、1月和2月结冰次数多,结冰持续时间长。（3）持续时间越长的结冰过程结冰次数越少。不同持续时间的结冰过程在各时段的结冰次数和累积结冰时长均存在明显差异,两者在各时刻的演变特征基本相同,峰值均出现在22：00,且在00：00之后呈明显减少趋势。（4）大部分站点慢车道比快车道更容易结冰,少部分站点快、慢车道结冰率近似,有的站点甚至快车道结冰率远高于慢车道。快、慢车道水或覆盖物厚度与交通站点所处位置的局地气象条件密切相关。（5）气温、道面温度和大地温度与水或覆盖物的厚度变化呈反相位,即温度越低,水或覆盖物的厚度越大。这可为北京地区道面结冰预报预警方法的开展提供着眼点和依据。     

陇中黄土高原LAS在非均匀下垫面的观测优势及其在区域陆面模式验证中的初步应用

目前通用的通量观测技术涡动相关仪（EC）在区域陆面模式验证中存在能量不闭合和空间代表性有限的问题,寻求改进EC观测热通量的新技术是提高陆面模式验证效果的关键环节,大孔径闪烁仪（LAS）的出现有效改善了这一现状。基于黄土高原定西站2010年1月和6月的同步综合观测资料以及目前比较有代表性的陆面过程模式CLM的模拟数据,分析研究了LAS对EC观测地表能量不平衡问题的改进以及LAS对EC在区域陆面模式验证中的提高,结果表明：LAS可有效解决EC观测存在的地表能量不平衡问题,提高EC的地表能量闭合度,在非均匀下垫面LAS观测优势突出;利用LAS观测的感热通量进行区域陆面模式的验证,能够很大程度地避免EC能量不闭合和空间尺度不匹配在验证中造成的偏差,LAS观测更适合于大尺度模拟的验证,验证效果更好。     

滇西腾冲更新世粗面安山岩Ar-Ar年代学、地球化学特征及其构造意义

腾冲芒棒地区广泛发育新近纪以来的火山岩,本文首次对该区粗面安山岩进行的高精度40Ar/39Ar同位素定年,其坪年龄分别为2.1±0.1Ma和2.0±0.1Ma,代表其喷溢时代为早更新世。岩石地球化学特征研究显示,粗面安山岩属高钾钙碱性系列,富集大离子亲石元素、亏损高场强元素,具有弧火山岩与大陆板内火山岩的双重属性,推测其形成于后造山环境。综合研究表明,腾冲地块更新世火山岩的形成与印度-欧亚板块俯冲-碰撞造山隆升后的伸展垮塌、富集地幔部分熔融作用密切相关,而印支陆块向东南逃逸以及腾冲弧形走滑断裂带(龙川江断裂带、大盈江断裂带)的活化促使了该期火山岩的喷发与侵位。     

基于修正Drucker-Prager准则的岩石次加载面模型

为模拟岩石在循环加、卸载作用下的变形特性,结合Drucker-Prager(简称D-P)屈服准则与次加载面理论,建立了能考虑岩石材料实际抗拉强度和应力角效应的循环加、卸载模型。首先考虑岩石在三轴压缩和在三轴拉伸状态下不同的强度特性,在传统的D-P屈服准则中引入角隅模型,并借鉴殷有泉提出的D-P-Y准则,形成了修正的D-P屈服准则;进一步结合修正屈服准则和次加载面理论,提出了适用于岩石循环加、卸载的次加载面模型,并通过自编程序,实现了循环加、卸载模型的数值表达与岩石循环加、卸载试验数值仿真。计算结果表明,基于修正D-P屈服准则的次加载面模型能较好地描述岩石的曼辛效应和棘轮效应;通过改变抗拉强度的大小,可以发现随着抗拉强度的减小,岩石的变形量相应增加,且增加的幅度逐渐提高,说明修正模型可准确反映抗拉强度对岩石的循环加、卸载变形的影响。     

基于GRASS GIS与TIN滑动面的边坡三维极限平衡方法研究

基于开源地理信息系统(GRASS GIS),建立三维极限平衡Hovland模型,研究了边坡稳定性计算方法。采用空间平面投影方法对椭球滑动面中各组成单元面积计算方法进行改进,提高了计算效率。考虑边坡主滑方向与各单元倾向之间的关系,给出了能反映正负变化的滑动面各单元安全系数表达式。通过引入约束Delaunay三角形剖分方法与网格优化算法来生成滑面不规则三角形格网(简称TIN)模型,弥补了Grid模型在滑动面边界准确表达方面的不足。相比Grid模型,在相同网格分辨率条件下,TIN滑动面获得的边坡整体安全系数更加接近正确解,能够有效提升计算精度。通过考题与工程实例,对所提方法与程序模块的正确性进行了验证。结果表明,通过分析滑动面上各单元的安全系数,可以进一步划分不稳定区域与阻滑区,能够为滑坡治理措施的制定提供重要参考。     

隧道锚围岩抗拔机制及抗拔力计算模式初步研究

在普立特大桥隧道锚现场模型试验的基础上,采用数值模拟技术揭示了隧道锚围岩变形破坏过程:围岩破坏面从锚体底部与围岩接触面附近启裂,并逐渐向外呈圆台状扩散,破坏形式为拉剪破坏。并且,锚体前部临空岩体被拱出而发生拉破坏。破坏面上的应力分布随着拉拔荷载增大而发生复杂变化。基于此,通过在破坏面上建立力的平衡关系,提出了隧道锚围岩抗拔力计算模式。该计算模式与现有文献不同,体现了夹持效应以及破坏面上的复杂应力变化。破坏面上的应力分布需要通过模型试验和数值模拟论证得到。今后,在针对不同强度、不同结构特征的岩体进行全面试验分析的基础上,可以对破坏面形态和应力大小进行取值建议。采用该模式验证了试验结果,估算得到大桥原型锚碇的极限抗拔力非常大。目前隧道锚设计普遍偏于保守,隧道锚在中、软岩中仍然可以使用。讨论了破坏面形态特征可能的变化、岩体结构特征对抗拔力的影响等问题。     

自适应差分进化算法在边坡滑面搜索中的应用

边坡三维临界滑裂面的确定及其安全系数的求解是边坡稳定分析中的关键问题,而现有分析方法尚难以方便、可靠地获得三维滑裂面。在有限元数值应力场的基础上,计算边坡的安全系数。通过构造扩展椭球形滑裂面,利用自适应差分进化(SaDE)算法结合反向移动操作,实现了复杂边坡三维临界滑裂面的搜索。扩展椭球滑裂面构造方便,搜索灵活,比传统的球面或椭球面更能反映实际工程情况。整个求解过程为显式求解,不存在收敛性问题。算例考核结果表明,提出的方法具有良好的搜索效率与全局寻优能力,结果与刚体极限平衡法、有限元强度折减法吻合,计算结果合理可靠。通过实际工程边坡,进一步验证了该方法可方便、有效地应用于具有复杂岩性及环境条件的边坡三维稳定性分析。