• Journal of the Korean Geotechnical Society
• /
• v.17 no.4
• /
• pp.7-14
• /
• 2001

구조물 하중에 의한 점성토 지반의 침하량을 보다 정확하게 평가하기 위해서는 지반 내의 흙요소가 경험하는 실제적인 응력경로와 이에 따른 변형양상이 적절하게 고려되어야만 한다. 따라서 본 연구에서는 축대칭 조건의 다양한 응력경로를 따라 발생하는 정규압밀 점성토의 변형 거동을 고찰한 기존의 실험적 연구결과를 바탕으로 응력경로법에 근거한 보다 간편하고 합리적인 침하량 평가기법을 제시하였다. 또한 본 연구에서는 제시된 평가기법을 기존의 1차원적인 침하량 평가기법들과 함께 실제와 유사한 조건을 가지는 가상지반의 침하량 산정에 적용해 보았으며, 동일한 조건에 대해 소성모델(MCC 모델)과 혼합압밀이론에 바탕을 둔 유한요소해석을 실시하였다. 그리고 이를 통해 얻어진 결과들을 비교.분석함으로써 기존 평가기법들의 문제점과 한계를 명확히 제시하였으며, 응력증분 평가방법이 침하량 평가에 미치는 영향을 분석하였다.

• Composites Research
• /
• v.29 no.1
• /
• pp.16-23
• /
• 2016

In this paper, numerical results of sectional analysis, stress recovery and energy release rate were compared with the results of VABS, 3-D FEM through the blade analysis library. The result of stress recovery analysis for one-dimensional model including the stiffness matrix is compared with stress results of three-dimensional FEM. We discuss the configuration of the blade analysis library and compare verifications of numerical analysis results of VABS. Blade analysis library through dimensional reduction and stress recovery is intended to be utilized in conjunction with pre- and post-processing of the analysis program of the composite blade, high-altitude uav's wing, wind blades and tilt rotor blade.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.38 no.7
• /
• pp.727-734
• /
• 2014

This paper describes multi-scale based ductile fracture simulation using finite element (FE) damage analysis. The maximum and crack initiation loads of cracked components were predicted using proposed virtual testing method. To apply the local approach criteria for ductile fracture, stress-modified fracture strain model was adopted as the damage criteria with modified calibration technique that only requires tensile and fracture toughness test data. Element-size-dependent critical damage model is also introduced to apply the proposed ductile fracture simulation to large-scale components. The results of the simulation were compared with those of the tests on SA333 Gr. 6 full-scale pipes at $288^{\circ}C$, performed by the Battelle Memorial Institute.

• Computational Structural Engineering
• /
• v.10 no.1
• /
• pp.201-211
• /
• 1997

A clearly consistent finite element formulation for geometrically non-linear analysis of space frames is presented by applying incremental equilibrium equations based on the updated Lagrangian formulation and introducing Vlasov's assumption. The improved displacement field for symmetric cross sections is introduced based on inclusion of second order terms of finite rotations, and the potential energy corresponding to the semitangential rotations and moments is consistently derived. For finite element analysis, elastic and geometric stiffness matrices of the space frame element are derived by using the Hermitian polynomials as shape functions. A co-rotational formulation in order to evaluate the unbalanced loads is presented by separating the rigid body rotations and pure deformations from incremental displacements and evaluating the updated direction cosines of the frame element due to rigid body rotations and incremental member forces from pure deformaions. Finite element solutions for the spatial buckling and post-buckling analysis of space frames are compared with available solutions and other researcher's results.

• Computational Structural Engineering
• /
• v.10 no.4
• /
• pp.205-216
• /
• 1997

In order to trace the lateral post-buckling behaviors of thin-wafled space frames, a geometrically nonlinear finite element formulation is presented by applying incremental equilibrium equations based on the updated Lagrangian formulation and introducing Vlasov's assumption. The improved displacement field for symmetric thin-walled cross sections is introduced based on inclusion of second order terms of finite rotations, and the potential energy corresponding to the semitangential rotations and moments is consistently derived. For finite element analysis, tangent stiffness matrices of the thinwalled space frame element with 7 degrees of freedom including the restrained warping for each node are derived by using the Hermition polynomials as shape functions. A co-rotational formulation in order to evaluate the unbalanced loads is presented by separating the rigid body rotations and pure deformations from incremental displacements and evaluating the updated direction cosines of the frame element due to rigid body rotations and incremental member forces from pure deformations. Finite element solutions for the spatial buckling and post-buckling analysis of thin-walled space frames are presented and compared with available solutions and other researcher's results.

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

