• Proceedings of the Computational Structural Engineering Institute Conference
• /
• 2004.10a
• /
• pp.355-362
• /
• 2004

An efficient boundary-based technique is developed for addressing shape design sensitivity analysis in various problems. An analytical sensitivity formula in the form of a boundary integral is derived based on the continuum formulation for a general functional defined in problems. The formula, which is expressed in terms of the boundary solutions and shape variation vectors, can be conveniently used for gradient computation in a variety of shape design problems. While the sensitivity can be calculated independent of the analysis means, such as the finite element method (FEM) or the boundary element method (BEM), the FEM is used for the analysis in this study because of its popularity and easy-to-use features. The advantage of using a boundary-based method is that the shape variation vectors are needed only on the boundary, not over the whole domain. The boundary shape variation vectors are conveniently computed by using finite perturbations of the shape geometry instead of complex analytical differentiation of the geometry functions. The supercavitating flow problem and fillet problem are chosen to illustrate the efficiency of the proposed methodology. Implementation issues for the sensitivity analysis and optimization procedure are also addressed in these problems.

• Proceedings of the KSME Conference
• /
• 2004.04a
• /
• pp.1047-1052
• /
• 2004

An efficient boundary-based technique is developed for addressing shape design sensitivity analysis in supercavitating flow problem. An analytical sensitivity formula in the form of a boundary integral is derived based on the continuum formulation for a general functional defined in potential flow problems. The formula, which is expressed in terms of the boundary solutions and shape variation vectors, can be conveniently used for gradient computation in a variety of shape design in potential flow problems. While the sensitivity can be calculated independent of the analysis means, such as the finite element method (FEM) or the boundary element method (BEM), the FEM is used for the analysis in this study because of its popularity and easy-touse features. The advantage of using a boundary-based method is that the shape variation vectors are needed only on the boundary, not over the whole domain. The boundary shape variation vectors are conveniently computed by using finite perturbations of the shape geometry instead of complex analytical differentiation of the geometry functions. The supercavitating flow problem is chosen to illustrate the efficiency of the proposed methodology. Implementation issues for and optimization procedure are addressed in this flow problem.

• Proceedings of the Korean Society of Machine Tool Engineers Conference
• /
• 1996.03a
• /
• pp.61-66
• /
• 1996

An overall feature to simulate composite behavior and to predict closed solution has been performed for the application to the stress analysis in a discontinuous composite solid. To obtain the internal field quantities of composite, the micromechanics analysis and finite element analysis (FEA) were implemented. For the numerical illustration, an aligned axisymmetric single fiber model has been employed to assess field quantities. Further, a micromechanics model to describe the elastic behavior of fiber or whisker reinforced metal matrix composites has been developed and the stress concentrations between reinforcements were investigated using the modified shear lag model with the comparions between reinforcements were investigated using the modified shear lag model with the comparison of finite element analysis (FEA). The rationale is based on the replacement of the matrix between fiber ends with the fictitious fiber to maintain the compatibility of displacement and traction. It was found that the new model gives a good agreement with FEA results in the small fiber aspect ratio regime as well as that in the large fiber aspect ratio regime. It was found that the proposed simulation methodology for stress analysis is applicable to the complicated inhomogeneous solid for the investigation of micromechanical behavior.

• Structural Engineering and Mechanics
• /
• v.12 no.5
• /
• pp.541-558
• /
• 2001

The design specifications of guard fences in Japan were reexamined and the revised specifications were implemented from April 1999. Because of the huge consumption in time and cost to test the performances of full-scale guard fences in the field, some assumptions are adopted while modifying the design specifications, and numerical analyses are necessary to confirm the impact performance and safety level of new types of steel highway guard fences. In this study, the finite element models are developed for the heavy trucks and steel highway guard fences to reenact their behaviors, and the solution approach is carried out using nonlinear dynamic analysis software of structures in three dimensions (LS-DYNA). The numerical simulation results are compared with the full-scale on-site testing results to verify the proposed analysis procedure. The collision process is simulated and it is also made possible to visualize the movement of the truck and the performances of guard fences. In addition, the energy shift of the truck kinetic energy to the truck and guard fence Internal energy, and the energy absorption of each guard fence component are studied for the development of a new design methodology of steel highway guard fences based on the energy absorption capacity.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.31 no.6
• /
• pp.381-388
• /
• 2018

This paper presents a generalized finite element formulation for plastic hinge modeling based on lumped plasticity in the classical Euler-Bernoulli beam elements. In this approach, the plastic hinges are effectively modeled using proper enrichment functions describing weak discontinuities of the solution. The proposed methodology enables the insertion of plastic hinges at an arbitrary location without modifying the connectivity of elements. The formations of plastic hinges are instead achieved by hierarchically adding degrees of freedom to existing elements. Convergence analyses such as h- and p-extensions are performed to investigate the effectiveness of the proposed method. The analysis results indicate that the proposed generalized finite element method can achieve theoretical convergence rates for both cases where plastic hinges are located at nodes and within an element, thus demonstrating its accuracy.

