• Geophysics and Geophysical Exploration
• /
• 제22권2호
• /
• pp.56-61
• /
• 2019

In this paper, we studied efficient modeling strategies when we simulate the 3D time-domain acoustic wave propagation using a cell-based finite difference method which can handle the variations of both P-wave velocity and density. The standard finite difference method assigns physical properties such as velocities of elastic waves and density to grid points; on the other hand, the cell-based finite difference method assigns physical properties to cells between grid points. The cell-based finite difference method uses average physical properties of adjacent cells to calculate the finite difference equation centered at a grid point. This feature increases the computational cost of the cell-based finite difference method compared to the standard finite different method. In this study, we used additional memory to mitigate the computational overburden and thus reduced the calculation time by more than 30 %. Furthermore, we were able to enhance the performance of the modeling on several media with limited density variations by using the cell-based and standard finite difference methods together.

• The Journal of the Acoustical Society of Korea
• /
• 제42권5호
• /
• pp.377-387
• /
• 2023

This paper presents accuracy of a method that directly estimates equivalent properties of a 1-3 piezocomposite for modeling it into the single phase homogeneous piezomaterial. This direct estimation method finds individual components of a material property matrix based on the piezoelectric constitutive equations, which represent mechanical and electrical behaviors and their couplings. Equivalent properties on a single 1-3 piezocomposite hydrophone are derived, and their accuracy depending on pairing of the constitutive equations is investigated by comparing them with finite element analysis for the whole domain. The accuracy is related to elastic characteristics of a matrix polymer, and the error is analyzed so that some guidelines for correct estimation are suggested. Fidelity of estimated properties and equivalent modeling is shown in a stave scale including hydrophones and surrounding acoustic structures as well, and reduced computational cost is verified.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• 제14권3호
• /
• pp.403-414
• /
• 1994

Nonlinear analysis of RC containment structure under thermal load and pressure is presented to trace the behaviour after an assumed LOCA. The temperature distribution varying with time through the wall thickness is determined by transient finite element analysis with the two time level scheme in time domain. The layered shell finite elements are used to represent the containment structures in nuclear power plants. Both geometric and material nonlinearities are taken into account in the finite element formulation. The constitutive relation of concrete is modeled according to Drucker-Prager yield criteria in compression. Tension stiffening model is used to represent the tensile behaviour of concrete including bond effect. The reinforcing bars are modeled by smeared layer at the location of reinforcements accounting elasto-plastic axial behaviors. The steel liner model under Von Mises yield criteria is adopted to represent elastic-perfect plastic behaviour. Geometric nonlinearity is formulated to consider the large displacement effect. Thermal stress components are determined by the initial strain concept during each time step. The temperature differential between any two consecutive time steps is considered as a load incremental. The numerical results from this study reveal that nonlinear temperature gradient based on transient thermal analysis will produces excessive large displacement. Nonlinear behavior of containment structures up to ultimate stage can be traced reallistically. The present study allows more realistic analysis of concrete containment structures in nuclear power plants.

• Journal of the Korean Society of Fisheries and Ocean Technology
• /
• 제38권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 the Earthquake Engineering Society of Korea
• /
• 제11권1호
• /
• pp.89-97
• /
• 2007

Parameters including the seismic sources and the elastic wave propagation characteristics were analysed using the observed ground motions from 12 Fukuoka region earthquakes. The Levenberg-Marquardt algorithm was applied to invert all the variables non-linearly and simultaneously with S wave energy in fiequency domain. Average stress drop of 12 events and local attenuation parameter $\kappa$ under seismic stations were estimated to about 79.2-bar and 0.043 respectively. Regional attenuation parameter, Qo and ${\eta}$, were also estimated to be about 248.1 and 0.558 respectively. Low value of Qo seems to caused by inhomogeneous tectonic characteristics between Japan island and southern Korean peninsula. $\kappa$ values are much higher than that characterizing EUS (Eastern United States) region, and nearly similar to that of WUS (Western Waited States) region. If the informations on site specific amplification of all the seismic stations are known, $\kappa$ values can be estimated more precisely. All the values including the seismic sources and the site and crustal scale propagation characteristics can be used as seismic design parameters.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• 제26권6호
• /
• pp.423-430
• /
• 2013

In this research, a paramteric study to account for the effect of interfacial strength and nanotube agglomeration on the elastoplastic behavior of carbon nanotube reinforced polypropylene composites is performed. At first, the elastoplastic behavior of nanocomposites is predicted from molecular dynamics(MD) simulations. By combining the MD simulation results with the nonlinear micromechanics model based on the Mori-Tanaka model, a two-step domain decomposition method is applied to inversely identify the elastoplastic behavior of adsorption interphase zone inside nanocomposites. In nonlinear micromechanics model, the secant moduli method combined with field fluctuation method is used to predict the elastoplastic behavior of nanocomposites. To account for the imperfect material interface between nanotube and matrix polymer, displacement discontinuity condition is applied to the micromechanics model. Using the elastoplastic behavior of the adsorption interphase zone obtained from the present study, stress-strain relation of nanocomposites at various interfacial bonding condition and local nanotube agglomeration is predicted from nonlinear micromechanics model with and without the adsorption interphase zone. As a result, it has been found that local nanotube agglomeration is the most important design factor to maximize reinforcing effect of nanotube in elastic and plastic behavior.

