浅水方程数值计算方法的研究

求解以水位为变量的连续方程,并根据Navier-Stokes方程压力修正算法的基本思想,建立了浅水方程的水位修正算法,放宽了对离散时间步长的限制.通过对离散方程系数矩阵的重新构造,建立了高分辨率有限元格式,该格式既具有较高的离散精度又避免了数值解的伪振荡.对动量方程的阻力项做负坡线性化处理,提高了露滩计算的稳定性.数值模拟结果与解析解吻合良好,表明所建立的数值计算方法是正确的和可靠的.     

基于应变模态法识别刚架桥梁的损伤

基于应变模态方法，对刚架桥的损伤识别进行了研究。通过对某刚架桥在不同损伤工况下的数值仿真计算，探讨了应变模态方法用于刚架桥损伤识别的能力。计算结果表明：利用应变模态差曲线能比较准确地识别出刚架桥的损伤位置；应变模态差曲线在刚架桥损伤单元处的跳跃幅值随单元损伤程度的增加而增大，依此可定性地识别出刚架桥的损伤程度。     

锚杆动测问题的数值解

建立了锚杆动测问题的有限元模型，并对锚杆中的缺陷进行了数值模拟，获得了位移场的数值解。对该解在时域和频域范围内的讨论表明，锚杆中的缺陷将引起动测曲线一系列的变化，为现场动测曲线的识别提供了依据。     

深厚软弱地基上相邻桩箱基础高层建筑地震反应的数值模拟

采用二维有限元整体分析法对深厚软弱地基上相邻的双幢桩箱基础高层建筑的地震反应进行了数值模拟，探讨了SSI效应对双幢高层建筑地震反应的影响规律；并分析了深厚软弱地基上双幢高层建筑和单幢高层建筑地震反应之间的差异。一般地，SSI效应使深厚软弱地基上相邻的双幢高层建筑各楼层的绝对加速度反应要比单幢高层建筑的反应略小一些，但也可能使双幢高层建筑中部楼层的绝对加速度反应大于单幢高层建筑的反应。与此相反，SSI效应一般使相邻的双幢高层建筑各楼层的相对位移反应有所增大；但也可能使双幢高层建筑各楼层的相对位移反应有所减小，两者的差异与场地条件、输入地震动特性有关。此外，还探讨了双幢高层建筑的间距对其地震反应的影响。随着间距的增大，双幢高层建筑的地震反应趋近于单幢高层建筑的反应。输入加速度峰值越大，场地土越软，相邻的双幢高层建筑的间距对其地震反应的影响越大。     

长江三角洲稻田生态系统综合增温潜势源汇交替的数值分析

将生物地球化学C, N循环模式与植被生态模式耦合, 建立了生物地球化学-植被生态耦合动力模式. 验证表明模型可以准确模拟稻田生态系统中作物生长及CH₄, N₂O排放过程. 利用模式分析了不同施肥量对CO₂吸收及CH₄, N₂O排放的影响, 以及3种温室气体综合增温潜势的演变规律, 发现了生长期内稻田综合增温潜势源汇交替现象, 提出最大汇施肥量以及零排放施肥量等概念. 利用模式研制了长江三角洲地区不同施肥方案, 为该地区可持续发展提供了科学依据.     

运用三维变密度潮汐效应模型确定滨海含水系统的海底等效边界——山东烟台夹河中下游地区为例

确定咸淡水界面的位置是滨海地区海水入侵研究的主要任务之一。对于天然条件承压含水层而言，含水层顶板向海底延伸的距离直接影响了咸淡水界面的位置，它可以通过承压含水层中地下水的潮汐效应信息来确定。考虑到咸淡水之间密度的差异，建立了山东省夹河中下游地区滨海含水系统地下水三维变密度潮汐效应模型。通过反复对比潮汐效应观测中的地下水水头波动与模型计算出的水头波动，确定了滨海承压含水系统的海底边界。同时，也初步估计出海区与近海陆区含水层的水文地质参数。     

具有脆性形变和蠕动变形的复合单断层系模型中的应力、滑动与地震

