• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.15 no.3
• /
• pp.491-499
• /
• 2002

A p-version finite element model based on degenerate shell element is proposed tot the analysis of orthotropic laminated plates. In the nonlinear formulation of the model, the total Lagrangian formulation is adopted with large deflection and moderate rotation being accounted tot in the sense of yon Karman hypothesis. The material model is based on the Huber-Mises yield criterion and Prandtl-Reuss flow rule in accordance with the theory of strain hardening yield function, which is generalized lot anisotropic materials by introducing the parameters of anisotropy. The model is also based on extension of equivalent-single layer laminate theory(ESL theory) with shear deformation, leading to continuous shear strain at the interface of two layers. The integrals of Legendre polynomials are used for shape functions with p-level varying from 1 to 10. Gauss-Lobatto numerical quadrature is used to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed P-version finite element model is demonstrated through several comparative points of iew in terms of ultimate load, convergence characteristics, nonlinear effect, and shape of plastic tone.

• Journal of the Microelectronics and Packaging Society
• /
• v.17 no.1
• /
• pp.69-73
• /
• 2010

Finite-element analyses were conducted to investigate the thermal stress in 3-dimensional stacked wafers package containing through-silicon-via (TSV), which is being widely used for 3-Dimensional integration. With finite element method (FEM), thermal stress was analyzed with the variation of TSV diameter, bonding diameter, pitch and TSV height. It was revealed that the maximum von Mises stresses occurred at the edge of top interface between Cu TSV and Si and the Si to Si bonding site. As TSV diameter increased, the von Mises stress at the edge of TSV increased. As bonding diameter increased, the von Mises stress at Si to Si bonding site increased. As pitch increased, the von Mises stress at Si to Si bonding site increased. The TSV height did not affect the von Mises stress. Therefore, it is expected that smaller Cu TSV diameter and pitch will ensure mechanical reliability because of the smaller chance of plastic deformation and crack initiation.

• Tunnel and Underground Space
• /
• v.29 no.3
• /
• pp.197-213
• /
• 2019

We simulated the fault reactivation experiment conducted at 'Main Fault' intersecting the low permeability clay formations of Mont Terri Underground Research Laboratory in Switzerland using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. We formulate the hydro-mechanical coupling relation of hydraulic aperture to consider the elastic fracture opening and failure-induced dilation for reproducing the abrupt changes in injection flow rate and monitoring pressure at fracture opening pressure. A parametric study was conducted to examine the effects of in-situ stress condition and fault deformation and strength parameters and to find the optimal parameter set to reproduce the field observations. In the best matching simulation, the fracture opening pressure and variations of injection flow rate and monitoring pressure showed good agreement with field experiment results, which suggests the capability of the numerical model to reasonably capture the fracture opening and propagation process. The model overestimated the fault displacement in shear direction and the range of reactivated zone, which was attributed to the progressive shear failures along the fault at high injection pressure. In the field experiment results, however, fracture tensile opening seems the dominant mechanism affecting the hydraulic aperture increase.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.27 no.1
• /
• pp.25-35
• /
• 2023

Considering the non-linear behavior of structure and soil when evaluating a nuclear power plant's seismic safety under a beyond-design basis earthquake is essential. In order to obtain the nonlinear response of a nuclear power plant structure, a time-domain SSI analysis method that considers the nonlinearity of soil and structure and the nonlinear Soil-Structure Interaction (SSI) effect is necessary. The Boundary Reaction Method (BRM) is a time-domain SSI analysis method. The BRM can be applied effectively with a Perfectly Matched Layer (PML), which is an effective energy absorbing boundary condition. The BRM has a characteristic that the magnitude of the response in far-field soil increases as the boundary interface of the effective seismic load moves outward. In addition, the PML has poor absorption performance of low-frequency waves. For this reason, the accuracy of the low-frequency response may be degraded when analyzing the combination of the BRM and the PML. In this study, the accuracy of the analysis response was improved by adjusting the PML input parameters to improve this problem. The accuracy of the response was evaluated by using the analysis response using KIESSI-3D, a frequency domain SSI analysis program, as a reference solution. As a result of the analysis applying the optimal PML parameter, the average error rate of the acceleration response spectrum for 9 degrees of freedom of the structure was 3.40%, which was highly similar to the reference result. In addition, time-domain nonlinear SSI analysis was performed with the soil's nonlinearity to show this study's applicability. As a result of nonlinear SSI analysis, plastic deformation was concentrated in the soil around the foundation. The analysis results found that the analysis method combining BRM and PML can be effectively applied to the seismic response analysis of nuclear power plant structures.

