• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2005.04a
• /
• pp.309-316
• /
• 2005

In this paper, we develop continuum-based design sensitivity analysis (DSA) methods using both direct differential method (DDM) and adjoint variable method (AVM) for non-shape design problems. The developed DSA method is further utilized for the topology design optimization of 3-dimensional structures. In numerical examples, the analytical DSA results are verified using finite difference ones. The topology optimization method yields very reasonable results in physical point of view.

• Geophysics and Geophysical Exploration
• /
• v.12 no.4
• /
• pp.281-288
• /
• 2009

Using analytical method, intrinsic absorption and scattering attenuation was lately separated from total quality factor $(Q^{-1}_t)$ on the seismic data of Korean Peninsula. However, we should use numerical method rather than the analytical method to consider depth dependent structure of scattering. The direct simulation Monte Carlo (DSMC) method, as a kind of the numerical method, is good option due to its extended availability from 1 to 3-dimensional model; but there is few study to use it. In this paper, we introduced the analytical method and the DSMC method, and compared the results of the two analysis applied to the isotropic scattering model. While the scattering attenuation coefficients $(\eta_s)$ are identical, the intrinsic absorption coefficients $(\eta_i)$ for the analytical method are larger than those for DSMC method and have large errors. In addition, the $(Q^{-1}_t)$ by the previous studies show closer to DSMC method than analytical method.

• Structural Engineering and Mechanics
• /
• v.38 no.2
• /
• pp.211-229
• /
• 2011

In this paper, we present a new explicit procedure using periodic cubic B-spline interpolation polynomials to solve linear and nonlinear dynamic equation of motion governing single degree of freedom (SDOF) systems. In the proposed approach, a straightforward formulation was derived from the approximation of displacement with B-spline basis in a fluent manner. In this way, there is no need to use a special pre-starting procedure to commence solving the problem. Actually, this method lies in the case of conditionally stable methods. A simple step-by-step algorithm is implemented and presented to calculate dynamic response of SDOF systems. The validity and effectiveness of the proposed method is demonstrated with four examples. The results were compared with those from the numerical methods such as Duhamel integration, Linear Acceleration and also Exact method. The comparison shows that the proposed method is a fast and simple procedure with trivial computational effort and acceptable accuracy exactly like the Linear Acceleration method. But its power point is that its time consumption is notably less than the Linear Acceleration method especially in the nonlinear analysis.

• Structural Engineering and Mechanics
• /
• v.83 no.3
• /
• pp.353-361
• /
• 2022

Linear Elastic Fracture Mechanics (LEFM) has been developed by applying stress analysis to determine the stress intensity factor (SIF, K). The finite element method (FEM) is widely used as a standard tool for evaluating the SIF for various crack configurations. The prediction accuracy can be achieved by applying an adaptive Delaunay triangulation combined with a FEM. The solution can be solved using either direct or iterative solvers. This work adopts the element-by-element preconditioned conjugate gradient (EBE-PCG) iterative solver into an adaptive FEM to solve the solution to heal problem size constraints that exist when direct solution techniques are applied. It can avoid the formation of a global stiffness matrix of a finite element model. Several numerical experiments reveal that the present method is simple, fast, and efficient compared to conventional sparse direct solvers. The optimum convergence criterion for two-dimensional LEFM analysis is studied. In this paper, four sample problems of a two-edge cracked plate, a center cracked plate, a single-edge cracked plate, and a compact tension specimen is used to evaluate the accuracy of the prediction of the SIF values. Finally, the efficiency of the present iterative solver is summarized by comparing the computational time for all cases.

• Proceedings of the Earthquake Engineering Society of Korea Conference
• /
• 2002.03a
• /
• pp.125-132
• /
• 2002

In this study, a numerical method is developed for nonlinear analysis for soil-pile-structure interaction system in time domain. Finite elements considering material nonlinearity are used for the near field and boundary elements for the far field. In the near field, frame elements are used for modeling a pile and plane-strain elements for surrounding soil and superstructure. In. the far field, boundary element formulation using the dynamic fundamental solution is adopted and coupled with the near field. Transformation of stiffness matrices of boundary elements into time domain is performed by inverse FFT. Stiffness matrices in the near field and far field are coupled. Newmark direct time integration method is applied. Developed soil-pile-structure interaction analysis method is verified with available literature and commercial code. Also, parametric studies by developed numerical method are performed. And seismic response analysis is performed using actual earthquake records.

