波—流在摩阻地形上的综合传播

本文采用波浪缓坡方程抛物型近似方法，研究了缓变地形上波－流共存情况下的传播问题。文章首先提出了包含摩阻项在内的抛物型控制方程，然后结合波－流作用守恒方程和文献（８）的成果，给出了波能衰减率的简便形式。算例表明，本文提出的计算方法可用于大型水域实际工程设计中波浪要素的计算。     

水锤计算中摩阻的处理方法

总结瞬变流控制方程摩阻项积分的处理方法,指出这种在传统的拟恒定摩阻假设基础上的处理方法存在问题,着重阐述国外对瞬变流的非恒定摩阻问题的研究进展,分析评述了几类典型的非恒定摩阻模型,并指出一些新思路和研究方向.     

水锤计算中摩阻的处理方法综述

管道流体摩擦引起两个问题：数值计算上摩阻项的(显)隐式有限差分近似；(非恒定)摩阻引起的水锤波衰减。后者一直未得到合理解释。本文先总结传统的在拟恒定假设基础上对摩阻的处理．再着重阐述国外对瞬变流下的非恒定摩阻问题的研究进展，并指出一些新思路和研究方向。     

倾斜底坡上波浪的传播与破碎

本文根据TUMMAC方法探讨了孤立波、微幅波和二阶Stokes波在1/30底坡上的传播与破碎,并采用示踪法测得了相同底坡上濒临破碎波的速度场。与前人的实验结果相比较表明,计算所得的波浪外形变化是一致的,依据运动学标准所得的各项破碎指标也大致相同;与岩垣及本文实验结果比较表明,速度的分布趋势也是一致的。     

基于EIT摩阻模型与随机摩阻模型的水击随机分析

在以往水击随机分析中,摩阻按确定量考虑,其影响往往被忽略;而实际的水击过程中,摩阻为非恒定值,为一不确定量,文中研究非恒定摩阻对随机水击波动过程的影响.从水击的基本方程出发,分别建立了基于广义不可逆热力学理论的水击随机模型(EIT摩阻模型)与含有随机摩阻项的水击随机模型(随机摩阻模型),采用矩特征线法对水击过程进行数值模拟.算例表明,与采用恒定摩阻的水击随机分析比较,本文提出的两种随机模型在计算最大水击压力值的同时,可以模拟水击压力的波动过程,而且随机摩阻模型的数值计算更为简便.     

平面潮流数值模拟中底床摩阻系数的修正

本文从三维潮流运动控制方程出发，假定潮流垂直分布服从对数律，推得了一个较全面反映水深、水位变化梯度和流向等因素的摩阻系数修正表达式，理论分析和应用表明，运用修正后的摩阻系数，平面二维潮流数值模型更为完善，并能适用于陡变地形，从而使复演模拟和工程计算结果更为科学合理。     

沙波形态对波浪要素的响应

针对波浪作用下沙质底床上普遍存在的沙波现象,从沙波尺度角度出发,采用1957—2018年共385组野外实测资料和水槽试验资料进行拟合分析,得到基于近底波浪水质点运动强度的沙波波高、沙波波长计算公式。结果表明:与Ni81规则波、Ni81不规则波、VR89和GK04等4种典型模型对比,新的沙波波高和沙波波长计算公式具有更高的计算精度且表达形式上更具一般性;利用该沙波形态计算公式可改进床面粗糙度的计算方法,当应用于波浪底摩阻系数的求解时计算结果与实测数据较为吻合,即沙波形态对底边界层有效粗糙度的影响显著;新的沙波波高和沙波波长计算公式可较好地应用于波浪底摩阻系数的计算。     

考虑材料特性E,ν不确定的随机边界元法

本文把边界元的灵敏度分析与摄动法相结合，提出了随机边界元法，利用这种方法，可以计算结构响应的均值和标准差，估计由于材料特性Ｅ，ν的不确定性导致的响应误差。     

斜向波与长棱柱型浮式结构物相互作用的边界元模拟

采用第二类变形Bessel函数对斜向线性入射波与任意截面形状的长棱柱型浮式海洋结构物的相互作用问题进行了边界元法模拟,给出了常数元和线性元的奇异积分处理方法以及开边界的处理方式。对矩形和圆形结构物的水动力系数和波浪力进行了计算,并与解析解及有限元计算结果进行了比较。结果表明,本文方法能得到满意的结果,且未出现不规则频率现象。     

An iterative Rankine boundary element method for wave diffraction of a ship with forward speed

A 3-D time-domain seakeeping analysis tool has been newly developed by using a higher-order boundary element method with the Rankine source as the kernel function. An iterative time-marching scheme for updating both kinematic and dynamic free-surface boundary conditions is adopted for achieving numerical accuracy and stability. A rectangular computational domain moving with the mean speed of ship is introduced. A damping beach at the outer portion of the truncated free surface is installed for satisfying the radiation condition. After numerical convergence checked, the diffraction unsteady problem of a Wigley hull traveling with a constant forward speed in waves is studied. Extensive results including wave exciting forces, wave patterns and pressure distributions on the hull are presented to validate the efficiency and accuracy of the proposed 3-D time-domain iterative Rankine BEM approach. Computed results are compared to be in good agreement with the corresponding experimental data and other published numerical solutions.     

