• Journal of the Korean Geotechnical Society
• /
• v.27 no.5
• /
• pp.45-59
• /
• 2011

This study introduces the diagrammatic Upper and Lower Bound (UB and LB) methods theoretically in order to derive the bearing capacity factor, $N_c$ in undrained clay and to compare with Prandtl's exact solution (1921). As a result of the theoretical study, an exact solution comes out when the UB and LB solutions are the same. In addition, the finite element analyses show that the failure loads approach to the bearing capacity factor of 5.14. Results of the FEA significantly depend on the finite element type, a number of elements, and a number of increments. From this study the exact solution defines that solutions from UB and LB are the same. However, this situation is very difficult to process, so we can confirm the exact solution as a range between UB and LB solutions.

• KSCE Journal of Civil and Environmental Engineering Research
• /
• v.4 no.2
• /
• pp.55-64
• /
• 1984

Presented is a nonlinear constitutive model for progressive tensile fracturing of concrete. The model is incremental, path-dependent, and tensorialy invariant. The total strain tensor is assumed to be a sum of a purely elastic component and an inelastic component. The material is considered to contain weak planes of all directions which characterize the planes of the microcracks. A one-to-one functional dependence is assumed between the normal stress and the normal strain across each of the weak planes. The tangential stiffness of concrete is then derived form the principle of virtual work. The present theory can be applied to loading histories which are nonproportional or during which the principal directions rotate. Good agreement with the available direct tensile test data which cover strain-softening is demonstrated.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.17 no.4
• /
• pp.375-387
• /
• 2004

A new finite element model will be presented to analyze the nonlinear behavior of RC beams and slabs strengthened by a patch repair. The numerical approach is based on the p-version degenerate shell element including theory of anisotropic laminated composites, theory of materially and geometrically nonlinear plates. In the nonlinear formulation of this model, the total Lagrangian formulation is adopted with large deflections and moderate rotations being accounted for in the sense of von Karman hypothesis. The material model is based on hardening rule, crushing condition, plate-end debonding strength model and so on. The Gauss-Lobatto numerical quadrature is applied to calculate the stresses at the nodal points instead of Gauss points. The validity of the proposed p-version nonlinear finite element model is demonstrated through the load-deflection curves, the ultimate loads, and the failure modes of RC beams or slabs bonded with steel plates or FRP plates compared with available result of experiment and other numerical methods.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.39 no.7
• /
• pp.693-700
• /
• 2015

To estimate the fatigue crack propagation behavior of compact tension shear (CTS) specimen under mixed-mode loads, crack path prediction theories and Tanaka's equation were applied. The stress intensity factor at a newly created crack tip was calculated using a finite element method via ANSYS, and the crack path and crack increment were then obtained from the crack path prediction theories, Tanaka's equation, and the Paris' equation, which were preprogrammed in Microsoft Excel. A new method called the finite element crack tip updating method (FECTUM) was developed. In this method, the finite element method and Microsoft Excel are used to calculate the stress intensity factors and the crack path, respectively, at the crack tip per each crack increment. The developed FECTUM was applied to simulate the fatigue crack propagation of a single-edge notched bending (SENB) specimen under eccentric three-point bending loads. The results showed that the number of cycles to failure of the specimen obtained experimentally and numerically were in good agreement within an error range of less than 3%.

• Journal of the Korean Geotechnical Society
• /
• v.31 no.1
• /
• pp.5-14
• /
• 2015

Seabed may undergo large excess pore water pressure in the case of long duration of high wave loading. This excess pore water pressure may reduce effective stress and, consequently, the seabed may liquefy. Thus, it is necessary to develop a numerical technique which can precisely evaluate the dynamic response of seabed due to wave action. In this study, a new numerical technique named mixed model (2D NIT & FLIP models) was proposed. The dynamic wave pressure and water flow velocity acting on the boundary between seabed and the wave field was estimated using 2D-NIT model. This result was used as input data in FLIP program for investigation of dynamic response of seabed. To secure the reliability of the mixed model, the numerical analysis results of the mixed model were compared with Yamamoto's solution and Chang's experiment results. The comparison results indicated that there were some differences between them, but the general trend of the effective stress increment and the excess pore water pressure along the depth of seabed was similar to each other. Thus, this study clearly supports the plausibility of the numerical analysis of the mixed model.

• Journal of Korean Association for Spatial Structures
• /
• v.11 no.1
• /
• pp.121-130
• /
• 2011

