一种基于全局优化的交错网格有限差分法

在地震波场数值模拟中, 交错网格有限差分技术得到了广泛的应用, 但是在弹性模量变化较大时, 通常会因插值而导致模拟误差增大. 旋转交错网格可以很好地克服这个缺点, 因而适合于各向异性介质正演模拟. 但是对于同样大小的网格单元, 旋转交错网格需要的步长比常规交错网格要大, 这会使梯度和散度算子的误差增大因而更易产生空间数值频散. 针对这些问题, 本文提出了旋转交错网格与紧致有限差分相结合的方法, 并基于模拟退火算法进行全局优化, 压制数值频散, 拓宽波数范围. 数值模拟结果表明, 此方法可以有效地压制数值频散, 且具有较高的模拟精度.      

近海风机单桩基础桩周土应力特征分析

桩周土体的力学响应规律对近海风机桩基础的变形和稳定性有重要影响。采用有限元方法对一个概化的近海风机单桩基础进行数值模拟,着重分析动荷载作用下桩基础和桩周土的变形规律、桩周土应力作用特征及应力主轴的旋转模式。发现桩周土体各方向的应力分量随循环荷载的作用呈现出周期性的往返变化,第一周期的应力路径与继后周期的应力路径明显不同,浅层桩周土存在显著的应力主轴往复旋转,旋转幅度可达60°。     

由书架型断层作用容纳的冰岛北部火山裂谷带的运动

沿中央海岭,延伸的地壳以10~1 000km的尺度分段（MacDonald et al,1988）。在裂谷段相互错开的地方,段与段之间的运动由垂直于主裂谷轴方向的转换断层所容纳。在段叠加的地方,不同几何形状的非转换断层（Grindlay et al,1991）容纳剪切运动。这里,我们采用微震数据来分析冰岛北部出露的两个叠加裂谷段的断层几何形状。在这两个裂谷段之间,我们识别了沿主裂谷走向排列的一系列断层。这些断层为左旋走滑运动。而叠加裂谷段之间的运动则是右旋运动。这些运动一起产生了断层的顺时针旋转,致使地壳块体按书架型断层机制（类似于书架上倾斜放置的书的排列而得名）运动（Mandl,1987）。这些断层或许活化了已经存在的地壳软弱面,如原始走向平行于主裂谷并具有15°顺时针旋转的岩脉侵入体。先期存在、平行于裂谷的软弱面的活化与典型的中央海岭转换断层相反,说明非转换断层容纳着叠加段之间的剪切错动。     

扇束CT探测器偏转的重建图像伪影校正

对于CT扇束扫描,射线源焦点与被测工件旋转中心的连线不垂直于探测器所在直线时,采用标准的滤波反投影重建算法重建出的图像会存在伪影,降低图像质量。本文基于经典滤波反投影重建算法理论,分析了此类伪影形成的原因,给出了利用投影原始数据（正弦图）计算校正算法所必需的几何参数的方法。本方法无需专用模体,不必进行二次扫描,方便、快捷并有效。计算机仿真及实验结果证明本方法的有效性。     

一种用于中国年最高(低)气温区划的新的聚类方法

采用聚类分析和旋转主分量分析相结合的方案,对我国年最高（低）气温的年际变化型态进行地理区划。这种两者相结合的分区方案可以互相补充,使区划更具客观性。结果表明,中国年最高气温和最低气温年际变化分别可划为12和11个不同类型的区域。对前者（高温）来说,各区域的年际差异比较大,气温变化的阶段性特征有所不同;对后者（低温）来说,各区域具有较强的增温一致性,但增温的特点和幅度存在明显区域差异。最高（低）气温和平均温度的关系在各区域是不同的。总体来说,年代际变化的一致性要好于年际变化的一致性。     

铁伊利石近双晶的电子衍射图

本文研究的伊利石采自天津蓟县中元古界长城系串岭沟组。EDX数据显示该样品富Fe贫Mg,这与XRD分析一致。其化学分子式为K0.59(Al1.43Fe0.83)(Si3.47Al0.53)4O10(OH)2,为铁伊利石。由高分辨透射电镜H-8100和H-9000获得了奇特的衍射花样。该电子衍射花样是由相互旋转了26°的两个h     

黏弹TTI介质中旋转交错网格高阶有限差分数值模拟

以Carcione黏弹各向异性理论为基础,给出了适用于黏弹性具有任意倾斜对称轴横向各向同性介质(黏弹TTI介质)的二维三分量一阶速度-应力方程,采用旋转交错网格任意偶数阶精度有限差分格式求解该方程,并推导出了二维黏弹TTI介质完全匹配层(PML)吸收边界条件公式和相应的旋转交错网格任意偶数阶精度有限差分格式,实现了该类介质的地震波场数值模拟.数值模拟结果表明:该方法模拟精度高,边界吸收效果好,可以得到高精度的波场快照和合成记录;并且波场快照和合成记录能较好地反映地下介质的各向异性特征和黏弹性特征.     

斜压海洋能量传递对风场的响应

风是海水运动的重要动力因素,也是海洋内部的主要能量来源.本文在应用陆架海洋模式HAMSOM对东中国海海水运动进行数值模拟的基础上,通过傅里叶变换、旋转谱分析等研究方法,对风向海洋的能量输运进行研究.研究结果显示,风场向海洋输运的最有效途径是风杨扰动量与流场扰动量的相互作用;惯性能量主要来源于海洋表层,由风场提供,向下传递;而潮频率能量大部分来自海底的内潮与底地形相互作用,向上传递.     

旋转触探数学模型及其试验分析

通过对旋转触探用双螺旋探头的受力分析,建立了螺旋探头的贯入阻力和扭矩与贯入速度和转速,以及岩土的力学性能参数间的数学模型。结合数学模型和试验数据对旋转触探技术的影响因素进行了分析,为设计探头结构和进一步理解该技术的力学机理提供了理论依据。      

利用部分投影的二维CT旋转中心偏移快速校正

现有的CT重建算法都依赖于转台旋转中心的准确定位,否则会引起重建图像的伪影。目前已有的旋转中心估计方法的共同特点是利用180°或360°的投影数据。获取这些投影,需要比较长的扫描时间,占用比较多的资源。为此,在平行束扫描和扇形束扫描两种情形下,本文提出利用部分视角的投影来对转台旋转中心进行快速定位的方法。在平行束情形,利用物体空间质心与投影空间的质心的对应关系定理,由60°的投影即可准确测量出旋转中心的位置。在扇形束情形,通过将扇形束投影重排成平行束投影,仍可由扇形角再进行60°的投影准确测量出旋转中心的位置。