• Journal of Mechanical Science and Technology
• /
• v.15 no.4
• /
• pp.422-432
• /
• 2001

There are many analysis formulas for determining the resultant shear force in welds. However, there is no general procedure which is applicable to a joint with all six possible loadings exerted simultaneously. A numerical methodology and computer program for such a problem were developed, and they are capable of analyzing a weld of any shape composed of straight or circular line segments. The computer program developed in this study can also display the design procedures and results using computer graphics. The development of such a design procedure and an interactive computer program for weldments analysis will lead to lower cost.

• Nuclear Engineering and Technology
• /
• v.56 no.5
• /
• pp.1627-1637
• /
• 2024

This research paper proposes a novel methodology for predicting the optimal location of friction supports to effectively mitigate vibrational energy in piping systems. The incorporation of friction forces in the dynamic characteristics of the system introduces inherent nonlinearity, making its analysis challenging. Typically, numerical solutions in the time domain are employed to circumvent the complexities associated with finding analytic solutions for nonlinear systems. However, time domain analysis (TDA) can be computationally intensive and demand significant computational resources due to the intricate calculations stemming from nonlinearity. To address this computational burden, this study presents an efficient approach based on linear analysis to predict the ideal position for installing friction supports as a replacement for fixed supports. Furthermore, we investigate the relationship between the installation positions of friction supports and their effectiveness in absorbing vibrations using the harmonic balanced method (HBM). Both methodologies are validated by comparing the obtained results with those obtained through time domain analysis (TDA) using the finite element method (FEM).

• Nuclear Engineering and Technology
• /
• v.51 no.3
• /
• pp.884-893
• /
• 2019

Spent fuel racks are steel structures used in the storage of the spent fuel removed from the nuclear power reactor. Rack units are submerged in the depths of the spent fuel pool to keep the fuel cool. Their free-standing design isolates their bases from the pool floor reducing structural stresses in case of seismic event. However, these singular features complicate their seismic analysis which involves a transient dynamic response with geometrical nonlinearities and fluid-structure interactions. An accurate estimation of the response is essential to achieve a safe pool layout and a reliable structural design. An analysis methodology based on the hydrodynamic mass concept and implicit integration algorithms was developed ad-hoc, but some dispersion of results still remains. In order to validate the analysis methodology, vibration tests are carried out on a reduced scale mock-up of a 2-rack system. The two rack mockups are submerged in free-standing conditions inside a rigid pool tank loaded with fake fuel assemblies and subjected to accelerations on a unidirectional shaking table. This article compares the experimental data with the numerical outputs of a finite element model built in ANSYS Mechanical. The in-phase motion of both units is highlighted and the water coupling effect is detailed. Results show a good agreement validating the methodology.

• Journal of the Earthquake Engineering Society of Korea
• /
• v.7 no.6
• /
• pp.1-7
• /
• 2003

This paper presents a novel hybrid time-frequency-domain method for nonlinear soil-structure interaction(SSI) analysis. It employs, in a practical manner, a computer code for equivalent linear SSI analysis and a general-purpose nonlinear finite element program. The proposed method first (calculates dynamic responses on a truncated finite element boundary utilizing an equivalent linear SSI program in the frequency domain. Then, a general purpose nonlinear finite element program is employed to analyze the nonlinear SSI problem in the time domain, in which boundary conditions at the truncated boundary are imposed with the responses calculated in the previous frequency domain SSI analysis, In order to validate the proposed method, seismic response analyses are carried out for a 2-D underground subway station in a multi-layered half-space, For the analyses, a equivalent linear SSI code KIESSI-2D is coupled to ANSYS program. The numerical results indicate that the proposed methodology can be a viable solution for nonlinear SSI problems.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.11
• /
• pp.1986-1996
• /
• 2003

Methodology based on the elastic-plastic fracture mechanics has been widely accepted in predicting the critical crack length(CCL) of pressure tubes of CANDU nuclear plants. A conservative estimate of CCL is obtained by employing the J-resistance curves measured with the specimens satisfying plane strain condition as suggested in the ASTM standard. Due to limited thickness of the pressure tubes the curved compact tension(CT) specimens taken out from tile pressure tube have been used in obtaining J-resistance curves. The curved CT specimen inevitably introduce slant fatigue crack during precracking. Hence, effect of specimen geometry and slant crack on J-resistance curve should be explored. In this study, the difference of J integral values between the standard CT specimens satisfying plane strain condition and the nonstandard curved CT with limited thickness (4.2mm) is estimated using finite element analysis. The fracture resistance curves of Zr-2.5Nb obtained previously by other authors are critically discussed. Various finite element analysis were conducted such as 2D analysis under plane stress and plane strain conditions and 3D analysis for flat CT, curved CT with straight crack and curved CT with slant crack front. J-integral values were determined by local contour integration near the crack tip, which was considered as accurate J-values. J value was also determined from the load versus load line displacement curve and the J estimation equation in the ASTM standard. Discrepancies between the two values were shown and suggestion was made for obtaining accurate J values from the load line displacement curves obtained by the curved CT specimens.