不平水底浅水波问题的位移法

该文研究考虑不平水底的二维浅水波问题的位移法。利用水平位移与竖向坐标无关的假定,建立了基于位移的浅水波方程和相应的Hamilton变分原理。针对文中提出的基于位移的浅水波方程,利用有限元和保辛算法建立了相应的计算格式。最后给出一个具有抛物线型河床的浅水周期晃荡算例。数值算例表明该方法可以有效模拟带动边界和不平水底浅水波问题。     

三维裂隙网络渗流计算的边界元法及程序

岩体力学性质受结构面控制，表现出强烈的各向异性。同时工程地质师多年的实践认为岩石结构面是具有统计性的，因此岩体结构面网络模拟成为研究岩石力性质的新方向。其中裂隙网络渗流模拟，特别是三维裂隙网络的渗流计算是结构面网络模拟研究工作的重点和难点。又由于对于具有随机性裂隙网络渗流的计算，需要多次形成裂隙网络样本，并对每个样本进行多次计算，因此对渗流计算比确定性渗流提出更高的要求。本文基于将单个裂隙看作为一个不可压缩的二维各向同性多孔介质的假定，运用加权剩余法，推导了单裂隙渗流的边界单元法公式。通过将裂隙之间的交线看作为内部边界，推导了多裂隙相交时的边界单元法公式。利用FORTRAN 90语言编制了相应的计算程序BIEMNF。为验证程序的正确性，利用计算实例进行了考核。     

白石水库泄水排沙底孔三维有限元静力分析

按水库运用的3种工况，采用三维有限元法，分析了在各种工况下泄水排沙底孔的应力分布规律，为优化底孔布置和结构设计提供了科学依据，对深入了解底孔运行后的工作性态有重要意义。     

坝基渗流的边界元分析

采用边界元法分析了坝基渗流问题,并与有限元法和裂隙网络法的计算结果进行了比较。结果表明,用边界元法处理同类问题,具有方便、快速、精度高等优点。     

用守恒元和解元法数值模拟二维溃坝洪水波

本文采用守恒元和解元方法（简称CE／SE法）来求解二维浅水方程,该方法从守恒积分型方程出发,严格遵守物理意义上的守恒律,通过设立守恒元和解元,从而提出一种新的数值模拟方法。用该方法计算得到的全溃的溃坝洪水波的数值解和理论解与文献中的实验值符合的较好。在此基础上,进一步对二维局部溃坝洪水波的演进过程进行模拟,结果表明,CE／SE法精度较高,稳定性较好,其溃坝洪水波的间断具有很高的分辨率,能很好的模拟溃坝洪水波的运动过程。     

等效黏性边界单元及黏性边界波动输入方法

为更好的利用大型通用软件进行结构-地基动力相互作用分析,本文在黏性边界的基础上提出黏性边界单元的概念,通过构造相同阻尼矩阵的方式提出了能够实现黏性边界功效的特殊有限单元-等效黏性边界单元.等效黏性边界单元能够较好的模拟地基的辐射阻尼效应,使黏性边界的实现更加方便.同时,为简化大型通用软件中地震波的输入,提出了适合于黏性边界条件下水平成层地基中垂直入射地震波的输入方法.该方法不需要借助其它辅助计算就可直接在各种大型通用有限元软件中用一个数值模型实现地震波动输入,具有很强的通用性.此输入方法极大的减少了地震波动输入的数据准备工作,降低了操作的复杂程度.数值算例表明,等效黏性边界单元及文中提出的地震波动输入是简单可靠的.     

波浪力计算中的一个新边界元方法

本文就波浪与结构物作用问题，提出了一个新的高阶边界元方法。较之Ｅａｔｏｃｋ Ｔｏｙｌｏｒ和Ｃｈａｕ的方法，本文在体内附助积分方程中应用了满足刚性水面条件和弱刚性底面条件的简单格林函数。计算速度相应地有所提高；漂浮结构的积分方程更为简洁；计算结果与Ｅａｔｏｃｋ Ｔａｙｌｏｒ和Ｃｈａｕ的方法具有几乎相同的精度和收敛速度。     

三角门底枢的非线性接触有限元分析

底枢的疲劳破坏和磨损是船闸闸门经常维修的主要原因之一,底枢在设计过程中进行了非线性接触分析,检验其结构强度并针对结构设计提出改进措施。以滁河汊河集船闸下游三角门为研究对象,详细介绍了Workbench对底枢非线性接触有限元分析的思路和方法,对比分析底枢在正反向水压力作用下的变形和接触应力情况。结果表明:底枢整体结构满足规范要求;底枢接触面在反向水压力作用下更易失效;局部应力集中部位需要设置圆弧倒角或过渡圆角,避免结构上产生较大的尺寸突变部位。该研究为底枢结构设计提供参考依据。