床面上直立圆柱的三维湍流数值模拟

从数值预报桥墩等结构物床面局部冲刷的角度发展绕直立圆柱的三维湍流的数值模拟技术. 基于Wilcox的k-ω两方程湍流模式，采用基于有限体积法的压力修正SIMPLE算法, 计算了绕床面直立圆柱的三维湍流流场，分析了光滑和粗糙床面两种情况下的流动情况. 通过系列的验证计算，表明该计算模型能够比较准确地反映不同外来流条件下绕直立圆柱的流场. 计算结果揭示了床面粗糙度对绕圆柱的湍流流动的影响.     

船-冰作用中海冰典型破坏模式的离散元数值模拟

在破冰船破冰过程中,冰排主要表现为挤压与弯曲两种破坏模式.本文基于黏结离散单元法对船-冰作用中的这两种典型破坏模式进行数值模拟.海冰离散元数值试样采用随机排布方式生成,采用单轴压缩试验与三点弯曲试验相结合的方式标定模型中的细观参数.将船-冰碰撞中的挤压和弯曲作用方式简化为直立或倾斜 平板与海冰的作用模式,构建挤压与弯曲破坏的海冰离散元数值试样,分析了挤压破坏模式中不同加载方位以及冰厚和加载速率对破坏模式的影响,以及弯曲破坏模式中不同冰排夹角以及冰厚和加载速率对破坏模式的影响.计算结果表明,离散元海冰数值模型可以很好地对冰排挤压与弯曲破坏现象进行模拟,可揭示冰排在不同破坏模式下的破坏机理.     

高速钢球在水中贯穿过程的数值计算

本文采用二维有限元数值计算方法对水中高速钢球的侵彻贯穿作用及贯穿通道的形成和发展进行了数值模拟计算,给出了钢球侵彻速度水中压力波和空腔发展,这些结果与实验观察较为吻合。本文的工作说明,数值方法可以作为研究创伤力学的重要手段,有着广阔的应用前景。     

基于平均法的旋转薄壁圆柱壳非线性行波振动研究

应用Donnell's简化壳理论,在考虑阻尼和几何非线性的情况下,基于平均法对旋转的薄壁悬臂圆柱壳在法向激振力作用下的非线性行波振动进行了研究.在分析过程中,首先,引入考虑阻尼及几何非线性的薄壁圆柱壳非线性波动方程,进行降阶处理后,得到模态坐标下的振动方程;其次,对模态方程进行平均化处理,确定转换矩阵,进行变量的幅值相角化,从而得到自治的标准化方程组;最后,由系统谐波共振周期解对应平均方程稳态解的原理,得到幅频特性方程.根据上述所得结果,进行了系统参数振动及稳定性研究,并进一步将结果与谐波平衡法及数值解作了比较.     

水中冲击波与弹性薄板耦合作用的研究

对水中冲击波与弹性薄板的流固耦合作用进行了研究,建立了实验装置。利用全息干涉法定量测量了水中动态流场,并测出受水中冲击波作用的弹性薄板在不同时刻的变形,得出相应的全息干涉条纹图。同时利用有限元法对水中冲击波与弹性薄板的流固耦合问题进行了数值计算。研究结果表明,用实验方法能够较好地模拟水中冲击波与弹性结构的耦合作用,用有限元法计算水中冲击波传播问题是可行的。     

旋转圆环动态特性的行波解

基于固体弹性介质中的行波理论,对旋转圆环的固有频率进行了计算.以扩张弯曲梁微元体的运动方程为基础,分析了运动条件下该波导中行波的色散关系.通过分离正行波数和负行波数,解决了行波法计算旋转圆环的固有频率问题.同时给出了波数分离的判断依据,并结合简谐波动的相位封闭原理和传递矩阵建立了旋转圆环的特征方程.数值模拟研究了旋转体...     

船舶在进出船厢运动中的附加质量的计算

采用三维有限元方法，计算船在进出船厢运动中的附加质量。文中对单圆柱体在无限域和半无限域中运动的附加质量及双圆柱体同心放置时内圆柱运动的附加质量分别进行了计算，计算值与理论值吻合很好；对船在进出船厢运动中的附加质量做了大量的计算，结果表明：随着船与船厢侧壁间距的减小，船的附加质量系数增大；随着船厢中水深的减小，船的附加质量系数增大，随着船离船厢封闭端的距离的减小，附加质量系数增大。     

船行波

