考虑循环软化的非线性动力本构模型在FLAC3D中的实现

已有震害研究表明,震后边坡会因持续变形而破坏,且伴随着土体强度逐渐降低的现象,即土的循环软化行为。因此,有必要研究考虑循环软化的非线性动力本构模型以用于复杂条件下地震边坡稳定性分析。在已有的非线性动力本构模型基础上,提出了考虑循环软化的处理方法。同时,在FLAC3D平台上实现了本构模型二次开发,并通过了理论公式与已有文献中试验数据的验证。结果表明：计算出的骨干曲线与理论公式一致,且计算出的动剪切模量比及阻尼比与试验数据吻合较好,能够克服Hardin-Drnevich模型和Davidenkov模型在较大应变处（>0.01%）过高地估算阻尼比的缺陷;考虑了循环软化后,计算出的剪切强度有明显降低,且当遇到骨干曲线剪应力可以连续地过渡到软化后的主干曲线上,模型的收敛性较好。所开发的本构模型可为大应变条件下软土场地及边坡地震灾害评估提供支持。     

衡山花岗岩体西缘韧性剪切带的特征,成因及岩体定位机制

衡山花岗岩体由南岳岩体(早期)和白石峰岩体(晚期)组成,衡山岩体西缘形成有分带明显的韧性剪切带,向西倾斜,在XZ面上呈倾向滑动剪切运动,矿物显微变形构造显示为倾向简单剪切特征。经测定,韧性剪切带形成时温度下限为343.6℃,古差应力下限为52～62MPa,应变速率下限为2.05×10~(-14)/s-3.13×10~(-14)/s,推算总位移量在13.5km左右。在花岗岩体以外没有韧性剪切变形现象,说明该韧性剪切带是衡山岩体定位过程中由白石峰岩体沿南岳岩体边缘侵位的简单剪切作用所形成。     

二维圆柱坐标下FDTD法对多频电磁波测井的数值模拟

在介质围绕Z轴旋转对称的情况下,三维直角坐标下的麦克斯韦方程可变成两组分别对应TM波和TE波的二维柱坐标下的偏微分方程.采用线圈作为电磁波测井的发射器和接收器的时候,可采用TE波对应的方程,这时相当于磁极子的情况.而采用电偶极子作为发射和接收器的时候,则对应TM波的方程.采用时间和空间均为二阶精确的有限差分方法,将偏微分方程进行差分化.这样,空间的电磁场可由时间域有限差分法(FDTD)来求解.     

基于自适应加权广义逆矢量方向滤波估计地震同相轴倾角

在断层等不连续位置处估计局部地震同相轴的方向较为困难。在断层、地质边界等处,同相轴不连续,结构不稳定,导致倾角等地震属性变化大。为了克服这个问题,本文提出一种自适应加权广义逆矢量方向滤波方法,通过光滑经过逆转的梯度矢量场来估计局部地震倾角剖面。计算过程是：首先由改进的梯度算子计算梯度矢量场;然后设置垂直向下为参考方向,对梯度矢量场进行逆转;最后通过自适应加权广义逆矢量方向滤波方法估计倾角剖面。自适应加权广义逆矢量方向滤波方法在不连续位置处具有传统倾角估计方法不具有的高分辨率,能够稳健高效地估计地震数据的倾角剖面。理论模型和实际地震资料处理得到的局部地震倾角剖面表明,所提出的方法能够高分辨率地保持构造细节信息。     

高阶交错网格有限差分弹性波场模拟的精度分析

交错网格波场数值模拟是目前地震正演中广泛使用的方法,为对比分析不同阶数的差分格式下产生的计算效率和精度差异,重新推导了弹性波方程的4种时间4阶、空间2N阶的差分公式及系数,并计算了他们的稳定性条件。利用这4种差分格式进行弹性波场数值模拟,对比分析了波场快照、合成地震记录及CPU时间。结果表明：时间4阶、空间6+6阶精度的交错网格有限差分方法在进行地震波场数值模拟时具有较高的计算精度和计算效率。     

地勘单位经济指标考核体系探讨

地勘单位执行多年的经济指标考核体系一般由决算总收入、新增资本积累、净资产总额、应收账款与存货占经营收入比率、提取折旧额、提交煤炭资源量,人均年收入等七项指标组成。该指标体系经过多年运行,已经趋于成熟。但是该指标考核体系没有考虑到不同单位的综合经济实力和生产能力水平的差异,忽略了不同单位之间的主营业务性质和区域性差异,使考核结果公平性存在一定的缺陷。本文依据全面考核原则、科学考核原则、可操作性原则,探讨了如何改进经济指标考核体系,并提出了可供实施的企业经济效益综合指数考核指标体系。     

基于CFD的大通孔开闭式钻头流场研究

采用计算流体力学方法对大通孔开闭式钻头\[1\]的流场进行数值模拟。对不同工况下的速度场进行分析,获得了最佳钻进风量,并发现流场中存在“阻塞漩涡”。对最佳流量工况下的压力场进行分析,发现流场中还存在“抽吸压差”。获得了“阻塞漩涡”和“抽吸压差”对排渣和散热的影响机理,为钻头的优化设计提供了一些参考依据,为松软突出煤层全孔筛管护孔瓦斯抽采技术在煤矿中推广应用奠定了基础。     

考虑起始坡降双层地基的大应变非线性固结

天然软土成层分布特性及土中渗流存在起始水力坡降的现象已被人们熟知。但变荷载下能同时考虑黏土中起始水力坡降、软土非线性压缩渗透特性及大应变特性的双层地基固结理论还鲜见报道。在拉格朗日坐标系中建立以超静孔压为变量的双层软土地基大应变非线性固结模型并给出其有限差分解。通过与考虑起始水力坡降的单层地基大应变非线性固结数值计算结果对比,验证了差分解的可靠性。着重分析了上、下土层起始坡降无量纲参数R1、R2对双层地基固结性状的影响,分析在大应变与小应变假定下双层地基超静孔压消散及固结沉降变形的异同。结果表明：上层土无量纲参数R1对双层地基固结性状的影响程度较下层土无量纲参数R2显著;大应变假定下双层地基渗流前锋的下移速度要快于小应变假定下的移动速度;大应变假定下考虑起始水力坡降的双层软土地基超静孔压消散速率要比小应变假定下快,且大应变假定下考虑起始水力坡降的双层地基最终沉降量要比小应变假定下大。     

珠江三角洲软土微结构的定量研究

对用不同长度的PVC管在珠海和广州两地采得的淤泥样柱施于一定荷载，然后在扫描电镜下研究软土受压后在样柱垂深方向上的微结构变化，并将样品的微结构照片扫描，在计算机内进行量化处理。结果表明，样品的孔隙分布分维差随着离施压点距离的增加而逐渐减小，从而得到两者的关系方程。据此，可了解在一定荷载下应力在软土层中向下的衰减趋势。文章还对软土的天然含水量和围压对压缩变形的影响作了初步的讨论。     

一阶弹性波交错网格时间高阶差分格式及稳定性分析

弹性波模拟或逆时偏移时,对空间偏导数采用高阶差分格式可提高计算精度,但这种算法的稳定性条件过于严格,要求差分离散的时间步长必须足够小以确保算法稳定。在常规空间高阶差分格式的基础上,将速度(应力)对时间的高阶导数转化为不同精度的应力(速度)对空间的差分,得到了一种新的基于交错网格的时间高阶、空间高阶差分格式。通过对交错网格时间高阶差分格式稳定性的分析,认为该算法的稳定性条件较常规算法宽松,在弹性波场的求解过程中可以采用更大的时间步长。