• Restorative Dentistry and Endodontics
• /
• v.33 no.3
• /
• pp.184-197
• /
• 2008

This study was to investigate the influence of combining composite resins with different elastic modulus, and occlusal loading condition on the stress distribution of restored notch-shaped non-carious cervical lesion using 3D finite element (FE) analysis. The extracted maxillary second premolar was scanned serially with Micro-CT. The 3D images were processed by 3D-DOCTOR. ANSYS was used to mesh and analyze 3D FE model. A notch-shaped cavity was modeled and filled with hybrid, flowable resin or a combination of both. After restoration, a static load of 500N was applied in a point-load condition at buccal cusp and palatal cusp. The stress data were analyzed using analysis of principal stress. Results showed that combining method such that apex was restored by material with high elastic modulus and the occlusal and cervical cavosurface margin by small amount of material with low elastic modulus was the most profitable method in the view of tensile stress that was considered as the dominant factor jeopardizing the restoration durability and promoting the lesion progression.

• Computational Structural Engineering
• /
• v.3 no.1
• /
• pp.59-70
• /
• 1990

The p-version of the finite element method is a new approach to finite element analysis in which the partition of the domain is held fixed while the degree p of approximating piecewise polynomials is increased. In this paper, the focus is on computer implementation of a new hierarchic p-convergence shell model based on blend mapping functions. Its rigid-body modes, round-off error, and convergence characteristics are investigated.

• The Journal of Engineering Geology
• /
• v.19 no.2
• /
• pp.205-215
• /
• 2009

The conventional numerical models to analyze flow in subsurface porous media under the transient state usually generate numerical oscillation and unstability due to local flux domain for critical cases such as infiltration into initially dry soil during rainfall period. In this case, it is required refined mesh and small time step, but it decrease efficiency of computation. In this study, numerical unstability in discontinuity domain is removed by applying particle tracking algorithm to simulate unsteady subsurface flow with inflow boundary condition. Finally the hybrid LE FEM improving numerical stability is proposed. The hypothetical domains with unsteady uniform and nonuniform flow field were used to demonstrated algorithm verification. In comparison with analytic solution, we obtained reasonable results and conducted simulation of hypothetical 3-D recharge/pumping area. The proposed algorithm can simulate saturated/unsaturated porous media with more practical problems and will greatly contribute to accuracy and stability of numerical computation.

• Computational Structural Engineering
• /
• v.9 no.2
• /
• pp.133-142
• /
• 1996

In this paper, two different techniques for mixed-mode type engineering fracture mechanics are investigated to estimate the stress intensity factors by using p-version finite element model. These two techniques are displacement extrapolation with COD and CSD method and J-integral with decomposition method. By decomposing the displacement field obtained from p-version of finite element analysis into symmetric and antisymmetric displacement fields with respect to the crack line, Mode-I and Mode-II stress intensity factors can be determined using aforementioned techniques. The example problems for validating the proposed techniques are centrally and centrally oblique cracked panels under tension. The numerical results associated with the variation of oblique angle and the ratio of crack length and panel width (a /W ratio) are compared with those by theoretical values and empirical solutions in literatures. Very good agreements with the existing solutions are shown.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.41 no.2
• /
• pp.91-97
• /
• 2013

Pyrolysis and surface recession of charring composites are progressed primarily in the thickness direction. The unmixing-mixing scheme is applied to describe the in-plane and through-thickness behaviors of porous composites. The extended unmixing-mixing equations are based on transverse isotropy of unidirectionally fiber-reinforced composites. The strain components of gas pressure in pores, thermal expansion, and chemical shrinkage are included in the constitutive model. By analyzing micromechanical representative volume elements of porous composites, the validity of the derived equations are examined.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.34 no.3
• /
• pp.167-174
• /
• 2021

Recent investigation into the integrity of nuclear containment buildings has highlighted the importance of developing an elaborate diagnostic method to evaluate the distribution and size of cavities inside concrete walls. As part of developing such a method, this paper presents a finite element approach to modeling elastic waves propagating in the containment building walls of a nuclear power plant. We introduce a perfectly matched layer (PML) wave-absorbing boundary to limit the large-scale nuclear containment wall to the region of interest. The formulation results in a semi-discrete form with symmetric damping and stiffness matrices. The transient elastic wave equations for a mixed unsplit-field PML were solved for displacement and stresses in the time domain. Numerical results show that the sensitivity of displacement, velocity, acceleration, and stresses is large depending on the size and location of the cavity. The dynamic response of the wall slightly differs depending on the existence of the containment liner plate. The results of this study can be applied to a full-waveform inversion approach for characterizing cavities inside a containment wall.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.31 no.6B
• /
• pp.507-514
• /
• 2011

A numerical model has been developed by employing a finite element method to simulate the depth-averaged 2-D dispersion of the heat pollutant, which is an important pollutant material in natural streams. Among the finite element methods, the Streamline Upwind/Petrov Galerkin (SUPG) method was applied. Also both linear and quadratic elements can be applied so that irregular river boundaries can be easily represented. To show the movement of heat pollutants, the reaction term describing heat transfer was represented as an equation in which sink/source term is proportional to the difference between the equilibrium temperature and water surface temperature. The equation was expressed so that the water surface temperature changes according to the temperature transfer coefficient and the equilibrium temperature. For the calibration of the model developed, analytic and numerical results from a case of rectangular channel with full width continuous injection have been compared in a steady state. The comparisons showed that the numerical results were in good agreement with analytical solutions. The application site was selected from the downstream of Paldang dam to Jamsil submerged weir, and overall length of this site is about 22.5 km. The change of water temperature caused by the discharge from the Guri sewage treatment plant has been simulated, and results were similar to the observed data. Overall it is concluded that the developed model can represent the water temperature changes due to heat transport accurately. But the verification using observed data will further enhance the validity of the model.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.20 no.7
• /
• pp.2067-2077
• /
• 1996

In adaptive finite element analysis, r- and h-methods are generally used on the basis of a discretization error estimator. In this paper, an rh-method is proposed as a new adaptive method which can improve the adaptivity performance by using both of them. This suggested rh-method moves nodal coordinates of initially given model to adjust element discretization errors and thereafter performes the h-method tdo obtain the specified accuracy of finite element solutions. Numerical experiments for various plane problems were performed using 4-noded isoparametric quadrilateral elements. As a result, the rh-method has been shown to be an accurate and efficient adaptive analysis method to obtain as improved solution.

• Journal of the Korean Society for Nondestructive Testing
• /
• v.20 no.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.

• Composites Research
• /
• v.25 no.2
• /
• pp.34-39
• /
• 2012

A thermo-poro-elastic constitutive model of unidirectionally fiber-reinforced composite materials is suggested by extending the unmixing-mixing scheme which is based upon composite micromechanics. The strain components of thermal expansion due to a temperature change, gas pressure in pores, and chemical shrinkage are included in the constitutive model. On purpose to verify the derived constitutive relations, the representative volume element of two-dimensional lamina subject to various loading conditions is analyzed by the finite element method. The overall stress and strain responses are obtained, and compared with the predicted values by the unmixing-mixing scheme. The numerical results show the usefulness of the proposed model to predict the thermoelastic behavior of porous composites.