船舶航行时,由船体对水的作用产生压力变化而形成的表面水波,称为船行波。因为在船的首部和尾部水流压力变化最大,所以在那里产生的波浪也最明显.由于产生波浪常需要能量,所以船行波越明显,船舶兴波阻力也往往越大。船舶兴波阻力的大小不仅和船速有关,而且还和船型有密切关系。随着航速的提高,兴波阻力占船舶总阻力的比例不断上升,高速舰船的兴波阻力可达到总阻力的50％以上。因此,如何应用消波的方法,减弱或消除船行波,从而减小兴波阻力以进一步提高航速,就成为一个很有现实意义的问题。     

刚性角隅内圆柱壳结构声辐射特性研究

本文对流场中处于刚性角隅内圆柱壳结构振动-声辐射问题进行研究．基于双反射方法，推导了位于两个垂直刚性壁面角隅区内圆柱壳结构的振动-声辐射方程．以此为基础，开展受线环向激励力作用的圆柱壳结构的振动-声辐射特性的数值计算．研究了刚性角隅内圆柱壳结构布置位置及计算频率对声辐射功率、声指向性的影响．计算结果可为分析含复杂声学边界的结构声振问题提供技术支持．     

含轴向裂纹柱壳裂纹尖端应力应变场及应力强度因子计算

本文从Reissner圆柱壳理论出发,应用摄动法获得了含轴向裂纹圆柱壳裂纹尖端应力应变场(包括Ⅰ、Ⅱ、Ⅲ型),并进一步应用Local-Global方法对不同尺寸块壳的应力强度因子进行了计算分析,同时对工程中常用的鼓胀系数进行了计算和分析讨论。计算结果表明,对于a/h较大的情况,经典公式是适用的,若a/h不太大时,经典理论将带来较大误差,本文给出了考虑剪切刚度影响的鼓胀系数的一些数值范围。     

The interaction of surface waves with fixed semi immersed cylinders having symmetric cross sections

The reflection of water waves by a semi immersed cylinder having a symmetric cross section is studied for both Dirichlet and Neumann boundary conditions on the cylinder. The method of conformal transformations as utilized by Ursell and by Tasai for the radiation problem is adapted to the present diffraction problem. The problem is solved by expansions of the reflected wave potential using nonorthogonal functions (wave free potentials). These functions are not complete, and an additional source and a dipole are required. Infinite systems of linear equations are obtained for the unknown expansion coefficients and the unknown strengths of the source and the dipole terms. Numerical results are obtained for the reflection coefficient, transmission coefficient, horizontal force on cylinder, vertical force on cylinder. In the long wave region analytical approximations are obtained for these functions when the cross section is circular. The reflection and transmission coefficients are very different for the two boundary conditions in the long wave region, the Dirichlet reflection coefficient being much larger than the corresponding Neumann coefficient. This behavior is similar to acoustic and electromagnetic diffraction problems in two dimensions. On leave of absence from Itek Corporation, Lexington (Mass.), U.S.A.     

面元法求解有限水深船舶兴波及水底压力变化

应用势流理论中的格林函数方法计算了船舶定常运动的水动力参数,将有限水深Kelvin移动兴波源格林函数分解成三部分:简单Rankine源集合、局部扰动项和波函数项。在亚临界和超临界航速时,采用不同的积分顺序来消除被积函数的奇异性。利用面元法在船体表面上分布Kelvin源,计算了有限水深下船体表面的源强、压力分布及表面兴波,比较了有限与无限水深结果的区别和联系,进一步求解了船舶航行时引起的水底压力变化,计算结果与实验测量结果吻合良好。     

Prediction of wave pattern and wave resistance of surface piercing bodies by a boundary element method

An iterative boundary element method, which was originally developed for both two‐ and three‐dimensional cavitating hydrofoils moving steadily under a free surface, is modified and extended to predict the wave pattern and wave resistance of surface piercing bodies, such as ship hulls and vertical struts. The iterative nonlinear method, which is based on the Green theorem, allows the separation of the surface piercing body problem and the free‐surface problem. The free‐surface problem is also separated into two parts; namely, left and right (with respect to x axis) free‐surface problems. Those all (three) problems are solved separately, with the effects of one on the other being accounted for in an iterative manner. The wetted surface of the body (ship hull or strut, including cavity surface if exists) and the left and right parts with respect to x axis of free surface are modelled with constant strength dipole and constant strength source panels. In order to prevent upstream waves, the source strengths from some distance in front of the body to the end of the truncated upstream boundary are enforced to be zero. No radiation condition is enforced for downstream and transverse boundaries. A transverse wave cut technique is used for the calculation of wave resistance. The method is first applied to a point source and a three‐dimensional submerged cavitating hydrofoil to validate the method and a Wigley hull and a vertical strut to compare the results with those of experiments. Copyright © 2007 John Wiley & Sons, Ltd.     

Diffraction of regular waves by an arbitrarily oriented system of vertical circular cylinders

