• Journal of Power System Engineering
• /
• v.14 no.6
• /
• pp.75-82
• /
• 2010

In this study, the topics for free-surface wave simulation, nonlinear hydrodynamic force, and the critical resonance frequency of so-called ${\tau}=U{\omega}/g$=1/4 are discussed. A high-order spectral/boundary element method is newly adapted as an efficient numerical tool. This method is one of the most efficient numerical methods by which the nonlinear gravity waves can be simulated and hydrodynamic forces also can be calculated in time domain. In this method, the velocity potential is expressed as the sum of surface potential and body potential. Then, surface potential is solved by using the high-order spectral method and body potential is solved by using the high-order boundary element method. By the combination of these two methods, the wave-making problems by a submerged sphere oscillating with forward speed under the free-surface are solved in time domain.

• Computational Structural Engineering
• /
• v.10 no.1
• /
• pp.185-191
• /
• 1997

Beam - and arch-like structures are two-dimensional bodies characterized by the fact of small thickness compared to the length of structures. Owing to this geometric feature, linear displacement approximations through the thickness such as Kirchhoff and Reissner-Mindlin theories which are more accessible one dimensional problems have been used. However, for accurate analysis of the behavior in the regions where the state of stresses is complex, two-dimensional linear elasicity or relatively high order of thickness polynomials is required. This paper analyses accuracy according to the order of thickness polynomials and introduces a technique for model combination for which several different polynomial orders are mixed in a single structure.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.14 no.6
• /
• pp.1289-1298
• /
• 1994

The new p-version crack model is proposed to estimate the stress intensity factors of the thick cracked plate under flexure. The proposed model is based on high order theory and $C^{\circ}$-plate element including shear deformation. The displacements fields are defined by integrals of Legendre polynomials which can be classified into three groups such as basic mode, side mode and internal mode. The computer implementation allows arbitrary variations of p-level Up to a maximum value of 10. The stress intensity factors are computed by virtual crack extention approach. The effects of ratios of thickness to crack length(h/a), crack length to width(a/W) and boundary conditions are investigated. Very good agreement with the existing solution in the literature are shown for the uncracked plate as well as the cracked plate.

• Journal of the Society of Naval Architects of Korea
• /
• v.29 no.1
• /
• pp.135-142
• /
• 1992

In the vibration analysis or submerged of floating bodies such as ship and offshore structures, the coupled system between structure and fluid satisfying the compatibility conditions on the wetted surface should be considered. It is well known that the hydroelastic analysis of structural vibration in contact with fluid can be solved by applying the finite element method to structure and the boundary element method to fluid domain. However such an approach is impractical, because fluid added mass matrix is fully coupled on whole wetted surface. To overcome this shortcoming, an efficient approach based on reanalysis scheme is proposed in this paper. The proposed method can be applied for cases with higher modes lacking 3-D reduction factor J as well as beam-like modes of marine structures. It is well known the traditional method using 2-D added mass and J-factor is good only for beam-like modes with reliable J values. The validity and the calculation efficiency of the proposed method are proved with numerical examples.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.32 no.4
• /
• pp.223-231
• /
• 2019

A continuum-based design sensitivity analysis(DSA) method is developed for electrostatic problems. To consider high order objective functions, we use 9-node finite element basis functions for analysis and DSA methods. As the design variables are parameterized with B-spline functions, smooth boundary variations are naturally obtained. To solve mesh entanglement problems during the optimization process, a mesh regularization scheme is employed. By minimizing the Dirichlet energy functional, mesh uniformity can be automatically achieved. In numerical examples for maximizing dielectrophoresis forces, the numerical results are compared with well-known electrode geometries and the obtained characteristics are discussed.

• Journal of Wind Energy
• /
• v.14 no.2
• /
• pp.26-33
• /
• 2023

In this study, the global response characteristics of floater and mooring for floating offshore wind turbine with a detachable mooring system are performed. Global motion and structural response result extracted from the coupled motion analysis of 10MW DTU floating offshore wind turbine with detachable mooring system modeled by high-order boundary element model and finite element mesh, were used to study the characteristics of tension on mooring lines subjected to three different types of ocean loads. Breaking limit of mooring line characterized by wind, current and wave load has a major effect on the distribution of mooring tension found in time domain analysis. Based on the numerical results of coupled motion analysis, governing equation for calculating the motion response of a floater under ocean loads, and excitation force and surge motion and tension respectively are presented using excursion curve. It is found that the response of floater is reliable and accurate for calculating the tension distributions along the mooring lines under complex loadings. This means that the minimun breaking limit of mooring system satisfied a design criteria at ultimate ocean environmental loading condtions.

