利用有限元方法求解双曲型缓坡方程

本文提出了一种双曲型缓坡方程的有限元计算方法 ,在建立有限元积分方程时通过在造波线处加入脉动源项来实现内部造波 ,并在开边界处利用阻尼层吸波 ,减少了在边界处由于数值处理引起的误差。数值计算结果与实测值吻合良好。本方法可用于大区域波浪场的计算中     

应用比例边界有限元法求解多种二维浮体水动力特性

比例边界有限元法(SBFEM)是一种半解析的数值方法,比完全数值方法具有更高的精度,该方法结合了有限元和边界元的优点,采用相对少的剖分单元就可以得到较高精度的模拟结果。通过改变有限子域内部比例中心的位置,使这种方法可以应用到多种形式浮体在波浪作用下的水动力特性的计算中。同时还给出了各种形式浮体的波浪力及反射、透射系数的数值结果,并与边界元方法(BEM)计算结果和特征函数展开方法得到解析解进行了比较,均吻合良好。研究表明比例边界有限元不仅可以计算矩形的浮体结构,而且对于多种结构形式的浮体都可以计算,这为多种结构形式浮体的水动力分析提供了一个可行的方法。     

钢悬链线立管与海底相互作用弹性基础梁模拟方法研究

采用弹性基础梁模型,模拟钢悬链线立管(steel catenary riser,SCR)与海床土的相互作用,推导了弹性基础梁模型的有限元公式,基于该模型编制了钢悬链式立管的动力分析程序,在此基础上,分析了在浮体运动及环境载荷作用下,钢悬链线立管触地点区域的动力响应特征,并与弹簧支撑模型进行了比较。研究表明,弹性基础梁模型对单元长度的敏感性较低,可采用较长的单元,从而减少计算量。     

Finite Element Method for Design of Reinforced Concrete Offshore Platforms

A design method of reinforced concrete (R. C.) offshore platforms with nonlinear finite element analysis is proposed. According to the method, a computer program is developed. In this program nonlinear constitutive relationships and strength criteria of concrete and steel bars are included, and the progressive cracking and crushing of the concrete are taken into account. Based on the stress distribution obtained by the nonlinear finite element analysis, the amount of reinforcement in the control sections can be computed and adjusted automatically by the program to satisfy the requirement of the design. The amount of reinforcement required in the control sections, which are obtained with the nonlinear finite element analysis, is agreeable to that obtained in the experiment. This shows that the design method of R. C. offshore platform with the nonlinear finite element method proposed by the authors is reliable for practical use.     

Scaled Boundary Finite Element Analysis of Wave Passing A Submerged Breakwater

The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique combining the advantage of the finite element method （FEM） and the boundary element method （BEM） with its unique properties. In this paper, the SBFEM is used for computing wave passing submerged breakwaters, and the reflection coeffcient and transmission coefficient are given for the case of wave passing by a rectangular submerged breakwater, a rigid submerged barrier breakwater and a trapezium submerged breakwater in a constant water depth. The results are compared with the analytical solution and experimental results. Good agreement is obtained. Through comparison with the results using the dual boundary element method （DBEM）, it is found that the SBFEM can obtain higher accuracy with fewer elements. Many submerged breakwaters with different dimensions are computed by the SBFEM, and the changing character of the reflection coeffcient and the transmission coefficient are given in the current study.     

水下加筋圆柱壳结构声散射特性研究

采用有限元/边界元(FEM/BEM)方法对水下加筋圆柱壳结构的声散射特性进行分析。考虑入射声场激发弹性结构振动产生二次辐射声场,总声场为入射场与散射场之和。基于Mindlin理论建立结构受迫振动的有限元模型,考虑入射声场为简谐激励,采用模态叠加法求出结构振动表面速度。根据声学边界元相关理论,利用有限元方法计算得出的弹性结构表面动力相应,可求出弹性结构的散射声场。研究结果表明,在同一波数下,不加筋结构后向散射强于加筋结构;加筋结构的散射指向性图案分叉较为明显;同种结构在波数增加的时候后向散射强度逐渐增大,前向散射强度逐渐减弱。     

Time stepping solutions of the two-dimensional nonlinear wave radiation problem

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.     

预修正快速傅里叶变换高阶边界元方法在波浪对物体作用问题中的应用

边界元方法被广泛应用于波浪对海上婕筑物作用领域,但由于传统边界元方法的存储量和计算量均为未知量的平方量级,很难满足大范围多未知量问题的计算需要.本文基于高阶边界元方法,应用预修正快速傅里叶变换方法,使计算量与存储量分别降低至O(NlogN)和O(N)量级,并得到一个连续的压强分布以适应结构设计的要求,同时可以通过使用满...     