The mixed mode problem (I and II) of a peny-shaped interface cracks in remote tension loading on a bi-material cylinder is studied using finite element method. The energy release rates for the tip of the crack in the interface were calibrated for several different moduli combinations and crack ratios using the modified crack closure integral technique and J-integral method, with numerical results obtained from a commercial finite element program. Numerical results show that non-dimensional value of$\sqrt{G_{II}E^*}/\sqrt[p]{\pi a}$ increases as the crack size or moduli ratio increases. Meanwhile, non-dimensional value of$\sqrt{G_{I}E^*}/\sqrt[p]{\pi a}$ decreases as the moduli ratio increases, but above the moduli ratio of 3 its value decreases then increases again as the crack size increases. Reliability of the numerical analysis in this study was acquired with comparison to an analytical solution for the peny-shaped interface crack in an infinite medium.

• Computational Structural Engineering
• /
• v.6 no.4
• /
• pp.57-66
• /
• 1993

The p-version crack model based on integrals of Legendre polynomial and virtual crack extension method is proposed with its potential for application to stress intensity factor computations in linear elastic fracture mechanics. The main advantage of this model is that the data preparation effort is minimal because only a small number of elements are used and high accuracy and the rapid convergence can be achieved in the vicinity of crack tip. There are two important findings from this study. Firstly, the limit value, the strain energy of the exact solution, can be estimated with successive three p-version approximations by ascertaining that the approximations enter the asymptotic range. Secondly, the rate of convergence of p-version model is almost twice that of h-version model on the basis of uniform or quasiuniform mesh refinement for the cracked panel problem subjected to tension.

• Journal of the Korean Society of Safety
• /
• v.11 no.4
• /
• pp.143-150
• /
• 1996

This is an explicit-Implicit, finite element analysis for linear as well as nonlinear hygrothermal stress problems. Additional features, such as moisture diffusion equation, crack element and virtual crack extension(VCE ) method for evaluating J-integral are implemented in this program. The Linear Elastic Fracture Mechanics(LEFM) Theory is employed to estimate the crack driving force under the transient condition for and existing crack. Pores in materials are assumed to be saturated with moisture in the liquid form at the room temperature, which may vaporize as the temperature increases. The vaporization effects on the crack driving force are also studied. The Ideal gas equation is employed to estimate the thermodynamic pressure due to vaporization at each time step after solving basic nodal values. A set of field equations governing the time dependent response of porous media are derived from balance laws based on the mixture theory Darcy's law Is assumed for the fluid flow through the porous media. Perzyna's viscoplastic model incorporating the Von-Mises yield criterion are implemented. The Green-Naghdi stress rate is used for the invariant of stress tensor under superposed rigid body motion. Isotropic elements are used for the spatial discretization and an iterative scheme based on the full newton-Raphson method is used for solving the nonlinear governing equations.

• Journal of the Korean Geotechnical Society
• /
• v.16 no.1
• /
• pp.51-64
• /
• 2000

This paper presents a procedure of calculating a safety factor of the unsaturated slope suffering from the rainfall infiltration. The process of infiltration into a slope due to rainfall and its effect on the behavior of the soil slope are examined by using a two dimensional finite element flow-deformation coupled analysis. A factor of safety is calculated at various elapsed times after the commencement of rainfall as in the following procedure. First, stresses are estimated at each Gaussian point from the coupled finite element analysis. Then, the global stress smoothing method is applied to get a continuous stress field. Based on this stress field, a factor of safety is calculated for a specified slip surface by a stress integration scheme. Then, a search strategy is used to find out a critical slip surface which is associated with the minimum factor of safety. Some numerical examples are analyzed in order to study the effect of hydraulic conductivity on the slope stability during rain-induced infiltration. According to the results, local failure zone can be formed near the slope surface due to inhomogeneous distribution of hydraulic conductivity If the failure zone is once formed, then the region extends until a large amount of slide activates. Therefore the local failure can be neglected no longer in the stability analysis.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.13 no.2
• /
• pp.245-256
• /
• 2000

A moment-curvature relationship to simulate the behavior of reinforced concrete beam under cyclic loading is introduced. Unlike previous moment-curvature models and the layered section approach, the proposed model takes into consideration the bond-slip effect by using monotonic moment-curvature relationship constructed on the basis of the bond-slip relation and corresponding equilibrium equation at each nodal point. In addition, the use of curved unloading and reloading branches inferred from the stress-strain relation of steel gives more exact numerical result. The advantages of the proposed model, comparing to layered section approach, may be on the reduction in calculation time and memory space in case of its application to large structures. The modification of the moment-curvature relation to reflect the fixed-end rotation and pinching effect is also introduced. Finally, correlation studies between analytical results and experimental studies are conducted to establish the validity of the proposed model.