• Structural Engineering and Mechanics
• /
• v.53 no.4
• /
• pp.703-723
• /
• 2015

In this paper optimization of volume fraction distribution in a thick hollow cylinder with finite length made of two-dimensional functionally graded material (2D-FGM) and subjected to steady state thermal and mechanical loadings is considered. The finite element method with graded material properties within each element (graded finite elements) is used to model the structure. Volume fractions of constituent materials on a finite number of design points are taken as design variables and the volume fractions at any arbitrary point in the cylinder are obtained via cubic spline interpolation functions. The objective function selected as having the normalized effective stress equal to one at all points that leads to a uniform stress distribution in the structure. Genetic Algorithm jointed with interior penalty-function method for implementing constraints is effectively employed to find the global solution of the optimization problem. Obtained results indicates that by using the uniform distribution of normalized effective stress as objective function, considerably more efficient usage of materials can be achieved compared with the power law volume fraction distribution. Also considering uniform distribution of safety factor as design criteria instead of minimizing peak effective stress affects remarkably the optimum volume fractions.

• Journal of Ship and Ocean Technology
• /
• v.4 no.2
• /
• pp.1-12
• /
• 2000

The present study concentrates on the weak-scatterer hypothesis for the nonlinear wave-body interaction problems. In this method, the free surface boundary conditions are linearized on the incoming wave profile and the exact body motion is applied. The considered problems are the diffraction problem near a circular cylinder and the ship response in oblique waves. The numerical method of solution is a Rankine panel method. The Rankine panel method of this study adopts the higher-order B spline basis function for the approximation of physical variables. A modified Euler scheme is applied for the time stepping, which has neutral stability. The computational result shows some nonlinear behaviors of disturbance waves and wave forces. Moreover, the ship response shows very close results to experimental data.

• Journal of the Society of Naval Architects of Korea
• /
• v.32 no.1
• /
• pp.42-57
• /
• 1995

An efficient and accurate scheme has been constructed by taking advantages of the hi-quadratic spline scheme and the higher-order boundary element method selectively depending on computation domains. Boundary surfaces are represented by 8-node boundary elements to describe curved surfaces of a ship and its neighboring free surface more accurately. The variation of the velocity potential complies with the characteristics of the 8-node element on the body surface. But on the free surface, it is assumed to follow that of the hi-quadratic spline scheme. By which, the free surface solution is free from numerical damping and has better numerical dispersion property. As numerical examples, steady and unsteady Neumann-Kelvin problems are considered. Numerical results for a submerged spheroid, Series 60($C_B=0.6$) and a modified support the proposed method. Finally, a new upstream radiation condition is derived using a wave equation operator in order to deal with problems for subcritical reduced frequency. The relevance of this operator has been confirmed in the case of unsteady Kelvin source potential.

• Bulletin of the Society of Naval Architects of Korea
• /
• v.27 no.3
• /
• pp.38-46
• /
• 1990

A numerical solution method of the boundary integral equation for axisymmetric potential flows is presented. Those are represented by ring source and ring vorticity distribution. Strengths of ring source and ring vorticity are approximated by linear functions of a parameter $\zeta$ on a segment. The geometry of the body is represented by a cubic B-spline. Limiting integral expressions as the field point tends to the surface having ring source and ring vorticity distribution are derived upto the order of ${\zeta}ln{\zeta}$. In numerical calculations, the principal value integrals over the adjacent segments cancel each other exactly. Thus the singular part proportional to $\(\frac{1}{\zeta}\)$ can be subtracted off in the calculation of the induced velocity by singularities. And the terms proportional to $ln{\zeta}$ and ${\zeta}ln{\zeta}$ can be integrated analytically. Thus those are subtracted off in the numerical calculations and the numerical value obtained from the analytic integrations for $ln{\zeta}$ and ${\zeta}ln{\zeta}$ are added to the induced velocity. The four point Gaussian Quadrature formula was used to evaluate the higher order terms than ${\zeta}ln{\zeta}$ in the integration over the adjacent segments to the field points and the integral over the segments off the field points. The root mean square errors, $E_2$, are examined as a function of the number of nodes to determine convergence rates. The convergence rate of this method approaches 2.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• v.38 no.2
• /
• pp.110-118
• /
• 2002

The element free Galerkin method has been suggested by Belyschko, etc. it is applied for linear elastic analysis of solid problems and its convergence characteristics and stability are shown according to the weight function, influence domain and scattered points. The various numerical examples are performed to check the efficiency of 1D EFG and 2D EFG program by changing factors. As a result it have the best results when it used the cubic spline weight function and the scaling parameter d$_{max}$=2. These programs were developed by mixed language programming method using Visual Basic and the C language. so it is fast and efficient. and visually shown the result.t.

• Journal of Navigation and Port Research
• /
• v.46 no.6
• /
• pp.471-482
• /
• 2022

This study aims to investigate the characteristics of wave generation in a deep ocean engineering basin and to develop a meta-model of the transfer function of the wavemaker that reflects the geometric characteristics of the deep ocean engineering basin. To this end, the two-dimensional frequency domain boundary element method was applied to achieve an efficient analysis that reflects the geometric characteristics of the deep ocean engineering basin. The developed numerical method was validated through comparison with the analytical solution. Numerical analyses were conducted for the boundary value problem of the wavemaker according to various periods and the positions of the movable bottom. The numerical results were used to investigate the effect of the geometric characteristics of the deep ocean engineering basin on the transfer function of the wavemaker, and the effect of depth on wave generation was checked by changing the position of the movable bottom. To efficiently utilize the various results of the boundary element method, a meta-model, an approximate model of the transfer function of the wave maker, was developed using a thin plate spline interpolation model. The validity of the developed meta-model was confirmed through a comparison of the results of the model tests.

• 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.