• Computational Structural Engineering
• /
• v.8 no.1
• /
• pp.173-186
• /
• 1995

It is well known that structures constructed by bonding two or more materials and then subjected to temperature change experience thermal stress. This stress results from thermal expansion mismatch of materials. The present paper derives formulas for the stresses in a bimaterial axisymmetric disk which is subjected to a uniform temperature change. First, an approximate solution following strength-of-materials principles is developed. However, the strength-of-materials solution has difficulty in predicting both the peak value of interfacial stresses and its associated distribution. Next, a solution consistent with the theory of elasticity is developed by way of an eigenfunction expansion approach. The eigenfunction analysis is compared with finite element stress analysis results for a specific numerical example. Finite element analysis results show that the interfacial stresses are adequately predicted by eigenfunction solution. Therefore, the method developed in this paper will be useful in determination of the interfacial stress state.

• Geotechnical Engineering
• /
• v.14 no.4
• /
• pp.129-140
• /
• 1998

Elasto-plastic finite element analysis of geotechnical boundary value problems necessitate the stress integration for the known strain increments. For the elasto-plastic constitutive model, the stress integration is generally achieved by numerical schemes, because analytical integration is impossible for general strain path. In this case, the accuracy of numerical stress integration has an important role on the overall accuracy of nonlinear finite element solution. In this study, the Sloan's substepping method which is one of explicit integration methods has been adopted and iris applicability has been checked. The unstability and inaccuracy of ifs results initiated from initial stress level were revealed. So. a new modified numerical integration method which employs the basic concept of modified Euler scheme for error control is proposed and accuracy and stability of the solutions are confirmed by triaxial test simulation.

• Composites Research
• /
• v.29 no.6
• /
• pp.388-394
• /
• 2016

Fiber reinforced plastic (FRP) composites have drawn increasing attentions worldwide for decades due to its outstanding properties. Stress concentration factor (SCF) as an essential parameter in materials science are critically considered in structure design and application, strength assessment and failure prediction. However, investigation of stress concentration in FRP composites has been rarely reported so far. In this study, three resource-conserving analyses (Isotropic analysis, Orthotropic analysis and Finite element analysis) were introduced to plot the $K_T^A-d/W$ curve for E-glass/epoxy composite plate with the geometrical defect of circular hole placed centrally. The plates were loaded to uniaxial direction for simplification. Finite element analysis (FEA) was carried out via ACP (ANSYS composite prepost module). Based on the least squares method, a simple expression of fitting equation could be given based on the simulated results of a set of discrete points. Finally, all three achievable solutions were presented graphically for explicit comparison. In addition, the investigation into customized efficient SCFs has also been carried out for further reference.

• Journal of the Korean Society of Manufacturing Process Engineers
• /
• v.13 no.2
• /
• pp.73-80
• /
• 2014

A camera module that gathers visual information via aerial observation reconnaissance is equipped inside a gimbal structure. This gimbal structure system must reduce dynamic responses in order to obtain clear images under all circumstances. Among many design specifications for this system, there is MIL-STD-810G as a shock standard. This specification indicates a limitation of the acceleration of the camera module under a base shock excitation on the gimbal structure. The satisfaction of this condition can usually be proved by experiment, because it includes bearings and dynamic isolators made of rubber. Numerical analysis must be proposed for design improvement of the gimbal structure. To achieve this goal, transient response analysis for the base shock excitation was performed using the finite element method. Experimental results were compared with numerical solutions and it is shown that the present method is useful.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.33 no.3
• /
• pp.865-873
• /
• 2013

This paper proposes a theoretical solution of elastic critical buckling load of infinitely long pipelines with non-uniform thickness under external pressure. The non-uniform cross section of pipelines can be considered as corroded or stiffened pipelines so that this paper can be a fundamental research of pipelines that are essential technology for offshore industries. The theoretical solution of pipelines with non-uniform thickness is derived with an assumption that a cylindrical shell under external pressure can be considered as a simple ring. The eigenfunctions are derived to obtain the critical buckling load. The reduced thickness and the reduced range are considered as variables in parametric analysis. The finite element analysis is performed to verify the theoretical solutions and the results of the analytic method and the finite element method are in good agreement.