The linear theory of small-amplitude waves is used to construct a solution to the problem of the diffraction of surface gravity waves by a system of arbitrarily oriented vertical circular cylinders. Analytic expressions are obtained for the wave forces and overturning moments acting on each cylinder of the system. A system of two rows of cylinders with three cylinders in each row is considered as an example. It is shown that for certain relationships between the diameter of the cylinders, the distance between them, the angle of approach of the wave, and its wavelength, the maximal values of the hydrodynamic forces acting on the cylinders of the system with allowance for the interaction of the diffracted fields may be appreciably greater than when no allowance is made for it.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 177–181, January–February, 1984.     

时域Rankine源法求解有限水深船舶在规则波中的水底压力变化

应用势流理论中的Rankine源面元法和时域步进法,求解了有限水深船舶在规则波中运动的水底压力变化。将速度势分解成基本势、局部势和记忆势,以叠模解作为基本势对自由表面条件和物面条件进行了线性化,通过在水底布置面元来满足水底条件。利用研制的水底压力-水面波浪测量系统,测量了不同入射波船模表面波形与水底压力的时历曲线,理论计算与实验结果符合较好,验证了自编程序的正确性。通过对比二者的等高线图发现,水底压力与表面波形的峰谷有较好的一致性,并且压力较波形更为平滑。     

直立码头前船波浪力耦合计算模型

建立了外域用Boussinesq方程、内域用刚体运动方程的直立码头前二维船剖面波浪力的时 域计算耦合模型，内域与外域在交界面的匹配条件是流量连续和压力相等. 进行了相关模型 实验，并把计算结果与实验结果进行了对比. 推导了船体与水底和直立码头之间间隙内流体 运动的自振频率，研究了间隙内流体运动的共振现象.     

Surface characters of internal waves generated by Rankine ovoid

A linear theory on the internal waves generated in the stratified fluid with a pycnocline is presented in this paper. The internal wave fields such as the velocity fields in the stratified fluid and velocity gradient fields at the free surface are also investigated by means of the theoretical and numerical method. From the numerical results, it is shown that the internal wave generated by horizontally moving Rankine ovoid is a sort of trapped wave which propagates in a wave guide, and its waveform is a kind of Mach front-type internal wave in the pycnocline. Influence of the internal wave on the flow fields at the free surface is represented by the velocity gradient fields resulted from the internal waves generated by motion of the Rankine ovoid. At the same time, it is also shown that under the hypothesis of inviscid fluid, the synchronism between the surface velocity gradient fields at the free surface and the internal wave fields in the fluid is retained. This theory opens a possibility to study further the modulated spectrum of the Bragg waves at the free surface.The project supported by the National Natural Science Foundation of China (40576010). The English text was polished by Keren Wang.     

A second order stationary solution of the three dimensional wave diffraction potential

A theory on the second order wave diffraction by a three dimensional body fixed in a regular sea has been developed in the present paper. By regarding the sinusoidal disturb potential as a stationary solution of an initial value problem, and using Laplace transformation method and Tauberian theorem, the boundary value problems of stationary solution of the first and second order diffraction potential have been derived in this paper. Furthermore, the explicit solution of the second order stationary diffraction potential has been obtained with the method of the double Fourier transformation. It is found that the asymptotic behaviour of the second order stationary solution at far field is dependent on two wave systems, the first is “free wave”, travelling independently of the first order wave system, the other is “phase locked waves”, which accompany the first order waves. At the same time, the radiation conditions of the second order diffraction problems are derived. We also find that one can still pursue a steady state formulation with the inclusion of Rayleigh damping. Finally, as an example, the second order wave forces upon a fixed vertical circular cylinder have been calculated, and the numerical results agree well with the experimental data.     

变深度浅水域中非定常船波

以Green—Naghdi(G—N)方程为基础，采用波动方程／有限元法计算船舶经过变深度浅水域时非定常波浪特性．把运动船舶对水面的扰动作为移动压强直接加在Green-Naghdi方程里，以描述运动船体和水面的相互作用．以Series60 CB＝0．6船为算例，给出自由面坡高，波浪阻力在船舶经过一个水下凸包时变化规律，并与浅水方程的结果进行了比较．计算结果表明，当船舶经过凸包时，波浪阻力先增加，后减少，并逐渐趋于正常．同时发现，当船速小于临界速度时(Fr＝√gh＜1．0)，G—N方程给出的船后尾波比浅水方程的结果明显，波浪阻力也比浅水方程的结果有所提高，频率散射必须考虑．当船速大于临界速度时(Fr＝√gh＞1．0)，G—N方程的计算结果与浅水方程差别不大，频率散射的影响可以忽略．