用面元法预报船舶的兴波阻力及波型

采用基于开尔文源的面元法求解航行于均匀流场中的船舶兴波问题．在船体湿表面上布置开尔文源，使其严格满足船体表面条件．数值计算中，首先采用分部积分对被积函数进行简化处理，然后采用高斯积分实现面元和线元上的积分，避免了被积函数为高频振荡函数所带来的数值计算的复杂性和不确定性．计算结果表明，其波型及兴波阻力系数与试验结果和其他数值方法相比均吻合较好．     

基于波叠加原理的辐射声场的计算研究

描述了用波叠加原理来计算任意形状的辐射体的辐射声场的方法。即用叠加来代替Kirchhoff-Helmholtz积分的方法．给定辐射体表面的振动速度就可求解辐射源的强度。进而求解辐射体表面的声压和辐射的声功率．文中对刚性球面上振动活塞源和脉动球源的辐射问题进行了求解。并将之与分析(解析)解进行了对比，结果表明用该方法进行计算时，可以用简单的数字方法产生矩阵元素，从而提高计算速度。结果的精度高。并能克服边界元法声场计算内在的奇异性问题．     

大深度水下航行体引起的水底压力变化

水下航行体引起的水底压力变化是航行体的重要目标特性之一,可用于水下运动目标探测和水中兵器研制.基于船舶水动力学势流理论,建立了水下航行体水底压力计算模型.利用HessSmith方法,对大深度水下航行体的压力场分布进行了计算,分析了航行体水压场沿纵向和潜深的变化规律,计算结果和实验结果有较好的一致性.     

螺旋桨非定常空泡数值预报中的几个问题

用基于速度势的低阶面元法对螺旋桨非定常空泡进行了预报分析，并对数值分析中的几个重点问题进行了分析讨论．通过桨叶表面划分不同的网格数进行计算，比较不同网格划分计算得到的空泡形状，讨论了空泡计算对网格疏密程度的要求．针对非定常空泡计算中不容易收敛的问题，采取的数值方法既能够较好地反映空泡的非定常性，也能够较快收敛，使得空泡数值模拟更加符合实际情况．计算得到空泡形状后，利用计算得到的奇点强度分布计算了空泡表面的压力分布，验证了空泡表面的压力等于水的汽化压力．     

基于NURBS高阶面元法的浅水船舶兴波阻力计算

实现了一种船舶兴波阻力和兴波波形的数值计算方法，将CFD与CAD技术联系在一起，应用NURBS曲面造型技术将船体曲面表达为统一的数学模型，使性能计算的真实物面条件得以满足，应用高阶面元法数值模拟单、双体船在浅水航道中定常直航状态下的兴波阻力与兴波波形，并与瑞典ShipFlow软件的计算结果进行比较分析，表明该方法是可行的．     

Rankine源法的船舶下沉量数值计算初步研究

船舶下沉量是影响船舶安全富余水深的重要因素,具有一定的研究价值.本文对基于Rankine源法的船舶下沉量数值计算进行了初步研究,该方法是一种基于势流理论的边界元法.使用Rankine源法计算船舶下沉量时,首先在船体、自由面、岸壁等边界上分布一系列Rankine源面;进而对边界条件线性化,将源强带入边界条件,建立源强的线性方程组;用迭代法求解该线性方程组,获得源强值,根据源强值计算出船体水动力、力矩和下沉量.研究结果表明,该方法计算船舶下沉量精度较高.     

穿透自由面平板横荡和道水动力计算

介绍了一处计算穿透自由在垂直平板横这运动水动力的三维Ｒａｎｋｉｎｅ面元法，采用在平板平均上及其尾流中的涡强谐变的涡分布，同时，彩在围绕平板的自由面上 强谐变的源分布以计人自由面影响，将所介绍的方法应用于小展弦比垂直平板，给出了计算结果并对结果进行了讨论。     

非线性兴波阻力计算