• Proceedings of the KSME Conference
• /
• 2000.04a
• /
• pp.596-601
• /
• 2000

The paper presents free vibration and stability analyses of a non-uniform Timoshenko beam resting on a two-parameter elastic soil. The soil parameters can vary along the spat and is assumed to be two-parameter model including the effects of both transverse shear deformation and elastic foundation Governing equations related to the vibration and the stability of the beam are derived from Hamilton's principle, and the resulting eigen-value problems can be solved to give natural frequencies and critical force by finite element method. Numerical results for both vibration and stability of beams under an axial force are presented and compared with other available solutions. Finally, vibration frequencies, mode shapes and critical forces are investigated for various thickness ratios, shear foundation parameter, Winkler foundation parameter and boundary conditions of tapered Timoshenko beams.

• The Journal of the Korea Contents Association
• /
• v.13 no.2
• /
• pp.505-511
• /
• 2013

A new boundary treatment technique which can be applied to axi-symmetric topography with inclined bottom was developed. Although the finite element method is good for complex geometry, there is no proper boundary treatment when a boundary is not a vertical section because the water depth at the coastline becomes zero. In this study, we developed a new boundary treatment for inclined bottom using the analytical solution for long wave. To develope a model, the mild-slope equation was used and then, a computational domain is divided into an analytical region and a numerical region. By combining a numerical and an analytical solutions, a complete solution was obtained. The developed solution was validated by comparing with a previous analytical solution.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.27 no.4
• /
• pp.551-558
• /
• 2003

An efficient and reliable subspace iteration algorithm using the block algorithm is proposed. The block algorithm is the method dividing eigenpairs into several blocks when a lot of eigenpairs are required. One of the key for the faster convergence is carefully selected initial vectors. As the initial vectors, the proposed method uses the modified Ritz vectors for guaranteering all the required eigenpairs and the quasi-static Ritz vectors for accelerating convergency of high frequency eigenvectors. Applying the quasi-static Ritz vectors, a shift is always required, and the proper shift based on the geometric average is proposed. To maximize efficiency, this paper estimates the proper number of blocks based on the theoretical amount of calculation in the subspace iteration. And it also considers the problems generated in the process of combining various algorithms and the solutions to the problems. Several numerical experiments show that the proposed subspace iteration algorithm is very efficient, reliable ,and accurate.

• Journal of Korea Water Resources Association
• /
• v.30 no.1
• /
• pp.35-44
• /
• 1997

Petrov-Galerkin finite element model for analyzing dynamic wave equation is applied to gradually and rapidly varied unsteady flow. The model in verified by applying to hydraulic jump, nonlinear disturbance propagation in frictionless horizontal channel and dam-break analysis. It shows stable and accurate results compared with analytical solutions for various cases. The model in applied to a surge propagation in a frictionless horizontal channel. Three-dimensional water surface profiles show that the computed result converges to the analytical one with sharp discontinuity. The model is also applied to the Taehaw River to analyze unsteady floodwave propagation. The computed results have good agreements with those of DWOPER model in terms of discharge and stage hydrographs.

• Journal of Electrical Engineering and Technology
• /
• v.5 no.4
• /
• pp.631-636
• /
• 2010

Gas-insulated transmission lines (GITL) are valued as technological solutions in hydro-power stations due to their enormous power handling capabilities. The performance of GITL is a function of the size of metallic particles inside the gas-insulated chamber. Electrostatic field (E-field) enhancement is a common phenomenon in gas-insulated lines due to these metallic particles. In this study, the E-field enhancement factor is calculated by considering metallic particles at various locations in the gas-insulated line/bus section, such as high-tension (HT) conductor, high-voltage shields, support insulator, and inner surface of grounded enclosure. For this purpose, a two-dimensional model based on finite element (FE) method is developed. The length of the metallic particle is in the range of 1 to 10 mm while the diameter is between 1 to 3 mm. E-field enhancement is also computed for various particle configurations of the gas-insulated system, with focus on dielectric coating made of epoxy on HT conductor and inner surface of grounded enclosure.