• Journal of Korean Tunnelling and Underground Space Association
• /
• 제17권2호
• /
• pp.75-90
• /
• 2015

As soil nails support a ground by the friction between nails and soils being reinforced, the length of soil nails is important factor for a ground stability. Thus, the soil nail length has to be accurately evaluated in order to secure a ground stability. The goal of this study is to suggest the applicability of the non-destructive method as the basic research for the evaluation of the soil nail length. First, the elastic and electromagnetic waves are adopted to select an applicable method for the soil nails connected by the coupler. Test results show that while the ultrasonic waves are not detected due to the coupler, the electromagnetic waves are free for the influence of the coupler. Second, electromagnetic waves are measured for combined soil nails with the length of 1 m~15 m for the investigation of the characteristics of electromagnetic waves. The travel time of the electromagnetic wave increases with an increase in the soil nail length. In addition, the ground cable is used to apply the electromagnetic waves to pre-installed soil nails. Test results show that the travel time of the electromagnetic wave by using the ground cable increases with an increase in soil nail length. This study demonstrates that the electromagnetic wave may be a promising method for the evaluation of the soil nail length.

• Journal of the Korean Society for Nondestructive Testing
• /
• 제20권2호
• /
• pp.138-149
• /
• 2000

Various modeling techniques for ultrasonic wave propagation and scattering problems in finite solid media are presented. Elastodynamic boundary value problems in inhomogeneous multi-layered plate-like structures are set up for modal analysis of guided wave propagation and numerically solved to obtain dispersion curves which show propagation characteristics of guided waves. As a powerful modeling tool to overcome such numerical difficulties in wave scattering problems as the geometrical complexity and mode conversion, the Boundary Element Method(BEM) is introduced and is combined with the normal mode expansion technique to develop the hybrid BEM, an efficient technique for modeling multi mode conversion of guided wave scattering problems. Time dependent wave forms are obtained through the inverse Fourier transformation of the numerical solutions in the frequency domain. 3D BEM program development is underway to model more practical ultrasonic wave signals. Some encouraging numerical results have recently been obtained in comparison with the analytical solutions for wave propagation in a bar subjected to time harmonic longitudinal excitation. It is expected that the presented modeling techniques for elastic wave propagation and scattering can be applied to establish quantitative nondestructive evaluation techniques in various ways.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• 제19권4호
• /
• pp.429-440
• /
• 2006

This paper comprises two specific objectives. The first is to examine the applicability of ordinary kriging interpolation(OK) to the p-adaptivity of the finite element method that is based on variogram modeling. The second objective Is to present the adaptive procedure by the hierarchical p-refinement in conjunction with a posteriori error estimator using the modified S.P.R. (superconvergent patch recovery) method. The ordinary kriging method that is one of weighted interpolation techniques is applied to obtain the estimated exact solution from the stress data at the Gauss points. The weight factor is determined by experimental and theoretical variograms for interpolation of stress data apart from the conventional interpolation methods that use an equal weight factor. In the p-refinement, the analytical domain has to be refined automatically to obtain an acceptable level of accuracy by increasing the p-level non-uniformly or selectively. To verify the performance of the modified S.P.R. method, the new error estimator based on limit value has been proposed. The validity of the proposed approach has been tested with the help of some benchmark problems of linear elastic fracture mechanics such as a centrally cracked panel, a single edged crack, and a double edged crack.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• 제27권1호
• /
• pp.9-16
• /
• 2014

Using a level set method and topological derivatives, a topological shape optimization method that is independent of an initial design is developed for linearly elastic structures. In the level set method, the initial domain is kept fixed and its boundary is represented by an implicit moving boundary embedded in the level set function, which facilitates to handle complicated topological shape changes. The "Hamilton-Jacobi(H-J)" equation and computationally robust numerical technique of "up-wind scheme" lead the initial implicit boundary to an optimal one according to the normal velocity field while minimizing the objective function of compliance and satisfying the constraint of allowable volume. Based on the asymptotic regularization concept, the topological derivative is considered as the limit of shape derivative as the radius of hole approaches to zero. The required velocity field to update the H-J equation is determined from the descent direction of Lagrangian derived from optimality conditions. It turns out that the initial holes are not required to get the optimal result since the developed method can create holes whenever and wherever necessary using indicators obtained from the topological derivatives. It is demonstrated that the proper choice of control parameters for nucleation is crucial for efficient optimization process.