近岸沿岸流及污染物运动的数值模拟

基于双曲型缓坡方程和近岸浅水方程对近岸波浪斜向入射破碎所生成的沿岸流及污染物在沿岸波流作用下的运动进行了数值模拟,并对数值模拟结果进行了验证分析。数值模拟结果表明,在相近工况参数下,随着入射波高的增大,沿岸流流速和平均水面升高值均明显增大;随着岸坡坡度的增加,沿岸流流速和平均水面升高值明显增大;随着入射波浪周期的增大,平均水面升高值明显增大。在沿岸缓坡区域,由斜向入射波浪破碎所产生的沿岸流对污染物的运动起着不可忽略的影响。     

一种斜井VSP的正演模型

在垂直入射斜井ＶＳＰ的条件下，假定地下为倾斜层状平面结构，在各层中地震射线为直线，界面的反射和透射服从Ｓｎｅｌｌ定律时，本文提出一种射线追踪法正演模型。射线追踪是采用迭代法实现的。用理论模型试算可以得到较为满意的结果。     

优化系数的四阶吸收边界条件

在波动方程正演模拟或偏移计算中往往要在边界上加上吸收边界条件来消除人工边界的虚假反射.传统的吸收边界条件方程阶数较低,对大角度入射波的吸收效果不理想.本文基于Higdon提出的渐进吸收边界条件方法,推导出四阶吸收边界条件方程,并提出了一种新的吸收边界条件角度参数的优化选择方法.实验证明应用此优化系数的四阶吸收边界条件对全角度入射的入射波都有较好的吸收效果.     

岩石圈挠曲-弹性薄板小挠度弯曲的新方程

岩石圈挠曲研究采用的是弹性薄板小挠度弯曲方程(即克希霍夫方程),克希霍夫方程基于薄板的前提,忽略并假设薄板内垂向应力为0。本文在无需垂向应力为0的这一与地质事实不相符的假设的情况下,由弹性体几何方程、物理方程和静力平衡方程推导出岩石圈挠曲-弹性薄板小挠度弯曲的新中面方程,具有同等的数学简洁性。取泊松比为0.25时,有DFF/D=1.125,即新挠曲方程中的挠曲刚度DFF要比经典的克希霍夫挠曲方程的D值大12.5%。本文推导的新方程不仅可以在岩石圈动力学,也可以在弹性力学中获得应用。     

波浪谱形对不规则波数值模拟的影响

通过数值模拟分析了波浪谱形对不规则波浪数值模拟结果的影响.采用不同参数的JONSWAP谱模拟入射波要素,基于抛物型缓坡方程模拟不规则波浪的传播,分析了波浪谱形状对波浪数值模拟结果的影响.结果表明,采用抛物型缓坡方程模拟不规则波浪时,入射波浪谱形对模拟结果影响不明显;但由于模型中非线性项的影响,采用不规则波模拟的波高分布和采用规则波模拟的结果略有差别.     

海洋声传播二维PE算法对海底边界条件的改进处理

运用海洋声传播二维 PE算法对不规则海底边界条件的处理 ,作者提出 1种改进方法 ,即利用抛物方程和界面边界条件 ,对界面附近的抛物方程的差分形式进行修正 ,以考虑复杂海底边界的影响。通过实际的数值计算表明 ,运用该方法计算声场 ,在提高计算精度方面起到有益的作用。     

分层海洋的内孤立波

本文从流体力学基本方程组出发,在非地转条件下导得了分层海洋的内孤立波方程—Kbv和mKdv方程,证实了在非地转条件下,一类海洋非线性波动是可以严格满足内孤立波方程的。在地转条件下采用f平面近似导出了KdV方程的演化形式一有源KdV方程,地转的影响含于源项中。由初步的分析得出,f对KdV方程的影响是微弱的。由已得的KdV和mKdV方程的解可知,内孤立波与线性波有着本质差别。     

大亚湾长周期波在缓变流作用下传播计算模型的研究

以存存缓变背景流的二维缓坡方程为控制方程建立一个计算模型来研究大亚湾长周期波的传播规律.按入射波源形式不同、周期不l司和有无背景潮流设计了几个数学试验,分别通过这儿个试验研究长周期波传播规律.从模拟结果发现,在大哑湾水域,仪在开边界附近不同源形式的入射波对波浪传播的影响有所变化;总的来说,从入射点到湾内波高削弱较快,但是对于不同周期的入射波,模拟所得波高分布变化较大;对于长周期的入射波,背景流也对波浪传播有一定影响.     

双突堤波浪绕射的数学模型