叙述了计算船舶非线性兴波阻力的一种方法，采用面元法满足完全非线性边条的势流场。特点是采用三角形单元，源面自自由面上置和自湿水面下潜，每次迭代均适应新的波面重新作网格划分，所采用的迭代方式未发生不收敛问题，只需少数迭代次数即达到收敛。对Ｔｏｄｄ６０系列Ｃｂ＝０．６￣０．７５及不同Ｌ／Ｂ，Ｂ／Ｔ作了计算，所得计算结果与试验值符合，令人满意。     

交通分布-交通分配组合模型研究

在交通规划实践中，交通分布和交通分配问题作为“四阶段预测模型”各自独立，尽管两个问题紧密相关．因此，有必要建立交通分布-交通分配组合模型，来反映这种反馈．文中介绍了Evans模型，并说明了逐次平均法求解模型的步骤，最后用算例说明了具体算法．这种方法完全可以用Microsoft Excel计算，计算方法简单，可以为组合模型的求解提供新的思路．     

���������ij������������ģ�����㷨�о�

在对多配送中心车辆调度问题进行直观描述的基础上，建立了该问题的数学模型。提出了采用距离最近分配法将多配送中心车辆调度问题分解为多个单配送中心车辆调度问题进行求解的策略．基于求解单配送中心车辆调度问题的禁忌搜索算法，设计了求解多配送中心车辆调度问题的算法，并进行了实验计算．计算结果表明，用本文设计的算法求解多配送中心车辆调度问题，不仅可以取得很好的计算结果，而且算法的计算效率较高，收敛速度较快，计算结果也较稳定．     

Prediction of hydrodynamic forces on a moored ship induced by a passing ship in shallow water using a high-order panel method

A three-dimensional high-order panel method based on non-uniform rational B-spline (NURBS) is developed for predicting the hydrodynamic interaction forces on a moored ship induced by a passing ship in shallow water. An NURBS surface is used to precisely represent the hull geometry. Velocity potential on the hull surface is described by B-spline after the source density distribution on the boundary surface is determined. A collocation approach is applied to the boundary integral equation discretization. Under the assumption of low passing speed, the effect of free surface elevation is neglected in the numerical calculation, and infinite image method is used to deal with the finite water depth effect. The time stepping method is used to solve the velocity potential at each time step. Detailed convergence study with respect to time step, panel size and Green function is undertaken. The present results of hydrodynamic forces are compared with those obtained by slender-body theory to show the validity of the proposed numerical method. Calculations are conducted for different water depths and lateral distances between ships, and the detail results are presented to demonstrate the effects of these factors.     

移动荷载与土体中孔洞相互作用的2.5D边界元法分析

为了研究具有复杂表面的饱和土体与移动荷载相互作用机理,根据饱和土Biot理论,采用Fourier变换和势函数分解法,推导了饱和土体频域-波数内的Green函数;建立了移动荷载作用下,饱和土体中孔洞动力响应的频域-波数域内边界元法方程(2.5D boundary element methods);利用快速Fourier逆变换法,得到了时间-空间域内饱和土体的动力响应.研究结果表明:在频域内,利用移动荷载方向的一致性建立的频域-波数域内边界积分方程,可将3D空间问题转化为频域-波数域的2D平面问题, 3D边界元简化为2D边界元,使得计算面转换为计算线,减小了计算规模.      

水下悬浮隧道管段结构流阻特性分析

为探讨结构断面形状对悬浮隧道管段表面压力分布的影响,采用ANSYS流体模块RNG(renormalizationgroup)k-ε湍流模型,对圆形、多边形、曲边形和椭圆形4种断面的悬浮隧道管段海水绕流场进行了数值模拟,分析了绕流场及管段表面的压力分布特点.结果表明,同种断面悬浮隧道管段表面的压力分布规律相同,压力值随来流速度变化明显;管段曲率是影响绕流场的重要参数,海流在流线形断面上分离晚,表面压力下降较慢,降低了管段上的海流阻力.     