Dynamic modeling of cable towed body using nodal position finite element method

This paper analyses nonlinear dynamics of cable towed body system. The cable has been modeled and analyzed using a new nodal position finite element method, which calculates the position of the cable directly instead of the displacement by the existing finite element method. The newly derived nodal position finite element method eliminates the need of decoupling the rigid body motion from the total motion, where numerical errors arise in the existing nonlinear finite element method, and the limitation of small rotation in each time step in the existing nonlinear finite element method. The towed body is modeled as a rigid body with six degrees of freedom while the tow ship motion is treated as a moving boundary to the system. A special procedure has been developed to couple the cable element with the towed body. The current approach can be used as design tool for achieving improved directional stability, maneuverability, safety and control characteristics with the cable towed body. The analysis results show the elegance and robustness of the proposed approach by comparing with the sea trial data.     

Nonlinear response of membranes to ocean waves using boundary and finite elements

The behavior of a highly deformable membrane to ocean waves was studied by coupling a nonlinear boundary element model of the fluid domain to a nonlinear finite element model of the membrane. The hydrodynamic loadings induced by water waves are computed assuming large body hydrodynamics and ideal fluid flow and then solving the transient diffraction/radiation problem. Either linear waves or finite amplitude waves can be assumed in the model and thus the nonlinear kinematic and dynamic free surface boundary conditions are solved iteratively. The nonlinear nature of the boundary condition requires a time domain solution. To implicitly include time in the governing field equation, Volterra's method was used. The approach is the same as the typical boundary element method for a fluid domain where the governing field equation is the starting point. The difference is that in Volterra's method the time derivative of the governing field equation becomes the starting point.The boundary element model was then coupled through an iterative process to a finite element model of membrane structures. The coupled model predicts the nonlinear interaction of nonlinear water waves with highly deformable bodies. To verify the coupled model a large scale test was conducted in the OH Hinsdale wave Research Laboratory at Oregon State University on a 3-ft-diameter fabric cylinder submerged in the wave tank. The model data verified the numerical prediction of the structure displacements and of the changes in the wave field.The boundary element model is an ideal modeling technique for modeling the fluid domain when the governing field equations is the Laplace equation. In this case the nonlinear boundary element model was coupled with a finite element model of membrane structures, but the model could have been coupled with other finite element models of more rigid structures, such as a pontoon floating breakwater.     

近岸波浪折射-绕射-破波耗散联合模式的有限元数值研究

建立了近岸波浪折射-绕射-破波耗散的有限元数值模式。采用的有限元方法为改进的混合元法,其中外域开边界条件得到改进,内域有限元为伽廖金有限元。用理论解检验了所建立的数值模式,并将该模式应用到一个模型港湾。     

桶型基础结构与土壤相互作用的有限元—无限元—接触元耦合线性分析

桶型基础是一种新型的海洋基础结构型式，在负压作用下的下沉过程中，桶型基础结构与地基产生了相互作用，考虑结构与土壤的相互作用的强度分析为桶型基础的设计提供依据，本文将有限元－无限地－接触元相耦合的线性数值分析方法引入桶型基础的结构与土壤之间的相互作用的强度分析，采用接触单元考虑结构与土壤之间的错动滑移及拉裂，采用无限元更有效地反映地基无限域的远场效应，还推导出三节点，七节点，八节点三种新形式的单向无限地坐标映射函数，应用该方法的计算结果与模型试验结果比较，取得较好的效果。     

Numerical Investigation on Dynamics of the Tendon System of A TLP by Applying Absolute Nodal Coordinate Formulation

In the present study, the dynamics of the tendon system of a tension-leg platform(TLP) is investigated through the absolute nodal coordinate formulation(ANCF). Based on the energy conversion principle, the stiffness, generalized elastic force, external load and mass matrices of the element are deduced to perform the element assembling by using the finite element method. Then the motion equation of the tendon/riser is established. In this study, the TLP in the International Ship Structures Committee(ISSC) model under the first and second wave forces is considered as the case study. The simulation is performed in the MATLAB environment. Moreover, the accuracy and reliability of the programs are verified for cases of beam model with theoretical solutions. It is found that the motion response of tendons is affected by the TLP movement and environmental load, simultaneously. Then, the motion response is calculated using the SESAM software and exported as the boundary of ANCF tendons. Finally, the static and dynamic characteristics of the four tendons of ISSC TLP are analyzed systematically by the ANCF method.Performed analysis proves the effectiveness and feasibility of the ANCF method. It is concluded that the proposed method is a powerful scheme for calculating the dynamics of tendon/riser in the field of ocean engineering.     

