共查询到20条相似文献，搜索用时 234 毫秒
1.
The two-dimensional nonlinear time domain free surface flow problem is analysed using potential flow theory. The problem is solved by a time marching method. At each time step two numerical approaches are used. One is based on the boundary element method in the complex plane. The complex potential is assumed to vary linearly within each element and the solution is obtained by imposing the boundary conditions at the nodes of the elements. The other approach is based on the finite element formulation. Triangular elements and linear shape functions are used. The solution is obtained by the Galerkin method. Numerical results are obtained for the wave elevation generated by a vertical wave maker. Results are also provided for a circular cylinder oscillating below the free surface. For these cases the finite element method is found to provide substantially more efficient computations than the boundary element method using equivalent discretizations. 相似文献
2.
G. X. WuR. Eatock Taylor 《Applied Ocean Research》1994,16(6):363-372
The two-dimensional nonlinear time domain free surface flow problem is analyzed by the finite element method. Two approaches are used. One is based on the velocity potential which is approximated by means of shape functions. The solution is obtained through use of a variational statement, and the velocity is obtained subsequently by the Galerkin method. The other approach is to write both potential and velocity in terms of the shape functions at the same time. Their solutions are derived from the same equation by using another variational statement. Numerical results are given for the vertical wave maker problem and for a transient wave in a rectangular container. They are compared with analytical solutions, and very good agreement is found. 相似文献
3.
Mass loss caused by glacier calving is one of the direct contributors to global sea level rise. Reliable calving laws are required for accurate modelling of ice sheet mass balance. Both continuous and discontinuous methods have been used for glacial calving simulations. In this study, the discrete element method(DEM) based on dilated polyhedral elements is introduced to simulate the calving process of a tidewater glacier. Dilated polyhedrons can be obtained from the Minkowski sum of a sphere and a core polyhedron. These elements can be utilized to generate a continuum ice material, where the interaction force between adjacent elements is modeled by constructing bonds at the joints of the common faces. A hybrid fracture model considering fracture energy is introduced. The viscous creep behavior of glaciers on long-term scales is not considered. By applying buoyancy and gravity to the modelled glacier, DEM results show that the calving process is caused by cracks which are initialized at the top of the glacier and spread to the bottom. The results demonstrate the feasibility of using the dilated polyhedral DEM method in glacier simulations, additionally allowing the fragment size of the breaking fragments to be counted. The relationship between crack propagation and internal stress in the glacier is analyzed during calving process. Through the analysis of the Mises stress and the normal stress between the elements, it is found that geometric changes caused by the glacier calving lead to the redistribution of the stress. The tensile stress between the elements is the main influencing factor of glacier ice failure. In addition, the element shape,glacier base friction and buoyancy are studied, the results show that the glacier model based on the dilated polyhedral DEM is sensitive to the above conditions. 相似文献
4.
The hydrodynamic problem of 3D planing surface is studied by a finite element approach. The planing surface is represented by a number of pressure patches whose strengths are constant at each element. The unknown pressure strength is obtained by using the free surface elevation condition under the planing surface and Kutta condition at the transom stern. Previous studies indicate that, when the constant pressure distribution method is used, the number of buttocks should be less than five or six, otherwise the calculated pressure distribution will start to oscillate and even become divergent. In the present study, after careful examination of the influence coefficients, it is found that the accuracy of the influence coefficients matrix is very important to the convergence of the solution, especially when the number of elements is relatively high. The oscillation of the pressure distribution can be avoided by constant element method if the influence coefficients are sufficiently accurate. The predicted results of the present paper with more number of buttocks are in good agreement with other researchers'. It is concluded that the irregularity of the pressure distribution found in previous studies is most likely caused by the low accuracy in their calculation of the influence coefficients, not by the method itself. 相似文献
5.
比例边界有限元法(SBFEM)是一种半解析的数值方法,比完全数值方法具有更高的精度,该方法结合了有限元和边界元的优点,采用相对少的剖分单元就可以得到较高精度的模拟结果。通过改变有限子域内部比例中心的位置,使这种方法可以应用到多种形式浮体在波浪作用下的水动力特性的计算中。同时还给出了各种形式浮体的波浪力及反射、透射系数的数值结果,并与边界元方法(BEM)计算结果和特征函数展开方法得到解析解进行了比较,均吻合良好。研究表明比例边界有限元不仅可以计算矩形的浮体结构,而且对于多种结构形式的浮体都可以计算,这为多种结构形式浮体的水动力分析提供了一个可行的方法。 相似文献
6.
An obvious motivation of this paper is to examine the effectiveness of the lateral vibration control of a jacket type offshore platform with an AMD control device, in conjunction with H2 control algorithm, which is an optimal frequency domain control method based on minimization of H2 norm of the system transfer function. In this study, the offshore platform is modeled numerically by use of the finite element method, instead of a lumped mass model. This structural model is later simplified to be single-degree-of-freedom (SDOF) system by extracting the first vibration mode of the structure. The corresponding "generalized" wave force is determined based on an analytical approximation of the first mode shape function, the physical wave loading being calculated from the linearized Morison equation. This approach facilitates the filter design for the generalized force. Furthermore, the present paper also intends to make numerical comparison between H2 active control and the corresponding passive control using a T 相似文献
7.
Flexible marine risers are commonly used in deepwater floating systems.Bend stiffeners are designed to protect flexible risers against excessive bending at the connection with the hull.The structure is usually analyzed as a cantilever beam subjected to an inclined point load.As deflections are large and the bend stiffener material exhibits nonlinear stress-strain characteristics,geometric and material nonlinearities are important considerations.A new approach has been developed to solve this nonlinear problem.Its main advantage is its simplicity;in fact the present method can be easily implemented on a spreadsheet.Finite element analysis using ABAQUS is performed to validate the method.Solid elements are used for the bend stiffener and flexible pipe.To simulate the near inextensibility of flexible risers,a simple and original idea of using truss elements is proposed.Through a set of validation studies,the present method is found to be in a good agreement with the finite element analysis.Further,parametric studies are performed by using both methods to identify the key parameters and phenomena that are most critical in design.The most important finding is that the common practice of neglecting the internal steel sleeve in the bend stiffener analysis is non-conservative and therefore needs to be reassessed. 相似文献
8.
护岸工程在建造和运营中 ,经常会遇到各种风险。为了保护胜利油田陆上的油气开发设施 ,应用风险分析技术对黄河口附近护岸 (海堤 )工程的风险进行了探讨 ,并提出了该海区护岸工程防范风险的对策 相似文献
9.
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. 相似文献
10.
应用基于势流理论的时域高阶边界元方法,建立一个完全非线性的三维数值波浪水槽,通过实时模拟推板造波运动的方式产生波浪。通过混合欧拉-拉格朗日方法和四阶Runge-Kutta方法更新自由水面和造波板的瞬时位置。利用所建模型分别模拟了有限水深波和浅水波,与试验结果、相关文献结果和浅水理论结果吻合较好,且波浪能够稳定传播。系统地讨论造波板的运动圆频率、振幅和水深等对波浪传播和波浪特性的影响,并对波浪的非线性特性进行分析,研究发现造波板运动频率、运动振幅以及水深均将对波浪形态和波浪非线性产生显著影响。结果为真实水槽造波机的运动控制以及波浪生成试验提供了依据,便于实验室设置更合理的参数来准确模拟不同条件下的波浪。 相似文献
11.
An effective boundary element method (BEM) is presented for the interaction between oblique waves and long prismatic structures in water of finite depth. The Green's function used here is the basic Green's function that does not satisfy any boundary condition. Therefore, the discretized elements for the computation must be placed on all the boundaries. To improve the computational efficiency and accuracy, a modified method for treatment of the open boundary conditions and a direct analytical approach for the singularity integrals in the boundary integral equation are adopted. The present BEM method is applied to the calculation of hydrodynamic coefficients and wave exciting forces for long horizontal rectangular and circular structures. The performance of the present method is demonstrated by comparisons of results with those generated by other analytical and numerical methods. 相似文献
12.
Numerical solutions for the hydroelastic problems of bodies are studied directly in the time domain using Neumann–Kelvin formulation. In the hydrodynamic part of problem, the exact initial boundary value problem is linearized using the free stream as a basis flow, replaced by the boundary integral equation applying Green theorem over the transient free surface Green function. The resultant boundary integral equation is discretized using quadrilateral elements over which the value of the potential is assumed to be constant and solved using the trapezoidal rule to integrate the memory or convolution part in time. In the structure part of the problem, the finite element method is used to solve the hydroelastic problem. The Mindlin plate as a bending element, which includes transverse shear effect and rotary inertia effect are used. The present numerical results show acceptable agreement with experimental, analytical, and other published numerical results. 相似文献
13.
《Coastal Engineering》1999,37(2):97-122
In this paper, a numerical model based on the improved Boussinesq equations derived by Beji and Nadaoka [Beji, S., Nadaoka, K., 1996. A formal derivation and numerical modeling of the improved Boussinesq equations for varying depth. Ocean Eng. 23 (8), 691–704] is presented. The finite element method is used to discretize the spatial derivatives. Quadrilateral elements with linear interpolating functions are employed for the two horizontal velocity components and the water surface elevation. The time integration is performed using the Adams–Bashforth–Moulton predictor–corrector method. Five test cases for which either theoretical solutions or laboratory results are available are employed to test the proposed scheme. The model is capable of giving satisfactory predictions in all cases. 相似文献
14.
The finite element method(FEM) is employed to analyze the resonant oscillations of the liquid confined within multiple or an array of floating bodies with fully nonlinear boundary conditions on the free surface and the body surface in two dimensions.The velocity potentials at each time step are obtained through the FEM with 8-node quadratic shape functions.The finite element linear system is solved by the conjugate gradient(CG) method with a symmetric successive overelaxlation(SSOR) preconditioner.The waves at the open boundary are absorbed by the combination of the damping zone method and the Sommerfeld-Orlanski equation.Numerical examples are given by an array of floating wedgeshaped cylinders and rectangular cylinders.Results are provided for heave motions including wave elevations,profiles and hydrodynamic forces.Comparisons are made in several cases with the results obtained from the second order solution in the time domain.It is found that the wave amplitude in the middle region of the array is larger than those in other places,and the hydrodynamic force on a cylinder increases with the cylinder closing to the middle of the array. 相似文献
15.
In the present study, we simulated the reel-lay installation process of deepwater steel catenary risers（SCRs） using the finite element method and proposed multiaxial fatigue analysis for reeled SCRs. The reel-lay method is one of the most efficient and economical pipeline installation methods. However, material properties of reeled risers may change, especially in the weld zone, which can affect the fatigue performance. Applying finite element analysis（FEA）, we simulated an installation load history through the reel, aligner, and straightener and analyzed the property variations. The impact of weld defects during the installation process, lack of penetration and lack of fusion, was also discussed. Based on the FEA results, we used the Brown-Miller criterion combined with the critical plane approach to predict the fatigue life of reeled and non-reeled models. The results indicated that a weld defect has a significant influence on the material properties of a riser, and the reel-lay method can significantly reduce the fatigue life of SCRs. The analysis conclusion can help designers understand the mechanical performance of welds during reel-lay installation. 相似文献
16.
冰脊对极地船舶及海洋工程结构的冰载荷设计及冰区安全运行具有重要的影响,其几何形态是重要的考虑因素。为研究冰脊的几何形态,本文基于Voronoi切割算法构造扩展多面体海冰单元,并使其相互冻结形成平整冰;采用离散单元法模拟海冰在相对运动时挤压形成冰脊的动力过程,统计分析脊帆高度、龙骨宽度、龙骨深度、龙骨水平倾角等主要几何参数,并确定这些几何参数间的对应关系。将离散元模拟结果同冰脊现场观测资料对比分析以验证扩展多面体离散元方法模拟冰脊形成过程的可行性。在此基础上讨论冰厚、冰速、海冰强度等因素对冰脊形态的影响。本文采用离散元方法对冰脊形成过程及几何特性的数值分析可为深入开展冰脊的形成机理及其对船舶、海洋工程结构物的冰载荷分析提供参考依据。 相似文献
17.
近年来,有限元方法[1]已被许多学者用来解浅水波方程,如潮汐潮流、近海环流、污染扩散等[2,3],也有用来讨论风暴潮问题的.本文给出用于风暴潮分析的一种有限元方法. 我国是一个多风暴潮的国家,因此,沿海地区防御风暴潮是很重要的.至今已提出单站预报、动力统计、差分数值等风暴潮预报方法.采用有限元方法能灵活的近似海岸线形状和海底形状,所得出的有限元公式能给出一个通用的计算机程序,适用于各个海区的预报工作.为计算简便,本文采用九点四边形Lagrange等参数元[4]和分离时间的时间积分格式,从而使得程序简化,并大大节约内存和计算机时间. 相似文献
18.
应用比例边界有限元法求解狭缝对双箱水动力的影响 总被引：1，自引：1，他引：0
比例边界有限元法(SBFEM)是一种半解析数值分析的新方法,既融合了有限元法和边界元法的优点,又有其特有的优点。用该方法可求解有限水深下狭缝对双箱水动力作用的影响,为波浪与多浮体超大型结构的相互作用探索一些规律。整个计算域划分成2个无限子域和4个有限子域,并利用加权余量法在各个子域上推导了SBFEM的积分方程;计算了4个数值算例并与边界元等其它数值方法进行了比较,验证了该方法是一种用很少单元便能得到精确结果的高效方法。应用SBFEM对不同箱体宽度、不同狭缝宽度、不同吃水深度条件的双箱作了计算,得出了狭缝对双箱水动力干涉影响的一些规律,对超大型浮体水动力分析和结构设计具有一定的参考价值。 相似文献
19.
The finite element method (FEM) is employed to analyze the resonant oscillations of the liquid confined within multiple or an array of floating bodies with fully nonlinear boundary conditions on the free surface and the body surface in two dimensions. The velocity potentials at each time step are obtained through the FEM with 8-node quadratic shape functions. The finite element linear system is solved by the conjugate gradient (CG) method with a symmetric successive overelaxlation (SSOR) preconditioner. The waves at the open boundary are absorbed by the combination of the damping zone method and the Sommerfeld-Orlanski equation. Numerical examples are given by an array of floating wedge- shaped cylinders and rectangular cylinders. Results are provided for heave motions including wave elevations, profiles and hydrodynamic forces. Comparisons are made in several cases with the results obtained from the second order solution in the time domain. It is found that the wave amplitude in the middle region of the array is larger than those in other places, and the hydrodynamic force on a cylinder increases with the cylinder closing to the middle of the array. 相似文献
20.
A finite element method is developed for the stiffened multi-layered airfoil/hydrofoil structure for the large deformation and finite strain problem. The kinematics of the airfoil/hydrofoil is set up. The Consistent Orthogonal Basis Function Space is applied for the airfoil/hydrofoil structure. Given the airfoil/hydrofoil configuration and boundary conditions, the basis function space can be uniquely determined, such that the diagonal mass matrix is obtained accurately and the basis functions are very identical with the mode shape functions of the structure. In order to satisfy the displacement compatibility condition between adjacent layers of the airfoil/hydrofoil, the traction degree-of-freedom is also induced.The post buckling analysis is presented for the wing (hydrofoil) structure of the underwater glider. The water pressure is applied on the outer surface and the critical buckling pressure is calculated. The post buckling equilibrium path is also given. The results are verified with ANSYS. The present study of the buckling analysis of the airfoil/hydrofoil under water pressure is helpful to the design of underwater glider. 相似文献