• Bulletin of the Korean Mathematical Society
• /
• v.50 no.5
• /
• pp.1737-1751
• /
• 2013

An efficient and stable point collocation scheme based on a meshfree method is studied for the stationary incompressible Navier-Stokes equations. We describe the diffuse derivatives associated with the moving least square method. Using these diffuse derivatives, we propose a point collocation method to fit in solving the Navier-Stokes equations which improves the stability of the direct point collocation scheme. The convergence of the numerical solution is investigated from numerical examples. The driven cavity ow and the backward facing step ow are implemented for the reliability of the scheme. Also, the viscous ow on complicated geometry is successfully calculated such as the ow past a circular cylinder in duct.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.5
• /
• pp.1383-1391
• /
• 1996

The steady-state response for a coupled structure-acoustic-cavity systme has been investigated by numerical technique using a directly coupled finite element method(FEM) and Boundary Element Method(BEM) model. The Laplace tranformed matrix equations for the structure and the acoustic cavity are coupled directly satisfying the necessary equilibrium and compatibility conditions. The coupled FEM-BEM code is verified by comparing its prediction for an example with known analytical, numerical and experimental results. The example involves a coupled structure-acoustic-cavity system which is a box-type cavity with one end as experimentally excited pinned-pinned plate.

• Journal of Ocean Engineering and Technology
• /
• v.22 no.1
• /
• pp.70-74
• /
• 2008

It is important to reduce the vibration and noise of submarines and ships. For the purpose of noise reduction, various researches are actively being conducted on the employment of complex structures. However, in the case of numerical analysis for complex structures with damping materials, substantial errors can be generated by the absence of an exact damping model. Thus experimental model analysis is necessary for the verification of a numerical analysis for complex structures. In this research, vibration experiments are conducted in order to ascertain the vibration properties of cantilever beam attached damping materials. First, an initial value is obtained by using a direct linear method. Next, based on this initial value, the exact modal parameters of the cantilever beam are obtained by using the Newton-Raphson method.

• International Journal of Naval Architecture and Ocean Engineering
• /
• v.6 no.2
• /
• pp.297-306
• /
• 2014

The paper presents the numerical analysis of wave run-up over rubble-mound breakwaters covered by antifer units using a technique integrating Computer-Aided Design (CAD) and Computational Fluid Dynamics (CFD) software. Direct application of Navier-Stokes equations within armour blocks, is used to provide a more reliable approach to simulate wave run-up over breakwaters. A well-tested Reynolds-averaged Navier-Stokes (RANS) Volume of Fluid (VOF) code (Flow-3D) was adopted for CFD computations. The computed results were compared with experimental data to check the validity of the model. Numerical results showed that the direct three dimensional (3D) simulation method can deliver accurate results for wave run-up over rubble mound breakwaters. The results showed that the placement pattern of antifer units had a great impact on values of wave run-up so that by changing the placement pattern from regular to double pyramid can reduce the wave run-up by approximately 30%. Analysis was done to investigate the influences of surface roughness, energy dissipation in the pores of the armour layer and reduced wave run-up due to inflow into the armour and stone layer.

• Geomechanics and Engineering
• /
• v.35 no.2
• /
• pp.209-219
• /
• 2023

This paper presents numerical investigations into the particle crushing effect on the shear properties of gravel under direct shear condition. A novel particle crushing model was developed based on the octahedral shear stress criterion and fragment replacement method. A series of direct shear tests were carried out on unbreakable particles and breakable particles with different strengths. The evolutions of the particle crushing, shear strength, volumetric strain behavior, and contact force fabric during shearing were analyzed. It was observed that the number of crushed particles increased with the increase of the shear displacement and axial pressure and decreased with the particle strength increasing. Moreover, the shear strength and volume dilatancy were obviously decreased with particle crushing. The shear displacement of particles starting to crush was close to that corresponding to the peak shear stress got. Besides, the shear-hardening behavior was obviously affected by the number of crushed particles. A microanalysis showed that due to particle crushing, the contact forces and anisotropy decreased. The mechanism of the particle crushing effect on the shear strength was further clarified in terms of the particle friction and interlock.