• Journal of the Society of Naval Architects of Korea
• /
• v.38 no.4
• /
• pp.39-47
• /
• 2001

In this paper, the hydroelastic behavior of a mat-type large floating structure is analyzed in time domain by using mode superposition method. The time-memory function is estimated by Fourier transforming the wave damping coefficients, which are computed by a higher-order boundary element method based on potential theory. Meanwhile, the structural response is obtained by time integrating the eigenmodes of the structure. Numerical examples are made for three test cases on the scaled model of a mat-type large floating structure ; weight pull-up case, weight drop case and weight moving case. In all three cases, the numerical results coincide well with experimental data.

• Journal of the Society of Naval Architects of Korea
• /
• v.35 no.4
• /
• pp.27-37
• /
• 1998

The sloshing is very important in a safe transport of the liquid cargo by a ship. With the increasing number of supertanker and LNG carriers, this problem has become increasingly more important. In order to study the magnitude and characteristics of impact pressures due to sloshing, experiments ware performed with a rectangular tank and compared with numerical results. Structural responses of tank wall under impulsive pressures were measured. Structural vibrations induced by the sloshing load were analysed by including hydroelastic erects in terms of added mass and damping. To check the validity of the numerical model, the natural frequencies of plate in air and water were compared with measurements, and a good agreement was found. In the case that a plate vibrates under impulsive loads, the pressure on the flexible plate is larger than that on the rigid plate without hydroelastic effects, which was confirmed experimentally. The frequency of oscillatory pressure as well as accel%pion coincides with the natural frequency of plate in water.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.24 no.1
• /
• pp.1-7
• /
• 2011

In this paper, based on an isogeometric approach, we have developed a shape design optimization method for plane elasticity problems subjected to design-dependent loads. The conventional shape optimization using the finite element method has some difficulties in the parameterization of geometry. In an isogeometric analysis, however, the geometric properties are already embedded in the B-spline basis functions and control points so that it has potential capability to overcome the aforementioned difficulties. The solution space for the response analysis can be represented in terms of the same NURBS basis functions to represent the geometry, which yields a precise analysis model that exactly represents the normal and curvature depending on the applied loads. A continuum-based isogeometric adjoint sensitivity is extensively derived for the plane elasticity problems under the design-dependent loads. Through some numerical examples, the developed isogeometric sensitivity analysis method is verified to show excellent agreement with finite difference sensitivity.

• Journal of the Society of Naval Architects of Korea
• /
• v.30 no.3
• /
• pp.116-126
• /
• 1993

The finite element method(FEM) has been commonly used for structural dynamic analysis. However, the direct global application of FEM to large complex structures such as ships and offshore structures requires considerable computational efforts, and remarkably more in structural dynamic optimization problems. Adoption of the component-mode synthesis method is an efficient means to overcome the above difficulty. Among three classes of the component-mode synthesis method, the free-interface mode method is recognized to have the advantages of better computational efficiency and easier implementation of substructures' experimental results, but the disadvantage of lower accuracy in analytical results. In this paper, an advanced method to improve the accuracy in the application of the free-interface mode method for the vibration analysis of large complex structures is presented. In order to compensate the truncation effect of the higher modes of substructures in the synthesis process, both residual inertia and stiffness effects are taken into account and a frequency shifting technique is introduced in the formulation of the residual compliance of substructures. The introduction of the frequency shrift ins not only excludes cumbersome manipulation of singular matrices for semi-definite substructural systems but gives more accurate results around the specified shifting frequency. Numerical examples of typical structural models including a ship-like two dimensional finite element model show that the analysis results based on the presented method are well competitive in accuracy with those obtained by the direst global FEM analysis for the frequencies which are lower than the highest one employed in the synthesis with remarkably higher computational efficiency and that the presented method is more efficient and accurate than the fixed-interface mode method.