双突堤是常见的一种防波堤布置型式。因此研究双突堤的波浪绕射数学模型和数值计算方法是有理论意义和实用价值的。本文引用数学物理方程中的单源点法,提出了两种双突堤波浪绕射的数学模型。本模型采用线性简谐波理论,假设水深为定值。     

一种推导Timoshenko梁频率方程的新方法

在Timoshenko梁横向弯曲自由振动方程的基础上,首次推导建立两端受弹性约束梁的频率方程,通过此方程,研究两端受弹性约束梁的性质,并通过简支梁模型,初步估计剪切变形与转动惯量对Timoshenko梁频率的影响,估算Timoshenko梁与Euler梁频率之间的关系,给出估算公式,并用数值算例进行验证,结果表明数值法解Timoshenko频率方程的可靠性.     

Brillouin scattering study on elastic and piezoelectric properties of laser irradiated ZnO single crystals

Brillouin light scattering technique can be successfully used to determine the whole set of elastic and piezoelectric constants of a ZnO single crystal irradiated by different laser energy densities, into a micron range (radiation layer thickness). It is found that the scattering intensity, the linewidth and the Brillouin scattering shift of acoustic phonons are all strongly dependent on laser energy density. Based on the sound propagation equations and these results, the directional dependences of the compressional and shear moduli of the irradiated ZnO sample in the (001) plane are investigated. It is found that under an appropriate laser condition, 248 nm KrF excimer laser irradiation can significantly improve the surface quality and increase the elastic properties of ZnO single crystal. This procedure has potential applications in the fabrication of ZnO-based surface acoustic wave and optic-electronic devices.     

基于有限元模型和模态参数的损伤识别

本文从特征值问题出发 ,结合结构系统的有限单元模型并引入单元损伤因子 ,导出了求解单元损伤因子的方程 ,从而可以确定结构的损伤位置和损伤程度。利用简支梁进行了数值模拟 ,建立了简支梁的有限单元模型 ,用不同单元弹性模量的降低来模拟各种不同的损伤工况 ,并对每种工况进行了损伤位置和损伤程度识别 ,数值模拟结果同实际情况吻合很好。     

Experimental investigating on the reflected waves from the caisson-type vertical porous seawalls

The hydrodynamic efficiencies of caisson-type vertical porous seawalls used for protecting coastal areas were calculated in this study. Physical models were developed to compare the wave reflection from vertical plane, semi-porous, and porous seawalls caused by both regular and random waves. Tests were carried out for a wide range of wave heights, wave periods, and different water depths (d=0.165, 0.270 and 0.375 m). The performance regarding the reflected waves from porous and semi-porous seawalls showed improvement when compared with those from the plane seawall. The reflection coefficients of the porous and semi-porous seawalls were calculated as 0.6 and 0.75, respectively, while the coefficient for the fully reflecting plane vertical wall was significantly higher (0.9). It was also observed that the reflection coefficient decreases with increase in wave steepness and relative water depth. In addition, the reduction in the reflection coefficient of porous and semi-porous seawalls, as compared to that of a plane seawall, was observed for both regular and random waves. New equations were also proposed to calculate the reflection coefficient of different types of seawalls with the aid of laboratory experiments. By verifying the developed equations using some other experimental data, it was validated that the equations could be used for practical situations. The results of the present study can be applied to optimize the design of vertical seawalls and for coastal protecting schemes.     

Axisymmetric viscoelastic response of flexible pipes in time domain

The challenges for determining the mechanical behavior of flexible pipes mainly arise from highly non-linear geometrical and material properties and complex contact interaction conditions between and within layers components. This paper develops an innovative model to investigate the linear viscoelastic behavior of flexible pipes under axisymmetric loads in time domain. The model is derived from an equivalent linear elastic axisymmetric model by invoking the elastic-viscoelastic correspondence principle. Analytical formulations that describe the behavior of the metallic helical layers based on a combination of differential geometry concepts and Clebsch–Kirchhoff equilibrium equations for initially curved slender elastic rods are presented. The elastic response of the homogenous polymeric cylindrical layers is also presented. The assemblage of both types of governing algebraic equations that approximate analytical solutions for force and moment distributions, deformations in each layer, as well as contact pressure between near layers, taking time-dependent characteristics of polymeric layers into account are provided and it is clear that the relationship between axial force and elongation is non-linear and encompasses a hysteretic response. Besides, the creep behavior in axial direction can also be found. Some insights into the differences in the behavior for several loading conditions are discussed by considering variable frequencies.     

Dynamic factor analysis considering elastic boom effects in heavy lifting operations