Space frame structures have the advantage of constructing a large space structures without column and it may be considered as a shell structure. Nevertheless, with the characteristics of thin and long term of spacing, the unstable problem of space structure could not be set up clearly, and there is a huge difference between theory and experiment. Therefore, in this work, the tangential stiffness matrix of space frame structures is studied to solve the instability problem, and the nonlinear incremental analysis of the structures considering rise-span ratio(${\mu}$) and the ratio of load($R_L$) is performed for searching unstable points. Basing on the results of the example, global buckling can be happened by low rise-span ratio(${\mu}$), nodal buckling can be occurred by high rise-span ratio(${\mu}$). And in case of multi node space structure applying the ratio of load($R_L$), the nodal buckling phenomenon occur at low the ratio of load($R_L$), the global buckling occur a1 high the ratio of load($R_L$). In case of the global buckling, the load of bifurcation is about from 50% to 70% of perfect one's snap-through load.

• Transactions of the Korean Society of Mechanical Engineers A
• /
• v.35 no.9
• /
• pp.1041-1050
• /
• 2011

An elasto-plastic finite element method using the theory of strain gradient plasticity is proposed to evaluate the size dependency of structural plasticity that occurs when the configuration size decreases to micron scale. For this method, we suggest a low-order plane and three-dimensional displacement-based elements, eliminating the need for a high order, many degrees of freedom, a mixed element, or super elements, which have been considered necessary in previous researches. The proposed method can be performed in the framework of nonlinear incremental analysis in which plastic strains are calculated and averaged at nodes. These strains are then interpolated and differentiated for gradient calculation. We adopted a strain-gradient-hardening constitutive equation from the Taylor dislocation model, which requires the plastic strain gradient. The developed finite elements are tested numerically on the basis of typical size-effect problems such as micro-bending, micro-torsion, and micro-voids. With respect to the strain gradient plasticity, i.e., the size effects, the results obtained by using the proposed method, which are simple in their calculation, are in good agreement with the experimental results cited in previously published papers.

• Journal of the Korean Society for Aeronautical & Space Sciences
• /
• v.35 no.10
• /
• pp.891-898
• /
• 2007

Optimal aerodynamic design for the pitch-controlled horizontal axis wind turbine and its aerodynamic performance for various pitch angles are performed numerically by using the blade element momentum theory. The numerical calculation includes effects such as Prandtl‘s tip loss, airfoil distribution, and wake rotation. Six different airfoils are distributed along the blade span, and the special airfoil i.e. airfoil of 40% thickness ratio is adopted at the hub side to have structural integrity. The nonlinear chord obtained from the optimal design procedure is linearized to decrease the weight and to increase the productivity with very little change of the aerodynamic performance. From the comparisons of the power, thrust, and torque coefficients with corresponding values of different pitch angles, the aerodynamic performance shows delicate changes for just $3^{\circ}$ increase or decrease of the pitch angle. For precisive pitch control, it requires the pitch control algorithm and its drive mechanism below $3^{\circ}$ increment of pitch angle. The maximum torque is generated when the speed ratio is smaller than the designed one.

• Magazine of the Korea Concrete Institute
• /
• v.8 no.3
• /
• pp.187-195
• /
• 1996

A finite element formulation based on the CFT(Compression Field Theory), considering the effect of compression softening in cracked concrete, and macro-scopic and rotating crack models etc., was presented for the nonlinear behaviour of structural concrete. Considering the computational efficency and the ability of modelling the post-ultimate behaviour as major concerns, the Incremental displacement solution algorithm involving initial material stiffnesses and the relaxation procedure for fast convergence was adopted and formulated in a type of 8-noded quadrilateral isoparametric elements. The analysis program NASCOM(Non1inear Analysis of Structural Concrete by FEM : Monotonic Loading) developed in this way enables the predictions of strength and deformation capacities in a full range, crack patterns and their corresponding widths, and yield extents of reinforcement. As the verification purpose of NASCOM, the predictions were made for Bhide's Panel(PB21) and Leonhardt's deep beam tests. The predicted results shows somewhat stiff behaviour for the panel test, and vice versa for deep beam tests. More refining process would be necessary hereafter in terms of more accurately simulating the effects of tension-stiffening and compression softening in concrete.

• Journal of the Computational Structural Engineering Institute of Korea
• /
• v.22 no.5
• /
• pp.487-499
• /
• 2009

The postprocessing technology has been advanced diversely to accommodate the tendency of increasingly refined and complicated practices of finite element modeling in pace with enhanced capabilities of computers and improved algorithm of equation solvers. As a result of such progresses in both hardware and software, it became practically meaningful to inspect and analyze the elasto-plastic behavior using the intermediate results from the increasing number of incremental and iterative processes. This paper is concerned about the new methods of postprocessing with computer graphic visualization of elasto-plastic behavior on the basis of the theoretically reorganized analysis results. This paper proposes a new method of rendering the plastic zone, and new approaches of analyzing and interpreting the elasto-plastic behavior using the graphical information visualized in the form of the yield surface and the stress path, or in the form of the Mohr circles and the failure envelope.