Computationally efficient approaches to fatigue analysis of deepwater risers

Riser long-term fatigue performance is an important design consideration. Although extensive application of irregular sea analysis in time domain with the rainflow counting technique for post-processing is regarded as the most accurate of the approaches for fatigue analysis, it does suffer from some limitations, such as the computational effort. For this reason, two computationally efficient approaches are employed to perform the fatigue analysis of a deepwater top-tensioned riser, based on the Longuet-Higgins distribution and time domain scaling respectively. With Longuet-Higgins distribution irregular wave sea states are expanded into their individual wave bins. These regular wave simulations are of short duration and consequently run quickly. Using the time domain scaling technique, the number of irregular wave runs can be performed for a comparatively small number of load cases and hence reduces the calculation time. The results showed a reasonable accuracy and significant efficiency for both approaches, compared with those from the equivalent rainflow analysis. With much less computational effort and disk storage requirement, the approaches outlined in this paper can therefore be used for the fatigue assessment of deepwater risers in industry practice.     

基于NURBS的船体曲面自适应三角网格剖分

针对面向曲面的三维船体性能计算和真实感图形显示问题,应用NURBS曲线、曲面理论,提出一种新颖的船体NURBS曲面三角形网格自动生成算法,运用四角编码方法和改进的曲面片平坦性检验方法,保证在递归分割船体NURBS曲面时,能够快速有效地分割出四边形网格,在曲面片的高度方向和边界处同时满足给定的精度要求,在此基础上,应用割角剖分算法将一个四边形网格剖分成两个或多个三角形网格。应用结果表明,应用该算法生成的三角形平面片能够较好地逼近船体曲面,避免出现网格间的裂缝,与二叉树、四叉树方法相比,四角编码方法明显节省了时间和空间,提高了算法效率。     

海洋环境中复合桩基水平方向波浪荷载研究

为研究承台淹没深度及波浪入射角对桩-承台复合基础水平方向波浪荷载的影响,以平潭海峡公铁两用大桥某复合基础为研究对象,通过求解RANS方程（Reynolds-averaged Navier-Stokes equation）和k-ε湍流模型,借助FLOW-3D软件建立了波浪与复合基础相互作用的三维数值模型.研究了承台淹没深度及波浪入射角对复合基础水平方向波浪荷载的影响,并给出了不同波浪入射角下复合基础水平方向波浪荷载幅值随承台淹没深度的变化曲线.研究结果表明:承台位于波峰与波谷之间时,复合基础所受波浪荷载较大,承台淹没深度为1倍承台高度时波浪荷载达到峰值;当波浪入射角定义为承台长轴线与波浪传播方向之间夹角时,复合基础所受波浪荷载随着波浪入射角增大,当波浪入射角达到90°时,波浪荷载最大,其值约为波浪入射角为0°时的1.4倍左右.      

基于NURBS的三维船体曲面重构

针对重大件货仿真系统中船体曲面造型问题，给出了基于NURBS(非均匀有理B样条方法)技术重构船体曲面的主要步骤。对于给定的船舶型值点，用累加弦长参数化方法构造节点矢量，用非均匀有理B样条方法进行横剖线和水线的全局插值，并构造出规整的船体曲面的插值点网格，进而重构出以NURBS为统一数学表达式的船体曲面。仿真结果表明应用NURBS方法重构出的三维船体具有高度真实感。     

高速列车突入隧道引起的压缩波的理论研究

应用气动声学理论对高速列车突入隧道引起的复杂压力场进行了研究．根据实际情况，对气动声学的Ffowcs Williams-Hawkings方程(FW—H方程)进行简化，以节省计算资源和计算时间．应用伽利略变换对简化后的方程进行变量变换，使得曲面函数仅与空间相关；再直接对方程进行傅立叶变换，将其从时域转换到频域．最后，采用格林函数法求解FW—H方程，得到了高速列车突入隧道产生的压缩波的波形曲线，该曲线与既有模型试验结果一致．