Nonlinear Finite Element Analysis of Ocean Cables

This study has focused on developing numerical procedures for the dynamic nonlinear analysis of cable structures subjected to wave forces and ground motions in the ocean. A geometrically nonlinear finite element procedure using the isoparamnetrie curved cable element based on the Lagrangian folmulation is briefly summarized. A simple and accurate method to determine the initial equilibrium state of cable systems associated with self-weights, buoyancy and the motion of end points is presented using the load incremental method combined with penalty method. Also the Newmark method is used for dynamic nonlinear analysis of ocean cables. Numerical examples are presented to validate the present numerical method.     

Finite Element Analysis of Deformed Legs of Offshore Platform Structures

The element stiffness matrix of the equivalent beam or pipe element of the deformed leg of the platform is derived by the finite element method. The stresses and displacements of some damaged components are calculated, and the numeri-cal solutions agree well with those obtained by the fine mesh finite element method. Finally, as an application of this method, the stresses of some platform structures are calculated and analyzed.     

具有裂纹损伤桩腿的海洋石油平台有限元分析

采用有限元分析方法和断裂力学方法,将含裂纹构件的裂纹以线弹簧处理,建立了带有裂纹损伤桩腿的等效梁和管单元的单元刚度矩阵的计算方法。对检验方法的正确性,使用该方法对几种损伤构件进行了强度计算,并与有限元细化网格计算结果和实验结果进行比较,最后将等效单元计算程序接入大型有限元程序ＳｕｐｅｒＳＡＰ,对受损后的平台整体结构进行了应力计算和强度分析。     

基于天绘一号卫星多光谱影像的海岸线自动提取方法

针对常规人工海岸线测量方法耗费人力物力、效率低、工作环境相对危险,难以快速反映海岸线动态变化的问题,提出一种基于天绘一号卫星多光谱影像的海岸线自动提取与分类方法。所提方法首先结合NDWI指数和水平集模型实现全自动水边线提取,然后采用面向对象分割技术和空间拓扑关系提取临海地物基元,最后结合遥感影像海岸线解译标志建立决策树完成临海基元的分类与海岸线提取。实验表明,海岸线提取结果能够作为1∶2.5万海图测制与更新的参考数据。     

Transverse vibrations of a rectangular cantilever plate with thickness varying in a discontinuous fashion

A numerical solution is obtained for the title problem by means of the popular finite element method. An experimental investigation of the problem is also presented for plates of square planform. The agreement between theoretical and experimental values is very good in the case of the fundamental frequency, and reasonably good when the second and third normal modes of transverse vibration are considered. It is concluded that beneficial effects are obtained in the sense that one generates a lighter structural element with higher fundamental frequency of transverse vibration. Hence, dynamic stiffening is achieved.     

Lateral Bearing Capacity of Multi-Bucket Foundation in Soft Ground

The bucket foundation is a new type of foundation for offshore application to intermediate-depth waters. It has advantages over conventional ones. However, there is no consensus in the analysis and design of this type of foundation. In this paper, the lateral bearing capacity and the failure mechanism of multi-bucket foundations are studied with different connection stiffness and bucket spacing by use of a three-dimensional finite element method. Based on the numerical analysis results, a limit analysis method of plasticity for evaluating the lateral bearing capacity of large-spacing multi-bucket foundation with rigid connection in soft ground is proposed. This method provides a simple procedure that gives results comparable to those from the finite element analyses.     

Numerical Analyses of Caisson Breakwaters on Soft Foundations Under Wave Cyclic Loading

A caisson breakwater is built on soft foundations after replacing the upper soft layer with sand. This paper presents a dynamic finite element method to investigate the strength degradation and associated pore pressure development of the intercalated soft layer under wave cyclic loading. By combining the undrained shear strength with the empirical formula of overconsolidation clay produced by unloading and the development model of pore pressure, the dynamic degradation law that describes the undrained shear strength as a function of cycle number and stress level is derived. Based on the proposed dynamic degradation law and M-C yield criterion, a dynamic finite element method is numerically implemented to predict changes in undrained shear strength of the intercalated soft layer by using the general-purpose FEM software ABAQUS, and the accuracy of the method is verified. The effects of cycle number and amplitude of the wave force on the degradation of the undrained shear strength of the intercalated soft layer and the associated excess pore pressure response are investigated by analyzing an overall distribution and three typical sections underneath the breakwater. By comparing the undrained shear strength distributions obtained by the static method and the quasi-static method with the undrained shear strength distributions obtained by the dynamic finite element method in the three typical sections, the superiority of the dynamic finite element method in predicting changes in undrained shear strength is demonstrated.