三角网格谱元法地震波场数值模拟

谱元法结合了有限元法的灵活性和谱方法的指数收敛性,高效且高精度,是近年来发展的一种重要的地震波场数值模拟方法.经典的谱元法采用四边形(六面体)网格,利用一维Gauss-Legendre-Lobatto(GLL)积分的张量积得到对角的质量矩阵,以大大提高计算效率,但是四边形(六面体)网格不能够灵活地刻画复杂的几何模型的弯曲界面.为此,在谱元法中引入三角形(四面体)网格到二维(三维)是十分必要的.不同于经典的谱元法,在非结构化网格中不能使用GLL积分的张量积,使得非结构化网格的谱元法的实现存在着诸多的困难.目前,比较流行的三角网格谱元法,通过使用KoornwinderDubiner(KD)正交多项式,并正交化这些KD多项式构建基函数,同时利用重合的插值节点和积分节点以获取对角的质量矩阵;它所使用的积分点为优化的点集——Fekete点,且这些积分点能与四边形网格完全耦合.相比于四边形,三角网格谱元法能显著提高复杂模型的描述能力,对起伏地表模型有很大优势.本文引入高效的最佳匹配层(PML)吸收边界条件,并通过数值试验将三角网格谱元法与经典的谱元法进行对比研究.相比于经典的谱元法,三角网格谱元法显著缺点为较低的计算精度.对于7阶谱元,为了能够精确地模拟面波,三角网格谱元法需要在每个最短的面波波长内至少有11个采样点,然而经典的谱元法仅需4个采样点,并且前者所需的内存量约为后者的5.5倍.     

上边界条件对多年冻土地温场数值模拟结果的影响分析

以玛多地区多年冻土为背景,建立多年冻土地温场的数值计算模型,以不同的方式考虑近 60 a 来的气温变化构成不同的上边界条件,通过模型计算分析不同上边界条件下的不同时期温度场、未来冻土退化特征.结果表明：在上边界条件中采用气象站实测近60 a波动温度值和采用近60 a平均恒定值时,浅层冻土地温差异明显,且越浅层地温与越近时间的上边界条件相关.预测未来100 a冻土地温变化趋势发现,相同升温速率和升温初始温度条件下,上边界采用实测60 a波动温度值对冻土退化过程影响较小;升温初始温度值提高到与趋势线衔接后,冻土退化起始时间从约第45年提前到约第20年;60 a实测温度和升温初始温度值均提高到与其初始温度场上边界条件衔接后,冻土退化起始时间从约第20年提前到约第15年;冻土退化从开始到完全退化经历时间为25 a左右.     

车载作用下箱形截面梁振动特性的瞬态动力学分析

论文根据磁浮交通道岔梁的结构及受力特点,通过建立合理的瞬态动力学数值模型,把列车作为高速运行的荷载动态作用于道岔梁上。采用瞬态拟动力学方法模拟钢结构梁动力响应特性。研究了道岔梁上关键点位置相应的振幅、反应加速度及其时效特性。通过计算比较精确地获得工程实际中最为关注的梁的最大振动位移和加速度。研究成果可为高速磁浮实际工程设计提供借鉴。     

大口径射电望远镜轮轨滚动接触有限元分析

通过显式有限元法建立轮轨滚动接触模型,用于分析110 m口径射电望远镜滚轮经过焊缝时的应力变化.建模时不仅考虑焊缝处材料的差异,同时还考虑轨道焊缝处表面不平度和摩擦系数的影响,分析了滚轮经过轨道焊缝时轮轨间作用力的瞬态响应.采用层级细分技术对轮轨模型进行网格划分,静力学分析结果表明模型求解精度较高.将网格划分模型用于瞬...     

基于混合边界条件的有限单元法GPR正演模拟