The dynamic factor is the ratio of the maximum dynamic load to the static load acting on the wire ropes between the boom of a floating crane and a cargo. In this paper, the dynamic factor is analyzed based on dynamic simulations of a floating crane and a cargo, considering an elastic boom. For the simulation, we designed a multibody system that consists of a floating crane barge, an elastic boom, and a cargo connected to the boom through wire ropes. The dynamic equations of motion of the system are based on flexible multibody system dynamics. Six-degree-of-freedom motions are considered for the floating crane and for the cargo, and three-dimensional deformations for the elastic boom. The hydrostatic force, the hydrodynamic force, the gravitational force, and the wire rope forces are considered as external forces. The dynamic factor is obtained by numerically solving the equation. The effects of the elastic boom on heavy cargo lifting are discussed by comparing the simulation results of an elastic boom and a rigid boom.     

Interaction of ocean waves with floating and submerged horizontal flexible structures in three-dimensions

A three-dimensional general mathematical hydroelastic model dealing with the problem of wave interaction with a floating and a submerged flexible structure is developed based on small amplitude wave theory and linear structural response. The horizontal floating and submerged flexible structures are modelled with a thin plate theory. The linearized long wave equations based on shallow water approximations are derived and results are compared. Three-dimensional Green’s functions are derived using fundamental source potentials in water of finite and infinite depths. The expansion formulae associated with orthogonal mode-coupling relations are derived based on the application of Fourier transform in finite and infinite depths in case of finite width in three-dimensions. The usefulness of the expansion formula is demonstrated by analysing a physical problem of surface gravity wave interaction with a moored finite floating elastic plate in the presence of a finite submerged flexible membrane in three-dimensions. The numerical accuracy of the method is demonstrated by computing the complex values of reflected wave amplitudes for different modes of oscillation and mooring stiffness. Further, the effect of compressive force and modes of oscillations on a free oscillation hydroelastic waves in a closed channel of finite width and length for floating and submerged elastic plate system is analysed.     

Vibration analysis of rectangular Mindlin plates on elastic foundations and vertically in contact with stationary fluid by the Ritz method

In this study Free vibration analysis of vertical rectangular Mindlin plates resting on Pasternak elastic foundation and fully or partially in contact with fluid on their one side is investigated for different combinations of boundary conditions. The plate is assumed to be one of vertical rectangular walls of a container in contact with fluid. In order to analyze the interaction of the Mindlin plate with the elastic foundation and fluid system, three displacement components of the plate are expressed in the Ritz method by adopting a set of static Timoshenko beam functions satisfying geometric boundary conditions in a Cartesian co-ordinate system. The method of separation of variables and the method of Fourier series expansion is used to model fluid and to obtain the exact expression of the motion of fluid in the form of integral equations. The fluid domain is finite in depth and width but infinite in the length direction. To demonstrate the accuracy of the present solution, convergence study is first carried out and then a few comparison studies are carried out with the available data in the literature. Finally, natural frequencies of rectangular plates are presented in tabular and graphical forms for different fluid levels, foundation parameters, aspect ratios, thickness to width ratios and boundary conditions.     

Gravity-wave-induced pressure fluctuations in the deep ocean

We establish a mathematically consistent theory of the pseudo-sound pressure fluctuation in the deep ocean induced by nonlinearly interacting random plane waves on the surface. In the process, a new set of the second-order perturbation equations is derived and power-correlation coefficients between random plane waves are introduced. A phenomenological model is adopted for wind pressure which excites the surface waves consisting of wind-driven sea and swell. By solving the first-order- and the second-order-perturbation equations with this wind pressure as the excitation, we obtain an expression for the pressure fluctuation and its power spectral density in the gravity-wave regime. It is concluded that only the swell part of the surface waves generates the pressure fluctuation and the spectral density is modified by the power correlation coefficient.     

稳恒水波的Fourier近似解研究

A computational method for steady water waves is presented on the basis of potential theory in the physical plane with spatial variables as independent quantities. The finite Fourier series are applied to approximating the free surface and potential function. A set of nonlinear algebraic equations for the Fourier coefficients are derived from the free surface kinetic and dynamic boundary conditions. These algebraic equations are numerically solved through Newton＇s iterative method, and the iterative stability is further improved by a relaxation technology. The integral properties of steady water waves are numerically analyzed, showing that （1） the set-up and the set-down are both non-monotonic quantities with the wave steepness, and （2） the Fourier spectrum of the free surface is broader than that of the potential function. The latter further leads us to explore a modification for the present method by approximating the free surface and potential function through different Fourier series, with the truncation of the former higher than that of the latter. Numerical tests show that this modification is effective, and can notably reduce the errors of the free surface boundary conditions.