• Progress in Medical Physics
• /
• v.23 no.4
• /
• pp.269-278
• /
• 2012

Aquaplast Thermoplastic (AT) is a tissue-equivalent oral compensator that has been developed to improve dose uniformity at the common boundary and around the treated area during radiotherapy in patients with head and neck cancer. In order to assess the usefulness of AT, the degree of improvement in dose distribution and physical properties were compared to those of oral compensators made using paraffin, alginate, and putty, which are materials conventionally used in dental imprinting. To assess the physical properties, strength evaluations (compression and drop evaluations) and natural deformation evaluations (volume change over time) were performed; a Gafchromic EBT2 film and a glass dosimeter inserted into a developed phantom for dose verification were used to measure the common boundary dose and the beam profile to assess the dose delivery. When the natural deformation of the oral compensators was assessed over a two-month period, alginate exhibited a maximum of 80% change in volume from moisture evaporation, while the remaining tissue-equivalent properties, including those of AT, showed a change in volume that was less than 3%. In a free-fall test at a height of 1.5 m (repeated 5 times as a strength evaluation), paraffin was easily damaged by the impact, but AT exhibited no damage from the fall. In compressive strength testing, AT was not destroyed even at 8 times the force needed for paraffin. In dose verification using a glass dosimeter, the results showed that in a single test, the tissue-equivalent (about 80 Hounsfield Units [HU]) AT delivered about 4.9% lower surface dose in terms of delivery of an output coefficient (monitor unit), which was 4% lower than putty and exhibited a value of about 1,000 HU or higher during a dose delivery of the same formulation. In addition, when the incident direction of the beam was used as a reference, the uniformity of the dose, as assessed from the beam profile at the boundary after passing through the oral compensators, was 11.41, 3.98, and 4.30 for air, AT, and putty, respectively. The AT oral compensator had a higher strength and lower probability of material transformation than the oral compensators conventionally used as a tissue-equivalent material, and a uniform dose distribution was successfully formed at the boundary and surrounding area including the mouth. It was also possible to deliver a uniformly formulated dose and reduce the skin dose delivery.

• Tunnel and Underground Space
• /
• v.28 no.5
• /
• pp.400-425
• /
• 2018

This study presents the research results and current status of the DECOVALEX-2019 project Task B. Task B named 'Fault slip modelling' is aiming at developing a numerical method to simulate the coupled hydro-mechanical behavior of fault, including slip or reactivation, induced by water injection. The first research step of Task B is a benchmark simulation which is designed for the modelling teams to familiarize themselves with the problem and to set up their own codes to reproduce the hydro-mechanical coupling between the fault hydraulic transmissivity and the mechanically-induced displacement. We reproduced the coupled hydro-mechanical process of fault slip using TOUGH-FLAC simulator. The fluid flow along a fault was modelled with solid elements and governed by Darcy's law with the cubic law in TOUGH2, whereas the mechanical behavior of a single fault was represented by creating interface elements between two separating rock blocks in FLAC3D. A methodology to formulate the hydro-mechanical coupling relations of two different hydraulic aperture models and link the solid element of TOUGH2 and the interface element of FLAC3D was suggested. In addition, we developed a coupling module to update the changes in geometric features (mesh) and hydrological properties of fault caused by water injection at every calculation step for TOUGH-FLAC simulator. Then, the transient responses of the fault, including elastic deformation, reactivation, progressive evolutions of pathway, pressure distribution and water injection rate, to stepwise pressurization were examined during the simulations. The results of the simulations suggest that the developed model can provide a reasonable prediction of the hydro-mechanical behavior related to fault reactivation. The numerical model will be enhanced by continuing collaboration and interaction with other research teams of DECOLVAEX-2019 Task B and validated using the field data from fault activation experiments in a further study.