从Maxwell方程组出发,推导了探地雷达(GPR)有限元波动方程.阐述了透射边界条件和Sarma边界条件的原理,推导了这两种边界条件的理论公式;通过在衰减层内加入过渡带优化了Sarma边界条件的加载方法,压制了介质区和衰减层交界面处的人为反射.考虑到透射边界条件与Sarma边界条件不同的理论机制,提出了一种结合透射边界条件和Sarma边界条件的混合边界条件,它利用Sarma边界条件对到达边界区域的GPR波能量衰减功能和透射边界对GPR波能量的透射功能,使GPR波经过Sarma边界条件的衰减吸收后,再通过透射边界条件将剩余能量透射出去,集成了二者的优势.并以二维均匀模型中的中心脉冲激励源方式为例,通过Matlab程序实现,以GPR的全波场快照的直观方式,对比了有、无边界条件及不同边界条件对人工截断边界的处理效果,说明了该混合边界条件对到达截断边界处的GPR波的处理优于单一边界条件.最后,以基于混合边界条件的有限单元法对两个典型的GPR地电模型进行了正演模拟,指导了GPR数据处理与工程实践.     

用于地震观测的瞬态滤波研究

目前,地震数据采集器使用的最后两级数字滤波器相位有最小相位和线性相位两种,相位特性直接影响着记录的地震数据波形质量.最小相位滤波的优点是:地震波形没有“前缀”波形出现,不会干扰分析人员对震相的判读;其缺点是:通带边沿频率附近不同频率分量群时延差异较大,导致最小相位滤波会带来波形失真,对初动半周期的波形影响显著.     

基于GPU并行加速的叠前逆时偏移方法

为了提高复杂地下介质的成像精度和偏移算法的计算效率,提出可高效对地下复杂构造进行准确成像的GPU加速叠前逆时偏移方法.该方法采用双程声波方程进行波场延拓,突破倾角限制,借助于高阶有限差分方法实现叠前逆时偏移成像;利用GPU(Graphic Processing Unit)并行加速技术对波场延拓和成像进行计算,相比于传统算法,其计算效率有较大提高,可以解决叠前逆时偏移算法计算量过大问题;在获取波场信息过程中,也采用随机边界条件,实施以计算换存储策略,解决逆时偏移计算中的海量存储问题.模型测试结果表明,该方法能够高效和高精度地对地下复杂地质体成像.     

饱和土体瞬态响应有限元分析

给出基于Biot多孔介质理论分析饱和土体在动载荷作用下瞬态响应的有限元公式,数值计算部分采用本文有限元法分别计算一维饱和土柱在两种不同类型动载荷作用下的瞬态响应,并将数值计算结果与文献中的解析解进行比较,二者结果十分吻合,从而验证本文方法的可行性。     

浅析发动机冷却系统

一、现代发动机冷却系统的特点 传统冷却系统的作用是可靠地保护发动机,而还应具有改善燃料经济性和降低排放的作用。为此,现代冷却系统要综合考虑下面的因素：发动机内部的摩擦损失;冷却系统水泵的功率;燃烧边界条件,如燃烧室温度、充量密度、充量温度。     

多年冻土区路基调查中瞬态面波勘探方法应用研究

多年冻土地区铁路路基的振动与环境变化对路基稳定性的影响问题日益突出,为了评价这种影响引起的青藏铁路多年冻土区冻土的变化状况,首次应用了一种简便、高效、环保的瞬态面波勘探方法对其天然地表和典型结构路基及其下部多年冻土层进行了勘探。试验数据良好、频散特征显著,通过反演分析得到了青藏高原多年冻土地区天然地表和典型结构路基下部土体的剪切波速度剖面和多年冻土的分布。勘探结果表明,青藏铁路北麓河试验段勘探区内多年冻土上限一般为2 m,层厚为20 m左右,多年冻土层Vs平均为600 m/s以上。该结果与该区地质钻探编录和钻孔波速测试结果较为一致,说明瞬态面波勘探在自然环境严峻的青藏高原多年冻土地区有着良好的适用性。