共查询到20条相似文献，搜索用时 31 毫秒
1.
On the basis of the reproducing kernel particle method (RKPM), a new meshless method, which is called the complex variable reproducing kernel particle method (CVRKPM), for two-dimensional elastodynamics is presented in this paper. The advantages of the CVRKPM are that the correction function of a two-dimensional problem is formed with one-dimensional basis function when the shape function is obtained. The Galerkin weak form is employed to obtain the discretised system equations, and implicit time integration method, which is the Newmark method, is used for time history analysis. And the penalty method is employed to apply the essential boundary conditions. Then the corresponding formulae of the CVRKPM for two-dimensional elastodynamics are obtained. Three numerical examples of two-dimensional elastodynamics are presented, and the CVRKPM results are compared with the ones of the RKPM and analytical solutions. It is evident that the numerical results of the CVRKPM are in excellent agreement with the analytical solution, and that the CVRKPM has greater precision than the RKPM. 相似文献
2.
本文采用变域变分原理,建立了导热几何形状反演问题的变分原理,同时获得了该问题所需满足的边界条件和附加条件.该变分原理能将未知形状的几何变量及控制方程结合在一个变分泛函中,使得数学描述简洁、紧凑,且几何变量及控制方程的求解能耦合地进行.介绍了运用该变分原理并结合有限元方法进行数值计算的方法. 相似文献
3.
4.
Simulation of Fully Nonlinear 3-D Numerical Wave Tank 总被引：6，自引：0，他引：6
A fully nonlinear numerical wave tank (NWT) has been simulated by use of a three-dimensional higher order bouodary element method (HOBEM) in the time domain. Within the frame of potential flow and the adoption of simply Rankine source, the resulting boundary integral equation is repeatedly solved at each time step and the fully nonlinear free surface boundary conditions are integrated with time to update its position and boundary values. A smooth technique is also adopted in order to eliminate the possible saw-tooth numerical instabilities. The incident wave at the uptank is given as theoretical wave in this paper. The outgoing waves are absorbed inside a damping zone by spatially varying artificial damping on the free surface at the wave tank end. The numerical results show that the NWT developed by these approaches has a high accuracy and good numerical stability. 相似文献
5.
《Applied Ocean Research》2005,27(4-5):224-234
The modified scaled boundary finite-element method (SBFEM), keeping the advantages of the original SBFEM, eliminates the restriction of the scaling center location so that this approach can solve two-dimensional problems with parallel side-faces. In this paper, the modified SBFEM is applied to solutions of two types of problems—wave diffraction by a single and twin surface rectangular obstacles and wave radiation induced by an oscillating mono-hull and twin-hull structures in a finite depth of water. For wave diffraction problems, numerical results agree extremely well with the analytic solution for the single obstacle case and other numerical results of a different approach for the twin obstacle case. For wave radiation problems, the particular solutions to the scaled boundary finite-element equation are presented for cases of heave, sway and roll motions. The added mass and damping coefficients for heave, sway and roll motions of a two-dimensional rectangular container are computed and the numerical results are compared with those from independent analytical solution and numerical solution using the boundary element method (BEM). It is found that the SBFEM method achieves equivalent accuracy to the conventional BEM with only a few degrees of freedom. In the last example, wave radiation by a two-dimensional twin-hull structure is analyzed. Comparisons of the results with those obtained using conventional Green's function method (GFM) demonstrate that the method presented in this paper is free from the irregular frequency problems. 相似文献
6.
A three-dimensional numerical model for determination of the interaction between non-linear water waves and a structure is developed. The model is based on a boundary integral equation method for the spatial solution of a potential theory problem, combined with a time-stepping method based on the fully non-linear free surface conditions for temporal updating of moments on a structure in the fluid domain. Comparison with experimental results shows good agreement. The present model is considered to be one of the steps towards a three-dimensional numerical model in which the wave-structure interaction in a wave tank can be simulated. 相似文献
7.
Wave-induced transient response of seabeds is numerically analyzed through a radial point interpolation meshless method (radial PIM). The Biot’s consolidation theory is employed and incorporated with virtual boundary conditions to describe this wave-induced transient response of the seabed. Displacement and pore water pressure are spatially discretized by the radial PIM with the same shape function. Compactly supported basis functions are proposed to obtain a banded system equation. Because the radial PIM passes through all nodal points within an influence domain, essential boundary conditions as well as virtual boundary conditions can be easily implemented at local level. Fully implicit integration scheme is used in time domain to avoid spurious ripple effect. The proposed algorithm is assessed through the comparison of numerical results with closed-form solution or finite element solutions. 相似文献
8.
Based on the integral equation transformed from three dimensional Laplace equation and by the adoption of the division manner of sub- region boundary element method, the numerical computations of the velocity potential of each sub-region are given considering the continuity conditions of potential and normal derivatives at the interface of sub-regions. Therefore, computation of wave deformation in offshore flow field is realized. The present numerical model provides a good solution for the application of boundary element method to the calculation of wave deformation in large areas. 相似文献
9.
在曲线坐标系下,建立了缓变水深水域波浪传播的数值模拟模型.模型适宜于复杂变化的边界形状,克服了各种代数坐标变换的局限性.在建立模型时,将原始的椭圆型缓坡方程的近似型式——依赖时间变化的抛物型方程,作为控制方程,既克服了一般抛物近似方法的缺点,又便利了方程的求解;从开边界条件、不同反射特性的固壁边界条件相统一的表达式出发,对边界条件进行处理;用ADI法数值求解控制方程.对模型的验证表明,数值解与物模实验值吻合良好,模型对于具有复杂边界的工程实际有较强的适应性. 相似文献
10.
A three-dimensional finite difference transport model appropriate for the coastal environment is developed for the solution of the three-dimensional convection-diffusion equation. A higher order upwind scheme is used for the convective terms of the convection-diffusion equation, to minimise the numerical diffusion. The validity of the numerical model is verified through five test problems, whose exact solutions are known. 相似文献
11.
A. Hayir B. Seseogullari İ. Kilinc A. Erturk H.K. Cigizoglu M.S. Kabdasli O. Yagci K. Day 《Coastal Engineering》2008
In this paper the tsunamis resulting from a submarine mass failure such as slides and slumps triggered by earthquakes or other environmental effects, which is settled at the bottom of the north eastern Sea of Marmara are examined in one sample region. As the solution method, one hybrid method is developed. The main objective of this method is to combine an analytical solution presenting near-field tsunami amplitudes above the submarine mass failure with a numerical solution indicating the tsunami amplitudes in the coastal regions. For this purpose, one common linear boundary between analytical and numerical solution domains is defined. Movements of Submarine Mass Failures (SMF) are modeled using one simple kinematics source model and the amplitudes of the tsunamis at the region that are closer to the landslide are computed by using the analytical method. SMF is modeled approximately from the bottom geometry, and an average depth is used. Scenarios of SMF are established depending on the velocities and thicknesses of the failure, and near-field tsunami amplitudes are obtained in the open sea during the source time. After the source times, the solutions are found in the numerical region using TELEMAC-2D software system with the mentioned boundary above. In this boundary, the output of the analytical solutions is taken as the boundary conditions or the disturbances for the numerical method. With these disturbances, the numerical method is performed and the amplitudes are calculated in the coastal area. The generation, propagation and coastal amplifications of the tsunamis are illustrated at some certain points and regions both in the open sea and near the coast line. The results have been visualized and discussed. 相似文献
12.
The problem of static three-dimensional, nonlinear, large deformation response of a marine riser is formulated within small strain theory and solved numerically. This type of analysis is necessary, for the new generation of drilling and production risers. The mathematical model takes properly into account the effects of internal and external pressure and the complete nonlinear boundary conditions, without linearizing the follower forces. The extensibility or inextensibility condition is used as the constitutive relation in the tangential direction. Torsion and bending are coupled. The external load and the boundary conditions are deformation dependent. A solution method is developed based on an incremental finite element algorithm, which involves a prediction-correction scheme. In the correction phase deformation dependent quantities are updated. The extensibility or inextensibility condition is used to reduce the degrees of freedom of the system. The numerical results of the developed computer code compare very well with available semi-analytical and numerical solutions. Three numerical applications are used to demonstrate the importance of large deformation, nonlinear and three-dimensional analyses. 相似文献
13.
Three-dimensional liquid sloshing in a tank with baffles 总被引：1，自引：0，他引：1
Dongming Liu 《Ocean Engineering》2009,36(2):202-52
A numerical model has been developed to study three-dimensional (3D) liquid sloshing in a tank with baffles. The numerical model solves the spatially averaged Navier-Stokes equations, which are constructed on a non-inertial reference frame having six degree-of-freedom (DOF) of motions. The large-eddy-simulation (LES) approach is employed to model turbulence by using the Smagorinsky sub-grid scale (SGS) closure model. The two-step projection method is employed in the numerical solutions, aided by the Bi-CGSTAB technique to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate volume-of-fluid (VOF) method is used to track the distorted and broken free surface. The baffles in the tank are modeled by the concept of virtual boundary force (VBF) method. The numerical model is first validated against the available analytical solution and experimental data for two-dimensional (2D) liquid sloshing in a tank without baffles. The 2D liquid sloshing in tanks with baffles is then investigated. The numerical results are compared with other results from available literatures. Good agreement is obtained. Finally, the model is used to study 3D liquid sloshing in a tank with vertical baffles. The effect of the baffle is investigated and discussed. 相似文献
14.
15.
A numerical model for coastal water wave motion that includes an effective method for treatment of non-reflecting boundaries is presented. The second-order one-way wave equation to approximate the non-reflecting boundary condition is found to be excellent and it ensures a very low level of reflection for waves approaching the boundary at a fairly wide range of the incidence angle. If the Newman approximation is adopted, the resulting boundary condition has a unique property to allow the free propagation of wave components along the boundary. The study is also based on a newly derived mild-slope wave equation system that can be easily made compatible to the one-way wave equation. The equation system is theoretically more accurate than the previous equations in terms of the mild-slope assumption. The finite difference method defined on a staggered grid is employed to solve the basic equations and to implement the non-reflecting boundary condition. For verification, the numerical model is then applied to three coastal water wave problems including the classical problem of plane wave diffraction by a vertical circular cylinder, the problem of combined wave diffraction and refraction over a submerged hump in the open sea, and the wave deformation around a detached breakwater. In all cases, the numerical results are demonstrated to agree very well with the relevant analytical solutions or with experimental data. It is thus concluded that the numerical model proposed in this study is effective and advantageous. 相似文献
16.
Lin Mingchung Hsiao Sungshan Hsu Yungcheng
Professor Dept. of Naval Architecture Ocean Engineering National Taiwan University Taiwan. Doctor Course Student Dept. of Naval Architecture Ocean Engineering National Taiwan University Taiwan 《中国海洋工程》1994,(3)
-Wave refraction-diffraction due to a large ocean structure and topography in the presence of a 'current are studied numerically. The mathematical model is the mild-slope equation developed by Kirby (1984). This equation is solved using a finite and boundary element method. The physical domain is devid-ed into two regions: a slowly varying topography region and a constant water depth region. For waves propagating in the constant water depth region, without current interfering, the mild- slope equation is then reduced to the Helmholtz equation which is solved by boundary element method. In varying topography region, this equation will be solved by finite element method. Conservation of mass and energy flux of the fluid between these two regions is required for composition of these two numerical methods. The numerical scheme proposed here is capable of dealing with water wave problems of different water depths with the main characters of these two methods. 相似文献
17.
A. P. K. De Zoysa 《Ocean Engineering》1978,5(3):209-223
A numerical solution and computer program developed to analyse undersea flexible cable problems is presented. The cable is subject to hydrodynamic loads from a surrounding fluid medium. The object is to determine the cable profile and tensile forces developed. The cable equations and boundary conditions result in a classical two-point boundary value problem. An iteration ‘shooting method’ is used for solution. The program has been written to analyse a class of towing and anchoring problems. However the analysis is general and applicable to most steady-state cable problems. 相似文献
18.
HUANG Yuying 《中国海洋工程》2000,14(1):1-13
—As a further development of the authors'work(Huang and Qian,1993),in this paper a newnumerical method based on the time domain boundary element technique is proposed for solving fluid-sol-id coupling problems,in which a rigid body impacts normally on the calm surface of a half-space fluid.Afundamental solution to the half-space potential flow problem is first derived with the method of images.Then,an equivalent boundary integral equation in the Laplace transform domain is established by meansof Green's second identity.Through the inverse Laplace transform and discretization in both time andboundary of the fluid region,the numerical calculation for the problem under consideration has been car-ried out.Several examples demonstrate that the present method is more efficient than existing ones,fromwhich it is also seen that the shape of the impacting body has a considerable effect on the total impactforce. 相似文献
19.
The boundary integral equation method (BIEM) is developed as a tool for studying two-dimensional, nonlinear water wave problems, including the phenomena of wave generation, propagation and run-up. The wave motions are described by a potential flow theory. Nonlinear free-surface boundary conditions are incorporated in the numerical formulation. Examples are given for either a solitary wave or two successive solitary waves. Special treatment is developed to trace the run-up and run-down along a shoreline. The accuracy of the present scheme is verified by comparing numerical results with experimental data of maximum run-up. 相似文献
20.
The interaction of obliquely incident surface gravity waves with a vertical flexible permeable submerged membrane wave barrier is investigated in the context of three-dimensional linear water wave theory. From the general formulation of the submerged membrane barrier, the performance of bottom-standing, surface-piercing and fully extended membrane wave barriers are analyzed for various values of wave and structural parameters. The analytic solution of the physical problem is obtained using eigenfunction expansion method and a coupled boundary element-finite difference method has been used to get the numerical solution. In the boundary element method, since the boundary condition on the membrane barrier is not known a priori, the membrane response and velocity potentials are solved simultaneously using appropriate discretization with the help of finite difference scheme. The convergence of the analytic and numerical solution techniques is discussed. The study reveals that for suitable combination of wave and structural parameters, approximately (45–50)% incident wave energy can be dissipated irrespective of membrane barrier configurations. Further, in certain situations, nearly full wave reflection and zero transmission occur for all barrier configurations. The study will be useful in the design of flexible permeable membrane to act as an effective wave barrier for creation of tranquility zone in the marine environment. 相似文献