对弹性半空间中的几个含有不连续断层系的三维弹性连续介质模型，在沿单元垂直走滑断层发生的滑动过程进行了数值模拟，模型的几何形状和所使用的边界条件大致与圣安德烈斯断层中部的情形相符合，数值模拟综合考虑了脆性形变和蠕动变形，断层带的总形变率为蠕变和摩擦滑动率之和。脆性断层的性质由破裂段（数值单元）的不同地震应力降的分布来给定。所假设的分布描述了两种理想化的情形，它们分别对应于不同断层滑动的极限状态：（1）以大尺度范围为特征的强无序状态，表示不成熟的断层带及扩展空间域。（2）相对有序的状态，以小尺度范围为特征，表示成熟的完全滑动断层。假定断层的蠕变特性在所有的情形中都相同。并且是用蠕滑速率与应力的幂次律关系来刻画的。利用这种脆性蠕变过程和模型参数可以得到一个在12．5km深度含有“脆性－延性”过渡带的应力－深度剖面，以及一个在1857年地震破裂西北65km处存在有“脆性－蠕变”过渡带的沿走向的应力剖面。模拟的震源空间分布在统计意义上与观测数据相类似，结果表明，断层强烈不均匀性的特征尺度范围会在断层系统的地球地震响应中显著表现出来，若尺度范围小，则会使震级－频率分布接近于特征地震分布，在时间上会层致如地震空区说所预言的大地震的准周期分布，而大尺度范围的地震活动性符合古登堡－里克特的震级－频度关系以及大地震的随机或是成簇的时间分布特征。数值模拟的结果表明，各种形式的震级－频度关系以及地震时间分布的统计特征都可以通过表征断层带不均匀性的尺度范围统一起来，这可以通过借助一个给定的断层带或是岩石层的结构特征给出一个清晰的物理解释，并且得到地震及断层观测资料的支持。在一些模拟中，小地震的震级－频度关系曲线明显低于自的古登堡－里克特关系曲线。这说明在构造加载过程中，小地震通过平滑应力的长波成分，为大地震的发生创造了条件。这一过程是通过与小地震的大量断裂有关的短波应力的消失过程来完成的，大地震的大尺度破裂的发生正好导致与上述长，短波应力的平滑与消失变化相反的过程。     

哈尔滨市磨盘山水利枢纽区三维渗流数值模拟模型研究

哈尔滨市磨盘山水库防渗型式拟采用垂直砼防渗墙及帷幕灌浆方案。为了分析左岸不同长度防渗对水库渗漏量及坝基渗透稳定性的影响。对磨盘山水利枢纽区渗流场进行了三维数值模拟研究。结果表明，随着左岸防渗长度的增大。渗漏量及下游渗流出口渗透比降基本呈线性趋势减小。当左岸防渗长度为200m时。渗漏量及各渗透比降均小于允许值，为推荐的渗控方案。     

最优化的矿物扩散地质速率计:北大别黄土岭麻粒岩缓慢冷却的证据

对北大别黄土岭麻粒岩中石榴子石进行了详细的电子探针工作，并根据一颗石榴子石的成分环带，建立了扩散模型并进行了数值模拟，结合最优化的方法，得到了极其缓慢的冷却速率0.1℃/a。从扩散动力学数值模拟的角度对大别山造山带的冷却速率进行了探索性研究。该结果表明，黄土岭麻粒岩在受麻粒岩相改造后经历的是一个极其缓慢、持续时间很长的冷却过程。     

地下水及溶质运移数值模拟系统

介绍了地下水及溶质运移数值模拟的研究进展、数值模型的分类和数值模拟的步骤。根据数值模拟中涉及到的地质信息资源的特点，提出了基于地理信息系统(GIS)技术的地下水及溶质运移数值模拟系统的设计思路。该系统能够对模拟区域的基础地质资料、水文地质勘探资料以及数值模型的输入、输出资料进行动态管理，而且随着勘探、科研和生产的不断进行，通过对数据库的及时更新以及对数值模型的不断校正，可为不同的科研和生产方案的预测与评价提供信息资源。